CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Programming & Development

Help to implement a new code to cure kinetics

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 17, 2012, 08:17
Default Help to implement a new code to cure kinetics
  #1
New Member
 
Victor
Join Date: Jul 2012
Location: Madrid
Posts: 5
Rep Power: 3
vikthor is on a distinguished road
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!!
Attached Files
File Type: h createFields.H (1.9 KB, 8 views)
File Type: h degreeEqn.H (203 Bytes, 11 views)
File Type: c myLaplacianFoam.C (2.3 KB, 9 views)
vikthor is offline   Reply With Quote

Old   July 17, 2012, 11:12
Default
  #2
Senior Member
 
akidess's Avatar
 
Anton Kidess
Join Date: May 2009
Location: Delft, Netherlands
Posts: 912
Rep Power: 16
akidess will become famous soon enough
Quote:
[...] I have problems to get it [...]
I'm having troubles guessing what those problems were...
__________________
*On twitter @akidTwit
*Spend as much time formulating your questions as you expect people to spend on their answer.
*Check out the scientific computing exchange http://scicomp.stackexchange.com
akidess is offline   Reply With Quote

Old   July 17, 2012, 11:34
Default
  #3
New Member
 
Victor
Join Date: Jul 2012
Location: Madrid
Posts: 5
Rep Power: 3
vikthor is on a distinguished road
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??
vikthor is offline   Reply With Quote

Old   July 17, 2012, 12:09
Default
  #4
Senior Member
 
akidess's Avatar
 
Anton Kidess
Join Date: May 2009
Location: Delft, Netherlands
Posts: 912
Rep Power: 16
akidess will become famous soon enough
If your assumption is correct, it should be solved by using:
Code:
fvm::ddt(T) - fvm::laplacian(DT, T) - fvc::ddt(DC, curedegree)
__________________
*On twitter @akidTwit
*Spend as much time formulating your questions as you expect people to spend on their answer.
*Check out the scientific computing exchange http://scicomp.stackexchange.com
akidess is offline   Reply With Quote

Old   July 17, 2012, 13:47
Default
  #5
New Member
 
Victor
Join Date: Jul 2012
Location: Madrid
Posts: 5
Rep Power: 3
vikthor is on a distinguished road
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.
vikthor is offline   Reply With Quote

Old   July 17, 2012, 16:18
Default
  #6
New Member
 
Helmut Roth
Join Date: Mar 2009
Posts: 23
Rep Power: 7
helmut is on a distinguished road
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?
helmut is offline   Reply With Quote

Old   July 18, 2012, 07:35
Default
  #7
New Member
 
Victor
Join Date: Jul 2012
Location: Madrid
Posts: 5
Rep Power: 3
vikthor is on a distinguished road
Hi Helmut

I linearized the equation N(curedegree,T) and run, but doesn't solve the curedegree, only solve the temperature field.
vikthor is offline   Reply With Quote

Reply

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
The FOAM Documentation Project - SHUT-DOWN holger_marschall OpenFOAM 242 March 7, 2013 12:30
Debugging Unsteady 2-D Panel Method Code: Wake Modeling RajeshAero Main CFD Forum 5 November 10, 2011 05:48
Open Source Vs Commercial Software MechE OpenFOAM 28 May 16, 2011 11:02
Design Integration with CFD? John C. Chien Main CFD Forum 19 May 17, 2001 15:56
What is the Better Way to Do CFD? John C. Chien Main CFD Forum 54 April 23, 2001 08:10


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