|
[Sponsors] |
September 11, 2013, 07:22 |
chtMultiRegionSimpleFoam
|
#1 |
Senior Member
Samuele Z
Join Date: Oct 2009
Location: Mozzate - Co - Italy
Posts: 520
Rep Power: 19 |
Dear All,
I am trying to run chtMultiRegionSimpleFoam. I have prepared a case with 1 fluid region and 6 solid regions. I have set my case and when I lauch it, I get this error: Code:
Adding to radiations Radiation model not active: radiationProperties not found Selecting radiationModel none Adding fvOptions No finite volume options present *** Reading solid mesh thermophysical properties for region domain5 Adding to thermos Selecting thermodynamics package { type heSolidThermo; mixture pureMixture; transport constIso; thermo hConst; equationOfState rhoConst; specie specie; energy sensibleEnthalpy; } Adding to radiations Radiation model not active: radiationProperties not found Selecting radiationModel none Adding fvOptions No finite volume options present Time = 1 Solving for fluid region part_2-solid DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 0.00888973, No Iterations 1 DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 0.00350072, No Iterations 1 DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 0.00662661, No Iterations 1 DILUPBiCG: Solving for h, Initial residual = 1, Final residual = 0.000600295, No Iterations 1 Min/max T:270 300 GAMG: Solving for p_rgh, Initial residual = 0.94945, Final residual = 0.00676794, No Iterations 7 time step continuity errors : sum local = 0.316082, global = -0.025282, cumulative = -0.025282 Min/max rho:1.15862 1.28736 Solving for solid region pcm --> FOAM FATAL ERROR: Attempt to cast type zeroGradient to type compressible::turbulentTemperatureCoupledBaffleMixed From function refCast<To>(From&) in file /home/zampini/OpenFOAM/OpenFOAM-2.2.0/src/OpenFOAM/lnInclude/typeInfo.H at line 114. FOAM aborting #0 Foam::error::printStack(Foam::Ostream&) in "/home/zampini/OpenFOAM/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #1 Foam::error::abort() in "/home/zampini/OpenFOAM/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #2 Foam::compressible::turbulentTemperatureCoupledBaffleMixedFvPatchScalarField const& Foam::refCast<Foam::compressible::turbulentTemperatureCoupledBaffleMixedFvPatchScalarField const, Foam::fvPatchField<double> const>(Foam::fvPatchField<double> const&) in "/home/zampini/OpenFOAM/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libcompressibleTurbulenceModel.so" #3 Foam::compressible::turbulentTemperatureCoupledBaffleMixedFvPatchScalarField::updateCoeffs() in "/home/zampini/OpenFOAM/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libcompressibleTurbulenceModel.so" #4 Foam::mixedFvPatchField<double>::evaluate(Foam::UPstream::commsTypes) in "/home/zampini/OpenFOAM/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so" #5 Foam::mixedEnergyFvPatchScalarField::updateCoeffs() in "/home/zampini/OpenFOAM/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so" #6 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricBoundaryField::updateCoeffs() in "/home/zampini/OpenFOAM/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam" #7 Foam::fvMatrix<double>::fvMatrix(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensionSet const&) in "/home/zampini/OpenFOAM/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam" #8 in "/home/zampini/OpenFOAM/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam" #9 in "/home/zampini/OpenFOAM/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam" #10 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #11 in "/home/zampini/OpenFOAM/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam" Annullato zampini@pc-zampini:~/Documenti/personali/Epta/SCC/steady$ I can not understand what it means. Could you help, please? Thanks a lot, Samuele Last edited by samiam1000; September 11, 2013 at 09:43. |
|
September 11, 2013, 17:22 |
|
#2 |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,981
Blog Entries: 45
Rep Power: 128 |
Hi Samuele,
Can you attach the boundary conditions for the "pcm" region, as well as any regions that are in touch with it? Although, if I have to guess, I think that this sentence: Code:
Attempt to cast type zeroGradient to type compressible::turbulentTemperatureCoupledBaffleMixed Best regards, Bruno
__________________
|
|
September 12, 2013, 03:27 |
|
#3 |
Senior Member
Samuele Z
Join Date: Oct 2009
Location: Mozzate - Co - Italy
Posts: 520
Rep Power: 19 |
Dear Bruno,
thanks for answering, first. Also, about your suggestions, I completely agree. The point is that I chacked the boundary conditions twice and I can't find any error. I am attaching to this email all the changeLogDictionary I use. Could you kindly have a look and let me know what's wrong with it? Thanks a lot, Samuele |
|
September 12, 2013, 04:53 |
|
#4 |
Senior Member
|
Dear Samuele,
it seems that you put this B.C for a patch in a region (compressible::turbulentTemperatureCoupledBaffleMi xed). then you put zeroGradient to the same patch in the other region which is not applicable. the B.C must be the same for same patch in different regions. hope it helps. BR, |
|
September 12, 2013, 05:08 |
|
#5 |
Senior Member
Samuele Z
Join Date: Oct 2009
Location: Mozzate - Co - Italy
Posts: 520
Rep Power: 19 |
I will check it again, then.
Thanks a lot, Samuele |
|
September 12, 2013, 06:07 |
|
#6 |
Senior Member
Samuele Z
Join Date: Oct 2009
Location: Mozzate - Co - Italy
Posts: 520
Rep Power: 19 |
I have found a first error: I haven't set the right BC in the fluid region.
I ran my case and I get a different error: could you help in solving this, too? Thanks a lot, Samuele The error message is: Code:
/*---------------------------------------------------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.2.1 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ Build : 2.2.1-57f3c3617a2d Exec : chtMultiRegionSimpleFoam Date : Sep 12 2013 Time : 11:47:31 Host : "lab-laptop" PID : 7539 Case : /home/lab/Documenti/Ethics/FRISBEE/CFD/SCC/steady nProcs : 1 sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster allowSystemOperations : Disallowing user-supplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create fluid mesh for region part_2-solid for time = 0 Create solid mesh for region pcm for time = 0 Create solid mesh for region packs_1 for time = 0 Create solid mesh for region packs_2 for time = 0 Create solid mesh for region part_2-solid.1 for time = 0 Create solid mesh for region domain2 for time = 0 Create solid mesh for region domain5 for time = 0 *** Reading fluid mesh thermophysical properties for region part_2-solid Adding to thermoFluid Selecting thermodynamics package { type heRhoThermo; mixture pureMixture; transport const; thermo hConst; equationOfState perfectGas; specie specie; energy sensibleEnthalpy; } Adding to rhoFluid Adding to UFluid Adding to phiFluid Adding to gFluid Adding to turbulence Selecting turbulence model type laminar Adding to ghFluid Adding to ghfFluid Radiation model not active: radiationProperties not found Selecting radiationModel none Adding fvOptions No finite volume options present *** Reading solid mesh thermophysical properties for region pcm Adding to thermos Selecting thermodynamics package { type heSolidThermo; mixture pureMixture; transport constIso; thermo hConst; equationOfState rhoConst; specie specie; energy sensibleEnthalpy; } Adding to radiations Radiation model not active: radiationProperties not found Selecting radiationModel none Adding fvOptions No finite volume options present *** Reading solid mesh thermophysical properties for region packs_1 Adding to thermos Selecting thermodynamics package { type heSolidThermo; mixture pureMixture; transport constIso; thermo hConst; equationOfState rhoConst; specie specie; energy sensibleEnthalpy; } Adding to radiations Radiation model not active: radiationProperties not found Selecting radiationModel none Adding fvOptions No finite volume options present *** Reading solid mesh thermophysical properties for region packs_2 Adding to thermos Selecting thermodynamics package { type heSolidThermo; mixture pureMixture; transport constIso; thermo hConst; equationOfState rhoConst; specie specie; energy sensibleEnthalpy; } Adding to radiations Radiation model not active: radiationProperties not found Selecting radiationModel none Adding fvOptions No finite volume options present *** Reading solid mesh thermophysical properties for region part_2-solid.1 Adding to thermos Selecting thermodynamics package { type heSolidThermo; mixture pureMixture; transport constIso; thermo hConst; equationOfState rhoConst; specie specie; energy sensibleEnthalpy; } Adding to radiations Radiation model not active: radiationProperties not found Selecting radiationModel none Adding fvOptions No finite volume options present *** Reading solid mesh thermophysical properties for region domain2 Adding to thermos Selecting thermodynamics package { type heSolidThermo; mixture pureMixture; transport constIso; thermo hConst; equationOfState rhoConst; specie specie; energy sensibleEnthalpy; } Adding to radiations Radiation model not active: radiationProperties not found Selecting radiationModel none Adding fvOptions No finite volume options present *** Reading solid mesh thermophysical properties for region domain5 Adding to thermos Selecting thermodynamics package { type heSolidThermo; mixture pureMixture; transport constIso; thermo hConst; equationOfState rhoConst; specie specie; energy sensibleEnthalpy; } Adding to radiations Radiation model not active: radiationProperties not found Selecting radiationModel none Adding fvOptions No finite volume options present Time = 1 Solving for fluid region part_2-solid DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 0.00888973, No Iterations 1 DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 0.00350072, No Iterations 1 DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 0.00662661, No Iterations 1 DILUPBiCG: Solving for h, Initial residual = 1, Final residual = 0.00187789, No Iterations 1 Min/max T:270 300 GAMG: Solving for p_rgh, Initial residual = 0.949213, Final residual = 0.00675391, No Iterations 7 time step continuity errors : sum local = 0.316054, global = -0.02528, cumulative = -0.02528 Min/max rho:1.15862 1.28736 Solving for solid region pcm DICPCG: Solving for h, Initial residual = 1, Final residual = 0.0743357, No Iterations 1 Min/max T:min(T) [0 0 0 1 0 0 0] 270 max(T) [0 0 0 1 0 0 0] 300 Solving for solid region packs_1 DICPCG: Solving for h, Initial residual = 1, Final residual = 0.0768139, No Iterations 1 Min/max T:min(T) [0 0 0 1 0 0 0] 288.897 max(T) [0 0 0 1 0 0 0] 300 Solving for solid region packs_2 #0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #1 Foam::sigSegv::sigHandler(int) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #2 in "/lib/x86_64-linux-gnu/libc.so.6" #3 Foam::polyMeshTetDecomposition::findFaceBasePts(Foam::polyMesh const&, double, bool) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #4 Foam::polyMesh::tetBasePtIs() const in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #5 Foam::mappedPatchBase::facePoints(Foam::polyPatch const&) const in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libmeshTools.so" #6 Foam::mappedPatchBase::calcMapping() const in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libmeshTools.so" #7 void Foam::mappedPatchBase::distribute<double>(Foam::List<double>&) const in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfiniteVolume.so" #8 Foam::compressible::turbulentTemperatureCoupledBaffleMixedFvPatchScalarField::updateCoeffs() in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libcompressibleTurbulenceModel.so" #9 Foam::mixedFvPatchField<double>::evaluate(Foam::UPstream::commsTypes) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfiniteVolume.so" #10 Foam::mixedEnergyFvPatchScalarField::updateCoeffs() in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so" #11 Foam::fvMatrix<double>::fvMatrix(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensionSet const&) in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam" #12 in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam" #13 in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam" #14 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #15 in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam" Segmentation fault (core dumped) lab@lab-laptop:~/Documenti/Ethics/FRISBEE/CFD/SCC/steady$ |
|
September 12, 2013, 06:37 |
|
#7 |
Senior Member
|
Dear Samuele,
this error always appears when there is a problem with three factors: cell size, time step, velocity value. you must compromise the three factors to get smooth run. hint: adjust dimensions and cell size such that patches boundaries lays on cells boundaries not in middle of it. hope it helps, Best Regards, Ahmed |
|
September 12, 2013, 06:49 |
|
#8 |
Senior Member
Samuele Z
Join Date: Oct 2009
Location: Mozzate - Co - Italy
Posts: 520
Rep Power: 19 |
Do you mean that it could be necessary to remesh my geometry?
|
|
September 12, 2013, 08:19 |
|
#10 |
Senior Member
Samuele Z
Join Date: Oct 2009
Location: Mozzate - Co - Italy
Posts: 520
Rep Power: 19 |
I think that being a steady simulation, the time-step does not influence the solution: is this right?
|
|
September 14, 2013, 10:27 |
|
#11 | ||
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,981
Blog Entries: 45
Rep Power: 128 |
Greetings to all!
@Samuele: Quote:
The unusual example might be LTSInterFoam: http://www.openfoam.org/version2.0.0/steady-vof.php Quote:
Code:
checkMesh -allGeometry -allTopology Bruno
__________________
|
|||
September 15, 2013, 12:39 |
|
#12 |
Senior Member
Samuele Z
Join Date: Oct 2009
Location: Mozzate - Co - Italy
Posts: 520
Rep Power: 19 |
Dear Bruno,
thanks for answering and pardon for the late reply. This is the output of the command you suggested: Code:
zampini@pc-zampini:~/Documenti/personali/Epta/SCC/steady$ checkMesh -allGeometry -allTopology /*---------------------------------------------------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.2.0 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ Build : 2.2.0 Exec : checkMesh -allGeometry -allTopology Date : Sep 15 2013 Time : 18:37:51 Host : "pc-zampini" PID : 1944 Case : /home/zampini/Documenti/personali/Epta/SCC/steady nProcs : 1 sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster allowSystemOperations : Disallowing user-supplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create polyMesh for time = 0 Enabling all (cell, face, edge, point) topology checks. Enabling all geometry checks. Time = 0 Mesh stats points: 621150 faces: 1816951 internal faces: 1761449 cells: 596440 faces per cell: 5.9996 boundary patches: 13 point zones: 0 face zones: 14 cell zones: 5 Overall number of cells of each type: hexahedra: 596200 prisms: 240 wedges: 0 pyramids: 0 tet wedges: 0 tetrahedra: 0 polyhedra: 0 Checking topology... Boundary definition OK. Cell to face addressing OK. Point usage OK. Upper triangular ordering OK. Face vertices OK. Topological cell zip-up check OK. Number of identical duplicate faces (baffle faces): 3200 Face-face connectivity OK. <<Writing 6400 faces with non-standard edge connectivity to set edgeFaces <<Writing 4 cells with two non-boundary faces to set twoInternalFacesCells Number of regions: 1 (OK). Checking patch topology for multiply connected surfaces... Patch Faces Points Surface topology Bounding box deflettori 6400 3403 multiply connected (shared edge) (-0.1 -0.05 0.45) (1.3 0.3 0.5) foam-pcm 4200 4402 ok (non-closed singly connected) (-0.05 0 -0.05) (1.25 0.3 -0.05) foam-part_2-solid 14391 14810 ok (non-closed singly connected) (-0.1 -0.05 -0.05) (1.3 0.3 0.75) foam-part_2-solid.1 8400 8662 ok (non-closed singly connected) (-0.05 -0.05 -0.05) (1.25 0 0.45) glass 5600 5781 ok (non-closed singly connected) (-0.1 -0.05 0.75) (1.3 0.3 0.75) inlet_1 400 451 ok (non-closed singly connected) (0.55 -0.05 -0.05) (0.6 0.3 -0.05) inlet_2 400 451 ok (non-closed singly connected) (0.6 -0.05 -0.05) (0.65 0.3 -0.05) intake_1 400 451 ok (non-closed singly connected) (-0.1 -0.05 -0.05) (-0.05 0.3 -0.05) intake_2 400 451 ok (non-closed singly connected) (1.25 -0.05 -0.05) (1.3 0.3 -0.05) symmetry-packs_1 2250 2346 ok (non-closed singly connected) (0 0.3 0) (0.5 0.3 0.45) symmetry-pcm 2500 2832 ok (non-closed singly connected) (-0.05 0.3 -0.05) (1.25 0.3 0.45) symmetry-part_2-solid7911 8250 ok (non-closed singly connected) (-0.1 0.3 -0.05) (1.3 0.3 0.75) symmetry-packs_2 2250 2346 ok (non-closed singly connected) (0.7 0.3 0) (1.2 0.3 0.45) <<Writing 3391 conflicting points to set nonManifoldPoints Checking geometry... Overall domain bounding box (-0.1 -0.05 -0.05) (1.3 0.3 0.75) Mesh (non-empty, non-wedge) directions (1 1 1) Mesh (non-empty) directions (1 1 1) Boundary openness (3.24011e-16 3.56288e-17 1.47307e-15) OK. Max cell openness = 3.10576e-16 OK. Max aspect ratio = 13.8392 OK. Minimum face area = 7.67269e-06. Maximum face area = 0.000250165. Face area magnitudes OK. Min volume = 3.06526e-08. Max volume = 2.41892e-06. Total volume = 0.392. Cell volumes OK. Mesh non-orthogonality Max: 54.2303 average: 11.4654 Non-orthogonality check OK. Face pyramids OK. Max skewness = 1.05523 OK. Coupled point location match (average 0) OK. Face tets OK. Min/max edge length = 0.00290397 0.0199751 OK. All angles in faces OK. Face flatness (1 = flat, 0 = butterfly) : average = 1 min = 1 All face flatness OK. Cell determinant (wellposedness) : minimum: 0 average: 5.679 ***Cells with small determinant found, number of cells: 80 <<Writing 80 under-determined cells to set underdeterminedCells Concave cell check OK. Failed 1 mesh checks. End |
|
September 15, 2013, 14:58 |
|
#13 | |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,981
Blog Entries: 45
Rep Power: 128 |
Hi Samuele,
Quote:
paraFoam provides you with the ability to also see the sets. Turn on that option and choose to see "underdeterminedCells" that should appear in the same list as the patches. Then try to see where exactly where the problem cells are and try to re-do your mesh. Another possibility is to follow the example shown here: http://openfoamwiki.net/index.php/SetSet#Usage_example - more specifically, to only remove the cells associated to "underdeterminedCells". But keep in mind that this kind of cell removal strategy has certain limitations, such as possibly and wrongly removing some important cells. Best regards, Bruno
__________________
|
||
September 16, 2013, 03:04 |
|
#14 |
Senior Member
Samuele Z
Join Date: Oct 2009
Location: Mozzate - Co - Italy
Posts: 520
Rep Power: 19 |
Dear Bruno,
I am attaching a picture of the whole volume where it is evident where the underdeterminedCells are. I can't understand what's wrong with them. Do you have any idea? First of all, I will try your suggestions. Thanks a lot, Samuele |
|
September 18, 2013, 06:41 |
|
#15 |
Senior Member
Samuele Z
Join Date: Oct 2009
Location: Mozzate - Co - Italy
Posts: 520
Rep Power: 19 |
Dear Bruno, Dear All,
after having re-meshed my geometry, I get a very strange result. First of all, all the mesh checks are ok! Hance I thought that my simulation would have started immediately, but.. ..but I got this error: Code:
lab@lab-laptop:~/Documenti/Ethics/FRISBEE/CFD/SCC/steady$ chtMultiRegionSimpleFoam /*---------------------------------------------------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.2.1 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ Build : 2.2.1-57f3c3617a2d Exec : chtMultiRegionSimpleFoam Date : Sep 18 2013 Time : 12:36:14 Host : "lab-laptop" PID : 4908 Case : /home/lab/Documenti/Ethics/FRISBEE/CFD/SCC/steady nProcs : 1 sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster allowSystemOperations : Disallowing user-supplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create fluid mesh for region part_2-solid for time = 0 Create solid mesh for region pcm for time = 0 Create solid mesh for region packs_1 for time = 0 Create solid mesh for region packs_2 for time = 0 Create solid mesh for region part_2-solid.1 for time = 0 Create solid mesh for region domain2 for time = 0 Create solid mesh for region domain5 for time = 0 *** Reading fluid mesh thermophysical properties for region part_2-solid Adding to thermoFluid Selecting thermodynamics package { type heRhoThermo; mixture pureMixture; transport const; thermo hConst; equationOfState perfectGas; specie specie; energy sensibleEnthalpy; } #0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #1 Foam::sigSegv::sigHandler(int) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #2 in "/lib/x86_64-linux-gnu/libc.so.6" #3 Foam::tmp<Foam::Field<double> > Foam::fvPatch::patchInternalField<double>(Foam::UList<double> const&) const in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam" #4 Foam::fvPatchField<double>::patchInternalField() const in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam" #5 Foam::basicSymmetryFvPatchField<double>::evaluate(Foam::UPstream::commsTypes) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfiniteVolume.so" #6 Foam::symmetryFvPatchField<double>::symmetryFvPatchField(Foam::fvPatch const&, Foam::DimensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfiniteVolume.so" #7 Foam::fvPatchField<double>::adddictionaryConstructorToTable<Foam::symmetryFvPatchField<double> >::New(Foam::fvPatch const&, Foam::DimensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfiniteVolume.so" #8 Foam::fvPatchField<double>::New(Foam::fvPatch const&, Foam::DimensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam" #9 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricBoundaryField::readField(Foam::DimensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam" #10 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::readFields(Foam::dictionary const&) in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam" #11 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::readFields() in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam" #12 at basicThermo.C:0 #13 Foam::basicThermo::lookupOrConstruct(Foam::fvMesh const&, char const*) const in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so" #14 Foam::basicThermo::basicThermo(Foam::fvMesh const&, Foam::word const&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so" #15 Foam::fluidThermo::fluidThermo(Foam::fvMesh const&, Foam::word const&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so" #16 Foam::rhoThermo::rhoThermo(Foam::fvMesh const&, Foam::word const&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so" #17 Foam::heThermo<Foam::rhoThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > >::heThermo(Foam::fvMesh const&, Foam::word const&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so" #18 Foam::rhoThermo::addfvMeshConstructorToTable<Foam::heRhoThermo<Foam::rhoThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > > >::New(Foam::fvMesh const&, Foam::word const&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so" #19 Foam::autoPtr<Foam::rhoThermo> Foam::basicThermo::New<Foam::rhoThermo>(Foam::fvMesh const&, Foam::word const&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so" #20 Foam::rhoThermo::New(Foam::fvMesh const&, Foam::word const&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so" #21 in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam" #22 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #23 in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam" Segmentation fault (core dumped) lab@lab-laptop:~/Documenti/Ethics/FRISBEE/CFD/SCC/steady$ And this happens for each region. Do you have any idea? Could you help? |
|
September 21, 2013, 15:22 |
|
#16 | |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,981
Blog Entries: 45
Rep Power: 128 |
Hi Samuele,
There seems to be a problem with a patch that is defined to be a symmetry plane: Quote:
Best regards, Bruno
__________________
|
||
September 22, 2013, 16:43 |
|
#17 |
Senior Member
Samuele Z
Join Date: Oct 2009
Location: Mozzate - Co - Italy
Posts: 520
Rep Power: 19 |
Here is the case: https://www.dropbox.com/sh/tgdwuqkfodgdffk/zFNdkvUgjH
Could you have a look? Thanks a lot, Samuele |
|
September 23, 2013, 06:48 |
|
#18 |
Super Moderator
Tobias Holzmann
Join Date: Oct 2010
Location: Bad Wörishofen
Posts: 2,711
Blog Entries: 6
Rep Power: 51 |
Code:
Solving for solid region packs_2 #0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #1 Foam::sigSegv::sigHandler(int) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #2 in "/lib/x86_64-linux-gnu/libc.so.6" #3 Foam::polyMeshTetDecomposition::findFaceBasePts(Foam::polyMesh const&, double, bool) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #4 Foam::polyMesh::tetBasePtIs() const in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #5 Foam::mappedPatchBase::facePoints(Foam::polyPatch const&) const in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libmeshTools.so" #6 Foam::mappedPatchBase::calcMapping() const in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libmeshTools.so" #7 void Foam::mappedPatchBase::distribute<double>(Foam::List<double>&) const in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfiniteVolume.so" #8 Foam::compressible::turbulentTemperatureCoupledBaffleMixedFvPatchScalarField::updateCoeffs() in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libcompressibleTurbulenceModel.so" #9 Foam::mixedFvPatchField<double>::evaluate(Foam::UPstream::commsTypes) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfiniteVolume.so" #10 Foam::mixedEnergyFvPatchScalarField::updateCoeffs() in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so" #11 Foam::fvMatrix<double>::fvMatrix(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensionSet const&) in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam" #12 in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam" #13 in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam" #14 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #15 in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam" Segmentation fault (core dumped) Code:
this error always appears when there is a problem with three factors: cell size, time step, velocity value. you must compromise the three factors to get smooth run. If I am right you have a Problem in the Boundary Conditions of T in packs_2 or in your fluid Region. Maybe your boundary file is wrong (patch type). Maybe you have no value set.? Your mesh seems okay. Regards Tobi |
|
September 23, 2013, 08:02 |
|
#19 |
Senior Member
Samuele Z
Join Date: Oct 2009
Location: Mozzate - Co - Italy
Posts: 520
Rep Power: 19 |
Dear Tobi,
thanks for answering. Actually, I have solved this very problem (it was due to a bad definition of the boudary conditions) and the simulation's running. However, the temperature seems to be meaningless: I do have a max temperature of about 740000 K. Too much, I say. I am going to check this problem, too. Any idea to begin to investigate the issue? Thanks a lot, Samuele. |
|
September 24, 2013, 06:20 |
|
#20 |
Super Moderator
Tobias Holzmann
Join Date: Oct 2010
Location: Bad Wörishofen
Posts: 2,711
Blog Entries: 6
Rep Power: 51 |
Hi,
write out the first 10 or 20 integrations and have a look at your Domain. You will be able to see the regions where you get the high temperature values. It could be possible that this Problem occure due to a mesh Problem. otherwise you see if your BC are incorrect or your Settings are wrong. Good luck |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|