# rhoCentralFoam - Issues when changing setup from symmetryPlane to full 2D geometry

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

December 5, 2023, 12:11
rhoCentralFoam - Issues when changing setup from symmetryPlane to full 2D geometry
#1
New Member

Join Date: Aug 2022
Location: Mexico
Posts: 16
Rep Power: 3
Hello all!

I'm currently engaged in a supersonic simulation project using rhoCentralFoam. However, most tutorial cases utilize a symmetryPlane, and my simulation involves an irregular geometry that cannot be treated with symmetry.

While setting up boundary conditions and names, I encounter difficulties running the case. I consistently face errors or encounter floating-point issues.

Essentially, I'm aiming to conduct a simulation using a logic similar to aerofoil analysis. Could you provide guidance on how to address these challenges and make the simulation successful? Right now I'm using the forwardStep tutorial case as a base.

I'm attaching a picture of the setup for reference and also a result using sonicFoam.

My data is below:

Pressure
Code:
```internalField   uniform 1;

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

outlet
{
}

bottom
{
type            symmetryPlane;
}

top
{
type            symmetryPlane;
}

obstacle
{
}

defaultFaces
{
type            empty;
}
}```
Temperature
Code:
```dimensions      [0 0 0 1 0 0 0];

internalField   uniform 1;

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

outlet
{
type            inletOutlet;
inletValue      uniform 1;
value           uniform 1;
}

bottom
{
type            symmetryPlane;
}

top
{
type            symmetryPlane;
}

obstacle
{
}

defaultFaces
{
type            empty;
}
}```
Velocity
Code:
```dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (3 0 0);

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

outlet
{
type            inletOutlet;
inletValue      uniform (7 0 0);
value           uniform (7 0 0);
}

bottom
{
type            symmetryPlane;
}

top
{
type            symmetryPlane;
}

obstacle
{
type            slip;
}

defaultFaces
{
type            empty;
}
}```
controlDict
Code:
```application     rhoCentralFoam;

startFrom       startTime;

startTime       0;

stopAt          endTime;

endTime         4;

deltaT          0.002;

writeInterval   0.1;

purgeWrite      10;

writeFormat     ascii;

writePrecision  6;

writeCompression off;

timeFormat      general;

timePrecision   6;

runTimeModifiable true;

maxCo           0.2;

maxDeltaT       1;```
thermophysicalProperties (I've tried to change these to "normal" air conditions with no luck)
Code:
```thermoType
{
type            hePsiThermo;
mixture         pureMixture;
transport       const;
thermo          hConst;
equationOfState perfectGas;
specie          specie;
energy          sensibleInternalEnergy;
}

// Note: these are the properties for a "normalised" inviscid gas
//       for which the speed of sound is 1 m/s at a temperature of 1K
//       and gamma = 7/5
mixture
{
specie
{
molWeight       11640.3;
}
thermodynamics
{
Cp              2.5;
Hf              0;
}
transport
{
mu              0;
Pr              1;
}
}```
Attached Images
 Setup.png (14.4 KB, 5 views) sonicFoam success.png (119.9 KB, 6 views)