wallHeatFlux error
1 Attachment(s)
I'm trying to simulate a forced convection problem using the buoyantBoussinesqSimpleFoam solver.
When I run the wallHeatFlux utility, I get this error: Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // contents of the constant folder attached. Thanks |
I am confused: You use buoyantBoussinesqSimpleFoam but with compressible turbulence model/wall function? Have a look at the tutorials for this solver. If I remember correctly it only uses an expansion factor to calculate the density change by temperature but not a full thermophysical model.
|
Dear Raghav
I am having same error. Can you please advise how you solved it. I shall be grateful. Thanks |
Hi Can any one please help me to overcome it.. Thank you very much.
Overall domain bounding box (0 0 0) (1 0.4 0.1) Error: Time = 10 GAMG: Solving for Ux, Initial residual = 0.424802, Final residual = 0.017996, No Iterations 1 GAMG: Solving for Uy, Initial residual = 0.332272, Final residual = 0.0140346, No Iterations 1 GAMG: Solving for Uz, Initial residual = 0.559406, Final residual = 0.0285128, No Iterations 1 GAMG: Solving for e, Initial residual = 0.574963, Final residual = 0.0496275, No Iterations 1 #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 ? in "/lib/x86_64-linux-gnu/libc.so.6" #3 Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam:: species::thermo<Foam::hConstThermo<Foam::perfectGa s<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() at ??:? #4 Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam:: species::thermo<Foam::hConstThermo<Foam::perfectGa s<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct() at ??:? #5 ? in "/opt/openfoam6/platforms/linux64GccDPInt32Opt/bin/rhoSimpleFoam" #6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #7 ? in "/opt/openfoam6/platforms/linux64GccDPInt32Opt/bin/rhoSimpleFoam" Floating point exception (core dumped) /*---------------------------------------------------------------------------------*/ Model: rhoSimpleFoam /*---------------------------------------------------------------------------------*/ Boundary Condition: 1. Pressure frontandback { type symmetry; } inlet { type zeroGradient; //type mixed; //refValue uniform 100000; //refGradient uniform 0; //valueFraction uniform 0.3; } outlet { type fixedValue; value uniform 100000; } wall { type zeroGradient; } 2. Velocity frontandback { type symmetry; } inlet { type fixedValue; value uniform (5 0 0); rhoInlet 1.5; } outlet { type zeroGradient; } wall { type noSlip; } 3. Temperature frontandback { type symmetry; } inlet { type fixedValue; value uniform 300; } outlet { type zeroGradient; } wall { type zeroGradient; } 4. k frontandback { type symmetry; } inlet { type turbulentMixingLengthDissipationRateInlet; mixingLength 0.14; value uniform 0.1; } outlet { type inletOutlet; inletValue uniform 1; value uniform 1; } wall { type kqRWallFunction; value uniform 1; } 5. epsilon frontandback { type symmetry; } inlet { type turbulentMixingLengthDissipationRateInlet; mixingLength 0.14; value uniform 0.006; } outlet { type inletOutlet; inletValue uniform 0.006; value uniform 0.006; } wall { type epsilonWallFunction; Cmu 0.09; kappa 0.41; E 9.8; value uniform 0.006; } 6. nut frontandback { type symmetry; } inlet { type calculated; value uniform 0; } outlet { type calculated; value uniform 0; } wall { type nutkWallFunction; Cmu 0.09; kappa 0.41; E 9.8; value uniform 0; } 7. alphat frontandback { type symmetry; } inlet { type calculated; value uniform 0; } outlet { type calculated; value uniform 0; } wall { type compressible::alphatWallFunction; Prt 0.85; value uniform 0; } /*------------------------------------------------------------------------------*/ fvSolution ddtSchemes { default steadyState; } gradSchemes { default Gauss linear; } divSchemes { default none; div(phi,U) bounded Gauss upwind; div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear; div(phi,e) bounded Gauss upwind; div(phi,epsilon) bounded Gauss upwind; div(phi,k) bounded Gauss upwind; div(phid,p) Gauss upwind; div(phi,Ekp) bounded Gauss upwind; div((phi|interpolate(rho)),p) Gauss upwind; } laplacianSchemes { default Gauss linear corrected; } interpolationSchemes { default linear; } snGradSchemes { default corrected; } /*----------------------------------------------------------------------*/ fvSolution solvers { p { solver GAMG; tolerance 1e-08; relTol 0.1; smoother GaussSeidel; nCellsInCoarsestLevel 20; } "(U|e|k|epsilon)" { solver GAMG; tolerance 1e-08; relTol 0.1; smoother GaussSeidel; nCellsInCoarsestLevel 20; } } SIMPLE { nNonOrthogonalCorrectors 0; rhoMin 0.1; rhoMax 1.5; //pMaxFactor 2; //pMinFactor 0.1; transonic yes; consistent yes; residualControl { p 1e-3; U 1e-4; e 1e-3; // possibly check turbulence fields "(k|epsilon|omega)" 1e-3; } } relaxationFactors { fields { p 0.3; rho 0.7; } equations { p 0.7; U 0.7; e 0.7; k 0.7; epsilon 0.7; } } |
All times are GMT -4. The time now is 10:27. |