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...

Inheritance diagram for OwenRyleyModel:

[legend]

Collaboration diagram for OwenRyleyModel:

[legend]

Public Member Functions
	TypeName ("OwenRyley")
	Runtime type information. More...

	OwenRyleyModel (const regionModels::areaSurfaceFilmModels::liquidFilmBase &film, const dictionary &dict)
	Construct from the base film model and dictionary. More...

virtual	~OwenRyleyModel ()=default

scalar	minHbyR1 () const noexcept

scalar	minInvR1 () const noexcept

virtual tmp< scalarField >	separatedMassRatio () const
	Calculate the mass ratio of film separation. More...

Public Member Functions inherited from filmSeparationModel
	TypeName ("filmSeparationModel")
	Runtime type information. More...

	declareRunTimeSelectionTable (autoPtr, filmSeparationModel, dictionary,(const regionModels::areaSurfaceFilmModels::liquidFilmBase &film, const dictionary &dict),(film, dict))

	filmSeparationModel (const filmSeparationModel &)=delete
	No copy construct. More...

void	operator= (const filmSeparationModel &)=delete
	No copy assignment. More...

	filmSeparationModel (const regionModels::areaSurfaceFilmModels::liquidFilmBase &film, const dictionary &dict)
	Construct from the base film model and dictionary. More...

virtual	~filmSeparationModel ()=default
	Destructor. More...

const regionModels::areaSurfaceFilmModels::liquidFilmBase &	film () const
	Return const access to the film properties. More...

const faMesh &	mesh () const noexcept
	Return const access to the finite-area mesh. More...

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

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
	=	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
	=	Radius of curvature
	=	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 225 of file OwenRyleyModel.C.

References Foam::exit(), Foam::FatalError, FatalErrorInFunction, liquidFilmBase::g(), Foam::mag(), and dimensioned< Type >::value().

Here is the call graph for this function:

◆ ~OwenRyleyModel()

virtual ~OwenRyleyModel ( )

virtualdefault

Member Function Documentation

◆ TypeName()

TypeName ( "OwenRyley" )

Runtime type information.

◆ minHbyR1()

scalar minHbyR1 ( ) const

inlinenoexcept

Definition at line 298 of file OwenRyleyModel.H.

◆ minInvR1()

scalar minInvR1 ( ) const

inlinenoexcept

Definition at line 301 of file OwenRyleyModel.H.

◆ 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(), faMesh::nFaces(), regionFaModel::regionMesh(), and Foam::Zero.