CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Ddt schemes update

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

LinkBack Thread Tools Display Modes
Old   April 11, 2005, 07:55
Default Here's a set of messages regar
Senior Member
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,810
Rep Power: 25
hjasak will become famous soon enough
Here's a set of messages regarding the ddt scheme update:



Attached is an update for ddt schemes. As I've indicated, once you start switching ddt schemes on a moving mesh, the mesh motion flux is calculated differently for each ddt scheme. Therefore, for anything other than Euler implicit mesh.phi() will not do the right thing.

Zeljko has introduced a new fvc::meshPhi(...) which takes a field as an argument. The field is there to find which ddtScheme is currently switched on and it will return the appropriate mesh motion flux.

In solvers (I did not pack those, do you want them?) fvc::meshPhi(U) (typically, you want U) should be used instead of mesh.phi().

You will notice some strangemess in handling mesh motion data in the fvMesh.C - that's because mesh.phi() now carries old-tome motion fluxes as well and needs to be handled with care - I suspect the same kind of problem will show up when we get dimensionedFields which carry oldTime() for mesh.V().


I need to consider these changes carefully because they may be incompatible with my plans to introduce an extended substantive derivative to avoid the algorithm differences for moving-mesh problems. I also want to create a super-class for U which incorporates phi so that all the relative/absolute phi crap is handled automatically. All this is in progress but will not be in the next release. The other issue is if there changes are compatible with the new phi predictor step I implemented for the pressure equation to avoid spacial-staggering on refined meshes.

Me again :-)

I will carry on with the current stuff for the moment - if someone needs higher order temporal discretisation on moving meshes any time soon, please give me a shout.

Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting:
hjasak is offline   Reply With Quote


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
FixedGradient BC update stefan82 OpenFOAM Running, Solving & CFD 4 April 28, 2009 03:10
Update Fluxes luca OpenFOAM Running, Solving & CFD 2 March 29, 2006 10:47
IcoTopoFoam update hjasak OpenFOAM Running, Solving & CFD 0 March 28, 2005 21:32
UDF update profile ramesh FLUENT 0 June 29, 2003 14:14
MOUSE - any update? Pei-Ying Hsieh Main CFD Forum 1 March 19, 2001 06:04

All times are GMT -4. The time now is 18:21.