Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Optimal Zone”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

TaZr_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys tantalum-zirconium (Ta-Zr). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach \cite{PAW}. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Ta and Zr. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on NERSC-Perlmutter and OLCF-Summit using the VASP 6.3.2. The VASP calculations for every atomic structure have been performed in 2 main steps: 1. Starting from an ideal body-centered-cubic (BCC) structure, geometry optimization with low precision has been executed to perform a preliminary optimization of the atomic structure. The output for this calculations is available in the files 0.CONTCAR, 0.OUTCAR, rlx1.out. 2. Using the atomic structure resulting from the preliminary geometry optimization, a second geometry optimization has been performed using normal precision. The output for this calculations is available in the files CONTCAR, OUTCAR, rlx2.out, vaspout.h5, and vasprun.xml. Cases 1-10 have been run without generating the file 'vaspout.h5'. Every chemical composition sampled across the composition range in the dataset has its own directory. The convention used to name the directories for binary alloys is AXBY, where A and B refer to the constituents, whereas X and Y are positive integers that represent the number of atoms for each constituent and their values still sum up to 128. Each atomic structure associated with a specific chemical composition has its own sub-directory within the directory of the corresponding chemical composition. The sub-directories for each atomic structure for each chemical composition are named 'case-*', where * is a positive integer that spans all the values from 1 through 100, extremes included. The files contained in each sub-directory 'case-*' for each atomic structure are as follows: FILES contained in each subdirectory with name "case-N" where N ranges between 11 and 80, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. 0.POSCAR: input file that defines the atomic structure of a system 4. 0.CONTCAR: output file that provides the atomic positions and cell parameters after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 5. 0.OUTCAR: output file that contains detailed information about the progress of a calculation after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.out: file with diagnostic information about the execution of the first geometry optimization with precision variable set to PREC=Low in the INCAR file 7. POSCAR: input file that defines the atomic structure of a system after the first geometry optimization has been run at low precision. This represents the input for the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. CONTCAR: output file that provides the atomic positions and cell parameters after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 9. OUTCAR: output file that contains detailed information about the progress of a calculation after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.out: file with diagnostic information about the execution of the second geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system FILES contained in each subdirectory with name "case-N" where N ranges between 1 and 10 and between 81 and 100, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. {ID}.POSCAR: input file that defines the atomic structure of a system at the beginning of ID execution of geometry optimization with PREC=LOW 4. {ID}.CONTCAR: output file that provides the atomic positions and cell parameters at the end of ID execution of geometry optimization with PREC=LOW in the INCAR file 5. {ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of geometry optimization that has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.{ID}.out: file with diagnostic information about the execution of the ID execution of the geometry optimization with precision variable set to PREC=Low in the INCAR file 7. N{ID}.POSCAR: input file that defines the atomic structure of a system after the geometry optimization run at low precision. This represents the input for the ID execution of the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. N{ID}.CONTCAR: output file that provides the atomic positions and cell parameters after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 9. N{ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.{ID}.out: file with diagnostic information about the ID execution of geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system This research is sponsored by the Artificial Intelligence Initiative as part of the Laboratory Directed Research and Development (LDRD) Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the US Department of Energy under contract DE-AC05-00OR22725. This work used resources of the Oak Ridge Leadership Computing Facility, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC05-00OR22725, under Directorate Discretionary awards MAT025 (Materials Science) and LRN026 (Machine Learning), and INCITE award MAT201. This work also used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231, under award ERCAP0025216. REFERENCES (1) Kresse, G. & Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. & Hafner, J. Ab initio molecular-dynamics simulation of the liquid-metal–amorphous-semiconductor transition in germanium. Phys. Rev. B 49, 14251 (1994) (3) Kresse, G. & Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. materials science 6, 15–50 (1996) (4) Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. review B 54, 11169 (1996) (5) Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗
DOE OSTI2025

Hund's coupling assisted orbital-selective superconductivity in Ba1-xKxFe2As2

While the superconducting transition temperature of hole-doped Ba_{1-x}K_{x}Fe_{2}As_{2} decreases past optimal doping, superconductivity does not completely disappear even for the fully doped KFe_{2}As_{2} compound. In fact, superconductivity is robust through a Lifshitz transition where electron bands become hole-like around the zone corner at around x=0.7, thus challenging the conventional understanding of superconductivity in iron-based systems. High-resolution angle-resolved photoemission spectroscopy is used to investigate the superconducting gap structure, as well as the normal state electronic structure, around optimal doping and across the Lifshitz transition. Our findings reveal a largely orbital-dependent superconducting gap structure, where the more strongly correlated d_{xy} band has a vanishing superconducting gap at higher doping, aligning with the Hund's metal behavior observed in the normal state. Notably, the superconducting gap on the d_{xy} band disappears before the Lifshitz transition, suggesting that the Fermi surface topology may play a secondary role. We discuss how these results point to orbital-selective superconducting pairing and how strong correlations via Hund's coupling may shape superconducting gap structures in iron-based and other multiorbital superconductors.

FOS: Physical sciences↗
DOE OSTI2023

A General Method for Estimating Zonal Transmission Interface Limits from Nodal Network Data: Preprint

Capacity expansion models for the electric power system often employ zonal (rather than nodal) resolution, necessitating estimates of aggregate power transfer limits across the interfaces between model zones. Interface limits between planning areas are sometimes published, but they are not generalizable to arbitrary zone shapes. There is thus a need for a reproducible method for estimating interface transfer limits (ITLs) between user-defined zones directly from nodal transmission system data. Here, we present a simple method for estimating ITLs using a DC power flow approximation via the power transfer distribution factor (PTDF) matrix. Linear optimization is performed to identify the distribution of power flows that maximizes the total flow on interface-crossing lines, subject to individual line ratings, limits on bus injection/withdrawal, and the relationships among flows, injections, and withdrawals imposed by the PTDF matrix. We demonstrate the application of the method on a 134-zone ~65000-bus system, and we explore the influence of flow direction, contingency level, and zone size on the estimated ITLs. There is significant heterogeneity in the ratio of the ITL to the sum of interface-crossing line ratings, which highlights the importance of accounting for the physical constraints on power flows imposed by Kirchhoff's laws when estimating zonal ITLs.

17 WIND ENERGY↗
DOE OSTI2020

Rapid generation of optimal generalized Monkhorst-Pack grids

In this report computational modeling of the properties of crystalline materials has become an increasingly important aspect of materials research, consuming hundreds of millions of CPU-hours at scientific computing centers around the world each year, if not more. A routine operation in such calculations is the evaluation of integrals over the Brillouin zone. We have previously demonstrated that performing such integrals using generalized Monkhorst-Pack k-point grids can roughly double the speed of these calculations relative to the widely-used traditional Monkhorst-Pack grids. However the generation of optimal generalized Monkhorst-Pack grids is not implemented in most software packages due to the computational cost and difficulty of identifying the best grids. To address this problem, we present new algorithms that allow rapid generation of optimal generalized Monkhorst-Pack grids on the fly. We demonstrate that the grids generated by these algorithms are on average significantly more efficient than those generated using existing algorithms across a range of grid densities. For grids that correspond to a real-space supercell with at least 50 Å between lattice points, which is sufficient to converge density functional theory calculations within 1 meV/atom for nearly all materials, our algorithm finds optimized grids in an average of 0.19 s on a single processing core. To facilitate the widespread adoption of this approach, we present new open-source tools including a library designed for integration with third-party software packages.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

On The Incorporation of Both Function-Driven Requirements and Topology Optimization in The Development of Lunar Launch & Landing Pads

In the design of additively manufactured structures, there can be a constant battle between focusing on function-driven design and topology optimization in an attempt to improve a products ability to meet all requirements. Function-driven optimization is defined as the process of identifying the best design parameters that satisfy project functionality requirements. Function-driven optimization is typically implemented by a designer(s) who, instead of computational limitations, has human bias for the input or inputs determined. On the other hand, topology optimization is defined as a mathematical method that optimizes material layout within a given design space, for a given set of loads, boundary conditions, and constraints with the goal of maximizing the performance of the system. In topology optimization, the work is defined by the limited mathematical inputs and computationally available goals/restraints. If not evaluated simultaneously, these two approaches can result in fundamentally different design concepts meant to address the same requirements. NASA’s Artemis program includes the development of a base camp near the lunar South Pole. After site selection is made, one of the first infrastructure elements needed will be a launch and landing pad (LLP). The construction of such a pad will prevent damage to nearby structures from plume-surface interaction (PSI) while also providing a stable platform for the lander. These benefits reduce the required distance between the established landing zones and other settlement infrastructure. In this paper, the authors evaluate the requirements associated with lunar LLP development, application of both functional and topological optimization techniques, methodologies to combine the two techniques, and potential impacts on final design of such a pad.

In-situ↗
DOE OSTI2024

Sensitization and Mechanical Response of Cu‐Containing Steel Rods

The iron and steel manufacturing sector significantly adds to global greenhouse gas emissions, caused primarily by the carbothermic reduction of iron ore. Recycling scrap steel offers an effective decarbonization strategy but introduces impurities like copper (Cu) that can negatively impact mechanical properties. This study investigates the effects of Cu content and heat treatment on the mechanical performance and sensitization of steel wire rods for tire manufacturing. Steel rods with 0.04 and 0.21 wt% Cu are heated to 1050 or 1200 °C, then air quenched, or furnace cooled. Tensile testing coupled with microscopic analysis is used to evaluate mechanical properties and assess the sensitization effects. Higher Cu content leads to larger sensitized zones with increased Cu precipitation along grain boundaries. Ductility and toughness, crucial for wire drawability, are found to be reduced, despite higher ultimate strength. Slower furnace cooling is seen to result in smaller sensitized zones compared to air quenching, suggesting a pivotal role of cooling rate in sensitization control. The findings provide insights into optimize heat treatment parameters and Cu content limits, balancing mechanical performance and maintaining drawability for enhanced scrap steel recycling in tire production.

copper↗
DOE OSTI2025

Characterizing and Modeling the Influence of Geometry on the Performance of Superconducting Nanowire Cryotrons

The scaling of superconducting nanowire detectors to larger arrays is often limited by room-temperature-readout cabling. Cryogenic integrated circuits constructed from nanowire cryotrons, or nanocryotrons, can address this limitation by performing signal processing on chip. In this study, we characterize key performance metrics of the nanocryotron to elucidate its potential as a logical element in cryogenic integrated circuits and develop an electro-thermal model to connect material parameters with device performance. We find that the performance of the nanocryotron depends on the device geometry, and trade-offs are associated with optimizing the gain, jitter, and energy dissipation. Here, we demonstrate that nanocryotrons fabricated on niobium nitride can achieve a grey zone less than 210 nA wide for a 5 ns long input pulse corresponding to a maximum achievable gain of 48 dB, an energy dissipation of less than 20 aJ per operation, and a jitter of less than 60 ps.

Superconductor↗
NASA NTRS1990

Experimental approaches for improving SiGe thermoelectric efficiency at JPL

Efforts to improve the thermoelectric efficiency of silicon-germanium alloys are discussed. These efforts are essentially focused on substantially enhancing carrier concentration levels and thus increasing the power factor values of n-type materials. With the help of thermodynamic phase diagram calculations, a combination of various crystal growth and dopant saturation techniques is used for this purpose. Bulk polycrystal growth of double-doped materials is achieved by zone leveling while a new method (the traveling heater method) for high-quality single-crystal growth of SiGe is being developed. Thin films of homogeneous single-crystalline multidoped Si80Ge20 have also been successfully grown from low-temperature melts to investigate combined dopant solid solubility variations. Very high carrier concentration levels were obtained by (P+As)-combined diffusion in zone-leveled and hot-pressed materials. The good overall understanding of doping mechanisms which has been achieved should enable optimization of the thermoelectric properties over this extended carrier concentration range. Reductions in thermal conductivity are being sought by introduction of very small phonon scattering centers (inert material) in the initial powders-precursors of Si-Ge alloys. This approach should benefit both n- and p-type materials.

Borshchevsky, Alex↗
NASA NTRS1977

Implications of outer-zone radiations on operations in the geostationary region utilizing the AE4 environmental model

The radiation exposure in the region of geostationary orbits is examined in search for means of optimizing human performance. It is found that the use of slightly inclined circular orbits is one means by which exposure and spacesuit thickness requirements can be reduced. Another effective technique is to limit the extravehicular activity to those days when the short term fluctuations result in low exposure. Space-suit shielding approaching 1/2 sq cm or less may be possible by utilizing work stoppages and inclined orbits. If aluminum and other low-atomic-number materials are used to construct the habitat, then excessive wall thicknesses are required. If special bremsstrahlung shielding is used, then the habitat shield may be reduced to as low as 2 g/sq cm. Numerous tables and graphs are presented for future analysis of dose in the geostationary region.

Wilson, J. W.↗
DOE OSTI2024

In-vivo Raman microspectroscopy reveals differential nitrate concentration in different developmental zones in Arabidopsis roots

Abstract Background Nitrate (NO 3 − ) is one of the two major forms of inorganic nitrogen absorbed by plant roots, and the tissue nitrate concentration in roots is considered important for optimizing developmental programs. Technologies to quantify the expression levels of nitrate transporters and assimilating enzymes at the cellular level have improved drastically in the past decade. However, a technological gap remains for detecting nitrate at a high spatial resolution. Using extraction-based methods, it is challenging to reliably estimate nitrate concentration from a small volume of cells (i.e., with high spatial resolution), since targeting a small or specific group of cells is physically difficult. Alternatively, nitrate detection with microelectrodes offers subcellular resolution with high cell specificity, but this method has some limitations on cell accessibility and detection speed. Finally, optical nitrate biosensors have very good ( in-vivo ) sensitivity (below 1 mM) and cellular-level spatial resolution, but require plant transformation, limiting their applicability. In this work, we apply Raman microspectroscopy for high-dynamic range in-vivo mapping of nitrate in different developmental zones of Arabidopsis thaliana roots in-situ . Results As a proof of concept, we have used Raman microspectroscopy for in-vivo mapping of nitrate content in roots of Arabidopsis seedlings grown on agar media with different nitrate concentrations. Our results revealed that the root nitrate concentration increases gradually from the meristematic zone (~ 250 µm from the root cap) to the maturation zone (~ 3 mm from the root cap) in roots grown under typical growth conditions used for Arabidopsis, a trend that has not been previously reported. This trend was observed for plants grown in agar media with different nitrate concentrations (0.5–10 mM). These results were validated through destructive measurement of nitrate concentration. Conclusions We present a methodology based on Raman microspectroscopy for in-vivo label-free mapping of nitrate within small root tissue volumes in Arabidopsis. Measurements are done in-situ without additional sample preparation. Our measurements revealed nitrate concentration changes from lower to higher concentration from tip to mature root tissue. Accumulation of nitrate in the maturation zone tissue shows a saturation behavior. The presented Raman-based approach allows for in-situ non-destructive measurements of Raman-active compounds.

Fernández González, Alma↗
DOE OSTI2020

Optimizing PCM-integrated walls for potential energy savings in U.S. Buildings

Buildings in the United States account for nearly half of total U.S. energy use. The energy used for space conditioning can be reduced by utilizing thermal energy storage, such as phase change materials (PCMs), into building envelopes; however, the energy savings of PCM-integrated building envelopes reported in the literature vary widely. In the absence of established guidelines, thermophysical requirements of an optimal PCM, its method of application into the building envelope, and the corresponding energy savings under various climates remain unknown. Here, we perform an extensive numerical investigation on the integration of PCM into building walls to establish the key conditions required for effective utilization of PCM in reducing heat gains in the cooling season and heat losses in the heating season. We also determine the optimal transition temperature, optimal PCM location in the wall, and the energy-saving potential of the PCM-integrated building walls in five U.S. cities located in different International Energy Conservation Code climate zones. Results show that employing PCMs in building walls does not always lead to an improvement; in fact, incorrect applications of PCMs can substantially increase energy use in the buildings. In the climates we studied, PCMs were found effective in reducing heat gains during the cooling season while mostly ineffective in managing heat losses during heating season. Depending on the climate, optimized PCMs in U.S. building walls can provide reduction in the annual heat gain in the range of 3.5% to 47.2% and the annual heat loss in the range of -2.8% to 8.3%. Future consideration of buildings with substantial solar gains in winter may lead to more reduction in heat losses by PCMs.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1992

Mixing in the dome region of a staged gas turbine combustor

To lower NO(x) emissions from gas-turbine engines the effect of dome design and operational changes on the mixing quality in the fuel-rich region is studied. A statistical analysis is employed to establish the parametric sensitivity in this complex flow. A mixing-effectiveness index is defined and used to optimize the gas-species uniformity and the extent of reaction at the exit plane of the dome. Mixing effectiveness is tied to the fuel and air injection locations, the macroscale structure of the dome aerodynamics, and the level of turbulence. Increases in nozzle/air to fuel ratio, reference velocities, and the dome expansion angle increased the level of turbulence. The optimum configuration featured counter-swirling fuel and air streams and produced a strong torroidal recirculation zone, an effective spray angle of 45 degrees, and azimuthal velocities that decayed to zero inside of two duct diameters. The results underscore the system specific nature of mixing optimization.

Sowa, W. A.↗
DOE OSTI2022

A liquid hydrogen tube moderator arrangement for SNS second target station

Liquid hydrogen filled tubes arranged in a triangular shape surrounded by light-water premoderators were investigated as cold moderators coupled to a neutron production zone of a tungsten target fed by a 1.3 GeV proton beam in a short-pulse mode. A moderator concept optimized for the tube length, premoderator thickness, and target position was found to deliver the highest pulse-integrated neutron brightness suitable to send cold neutron beams for neutron scattering instrumentation geared toward medium resolution and small samples. Reducing the tube diameter from 50 to 10 mm boosts the cold neutron brightness (E < 10 meV) by a factor of 2.25, also with enhancements seen by downsizing the height of cylindrical moderators in earlier studies. Attempts were made to characterize this novel moderator system with regard to pulse shapes, local brightness distribution, and angular distribution of emission at the neutron emission port.

Gallmeier, Franz X. (ORCID:0000000171500764)↗
NASA NTRS1973

Application of ERTS-1 data to the protection and management of New Jersey's coastal environment

ERTS-1 imagery is being used by the New Jersey Department of Environmental Protection (NJDEP) to develop information products that will assist the state in optimally managing its coastal resources and in allocating funds. Interviews with NJDEP personnel have identified significant problem areas in the coastal zone, and the types of remote sensor derived information products that can be used in real-time decision making. Initial analyses of imagery from several successive ERTS-1 orbits have shown the extent, predominant drift, and dispersion characteristics of waste disposal in coastal New Jersey waters. Imagery (MSS Bands 4 and 5) for several orbits, shows that New-York Harbor tidal discharge extending as far south as Long Branch, New Jersey.

Mairs, R. L.↗
DOE OSTI2026

Bent Crystal Channeling for Optimized Beam Shadowing and Proton Extraction at Mu2e

The Mu2e experiment is designed to investigate the CLFV through the observation of a neutrinoless muon-to-electron conversion in the field of an Al nucleus. The observation of such a process would be clear evidence of physics beyond the standard model. Due to the rarity of this process, a cutting-edge, intense muon beam is required to achieve an improvement of the current single-event sensitivity by 4 orders of magnitude. To achieve this goal, a primary proton beam with 8 GeV is extracted from the Fermilab Delivery Ring using the slow resonant extraction technique. Mu2e requires ~3.6x10$^{20}$ protons-on-target to meet its goal; hence, it is crucial to minimize the extraction losses. An important source of such losses are the particles impacting on the electrostatic septum blade. A very promising solution to the problem lies in the beam shadowing scheme tested at CERN SPS. In this approach, a bent crystal is strategically placed upstream of the septum, deflecting particles from the blade at a precise angle via the phenomenon of channeling. As a result, a zone with reduced particle flux is created downstream of the crystal, safeguarding the septum anode by minimizing interactions with the beam. This work explores the optimization of beam shadowing design and the process in the manufacturing and characterization of the bent crystal sample. It emphasizes the promising role of channeling in bent crystals, and it underscores the significant potential of channeling in bent crystals to assist the Mu2e experiment.

Fedeli, Pierluigi [Ferrara U.; INFN, Ferrara] (ORC↗
DOE OSTI2021

Manifold learning for coarse-graining atomistic simulations: Application to amorphous solids

In this work, we introduce a generalized machine learning framework to probabilistically parameterize upper-scale models in the form of nonlinear PDEs consistent with a continuum theory, based on coarse-grained atomistic simulation data of mechanical deformation and flow processes. The proposed framework utilizes a hypothesized coarse-graining methodology with manifold learning and surrogate-based optimization techniques. Coarse-grained high-dimensional data describing quantities of interest of the multiscale models are projected onto a nonlinear manifold whose geometric and topological structure is exploited for measuring behavioral discrepancies in the form of manifold distances. A surrogate model is constructed using Gaussian process regression to identify a mapping between stochastic parameters and distances. Derivative-free optimization is employed to adaptively identify a unique set of parameters of the upper-scale model capable of rapidly reproducing the system's behavior while maintaining consistency with coarse-grained atomic-level simulations. The proposed method is applied to learn the parameters of the shear transformation zone (STZ) theory of plasticity that describes plastic deformation in amorphous solids as well as coarse-graining parameters needed to translate between atomistic and continuum representations. We show that the methodology is able to successfully link coarse-grained microscale simulations to macroscale observables and achieve a high-level of parity between the models across scales.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Optimising response to an introduction of African swine fever in wild pigs

African swine fever virus (ASFv) is a virulent pathogen that threatens domestic swine industries globally and persists in wild boar populations in some countries. Persistence in wild boar can challenge elimination and prevent disease-free status, making it necessary to address wild swine in proactive response plans. In the United States, invasive wild pigs are abundant and found across a wide range of ecological conditions that could drive different epidemiological dynamics among populations. Information on the size of the control areas required to rapidly eliminate the ASFv in wild pigs and how this area should change with management constraints and local ecology is needed to optimize response planning. We developed a spatially explicit disease transmission model contrasting wild pig movement and contact ecology in two ecosystems in Southeastern United States. We simulated ASFv spread and determined the optimal response area (reported as the radius of a circle) for eliminating ASFv rapidly over a range of detection times (when ASFv was detected relative to the true date of introduction), culling capacities (proportion of wild pigs in the culling zone removed weekly) and wild pig densities. Large radii for response areas (14 km) were needed under most conditions but could be shortened with early detection (≤ 8 weeks) and high culling capacities (≥ 15% weekly). Under most conditions, the ASFv was eliminated in less than 22 weeks using optimal control radii, although ecological conditions with high rates of wild pig movement required higher culling capacities (≥ 10% weekly) for elimination within 1 year. The results highlight the importance of adjusting response plans based on local ecology and show that wild pig movement is a better predictor of the optimal response area than the number of ASFv cases early in the outbreak trajectory. Furthermore, our framework provides a tool for determining optimal control plans in different areas, guiding expectations of response impacts, and planning resources needed for rapid elimination.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

An Optimized Parameterization of Sub‐Grid Scale Advection for Convection Permitting Models

Convection‐permitting models (CPMs) explicitly resolve deep convection yet under‐resolve the organized lateral exchanges among drafts and their environment that control entrainment/detrainment, precipitation efficiency, and mesoscale structure. In this work, we introduce the Optimized Advection Scheme (OAS), which introduces a small rotation of the Cartesian frame of reference for the horizontal winds relative to other variables used in advection that induces cross‐gradient transport to mimic under‐resolved convective mixing. The rotation angle is selected to minimize the Kullback–Leibler divergence between the simulated and satellite observed precipitation intensity distributions, yielding a physically consistent perturbation that is computationally inexpensive and portable. Optimized Advection Scheme is implemented in WRF and evaluated over Amazon (April 2014). It shifts precipitation–precipitable‐water joint distributions toward lighter rain, reduces overly intense rates, and improves mesoscale convective system (MCS) lifetime and propagation. Mechanistically, the added cross‐gradient transport promotes convective detrainment and environmental mixing, which cools and moistens the mid‐troposphere, weakens downward momentum transport, alleviates excessive downwelling shortwave biases, and warms the surface temperature. The optimized rotation angle yields comparable improvements at 4‐km and 1‐km grid spacing, demonstrating resolution‐independent benefits across the CPM gray zone. By targeting the dynamical root of under‐mixed convective circulations, rather than tuning model microphysics or closures, OAS delivers robust, scale‐aware improvements in precipitation statistics, cloud vertical structure, and characteristics of MCS (MCSs), offering a practical pathway to more reliable CPM simulations for weather and climate applications.

CPM↗