laplacianFoam giving weird results:HELP

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

 June 13, 2013, 11:47 laplacianFoam giving weird results:HELP #1 Member   shiv Join Date: Jun 2012 Location: Lucknow, IN Posts: 51 Rep Power: 6 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.

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post deeimproptus FLUENT 0 October 28, 2011 00:25 rahulrp OpenFOAM Pre-Processing 6 April 8, 2010 00:03 dakos OpenFOAM 1 November 14, 2009 07:15 maxime FLUENT 0 June 13, 2007 12:01 Sireesha FLUENT 1 July 10, 2004 17:54

All times are GMT -4. The time now is 21:23.