CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Programming & Development (http://www.cfd-online.com/Forums/openfoam-programming-development/)
-   -   Help to implement a new code to cure kinetics (http://www.cfd-online.com/Forums/openfoam-programming-development/104857-help-implement-new-code-cure-kinetics.html)

vikthor July 17, 2012 08:17

Help to implement a new code to cure kinetics
 
3 Attachment(s)
Hi foamers!

I'm modifying the code of laplacianFoam to calculate the cure kinetics, but I have problems to get it. I believe that the problem is that I unknow how openFoam solve the ecuations. If someone could help me, I would appreciate it. I attached the 3 archives I have modified.

Thanks a lot!!

akidess July 17, 2012 11:12

Quote:

[...] I have problems to get it [...]
I'm having troubles guessing what those problems were...

vikthor July 17, 2012 11:34

The error displayed is

--> FOAM FATAL ERROR:
incompatible fields for operation
[T] - [curedegree]

I believe the error appears because I'm solving one equation with two variables. Some ideas??

akidess July 17, 2012 12:09

If your assumption is correct, it should be solved by using:
Code:

fvm::ddt(T) - fvm::laplacian(DT, T) - fvc::ddt(DC, curedegree)

vikthor July 17, 2012 13:47

thanks Anton
now it runs and calculate correctly de temperature field, but it doesn´t solve the curedegree (the initial and final residual of all steps 0 and the number of iteracion in each step is 0). It seems that ignore the cure kinetics terms and solve like the laplacianFoam.

helmut July 17, 2012 16:18

Looks like your degree eqn has the form
ddt(curedegree) - N(cureDegree,T) = 0,
where N is a nonlinear function. Maybe some linearizing is required?

vikthor July 18, 2012 07:35

Hi Helmut

I linearized the equation N(curedegree,T) and run, but doesn't solve the curedegree, only solve the temperature field.


All times are GMT -4. The time now is 15:37.