- **OpenFOAM Programming & Development**
(*https://www.cfd-online.com/Forums/openfoam-programming-development/*)

- - **develope a solver**
(*https://www.cfd-online.com/Forums/openfoam-programming-development/123797-develope-solver.html*)

develope a solverHi everyone
I want to develop a solver by adding the following Eq to simpleFoam solver: http://latex.codecogs.com/gif.latex?...wn%20T%29/T%5D http://latex.codecogs.com/gif.latex?%5Crho,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` http://latex.codecogs.com/gif.latex?%5Crho [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:
can anybody help? Thanks Arjang |

Quote:
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. |

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. |

All times are GMT -4. The time now is 10:47. |