Hallo Foamers,
Ich habe folgendes Problem: bei der Simulation von zwei Region (KK für Kühlkörper und bottomAir) nach Konvektion bekomme ich ständig die folgende Fehlermeldung:
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / 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 : chtMultiRegionSimpleFoam
Date : Aug 09 2013
Time : 13:30:00
Host : "gpi12175.iavgroup.local"
PID : 13403
Case : /usr2/bottomAir_KK/case
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 bottomAir for time = 0
Create solid mesh for region KK for time = 0
*** Reading fluid mesh thermophysical properties for region bottomAir
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 KK
Adding to thermos
Selecting thermodynamics package
{
type heSolidThermo;
mixture pureMixture;
transport constIso;
thermo hConst;
equationOfState rhoConst;
specie specie;
energy sensibleEnthalpy;
}
Adding to radiations
Selecting radiationModel opaqueSolid
Selecting absorptionEmissionModel constantAbsorptionEmission
Selecting scatterModel none
Adding fvOptions
No finite volume options present
Time = 1
Solving for fluid region bottomAir
DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 0.024221, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 0.0239137, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 0.0235421, No Iterations 2
DILUPBiCG: Solving for h, Initial residual = 1, Final residual = 0.0182296, No Iterations 3
Min/max T:294.282 300.428
GAMG: Solving for p_rgh, Initial residual = 0.841536, Final residual = 0.00288963, No Iterations 6
time step continuity errors : sum local = 0.0554715, global = 2.33029e-18, cumulative = 2.33029e-18
Min/max rho:1.15697 1.18114
Solving for solid region KK
DICPCG: Solving for h, Initial residual = 1, Final residual = 0.0205543, No Iterations 2
Min/max T:min(T) [0 0 0 1 0 0 0] 296.988 max(T) [0 0 0 1 0 0 0] 360
ExecutionTime = 1.8 s ClockTime = 2 s
Time = 2
Solving for fluid region bottomAir
DILUPBiCG: Solving for Ux, Initial residual = 0.492912, Final residual = 0.0174365, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.197266, Final residual = 0.01101, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.488943, Final residual = 0.0202762, No Iterations 2
DILUPBiCG: Solving for h, Initial residual = 0.426496, Final residual = 0.0100866, No Iterations 3
Min/max T:294.424 307.628
GAMG: Solving for p_rgh, Initial residual = 0.970015, Final residual = 0.00411204, No Iterations 7
time step continuity errors : sum local = 0.0211098, global = -2.13654e-18, cumulative = 1.93753e-19
Min/max rho:1.1284 1.18805
Solving for solid region KK
DICPCG: Solving for h, Initial residual = 0.209534, Final residual = 0.00386809, No Iterations 2
Min/max T:min(T) [0 0 0 1 0 0 0] 296.966 max(T) [0 0 0 1 0 0 0] 360
ExecutionTime = 2.47 s ClockTime = 3 s
Time = 3
Solving for fluid region bottomAir
DILUPBiCG: Solving for Ux, Initial residual = 0.471824, Final residual = 0.0371552, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.521667, Final residual = 0.0517465, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.470758, Final residual = 0.0320942, No Iterations 2
DILUPBiCG: Solving for h, Initial residual = 0.479107, Final residual = 0.0243738, No Iterations 2
Min/max T:248.206 316.052
GAMG: Solving for p_rgh, Initial residual = 0.717735, Final residual = 0.00578444, No Iterations 6
time step continuity errors : sum local = 0.0474212, global = 1.31978e-18, cumulative = 1.51353e-18
Min/max rho:0.859353 1.66003
Solving for solid region KK
DICPCG: Solving for h, Initial residual = 0.109024, Final residual = 0.00190127, No Iterations 2
Min/max T:min(T) [0 0 0 1 0 0 0] 296.863 max(T) [0 0 0 1 0 0 0] 360
ExecutionTime = 3.05 s ClockTime = 3 s
Time = 4
Solving for fluid region bottomAir
DILUPBiCG: Solving for Ux, Initial residual = 0.423113, Final residual = 0.0375891, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.725986, Final residual = 0.0713627, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.420529, Final residual = 0.0279138, No Iterations 2
DILUPBiCG: Solving for h, Initial residual = 0.710462, Final residual = 0.0287826, No Iterations 2
Min/max T:-687.938 321.57
GAMG: Solving for p_rgh, Initial residual = 0.791067, Final residual = 0.00474987, No Iterations 5
time step continuity errors : sum local = 0.113992, global = -7.09645e-19, cumulative = 8.03888e-19
Min/max rho:0.2 2
Solving for solid region KK
DICPCG: Solving for h, Initial residual = 0.0719346, Final residual = 0.00125857, No Iterations 2
Min/max T:min(T) [0 0 0 1 0 0 0] 295.145 max(T) [0 0 0 1 0 0 0] 360
ExecutionTime = 3.64 s ClockTime = 4 s
Time = 5
Solving for fluid region bottomAir
DILUPBiCG: Solving for Ux, Initial residual = 0.348038, Final residual = 0.0172571, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.480059, Final residual = 0.0462242, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.363602, Final residual = 0.0222914, No Iterations 2
DILUPBiCG: Solving for h, Initial residual = 0.446784, Final residual = 0.0196146, No Iterations 2
--> FOAM FATAL ERROR:
Maximum number of iterations exceeded
From function thermo<Thermo, Type>::T(scalar f, scalar T0, scalar (thermo<Thermo, Type>::*F)(const scalar) const, scalar (thermo<Thermo, Type>::*dFdT)(const scalar) const, scalar (thermo<Thermo, Type>::*limit)(const scalar) const) const
in file /usr2/sw/OpenFOAM//OpenFOAM-2.2.0/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 76.
FOAM aborting
#0 Foam::error: rintStack(Foam::Ostream&) in "/cax/sw/OPENFOAM/LINUX_x86_64/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::error::abort() in "/cax/sw/OPENFOAM/LINUX_x86_64/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 Foam::species::thermo<Foam::hConstThermo<Foam: erfectGas<Foam::specie> >, Foam::sensibleEnthalpy>::T(double, double, double, double (Foam::species::thermo<Foam::hConstThermo<Foam: erfectGas<Foam::specie> >, Foam::sensibleEnthalpy>::*)(double, double) const, double (Foam::species::thermo<Foam::hConstThermo<Foam: erfectGas<Foam::specie> >, Foam::sensibleEnthalpy>::*)(double, double) const, double (Foam::species::thermo<Foam::hConstThermo<Foam: erfectGas<Foam::specie> >, Foam::sensibleEnthalpy>::*)(double) const) const in "/cax/sw/OPENFOAM/LINUX_x86_64/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#3 Foam::heRhoThermo<Foam::rhoThermo, Foam: ureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam: erfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > >::calculate() in "/cax/sw/OPENFOAM/LINUX_x86_64/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#4 Foam::heRhoThermo<Foam::rhoThermo, Foam: ureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam: erfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > >::correct() in "/cax/sw/OPENFOAM/LINUX_x86_64/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#5
in "/cax/sw/OPENFOAM/LINUX_x86_64/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam"
#6 __libc_start_main in "/lib64/libc.so.6"
#7
in "/cax/sw/OPENFOAM/LINUX_x86_64/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam"
____________________________________________________________________________________________________________________
Dabei ist es klar zu sehen, dass der Konvergenz von h in der Luft instabil ist(Sehe Bild). Der Grund dafür habe ich nach mehrere Versuche indem ich viel BC benutzt, Relaxationfaktor geändert, div Scheme geändert, andere Netz benutzt... immer noch nicht gefunden.
____________________________________________________________________________________________________________________
Hier sind die 0-Ordner und system-Ordner Dateien:
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object epsilon;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -3 0 0 0 0];
internalField uniform 8.24e-7;
boundaryField
{
bottomMinX
{
type compressible::epsilonWallFunction;
value uniform 8.24e-7;
}
bottomMinY
{
type compressible::epsilonWallFunction;
value uniform 8.24e-7;
}
bottomMinZ
{
type compressible::epsilonWallFunction;
value uniform 8.24e-7;
}
bottomMaxX
{
type compressible::epsilonWallFunction;
value uniform 8.24e-7;
}
bottomMaxY
{
type compressible::epsilonWallFunction;
value uniform 8.24e-7;
}
bottomMaxZ
{
type compressible::epsilonWallFunction;
value uniform 8.24e-7;
}
bottomAir_to_KK
{
type compressible::epsilonWallFunction;
value uniform 8.24e-7;
}
}
// ************************************************************************* //
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -2 0 0 0 0];
internalField uniform 3.75e-5;
boundaryField
{
bottomMinX
{
type compressible::kqRWallFunction;
value uniform 3.75e-5;
}
bottomMinY
{
type compressible::kqRWallFunction;
value uniform 3.75e-5;
}
bottomMinZ
{
type compressible::kqRWallFunction;
value uniform 3.75e-5;
}
bottomMaxY
{
type compressible::kqRWallFunction;
value uniform 3.75e-5;
}
bottomMaxX
{
type compressible::kqRWallFunction;
value uniform 3.75e-5;
}
bottomMaxZ
{
type compressible::kqRWallFunction;
value uniform 3.75e-5;
}
bottomAir_to_KK
{
type compressible::kqRWallFunction;
value uniform 3.75e-5;
}
}
// ************************************************************************* //
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [1 -1 -2 0 0 0 0];
internalField uniform 1e5;
boundaryField
{
bottomMinX
{
type calculated;
value uniform 1e5;
}
bottomMinY
{
type calculated;
value uniform 1e5;
}
bottomMinZ
{
type calculated;
value uniform 1e5;
}
bottomMaxY
{
type calculated;
value uniform 1e5;
}
bottomMaxX
{
type calculated;
value uniform 1e5;
}
bottomMaxZ
{
type calculated;
value uniform 1e5;
}
bottomAir_to_KK
{
type calculated;
value uniform 1e5;
}
}
// ************************************************************************* //
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p_rgh;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [1 -1 -2 0 0 0 0];
internalField uniform 1e5;
boundaryField
{
bottomMinX
{
type fixedFluxPressure;
value uniform 1e5;
}
bottomMinY
{
type fixedFluxPressure;
value uniform 1e5;
}
bottomMinZ
{
type fixedFluxPressure;
value uniform 1e5;
}
bottomMaxX
{
type fixedFluxPressure;
value uniform 1e5;
}
bottomMaxY
{
type fixedFluxPressure;
value uniform 1e5;
}
bottomMaxZ
{
type fixedFluxPressure;
value uniform 1e5;
}
bottomAir_to_KK
{
type fixedFluxPressure;
value uniform 1e5;
}
}
// ************************************************************************* //
FoamFile
{
version 2.0;
format binary;
class volScalarField;
location "1";
object T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 0 1 0 0 0];
internalField uniform 297;
boundaryField
{
bottomMinX
{
type zeroGradient;
}
bottomMinY
{
type zeroGradient;
}
bottomMinZ
{
type zeroGradient;
}
bottomMaxX
{
type zeroGradient;
}
bottomMaxY
{
type zeroGradient;
}
bottomMaxZ
{
type zeroGradient;
}
bottomAir_to_KK
{
type compressible::turbulentTemperatureCoupledBaffleMixed;
value uniform 297;
neighbourFieldName T;
kappa fluidThermo;
kappaName none;
}
}
// ************************************************************************* //
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform ( 0 -0.1 0 );
boundaryField
{
bottomMinX
{
type fixedValue;
value uniform (0 0 0);
}
bottomMinY
{
type fixedValue;
value uniform (0 0 0);
}
bottomMinZ
{
type fixedValue;
value uniform (0 0 0);
}
bottomMaxX
{
type fixedValue;
value uniform (0 0 0);
}
bottomMaxY
{
type fixedValue;
value uniform (0 0 0);
}
bottomMaxZ
{
type fixedValue;
value uniform (0 0 0);
}
bottomAir_to_KK
{
type fixedValue;
value uniform (0 0 0);
}
}
// ************************************************************************* //
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
ddtSchemes
{
default steadyState;
}
gradSchemes
{
default Gauss linear;
}
divSchemes
{
default none;
div(phi,U) bounded Gauss upwind;
div(phi,K) bounded Gauss upwind;
div(phi,h) bounded Gauss upwind;
div(phi,k) bounded Gauss upwind;
div(phi,K) bounded Gauss upwind;
div(phi,epsilon) bounded Gauss upwind;
div(phi,R) bounded Gauss upwind;
div(R) Gauss linear;
div((muEff*dev2(T(grad(U))))) Gauss linear;
}
laplacianSchemes
{
default none;
laplacian(muEff,U) Gauss linear uncorrected;
laplacian(Dp,p_rgh) Gauss linear uncorrected;
laplacian(alphaEff,h) Gauss linear uncorrected;
laplacian(DkEff,k) Gauss linear uncorrected;
laplacian(DepsilonEff,epsilon) Gauss linear uncorrected;
laplacian(DREff,R) Gauss linear uncorrected;
}
interpolationSchemes
{
default linear;
}
snGradSchemes
{
default uncorrected;
}
fluxRequired
{
default no;
p_rgh;
}
// ************************************************************************* //
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
solvers
{
rho
{
solver PCG
preconditioner DIC;
tolerance 1e-7;
relTol 0;
}
p_rgh
{
solver GAMG;
tolerance 1e-7;
relTol 0.01;
smoother GaussSeidel;
cacheAgglomeration true;
nCellsInCoarsestLevel 10;
agglomerator faceAreaPair;
mergeLevels 1;
}
"(U|h|k|epsilon)"
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-7;
relTol 0.1;
}
}
SIMPLE
{
momentumPredictor on;
nNonOrthogonalCorrectors 0;
pRefCell 0;
pRefValue 100000;
rhoMin rhoMin [1 -3 0 0 0] 0.2;
rhoMax rhoMax [1 -3 0 0 0] 2;
}
relaxationFactors
{
fields
{
rho 1;
p_rgh 0.3;
}
equations
{
U 0.7;
h 0.7;
nuTilda 0.7;
k 0.3;
epsilon 0.3;
omega 0.7;
"ILambda.*" 0.7;
}
}
// ************************************************************************* //
Falls dass jemand eine Idee hat woran könnte es liegen, dass meine Berechnung immerwieder abstürtzt. Bitte Bitte mitteilen.
Ich bedanke mich im voraus.
MfG
Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP