Skip to content
Snippets Groups Projects
Commit 16cf8eaf authored by mattijs's avatar mattijs
Browse files

ENH: Changed flux condition into gradient bc

parent fd952687
Branches
Tags
No related merge requests found
...@@ -37,7 +37,7 @@ solidWallHeatFluxTemperatureFvPatchScalarField ...@@ -37,7 +37,7 @@ solidWallHeatFluxTemperatureFvPatchScalarField
const DimensionedField<scalar, volMesh>& iF const DimensionedField<scalar, volMesh>& iF
) )
: :
fixedValueFvPatchScalarField(p, iF), fixedGradientFvPatchScalarField(p, iF),
q_(p.size(), 0.0), q_(p.size(), 0.0),
KName_("undefined-K") KName_("undefined-K")
{} {}
...@@ -52,7 +52,7 @@ solidWallHeatFluxTemperatureFvPatchScalarField ...@@ -52,7 +52,7 @@ solidWallHeatFluxTemperatureFvPatchScalarField
const fvPatchFieldMapper& mapper const fvPatchFieldMapper& mapper
) )
: :
fixedValueFvPatchScalarField(ptf, p, iF, mapper), fixedGradientFvPatchScalarField(ptf, p, iF, mapper),
q_(ptf.q_, mapper), q_(ptf.q_, mapper),
KName_(ptf.KName_) KName_(ptf.KName_)
{} {}
...@@ -66,7 +66,7 @@ solidWallHeatFluxTemperatureFvPatchScalarField ...@@ -66,7 +66,7 @@ solidWallHeatFluxTemperatureFvPatchScalarField
const dictionary& dict const dictionary& dict
) )
: :
fixedValueFvPatchScalarField(p, iF, dict), fixedGradientFvPatchScalarField(p, iF, dict),
q_("q", dict, p.size()), q_("q", dict, p.size()),
KName_(dict.lookup("K")) KName_(dict.lookup("K"))
{} {}
...@@ -78,7 +78,7 @@ solidWallHeatFluxTemperatureFvPatchScalarField ...@@ -78,7 +78,7 @@ solidWallHeatFluxTemperatureFvPatchScalarField
const solidWallHeatFluxTemperatureFvPatchScalarField& tppsf const solidWallHeatFluxTemperatureFvPatchScalarField& tppsf
) )
: :
fixedValueFvPatchScalarField(tppsf), fixedGradientFvPatchScalarField(tppsf),
q_(tppsf.q_), q_(tppsf.q_),
KName_(tppsf.KName_) KName_(tppsf.KName_)
{} {}
...@@ -91,7 +91,7 @@ solidWallHeatFluxTemperatureFvPatchScalarField ...@@ -91,7 +91,7 @@ solidWallHeatFluxTemperatureFvPatchScalarField
const DimensionedField<scalar, volMesh>& iF const DimensionedField<scalar, volMesh>& iF
) )
: :
fixedValueFvPatchScalarField(tppsf, iF), fixedGradientFvPatchScalarField(tppsf, iF),
q_(tppsf.q_), q_(tppsf.q_),
KName_(tppsf.KName_) KName_(tppsf.KName_)
{} {}
...@@ -104,7 +104,7 @@ void Foam::solidWallHeatFluxTemperatureFvPatchScalarField::autoMap ...@@ -104,7 +104,7 @@ void Foam::solidWallHeatFluxTemperatureFvPatchScalarField::autoMap
const fvPatchFieldMapper& m const fvPatchFieldMapper& m
) )
{ {
fixedValueFvPatchScalarField::autoMap(m); fixedGradientFvPatchScalarField::autoMap(m);
q_.autoMap(m); q_.autoMap(m);
} }
...@@ -115,7 +115,7 @@ void Foam::solidWallHeatFluxTemperatureFvPatchScalarField::rmap ...@@ -115,7 +115,7 @@ void Foam::solidWallHeatFluxTemperatureFvPatchScalarField::rmap
const labelList& addr const labelList& addr
) )
{ {
fixedValueFvPatchScalarField::rmap(ptf, addr); fixedGradientFvPatchScalarField::rmap(ptf, addr);
const solidWallHeatFluxTemperatureFvPatchScalarField& hfptf = const solidWallHeatFluxTemperatureFvPatchScalarField& hfptf =
refCast<const solidWallHeatFluxTemperatureFvPatchScalarField>(ptf); refCast<const solidWallHeatFluxTemperatureFvPatchScalarField>(ptf);
...@@ -131,14 +131,14 @@ void Foam::solidWallHeatFluxTemperatureFvPatchScalarField::updateCoeffs() ...@@ -131,14 +131,14 @@ void Foam::solidWallHeatFluxTemperatureFvPatchScalarField::updateCoeffs()
return; return;
} }
const scalarField& Kw = const scalarField& Kw = patch().lookupPatchField<volScalarField, scalar>
patch().lookupPatchField<volScalarField, scalar>(KName_); (
KName_
const fvPatchScalarField& Tw = *this; );
operator==(q_/(patch().deltaCoeffs()*Kw) + Tw.patchInternalField()); gradient() = q_/Kw;
fixedValueFvPatchScalarField::updateCoeffs(); fixedGradientFvPatchScalarField::updateCoeffs();
} }
...@@ -147,9 +147,10 @@ void Foam::solidWallHeatFluxTemperatureFvPatchScalarField::write ...@@ -147,9 +147,10 @@ void Foam::solidWallHeatFluxTemperatureFvPatchScalarField::write
Ostream& os Ostream& os
) const ) const
{ {
fixedValueFvPatchScalarField::write(os); fixedGradientFvPatchScalarField::write(os);
q_.writeEntry("q", os); q_.writeEntry("q", os);
os.writeKeyword("K") << KName_ << token::END_STATEMENT << nl; os.writeKeyword("K") << KName_ << token::END_STATEMENT << nl;
this->writeEntry("value", os);
} }
......
...@@ -45,7 +45,7 @@ SourceFiles ...@@ -45,7 +45,7 @@ SourceFiles
#ifndef solidWallHeatFluxTemperatureFvPatchScalarField_H #ifndef solidWallHeatFluxTemperatureFvPatchScalarField_H
#define solidWallHeatFluxTemperatureFvPatchScalarField_H #define solidWallHeatFluxTemperatureFvPatchScalarField_H
#include "fixedValueFvPatchFields.H" #include "fixedGradientFvPatchFields.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
...@@ -58,7 +58,7 @@ namespace Foam ...@@ -58,7 +58,7 @@ namespace Foam
class solidWallHeatFluxTemperatureFvPatchScalarField class solidWallHeatFluxTemperatureFvPatchScalarField
: :
public fixedValueFvPatchScalarField public fixedGradientFvPatchScalarField
{ {
// Private data // Private data
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment