Results for “ground deformation” — Page 15

DOE OSTI2023

Recent Developments in Hydrogeologic Applications for Strain Tensor Analyses

Changes in fluid pressure deform porous media and this effect occurs in a variety of hydrogeologic processes, from the change in storage during pumping or injection to fluctuations in water levels caused by barometric pressure. We have developed instruments for measuring small strains in porous media, and we have used the resulting strain data to evaluate well testing, hydraulic fracturing, manual loading at the ground surface, and ambient hydrologic processes, like rainfall and evaporation. A particularly important application is the use of strain tensor data measured at shallow depths to analyze well tests or hydraulic fractures conducted at much greater depths. An early demonstration of this technique was conducted at the North Avant Field north of Tulsa, Oklahoma, where Pennsylvanian sandstone creates a confined aquifer and oil reservoir at a depth of 530m. We have showed that the strains caused by injecting into the aquifer could be measured at a depth of 30m and used to evaluate the properties of the aquifer. We recently expanded the array of strainmeters at the North Avant Field by deploying three more instruments at shallow depth (30m) along with a deep instrument at 520m depth in the winter, 2021. The deep instrument is deployed in shale caprock slightly above the aquifer. To our knowledge, the deep strainmeter at the North Avant Field is the deepest strainmeter ever deployed and this required refining methods originally developed for shallow deployments. The instrument was lowered to depth on oil field tubing and cemented in place using techniques and materials developed for use in oil wells. Optical fiber used to communicate with the instrument was cemented in the annulus on the outside of tubing. This is significant because the techniques we used could readily be extended to greater depth, suggesting that strainmeters can be deployed over a wide range of depths for monitoring critical subsurface processes. For example, it suggests that strainmeters could be deployed through the caprock to monitor for leaks from underlying CO2 storage reservoirs. The strainmeter array at the North Avant Field has been used to characterize deformation during a series of injection tests in the spring and summer, 2021. All the new strainmeters respond to pumping, and the strainmeters we installed earlier also responded. To our knowledge, this is the first application of well testing in a deep aquifer that was monitored by an array of strainmeters—our earlier work used strainmeters at a single location. We are currently analyzing the strain data using an analytical solution, a proxy-based Bayesian inversion algorithm, and other methods. Strainmeter data has also been used to characterize periodic pumping tests by us and Riley Blais. A periodically varying pumping rate causes both hydraulic head and strain signals that vary with the same period as the pumping. The peaks and troughs of the head in monitoring wells lag behind the peaks and troughs of the head in the pumping well, and this lag time increases with distance from the pumping well. The lag time of the pressure and the distance to the monitoring well can be used in a simple analysis to estimate the hydraulic diffusivity of the aquifer. The lag time determined from strain data can be used to estimate aquifer properties using the same analysis that works for the pressure only for strain data measured at particular locations. That is because the strain field in a confining unit advances upward, laterally and then downward even though the pressure in the underlying aquifer only advances laterally, according to our recent simulations. We have field data showing that a small periodic signal superimposed on an injection rate at the North Avant Field will create a periodic strain signal at shallow strainmeters. The field data and the recent simulations suggest that including a periodic component to injection or pumping and then monitoring the resulting strain signal could be a way to monitor the subsurface.

Murdoch, Larry↗

DOE OSTI2023

Recent Developments in Hydrogeologic Applications for Strain Tensor Analyses

Changes in fluid pressure deform porous media and this effect occurs in a variety of hydrogeologic processes, from the change in storage during pumping or injection to fluctuations in water levels caused by barometric pressure. We have developed instruments for measuring small strains in porous media, and we have used the resulting strain data to evaluate well testing, hydraulic fracturing, manual loading at the ground surface, and ambient hydrologic processes, like rainfall and evaporation. A particularly important application is the use of strain tensor data measured at shallow depths to analyze well tests or hydraulic fractures conducted at much greater depths. An early demonstration of this technique was conducted at the North Avant Field north of Tulsa, Oklahoma, where Pennsylvanian sandstone creates a confined aquifer and oil reservoir at a depth of 530m. We have showed that the strains caused by injecting into the aquifer could be measured at a depth of 30m and used to evaluate the properties of the aquifer. We recently expanded the array of strainmeters at the North Avant Field by deploying three more instruments at shallow depth (30m) along with a deep instrument at 520m depth in the winter, 2021. The deep instrument is deployed in shale caprock slightly above the aquifer. To our knowledge, the deep strainmeter at the North Avant Field is the deepest strainmeter ever deployed and this required refining methods originally developed for shallow deployments. The instrument was lowered to depth on oil field tubing and cemented in place using techniques and materials developed for use in oil wells. Optical fiber used to communicate with the instrument was cemented in the annulus on the outside of tubing. This is significant because the techniques we used could readily be extended to greater depth, suggesting that strainmeters can be deployed over a wide range of depths for monitoring critical subsurface processes. For example, it suggests that strainmeters could be deployed through the caprock to monitor for leaks from underlying CO2 storage reservoirs. The strainmeter array at the North Avant Field has been used to characterize deformation during a series of injection tests in the spring and summer, 2021. All the new strainmeters respond to pumping, and the strainmeters we installed earlier also responded. To our knowledge, this is the first application of well testing in a deep aquifer that was monitored by an array of strainmeters—our earlier work used strainmeters at a single location. We are currently analyzing the strain data using an analytical solution, a proxy-based Bayesian inversion algorithm, and other methods. Strainmeter data has also been used to characterize periodic pumping tests by us and Riley Blais. A periodically varying pumping rate causes both hydraulic head and strain signals that vary with the same period as the pumping. The peaks and troughs of the head in monitoring wells lag behind the peaks and troughs of the head in the pumping well, and this lag time increases with distance from the pumping well. The lag time of the pressure and the distance to the monitoring well can be used in a simple analysis to estimate the hydraulic diffusivity of the aquifer. The lag time determined from strain data can be used to estimate aquifer properties using the same analysis that works for the pressure only for strain data measured at particular locations. That is because the strain field in a confining unit advances upward, laterally and then downward even though the pressure in the underlying aquifer only advances laterally, according to our recent simulations. We have field data showing that a small periodic signal superimposed on an injection rate at the North Avant Field will create a periodic strain signal at shallow strainmeters. The field data and the recent simulations suggest that including a periodic component to injection or pumping and then monitoring the resulting strain signal could be a way to monitor the subsurface.

Murdoch, Larry↗

DOE OSTI2021

Regional Subsidence Analysis Through a Multi-Scale Modeling Framework Based on Breakage Mechanics

Although poroelastic models are often used to explain the delay between subsurface fluid depletion and ground subsidence, inelastic compaction involving permanent changes of rock microstructure may exacerbate hydro-mechanical coupling, thus influencing the interpretation of measurements and long-term forecasts. Here, a multi-scale modeling approach is discussed, which accounts for the inherent connection between rock microstructure, hydraulic conductivity, and pore compaction. A constitutive model built within the framework of breakage mechanics is proposed to link the hydraulic conductivity of granular rocks with inelastic deformations and changes in grading caused by injection-depletion cycles at stress levels far from yielding. The proposed model has been incorporated into large-scale simulation frameworks, thus enabling the spatiotemporal mapping of regional subsidence through a hybrid, semi-analytical approach. Numerical results based on this strategy show that the model allows isolating near-field and far-field effects into the computation of land subsidence and can generate forecasts for different modeling scenarios (e.g., elastic and inelastic compaction, constant permeability, and concurrent change of compressibility and permeability). In particular, examples of simulations for the case of the Groningen gas field are discussed, showing the model capabilities to use both field measurements and laboratory tests for the generation of reasonable subsidence maps, without expensive computational costs. Results indicate that ignoring coupled inelastic effects has major consequences on the predicted timescale of subsidence. Specifically, while all the model scenarios produced similar long-term ground settlements, those ignoring breakage-dependent permeability changes result in a variation of the temporal window of residual subsidence of the order of several decades.

58 GEOSCIENCES↗

DOE OSTI2020

Arctic Coastal Erosion: Modeling and Experimentation

Increasing Arctic coastal erosion rates have put critical infrastructure and native communities at risk while also mobilizing ancient organic carbon into modern carbon cycles. Although the Arctic comprises one-third of the global coastline and has some of the fastest eroding coasts, current tools for quantifying permafrost erosion are unable to explain the episodic, storm-driven erosion events. Our approach, mechanistically coupling oceanographic predictions with a terrestrial model to capture the thermo-mechanical dynamics of erosion, enables this much needed treatment of transient erosion events. The Arctic Coastal Erosion Model consists of oceanographic and atmospheric boundary conditions that force a coastal terrestrial permafrost environment in Albany (a multi-physics based finite element model). An oceanographic modeling suite (consisting of WAVEWATCH III, Delft3D-FLOW, and Delft3D-WAVE) produced time-dependent surge and run-up boundary conditions for the terrestrial model. In the terrestrial model, a coupling framework unites the mechanical and thermal aspects of erosion. 3D stress/strain fields develop in response to a plasticity model of the permafrost that is controlled by the frozen water content determined by modeling 3D heat conduction and solid-liquid phase change. This modeling approach enables failure from any allowable deformation (block failure, slumping, etc.). Extensive experimental work has underpinned the ACE Model development including field campaigns to measure in situ ocean and erosion processes, strength properties derived from thermally driven geomechanical experiments, as well as extensive physical composition and geochemical analyses. Combined, this work offers the most comprehensive and physically grounded treatment of Arctic coastal erosion available in the literature. The ACE model and experimental results can be used to inform scientific understanding of coastal erosion processes, contribute to estimates of geochemical and sediment land-to-ocean fluxes, and facilitate infrastructure susceptibility assessments.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Direct observation of ultrafast symmetry reduction during internal conversion of 2-thiouracil using Coulomb explosion imaging

The photochemistry of heterocyclic molecules plays a decisive role for processes and applications like DNA photo-protection from UV damage and organic photocatalysis. The photochemical reactivity of heterocycles is determined by the redistribution of photoenergy into electronic and nuclear degrees of freedom, initially involving ultrafast internal conversion. Most heterocycles are planar in their ground state and internal conversion requires symmetry breaking. To lower the symmetry, the molecule must undergo an out-of-plane motion, which has not yet been observed directly. Here we show using the example of 2-thiouracil, how Coulomb explosion imaging can be utilized to extract comprehensive information on this molecular deformation, linking the extracted deplanarization of the molecular geometry to the previously studied temporal evolution of its electronic properties. Particularly, the protons of the exploded molecule are well-suited messengers carrying rich information on its geometry at distinct times after electronic excitation. We expect that our new analysis approach centered on these peripheral protons can be adapted as a general concept for future time-resolved studies of complex molecules in the gas phase.

Chemical physics↗

DOE OSTI2025

Twisting the Hubbard model into the momentum-mixing Hatsugai–Kohmoto model

The Hubbard model is a standard theoretical tool for studying materials with strong electron–electron interactions, such as cuprate superconductors. Unfortunately, interaction-driven phenomena, such as a transition into the strongly correlated Mott insulator phase, are difficult to treat with established theoretical techniques. However, the exactly solvable Hatsugai–Kohmoto model displays similar Mott physics. In this work we show how the Hatsugai–Kohmoto model can be deformed continuously into the Hubbard model. The trick is to systematically reintroduce all the momentum mixing that the original Hatsugai–Kohmoto model omits. This can be accomplished by grouping n momenta into a cell and hybridizing them, resulting in the momentum-mixing Hatsugai–Kohmoto model. We recover the Bethe ansatz ground-state energy of the one-dimensional Hubbard model to within 1% from only ten mixed momenta. Overall, the convergence scales as 1/n2 as opposed to the inverse linear behaviour of standard finite-cluster techniques. Our results for a square lattice reproduce all the known features from state-of-the-art simulations also with only a few mixed momenta. Consequently, we believe that the momentum-mixing Hatsugai–Kohmoto model offers an alternative tool for strongly correlated quantum matter.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Core of 25 F in the rotational model

In a recent experiment, carried out at the Radioactive Isotope Beam Factory/RIKEN, the 25 F(p,2p) 24 O reaction was studied at 270 MeV/A in inverse kinematics. Derived spectroscopic factors suggest that the effective core of 25 F significantly differs from a free 24 O nucleus. We interpret these results within the particle-rotor model and show that the experimental level scheme of 25 F can be understood in the rotation-aligned coupling scheme with its 5/2$^{+}_{1}$ ground state as the bandhead of a decoupled band. Here, the excitation energies of the observed 1/2$^{+}_{1}$ and 9/2$^{+}_{1}$ states correlate strongly with the rotational energy of the effective core, seen by the odd proton, and allow us to estimate its 2 + energy at ≈3.2MeV and a moderate quadrupole deformation ϵ 2 ≈0.15. The measured fragmentation of the πd 5/2 single-particle strength is discussed, and some further experiments suggested.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Induced seismicity and surface deformation associated with long-term and abrupt geothermal operations in Blue Mountain, Nevada

Not Available

58 GEOSCIENCES↗

NASA NTRS2016

Asteroid Fragmentation Approaches for Modeling Atmospheric Energy Deposition

During asteroid entry, energy is deposited in the atmosphere through thermal ablation and momentum-loss due to aerodynamic drag. Analytic models of asteroid entry and breakup physics are used to compute the energy deposition, which can then be compared against measured light curves and used to estimate ground damage due to airburst events. This work assesses and compares energy deposition results from four existing approaches to asteroid breakup modeling, and presents a new model that combines key elements of those approaches. The existing approaches considered include a liquid drop or "pancake" model where the object is treated as a single deforming body, and a set of discrete fragment models where the object breaks progressively into individual fragments. The new model incorporates both independent fragments and aggregate debris clouds to represent a broader range of fragmentation behaviors and reproduce more detailed light curve features. All five models are used to estimate the energy deposition rate versus altitude for the Chelyabinsk meteor impact, and results are compared with an observationally derived energy deposition curve. Comparisons show that four of the five approaches are able to match the overall observed energy deposition profile, but the features of the combined model are needed to better replicate both the primary and secondary peaks of the Chelyabinsk curve.

Modeling↗

DOE OSTI2020

Large classes of quantum scarred Hamiltonians from matrix product states

Motivated by the existence of exact many-body quantum scars in the Affleck-Kennedy-Lieb-Tasaki (AKLT) chain, in this work we explore the connection between matrix product state (MPS) wave functions and many-body quantum scarred Hamiltonians. We provide a method to systematically search for and construct parent Hamiltonians with towers of exact eigenstates composed of quasiparticles on top of an MPS wave function. These exact eigenstates have low entanglement in spite of being in the middle of the spectrum, thus violating the strong eigenstate thermalization hypothesis. Using our approach, we recover the AKLT chain starting from the MPS of its ground state, and we derive the most general nearest-neighbor Hamiltonian that shares the AKLT quasiparticle tower of exact eigenstates. We further apply this formalism to other simple MPS wave functions, and derive families of Hamiltonians that exhibit AKLT-like quantum scars. As a consequence, we also construct a scar-preserving deformation that connects the AKLT chain to the integrable spin-1 pure biquadratic model. Finally, we also derive other families of Hamiltonians that exhibit types of exact quantum scars, including a $\textit{U}$(1)-invariant perturbed Potts model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1986

The Geoscience Laser Altimetry/Ranging System (GLARS)

The Geoscience Laser Altimetry Ranging System (GLARS) is a highly precise distance measurement system to be used for making extremely accurate geodetic observations from a space platform. It combines the attributes of a pointable laser ranging system making observations to cube corner retroreflectors placed on the ground with those of a nadir looking laser altimeter making height observations to ground, ice sheet, and oceanic surfaces. In the ranging mode, centimeter-level precise baseline and station coordinate determinations will be made on grids consisting of 100 to 200 targets separated by distances from a few tens of kilometers to about 1000 km. These measurements will be used for studies of seismic zone crustal deformations and tectonic plate motions. Ranging measurements will also be made to a coarser, but globally distributed array of retroreflectors for both precise geodetic and orbit determination applications. In the altimetric mode, relative height determinations will be obtained with approximately decimeter vertical precision and 70 to 100 meter horizontal resolution. The height data will be used to study surface topography and roughness, ice sheet and lava flow thickness, and ocean dynamics. Waveform digitization will provide a measure of the vertical extent of topography within each footprint. The planned Earth Observing System is an attractive candidate platform for GLARS since the GLAR data can be used both for direct analyses and for highly precise orbit determination needed in the reduction of data from other sensors on the multi-instrument platform. (1064, 532, and 355 nm)Nd:YAG laser meets the performance specifications for the system.

Cohen, S. C.↗

NASA NTRS2016

The Adaptive Optics and Transmit System for NASA's Laser Communications Relay Demonstration Project

The Laser Communication Relay Demonstration is NASA’s multi-year demonstration of laser communication to a geosynchronous satellite. We are currently assembling the optical system for the first of the two baseline ground stations. The optical system consists of an adaptive optics system, the transmit system and a camera for target acquisition. The adaptive optics system is responsible for compensating the downlink beam for atmospheric turbulence and coupling it into the modem’s single mode fiber. The adaptive optics system is a woofer/tweeter design, with one deformable mirror correcting for low spatial frequencies with large amplitude and a second deformable mirror correcting for high spatial frequencies with small amplitude. The system uses a Shack-Hartmann wavefront sensor. The transmit system relays four beacon beams and one communication laser to the telescope for propagation to the space terminal. Both the uplink and downlink beams are centered at 1.55 microns. We present an overview of the design of the system as well as performance predictions including time series of coupling efficiency and expected uplink beam quality.

Laser Communication↗

NASA NTRS2012

Photogrammetry Measurements During a Tanking Test on the Space Shuttle External Tank, ET-137

On November 5, 2010, a significant foam liberation threat was observed as the Space Shuttle STS-133 launch effort was scrubbed because of a hydrogen leak at the ground umbilical carrier plate. Further investigation revealed the presence of multiple cracks at the tops of stringers in the intertank region of the Space Shuttle External Tank. As part of an instrumented tanking test conducted on December 17, 2010, a three dimensional digital image correlation photogrammetry system was used to measure radial deflections and overall deformations of a section of the intertank region. This paper will describe the experimental challenges that were overcome in order to implement the photogrammetry measurements for the tanking test in support of STS-133. The technique consisted of configuring and installing two pairs of custom stereo camera bars containing calibrated cameras on the 215-ft level of the fixed service structure of Launch Pad 39-A. The cameras were remotely operated from the Launch Control Center 3.5 miles away during the 8 hour duration test, which began before sunrise and lasted through sunset. The complete deformation time history was successfully computed from the acquired images and would prove to play a crucial role in the computer modeling validation efforts supporting the successful completion of the root cause analysis of the cracked stringer problem by the Space Shuttle Program. The resulting data generated included full field fringe plots, data extraction time history analysis, section line spatial analyses and differential stringer peak ]valley motion. Some of the sample results are included with discussion. The resulting data showed that new stringer crack formation did not occur for the panel examined, and that large amounts of displacement in the external tank occurred because of the loads derived from its filling. The measurements acquired were also used to validate computer modeling efforts completed by NASA Marshall Space Flight Center (MSFC).

Littell, Justin D.↗

NASA NTRSDate not supplied

Understanding the Vertical Slope of Mesoscale Snow Band Reflectivity Signatures

Mesoscale snow bands are frequently responsible for dumping large amounts of snow over a relatively small region and pose a considerable forecasting challenge. The ongoing NASA Investigation of Microphysics and Precipitation for Atlantic Coast-Threatening Snowstorms (IMPACTS) field campaign observed several instances of mesoscale snow bands during January and February of 2020 and 2022. These snow bands were sampled by a combination of airborne (NASA ER-2 and P-3) and ground-based remote sensing and in-situ instruments. Along-band cross sections of these snow bands reveal that many of them are vertically sloped. Our preliminary analysis into one of these snow bands suggests that the vertical slope is not maintained by isentropic lift and does not appear to be tied to horizontal deformation. A combined analysis of the reflectivity and wind retrievals from the ER-2 radars raises the possibility that these features could be produced by vertical wind shear layers deforming fall streaks into bands that extend hundreds of kilometers. The present study will focus on extending our preliminary analysis to additional cases and understanding the mechanisms that support these features using a combination of radar data, in-situ aircraft measurements, and numerical modeling.

Charles Helms↗

DOE OSTI2022

Machine learning approach to pattern recognition in nuclear dynamics from the ab initio symmetry-adapted no-core shell model

A novel machine learning approach is used to provide further insight into atomic nuclei and to detect orderly patterns amid a vast data of large-scale calculations. The method utilizes a neural network that is trained on ab initio results from the symmetry-adapted no-core shell model (SA-NCSM) for light nuclei. We show that the SA-NCSM, which expands ab initio applications up to medium-mass nuclei by using dominant symmetries of nuclear dynamics, can reach heavier nuclei when coupled with the machine learning approach. In particular, we find that a neural network trained on probability amplitudes for s- and p-shell nuclear wave functions not only predicts dominant configurations for heavier nuclei but in addition, when tested for the 20 Ne ground state, accurately reproduces the probability distribution. The non-negligible configurations predicted by the network provide an important input to the SA-NCSM for reducing ultralarge model spaces to manageable sizes that can be, in turn, utilized in SA-NCSM calculations to obtain accurate observables. The neural network is capable of describing nuclear deformation and is used to track the shape evolution along the 20-42 Mg isotopic chain, suggesting a shape coexistence that is more pronounced toward the very neutron-rich isotopes. We provide first descriptions of the structure and deformation of 24 Si and 40 Mg of interest to x-ray burst nucleosynthesis, and even of the extremely heavy nuclei such as 166,168 Er and 236 U, that build on first-principles considerations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Symmetry and Shape Coexistence in 10 Be

Within the low-lying spectrum of 10 Be, multiple rotational bands are found, with strikingly different moments of inertia. A proposed interpretation has been that these bands variously represent triaxial rotation and prolate axially-deformed rotation. The bands are well-reproduced in ab initio no-core configuration interaction (NCCI) calculations. Here, we use the calculated wave functions to elucidate the nuclear shapes underlying these bands, by examining the Elliott SU(3) symmetry content of these wave functions. The ab initio results support an interpretation in which the ground-state band, along with an accompanying K = 2 side band, represent a triaxial rotor, arising from an SU(3) irreducible representation in the 0hw space. Then, the lowest excited K = 0 band represents a prolate rotor, arising from an SU(3) irreducible representation in the 2hw space.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1962

Linearized Theory of Wind-tunnel Jet-boundary Corrections and Ground Effect for VTOL-STOL Aircraft

A linearized theory is developed to obtain interference factors for wind tunnels and ground effect. The calculated results, presented in tabular form in NASA Technical Notes D-933, D-934, D-935, and D-936, indicate that the degree to which the wake is deflected downward has a primary effect on the magnitude of the interference. When the wake is undeflected the corrections are essentially the same as those of classical theory. When the wake is deflected severely downward, the corrections are primarily determined by the floor of the wind tunnel. Under these conditions, the corrections are, in general, much larger than those given by previously available results. Because of deformations at the lower boundary of the jet, the corrections for an open wind tunnel are uncertain. Testing at low speeds and high lift coefficients in open wind tunnels, therefore, is not recommended. Although the theory is developed expressly for single-element, vanishingly small models, methods, with sample calculations, are given for extending the results to multielement and finite span models. The theory is at least partially verified by available test data; however, complete verification would require a substantially more detailed experimental study.

V/STOL aircraft↗

DOE OSTI2021

Possible quenching of static neutron pairing near the N=98 deformed shell gap: Rotational structures in 160,161 Gd

A 160 Gd beam was accelerated to an energy of 1000 MeV and, separately, bombarded thick targets of 154 Sm and 164 Dy in order to observe neutron-rich, rare-earth nuclei via deep-inelastic collision processes. Gammasphere was utilized to observe g-ray emissions. Here, many new states and transitions were observed in 160 Gd as a result of so-called "unsafe" Coulomb excitation. The ground-state band in 160 Gd has been extended to I π = 20 + and a rotational band based on the K π = 4 + state, previously associated with a hexadecapole vibration, was observed up to 18 + . The quasiparticle configuration of the K π = 4 + band has been determined, and its unusual alignment behavior may result from a possible quenching of static neutron pairing. In addition, the band based on the [523]5/2 quasineutron orbital in 161 Gd was extended from 11/2 – to 33/2 – , and also displays the same unusual alignment behavior.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗