[sino75]

Hello,

I'm trying to use waveTransmissive for my wall, but I get following error:

--> FOAM FATAL ERROR:
Unsupported temporal differencing scheme : steadyState
on patch WALL of field p in file "/home/.../OpenFOAM/.../0/p"

From function advectiveFvPatchField<Type>::updateCoeffs()
in file fields/fvPatchFields/derived/advective/advectiveFvPatchField.C at line 255.

FOAM exiting

I tried to read the function, but I had no success. Can anyone help me please?

WALL
{
type waveTransmissive;
value uniform 80000;
field p;
gamma 1.4;
phi phi;
rho rho;
psi thermosi;
lInf 0.5;
fieldInf 80000;
}

Thanks in regard.

[alexeym]

Hi,

Well, as waveTransmissive is a child of advective fvPatchField (and in fact the error happens in advectiveFvPatchField constructor), here the list of supported time differencing schemes:

if (lInf_ > 0)
    {
        // Calculate the field relaxation coefficient k (See notes)
        const scalarField k(w*deltaT/lInf_);

        if
        (
            ddtScheme == fv::EulerDdtScheme<scalar>::typeName
         || ddtScheme == fv::CrankNicolsonDdtScheme<scalar>::typeName
        )
        {
            this->refValue() =
            (
                field.oldTime().boundaryField()[patchi] + k*fieldInf_
            )/(1.0 + k);

            this->valueFraction() = (1.0 + k)/(1.0 + alpha + k);
        }
        else if (ddtScheme == fv::backwardDdtScheme<scalar>::typeName)
        {
            this->refValue() =
            (
                2.0*field.oldTime().boundaryField()[patchi]
              - 0.5*field.oldTime().oldTime().boundaryField()[patchi]
              + k*fieldInf_
            )/(1.5 + k);

            this->valueFraction() = (1.5 + k)/(1.5 + alpha + k);
        }
        else
        {
             // Error message here
        }

There is no oldTime in steadyState calculations but it is necessary for calculation of refValue.

[sino75]

Thanks Alex,

sorry for my late response. I had to check some points. I have changed ddt Scheme from steadyState to Euler. But I get the same failure. I'm quite new to OF, so could you please help?

[alexeym]

Well,

1. What solver do you use?
2. Please, post your fvSchemes.
3. "same failure" is rather vague description. Could you post the error?

[sino75]

Hi,
I use the flameletSimpleFoam solver. I'm just running something, so I will paste the error message as soon as it's finished. It's the same which is given in the code you posted ('Error message here').

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.2.x |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.2;
format ascii;
class dictionary;
location "system";
object fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
default Euler;
}

gradSchemes
{
default cellLimited Gauss linear 1;
}

divSchemes
{
div(phi,U) bounded Gauss limitedLinear 1;
div(phi,epsilon) bounded Gauss upwind;
div(phi,k) bounded Gauss upwind;

div(phi,H) bounded Gauss limitedLinear 1;
div(phi,Z) bounded Gauss limitedLimitedLinear 1 0 1;
div(phi,Zvar) bounded Gauss limitedLimitedLinear 1 0 0.25;

div((muEff*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{

laplacian(muEff,U) Gauss linear corrected;
laplacian((rho*(1|A(U))),p) Gauss linear corrected;
laplacian(DkEff,k) Gauss linear corrected;

laplacian((rho|A(U)),p) Gauss linear corrected;

laplacian(DepsilonEff,epsilon) Gauss linear corrected;
laplacian((muEff|sigmat),Z) Gauss linear corrected;
laplacian((muEff|sigmat),H) Gauss linear corrected;
laplacian((mut|sigmat),Zvar) Gauss linear corrected;
}

interpolationSchemes
{
default linear;
}

snGradSchemes
{
default corrected;
}

fluxRequired
{
default no;
p ;
}


// ************************************************** *********************** //

[alexeym]

Guessing from the name flameletSimpleFoam is steady state solver, waveTransmissive BC can't be used with steady-state solvers, as there are no fvm::ddt(...) calls, so no oldTime values of the field.

[sino75]

Now I understand, thanks a lot.