Simulation of Gaussian Laser Beam

[LYS]

Dear all !

I am new to Elmer. My Elmer test case is a heat conduction (convection and radiation) case, and I have a some problem to simulate circular pulsed laser beam on CFRP material (see maillageGMSH_1Spot.gif)
The pulse is generated in my Elmer program by introducing a new variable "DiracFonction" in Solver 1. This variable is used in Boundary Condition 1 (see attached file case_19-10-10.sif) first to create the dirac function and after to the heat flux.

My problem is the error message (unlisted keyword : [diracfonction] in section [Boundary condition 1]) when I simulate (cf. Pb_Déclaration_Variable.jpg).

What am I doing wrong ?

Thanks a lot in advance, for Elmer and this helpboard.

Regards
Elisabeth

[raback]

Hi

The property "DiracFonction" that you introduce is not a variable that could be used in MATC expression. Why not make use fully algebraic expression for the heat flux. I mean, you can estimate Dirac delta by a narrow gaussian distribution, for example.

-Peter

[LYS]

Thank you very much Peter, now it's working!!

I have an other problem : I want to simulate a material with a Variable Heat Conductivity as a function of direction (x, y and z) and I don't now how can I do it.

Ithink about : Heat Conductivity(3) = Variable Coordinate(3)
Real
Heat conductivity 1 = 7
Heat conductivity 1 = 7
Heat conductivity 1 = 0.2
End
But it doesn't work
Have you an idea ??

Thanks a lot in advance

[raback]

Hi

So do you want a scalar or a tensor valued heat conductivity? The first one goes, for example, as

Code: Select all

Heat Conductivity = Variable Coordinate
  Real MATC "1.0+1.2*tx(0)+3.4*tx(1)+5.6*tx(2)",

Where the tx refers to the three components of the Coordinate in C conventin, i.e. indexes start from 0.

-Peter

[LYS]

Hi Peter,
I search to simulate the heat equation with conductivity kx, ky and kz which are respectively the conductivity in the direction vector x, y and z of the 3D cartesian system.
(i.e. vector kx =( kx * vector x) with kx = 7 for example ; vector ky =( ky * vector y) with ky = 3 for example; vector kz =( kz * vector z) with kz = 0.2 for example)

[Juha]

Hi,

if i understood correctly, this would read

Heat Conductivity(3) = 7 3 0.2

BR, Juha

[LYS]

Hi Juha,

Thank you very much for your answer, I tested, but it doesn't work

[Juha]

Hi,

how about giving some defintion of the "not working"?

Juha

[LYS]

Hi Juha,
Now it work with a space between Conductivity and (3)
Thank you !