|
[Sponsors] |
turn off viscus heating in thermophysicalProperties |
![]() |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
![]() |
![]() |
#1 |
New Member
yahya rezvani
Join Date: Aug 2019
Posts: 18
Rep Power: 7 ![]() |
Hello. i want use openfoam to solve laminar viscos fluid with rotating wall. my problem solved in fluent by turn off viscous heating and use simple scheme for Coupling scheme & use heat source in rotating wall instead of viscus heating. but in openfoam i get error of inital temprarure.
this is my setting: Code:
ddtSchemes { default Euler; } gradSchemes { default Gauss linear; grad(U) cellLimited Gauss linear 1; } divSchemes { div(phi,alpha) Gauss interfaceCompression vanLeer 1; div((((alpha.air*thermo:rho.air)*nuEff.air)*dev2(T(grad(U))))) Gauss linear; div((((alpha.polymer*thermo:rho.polymer)*(thermo:mu.polymer|thermo:rho.polymer))*dev2(T(grad(U))))) Gauss linear; div(((alpha.polymer*thermo:rho.polymer)*sigma.polymer)) Gauss linear; div((((alpha.polymer*thermo:rho.polymer)*nuM)*grad(U))) Gauss linear; div(rhoPhi,U) Gauss linearUpwindV grad(U); div(rhoPhi,K) Gauss upwind; div((phi+meshPhi),p) Gauss upwind; div(rhoPhi,T) Gauss linearUpwind grad(T); div(alphaRhoPhi.polymer,sigma.polymer) Gauss linearUpwind grad(sigma); div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear; } laplacianSchemes { default Gauss linear limited corrected 0.33; } interpolationSchemes { default linear; } snGradSchemes { default limited corrected 0.33; } Code:
thermoType { type heRhoThermo; mixture pureMixture; transport WLF; thermo hConst; equationOfState rPolynomial; specie specie; energy sensibleEnthalpy; } mixture { specie { molWeight 291.3; } equationOfState { C (0.005218 -3.1055e-06 4.8689e-09 5.7772e-13 -3.2225e-16); } thermodynamics { Cp 1800; Hf 0; } transport { mu0 53763179097.64641453; Tr 115; C1 31.24; C2 107.53; Pr 6.15; } } Code:
Create time Create mesh for time = 0 Selecting dynamicFvMesh dynamicMotionSolverFvMesh Selecting motion solver: solidBody Selecting solid-body motion function rotatingMotion Applying solid body motion to cellZone propellerStem PIMPLE: No convergence criteria found PIMPLE: No corrector convergence criteria found Calculations will do 2 corrections PIMPLE: Operating solver in transient mode with 2 outer correctors Reading field p_rgh Reading field U Reading/calculating face flux field phi Constructing twoPhaseMixtureThermo Selecting thermodynamics package { type heRhoThermo; mixture pureMixture; transport WLF; thermo hConst; equationOfState rPolynomial; specie specie; energy sensibleEnthalpy; } Selecting thermodynamics package { type heRhoThermo; mixture pureMixture; transport const; thermo hConst; equationOfState rhoConst; specie specie; energy sensibleInternalEnergy; } Reading thermophysical properties Reading g Reading hRef Calculating field g.h Selecting turbulence model type laminar Selecting laminar stress model Stokes Creating field kinetic energy K No MRF models present Creating finite volume options from "system/fvOptions" Selecting finite volume options model type limitTemperature Source: limitT - selecting all cells - selected 1131591 cell(s) with volume 0.0005971852225 Courant Number mean: 0.001553287087 max: 0.06698099517 Constructing face velocity Uf Starting time loop Courant Number mean: 0.001553287087 max: 0.06698099517 Interface Courant Number mean: 0 max: 0 deltaT = 0.001199040767 Time = 0.00119904 PIMPLE: Iteration 1 smoothSolver: Solving for alpha.polymer, Initial residual = 4.148307397e-07, Final residual = 4.148307397e-07, No Iterations 0 Phase-1 volume fraction = 0.2935706537 Min(alpha.polymer) = 0 Max(alpha.polymer) = 1 MULES: Correcting alpha.polymer MULES: Correcting alpha.polymer Phase-1 volume fraction = 0.2935706537 Min(alpha.polymer) = 0 Max(alpha.polymer) = 1 DILUPBiCGStab: Solving for T, Initial residual = 0.7538694576, Final residual = 7.468251956e-07, No Iterations 59 DICPCG: Solving for p_rgh, Initial residual = 0.9999999577, Final residual = 14.02207328, No Iterations 1000 PIMPLE: Iteration 2 smoothSolver: Solving for alpha.polymer, Initial residual = 0.0001202851214, Final residual = 9.627369821e-07, No Iterations 18 Phase-1 volume fraction = 0.2935690002 Min(alpha.polymer) = 0 Max(alpha.polymer) = 1 MULES: Correcting alpha.polymer MULES: Correcting alpha.polymer Phase-1 volume fraction = 0.2935690002 Min(alpha.polymer) = -0.1114075347 Max(alpha.polymer) = 1.245386626 DILUPBiCGStab: Solving for T, Initial residual = 0.2150890754, Final residual = 3.474027855e-07, No Iterations 169 #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 ? in "/lib/x86_64-linux-gnu/libc.so.6" #3 ? in "/lib/x86_64-linux-gnu/libm.so.6" #4 exp in "/lib/x86_64-linux-gnu/libm.so.6" #5 Foam::heRhoThermo<Foam::rhoThermo, Foam::pureMixture<Foam::WLFTransport<Foam::species::thermo<Foam::hConstThermo<Foam::rPolynomial<Foam::specie> >, Foam::sensibleEnthalpy> > > >::calculate() at ??:? #6 Foam::heRhoThermo<Foam::rhoThermo, Foam::pureMixture<Foam::WLFTransport<Foam::species::thermo<Foam::hConstThermo<Foam::rPolynomial<Foam::specie> >, Foam::sensibleEnthalpy> > > >::correct() at ??:? #7 Foam::twoPhaseMixtureThermo::correctThermo() at ??:? #8 ? in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/compressibleInterFoam" #9 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #10 ? in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/compressibleInterFoam" Floating point exception (core dumped) best regards. |
|
![]() |
![]() |
![]() |
![]() |
#2 |
New Member
atilla altintas
Join Date: Oct 2010
Posts: 7
Rep Power: 16 ![]() |
Hi,
It could be a bit late but I had a similar problem and solved it with below change in the thermophysicalProperties dictionary: thermoType { type heSolidThermo; mixture pureMixture; transport constIso; thermo hConst; equationOfState rhoConst; specie specie; energy sensibleEnthalpy; } changed to: thermoType { type heRhoThermo; mixture pureMixture; transport const; thermo hConst; equationOfState rhoConst; specie specie; energy sensibleEnthalpy; } |
|
![]() |
![]() |
![]() |
Thread Tools | Search this Thread |
Display Modes | |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
Sudden stop of heating - convergence problem? | PSander | OpenFOAM Running, Solving & CFD | 0 | February 2, 2021 05:34 |
Comparison of under-floor heating and radiator heating | Mountainbike | OpenFOAM Running, Solving & CFD | 1 | November 27, 2020 07:55 |
Viscous Heating | aeroramji@gmail.com | FLUENT | 0 | January 10, 2017 02:20 |
is it possible to modelize induction heating in a frequency domain with OpenFOAM? | julien23 | OpenFOAM | 0 | January 12, 2015 06:47 |
Heating value of fuel in non-premixed model? | Jeronimo97 | FLUENT | 0 | February 15, 2013 15:27 |