CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Programming & Development

U, phi

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

LinkBack Thread Tools Display Modes
Old   March 15, 2013, 05:44
Default U, phi
New Member
Join Date: Nov 2012
Posts: 5
Rep Power: 7
ariaballack is on a distinguished road
hi everybody,
I am solving concentration using diffusion equation and for the convection term I need the velocity (U) . but after 2 seconds, the velocity remains constant (steady) and solving it anymore is just wasting time. Does anyone know how can I read phi after 2 second and continue without solving U anymore?
ariaballack is offline   Reply With Quote

Old   March 15, 2013, 06:29
Senior Member
Join Date: Sep 2009
Location: Delft
Posts: 790
Rep Power: 15
Bernhard is on a distinguished road
I don't know how important your start-up phase is (I assume it isn't important), but you can

*write a solver that first calculates the velocity field, and afterwards calculate the diffusion/convection.
* update the present solver with an if( t>2 ) around the solution step for U
* use a solver for flow first, then switch to a seperate solver for convection/diffusion.
Bernhard is offline   Reply With Quote

Old   March 21, 2013, 16:32
Jim Knopf
Join Date: Dec 2010
Posts: 60
Rep Power: 9
JimKnopf is on a distinguished road
A not so elegant but very fast way would be, to write a if statement around the #include "uEqn.H" or the code where you solve the momentum equation. You should also stopp solving the pressure correction.

This approch increases the computational work by one if statement every iteration, wouldn't be very flexible and only valid if U and p are constant after some time. But it is done within a minute.

JimKnopf 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
phi -= pEqn.flux() vs. linearInterpolate(U) & mesh.Sf() santiagomarquezd OpenFOAM Programming & Development 32 June 12, 2014 01:50
Recreation of phi file, icoFoam Madeleine P. Vincent OpenFOAM Running, Solving & CFD 0 May 7, 2012 11:17
Turbulence Model phi vs phi_ doug OpenFOAM Running, Solving & CFD 4 November 10, 2009 05:33
Another phi question ehsan_vaghefi OpenFOAM Running, Solving & CFD 0 October 24, 2008 19:56

All times are GMT -4. The time now is 02:44.