Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “field performance”

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 523 records · Page 29

Automation of Active Space Selection for Multireference Methods via Machine Learning on Chemical Bond Dissociation

Predicting and understanding the chemical bond is one of the major challenges of computational quantum chemistry. Kohn–Sham density functional theory (KS-DFT) is the most common method, but approximate density functionals may not be able to describe systems where multiple electronic configurations are equally important. Multiconfigurational wave functions, on the other hand, can provide a detailed understanding of the electronic structures and chemical bonds of such systems. In the complete active space self-consistent field (CASSCF) method, one performs a full configuration interaction calculation in an active space consisting of active electrons and active orbitals. However, CASSCF and its variants require the selection of these active spaces. This choice is not black box; it requires significant experience and testing by the user, and thus active space methods are not considered particularly user-friendly and are employed only by a minority of quantum chemists. Our goal is to popularize these methods by making it easier to make good active space choices. Here, we present a machine learning protocol that performs an automated selection of active spaces for chemical bond dissociation calculations of main group diatomic molecules. The protocol shows high prediction performance for a given target system as long as a properly correlated system is chosen for training. Good active spaces are correctly predicted with a considerably better success rate than random guess (larger than 80% precision for most systems studied). Our automated machine learning protocol shows that a “black-box” mode is possible for facilitating and accelerating the large-scale calculations on multireference systems where single-reference methods such as KS-DFT cannot be applied.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling Equilibrium Solid–Liquid Interfaces under Effective Constant Chemical Potential Using Machine Learning Interatomic Potentials

The chemical potential (μ) of species in solution is essential for understanding various chemical processes at interfaces. Molecular dynamics (MD) simulations, constrained by fixed compositions, cannot maintain constant chemical potential with reference to a targeted concentration or chemical potential under nonequilibrium or dynamic conditions, as solute species can migrate to the interface and deplete (or enrich) the bulk due to solute-interface interactions. In this study, we introduce a simple and computationally efficient approach named iterative quasi-constant chemical potential molecular dynamics (iqCμMD) simulation, which helps simulate targeted molar concentrations of species in solution. iqCμMD overcomes the limitations of conventional MD by adjusting the number of species in the solution to reach a target bulk concentration (chemical potential), which allows simulation of the interface under the bulk conditions comparable to experiment. We demonstrate our approach using machine learning interatomic potential (MLIP)-based MD simulations of the Na 2 SO 4,aq –graphene interface, and to show the transferability of our approach, we also perform classical force field-based MD simulations of NaCl aq –air and NaCl aq –graphite interfaces, which produce comparable results to previous CμMD simulations. Our results also show that the iqCμMD approach efficiently achieves the desired bulk ion concentration within two iterations, and by utilizing MLIPs, we can achieve converged results using relatively small-scale simulations compared to previous CμMD simulations. By combining iqCμMD with MLIP-driven simulations, solid–liquid interfaces can be modeled under an effective constant chemical potential with DFT-level accuracy. Here, we show that iqCμMD offers a robust and simple computational framework for constant chemical potential simulations, as its only requirement is to be able to converge interfacial simulations with a measurable bulk region.

Chemical structure↗

Spherical compression of an applied magnetic field in inertial confinement fusion

Applying an external magnetic field to laser-driven inertial confinement fusion implosions is a promising approach for enhancing fusion yield. The field is compressed with the plasma, producing a magnetized hotspot that anisotropically suppresses thermal losses and traps alpha particles, making performance sensitive to the compressed field orientation. Here, we derive a simple, readily applicable analytic model that enables rapid evaluation of the compressed field topology and show that ablation into the hotspot amplifies the central field, while the ablated ice near the hotspot edge develops a decaying, radially bent field, with a discontinuity in the field direction. The radially bent field renders thermal insulation at the hotspot edge negligible and largely independent of the applied field strength, whereas insulation in the hotspot core still depends strongly on the applied field. Applying the model to non-axial initial field configurations, we find that an initially applied mirror field provides the greatest suppression, followed by the standard axial field.

Physics - Plasma physics↗

A Code-Agnostic Driver Application for Coupled Neutronics and Thermal-Hydraulic Simulations

While the literature has numerous examples of Monte Carlo and computational fluid dynamics (CFD) coupling, most are hard-wired codes intended primarily for research rather than as standalone, general-purpose applications. In this work, we describe an open source application, ENRICO, that enables coupled neutronic and thermal-hydraulic simulations between multiple codes that can be chosen at runtime (as opposed to a coupling between two specific codes). The application has been designed such that the control flow logic, domain mapping, nonlinear fixed-point iteration, solution transfers, and convergence checks are all agnostic to the underlying physics solvers used. Special emphasis has also been placed on enabling efficient execution on distributed-memory computing environments. The transfer of solution fields between solvers is performed in memory rather than through filesystem I/O. Additionally, solvers can be configured to run on overlapping or disjoint sets of processes. To date, coupling with the OpenMC and Shift Monte Carlo codes, the Nek5000 CFD code, and a simplified heat diffusion and subchannel solver has been implemented in ENRICO. We present results for coupled simulations of a single light-water reactor fuel assembly based on the NuScale reactor using various combinations of the physics solvers. For this problem, the coupled simulations are shown to converge in about four Picard iterations. A comparison of the heat source and temperature distributions computed by ENRICO using OpenMC coupled with Nek5000 and Shift coupled with Nek5000 illustrates remarkable agreement between the codes.

42 ENGINEERING↗

Detection of Diversion in a Realistic Heat Pipe Microreactor Using Supervised Machine Learning

Microreactors (MRs) pose new challenges for international safeguards. Here, their small size and mass reproducibility make them ideal for deployment in greater numbers and in remote locations, making the job of safeguards inspectors more challenging. Machine learning (ML) is currently being applied to many fields to augment human performance and increase automation; in particular, ML could be used to provide insight for international inspectors to help detect the diversion of nuclear fuel from MR cores. Four ML model types (k-nearest neighbors, decision tree, random forest, and histogram-based gradient boosted ensemble) were trained on integrated flux and critical control drum angle data generated with Serpent 2 for a realistic heat pipe MR design, achieving nearly 100% binary classification accuracy of nominal and diversion core configurations by the end of 1 full power year for three of the four model types. Regression model variants were also trained, using the same input data, for predicting the number of fuel pins diverted. Root-mean-square errors below 5% of the total number of fuel pins were achieved by the 1 full power year mark for all models.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Trapped flux in a small crystal of CaKFe 4 As 4 at ambient pressure and in a diamond anvil pressure cell

In an extension of our previous work, the measurements of temperature dependent magnetization associated with trapped magnetic flux in a small single crystal of CaKFe 4 As 4 , using zero - field - cooled and field - cooled protocols were performed, on the same crystal, at ambient pressure without a pressure cell and at 2.2 GPa in a commercial diamond anvil cell (DAC), showing comparable results. The data show that with a proper care and understanding, trapped flux measurements in superconductors indeed can be performed on samples in DACs under pressure, as was done on superhydrides.

36 MATERIALS SCIENCE↗

Behavior and mechanisms of Doppler wind lidar error in complex terrain: stable flow case study at Perdigão

A numerical experiment is carried out investigating the magnitude of biases in ground-based lidar measurements in complex flow conditions. Biases assessed include those arising from flow curvature and from the interaction of turbulence with the wind field reconstruction (WFR) algorithms used by a WindCube lidars and anemometers. RANS-CFD and WRF-LES simulations were performed for the Perdig˜ao Field Experiment site for a range of atmospheric conditions. Virtual anemometer and lidar data were generated for four locations: two near exposed ridge tops and two in low-speed regions in the valley. The LES data at these four locations show that the scalar inflation terms (the relation between scalar and vector averaged wind speed) for virtual lidar and virtual cups agree very well with predictions using perturbation theory. While the lidar errors vary greatly with location and height, the contribution from the flow curvature tends to be larger than the differences arising from scalar inflation. For one lidar/mast pair near the ridge top, comparisons between simulations and measurements are carried out for a resonant mountain wave event on June 14th, 2017, and for the whole duration of the Perdigão campaign for winds perpendicular to the ridges. The lidar error during the mountain wave, a period of strong stability and low inversion height, is significantly larger than the campaign average. The sensitivity of the lidar error to atmospheric stability is confirmed by the RANS simulations, which suggests strong sensitivity of flow curvature error to stability conditions and to the shape of the wind speed profile near the top of the boundary layer.

17 WIND ENERGY↗

Monte Carlo simulation methodology for characterizing the energy response function and resolving power of a digital Flat Panel Detector (FPD) assembly

The use of high-resolution, digital flat panel detectors (FPDs) makes quantitative assessments of radiographic and 3D tomographic representations feasible. However, the energy response function and spatial resolving power of the detector must be known. Here, we present a methodology for characterizing an FPD assembly using Monte Carlo simulations. The detector assembly modeled here was fielded in the 3-Ring Flash X-ray Imaging System designed and built at Lawrence Livermore National Laboratory (LLNL), which employs fifteen Carestream ® HPX-DR3543 FPDs enclosed in protective housings. Detailed Monte Carlo (MCNP6) simulations of an as-fielded detector assembly were performed to estimate its energy response function in terms of analog-to-digital unit (ADU) “counts” per pixel per incident x-ray as a function of energy up to ∼ 500 keV and its theoretical (idealized) resolving power. The practical resolving power of the FPD assembly in a realistic scenario was then estimated by simulating radiographic images of a variant of the United Kingdom (U.K.) “Kaleidoscope” resolution target.

Physics↗

Magnetism and spin dynamics in room-temperature van der Waals magnet Fe 5 GeTe 2

Two-dimensional van der Waals (vdWs) materials have gathered a lot of attention recently. However, the majority of these materials have Curie temperatures that are well below room temperature, making it challenging to incorporate them into device applications. In this work, we synthesized a room-temperature vdW magnetic crystal Fe 5 GeTe 2 with a Curie temperature T c = 332 K, and studied its magnetic properties by vibrating sample magnetometry (VSM) and broadband ferromagnetic resonance (FMR) spectroscopy. The experiments were performed with external magnetic fields applied along the c -axis (H∥ c ) and the ab -plane (H∥ ab ), with temperatures ranging from 300 to 10 K. We have found a sizable Landé g -factor difference between the H∥ c and H∥ ab cases. In both cases, the Landé g -factor values deviated from g = 2. This indicates contribution of orbital angular momentum to the magnetic moment. The FMR measurements reveal that Fe 5 GeTe 2 has a damping constant comparable to Permalloy. With reducing temperature, the linewidth was broadened. Together with the VSM data, our measurements indicate that Fe 5 GeTe 2 transitions from ferromagnetic to ferrimagnetic at lower temperatures. Our experiments highlight key information regarding the magnetic state and spin scattering processes in Fe 5 GeTe 2 , which promote the understanding of magnetism in Fe 5 GeTe 2 , leading to implementations of Fe 5 GeTe 2 based room-temperature spintronic devices.

36 MATERIALS SCIENCE↗

Atomistic modeling of lanthanide diffusion in refractory body-centered cubic molybdenum

Lanthanide fission products can strongly interact with candidate cladding alloys, but their transport properties in refractory metals remain poorly understood. Here, in this work, we investigate the atomic-scale diffusion behavior of La, Ce, Pr, and Nd in body-centered cubic (bcc) molybdenum, a potential candidate for advanced nuclear cladding. Self-consistent mean-field transport modeling is performed to evaluate the fission product transport and vacancy mobility, informed by first-principles and nudged elastic band calculations of vacancy formation energies, migration barriers, and solute–vacancy binding characteristics. Compared with bcc Fe, lanthanide solutes in bcc Mo exhibit slower tracer diffusion due to higher vacancy formation and migration energies. Furthermore, the calculations reveal that the influence of fission products on migration barriers in bcc Mo are not as extensive in range compared to bcc Fe. Among the studied lanthanides, La exhibits the strongest vacancy binding while also being the fastest diffuser in Mo. These findings highlight how refractory bcc alloys can reduce fission product infiltration, offering valuable insight into the development of durable cladding systems for advanced reactors.

36 - MATERIALS SCIENCE↗

Effect of localization on photoluminescence and zero-field splitting of silicon color centers

The study of defect centers in silicon has been recently reinvigorated by their potential applications in optical quantum information processing. A number of silicon defect centers emit single photons in the telecommunication O-band, making them promising building blocks for quantum networks between computing nodes. The two-carbon G-center, self-interstitial W-center, and spin-1/2 T-center are the most intensively studied silicon defect centers, yet despite this, there is no consensus on the precise configurations of defect atoms in these centers, and their electronic structures remain ambiguous. Here, in this work, we employ ab initio density functional theory to characterize these defect centers, providing insight into the relaxed structures, band structures, and photoluminescence spectra, which are compared to experimental results. Motivation is provided for how these properties are intimately related to the localization of electronic states in the defect centers. In particular, we present the calculation of the zero-field splitting for the excited triplet state of the G-center defect as the structure is linearly interpolated from the A-configuration to the B-configuration, showing a sudden increase in the magnitude of the D zz component of the zero-field-splitting tensor. By performing projections onto the local orbital states of the defect, we analyze this transition in terms of the symmetry and bonding character of the G-center defect, which sheds light on its potential application as a spin-photon interface.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Topological and magnetic phase transitions in the bilayer Kitaev-Ising model

We investigate the phase diagram of a bilayer Kitaev honeycomb model with Ising interlayer interactions, deriving effective models via perturbation theory and performing Majorana mean-field theory calculations. We show that a diverse array of magnetic and topological phase transitions occur, depending on the direction of the interlayer Ising interaction and the relative sign of Kitaev interactions. When two layers have the same sign of the Kitaev interaction, a first-order transition from a Kitaev spin liquid to a magnetically ordered state takes place. The magnetic order points along the Ising axis and it is (anti)ferromagnetic for (anti)ferromagnetic Kitaev interactions. However, when two layers have opposite signs of the Kitaev interaction, we observe a notable weakening of magnetic ordering tendencies and the Kitaev spin liquid survives up to a remarkably larger interlayer exchange. Our mean-field analysis suggests the emergence of an intermediate gapped ℤ 2 spin-liquid state, which eventually becomes unstable upon vison condensation. The confined phase is described by a highly frustrated 120° compass model. We furthermore use perturbation theory to study the model with the Ising axis pointing along the $\hat{z}$ axis or lying in the 𝑥⁢𝑦 plane. In both cases, our analysis reveals the formation of one-dimensional Ising chains, which remain decoupled in perturbation theory, resulting in a subextensive ground-state degeneracy. Our results highlight the interplay between topological order and magnetic ordering tendencies in bilayer quantum spin liquids.

Kitaev model↗

From Trivial Kondo Insulator Ce 3 Pt 3 Bi 4 to Topological Nodal-Line Semimetal Ce 3 Pd 3 Bi 4

In this work, using the density functional theory combined with dynamical mean-field theory, we have performed systematic study of the electronic structure and its band topology properties of Ce 3 Pt 3 Bi 4 and Ce 3 Pd 3 Bi 4 . At high temperatures ( ~ 290 K ), the electronic structures of both compounds resemble the open-core 4 f density functional calculation results. For Ce 3 Pt 3 Bi 4 , clear hybridization gap can be observed below 72 K, and its coherent momentum-resolved spectral function below 18 K exhibits an topologically trivial indirect gap of ~ 6 meV and resembles density functional band structure with itinerant 4 f state. For Ce 3 Pd 3 Bi 4 , no clear hybridization gap can be observed down to 4 K, and its momentum-resolved spectral function resembles electron-doped open-core 4 f density functional calculations. The band nodal points of Ce 3 Pd 3 Bi 4 at 4 K are protected by the gliding-mirror symmetry and form ringlike structure. Therefore, the Ce 3 Pt 3 Bi 4 compound is topologically trivial Kondo insulator while the Ce 3 Pd 3 Bi 4 compound is topological nodal-line semimetal.

36 MATERIALS SCIENCE↗

Role of the isovector spin-orbit potential in mitigating the CREX-PREX dilemma

Pioneering electroweak measurements of the neutron skin thickness in lead-208 and calcium-48 are challenging our understanding of nuclear dynamics. Many theoretical models suggest that the slope of the symmetry energy controls the development of a neutron skin in neutron-rich nuclei. This led to the expectation that if lead-208 exhibits a large neutron skin, calcium-48 should as well. Given that the PREX Collaboration reported a relatively thick neutron skin in lead, we anticipated that calcium would also have a significant neutron skin. Instead, the CREX Collaboration reported a thin neutron skin in calcium. Although many suggestions have been proposed, the “CREX-PREX dilemma” remains unsolved. Recently, an intriguing scenario has emerged, suggesting that an enhanced isovector spin-orbit interaction could simultaneously account for both results. Following this approach, we performed relativistic mean-field calculations with an increased isovector spin-orbit potential. Our findings indicate that, while this modification significantly affects the structure of calcium-48, it has only a marginal impact on lead-208, thereby bringing the results into better agreement with experiment. Furthermore, the strong enhancement required to mitigate the CREX-PREX dilemma destroys the agreement with a successful spin-orbit phenomenology, primarily by modifying the well-known ordering of spin-orbit partners.

Neutron skin thickness↗

Effects of magnetic and non-magnetic doping on the vortex lattice in MgB 2

Using small-angle neutron scattering we have studied the vortex lattice in superconducting MgB 2 , doped with either manganese or carbon to achieve a similar suppression of the critical temperature. Measurements were performed with the magnetic field applied along the c-axis, where the vortex lattice in pure MgB 2 is known to undergo a field-and temperature-driven 30° rotation transition. For Mn-doping, the vortex lat- tice phase diagram remains qualitatively similar to that of pure MgB 2 , indicating only a modest effect on the vortex-vortex interaction. In contrast, the vortex lattice rotation transition is completely suppressed in the C-doped case, likely due to a change in the electronic structure which affects the two-band/two-gap nature of superconductivity in MgB 2 . The vortex lattice longitudinal correlation length shows the opposite behavior, remaining roughly unchanged between pure and C-doped MgB 2 while it is significantly reduced in the Mn-doped case. However, the extensive vortex lattice metastability and related activated behavior, observed in conjunction with the vortex lattice transition in pure MgB 2 , is also seen in the Mn doped sample. This shows that the vortex lattice disordering is not associated with a substantially increased vortex pinning.

36 MATERIALS SCIENCE↗

Geomagnetically Induced Current Field Test on Large Grid-Connected Power Transformers: Analysis, Model Development, and Simulations

Geomagnetic-induced current (GIC) flow in power grids can cause undesirable effects such as transformer overheating, harmonics, higher reactive power demand, etc. Many simulation models have been developed to study these effects, but real-world verification on modern transformer designs is rare. Here, this paper presents the first long-duration GIC field test in the U.S. performed on high-voltage, grid-connected transformers featuring winding clamps and tie rods instead of conventional tie bars. Field measurements were taken to evaluate GIC effects. These measurements also aided in developing and validating thermal and electromagnetic transient (EMT) models of the transformers. During the test, significant current and voltage distortions were observed along with considerable transformer reactive power losses. Analysis of the field measurements showed that the transformers’ hottest spot was at the inner windings, and their k-factors were close to factory test and software default values. Thermal simulations indicated that the transformers would not violate their thermal limits even for a GIC waveform that peaks at about 200 A/phase. EMT simulations revealed that increased transformer loading may reduce GIC-induced reactive power demand and harmonics in certain scenarios. The study also highlighted potential inaccuracies in using the k-factor method to calculate transformer reactive power losses.

EMTDC↗

LatticeAnalytics: Strut-Level Visualization and Inspection of Additively Manufactured Lattice Structures

Additive manufacturing (AM) is revolutionizing the production of custom components with complex internal geometries, essential for high-performance applications in diverse fields such as medicine and defense. These AM parts optimize strength while minimizing weight by utilizing internal lattice structures consisting of large quantities of small interconnected struts. However, the complexity of these structures, combined with the challenges of using X-ray Computed Tomography (XCT) data, makes validation of part reliability difficult. This ultimately inhibits the development of novel parts for our collaborating material scientists. Here, we introduce LatticeAnalytics, a novel framework specifically designed for visual inspection of defects in these lattice structures. Our framework offers an end-to-end solution that includes the data management of XCT scans, enables remote access for geographically dispersed teams through a web-based dashboard, and incorporates novel visualizations. Our analysis is facilitated by a coarse alignment between the lattice’s nominal model, a spatial graph, and the XCT data. We employ a simple VR-based approach for fast and rough alignment, followed by an offline registration and identification of the struts. With the nodes and struts aligned and identified in the volume, our framework allows querying of subvolumes containing a single strut at multiple resolutions. This avoids computation over the entire lattice and also allow for easy parallelization of down-stream computations, such as strut-specific metrics. To depict a fast overview of the strut quality, we introduce two innovative visual encodings, crucial for our collaborators’ research in creating novel AM parts: the Contour View and the Roughness Map, which depict critical geometrical and surface features of individual struts in standardized two 2D views. We evaluated the integrated system through expert interviews. The feedback confirms the framework’s practicality and its effectiveness in enhancing current inspection workflows. It solves major bottlenecks for our collaborators, ultimately helping them create novel parts with advanced properties.

Miao, Haichao [Lawrence Livermore National Laborat↗

Zinc–air batteries in neutral/near-neutral electrolytes

The need to fully utilize intermittent renewable energy and have energy security, large-scale batteries are necessary. The aqueous zinc–air battery (ZAB) is a promising potential candidate for its safety, low-cost, and theoretical capacity. This research on ZAB mainly focused on alkaline electrolytes. These are favored for their conductivity, but greatly reduce the stability of the battery by zinc corrosion from the hydrogen evolution reaction (HER), dendrite formation, and carbonate formation. To address these issues, recently neutral and near-neutral electrolytes have been applied to aqueous ZAB to suppress HER, dendrite formation, and minimize carbonate formation. These include the use of chloride-based, potassium nitrate, aqueous organic, solid-state electrolytes supplemented with additives to improve the performance and stability. This field is still young with significant opportunities available for research.

25 ENERGY STORAGE↗