Hello(adsbygoogle = window.adsbygoogle || []).push({});

I'm a french student, I'm actually not sure this is the good place to ask my question but as it deals with the nuclear fission I try here... don't hesitate to tell me if there is a better forum... thx..

well, I'm trying to solve numerically the Langevin equation, initially for the brownian motion, it is also used in nuclear physics as you probably know, to describe the fission of nuclei...

the Langevin equation is :

[tex]\ddot{x}+\beta \dot{x} + \frac{1}{m}\frac{\partial U}{\partial x} = \Gamma(t)[/tex]

where [tex]\beta=\frac{\gamma}{m}[/tex] is a friction term (from a stocke's law) divided by the mass of the brownian particle. (so [tex]\beta[/tex] is [tex]s^{-1}[/tex])

[tex]U(x)[/tex] is an extern potential

and [tex]\Gamma(t)[/tex] is a stochastic force divided my the mass.

well, now I decide, as Kramers did, to tell that my potential U is represented by two parabolas as you can see on http://nicolas.aunai.free.fr/courb.htm [Broken]

I call B the height of the maximum, and [tex]x_b[/tex] the correspondant X-coordonnate (assuming that the minimum will be at x=0)

I decide that B will be my caracteristic energy for the problem, and [tex]x_b[/tex] my caracteristc length.

I want, in order to solve my equation numerically, to write it without dimension, so I introduce the following new variables :

[tex]Q = \frac{x}{x_b}[/tex]

[tex]\Pi = \frac{U}{B}[/tex]

[tex]t'=\omega t[/tex] where [tex]\omega[/tex] is the caracteristic pulsation of my parabolas

i.e. the equations of the parabolas are :

[tex]U_1(x)=\frac{1}{2}m\omega^2 x^2[/tex] for [tex]x<x_0[/tex]

[tex]U_2(x)=-\frac{1}{2}m\omega^2(x-x_b)^2 + B[/tex] for [tex]x>x_0[/tex]

assuming the continuity of the potential, the functions and the derivatives shoud be equal between them at [tex]x_0[/tex] the junction point. these conditions give us the height B which is equals to :

[tex]\frac{1}{4}m\omega^2 x_b^2[/tex]

well now we have to re-write the langevin equation, replacing by the new variables... and here began my problems....

[tex]\ddot{x}+\beta \dot{x} + \frac{1}{m}\frac{\partial U}{\partial x} = \Gamma(t)[/tex]

[tex]\dot{x}[/tex] means [tex]\frac{dx}{dt}[/tex], so since we have [tex]t'=\omega t[/tex], the derivation by t' gives us an [tex]\omega[/tex]

[tex]\omega^2 x_b\ddot{Q}+\beta \omega x_b \dot{Q} + \frac{B}{m}\frac{\partial \Pi}{\partial Q} \frac{\partial Q}{\partial x} = \Gamma(t)[/tex]

we know that [tex]\frac{\partial Q}{\partial x} = \frac{1}{x_b}[/tex]

so we have :

[tex]\omega^2 x_b\ddot{Q}+\beta \omega x_b \dot{Q} + \frac{B}{mx_b}\frac{\partial \Pi}{\partial Q}= \Gamma(t)[/tex]

and with the B value we have :

[tex]\omega^2 x_b\ddot{Q}+\beta \omega x_b \dot{Q} + \frac{\omega^2 x_b}{4}\frac{\partial \Pi}{\partial Q}= \Gamma(t)[/tex]

finally, dividing by [tex]\omega^2 x_b[/tex] we obtain :

[tex]\ddot{Q}+\frac{\beta}{\omega}\dot{Q} + \frac{1}{4}\frac{\partial \Pi}{\partial Q}= \frac{1}{\omega^2 x_b}\Gamma(t)[/tex]

which I think is almost the result that I'm looking for.. since the first term is without dimention, [tex]\frac{\beta}{\omega}[/tex] is without dimension too, and the third term too..

my problem is the stochastic term...

now I don't know how to do with it.... if you have an idea...

I now, that the autocorrelation function of the Langevin Force is :

[tex]<\Gamma(t)\Gamma(t')> = 2\beta T \delta(t_1-t_2)[/tex]

since [tex]t_i = \frac{t_i'}{\omega}[/tex], can I write the following equation :

[tex]<\Gamma(t)\Gamma(t')> = 2\beta T \delta(\frac{t_1'}{\omega}-\frac{t_2'}{\omega})[/tex]

and with [tex]\delta(a x)=\frac{1}{a}\delta(x)[/tex]

I find :

[tex]<\Gamma(t)\Gamma(t')> = 2\omega \beta T \delta(t_1'-t_2')[/tex]

so [tex]\Gamma(t)[/tex] should be proportionnal to [tex]\sqrt{(2\omega\beta T \delta(t_1'-t_2'))}[/tex]

in the nuclear system, k (the Boltzmann constant) is equal to 1, so we can also introduce the new variable :

[tex]\Theta=\frac{T}{B}[/tex]

so we can put it in our relation :

[tex]\Gamma(t)\ \ \ \alpha\ \ \ \sqrt{(2\omega\beta B\Theta\delta(t_1'-t_2'))}[/tex]

well... I'm not sure where I'm going with this... can someone help me ?

thanks !

**Physics Forums - The Fusion of Science and Community**

Join Physics Forums Today!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# Simulation of the fission as a stochastic process

Can you offer guidance or do you also need help?

Draft saved
Draft deleted

Loading...

Similar Threads for Simulation fission stochastic | Date |
---|---|

I Why Monte Carlo Simulations? What do they tell us? | Feb 22, 2017 |

I Where are the four muons in this picture | Jul 31, 2016 |

A How reliable are MC simulations when it comes to jets? | Jul 16, 2016 |

How to simulate protons on a computer program to observe them? | Jan 12, 2016 |

Simulation software for modeling a dusty fission reactor variant. | Mar 21, 2014 |

**Physics Forums - The Fusion of Science and Community**