Turbulence Model with Elliptic Relaxation
I've been working on implementing a damping function, fmu, defined by an elliptic relaxation equation. This is the first modification I've done that has required defining my own discretized equation so I'm sure there are errors in my implementation. Also I believe my method of setting the initial and boundary conditions could use improvement. Currently my simpleFoam case compiles and runs but the results are incorrect.
fmu must be between 0 and 1 and equal 0 on wall boundaries. To initialize this I added the code below to the constructor.
In /0/fmu I set a value of 0 to my wall patch and 0.5 to the internalField and all other patches.
label patchWall = mesh_.boundaryMesh().findPatchID("wall")
The fmu equation is
Does this procedure look correct? Is using the boudaryField()[patchWall]=0.0 the proper way to correct the wall boundary condition?
Thanks in advance for any input,
You can do something like this so that you do not need to create Initial Condition files for the fields you want to get printed out:
dimensionedScalar("LPrint", dimensionSet(0, 1, 0, 0, 0, 0, 0), 0.0),
Thanks for the input.
Also, as an update, I use the following to solve the elliptic equation. It seems to be working, but I'm still testing.
|All times are GMT -4. The time now is 22:48.|