OpenFOAM® v1712: New boundary conditions

31/12/2017

New wave modelling: stream function

Integration of functionality produced by The Environmental Hydraulics Institute IHCantabria to model wave inlet conditions based on stream function theory.

The new boundary condition is specified using, e.g.:

inlet
{
    alpha           alpha.water;
    waveModel       streamFunction;
    nPaddle         1;
    waveHeight      0.1517;
    waveAngle       0.0;
    rampTime        6.034;
    activeAbsorption yes;
    wavePeriod      3.017;
    uMean           2.0825;
    waveLength      6.2832;

    Bjs
    (
        8.6669014e-002
        2.4849799e-002
        7.7446850e-003
        2.3355420e-003
        6.4497731e-004
        1.5205114e-004
        2.5433769e-005
       -2.2045436e-007
       -2.8711504e-006
       -1.2287334e-006
    );

    Ejs
    (
    5.6009609e-002
    3.1638171e-002
    1.5375952e-002
    7.1743178e-003
    3.3737077e-003
    1.6324880e-003
    8.2331980e-004
    4.4403497e-004
    2.7580059e-004
    2.2810557e-004
    );
}

These settings lead to wave profiles such as:

[Picture]

Source code
$FOAM_SRC/waveModels/waveGenerationModels/derived/streamFunction
Examples
$FOAM_TUTORIALS/multiphase/interFoam/laminar/waveExampleStreamFunction

Attribution
This boundary condition was supplied by The Environmental Hydraulics Institute IHCantabria - see commit 485ac6
Authors: Gabriel Barajas
Integration
The code has been updated by OpenCFD and added to the $FOAM_SRC/waveModels library available to the interFoam family of solvers - see commit 310756

New blended turbulence viscosity wall function

The new nutUBlendedWallFunction is based on the automatic wall treatment method described in the publication:

Menter, F., Carregal Ferreira, J., Esch, T., Konno, B. (2003). The SST Turbulence Model with Improved Wall Treatment for Heat Transfer Predictions in Gas Turbines. Proceedings of the International Gas Turbine Congress 2003 Tokyo

The turbulence viscosity is derived from the friction velocity, defined as a blend of contributions from the viscous sub-layer and logarithmic regions.

pict\relax \special {t4ht=

where the coefficient n  \relax \special {t4ht= is set to 4  \relax \special {t4ht= by default. The wall function is specified using, e.g.

myWallPatch
{
    type        nutUBlendedWallFunction;
}

Source code
$FOAM_SRC//TurbulenceModels/turbulenceModels/derivedFvPatchFields/wallFunctions/nutWallFunctions/nutUBlendedWallFunction