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

Help with chtMultiRegionFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   December 21, 2010, 13:07
Default Help with chtMultiRegionFoam
  #1
New Member
 
Jonathan Brodie
Join Date: Dec 2009
Posts: 4
Rep Power: 6
jbvw96 is on a distinguished road
Hello,

I am having trouble diagnosing the issue I am having getting CHTmultirgeion foam running. I have had a few of these Floating Point exception errors now and I have been able to find mistakes in the boundary conditions. But now I'm stuck, here is the output: where should I be looking?

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6.x |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 1.6.x-05e66c6c0080
Exec : chtMultiRegionFoam
Date : Dec 21 2010
Time : 08:04:48
Host : jbrodie
PID : 4193
Case : /home/jonathan/OpenFOAM/OpenFOAM-1.6.x/run/snappyShellTube
nProcs : 1
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).

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

Create fluid mesh for region domain0 for time = 0.001

Create fluid mesh for region internal for time = 0.001

Create solid mesh for region tube for time = 0.001

*** Reading fluid mesh thermophysical properties for region domain0

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 laminar
Adding to DpDtFluid

*** Reading fluid mesh thermophysical properties for region internal

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 laminar
Adding to DpDtFluid

*** Reading solid mesh thermophysical properties for region tube

Adding to rhos

Adding to cps

Adding to Ks

Adding to Ts

Region: domain0 Courant Number mean: 0 max: 129.5818
Region: internal Courant Number mean: 0 max: 111.2308
#0 Foam::error:rintStack(Foam::Ostream&) in "/home/jonathan/OpenFOAM/OpenFOAM-1.6.x/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/jonathan/OpenFOAM/OpenFOAM-1.6.x/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Uninterpreted:
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/home/jonathan/OpenFOAM/OpenFOAM-1.6.x/lib/linuxGccDPOpt/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 "/home/jonathan/OpenFOAM/OpenFOAM-1.6.x/applications/bin/linuxGccDPOpt/chtMultiRegionFoam"
#5 Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > Foam:perator/<Foam::fvsPatchField, Foam::surfaceMesh>(Foam::tmp<Foam::GeometricField< double, Foam::fvsPatchField, Foam::surfaceMesh> > const&, Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > const&) in "/home/jonathan/OpenFOAM/OpenFOAM-1.6.x/applications/bin/linuxGccDPOpt/chtMultiRegionFoam"
#6 Foam::solidRegionDiffNo(Foam::fvMesh const&, Foam::Time const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/home/jonathan/OpenFOAM/OpenFOAM-1.6.x/applications/bin/linuxGccDPOpt/chtMultiRegionFoam"
#7 main in "/home/jonathan/OpenFOAM/OpenFOAM-1.6.x/applications/bin/linuxGccDPOpt/chtMultiRegionFoam"
#8 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#9 _start at /build/buildd/eglibc-2.10.1/csu/../sysdeps/i386/elf/start.S:122
Floating point exception



Thanks for any ideas you could contribute
jbvw96 is offline   Reply With Quote

Old   December 22, 2010, 12:43
Default check time step
  #2
New Member
 
Jean El-Hajal
Join Date: Jun 2010
Location: Ulm
Posts: 16
Rep Power: 6
Jean El-Hajal is on a distinguished road
Hi Jonathan,

You probably have to check the Boundary Condition and/or the time step. The courant Number should be < 0.5.

The error massage say Courant number > 100 !

Hope this can help you.

Jean
Jean El-Hajal is offline   Reply With Quote

Old   December 26, 2010, 17:16
Default
  #3
Member
 
Paulo Bufacchi
Join Date: Jun 2010
Posts: 50
Rep Power: 6
Bufacchi is on a distinguished road
Hi Jonathan,
It seems to me you have a problem in your region definition. I see a domain0 in you log file, which shouldn't exist. I understand your regions are internal and tube, so you should use ./Allclean in your case directory to clean it up and then make sure all files and directories reflect only these two regions. Delete all directories you may find in 0, constants and system directories.

Regards,

Paulo
Bufacchi is offline   Reply With Quote

Reply

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
Error in chtMultiRegionFoam kirankarki OpenFOAM 5 October 24, 2012 17:20
Problem with running chtMultiRegionFoam after using setSet utility Victor OpenFOAM 10 September 10, 2012 09:30
Using a Different Thermodynamics Package with the chtMultiRegionFoam Solver m.nichols19 OpenFOAM 7 March 17, 2011 17:26
chtMultiRegionFoam with liquid phsieh2005 OpenFOAM Running, Solving & CFD 0 October 9, 2010 06:07
How to add 3rd region to chtMultiRegionFoam benk OpenFOAM 1 April 3, 2010 17:22


All times are GMT -4. The time now is 19:33.