[Tegdlemat]

Currently I am working on a project about the stability of hypersonic flow across a double wedge, shown as below:
image_2021-07-08_133636.png
I modified the geometry according to the wedge tutorial blockMesh file by OpenFOAM. However since the bc of the tutorial doesnt match the condition I wish to test obviously, I was trying to modify the files

U:
[QUOTE]/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v2106 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 1 -1 0 0 0 0];

internalField uniform (2382.03 0 0);

boundaryField
{
inlet
{
type fixedValue;
value uniform (2382.03 0 0);
}

outlet
{
type zeroGradient;
}

bottom
{
type symmetryPlane;
}

top
{
type symmetryPlane;
}

obstacle
{
type slip;
}

defaultFaces
{
type empty;
}
}


// ************************************************** *********************** //
p:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v2106 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [1 -1 -2 0 0 0 0];

internalField uniform 391;

boundaryField
{
inlet
{
type fixedValue;
value uniform 391;
}

outlet
{
type zeroGradient;
}

bottom
{
type symmetryPlane;
}

top
{
type symmetryPlane;
}

obstacle
{
type zeroGradient;
}

defaultFaces
{
type empty;
}
}


// ************************************************** *********************** //
T:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v2106 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 0 0 1 0 0 0];

internalField uniform 300;

boundaryField
{
inlet
{
type fixedValue;
value uniform 300;
}

outlet
{
type zeroGradient;
}

bottom
{
type symmetryPlane;
}

top
{
type symmetryPlane;
}

obstacle
{
type zeroGradient;
}

defaultFaces
{
type empty;
}
}


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

I read the manual on OpenFOAM which stated that the units of U is m/s, T is K, and p is Pa. And when I tried to input the p as 391 Pa, U as 2382 m/s and T as 300K, then decomposePar and run mpiexec, all I get is:

Primary job terminated normally, but 1 process returned
a non-zero exit code. Per user-direction, the job has been aborted.
--------------------------------------------------------------------------
--------------------------------------------------------------------------
mpiexec detected that one or more processes exited with non-zero status, thus causing
the job to be terminated. The first process to do so was:

Process name: [[2557,1],0]
Exit code: 144
--------------------------------------------------------------------------

I checked my mesh already with the tutorial case which is fine, so I suspect its about the bc setups. Since I am new to hypersonic and rhoCentralFoam, I am just simulating laminar flow atm. I would like to know how to correct my bc setups to get rid of the error.

Thanks guys!

[janibo]

Hi Matthew,

I'm only an OF beginner myself, but what i have noticed is that
you have overspecified your problem, since you have given 3 of 3
bc's at the inlet. You should only specify 2 of 3 bc's and have the
third calculated e.g. for U using the pressureInletOutletVelocity bc.

[fuad]

If you have solved this issue please share it here. I am having exactly the same issue.
The compressible/rhoCentralFoam/forwardStep case in OF-2106 just fine until I change the 0/U, 0/p, 0/T, constant/thermoPhysicalProperties files. I want to use S.I. units where the tutorial uses normalised parameters. Is it possible to use values in S.I. units of parameters for this solver? If not, how can I normalize boundary conditions values and thermophysical properties.