Open access

New Methods in Doppler Broadening Function Calculation

Written By

Daniel Artur P. Palma, Alessandro da C. Gonçalves, Aquilino Senra Martinez and Amir Zacarias Mesquita

Submitted: 16 May 2012 Published: 06 February 2013

DOI: 10.5772/52464

From the Edited Volume

Current Research in Nuclear Reactor Technology in Brazil and Worldwide

Edited by Amir Zacarias Mesquita

Chapter metrics overview

2,105 Chapter Downloads

View Full Metrics

1. Introduction

In all nuclear reactors some neutrons can be absorbed in the resonance region and, in the design of these reactors, an accurate treatment of the resonant absorptions is essential. Apart from that, the resonant absorption varies with fuel temperature, due to the Doppler broadening of the resonances (Stacey, 2001). The thermal agitation movement of the reactor core is adequately represented in microscopic cross-section of the neutron-core interaction through the Doppler Broadening function. This function is calculated numerically in modern systems for the calculation of macro-group constants, necessary to determine the power distribution in a nuclear reactor. This function has also been used for the approximate calculations of the resonance integrals in heterogeneous fuel cells (Campos and Martinez, 1989). It can also be applied to the calculation of self-shielding factors to correct the measurements of the microscopic cross-sections through the activation technique (Shcherbakov and Harada, 2002). In these types of application we can point out the need to develop precise analytical approximations for the Doppler broadening function to be used in the codes that calculates the values of this function. Tables generated from such codes are not convenient for some applications and experimental data processing.

This chapter will present a brief retrospective look at the calculation methodologies for the Doppler broadening function as well as the recent advances in the development of simple and precise analytical expressions based on the approximations of Beth-Plackzec according to the formalism of Briet-Wigner.

Advertisement

2. The Doppler broadening function

Let us consider a medium with a temperature where the target nuclei are in thermal movement. In a state of thermal equilibrium for a temperature T, the velocities are distributed according to Maxwell-Boltzmann distribution (Duderstadt and Hamilton, 1976),

f(V)=N(M2πkT)32eMV22kT,E1

where N is the total number of nucleus, M is the mass of the nucleus and k is Boltzmann’s constant.

Considering the neutrons as an ideal gas in thermal equilibrium, it is possible to write the average cross-section for neutron-nucleus interaction taking into consideration the movement of the neutrons and of the nucleus as:

σ(v,T)=1vNd3V(|vV|)σ(|vV|)f(V),E2

where f(V) is the distribution function of Maxwell-Boltzmann as given by equation (1) and V=VΩ^ is the velocity of the target nuclei. Denoting vr=vV the relative velocity between the movement of the neutron and the movement of the target nucleus and considering the isotropic case, that is, with no privileged direction, it is possible to separate the integration contained in equation (2) in the double integral:

σ(v,T)=1vN0dVV2f(V)4πvrσ(vr)dΩ^.E3

It is possible to see clearly in equation (3) that the cross-section depends of the relative velocity between the neutrons and the target nuclei. As the nuclei are in thermal movement, the relative velocity can increase or decrease. This difference between relative velocities rises to the Doppler deviation effect in cross-section behaviour. After integrating equation (3) in relation to the azimuthal angle (ϕ) the average cross-section for neutron-nucleus interaction can be written thus:

σ(v,T)=2πvN0dVV2f(V)0πvrσ(vr)sinθdθ.E4

Denoting μ=cosθ so that dμ=sinθdθ, equation (4) takes the form of:

σ(v,T)=2πvN0dVV2f(V)11vrσ(vr)dμ.E5

From the definition of the relative velocity one has the relation,

vr2=v2+V22vVμ,E6

and, as a result,

dμ=vrdvrvV.E7

With the aid of a simple substitution, using relations (6) and (27), equation (5) is thus written as:

σ(v,T)=2πv2N0dVVf(V)|vV||v+V|vr2σ(vr)dvr.E8

In equation (8), the limits of integration are always positive due to the presence of the module. As a result, one should separate the integral found in equation (8) into two separate integrals, as follows,

σ(v,T)=2πv2N[0vdVVf(V)vVv+Vvr2σ(vr)dvr+vdVVf(V)Vvv+Vvr2σ(vr)dvr].E9

It is possible to modify the limits of integration for equation (9) taking into account that the mass of the target nucleus is much larger than the mass of the incident neutron. In terms of relative velocity, equation (9) can be written as:

σ(v,T)=2πv2N[vvrv+vrvr2σ(vr)dvr0vdVVf(V)+vrvv+vrvr2σ(vr)dvrvdVVf(V)].E10

In replacing the expression of the Boltzmann distribution function, equation (1), in equation (10) one has:

σ(v,T)=2β3v2π[vvrv+vrvr2σ(vr)dvr0vdVVeβ2V2+vrvv+vrvr2σ(vr)dvrvdVVeβ2V2],E11

where it was defined β2M2kT. Introducing the variables for reduced velocities ϖr=βvr and ϖ=βv, equation (11) is written by:

σ(v,T)=2β2ϖ2π×[0ϖ/βϖr2σ(ϖr)dϖr(ϖϖr)/β(ϖ+ϖr)/βdVVeβ2V2+ϖ/βϖr2σ(ϖr)dϖr(ϖrϖ)/β(ϖ+ϖr)/βdVVeβ2V2].E12

Integrating equation (12) in relation to V one gets to the expression:

σ(v,T)=1ϖ2π0ϖr2σ(ϖr)[e(ϖϖr)2e(ϖ+ϖr)2]dϖr.E13

For resonances (that is, for the energy levels of the composed nucleus) it is possible to describe the energy dependence of the absorption cross-section by a simple formula, valid for T=0K, known as Breit-Wigner formula for resonant capture, expressed in function of the energy of the centre-of-mass by,

σγ(ECM)=σ0ΓγΓ(E0ECM)1/211+4Γ2(ECME0)2,E14

where E0 is the energy where the resonance occurs and ECM is the energy of the centre-of-mass of the neutron–nucleus system. Apart from that, we find in equation (14) the term σ0, that is the value of the total cross-section σtotal(E) in resonance energy E0 that can be written in terms of the reduced wavelength ƛ0 by:

σ0=4πƛ02ΓnΓg=2.608×106(A+1)2A2E(eV)ΓnΓg,E15

where the statistical spin factor g is given by the expression:

g=2J+12(2I+1),E16

where I is the nuclear spin and J is the total spin (Bell and Glasstone, 1970).

In replacing the expression (14) in equation (13) one finds an exact expression for the average cross-section, valid for any temperature:

σγ(v,T)=σ0ΓγΓβ2πv2×0dvr(E0ECM)1/2vr21+4Γ2(ECME0)2[eβ2(vvr)2eβ2(v+vr)2],E17

In a system with two bodies it is possible to write the kinetic energy in the centre-of-mass system, by

ECM=MRvr22,E18

where MR=mMm+M is the reduced mass of the system.

For the problem at hand, of a neutron that is incident in a thermal equilibrium system with a temperature T, it is a good approximation to assume that vvr. Thus, the ratio between the kinetic energy of the incident neutron and the kinetic energy of the centre-of-mass system is thus written

ECME0=A+1A,E19

where A is the atomic mass of the target core. Resulting:

y=2Γ(ECME0)(a)x=2Γ(EE0),(b)E20

and denoting β2=12vth2 one finally obtains the expression for the cross-section of radioactive capture near any isolated resonance with an energy peak E0, as written by:

σγ(E,T)=σ0ΓγΓ(E0E)1/2Ψ(x,ξ),E21

where

Ψ(x,ξ)=ξ22E/Γ+dy1+y2[exp((vvr)22v2th)exp((v+vr)22v2th)],E22

where vr is the module for relative neutron-nucleus velocity, v is the module for neutron velocity, and

ξΓΓD.E23

The Doppler width for resonance ΓD is expressed by:

ΓD=(4E0kT/A)1/2.E24

All the other nuclear parameters listed below are well established in the literature,

  • A= mass number;

  • T= absolute temperature;

  • E= energy of incident neutron;

  • ECM= energy of centre-of-mass;

  • E0= energy where the resonance occurs;

  • Γ= total width of the resonance as measured in the lab coordinates;

  • ΓD=(4E0kT/A)1/2= Doppler width of resonance;

  • v= neutron velocity module;

  • vr=|vV|= module of the relative velocity between neutron movement and nucleus movement;

  • vth=2kTM= module of the velocity for each target nucleus.

Advertisement

3. The Bethe and Placzek approximations

The expression proposed by Bethe and Placzek for the Doppler broadening function ψ(x,ξ) is obtained from some approximations, as follows:

  1. one neglects the second exponential in equation (22), given that it decreases exponentially and is negligible in relation to first integral in equation (22) given that (v+vr)2>>(vvr)2.

  2. it is a good approximation to extend the lower limit for integration down to in equation (22), given that the ratio between the energy of neutron incidence and the practical width is big.

  3. being ECM the energy of the system in the centre-of-mass system and E the energy of the incident neutron, the following relation is always met:

ECM1/2=E1/2(1+ECMEE)1/2=E1/2(1+η)1/2,E25

where it was denoted that η=ECMEE. Equation (25) can be expanded in a Taylor series and, to the first order, is written by

ECM1/2=E1/2(1+η2η24+...)E1/2(1+ECME2E)E26

In terms of the masses and velocities, equation (26) is written as follows:

(MRvr22)1/2=(mv22)1/2(1+MRvr22mv22mv2),E27

where MR is the reduced mass of the system. For heavy nucleus MRm and equation (27) can be written as:

vr=vr2+v22v,E28

so that,

vvr=vvr2+v22v=v2vr22v.E29

In replacing approximation equation (29) in the remaining exponential of equation (22) one finally obtains the Doppler broadening function that will be approached in this chapter,

ψ(x,ξ)ξ2π+dy1+y2exp[ξ24(xy)2].E30

The approximations made in this section apply in almost all the practical cases, and are not applicable only in situations of low resonance energies (E<1eV) and very high temperatures.

Advertisement

4. Properties of the Doppler broadening function ψ(x,ξ)

The function ψ(x,ξ) as proposed by the approximation of Bethe and Placzek has an even parity, is strictly positive and undergoes a broadening as that variable ξ diminishes, that is, varies inversely with the absolute temperature of the medium. For low temperatures, that is, when temperature in the medium tend to zero, the Doppler broadening function can be represented as shown below:

limT0ψ(x,ξ)=ξ2πlimT0+dy1+y2exp[ξ24(xy)2]=11+x2E31

Equation (31) is known as an asymptotic approximation of the Doppler broadening function. For high temperatures, that is, when the temperature of the medium tends to infinite, the Doppler broadening function can be represented through the Gaussian Function, given that:

limTψ(x,ξ)=ξ2πlimT+dy1+y2exp[ξ24(xy)2]=ξ2πexp(ξ24x2).E32

Figure 1.

The Doppler broadening function for ξ=0.05, 0.15 and 0.5.

The area over the curve of the Doppler Broadening function is written as below and, as it consists of separable and known integers it is possible to write:

+ψ(x,ξ)dx=ξ2π+dy1+y2+eξ24(xy)2dx=ξ2π(π)(π2ξ)=π.E33

From equation (33) once concludes that the area over the curve of the Doppler Broadening function is constant for the intervals of temperature and energy of interest in thermal reactors. This property is valid even for broadened resonances as shown in Figure 1, considering the different values of variable ξ.

Advertisement

5. Analytical approximations for the Doppler broadening function

This section describes the main approximation methods for the Doppler broadening function, according to the approximation of Bethe and Placzek, equation (30).

5.1. Asymptotic expansion

A practical choice to calculate the Doppler broadening function is its asymptotic expression resulting from the expansion of the term 11+y2 in equation (30) in a Taylor series around y=x.

11+y2=11+x22x(1+x2)2(yx)+1+3x2(1+x2)2(yx)24x(1+x2)(1+x2)4(yx)3+E34

In replacing equation (34) in equation (30) and integrating term by term, one obtains the following the asymptotic expansion:

ψ(x,ξ)=11+x2{1+2ξ2(3x21)(1+x2)2+12ξ4(5x410x21)(1+x2)4+...}.E35

Despite equation (35) being valid only for |x.ξ|>6, it is quite useful to determine the behaviour of the Doppler Broadening function in specific conditions. For high values of x, it is possible to observe that function ψ(x,ξ) presents the following asymptotic form:

ψ(x,ξ)11+x2 E36

5.2. Method of Beynon and Grant

Beynon and Grant (Beynon and Grant, 1963) proposed a calculation method for the Doppler broadening function that consists of expanding the exponential part of the integrand of the Doppler broadening function ψ(x,ξ) in the Chebyshev polynomials and integrate, term by term, the resulting expression, which allows writing:

ψ(a,b)=1a{πcos(ab)[1- E2(a)]ea2+J(a,b)}e-14b2,E37

where a=12ξ and b=ξx and still,

J(a,b)=1a{12!(ab)214!(ab)4+16!(ab)6}+12a3{14!(ab)416!(ab)6}++1πa2n+1Γ(2n+12){1(2n+2)!(ab)2(n+1)}+E38

and,

E2(a)=2π0aey2 dy.E39

For values where the condition |x.ξ|>6 was met (Beynon and Grant, 1963) is recommended the use of the asymptotic expression of the function ψ(x,ξ), equation (36). It should be pointed that the results obtained by this method have become a reference in several works on the Doppler broadening function.

5.3. Method of Campos and Martinez

The core idea of the method proposed by (Campos and Martinez, 1987) is to transform the Doppler broadening function from its integral form into a differential partial equation subjected to the initial conditions. Differentiating equation (30) in relation to x one obtains:

ψ(x,ξ)x=ξ22×{ξx2πdy1+y2exp[ξ24(xy)2]+ξ2πydy1+y2exp[ξ24(xy)2]}.E40

Acknowledging in equation (40) the very Doppler broadening function and the term of interference as defined by the integral:

χ(x,ξ)=ξπ+ydy1+y2exp[ξ24(xy)2],E41

it is possible to write:

ψ(x,ξ)x=ξ22[xψ(x,ξ)+χ(x,ξ)2].E42

Deriving equation (42) again in relation to x, after expliciting function χ(x,ξ) in the same equation, one has:

2ψ(x,ξ)x2+xξ2ψ(x,ξ)x+ξ24(x2ξ2+2)ψ(x,ξ)=ξ58πy2dy1+y2exp[ξ24(xy)2].E43

The right side of equation (43) can be written in another way, given that y21+y2=111+y2,

y2dy1+y2exp[ξ24(xy)2]=2πξ[1ψ(x,ξ)].E44

In replacing the result obtained in (44) in equation (43) one obtains the differential equation where Campos and Martinez based themselves to obtain an analytical approximation for the broadening function ψ(x,ξ):

4ξ22ψ(x,ξ)x2+4xψ(x,ξ)x+(ξ2x2+ξ2+2)ψ(x,ξ)=ξ2,E45

subjected to the initial conditions:

ψ(x,ξ)|x=0ψ0=ξπ2exp(ξ24)[1erf(ξ2)](a)ψ(x,ξ)x|x=0=0.(b)E46

Admitting that function ψ(x,ξ) may be expanded in series,

ψ(x,ξ)=n=0cn(ξ) xnE47

and in replacing-se equation (47) in the differential equation as given by equation (45), one obtains after some algebraic manipulation the following polynomial equation:

[8ξ2c2+(ξ2+2)c0]+[24ξ2c3+(ξ2+6)c1]x++n=2[4ξ2(n+2)(n+1)cn+2+(4n+ξ2+2)cn+ξ2cn2]xn=ξ2,E48

where:

c0=ψ0,

c1=ξ28[ξ2(ξ2+2)ψ0],

and all the other terms are calculated from the following relation of recurrence:

cn+1=ξ24(4n+ξ2+2)cn+ξ2cn1(n+2)(n+1)

The representation in series for the Doppler broadening function, as given by equation (47), is valid only for |x.ξ|<6. For the cases where |x.ξ|>6, (Campos and Martinez, 1987) used the asymptotic form as given by equation (35), as well as proposed by Beynon and Grant.

5.4. Four order method of Padé

The Padé approximation is one of the most frequently used approximations for the calculation of the Doppler broadening function and its applications and can efficiently represent functions, through a rational approximation, that is, a ratio between polynomials. For the four-order Padé approximation (Keshavamurthy& Harish, 1993) they proposed the following polynomial ratio:

ψ(x,ξ)=ha0+a2(hx)2+a4(hx)4+a6(hx)6b0+b2(hx)2+b4(hx)4+b6(hx)6+b8(hx)8,E49

whose coefficients are given in Tables 1 and 2.

p0=π q1=π(9π+28)2(6π229π+32)
p1=15π2+88π1282(6π229π+32) q2=36π2195π+2566(6π229π+32)
p2=π(33π104)6(6π229π+32) q3=π(33π+104)6(6π229π+32)
p3=9π2+69π1283(6π229π+32) q4=9π269π+1283(6π229π+32)

Table 1.

Coefficients p and q of the four-order Padé Approximation

h=ξ2
a0=(p0+p1hp2h2p3h3)(1q1hq2h2+q3h3+q4h4)
a2=(p2+3p3h)(1q1hq2h2+q3h3+q4h4)+(p0+p1hp2h2p3h3)(q23q3h6q4h2)+(p1+2p2h+3p3h2)(q1+2q2h3q3h24q4h3)
a4=q4(p0+p1hp2h2p3h3)+(p2+3p3h)(q23q3h6q4h2)p3(q1+2q2h3q3h24q4h3)+(p1+2p2h+3p3h2)(q3+4qh4)
a6=q4(p2+3p3h)p3(q3+4qh4)
b0=(1q1hq2h2+q3h3+q4h4)2
b2=2(1q1hq2h2+q3h3+q4h4)(q23q3h6q4h2)+(q1+2q2h3q3h24q4h3)2
b4=(q23q3h6q4h2)2+2q4(1q1hq2h2+q3h3+q4h4)+2(q1+2q2h3q3h24q4h3)(q3+4qh4)
b6=2q4(q23q3h6q4h2)+(q3+4qh4)2
b8=q42

Table 2.

Coefficients h, a and b of the four-order Padé Approximation

From the coefficients of Tables 1 and 2, and of equation (49), one obtains in the end the following analytical approximation for function ψ(x,ξ), according to the four-order Padé approximation:

ψ(x,ξ)=η(x,ξ)ω(x,ξ),E50

where η(x,ξ) and ω(x,ξ) are the following polynomials:

η(x,ξ)=2ξ(7,0898154041022+1,1467508441023ξ+8,3997250591022ξ2+3,6222070531022ξ3+9,9577517401021ξ4+1,7490672581021ξ5+1,8351652131020ξ6+8,9400726991018ξ72,5397366571021ξ2x2+2,0694839911021ξ3x2+3,9723935481021ξ4x2+1,9193195601021ξ5x2+3,6703304261020ξ6x2+2,6820218081019ξ7x2+1,0487480261019ξ4x4+1,7025230081020ξ5x4+1,8351652091020ξ6x4+2,6820218061019ξ7x4+8,9400726881018ξ7x6),E51

and

η(x,ξ)=(3,4906429251023ξ+3,4649993811023ξ2+2,0501509911023ξ3+7,9337711181022ξ4+3,6703304271020ξ7x6+1,7880145391019ξ8x8+3,6703304261020ξ7+3,5338948061021ξ6+1,7880145411019ξ8+2,0628594601022ξ5+3,4268437961022ξ2x2+5,5866136301022ξ4x2+2,6497033231022ξ5x2+6,6135126251022ξ3x2+1,1010991291021ξ7x2+7,3010133531021ξ6x2+3,5907744131021ξ4x4+1,1010991251021ξ7x4+5,8684385811021ξ5x4+4,0003422611021ξ6x4+7,1520581561019ξ8x2+2,3322373051020ξ6x6+1,0728087211020ξ8x4+7,1520581521019ξ8x6.E52

5.5. Frobenius method

In this method the homogeneous part of the differential equation that rules the Doppler broadening function, equation (45), is solved using the Frobenius Method (Palma et. al., 2005) that consists fundamentally of seeking a solution of the differential equation in the form of series around the point x=x0, with a free parameter, that is, as follows:

ψ(x,ξ)=xsn=0cn(ξ) xn=n=0cn(ξ) xn+s,E53

with c00 and where s is the parameter that grants the method flexibility.

Deriving equation (53) and replacing it in the homogeneous equation associated to equation (45) one obtains, after grouping the similar terms:

n=0cn(n+s)(n+s1) xn+s2+n=0cnξ2[(n+s)+ξ2+24] xn+s+ξ24n=0cn xn+s2=0.E54

The initial equation of the problem, obtained when n=0, remembering that c00 is

c0s(s1)=0.E55

From equation (55), as c00, one obtains that s=0 or s=1. Using first s=0 and c00 one obtains the following relations of recurrence:

cn=ξ2(4n+ξ26)4n(n+1)cn2, valid forn=2orn=3E56
cn=ξ2[cn2(4n+ξ26)+cn4ξ2]4n(n+1),valid for  n4.         E57

Considering the case where s=1, one obtains the other series linearly independent with the first term, not null, denoted by c˜0:

cn=ξ2(4n+ξ22)4n(n+1)cn2, valid for  n4E58
cn=ξ2[cn2(4n+ξ22)+cn4ξ2]4n(n+1), valid for  n4E59

With this the homogeneous solution assumes the following form:

ψh(x,ξ)=(c0+c2x2+c4x4+...)+(c˜0x+c˜2x3+c˜4x5+...),E60

where the coefficients are all known from the relations of recurrence, equations (56) to (59). In writing function ψh(x,ξ) as:

ψh(x,ξ)=exp(ξ2x24)n=0Anxn=(1ξ2x24+ξ4x432+...)n=0Anxn=A0+A1x+(A2ξ24A0)x2+(A3ξ24A1)x3+(A4ξ24A2+ξ432A0)x4+(A5ξ24A3+ξ432A1)x5+(A6ξ24A4+ξ432A2ξ6384A0)+(A7ξ24A5+ξ432A3ξ6384A1)+.E61

it is possible to determine all the coefficients An equalling, term by term, equations (60) and (61) so to write:

ψh(x,ξ)=k1exp(ξ2x24)[112(ξ2x2)2+124(ξ2x2)4+1720(ξ2x2)6+...]+k2exp(ξ2x24)[ξ2x216(ξ2x2)3+1120(ξ2x2)5+15040(ξ2x2)7+...],E62

Acknowledging the expansion of the cosine and sine functions, one obtains an analytical form to solve the homogeneous part of the differential equations that rule the Doppler broadening function:

ψh(x,ξ)=exp(ξ2x24)[k1cos(ξ2x2)+k2sin(ξ2x2)].E63

In order to obtain the particular solutions of equation (45), and consequently its general solution, it is possible to apply the method of parameter variation from the linearly independent solutions:

ψ1(x,ξ)=exp(ξ2x24)cos(ξ2x2)E64
ψ2(x,ξ)=exp(ξ2x24)sin(ξ2x2),E65

Supposing a solution thus,

ψp(x,ξ)=u1ψ1(x,ξ)+u2ψ2(x,ξ)E66

where functions u1(x) and u2(x) are determined after the imposition of the initial conditions expressed by equations (46a) and (46b) and of the imposition of the nullity of the expression:

u'1(x)ψ1(x,ξ)+u'2(x)ψ2(x,ξ)=0,E67

That, along with the condition,

u'1(x)ψ'1(x,ξ)+u'2(x)ψ'2(x,ξ)=ξ44,E68

Which results from the very equation (45), form a linear system whose solution is given by the equations:

u1'(x)=ξ22exp(ξ2x24)sin(ξ2x2)u1(x)=ξ220xdx'exp(ξ2x'24)sin(ξ2x'2)E69
u2'(x)=ξ22exp(ξ2x24)cos(ξ2x2)u2(x)=ξ220xdx'exp(ξ2x'24)cos(ξ2x'2).E70

Integrating equations (69) and (70),

u1(x)=ξπ4exp(ξ24)[erf(iξxξ2)erf(iξx+ξ2)+2erf(ξ2)]E71
u2(x)=iξπ4exp(ξ24)[erf(iξxξ2)+erf(iξx+ξ2)],E72

it is possible to write the solution particular of equation (45) as follows:

ψp(x,ξ)=iξπ4sin(ξ2x2)exp[ξ24(x21)][erf(iξxξ2)+erf(iξx+ξ2)]+ξπ4cos(ξ2x2)exp[ξ24(x21)][erf(iξxξ2)erf(iξx+ξ2)+2erf(ξ2)].E73

As the general solution of differential equation (45) is the sum of the solution of the homogeneous and particular equations, the initial conditions are imposed, as expressed by equations (46a) and (46b), to determine the constants:

k1=ξπ2exp(ξ24)[1erf(ξ2)]E74
k2=0.E75

Finally, according to the Frobenius Method, the Doppler Broadening function can be written thus:

ψ(x,ξ)=ξπ2exp[14ξ2(x21)]cos(ξ2x2)×{1+Reϕ(x,ξ)+tan(ξ2x2)Imϕ(x,ξ)},E76

where ϕ(x,ξ)=erf(iξxξ2).

5.6. Fourier transform method

In doing the transformation of variables u=ξ2(xy) in the full representation of the Doppler broadening function, equation (30), one obtains the expression

ψ(ξ,x)=1π+eu2du1+(x2uξ)2 ,E77

that can be mathematically interpreted as the convolution of the Lorentzian function with a gaussian function, as exemplified by the equation below:

ψ(ξ,x)=f*g+g(u) f(x-u) du ,E78

where f(x-u)=11+(x2uξ)2 is the lorentzian function and g(u)=1πeu2 the gaussian function. Function f(x-u) admits a full representation through the Fourier cosine transform (Polyanin and Manzhirov, 1998), as being

f(x-u)=0ewcos[(x2uξ)w] dw.E79

In replacing-se equation (79) in the integer of convolution, as given by equation (78), applying the properties of the integrals of convolution one gets to the following expression:

ψ(ξ,x)=f*g 0e-w+g(u) cos[(x-2uξ )w] du dw =0e-w I(w) dw ,E80

where,

I(w)  1π+eu2 cos[(x-2uξ) w] du=1πcos(xw)+eu2cos[(2uξ) w] du=e-wξ22cos(xw).E81

In replacing equation (81) in the equation (80), one obtains a new full representation of the Doppler broadening function, interpreted as a Fourier cosine transform (Gonçalves et. al., 2008):

ψ(ξ,x) = 0e-w2ξ2wcos(wx) dw= 12[ 0ew2ξ22wadw +  0ew2ξ22wbdw ],E82

where a(1ix)2 and (1+ix)2.

The integrals on the right side of equations (3.25) and (3.26) are known as complementary error functions, in which case one can conclude that:

0ew2ξ22 wadw =ξπ2e(xi-1)24ξ2erfc (ξiξx2)E83
0ew2ξ22 wbdw =ξπ2e(xi+1)24ξ2erfc (ξ+iξx2).E84

In replacing equations (83) and (84) in equation (82) it is possible to write the following expression for the Doppler broadening function:

ψ(x,ξ)= ξπ4e(xi1)2ξ42[1+erf (iξxξ2)]+ξπ4e(xi+1)2ξ42[1erf (iξx+ξ2)].E85

With some algebraic manipulation it is easy to prove that the Fourier transform method and the Frobenius method, equations (85) and (76) respectively, provide identical results.

5.7. Fourier series method

From the representation of the Doppler broadening function in a Fourier cosine transform, equation (82), it is possible to write

ψ(ξ,x) = 0e-w2ξ2wcos(wx) dw=0G(w)ewcos(wx) dw,E86

where function G(w)=e-w2ξ2 is even and can be expanded into a Fourier series in cosines:

G(w)=a02+n=1ancos(nπwL),E87

where

a0=ξπLerf(Lξ)E88
an=ξπ2Le(nπξ2L)2[erf (2L+nπξ2i2ξL)+erf (2Lnπξ2i2ξL)].E89

In replacing equation (87) in equation (86) and integrand, it is possible to write the following expression for the Doppler broadening function in the form of Fourier series:

ψ(x,ξ)=ξπ2L(1+x2)erf(Lξ)+ξπLn=1Fn(x,ξ,L) Re [Z(ξ,L)],E90

where

Fn(x,ξ,L)[(nπ)2+L2(1+x2)]e (n π ξ2 L)2L2(1+x2)2+(nπ)2(22x2+(nπ/L)2),E91
Z(n,ξ,L)=erf (nπξ2i+2L22ξL).E92

5.8. Representation of function ψ(x,ξ) using Salzer expansions

Although the formulations obtained for function ψ(x,ξ) from the Frobenius method, Fourier transform and Fourier series methods only contain functions that are well-known in literature, it can be inconvenient to work with error functions that contain an imaginary argument. One of the ways to overcome this situation is to calculate the real and imaginary parts of function ϕ(x,ξ) using the expansions proposed by Salzer (Palma and Martinez, 2009)

ϕ(x,ξ)=erf(iξxξ2)=Reϕ(x,ξ)+Imϕ(x,ξ)i,E93

where:

Reϕ(x,ξ)erf(ξ2)+exp(ξ24)×{1πξ[cos(ξ2x2)1]+2πn=1nmaxexp(n2/4)n2+ξ2 fn(x,ξ)}E94
Imϕ(x,ξ)exp(ξ24){1πξsin(ξ2x2)+2πn=1nmaxexp(n2/4)n2+ξ2 gn(x,ξ)},E95

where auxiliary functions fn(x,ξ) and gn(x,ξ) are written by:

fn(x,ξ)=ξ+ξcosh(nξx2)cos(ξ2x2)nsinh(nξx2)sin(ξ2x2),E96
gn(x,ξ)=ξcosh(nξx2)sin(ξ2x2)+nsinh(nξx2)cos(ξ2x2).E97

5.9. The Mamedov method

Mamedov (Mamedov, 2009) put forward an analytical formulation to calculate function ψ(x,ξ), based on its representation in the form of a Fourier transform, equation (82). Using the expansions in series of the exponential and cosine functions,

ex=limNi=1N(1)ixii!,E98
cosx=limNi=0N(1)ix2i(2i)!,E99

and the well-known binomial expansion

(x±y)n=limNm=0N(±1)mFm(n)xnmym,E100

Mamedov proposed the following expressions for the Doppler broadening function:

for ξ>1 and x>1

ψ(x,ξ)=ξ2πlimLi=0LFi(1){1(x2+1)i+1j=0iFj(i)[1+(1)j]×22ixjξ2ij+1γ(2ij+12,ξ2(x2+1)24)+(x2+1)ilimMk=0MFk(1i)×[1(1)k]xkξ2i+k+122i+2Γ(2i+k+12,ξ2(x2+1)24)},E101

for ξ1 and x15

ψ(x,ξ)=ξ2πlimLi=0Lj=0M(1)i+jx2jξ2j+i+12i!(2j!)Γ(i+2j+12),E102

for ξ1 and x=0

ψ(x,ξ)=ξ2πlimLi=0L(1)iξi+12i!Γ(i+12),E103

where Γ(x,ξ), γ(x,ξ) and Γ(x) are the well-know incomplete Gamma functions and Fm(n) are binomials coefficients defined by:

Fm(n)={n(n1)...(nm+1)m!, for integer n(1)mΓ(mn)m!Γ(n), for non-integer nE104
Advertisement

6. Numerical calculation of function ψ(x,ξ)

The numerical calculation of the Doppler broadening function consists of calculating a defined integral. There are many methods in the literature for this calculation, but in this chapter we will describe a numerical reference method based on the Gauss-Legendre quadrature. In basic terms, the Gauss-Legendre quadrature method consists of approximating a defined integer through the following expression:

11f(ba2η+b+a2)dxba2i=1Nwif(ba2ηi+b+a2),E105

where N is the order of the quadrature, ηiis the point of the quadrature and wi the weight corresponding to the point of quadrature. The points of the Gauss-Legendre quadrature are the roots of the polynomials of Legendre (Arfken, 1985) in the interval [1,1], as generated from the Rodrigues‘ formula,

Pn(x)=12nn!dndxn{(x21)n}.E106

for an isotope at a given temperature, that is, for a fixed value for variable ξ, the function ψ(x,ξ) decreases rapidly and a very high value is not necessary for what we will consider our numerical infinite. This fact can be evidenced at Figure 1.

For that an adequate numerical infinite (x=5000) was considered, as well as a high-order quadrature (N=15), whose points of Legendre and respective weights are found in Table 3. The results obtained with this method, whose handicap is the high computing cost, can be seen in Table 4.

i xi ω(xi)
1 0.9879925 0.0307532
2 0.9372734 0.0703660
3 0.8482066 0.1071592
4 0.7244177 0.1395707
5 0.5709722 0.1662692
6 0.3941513 0.1861610
7 0.2011941 0.1984315
8 0.0000000 0.2025782
9 -0.2011941 0.1984315
10 -0.3941513 0.1861610
11 -0.5709722 0.1662692
12 -0.7244177 0.1395707
13 -0.8482066 0.1071592
14 -0.9372734 0.0703660
15 -0.9879925 0.0307532

Table 3.

Points of Legendre ηi and respective wi weights.

ξ/x 0 0.5 1 2 4 6 8 10 20 40
0.01 0.00881 0.00881 0.00881 0.00881 0.00881 0.00880 0.00880 0.00879 0.00873 0.00847
0.02 0.01753 0.01753 0.01752 0.01752 0.01750 0.01746 0.01742 0.01735 0.01685 0.01496
0.03 0.02614 0.02614 0.02614 0.02612 0.02605 0.02594 0.02578 0.02557 0.02393 0.01836
0.04 0.03466 0.03466 0.03465 0.03461 0.03445 0.03418 0.03381 0.03333 0.02965 0.01857
0.05 0.04309 0.04308 0.04306 0.04298 0.04267 0.04216 0.04145 0.04055 0.03380 0.01639
0.10 0.08384 0.08379 0.08364 0.08305 0.08073 0.07700 0.07208 0.06623 0.03291 0.00262
0.15 0.12239 0.12223 0.12176 0.11989 0.11268 0.10165 0.08805 0.07328 0.01695 0.00080
0.20 0.15889 0.15854 0.15748 0.15331 0.13777 0.11540 0.09027 0.06614 0.00713 0.00070
0.25 0.19347 0.19281 0.19086 0.18325 0.15584 0.11934 0.08277 0.05253 0.00394 0.00067
0.30 0.22624 0.22516 0.22197 0.20968 0.16729 0.11571 0.07043 0.03881 0.00314 0.00065
0.35 0.25731 0.25569 0.25091 0.23271 0.17288 0.10713 0.05726 0.02816 0.00289 0.00064
0.40 0.28679 0.28450 0.27776 0.25245 0.17360 0.09604 0.04569 0.02110 0.00277 0.00064
0.45 0.31477 0.31168 0.30261 0.26909 0.17052 0.08439 0.03670 0.01687 0.00270 0.00064
0.50 0.34135 0.33733 0.32557 0.28286 0.16469 0.07346 0.03025 0.01446 0.00266 0.00063

Table 4.

Reference values for Doppler Broadening Function ψ(x,ξ).

Advertisement

7. Conclusion

A brief retrospective look at the calculation methodologies for the Doppler broadening function considering the approximations of Beth-Plackzec according to the formalism ofBriet-Wigner was presented in this chapter.

Advertisement

Acknowledgments

This research project is supported by Instituto Alberto Luiz Coimbra de Pós-graduação e PesquisaemEngenharia/Universidade Federal do Rio de Janeiro (COPPE/UFRJ), ComissãoNacional de Energia Nuclear (CNEN), Centro de Desenvolvimento de Tecnologia Nuclear (CDTN). It was financially supported by through National Institute of Science and Technology of Innovative Nuclear Reactors,Brazilian Council for Scientific and Technological Development (CNPq,) and Research Support Foundation of the State of Minas Gerais (FAPEMIG).

References

  1. 1. Arfken, G. (1985). Mathematical Method for Physicists. Academic Press Inc, London.
  2. 2. Bell, G.I., Glasstone, S. (1970).Nuclear Reactor Theory. Van Nostrand Reinhold Co., New York.
  3. 3. Campos, T.P.R. & Martinez, A.S. (1989). Approximate Calculation of the Resonance Integral for Isolated Resonances. Journal of Nuclear Science and Technology, Vol. 102, No. 3, (July 1989), pp. 211-218, ISSN 0029-5639.
  4. 4. Campos, T.P.R., Martinez, A.S. (1987). The dependence of pratical width on temperature. Annals of Nuclear Energy, Vol. 34, No. 1-2, (May 1987), pp. 68-82, ISSN 0306-4549.
  5. 5. Duderstadt, J. J, Hamilton, L.J. (1976). Nuclear Reactor Analysis. John Wiley and Sons, New York.
  6. 6. Gonçalves, A.C., Martinez, A.S., Silva, F.C. (2008). Solution of the Doppler broadening function based on the Fourier cosine transform. Annals of Nuclear Energy, Vol. 35, No. 10, (June 2008), pp. 1871-1881, ISSN 0306-4549.
  7. 7. Keshavamurthy, R.S., Harish, R. (1993). Use of Padé Approximations in the Analytical Evaluation of the J(θ,β) Function and its Temperature Derivative. Nuclear Science and Engineering, Vol. 115, No. 1, (September 1993), pp. 81-88, ISSN 0029-5639.
  8. 8. Mamedov, B.A. (2009). Analytical evaluation of Doppler functions arising from resonance effects in nuclear processes. Nuclear Instruments & Methods in Physics Research. Section A, Accelerators, Spectrometers, Detectors and Associated, Vol. 608, (July 2009), pp. 336-338, ISSN 0306-4549.
  9. 9. Palma, D.A.P., Martinez, A. S (2009). A faster procedure for the calculation of the J(ζ,β)). Annals of Nuclear Energy, Vol. 36, N0.10, (October 2009), pp. 1516-1520, ISSN 0168-9002.
  10. 10. Palma, D.A.P., Martinez, A.S., Silva, F.C. (2005). The Derivation of the Doppler Broadening Function using Frobenius Method. Journal of Nuclear Science and Technology, Vol. 43, No. 6, (December 2005), pp. 617-622, ISSN 0022-3131.
  11. 11. Polyanin, A.D &Manzhirov, A.V. (1998). Handbook of Integral Equation. CRC Press, New York.
  12. 12. Shcherbakov, O., Harada, H. (2002). Resonance Self-Shielding Corrections for Activations Cross Section Measurements. Journal of Nuclear Science and Technology, Vol. 39, No. 5, (May 2002), pp. 548-553, ISSN 0022-3131.
  13. 13. Stacey, W.M. (2001). Nuclear Reactor Physics. John Wiley and Sons, New York.

Written By

Daniel Artur P. Palma, Alessandro da C. Gonçalves, Aquilino Senra Martinez and Amir Zacarias Mesquita

Submitted: 16 May 2012 Published: 06 February 2013