|
[Sponsors] |
May 12, 2009, 04:51 |
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)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.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 system.tar.gz).If you have any solutions please don't esitate. Thank you. Antonin |
|
May 16, 2009, 11:38 |
|
#2 |
Senior Member
Aram Amouzandeh
Join Date: Mar 2009
Location: Vienna, Vienna, Austria
Posts: 190
Rep Power: 17 |
Antonin!
One needs to know more about your set-up (e.g. BCs). Unfortunately I was not able to download the system files. Aram |
|
May 18, 2009, 12:28 |
|
#4 |
New Member
Antonin
Join Date: Apr 2009
Location: France
Posts: 7
Rep Power: 17 |
chtMultiregion.tar.gz
I add this files for more comprehension |
|
April 12, 2011, 01:50 |
|
#5 |
Senior Member
Join Date: Mar 2011
Posts: 158
Rep Power: 15 |
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.486518e-05 Time = 4.48652e-05 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.5910077e-08, No Iterations 6 DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 5.0421213e-08, No Iterations 4 DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 8.4322666e-08, 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.7740028e-08, 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.3361364e-07, global = 9.2719095e-09, cumulative = 9.2719095e-09 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.9117031e-08, global = -1.0869296e-09, cumulative = 8.1849799e-09 DILUPBiCG: Solving for epsilon, Initial residual = 0.032800783, Final residual = 9.1607404e-08, No Iterations 10 DILUPBiCG: Solving for k, Initial residual = 1, Final residual = 6.8377323e-08, 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.9958435e-05 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.4983523e-07, No Iterations 96 DICPCG: Solving for T, Initial residual = 0.21773582, Final residual = 7.8090756e-07, 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.6333806e-05 max: 0.29998597 deltaT = 4.4867193e-05 Time = 8.97324e-05 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.6027038e-08, No Iterations 10 DILUPBiCG: Solving for Uy, Initial residual = 0.48669134, Final residual = 4.2404301e-08, No Iterations 13 DILUPBiCG: Solving for Uz, Initial residual = 0.59158917, Final residual = 8.5674749e-08, 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.6251018e-08, 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.5623605e-08, global = -9.3885069e-09, cumulative = -1.2035271e-09 GAMG: Solving for p, Initial residual = 0.060360581, Final residual = 5.0485323e-05, No Iterations 2 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors (air): sum local = 1.4832474e-08, global = 3.1777063e-10, cumulative = -8.8575644e-10 DILUPBiCG: Solving for epsilon, Initial residual = 0.0034727837, Final residual = 3.0269126e-08, No Iterations 8 DILUPBiCG: Solving for k, Initial residual = 0.40375756, Final residual = 4.1736344e-08, 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.9972721e-05 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.7977311e-07, No Iterations 72 DICPCG: Solving for T, Initial residual = 0.079959901, Final residual = 9.5937062e-07, 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.9389894e-05 max: 0.30171661 deltaT = 4.4611786e-05 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 steady-state 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 |
|
April 12, 2011, 02:35 |
|
#6 |
Senior Member
Aram Amouzandeh
Join Date: Mar 2009
Location: Vienna, Vienna, Austria
Posts: 190
Rep Power: 17 |
hi!
you might have set Code:
adjustTimeStep yes; maxCo 0.3; Code:
adjustTimeStep no; cheers, aram |
|
April 12, 2011, 06:53 |
|
#7 |
Senior Member
Join Date: Mar 2011
Posts: 158
Rep Power: 15 |
Perfect!
That was exactly the point. :-) Thank you! |
|
April 18, 2012, 05:36 |
chtMultiRegionFoam
|
#8 |
New Member
namdar
Join Date: Mar 2012
Posts: 17
Rep Power: 14 |
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 |
|
April 18, 2012, 05:38 |
chtMultiRegionFoam
|
#9 |
New Member
namdar
Join Date: Mar 2012
Posts: 17
Rep Power: 14 |
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 |
|
April 18, 2012, 06:42 |
|
#10 | |
Senior Member
mauricio
Join Date: Jun 2011
Posts: 172
Rep Power: 18 |
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 pre-processing tools alone, but it is feasible. gl
__________________
Best Regards /calim "Elune will grant us the strength" |
||
April 24, 2012, 09:50 |
|
#11 |
New Member
namdar
Join Date: Mar 2012
Posts: 17
Rep Power: 14 |
Thanks a lot "Calim"
|
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Problem setting C_MU_T with DEFINE_ADJUST | Andy | FLUENT | 3 | June 22, 2018 02:53 |
HELP!: swimming pool flow problem | Tee | FLUENT | 1 | September 8, 2005 08:29 |
Problem setting AV | Marko | CFX | 0 | March 18, 2005 08:06 |
Problem setting C_MU_T with DEFINE_ADJUST | sarah_ron | FLUENT | 1 | February 24, 2005 12:11 |
Help setting up combustor problem | Tom | FLUENT | 0 | February 5, 2005 01:20 |