Inertial Confinement Fusion (ICF) is a process by which energy is produced by nuclear fusion reactions. The fuel pellet, generally called the target, is conventionally a spherical device which contains fuel for the fusion process. Lasers and ion beams have been considered for energizing and imploding the target. If the fuel contained in the target is sufficiently heated and compressed, a self-sustaining fusion reaction can occur.
The conventional approach to ICF target design is exemplified by the Department of Energy's program, National Ignition Facility (“NIF”). NIF target designs, as described in Lindl (Physics of Plasmas v11, number 2), consists of a mostly plastic or beryllium ablator region which surrounds a cryogenic DT ice, and a central void which is filled with very low density DT gas. The target is then placed in a cylindrical hohlraum. The entire target assembly (hohlraum and target) are then placed in the target chamber, and the hohlraum illuminated with a number (192) of discrete beams of laser light. The hohlraum then converts the energy to x-rays, which then ablate the ablator region, and by the reactive force, drives the DT inward. The combination of a non-spherical hohlraum and illumination leads directly to spatial non-uniformities in the target absorbed energy. Even in configurations with spherical hohlraums, the illumination is never spherically symmetric because entrance holes are required to admit the beams. This unsymmetrical illumination leads to unsymmetrical energy absorption by the target which in turn seeds instabilities that can prevent the temperature and density from achieving the necessary values to initiate a useful or self-sustaining fusion reaction.
In ICF targets, the energy absorbed by the generally spherical target needs to be uniform at the σabsE/Eabs≤1% level (RMS of absorbed energy divided by the absorbed energy) and in some arrangements «1%. In a hohlraum arrangement this is achievable with multiple incoming beams but adds additional requirements on laser or particle beam uniformity and further complicates the beam delivery system. It is the purpose of this invention to reduce σabsE/Eabs by trading off a multiplicity of beams for simpler illumination geometries by modifying some measureable properties of the target or hohlraum wall. It is recognized that feasible arrangements for production of power by ICF require energizing the target assembly over a relatively small solid angle (Tabak, 1997, UCRL-JC-128366 Rev 2, Design of a Distributed Radiator Target for Inertial Fusion Driven from Two Sides with Heavy Ion Beams). In commercial arrangements utilizing fusion power, the solid angle occupied by input beams (laser or particle) must necessarily be small and non-uniform so that the energy generated by a self-sustaining fusion reaction of the target can be efficiently absorbed by the reactor structure and because the target output (energy yielded by the fusion process) may be tailored to project energy into certain angular sectors that must be substantially free from beam entrance apertures.
Eabs
when the target has constant reflectivity Rdi.
Eabs
≤1%.
Eabs
≤1% overlaid with curves of constant absorption efficiency.
Eabs
≤1% overlaid with curves of constant absorption efficiency.
Eabs
≤1% overlaid with curves of constant absorption efficiency for an annular beam incident on the equator of a spherical target inside a spherical hohlraum.
Eabs
≤1% overlaid with curves of constant absorption efficiency.
Eabs
when target of
Eabs
.
In a target assembly 100 (
E
abs(x)=(1−Rdi)*Ti←o/(1−DS)Io (1)
where Rdi is the diffuse reflectivity of target surface 104, Ti←o is the radiation transfer operator from inner wall 112 of hohlraum 110 to the outer surface 104 of target 102, DS is the round trip scattering operator, and Io the externally incident intensity on inner hohlraum wall 112. The combined operator Ti←o/(1−DS) represents multiple scattering or rattling around of radiation between target 102 and hohlraum 110 and tends to spatially homogenize the radiation field. Consequently, varying the reflectivity of the hohlraum wall has less effect on the spatial uniformity of Eabs({right arrow over (x)}) than varying the reflectivity of target surface 104. For the same reason, spatial variation of Rdo (diffuse reflectivity of hohlraum surface 112) does not substantially degrade the homogenization effect of operator Ti←o/(1−DS). Based on this observation, to leading order the change in Rdi (ΔRdi) required to null out any non-uniformity in Eabs(x), ΔEabs(x)=Eabs(x)−Eabs
is
ΔRdi=(1−Rdi)*ΔEabs/Eabs
(2)
where Eabs
is the average energy/unit area absorbed by target 102. Equation 2 can be iterated to solution. Alternatively, and in conjunction with this relationship, ΔRi is expressed as a sum of spherical harmonics (Ylm′s) and the coefficients varied so as to minimize σabsE/
Eabs
(RMS of absorbed energy divided by the absorbed energy). During this computation, the reflectivity at all points on target surface 104 is held within physically realizable limits Rdi 0.4:0.8-0.9 for instance.
As a first example, beam 130 is incident on target assembly 100, enters hohlraum 110 through entrance aperture 114, is absorbed and reradiated by a volume absorber/radiator 120 resulting in an initial radiation pattern which in this example is approximated by a point source with angular distribution equal to the diffuse (Lambertian) reflection of a narrow beam at location of entrance aperture 114. Hohlraum inner wall surface 112 has a high constant diffuse reflectivity Rdo=0.9 and spherical shape with radius router=1 while target outer surface 104 has radius rinner=0.4 and constant diffuse reflectivity (Rdi=0.6). A polar plot (Eabs
shows variations as large as 3:1 and RMS nonuniformity (σabsE/
Eabs
) of 40.6%. After optimizing the reflectivity of 104 by varying the reflectivity of the first 10 (1=1:10) Legendre modes while keeping the average (Rdi(l=0)=0.6) constant the nonuniformity σabsE/
Eabs
is reduced to 0.6%, a 66-fold reduction. The resulting reflectivity profile (
Eabs
can be reduced to 1% or less by varying the first 10 (1=1:10) Legendre modes. The solid lines labeled 0.2, 0.3 etc. are contours of constant target absorption efficiency (
Eabs
/Etotal where Etotal=total energy initially reradiated by volume absorbers 120). As we see from the figure, net absorption efficiencies as high as 45% with absorption uniformities <1% are achievable with a single narrow incident beam. This is probably the most asymmetric geometric situation involving a concentric spherical target and hohlraum. More symmetric initial configurations lead to better results (wider operating window and greater energy deposition uniformity).
Another single sided illumination geometry (Eabs
<1% (triangles). Curves of constant (
Eabs
/Etotal overlay the operating envelope. The resulting operating envelope goes out to absorption efficiencies as high as ≈55%.
Yet another single sided illumination geometry (not shown) is similar to the annular illumination arrangement of
The simplest two sided illumination geometry (
Eabs
<1%.
In another embodiment (
As a specific example of optimizing the reflectivity of the target LH/Dtar=DH/Dtar=2, Eabs
for an uncompensated target. In this case, σabsE/
Eabs
=34%. After adjusting the reflectivity of target 102 for low order modes 1=1:10 (
Eabs
=0.13%, a 260× reduction. More complex illumination geometries involving narrow beams, multiple annular beams and combinations thereof can also have their target reflectivity optimized. So far the specific calculations displayed have cylindrical symmetry and we only needed to vary the spherical harmonics Yl,m=0(θ, ϕ)˜Pl(θ) along the axis of symmetry. Indeed, for the sphere in cylinder calculation (illustrated in
Beyond cylindrical and spherical hohlraums, this technique also applies to elliptical (
As described above in detail for some simple geometries we can also greatly decrease σabsE/Eabs
in more general target assembly geometries, including asymmetric hohlraums and targets. Also, any incident radiation condition whose spatial and angular radiation characteristics are known in advance can be compensated for to greatly decrease σabsE/
Eabs
. This includes finite area, unequal intensity beams.
In another aspect of this invention (generally referring to
In still another aspect of this invention, the shape of hohlraum inner surface 112 can be varied and this can be helpful in reducing the total required variation (max|θ,φ(Rdi(θ, ϕ))−min|θ,φ(Rdi(θ, ϕ)) of the target surface reflectivity.
The calculations presented so far utilized a time independent view factor calculation assuming a very short (delta function) beam temporal pulse length. The validity of these calculations was verified using time dependent view factor calculations with a non-zero pulse length.
Spatial variation of Rd (on target or hohlraum surface) requires a corresponding variation in materials.
If we cannot vary the target surface reflectivity (Rdi(θ, ϕ)) using discrete tessellation, then thin lamina consisting of two or more materials are deposited in alternating layers whose thickness varies continuously with position (θ, ϕ) on the target surface.
In the foregoing we have referred to volume absorbers or volume radiators 120.
In yet another aspect of this invention, volume absorbers embedded within hohlraum wall 112 (as in
As an example, in a target assembly (
Target assembly 1300 comprises a central spherical fuel region, the inner fuel region 1302. Surrounding the inner fuel region 1302 is an inner shell 1304 and outer shell 1308. In the space between the inner shell 1304 and outer shell 1308 is an outer fuel region 1306. Surrounding the outer shell 1308 is a propellant region 1312. A plurality of gold foam radiators 1314 are arranged in a one-to-one correspondence with the cylindrical beam channels 1320 located in case 1316. The cylindrical beam channels 1320 completely penetrate through case 1316. Relative to ICF target assembly previously discussed, case 1316 may be considered analogous to the hohlraum (110,
In another variation of this invention,
Additionally, the embodiments discussed in this application are exemplary and not an exhaustive enumeration of variants. Features discussed as part of separate embodiments may be combined into a single embodiment. Further, embodiments may make use of other features known in the art but not explicitly cited in this application.
This application claims the benefit of U.S. Provisional Application No. 62/546,440 filed on Aug. 16, 2017, which is incorporated herein by reference.
Number | Date | Country | |
---|---|---|---|
62546440 | Aug 2017 | US |