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

laplacianFoam giving weird results:HELP

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 13, 2013, 11:47
Default laplacianFoam giving weird results:HELP
  #1
Member
 
shiv
Join Date: Jun 2012
Location: Lucknow, IN
Posts: 51
Rep Power: 5
shash is on a distinguished road
hi all,

I am using laplacianfoam to simulate the bratu problem laplacian(b)+epsilon*exp(b/(1+mu*b))=0, thou i am getting expected variation on applying BCs: b=100 on both sides of my rectangular domain but on applying b=0 i am getting no changes(why??).
In my new solver i have linearised the exponential term by taylor exp. and have not removed d/dt term of laplacianfoam instead specified steadyState in fvSchemes.

0/b :

FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object b;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 0 0 0 0 0];
internalField uniform 100;
boundaryField
{
patch_1
{
type zeroGradient;
}
patch_2
{
type fixedValue;
value uniform 100;
}
patch_3
{
type zeroGradient;
}
patch_4
{
type fixedValue;
value uniform 100;
}}

and my solver

dimensionedScalar one
(
"one",
dimensionSet(0, 2, -1, 0, 0, 0 ,0),
1.0
);

simpleControl simple(mesh);

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

Info<< "\nCalculating distribution of b\n" << endl;

while (simple.loop())
{
Info<< "Time = " << runTime.timeName() << nl << endl;

while (simple.correctNonOrthogonal())
{
solve
(
fvm::ddt(b) - fvm::laplacian(one, b) == fvm::Sp(epsilon*((2 + fvc::Sp(1 +2*mu, b)))/(2*(1+fvc::Sp(mu, b))), b)
);
}

#include "write.H"

Last edited by shash; June 13, 2013 at 13:59.
shash is offline   Reply With Quote

Reply

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
Solution in FLuent is giving weird results deeimproptus FLUENT 0 October 28, 2011 00:25
create own case for laplacianFoam rahulrp OpenFOAM Pre-Processing 6 April 8, 2010 00:03
laplacianfoam with solidification dakos OpenFOAM 1 November 14, 2009 07:15
Weird pathlines on periodic sides maxime FLUENT 0 June 13, 2007 12:01
Problem in giving heat source Sireesha FLUENT 1 July 10, 2004 17:54


All times are GMT -4. The time now is 00:15.