[driu]

Dear Foamers,

I have recently started using OpenFOAM and gradually ran into some convergence issues that I cannot tackle myself.

What I am trying to do is to compare my OpenFOAM simulation with a simulation done in ANSYS CFX.
The case is a swirl induced by turning vanes. Please find attached the geometry domain screenshot.

I am running the case with simpleFoam using kOmegaSST turbulence model.
Please see the case files attached or let me know if you require any specific details.


Now the problem is convergence in my OpenFOAM simulation. As soon as I try using second order divergence schemes, I get problems with convergence. I am quite lost where the problem is.

Please see the residual plot attached.

For comparison, the residuals in CFX go down below 1e-6 (RMS) or 1e-3 (max). I am using ‘High resolution’ Advection Scheme setting there.

I would really appreciate any guidance on this and apologies if it’s just another beginner level question.

Is there anything obvious I should check first or anything obviously wrong in the setup?

Many thanks for your time!

Kind regards,
Audrius

[driu]

Hello Everyone,

Recently started learning OpenFOAM. I am trying to define the inlet turbulence exactly how I have it in ANSYS CFX.

In CFX I have used Medium, which means 5% intensity and a viscosity ratio μ_t/μ equal to 10.

In OpenFOAM I have used Intensity and Length Scale:

For k

inlet
{
type turbulentIntensityKineticEnergyInlet;
intensity 0.05;
value $internalField;
}

For omega

inlet
{
type turbulentMixingLengthFrequencyInlet;
mixingLength 0.021; //0.07 x Hydraulic Diameter
value $internalField;
}

It's kOmegaSST turbulence model BTW.

I know the problem is in the k and omega specification because if I specify them as fixed values, OpenFOAM and CFX give identical results.

So my questions are:

1) Is it possible in OpenFOAM to specify the intensity and the viscosity ratio instead of length scale?
2) If not, can I calculate turbulent length scale from the viscosity ratio?
3) Maybe I should'n stick to the CFX's default turbulence setting and there is a better way of estimating the turbulence values for the inlet?

Many thanks,
Audrius


----------------------------


Hello Everyone,

I have just started using OpenFOAM and got into some issues trying to specify identical inlet conditions to ANSYS CFX.

In CFX I was using 'Medium' turbulence, which is 5% intensity and viscosity ratio μ_t/μ equal to 10.

In OpenFOAM I have found a way to specify the Intensity in k using turbulentIntensityKineticEnergyInlet and length scale in omega using turbulentMixingLengthFrequencyInlet. For the length scale I have specified 0.07 x Hydraulic Diameter.

Now this does not give me identical results to CFX. I know this discrepancy is because of k and omega values at the inlet. When I specify k and omega as fixed values, I get identical results in OpenFOAM and in CFX.

So my questions are:

1) Is there a way to specify viscosity ration in OpenFOAM?
2) Is there a way of estimating the length-scale from the viscosity ratio?

Many thanks,
Audrius

[Mash]

hello foamers, I'm a bit of a newbie to the software. ive been trying to use the k-omega SST model without a wall function to model the effect of a ducted tip hydrofoil, and it effects on the CL/CD. i have my Y+~ 1, but i am still wondering what values for K and Omega to use in the initial conditions.

i have found several different sources to giving different equations to calculate the these values.

currently using:
k=10e-6*U^2
Omega = 5U/L

l=chord length
this is to model the tip-vortex of a hydrofoil.
Reynolds number ~ 1.5e6

anyadvice is much appreciated

[tareqkh]

Hi mash,

If you are not using wall function i.e y+<6, the following initial condition for k and omega should be used.

k@omega.png

Note: for omega use 1e+10 instead of minus (close to the infinity) but you have the freedom to try them both.

Best,

[fredicenci]

I am trying to reproduce a transient flow around a circular cylinder. I am setting the values of k, omega and nut as the equations on this link https://www.cfd-online.com/Tools/turbulence.php, it might help you.

I normally give the following entries: Turbulent intensity, eddy viscosity ratio, free stream velocity and laminar viscosity.

I have to point out that, in my case, i am using nut/nu = 0.01, however my simulation doesnt give this value back in my mesh. I am trying to investigate why.


Bests Regards.

[Mash]

im having a similar issue, i have seen a ridiculous number of different methods for estimating these values, all of which get different answers,
im working with a 3d hydrofoil with chord 0.1m and inlet velocity of 16m/s
also finding many different formula for K and Omega.
if anyone has a method they are sure of it would be much appreciated
cheers matteo

[wyldckat]

Greetings to all!

I've merged a few threads that were all on the same topic into this single thread, given that pretty much the same answers can be given here. Therefore, for those who asked the questions, please re-read the whole thread to see if there is already an answer for your question.


As for the questions about "eddy viscosity ratio", I can't find any public implementation for OpenFOAM of a boundary condition for the "omega" field that takes as a parameter the "eddy viscosity ratio". I wasn't even aware that this was sometimes a better alternative to using mixing length.

I've done a bit of research and the following links are relevant to this topic:

• What "eddy viscosity ratio" actually is: https://www.cfd-online.com/Wiki/Eddy_viscosity_ratio
• How it's used for epsilon and omega fields: https://www.cfd-online.com/Wiki/Turb...ary_conditions

This would be a great addition for OpenFOAM, but I'm intrigued about a particular detail: Is this still also commonly used for compressible and multiphase flow?




Right now I'm unable to try and code this myself, but if anyone can provide a simple test case in which this can be verified that it's easier to use this type of boundary condition instead of using mixing length, I can try and code the boundary conditions some weekend in May 2017...


Best regards,
Bruno

[arvindpj]

Attached is the calculator, I have used for my cases.
Cheers,
J

[mahtin360]

Hello,

I am trying different settings for the k-w SST model in OpenFOAM 2.3.0 at the moment.
I have come across a number of different settings, espcially the internalField set up and started to wonder what is the best/right way?

are there any differences between using

internalField uniform $turbulentOmega;
and
internalField uniform [value];

does the second mean that omega is fixed?

Does anyone have any experience with this?

[arvindpj]

Originally Posted by mahtin360

Hello,

I am trying different settings for the k-w SST model in OpenFOAM 2.3.0 at the moment.
I have come across a number of different settings, espcially the internalField set up and started to wonder what is the best/right way?

are there any differences between using

internalField uniform $turbulentOmega;
and
internalField uniform [value];

does the second mean that omega is fixed?

Does anyone have any experience with this?

They are both the same. $NAME is just the variable name. In this case, the value is stored in $NAME.

Cheers,
Jay

[vava10]

Hey,

I am simulating a kayak for velocity 4.5 and 5.5 m/s and I am using K-omega SST. I found the formulas for finding the values as shown in the image. I am not an expert in CFD. I am getting really high value for omega. I don't know if I am making any mistakes consept.
But when I try these formula on some on the examples cases, I am not getting the correct answer

length 5.19

Following are the values

Re = 23355000
I = 0.0191
k = 0.0111
L/sqr(Re) = 1.073e-3
l =4.29e-4
omega =488

in some of the papers which I read they just assume the value of intensity I as 0.5%


can some one help?

thanks in advance
vava10