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

Like Tree1Likes
  • 1 Post By Achuth

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: 7
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, 21 views)
File Type: h degreeEqn.H (203 Bytes, 26 views)
File Type: c myLaplacianFoam.C (2.3 KB, 28 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: Germany
Posts: 1,262
Rep Power: 23
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.
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: 7
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: Germany
Posts: 1,262
Rep Power: 23
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.
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: 7
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: 10
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: 7
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

Old   May 16, 2017, 10:32
Default
  #8
New Member
 
Achuth Nair Balachandran Nair
Join Date: Apr 2017
Location: Austria
Posts: 3
Rep Power: 2
Achuth is on a distinguished road
You can solve for the degree of cure by means of an explicit finite difference scheme. If you know the parameters of the equation, you can just put them as a separate header file to solve for the rate of cure and this can be implemented in the solver.

So, in general the equation will look like this in OpenFOAM.

fvc::ddt(c) == A*exp(-Ea/R*T)*pow(c,m)*pow((1-c),n)

With A, Ea, R, m, n known and with an initial conversion c0 (say 10^-20) you could solve the rate of cure explicitly.

fvc - stands for explicitly solving the equation.
Achuth is offline   Reply With Quote

Old   August 3, 2017, 14:18
Default Composite curring simulation in Openfoam
  #9
New Member
 
Nitish Kumar
Join Date: Aug 2017
Posts: 3
Rep Power: 2
nkv1990 is on a distinguished road
I am new to the openfoam and i want to simulate in curring of composite but dont know which solver i have to start and how to incorporate the cure kinetics in the solver for solving it implicitly. Can any help me in this regard.
nkv1990 is offline   Reply With Quote

Old   August 3, 2017, 14:19
Default Composite curring simulation in Openfoam
  #10
New Member
 
Nitish Kumar
Join Date: Aug 2017
Posts: 3
Rep Power: 2
nkv1990 is on a distinguished road
I am new to the openfoam and i want to simulate curring of composite but dont know which solver i have to start with and how to incorporate the cure kinetics in the solver for solving it implicitly. Can anybody help me in this regard.
nkv1990 is offline   Reply With Quote

Old   August 3, 2017, 15:54
Default
  #11
New Member
 
Achuth Nair Balachandran Nair
Join Date: Apr 2017
Location: Austria
Posts: 3
Rep Power: 2
Achuth is on a distinguished road
I have solved the problem I had mentioned earlier in this thread.

It depends on which kind of problem you want to solve.

If you want to just solve the static cure kinetics problem, then take laplacianFoam and work with it.

If you want to work with mold filling problems, take the interFoam solver and take on of the phases as resin and the other as air. The alphaEqn.H solves the phase fraction equation. Then you will have to implement the conservation of energy equation in terms of temperature. Then, integrate the transport equation for the degree of cure.

Its a long process to understand the physics behind the mold filling problem. I am currently working at a research facility and made the solver for mold filling process for LCM process. But I cannot provide the solver. Its strictly for our research facility.

If you want any help I can help you.
nkv1990 likes this.
Achuth is offline   Reply With Quote

Old   August 4, 2017, 08:46
Default
  #12
New Member
 
Nitish Kumar
Join Date: Aug 2017
Posts: 3
Rep Power: 2
nkv1990 is on a distinguished road
Thanks for reply .... let me search some more details on it
nkv1990 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 13:30
Debugging Unsteady 2-D Panel Method Code: Wake Modeling RajeshAero Main CFD Forum 5 November 10, 2011 06: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 19:44.