LASER Log Out | Topics | Search
Moderators | Register | Edit Profile

FlexPDE User's Forum » User Postings » LASER « Previous Next »

Author Message
Top of pagePrevious messageNext messageBottom of page Link to this 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.
Top of pagePrevious messageNext messageBottom of page Link to this message

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.
Top of pagePrevious messageNext messageBottom of page Link to this message

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.
Top of pagePrevious messageNext messageBottom of page Link to this message

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.

Add Your Message Here
Post:
Username: Posting Information:
This is a private posting area. Only registered users and moderators may post messages here.
Password:
Options: Enable HTML code in message
Automatically activate URLs in message
Action:

Topics | Last Day | Last Week | Tree View | Search | Help/Instructions | Program Credits Administration