max time step related to diffusion time
I am interested in adjusting the time step for a transient solver on both the Courant number (already implemented in Openfoam) and on the time required for diffusion to be captured. Right now Im starting simple, using the icoFoam solver with species transport. Since I would like to add this to a turbulent solver later (diffusivity depends on turbulent schmidt number), I have made diffusivity a field. I want to find the average diffusion time and the minimum diffusion time based on:
t_diffusion = dx^2/D
I looked at the CourantNo.H file and found:
surfaceScalarField SfUfbyDelta =
CoNum = max(SfUfbyDelta/mesh.magSf())
meanCoNum = (sum(SfUfbyDelta)/sum(mesh.magSf()))
velMag = max(mag(phi)/mesh.magSf()).value();
but I need to extract the delta X values, square them, and then divide them by the cell specific diffusivity. How does one do this exactly? Any help is much appreciated.
What do you mean exactly by delta X values?
dalta x as in h, spacing, or cell-cell distance
Ok then you can get it from
This is a surfaceScalarField with value 1/delta - the cell-centre to cell-centre distances. (Or face centre to cell centre distance at the boundary.)
|All times are GMT -4. The time now is 15:04.|