CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General Forums > Main CFD Forum

three coupled equations

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 20, 2016, 14:32
Red face three coupled equations
  #1
New Member
 
Join Date: Dec 2015
Posts: 6
Rep Power: 10
ptub is on a distinguished road
Hi,

I am trying to solve three coupled equations in 1D computational domain. For example, if I want to write my first equations, there are six unknowns a(i+1), a(i), a(i-1) ,c(i-1) c(i),c(i+1). and so for second and third equation. However, I couldn't get tri-diagonal matrix I didn't get how to set my unknowns into matrix. Can somebody help? Thanks.
ptub is offline   Reply With Quote

Old   February 20, 2016, 14:59
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,760
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
But it is a problem for you if the matrix is not tri-diagonal?
your unknown vector x is something like this?

A.x=q

x=
[a(1),...a(N), b(1), ....b(N), c(1),....c(N)]^T
FMDenaro is offline   Reply With Quote

Old   February 20, 2016, 15:04
Default
  #3
New Member
 
Join Date: Dec 2015
Posts: 6
Rep Power: 10
ptub is on a distinguished road
Yeap. My first equation for a(1) includes unknowns C(1) and C(2). Also b(1) includes C(1) and C(2). C(1) includes A(1) and B(1). At the end, I didn't get tri-diagonal matrix. I had tri-diagonal matrix solver. I don't have sparse matrix solver. I am not sure how to proceed.
ptub is offline   Reply With Quote

Old   February 20, 2016, 15:29
Default
  #4
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,760
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by ptub View Post
Yeap. My first equation for a(1) includes unknowns C(1) and C(2). Also b(1) includes C(1) and C(2). C(1) includes A(1) and B(1). At the end, I didn't get tri-diagonal matrix. I had tri-diagonal matrix solver. I don't have sparse matrix solver. I am not sure how to proceed.
Have you tried to organize a block-matrix form?
FMDenaro is offline   Reply With Quote

Old   February 20, 2016, 15:38
Default
  #5
New Member
 
Join Date: Dec 2015
Posts: 6
Rep Power: 10
ptub is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
But it is a problem for you if the matrix is not tri-diagonal?
your unknown vector x is something like this?

A.x=q

x=
[a(1),...a(N), b(1), ....b(N), c(1),....c(N)]^T
I am not sure how to do that. If I do that, my unknown matrix won't be like this x.
ptub is offline   Reply With Quote

Old   February 21, 2016, 13:33
Default
  #6
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,760
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Check for an interative splitting method so that you can work with iterated Thomas solver
FMDenaro is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
coupled source term in 3 transport equations mhassani OpenFOAM Running, Solving & CFD 1 September 10, 2018 10:35
Coupled Solution of the RANS turbulence equations sfigato OpenFOAM Programming & Development 0 September 14, 2014 09:59
Coupled Equations ? T.D. OpenFOAM 0 May 4, 2011 07:28
problem with gradient as source term in coupled equations Cyp OpenFOAM Programming & Development 0 August 16, 2010 14:12
coupled Navier-Stokes and convection-diffusion equations using in cfx rystokes CFX 0 November 2, 2009 06:14


All times are GMT -4. The time now is 20:07.