Author |
Message |
Milton Lima (milton)
New member Username: milton
Post Number: 1 Registered: 06-2004
| Posted on Tuesday, June 01, 2004 - 10:46 am: | |
Hi: I would like to use the FlexPDE (student version) for the study of pulsed laser heating. The problem is to use a small diameter laser spot at a given material surface (e.g. steel) for a given time. The laser spot has gaussian type distribution and about 100 nanoseconds pulse length. All I need is the temperature evolution, heat fluxes and gradients, before and after the laser pulse. All my attempts to created a code failed. I would be most grateful if someone could help me with this project. The collaborator will be cited, together with PDE solutions, in my report. Thank you, Milton.
|
Robert G. Nelson (rgnelson)
Moderator Username: rgnelson
Post Number: 150 Registered: 06-2003
| Posted on Tuesday, June 01, 2004 - 04:29 pm: | |
You may be encountering trouble resolving the small spot, because of the mesh size constraints of the student version. Send me your descriptor, and I will try to determine the cause of the trouble. |
Milton Lima (milton)
New member Username: milton
Post Number: 2 Registered: 06-2004
| Posted on Friday, June 04, 2004 - 07:05 am: | |
Dear Robert: Thank you for your reply. I am trying to use the Moving Heat Source problem of the Fields of Physics workbook. The original code looks like this (example 203, page 211): TITLE' Steel Bar Heated by Moving Torch' SELECT errlim=1e-3 spectral_colors VARIABLES temp(1000) DEFINITIONS Lx=1.0 Ly=0.2 d0=0.1 vx= 0e-3 { Torch velocity } heat=1000 tempi=300 k=82 rcp=7.87e3*449 { Steel } fluxd_x=-k*dx(temp) fluxd_y=-k*dy(temp) fluxd=vector( fluxd_x, fluxd_y) fluxdm=magnitude( fluxd) fluxd0= -1e6* [ ustep( vx*t+d0- x)- ustep( vx*t- x) ] INITIAL VALUES temp= tempi EQUATIONS div( fluxd)+ rcp*dt( temp)=heat BOUNDARIES region 'domain' start (0,0) value(temp)=tempi line to (Lx,0) natural(temp)=0 line to (Lx,Ly) natural(temp)=fluxd0 line to (0,Ly) natural(temp)=0 line to finish TIME from 0 to 800 PLOTS for t=50 by 50 to 800 elevation( temp) from (0,Ly) to (Lx,Ly) contour( temp) painted vector( fluxd) norm END 3907 Any attempts to change this code to the pulsed heat source did not work. Could you please give me help? Thank you, Milton.
|
Robert G. Nelson (rgnelson)
Moderator Username: rgnelson
Post Number: 157 Registered: 06-2003
| Posted on Friday, June 04, 2004 - 03:58 pm: | |
You showed me what Backstrom did, but you didn't show me what YOU did, or define what "didn't work" means. Send me your descriptor. |
|