|November 15, 2011, 10:05||
how to define this PDE:Uxx-a^2Ut+A＊exp(-a*x)＝0
Join Date: Aug 2011
Posts: 15Rep Power: 6
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.
|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|