Trouble with simple case for a new OpenFOAM user
 

Hello all,

I am a long time lurker on these forums, and I've joined now hoping to get some help with OpenFOAM. I've recently started trying to learn the software in order to repeat some work for graduate school, and I am lost setting up the problem.

I have a free shear layer, with high- and low-speed inlets on the left side, an outlet on the right side, the low-speed freestream at the bottom boundary, and the high-speed freestream at the top boundary. All I want is to run LES and DES simulations of this problem. I've done so in Fluent and Cobalt, but OpenFOAM has me baffled.

At each inlet, I'm trying to define a velocity and static temperature profile with approximately a 5mm boundary layer. From what I have found in sources online, I've written .csv files for these profiles, and I've used a fixedProfile BC. However, when following online tutorials, I get many errors, and a lot of examples online don't work.

Is there someone willing to take a look at my case and help me to understand the simple things I must be missing in order to get this case going? I would greatly appreciate any help. Attached are my 0/U and 0/T files (uploaded as txt files)

Hello,

What you are doing seems correct but I may suggest to you to check the U and T via ParaView by de-activating "Skip zero folder". So, you may see your mistake (if any). Otherwise, write the error that you get, so it might be easier to spot the problem for us, I guess.

Thank you for the reply.

I have tried to view the BC's in paraFoam, and I've found the "skip zero folder" option, but even when I uncheck that box and click apply, no options to view T or U are present in the drop down box at the top. I'm not sure why it isn't showing up, apart from a syntax error that I'm not spotting.

When I try to run the case using rhoPimpleFoam, I get the following output:

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


PIMPLE: no residual control data found. Calculations will employ 3 corrector loops

Reading thermophysical properties

Selecting thermodynamics package
{
type hePsiThermo;
mixture pureMixture;
transport sutherland;
thermo eConst;
equationOfState perfectGas;
specie specie;
energy absoluteEnthalpy;
}

#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigSegv::sigHandler(int) at ??:?
#2 ? in "/lib64/libpthread.so.0"
#3 std::string::assign(std::string const&) in "/lib64/libstdc++.so.6"
#4 Foam::basicThermo::splitThermoName(Foam::word const&, int) at ??:?
#5 Foam::HashTable<Foam::autoPtr<Foam::fluidThermo> (*)(Foam::fvMesh const&, Foam::word const&), Foam::word, Foam::string::hash>::iterator Foam::basicThermo::lookupThermo<Foam::fluidThermo, Foam::HashTable<Foam::autoPtr<Foam::fluidThermo> (*)(Foam::fvMesh const&, Foam::word const&), Foam::word, Foam::string::hash> >(Foam::dictionary const&, Foam::HashTable<Foam::autoPtr<Foam::fluidThermo> (*)(Foam::fvMesh const&, Foam::word const&), Foam::word, Foam::string::hash>*, int, char const**, Foam::word const&) at ??:?
#6 Foam::HashTable<Foam::autoPtr<Foam::fluidThermo> (*)(Foam::fvMesh const&, Foam::word const&), Foam::word, Foam::string::hash>::iterator Foam::basicThermo::lookupThermo<Foam::fluidThermo, Foam::HashTable<Foam::autoPtr<Foam::fluidThermo> (*)(Foam::fvMesh const&, Foam::word const&), Foam::word, Foam::string::hash> >(Foam::dictionary const&, Foam::HashTable<Foam::autoPtr<Foam::fluidThermo> (*)(Foam::fvMesh const&, Foam::word const&), Foam::word, Foam::string::hash>*) at ??:?
#7 Foam::fluidThermo::New(Foam::fvMesh const&, Foam::word const&) at ??:?
#8 ? at ??:?
#9 __libc_start_main in "/lib64/libc.so.6"
#10 ? at ??:?
Segmentation fault