CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Inexplicable behaviour of Crank-Nicholson Scheme

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 8, 2009, 06:18
Default Inexplicable behaviour of Crank-Nicholson Scheme
  #1
New Member
 
Takeshi Omori
Join Date: Apr 2009
Posts: 2
Rep Power: 0
tomori is on a distinguished road
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
Attached Files
File Type: gz code_with_case_files.tar.gz (6.3 KB, 2 views)
tomori is offline   Reply With Quote

Old   June 14, 2009, 20:00
Default My applogies
  #2
New Member
 
Takeshi Omori
Join Date: Apr 2009
Posts: 2
Rep Power: 0
tomori is on a distinguished road
Starting from some reasonable initial field and avoiding a "rocket-start" solution, Crank-Nicholson was all right. My deepest appologies.
tomori is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Crank Nicolson scheme msrinath80 OpenFOAM Running, Solving & CFD 6 November 14, 2010 14:59
Godunov scheme Peter Main CFD Forum 6 November 27, 2009 11:51
extrapolation in MUSCL scheme Chandra Main CFD Forum 6 February 14, 2007 12:21
Crank Nicolson scheme contd msrinath80 OpenFOAM Running, Solving & CFD 0 December 12, 2006 23:41
Crank Nicholson noName Main CFD Forum 3 September 15, 2004 09:10


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