CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   Problem setting with chtmultiregionFoam (http://www.cfd-online.com/Forums/openfoam/64455-problem-setting-chtmultiregionfoam.html)

Antonin May 12, 2009 04:51

Problem setting with chtmultiregionFoam
 
1 Attachment(s)
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

mabinty May 16, 2009 11:38

Antonin!

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

Aram

Antonin May 18, 2009 04:09

1 Attachment(s)
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

Antonin May 18, 2009 12:28

1 Attachment(s)
Attachment 312
I add this files for more comprehension

tH3f0rC3 April 12, 2011 01:50

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

mabinty April 12, 2011 02:35

hi!

you might have set

Code:

adjustTimeStep  yes;
maxCo          0.3;

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

Code:

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

cheers,
aram

tH3f0rC3 April 12, 2011 06:53

Perfect!
That was exactly the point.

:-)

Thank you!

Budlo April 18, 2012 05:36

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

Budlo April 18, 2012 05:38

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

calim_cfd April 18, 2012 06:42

Quote:

Originally Posted by Budlo (Post 355330)
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

Budlo April 24, 2012 09:50

Thanks a lot "Calim"


All times are GMT -4. The time now is 20:53.