Problem setting with chtmultiregionFoam
 

Hello,

I am using openfoam 1.5 on ubuntu (installed from ppa repo)and I'm novice in OpenFoam
To get a quick summary: I want to simulate an hot fluid flow in a steel tube.
I am trying to setup a 3D case with chtmultiregionFoam, but I'm getting some trouble with it.I am starting from the tutorial one
(see http://www.cfd-online.com/OpenFOAM_D...eater-8318.unk ).With my new case, i get a segmentation faults.

My problem is, when I run chtmultiregion this message apair:

Create time

Create fluid mesh for region gaz for time = 0.001

Create solid mesh for region Solid for time = 0.001

*** Reading fluid mesh thermophysical properties for region gaz

Adding to pdf

Adding to thermof

Selecting thermodynamics package hThermo<pureMixture<constTransport<specieThermo<hC onstThermo<perfectGas>>>>>
Adding to rhof

Adding to Kf

Adding to Uf

Adding to phif

Adding to turb

Selecting RAS turbulence model kEpsilon
Adding to DpDtf

Adding to ghf

Updating p from pd

*** Reading solid mesh thermophysical properties for region Solid

Adding to rhos

Adding to cps

Adding to Ks

Adding to Ts

Region: gaz Courant Number mean: 0 max: 0.01
Region: gaz Courant Number mean: 0 max: 0.2994012
deltaT = 0.02994012
Time = 0.0309401


Solving for fluid region gaz
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
#0 Foam::error::printStack(Foam::Ostream&) in "/usr/lib/OpenFOAM-1.5/lib/libOpenFOAM.so"
#1 Foam::sigSegv::sigSegvHandler(int) in "/usr/lib/OpenFOAM-1.5/lib/libOpenFOAM.so"
#2 ?? in "/lib/libc.so.6"
#3 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::PatchFieldType const& Foam::fvPatch::lookupPatchField<Foam::GeometricFie ld<double, Foam::fvPatchField, Foam::volMesh>, double>(Foam::word const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const*, double const*) const in "/usr/lib/OpenFOAM-1.5/applications/bin/chtMultiRegionFoam"
#4 Foam::fvPatchField<double>::adddictionaryConstruct orToTable<Foam::solidWallHeatFluxTemperatureFvPatc hScalarField>::New(Foam::fvPatch const&, Foam::DimensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) in "/usr/lib/OpenFOAM-1.5/applications/bin/chtMultiRegionFoam"
#5 Foam::fixedGradientFvPatchField<double>::evaluate( Foam::Pstream::commsTypes) in "/usr/lib/OpenFOAM-1.5/applications/bin/chtMultiRegionFoam"
#6 Foam::gradientEnthalpyFvPatchScalarField::updateCo effs() in "/usr/lib/OpenFOAM-1.5/lib/libbasicThermophysicalModels.so"
#7 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricBoundaryField::updateCoef fs() in "/usr/lib/OpenFOAM-1.5/lib/libfiniteVolume.so"
#8 Foam::fvMatrix<double>::fvMatrix(Foam::GeometricFi eld<double, Foam::fvPatchField, Foam::volMesh>&, Foam::dimensionSet const&) in "/usr/lib/OpenFOAM-1.5/lib/libfiniteVolume.so"
#9 Foam::fv::gaussLaplacianScheme<double, double>::fvmLaplacianUncorrected(Foam::GeometricFi eld<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&) in "/usr/lib/OpenFOAM-1.5/lib/libfiniteVolume.so"
#10 Foam::fv::gaussLaplacianScheme<double, double>::fvmLaplacian(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&) in "/usr/lib/OpenFOAM-1.5/lib/libfiniteVolume.so"
#11 Foam::fv::laplacianScheme<double, double>::fvmLaplacian(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&) in "/usr/lib/OpenFOAM-1.5/lib/libfiniteVolume.so"
#12 Foam::tmp<Foam::fvMatrix<double> > Foam::fvm::laplacian<double, double>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::word const&) in "/usr/lib/OpenFOAM-1.5/applications/bin/chtMultiRegionFoam"
#13 Foam::tmp<Foam::fvMatrix<double> > Foam::fvm::laplacian<double, double>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&) in "/usr/lib/OpenFOAM-1.5/applications/bin/chtMultiRegionFoam"
#14 Foam::tmp<Foam::fvMatrix<double> > Foam::fvm::laplacian<double, double>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&) in "/usr/lib/OpenFOAM-1.5/applications/bin/chtMultiRegionFoam"
#15 Foam::fvPatchField<double>::adddictionaryConstruct orToTable<Foam::solidWallTemperatureCoupledFvPatch ScalarField>::New(Foam::fvPatch const&, Foam::DimensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) in "/usr/lib/OpenFOAM-1.5/applications/bin/chtMultiRegionFoam"
#16 Foam::fvPatchField<double>::adddictionaryConstruct orToTable<Foam::solidWallTemperatureCoupledFvPatch ScalarField>::New(Foam::fvPatch const&, Foam::DimensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) in "/usr/lib/OpenFOAM-1.5/applications/bin/chtMultiRegionFoam"
#17 __libc_start_main in "/lib/libc.so.6"
#18 ?? in "/usr/lib/OpenFOAM-1.5/applications/bin/chtMultiRegionFoam"
Segmentation fault


I do not know where this error come,I think I misunderstood creatPatchDict set in the file system, but I'm not sure because I did not understand its usefulness.
For more information I have joind to my post my file system (system.tar.gz Attachment 273).If you have any solutions please don't esitate.
Thank you.

Antonin

Antonin!

One needs to know more about your set-up (e.g. BCs). Unfortunately I was not able to download the system files.

Aram

Aram

Thank you, for the speed of the reponce.
I joined my post my file 0 Attachment 310 , if you want more information, tell me.
I hope that this, could shed light on things.

Antonin

Attachment 312
I add this files for more comprehension

Hi,

I want to know if it is right that a transient solver solves like this:

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

Create fluid mesh for region air for time = 0

Create solid mesh for region halterungsgestell for time = 0

*** Reading fluid mesh thermophysical properties for region air

Adding to thermoFluid

Selecting thermodynamics package hPsiThermo<pureMixture<constTransport<specieThermo <hConstThermo<perfectGas>>>>>
Adding to rhoFluid

Adding to KFluid

Adding to UFluid

Adding to phiFluid

Adding to gFluid

Adding to turbulence

Selecting turbulence model type RASModel
Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
Cmu 0.09;
C1 1.44;
C2 1.92;
C3 -0.33;
sigmak 1;
sigmaEps 1.3;
Prt 1;
}

Adding to radiation

Selecting radiationModel P1
Selecting absorptionEmissionModel constantAbsorptionEmission
Selecting scatterModel constantScatter
Adding to DpDtFluid

*** Reading solid mesh thermophysical properties for region halterungsgestell

Adding to rhos

Adding to cps

Adding to Ks

Adding to Ts

Region: air Courant Number mean: 0 max: 6686.693
Region: air Courant Number mean: 0 max: 0.29999969
deltaT = 4.486518e-05
Time = 4.48652e-05


Solving for fluid region air
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 4.5910077e-08, No Iterations 6
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 5.0421213e-08, No Iterations 4
DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 8.4322666e-08, No Iterations 10
DICPCG: Solving for G, Initial residual = 1, Final residual = 0.090942605, No Iterations 37
DILUPBiCG: Solving for h, Initial residual = 0.99999999, Final residual = 9.7740028e-08, No Iterations 12
Min/max T:293.02624 973.15
GAMG: Solving for p, Initial residual = 1, Final residual = 0.0021243927, No Iterations 2
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (air): sum local = 1.3361364e-07, global = 9.2719095e-09, cumulative = 9.2719095e-09
GAMG: Solving for p, Initial residual = 0.16696896, Final residual = 0.0001442381, No Iterations 2
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (air): sum local = 1.9117031e-08, global = -1.0869296e-09, cumulative = 8.1849799e-09
DILUPBiCG: Solving for epsilon, Initial residual = 0.032800783, Final residual = 9.1607404e-08, No Iterations 10
DILUPBiCG: Solving for k, Initial residual = 1, Final residual = 6.8377323e-08, No Iterations 8
Building global boundary list
Mass flux at Kettengetriebe = 0
Mass flux at Tuer = 0
Mass flux at Inlet_Outlet = 3.2890796
Mass flux at Einduesungsanlage = 0
Mass flux at Dueseneintritt = -3.0569587
Mass flux at Aussenwand = 0
Mass flux at untere_Einduesungsaustritte = 0
Mass flux at Abgasoutlet = 2.9958435e-05
Mass flux at Querverstrebungen = 0
Mass flux at Halterungsgestell = 0
Mass flux at Brenner_Einlass = 0
Net mass flux = 0.23215087


Solving for solid region halterungsgestell
DICPCG: Solving for T, Initial residual = 1, Final residual = 8.4983523e-07, No Iterations 96
DICPCG: Solving for T, Initial residual = 0.21773582, Final residual = 7.8090756e-07, No Iterations 85
Min/max T:min(T) [0 0 0 1 0 0 0] 572.93858 max(T) [0 0 0 1 0 0 0] 572.97669
ExecutionTime = 14.96 s ClockTime = 16 s

Region: air Courant Number mean: 4.6333806e-05 max: 0.29998597
deltaT = 4.4867193e-05
Time = 8.97324e-05


Solving for fluid region air
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 0.37506564, Final residual = 6.6027038e-08, No Iterations 10
DILUPBiCG: Solving for Uy, Initial residual = 0.48669134, Final residual = 4.2404301e-08, No Iterations 13
DILUPBiCG: Solving for Uz, Initial residual = 0.59158917, Final residual = 8.5674749e-08, No Iterations 12
DICPCG: Solving for G, Initial residual = 0.42589132, Final residual = 0.041371532, No Iterations 3
DILUPBiCG: Solving for h, Initial residual = 0.64209013, Final residual = 9.6251018e-08, No Iterations 13
Min/max T:561.21 973.15
GAMG: Solving for p, Initial residual = 0.45521647, Final residual = 0.00046735382, No Iterations 2
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (air): sum local = 9.5623605e-08, global = -9.3885069e-09, cumulative = -1.2035271e-09
GAMG: Solving for p, Initial residual = 0.060360581, Final residual = 5.0485323e-05, No Iterations 2
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (air): sum local = 1.4832474e-08, global = 3.1777063e-10, cumulative = -8.8575644e-10
DILUPBiCG: Solving for epsilon, Initial residual = 0.0034727837, Final residual = 3.0269126e-08, No Iterations 8
DILUPBiCG: Solving for k, Initial residual = 0.40375756, Final residual = 4.1736344e-08, No Iterations 13
Building global boundary list
Mass flux at Kettengetriebe = 0
Mass flux at Tuer = 0
Mass flux at Inlet_Outlet = 3.2784506
Mass flux at Einduesungsanlage = 0
Mass flux at Dueseneintritt = -3.1253266
Mass flux at Aussenwand = 0
Mass flux at untere_Einduesungsaustritte = 0
Mass flux at Abgasoutlet = 2.9972721e-05
Mass flux at Querverstrebungen = 0
Mass flux at Halterungsgestell = 0
Mass flux at Brenner_Einlass = 0
Net mass flux = 0.15315391


Solving for solid region halterungsgestell
DICPCG: Solving for T, Initial residual = 0.24529809, Final residual = 9.7977311e-07, No Iterations 72
DICPCG: Solving for T, Initial residual = 0.079959901, Final residual = 9.5937062e-07, No Iterations 61
Min/max T:min(T) [0 0 0 1 0 0 0] 572.92759 max(T) [0 0 0 1 0 0 0] 572.97414
ExecutionTime = 27.23 s ClockTime = 29 s

Region: air Courant Number mean: 7.9389894e-05 max: 0.30171661
deltaT = 4.4611786e-05
Time = 0.000134344


Solving for fluid region air
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
__________________________________________________ ____________________
I miss the more clear structure of the steady-state solvers like chtMultiSimpleFoam.
I have entered into controlDict as a time step 1. But the chtMultiRadFoam doesn't solve as instructed.

Does someone can explain that to me?
Many thanks!

Best regards,
tH3f0rC3

hi!

you might have set

in the controclDict, i.e. the time step is controlled by a constant Courant number of 0.3. set

then the prescribed time step "deltaT" is used by the solver.

cheers,
aram

Perfect!
That was exactly the point.

:-)

Thank you!

chtMultiRegionFoam
 

Hi
I want take a simple run from example "chtMultiRegionFoam" but when It give error. Error is " Can not read file "pints" at time 0 ".
How do you take a simple run from This solver (after blockMesh)?
Thanks

chtMultiRegionFoam
 

Hi
I want take a simple run from example "chtMultiRegionFoam" but when It give error. Error is " Can not read file "pints" at time 0 ".
How do you take a simple run from This solver (after blockMesh)?
Thanks

Hi!
take a better look at the case files and folders at
/tutorials/heatTransfer/chtMultiRegion*

but run these cases b4 taking any conclusions. The Allrun script will do a few things. Once the case is run, take a look at the file/folder structure you got and the log files.

cht makes sense with more than 1 volume region, (solid and air,... etc) that you'll have to create using blockmesh or splitMeshregions app. Unless your case is simple it might be a little hard to set the case using OF pre-processing tools alone, but it is feasible.
gl

Thanks a lot "Calim"