[andyreagan]

Hi foamers,

I've been working to simulate a thermal convection loop, in other words a "hula hoop filled with water, heated on the bottom half and cooled on top half" for a little while now.

I have a lot of questions, and have been reading the forums for a few days trying to answer most of them, but I'm stuck with buoyantPimpleFoam not starting.

The case is pretty basic, in trying to get it to run I've narrowed down to the only difference from the hotRoom tutorial being my 2D meshed loop. The mesh works well in both Boussinesq solvers.

I've attached my case, and here is the error that I'm getting out of the solver (if I switch the RASModel to "laminar") the error delays until beginning the run.

Creating turbulence model

Selecting turbulence model type RASModel
Selecting RAS turbulence model kEpsilon
#0  Foam::error::printStack(Foam::Ostream&) in "/users/a/r/areagan/OpenFOAM/OpenFOAM-2.2.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/users/a/r/areagan/OpenFOAM/OpenFOAM-2.2.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2  ?? in "/lib64/libc.so.6"
#3  Foam::compressible::mutkWallFunctionFvPatchScalarField::calcMut() const in "/users/a/r/areagan/OpenFOAM/OpenFOAM-2.2.1/platforms/linux64GccDPOpt/lib/libcompressibleRASModels.so"
#4  Foam::compressible::mutWallFunctionFvPatchScalarField::updateCoeffs() in "/users/a/r/areagan/OpenFOAM/OpenFOAM-2.2.1/platforms/linux64GccDPOpt/lib/libcompressibleRASModels.so"
#5  Foam::fvPatchField<double>::evaluate(Foam::UPstream::commsTypes) in "/users/a/r/areagan/OpenFOAM/OpenFOAM-2.2.1/platforms/linux64GccDPOpt/bin/buoyantPimpleFoam"
#6  Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricBoundaryField::evaluate() in "/users/a/r/areagan/OpenFOAM/OpenFOAM-2.2.1/platforms/linux64GccDPOpt/bin/buoyantPimpleFoam"
#7  Foam::compressible::RASModels::kEpsilon::kEpsilon(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::fluidThermo const&, Foam::word const&, Foam::word const&) in "/users/a/r/areagan/OpenFOAM/OpenFOAM-2.2.1/platforms/linux64GccDPOpt/lib/libcompressibleRASModels.so"
#8  Foam::compressible::RASModel::adddictionaryConstructorToTable<Foam::compressible::RASModels::kEpsilon>::New(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::fluidThermo const&, Foam::word const&) in "/users/a/r/areagan/OpenFOAM/OpenFOAM-2.2.1/platforms/linux64GccDPOpt/lib/libcompressibleRASModels.so"
#9  Foam::compressible::RASModel::New(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::fluidThermo const&, Foam::word const&) in "/users/a/r/areagan/OpenFOAM/OpenFOAM-2.2.1/platforms/linux64GccDPOpt/lib/libcompressibleRASModels.so"
#10  Foam::compressible::turbulenceModel::addturbulenceModelConstructorToTable<Foam::compressible::RASModel>::NewturbulenceModel(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::fluidThermo const&, Foam::word const&) in "/users/a/r/areagan/OpenFOAM/OpenFOAM-2.2.1/platforms/linux64GccDPOpt/lib/libcompressibleRASModels.so"
#11  Foam::compressible::turbulenceModel::New(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::fluidThermo const&, Foam::word const&) in "/users/a/r/areagan/OpenFOAM/OpenFOAM-2.2.1/platforms/linux64GccDPOpt/lib/libcompressibleTurbulenceModel.so"
#12  main in "/users/a/r/areagan/OpenFOAM/OpenFOAM-2.2.1/platforms/linux64GccDPOpt/bin/buoyantPimpleFoam"
#13  __libc_start_main in "/lib64/libc.so.6"
#14  __gxx_personality_v0 in "/users/a/r/areagan/OpenFOAM/OpenFOAM-2.2.1/platforms/linux64GccDPOpt/bin/buoyantPimpleFoam"

[andyreagan]

Okay,

So after re-copying over the 0/* files from the hotRoom tutorial, and manually adding my patches to each file, the case now runs.

However, the Co number quickly explodes and the simulation stops after time 1 or 2. I've reduced the time step, and now I'm looking at improving my mesh.

I'll report back with progress.