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.

Time dependant heat flux at the boundary

Please login with a confirmed email address before reporting spam

Hello,

I am trying to model 3D tissue. I wish to apply a time-dependent heat flux for one boundary. I have introduced an interpolation function under global definition. Basically, the function is a sin wave,i.e., f(t)=sin(wt). I am getting a smooth plot of the sine wave. The problem is that, I am unable to figure out the way to implement that heat flux to behave like a sin wave.

Thanking You

Arka

11 Replies Last Post 5 nov. 2016, 17:34 UTC−4
Ivar KJELBERG COMSOL Multiphysics(r) fan, retired, former "Senior Expert" at CSEM SA (CH)

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 15 sept. 2012, 05:20 UTC−4
Hi

I do not really understand what you are missing, if you have a boundary heat flux condition and you write
P0[W/m^2]*sin(omega[1/s]*t) with a time dependent solver and omega some frequency value then you should see your system being influenced by the sinus excitation.

But do not forget that a time series solver is a transient type analysis, normally for a few periods, if you are interested in the steady state sinus exciatiaion response you should rather consider the frequency domain solvers (harmonic developments of the equations, see the Equation view

--
Good luck
Ivar
Hi I do not really understand what you are missing, if you have a boundary heat flux condition and you write P0[W/m^2]*sin(omega[1/s]*t) with a time dependent solver and omega some frequency value then you should see your system being influenced by the sinus excitation. But do not forget that a time series solver is a transient type analysis, normally for a few periods, if you are interested in the steady state sinus exciatiaion response you should rather consider the frequency domain solvers (harmonic developments of the equations, see the Equation view -- Good luck Ivar

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 15 sept. 2012, 05:53 UTC−4
Thank you Mr. Ivar

Actually I understand that the problem was with unit. I have considered omega a constant and comsol understand time,t in sec. But inside Sin function t should be [1/s], and that solves problem. Now it is working fine.

Kind Regards

Arka
Thank you Mr. Ivar Actually I understand that the problem was with unit. I have considered omega a constant and comsol understand time,t in sec. But inside Sin function t should be [1/s], and that solves problem. Now it is working fine. Kind Regards Arka

Ivar KJELBERG COMSOL Multiphysics(r) fan, retired, former "Senior Expert" at CSEM SA (CH)

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 15 sept. 2012, 16:32 UTC−4
Hi

It's always useful to check the units, as "t" is in seconds you can either write omega[1]*t[1/s] or omega[1/s]*t to get a unitless call for sinus(), but do not forget an orange unit "error" is a warning, COMSOl calculate anyhow, while red is an error and COMSOL will flag it and halt

--
Good luck
Ivar
Hi It's always useful to check the units, as "t" is in seconds you can either write omega[1]*t[1/s] or omega[1/s]*t to get a unitless call for sinus(), but do not forget an orange unit "error" is a warning, COMSOl calculate anyhow, while red is an error and COMSOL will flag it and halt -- Good luck Ivar

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 17 sept. 2012, 05:44 UTC−4
Hello Sir,

I have another query. Regarding the temperature profile? Since I am irradiating the boundary using a sin wave. So, I am assuming that the temperature profile at the irradiated boundary should behave a similar nature like a wave. But I am getting temperature curve at a point on the boundary straight lines. It is not following a sin nature. I am using time step of 0.1 for time marching procedure. But I think the temperature profile is not taking all the time steps it is skipping certain points. It is showing irregular straight lines like it has skipped few points. I am not sure what is the problem.

Thanking You

Arka
Hello Sir, I have another query. Regarding the temperature profile? Since I am irradiating the boundary using a sin wave. So, I am assuming that the temperature profile at the irradiated boundary should behave a similar nature like a wave. But I am getting temperature curve at a point on the boundary straight lines. It is not following a sin nature. I am using time step of 0.1 for time marching procedure. But I think the temperature profile is not taking all the time steps it is skipping certain points. It is showing irregular straight lines like it has skipped few points. I am not sure what is the problem. Thanking You Arka


Ivar KJELBERG COMSOL Multiphysics(r) fan, retired, former "Senior Expert" at CSEM SA (CH)

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 17 sept. 2012, 09:10 UTC−4
Hi

indeed it looks like the time stepping is being too coarse, the thing is that by default COMSOL assumes a diffusion type model decay, so when you have a periodic BC you must tweak a little the time stepping settings.

Use rather an intermediate or a strict time stepping (and not the default automatic for periodic BC conditions) then be sure you define enough time steps per periode, 5 as absolute mnimum 10 or more is better (it all depends on the time you ant to wait for the solution ...

Thne you can also consider a gen-alpha, in case of converging issues

--
Good luck
Ivar
Hi indeed it looks like the time stepping is being too coarse, the thing is that by default COMSOL assumes a diffusion type model decay, so when you have a periodic BC you must tweak a little the time stepping settings. Use rather an intermediate or a strict time stepping (and not the default automatic for periodic BC conditions) then be sure you define enough time steps per periode, 5 as absolute mnimum 10 or more is better (it all depends on the time you ant to wait for the solution ... Thne you can also consider a gen-alpha, in case of converging issues -- Good luck Ivar

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 18 févr. 2013, 14:48 UTC−5
Dear Arka;
I have the same problem as you. I want to define a SIN Temperature for the comsol but i can not define it very well. Would you please help me and tell me that how did you define a Time dependant heat for the software?

Best Regards
Dear Arka; I have the same problem as you. I want to define a SIN Temperature for the comsol but i can not define it very well. Would you please help me and tell me that how did you define a Time dependant heat for the software? Best Regards

Please login with a confirmed email address before reporting spam

Posted: 8 years ago 3 nov. 2016, 08:05 UTC−4
Hi,

I have a similar problem. I have to set my convective heat flux depending how much time as passed.
For example;
Time (s) 0-500 set h=10
Time 501-1000 set h=20
Time 1001-1500 set h=35
Time 1501-2000 set h=15

How do I covert this into a IF statement to insert into the "Convective heat flux" field ?

Thanks,
Hi, I have a similar problem. I have to set my convective heat flux depending how much time as passed. For example; Time (s) 0-500 set h=10 Time 501-1000 set h=20 Time 1001-1500 set h=35 Time 1501-2000 set h=15 How do I covert this into a IF statement to insert into the "Convective heat flux" field ? Thanks,

Please login with a confirmed email address before reporting spam

Posted: 8 years ago 3 nov. 2016, 08:32 UTC−4
Hi

Make a piecewise function, say pw1, and call it then pw1(t). See attachment.

Lasse
Hi Make a piecewise function, say pw1, and call it then pw1(t). See attachment. Lasse


Please login with a confirmed email address before reporting spam

Posted: 8 years ago 3 nov. 2016, 11:51 UTC−4
Thank you,
I have done the piecewise function, but how are you plotting that graph with h as a function of t?
When I plot the graph I don't get any labels at all.

Thank you,


Hi

Make a piecewise function, say pw1, and call it then pw1(t). See attachment.

Lasse


Thank you, I have done the piecewise function, but how are you plotting that graph with h as a function of t? When I plot the graph I don't get any labels at all. Thank you, [QUOTE] Hi Make a piecewise function, say pw1, and call it then pw1(t). See attachment. Lasse [/QUOTE]

Please login with a confirmed email address before reporting spam

Posted: 8 years ago 3 nov. 2016, 17:30 UTC−4
The function set-up window (my attachment) asks for the units of the argument and the function. You must fill in them properly.
The function set-up window (my attachment) asks for the units of the argument and the function. You must fill in them properly.

Please login with a confirmed email address before reporting spam

Posted: 8 years ago 5 nov. 2016, 17:34 UTC−4
Ok, I got it now, thanks for all the help.
Ok, I got it now, thanks for all the help.

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.