[ANacc]

I am interested in applying a Johnson-Jackson boundary condition to a fluidized bed I am looking at. I am using twoPhaseEulerFoam in OF 2.3.0.

In searching I have seen instances of a JohnsonJacksonParticleSlip boundary condition being used but cannot find any specifics. Can anyone clarify if this is readily available in OF or if I need to create a new boundary condition to do this?

Thank you in advance for any help you can provide

[cesarvenier]

Hi,

The main discussion of this topic may be found here:
www.cfd-online.com/Forums/openfoam-solving/60325-new-boundary-condition.html

In case you have already seen it and want to know how to set your fields with J&J BC, I recommend something like this:

For the particles velocity:

boundaryField
{
...
walls
{
type particleSlipJohnsonJackson;
specularityCoefficient 0.6;
value $internalField;
}
...
}

For the granular temperature:

boundaryField
{
...
walls
{
type particleThetaJohnsonJackson;
specularityCoefficient 0.6;
restitutionCoefficient 0.8;
value $internalField;
}
...
}

Make sure to specify the specularity and restitution coefficients as they are needed by the model.

I'm using OF2.2.0 so I'm not sure if J&J BC are currently implemented in the latest release of OF. You can check it by running your case with these BC for the Usolid and Theta. If the J&J lib is not compiled, OF will tell you.

I've wrote my J&J BC patches based on the work of Alberto Passalacqua.

Hope it helps!...tell me how it went.

[ANacc]

Thank you for the link to that thread.

Shortly after my original post I found the J-J boundary conditions implemented in 2.3.x and have been using them since.

Anthony

[mwaqas]

Hello everyone

The above discussion is very useful.
I am also working with fluidizedBed/RAS (twoPhaseEuler) model. I am using 2.3.0. In my simulation I want to apply partially slip Johanson Jackson boundary condition. my theta.Particles and U.particles boundary conditions (in 0 folder) are as follows

for theeta
walls
{
type zeroGradient;
}

for U.particles
walls
{
type fixedValue;
value uniform (0 0 0);
}

does it mean no slip boundary condition ? If it is so, how can I define partially slip boundary condition. Also where to define "Specularity coefficient"

PS: I have already tried the boundary conditions mentioned by cesarvenier in above comment, but it is not working

Thank you

[cesarvenier]

Hi,

Yes, it means that the particles phase behave like a fluid near walls (no-slip).

What kind of error do you get when you try

type particleSlipJohnsonJackson;

for the U1 field?

[mwaqas]

I got this error when I applied boundary conditions as you mentioned

Unknown patchField type particleSlipJohnsonJackson for patch type wall

Valid patchField types are :

75
(
SRFFreestreamVelocity
SRFVelocity
activeBaffleVelocity
activePressureForceBaffleVelocity
advective
atmBoundaryLayerInletVelocity
calculated
codedFixedValue
codedMixed
cyclic
cyclicACMI
cyclicAMI
cyclicSlip
cylindricalInletVelocity
directionMixed
empty
externalCoupled
fixedGradient
fixedInternalValue
fixedJump
fixedJumpAMI
fixedMean
fixedNormalSlip
fixedShearStress
fixedValue
flowRateInletVelocity
fluxCorrectedVelocity
freestream
inletOutlet
interstitialInletVelocity
kqRWallFunction
mapped
mappedField
mappedFixedInternalValue
mappedFixedPushedInternalValue
mappedFlowRate
mappedVelocityFlux
mixed
movingWallVelocity
nonuniformTransformCyclic
oscillatingFixedValue
outletInlet
outletMappedUniformInlet
outletPhaseMeanVelocity
partialSlip
pressureDirectedInletOutletVelocity
pressureDirectedInletVelocity
pressureInletOutletParSlipVelocity
pressureInletOutletVelocity
pressureInletUniformVelocity
pressureInletVelocity
pressureNormalInletOutletVelocity
processor
processorCyclic
rotatingPressureInletOutletVelocity
rotatingWallVelocity
sliced
slip
supersonicFreestream
surfaceNormalFixedValue
swirlFlowRateInletVelocity
symmetry
symmetryPlane
timeVaryingMappedFixedValue
translatingWallVelocity
turbulentInlet
uniformFixedGradient
uniformFixedValue
uniformInletOutlet
uniformJump
uniformJumpAMI
variableHeightFlowRateInletVelocity
waveTransmissive
wedge
zeroGradient
)


file: /home/waqas/OpenFOAM/waqas-2.3.0/run/standard_tutorial/fluidisedBed_RAS6/0/U.particles.boundaryField.walls from line 37 to line 19.

From function fvPatchField<Type>::New(const fvPatch&, const DimensionedField<Type, volMesh>&, const dictionary&)
in file /home/opencfd/OpenFOAM/OpenFOAM-2.3.0/src/finiteVolume/lnInclude/fvPatchFieldNew.C at line 143.

FOAM exiting


These are my boundary conditons

U.Particles

walls
{
type particleSlipJohnsonJackson;
specularityCoefficient 0.6;
value $internalField;
}

Theeta
walls
{
type particleThetaJohnsonJackson;
specularityCoefficient 0.6;
restitutionCoefficient 0.8;
value $internalField;
}

[cesarvenier]

Ok, that means that these BC patches are not included in your compilation. The patches should be located here:
twoPhaseEulerFoam/phaseCompressibleTurbulenceModels/kineticTheoryModels/derivedFvPatchFields

You need to include those when you compile your solver. Here's an example:

http://www.tfd.chalmers.se/~hani/kur...yCondition.pdf

Good luck!

[mwaqas]

Thank you very much for quick response. I go through this pdf and try to implement.

mwaqas

[spf521]

Hi Cesar,

Are you still working on OpenFOAM? Now I want to modify the JohnsonJacksonParticleSlip BC to account for the frictional strength effect of particle on BC based on the original paper of Jackson(published on 1987). Do you some idea on how to impletment this modification?

Thanks for your help!
Bryant