Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical stability”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Faster ablative Kelvin–Helmholtz instability growth in a magnetic field

Shear flows along a plasma interface will quickly grow unstable due to the Kelvin–Helmholtz instability. If there is a concurrent temperature gradient across the interface, higher modes are stabilized by the thermal diffusion. These ablative effects must be considered in, for example, jet features in inertial confinement fusion hot-spots, or plasma plumes in young supernovae. We show that magnetization of the plasma can greatly affect the instability, even if magnetic pressure is small. This is because electrons are localized by their gyromotion, reducing the heat flux and material ablation. We use a two-dimensional numerical extended-magnetohydrodynamics approach to assess this effect for dense fusion conditions. In comparison with the unmagnetized case, self-generated Biermann fields make only a minor difference to growth rates. However, simulations with a large 50 kT external field found that the growth rate of the least stable mode increased by 40%. In conclusion, this has implications for mix processes in Z-pinches and magnetized inertial confinement fusion concepts.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Efficient Simulation of Logical Magic State Preparation Protocols

Developing space- and time-efficient logical magic state preparation (MSP) protocols will likely be an essential step toward building a large-scale fault-tolerant quantum computer. Motivated by this need, we introduce a scalable method for simulating logical MSP protocols under the standard circuit-level noise model. When applied to protocols based on code-switching, magic state cultivation, and magic state distillation, our method yields a complexity polynomial in (i) the number of qubits and (ii) the nonstabilizerness, e.g., stabilizer rank or Pauli rank, of the target encoded magic state. The efficiency of our simulation method is rooted in a curious fact: every circuit-level Pauli error in these protocols propagates to a Clifford error at the end. This property is satisfied by a large family of protocols, including those that repeatedly measure a transversal Clifford that squares to a Pauli. We provide a proof-of-principle numerical simulation that prepares a magic state using such logical Clifford measurements. Our work enables practical simulation of logical MSP protocols without resorting to approximations or resource-intensive state-vector simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

New Safety Feedback Control Design to Guarantee Adequate Frequency Performance in Microgrids

Safety analysis of power systems is concerned with the system's ability to maintain critical variables within specified limits following a disturbance. Frequency control adequacy has become increasingly important as the system inertia decreases due to the increase in renewable energy penetration. Various controllers for inverters have been proposed to improve the system frequency response and few are capable to ensure the safety of the response. In this article, a diesel-wind energy system is considered and modeled as a switching system between normal, faulted, and post-fault modes. Further, a safety feedback controller is designed as a supplementary signal for a wind turbine generator such that the speed of the diesel generator stays within a permissible range in the presence of a finite energy disturbance. Numerical results on the modified 33-bus microgrid system obtained of the proposed novel approach indicate that the suggested control configuration can guarantee adequate frequency response without excessive conservativeness.

barrier function↗
DOE OSTI2025

Operation of helium sub-atmospheric multistage cryogenic centrifugal compressor trains: Part 2 – Transient modeling and pump-down path selection

Low-pressure conditions required for operation of helium cryogenic systems below the normal boiling point (i.e. 4.2 K) are established through a transient process, commonly referred to as ‘pump-down’. This process is defined as the transition from pressures above atmospheric conditions to the saturation pressure which corresponds to a specified operational temperature. The FRIB 2 K system consists of five cryogenic centrifugal compressors which are operated in series. Historically, the pump-down process path has been established through empirical methods and system operator experience. Investigation into the pump-down process at FRIB aimed to develop a pump-down methodology which relies on theoretical model predictions rather than empirically developed process paths. Ensuring stable operation during the pump-down process involved application of a centrifugal compressor performance prediction model, which is described in Part 1 of this paper. Compressor performance maps can be directly used to evaluate the stability of a selected pump-down path and anticipate the overall reliability of the selected path. In conjunction with the compressor performance maps, a system pressure model was developed to estimate the transient pressure response during the pump-down process. Lastly, an explicit equation was developed to establish a mass flow rate profile for the pump-down process. Implementation of the presented methodology (including the developed models) allows for the system operator to determine a continuous pump-down path which maintains compressor stability while conforming to overall system capabilities. Altogether, the methodology presented has resulted in simplification of transient pump-down operations and increased the reliability, stability and efficiency of the pump-down process.

Compressor train control↗
DOE OSTI2021

Identification of opportunities for integrating chemical processes for carbon (dioxide) utilization to nuclear power plants

The energy and industrial sectors are responsible for more than 75% of global carbon dioxide emissions, which are primary contributors to climate change. Decarbonization strategies to reduce overall carbon dioxide emissions are urgently needed. Among the numerous decarbonization strategies that are being defined and implemented, a particular strategy of interest for the present work is “the Re-X or zero waste” strategy, which involves recycling, reuse, repurposing, and remanufacturing of industrial emissions, by-products, and waste in general. A circular value chain or Circular Economy is embedded within this strategy. In a circular carbon economy, efficient reuse of emitted carbon dioxide is considered. However, the stability of carbon dioxide implies the requirement of significant amounts of energy for its transformation into value-added chemicals or products, which must be supplied by low carbon emitting energy sources. Nuclear power plants are low carbon energy sources that, additionally to electricity, could also supply heat and radiation for chemical transformations. This work identifies the opportunities and challenges for the development of integrated energy systems to upgrade and transform carbon dioxide, involving chemical and nuclear energy. The analysis encompasses possible use of the different forms of energy that can be obtain from nuclear reactors (i.e., radiation, electricity, and heat), by reviewing the published literature on potential routes for its conversion. Overall, our review indicates that an universal technology of the one-fits-all solution-type is an utopic dream. Instead, a suite of (contaminants tolerant) technologies for processing various concentrations CO 2 and producing a variety of products is currently needed.

08 HYDROGEN↗
DOE OSTI2021

Initial conditions for simulations of beam physics in linear induction accelerators

Flash radiography of hydrodynamic experiments driven by high explosives is a well-known diagnostic technique in use at many laboratories. At Los Alamos, the DualAxis Radiographic Hydrodynamic Test (DARHT) facility two linear induction electron accelerators (LIAs) make the bremsstrahlung radiographic source spots for point projection radiographs from orthogonal views. A new LIA, called Scorpius, is presently under development to advance this technology. To better understand electron-beam physics in these LIAs, numerical simulations are frequently performed with the objective of improving the radiography. At Los Alamos we frequently use the TRAK ray-trace and LSP particle-in-cell (PIC) codes to simulate the injector, and the XTR and LAMDA envelope/centroid codes along with LSP to simulate transport of the accelerated beam through the LIAs. The LIA simulations need the injected beam parameters as initial conditions for calculating beam transport and stability. The determination of these initial conditions is the topic of this note.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Correct Initial Conditions for Simulations of Beam Physics in Linear Induction Accelerators

Flash radiography of hydrodynamic experiments driven by high explosives is a well-known diagnostic technique in use at many laboratories. At Los Alamos, the Dual Axis Radiographic Hydrodynamic Test (DARHT) facility two linear induction electron accelerators (LIAs) make the bremsstrahlung radiographic source spots for point projection radiographs from orthogonal views. A new LIA, called Scorpius, is presently under development to advance this technology. To better understand electron-beam physics in these LIAs, numerical simulations are frequently performed with the objective of improving the radiography. At Los Alamos we frequently use the TRAK ray-trace and LSP particle-in-cell (PIC) codes to simulate the injector, and the XTR and LAMDA envelope/centroid codes along with LSP to simulate transport of the accelerated beam through the LIAs. The LIA simulations need the injected beam parameters as initial conditions for calculating beam transport and stability. The determination of these initial conditions is the topic of this note.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

A Phase-Field Model for In-Space Manufacturing of Binary Alloys

The integrity of the final printed components is mostly dictated by the adhesion between the particles and phases that form upon solidification, which is a major problem in printing metallic parts using available In-Space Manufacturing (ISM) technologies based on the Fused Deposition Modeling (FDM) methodology. Understanding the melting/solidification process helps increase particle adherence and allows to produce components with greater mechanical integrity. We developed a phase-field model of solidification for binary alloys. The phase-field approach is unique in capturing the microstructure with computationally tractable costs. The developed phase-field model of solidification of binary alloys satisfies the stability conditions at all temperatures. The suggested model is tuned for Ni-Cu alloy feedstocks. We derived the Ginzburg-Landau equations governing the phase transformation kinetics and solved them analytically for the dilute solution. We calculated the concentration profile as a function of interface velocity for a one-dimensional steady-state diffuse interface neglecting elasticity and obtained the partition coefficient, k, as a function of interface velocity. Numerical simulations for the diluted solution are used to study the interface velocity as a function of undercooling for the classic sharp interface model, partitionless solidification, and thin interface.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Potent Antiviral Activity against HSV-1 and SARS-CoV-2 by Antimicrobial Peptoids

Viral infections, such as those caused by Herpes Simplex Virus-1 (HSV-1) and SARS-CoV-2, affect millions of people each year. However, there are few antiviral drugs that can effectively treat these infections. The standard approach in the development of antiviral drugs involves the identification of a unique viral target, followed by the design of an agent that addresses that target. Antimicrobial peptides (AMPs) represent a novel source of potential antiviral drugs. AMPs have been shown to inactivate numerous different enveloped viruses through the disruption of their viral envelopes. However, the clinical development of AMPs as antimicrobial therapeutics has been hampered by a number of factors, especially their enzymatically labile structure as peptides. We have examined the antiviral potential of peptoid mimics of AMPs (sequence-specific N-substituted glycine oligomers). These peptoids have the distinct advantage of being insensitive to proteases, and also exhibit increased bioavailability and stability. Our results demonstrate that several peptoids exhibit potent in vitro antiviral activity against both HSV-1 and SARS-CoV-2 when incubated prior to infection. In other words, they have a direct effect on the viral structure, which appears to render the viral particles non-infective. Visualization by cryo-EM shows viral envelope disruption similar to what has been observed with AMP activity against other viruses. Furthermore, we observed no cytotoxicity against primary cultures of oral epithelial cells. These results suggest a common or biomimetic mechanism, possibly due to the differences between the phospholipid head group makeup of viral envelopes and host cell membranes, thus underscoring the potential of this class of molecules as safe and effective broad-spectrum antiviral agents. We discuss how and why differing molecular features between 10 peptoid candidates may affect both antiviral activity and selectivity.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Heterogeneous Nucleating Agent for High‐Boiling‐Point Nonhalogenated Solvent‐Processed Organic Solar Cells and Modules

Abstract High‐boiling‐point nonhalogenated solvents are superior solvents to produce large‐area organic solar cells (OSCs) in industry because of their wide processing window and low toxicity; while, these solvents with slow evaporation kinetics will lead excessive aggregation of state‐of‐the‐art small molecule acceptors (e.g. L8‐BO), delivering serious efficiency losses. Here, a heterogeneous nucleating agent strategy is developed by grafting oligo (ethylene glycol) side‐chains on L8‐BO (BTO‐BO). The formation energy of the obtained BTO‐BO; while, changing from liquid in a solvent to a crystalline phase, is lower than that of L8‐BO irrespective of the solvent type. When BTO‐BO is added as the third component into the active layer (e.g. PM6:L8‐BO), it easily assembles to form numerous seed crystals, which serve as nucleation sites to trigger heterogeneous nucleation and increase nucleation density of L8‐BO through strong hydrogen bonding interactions even in high‐boiling‐point nonhalogenated solvents. Therefore, it can effectively suppress excessive aggregation during growth, achieving ideal phase‐separation active layer with small domain sizes and high crystallinity. The resultant toluene‐processed OSCs exhibit a record power conversion efficiency (PCE) of 19.42% (certificated 19.12%) with excellent operational stability. The strategy also has superior advantages in large‐scale devices, showing a 15.03‐cm 2 module with a record PCE of 16.35% (certificated 15.97%).

Chemistry↗
DOE OSTI2022

Entropy–Preserving and Entropy–Stable Relaxation IMEX and Multirate Time–Stepping Methods

In this work, we propose entropy-preserving and entropy-stable partitioned Runge–Kutta (RK) methods. In particular, we extend the explicit relaxation Runge–Kutta methods to IMEX–RK methods and a class of explicit second-order multirate methods for stiff problems arising from scale-separable or grid-induced stiffness in a system. The proposed approaches not only mitigate system stiffness but also fully support entropy-preserving and entropy-stability properties at a discrete level. The key idea of the relaxation approach is to adjust the step completion with a relaxation parameter so that the time-adjusted solution satisfies the entropy condition at a discrete level. The relaxation parameter is computed by solving a scalar nonlinear equation at each timestep in general; however, as for a quadratic entropy function, we theoretically derive the explicit form of the relaxation parameter and numerically confirm that the relaxation parameter works the Burgers equation. Several numerical results for ordinary differential equations and the Burgers equation are presented to demonstrate the entropy-conserving/stable behavior of these methods. We also compare the relaxation approach and the incremental direction technique for the Burgers equation with and without a limiter in the presence of shocks.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Formation of 6H-Ba 3 Ce 0.75 Mn 2.25 O 9 during Thermochemical Reduction of 12R-Ba 4 CeMn 3 O 12 : Identification of a Polytype in the Ba(Ce,Mn)O 3 Family

The resurgence of interest in a hydrogen economy and the development of hydrogen-related technologies has initiated numerous research and development efforts aimed at making the generation, storage, and transportation of hydrogen more efficient and affordable. Solar thermochemical hydrogen production (STCH) is a process that potentially exhibits numerous benefits such as high reaction efficiencies, tunable thermodynamics, and continued performance over extended cycling. Although CeO 2 has been the de facto standard STCH material for many years, more recently 12R-Ba 4 CeMn 3 O 12 (BCM) has demonstrated enhanced hydrogen production at intermediate H 2 /H 2 O conditions compared to CeO 2 , making it a contender for large-scale hydrogen production. However, the thermo-reduction stability of 12R-BCM dictates the oxygen partial pressure (pO 2 ) and temperature conditions optimal for cycling. In this study, we identify the formation of a 6H-BCM polytype at high temperature and reducing conditions, experimentally and computationally, as a mechanism and pathway for 12R-BCM decomposition. 12R-BCM was synthesized with high purity and then controllably reduced using thermogravimetric analysis (TGA). Synchrotron X-ray diffraction (XRD) data is used to identify the formation of a 6H-Ba 3 Ce 0.75 Mn 2.25 O 9 (6H-BCM) polytype that is formed at 1350 degrees C under strongly reducing pO 2 . Density functional theory (DFT) total energy and defect calculations show a window of thermodynamic stability for the 6H-polytype consistent with the XRD results. These data provide the first evidence of the 6H-BCM polytype and could provide a mechanistic explanation for the superior water-splitting behaviors of 12R-BCM.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

A patterned phase-changing vanadium dioxide film stacking with VO 2 nanoparticle matrix for high performance energy-efficient smart window applications

A vanadium dioxide (VO 2 ) based solid-to-solid phased changing material has been attracting great interest in smart window applications. However, achieving high solar modulation and high transparency simultaneously in visible light is the major challenge for the practical application of this smart material. To resolve this issue, in this paper, a smart film composed of a VO 2 nanoparticle matrix and a patterned VO 2 film is presented. Furthermore, numerical modeling and electromagnetic simulation are carried out to characterize the performance in terms of solar modulation and luminous transmittance, and a parametric study is carried out to optimize the proposed smart window film. Compared with the VO 2 nanoparticle matrix, the proposed structure can obtain 23% solar modulation and 57% luminous transmittance but with a much thinner thickness, which will significantly reduce the cost and fabrication complexity and extend the environment stability.

42 ENGINEERING↗
DOE OSTI2022

Revisiting the link between magnetic properties and chemisorption at graphene nanoribbon zigzag edge

Graphene has received tremendous interest in both chemical and physical fields. Among different edges of the graphene system, the zigzag edge terminated graphene nanoribbons (ZGNRs) show unique magnetic properties in the antiferromagnetic (AFM) state. However, to date, the understanding of ZGNR chemical properties is mainly based on the partial radical concept, and in previous studies, the energy differences between the ferromagnetic (FM) and AFM states are smaller than experimental evidence. Here, we report that the strongly constrained and appropriately normed functional gives a significantly larger energy difference, which matches the experimental observation. Furthermore, utilizing the energetics in the large difference case, we propose a conceptual supplement to the previous partial radical concept: the overall stabilization of the AFM state compared to the nonmagnetic (NM) state consists of two parts that affect the adsorption energy conversely. The NM-FM energy differences will strengthen the adsorption, being in line with the previous partial radical concept. The FM-AFM energy differences will instead weaken the adsorption. Here, we perform calculations of H, OH, and LiS radical adsorption energies on ZGNRs to show that this weakening effect is numerically non-negligible: at least a ~0.2 eV difference in the adsorption energies is found. We expect that this refinement of the partial radical concept can provide a more comprehensive understanding of the chemical properties of ZGNRs. The differences in adsorption energies for the H, OH, and LiS radicals found here lead to significant changes in the predicted reactivity of the ZGNR models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Structure of native four-repeat satellite III sequence with non-canonical base interactions

Abstract Tandem-repetitive DNA (where two or more DNA bases are repeated numerous times) can adopt non-canonical secondary structures. Many of these structures are implicated in important biological processes. Human Satellite III (HSat3) is enriched for tandem repeats of the sequence ATGGA and is located in pericentromeric heterochromatin in many human chromosomes. Here, we investigate the secondary structure of the four-repeat HSat3 sequence 5′-ATGGA ATGGA ATGGA ATGGA-3′ using X-ray crystallography, NMR, and biophysical methods. Circular dichroism spectroscopy, thermal stability, native PAGE, and analytical ultracentrifugation indicate that this sequence folds into a monomolecular hairpin with non-canonical base pairing and B-DNA characteristics at concentrations below 0.9 mM. NMR studies at 0.05–0.5 mM indicate that the hairpin is likely folded-over into a compact structure with high dynamics. Crystallographic studies at 2.5 mM reveal an antiparallel self-complementary duplex with the same base pairing as in the hairpin, extended into an infinite polymer. The non-canonical base pairing includes a G–G intercalation sandwiched by sheared A–G base pairs, leading to a cross-strand four guanine stack, so called guanine zipper. The guanine zippers are spaced throughout the structure by A–T/T–A base pairs. Our findings lend further insight into recurring structural motifs associated with the HSat3 and their potential biological functions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Understanding electronic peculiarities in tetragonal FeSe as local structural symmetry breaking

Traditional band theory of perfect crystalline solids often uses as input the structure deduced from diffraction experiments; when modeled by the minimal unit cell this often produces a spatially averaged model. The present study illustrates that this is not always a safe practice unless one examines if the intrinsic bonding mechanism is capable of benefiting from the formation of a distribution of lower symmetry local environments that differ from the macroscopically averaged structure. This can happen either due to positional or to magnetic symmetry breaking. By removing the constraint of a small crystallographic cell, the energy minimization in the density functional theory finds atomic and spin symmetry breaking, not evident in conventional diffraction experiments but being found by local probes such as atomic pair distribution function analysis. Here in this paper we report that large atomic and electronic anomalies in bulk tetragonal FeSe emerge from the existence of distributions of local positional and magnetic moment motifs. The found symmetry-broken motifs obtained by minimization of the internal energy represent what chemical bonding in the tetragonal phase prefers as intrinsic energy lowering (stabilizing) static distortions. This explains observations of band renormalization, predicts orbital order and enhanced nematicity, and provides unprecedented close agreement with spectral function measured by photoemission and local atomic environment revealed by the pair distribution function. While the symmetry-restricted strong correlation approach has been argued previously to be the exclusive theory needed for describing the main peculiarities of FeSe, we show here that the symmetry-broken mean-field approach addresses numerous aspects of the problem, provides intuitive insight into the electronic structure, and opens the door for large-scale mean-field calculations for similar d -electron quantum materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Accelerating the Discovery of New, Single Phase High Entropy Ceramics via Active Learning

High-entropy ceramics have garnered interest due to their remarkable hardness, compressive strength, thermal stability, and fracture toughness; yet the discovery of new high-entropy ceramics (out of a tremendous number of possible elemental permutations) still largely requires costly, inefficient, trial-and-error experimental and computational approaches. The entropy forming ability (EFA) factor was recently proposed as a computational descriptor that positively correlates with the likelihood that a 5-metal high-entropy carbide (HECs) will form the desired single phase, homogeneous solid solution; however, discovery of new compositions is computationally expensive. If you consider 8 candidate metals, the HEC EFA approach uses 49 optimizations for each of the 56 unique 5-metal carbides, requiring a total of 2744 costly density functional theory calculations. Here, we describe an orders-of-magnitude more efficient active learning (AL) approach for identifying novel HECs. To begin, we compared numerous methods for generating composition-based feature vectors (e.g., magpie and mat2vec), deployed an ensemble of machine learning (ML) models to generate an average and distribution of predictions, and then utilized the distribution as an uncertainty. Here we then deployed an AL approach to extract new training data points where the ensemble of ML models predicted a high EFA value or was uncertain of the prediction. Our approach has the combined benefit of decreasing the amount of training data required to reach acceptable prediction qualities and biases the predictions toward identifying HECs with the desired high EFA values, which are tentatively correlated with the formation of single phase HECs. Using this approach, we increased the number of 5-metal carbides screened from 56 to 15,504, revealing 4 compositions with record-high EFA values that were previously unreported in the literature. Our AL framework is also generalizable and could be modified to rationally predict optimized candidate materials/combinations with a wide range of desired properties (e.g., mechanical stability, thermal conductivity).

36 MATERIALS SCIENCE↗
DOE OSTI2023

MAS: A versatile Landau-fluid eigenvalue code for plasma stability analysis in general geometry

We have developed a new global eigenvalue code, Multiscale Analysis for plasma Stabilities (MAS), for studying plasma problems with wave toroidal mode number (n) and frequency (ω) in a broad range of interest in general tokamak geometry, based on a five-field Landau-fluid description of thermal plasmas. Beyond keeping the necessary plasma fluid response, we further retain the important kinetic effects including diamagnetic drift, ion finite Larmor radius, finite parallel electric field (E||), ion and electron Landau resonances in a self-consistent and non-perturbative manner without sacrificing the attractive efficiency in computation. The physical capabilities of the code are evaluated and examined in the aspects of both theory and simulation. In theory, the comprehensive Landau-fluid model implemented in MAS can be reduced to the well-known ideal MHD model, electrostatic ion-fluid model, and drift-kinetic model in various limits, which clearly delineates the physics validity regime. In simulation, MAS has been well benchmarked with theory and other gyrokinetic and kinetic-MHD hybrid codes in a manner of adopting the unified physical and numerical framework, which covers the kinetic Alfv\'en wave (KAW), ion sound wave (ISW), low-n kink, high-n ion temperature gradient mode (ITG) and kinetic ballooning mode (KBM). Moreover, MAS is successfully applied to model the Alfv\'en eigenmode (AE) activities in DIII-D discharge #159243, which faithfully captures the frequency sweeping of reversed shear Alfv\'en eigenmode (RSAE), the tunneling damping of toroidal Alfv\'en eigenmode (TAE), as well as the polarization characteristics of kinetic beta-induced Alfv\'en eigenmode (KBAE) and beta-induced Alfv\'en-acoustic eigenmode (BAAE) being consistent with former gyrokinetic theory and simulation. With respect to the key progress contributed to the community, MAS has the advantage of combining rich physics ingredients, realistic global geometry and high computation efficiency together for plasma stability analysis in linear regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗