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.

SourceResearch connectionFrom yearThrough year

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

At least 487 records · Page 27

All sources NASA NTRS DOE OSTI

NASA NTRS2017

The 2017 Meteor Shower Activity Forecast for Earth Orbit

Most meteor showers will display typical activity levels in 2017. Perseid activity is expected to be higher than normal but less than in 2016; rates may reach 80% of the peak ZHR in 2016. Despite this enhancement, the Perseids rank 4th in flux for 0.04-cm-equivalent meteoroids: the Geminids (GEM), Daytime Arietids (ARI), and Southern delta Aquariids (SDA) all produce higher fluxes. Aside from heightened Perseid activity, the 2017 forecast includes a number of changes. In 2016, the Meteoroid Environment Office used 14 years of shower flux data to revisit the activity profiles of meteor showers included in the annual forecast. Both the list of showers and the shape of certain major showers have been revised. The names and three-letter shower codes were updated to match those in the International Astronomical Union (IAU) Meteor Data Center, and a number of defunct or insignificant showers were removed. The most significant of these changes are the increased durations of the Daytime Arietid (ARI) and Geminid (GEM) meteor showers. This document is designed to supplement spacecraft risk assessments that incorporate an annual averaged meteor shower flux (as is the case with all NASA meteor models). Results are presented relative to this baseline and are weighted to a constant kinetic energy. Two showers - the Daytime Arietids (ARI) and the Geminids (GEM) - attain flux levels approaching that of the baseline meteoroid environment for 0.1-cm-equivalent meteoroids. This size is the threshold for structural damage. These two showers, along with the Quadrantids (QUA) and Perseids (PER), exceed the baseline flux for 0.3-cm-equivalent particles, which is near the limit for pressure vessel penetration. Please note, however, that meteor shower fluxes drop dramatically with increasing particle size. As an example, the Arietids contribute a flux of about 5x10(exp -6) meteoroids m(exp -2) hr-1 in the 0.04-cm-equivalent range, but only 1x10(exp -8) meteoroids m(sub -2) hr-1 for the 0.3-cmequivalent and larger size regime. Thus, a PNP risk assessment should use the flux and flux enhancements corresponding to the smallest particle capable of penetrating a component, because the flux at this size will be the dominant contributor to the risk.

Moorhead, Althea↗

DOE OSTI2026

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 (↗

DOE OSTI2018

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↗

NASA NTRS2005

Buoyancy Suppression in Gases at High Temperatures

The computational fluid dynamics code FLUENT was used to study Rayleigh instability at large temperature differences in a sealed gas-filled enclosure with a cold top surface and a heated bottom wall (Benard problem). Both steady state and transient calculations were performed. The results define the boundaries of instability in a system depending on the geometry, temperature and pressure. It is shown that regardless of how fast the bottom-wall temperature can be ramped up to minimize the time spent in the unstable region of fluid motion, the eventual stability of the system depends on the prevailing final pressure after steady state has been reached. Calculations also show that the final state of the system can be different depending on whether the result is obtained via a steady-state solution or is reached by transient calculations. Changes in the slope of the pressure-versus-time curve are found to be a very good indicator of changes in the flow patterns in the system.

Kuczmarski, Maria A.↗

NASA NTRS2022

Bringing Planetary Science Mission Outreach to the Deaf and Blind Communities

Introduction: Technology for enhancing outreach, like 3D printing, and science communication products, such as videos and podcasts, can be utilized within the planetary science community, especially for the engagement and excitement of current or upcoming planetary exploration missions. However, these communication products can also be further enhanced for the benefit of the blind and deaf communities. While such products may already be readily available, such projects are not easily accessible to blind and/or deaf certified educators, which often rely on making their own resources or do not have the funds to provide such resources (e.g., cost of 3D printers or cost of braille books). The planetary science community can have better practices to reach these broader audiences. Best practices can include transcripts from podcasts, transcripts in videos, and large-font captions. Images on websites and social media accounts should also include alt-text descriptive captions. 3D printing can also enhance planetary science for the blind community, through tactile posters, maps, and pamphlets. Planetary data can also be augmented by providing different tactile geological maps (e.g., topography or various datasets), and audio-visual videos freely available for educators. Visual Engagement: Visual engagement consists of several avenues to consider, the three main themes includes: 1) swag; 2) videos; 3) interactive exploration. Swag can include the fun visual take-home materials, such as stickers, posters, bookmarks, etc. Videos can include educational-specific videos (available freely via YouTube or by other educational-specific streaming avenues, such as Nebula or Curiosity Stream), provided they have Closed Captioning (CC). The use of QR codes to such videos or websites can also benefit to being added on swag. Interactive exploration can also be sub-divided by different types of engagement. A popular and still fairly new technology for public engagement is the use of virtual reality (VR). While this has been mainly for martian and lunar surface exploration [1], the deaf communities can benefit from VR through a more extensive look at our solar system and beyond (for example, a VR experience of the flight path, or visual map of the heliosphere/dynamics of our Sun). Audio Engagement: Audio tools can also be a useful avenue of communication, especially for the blind communities. Audio archiving can certainly be transcripts from the video engagements, but also the use of podcasts can also be a benefit. Podcasting can take on two forms: 1) interview engagement; and 2) update engagement. For interviews, scientists can communicate with STEM-specific podcast platforms to make other listener-bases aware of what is going on with a specific mission. For update-type communication, missions may opt to have an archived podcast of news, updates, and the teams involved. The most important aspect of such podcasts would be for the need of complimentary transcripts (including descriptive transcripts if sounds are included), and the limited use of jargon. Research Engagement: There have been several examples of involving the blind and low-vision communities in citizen science, such as through the NASA Heliophysics division. Examples include the NASA PUNCH (Polarimeter to Unify the Corona and Heliosphere) mission led by the Southwest Research Institute [2], which include blind and visually-impaired citizens to assist in the Sun’s coronal rhythms.” Another example is the Eclipse Soundscapes: Citizen Science Project (ES:CSP), which documents observations of acoustical changes of nature and ecosystems during solar eclipse events [3]. Inclusivity: A major theme that is necessary for public engagement is inclusivity and the awareness of reaching broader audiences. Outreach to include hearing/seeing impaired communities are still lacking in the sciences. There are several opportunities that the planetary sciences could take. Other projects that have emerged from the space sciences include adding transcripts to visual engagement [4], and the use of 3D printing for the visually-impaired [5]. References: [1] Olgin, J. (2020) 51st LPSC, Abstract 2137. [2] https://scitechdaily.com/outreach-for-nasa-punch-mission-embraces-ancient-and-modern-sun-watching-theme/ [3] https://science.nasa.gov/science-activation-team/eclipse-soundscapes [4] NASA International Observe the Moon Night (Blind and Deaf Accessible), Youtube Video ( https://www.youtube.com/watch?v=neHCfg0S3-Q) [5] Richardson, J., et al. (2018) AGU Fall Meeting, Abstract ED23F-0963.

C J Ahrens↗

NASA NTRS1995

Electron distributions in quasi-perpendicular shocks: Simulations

One of the outstanding problems of collisionless shock physics is the microphysical description of electron thermalization and associated wave processes. Observations at the Earth's bow shock have shown various features in the electron velocity distribution such as a loss-cone, temperature anisotropies, and beams before a thermalized, flat-topped distribution is reached downstream. It is generally believed that in addition to the macroscopic fields, which determine the overall mapping of the distribution function, electromagnetic whistler waves and lower hybrid waves are responsible for the pitch-angle scattering and thermalization of the electrons. However, a genuine physical explanation requires analysis of the evolution of the electron distribution in the self-consistent wave fields. Particle simulations are ideally suited to do this, and include processes in which both electrons and ions are involved, such as cross-field streaming instabilities. We show results from l-D and 2-D simulations using an implicit code, which allows time steps large compared to the inverse electron plasma frequency. It is shown that in order to capture the correct propagation characteristics of the waves, a large mass ratio between the protons and electrons is generally necessary. Detailed diagnostics of the waves as well as the evolution of the electron velocity distribution will be presented and will be discussed in the frame work of satellite observations.

Pantellini, F.↗

DOE OSTI2019

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↗

DOE OSTI2021

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↗

DOE OSTI2022

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↗

NASA NTRS1996

CFD Study of Turbo-Ramjet Interactions in Hypersonic Airbreathing Propulsion System

Advanced airbreathing propulsion systems used in Mach 4-6 mission scenarios, usually involve turbo-ramjet configurations. As the engines transition from turbojet to ramjet, there is an operational envelope where both engines operate simultaneously. In the first phase of our study, an over/under nozzle configuration was analyzed. The two plumes from the turbojet and ramjet interact at the end of a common 2-D cowl, where they both reach an approximate Mach 3.0 condition and then jointly expand to Mach 3.6 at the common nozzle exit plane. For the problem analyzed, the turbojet engine operates at a higher nozzle pressure ratio than the ramjet, causes the turbojet plume overpowers the ramjet plume, deflecting it approximately 12 degrees downward and in turn the turbojet plume is deflected 6 degrees upward. In the process, shocks were formed at the deflections and a shear layer formed at the confluence of the two jets. This particular case was experimentally tested and the data were used to compare with a computational fluid dynamics (CFD) study using the PARC2D code. The CFD results were in good agreement with both static pressure distributions on the cowl separator and on nozzle walls. The thrust coefficients were also in reasonable agreement. In addition, inviscid relationships were developed around the confluence point, where the two exhaust jets meet, and these results compared favorably with the CFD results. In the second phase of our study, a 3-D CFD solution was generated to compare with the 2-D solution. The major difference between the 2-D and 3-D solutions was the interaction of the shock waves, generated by the plume interactions, on the sidewall. When a shock wave interacts with a sidewall and sidewall boundary layer, it is called a glancing shock sidewall interaction. These interactions entrain boundary layer flow down the shockline into a vortical flow pattern. The 3-D plots show the streamlines being entrained down the shockline. The pressure of the flow also decreases slightly as the sidewall is approached. Other difference between the 2-D and 3-D solutions were a lowering of the nozzle thrust coefficient value from 0.9850 (2-D) to 0.9807 (3-D), where the experimental value was 0.9790. In the third phase of our study, a different turbo-ramjet configuration was analyzed. The confluence of a supersonic turbojet and a subsonic ramjet in the turbine based combined-cycle (TBCC) propulsion system was studied by a 2-D CFD code. In the analysis, Mach 1.4 primary turbojet was mixed with the subsonic ramjet secondary flow in an ejector mode operation. Reasonable agreements were obtained with the supplied I-D TBCC solutions. For low downstream backpressure, the Fabri choke condition (Break-Point condition) was observed in the secondary flow within mixing zone. For sufficient high downstream backpressure, the Fabri choke no longer exist, the ramjet flow was reduced and the ejector flow became backpressure dependent. Highly non-uniform flow at ejector exit were observed, indicated that for smooth downstream combustion, the mixing of the two streams probably required some physical devices.

Chang, Ing↗

NASA NTRS2011

Prediction of Acoustic Loads Generated by Propulsion Systems

NASA Stennis Space Center is one of the nation's premier facilities for conducting large-scale rocket engine testing. As liquid rocket engines vary in size, so do the acoustic loads that they produce. When these acoustic loads reach very high levels they may cause damages both to humans and to actual structures surrounding the testing area. To prevent these damages, prediction tools are used to estimate the spectral content and levels of the acoustics being generated by the rocket engine plumes and model their propagation through the surrounding atmosphere. Prior to the current work, two different acoustic prediction tools were being implemented at Stennis Space Center, each having their own advantages and disadvantages depending on the application. Therefore, a new prediction tool was created, using NASA SP-8072 handbook as a guide, which would replicate the same prediction methods as the previous codes, but eliminate any of the drawbacks the individual codes had. Aside from replicating the previous modeling capability in a single framework, additional modeling functions were added thereby expanding the current modeling capability. To verify that the new code could reproduce the same predictions as the previous codes, two verification test cases were defined. These verification test cases also served as validation cases as the predicted results were compared to actual test data.

Perez, Linamaria↗

DOE OSTI2023

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 ↗

DOE OSTI2023

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↗

DOE OSTI2018

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↗

NASA NTRS1988

Driven, steady-state RFP computations

The pseudospectral three-dimensional MHD code of Dahlburg et al. (1986 and 1987) is used to compute the dynamical behavior of a channel of magnetofluid carrying an axial current and magnetic flux. This situation contains the essential MHD behavior of the reversed-field pinch (RFP). An externally imposed electric field is applied to an initially current-free magnetofluid and drives currents that rise and eventually fluctuate about values corresponding to pinch ratios Theta of about 1.3, 2.2, and 4.5. A period of violent turbulence leads to an approximately force-free core, surrounded by an active MHD boundary layer that is not force-free. A steady state is reached that can apparently be sustained indefinitely (for several hundred Alfven transit times or longer). The turbulence level and time variability in the steady state increase with increasing Theta. The average toroidal magnetic field at the wall reverses for Theta = 2.2 and 4.5, but not for Theta = 1.3. Negative toroidal current filaments are observed. The Lundquist numbers are of the order of a few hundred.

Dahlburg, J. P.↗

NASA NTRS2001

Jet Stability and the Generation of Superluminal and Stationary Components

We present a numerical simulation of the response of an expanding relativistic jet to the ejection of a superluminal component. The simulation has been performed with a relativistic time-dependent hydrodynamical code from which simulated radio maps are computed by integrating the transfer equations for synchrotron radiation. The interaction of the superluminal component with the underlying jet results in the formation of multiple conical shocks behind the main perturbation. These trailing components can be easily distinguished because they appear to be released from the primary superluminal component instead of being ejected from the core. Their oblique nature should also result in distinct polarization properties. Those appearing closer to the core show small apparent motions and a very slow secular decrease in brightness and could be identified as stationary components. Those appearing farther downstream are weaker and can reach superluminal apparent motions. The existence of these trailing components indicates that not all observed components necessarily represent major perturbations at the jet inlet; rather, multiple emission components can be generated by a single disturbance in the jet. While the superluminal component associated with the primary perturbation exhibits a rather stable pattern speed, trailing components have velocities that increase with distance from the core but move at less than the jet speed. The trailing components exhibit motion and structure consistent with the triggering of pinch modes by the superluminal component. The increase in velocity of the trailing components is an indirect consequence of the acceleration of the expanding fluid, which is assumed to be relativistically hot; if observed, such accelerations would therefore favor an electron-positron (as opposed to proton rest mass) dominated jet.

Agudo, Ivan↗

NASA NTRSDate not supplied

Open-source Numerical Modeling of Solidification Cracking Susceptibility: Application to Refractory Alloy Systems

Introduction. Alloys such as aluminum, nickel-base, and austenitic stainless steels are susceptible to solidification cracking during welding and 3D printing. Compositional optimization is one method used to effectively mitigate solidification cracking of those alloy systems. With the surge in hypersonic and in-space propulsion activities, refractory metals (Nb, Mo, Ta, W, and Re) and their alloy derivatives are increasing in importance due to their extreme high melting point and retention of high-temperature strength; however, their chemistry was most typically optimized to promote ductility during mechanical operations such as drawing and forming. Welding of such alloys has been a challenge due to a number of issues including solidification cracking, atmospheric contamination (O, C, and N), as well as a shift in ductile-to-brittle transition to higher temperature following grain growth induced by welding. Compositional optimization of refractory alloys for solidification cracking resistance in particular is desirable as their usage increases with the advent of advanced manufacturing methods such as 3D printing. This work evaluates the effect of compositional variation in refractory metal systems on the solidification cracking susceptibility with the goals of optimizing existing alloys and joining process techniques, and formulating new alloys with increased solidification cracking resistance. Experimental Procedures. A python code was developed in a Jupyter notebook environment (Michael and Sowards, 2023) to facilitate the calculation of crack susceptibility index proposed by Kou (2015). Composition is entered as a single point, or as a 1-D or 2-D array. The notebook calls pycalphad (Otis and Liu, 2017 and Bocklund et al, 2020) to calculate the evolution of fraction solid as a function of temperature (under either Scheil or equilibrium assumptions) and then evaluates steepness of the fraction solid curve near the terminal stage of solidification to predict solidification cracking resistance. Open source thermodynamic databases available at online repositories are used (van de Walle). The process is setup in an automated fashion to generate plots that show variation in solidification cracking susceptibility according to composition on 1-D line plots or 2-D contour plots. The Jupyter notebook and crack susceptibility algorithm was also integrated with a widely used commercial CALPHAD code for validation and alloy exploration. Results and Discussion. The crack susceptibility model was first validated against a series of refractory alloy compositions evaluated in past work which utilized a specialized Varestraint test built inside a vacuum chamber environment (Lessman and Gold, 1971). The alloys tested in the Varestraint apparatus included T-111 (Ta-8W-2Hf), ASTAR-811C (Ta-8W-1Re-0.7Hf-0.025C), FS-85 (Nb-27Ta-10W-1Zr), T-222 (Ta-9.6W-2.4Hf-0.01C), Ta-10W, B-66 (Nb-5Mo-5V-1Zr), and SCb-291 (Nb-10W-10Ta). The initial test of the model showed a strong correlation with empirical Varestraint data, i.e., a Spearman rank correlation between model predictions and hot cracking measurements was observed to be greater than 0.8. Following the validation, a set of refractory metal binary mixtures was investigated to evaluate sensitivity of Nb, Mo, W, and Ta to C, N, and O content. A series of plots were produced that suggest ppmw ranges of C, N, and O where solidification cracking increases significantly and reaches a maximum. Also comparative ranking of each primary refractory metal to each interstitial was produced. For example C produces greater cracking response in Mo whereas O produces greater cracking response in Ta and Nb. Such compositional values have utility in setting limits on pickup of these interstitial elements during welding and printing rather than using a one-size-fits-all approach. Furthermore, the results have use in determining additive powder recycling requirements, which is especially pertinent for refractory metal powders due to their high cost compared to conventional alloys. Another application created thousands of hypothetical alloys within the nominal specified composition range of two widely used refractory alloys C103 (Nb-10Hf-1Ti) and TZM (Mo-0.5Ti-0.1Zr). The cracking index was calculated for the alloys and results were fed into machine learning regression techniques including Multiple Linear Regression, Ridge Regression, and Lasso Regression to determine relative potency each alloying element had on computed solidification cracking index. A series of linear equations were produced that relate composition of C103 and TZM to solidification cracking index. The crack susceptibility of C103 for example is described by an equation of the form: cracking index ~ O + 0.667*C + 0.635*N + 0.00037*Ta – 0.0008*Hf (in wt.%) From that equation, it is clear that O has strong propensity to induce solidification cracking. Interestingly, Hf is shown to reduce calculated cracking response. Finally, realizing the potential of this method to discover new refractory alloy formulations across the period table that have low solidification cracking sensitivity, the code was applied to new untested alloy systems including W-Zr-C, W-Ta-C, and others. Conclusions. In summary, an open source numerical method has been developed using Python code to calculate Kou’s crack susceptibility index. The method was applied to refractory metals which are inherently difficult to study from a weldability testing standpoint since inert shielding gas is not sufficient and welding is typically done in vacuum, especially in light of findings presented here where oxygen has profound influence on solidification cracking. This work revealed the effect of compositional variations on a series of refractory metals and showed the framework defined here will be useful in 1) the development of new alloys that have improved weldability and 3D printability, 2) placing compositional limits on existing alloys, and 3) ensuring adequate controls of manufacturing processes such as 3D printing where powder reuse is critical. Keywords. pycalphad; Python; refractory metals; solidification cracking. References. B. Bocklund et. al. (2020) http://doi.org/10.5281/zenodo.3630657. S. Kou. (2015) https://doi.org/10.1016/j.actamat.2015.01.034. G.G. Lessmann and R.E. Gold. Welding Journal, issue 1, pp. 1-s – 8-s (1971). F.N. Michael and J.W. Sowards. NASA/TM-20230002218 (2023). R. Otis and Z.-K. Liu. (2017) http://doi.org/10.5334/jors.140. A. Van de Wallle et. al. (2018) https://doi.org/10.1016/j.calphad.2018.04.003.

pycalphad↗

DOE OSTI2018

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↗