CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Problem setting with chtmultiregionFoam (http://www.cfd-online.com/Forums/openfoam-solving/63860-problem-setting-chtmultiregionfoam.html)

Antonin April 22, 2009 05:26

Problem setting with chtmultiregionFoam
 
Hello,

I am using openfoam 1.5 on ubuntu (installed from ppa repo)
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).If you have any solutions please don't esitate.
Thank you.

Antonin

Khelian973 July 2, 2009 05:09

Hi Antonin, Foamers

I have also a error message with chtMultiRegionFoam. It's not a segmentation fault, but maybe someone can help me or has an idea.

Create time
Create fluid mesh for region water for time = 0
Create solid mesh for region alu for time = 0
*** Reading fluid mesh thermophysical properties for region water
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 laminar
Adding to DpDtf

#0 Foam::error::printStack(Foam::Ostream&) in "/DATA/Softs/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"

#1 Foam::sigFpe::sigFpeHandler(int) in "/DATA/Softs/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"

#2 ?? in "/lib/libc.so.6"

#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/DATA/Softs/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"

#4 void Foam::divide<Foam::fvsPatchField, Foam::surfaceMesh>(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/DATA/Softs/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/chtMultiRegionFoam"

#5 Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > Foam::operator/<Foam::fvsPatchField, Foam::surfaceMesh>(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > const&) in "/DATA/Softs/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/chtMultiRegionFoam"

#6 main in "/DATA/Softs/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/chtMultiRegionFoam"

#7 __libc_start_main in "/lib/libc.so.6"

#8 Foam::regIOobject::readIfModified() in "/DATA/Softs/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/chtMultiRegionFoam"

Floating point exception


Thanks for any Help
Kyian

Khelian973 July 20, 2009 05:12

Error message chtMultiRegionFoam
 
Hi Foamers,

I succeed to solve my problem to build the case and the different regions. But now i have another problem concerning the solver. Indeed, the calculations stop after around 10 iterations and I have this message

Solving for fluid region gn2
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for h, Initial residual = 0.9297533, Final residual = 8.310485e-12, No Iterations 1

#0 Foam::error::printStack(Foam::Ostream&) in "/DATA/Softs/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"

#1 Foam::sigFpe::sigFpeHandler(int) in "/DATA/Softs/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"

#2 ?? in "/lib/libc.so.6"

#3 Foam::hThermo<Foam::pureMixture<Foam::sutherlandTr ansport<Foam::specieThermo<Foam::hConstThermo<Foam ::perfectGas> > > > >::calculate() in "/DATA/Softs/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libbasicThermophysicalModels.so"

#4 Foam::hThermo<Foam::pureMixture<Foam::sutherlandTr ansport<Foam::specieThermo<Foam::hConstThermo<Foam ::perfectGas> > > > >::correct() in "/DATA/Softs/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libbasicThermophysicalModels.so"

#5 solveEnthalpyEquation(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::compressible::RASModel const&, Foam::basicThermo&) in "/DATA/Softs/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/chtMultiRegionFoam"

#6 main in "/DATA/Softs/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/chtMultiRegionFoam"

#7 __libc_start_main in "/lib/libc.so.6"

#8 Foam::regIOobject::readIfModified() in "/DATA/Softs/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/chtMultiRegionFoam"

Floating point exception


I think it's probably due to the thermophysical properties, particularly sutherland Transport coefficient. But i double checked my figures and they are good.
If you have an idea, please


Kyian

henrik July 20, 2009 06:52

Dear Khelian973,

my bet is that the temperature is going below zero and then the sqrt(T) in the sutherlandTransport blows up.

Check your boundary and initial conditions and if correct try a smaller time step, more coupling until you you get a first time dump which you can inspect. Note that the regions are explicitly coupled - so please be gentle.

Henrik


All times are GMT -4. The time now is 02:01.