Loading...
Searching...
No Matches
OwenRyleyModel Class Reference

Computes film-separation properties from round edges for full separation (Owen & Ryley, 1983). The model assesses the film curvature based on the mesh geometry and flow field, and then calculates a force balance of the following form: More...

#include <OwenRyleyModel.H>

Inheritance diagram for OwenRyleyModel:
Collaboration diagram for OwenRyleyModel:

Public Member Functions

 TypeName ("OwenRyley")
 Runtime type information.
 OwenRyleyModel (const regionModels::areaSurfaceFilmModels::liquidFilmBase &film, const dictionary &dict)
 Construct from the base film model and dictionary.
virtual ~OwenRyleyModel ()=default
scalar minHbyR1 () const noexcept
scalar minInvR1 () const noexcept
virtual tmp< scalarFieldseparatedMassRatio () const
 Calculate the mass ratio of film separation.
Public Member Functions inherited from filmSeparationModel
 TypeName ("filmSeparationModel")
 Runtime type information.
 declareRunTimeSelectionTable (autoPtr, filmSeparationModel, dictionary,(const regionModels::areaSurfaceFilmModels::liquidFilmBase &film, const dictionary &dict),(film, dict))
 filmSeparationModel (const filmSeparationModel &)=delete
 No copy construct.
void operator= (const filmSeparationModel &)=delete
 No copy assignment.
 filmSeparationModel (const regionModels::areaSurfaceFilmModels::liquidFilmBase &film, const dictionary &dict)
 Construct from the base film model and dictionary.
virtual ~filmSeparationModel ()=default
 Destructor.
const regionModels::areaSurfaceFilmModels::liquidFilmBasefilm () const
 Return const access to the film properties.
const faMeshmesh () const noexcept
 Return const access to the finite-area mesh.

Additional Inherited Members

Static Public Member Functions inherited from filmSeparationModel
static autoPtr< filmSeparationModelNew (const regionModels::areaSurfaceFilmModels::liquidFilmBase &film, const dictionary &dict)
 Return a reference to the selected filmSeparation model.

Detailed Description

Computes film-separation properties from round edges for full separation (Owen & Ryley, 1983). The model assesses the film curvature based on the mesh geometry and flow field, and then calculates a force balance of the following form:

\‍[    F_{net} = F_{inertial} + F_{body} + F_{surface}
\‍]

with:

\‍[    F_{inertial} = -\frac{4}{3} \frac{h \, \rho \, |\mathbf{u}|^2}{R_1}
\‍]

\‍[    F_{body} =
        -\frac{\rho |\mathbf{g}| (R_1^2 - R_2^2) cos(\theta)}{2 \, R_1}
\‍]

\‍[    F_{surface} = \frac{\sigma}{R_2}
\‍]

where:

$      F_{net}       $=Magnitude of net force on the film
$      F_{inertial}  $=Magnitude of inertial forces on the film
$      F_{body}      $=Magnitude of body forces on the film
$      F_{surface}   $=Magnitude of surface forces on the film
$      h             $=Film thickness
$      \rho          $=Film density
$      \sigma        $=Film surface tension
$      \mathbf{u}    $=Film velocity
$      \mathbf{g}    $=Gravitational accelaration
$      \theta $=Approximate angle between gravity vector and outflow direction
$      R_1           $=Radius of curvature
$      R_2           $=Sum of the radius of curvature and film thickness

The film is considered separated when $F_{net}<0$; otherwise, it remains attached.

Reference:

    Governing equations (tag:OR):
        Owen, I., & Ryley, D. J. (1983).
        The flow of thin liquid films around corners.
        International journal of multiphase flow, 11(1), 51-62.
        DOI:10.1016/0301-9322(85)90005-9
Usage
Minimal example in boundary-condition files:
filmSeparationCoeffs
{
    // Mandatory entries
    model              OwenRyley;

    // Optional entries
    fThreshold         <scalar>;
    definedPatchRadii  <scalar>;
    deltaByR1Min       <scalar>;
    minInvR1           <scalar>;
}

where the entries mean:

Property Description Type Reqd Deflt
model Model name: OwenRyley word yes -
fThreshold Threshold force for separation scalar no 1e-8
definedPatchRadii Patch radius scalar no 0
deltaByR1Min Minimum gravity driven film thickness scalar no 0
minInvR1 Minimum inv R1 for separation scalar no 5
Note
  • The assumptions underlying the derivation of the model indicate that the model is better suited for round corners than for sharp corners.
Source files

Definition at line 214 of file OwenRyleyModel.H.

Constructor & Destructor Documentation

◆ OwenRyleyModel()

OwenRyleyModel ( const regionModels::areaSurfaceFilmModels::liquidFilmBase & film,
const dictionary & dict )

Construct from the base film model and dictionary.

Definition at line 224 of file OwenRyleyModel.C.

References dict, e, filmSeparationModel::film(), filmSeparationModel::filmSeparationModel(), and Foam::mag().

Here is the call graph for this function:

◆ ~OwenRyleyModel()

virtual ~OwenRyleyModel ( )
virtualdefault

Member Function Documentation

◆ TypeName()

TypeName ( "OwenRyley" )

Runtime type information.

References dict, and filmSeparationModel::film().

Here is the call graph for this function:

◆ minHbyR1()

scalar minHbyR1 ( ) const
inlinenoexcept

Definition at line 298 of file OwenRyleyModel.H.

References Foam::noexcept.

◆ minInvR1()

scalar minInvR1 ( ) const
inlinenoexcept

Definition at line 301 of file OwenRyleyModel.H.

References Foam::noexcept.

◆ separatedMassRatio()

tmp< scalarField > separatedMassRatio ( ) const
virtual

Calculate the mass ratio of film separation.

Implements filmSeparationModel.

Definition at line 248 of file OwenRyleyModel.C.

References tmp< T >::cref(), filmSeparationModel::film(), forAll, filmSeparationModel::mesh(), tmp< T >::New(), regionFaModel::regionMesh(), and Foam::Zero.

Here is the call graph for this function:

The documentation for this class was generated from the following files:
  • src/regionFaModels/liquidFilm/subModels/kinematic/injectionModel/filmSeparation/filmSeparationModels/OwenRyleyModel/OwenRyleyModel.H
  • src/regionFaModels/liquidFilm/subModels/kinematic/injectionModel/filmSeparation/filmSeparationModels/OwenRyleyModel/OwenRyleyModel.C