|
[Sponsors] |
Inexplicable behaviour of Crank-Nicholson Scheme |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
April 8, 2009, 06:18 |
Inexplicable behaviour of Crank-Nicholson Scheme
|
#1 |
New Member
Takeshi Omori
Join Date: Apr 2009
Posts: 2
Rep Power: 0 |
Dear all,
when I run transientSimpleFoam (without any turbulence models, codes are attached) on a very trivial case where the flow is just going through a 2D friction-free channel (blockMeshDict and boundary condition files etc. are attached too), I observe unphysical uniform pressure gradient in the direction of the flow. As the calculation proceeds in time, the pressure gradient tends to vanish, but even after 100 Flow-Through-Time it still stays there. Another thing is that when I edit icoFoam as the following, phi = (fvc::interpolate(U) & mesh.Sf()); // + fvc::ddtPhiCorr(rUA, U, phi); I obtain the same unphysical pressure gradient. As far as I know, the term I omitted is for the calculation with a very small time step width and not supposed to behave like this. Let me also note that the pressure gradient appears even when the calculation starts from the uniform velocity field with the inlet velocity magnitude which can be obtained by another solver or whatever. When I run this icoFoam with the implicit Euler scheme, the unphysical pressure gradient does not appear. Unstability is, however, no excuse on this very trivial case (Courant number was set less than 1, just for info). If anyone would tell me why this happens, I deeply appreciate it. Thank you, tomori |
|
June 14, 2009, 20:00 |
My applogies
|
#2 |
New Member
Takeshi Omori
Join Date: Apr 2009
Posts: 2
Rep Power: 0 |
Starting from some reasonable initial field and avoiding a "rocket-start" solution, Crank-Nicholson was all right. My deepest appologies.
|
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Crank Nicolson scheme | msrinath80 | OpenFOAM Running, Solving & CFD | 6 | November 14, 2010 13:59 |
Godunov scheme | Peter | Main CFD Forum | 6 | November 27, 2009 10:51 |
extrapolation in MUSCL scheme | Chandra | Main CFD Forum | 6 | February 14, 2007 11:21 |
Crank Nicolson scheme contd | msrinath80 | OpenFOAM Running, Solving & CFD | 0 | December 12, 2006 22:41 |
Crank Nicholson | noName | Main CFD Forum | 3 | September 15, 2004 09:10 |