CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Programming & Development

Boundary Condition for granular temperature Theta

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 1, 2012, 21:08
Default Boundary Condition for granular temperature Theta
  #1
Member
 
Charlie
Join Date: Dec 2010
Location: 415 Kinross Dr. Newark, DE 19711
Posts: 78
Rep Power: 6
cheng1988sjtu is on a distinguished road
Hi,

I'm using twoPhaseEulerFoam, including the kinetic theory. I've decided to create a file to solve the PDE for Theta, i.e the granular temperature, that is, I don't want to use the existing kinetic theory in OpenFOAM, basically, the work is done, and the compilation is OK, however, when I tried to run a case, the error message came out.

The following is how I did that, it's simple, I just create a *.H file, which contain the fvScalarMatrix for Theta:

/**************************/
fvScalarMatrix ThetaEqn
(
fvm::ddt(ss,Theta)
+ ...
==
...
);

ThetaEqn.relax();
ThetaEqn.solve();
/****************************/

The compilation is OK, and the initial file for Theta is defined as:
walls
{
type fixedValue;
uniform 0;
}
inandouthalf11
{
type cyclic;
}
inandouthalf12
{
type cyclic;
}
inandouthalf21
{
type cyclic;
}
inandouthalf22
{
type cyclic;
}
top
{
type fixedValue;
uniform 0;
}
frontAndBackPlanes
{
type empty;
}
/************************************************** ******/

However, When I tried to run the case, the error message is that :

--> FOAM FATAL ERROR:

gradientInternalCoeffs cannot be called for a calculatedFvPatchField
on patch top of field Theta in file "/*/*/*/*/0/Theta"
You are probably trying to solve for a field with a default boundary condition.

From function calculatedFvPatchField<Type>::gradientInternalCoef fs() const
in file fields/fvPatchFields/basic/calculated/calculatedFvPatchField.C at line 186.

FOAM exiting

It seems that the problem is that I didn't define the boundary condition, however, I'm basically following Alberto's work on the boundary condition, that is JohnsonJackson's boundary condition, should I include the boundary condition explicitly in the theta equation file? Does anyone has that kind of experience?

Thanks!
cheng1988sjtu is offline   Reply With Quote

Old   June 1, 2012, 23:15
Default Problem solved
  #2
Member
 
Charlie
Join Date: Dec 2010
Location: 415 Kinross Dr. Newark, DE 19711
Posts: 78
Rep Power: 6
cheng1988sjtu is on a distinguished road
Never mind, it's just due to a stupid mistake,
in createField.H, I initiate Theta in a wrong way, it should be:

volScalarField Theta
(
IOobject
(
"Theta",
runTime.timeName(),
mesh,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
mesh
// mesh,
// dimensionedScalar("zero", dimensionSet(0, 2, -2, 0, 0), 0.0)
);

cheers!
cheng1988sjtu is offline   Reply With Quote

Reply

Tags
boundary conditions, kinetic theory, two fluid model

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
Temperature dependant mixed Boundary condition argonaut OpenFOAM Pre-Processing 2 February 15, 2011 12:02
External Radiation Boundary Condition (Two sided wall), Grid Interface CFD XUE FLUENT 0 July 8, 2010 06:49
How to apply the temperature boundary condition on the outside of the mould? alefem FLOW-3D 5 May 25, 2010 05:32
Pressure Boundary Condition Matt Umbel Main CFD Forum 0 January 11, 2002 11:06
boundary condition : temperature J.D.Yoon FLUENT 1 August 29, 2000 04:08


All times are GMT -4. The time now is 01:50.