CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Bugs

Bug caused by CrankNicholson scheme

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

Like Tree1Likes
  • 1 Post By alberto

Reply
 
LinkBack Thread Tools Display Modes
Old   September 3, 2008, 17:56
Default Using the Crank-Nicholson Sche
  #1
New Member
 
Chris Beck
Join Date: Mar 2009
Posts: 4
Rep Power: 8
cbeck is on a distinguished road
Using the Crank-Nicholson Scheme in compressible transient solvers causes a unit problem in the pressure correction.

I can reproduce the problem with the coodles tutorial. Can anybody else confim this?

Version:
/openfoam-extend/branches/OpenCFD_Release/OpenFOAM-1.5_patch/ (a couple of weeks old)

kind regards,
Christian Beck
cbeck is offline   Reply With Quote

Old   September 4, 2008, 03:41
Default Thanks for the bug report. I
  #2
Senior Member
 
Join Date: Mar 2009
Posts: 854
Rep Power: 13
henry is on a distinguished road
Thanks for the bug report. I have fixed the problem and running the coodles tutorial at the moment. If all is well I will upload the patch to our 1.5.x git repository (see http://www.opencfd.co.uk/openfoam/download.html for download details) later today.

H
henry is offline   Reply With Quote

Old   September 12, 2008, 09:41
Default I observed the same problem wi
  #3
Member
 
Carsten Thorenz
Join Date: Mar 2009
Location: Germany
Posts: 32
Rep Power: 8
carsten is on a distinguished road
I observed the same problem with lesInterFoam: If I switch to Crank-Nicholson, a dimension error appears.

Bye,

Carsten
carsten is offline   Reply With Quote

Old   September 12, 2008, 11:07
Default Even when using the version I
  #4
Senior Member
 
Join Date: Mar 2009
Posts: 854
Rep Power: 13
henry is on a distinguished road
Even when using the version I fixed and uploaded to our git repository?

H
henry is offline   Reply With Quote

Old   September 15, 2008, 05:32
Default Uhhhh. Difficult to say, sorry
  #5
Member
 
Carsten Thorenz
Join Date: Mar 2009
Location: Germany
Posts: 32
Rep Power: 8
carsten is on a distinguished road
Uhhhh. Difficult to say, sorry. I can't reconstruct exactly when I did check-out and I fear to check it out again, because I don't want to loose my settings/changes, which were necessary for compilation on our machine.

Anyhow, I will try to recompile in another directory and will report my findings.

Bye,

Carsten
carsten is offline   Reply With Quote

Old   September 15, 2008, 11:02
Default O.k., problem still exists. I
  #6
Member
 
Carsten Thorenz
Join Date: Mar 2009
Location: Germany
Posts: 32
Rep Power: 8
carsten is on a distinguished road
O.k., problem still exists. I checked out a new tree this morning and still have the problem.

You can recreate it with the lesInterFoam-tutorial, simply switch to CrankNicholson and it will die.

Bye,

Carsten
carsten is offline   Reply With Quote

Old   September 15, 2008, 11:48
Default OK I will look into it... H
  #7
Senior Member
 
Join Date: Mar 2009
Posts: 854
Rep Power: 13
henry is on a distinguished road
OK I will look into it...

H
henry is offline   Reply With Quote

Old   March 27, 2009, 05:43
Default
  #8
Member
 
Virginie Ehrlacher
Join Date: Mar 2009
Posts: 52
Rep Power: 8
virginie_e is on a distinguished road
Hello,

I also have a bug when trying to use the crankNicholson scheme with interTrackFoam. However, it is not a dimension problem. The error message is rather:

attempt to read beyond EOF
file: /users/V1117324/OpenFOAM/v1117324-1.5-dev/run/tutorials/interTrackFoam/hydrofoil/system/fvSchemes::ddt(rho,U) at line 21.
From function ITstream::read(token& t)
in file db/IOstreams/Tstreams/ITread.C at line 64.
FOAM exiting


Do you know what this means? Thank you.

Virginie
virginie_e is offline   Reply With Quote

Old   March 28, 2009, 02:48
Default
  #9
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
Hi Virginie_e,
that's not a bug but a typo or an error in the specification of the scheme. Did you specify a number between 0 (Euler) and 1 (pure Crank Nicholson) after the scheme name?
If this is not the case, please post your fvScheme file in the appropriate forum. This one is specific for bugs

Best,
mm.abdollahzadeh likes this.
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto 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
CrankNicholson errors lr103476 OpenFOAM Running, Solving & CFD 12 October 20, 2013 06:21
CrankNicholson returns wrong dimensions and sound speed when mesh is moving lillberg OpenFOAM Bugs 6 May 4, 2009 17:53
Forces caused by fluids anja OpenFOAM Running, Solving & CFD 20 May 8, 2006 06:37
Adiabatic compression caused by gas ram Chris FLUENT 0 December 13, 2005 05:42
Probelm with CrankNicholson scheme in channelOodles maka OpenFOAM Pre-Processing 2 November 2, 2005 11:09


All times are GMT -4. The time now is 09:58.