# develope a solver

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

 September 21, 2013, 07:21 develope a solver #1 Member   Arjang Behnoud Join Date: Oct 2012 Posts: 54 Rep Power: 3 Hi everyone I want to develop a solver by adding the following Eq to simpleFoam solver: $\rho cV.\triangledown T=\triangledown .k\triangledown T+\rho c[D\triangledown L.\triangledown T+S(\triangledown T.\triangledown T)/T]$ $\rho$,c,k,D and S are constant. V:velocity vector T and L are scalars. I've considered these assumptions: m=fvc::grad(L) n=fvc::grad(T)/T and define m and n in createfields and finally write TEqn as bellow: Code:  fvScalarMatrix TEqn ( rho*c*fvm::div(phi,T) -fvm::laplacian(k,T) ); solve(TEqn == rho*c*T*(D*(m&n)+S*(n&n))); m=fvc::grad(V); n=fvc::grad(T)/T; dimensions of variable and constants are listed below: $\rho$ [1 -3 0 0 0 0 0 ] c [0 2 -2 -1 0 0 0 ] D&S [0 2 -1 0 0 0 0] k [1 1 -3 -1 0 0 0] variable L is volume fraction and has no dimension. variable T is temperature:[0 0 0 1 0 0 0] I've made this new solver with wmake command successfully but when I run the case and want to complete the solver I confronted the following Fatal Error: Code: --> FOAM Warning : From function dlOpen(const fileName&, const bool) in file POSIX.C at line 1175 dlopen error : libsimpleSwakFunctionObjects.so: cannot open shared object file: No such file or directory --> FOAM Warning : From function dlLibraryTable::open(const fileName&, const bool) in file db/dynamicLibrary/dlLibraryTable/dlLibraryTable.C at line 96 could not load "libsimpleSwakFunctionObjects.so" --> FOAM FATAL ERROR: incompatible dimensions for operation [T[1 -1 -3 0 0 0 0] ] - [T[1 -1 -3 2 0 0 0] ] From function checkMethod(const fvMatrix&, const fvMatrix&) in file /opt/openfoam211/src/finiteVolume/lnInclude/fvMatrix.C at line 1316. FOAM aborting #0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #1 Foam::error::abort() in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #2 void Foam::checkMethod(Foam::fvMatrix const&, Foam::fvMatrix const&, char const*) in "/home/amir/OpenFOAM/amir-2.1.1/platforms/linux64GccDPOpt/bin/nanoHeatSimpleFoam" #3 at nanoHeatSimpleFoam.C:0 #4 in "/home/amir/OpenFOAM/amir-2.1.1/platforms/linux64GccDPOpt/bin/nanoHeatSimpleFoam" #5 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #6 in "/home/amir/OpenFOAM/amir-2.1.1/platforms/linux64GccDPOpt/bin/nanoHeatSimpleFoam" Aborted (core dumped) can anybody help? Thanks Arjang

September 22, 2013, 20:11
#2
Senior Member

Dongyue Li
Join Date: Jun 2012
Location: Beijing,China
Posts: 429
Rep Power: 5
Quote:
 Originally Posted by AmirBaqa1987 Hi everyone I want to develop a solver by adding the following Eq to simpleFoam solver: $\rho cV.\triangledown T=\triangledown .k\triangledown T+\rho c[D\triangledown L.\triangledown T+S(\triangledown T.\triangledown T)/T]$ $\rho$,c,k,D and S are constant. V:velocity vector T and L are scalars.
Hi Arjang,

Im not an expert in this field. but are you sure your equation's dimension is compatible? In the following equation, I think the left's dimension is not the same with the right one.

 September 23, 2013, 01:17 #3 Member   Arjang Behnoud Join Date: Oct 2012 Posts: 54 Rep Power: 3 Hi Forrest I've found the problem. I've made a miss and set the dimension of k in the Transport Properties Dictionary wrong, I mean: K [1 1 -3 1 0 0 0] instead of k [1 1 -3 -1 0 0 0]. I modify this dictionary and my problem is solved.

 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 Peter_600 OpenFOAM 3 May 4, 2013 08:19 Smaras FLUENT 2 February 19, 2013 06:58 danvica OpenFOAM Running, Solving & CFD 16 December 22, 2012 02:09 bearcat CFX 6 April 28, 2008 14:08 ztdep OpenFOAM Running, Solving & CFD 0 September 15, 2005 10:41

All times are GMT -4. The time now is 07:33.