CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Running, Solving & CFD

how to define this PDE:Uxx-a^2Ut+A*exp(-a*x)=0

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

LinkBack Thread Tools Display Modes
Old   November 15, 2011, 10:05
Default how to define this PDE:Uxx-a^2Ut+A*exp(-a*x)=0
New Member
Join Date: Aug 2011
Posts: 15
Rep Power: 8 is on a distinguished road
hello every one, i need your help now.
i'm tring to solve this PDE:Uxx-a^2Ut+A*exp(-a*x)=0 with openFoam,

after copying the solver "laplacianFoam"to my own directory and changing it to myLaplacianFoam, i add two lines in the "main" function in the myLaplacianFoam.c like this:

volVectorField centers = T.mesh().C();
volScalarField x = (centers&vector(1,0,0));

then define the equation :

-fvm::ddt(sqr(a),T) - fvm::laplacian(T) == -A*exp(-a*x);

the parameters A and a are defined in creatFiled.H and in transportProperties dictionary file within the case directory.

when i try to use this solver, the error echoed :

Argument of trancendental function not dimensionless

how could i deal with this problem?
shoul i make the argument in the exp() function dimensioniess or define the A*exp(-a*x) as a volScalarField "f" in createFiled.H,and give "f" the original and boundary value in the case directory when i need to run the solver?

by the way ,if i want to utilize openFoam to solve some typical PDEs,where can i find the reference source?

thanks for your attention. is offline   Reply With Quote


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
Installing OF 1.6 on Mac OS X gschaider OpenFOAM Installation on Windows, Mac and other Unsupported Platforms 129 June 19, 2010 09:23
HELP----Surface Reaction UDF Ashi Fluent UDF and Scheme Programming 0 May 25, 2009 09:39
Missing math.h header Travis FLUENT 4 January 15, 2009 12:48
Free surface boudary conditions with SOLA-VOF Fan Main CFD Forum 10 September 9, 2006 12:24
UDF FOR UNSTEADY TIME STEP mayur FLUENT 3 August 9, 2006 10:19

All times are GMT -4. The time now is 06:12.