k-Epsilon model for compressible air

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 June 7, 2021, 13:23 k-Epsilon model for compressible air #1 New Member   Carlos Join Date: May 2017 Posts: 15 Rep Power: 7 Hello all, I hope that you may help me. I am building a venturing nozzle based on the "Supersonic flow over a forward-facing step "([URL="https://www.openfoam.com/documentation/tutorial-guide/3-compressible-flow/3.2-supersonic-flow-over-a-forward-facing-step"]). The code runs well for my geometry. I have added a sketch of my geometry and the names of the patches. Now, I am trying to add a k-Epsilon turbulent model to account for turbulence. However, I am struggling a lot with this. The case crushes after 4 time steps: Courant Number mean: 0.00183398 max: 0.296916 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCGStab: Solving for Ux, Initial residual = 0.848478, Final residual = 7.41534e-07, No Iterations 1000 DILUPBiCGStab: Solving for Uy, Initial residual = 0.905086, Final residual = 9.39327e-06, No Iterations 1000 DILUPBiCGStab: Solving for e, Initial residual = 0.357641, Final residual = 8.74872e-10, No Iterations 828 DILUPBiCGStab: Solving for p, Initial residual = 0.999241, Final residual = 1177.16, No Iterations 1000 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 199.308, global = 4.46821, cumulative = 4.46821 DILUPBiCGStab: Solving for p, Initial residual = 0.979981, Final residual = 7.17817e+07, No Iterations 1000 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 88132.4, global = -30178.3, cumulative = -30173.8 DILUPBiCGStab: Solving for epsilon, Initial residual = 1, Final residual = 9.78935e-06, No Iterations 1000 bounding epsilon, min: -6.89085e+36 max: 3.41706e+36 average: 4.37867e+32 DILUPBiCGStab: Solving for k, Initial residual = 0.999482, Final residual = 4.61607e-11, No Iterations 6 bounding k, min: -5.46211e+31 max: 6.0254e+29 average: -3.7871e+26 #0 Foam::error:rintStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 ? in "/lib64/libc.so.6" #3 Foam::nutUSpaldingWallFunctionFvPatchScalarField:: calcUTau(Foam::Field const&) const at ??:? #4 Foam::nutUSpaldingWallFunctionFvPatchScalarField:: calcNut() const at ??:? #5 Foam::nutWallFunctionFvPatchScalarField::updateCoe ffs() at ??:? #6 Foam::fvPatchField::evaluate(Foam::UPstrea m::commsTypes) at ??:? #7 Foam::RASModels::kEpsilon > > >::correctNut() at ??:? #8 Foam::RASModels::kEpsilon > > >::correct() at ??:? #9 ? at ??:? #10 __libc_start_main in "/lib64/libc.so.6" #11 ? at ??:? Floating point exception (core dumped) I think that my problem is in the BC for the turbulence. I have attached all the files in this message. I am wondering if someone could help me. I am not sure what I am doing wrong. Last edited by Ricardo1992; June 9, 2021 at 09:51.

 Tags compressible air, k-epsilon model, openfoam, turbulence

 Thread Tools Search this Thread Search this Thread: Advanced Search Display Modes Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are Off Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Guille1811 CFX 25 November 12, 2017 17:38 ashtonJ CFX 6 August 11, 2014 14:32 malaboss OpenFOAM Verification & Validation 1 February 1, 2013 16:36 Maicol Main CFD Forum 0 November 30, 2012 16:25 Saba FLUENT 0 February 10, 2009 12:05

All times are GMT -4. The time now is 00:47.

 Contact Us - CFD Online - Privacy Statement - Top