|
[Sponsors] |
Initial Boundary Conditions setup for hypersonic flow of a double wedge |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
July 8, 2021, 02:43 |
Initial Boundary Conditions setup for hypersonic flow of a double wedge
|
#1 |
New Member
Matthew
Join Date: May 2021
Posts: 10
Rep Power: 5 |
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! Last edited by Tegdlemat; July 8, 2021 at 02:45. Reason: typo |
|
August 4, 2021, 11:00 |
|
#2 |
New Member
Niklas Boehnke
Join Date: Jul 2021
Location: Lübeck
Posts: 12
Rep Power: 5 |
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. |
|
May 13, 2022, 15:17 |
Getting exactly the same problem
|
#3 |
New Member
Fuad Hasibul Hasan
Join Date: Feb 2022
Location: Dhaka, Bangladesh
Posts: 7
Rep Power: 4 |
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. |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Table bounds warnings at: END OF TIME STEP | CFXer | CFX | 4 | July 17, 2020 00:44 |
Maximum number of iterations exceeded chtmultiregionsimpleFoam | Moncef | OpenFOAM Running, Solving & CFD | 28 | July 13, 2020 15:26 |
chtMultiRegionSimpleFoam turbulent case | Aditya Patil | OpenFOAM Running, Solving & CFD | 6 | April 24, 2017 23:13 |
Wrong flow in ratating domain problem | Sanyo | CFX | 17 | August 15, 2015 07:20 |
Error finding variable "THERMX" | sunilpatil | CFX | 8 | April 26, 2013 08:00 |