[Henning86]

hi,



i made some modification on the wavetransmissive/advective Bc to include the localEuler ddt scheme:

it is based on the thread( and 2 links for local Euler for rhoCentralFoam):


http://www.cfd-online.com/Forums/ope...-openfoam.html

http://www.cfd-online.com/Forums/ope...nge-kutta.html

--- ../OpenFOAM-2.1.1/src/finiteVolume/fields/fvPatchFields/derived/advective/advectiveFvPatchField.C +++ /opt/openfoam211/src/finiteVolume/fields/fvPatchFields/derived/advective/advectiveFvPatchField.C @@ -30,6 +30,7 @@
#include "EulerDdtScheme.H"
#include "CrankNicholsonDdtScheme.H"
#include "backwardDdtScheme.H"
+#include "localEulerDdtScheme.H" //add
 
 
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
@@ line 215 ++@@
 
scalarField alpha(w*deltaT*this->patch().deltaCoeffs()); // removed the const
label patchi = this->patch().index();
// Non-reflecting outflow boundary
// If lInf_ defined setup relaxation to the value fieldInf_.
if (lInf_ > SMALL)
{
// Calculate the field relaxation coefficient k (See notes)
scalarField k(w*deltaT/lInf_); // removed the const
 
 
@@ line 251++@@
//insert this:
else if (ddtScheme == fv::localEulerDdtScheme<scalar>::typeName)
{
// only works with rDeltaT hardcoded
const volScalarField& vsfrDeltaT =
this->db().objectRegistry::lookupObject<volScalarField>("rDeltaT");
const scalarField& prDeltaT = vsfrDeltaT.boundaryField()[patchi];
scalarField alpha(w/prDeltaT*this->patch().deltaCoeffs());
scalarField k(w/prDeltaT/lInf_); 
 
this->refValue() =
(
field.oldTime().boundaryField()[patchi] + k*fieldInf_
)/(1.0 + k);
this->valueFraction() = (1.0 + k)/(1.0 + alpha + k);
}
@@ line 295++@@
//insert this:
 
else if (ddtScheme == fv::localEulerDdtScheme<scalar>::typeName)
{
// only works with rDeltaT hardcoded
const volScalarField& vsfrDeltaT =
this->db().objectRegistry::lookupObject<volScalarField>("rDeltaT");
const scalarField& prDeltaT = vsfrDeltaT.boundaryField()[patchi];
scalarField alpha(w/prDeltaT*this->patch().deltaCoeffs());
this->refValue() = field.oldTime().boundaryField()[patchi];
this->valueFraction() = 1.0/(1.0 + alpha);
}

if you want to use the wave transmissive BC with LTS you need to recompile the libary "finiteVolume". i would recommend it to recompile it in at the $(LIB_SRC) (main folder). If you do it in your local folder (FOAM_USER_LIBBIN) you get a double free or corruption error at the end of the run because this libary is pretty much used by every other libary


dont forget to make a backup before before compiling

Best regards

Henning