Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reach code”

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 451 records · Page 25

LOW MACH NUMBER MODELING OF CORE CONVECTION IN MASSIVE STARS

This work presents three-dimensional simulations of core convection in a 15 M{sub Sun} star halfway through its main sequence lifetime. To perform the necessary long-time calculations, we use the low Mach number code MAESTRO, with initial conditions taken from a one-dimensional stellar model. We first identify several key factors that the one-dimensional initial model must satisfy to ensure efficient simulation of the convection process. We then use the three-dimensional simulations to examine the effects of two common modeling choices on the resulting convective flow: using a fixed composition approximation and using a reduced domain size. We find that using a fixed composition model actually increases the computational cost relative to using the full multi-species model because the fixed composition system takes longer to reach convection that is in a quasi-static state. Using a reduced (octant rather than full sphere) simulation domain yields flow with statistical properties that are within a factor of two of the full sphere simulation values. Both the octant and full sphere simulations show similar mixing across the convection zone boundary that is consistent with the turbulent entrainment model. However, the global character of the flow is distinctly different in the octant simulation, showing more rapid changes in the large-scale structure of the flow and thus a more isotropic flow on average.

79 ASTRONOMY AND ASTROPHYSICS↗

Results and Lessons Learned From the DOE Commercial Building Partnerships: Preprint

Over the course of 5 years, NREL worked with commercial building owners and their design teams in the DOE Commercial Building Partnerships (CBP) to cut energy consumption by 50% in new construction (versus code) and by 30% in existing building pilot projects (versus code or pre-retrofit operational energy use depending on the preference of the Partner) using strategies that could be replicated across their building portfolios. A number of different building types were addressed, including supermarket, retail merchandise, combination big box (general merchandise and food sales), high rise office space, and warehouse. The projects began in pre-design and included a year of measurement data to evaluate performance against design expectations. Focused attention was required throughout the entire process to achieve a design with the potential to hit the energy performance target and to operate the resulting building to reach this potential. This paper will report quantitative results and cover both the technical and the human sides of CBP, including the elements that were required to succeed and where stumbling blocks were encountered. It will also address the impact of energy performance goals and intensive energy modeling on the design process innovations and best practices.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Clad to coolant heat transfer during RIA in PWR conditions

The power pulse during an insertion of reactivity accident (RIA) leads to a transient heating of the fuel. Heat transfer toward the coolant through the clad and associated temperature transient of the clad is of safety significance since it may largely impact the risk of clad failure. For sufficiently large enthalpy deposited within the fuel, clad temperature increases rapidly and this can lead to coolant boiling. This paper concerns the critical review of the present understanding and of the models for the clad to coolant heat transfer in RIA for PWR conditions. The results of past experiments are reviewed and the reliability of current models dedicated to these conditions is analyzed. It appears that the use of classical models for wall to fluid heat transfer does not allow predicting the clad temperature transient as measured in RIA-related experiments. The deviation is all the more large when film boiling conditions are reached. To reproduce the experimental data, some correlations of heat transfer can be fitted but the reliability of such model out of the range of conditions of the experimental data is low. It is mainly due to a rather low level of understanding of the heat transfer in very transient conditions (i.e. when the time scale of variation of the main heat transfer parameters, like a wall temperature, is less than the typical internal time scale of the heat transfer process itself). It results in a very large scattering in the numerical simulation of clad temperature when several fuel codes are considered. To improve the reliability of these calculations, it is believed that the study of transient boiling heat transfer thanks to devoted experimental facilities is required.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Augmenting Molecular Graphs with Geometries via Machine Learning Interatomic Potentials

Accurate molecular property predictions require 3D geometries, which are typically obtained using expensive methods such as density functional theory (DFT). Here, we attempt to obtain molecular geometries by relying solely on machine learning interatomic potential (MLIP) models. To this end, we first curate a large-scale molecular relaxation dataset comprising 3.5 million molecules and 300 million snapshots. Then MLIP pre-trained models are trained with supervised learning to predict energy and forces given 3D molecular structures. Once trained, we show that the pre-trained models can be used in different ways to obtain geometries either explicitly or implicitly. First, it can be used to obtain approximate low-energy 3D geometries via geometry optimization. While these geometries do not consistently reach DFT-level chemical accuracy or convergence, they can still improve downstream performance compared to non-relaxed structures. To mitigate potential biases and enhance downstream predictions, we introduce geometry fine-tuning based on the relaxed 3D geometries. Second, the pre-trained models can be directly fine-tuned for property prediction when ground truth 3D geometries are available. Our results demonstrate that MLIP pre-trained models trained on relaxation data can learn transferable molecular representations to improve downstream molecular property prediction and can provide practically valuable but approximate molecular geometries that benefit property predictions. Our code is publicly available at: https://github.com/divelab/AIRS/.

Fu, Cong [Texas A & M Univ., College Station, TX (↗

Stop coannihilation in the CMSSM and SubGUT models

Stop coannihilation may bring the relic density of heavy supersymmetric dark matter particles into the range allowed by cosmology. The efficiency of this pro-cess is enhanced by stop-antistop annihilations into the lon-gitudinal (Goldstone) modes of the W and Z bosons, as well as by Sommerfeld enhancement of stop annihilations and the effects of bound states. Since the couplings of the stops to the Goldstone modes are proportional to the trilin-ear soft supersymmetry-breaking A -terms, these annihila-tions are enhanced when the A-terms are large. Yet, the Higgs mass may be reduced below the measured value if the A-terms are too large. Unfortunately, the interpretation of this constraint on the stop coannihilation strip is clouded by differences between the available Higgs mass calculators. For our study, we use as our default calcula-tor FeynHiggs 2.13.0, the most recent publicly avail-able version of this code. Exploring the CMSSM parameter space, we find that along the stop coannihilation strip the masses of the stops are severely split by the large A-terms. This suppresses the Higgs mass drastically for μ and A 0 > 0, whilst the extent of the stop coannihilation strip is limited for A 0 < 0 and either sign of μ. However, in sub-GUT models, reduced renormalization-group running mitigates the effect of the large A-terms, allowing larger LSP masses to be consistent with the Higgs mass calculation. We give examples where the dark matter particle mass may reach ≳ 8 TeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Kinetic simulation of the O-X conversion process in dense magnetized plasmas

One scheme for heating a dense magnetized plasma core, such as in a tokamak, involves launching an ordinary (O) electromagnetic wave at the low density edge. It is converted to a reflected extraordinary (X) electromagnetic wave under certain conditions, and then transformed into an electron Bernstein wave able to reach high density regions inaccessible to most other waves. The O-X mode conversion is important in heating and diagnostic processes in different devices such as tokamaks, stellarators, and some types of pinches. The goal of this study has been to demonstrate that the kinetic particle-in-cell (PIC) scheme is suitable for modeling the O-X conversion process as the first step toward a more complete simulation of O-X-B heating. The O-X process is considered and simulated with a kinetic particle model for parameters of the TJ-II stellarator using the PIC code, XOOPIC. This code is able to model the non-monotonic density and the magnetic profile of the TJ-II stellarator. It can also statistically represent the self-consistent distribution function of the plasma, which has not been possible in previous fluid models. By considering the electric and magnetic components of launched and reflected waves, the O-mode and X-mode waves can be detected, and the O-X conversion can be demonstrated. In this work, the optimum angle for conversion efficiency, as predicted by the previous theory and experimentally confirmed, is used. Via considering the power of the launched O-mode wave and the converted X-mode wave, the efficiency of 63% for O-X conversion for the optimum theoretical launch angle of 47{sup ∘} is obtained, which is in good agreement with efficiencies computed via full-wave simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

EXTENDED Ly{alpha} EMISSION FROM INTERACTING GALAXIES AT HIGH REDSHIFTS

Recent observations have discovered a population of extended Ly{alpha} sources, dubbed Ly{alpha} blobs (LABs), at high redshift z {approx} 2-6.6. These LABs typically have a luminosity of L {approx} 10{sup 42}-10{sup 44} erg s{sup -1}, and a size of tens of kiloparsecs, with some giant ones reaching up to D {approx} 100 kpc. However, the origin of these LABs is not well understood. In this paper, we investigate a merger model for the formation of LABs by studying Ly{alpha} emission from interacting galaxies at high redshifts by means of a combination of hydrodynamics simulations with three-dimensional radiative transfer calculations. Our galaxy simulations focus on a set of binary major mergers of galaxies with a mass range of 3-7 Multiplication-Sign 10{sup 12} M{sub Sun} in the redshift range z {approx} 3-7, and we use the newly improved ART{sup 2} code to perform the radiative transfer calculations, which couple multi-wavelength continuum, ionization of hydrogen, and Ly{alpha} line emission. We find that intense star formation and enhanced cooling induced by gravitational interaction produce strong Ly{alpha} emission from these merging galaxies. The Ly{alpha} emission appears to be extended due to the extended distribution of sources and gas. During the close encounter of galaxy progenitors when the star formation rate peaks at {approx}10{sup 3} M{sub Sun} yr{sup -1}, our model produces LABs with luminosity of L {approx} 10{sup 42}-10{sup 44} erg s{sup -1}, and size of D {approx} 10-20 kpc at z > 6 and D {approx} 20-50 kpc at z {approx} 3, in broad agreement with observations in the same redshift range. Our results suggest that merging galaxies may produce some typical LABs as observed, but the giant ones may be produced by mergers more massive than those in our model, or a combination of mergers and cold accretion from filaments on a large scale.

79 ASTRONOMY AND ASTROPHYSICS↗

24.77 Pflops on a Gravitational Tree-Code to Simulate the Milky Way Galaxy with 18600 GPUs, In: SC '14: Proceedings of the International Conference for High Performance Computing, Networking, Storage and Analysis

We have simulated, for the first time, the long term evolution of the Milky Way Galaxy using 51 billion particles on the Swiss Piz Daint supercomputer with our N-body gravitational tree-code Bonsai. Herein, we describe the scientific motivation and numerical algorithms. The Milky Way model was simulated for 6 billion years, during which the bar structure and spiral arms were fully formed. This improves upon previous simulations by using 1000 times more particles, and provides a wealth of new data that can be directly compared with observations. We also report the scalability on both the Swiss Piz Daint and the US ORNL Titan. On Piz Daint the parallel efficiency of Bonsai was above 95%. The highest performance was achieved with a 242 billion particle Milky Way model using 18600 GPUs on Titan, thereby reaching a sustained GPU and application performance of 33.49 Pflops and 24.77 Pflops respectively.

Computer Science↗

Simulations of high-velocity impacts on metal in preparation for the Psyche mission

In 2026, the NASA Discovery mission Psyche will orbit the asteroid (16) Psyche, the largest known metal-rich asteroid in the main belt. To estimate relative ages of the surface, identify re-surfacing events and better constrain Psyche's history, impact craters will be counted and characterized. No spacecraft has ever visited a metal-rich small body; therefore, laboratory-scale impact experiments and numerical simulations will play an important role in the interpretation of the mission's data. However, the planetary applications of high-velocity impacts have so far mostly been studied for silicate targets. Limited attention has been drawn to planetary objects predominantly made of metal, and more laboratory experiments and numerical calibrations are needed. As part of this effort, we present a suite of numerical simulations using an adaptative smoothed particles hydrodynamics numerical code (ASPH) reproducing a high-velocity impact experiment conducted on a steel target. This work primarily focuses on the influence of the chosen equation of state and initial distribution of flaws in the material on the estimated crater dimensions, damage and temperature. Additionally, we find that changing the EOS and initial flaw distribution affects the crater dimensions, though for a vast majority of cases the dimensions remain within 20% of the experimental values. The target is in most cases only locally weakened but not fully damaged, independently from the EOS chosen. Finally, temperatures at the impact point and around the forming crater can reach values above the melting point of iron at <100 GPa, which is in agreement with experimental observations. These results allow us to speculate on the differences expected between the surfaces of visited silicate-rich asteroids and that of the metal-rich target of the Psyche mission.

36 MATERIALS SCIENCE↗

Fast ion transport by sawtooth instability in the presence of ICRF-NBI synergy in JET plasmas

JET experiments have shown that the three-ion scenarios using waves in the ion cyclotron range of frequencies (ICRF) is an efficient way to build fast ion population through beam ion acceleration by radio frequency (RF) waves. Here, such a heating scheme is applied to plasmas with at least two thermal ion species. Analysis of mixed discharges with complex heating schemes requires a workflow that allows to model thermal and fast ion transport consistently. This paper is dedicated to modelling of a mixed plasma discharge with significant fraction of fast ions and contributes to development of fast ion transport models. For interpretive analysis with the TRANSP code a JET hydrogen-deuterium (HD) plasma discharge with neutral beam injection (NBI) and ICRF heating has been chosen. The task is complicated by NBI-ICRF synergy and plasma magnetohydrodynamic activity, like sawtooth crashes. D beam ions accelerated by RF waves form a high energy tail in fast ion distribution. Significant difference between the neutron rate computed by TRANSP and measured one is observed if the same diffusivity for electrons and ions is assumed. Sensitivity studies show that uncertainties in input plasma parameters and thermal ion transport models are crucial for modelling mixed plasma discharges and increased D transport is required to reach the plasma composition consistent with diagnostic measurements at the plasma edge. Fast ion redistribution by a sawtooth instability is characterised by non-resonant transport due to reconnection of magnetic field lines and resonant transport caused by resonance interaction between the instability and fast ions. With ORBIT simulations it has been shown that resonant interaction strongly affects fast ions of high energies, like beam ions accelerated by RF waves and fusion products. For the considered case, fast ion profiles simulated by ORBIT remain peaked after the sawtooth crashes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Strategic Petroleum Reserve Cavern Leaching Monitoring CY21

Th e U.S. Strategic Petroleum Reserve (SPR) is a crude oil storage system administered by the U.S. Department of Energy. The reserve consists of 60 active storage caverns located in underground salt domes spread across four sites in Louisiana and Texas, near the Gulf of Mexico. Beginning in 2016, the SPR started executing C ongressionally mandated oil sales. The configuration of the reserve, with a total capacity of greater than 700 million barrels ( MMB ) , re quires that unsaturated water (referred to herein as ?raw? water) is injected into the storage caverns to displace oil for sales , exchanges, and drawdowns . As such, oil sales will produce cavern growth to the extent that raw water contacts the salt cavern walls and dissolves (leaches) the surrounding salt before reaching brine saturation. SPR injected a total of over 45 MMB of raw water into twenty - six caverns as part of oil sales in CY21 . Leaching effects were monitored in these caverns to understand how the sales operations may impact the long - term integrity of the caverns. While frequent sonars are the most direct means to monitor changes in cavern shape, they can be resource intensive for the number of caverns involved in sales and exchanges. An interm ediate option is to model the leaching effects and see if any concerning features develop. The leaching effects were modeled here using the Sandia Solution Mining Code , SANSMIC . The modeling results indicate that leaching - induced features do not raise co ncern for the majority of the caverns, 15 of 26. Eleven caverns, BH - 107, BH - 110, BH - 112, BH - 113, BM - 109, WH - 11, WH - 112, WH - 114, BC - 17, BC - 18, and BC - 19 have features that may grow with additional leaching and should be monitored as leaching continues in th ose caverns. Additionally, BH - 114, BM - 4, and BM - 106 were identified in previous leaching reports for recommendation of monitoring. Nine caverns had pre - and post - leach sonars that were compared with SANSMIC results. Overall, SANSMIC was able to capture the leaching well. A deviation in the SANSMIC and sonar cavern shapes was observed near the cavern floor in caverns with significant floor rise, a process not captured by SANSMIC. These results validate that SANSMIC continues to serve as a useful tool for mon itoring changes in cavern shape due to leaching effects related to sales and exchanges.

02 PETROLEUM↗

Optical and ultraviolet observations of the very young Type IIP SN 2014cx in NGC 337

Extensive photometric and spectroscopic observations are presented for SN 2014cx, a Type IIP supernova (SN) exploding in the nearby galaxy NGC 337. The observations are performed in optical and ultraviolet bands, covering from −20 to +400 days from the peak light. The stringent detection limit from prediscovery images suggests that this supernova was actually detected within about one day after explosion. Evolution of the very early time light curve of SN 2014cx is similar to that predicted from a shock breakout and post-shock cooling decline before reaching the optical peak. Our photometric observations show that SN 2014cx has a plateau duration of ∼100 days, an absolute V-band magnitude of ∼−16.5 mag at t≈50 days, and a nickel mass of 0.056 ± 0.008 M{sub ⊙}. The spectral evolution of SN 2014cx resembles that of normal SNe IIP like SN 1999em and SN 2004et, except that it has a slightly higher expansion velocity (∼4200 km s{sup −1} at 50 days). From the cooling curve of photospheric temperature, we derive that the progenitor has a pre-explosion radius of ∼640 R {sub ⊙}, consistent with those obtained from SuperNova Explosion Code modeling (∼620 R {sub ⊙}) and hydrodynamical modeling of the observables (∼570 R {sub ⊙}). Moreover, the hydrodynamical simulations yield a total explosion energy of ∼0.4×10{sup 51} erg, and an ejected mass of ∼8 M{sub ⊙}. These results indicate that the immediate progenitor of SN 2014cx is likely a red supergiant star with a mass of ∼10 M{sub ⊙}.

79 ASTRONOMY AND ASTROPHYSICS↗

SATURATION OF THE MAGNETO-ROTATIONAL INSTABILITY IN STRONGLY RADIATION-DOMINATED ACCRETION DISKS

The saturation level of the magneto-rotational instability (MRI) in a strongly radiation-dominated accretion disk is studied using a new Godunov radiation MHD code in the unstratified shearing box approximation. Since vertical gravity is neglected in this work, our focus is on how the MRI saturates in the optically thick mid-plane of the disk. We confirm that turbulence generated by the MRI is very compressible in the radiation-dominated regime, as found by previous calculations using the flux-limited diffusion approximation. We also find little difference in the saturation properties in calculations that use a larger horizontal domain (up to four times the vertical scale height in the radial direction). However, in strongly radiation pressure dominated disks (one in which the radiation energy density reaches {approx}1% of the rest mass energy density of the gas), we find that Maxwell stress from the MRI turbulence is larger than the value produced when radiation pressure is replaced with the same amount of gas pressure. At the same time, the ratio between Maxwell stress and Reynolds stress is increased by almost a factor of eight compared with the gas pressure dominated case. We suggest that this effect is caused by radiation drag, which acts like bulk viscosity and changes the effective magnetic Prandtl number of the fluid. Radiation viscosity significantly exceeds both the microscopic plasma viscosity and resistivity, ensuring that radiation-dominated systems occupy the high magnetic Prandtl number regime. Nevertheless, we find that radiative shear viscosity is negligible compared to the Maxwell stress and Reynolds stress in the flow. This may have important implications for the structure of radiation-dominated accretion disks.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

{sup 6}LiF:ZnS(Ag) Neutrons Scintillator Detector Configuration for Optimal Readout

A Chromatic Analysis Neutron Diffractometer Or Reflectometer (CANDOR) is under development at the NIST Center for Neutron Research (NCNR). The CANDOR neutron sensor will rely on scintillator material for detecting the neutrons scattered by the sample under test. It consists of {sup 6}LiF:ZnS(Ag) scintillator material into which wavelength shifting (WLS) fibers have been embedded. Solid state photo-sensors (silicon photomultipliers) coupled to the WLS fibers are used to detect the light produced by the neutron capture event ({sup 6}Li (n,α) {sup 3}H reaction) and ionization of the ZnS(Ag). This detector configuration has the potential to accomplish the CANDOR performance requirements for efficiency of 90% for 5 A (3.35 meV) neutrons with high gamma rejection (10{sup 7}) along with compact design, affordable cost and materials availability. However this novel design includes challenges for precise neutron detection. The recognizing of the neutron signature versus the noise event produce by gamma event cannot be easy overcome by pulse height discrimination obstacle as can be achieved with {sup 3}He gas tube. Furthermore the selection of silicon photomultipliers (SiPM) as the light sensor maintains the obstacle of dark noise that does not exist when a photomultiplier tube is coupled to the scintillator. A proper selection of SiPM should focus on increasing the output signal and reducing the dark noise in order to optimize the detection sensitivity and to provide a clean signal pulse shape discrimination. The main parameters for evaluation are: - Quantum Efficiency (QE) - matching the SiPM peak QE with the peak transmission wavelength emission of the WLS. - Recovery time - a short recovery time is preferred to minimize the pulse width beyond the intrinsic decay time of the scintillator crystal (improves the gamma rejection based output pulse shape (time)). - Diode dimensions -The dark noise is proportional to the diode active area while the signal is provided by the WLS fibers; therefore the diode area should ideally be only minimally larger than fiber bundle area. - Low dark noise - it is desirable to minimize the dark noise during the pulse integration period so as to minimize the background for pulse shape discrimination. - Photon Detection Efficiency - it is desirable to increase the SiPM PDE in order to enhance light collection. This will increase the likelihood of detecting neutron events with lower light production and will present a cleaner raw signal for pulse shape discrimination. We will present the SiPM optimization process and studies of dark noise and gamma and neutron sensitivity as a function of bias voltage and operating temperature that have enabled us to optimize the detector sensitivity and gamma rejection. The gamma rejection performance goal requires to overcome the challenge of discriminating between the light signature accepted by neutron event to the one received by the noise. In addition there is a huge variation between the number of light photons that reaching the WLS fibers for different neutron events caused by the heavy ions energy losses prior to ionizing the ZnS(Ag) and the high light attenuation of the scintillation mixture. This variation in the light signal along with the long decay time of the ZnS(Ag) (tens of microseconds) can cause double counting of the same neutron event in the case of high light output signature or preventing the detection of low sequential light output signature neutron event. We will presents the algorithm developed for {sup 6}LiF:ZnS(Ag) sensor readout and the results achieved by an off-line analysis by Matlab software code that successfully achieved both the high gamma rejection with a sensitive and accurate neutron event detection. (authors)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Capability Demonstration of a 3D CdZnTe Detector on a High-Altitude Balloon Flight

In collaboration between the University of Michigan and Los Alamos National Laboratory, a 3D position-sensing CdZnTe (CZT) detector prototype was built and integrated into a high-altitude balloon platform to evaluate its performance in a space-like mixed-radiation environment. The detector prototype, Orion Eagle, was designed to operate in near-vacuum environments without any temperature regulation. Orion Eagle was hand-launched from NASA’s Columbia Scientific Balloon Facility (CSBF) at Fort Sumner, NM on September 26, 2021, and successfully operated throughout a 9-hour flight, which reached 38.5 km in altitude. The flight met its objectives, successfully detecting atmospheric gamma rays and galactic cosmic rays, and raising the Technical Readiness Level from 4 to 6 for large-volume 3D CZT detector technology for space applications. Ionization tracks produced by charged particles create spatial signatures in the detector that are distinguishable from discrete gamma-ray interactions. Therefore, the 3D position-sensing capabilities using pixelated electrodes on a CZT detector can help enable discrimination of background charged particles from gamma-ray events without an anticoincidence shield. The potential for background rejection capability, ambient-temperature operation, gamma-ray coded-aperture and Compton imaging, and near High Purity Germanium (HPGe) energy resolution motivate the use of large-volume 3D CZT imaging spectrometers in future space missions.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Turbulent dynamo action and its effects on the mixing at the convective boundary of an idealized oxygen-burning shell

Convection is one of the most important mixing processes in stellar interiors. Hydrodynamic mass entrainment can bring fresh fuel from neighboring stable layers into a convection zone, modifying the structure and evolution of the star. Because flows in stellar convection zones are highly turbulent, multidimensional hydrodynamic simulations are fundamental to accurately capture the physics of mixing processes. Under some conditions, strong magnetic fields can be sustained by the action of a turbulent dynamo, adding another layer of complexity and possibly altering the dynamics in the convection zone and at its boundaries. In this study, we used our fully compressible S EVEN -L EAGUE H YDRO code to run detailed and highly resolved three-dimensional magnetohydrodynamic simulations of turbulent convection, dynamo amplification, and convective boundary mixing in a simplified setup whose stratification is similar to that of an oxygen-burning shell in a star with an initial mass of 25 M ⊙ . We find that the random stretching of magnetic field lines by fluid motions in the inertial range of the turbulent spectrum (i.e., a small-scale dynamo) naturally amplifies the seed field by several orders of magnitude in a few convective turnover timescales. During the subsequent saturated regime, the magnetic-to-kinetic energy ratio inside the convective shell reaches values as high as 0.33, and the average magnetic field strength is ∼10 10 G. Such strong fields efficiently suppress shear instabilities, which feed the turbulent cascade of kinetic energy, on a wide range of spatial scales. The resulting convective flows are characterized by thread-like structures that extend over a large fraction of the convective shell. The reduced flow speeds and the presence of magnetic fields with strengths up to 60% of the equipartition value at the upper convective boundary diminish the rate of mass entrainment from the stable layer by ≈20% as compared to the purely hydrodynamic case.

79 ASTRONOMY AND ASTROPHYSICS↗

Integrated pre-feasibility study for CO 2 geological storage in the Cascadia Basin, offshore Washington State, British Columbia

The Cascadia CarbonSAFE project conducted a pre-feasibility study to evaluate technical and nontechnical aspects of storing 50 million metric tons (MMT) of carbon dioxide in a safe, ocean basalt reservoir offshore Washington State and British Columbia. Sub-seafloor basalts are very common on Earth and enable geological mineralization as a long-term storage mechanism, permanently sequestering the carbon in solid rock form. This project evaluated the offshore storage complex, developed potential industrial source/transport scenarios, built an inventory of existing geophysical/geological data and environmental monitoring capabilities, conducted laboratory studies and reservoir modeling to determine storage capacity, and analyzed economic factors, regulatory requirements, and project management risks. Our team included researchers at Columbia University in New York, Pacific Northwest National Laboratory (PNNL) in Washington, University of Victoria in British Columbia, University of California at Santa Cruz, and GHG Underground in Maine, as well as collaborators at University of Iceland in Reykjavik and University of Washington in Saint Louis. In the study region, large emitters (>100,000 MT/year) generate a total of approximately 40 MMT annually from stationary sources, such as power plants, refineries, ammonia production, and mineral processing plants. Five potential industry source/transport scenarios were identified (three in U.SA S. and two in Canada) to potentially provide 50 MMT of carbon dioxide to the offshore reservoir. One scenario included net carbon-negative sources that would reduce atmospheric levels. The inventory of existing data illustrated that basalt properties in the region are potentially beneficial for long-term storage. Permeability is on the order of 0.1 to 1 Darcies within the uppermost 600 m of ocean crust and density-derived porosity reaches 10–20% in thin layers that provide act as flow channels and provide access to porous and permeable basalt. Fine-grained sediments overlying basalt provide a low-permeability barrier separating the storage reservoir from the overlying ocean. Subsurface monitoring of the storage complex benefits from a unique opportunity close proximity of to the nearby NEPTUNE cabled seafloor observatory, which actively transmits data to shore and could be used to monitor reservoir properties. Injection pressure changes, reservoir chemistry, seismicity, well head corrosion, and potential leakage monitoring would be feasible in real time. Using existing data and published models, the STOMP-CO 2 modeling code developed at PNNL was used to simulate carbon dioxide injection, storage, and mineralization in the reservoir. Preliminary simulations indicate that a 50 MMT carbon dioxide plume injected over a 20-year period remains well within the reservoir area, both laterally and vertically, for at least a 50-year period of post-injection monitoring. Laboratory studies on core samples from Cascadia measured dissolution rates in CO 2 –basalt–seawater mixtures and indicated large variability in rates with the potential for complete mineralization. Compared to published data, the dissolution rate measured in these samples is faster, decreasing the time until full conversion to carbonate minerals would occur. Adding geochemical reactivity in the STOMP-CO 2 model indicates that carbon dioxide would be fully converted to carbonate minerals in 135 years or less after injection ceases. The project assessed the viability and costs for the five scenarios using either dedicated offshore pipelines or shipping vessels. Carbon dioxide capture costs were not considered. Shipping incurred lower upfront CAPEX costs, and because most viable sources were located near the shore and waterways, shipping is preferred; no dedicated carbon dioxide pipelines exist in the region. The project also considered potential risks from mixed-use impacts on surrounding lands and ocean areas, such as active use by commercial and recreational fisherman, and other community considerations in planning for public engagement in the future. Other project risks were identified, including project management, permitting and compliance, and financing, as well as technical issues such as long-term instrumentation, plume migration, and model development. Geological concerns like reservoir permeability, seal, and induced seismicity held relatively low risk. The regulatory framework for offshore storage in the U.S. and Canada may be covered by general legal/regulatory programs, but is not yet fully defined. Existing programs appear to restrict the transportation of carbon dioxide across national boundaries (e.g., Canada to U.S.) and its sub-seabed disposal in some circumstances. Our assessment suggests that developing specific and separate pilot projects will help to advance the technical feasibility and regulatory requirements in each country, and lessons learned from such activities may be transferable elsewhere around the globe.

42 ENGINEERING↗

Late-time mixing and turbulent behavior in high-energy-density shear experiments at high Atwood numbers

The LANL Shear Campaign uses millimeter-scale initially solid shock tubes on the National Ignition Facility to conduct high-energy-density hydrodynamic plasma experiments, capable of reaching energy densities exceeding 100 kJ/cm3. These shock-tube experiments have for the first time reproduced spontaneously emergent coherent structures due to shear-based fluid instabilities [i.e., Kelvin-Helmholtz (KH)], demonstrating hydrodynamic scaling over 8 orders of magnitude in time and velocity. The KH vortices, referred to as “rollers,” and the secondary instabilities, referred to as “ribs,” are used to understand the turbulent kinetic energy contained in the system. Their evolution is used to understand the transition to turbulence and that transition's dependence on initial conditions. Experimental results from these studies are well modeled by the RAGE (Radiation Adaptive Grid Eulerian) hydro-code using the Besnard-Harlow-Rauenzahn turbulent mix model. Information inferred from both the experimental data and the mix model allows us to demonstrate that the specific Turbulent Kinetic Energy (sTKE) in the layer, as calculated from the plan-view structure data, is consistent with the mixing width growth and the RAGE simulations of sTKE.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗