Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complex simulations”

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

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

At least 1,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Lefschetz thimble quantum Monte Carlo for spin systems

Monte Carlo simulations are useful tools for modeling quantum systems, but in some cases they suffer from a sign problem, leading to an exponential slow down in their convergence to a value. While solving the sign problem is generically NP hard, many techniques exist for mitigating the sign problem in specific cases; in particular, the technique of deforming the Monte Carlo simulation's plane of integration onto Lefschetz thimbles (complex hypersurfaces of stationary phase) has seen significant success in the context of quantum field theories. We extend this methodology to spin systems by utilizing spin coherent state path integrals to reexpress the spin system's partition function in terms of continuous variables. Using some toy systems, we demonstrate its effectiveness at lessening the sign problem in this setting, despite the fact that the initial mapping to spin coherent states introduces its own sign problem. The standard formulation of the spin coherent path integral is known to make use of uncontrolled approximations; despite this, for large spins they are typically considered to yield accurate results, so it is somewhat surprising that our results show significant systematic errors. Furthermore, possibly of independent interest, our use of Lefschetz thimbles to overcome the intrinsic sign problem in spin coherent state path integral Monte Carlo enables a novel numerical demonstration of a breakdown in the spin coherent path integral.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Quantum criticality using a superconducting quantum processor

Quantum criticality emerges from the collective behavior of many interacting quantum particles, often at the transition between different phases of matter. It is one of the cornerstones of condensed matter physics, which we access on noisy intermediate-scale (NISQ) quantum devices by leveraging a dynamically driven phenomenon. Here we probe the critical properties of the one-dimensional quantum Ising model on a programmable superconducting quantum chip via a Kibble-Zurek process, obtain scaling laws, and estimate critical exponents despite inherent sources of errors on the hardware. In addition, we investigate how the improvement of NISQ computers (more qubits, less noise) will consolidate the computation of those universal physical properties. A one-parameter noise model captures the effect of imperfections and reproduces the experimental data. Its systematic study reveals that the noise, analogously to temperature, induces a new length scale in the system. We introduce and successfully verify modified scaling laws, directly accounting for the noise without any prior knowledge. It makes data analyses for extracting physical properties transparent to noise. By understanding how imperfect quantum hardware modifies the genuine properties of quantum states of matter, we enhance the power of NISQ processors considerably for addressing quantum criticality and potentially other phenomena and algorithms.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Model emulation and closure tests for (3+1)D relativistic heavy-ion collisions

In nuclear and particle physics, reconciling sophisticated simulations with experimental data is vital for understanding complex systems like the Quark Gluon Plasma (QGP) generated in heavy ion collisions. However, computational demands pose challenges, motivating using Gaussian Process emulators for efficient parameter extraction via Bayesian calibration. We conduct a comparative analysis of Gaussian Process emulators in heavy-ion physics to identify the most adept emulator for parameter extraction with minimal uncertainty. Furthermore, our study contributes to advancing computational techniques in heavy-ion physics, enhancing our ability to interpret experimental data and understand QGP properties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Chaos in inhomogeneous neutrino fast flavor instability

In dense neutrino gases, the neutrino-neutrino coherent forward scattering gives rise to a complex flavor oscillation phenomenon not fully incorporated in simulations of neutron star mergers (NSM) and core collapse supernovae (CCSNe). Moreover, it has been proposed to be chaotic, potentially limiting our ability to predict neutrino flavor transformations in simulations. To address this issue, we explore how small flavor perturbations evolve in the nonlinear regime of the neutrino quantum kinetic equation within a narrow centimeter-scale region inside a NSM and a toy neutrino distribution. Our findings reveal that paths in the flavor state space of solutions with similar initial conditions diverge exponentially, exhibiting chaos. This inherent chaos makes the microscopic scales of neutrino flavor transformations unpredictable. However, the domain-averaged neutrino density matrix remains relatively stable, with chaos minimally affecting it. This particular property suggests that domain-averaged quantities remain reliable despite the exponential amplification of errors. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2025

Quantum algorithm to simulate Lindblad master equations

We present a quantum algorithm for simulating a family of Markovian master equations that can be realized through a probabilistic application of unitary channels and state preparation. Our approach employs a second-order product formula for the Lindblad master equation, achieved by decomposing the dynamics into dissipative and Hamiltonian components and replacing the dissipative segments with randomly compiled, easily implementable elements. The sampling approach eliminates the need for ancillary qubits to simulate the dissipation process and reduces the gate complexity in terms of the number of jump operators. We provide a rigorous performance analysis of the algorithm. We also extend the algorithm to time-dependent Lindblad equations, generalize the family of Markovian master equations it can be applied to, and explore applications beyond the Markovian noise model. A new error bound, in terms of the diamond norm, for second-order product formulas for time-dependent Liouvillians is provided that might be of independent interest. Published by the American Physical Society 2025

Borras, Evan (ORCID:000900017709037X)↗
DOE OSTI2022

C6 Hydroxymethyl-Substituted Carbapenem MA-1-206 Inhibits the Major Acinetobacter baumannii Carbapenemase OXA-23 by Impeding Deacylation

Acinetobacter baumannii has become a major nosocomial pathogen, as it is often multidrug-resistant, which results in infections characterized by high mortality rates. The bacterium achieves high levels of resistance to β-lactam antibiotics by producing β-lactamases, enzymes which destroy these valuable agents. Historically, the carbapenem family of b-lactam antibiotics have been the drugs of choice for treating A. baumannii infections. However, their effectiveness has been significantly diminished due to the pathogen’s production of carbapenem-hydrolyzing class D β-lactamases (CHDLs); thus, new antibiotics and inhibitors of these enzymes are urgently needed. Here, we describe a new carbapenem antibiotic, MA-1-206, in which the canonical C6 hydroxyethyl group has been replaced with hydroxymethyl. The antimicrobial susceptibility studies presented here demonstrated that this compound is more potent than meropenem and imipenem against A. baumannii producing OXA-23, the most prevalent CHDL of this pathogen, and also against strains producing the CHDL OXA-24/40 and the class B metallo-β-lactamase VIM-2. Our kinetic and mass spectrometry studies revealed that this drug is a reversible inhibitor of OXA-23, where inhibition takes place through a branched pathway. X-ray crystallographic studies, molecular docking, and molecular dynamics simulations of the OXA-23-MA-1-206 complex show that the C6 hydroxymethyl group forms a hydrogen bond with the carboxylated catalytic lysine of OXA-23, effectively preventing deacylation. These results provide a promising strategy for designing a new generation of CHDL-resistant carbapenems to restore their efficacy against deadly A. baumannii infections. Carbapenem antibiotics are the drugs of choice for treatment of deadly infections caused by Gram-negative bacteria. However, their efficacy is severely compromised by the wide spread of carbapenem-hydrolyzing class D β-lactamases (CHDLs). The importance of this research is the discovery that substitution of the canonical hydroxyethyl group of carbapenems by a hydroxymethyl significantly enhances stability against inactivation by the major CHDL of Acinetobacter baumannii, OXA-23. These results provide a novel strategy for designing next-generation, carbapenemase-stable carbapenems to fight multidrug-resistant infections caused by Gram-negative pathogens

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Nucleon-induced inelastic scattering with statistical strength functions and the ECIS direct reaction code

Modern theoretical descriptions of inelastic scattering make use of multi-step direct reaction approaches together with transition potentials obtained from sophisticated nuclear structure models. Here we demonstrate how the complexity of such calculations can be reduced to permit simpler ones, also using the ECIS code, but providing an almost equally precise alternative to a much more detailed calculation. We have studied the transition form factors within the random phase approximation (RPA), where these are obtained as linear combinations of particle–hole states. At moderate to high excitation energies, where interference effects tend to disappear, we have proposed an independent particle–hole formalism in which particle–hole states are spread in energy with an appropriate strength function obtained from the RPA. The effects of more complex modes, such as 2p–2h ones, are simulated with widths calculated in a semi-classical context. In this work, we verify the validity of our approximations for pre-equilibrium proton-induced reactions on 90 Zr target. Our calculations provide a good description of the reaction data and point toward a simplification of the description of nucleon-induced reactions based on averages of microscopic details of the projectile–target interaction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Data-Driven Model Reduction of Multiphase Flow in a Single-Hole Automotive Injector

Fuel injector design has a substantial influence on the performance and emissions of direct injection engines. To date, large eddy simulations coupled with a single-fluidmixture modeling approach have shown great success in evaluating the complex interplay among injector design, fuel properties, and operating conditions on the injector performance. However, this simulation approach is too computationally expensive to be used by industry routinely for injector design due to the fine temporal and spatial resolution required to resolve wall-bounded flow within the injector. The work presented in this paper highlights a potential pathway to addressing this issue. To study the influence of injector design, fuel properties, and operating conditions on injector performance, large eddy simulations were performed to model the turbulent multiphase flow development through a side-oriented singlehole diesel injector. Using Latin hypercube sampling, the design space spanning a range of fuel properties, operating conditions, and needle lifts were explored. Two techniques for dimensionality reduction, namely proper orthogonal decomposition and autoencoders, were compared to evaluate their accuracy in representing the flow in a reduced dimensional space. Finally, based on the findings from this work, recommendations are provided in using machine learning approaches within the context of emulation to construct reduced-ordermodels for internal flow development relevant to automotive applications.

42 ENGINEERING↗
DOE OSTI2024

The LANL Common Modeling Framework (CMF) [Slides]

The Common Modeling Framework is a modern software repository built and maintained by the LANL Weapons Physics Directorate (ALDX). ALDX’s mission space is to develop and apply cutting-edge theory, computational models, and multi-physics simulation codes, and to design, execute, and analyze complex experiments

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Design Load Basis Guidance for Distributed Wind Turbines

Aeroelastic modeling (AM) is the primary methodology for structural and performance assessment of any wind turbine. Nonetheless, the use of AM in the distributed wind (DW) industry sector is limited due to several challenges (Damiani, Davis, & Summerville, 2022). One of these challenges lies in the perceived complexity of generating a proper set of numerical simulations to extract and process the key outputs for component design and verification, and, ultimately, achieve certification. This makes it difficult to reliably predict the structural and performance response of small wind turbines. From the investigation carried out in (Damiani & Davis, 2022), it is apparent that many stakeholders in this sector believe that a comprehensive guide for developing a design load basis (DLB) for distributed wind turbines (DWTs) is necessary.

17 WIND ENERGY↗
DOE OSTI2023

Advances in Computational Approaches for Estimating Passive Permeability in Drug Discovery

Passive permeation of cellular membranes is a key feature of many therapeutics. The relevance of passive permeability spans all biological systems as they all employ biomembranes for compartmentalization. A variety of computational techniques are currently utilized and under active development to facilitate the characterization of passive permeability. These methods include lipophilicity relations, molecular dynamics simulations, and machine learning, which vary in accuracy, complexity, and computational cost. This review briefly introduces the underlying theories, such as the prominent inhomogeneous solubility diffusion model, and covers a number of recent applications. Various machine-learning applications, which have demonstrated good potential for high-volume, data-driven permeability predictions, are also discussed. Due to the confluence of novel computational methods and next-generation exascale computers, we anticipate an exciting future for computationally driven permeability predictions.

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

Using the Monte Carlo Method to Evaluate the Reliability of Screening Multifamily Housing for Radon

When screening for radon in a multifamily housing complex using a fixed sample density (e.g., testing 1 in 10 (10%) or 1 in 4 units (25%)), the statistical confidence is dependent upon the assumed elevated radon frequency. For example, testing 25% of units in a complex estimated to have eight units with elevated radon levels will provide 90% confidence. However, if it is assumed that, in the same complex, there are only three units with elevated radon levels, the confidence drops to around 58% for the same sample density. Furthermore, in the previous example, if elevated radon levels are not found during the screening, all that can be stated is that the screening provides 90% confidence that there are no more than seven units in the complex with elevated radon levels. To more fully illustrate this uncertainty, ten separate multifamily housing radon data sets with 1 to 10 units with radon levels ≥4 pCi/L (Table 1) were selected for analysis using the Monte Carlo statistical method. Unlike other mathematically based statistical approaches, the Monte Carlo statistical method relies on repeated analysis of randomly selected data. from a 100% sampled complex at various sample densities. Success for each simulation is defined as finding at least one unit with elevated radon levels. In this statistical method, confidence in the overarching conclusion can be greatly enhanced by repeating the simulated screening hundreds or even thousands of times. For this study, each of the ten data sets was simulated 1,000 times. For each simulation, the data set was randomized three times before the fixed percentage of data was selected.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2021

Control-Theory-Informed Feature Selection for Detecting Malicious Tampering in Additive Layer Manufacturing Processes

Additive layer manufacturing (ALM) is rapidly becoming an appealing solution to the low-volume manufacturing of metal, polymer, or composite parts. However, ALM’s reliance on digital part specifications, microcontrollers, and modern networking makes these devices vulnerable to malicious tampering by cyber attackers, which can negatively affect part performance and even result in catastrophic failure. We present a hybrid analytic approach to feature discovery using control theoretic techniques and linear modelling on input-output data collected from a representative controller system. Employing this approach, we design, train, and test an anomaly detection system. The preliminary results show that the proposed approach effectively discovers useful input-output relationships for anomaly detection in a simulated ALM process. Application to larger and more complex systems are discussed.

Dawson, Joel↗
DOE OSTI2025

A Refined BISON-Griffin Workflow For Designing Targeted TRISO Transient Experiments In The Transient Reactor Test Facility

U.S. industry is expected to deploy tristructural isotropic (TRISO) particle fuel technologies for commercial reactors within the next decade. In our previous work, we defined a preliminary transient design space for TRISO fuels, identified potential gaps in the available data, and began to develop multiphysics modeling tools that could be applied to design targeted Transient Reactor Test Facility (TREAT) experiments to fill these gaps. This work builds on that foundation by (1) updating BISON fuel performance and Griffin reactor physics models to reflect the current TREAT experiment tube and capsule designs,(2) coupling the codes to improve the accuracy and usability of the transient design analyses, and (3) demonstrating their use over an expanded design space that includes fuel burnup. The simulated mechanical responses of the TRISO particles were complex functions of fission product accumulation, fission gas release, and irradiation-induced dimensional change in the pyrolytic carbon layers. Predicted tangential stresses in particle silicon carbide layers were least compressive for preheated tests involving fresh fuels but remained compressive throughout the ranges of temperature, heat rate, and burnup considered in this work. Finally, comparisons between the potential TREAT transients and historical test reactor irradiations showed that the TREAT tests would produce significantly lower average energy deposition rates, yielding less severe transients with greater relevance to near-term commercial applications. Use of these predictive capabilities has the potential to increase the value of each test, improving the overall efficiency and cost effectiveness of transient testing for TRISO and other advanced fuels.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

5-minute Wind Power Data based on WFIP2 WRF Simulation

The second Wind Forecast Improvement Project (WFIP2) was a public-private partnership funded by the U.S. Department of Energy and NOAA, aimed at enhancing the forecast skill of numerical weather prediction models for turbine-height winds in regions with complex terrain. An 18-month Weather Research and Forecasting (WRF) model simulation was conducted over the Pacific Northwest, with model outputs validated against observational data collected during WFIP2. Simulated wind speeds were used to estimate wind power generation using reV (the Renewable Energy Potential Model developed by NREL) at ten wind project sites. Two sets of results were produced: one using wind speeds extracted from the model grid cell at the project centroid, and another using wind speeds from the actual turbine locations. For each dataset, power output was calculated using both actual turbine-specific power curves and nine generic power curves to convert wind speed into power.

17 WIND ENERGY↗
DOE OSTI2024

Performance of explicit and IMEX MRI multirate methods on complex reactive flow problems within modern parallel adaptive structured grid frameworks

Large-scale multiphysics simulations are computationally challenging due to the coupling of multiple processes with widely disparate time scales. The advent of exascale computing systems exacerbates these challenges since these systems enable ever-increasing size and complexity. In recent years, there has been renewed interest in developing multirate methods as a means to handle the large range of time scales, as these methods may afford greater accuracy and efficiency than more traditional approaches of using implicit-explicit (IMEX) and low-order operator splitting schemes. However, to date there have been few performance studies that compare different classes of multirate integrators on complex application problems. In this work, we study the performance of several newly developed multirate infinitesimal (MRI) methods, implemented in the SUNDIALS solver package, on two reacting flow model problems built on structured mesh frameworks. The first model revisits prior work on a compressible reacting flow problem with complex chemistry that is implemented using BoxLib but where we now include comparisons between a new explicit MRI scheme with the multirate spectral deferred correction (SDC) methods in the original paper. The second problem uses the same complex chemistry as the first problem, combined with a simplified flow model, but runs at a large spatial scale where explicit methods become infeasible due to stability constraints. Two recently developed IMEX MRI multirate methods are tested. These methods rely on advanced features of the AMReX framework on which the model is built, such as multilevel grids and multilevel preconditioners. The results from these two problems show that MRI multirate methods can offer significant performance benefits on complex multiphysics application problems and that these methods may be combined with advanced spatial discretization to compound the advantages of both.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A practical extension of the recursive multi-fidelity model for the emulation of hole closure experiments

We report in regimes of high strain rate, the strength of materials often cannot be measured directly in experiments. Instead, the strength is inferred based on an experimental observable, such as a change in shape, that is matched by simulations supported by a known strength model. In hole closure experiments, the rate and degree to which a central hole in a plate of material closes during a dynamic loading event are used to infer material strength parameters. Due to the complexity of the experiment, many computationally expensive, three-dimensional simulations are necessary to train an emulator for calibration or other analyses. These simulations can be run at multiple grid resolutions, where dense grids are slower but more accurate. In an effort to reduce the computational cost, a combination of simulations with different resolutions can be combined to develop an accurate emulator within a limited training time. We explore the novel design and construction of an appropriate functional recursive multi-fidelity emulator of a strength model for tantalum in hole closure experiments that can be applied to arbitrarily large training data. Hence, by formulating a multi-fidelity model to employ low-fidelity simulations, we were able to reduce the error of our emulator by approximately 81% with only an approximately 1.6% increase in computing resource utilization.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Changes in polarization dictate necessary approximations for modeling electronic deexcitation intensity: Application to x-ray emission

Accurate simulation of electronic excitations and deexcitations are critical for complementing complex spectroscopic experiments and can provide validation to theoretical approaches. Here, using a generalized framework, we contrast the accuracy and validity of orbital-constrained and linear-response approaches that build upon Kohn-Sham density functional theory (DFT) to simulate emission spectra of electronic origin and propose an efficient approximation, named many-body x-ray emission spectroscopy or MBXES, for simulating such processes. We show analytically as well as with computed examples that for electronic (de)excitation leading to an appreciable change in polarization (i.e., density rearrangement), the adiabatic approximation in a response-based formalism will be inadequate for the calculation of oscillator strength. Thus, such a change (e.g., in the net electrostatic dipole moment of a finite system) can be used as a metric for evaluating the applicability of the adiabatic response-based approach and can be particularly valuable in x-ray emission spectroscopy. On the other hand, MBXES, the flexible method introduced in this paper, can compute oscillator strengths accurately at a much lower computational expense on the basis of two DFT-based self-consistent field calculations. Using illustrative examples of emission spectra, the efficacy of the MBXES method is demonstrated by comparison with its parent theory, orbital-optimized DFT, and with experiments.

74 ATOMIC AND MOLECULAR PHYSICS↗