|
[Sponsors] |
April 22, 2009, 06:26 |
Problem setting with chtmultiregionFoam
|
#1 |
New Member
Antonin
Join Date: Apr 2009
Location: France
Posts: 7
Rep Power: 17 |
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 13:16. |
|
July 2, 2009, 06:09 |
|
#2 |
New Member
Kyian Barrat
Join Date: Apr 2009
Posts: 25
Rep Power: 17 |
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 06:07. |
|
July 20, 2009, 06:12 |
Error message chtMultiRegionFoam
|
#3 |
New Member
Kyian Barrat
Join Date: Apr 2009
Posts: 25
Rep Power: 17 |
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 |
|
July 20, 2009, 07:52 |
|
#4 |
Senior Member
Henrik Rusche
Join Date: Mar 2009
Location: Wernigerode, Sachsen-Anhalt, Germany
Posts: 281
Rep Power: 18 |
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 |
|
Tags |
antonin |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Problem setting C_MU_T with DEFINE_ADJUST | Andy | FLUENT | 3 | June 22, 2018 03:53 |
HELP!: swimming pool flow problem | Tee | FLUENT | 1 | September 8, 2005 09: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 |