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 |
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 |
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 |
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 00:16. |