Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “EXPLOSIVE”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,333 records · Page 74

On the Symmetry of Blast Waves

This paper presents a brief historical review of G. I. Taylor’s solution of the point blast wave problem which was applied to the Trinity test of the first atomic bomb. Lie group symmetry techniques (also referred to throughout this paper as geometric techniques) are used to derive Taylor’s famous two-fifths law that relates the position of a blast wave to the time after the explosion and the total energy released. The theory of exterior differential systems is combined with the method of characteristics to demonstrate that the solution of the blast wave problem is directly related to the basic relationships that exist between the symmetry (or geometry) and the physics of wave propagation through the equations of motion. The point blast wave model is cast in terms of two exterior differential systems, and both systems are shown to be integrable with local solutions for the velocity, pressure, and density along curves in space and time behind the blast wave. This work is dedicated to the memory of Professor Roy Axford, who introduced many of his students to the topic of symmetry analysis of differential equations.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

CADIS and FW-CADIS Variance Reduction in Gamma Transport for Predicting Prompt Forensics Signatures

The goal of prompt nuclear forensics is to determine the characteristics of a nuclear detonation based on the signatures available almost immediately after the explosion. An important characteristic is the reaction time history (RTH), a measure of the device’s rate of neutron multiplication. The RTH can be estimated by observation of the gamma radiation emitted from the detonation, which can be detected directly or observed indirectly as Teller light. Gamma transport simulations used to predict these radiation fields are often modeled stochastically using the Monte Carlo N-Particle (MCNP) code, which can be a computationally demanding task due to the number of particle histories needed to achieve statistical convergence. In an attempt to improve the efficiency of these calculations, we evaluate two variance reduction techniques: Consistent Adjoint-Driven Importance Sampling (CADIS) and Forward-Weighted Consistent Adjoint-Driven Importance Sampling (FW-CADIS). These methods use a deterministically calculated adjoint flux to create weight windows and source biasing that guide MCNP sampling. We study the utility of CADIS and FW-CADIS for their use in MCNP gamma transport for nuclear forensics prediction simulations. Furthermore, the results demonstrate that both CADIS and FW-CADIS improve the accuracy for forensics-focused simulations, with CADIS being most beneficial in direct detection and FW-CADIS being ideal for computing a global Teller light source.

CADIS↗

Co-Active Subspace Methods for the Joint Analysis of Adjacent Computer Models

Active subspace (AS) methods are a valuable tool for understanding the relationship between the inputs and outputs of a Physics simulation. In this article, an elegant generalization of the traditional ASM is developed to assess the co-activity of two computer models. This generalization, which we refer to as a Co-Active Subspace (Co-AS) Method, allows for the joint analysis of two or more computer models allowing for thorough exploration of the alignment (or non-alignment) of the respective gradient spaces. We define co-active directions, co-sensitivity indices, and a scalar “concordance” metric (and complementary “discordance” pseudo-metric) and we demonstrate that these are powerful tools for understanding the behavior of a class of computer models, especially when used to supplement traditional AS analysis. Details for efficient estimation of the Co-AS and an accompanying R package (concordance) are provided. Practical application is demonstrated through analyzing a set of simulated rate stick experiments for PBX 9501, a high explosive, offering insights into complex model dynamics.

97 MATHEMATICS AND COMPUTING↗

PETN spark-gap detonators

The well-developed theory of Lorentz plasma that is dominated by electron–ion interactions is used to calculate the PETN arc characteristics. The spark-gap discharge current is a ramp with 10 to 25 ns rise time to peak and remaining constant subsequently. The approximate formulas for the arc channel conductivity, arc temperature, arc radius, and shock pressure from the arc are obtained from a system of nonlinear ordinary differential equations, which is the similarity solution of hydrodynamic equations similar to the Braginskii approximation. These arc parameters are given for the peak current ranging from 100 A to 1000 A and with different rise times. Representative cases are compared to the nonlinear ordinary differential equation code results. The shock pressures at the peak current are comparable to those from a typical commercial EBW bridgewire burst reported in the literature; the arc radius at the peak current is comparable to a typical bridgewire diameter of 0.0375 mm (e.g., RISI detonators, RP-1, and RP-80). The relevant Pop-Plot for low-density PETN is converted into an empirical detonation criterion, which is applicable to explosives subject to shocks of variable pressure. Finally, this criterion is then used to determine the detonation thresholds, which are comparable with test data obtained by Tucker, et al.

42 ENGINEERING↗

Quantifying CTFE content in FK-800 using ATR-FTIR and time to peak crystallization

The chlorotrifluoroethylene (CTFE) content in vinylidene fluoride (VDF)-CTFE copolymers greatly influences the chemical and physical properties and ultimately determines the application. The copolymer known as FK-800 has a VDF:CTFE ratio of ~1:3 and is used almost exclusively as the binder in the insensitive high explosive (HE) formulation PBX 9502. Due to the dangerous nature of HE work and the effect CTFE content variations have, the ability to quickly and easily quantify CTFE content is of interest. Two techniques, attenuated total reflectance-Fourier transform infrared spectroscopy (ATR-FTIR) and time to peak crystallization measured by heat flow calorimetry, were evaluated here as potential means for quantifying CTFE content. Based on the certificate of analyses (COAs), FK-800 lots with varied CTFE content encompassing the range historically used in HE applications were used as standards for both methodologies. For FTIR, a linear regression was performed on the peak area of the C-Cl IR stretch and the CTFE content; the methodology was then demonstrated on four samples, two lots of unadulterated material and material recovered from two HE samples. The calculated CTFE content and COAs were in good agreement for all four samples. Heat flow calorimetry revealed the relationship between time to peak crystallization and CTFE content was best fit by a power regression; the methodology was then tested on two different FK-800 lots. Only one sample showed good agreement between the calculated CTFE content and the COA. This discrepancy indicates the method developed is not suitable for quantifying CTFE, but provides valuable insight regarding the crystallization behavior of VDF-CTFE copolymers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Multi-Analysis Approach for Estimating Regional Health Impacts from the 2017 Northern California Wildfires

In the evening of October 8 and early hours of October 9, 2017, high winds in Northern California downed trees and power lines, igniting some of the most devastating wildfires the state had seen, and within hours unhealthy air quality impacted millions of people. We simulated these air quality conditions using fire detection information from the MODIS, VIIRS, and GOES-16 ABI satellite instruments, and applying a set of three WRF–CMAQ simulations, one data fusion, and three machine learning methods. We investigated using the 5-min available GOES-16 fire detection data to simulate timing of fire activity to allocate emissions hourly for the WRF-CMAQ air quality modeling system. Interestingly, this approach did not necessarily improve results compared to the baseline case, which used a default time profile. However, this approach was key to simulating the initial 12-hr explosive fire activity and smoke impacts. The WRF-CMAQ simulations compared well with observational data for the October 8-15 time period and tended to overestimate concentrations October 16-20. To improve these results, we applied three machine learning algorithms. We also had a unique opportunity to evaluate results with temporary monitors deployed specifically for wildfires, and performance was markedly different. For example, at the permanent monitoring locations, the WRF-CMAQ simulations had a Pearson correlation of 0.65, and the data fusion approach improved this (Pearson correlation = 0.95), while at the temporary monitor locations across the WRF-CMAQ, data fusion, and machine learning datasets, the best Pearson correlation was 0.5. The data fusion and machine learning results were biased low and WRF-CMAQ results were biased high. Finally, we applied the optimized PM2.5 exposure estimate in a short-term exposure-response function. Total estimated mortality attributable to PM2.5 exposure during the smoke episode was 83 (95% confidence interval: 0, 196) with 47% of these deaths attributable to wildland fire smoke.

O'Neill, Susan↗

Material Plasma Exposure eXperiment High Heat Flux Microwave Absorber Design, Manufacture, and Articles Test

The Material Plasma Exposure eXperiment (MPEX) at Oak Ridge National Laboratory is in the final design phase. MPEX will be capable of exposing neutron-irradiated materials to plasmas for the study of plasma-material interaction. This facility will provide information about the complex effects of plasmas on materials and contribute to examining new materials that can withstand high heat fluxes and high ion fluences for future fusion devices. MPEX plasma is heated by 70-GHz or 105-GHz electron Bernstein wave/electron cyclotron heating (ECH), and the high-frequency microwaves are prone to scattering microwave power, which can have detrimental effects, especially on diagnostic components. A large portion of the injected ECH power is expected to be absorbed by plasma, but the remainder requires that microwave absorbers be placed immediately upstream and downstream of the ECH launcher to minimize stray microwaves leaving the ECH region. These microwaves can inadvertently heat components that cannot be shielded or otherwise protected. The microwave absorber design is based on an array of pyramid-shaped ceramic tiles brazed to a water-cooled explosion-bonded heat sink and a stainless steel plate to produce one tile module. Computational fluid dynamics and structural analyses were performed to optimize and validate the design. Multiple test coupons were produced to validate the process for brazing the two different tile materials to the Glidcop AL-15 baseplate. The articles were tested to evaluate the reliability and thermal performance through exposure to an electron beam with a heat flux of up to 1.5 MW/m2. Nondestructive testing was performed before and after testing to identify voids or separations that may have been introduced by the high heat flux. This paper discusses the details of high heat flux microwave absorber design, manufacturing details and associated challenges, and test results, demonstrating the effectiveness of the proposed design.

Hussain, Aftab↗

Separation Process of Plant Fibers for Textile and Composite Application: A Review of Recent Advances

Plant fiber resources have gained significant attention for value-added utilization due to their renewability, sustainability, abundance, and widely acknowledged physical properties. The efficient and pragmatic separation of plant fibers is a critical process for their efficient utilization, yet a substantial gap persists between laboratory research advancements and their commercialization. To increase the possibility of research advancements for industrial application, this review summarizes the recent advances in different extraction methodologies of plant fiber research in textile and composite fields. It systematically outlines, compares, and contrasts physical (cryogenic, supercritical carbon dioxide, ultrasonic, steam explosion and microwave heating treatment), chemical (alkali, oxidation, organic solvents and deep eutectic solvents methods), and biological (natural retting, enzymatic and microorganism approaches) methods, addressing their respective mechanisms, strengths, limitations, research progress, and future prospects. In general, traditional chemical approaches have proven significantly effective but are accompanied by high pollution. Conversely, novel chemical treatments such as deep eutectic solvents and organic solvents offer a promising blend of efficiency and environmental friendliness but require deeper studies currently. Meanwhile, physical and biological treatments, though largely eco-friendly, tend to suffer from lower separation efficiencies. The research needs and future direction are also addressed to bridge the gap between scientific advancements and their widespread industrial application.

60 APPLIED LIFE SCIENCES↗

Time-resolved site-selective imaging of predissociation and charge transfer dynamics: the CH 3 I B-band

The predissociation dynamics of the 6s (B 2 E) Rydberg state of gas-phase CH 3 I were investigated by time-resolved Coulomb-explosion imaging using extreme ultraviolet (XUV) free-electron laser pulses. Inner-shell ionization at the iodine 4d edge was utilized to provide a site-specific probe of the ensuing dynamics. The combination of a velocity-map imaging (VMI) spectrometer coupled with the pixel imaging mass spectrometry (PImMS) camera permitted three-dimensional ionic fragment momenta to be recorded simultaneously for a wide range of iodine charge states. In accord with previous studies, initial excitation at 201.2 nm results in internal conversion and subsequent dissociation on the lower-lying A-state surface on a picosecond time scale. Examination of the time-dependent yield of low kinetic energy iodine fragments yields mechanistic insights into the predissociation and subsequent charge transfer following multiple ionization of the iodine products. The effect of charge transfer was observed through differing delay-dependencies of the various iodine charge states, from which critical internuclear distances for charge transfer could be inferred and compared to a classical over-the-barrier model. Time-dependent photofragment angular anisotropy parameters were extracted from the central slice of the Newton sphere, without Abel inversion, and highlight the effect of rotation of the parent molecule before dissociation, as observed in previous works. Our results demonstrate the ability to perform three-dimensional ion imaging at high event rates and showcase the potential benefits of this approach, particularly in relation to further time-resolved studies at free-electron laser facilities.

74 ATOMIC AND MOLECULAR PHYSICS↗

Imaging multiphoton ionization dynamics of CH 3 I at a high repetition rate XUV free-electron laser

XUV multiphoton ionization of molecules is commonly used in free-electron laser experiments to study charge transfer dynamics. However, molecular dissociation and electron dynamics, such as multiple photon absorption, Auger decay, and charge transfer, often happen on competing time scales, and the contributions of individual processes can be difficult to unravel. We experimentally investigate the Coulomb explosion dynamics of methyl iodide upon core–hole ionization of the shallow inner-shell of iodine (4d) and classically simulate the fragmentation by phenomenologically introducing ionization dynamics and charge transfer. Under our experimental conditions with medium fluence and relatively long XUV pulses (~75 fs), we find that fast Auger decay prior to charge transfer significantly contributes to the charging mechanism, leading to a yield enhancement of higher carbon charge states upon molecular dissociation. Furthermore, we argue for the existence of another charging mechanism for the weak fragmentation channels leading to triply charged carbon atoms. This study shows that classical simulations can be a useful tool to guide the quantum mechanical description of the femtosecond dynamics upon multiphoton absorption in molecular systems.

Optics↗

IceCube-Gen2: the window to the extreme Universe

The observation of electromagnetic radiation from radio to γ-ray wavelengths has provided a wealth of information about the Universe. However, at PeV (10 15 eV) energies and above, most of the Universe is impenetrable to photons. New messengers, namely cosmic neutrinos, are needed to explore the most extreme environments of the Universe where black holes, neutron stars, and stellar explosions transform gravitational energy into non-thermal cosmic rays. These energetic particles havemillions of times higher energies than those produced in the most powerful particle accelerators on Earth. As neutrinos can escape from regions otherwise opaque to radiation, they allow an unique view deep into exploding stars and the vicinity of the event horizons of black holes. The discovery of cosmic neutrinos with IceCube has opened this new window on the Universe. IceCube has been successful in finding first evidence for cosmic particle acceleration in the jet of an active galactic nucleus. Yet, ultimately, its sensitivity is too limited to detect even the brightest neutrino sources with high significance, or to detect populations of less luminous sources. In thiswhite paper, we present an overview of a next-generation instrument, IceCube-Gen2, which will sharpen our understanding of the processes and environments that govern the Universe at the highest energies. IceCube-Gen2 is designed to: (a) Resolve the high-energy neutrino sky from TeV to EeV energies (b) Investigate cosmic particle acceleration through multi-messenger observations (c) Reveal the sources and propagation of the highest energy particles in the Universe (d) Probe fundamental physics with high-energy neutrinos IceCube-Gen2 will enhance the existing IceCube detector at the South Pole. It will increase the annual rate of observed cosmic neutrinos by a factor of ten compared to IceCube, and will be able to detect sources five times fainter than its predecessor. Furthermore, through the addition of a radio array, IceCube- Gen2 will extend the energy range by several orders of magnitude compared to IceCube. Construction will take 8 years and cost about $350M. The goal is to have IceCube-Gen2 fully operational by 2033. IceCube-Gen2 will play an essential role in shaping the new era of multimessenger astronomy, fundamentally advancing our knowledge of the highenergy Universe. Overall, this challenging mission can be fully addressed only through the combination of the information from the neutrino, electromagnetic, and gravitational wave emission of high-energy sources, in concert with the new survey instruments across the electromagnetic spectrum and gravitational wave detectors which will be available in the coming years.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

16 O + 16 O fusion at sub-coulomb energies

A simple optical model (OPM) method using non-monotonic (NM) potentials characterized by a repulsive core from the microscopic theory of the Pauli-led energy-density functional (EDF) has been developed to investigate the 16 O + 16 O fusion at sub-Coulomb energies relevant to the oxygen burning. Here, the study involves the analysis of the experimental fusion cross-section (FCS) data in energy range 6.92 ≤ E cm ≤ 13.83 MeV, which includes the Coulomb barrier region E cm = 10.0–11.92 MeV. Apart from an excellent description of the existing FCS data in the studied energy range, the associated hindrance, characterized by the S-factor reaching a maximum and then gradually falling off at lower energies, so far suggested empirically for the system, is reproduced down to 4 MeV for the first time in the simple OPM. An S-factor maximum of S 0 = 3.15 × 10 25 MeV.b at E 0 = 7.14 MeV is observed with T ≃ 2.6 GK, which conforms to the values reported for quiescent and explosive burning. Our reaction rate, deduced from the NM potential, compares well with the Caughlan and Fowler data. Dominant partial waves implicit in the observed maximum S-factor in the studied Gamow range are also explored. Our present findings, with the success of NM potentials, suggest that the nucleus–nucleus potential is non-monotonic.

16O+16O↗

Developing solid-surface plasma facing components for pilot plants and reactors with replenishable wall claddings and continuous surface conditioning. Part A: concepts and questions

Abstract It is estimated that pilot plants and reactors may experience rates of net erosion and deposition of solid plasma facing component (PFC) material of 10 3 –10 5 kg yr −1 . Even if the net erosion (wear) problem can be solved, the redeposition of so much material has the potential for major interference with operation, including disruptions due to so-called ‘unidentified flying objects (UFOs)’ and unsafe dust levels. The potential implications appear to be no less serious than for plasma contact with the divertor target: a dust explosion or a major UFO-disruption could be as damaging for an actively-cooled deuterium-tritium (DT) tokamak as target failure. It will therefore be necessary to manage material deposits to prevent their fouling operation. This situation appears to require a fundamental paradigm shift with regard to meeting the challenge of taming the plasma–material interface: it appears that any acceptable solid PFC material will in effect be flow-through , like liquid–metal PFCs, although at far lower mass flow rates. Solid PFC material will have to be treated as a consumable , like brake pads in cars. ITER will use high-Z (tungsten) armor on the divertor targets and low-Z (beryllium) on the main walls. The ARIES-AT reactor design calls for a similar arrangement, but with SiC cladding on the main walls. Non-metallic low-Z refractory materials such as ceramics (graphite, SiC, etc) used as in situ replenishable, relatively thin—of order mm—claddings on a substrate which is resistant to neutron damage could provide a potential solution for the main walls, while reducing the risk of degrading the confined plasma. Separately, wall conditioning has proven essential for achieving high performance. For DT devices, however, standard methods appear to be unworkable, but recently powder droppers injecting low-Z material ∼continuously into discharges have been quite effective and may be usable in DT devices as well. The resulting massive generation of low-Z debris, however, has the same potential to seriously disrupt operation as noted above. Powder droppers provide a unique opportunity to carry out controlled studies on the management of low-Z slag in all current tokamaks, independent of whether their protection tiles use low-Z or high-Z material.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Present understanding of ignition and gain using indirect-drive inertial confinement fusion target designs on the U.S. National Ignition Facility

For many decades, the running joke in fusion research has been that 'fusion' is thirty years away and always will be. Yet, these past few years we find ourselves in a position where we can now talk about the milestones of burning plasmas, fusion ignition, and target energy gain greater than unity (scientific breakeven) in the past tense. Fusion is no longer a joke! Yet getting to fusion ignition, the tipping-point of thermonuclear instability resulting in an explosive increase in ion thermal temperature and fusion reaction-rate, and scientific breakeven (target gain, $G_{target} =$ fusion yield/deposited laser energy >1, in the laser-driven inertial confinement fusion context) has not been easy. Here, in this publication, we discuss our present understanding of the physics and technological challenges surrounding ignition and Gain as well as highlight some outstanding problems that still need resolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Neutrino halo effect on collective neutrino oscillation in iron core-collapse supernova model of a 9.6 M ⊙ star

We extend the multi-angle computational framework and investigate the time evolution of the neutrino halo on collective neutrino oscillation in the core collapse of an iron core progenitor. We find that in the case of the 9.6 M progenitor adopted in this work, there are windows of time when the effects of neutrino halo and collective neutrino oscillation are not simultaneously large. Inside the shock, the impact of the inward-scattered halo neutrino cannot in general be neglected compared to the outward-propagating neutrino flux. However, during early epochs, collective neutrino oscillation is effectively shut down by multi-angle matter suppression. During the intermediate epoch, collective neutrino oscillation is not suppressed, but its onset radius is beyond the still relatively small explosion shock front where the halo is prominent. We also find in the case of the 9.6 M progenitor the halo neutrinos induce a delay in the onset of collective neutrino oscillations. This causes novel flavor conversions which sharpen collective neutrino oscillation spectral features. We predict that the inclusion of neutrino halo effects makes neutrino signals that are more clearly distinct from thermal emission that when halo neutrinos are omitted.

79 ASTRONOMY AND ASTROPHYSICS↗

Flavor triangle of the diffuse supernova neutrino background

Although Galactic core-collapse supernovae (SNe) only happen a few times per century, every hour a vast number of explosions happen in the whole universe, emitting energy in the form of neutrinos, resulting in the diffuse supernova neutrino background (DSNB). The DSNB has not yet been detected, but Super-Kamiokande doped with gadolinium is expected to yield the first statistically significant observation within the next several years. Since the neutrinos produced at the core collapse undergo mixing during their propagation to Earth, the flavor content at detection is a test of oscillation physics. In this paper, we estimate the expected DSNB data at the DUNE, Hyper-K and JUNO experiments which when combined are sensitive to all different neutrino flavors. Here, we determine how well the flavor content of the DSNB will be reconstructed in the future, for a Mikheyev-Smirnov-Wolfenstein (MSW) scenario as well as a neutrino decay scenario. A large fraction of the flavor space will be excluded, but the heavy-lepton neutrino flux remains a challenge.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Exploring the MeV sky with a combined coded mask and Compton telescope: the Galactic Explorer with a Coded aperture mask Compton telescope (GECCO)

The sky at MeV energies is currently poorly explored. Here we present an innovative mission concept that builds upon the heritage of past and current missions improving the sensitivity and, very importantly, the angular resolution. This consists in combining a Compton telescope and a coded-mask telescope. We delineate the motivation for such a concept and we define the scientific goals for such a mission. The Galactic Explorer with a Coded Aperture Mask Compton Telescope (GECCO) is a novel concept for a next-generation telescope covering hard X-ray and soft gamma-ray energies. The potential and importance of this approach that bridges the observational gap in the MeV energy range are presented. With the unprecedented angular resolution of the coded mask telescope combined with the sensitive Compton telescope, a mission such as GECCO can disentangle the discrete sources from the truly diffuse emission. Individual Galactic and extragalactic sources are detected. This also allows to understand the gamma-ray Galactic center excess and the Fermi Bubbles, and to trace the low-energy cosmic rays, and their propagation in the Galaxy. Nuclear and annihilation lines are spatially and spectrally resolved from the continuum emission and from sources, addressing the role of low-energy cosmic rays in star formation and galaxy evolution, the origin of the 511 keV positron line, fundamental physics, and the chemical enrichment in the Galaxy. Such an instrument also detects explosive transient gamma-ray sources, which, in turn, enables identifying and studying the astrophysical objects that produce gravitational waves and neutrinos in a multi-messenger context. By looking at a poorly explored energy band it also allows discoveries of new astrophysical phenomena.

79 ASTRONOMY AND ASTROPHYSICS↗

Successful νp -process in neutrino-driven outflows in core-collapse supernovae

The origin of the solar system abundances of several proton-rich isotopes, especially 92,94 Mo and 96,98 Ru, has been an enduring mystery in nuclear astrophysics. An attractive proposal to solve this problem is the νp-process, which can operate in neutrino-driven outflows in a core-collapse supernova after the shock is launched. Years of detailed studies, however, have cast doubt over the ability of this process to generate sufficiently high absolute and relative amounts of various p-nuclei. The νp-process is also thought to be excluded by arguments based on the long-lived radionuclide 92 Nb. Here, we present explicit calculations, in which both the abundance ratios and the absolute yields of the p-nuclei up to A ≲ 105 are successfully reproduced, even when using the modern (medium enhanced) triple-α reaction rates. The process is also shown to produce the necessary amounts of 92 Nb. The models are characterized by subsonic outflows and by the protoneutron star masses in the ≳ 1.7 M ⊙ range. In conclusion, this suggests that the Mo and Ru p-nuclides observed in the Solar System were made in CCSN explosions characterized by an extended accretion stage.

79 ASTRONOMY AND ASTROPHYSICS↗