particleSlipJohnsonJackson
 

hi Alex,

I don not have in the case folder but I added in the system/controlDict file like:

libs ("libJohnsonJackson.so");

To solve it the solver which you use should read kineticTheoryProperties file, as it reads for example transportProperties file.
Check twoPhaseEulerFoam solver if it reads kineticTheoryProperties, and if not add reading of kineticTheoryPropertie.

No, you miss-understood:
Do you have constant/kineticTheoryProperties file in your case folder?

particleSlipJohnsonJackson
 

yes I have it in bad2/constant/kineticTheoryProperties

so, then check if twoPhaseEulerFoam reads it?

particleSlipJohnsonJackson
 

I have really stuid question is there any way to check it or I should just change the variable there and it is reading when I am using classical BC like noSlip

stupid answer: "look to the source code of twoPhaseEulerFoam")))

nice answer :) there is kinetictheoryModels only included and they are connected to the kineticTheoryProperties file

particleSlipJohnsonJackson
 

normal tutorial case(case bed2 and solver twoPhaseEulerFoam) it is running so it means it reading the kineticTheoryProperties file:

{
version 2.0;
format ascii;
class dictionary;
location "constant";
object transportProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

phasea
{
rho rho [ 1 -3 0 0 0 ] 2500;
nu nu [ 0 2 -1 0 0 ] 1e-06;
d d [ 0 1 0 0 0 0 0 ] 0.0003;
}

phaseb
{
rho rho [ 1 -3 0 0 0 ] 1.2;
nu nu [ 0 2 -1 0 0 ] 1.5e-05;
d d [ 0 1 0 0 0 0 0 ] 1;
}

Cvm Cvm [ 0 0 0 0 0 ] 0;

Cl Cl [ 0 0 0 0 0 ] 0;

Ct Ct [ 0 0 0 0 0 ] 0;

alphaAlpha alphaAlpha [ 0 0 0 0 0 ] 0;

it seems everything is ok but actually not :(

this is kineticTheoryProperties :
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "constant";
object kineticTheoryProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

kineticTheory on;

equilibrium on;

e e [ 0 0 0 0 0 0 0 ] 0.8;

alphaMax alphaMax [ 0 0 0 0 0 0 0 ] 0.62;

alphaMinFriction alphaMinFriction [ 0 0 0 0 0 0 0 ] 0.65;

Fr Fr [ 1 -1 -2 0 0 0 0 ] 0.05;

eta eta [ 0 0 0 0 0 0 0 ] 2;

p p [ 0 0 0 0 0 0 0 ] 5;

phi phi [ 0 0 0 0 0 0 0 ] 0;

viscosityModel Gidaspow;

conductivityModel Gidaspow;

granularPressureModel Lun;

frictionalStressModel JohnsonJackson;

radialModel Gidaspow;

HrenyaSinclairCoeffs
{
L L [ 0 1 0 0 0 0 0 ] 0.0005;
}

Ok, it seems that it is the problem. BC looks alphaMax in kinetic dictionary and it is renamed to transportProperties or something similar.

Try updates in particleSlipJohnsonJacksonFvPatchVectorField.C which I marked with bold, recompile library and try again to run a solver.

particleSlipJohnsonJackson
 

hi Alex,

I got the same error...:( again. this boundary conditions are made by Alberto and Juho..

best!
P.S: while compiling I got some message "statement has no effect":
SOURCE=particleSlipJohnsonJackson/particleSlipJohnsonJacksonFvPatchVectorField.C ; g++ -m64 -Dlinux64 -DWM_DP -Wall -Wextra -Wno-unused-parameter -Wold-style-cast -Wnon-virtual-dtor -O3 -DNoRepository -ftemplate-depth-100 -I/opt/openfoam210/src/triSurface/lnInclude -I/opt/openfoam210/src/meshTools/lnInclude -I/opt/openfoam210/src/finiteVolume/lnInclude -IlnInclude -I. -I/opt/openfoam210/src/OpenFOAM/lnInclude -I/opt/openfoam210/src/OSspecific/POSIX/lnInclude -fPIC -c $SOURCE -o Make/linux64GccDPOpt/particleSlipJohnsonJacksonFvPatchVectorField.o
particleSlipJohnsonJackson/particleSlipJohnsonJacksonFvPatchVectorField.C: In member function ‘virtual void Foam::particleSlipJohnsonJacksonFvPatchVectorField ::updateCoeffs()’:
particleSlipJohnsonJackson/particleSlipJohnsonJacksonFvPatchVectorField.C:125 :6: warning: statement has no effect [-Wunused-value]
'/opt/openfoam210/platforms/linux64GccDPOpt/lib/libJohnsonJackson.so' is up to date.
Recep

OpenFoam for Multiphase Solid-Gas Model
 

hi Alex
anything else you can suggest because it is really important for me to have kind of BC

thanks a lot for your help

I think you make some thing not exactly the same way I told you. Please upload here archived library folder (source files and Make directory), archived case folder you use and clarify which solver in which OF version you are trying to use. Tomorrow I will check it myself. Too busy now with my project.

Cheers!

particleSlipJohnsonJackson
 

Hi Alex,

Sorry for disturbing you. Thanks a lot...I am using version of OpenFoam210 and solver twoPhaseEulerFoam. I will send the case and particleSlipJohnsonJackson BC file as well in your email address because they are too big for here.

Thanks a lot and sorry again!!!

Best!

Recep

SpoutBed with here presented BC's
 

Hello,

I'm new at OpenFOAM and I've tried to simulate the spoutBed presented by Buijtenen. "Numerical and experimental study on multiple-spout fluidized beds"

After adapting the twoPhaseEulerFoam solver (defining a mapping field in which region the velocity of granular phase and temperature should be solved, otherwise the bed collapses allways in OF 1.7.1) and using fixedValue (0 0 0) at the wall for the granular phase I got nice and useable results, but with zero velocity at the wall.

I want to improve them by adding the boundary conditions of Schneiderbauer et. al.
("CFD study of a single-spout pseudo-2D bed: the impact of solids wall boundary conditions", 2012)
and for comparison Johnson and Joackson with spec. coeff = 0.15. Schneiderbauer et. al. BC's have the benefit that the BC's only use material parameters and are more accurate as Johnson and Jackson BC's with adaptive Specularity Coefficient calculation.

So I saw the discussion here and used the way presented by alberto and started with partialSlip boundary Condition

valueFraction*U + (1-valueFraction)*gradU*delta n = 0

The results of both BC's were nearly identically, but I also had in my bed allways velocity zero (or nearly zero) at the wall. This isn't the case in reality and measurements shown in the paper, also our CFDEM results show other tendencies, so for my opinion this useage is maybe mathematical correct, but not physically and we need physically correct results, but it give the best result for this simulation I have got now.

Now i'm searching for other ideas to implement the two boundary conditions in OpenFOAM.

For example with a mixedBoundary condition

valueFraction*(U-Uref)+ (1-valueFraction)*(gradU*delta n-gradUref*delta n) = 0,

or a calculated one. Does anybody has some knowledge about calculated boundary conditions and how the work?

David, your question is too general. There are a lot of discussions here about different BC implementations. You can create your own copying e.g. mixed typed BC from OpenFOAM and update it according to your formulation. You just have to know how your BC influences diagonal and source parts of the matrix. Another opportunity is to use swak4foam... So choice is yours...

Regarding calculated BC: boundary values are obtained from field calculations, e.g. if Field_c = Field_a <operation > Field_b, boundary values for Field_c are already calculated from previous relation and are not updated in some special way.

Thank you for the quick answer. I will have a look at swak4foam.

particleSlipJohnsonJackson
 

Hi,

I am implemeting JohnsonJackson BC , I have some problems:

walls
{
type particleThetaJohnsonJackson;
specularityCoefficient 0.5;
}

Is it correct to write like this?

thanks ...

Best!

Recep

Don't forget
 

to inlcude the value !!

walls
{
type particleSlipJohnsonJackson;
specularityCoefficient 0.5;
value uniform(0 0 0);
}

I think, this was already discussed here in the CFD-online forum.

regards,
David