CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Pre-Processing

waveTransmissive steadyState error

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 6, 2015, 05:59
Default waveTransmissive steadyState error
  #1
New Member
 
Sinasi Bahceci
Join Date: Jan 2015
Posts: 13
Rep Power: 11
sino75 is on a distinguished road
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.
sino75 is offline   Reply With Quote

Old   March 6, 2015, 07:10
Default
  #2
Senior Member
 
Alexey Matveichev
Join Date: Aug 2011
Location: Nancy, France
Posts: 1,938
Rep Power: 38
alexeym has a spectacular aura aboutalexeym has a spectacular aura about
Send a message via Skype™ to 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:

Code:
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.
alexeym is offline   Reply With Quote

Old   March 9, 2015, 05:09
Default
  #3
New Member
 
Sinasi Bahceci
Join Date: Jan 2015
Posts: 13
Rep Power: 11
sino75 is on a distinguished road
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?
sino75 is offline   Reply With Quote

Old   March 9, 2015, 05:14
Default
  #4
Senior Member
 
Alexey Matveichev
Join Date: Aug 2011
Location: Nancy, France
Posts: 1,938
Rep Power: 38
alexeym has a spectacular aura aboutalexeym has a spectacular aura about
Send a message via Skype™ to 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?
alexeym is offline   Reply With Quote

Old   March 9, 2015, 07:28
Default
  #5
New Member
 
Sinasi Bahceci
Join Date: Jan 2015
Posts: 13
Rep Power: 11
sino75 is on a distinguished road
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 ;
}


// ************************************************** *********************** //
sino75 is offline   Reply With Quote

Old   March 9, 2015, 07:35
Default
  #6
Senior Member
 
Alexey Matveichev
Join Date: Aug 2011
Location: Nancy, France
Posts: 1,938
Rep Power: 38
alexeym has a spectacular aura aboutalexeym has a spectacular aura about
Send a message via Skype™ to 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.
alexeym is offline   Reply With Quote

Old   March 9, 2015, 08:01
Default
  #7
New Member
 
Sinasi Bahceci
Join Date: Jan 2015
Posts: 13
Rep Power: 11
sino75 is on a distinguished road
Now I understand, thanks a lot.
sino75 is offline   Reply With Quote

Reply

Tags
differencing scheme, openfoam, steadystate, wavetransmissive

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
[OpenFOAM.org] compile error in dynamicMesh and thermophysicalModels libraries NickG OpenFOAM Installation 3 December 30, 2019 00:21
[blockMesh] blockMesh with double grading. spwater OpenFOAM Meshing & Mesh Conversion 92 January 12, 2019 09:00
[swak4Foam] GroovyBC the dynamic cousin of funkySetFields that lives on the suburb of the mesh gschaider OpenFOAM Community Contributions 300 October 29, 2014 18:00
OpenFOAM without MPI kokizzu OpenFOAM Installation 4 May 26, 2014 09:17
DecomposePar links against liblamso0 with OpenMPI jens_klostermann OpenFOAM Bugs 11 June 28, 2007 17:51


All times are GMT -4. The time now is 08:50.