CFD Online Logo CFD Online URL
Home > Forums > General Forums > Main CFD Forum

Differences in 1D and 2D flow results using Gauss-Seidal/Jacobi

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

LinkBack Thread Tools Search this Thread Display Modes
Old   November 19, 2016, 17:21
Default Differences in 1D and 2D flow results using Gauss-Seidal/Jacobi
New Member
Benjamin Harris
Join Date: Nov 2015
Posts: 2
Rep Power: 0
u05harrisb is on a distinguished road
Dear All,

Thank you in advance for reading and contributing! I've had a good search around and haven't found a previous thread, so thought I'd start a new one.

I am using Gauss-Seidal and Jacobi iterative methods to calculate a simple square duct flow in one and two dimensions. In theory, in my mind, the 1D case should result in precisely the same flow characteristics as the 2D does down its symmetrical center. In practice I have found that the max velocities found differ quite significantly from one another, and I do not understand why.

I am coding these problems using a Fortran code given by my supervisor. As such, I believe the code to be accurate having checked it very closely. Which leads me to the conclusion that there is a known difference between calculating 1D and 2D flows using these methods?

If any of you can shine some light on this matter it would be terrifically appreciated. Is there indeed a reason as to why this is happening, and if so can you elaborate or point me to material that can explain further?

Many thanks,
u05harrisb is offline   Reply With Quote

Old   November 20, 2016, 04:01
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,061
Rep Power: 64
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
To tell the true, I do not understand what you want to compare... 1D and 2D model produce different physics... The 1D duct is a model where the tangential stress can only be "modelled" in the sense of Fanno theory. The 2D model produce a flow rate (and a velocity at centerline) that depends on the effect of the stress.

What kind of physical models are you using? For sure, the problem is not the iterative solver.
FMDenaro is offline   Reply With Quote


coding, duct flow, fortran, gauss-seidel, jacobian matrix

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

All times are GMT -4. The time now is 01:49.