|
[Sponsors] |
May 25, 2023, 10:34 |
Error in chtMultiRegionFoam
|
#1 |
New Member
TGS
Join Date: May 2023
Location: United Kingdom
Posts: 29
Rep Power: 3 |
I have three regions in my OpenFoam simulation. My geometry contains a solid fin that was covered by air at room temperature. The air flow is required to go through the fin as well. The one side wall of a fin has a temperature of 400 K. This is simply an air-cooled fin scenario. To solve this problem, I'm using the chtMultiRegionFoam solver, and I've successfully done blockMesh, TopoSet, splitMultiRegions, and snappyHexMesh. However, while I am running the chtMultiRegionsFoam, there is an error generated, and the error looks like this: I am new to OpenFoam and please help me to solve this error.
Thanks. TGS /*---------------------------------------------------------------------------*\ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | Website: https://openfoam.org \\ / A nd | Version: 10 \\/ M anipulation | \*---------------------------------------------------------------------------*/ Build : 10-c4cf895ad8fa Exec : chtMultiRegionFoam Date : May 25 2023 Time : 12:20:26 Host : "AIE-OpenFOAM" PID : 129243 I/O : uncollated Case : /home/thamasha/OpenFOAM/thamasha-10/run/NewFinT nProcs : 1 sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10) allowSystemOperations : Allowing user-supplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create fluid mesh for region air for time = 0 Create solid mesh for region fin for time = 0 Create solid mesh for region heatedBoundary for time = 0 PIMPLE: Region air PIMPLE: No convergence criteria found PIMPLE: Region fin PIMPLE: No convergence criteria found PIMPLE: Region heatedBoundary PIMPLE: No convergence criteria found PIMPLE: Operating solver in steady-state mode with 1 outer corrector PIMPLE: Operating solver in SIMPLE mode *** Reading fluid mesh thermophysical properties for region air Adding to thermoFluid Selecting thermodynamics package { type heRhoThermo; mixture pureMixture; transport const; thermo hConst; equationOfState rhoConst; specie specie; energy sensibleEnthalpy; } Adding to rhoFluid Adding to UFluid Adding to phiFluid Adding to gFluid Adding to hRefFluid Adding to pRefFluid Adding to ghFluid Adding to ghfFluid Adding to turbulenceFluid Selecting turbulence model type laminar Selecting laminar stress model Stokes Adding to thermophysicalTransport Selecting thermophysical transport type laminar Selecting default laminar thermophysical transport model unityLewisFourier Adding to reactionFluid Combustion model not active: combustionProperties not found Selecting combustion model none Adding to KFluid Adding to dpdtFluid Adding to fieldsFluid Adding MRF No MRF models present Adding fvModelsFluid Creating fvModels from "constant/fvModels" Selecting finite volume model type interRegionHeatTransfer Name: air_to_heatedBoundary Selecting heatTransferModel constant Adding fvConstraintsFluid --> FOAM Warning : Creating fvConstraints from "system/fvOptions" Selecting finite volume constraint type limitTemperature Name: limitT - selecting all cells - selected 1057 cell(s) with volume 79275 *** Reading solid mesh thermophysical properties for region fin Adding to thermoSolid Selecting thermodynamics package { type heSolidThermo; mixture pureMixture; transport constIsoSolid; thermo eConst; equationOfState rhoConst; specie specie; energy sensibleInternalEnergy; } Adding to fvModelsSolid --> FOAM Warning : Creating fvModels from "system/fvOptions" rather than constant/fvModels Adding fvConstraintsSolid --> FOAM Warning : Creating fvConstraints from "system/fvOptions" Selecting finite volume constraint type limitTemperature Name: limitT - selecting all cells - selected 62 cell(s) with volume 4650 *** Reading solid mesh thermophysical properties for region heatedBoundary Adding to thermoSolid Selecting thermodynamics package { type heSolidThermo; mixture pureMixture; transport constIsoSolid; thermo eConst; equationOfState rhoConst; specie specie; energy sensibleInternalEnergy; } Adding to fvModelsSolid Creating fvModels from "constant/fvModels" Selecting finite volume model type interRegionHeatTransfer Name: heatedBoundary_to_air Selecting heatTransferModel constant Adding fvConstraintsSolid Creating fvConstraints from "system/fvConstraints" Selecting finite volume constraint type fixedTemperatureConstraint Name: fixedTemperature - selecting all cells - selected 1 cell(s) with volume 75 Region: air Courant Number mean: 1.35289e-05 max: 5e-05 Region: fin Diffusion Number mean: 4.88372e-10 max: 3.65798e-09 --> FOAM Warning : From function Type Foam::gAverage(const Foam::UList<T>&, Foam::label) [with Type = double; Foam::label = int] in file /home/ubuntu/OpenFOAM/OpenFOAM-10/src/OpenFOAM/lnInclude/FieldFunctions.C at line 597 empty field, returning zero. Region: heatedBoundary Diffusion Number mean: 0 max: 3.65798e-09 --> FOAM Warning : From function virtual void Foam:robes::findElements(const Foam::fvMesh&) in file probes/probes.C at line 117 Did not find location (-0.054 0.018 0.017) in any cell. Skipping location. --> FOAM Warning : From function virtual void Foam:robes::findElements(const Foam::fvMesh&) in file probes/probes.C at line 117 Did not find location (-0.054 0.018 -0.017) in any cell. Skipping location. --> FOAM Warning : From function virtual void Foam:robes::findElements(const Foam::fvMesh&) in file probes/probes.C at line 117 Did not find location (-0.054 0.018 0) in any cell. Skipping location. --> FOAM Warning : From function virtual void Foam:robes::findElements(const Foam::fvMesh&) in file probes/probes.C at line 117 Did not find location (-0.07 0.018 0) in any cell. Skipping location. fieldAverage volAvgFinT: Starting averaging at time 0 fieldAverage volAvgAirT: Starting averaging at time 0 deltaT = 0.0001 Region: air Courant Number mean: 1.35289e-05 max: 5e-05 Region: fin Diffusion Number mean: 4.88372e-10 max: 3.65798e-09 --> FOAM Warning : From function Type Foam::gAverage(const Foam::UList<T>&, Foam::label) [with Type = double; Foam::label = int] in file /home/ubuntu/OpenFOAM/OpenFOAM-10/src/OpenFOAM/lnInclude/FieldFunctions.C at line 597 empty field, returning zero. Region: heatedBoundary Diffusion Number mean: 0 max: 3.65798e-09 deltaT = 0.000119048 Time = 0.000119048s Solving for solid region fin GAMG: Solving for e, Initial residual = 0.465165, Final residual = 0.401726, No Iterations 1000 Min/max T:300 300 Solving for solid region heatedBoundary - selecting inter region mapping Creating mesh-to-mesh addressing for heatedBoundary and air regions using cellVolumeWeight Overlap volume: 0 diagonal: Solving for e, Initial residual = 0, Final residual = 0, No Iterations 0 Min/max T:300 399.987 Solving for fluid region air DILUPBiCGStab: Solving for Ux, Initial residual = 0.998985, Final residual = 0.00181951, No Iterations 1 DILUPBiCGStab: Solving for Uy, Initial residual = 0.997048, Final residual = 0.00222112, No Iterations 1 DILUPBiCGStab: Solving for Uz, Initial residual = 1, Final residual = 0.00429678, No Iterations 1 --> FOAM FATAL ERROR: request for polyMesh region0 from objectRegistry NewFinT failed available objects of type polyMesh are 3 ( air fin heatedBoundary ) From function const Type& Foam:bjectRegistry::lookupObject(const Foam::word&) const [with Type = Foam:olyMesh] in file /home/ubuntu/OpenFOAM/OpenFOAM-10/src/OpenFOAM/lnInclude/objectRegistryTemplates.C at line 211. FOAM aborting #0 Foam::error:rintStack(Foam::Ostream&) at ??:? #1 Foam::error::abort() at ??:? #2 Foam:olyMesh const& Foam:bjectRegistry::lookupObject<Foam:olyMesh> (Foam::word const&) const at ??:? #3 Foam::compressible::turbulentTemperatureCoupledBaf fleMixedFvPatchScalarField::updateCoeffs() at ??:? #4 Foam::mixedFvPatchField<double>::evaluate(Foam::UP stream::commsTypes) at ??:? #5 Foam::mixedEnergyFvPatchScalarField::updateCoeffs( ) at ??:? #6 Foam::fvMatrix<double>::fvMatrix(Foam::GeometricFi eld<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensionSet const&) in "/opt/openfoam10/platforms/linux64GccDPInt32Opt/bin/chtMultiRegionFoam" #7 Foam::tmp<Foam::fvMatrix<double> > Foam::fvModels::source<double, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> >(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&, Foam::dimensionSet const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/opt/openfoam10/platforms/linux64GccDPInt32Opt/bin/chtMultiRegionFoam" #8 ? in "/opt/openfoam10/platforms/linux64GccDPInt32Opt/bin/chtMultiRegionFoam" #9 ? in "/lib/x86_64-linux-gnu/libc.so.6" #10 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #11 ? in "/opt/openfoam10/platforms/linux64GccDPInt32Opt/bin/chtMultiRegionFoam" Aborted (core dumped) |
|
May 25, 2023, 15:53 |
|
#2 |
Senior Member
Lucky
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 5,762
Rep Power: 66 |
You have 3 meshed regions named air, fin, and heated boundary. Your build is trying to find a meshed region that has the name "region0" and it doesn't exist.
Search all your dicts and find region0 and rename it to the correct region. Since the solver stops right after completing one step through the air region but before printing the Courant number, the most likely culprit is probably in your controlDict in the CourantNo object. You most likely yoinked it from a tutorial and forgot to specify the region. Change region0 to air. Code:
Co1 { // Mandatory entries (unmodifiable) type CourantNo; libs (fieldFunctionObjects); // Optional entries (runtime modifiable) rho rho; // Optional (inherited) entries field <phi>; result <fieldResult>; region region0; enabled true; log true; timeStart 0; timeEnd 1000; executeControl timeStep; executeInterval 1; writeControl timeStep; writeInterval 1; } |
|
May 26, 2023, 04:50 |
|
#3 |
New Member
TGS
Join Date: May 2023
Location: United Kingdom
Posts: 29
Rep Power: 3 |
Dear Lucky,
Thank you very much for the detailed answer and your time of helping me. I will try according to your advice and let you know. Many thanks again. TGS |
|
May 26, 2023, 05:56 |
|
#4 |
New Member
TGS
Join Date: May 2023
Location: United Kingdom
Posts: 29
Rep Power: 3 |
Dear Lucky,
I tried by using the function that you mentioned in the thread to Courant number and changed the region to air. However, still I am getting the same error . Here I attached my controlDict file and the error message as well. If you have any free time ,can you please check it. many thanks in advance. TGS. ControlDict file. /*--------------------------------*- C++ -*----------------------------------*\ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | Website: https://openfoam.org \\ / A nd | Version: 8 \\/ M anipulation | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; location "system"; object controlDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // application chtMultiRegionFoam; startFrom startTime; startTime 0; stopAt endTime; endTime 10; deltaT 0.0001; writeControl adjustableRunTime; writeInterval 0.01; purgeWrite 0; writeFormat binary; writePrecision 6; writeCompression off; timeFormat general; timePrecision 6; runTimeModifiable on; adjustTimeStep on; functions { probes { libs ("libsampling.so"); type probes; // Name of the directory for probe data. name probes; // Write at same frequency as fields. region air ; enabled true; log true; timeStart 0; timeEnd 10; executeControl timeStep; executeInterval 1; writeControl timeStep; writeInterval 1; // Fields to be probed. fields ( p U T ); probeLocations ( ( 26.25 18.75 1.25 ) // Exit from fin. ( 27.00 10.00 35.00 ) // Entrance to fin ( 28.50 20.00 40.00 ) // Above fin. ); } Co1 { // Mandatory entries (unmodifiable) type CourantNo; libs ("libfieldFunctionObjects.so"); // Optional entries (runtime modifiable) // Optional (inherited) entries region air; enabled true; log true; timeStart 0; timeEnd 10; executeControl timeStep; executeInterval 1; writeControl timeStep; writeInterval 1; } } // ************************************************** *********************** // Error Message *---------------------------------------------------------------------------*\ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | Website: https://openfoam.org \\ / A nd | Version: 10 \\/ M anipulation | \*---------------------------------------------------------------------------*/ Build : 10-c4cf895ad8fa Exec : chtMultiRegionFoam Date : May 26 2023 Time : 09:48:52 Host : "AIE-OpenFOAM" PID : 144895 I/O : uncollated Case : /home/thamasha/OpenFOAM/thamasha-10/run/NewFinT nProcs : 1 sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10) allowSystemOperations : Allowing user-supplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create fluid mesh for region air for time = 0 Create solid mesh for region fin for time = 0 Create solid mesh for region heatedBoundary for time = 0 PIMPLE: Region air PIMPLE: No convergence criteria found PIMPLE: Region fin PIMPLE: No convergence criteria found PIMPLE: Region heatedBoundary PIMPLE: No convergence criteria found PIMPLE: Operating solver in steady-state mode with 1 outer corrector PIMPLE: Operating solver in SIMPLE mode *** Reading fluid mesh thermophysical properties for region air Adding to thermoFluid Selecting thermodynamics package { type heRhoThermo; mixture pureMixture; transport const; thermo hConst; equationOfState rhoConst; specie specie; energy sensibleEnthalpy; } Adding to rhoFluid Adding to UFluid Adding to phiFluid Adding to gFluid Adding to hRefFluid Adding to pRefFluid Adding to ghFluid Adding to ghfFluid Adding to turbulenceFluid Selecting turbulence model type laminar Selecting laminar stress model Stokes Adding to thermophysicalTransport Selecting thermophysical transport type laminar Selecting default laminar thermophysical transport model unityLewisFourier Adding to reactionFluid Combustion model not active: combustionProperties not found Selecting combustion model none Adding to KFluid Adding to dpdtFluid Adding to fieldsFluid Adding MRF No MRF models present Adding fvModelsFluid Creating fvModels from "constant/fvModels" Selecting finite volume model type interRegionHeatTransfer Name: air_to_heatedBoundary Selecting heatTransferModel constant Adding fvConstraintsFluid --> FOAM Warning : Creating fvConstraints from "system/fvOptions" Selecting finite volume constraint type limitTemperature Name: limitT - selecting all cells - selected 1057 cell(s) with volume 79275 *** Reading solid mesh thermophysical properties for region fin Adding to thermoSolid Selecting thermodynamics package { type heSolidThermo; mixture pureMixture; transport constIsoSolid; thermo eConst; equationOfState rhoConst; specie specie; energy sensibleInternalEnergy; } Adding to fvModelsSolid --> FOAM Warning : Creating fvModels from "system/fvOptions" rather than constant/fvModels Adding fvConstraintsSolid --> FOAM Warning : Creating fvConstraints from "system/fvOptions" Selecting finite volume constraint type limitTemperature Name: limitT - selecting all cells - selected 62 cell(s) with volume 4650 *** Reading solid mesh thermophysical properties for region heatedBoundary Adding to thermoSolid Selecting thermodynamics package { type heSolidThermo; mixture pureMixture; transport constIsoSolid; thermo eConst; equationOfState rhoConst; specie specie; energy sensibleInternalEnergy; } Adding to fvModelsSolid Creating fvModels from "constant/fvModels" Selecting finite volume model type interRegionHeatTransfer Name: heatedBoundary_to_air Selecting heatTransferModel constant Adding fvConstraintsSolid Creating fvConstraints from "system/fvConstraints" Selecting finite volume constraint type fixedTemperatureConstraint Name: fixedTemperature - selecting all cells - selected 1 cell(s) with volume 75 Region: air Courant Number mean: 1.35289e-05 max: 5e-05 Region: fin Diffusion Number mean: 4.88372e-10 max: 3.65798e-09 --> FOAM Warning : From function Type Foam::gAverage(const Foam::UList<T>&, Foam::label) [with Type = double; Foam::label = int] in file /home/ubuntu/OpenFOAM/OpenFOAM-10/src/OpenFOAM/lnInclude/FieldFunctions.C at line 597 empty field, returning zero. Region: heatedBoundary Diffusion Number mean: 0 max: 3.65798e-09 deltaT = 0.0001 functionObjects::CourantNo Co1 writing field: Co Region: air Courant Number mean: 1.35289e-05 max: 5e-05 Region: fin Diffusion Number mean: 4.88372e-10 max: 3.65798e-09 --> FOAM Warning : From function Type Foam::gAverage(const Foam::UList<T>&, Foam::label) [with Type = double; Foam::label = int] in file /home/ubuntu/OpenFOAM/OpenFOAM-10/src/OpenFOAM/lnInclude/FieldFunctions.C at line 597 empty field, returning zero. Region: heatedBoundary Diffusion Number mean: 0 max: 3.65798e-09 deltaT = 0.000119048 Time = 0.000119048s Solving for solid region fin GAMG: Solving for e, Initial residual = 0.465165, Final residual = 0.401726, No Iterations 1000 Min/max T:300 300 Solving for solid region heatedBoundary - selecting inter region mapping Creating mesh-to-mesh addressing for heatedBoundary and air regions using cellVolumeWeight Overlap volume: 0 diagonal: Solving for e, Initial residual = 0, Final residual = 0, No Iterations 0 Min/max T:300 399.987 Solving for fluid region air DILUPBiCGStab: Solving for Ux, Initial residual = 0.998985, Final residual = 0.00181951, No Iterations 1 DILUPBiCGStab: Solving for Uy, Initial residual = 0.997048, Final residual = 0.00222112, No Iterations 1 DILUPBiCGStab: Solving for Uz, Initial residual = 1, Final residual = 0.00429678, No Iterations 1 --> FOAM FATAL ERROR: request for polyMesh region0 from objectRegistry NewFinT failed available objects of type polyMesh are 3 ( air fin heatedBoundary ) From function const Type& Foam:bjectRegistry::lookupObject(const Foam::word&) const [with Type = Foam:olyMesh] in file /home/ubuntu/OpenFOAM/OpenFOAM-10/src/OpenFOAM/lnInclude/objectRegistryTemplates.C at line 211. FOAM aborting #0 Foam::error:rintStack(Foam::Ostream&) at ??:? #1 Foam::error::abort() at ??:? #2 Foam:olyMesh const& Foam:bjectRegistry::lookupObject<Foam:olyMesh> (Foam::word const&) const at ??:? #3 Foam::compressible::turbulentTemperatureCoupledBaf fleMixedFvPatchScalarField::updateCoeffs() at ??:? #4 Foam::mixedFvPatchField<double>::evaluate(Foam::UP stream::commsTypes) at ??:? #5 Foam::mixedEnergyFvPatchScalarField::updateCoeffs( ) at ??:? #6 Foam::fvMatrix<double>::fvMatrix(Foam::GeometricFi eld<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensionSet const&) in "/opt/openfoam10/platforms/linux64GccDPInt32Opt/bin/chtMultiRegionFoam" #7 Foam::tmp<Foam::fvMatrix<double> > Foam::fvModels::source<double, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> >(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&, Foam::dimensionSet const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/opt/openfoam10/platforms/linux64GccDPInt32Opt/bin/chtMultiRegionFoam" #8 ? in "/opt/openfoam10/platforms/linux64GccDPInt32Opt/bin/chtMultiRegionFoam" #9 ? in "/lib/x86_64-linux-gnu/libc.so.6" #10 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #11 ? in "/opt/openfoam10/platforms/linux64GccDPInt32Opt/bin/chtMultiRegionFoam" Aborted (core dumped) |
|
May 27, 2023, 02:33 |
|
#5 |
Senior Member
Lucky
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 5,762
Rep Power: 66 |
It still means that you have region region0 declared somewhere in one of your dicts. Ctrl + F is your friend
|
|
May 31, 2023, 04:38 |
|
#6 |
New Member
TGS
Join Date: May 2023
Location: United Kingdom
Posts: 29
Rep Power: 3 |
Dear Lucky,
Okay. I will search thank you so much. TGS |
|
Tags |
chtmultiregion, openfaom, region0 |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
chtMultiRegionFoam solver stops without any error | amol_patel | OpenFOAM Running, Solving & CFD | 4 | July 5, 2024 02:41 |
Help with PIMPLE algorithm in chtMultiregionFoam | Chris T | OpenFOAM Running, Solving & CFD | 0 | August 30, 2022 09:49 |
Error in thermophysical properties (chtMultiRegionFoam) | mukut | OpenFOAM Pre-Processing | 28 | November 23, 2021 07:34 |
Error in chtMultiRegionFoam | kirankarki | OpenFOAM | 6 | August 21, 2018 09:00 |
Simulation of a sample in a furnace w/ chtMultiRegionFoam | sergimart7 | OpenFOAM Running, Solving & CFD | 7 | August 12, 2015 07:48 |