CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

New boundary condition

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   February 28, 2012, 07:54
Default particleSlipJohnsonJackson
  #41
Senior Member
 
rkhr
Join Date: May 2011
Posts: 222
Rep Power: 7
Kanarya is on a distinguished road
hi Alex,

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

libs ("libJohnsonJackson.so");
Kanarya is offline   Reply With Quote

Old   February 28, 2012, 07:55
Default
  #42
Senior Member
 
Dr. Alexander Vakhrushev
Join Date: Mar 2009
Posts: 213
Rep Power: 10
makaveli_lcf is on a distinguished road
Send a message via ICQ to makaveli_lcf
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.
__________________
Best regards,

Dr. Alexander VAKHRUSHEV

Christian Doppler Laboratory for "Advanced Process Simulation of
Solidification and Melting"

Simulation and Modelling of Metallurgical Processes
Department of Metallurgy
University of Leoben

Franz-Josef-Str. 18
A - 8700 Leoben
Österreich / Austria
Tel.: +43 3842 - 402 - 3125
http://smmp.unileoben.ac.at
makaveli_lcf is offline   Reply With Quote

Old   February 28, 2012, 08:00
Default
  #43
Senior Member
 
Dr. Alexander Vakhrushev
Join Date: Mar 2009
Posts: 213
Rep Power: 10
makaveli_lcf is on a distinguished road
Send a message via ICQ to makaveli_lcf
No, you miss-understood:
Do you have constant/kineticTheoryProperties file in your case folder?
__________________
Best regards,

Dr. Alexander VAKHRUSHEV

Christian Doppler Laboratory for "Advanced Process Simulation of
Solidification and Melting"

Simulation and Modelling of Metallurgical Processes
Department of Metallurgy
University of Leoben

Franz-Josef-Str. 18
A - 8700 Leoben
Österreich / Austria
Tel.: +43 3842 - 402 - 3125
http://smmp.unileoben.ac.at
makaveli_lcf is offline   Reply With Quote

Old   February 28, 2012, 08:03
Default particleSlipJohnsonJackson
  #44
Senior Member
 
rkhr
Join Date: May 2011
Posts: 222
Rep Power: 7
Kanarya is on a distinguished road
yes I have it in bad2/constant/kineticTheoryProperties
Kanarya is offline   Reply With Quote

Old   February 28, 2012, 08:04
Default
  #45
Senior Member
 
Dr. Alexander Vakhrushev
Join Date: Mar 2009
Posts: 213
Rep Power: 10
makaveli_lcf is on a distinguished road
Send a message via ICQ to makaveli_lcf
so, then check if twoPhaseEulerFoam reads it?
__________________
Best regards,

Dr. Alexander VAKHRUSHEV

Christian Doppler Laboratory for "Advanced Process Simulation of
Solidification and Melting"

Simulation and Modelling of Metallurgical Processes
Department of Metallurgy
University of Leoben

Franz-Josef-Str. 18
A - 8700 Leoben
Österreich / Austria
Tel.: +43 3842 - 402 - 3125
http://smmp.unileoben.ac.at
makaveli_lcf is offline   Reply With Quote

Old   February 28, 2012, 08:08
Default particleSlipJohnsonJackson
  #46
Senior Member
 
rkhr
Join Date: May 2011
Posts: 222
Rep Power: 7
Kanarya is on a distinguished road
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
Kanarya is offline   Reply With Quote

Old   February 28, 2012, 08:15
Default
  #47
Senior Member
 
Dr. Alexander Vakhrushev
Join Date: Mar 2009
Posts: 213
Rep Power: 10
makaveli_lcf is on a distinguished road
Send a message via ICQ to makaveli_lcf
stupid answer: "look to the source code of twoPhaseEulerFoam")))
__________________
Best regards,

Dr. Alexander VAKHRUSHEV

Christian Doppler Laboratory for "Advanced Process Simulation of
Solidification and Melting"

Simulation and Modelling of Metallurgical Processes
Department of Metallurgy
University of Leoben

Franz-Josef-Str. 18
A - 8700 Leoben
Österreich / Austria
Tel.: +43 3842 - 402 - 3125
http://smmp.unileoben.ac.at
makaveli_lcf is offline   Reply With Quote

Old   February 28, 2012, 08:34
Default
  #48
Senior Member
 
rkhr
Join Date: May 2011
Posts: 222
Rep Power: 7
Kanarya is on a distinguished road
nice answer there is kinetictheoryModels only included and they are connected to the kineticTheoryProperties file
Kanarya is offline   Reply With Quote

Old   February 28, 2012, 08:47
Default particleSlipJohnsonJackson
  #49
Senior Member
 
rkhr
Join Date: May 2011
Posts: 222
Rep Power: 7
Kanarya is on a distinguished road
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
Kanarya is offline   Reply With Quote

Old   February 28, 2012, 08:49
Default
  #50
Senior Member
 
rkhr
Join Date: May 2011
Posts: 222
Rep Power: 7
Kanarya is on a distinguished road
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;
}
Kanarya is offline   Reply With Quote

Old   February 28, 2012, 09:01
Default
  #51
Senior Member
 
Dr. Alexander Vakhrushev
Join Date: Mar 2009
Posts: 213
Rep Power: 10
makaveli_lcf is on a distinguished road
Send a message via ICQ to makaveli_lcf
Ok, it seems that it is the problem. BC looks alphaMax in kinetic dictionary and it is renamed to transportProperties or something similar.

Quote:
// Update the coefficients associated with the patch field
void particleSlipJohnsonJacksonFvPatchVectorField::upda teCoeffs()
{
if (updated())
{
return;
}

if ((specularityCoefficient_ < 0) || (specularityCoefficient_ > 1))
{
FatalErrorIn
(
"particleSlipJohnsonJacksonFvPatchScalarField: :"
"updateCoeffs()"
) << "The value of the specularity coefficient has to be between 0 and 1."
<< abort(FatalError);
}

const dictionary& transportProperties = db().lookupObject<IOdictionary>
(
"transportProperties"
);

// const dictionary& kineticTheoryProperties = db().lookupObject<IOdictionary> ("kineticTheoryProperties");

dictionary phaseaDictionary
(
transportProperties.subDict("phasea")
);

dimensionedScalar rhoa(phaseaDictionary.lookup("rho"));
// dimensionedScalar alphaMax(kineticTheoryProperties.lookup("alphaMax" ));
dimensionedScalar alphaMax(transportProperties.lookup("alphaMax"));


const fvPatchScalarField& alpha =
patch().lookupPatchField<volScalarField, scalar>("alpha");

const fvPatchScalarField& g0 =
patch().lookupPatchField<volScalarField, scalar>("gs0");

const fvPatchScalarField& mua =
patch().lookupPatchField<volScalarField, scalar>("mua");

scalarField alphaPatch = alpha.patchInternalField() + 1.0e-6;
scalarField ThetaPatch = max(alphaPatch, 1.0e-6);
scalarField g0Patch = g0.patchInternalField();
scalarField muaPatch = mua.patchInternalField();

if (db().foundObject<volScalarField>("Theta"))
{
const fvPatchScalarField& Theta =
patch().lookupPatchField<volScalarField, scalar>("Theta");

ThetaPatch = Theta.patchInternalField();
}

// The partial slip BC in OpenFOAM is implemented as
//
// valueFraction*U + (1-valueFraction)*grad(U) = 0
//
// To find valueFraction, we re-write Johnson and Jackson BC as
//
// c*U + grad(U) = 0
//
// where
//
// c = valueFraction/(1 - valueFraction)
//
// As a consequence
//
// valueFraction = 1/(c + 1)

scalarField c = (6.0*muaPatch*alphaMax.value())/
(M_PI*rhoa.value()*alphaPatch*g0Patch*
specularityCoefficient_*sqrt(3.0*ThetaPatch));

this->valueFraction() = scalar(1)/(c + scalar(1));

partialSlipFvPatchVectorField::updateCoeffs();
}
Try updates in particleSlipJohnsonJacksonFvPatchVectorField.C which I marked with bold, recompile library and try again to run a solver.
__________________
Best regards,

Dr. Alexander VAKHRUSHEV

Christian Doppler Laboratory for "Advanced Process Simulation of
Solidification and Melting"

Simulation and Modelling of Metallurgical Processes
Department of Metallurgy
University of Leoben

Franz-Josef-Str. 18
A - 8700 Leoben
Österreich / Austria
Tel.: +43 3842 - 402 - 3125
http://smmp.unileoben.ac.at

Last edited by makaveli_lcf; February 28, 2012 at 09:02. Reason: typo
makaveli_lcf is offline   Reply With Quote

Old   February 28, 2012, 09:52
Default particleSlipJohnsonJackson
  #52
Senior Member
 
rkhr
Join Date: May 2011
Posts: 222
Rep Power: 7
Kanarya is on a distinguished road
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:articleSlipJohnsonJacksonFvPatchVectorField ::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
Kanarya is offline   Reply With Quote

Old   February 28, 2012, 12:52
Unhappy OpenFoam for Multiphase Solid-Gas Model
  #53
Senior Member
 
rkhr
Join Date: May 2011
Posts: 222
Rep Power: 7
Kanarya is on a distinguished road
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
Kanarya is offline   Reply With Quote

Old   February 28, 2012, 13:46
Default
  #54
Senior Member
 
Dr. Alexander Vakhrushev
Join Date: Mar 2009
Posts: 213
Rep Power: 10
makaveli_lcf is on a distinguished road
Send a message via ICQ to makaveli_lcf
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!
__________________
Best regards,

Dr. Alexander VAKHRUSHEV

Christian Doppler Laboratory for "Advanced Process Simulation of
Solidification and Melting"

Simulation and Modelling of Metallurgical Processes
Department of Metallurgy
University of Leoben

Franz-Josef-Str. 18
A - 8700 Leoben
Österreich / Austria
Tel.: +43 3842 - 402 - 3125
http://smmp.unileoben.ac.at
makaveli_lcf is offline   Reply With Quote

Old   February 28, 2012, 14:02
Smile particleSlipJohnsonJackson
  #55
Senior Member
 
rkhr
Join Date: May 2011
Posts: 222
Rep Power: 7
Kanarya is on a distinguished road
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
Kanarya is offline   Reply With Quote

Old   March 20, 2012, 02:42
Default SpoutBed with here presented BC's
  #56
New Member
 
David Schellander
Join Date: Jan 2010
Posts: 5
Rep Power: 7
Lindemann1712 is on a distinguished road
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?
Lindemann1712 is offline   Reply With Quote

Old   March 20, 2012, 03:44
Default
  #57
Senior Member
 
Dr. Alexander Vakhrushev
Join Date: Mar 2009
Posts: 213
Rep Power: 10
makaveli_lcf is on a distinguished road
Send a message via ICQ to makaveli_lcf
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.
__________________
Best regards,

Dr. Alexander VAKHRUSHEV

Christian Doppler Laboratory for "Advanced Process Simulation of
Solidification and Melting"

Simulation and Modelling of Metallurgical Processes
Department of Metallurgy
University of Leoben

Franz-Josef-Str. 18
A - 8700 Leoben
Österreich / Austria
Tel.: +43 3842 - 402 - 3125
http://smmp.unileoben.ac.at
makaveli_lcf is offline   Reply With Quote

Old   March 20, 2012, 03:50
Default
  #58
New Member
 
David Schellander
Join Date: Jan 2010
Posts: 5
Rep Power: 7
Lindemann1712 is on a distinguished road
Thank you for the quick answer. I will have a look at swak4foam.
Lindemann1712 is offline   Reply With Quote

Old   March 26, 2012, 10:14
Default particleSlipJohnsonJackson
  #59
Senior Member
 
rkhr
Join Date: May 2011
Posts: 222
Rep Power: 7
Kanarya is on a distinguished road
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
Kanarya is offline   Reply With Quote

Old   March 27, 2012, 03:35
Default Don't forget
  #60
New Member
 
David Schellander
Join Date: Jan 2010
Posts: 5
Rep Power: 7
Lindemann1712 is on a distinguished road
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
Lindemann1712 is offline   Reply With Quote

Reply

Tags
heat transfer, new bc

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Boundary condition for UDS Tomik FLUENT 0 December 5, 2006 18:37
Boundary condition of the third kind or Danckwertz boundary condition plage OpenFOAM Running, Solving & CFD 4 October 3, 2006 12:21
Slip Boundary Condition for Moving Boundary Shukla Main CFD Forum 3 November 11, 2005 16:02
UDF boundary condition Jeff FLUENT 2 November 20, 2003 18:15
Boundary Condition in LES Zhang Tsiang Main CFD Forum 3 February 5, 2002 21:15


All times are GMT -4. The time now is 13:25.