Discussion Closed This discussion was created more than 6 months ago and has been closed. To start a new discussion with a link back to this one, click here.
heating a block of material by a moving laser
Posted 11.08.2014, 07:39 GMT-4 Wave Optics, Heat Transfer & Phase Change, Materials, Modeling Tools & Definitions, Parameters, Variables, & Functions Version 5.2 10 Replies
Please login with a confirmed email address before reporting spam
I'm a beginner at COMSOL, so probably my question is quite stupid to the frequent users of COMSOL.
I want to simulate a laser hardening process. At this time only the temperature field (not yet the phase thansformations etc.). Therefore I want to simulate a block of material (160 mm long, 40 mm wide, 8 mm high), and a laser power of 450W, with a spotsize of 2.7 mm, moving at a speed of 100 mm/min.
The laser is comming out of a fiber thus with a tophat profile (not gaussian).
I started from the laser heating tutorial (silicon wafer).
I tried to make instead of the gaussian profile as in the tutorial, a top hat distribution. I created a rectangular function, with limits -r_spot/2 and +r_spot/2. The analytic function used for the power distribution is than: p_laser*rect1(x)*rect1(y)
First problem I have: this is not a top hat profile!
www.holoor.co.il/images/image1.jpg
Left: the gaussian profile, right upper image: what I have, right lower image: what I want...
So the analytic function should be different, using a circular formula. Any idea how to do this?
For the heat transfer in solids (ht), I have added the translational motion, with v_feed in x direction.
Attached you can find the model.
I have 2 problems: to small temperature differences, and no moving laser source...
Can anyone help me out please?
Best regards,
Jan
Attachments:
Please login with a confirmed email address before reporting spam
-I removed the rectangle and used the gaussian pulse (suppose a gaussian beam instead of the wanted tophat), as in the example. Change the expression for the analytic component into:
p_laser*gp(x-v_feed*t)*gp(y) , with gp1: the gaussian pulse.
So far it's looking promising. The results (about 1100°C) are realistic.
However, when looking at the attached image comsol2, I believe it is strange there is this variation in the max. temperature. Why is this?
So, the moving problem is fixed. But what about how to create a Top Hat distribution. Anyone an idea?
Attached the current model.
Attachments:
Please login with a confirmed email address before reporting spam
About the variation of the maximum: I have not counted but it seems like the amount of peaks in the t-Tmax graph is roughly equal to the amount of mesh elements encountered along the length of the bar. Try a finer mesh and see if it makes any difference.
Please login with a confirmed email address before reporting spam
I would like to create a model in which heat transfer to the double layered material is done using a laser moving in meander shape. Can any one please help me out in creating a meander shape moment of laser for heat transfer. any ideas how to do??
with regards;
meera
Please login with a confirmed email address before reporting spam
Please login with a confirmed email address before reporting spam
Thus moving in Y+ direction, at for example 200 mm/min, and oscillating in the X+ and - direction at 100 Hz. The laser in my case has a tophat power distribution. But a gaussian is of course also possible.
Is this what you mean, or what do you mean by meander?
Please login with a confirmed email address before reporting spam
for example my material is of length 1cm and width 1 cm, and my laser has to move in straight line in y-direncton example 0.75cm/min and then move straightly in x-direction for example 0.25cm/min, then again move in a straight line y-direction same 0.75cm/min, then again move in x-direction for 0.25cm/min... like wise.. can you please suggest me some ideas
Please login with a confirmed email address before reporting spam
for example x: x-v_feed*t for [x:... to ...] etc.?
Please login with a confirmed email address before reporting spam
I have seen few examples like , "Heating with a moving laser", here the laser is moving in circular direction and thez draw that circle on the wafer, but i didnt understand how they represented in the form of equatio in comsol. And another example "laser heating of silicon wafer", here the wafer is rotating and the laser is also moving along its diameter, but here also they didnt explain the moment of laser in the form of equation in comsol...
Is their any better to understand and represent the equations. I am stuck at the begining only..
Please login with a confirmed email address before reporting spam
Your laser light will be absorbed in the material on which it shines, acting as a heat source. The amount of heat generated is proportional to the light absorption, which is proportional to the light intensity. In the horizontal direction the intensity of your laser beam can have any distribution, although gaussian is quite common. In the vertical direction (the light propagation direction) the light will decay because of absorption. In formulas now the strenght of your heat source can be described as
A*exp(B*(-x^2-y^2))*exp(C*z), where A is related to the light intensity, B defines the width of your beam and C determines how fast the light is absorbed. exp(B*(-x^2-y^2)) is the gaussian distribution and exp(C*z) is the exponential decay.
If it is described like this the center of your beam will be where x=0 and y=0. If we now want to shift it a distance d in the positive x direction we can write it as A*exp(B*(-(x-d)^2-y^2))*exp(C*z).
If we want to shift it a distance d per second we can write it as A*exp(B*(-(x-d*t)^2-y^2))*exp(C*z).
x, y and z are defined in comsol, and t also if you do a time dependent study.
Now if you want to want the beam to move in a square you can write it as A*exp(B*(-(x-fun1(t))^2-(y-fun2(t)^2))*exp(C*z) and you have to define the functions fun1 and fun2 correctly (so first fun1 increases gradually and fun2 remains the same, then the other way around, then fun1 decreases gradually and fun2 remains the same and then again the other way around). There should be something in the comsol help/manual about defining your own functions but if you do not manage please ask.
Does this help you?
Please login with a confirmed email address before reporting spam
I am trying to model a laser pulse which has a rectangular shape in 2D axisymmetric domain.
I have already modeled a Guassian shape with the equation: I0*exp(-2*(r/W)^2) as the boundary condition, where: I0: peak light intensity
r radial direction (i.e. the x axis) (Y axis is light falling on the surface)
W: beam width.
How do I express rectangular shape in mathematical expression using either step function or rectangular function?
Please suggest some ideas.
Thank you,
Ishan
Note that while COMSOL employees may participate in the discussion forum, COMSOL® software users who are on-subscription should submit their questions via the Support Center for a more comprehensive response from the Technical Support team.
Suggested Content
- FORUM Coil heating simulation
- FORUM surface laser
- BLOG Using the General Extrusion Coupling Operator in COMSOL: Dynamic Probe
- KNOWLEDGE BASE How does COMSOL Multiphysics handle complex-valued numbers and problems in the frequency domain?
- KNOWLEDGE BASE Solving electromagnetic heating models with high-frequency loads