[shiromaniac]

Hey Foamers,
I am an amateur in CFD and OpenFoam and I have been assigned to do a relatively difficult task( as far as my standards are concerned).
For starters, I am supposed to perform incompressible flow past a heated(having internal heat generation) rectangular cylinder using chtMultiRegionSimpleFoam.
Although, it seemed like a herculean task at first but I somehow managed to mesh both my regions (Domain and Cell) using salome and used ideasUnvToFoam to convert the mesh into OpenFoam readable format.
The issue which I am facing is I just cant get my solver to run. I am attaching my case as well as the errors coming when I run simulation. Please, please please somebody help!


On running Allrun script the following output comes:


/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 5.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 5.x-197d9d3bf20a
Exec : chtMultiRegionSimpleFoam
Date : Jan 31 2019
Time : 10:08:42
Host : "openfoam-VirtualBox"
PID : 18147
I/O : uncollated
Case : /home/openfoam/Desktop/planewall/planeWall2Docomp
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

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

Create fluid mesh for region Domain for time = 0

Create solid mesh for region Cell for time = 0

*** Reading fluid mesh thermophysical properties for region Domain

Adding to thermoFluid

Selecting thermodynamics package
{
type heRhoThermo;
mixture pureMixture;
transport const;
thermo hConst;
equationOfState perfectGas;
specie specie;
energy sensibleEnthalpy;
}

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::heRhoThermo<Foam::rhoThermo, Foam:ureMixture<Foam::constTransport<Foam::speci es::thermo<Foam::hConstThermo<Foam:erfectGas<Foa m::specie> >, Foam::sensibleEnthalpy> > > >::calculate() at ??:?
#4 Foam::rhoThermo::addfvMeshConstructorToTable<Foam: :heRhoThermo<Foam::rhoThermo, Foam:ureMixture<Foam::constTransport<Foam::speci es::thermo<Foam::hConstThermo<Foam:erfectGas<Foa m::specie> >, Foam::sensibleEnthalpy> > > > >::New(Foam::fvMesh const&, Foam::word const&) at ??:?
#5 Foam::autoPtr<Foam::rhoThermo> Foam::basicThermo::New<Foam::rhoThermo>(Foam::fvMe sh const&, Foam::word const&) at ??:?
#6 Foam::rhoThermo::New(Foam::fvMesh const&, Foam::word const&) at ??:?
#7 ? at ??:?
#8 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9 ? at ??:?
Floating point exception (core dumped)



I would have shared my case but its size is exceeding the allowable size limit.
Please Help!

[Bloerb]

Well this sadly tells us nothing. It crashes during setup hence something fundamental is wrong. Delete everything inside all of the polyMesh folders except the boundary file and upload it.

[shiromaniac]

Thanks for your quick reply. I did what you said but it turned out that it is my mesh which is large and even after compressing it comes down to 1.4MB. I am attaching my case here. As far as mesh is concerned could you please give me your email ID so that I can mail it to you. Thanks in advance.

[Bloerb]

There are regions in there which are not solved. Files which set up the regions without them being used etc. Please clean this up, that makes it difficult to tell what is used and what isn't. The only files which should add to the file size are
A the files inside polyMesh except the boundary file and the cellToRegion files. Please add the boundary files since without them it is difficult to tell what your mesh boundary coupling looks like

[shiromaniac]

Hey Bloerb,
Sorry for my late reply as I could not log in due to some personal issues. I am attaching the updated case folder( very sorry as I dint post the case properly last time). Kindly help.





When I run the case following error comes:


bash: source: /opt/openfoam5/bin/: is a directory
openfoam@openfoam-VirtualBox:~/Desktop/planewall/FlowPastCell1$ ./Allclean
Cleaning /home/openfoam/Desktop/planewall/FlowPastCell1 case
openfoam@openfoam-VirtualBox:~/Desktop/planewall/FlowPastCell1$ ./Allrun
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 5.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 5.x-197d9d3bf20a
Exec : ideasUnvToFoam Mesh_nointface1.unv
Date : Feb 04 2019
Time : 16:50:35
Host : "openfoam-VirtualBox"
PID : 3416
I/O : uncollated
Case : /home/openfoam/Desktop/planewall/FlowPastCell1
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

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

Processing tag:164
Starting reading units at line 3.
l:1
units:" SI: Meter (newton)"
unitType:2
Unit factors:
Length scale : 1
Force scale : 1
Temperature scale : 1
Temperature offset : 273.15


Processing tag:2420
Skipping tag 2420 on line 9
Skipping section at line 9.

Processing tag:2411
Starting reading points at line 20.
Read 10421 points.

Processing tag:2412
Starting reading cells at line 20865.
First occurrence of element type 11 for cell 1 at line 20866
First occurrence of element type 41 for cell 265 at line 21658
First occurrence of element type 111 for cell 5325 at line 31778
Read 60088 cells and 5060 boundary faces.

Processing tag:2467
Starting reading patches at line 151956.
For group 1 named Cell trying to read 17520 patch face indices.
For group 2 named Inlet trying to read 84 patch face indices.
For group 3 named Domain trying to read 42568 patch face indices.
For group 4 named Outlet trying to read 84 patch face indices.
For group 5 named Walls trying to read 416 patch face indices.
For group 6 named CellSurf trying to read 4476 patch face indices.

Of 5060 so-called boundary faces 4476 belong to two cells and are therefore internal
Sorting boundary faces according to group (patch)
0: Cell is cellZone
1: Inlet is patch
2: Domain is cellZone
3: Outlet is patch
4: Walls is patch
5: CellSurf is faceZone

Constructing mesh with non-default patches of size:
Inlet 84
Outlet 84
Walls 416

Adding cell and face zones
Cell Zone Domain 42568
Cell Zone Cell 17520
Face Zone CellSurf 4476

End

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 5.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 5.x-197d9d3bf20a
Exec : splitMeshRegions -cellZonesOnly -overwrite
Date : Feb 04 2019
Time : 16:50:41
Host : "openfoam-VirtualBox"
PID : 3418
I/O : uncollated
Case : /home/openfoam/Desktop/planewall/FlowPastCell1
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

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

Create mesh for time = 0

Creating single patch per inter-region interface.

Using current cellZones to split mesh into regions. This requires all cells to be in one and only one cellZone.


Number of regions:2

Writing region per cell file (for manual decomposition) to "/home/openfoam/Desktop/planewall/FlowPastCell1/constant/cellToRegion"

Writing region per cell as volScalarField to "/home/openfoam/Desktop/planewall/FlowPastCell1/0/cellToRegion"

Region Cells
------ -----
0 42568
1 17520

Region Zone Name
------ ---- ----
0 0 Domain
1 1 Cell

Sizes of interfaces between regions:

Interface Region Region Faces
--------- ------ ------ -----
0 0 1 4476

Reading volScalarField cellToRegion
Reading volScalarField p
Reading volScalarField T
Reading volScalarField k
Reading volScalarField epsilon
Reading volScalarField p_rgh
Reading volVectorField U


Adding patches


Adding patches

For interface between region Domain and Cell added patches
3 Domain_to_Cell
4 Cell_to_Domain

Region 0
--------
Creating mesh for region 0 Domain
Testing:"/home/openfoam/Desktop/planewall/FlowPastCell1/system/Domain/fvSchemes"
Mapping fields
Mapping field cellToRegion
Mapping field k
Mapping field p_rgh
Mapping field p
Mapping field T
Mapping field epsilon
Mapping field U
Deleting empty patches
Writing new mesh
Writing addressing to base mesh
Writing map pointRegionAddressing from region0 points back to base mesh.
Writing map faceRegionAddressing from region0 faces back to base mesh.
Writing map cellRegionAddressing from region0 cells back to base mesh.
Writing map boundaryRegionAddressing from region0 boundary back to base mesh.

Region 1
--------
Creating mesh for region 1 Cell
Testing:"/home/openfoam/Desktop/planewall/FlowPastCell1/system/Cell/fvSchemes"
Mapping fields
Mapping field cellToRegion
Mapping field k
Mapping field p_rgh
Mapping field p
Mapping field T
Mapping field epsilon
Mapping field U
Deleting empty patches
Writing new mesh
Writing addressing to base mesh
Writing map pointRegionAddressing from region1 points back to base mesh.
Writing map faceRegionAddressing from region1 faces back to base mesh.
Writing map cellRegionAddressing from region1 cells back to base mesh.
Writing map boundaryRegionAddressing from region1 boundary back to base mesh.
End

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 5.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 5.x-197d9d3bf20a
Exec : chtMultiRegionSimpleFoam
Date : Feb 04 2019
Time : 16:50:44
Host : "openfoam-VirtualBox"
PID : 3423
I/O : uncollated
Case : /home/openfoam/Desktop/planewall/FlowPastCell1
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

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

Create fluid mesh for region Domain for time = 0

Create solid mesh for region Cell for time = 0

*** Reading fluid mesh thermophysical properties for region Domain

Adding to thermoFluid

Selecting thermodynamics package
{
type heRhoThermo;
mixture pureMixture;
transport const;
thermo hConst;
equationOfState rhoConst;
specie specie;
energy sensibleEnthalpy;
}

Adding to rhoFluid

Adding to UFluid

Adding to phiFluid

Adding to gFluid

Adding to hRefFluid

Adding to ghFluid

Adding to ghfFluid

Adding to turbulence

Selecting turbulence model type laminar
Selecting laminar stress model Stokes
Selecting radiationModel none
Adding MRF

No MRF models present

Adding fvOptions

No finite volume options present

*** Reading solid mesh thermophysical properties for region Cell

Adding to thermos

Selecting thermodynamics package
{
type heSolidThermo;
mixture pureMixture;
transport constIso;
thermo hConst;
equationOfState rhoConst;
specie specie;
energy sensibleEnthalpy;
}

Adding to radiations

Selecting radiationModel none
Adding fvOptions

No finite volume options present

Time = 1


Solving for fluid region Domain


--> FOAM FATAL ERROR:
cannot be called for a calculatedFvPatchField
on patch Walls of field U in file "/home/openfoam/Desktop/planewall/FlowPastCell1/0/Domain/U"
You are probably trying to solve for a field with a default boundary condition.

From function Foam::tmp<Foam::Field<Type> > Foam::calculatedFvPatchField<Type>::gradientIntern alCoeffs() const [with Type = Foam::Vector<double>]
in file fields/fvPatchFields/basic/calculated/calculatedFvPatchField.C at line 187.

FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::error::abort() at ??:?
#2 Foam::calculatedFvPatchField<Foam::Vector<double> >::gradientInternalCoeffs() const at ??:?
#3 Foam::fv::gaussLaplacianScheme<Foam::Vector<double >, Foam::SymmTensor<double> >::fvmLaplacianUncorrected(Foam::GeometricField<do uble, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#4 Foam::fv::gaussLaplacianScheme<Foam::Vector<double >, double>::fvmLaplacian(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#5 Foam::fv::laplacianScheme<Foam::Vector<double>, double>::fvmLaplacian(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#6 Foam::tmp<Foam::fvMatrix<Foam::Vector<double> > > Foam::fvm::laplacian<Foam::Vector<double>, double>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) at ??:?
#7 Foam::tmp<Foam::fvMatrix<Foam::Vector<double> > > Foam::fvm::laplacian<Foam::Vector<double>, double>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#8 Foam::tmp<Foam::fvMatrix<Foam::Vector<double> > > Foam::fvm::laplacian<Foam::Vector<double>, double>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#9 Foam::linearViscousStress<Foam::laminarModel<Foam: :ThermalDiffusivity<Foam::CompressibleTurbulenceMo del<Foam::fluidThermo> > > >::divDevRhoReff(Foam::GeometricField<Foam::Vector <double>, Foam::fvPatchField, Foam::volMesh>&) const at ??:?
#10 ? at ??:?
#11 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#12 ? at ??:?
Aborted (core dumped)
openfoam@openfoam-VirtualBox:~/Desktop/planewall/FlowPastCell1$

[Bloerb]

Nearly all your boundary conditions are wrong.

1. You can not impose the velocity on inlet AND outlet. Set the outlet to zeroGradient!
2. Only specify p_rgh and set p to calculated everywhere
3. Change the inlet and outlet boundary conditions for k and epsilon
4. and you need to use changeDictionary afer splitting the mesh.

[shiromaniac]

Hey Bloerb,
I am really thankful to you for having a look at my case. Because of you I was finally able to run it by keeping the Cell(the Solid) at a constant temperature of 400K. Now that i have run the case, my supervisor has asked me to add a constant heat generation source inside the cell. Could you please help me out with that??

Regards,
Shiromani

[Bloerb]

Add an fvOption file to your region folder. Either constant/cell oder system/cell depending on your OF version.

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  dev                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system/cell";
    object      fvOptions;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

options
{
    energySource
    {
        type            scalarSemiImplicitSource;

        selectionMode   all; // or only part of the domain via a cellZone e.g defined with topoSet

        volumeMode      specific; // source term per volume or total
        injectionRateSuSp
        {
            h      (1.5e6 0);   // explicit and implicit part of the source term
        }
    }
}

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

This adds a temperature independent source term in either W or W/m³ depending on the volumeMode. The second number in the injectionRateSuSp(here 0) would be the implicit part and hence the temperature dependent one.

[shiromaniac]

Hey Bloerb,
Thank you so much for all the knowledge you shared with me. I have one more doubt. What should I keep my interface boundary conditions for the region Cell? As of now I have the following interface conditions:
*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.x |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object changeDictionaryDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

T
{
internalField uniform 300;

boundaryField
{

Cell_to_Domain
{
type compressible::turbulentTemperatureCoupledBaffleMix ed;
Tnbr T;
kappaMethod solidThermo;
kappaName none;
value uniform 300;
}

}
}

[Anushaa]

Hi , I am trying to solve my simulation and I am stuck with this ...pls let me know what is the issue

Create time

Create fluid mesh for region waterchannel for time = 0

Create solid mesh for region baseplate for time = 0

*** Reading fluid mesh thermophysical properties for region waterchannel

Adding to thermoFluid

Selecting thermodynamics package
{
type heRhoThermo;
mixture pureMixture;
transport const;
thermo hConst;
equationOfState perfectGas;
specie specie;
energy sensibleEnthalpy;
}

Adding to rhoFluid

Adding to UFluid

Adding to phiFluid

Adding to hRefFluid

Adding to ghFluid

Adding to ghfFluid

Adding to turbulence

Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
Selecting radiationModel none
Adding MRF

No MRF models present

Adding fvOptions

*** Reading solid mesh thermophysical properties for region baseplate

Adding to thermos

Selecting thermodynamics package
{
type heSolidThermo;
mixture pureMixture;
transport constIso;
thermo hConst;
equationOfState rhoConst;
specie specie;
energy sensibleEnthalpy;
}

Adding to radiations

Selecting radiationModel opaqueSolid
Selecting absorptionEmissionModel constantAbsorptionEmission
Selecting scatterModel none
Adding fvOptions

Time = 1


Solving for fluid region waterchannel

[lonewanderer]

Hi Anusha,
what are your boundary conditions?

[amol_patel]

Hi lonewanderer,

I have a similar problem while solving my case of shell and tube heat exchanger
Can you kindly check my thread "chtMultiRegionFoam solver stops without any error"
And if possible, help me please.

Thanks ,
amol.