|July 5, 2013, 04:11||
Sub-iterations for transport equation in drift-flux
Join Date: Jul 2013
Posts: 3Rep Power: 2
I am currently working on a drift-flux (or algebraic-slip) model to simulate the sedimentation of sand in water.
After solving the transport equation for sand, I use the packingLimiter.H like in twoPhaseEulerFoam in order to maintain the maximum packing concentration of sand (say 0.6).
If I take the time-step dt too large, a certain computational cell can be filled with an amount of sand such that the packingLimiter cannot reduce the sand to the maximum packing concentration. If the packing concentration exceeds a value of 1.0, the solving process is aborted by the solver.
I am wondering if I could decrease the timestep dt only for solving the transport equation. Say, 10 sub-iterations on the transport equation are chosen with for each sub-iteration the timestep 0.1*dt. Each time the packingLimiter is applied after a sub-iteration.
Does anybody know how I can call the magnitude of the timestep in the source code, such that I can modify it?
Thanks in advance!
|July 5, 2013, 05:03||
Join Date: Oct 2011
Location: Göteborg, Sweden
Posts: 10Rep Power: 4
For the twoPhaseEulerFoam solver the time step is set from the include of:
For subiterations, perhaps you would find some interest in: how to program two independent time loops.
I guess you would like to use subCycles to achieve subiterations for some part of the code. Have a look at e.g the alphaEqn.H in twoPhaseEulerFoam. There a subCycle is used to solve the dispersed phase continuity equation. Inside the loop the deltaT is split according to how many sub-iterations you have specified. The code, subCycle.H and subCycleTime.H, will help you with the time step splitting.
|drift-flux, sub-iterations, transport equation|
|Thread||Thread Starter||Forum||Replies||Last Post|
|pimpleFoam: turbulence->correct(); is not executed when using residualControl||hfs||OpenFOAM Running, Solving & CFD||3||October 29, 2013 10:35|
|Moving mesh||Niklas Wikstrom (Wikstrom)||OpenFOAM Running, Solving & CFD||121||March 7, 2013 18:21|
|Interfoam blows on parallel run||danvica||OpenFOAM Running, Solving & CFD||16||December 22, 2012 03:09|
|Differences between serial and parallel runs||carsten||OpenFOAM Bugs||11||September 12, 2008 12:16|
|Convergence moving mesh||lr103476||OpenFOAM Running, Solving & CFD||30||November 19, 2007 15:09|