|
[Sponsors] |
Possibility of local mechanical turbulence implementation |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
March 12, 2013, 18:23 |
Possibility of local mechanical turbulence implementation
|
#1 |
Senior Member
Canakkale Dardanelspor
Join Date: Aug 2012
Posts: 135
Rep Power: 13 |
Hi,
Is there any possibility to add turbulent kinetic energy source term, k, for only specific part of any domain different from the turbulent kinetic energy value of the rest of the domain? == 1 == Turbulent kinetic energy equation in kOmegaSST model of OpenFOAM 2.1.1 under kOmegaSST.C: Code:
// Turbulent kinetic energy equation tmp<fvScalarMatrix> kEqn ( fvm::ddt(k_) + fvm::div(phi_, k_) - fvm::Sp(fvc::div(phi_), k_) - fvm::laplacian(DkEff(F1), k_) == min(G, c1_*betaStar_*k_*omega_) - fvm::Sp(betaStar_*omega_, k_) ); An arbitrary specified circular region in the computational domain declaration is defined under Member functions: of kOmegaSST.C: Code:
const scalar xo = 0; const scalar yo = 0; const scalar zmin = -205.5; const scalar zmax = -204.5; const scalar radius1 = 10; const scalar radius2 = 5; scalar addedTurbulence; // Assumed in units m^2 / s^3. forAll(mesh_.cells(),celll) { const vector& cellCenter=mesh_.C()[celll]; const scalar dx = cellCenter[0] - xo; const scalar dy = cellCenter[1] -yo; const vector r(dx, dy, 0); if ((cellCenter[2] < zmax) && (cellCenter[2] > zmin) && (mag(r) < tipRadius) && (mag(r) > hubRadius)) { addedTurbulence== 10; }else{ addedTurbulence== 1; } I thought to add the addedTurbulence term would be adequate which is obviously a wrong approach: Code:
// Turbulent kinetic energy equation tmp<fvScalarMatrix> kEqn ( fvm::ddt(k_) + fvm::div(phi_, k_) - fvm::Sp(fvc::div(phi_), k_) - fvm::laplacian(DkEff(F1), k_) == min(G, c1_*betaStar_*k_*omega_) - fvm::Sp(betaStar_*omega_, k_) + addedTurbulence addedTurbulence must be defined as volScalarField, I appreciate that. Nevertheless, I couldn't make it work with it; therefore, I tried the above approach which considers addedTurbulence term is scalar. == 5 == Appreciate any piece of advice or suggestion. Many thanks in advance. |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Multiple floating objects | CKH | OpenFOAM Running, Solving & CFD | 14 | February 20, 2019 09:08 |
Discussion: Reason of Turbulence!! | Wen Long | Main CFD Forum | 3 | May 15, 2009 09:52 |
[blockMesh] BlockMeshmergePatchPairs | hjasak | OpenFOAM Meshing & Mesh Conversion | 11 | August 15, 2008 07:36 |
A reference on implementation of Spalart-Alam. Turbulence Model? | Mohammad Kermani | Main CFD Forum | 2 | December 26, 1999 02:56 |
turbulence modeling questions | llowen | Main CFD Forum | 3 | September 11, 1998 04:24 |