- **OpenFOAM Programming & Development**
(*http://www.cfd-online.com/Forums/openfoam-programming-development/*)

- - **A new turbulent model without empirical coefficients**
(*http://www.cfd-online.com/Forums/openfoam-programming-development/119879-new-turbulent-model-without-empirical-coefficients.html*)

A new turbulent model without empirical coefficients1 Attachment(s)
Hi there,
This is my first post. I read a academic paper "A Partial Average Based Study of Compressible Turbulent Flows" http://www.jomse.org/paperInfo.aspx?ID=88 For incompressible flow, the governing equations are shown in the figure. Attachment 22973 Equation (1) and (2) are averaged continuity and momentum equations. Equation (3) and (4) are fluctuation continuity and momentum equations. Equation (5) is the turbulent energy dissipation equation. The equations describe a new turbulence model without any empirical coefficients, which can simulate statistical mean behaviors and coherent structures of various benchmark turbulent flows. The chain of momentum transference contains two stages.The first stage is from the mean flow to the fluctuation flow as the viscous dissipation term of the mean flow equals the source term of the fluctuation flow. The second stage is from the fluctuation flow to molecular motion as the viscous dissipation term of the fluctuation flow plays its role. I have followed the tutorial that added temperature to icoFoam, now I plan to modify pisoFoam solver to solve the equations in OpenFoam. Any suggestion is really appreciated. Sam Chio |

As I known, the equations used by pisoFOAM solver are similar to equation (1) and (2). So I can keep the velocity predictor UEqn using the following code:
Code:
`fvVectorMatrix UEqn` |

Dear all,
I consider to build up the new turbulence solver by modifying the RASModel. As I mentioned in foregoing thread, OpenFOAM calculates the turbulent shear stresses in incompressible RASModel using the divDevReff term: Code:
`tmp<fvVectorMatrix> kEpsilon::divDevReff(volVectorField& U) const` Noting that 'nuEff' is the so-called effective viscosity that is the sum of the laminar molecular viscosity and turbulent eddy viscosity. The expression is in the form: Code:
`virtual tmp<volScalarField> nuEff() const` Sam |

Next it is required to solve drift velocity udrift, drift pressure pdrift and turbulence length scale l. Header file createField.h may take the form like:
Code:
` udrift_` |

Considering the C++ code of equation (4), we have the expression in the form
Code:
`fvVectorMatrix UdriftEqn` |

It has been a while. how is it going dude?
AFAIK, Prof. Gao's student is doing the same thing and their works had been published in Chinese Journal. |

Quote:
Code:
`+ fvc::div(phidrift, U)` |

All times are GMT -4. The time now is 13:23. |