Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “M CODES”

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

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

At least 1,153 records · Page 64

Risk-based area of review estimation in overpressured reservoirs to support injection well storage facility permit requirements for CO 2 storage projects

This paper by the Energy & Environmental Research Center presents a workflow and modeling approach for delineating a risk-based area of review (AOR) to support a U.S. Environmental Protection Agency (EPA) Class VI permit for a carbon dioxide (CO 2 ) storage project. The approach combines semianalytical solutions for estimating formation fluid leakage through a hypothetical leaky wellbore with the results of numerical reservoir simulations to define the AOR. The modeling utilizes 1) semianalytical solutions from the peer-reviewed literature for formation fluid leakage through abandoned wellbores by Raven (1990) and Avci (1994), 2) a FORTRAN model compiled and described in Cihan et al. (2011, 2012) called ASLMA (Analytical Solution for Leakage in Multilayered Aquifers), and 3) a computational framework for estimating a risk-based AOR first proposed by Oldenburg et al. (2014, 2016). Therefore, the approach builds upon well-established research and underlying hydrogeological principles that have been upheld for nearly three decades. Moreover, the ASLMA model has been broadly applied to an array of storage projects. The work presented herein extends these earlier works using a custom wrapper written in the software environment, R (R Core Team, 2020), which was developed to perform multiple runs of the ASLMA model using given ranges for one or more input parameters. In addition, the current work simulates the pressure buildup within the storage reservoir in response to CO 2 injection using a compositional simulator to better accommodate the temporospatial evolution of pressure buildup within the storage reservoir that is more accurately modeled using a heterogeneous geologic model and a compositional simulator that accounts for the multiphase interactions. The workflow is demonstrated using a case study for a 180,000-metric-ton-per-year storage project located in the PCOR (Plains CO 2 Reduction) Partnership region. For the storage project evaluated here, under the scenario where the leaky wellbore is open to a saline aquifer (thief zone) between the overlying seal (cap rock) and the underground sources of drinking water (USDW), the risk-based AOR essentially collapses to the areal extent of the CO 2 plume in the storage reservoir because the pressure buildup in the storage reservoir beyond the CO 2 plume is insufficient to drive formation fluids up a hypothetical leaky wellbore into the USDW. However, even under the conservative assumption that the leaky wellbore is not open to a thief zone, beyond the areal extent of the CO 2 plume, the incremental leakage is less than 400 m 3 over 20 years, which represents ~0.0001% or less of the total volume of water contained within the USDW rock volume. As discussed in the text, the threshold criterion for defining the risk-based AOR is site-specific and should be informed by the results of the sensitivity analysis and available site characterization data. The approach outlined in this paper is designed to be protective of USDWs and, therefore, comply with the Safe Drinking Water Act requirements and provisions for the U.S. EPA Class VI Underground Injection Control (UIC) Program (Class VI Rule) and North Dakota Administrative Code Chapter 43-05-01.

54 ENVIRONMENTAL SCIENCES↗

Application and assessment of deep learning for the generation of potential NMDA receptor antagonists

Uncompetitive antagonists of the N-methyl D-aspartate receptor (NMDAR) have demonstrated therapeutic benefit in the treatment of neurological diseases such as Parkinson's and Alzheimer's, but some also cause dissociative effects that have led to the synthesis of illicit drugs. The ability to generate NMDAR antagonists in silico is therefore desirable for both new medication development and preempting and identifying new designer drugs. Recently, generative deep learning models have been applied to de novo drug design as a means to expand the amount of chemical space that can be explored for potential drug-like compounds. In this study, we assess the application of a generative model to the NMDAR to achieve two primary objectives: (i) the creation and release of a comprehensive library of experimentally validated NMDAR phencyclidine (PCP) site antagonists to assist the drug discovery community and (ii) an analysis of both the advantages conferred by applying such generative artificial intelligence models to drug design and the current limitations of the approach. We apply, and provide source code for, a variety of ligand- and structure-based assessment techniques used in standard drug discovery analyses to the deep learning-generated compounds. Overall, we present twelve candidate antagonists that are not available in existing chemical databases to provide an example of what this type of workflow can achieve, though synthesis and experimental validation of these compounds are still required.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Charged-particle transport in high energy density plasmas

This Special Topic Collection grew out of two gatherings of researchers active in the high energy density (HED) physics community: a mini-conference on charged-particle transport in HED plasma held during the 64th annual meeting of the American Physical Society's Division of Plasma Physics (Spokane, WA, November 2022) and a dedicated charged-particle transport coefficient code comparison workshop (Livermore, CA, July 2023). These gatherings provided opportunities for theoretical, computational, and experimental researchers to discuss the state of the field, including current capabilities and methods, needs of hydrodynamic simulations, and frontiers for future research. Finally, this special issue collects a total of 13 research and review articles on charged-particle transport in HED plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Non-local transport in radiation-hydrodynamics codes for ICF by efficient coupling to an external Vlasov–Fokker–Planck code

Accurately incorporating non-local transport into radiation-hydrodynamics codes, and indeed any fluid system, has long been elusive. To date, a simplified and accurate theory that can be easily integrated has not been available. This limitation affects modeling in inertial confinement fusion (ICF) and magnetic confinement fusion systems, among others, where non-local transport is well-known to be present. Here, we present a coupling methodology between a full Vlasov–Fokker–Planck (VFP) electron kinetic code and radiation-hydrodynamics (rad-hydro) codes. The VFP code is used to adjust native electron transport in the rad-hydro code, thus enabling improved transport without the need to integrate a full electron VFP solver into the rad-hydro code. This approach necessitates only occasional invocation of the VFP code, reducing computational intensity compared to following the dynamic evolution entirely with the VFP code on fluid time scales. We illustrate that the methodology is more accurate than other simplified methods in thermal decay systems relevant to ICF and can replicate standard theoretical results with high accuracy.

Electronic transport↗

A Framework for International Collaboration on ITER Using Large-Scale Data Transfer to Enable Near-Real-Time Analysis

The global nature of the ITER project along with its projected ~ petabyte per day data generation presents a unique challenge, but also an opportunity for the fusion community to rethink, optimize, and enhance our scientific discovery process. Recognizing this, collaborative research with computational scientists was undertaken over the past several years to create a framework for large-scale data movement across wide-area networks (WANs), to enable global near-real time analysis of fusion data. This would broaden the available computational resources for analysis/simulation, and increase the number of researchers actively participating in experiments. An official demonstration of this framework for fast, large data transfer and real-time analysis was carried out between the KSTAR tokamak in Daejeon, Korea and PPPL in Princeton, USA. Streaming large data transfer, with near real-time movie creation and analysis of the KSTAR Electron Cyclotron Emission imaging (ECEI) data, was performed using the I/O framework ADIOS, and comparisons made at PPPL with simulation results from the XGC1 code. These demonstrations were made possible utilizing an optimized network configuration at PPPL, which achieved over 8.8 Gbps (88% utilization) in throughput tests from NFRI to PPPL. This demonstration showed the feasibility for large-scale data analysis of KSTAR data, and provides a nascent framework to enable use of globally distributed computational and personnel resources in pursuit of scientific knowledge from the ITER experiment.

43 PARTICLE ACCELERATORS↗

Pedestal particle balance studies in JET-ILW H-mode plasmas

JET-ILW type I ELMy H-modes at 2.5 MA/2.8 T with constant NBI heating (23 MW) and gas fuelling rate were performed, utilising edge localised mode (ELM) pacing by vertical kicks and plasma shaping (triangularity, δ) as tools to disentangle the effects of ELMs, inter-ELM transport and edge stability on the pedestal particle balance. In agreement with previous studies, the pedestal confinement improves with increasing δ, mostly due to a significant increase in pedestal density while the ELM frequency ($f_{\mathrm{ELM}}$) is decreased. Improved pedestal confinement with increasing δ was observed even when the pedestal MHD stability was degraded artificially by vertical kicks, implying that increased triangularity may favourably affect the inter-ELM pedestal recovery. The workflow developed to quantify the pedestal particle balance uses high time-resolution profile reflectometry to characterise the inter-ELM evolution of the plasma particle content ($\mathrm{d} N/\mathrm{d}t$), the NEO drift-kinetic solver to evaluate the neoclassical fluxes and interpretative EDGE2D-EIRENE simulations to estimate the edge particle source. The edge particle source is then constrained by deuterium Balmer-α line intensity measurements in the main chamber, which are, however, strongly affected by reflections from the metal walls. The reflections are accounted for by the CHERAB code taking the divertor emission (the brightest light source in the torus) distribution from imaging spectroscopy measurements as input. Our analysis shows that in the second half of the ELM cycle, the volume-integrated particle source is larger than $\mathrm{d} N/\mathrm{d}t$, indicating that transport plays a key role in the inter-ELM pedestal recovery.

fuelling↗

Scenario optimization for the tokamak ramp-down phase in RAPTOR: Part B. safe termination of DEMO plasmas

An optimized plasma current ramp-down strategy is critical for safe and fast termination of plasma discharges in a tokamak demonstration fusion reactor (DEMO), both in planned and emergency scenarios, avoiding plasma disruptions and excessive heat loads to the first wall. Plasma stability limits and machine-specific technical requirements constrain the stable envelope through which the plasma must be navigated. Large amounts of auxiliary heating are required throughout the ramp-down phase, to avoid a radiative collapse in the presence of intrinsic tungsten and seeded xenon impurities, as quantitatively estimated in this work. As the plasma current is reduced, the current density becomes increasingly peaked, reflected by a growing value of the internal inductance $\ell$ $i3$ , resulting in reduced controllability of the vertical position of the plasma. The feasibility of different plasma current ramp-down rates is tested by applying an automated optimization framework embedding the RAPTOR core transport solver. Optimal time traces for plasma current $I$ $p$ $(t)$ and plasma elongation $κ(t)$ are proposed, to satisfy an $I$ $p$ -dependent upper limit on the plasma internal inductance, as obtained from vertical stability studies using the CREATE-NL code, as well as a constraint on the time evolution of $q$ 95 , to avoid an ideal MHD mode. A negative current density near the plasma edge is observed in our simulations, even for the most conservative $I$ $p$ ramp-down rate, indicating significant transient dynamics due to a large resistive time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Overview of T and D–T results in JET with ITER-like wall

In 2021 JET exploited its unique capabilities to operate with T and D–T fuel with an ITER-like Be/W wall (JET-ILW). This second major JET D–T campaign (DTE2), after DTE1 in 1997, represented the culmination of a series of JET enhancements—new fusion diagnostics, new T injection capabilities, refurbishment of the T plant, increased auxiliary heating, in-vessel calibration of 14 MeV neutron yield monitors—as well as significant advances in plasma theory and modelling in the fusion community. DTE2 was complemented by a sequence of isotope physics campaigns encompassing operation in pure tritium at high T-NBI power. Carefully conducted for safe operation with tritium, the new T and D–T experiments used 1 kg of T (vs 100 g in DTE1), yielding the most fusion reactor relevant D–T plasmas to date and expanding our understanding of isotopes and D–T mixture physics. Furthermore, since the JET T and DTE2 campaigns occurred almost 25 years after the last major D–T tokamak experiment, it was also a strategic goal of the European fusion programme to refresh operational experience of a nuclear tokamak to prepare staff for ITER operation. The key physics results of the JET T and DTE2 experiments, carried out within the EUROfusion JET1 work package, are reported in this paper. Progress in the technological exploitation of JET D–T operations, development and validation of nuclear codes, neutronic tools and techniques for ITER operations carried out by EUROfusion (started within the Horizon 2020 Framework Programme and continuing under the Horizon Europe FP) are reported in (Litaudon et al Nucl. Fusion accepted), while JET experience on T and D–T operations is presented in (King et al Nucl. Fusion submitted).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Lifetimes of two- and four-neutron-hole states in 204 Pb and 206 Pb

The centroid-shift method was used to determine half-lives of excited levels in 204 Pb and 206 Pb. For the states with 𝐼 𝜋 = 12 + , 17 − , and 19 − in 204 Pb, values of 𝑇 1/2 = 4.9⁢(3), 3.9(3), and 0.9(3) ns, respectively, were determined. The half-life of the 𝐼 𝜋 = 9 − level in 206 Pb was measured to be 1.1(3) ns. In 204 Pb, the 12 + state arises from two neutron holes in the 𝑖 13/2 subshell, while the 17 − and 19 − levels are built on 𝜈⁡(𝑖$^{−3}_{13/2}$, 𝑝$^{−1}_{1/2}$) and 𝜈⁡(𝑖$^{−3}_{13/2}$, 𝑓$^{−1}_{5/2}$) four-neutron-hole configurations, respectively. The 9 − state in 206 Pb results from a 𝜈⁡(𝑖$^{−1}_{13/2}$, 𝑓${−1}_{5/2}$) excitation. As a result, the transition probabilities for the decay of these isomeric levels are reasonably consistent with the expected single-particle values, and these states are described well by shell-model calculations using the kshell code with the KHHE interaction.

190 ≤ A ≤ 219↗

Transfer reactions in 206 $\mathrm{Pb}$ + 118 $\mathrm{Sn}$: From quasielastic to deep-inelastic processes

Here, we measured multinucleon transfer reactions for the 206 Pb+ 118 Sn system at E lab =1200 MeV by employing the large solid angle magnetic spectrometer PRISMA. Differential and total cross sections and Q-value distributions have been obtained for a variety of neutron and proton pick-up and stripping channels. The Q-value distributions show how the quasi-elastic and deep inelastic processes depend on the mass and charge of the transfer products. The corresponding cross sections have been compared with calculations performed with the GRAZING code. An overall good agreement is found for most of the few nucleon transfer channels. The underestimation of the data for channels involving a large number of transferred nucleons indicates that more complicated processes populate the given isotopes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Cross sections for the formation of Rb84m,g, Rb83, and Rb82m in Sr86(d,x) reactions up to deuteron energies of 49 MeV: Competition between α-particle and multinucleon emission processes

Cross sections of Sr86(d,x) reactions leading to the products Rb84m,g, Rb83, and Rb82m were measured by the stacked-sample activation technique up to deuteron energies of 49 MeV. Nuclear model calculations were performed using the codes talys and empire, which combine the statistical, precompound, and direct interaction components. In all cases, the empire results were much higher than the talys calculation. Fairly good agreement was obtained between measured data and the talys calculation after some optimization of the input model parameters. Insight into competition between α-particle and multinucleon emission in the Y88 compound-nucleus system was also gained.

59 ≤ A ≤ 89↗

MatRIS: Multi-level Math Library Abstraction for Heterogeneity and Performance Portability using IRIS Runtime

Vendor libraries are tuned for a specific architecture and are not portable to others. Moreover, they lack support for heterogeneity and multi-device orchestration, which is required for efficient use of contemporary HPC and cloud resources. To address these challenges, we introduce MatRIS—a multilevel math library abstraction for scalable and performance-portable sparse/dense BLAS/LAPACK operations using IRIS runtime. The MatRIS-IRIS co-design introduces three levels of abstraction to make the implementation completely architecture agnostic and provide highly productive programming. We demonstrate that MatRIS is portable without any change in source code and can fully utilize multi-device heterogeneous systems by achieving high performance and scalability on Summit, Frontier, and a CADES cloud node equipped with four NVIDIA A100 GPUs and four AMD MI100 GPUs. A detailed performance study is presented in which MatRIS demonstrates multi-device scalability. When compared, MatRIS provides competitive and even better performance than libraries from vendors and other third parties.

Monil, M. A. H.↗

Study of recirculating liquid fuel in a 1D critical stationary system

Several studies have been conducted to investigate the physics of liquid fuel reactors, showing also applications with molten salts. A liquid nuclear fuel implies changes in the neutron balance equation to take into account the precursors' displacement and the emission of delayed neutrons in a different position than at the original fission. This requires to upgrade the computer codes normally used to calculate nuclear reactors using only solid fuel. In this work, we revisit a simple problem with liquid fuel, which is proposed for the verification of the numerical solutions obtained by advanced computer codes. This problem studies criticality with constant coefficients, thus neglecting thermal feedback. We elaborate on the analytical solution of the problem, deriving also a generalized eigenvalue problem by finite-volume integration over the cells of a discretized mesh to study the evolution of the dominance ratio with fuel velocity. Finally, we investigate the influence of the fuel velocity on the reactivity of the system. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A pixel detector ROIC with a 3-bit, 40 Msps synchronous ADC per pixel for HL LHC

The prototype Fermi Compact Muon Solenoid (CMS) Pixel (FCP130), manufactured in a 130 nm process, is a readout integrated circuit (ROIC) designed for the High Luminosity upgrade of the Large Hadron Collider. It employs a synchronous analog to digital converter (ADC) for the front-end design, with signal processing and data conversion within a single bunch crossing of 25 ns. It is therefore capable of accurately detecting hits occurring in consecutive bunch crossings, without off-time registration of events, making it particularly suitable for the inner most layers of the vertex detector. The ROIC consists of a matrix of 48 × 160 pixels, each 100 × 30 μm2 in size. Each pixel contains a charge sensitive preamplifier with self-biasing leakage current compensation, eight auto-zero comparators and a Gray code encoder for a 3-bit flash-type ADC. The total power consumption is approximately 24 μW per pixel. The measured noise at the output of all the hit comparators across the ROIC is < 80 e- with threshold dispersion < 110 e-RMS, which allows an in-time threshold setting of ~ 1000 e-.

43 PARTICLE ACCELERATORS↗

Method of regulating gene expression in a cell

The disclosure provides methods of making a tetracycline inducible expression system in a cell. The methods include providing the cell with a first nucleic acid sequence comprising a first promoter operably linked to a tetracycline repressor gene coding sequence, providing the cell with a second nucleic acid sequence comprising a second promoter operably linked to a coding sequence of a gene of interest wherein the second promoter is modified to include one or more tetracycline repressor protein binding sites, and determining the expression of the gene of interest in the presence or absence of tetracycline. The disclosure further provides nucleic acid sequences, vectors and cells including the tetracycline inducible modified promoter.

59 BASIC BIOLOGICAL SCIENCES↗

Expanded Validation of Uranium Systems with the KENO Monte Carlo Codes and SCALE 6.2.4 [Slides]

There is likely a trend in bias for uranium-fueled systems as a function of spectrum, but probably not as a function of enrichment. Some enrichment trends appear, but are potentially the result of outlier experimental results; further investigation is needed where possible. There are good examples in this data of markedly different biases from different systems. Future work includes review and vetting for inclusion in the VALID library and analysis with TSUNAMI.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Solidification and crystallographic texture modeling of laser powder bed fusion Ti-6Al-4V using finite difference-monte carlo method

Laser powder bed fusion (LPBF) additive manufacturing makes near-net-shaped parts with reduced material cost and time, rising as a promising technology to fabricate Ti-6Al-4V, a widely used titanium alloy in aerospace and medical industries. However, LPBF Ti-6Al-4V parts produced with 67° rotation between layers, a scan strategy commonly used to reduce microstructure and property inhomogeneity, have varying grain morphologies and weak crystallographic textures that change depending on processing parameters. Here, this study predicts LPBF Ti-6Al-4V solidification at three energy levels using a finite difference-Monte Carlo method and validates the simulations with large-area electron backscatter diffraction (EBSD) scans. The developed model accurately shows that a <001> texture forms at low energy and a <111> texture occurs at higher energies parallel to the build direction but with a lower strength than the textures observed from EBSD. A validated and well-established method of combining spatial correlation and general spherical harmonics representation of texture is developed to calculate a difference score between simulations and experiments. The quantitative comparison enables effective fine-tuning of nucleation density (N 0 ) input, which shows a nonlinear relationship with increasing energy level. Future improvements in texture prediction code and a more comprehensive study of N 0 with different energy levels will further advance the optimization of LPBF Ti-6Al-4V components. These developments contribute a novel understanding of crystallographic texture formation in LPBF Ti-6Al-4V, the development of robust model validation and calibration pipeline methodologies, and provide a platform for mechanical property prediction and process parameter optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The General Atomic and Molecular Electronic Structure System (GAMESS): Novel Methods on Novel Architectures

The primary focus of GAMESS over the last 5 years has been the development of new high-performance codes that are able to take effective and efficient advantage of the most advanced computer architectures, both CPU and accelerators. These efforts include employing density fitting and fragmentation methods to reduce the high scaling of well-correlated (e.g., coupled-cluster) methods as well as developing novel codes that can take optimal advantage of graphical processing units and other modern accelerators. Because accurate wave functions can be very complex, an important new functionality in GAMESS is the quasi-atomic orbital analysis, an unbiased approach to the understanding of covalent bonds embedded in the wave function. Finally, best practices for the maintenance and distribution of GAMESS are also discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗