Note: This discussion is about an older version of the COMSOL Multiphysics® software. The information provided may be out of date.

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.

How Can A Temperature Dependent Variable Changes Permanently?

Please login with a confirmed email address before reporting spam

Hello, I'm using COMSOL 4.3b for simulating a laser-material interaction, and I have predefined several parameters for specifying material's physical properties such as material's density, material's thermal conductivity, etc. with temperature dependence (interpolation function), which I run no problems at all.

However I've found out, according to the simulation result, these properties couldn't be changed permanently. Take the material density for an example*, the material's density changes with increased temperature as a result of absorption to laser induced energy, which is correct, however it changes back to original value when material cools down after laser beam moves away, which is physically incorrect as the material's composition has changed, which led a permanent density change, so the density value shouldn't be supposed to return back.

* rho v.s T is plotted as: farm8.staticflickr.com/7316/9580528024_451830a09c.jpg

Therefore may anyone could kindly offer any advice(s) for suggesting how I can specify a one-way directional relationship for making sure that a predefined physical property ONLY changes with INCREASING temperature please? In other words, reverse change is forbidden. (in this case, density reduction with decreasing temperature is not allowed)

Thanks!

P.S. I have thought about an alternative, which I can convert these physical properties to time-dependent values, based on a good assumption that could fix the model.

2 Replies Last Post 26.08.2013, 07:56 GMT-4

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 24.08.2013, 16:41 GMT-4
Hi,


Create a new time dependent "Domain ODEs and DAEs" equation to represent the variable material density (or any) and let us call it "max_rho_mat". If you did not know, this is available from "Mathematics > ODE and DAE interfaces > Domain ODEs and DAEs (dode)".

Once you have added the model, define the following coefficients:

f = ((nojac(rho_mat)-max_rho_mat)/timestep)*(rho_mat>max_rho_mat)
da = 1
ea = 0
under the sub-node "initial values" give the initial value of your material density, don't forget this aspect
rho_mat is your material density that is a function of temperature.

This will remember the largest value of the material density ever reached at any location at any give time (because this is what you asked for). So just use max_rho_mat instead of rho_mat in your heat transfer equation.


Suresh
Hi, Create a new time dependent "Domain ODEs and DAEs" equation to represent the variable material density (or any) and let us call it "max_rho_mat". If you did not know, this is available from "Mathematics > ODE and DAE interfaces > Domain ODEs and DAEs (dode)". Once you have added the model, define the following coefficients: f = ((nojac(rho_mat)-max_rho_mat)/timestep)*(rho_mat>max_rho_mat) da = 1 ea = 0 under the sub-node "initial values" give the initial value of your material density, don't forget this aspect rho_mat is your material density that is a function of temperature. This will remember the largest value of the material density ever reached at any location at any give time (because this is what you asked for). So just use max_rho_mat instead of rho_mat in your heat transfer equation. Suresh

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 26.08.2013, 07:56 GMT-4
Hi Suresh,

Thanks so much for your suggestion! I will have a try later today.


Best regards,

Liwei



Hi,


Create a new time dependent "Domain ODEs and DAEs" equation to represent the variable material density (or any) and let us call it "max_rho_mat". If you did not know, this is available from "Mathematics > ODE and DAE interfaces > Domain ODEs and DAEs (dode)".

Once you have added the model, define the following coefficients:

f = ((nojac(rho_mat)-max_rho_mat)/timestep)*(rho_mat>max_rho_mat)
da = 1
ea = 0
under the sub-node "initial values" give the initial value of your material density, don't forget this aspect
rho_mat is your material density that is a function of temperature.

This will remember the largest value of the material density ever reached at any location at any give time (because this is what you asked for). So just use max_rho_mat instead of rho_mat in your heat transfer equation.


Suresh


Hi Suresh, Thanks so much for your suggestion! I will have a try later today. Best regards, Liwei [QUOTE] Hi, Create a new time dependent "Domain ODEs and DAEs" equation to represent the variable material density (or any) and let us call it "max_rho_mat". If you did not know, this is available from "Mathematics > ODE and DAE interfaces > Domain ODEs and DAEs (dode)". Once you have added the model, define the following coefficients: f = ((nojac(rho_mat)-max_rho_mat)/timestep)*(rho_mat>max_rho_mat) da = 1 ea = 0 under the sub-node "initial values" give the initial value of your material density, don't forget this aspect rho_mat is your material density that is a function of temperature. This will remember the largest value of the material density ever reached at any location at any give time (because this is what you asked for). So just use max_rho_mat instead of rho_mat in your heat transfer equation. Suresh [/QUOTE]

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.