Thank you so much for your reply.
I have one more question:
I try to explain my problem.
I am using chtMultiRegionSimpleFoam. My Openfoam version is 4.1.
I made my geometry in salome, and then imported to OpenFoam as a UNV file.
I have multiple regions in my geometry, (heaters, box, and fluid).
Both heaters and fluid are in the box.
Now, I want to simulate the turbulent flow in the fluid region. I am using kEpsilon Model.
I am trying to model incompressible flow.
I put boundary conditions for alphat, k, epsilon, nut etc in system/fluid/changeDictionaryDict file. This file is shown below:
Code:
boundary
{
inlet
{
type patch;
}
outlet
{
type patch;
}
}
T
{
internalField uniform 300;
boundaryField
{
inlet
{
type fixedValue;
value $internalField;
}
outlet
{
type zeroGradient;
value $internalField;
}
"fluid_to_box"
{
type zeroGradient;
}
}
}
U
{
internalField uniform (0 0 0);
boundaryField
{
inlet
{
type fixedValue;
value uniform (0 0 -1.33);
}
outlet
{
type zeroGradient;
}
"fluid_to_box"
{
type noSlip;
}
}
}
epsilon
{
internalField uniform 0.01;
boundaryField
{
inlet
{
type fixedValue;
value uniform 0.01;
}
outlet
{
type zeroGradient;
//value uniform 0;
}
".*"
{
type epsilonWallFunction;
value uniform 0.01;
}
}
}
k
{
internalField uniform 0.1;
boundaryField
{
inlet
{
type fixedValue;
value uniform 0.1;
}
outlet
{
type zeroGradient;
//value uniform 0;
}
"fluid_to_box"
{
type kqRWallFunction;
value uniform 0.1;
}
}
}
p_rgh
{
internalField uniform 0;
boundaryField
{
inlet
{
type zeroGradient;
value uniform 0;
}
outlet
{
type fixedValue;
value uniform 0;
}
".*"
{
type fixedFluxPressure;
value uniform 0;
}
}
}
p
{
internalField uniform 0;
boundaryField
{
".*"
{
type calculated;
value uniform 0;
}
}
}
alphat
{
internalField uniform 0;
boundaryField
{
inlet
{
type calculated;//fixedValue;
value uniform 0;
}
outlet
{
type calculated;
value uniform 0;
}
"fluid_to_box"
{
type alphatJayatillekeWallFunction;
value uniform 0;
}
}
}
nut
{
internalField uniform 0;
boundaryField
{
inlet
{
type calculated;
value uniform 0;
}
outlet
{
type calculated;
value uniform 0;
}
"fluid_to_box"
{
type nutkWallFunction;
value uniform 0;
}
}
}
// ************************************************************************* //
The error I posted before, I started getting when I shifted from laminar to Turbulent.
Now, I have the following doubts:
1. How and where actually we define that out flow would be incompressible?
2. As you can see in the file above, in the
alphat portion, for the boundary
fluid_to_box, I have put the wall function. before it was compressible written before the wall function, but I removed it because I am trying to simulate incompressible. But it gave error that:
Code:
--> FOAM FATAL IO ERROR:
Unknown patchField type alphatJayatillekeWallFunction for patch type mappedWall
Valid patchField types are :
112
(
MarshakRadiation
MarshakRadiationFixedTemperature
advective
atmBoundaryLayerInletEpsilon
atmBoundaryLayerInletK
calculated
codedFixedValue
codedMixed
compressible::alphatJayatillekeWallFunction
compressible::alphatWallFunction
compressible::thermalBaffle1D<hConstSolidThermoPhysics>
compressible::thermalBaffle1D<hPowerSolidThermoPhysics>
compressible::turbulentHeatFluxTemperature
compressible::turbulentTemperatureCoupledBaffleMixed
compressible::turbulentTemperatureRadCoupledMixed
convectiveHeatTransfer
cyclic
cyclicACMI
cyclicAMI
cyclicSlip
directionMixed
empty
energyJump
energyJumpAMI
epsilonLowReWallFunction
epsilonWallFunction
externalCoupled
externalCoupledTemperature
externalWallHeatFluxTemperature
extrapolatedCalculated
fWallFunction
fan
fanPressure
fixedEnergy
fixedFluxExtrapolatedPressure
fixedFluxPressure
fixedGradient
fixedInternalValue
fixedJump
fixedJumpAMI
fixedMean
fixedPressureCompressibleDensity
fixedProfile
fixedUnburntEnthalpy
fixedValue
freestream
freestreamPressure
gradientEnergy
gradientUnburntEnthalpy
greyDiffusiveRadiation
greyDiffusiveRadiationViewFactor
inletOutlet
inletOutletTotalTemperature
kLowReWallFunction
kqRWallFunction
mapped
mappedField
mappedFixedInternalValue
mappedFixedPushedInternalValue
mixed
mixedEnergy
mixedUnburntEnthalpy
nonuniformTransformCyclic
nutLowReWallFunction
nutTabulatedWallFunction
nutURoughWallFunction
nutUSpaldingWallFunction
nutUWallFunction
nutkAtmRoughWallFunction
nutkRoughWallFunction
nutkWallFunction
omegaWallFunction
outletInlet
outletMappedUniformInlet
partialSlip
phaseHydrostaticPressure
plenumPressure
porousBafflePressure
prghPressure
prghTotalHydrostaticPressure
prghTotalPressure
processor
processorCyclic
rotatingTotalPressure
sliced
slip
symmetry
symmetryPlane
syringePressure
timeVaryingMappedFixedValue
totalFlowRateAdvectiveDiffusive
totalPressure
totalTemperature
turbulentInlet
turbulentIntensityKineticEnergyInlet
turbulentMixingLengthDissipationRateInlet
turbulentMixingLengthFrequencyInlet
uniformDensityHydrostaticPressure
uniformFixedGradient
uniformFixedValue
uniformInletOutlet
uniformJump
uniformJumpAMI
uniformTotalPressure
v2WallFunction
variableHeightFlowRate
wallHeatTransfer
waveSurfacePressure
waveTransmissive
wedge
wideBandDiffusiveRadiation
zeroGradient
)
3. In case of incompressible, which wall function shall I put here on the boundary
fluid_to_box(mappedwall)?
4. Can I use other boundary conditions, other than wallfunctions on this boundary
fluid_to_box(mappedwall)?
I tried to explain my problem, but if you need any other information, I would be happy to give.
Thank you