mirror of
https://develop.openfoam.com/Development/openfoam.git
synced 2025-11-28 03:28:01 +00:00
ENH: electrostaticDeposition: new finiteVolume boundary condition
This commit is contained in:
Kutalmis Bercin
committed by
Andrew Heather
Andrew Heather
parent
aaeddba466
commit
19f7825a04
@ -241,6 +241,7 @@ $(derivedFvPatchFields)/swirlFanVelocity/swirlFanVelocityFvPatchField.C
|
||||
$(derivedFvPatchFields)/acousticWaveTransmissive/acousticWaveTransmissiveFvPatchFields.C
|
||||
$(derivedFvPatchFields)/prghPermeableAlphaTotalPressure/prghPermeableAlphaTotalPressureFvPatchScalarField.C
|
||||
$(derivedFvPatchFields)/pressurePermeableAlphaInletOutletVelocity/pressurePermeableAlphaInletOutletVelocityFvPatchVectorField.C
|
||||
$(derivedFvPatchFields)/electrostaticDeposition/electrostaticDepositionFvPatchScalarField.C
|
||||
|
||||
$(derivedFvPatchFields)/mappedMixed/mappedMixedFvPatchFields.C
|
||||
$(derivedFvPatchFields)/mappedField/Sampled/makeSampledPatchFunction1s.C
|
||||
|
||||
@ -0,0 +1,593 @@
|
||||
/*---------------------------------------------------------------------------*\
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | www.openfoam.com
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
Copyright (C) 2021 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 "electrostaticDepositionFvPatchScalarField.H"
|
||||
#include "addToRunTimeSelectionTable.H"
|
||||
#include "fvPatchFieldMapper.H"
|
||||
#include "volFields.H"
|
||||
#include "surfaceFields.H"
|
||||
|
||||
// * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * //
|
||||
|
||||
Foam::electrostaticDepositionFvPatchScalarField&
|
||||
Foam::electrostaticDepositionFvPatchScalarField::eVPatch
|
||||
(
|
||||
const label patchi
|
||||
) const
|
||||
{
|
||||
const auto& eV =
|
||||
db().lookupObject<volScalarField>(this->internalField().name());
|
||||
|
||||
const volScalarField::Boundary& bf = eV.boundaryField();
|
||||
|
||||
const auto& eVpf =
|
||||
refCast<const electrostaticDepositionFvPatchScalarField>(bf[patchi]);
|
||||
|
||||
return const_cast<electrostaticDepositionFvPatchScalarField&>(eVpf);
|
||||
}
|
||||
|
||||
|
||||
void Foam::electrostaticDepositionFvPatchScalarField::setMaster() const
|
||||
{
|
||||
if (master_ != -1)
|
||||
{
|
||||
return;
|
||||
}
|
||||
|
||||
const auto& eV =
|
||||
db().lookupObject<volScalarField>(this->internalField().name());
|
||||
|
||||
const volScalarField::Boundary& bf = eV.boundaryField();
|
||||
|
||||
label master = -1;
|
||||
forAll(bf, patchi)
|
||||
{
|
||||
if (isA<electrostaticDepositionFvPatchScalarField>(bf[patchi]))
|
||||
{
|
||||
electrostaticDepositionFvPatchScalarField& eVpf = eVPatch(patchi);
|
||||
|
||||
if (master == -1)
|
||||
{
|
||||
master = patchi;
|
||||
}
|
||||
|
||||
eVpf.master() = master;
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
void Foam::electrostaticDepositionFvPatchScalarField::round
|
||||
(
|
||||
scalarField& fld,
|
||||
const scalar dcml
|
||||
) const
|
||||
{
|
||||
for (auto& f : fld)
|
||||
{
|
||||
f = std::round(f*dcml)/dcml;
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
void Foam::electrostaticDepositionFvPatchScalarField::writeFilmFields() const
|
||||
{
|
||||
const auto& eV =
|
||||
db().lookupObject<volScalarField>(this->internalField().name());
|
||||
|
||||
const volScalarField::Boundary& bf = eV.boundaryField();
|
||||
|
||||
const fvMesh& mesh = eV.mesh();
|
||||
|
||||
volScalarField h
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
IOobject::scopedName(word("electrostaticDeposition"), "h"),
|
||||
mesh.time().timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::NO_WRITE,
|
||||
false // do not register
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar(dimLength)
|
||||
);
|
||||
|
||||
forAll(bf, patchi)
|
||||
{
|
||||
if (isA<electrostaticDepositionFvPatchScalarField>(bf[patchi]))
|
||||
{
|
||||
electrostaticDepositionFvPatchScalarField& eVpf = eVPatch(patchi);
|
||||
|
||||
auto& hp = h.boundaryFieldRef()[patchi];
|
||||
|
||||
hp = eVpf.h();
|
||||
}
|
||||
}
|
||||
|
||||
h.write();
|
||||
}
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
|
||||
|
||||
Foam::electrostaticDepositionFvPatchScalarField::
|
||||
electrostaticDepositionFvPatchScalarField
|
||||
(
|
||||
const fvPatch& p,
|
||||
const DimensionedField<scalar, volMesh>& iF
|
||||
)
|
||||
:
|
||||
fixedValueFvPatchScalarField(p, iF),
|
||||
h_(p.size(), 0),
|
||||
qcum_(p.size(), 0),
|
||||
Vfilm_(p.size(), 0),
|
||||
Ceffptr_(nullptr),
|
||||
rptr_(nullptr),
|
||||
jMin_(0),
|
||||
qMin_(0),
|
||||
Rbody_(0),
|
||||
Vi_(0),
|
||||
Vanode_(GREAT),
|
||||
phasesDict_(),
|
||||
phaseNames_(),
|
||||
phases_(),
|
||||
sigmas_(),
|
||||
sigma_(sqr(dimCurrent)*pow3(dimTime)/(dimMass*pow3(dimLength)), scalar(1)),
|
||||
timei_(-1),
|
||||
master_(-1)
|
||||
{}
|
||||
|
||||
|
||||
Foam::electrostaticDepositionFvPatchScalarField::
|
||||
electrostaticDepositionFvPatchScalarField
|
||||
(
|
||||
const fvPatch& p,
|
||||
const DimensionedField<scalar, volMesh>& iF,
|
||||
const dictionary& dict
|
||||
)
|
||||
:
|
||||
fixedValueFvPatchScalarField(p, iF, dict, false),
|
||||
h_("h", dict, p.size()),
|
||||
qcum_
|
||||
(
|
||||
dict.found("qCumulative")
|
||||
? scalarField("qCumulative", dict, p.size())
|
||||
: scalarField(p.size(), 0)
|
||||
),
|
||||
Vfilm_
|
||||
(
|
||||
dict.found("Vfilm")
|
||||
? scalarField("Vfilm", dict, p.size())
|
||||
: scalarField(p.size(), 0)
|
||||
),
|
||||
Ceffptr_
|
||||
(
|
||||
PatchFunction1<scalar>::New(p.patch(), "CoulombicEfficiency", dict)
|
||||
),
|
||||
rptr_(PatchFunction1<scalar>::New(p.patch(), "resistivity", dict)),
|
||||
jMin_(dict.getCheckOrDefault<scalar>("jMin", 0, scalarMinMax::ge(0))),
|
||||
qMin_(dict.getCheckOrDefault<scalar>("qMin", 0, scalarMinMax::ge(0))),
|
||||
Rbody_(dict.getCheckOrDefault<scalar>("Rbody", 0, scalarMinMax::ge(0))),
|
||||
Vi_(dict.getOrDefault<scalar>("Vi", 0)),
|
||||
Vanode_(dict.getOrDefault<scalar>("Vanode", GREAT)),
|
||||
phasesDict_(dict.subOrEmptyDict("phases")),
|
||||
phaseNames_(),
|
||||
phases_(),
|
||||
sigmas_(),
|
||||
sigma_
|
||||
(
|
||||
dimensionedScalar
|
||||
(
|
||||
sqr(dimCurrent)*pow3(dimTime)/(dimMass*pow3(dimLength)),
|
||||
dict.getCheckOrDefault<scalar>
|
||||
(
|
||||
"sigma",
|
||||
scalar(1),
|
||||
scalarMinMax::ge(SMALL)
|
||||
)
|
||||
)
|
||||
),
|
||||
timei_(-1),
|
||||
master_(-1)
|
||||
{
|
||||
if (dict.found("value"))
|
||||
{
|
||||
fvPatchScalarField::operator=
|
||||
(
|
||||
scalarField("value", dict, p.size())
|
||||
);
|
||||
}
|
||||
else
|
||||
{
|
||||
fvPatchScalarField::operator=(patchInternalField());
|
||||
}
|
||||
|
||||
// If flow is multiphase
|
||||
if (!phasesDict_.empty())
|
||||
{
|
||||
phaseNames_.setSize(phasesDict_.size());
|
||||
phases_.setSize(phasesDict_.size());
|
||||
sigmas_.setSize(phasesDict_.size());
|
||||
|
||||
label phasei = 0;
|
||||
forAllConstIters(phasesDict_, iter)
|
||||
{
|
||||
const word& key = iter().keyword();
|
||||
|
||||
if (!phasesDict_.isDict(key))
|
||||
{
|
||||
FatalErrorInFunction
|
||||
<< "Found non-dictionary entry " << iter()
|
||||
<< " in top-level dictionary " << phasesDict_
|
||||
<< exit(FatalError);
|
||||
}
|
||||
|
||||
const dictionary& subDict = phasesDict_.subDict(key);
|
||||
|
||||
phaseNames_[phasei] = key;
|
||||
|
||||
sigmas_.set
|
||||
(
|
||||
phasei,
|
||||
new dimensionedScalar
|
||||
(
|
||||
sqr(dimCurrent)*pow3(dimTime)/(dimMass*pow3(dimLength)),
|
||||
subDict.getCheck<scalar>
|
||||
(
|
||||
"sigma",
|
||||
scalarMinMax::ge(SMALL)
|
||||
)
|
||||
)
|
||||
);
|
||||
|
||||
++phasei;
|
||||
}
|
||||
|
||||
forAll(phaseNames_, i)
|
||||
{
|
||||
phases_.set
|
||||
(
|
||||
i,
|
||||
db().getObjectPtr<volScalarField>(phaseNames_[i])
|
||||
);
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
Foam::electrostaticDepositionFvPatchScalarField::
|
||||
electrostaticDepositionFvPatchScalarField
|
||||
(
|
||||
const electrostaticDepositionFvPatchScalarField& ptf,
|
||||
const fvPatch& p,
|
||||
const DimensionedField<scalar, volMesh>& iF,
|
||||
const fvPatchFieldMapper& mapper
|
||||
)
|
||||
:
|
||||
fixedValueFvPatchScalarField(ptf, p, iF, mapper),
|
||||
h_(ptf.h_, mapper),
|
||||
qcum_(ptf.qcum_, mapper),
|
||||
Vfilm_(ptf.Vfilm_, mapper),
|
||||
Ceffptr_(ptf.Ceffptr_.clone(p.patch())),
|
||||
rptr_(ptf.rptr_.clone(p.patch())),
|
||||
jMin_(ptf.jMin_),
|
||||
qMin_(ptf.qMin_),
|
||||
Rbody_(ptf.Rbody_),
|
||||
Vi_(ptf.Vi_),
|
||||
Vanode_(ptf.Vanode_),
|
||||
phasesDict_(ptf.phasesDict_),
|
||||
phaseNames_(ptf.phaseNames_),
|
||||
phases_(ptf.phases_),
|
||||
sigmas_(),
|
||||
sigma_(ptf.sigma_),
|
||||
timei_(ptf.timei_),
|
||||
master_(-1)
|
||||
{}
|
||||
|
||||
|
||||
Foam::electrostaticDepositionFvPatchScalarField::
|
||||
electrostaticDepositionFvPatchScalarField
|
||||
(
|
||||
const electrostaticDepositionFvPatchScalarField& ptf
|
||||
)
|
||||
:
|
||||
fixedValueFvPatchScalarField(ptf),
|
||||
h_(ptf.h_),
|
||||
qcum_(ptf.qcum_),
|
||||
Vfilm_(ptf.Vfilm_),
|
||||
Ceffptr_(ptf.Ceffptr_.clone(patch().patch())),
|
||||
rptr_(ptf.rptr_.clone(patch().patch())),
|
||||
jMin_(ptf.jMin_),
|
||||
qMin_(ptf.qMin_),
|
||||
Rbody_(ptf.Rbody_),
|
||||
Vi_(ptf.Vi_),
|
||||
Vanode_(ptf.Vanode_),
|
||||
phasesDict_(ptf.phasesDict_),
|
||||
phaseNames_(ptf.phaseNames_),
|
||||
phases_(ptf.phases_),
|
||||
sigmas_(),
|
||||
sigma_(ptf.sigma_),
|
||||
timei_(ptf.timei_),
|
||||
master_(-1)
|
||||
{}
|
||||
|
||||
|
||||
Foam::electrostaticDepositionFvPatchScalarField::
|
||||
electrostaticDepositionFvPatchScalarField
|
||||
(
|
||||
const electrostaticDepositionFvPatchScalarField& ptf,
|
||||
const DimensionedField<scalar, volMesh>& iF
|
||||
)
|
||||
:
|
||||
fixedValueFvPatchScalarField(ptf, iF),
|
||||
h_(ptf.h_),
|
||||
qcum_(ptf.qcum_),
|
||||
Vfilm_(ptf.Vfilm_),
|
||||
Ceffptr_(ptf.Ceffptr_.clone(patch().patch())),
|
||||
rptr_(ptf.rptr_.clone(patch().patch())),
|
||||
jMin_(ptf.jMin_),
|
||||
qMin_(ptf.qMin_),
|
||||
Rbody_(ptf.Rbody_),
|
||||
Vi_(ptf.Vi_),
|
||||
Vanode_(ptf.Vanode_),
|
||||
phasesDict_(ptf.phasesDict_),
|
||||
phaseNames_(ptf.phaseNames_),
|
||||
phases_(ptf.phases_),
|
||||
sigmas_(),
|
||||
sigma_(ptf.sigma_),
|
||||
timei_(ptf.timei_),
|
||||
master_(-1)
|
||||
{}
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
|
||||
|
||||
void Foam::electrostaticDepositionFvPatchScalarField::autoMap
|
||||
(
|
||||
const fvPatchFieldMapper& m
|
||||
)
|
||||
{
|
||||
fixedValueFvPatchScalarField::autoMap(m);
|
||||
|
||||
h_.autoMap(m);
|
||||
qcum_.autoMap(m);
|
||||
Vfilm_.autoMap(m);
|
||||
|
||||
if (Ceffptr_)
|
||||
{
|
||||
Ceffptr_->autoMap(m);
|
||||
}
|
||||
|
||||
if (rptr_)
|
||||
{
|
||||
rptr_->autoMap(m);
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
void Foam::electrostaticDepositionFvPatchScalarField::rmap
|
||||
(
|
||||
const fvPatchScalarField& ptf,
|
||||
const labelList& addr
|
||||
)
|
||||
{
|
||||
fixedValueFvPatchScalarField::rmap(ptf, addr);
|
||||
|
||||
const auto& tiptf =
|
||||
refCast<const electrostaticDepositionFvPatchScalarField>(ptf);
|
||||
|
||||
h_.rmap(tiptf.h_, addr);
|
||||
qcum_.rmap(tiptf.qcum_, addr);
|
||||
Vfilm_.rmap(tiptf.Vfilm_, addr);
|
||||
|
||||
if (Ceffptr_)
|
||||
{
|
||||
Ceffptr_->rmap(tiptf.Ceffptr_(), addr);
|
||||
}
|
||||
|
||||
if (rptr_)
|
||||
{
|
||||
rptr_->rmap(tiptf.rptr_(), addr);
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
Foam::tmp<Foam::scalarField>
|
||||
Foam::electrostaticDepositionFvPatchScalarField::sigma() const
|
||||
{
|
||||
const label patchi = patch().index();
|
||||
|
||||
if (phases_.size())
|
||||
{
|
||||
tmp<scalarField> tsigma =
|
||||
phases_[0].boundaryField()[patchi]*sigmas_[0].value();
|
||||
|
||||
for (label i = 1; i < phases_.size(); ++i)
|
||||
{
|
||||
tsigma.ref() +=
|
||||
phases_[i].boundaryField()[patchi]*sigmas_[i].value();
|
||||
}
|
||||
|
||||
return tsigma;
|
||||
}
|
||||
|
||||
return tmp<scalarField>::New(patch().size(), sigma_.value());
|
||||
}
|
||||
|
||||
|
||||
void Foam::electrostaticDepositionFvPatchScalarField::updateCoeffs()
|
||||
{
|
||||
if (updated())
|
||||
{
|
||||
return;
|
||||
}
|
||||
|
||||
if (timei_ == db().time().timeIndex())
|
||||
{
|
||||
return;
|
||||
}
|
||||
|
||||
const scalar t = db().time().timeOutputValue();
|
||||
const scalar dt = db().time().deltaTValue();
|
||||
const label patchi = patch().index();
|
||||
|
||||
const auto& eV =
|
||||
db().lookupObject<volScalarField>(this->internalField().name());
|
||||
|
||||
// Current density on film interface
|
||||
tmp<scalarField> tjnp = -this->sigma()*eV.boundaryField()[patchi].snGrad();
|
||||
scalarField& jnp = tjnp.ref();
|
||||
jnp = max(jnp, scalar(0)); // experimental - do not allow any negative jnp
|
||||
// experimental - avoid micro/nano currents/volts
|
||||
// to reduce snowballing effects of lateral gradients on the patch
|
||||
round(jnp);
|
||||
|
||||
|
||||
// Calculate film-thickness finite increments
|
||||
tmp<scalarField> tCoulombicEfficiency = Ceffptr_->value(t);
|
||||
tmp<scalarField> tdh = tCoulombicEfficiency*(jnp - jMin_)*dt;
|
||||
scalarField& dh = tdh.ref();
|
||||
|
||||
// Do not allow any depletion or abrasion of deposition
|
||||
dh = max(dh, scalar(0));
|
||||
|
||||
// Do not allow any deposition when accumulative specific
|
||||
// charge is less than minimum accumulative specific charge
|
||||
qcum_ += jnp*dt;
|
||||
|
||||
forAll(dh, i)
|
||||
{
|
||||
if (qcum_[i] < qMin_)
|
||||
{
|
||||
dh[i] = 0;
|
||||
}
|
||||
}
|
||||
|
||||
// Add finite increments of film thickness to total film thickness
|
||||
h_ += dh;
|
||||
|
||||
|
||||
// Calculate incremental electric potential due to film resistance
|
||||
tmp<scalarField> tresistivity = rptr_->value(t);
|
||||
tmp<scalarField> tRfilm = tresistivity*tdh;
|
||||
tmp<scalarField> tdV = jnp*tRfilm;
|
||||
Vfilm_ += tdV;
|
||||
Vfilm_ = min(Vfilm_, Vanode_);
|
||||
|
||||
|
||||
// Calculate electric potential due to body resistance
|
||||
tmp<scalarField> tVbody = tjnp*Rbody_;
|
||||
|
||||
|
||||
// Add all electric potential contributions
|
||||
operator==(min(Vi_ + Vfilm_ + tVbody, Vanode_));
|
||||
|
||||
|
||||
fixedValueFvPatchScalarField::updateCoeffs();
|
||||
|
||||
timei_ = db().time().timeIndex();
|
||||
|
||||
{
|
||||
const scalar hMin = gMin(h_);
|
||||
const scalar hMax = gMax(h_);
|
||||
const scalar hAvg = gAverage(h_);
|
||||
|
||||
if (Pstream::master())
|
||||
{
|
||||
Info<< " patch: " << patch().name()
|
||||
<< ", h: min = " << hMin
|
||||
<< ", max = " << hMax
|
||||
<< ", average = " << hAvg << nl
|
||||
<< endl;
|
||||
}
|
||||
}
|
||||
|
||||
// Write here to avoid any upset to redistributePar-decompose
|
||||
if (db().time().writeTime())
|
||||
{
|
||||
// Write film thickness fields as patch fields of a volScalarField
|
||||
setMaster();
|
||||
|
||||
if (patch().index() == master_)
|
||||
{
|
||||
writeFilmFields();
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
void Foam::electrostaticDepositionFvPatchScalarField::write(Ostream& os) const
|
||||
{
|
||||
fvPatchScalarField::write(os);
|
||||
|
||||
h_.writeEntry("h", os);
|
||||
|
||||
if (Ceffptr_)
|
||||
{
|
||||
Ceffptr_->writeData(os);
|
||||
}
|
||||
|
||||
if (rptr_)
|
||||
{
|
||||
rptr_->writeData(os);
|
||||
}
|
||||
|
||||
if (!phasesDict_.empty())
|
||||
{
|
||||
phasesDict_.writeEntry(phasesDict_.dictName(), os);
|
||||
}
|
||||
else
|
||||
{
|
||||
sigma_.writeEntry("sigma", os);
|
||||
}
|
||||
|
||||
os.writeEntryIfDifferent<scalar>("jMin", 0, jMin_);
|
||||
os.writeEntryIfDifferent<scalar>("qMin", 0, qMin_);
|
||||
os.writeEntryIfDifferent<scalar>("Rbody", 0, Rbody_);
|
||||
os.writeEntryIfDifferent<scalar>("Vi", 0, Vi_);
|
||||
os.writeEntryIfDifferent<scalar>("Vanode", GREAT, Vanode_);
|
||||
qcum_.writeEntry("qCumulative", os);
|
||||
Vfilm_.writeEntry("Vfilm", os);
|
||||
|
||||
writeEntry("value", os);
|
||||
}
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
namespace Foam
|
||||
{
|
||||
makePatchTypeField
|
||||
(
|
||||
fvPatchScalarField,
|
||||
electrostaticDepositionFvPatchScalarField
|
||||
);
|
||||
}
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -0,0 +1,385 @@
|
||||
/*---------------------------------------------------------------------------*\
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | www.openfoam.com
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
Copyright (C) 2021 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::electrostaticDepositionFvPatchScalarField
|
||||
|
||||
Group
|
||||
grpGenericBoundaryConditions
|
||||
|
||||
Description
|
||||
The \c electrostaticDeposition is a boundary condition to
|
||||
calculate electric potential (\c V) on a given boundary
|
||||
based on film thickness (\c h) and film resistance (\c R) fields
|
||||
which are updated based on a given patch-normal current density
|
||||
field (\c jn), Coulombic efficiency and film resistivity.
|
||||
|
||||
\f[
|
||||
j_n = - \sigma \nabla^\perp_p V = - \sigma (\vec{n}\cdot(\nabla V)_p)
|
||||
\f]
|
||||
|
||||
\f[
|
||||
\frac{dh}{dt} = C_{eff} (j_n - j_{min})
|
||||
\f]
|
||||
|
||||
\f[
|
||||
\frac{dR}{dt} = \rho \frac{dh}{dt} = \rho C_{eff} (j_n - j_{min})
|
||||
\f]
|
||||
|
||||
\f[
|
||||
V_{film}^n = V_{film}^o + j_n R_\Delta
|
||||
\f]
|
||||
|
||||
\f[
|
||||
V_{body} = j_n R_{body}
|
||||
\f]
|
||||
|
||||
\f[
|
||||
V_p^n = V_i + V_{body} + V_{film}^n
|
||||
\f]
|
||||
|
||||
where
|
||||
\vartable
|
||||
j_n | Patch-normal current density [A/m^2]
|
||||
V_p^n | Electric potential on film-fluid interface [volt = kg m^2/(A s^3)]
|
||||
V_p^o | Previous time-step electric potential on the interface [volt]
|
||||
V_{film} | Electric potential due to film resistance [volt]
|
||||
V_{body} | Electric potential due to body resistance [volt]
|
||||
V_i | Initial electric potential [volt]
|
||||
R_\Delta| Film resistance (finite increment) [ohm m^2 = kg m^4/(A^2 s^3)]
|
||||
R_{body} | Body resistance [ohm m^2 = kg m^4/(A^2 s^3)]
|
||||
\rho | Isotropic film resistivity [ohm m = kg m^3/(A^2 s^3)]
|
||||
h | Film thickness [m]
|
||||
C_{eff} | Volumetric Coulombic efficiency [m^3/(A s)]
|
||||
j_{min} | Minimum current density for deposition onset [A/m^2]
|
||||
\sigma | Isotropic conductivity of mixture [S/m = A^2 s^3/(kg m^3)]
|
||||
\vec{n} | Patch-normal unit vector [-]
|
||||
\endvartable
|
||||
|
||||
Usage
|
||||
Example of the boundary condition specification:
|
||||
\verbatim
|
||||
<patchName>
|
||||
{
|
||||
// Mandatory entries
|
||||
type electrostaticDeposition;
|
||||
h <scalarField>;
|
||||
CoulombicEfficiency <PatchFunction1>;
|
||||
resistivity <PatchFunction1>;
|
||||
|
||||
// Conditional mandatory entries
|
||||
|
||||
// Option-1: single-phase
|
||||
sigma <scalar>;
|
||||
|
||||
// Option-2: multiphase
|
||||
phases
|
||||
{
|
||||
alpha.air
|
||||
{
|
||||
sigma <scalar>;
|
||||
}
|
||||
alpha.water
|
||||
{
|
||||
sigma <scalar>;
|
||||
}
|
||||
alpha.mercury
|
||||
{
|
||||
sigma <scalar>;
|
||||
}
|
||||
...
|
||||
}
|
||||
|
||||
// Optional entries
|
||||
jMin <scalar>;
|
||||
qMin <scalar>;
|
||||
Rbody <scalar>;
|
||||
Vi <scalar>;
|
||||
Vanode <scalar>;
|
||||
qCumulative <scalarField>;
|
||||
|
||||
// Inherited entries
|
||||
...
|
||||
}
|
||||
\endverbatim
|
||||
|
||||
where the entries mean:
|
||||
\table
|
||||
Property | Description | Type | Reqd | Deflt
|
||||
type | Type name: electrostaticDeposition | word | yes | -
|
||||
h | Film thickness | scalarField | yes | -
|
||||
CoulombicEfficiency | Coulombic efficiency <!--
|
||||
--> | PatchFunction1\<scalar\> | yes | -
|
||||
resistivity | Isotropic film resistivity <!--
|
||||
--> | PatchFunction1\<scalar\> | yes | -
|
||||
sigma | Isotropic electrical conductivity of phase | scalar | yes | -
|
||||
jMin | Minimum current density for deposition onset | scalar | no | 0
|
||||
qMin | Minimum accumulative specific charge for deposition onset <!--
|
||||
--> | scalar | no | 0
|
||||
Rbody | Resistance due to main body and/or <!--
|
||||
--> pretreatment layers | scalar | no | 0
|
||||
Vi | Initial electric potential | scalar | no | 0
|
||||
Vanode | Anode electric potential | scalar | no | GREAT
|
||||
qCumulative | Accumulative specific charge [A s/m^2] <!--
|
||||
--> | scalarField | no | 0
|
||||
\endtable
|
||||
|
||||
The inherited entries are elaborated in:
|
||||
- \link fixedValueFvPatchFields.H \endlink
|
||||
- \link PatchFunction1.H \endlink
|
||||
|
||||
Note
|
||||
- Depletion or abrasion of material due to negative current is not allowed.
|
||||
- When accumulative specific charge (\c qCumulative) is less than minimum
|
||||
accumulative specific charge (\c qMin), no deposition occurs.
|
||||
- Boundary-condition updates are not allowed during outer corrections
|
||||
to prevent spurious accumulation of film thickness.
|
||||
- \c resistivity, \c jMin, \c qMin and \c Rbody are always non-negative.
|
||||
|
||||
SourceFiles
|
||||
electrostaticDepositionFvPatchScalarField.C
|
||||
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
#ifndef electrostaticDepositionFvPatchScalarField_H
|
||||
#define electrostaticDepositionFvPatchScalarField_H
|
||||
|
||||
#include "fixedValueFvPatchFields.H"
|
||||
#include "Enum.H"
|
||||
#include "PatchFunction1.H"
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
namespace Foam
|
||||
{
|
||||
|
||||
/*---------------------------------------------------------------------------*\
|
||||
Class electrostaticDepositionFvPatchScalarField Declaration
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
class electrostaticDepositionFvPatchScalarField
|
||||
:
|
||||
public fixedValueFvPatchScalarField
|
||||
{
|
||||
// Private Data
|
||||
|
||||
//- Film thickness
|
||||
scalarField h_;
|
||||
|
||||
//- Accumulative specific charge
|
||||
scalarField qcum_;
|
||||
|
||||
//- Electric potential due to film resistance
|
||||
scalarField Vfilm_;
|
||||
|
||||
//- Coulombic efficiency
|
||||
autoPtr<PatchFunction1<scalar>> Ceffptr_;
|
||||
|
||||
//- Isotropic film resistivity
|
||||
autoPtr<PatchFunction1<scalar>> rptr_;
|
||||
|
||||
//- Minimum current density for the deposition onset
|
||||
scalar jMin_;
|
||||
|
||||
//- Minimum accumulative specific charge for the deposition onset
|
||||
scalar qMin_;
|
||||
|
||||
//- Resistance due to main body and/or pretreatment layers
|
||||
scalar Rbody_;
|
||||
|
||||
//- Initial electric potential
|
||||
scalar Vi_;
|
||||
|
||||
//- Anode electric potential
|
||||
scalar Vanode_;
|
||||
|
||||
// Isotropic electrical conductivitiy properties
|
||||
|
||||
//- Dictionary of phase data
|
||||
dictionary phasesDict_;
|
||||
|
||||
//- List of phase names
|
||||
wordList phaseNames_;
|
||||
|
||||
//- Unallocated list of phase fields
|
||||
UPtrList<volScalarField> phases_;
|
||||
|
||||
//- List of isotropic electrical conductivity of phases
|
||||
PtrList<dimensionedScalar> sigmas_;
|
||||
|
||||
//- Isotropic electrical conductivity of a single phase
|
||||
dimensionedScalar sigma_;
|
||||
|
||||
//- Time index - used to prevent film accumulation during outer iters
|
||||
label timei_;
|
||||
|
||||
//- Master patch ID
|
||||
mutable label master_;
|
||||
|
||||
|
||||
// Private Member Functions
|
||||
|
||||
//- Return non-const access to an electrostaticDeposition patch
|
||||
electrostaticDepositionFvPatchScalarField& eVPatch
|
||||
(
|
||||
const label patchi
|
||||
) const;
|
||||
|
||||
//- Return non-const access to the master patch ID
|
||||
label& master() noexcept
|
||||
{
|
||||
return master_;
|
||||
}
|
||||
|
||||
//- Set master patch ID
|
||||
void setMaster() const;
|
||||
|
||||
//- Round scalars of a given scalar field to dcml points decimal
|
||||
void round(scalarField& fld, const scalar dcml=1e8) const;
|
||||
|
||||
//- Write film thickness field to facilitate postprocessing
|
||||
void writeFilmFields() const;
|
||||
|
||||
|
||||
public:
|
||||
|
||||
//- Runtime type information
|
||||
TypeName("electrostaticDeposition");
|
||||
|
||||
|
||||
// Constructors
|
||||
|
||||
//- Construct from patch and internal field
|
||||
electrostaticDepositionFvPatchScalarField
|
||||
(
|
||||
const fvPatch&,
|
||||
const DimensionedField<scalar, volMesh>&
|
||||
);
|
||||
|
||||
//- Construct from patch, internal field and dictionary
|
||||
electrostaticDepositionFvPatchScalarField
|
||||
(
|
||||
const fvPatch&,
|
||||
const DimensionedField<scalar, volMesh>&,
|
||||
const dictionary&
|
||||
);
|
||||
|
||||
//- Construct by mapping given
|
||||
//- electrostaticDepositionFvPatchScalarField onto a new patch
|
||||
electrostaticDepositionFvPatchScalarField
|
||||
(
|
||||
const electrostaticDepositionFvPatchScalarField&,
|
||||
const fvPatch&,
|
||||
const DimensionedField<scalar, volMesh>&,
|
||||
const fvPatchFieldMapper&
|
||||
);
|
||||
|
||||
//- Construct as copy
|
||||
electrostaticDepositionFvPatchScalarField
|
||||
(
|
||||
const electrostaticDepositionFvPatchScalarField&
|
||||
);
|
||||
|
||||
//- Construct and return a clone
|
||||
virtual tmp<fvPatchScalarField> clone() const
|
||||
{
|
||||
return tmp<fvPatchScalarField>
|
||||
(
|
||||
new electrostaticDepositionFvPatchScalarField(*this)
|
||||
);
|
||||
}
|
||||
|
||||
//- Construct as copy setting internal field reference
|
||||
electrostaticDepositionFvPatchScalarField
|
||||
(
|
||||
const electrostaticDepositionFvPatchScalarField&,
|
||||
const DimensionedField<scalar, volMesh>&
|
||||
);
|
||||
|
||||
//- Construct and return a clone setting internal field reference
|
||||
virtual tmp<fvPatchScalarField> clone
|
||||
(
|
||||
const DimensionedField<scalar, volMesh>& iF
|
||||
) const
|
||||
{
|
||||
return tmp<fvPatchScalarField>
|
||||
(
|
||||
new electrostaticDepositionFvPatchScalarField(*this, iF)
|
||||
);
|
||||
}
|
||||
|
||||
|
||||
// Member Functions
|
||||
|
||||
// Access
|
||||
|
||||
//- Return const access to film thickness patch field
|
||||
const scalarField& h() const noexcept
|
||||
{
|
||||
return h_;
|
||||
}
|
||||
|
||||
|
||||
// Mapping
|
||||
|
||||
//- Map (and resize as needed) from self given a mapping object
|
||||
virtual void autoMap
|
||||
(
|
||||
const fvPatchFieldMapper&
|
||||
);
|
||||
|
||||
//- Reverse map the given fvPatchField onto this fvPatchField
|
||||
virtual void rmap
|
||||
(
|
||||
const fvPatchScalarField&,
|
||||
const labelList&
|
||||
);
|
||||
|
||||
|
||||
// Evaluation
|
||||
|
||||
//- Return the isotropic electrical conductivity field of mixture
|
||||
tmp<scalarField> sigma() const;
|
||||
|
||||
//- Update the coefficients associated with the patch field
|
||||
virtual void updateCoeffs();
|
||||
|
||||
|
||||
// I-O
|
||||
|
||||
//- Write
|
||||
virtual void write(Ostream&) const;
|
||||
};
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
} // End namespace Foam
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
#endif
|
||||
|
||||
// ************************************************************************* //
|
||||
Reference in New Issue
Block a user