CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Problem setting with chtmultiregionFoam

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   April 22, 2009, 05:26
Lightbulb Problem setting with chtmultiregionFoam
  #1
New Member
 
Antonin
Join Date: Apr 2009
Location: France
Posts: 7
Rep Power: 8
Antonin is on a distinguished road
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:rintStack(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:imensionedField<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:imensionedField<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:imensionedField<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

Last edited by Antonin; May 18, 2009 at 12:16.
Antonin is offline   Reply With Quote

Old   July 2, 2009, 05:09
Default
  #2
New Member
 
Kyian Barrat
Join Date: Apr 2009
Posts: 25
Rep Power: 8
Khelian973 is on a distinguished road
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:rintStack(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:perator/<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

Last edited by Khelian973; July 20, 2009 at 05:07.
Khelian973 is offline   Reply With Quote

Old   July 20, 2009, 05:12
Default Error message chtMultiRegionFoam
  #3
New Member
 
Kyian Barrat
Join Date: Apr 2009
Posts: 25
Rep Power: 8
Khelian973 is on a distinguished road
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:rintStack(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:ureMixture<Foam::sutherlandTr ansport<Foam::specieThermo<Foam::hConstThermo<Foam :erfectGas> > > > >::calculate() in "/DATA/Softs/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libbasicThermophysicalModels.so"

#4 Foam::hThermo<Foam:ureMixture<Foam::sutherlandTr ansport<Foam::specieThermo<Foam::hConstThermo<Foam :erfectGas> > > > >::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
Khelian973 is offline   Reply With Quote

Old   July 20, 2009, 06:52
Default
  #4
Senior Member
 
Henrik Rusche
Join Date: Mar 2009
Location: Braunschweig, Niedersachsen, Germany
Posts: 275
Rep Power: 9
henrik is on a distinguished road
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
henrik is offline   Reply With Quote

Reply

Tags
antonin

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Problem setting C_MU_T with DEFINE_ADJUST Andy FLUENT 2 July 30, 2014 05:05
HELP!: swimming pool flow problem Tee FLUENT 1 September 8, 2005 08:29
Problem setting AV Marko CFX 0 March 18, 2005 09:06
Problem setting C_MU_T with DEFINE_ADJUST sarah_ron FLUENT 1 February 24, 2005 13:11
Help setting up combustor problem Tom FLUENT 0 February 5, 2005 02:20


All times are GMT -4. The time now is 14:13.