CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Turbulence for natural external convection

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Search this Thread Display Modes
Old   April 17, 2024, 07:55
Default Turbulence for natural external convection
New Member
Join Date: Mar 2024
Location: Brazil
Posts: 7
Rep Power: 2
gushrqcosta is on a distinguished road
Hello. I'm trying to set up a case of natural external convection where a fluid flows around a hot cylinder due only to temperature differences. I have done laminar tests and it worked fine with great results.

Now I wanna add turbulence (by modifying Rayleigh parameters) but I can't figure out how to set up the initial BCs for turbulence.

My professor told me to use k-omega-SST model, but I'm not able to judge the best BCs for the files of k/omega/alphat/nut.

In my case the fluid is initially quiescent (zero velocity) and starts moving only due to the natural convection, but because of this I can't calculate what the initial values for k/omega should be, since the fluid isn't moving at the beggining.

Also, my mesh consists of:

inlet/outlet = patches
superiorPlane/inferiorPlane = symmetryPlane
cylinderWall = wall
frontAndBack = empty

I'm a begginer at OpenFOAM, so I'm sorry for stupid mistakes.

No fvModels present
No fvConstraints present
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/"
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#4 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam:perator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<doub le, Foam::fvPatchField, Foam::volMesh> > const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
#5 Foam::fluidThermo::nu() const at ??:?
#6 Foam::RASModel<Foam::compressibleMomentumTransport Model>::nu() const at ??:?
#7 Foam::kOmegaSST<Foam::eddyViscosity<Foam::RASModel <Foam::compressibleMomentumTransportModel> >, Foam::compressibleMomentumTransportModel>::F2() const at ??:?
#8 Foam::kOmegaSST<Foam::eddyViscosity<Foam::RASModel <Foam::compressibleMomentumTransportModel> >, Foam::compressibleMomentumTransportModel>::F23() const at ??:?
#9 Foam::kOmegaSST<Foam::eddyViscosity<Foam::RASModel <Foam::compressibleMomentumTransportModel> >, Foam::compressibleMomentumTransportModel>::correct Nut() at ??:?
#10 ? in "/home/workstation1/OpenFOAM-10/OpenFOAM-10/platforms/linux64GccDPInt32Opt/bin/buoyantFoam"
#11 __libc_start_main in "/lib/x86_64-linux-gnu/"
#12 ? in "/home/workstation1/OpenFOAM-10/OpenFOAM-10/platforms/linux64GccDPInt32Opt/bin/buoyantFoam"
Floating point exception (core dumped)
Attached Files
File Type: zip (5.5 KB, 1 views)
File Type: zip (2.9 KB, 0 views)
File Type: zip (2.3 KB, 0 views)

Last edited by gushrqcosta; April 17, 2024 at 10:01.
gushrqcosta is offline   Reply With Quote

Old   April 17, 2024, 10:40
Senior Member
Join Date: Dec 2021
Posts: 211
Rep Power: 5
Alczem is on a distinguished road

It looks like a division by 0 occured somewhere!

Try to avoir 0 values in your turbulent boundary conditions (you can use 1e-8 for instance) and in your fluid properties.

Concerning initial values for turbulence, I guess you could estimate the max velocity reached once your convective pattern is established (a few m/s?) and get approximative values for k and omega. Initial values are just that, for the initial steps of the simulation. The turbulent quantities should converge towards realistic values afterwars.
Alczem is offline   Reply With Quote


buoyantpimplefoam, komegasst, natural convection, turbulence

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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 Off
Trackbacks are Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Rotating body with internal and external convection RomuloSc FLUENT 2 October 27, 2017 03:08
Natural convection with or without turbulence model chrisf90 FLUENT 0 June 30, 2016 03:26
Natural convection - Inlet boundary condition max91 CFX 1 July 29, 2008 20:28
Appropiate turbulence model for natural convection Diego Peinado FLUENT 0 August 10, 2005 07:24
turbulence modeling for external aerodynamics Charles Renard Main CFD Forum 3 June 27, 2001 07:55

All times are GMT -4. The time now is 07:43.