Skip to content
Snippets Groups Projects
Commit 4c6b7217 authored by Kutalmış Berçin's avatar Kutalmış Berçin Committed by Andrew Heather
Browse files

ENH: KinematicWeberNumber: add new cloud function object

parent 9c8c69a5
No related branches found
No related tags found
1 merge request!648ENH: KinematicWeberNumber: New cloud function object
...@@ -45,6 +45,7 @@ License ...@@ -45,6 +45,7 @@ License
#include "RemoveParcels.H" #include "RemoveParcels.H"
#include "VoidFraction.H" #include "VoidFraction.H"
#include "KinematicReynoldsNumber.H" #include "KinematicReynoldsNumber.H"
#include "KinematicWeberNumber.H"
#include "ParticleDose.H" #include "ParticleDose.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
...@@ -67,6 +68,7 @@ License ...@@ -67,6 +68,7 @@ License
makeCloudFunctionObjectType(RemoveParcels, CloudType); \ makeCloudFunctionObjectType(RemoveParcels, CloudType); \
makeCloudFunctionObjectType(VoidFraction, CloudType); \ makeCloudFunctionObjectType(VoidFraction, CloudType); \
makeCloudFunctionObjectType(KinematicReynoldsNumber, CloudType); \ makeCloudFunctionObjectType(KinematicReynoldsNumber, CloudType); \
makeCloudFunctionObjectType(KinematicWeberNumber, CloudType); \
makeCloudFunctionObjectType(ParticleDose, CloudType); makeCloudFunctionObjectType(ParticleDose, CloudType);
......
/*---------------------------------------------------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | www.openfoam.com
\\/ M anipulation |
-------------------------------------------------------------------------------
Copyright (C) 2023 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
OpenFOAM is free software: you can redistribute it and/or modify it
under the terms of the GNU General Public License as published by
the Free Software Foundation, either version 3 of the License, or
(at your option) any later version.
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
for more details.
You should have received a copy of the GNU General Public License
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
\*---------------------------------------------------------------------------*/
#include "KinematicWeberNumber.H"
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
template<class CloudType>
Foam::KinematicWeberNumber<CloudType>::KinematicWeberNumber
(
const dictionary& dict,
CloudType& owner,
const word& modelName
)
:
CloudFunctionObject<CloudType>(dict, owner, modelName, typeName),
sigma_(dict.getScalar("sigma"))
{}
template<class CloudType>
Foam::KinematicWeberNumber<CloudType>::KinematicWeberNumber
(
const KinematicWeberNumber<CloudType>& we
)
:
CloudFunctionObject<CloudType>(we),
sigma_(we.sigma_)
{}
// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
template<class CloudType>
void Foam::KinematicWeberNumber<CloudType>::postEvolve
(
const typename parcelType::trackingData& td
)
{
auto& c = this->owner();
auto* resultPtr = c.template getObjectPtr<IOField<scalar>>("We");
if (!resultPtr)
{
resultPtr = new IOField<scalar>
(
IOobject
(
"We",
c.time().timeName(),
c,
IOobject::NO_READ,
IOobject::NO_WRITE,
IOobject::REGISTER
)
);
resultPtr->store();
}
auto& We = *resultPtr;
We.resize(c.size());
label parceli = 0;
forAllConstIters(c, parcelIter)
{
const parcelType& p = parcelIter();
We[parceli++] = p.We(td, sigma_);
}
const bool haveParticles = c.size();
if (c.time().writeTime() && returnReduceOr(haveParticles))
{
We.write(haveParticles);
}
}
// ************************************************************************* //
/*---------------------------------------------------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | www.openfoam.com
\\/ M anipulation |
-------------------------------------------------------------------------------
Copyright (C) 2023 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
OpenFOAM is free software: you can redistribute it and/or modify it
under the terms of the GNU General Public License as published by
the Free Software Foundation, either version 3 of the License, or
(at your option) any later version.
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
for more details.
You should have received a copy of the GNU General Public License
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
Class
Foam::KinematicWeberNumber
Group
grpLagrangianIntermediateFunctionObjects
Description
Calculates and writes particle Weber number field on the cloud.
\f[
\mathrm{We}_p =
\frac{\rho_c \, | \mathbf{u}_\mathrm{rel} |^2 \, d_p }{\sigma}
\f]
\vartable
\mathrm{We}_p | Particle Weber number
\rho_c | Density of carrier
d_p | Particle diameter
\mathbf{u}_\mathrm{rel} | Relative velocity between particle and carrier
\endvartable
Operands:
\table
Operand | Type | Location
input | - | -
output file | - | -
output field | scalarField | \<time\>/lagrangian/\<cloud\>/We
\endtable
Usage
Minimal example by using \c constant/\<CloudProperties\>:
\verbatim
cloudFunctions
{
KinematicWeberNumber1
{
// Mandatory entries
type WeberNumber;
sigma <scalar>;
}
}
\endverbatim
where the entries mean:
\table
Property | Description | Type | Reqd | Deflt
type | Type name: WeberNumber | word | yes | -
sigma | Surface tension [N/m] | scalar | yes | -
\endtable
See also
- Foam::ReactingWeberNumber
SourceFiles
KinematicWeberNumber.C
\*---------------------------------------------------------------------------*/
#ifndef KinematicWeberNumber_H
#define KinematicWeberNumber_H
#include "CloudFunctionObject.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
namespace Foam
{
/*---------------------------------------------------------------------------*\
Class KinematicWeberNumber Declaration
\*---------------------------------------------------------------------------*/
template<class CloudType>
class KinematicWeberNumber
:
public CloudFunctionObject<CloudType>
{
// Private Typedefs
//- Convenience typedef for parcel type
typedef typename CloudType::parcelType parcelType;
// Private Data
//- Surface tension [N/m]
scalar sigma_;
public:
//- Runtime type information
TypeName("WeberNumber");
// Generated Methods
//- No copy assignment
void operator=(const KinematicWeberNumber<CloudType>&) = delete;
// Constructors
//- Construct from dictionary
KinematicWeberNumber
(
const dictionary& dict,
CloudType& owner,
const word& modelName
);
//- Copy construct
KinematicWeberNumber(const KinematicWeberNumber<CloudType>& vf);
//- Construct and return a clone
virtual autoPtr<CloudFunctionObject<CloudType>> clone() const
{
return autoPtr<CloudFunctionObject<CloudType>>
(
new KinematicWeberNumber<CloudType>(*this)
);
}
//- Destructor
virtual ~KinematicWeberNumber() = default;
// Member Functions
//- Post-evolve hook
virtual void postEvolve(const typename parcelType::trackingData& td);
};
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
} // End namespace Foam
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
#ifdef NoRepository
#include "KinematicWeberNumber.C"
#endif
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
#endif
// ************************************************************************* //
...@@ -102,7 +102,13 @@ subModels ...@@ -102,7 +102,13 @@ subModels
cloudFunctions cloudFunctions
{} {
KinematicWeberNumber1
{
type WeberNumber;
sigma 1.0;
}
}
// ************************************************************************* // // ************************************************************************* //
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment