Set Parameters of a Turbulence Model

alinaghedifar · October 29, 2019, 08:07

Hi everyone
I need to change some parameters of turbulence model (e.g. TurbSA_E_COMP) which are not normally available in config file. For example, I want change "cw3_6" parameter in "TurbSA_E_COMP". Could you please help me?
The code of "TurbSA_E_COMP" is located in "src\numerics_direct_turbulent.cpp":

CSourcePieceWise_TurbSA_E_COMP::CSourcePieceWise_T urbSA_E_COMP(unsigned short val_nDim, unsigned short val_nVar,
CConfig *config) : CNumerics(val_nDim, val_nVar, config) {

incompressible = (config->GetKind_Regime() == INCOMPRESSIBLE);
rotating_frame = config->GetRotating_Frame();

/*--- Spalart-Allmaras closure constants ---*/

cv1_3 = pow(7.1, 3.0);
k2 = pow(0.41, 2.0);
cb1 = 0.1355;
cw2 = 0.3;
ct3 = 1.2;
ct4 = 0.5;
cw3_6 = pow(2.0, 6.0);
sigma = 2./3.;
cb2 = 0.622;
cb2_sigma = cb2/sigma;
cw1 = cb1/k2+(1.0+cb2)/sigma;

}

pcg · October 30, 2019, 04:18

Hi,
Have a look at this thread on the SU2 github: https://github.com/su2code/SU2/issues/679
You might be able to use the branch @amirkb91 was working on.

October 29, 2019, 08:07	Set Parameters of a Turbulence Model	#1
alinaghedifar New Member S. Ali Naghedifar Join Date: Feb 2014 Posts: 7 Rep Power: 12	Hi everyone I need to change some parameters of turbulence model (e.g. TurbSA_E_COMP) which are not normally available in config file. For example, I want change "cw3_6" parameter in "TurbSA_E_COMP". Could you please help me? The code of "TurbSA_E_COMP" is located in "src\numerics_direct_turbulent.cpp": CSourcePieceWise_TurbSA_E_COMP::CSourcePieceWise_T urbSA_E_COMP(unsigned short val_nDim, unsigned short val_nVar, CConfig config) : CNumerics(val_nDim, val_nVar, config) { incompressible = (config->GetKind_Regime() == INCOMPRESSIBLE); rotating_frame = config->GetRotating_Frame(); /--- Spalart-Allmaras closure constants ---*/ cv1_3 = pow(7.1, 3.0); k2 = pow(0.41, 2.0); cb1 = 0.1355; cw2 = 0.3; ct3 = 1.2; ct4 = 0.5; cw3_6 = pow(2.0, 6.0); sigma = 2./3.; cb2 = 0.622; cb2_sigma = cb2/sigma; cw1 = cb1/k2+(1.0+cb2)/sigma; } Kaanytk likes this.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
interFoam wave propagation and explosion of Courant number and residuals	ChiaraViola	OpenFOAM Running, Solving & CFD	1	June 26, 2019 05:36
Pro/E to ANSYS Parameterization Guide	Trues	ANSYS	4	April 18, 2018 05:52
[snappyHexMesh] determining displacement for added points	CFDnewbie147	OpenFOAM Meshing & Mesh Conversion	1	October 22, 2013 09:53
How to decide to Turbulence model	shipman	OpenFOAM	2	August 18, 2013 03:00
Centrifugal Pump and Turbulence Model	Michiel	CFX	12	January 25, 2010 03:20

October 30, 2019, 04:18		#2
pcg Senior Member Pedro Gomes Join Date: Dec 2017 Posts: 465 Rep Power: 13	Hi, Have a look at this thread on the SU2 github: https://github.com/su2code/SU2/issues/679 You might be able to use the branch @amirkb91 was working on.