View Issue Details
| ID | Project | Category | View Status | Date Submitted | Last Update |
|---|---|---|---|---|---|
| 0004131 | OpenFOAM | Patch | public | 2024-08-12 07:11 | 2024-08-12 10:45 |
| Reporter | tniemi | Assigned To | henry | ||
| Priority | normal | Severity | minor | Reproducibility | always |
| Status | resolved | Resolution | fixed | ||
| Product Version | dev | ||||
| Fixed in Version | dev | ||||
| Summary | 0004131: Patch for multiphaseEulerFoam interfaceTurbulenceDamping | ||||
| Description | Hi, I have attached a patch which applies the same changes to multiphaseEulerFoam interfaceTurbulenceDamping as was applied to VOF-variants in this commit: https://github.com/OpenFOAM/OpenFOAM-dev/commit/7b6aaf1b58975c73692f8ce6168e757fd732c470 This allows to run a case without having to explicitly specify turbulence coeffs. | ||||
| Steps To Reproduce | Take eg. tutorials/multiphaseEuler/bubblePipe and add interfaceTurbulenceDamping { type interfaceTurbulenceDamping; libs ("libmultiphaseEulerFvModels.so"); phase water; delta 1e-4; } to fvModels. | ||||
| Tags | No tags attached. | ||||
|
|
damping.diff (1,759 bytes)
diff --git a/applications/modules/multiphaseEuler/fvModels/interfaceTurbulenceDamping/interfaceTurbulenceDamping.C b/applications/modules/multiphaseEuler/fvModels/interfaceTurbulenceDamping/interfaceTurbulenceDamping.C
index 21272fb4c5..882fe81c71 100644
--- a/applications/modules/multiphaseEuler/fvModels/interfaceTurbulenceDamping/interfaceTurbulenceDamping.C
+++ b/applications/modules/multiphaseEuler/fvModels/interfaceTurbulenceDamping/interfaceTurbulenceDamping.C
@@ -195,9 +195,9 @@ Foam::fv::interfaceTurbulenceDamping::interfaceTurbulenceDamping
(
mesh.lookupType<phaseCompressible::momentumTransportModel>(phaseName_)
),
- C2_("C2", dimless, 0),
- betaStar_("betaStar", dimless, 0),
- beta_("beta", dimless, 0)
+ C2_("C2", dimless, dict, 1.92),
+ betaStar_("betaStar", dimless, dict, 0.09),
+ beta_("beta", dimless, dict, 0.072)
{
const word epsilonName(IOobject::groupName("epsilon", phaseName_));
const word omegaName(IOobject::groupName("omega", phaseName_));
@@ -205,23 +205,10 @@ Foam::fv::interfaceTurbulenceDamping::interfaceTurbulenceDamping
if (mesh.foundObject<volScalarField>(epsilonName))
{
fieldName_ = epsilonName;
- C2_.read(turbulence_.coeffDict());
}
else if (mesh.foundObject<volScalarField>(omegaName))
{
fieldName_ = omegaName;
- betaStar_.read(turbulence_.coeffDict());
-
- // Read beta for k-omega models or beta1 for k-omega SST
- if (turbulence_.coeffDict().found("beta"))
- {
- beta_.read(turbulence_.coeffDict());
- }
- else
- {
- beta_ =
- dimensionedScalar("beta1", dimless, turbulence_.coeffDict());
- }
}
else
{
|
|
|
Resolved by commit c92b91f63187bef2509eb9530922271695bc2ff4 |
| Date Modified | Username | Field | Change |
|---|---|---|---|
| 2024-08-12 07:11 | tniemi | New Issue | |
| 2024-08-12 07:11 | tniemi | File Added: damping.diff | |
| 2024-08-12 10:45 | henry | Assigned To | => henry |
| 2024-08-12 10:45 | henry | Status | new => resolved |
| 2024-08-12 10:45 | henry | Resolution | open => fixed |
| 2024-08-12 10:45 | henry | Fixed in Version | => dev |
| 2024-08-12 10:45 | henry | Note Added: 0013360 |
