Problem setting with chtmultiregionFoam
1 Attachment(s)
Hello,
I am using openfoam 1.5 on ubuntu (installed from ppa repo)and I'm novice in OpenFoam 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.cfdonline.com/OpenFOAM_D...eater8318.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/OpenFOAM1.5/lib/libOpenFOAM.so" #1 Foam::sigSegv::sigSegvHandler(int) in "/usr/lib/OpenFOAM1.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/OpenFOAM1.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/OpenFOAM1.5/applications/bin/chtMultiRegionFoam" #5 Foam::fixedGradientFvPatchField<double>::evaluate( Foam::Pstream::commsTypes) in "/usr/lib/OpenFOAM1.5/applications/bin/chtMultiRegionFoam" #6 Foam::gradientEnthalpyFvPatchScalarField::updateCo effs() in "/usr/lib/OpenFOAM1.5/lib/libbasicThermophysicalModels.so" #7 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricBoundaryField::updateCoef fs() in "/usr/lib/OpenFOAM1.5/lib/libfiniteVolume.so" #8 Foam::fvMatrix<double>::fvMatrix(Foam::GeometricFi eld<double, Foam::fvPatchField, Foam::volMesh>&, Foam::dimensionSet const&) in "/usr/lib/OpenFOAM1.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/OpenFOAM1.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/OpenFOAM1.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/OpenFOAM1.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/OpenFOAM1.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/OpenFOAM1.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/OpenFOAM1.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/OpenFOAM1.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/OpenFOAM1.5/applications/bin/chtMultiRegionFoam" #17 __libc_start_main in "/lib/libc.so.6" #18 ?? in "/usr/lib/OpenFOAM1.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 Attachment 273).If you have any solutions please don't esitate. Thank you. Antonin 
Antonin!
One needs to know more about your setup (e.g. BCs). Unfortunately I was not able to download the system files. Aram 
1 Attachment(s)
Aram
Thank you, for the speed of the reponce. I joined my post my file 0 Attachment 310 , if you want more information, tell me. I hope that this, could shed light on things. Antonin 
1 Attachment(s)
Attachment 312
I add this files for more comprehension 
Hi,
I want to know if it is right that a transient solver solves like this: // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create fluid mesh for region air for time = 0 Create solid mesh for region halterungsgestell for time = 0 *** Reading fluid mesh thermophysical properties for region air 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 RASModel Selecting RAS turbulence model kEpsilon kEpsilonCoeffs { Cmu 0.09; C1 1.44; C2 1.92; C3 0.33; sigmak 1; sigmaEps 1.3; Prt 1; } Adding to radiation Selecting radiationModel P1 Selecting absorptionEmissionModel constantAbsorptionEmission Selecting scatterModel constantScatter Adding to DpDtFluid *** Reading solid mesh thermophysical properties for region halterungsgestell Adding to rhos Adding to cps Adding to Ks Adding to Ts Region: air Courant Number mean: 0 max: 6686.693 Region: air Courant Number mean: 0 max: 0.29999969 deltaT = 4.486518e05 Time = 4.48652e05 Solving for fluid region air diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 4.5910077e08, No Iterations 6 DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 5.0421213e08, No Iterations 4 DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 8.4322666e08, No Iterations 10 DICPCG: Solving for G, Initial residual = 1, Final residual = 0.090942605, No Iterations 37 DILUPBiCG: Solving for h, Initial residual = 0.99999999, Final residual = 9.7740028e08, No Iterations 12 Min/max T:293.02624 973.15 GAMG: Solving for p, Initial residual = 1, Final residual = 0.0021243927, No Iterations 2 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors (air): sum local = 1.3361364e07, global = 9.2719095e09, cumulative = 9.2719095e09 GAMG: Solving for p, Initial residual = 0.16696896, Final residual = 0.0001442381, No Iterations 2 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors (air): sum local = 1.9117031e08, global = 1.0869296e09, cumulative = 8.1849799e09 DILUPBiCG: Solving for epsilon, Initial residual = 0.032800783, Final residual = 9.1607404e08, No Iterations 10 DILUPBiCG: Solving for k, Initial residual = 1, Final residual = 6.8377323e08, No Iterations 8 Building global boundary list Mass flux at Kettengetriebe = 0 Mass flux at Tuer = 0 Mass flux at Inlet_Outlet = 3.2890796 Mass flux at Einduesungsanlage = 0 Mass flux at Dueseneintritt = 3.0569587 Mass flux at Aussenwand = 0 Mass flux at untere_Einduesungsaustritte = 0 Mass flux at Abgasoutlet = 2.9958435e05 Mass flux at Querverstrebungen = 0 Mass flux at Halterungsgestell = 0 Mass flux at Brenner_Einlass = 0 Net mass flux = 0.23215087 Solving for solid region halterungsgestell DICPCG: Solving for T, Initial residual = 1, Final residual = 8.4983523e07, No Iterations 96 DICPCG: Solving for T, Initial residual = 0.21773582, Final residual = 7.8090756e07, No Iterations 85 Min/max T:min(T) [0 0 0 1 0 0 0] 572.93858 max(T) [0 0 0 1 0 0 0] 572.97669 ExecutionTime = 14.96 s ClockTime = 16 s Region: air Courant Number mean: 4.6333806e05 max: 0.29998597 deltaT = 4.4867193e05 Time = 8.97324e05 Solving for fluid region air diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for Ux, Initial residual = 0.37506564, Final residual = 6.6027038e08, No Iterations 10 DILUPBiCG: Solving for Uy, Initial residual = 0.48669134, Final residual = 4.2404301e08, No Iterations 13 DILUPBiCG: Solving for Uz, Initial residual = 0.59158917, Final residual = 8.5674749e08, No Iterations 12 DICPCG: Solving for G, Initial residual = 0.42589132, Final residual = 0.041371532, No Iterations 3 DILUPBiCG: Solving for h, Initial residual = 0.64209013, Final residual = 9.6251018e08, No Iterations 13 Min/max T:561.21 973.15 GAMG: Solving for p, Initial residual = 0.45521647, Final residual = 0.00046735382, No Iterations 2 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors (air): sum local = 9.5623605e08, global = 9.3885069e09, cumulative = 1.2035271e09 GAMG: Solving for p, Initial residual = 0.060360581, Final residual = 5.0485323e05, No Iterations 2 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors (air): sum local = 1.4832474e08, global = 3.1777063e10, cumulative = 8.8575644e10 DILUPBiCG: Solving for epsilon, Initial residual = 0.0034727837, Final residual = 3.0269126e08, No Iterations 8 DILUPBiCG: Solving for k, Initial residual = 0.40375756, Final residual = 4.1736344e08, No Iterations 13 Building global boundary list Mass flux at Kettengetriebe = 0 Mass flux at Tuer = 0 Mass flux at Inlet_Outlet = 3.2784506 Mass flux at Einduesungsanlage = 0 Mass flux at Dueseneintritt = 3.1253266 Mass flux at Aussenwand = 0 Mass flux at untere_Einduesungsaustritte = 0 Mass flux at Abgasoutlet = 2.9972721e05 Mass flux at Querverstrebungen = 0 Mass flux at Halterungsgestell = 0 Mass flux at Brenner_Einlass = 0 Net mass flux = 0.15315391 Solving for solid region halterungsgestell DICPCG: Solving for T, Initial residual = 0.24529809, Final residual = 9.7977311e07, No Iterations 72 DICPCG: Solving for T, Initial residual = 0.079959901, Final residual = 9.5937062e07, No Iterations 61 Min/max T:min(T) [0 0 0 1 0 0 0] 572.92759 max(T) [0 0 0 1 0 0 0] 572.97414 ExecutionTime = 27.23 s ClockTime = 29 s Region: air Courant Number mean: 7.9389894e05 max: 0.30171661 deltaT = 4.4611786e05 Time = 0.000134344 Solving for fluid region air diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 __________________________________________________ ____________________ I miss the more clear structure of the steadystate solvers like chtMultiSimpleFoam. I have entered into controlDict as a time step 1. But the chtMultiRadFoam doesn't solve as instructed. Does someone can explain that to me? Many thanks! Best regards, tH3f0rC3 
hi!
you might have set Code:
adjustTimeStep yes; Code:
adjustTimeStep no; cheers, aram 
Perfect!
That was exactly the point. :) Thank you! 
chtMultiRegionFoam
Hi
I want take a simple run from example "chtMultiRegionFoam" but when It give error. Error is " Can not read file "pints" at time 0 ". How do you take a simple run from This solver (after blockMesh)? Thanks 
chtMultiRegionFoam
Hi
I want take a simple run from example "chtMultiRegionFoam" but when It give error. Error is " Can not read file "pints" at time 0 ". How do you take a simple run from This solver (after blockMesh)? Thanks 
Quote:
take a better look at the case files and folders at /tutorials/heatTransfer/chtMultiRegion* but run these cases b4 taking any conclusions. The Allrun script will do a few things. Once the case is run, take a look at the file/folder structure you got and the log files. cht makes sense with more than 1 volume region, (solid and air,... etc) that you'll have to create using blockmesh or splitMeshregions app. Unless your case is simple it might be a little hard to set the case using OF preprocessing tools alone, but it is feasible. gl 
Thanks a lot "Calim"

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