[mehtab]

Hi

I would like to implement some boundary conditions using groovyBC and need some help in this regard.

1> Temperature boundary at Inlet

T_s = T_0 + 1/C_p(Q_f/m" - Hvap)


Here Q_f and m ̇ " needs to be computed for inlet patch at previous time step. Also C_p needs to be updated from thermo



Is it possible to do it using groovyBC? Which option is better groovyBC or codedFixedValue in OpenFOAM.

Thanks in advance.....

[mehtab]

Hi,

I tried implementing temperature boundary condition at Inlet as:

T_s = T_0 + 1/C_p(Q_f/m" - Hv)


Here is my files

boundaryField
{
outlet
{
type inletOutlet;
inletValue uniform 298.15;
value uniform 298.15;
}
sides
{
type inletOutlet;
inletValue uniform 298.15;
value uniform 298.15;
}
base
{
type zeroGradient;
}
inlet
{
type groovyBC;
valueExpression "T0+((1/cpi)*((Qf/m)-Hv))";
variables "T0=300;m{inlet}=phi;Hv=8500;Qf{inlet}=Qr;cpi{inle t}=cp;";
value 300;

}

controlDict
application fireFoam;

startFrom startTime;

startTime 0.0;

stopAt endTime;

endTime 1.0;

deltaT 1.0e-2;

writeControl adjustableRunTime;

writeInterval 0.01;

purgeWrite 0;

writeFormat ascii;

writePrecision 6;

writeCompression uncompressed;

timeFormat general;

timePrecision 6;

graphFormat raw;

runTimeModifiable yes;

adjustTimeStep yes;

maxCo 0.6;

maxDeltaT 0.1;

functions
{
fieldAverage
{
type fieldAverage;
functionObjectLibs ( "libfieldFunctionObjects.so" );
enabled true;
outputControl outputTime;

fields
(
U
{
mean on;
prime2Mean on;
base time;
}
T
{
mean on;
prime2Mean on;
base time;
}
CH4
{
mean on;
prime2Mean on;
base time;
}
CO2
{
mean on;
prime2Mean on;
base time;
}
H2O
{
mean on;
prime2Mean on;
base time;
}
O2
{
mean on;
prime2Mean on;
base time;
}
);
}

libs (
"libOpenFOAM.so"
"libgroovyBC.so"
"libgroovyStandardBCs.so"
"libsimpleFunctionObjects.so"
"libsimpleSwakFunctionObjects.so"
"libswakFunctionObjects.so"
"libswakPythonIntegration.so"
"libswakTransportTurbFunctionPlugin.so"
) ;

I also checked my solver fireFoam createFields.H to check how Cp is being used.

createFields.H

volScalarField cp
(
IOobject
(
"cp",
runTime.timeName(),
mesh,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
thermo.Cp()
);

I am getting this error:

swak4Foam: Setting default mesh
swak4Foam: Allocating new repository for sampledGlobalVariables
--> FOAM Warning :
From function ConcretePluginFunction<DriverType>::exists
in file lnInclude/ConcretePluginFunction.C at line 121
Constructor table of plugin functions for PatchValueExpressionDriver is not initialized


--> FOAM FATAL ERROR:
Parser Error for driver PatchValueExpressionDriver at "1.1-2" :"field Qr not existing or of wrong type"
"Qr"
^^
--|

Context of the error:


- Driver constructed from scratch
Evaluating expression "Qr"


From function parsingValue
in file lnInclude/CommonValueExpressionDriverI.H at line 1189.

FOAM exiting


What should I modify in me expressions? Have I defined cp correctly or should I change to thermo.Cp? Here Qr is totat radiative flux being calculated through radiation model. I am not able to understand the error. Please help

Thanks in advance....