Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Simulation–Extrapolation for Bias Correction with Exposure Uncertainty in Radiation Risk Analysis Utilizing Grouped Data

Summary In observational epidemiological studies, the exposure that is received by an individual often cannot be precisely observed, resulting in measurement error, and a common approach to dealing with measurement error is regression calibration (RC). Use of RC, which requires assumptions about the distribution of unknown error-free (true) variables, leads to concern about the possibility of bias due to misspecification of that distribution. The simulation–extrapolation (SIMEX) method, in contrast, does not require a distributional assumption. However, analyses of large cohorts may be performed by using grouped or person-year data, and application of SIMEX to grouped data is not straightforward, particularly when there is a mixture of classical and Berkson measurement errors. We compared RC and SIMEX with grouped data analyses to assess robustness of the RC method to misspecification of the true dose distribution. We also applied SIMEX assuming mixtures of classical and Berkson errors and compared the results with those obtained by using RC for classical error only. SIMEX had less bias than RC and performed well regardless of the true dose distribution, whereas RC based on a misspecified true dose distribution showed greater bias than when based on the correctly specified true dose distribution.

Misumi, Munechika↗
DOE OSTI2022

Parallel simulated annealing with embedded machine learning and multifidelity models for reactor core design

This paper presents extensions to a penalty-free, parallel simulated annealing (SA) algorithm for multi-constrained combinatorial optimization with the aim of embedding multi-fidelity physics models into the annealing procedure. The method uses a low-fidelity, quickly executing model for rapid design space exploration and a high-fidelity model for detailed constraint resolution and on-the-fly bias correction. Machine learning models updated within the annealing procedure were used to bridge the gap between the multi-fidelity models, which led to accurate rapid exploration and efficient detailed constraint resolution. A software implementation of the new multi-fidelity optimization methods, called ML-PSA, was demonstrated on a continuous multi-fidelity optimization problem and a constrained combinatorial PWR lattice design problem. These problems demonstrate some of the features, parallel performance characteristics, and extensible nature of the multi-fidelity SA methods. This paper shows that the developed software and procedure are a general optimization tool that can be applied to a wide variety of scientific and engineering design optimization applications. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2018

Risk Reduction of CO 2 Storage with Stochastic Simulations

The purpose of this project was to create an efficient and flexible tool based on principal component analysis (PCA) for a generation of differentiable realizations of porosity and permeability fields. Of particular interest were data which had significant connectivity between patches of the same rock type, which is called “binary”. The efficiency requirement needed a so called kernel PCA, and the binary images were found with the method called optimization based PCA (OPCA) for which interpretation is provided. The ability to honor given data at particular locations is also incorporated. The tool worked very well and efficiently for two dimensional (2-D) simulations. Furthermore, the tool was applied to a three-dimensional (3-D) dataset and determined that conditioning can be used to maintain the connectivity between vertical layers.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Simulation of Chilled-water Injection at EGS Collab Testbed 2 using the GEOS simulation framework

This work seeks to model the chilled-water injection that occurred at EGS Collab Testbed 2 as part of experiment 3, using the open-source GEOS simulation framework. Modeling of this process, previously conducted using STOMP-GT, confirmed that thermal breakthrough was expected to occur during the chill water injection and may have been masked by strong Joule-Thompsons effects induced by large pressure drops between the fractures and the production wells. Building a GEOS numerical model of this process will allow to include additional physics. In fact, GEOS can model fully coupled thermo-poromechanical processes in fractured porous media and can be employed to model the process of artificially enhancing rock permeability by hydraulically fracturing the rocks. For example, the GEOS model could help understanding the dynamic nature of the EGS collab experiment 3 flow system. Since the limited time and resources do not allow for an extensive study, the focus is on building a baseline model of this process. As such, the main outcomes of this modelling efforts are listed below.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Simulating Superconducting Devices in Novel Materials Using the G4CMP Simulation Toolkit

Understanding the kinematics of phonons and charge propagation in superconducting devices is crucial for conducting low-threshold dark matter searches and minimizing correlated errors in superconducting qubits. For nearly a decade, the Geant4 Condensed Matter Physics (G4CMP) package, originally developed for the Cryogenic Dark Matter Search (CDMS) experiment, has been limited to simulating charge and phonon transport in silicon and germanium materials. In this work, we have expanded the capabilities of G4CMP to include novel substrate materials such as sapphire (Al2O3), Gallium Arsenide (GaAs), Lithium Fluoride (LiF), Calcium Tungstate (CaWO4), and Calcium Fluoride (CaF2). We demonstrate the use of this toolkit in generating phonon transport properties for these materials and compare the results with experimentally determined values where available. Additionally, we explore the performance of superconducting devices utilizing these materials.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Simulation of 15% and 50% Thermal Power Dispatch to an Industrial Facility Using a Flexible Generic Full-Scope Pressurized Water Reactor Plant Simulator

Nuclear power plants in the United States are increasingly challenged to compete in wholesale electricity markets due to the low electricity costs and increasingly dynamic grid conditions from competing generation sources. An alternative market for nuclear power is industrial facilities that can use the thermal and/or electrical power generated by a nuclear power plant to offset the economic losses incurred by electricity market challenges. A generic pressurized water reactor (PWR) simulator was used to show the results of a basic design for a generic thermal power extraction system and tests were run using a set of procedures to show what happens when a nuclear power plant transitions from full electrical power dispatch to 15% and 50% thermal power dispatch. This type of operation leads to losses in turbine performance efficiency due to the deviation from the design operating point, but because the thermal power is also used by the industry load without conversion losses, the combined thermal efficiency of the PWR increases. For the 15% case, the thermal efficiency increased from 32% to 41.9%, while for the 50% case, the efficiency increased up to 60.1%. In addition, for 50% thermal power dispatch, the power dissipated by the condenser decreased from approximately 2000 to approximately 1300 MW (thermal), indicating a substantially diminished impact on the environment in terms of releasing heat into the cooling water reservoir.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Molecular Dynamics Free Energy Simulations Reveal the Mechanism for the Antiviral Resistance of the M66I HIV-1 Capsid Mutation

While drug resistance mutations can often be attributed to the loss of direct or solvent-mediated protein−ligand interactions in the drug-mutant complex, in this study we show that a resistance mutation for the picomolar HIV-1 capsid (CA)-targeting antiviral (GS-6207) is mainly due to the free energy cost of the drug-induced protein side chain reorganization in the mutant protein. Among several mutations, M66I causes the most suppression of the GS-6207 antiviral activity (up to ~84,000-fold), and only 83- and 68-fold reductions for PF74 and ZW-1261, respectively. To understand the molecular basis of this drug resistance, we conducted molecular dynamics free energy simulations to study the structures, energetics, and conformational free energy landscapes involved in the inhibitors binding at the interface of two CA monomers. To minimize the protein−ligand steric clash, the I66 side chain in the M66I−GS-6207 complex switches to a higher free energy conformation from the one adopted in the apo M66I. In contrast, the binding of GS-6207 to the wild-type CA does not lead to any significant M66 conformational change. Based on an analysis that decomposes the absolute binding free energy into contributions from two receptor conformational states, it appears that it is the free energy cost of side chain reorganization rather than the reduced protein−ligand interaction that is largely responsible for the drug resistance against GS-6207.

Sun, Qinfang (ORCID:0000000173621543)↗
DOE OSTI2024

General circulation models simulate negative liquid water path–droplet number correlations, but anthropogenic aerosols still increase simulated liquid water path

General circulation models' (GCMs) estimates of the liquid water path adjustment to anthropogenic aerosol emissions differ in sign from other lines of evidence. This reduces confidence in estimates of the effective radiative forcing of the climate by aerosol–cloud interactions (ERFaci). The discrepancy is thought to stem in part from GCMs' inability to represent the turbulence–microphysics interactions in cloud-top entrainment, a mechanism that leads to a reduction in liquid water in response to an anthropogenic increase in aerosols. In the real atmosphere, enhanced cloud-top entrainment is thought to be the dominant adjustment mechanism for liquid water path, weakening the overall ERFaci. We show that the latest generation of GCMs includes models that produce a negative correlation between the present-day cloud droplet number and liquid water path, a key piece of observational evidence supporting liquid water path reduction by anthropogenic aerosols and one that earlier-generation GCMs could not reproduce. However, even in GCMs with this negative correlation, the increase in anthropogenic aerosols from preindustrial to present-day values still leads to an increase in the simulated liquid water path due to the parameterized precipitation suppression mechanism. This adds to the evidence that correlations in the present-day climate are not necessarily causal. We investigate sources of confounding to explain the noncausal correlation between liquid water path and droplet number. These results are a reminder that assessments of climate parameters based on multiple lines of evidence must carefully consider the complementary strengths of different lines when the lines disagree.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

High-frequency, high-power performance of AlGaN-channel high-electron-mobility transistors: an RF simulation study [High-frequency high power performance of AlGaN channel HEMTs: an RF simulation study]

The emerging Al-rich AlGaN-channel Al x Ga 1–x N/Al y Ga 1–y N high-electron-mobility transistors (HEMTs) with 0.7 ≤ y < x ≤ 1.0 have the potential to greatly exceed the power handling capabilities of today's GaN HEMTs, possibly by five times. This projection is based on the expected 4× enhancement of the critical electric field, the 2× enhancement of sheet carrier density, and the parity of the electron saturation velocity for Al-rich AlGaN-channel HEMTs relative to GaN-channel HEMTs. In this paper, the expected increased RF power density in Al-rich AlGaN-channel HEMTs is calculated by theoretical analysis and computer simulations, based on existing data on long-channel AlGaN devices. It is shown that a saturated power density of 18 W mm –1 , a power-added efficiency of 55% and an output third-order intercept point over 40 dB can be achieved for this technology. Furthermore, the method for large-signal RF performance estimation presented in this paper is generic and can be applied to other novel high-power device technologies at the early stages of development.

42 ENGINEERING↗
DOE OSTI2021

Suitable low Global Warming Potential (GWP) refrigerants for two-speed Heat Pumps for residential applications based on speed Heat Pumps for residential applications based on simulated performance simulated performance

Cold climate Heat Pumps (CCHPs) using single-speed compressors in tandem have been successfully built and tested to reach its performance metrics with outdoor temperatures as low as -25°C. CCHPs must be designed to alleviate the problems of too high discharge temperatures, low suction pressure and high-pressure ratio at low ambient conditions and insufficient heating capacity relative to the heating capacity at 8.3°C. In this paper we recognize those concerns. Low-GWP refrigerants are screened based on the shape of their T-S saturation boundary. Simulations of the performance of low-GWP refrigerants (R32, R454B, R466A, and R452B) relative to R410A is done using the DOE/ORNL Heat Pump Design Model at the component level. Systemic inefficiency, the root cause of lost work, is addressed by component-level exergy analysis to refine design options. CCHPs using low-GWP refrigerants addresses issues of reducing energy consumption, lowering carbon footprint, and enabling environmentally sustainability. The CCHP was designed to meet a HSPF > 11.0 and SEER > 15.0 having adequate heat pump heating capacity down to -25°C. It is estimated that at least 2.6 million U.S homes using electric furnaces and heat pumps in cold or very old climates with an annual energy consumption of 0.17EJ would benefit by CCHPs.

Ally, Moonis↗
DOE OSTI2016

Stabilized FE simulation of prototype thermal-hydraulics problems with integrated adjoint-based capabilities

A critical aspect of applying modern computational solution methods to complex multiphysics systems of relevance to nuclear reactor modeling, is the assessment of the predictive capability of specific proposed mathematical models. In this respect the understanding of numerical error, the sensitivity of the solution to parameters associated with input data, boundary condition uncertainty, and mathematical models is critical. Additionally, the ability to evaluate and or approximate the model efficiently, to allow development of a reasonable level of statistical diagnostics of the mathematical model and the physical system, is of central importance. In this study we report on initial efforts to apply integrated adjoint-based computational analysis and automatic differentiation tools to begin to address these issues. The study is carried out in the context of a Reynolds averaged Navier–Stokes approximation to turbulent fluid flow and heat transfer using a particular spatial discretization based on implicit fully-coupled stabilized FE methods. Initial results are presented that show the promise of these computational techniques in the context of nuclear reactor relevant prototype thermal-hydraulics problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

A mass and momentum conserving unsplit semi-Lagrangian framework for simulating multiphase flows

In this work, we present a computational methodology for convection and advection that handles discontinuities with second order accuracy and maintains conservation to machine precision. This method can transport a variety of discontinuous quantities and is used in the context of an incompressible gas–liquid flow to transport the phase interface, momentum, and scalars. The proposed method provides a modification to the three-dimensional, unsplit, second-order semi-Lagrangian flux method of Owkes & Desjardins (JCP, 2014). The modification adds a refined grid that provides consistent fluxes of mass and momentum defined on a staggered grid and discrete conservation of mass and momentum, even for flows with large density ratios. Additionally, the refined grid doubles the resolution of the interface without significantly increasing the computational cost over previous non-conservative schemes. This is possible due to a novel partitioning of the semi-Lagrangian fluxes into a small number of simplices. The proposed scheme is tested using canonical verification tests, rising bubbles, and an atomizing liquid jet.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Frustration- and doping-induced magnetism in a Fermi–Hubbard simulator

Geometrical frustration in strongly correlated systems can give rise to a plethora of novel ordered states and intriguing magnetic phases, such as quantum spin liquids. Promising candidate materials for such phases can be described by the Hubbard model on an anisotropic triangular lattice, a paradigmatic model capturing the interplay between strong correlations and magnetic frustration. However, the fate of frustrated magnetism in the presence of itinerant dopants remains unclear, as well as its connection to the doped phases of the square Hubbard model. Here we investigate the local spin order of a Hubbard model with controllable frustration and doping, using ultracold fermions in anisotropic optical lattices continuously tunable from a square to a triangular geometry. At half-filling and strong interactions U/t ≈ 9, we observe at the single-site level how frustration reduces the range of magnetic correlations and drives a transition from a collinear Néel antiferromagnet to a short-range correlated 120° spiral phase. Away from half-filling, the triangular limit shows enhanced antiferromagnetic correlations on the hole-doped side and a reversal to ferromagnetic correlations at particle dopings above 20%, hinting at the role of kinetic magnetism in frustrated systems. Furthermore, this work paves the way towards exploring possible chiral ordered or superconducting phases in triangular lattices and realizing t–t' square lattice Hubbard models that may be essential to describe superconductivity in cuprate materials.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2016

Initial performances of first undulator-based hard x-ray beamlines of NSLS-II compared to simulations

Commissioning of the first X-ray beamlines of NSLS-II included detailed measurements of spectral and spatial distributions of the radiation at different locations of the beamlines, from front-ends to sample positions. Additionally, comparison of some of these measurement results with high-accuracy calculations of synchrotron (undulator) emission and wavefront propagation through X-ray transport optics, performed using SRW code, is presented.

47 OTHER INSTRUMENTATION↗
DOE OSTI2016

Time-dependent analysis of visible helium line-ratios for electron temperature and density diagnostic using synthetic simulations on NSTX-U

Helium line-ratios for electron temperature (T{sub e}) and density (n{sub e}) plasma diagnostic in the Scrape-Off-Layer (SOL) and edge regions of tokamaks are widely used. Due to their intensities and proximity of wavelengths, the singlet, 667.8 and 728.1 nm, and triplet, 706.5 nm, visible lines have been typically preferred. Time-dependency of the triplet line (706.5 nm) has been previously analyzed in detail by including transient effects on line-ratios during gas-puff diagnostic applications. In this work, several line-ratio combinations within each of the two spin systems are analyzed with the purpose of eliminating transient effects to extend the application of this powerful diagnostic to high temporal resolution characterization of plasmas. The analysis is done using synthetic emission modeling and diagnostic for low electron density NSTX SOL plasma conditions by several visible lines. Quasi-static equilibrium and time-dependent models are employed to evaluate transient effects of the atomic population levels that may affect the derived electron temperatures and densities as the helium gas-puff penetrates the plasma. The analysis of a wider range of spectral lines will help to extend this powerful diagnostic to experiments where the wavelength range of the measured spectra may be constrained either by limitations of the spectrometer or by other conflicting lines from different ions.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗