SEARCH · Engineering Papers
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.
Kynema-driver (FKA: exawind-driver) is a C++ code that is part of the kynema software stack. It is a driver for coupled Kynema-SGF/UGF simulations. It was designed to couple and drive hybrid-solver computational fluid dynamics (CFD) simulations where NLR'S kynema-sgf (SWR-20-85) software and NLR'S kynema-ugf (SWR-20-27) CFD codes are run simultaneously and are two-way coupled via overset meshes and the TIOGA overset-mesh library.
The Superconducting Quantum Materials and Systems (SQMS) Center, a Fermilab lead NQI center, will develop and deploy a beyond state of the art quantum computers. These systems will exhibit long coherent times, and a high degree of connectivity. We are interested in exploring the utility of these devices to address computational challenges in HEP, such as the simulation of strongly coupled quantum field theories. A progression of toy models and simulations can begin to chart a course towards simulating QCD-like systems.
A primary concern of CO2-induced warming is the associated rise of tropical (10S-10N) seasurface temperatures (SSTs). GISS Model-E was used to produce two sets of simulations-one with the present-day and one with doubled CO2 in the atmosphere. The intrinsic usefulness of model guidance in the tropics was confirmed when the model simulated realistic convective coupling between SSTs and atmospheric soundings and that the simulated-data correlations between SSTs and 300 hPa moiststatic energies were found to be similar to the observed. Model predicted SST limits: (i) one for the onset of deep convection and (ii) one for maximum SST, increased in the doubled C02 case. Changes in cloud heights, cloud frequencies, and cloud mass-fractions showed that convective-cloud changes increased the SSTs, while warmer mixed-layer of the doubled CO2 contained approximately 10% more water vapor; clearly that would be conducive to more intense storms and hurricanes.
Water availability is fundamental to societies and ecosystems, but our understanding of variations in hydroclimate (including extreme events, flooding, and decadal periods of drought) is limited because of a paucity of modern instrumental observations that are distributed unevenly across the globe and only span parts of the 20th and 21st centuries. Such data coverage is insufficient for characterizing hydroclimate and its associated dynamics because of its multidecadal to centennial variability and highly regionalized spatial signature. High-resolution (seasonal to decadal) hydroclimatic proxies that span all or parts of the Common Era (CE) and paleoclimate simulations from climate models are therefore important tools for augmenting our understanding of hydroclimate variability. In particular, the comparison of the two sources of information is critical for addressing the uncertainties and limitations of both while enriching each of their interpretations. We review the principal proxy data available for hydroclimatic reconstructions over the CE and highlight the contemporary understanding of how these proxies are interpreted as hydroclimate indicators. We also review the available last-millennium simulations from fully coupled climate models and discuss several outstanding challenges associated with simulating hydroclimate variability and change over the CE. A specific review of simulated hydroclimatic changes forced by volcanic events is provided, as is a discussion of expected improvements in estimated radiative forcings, models, and their implementation in the future. Our review of hydroclimatic proxies and last-millennium model simulations is used as the basis for articulating a variety of considerations and best practices for how to perform proxy-model comparisons of CE hydroclimate. This discussion provides a framework for how best to evaluate hydroclimate variability and its associated dynamics using these comparisons and how they can better inform interpretations of both proxy data and model simulations.We subsequently explore means of using proxy-model comparisons to better constrain and characterize future hydroclimate risks. This is explored specifically in the context of several examples that demonstrate how proxy-model comparisons can be used to quantitatively constrain future hydroclimatic risks as estimated from climate model projections.
Estimate the mass of the Thermal Protection System (TPS) for a single design construct of an atmospheric entry probe with a rigid aeroshell, which could be used at five destinations, i.e. Venus, Saturn, Uranus, Neptune, and perhaps, Jupiter. The entry mass of the probe is 400 kg with a ballistic coefficient of 216 kg/m2. Process: The 3DoF trajectory simulation program Traj, coupled with the TPS response program FIAT was used for simulation and design. The assumed atmospheric models were VIRA (Venus-GRAM) for Venus, the Julianne Moses' model for Saturn, a NASA Ames engineering model for Uranus, Neptune-GRAM for Neptune, and Galileo Probe (Al Seiff's) result for Jupiter.
In this Letter, we explore a non-degenerate phase independent parametric quality factor (Q)-enhancement technique for aluminum nitride (AlN) Lamb wave resonators. Unlike other active Q-enhancement techniques which require precise phase control of the electronic feedback loop, this technique is implemented by parametrically pumping AlN material stiffness to realize a negative resistance seen at the signal path. The negative resistance is dependent on the nonlinear material modulation and multi-resonance coupling in the device. A nonlinear circuit model is developed to simulate the parametric coupling of each resonance and extract the nonlinearity of AlN from experimental data. With proper pump frequency and pump power, the device quality factor is boosted in both simulation and experiment. The demonstrated Q-enhancement method is simple to implement and can be applied to other types of resonators that have nonlinear behavior and support multi-resonance operation.
Abstract Atom‐ and energy‐efficient chemical separations are urgently needed to meet the surging demand for critical materials that has strained supply chains and threatened environmental damage. In this study, we used reaction‐diffusion coupling to separate iron, neodymium, and dysprosium ions from model feedstocks of permanent magnets, which are typically found in electronic wastes. Feedstock solutions were placed in contact with a hydrogel loaded with potassium hydroxide and/or dibutyl phosphate, resulting in complex precipitation patterns as the various metal ions diffused into the reaction medium. Specifically, we observed the precipitation of up to 40 mM of iron from the feedstock, followed by the enrichment of 73 % dysprosium, and the extraction of >95 % neodymium product at a further distance from the solution‐gel interface. We designed a series of experiments and simulations to determine the relevant ion diffusivities, D Nd =5.4×10 −10 and D Dy =5.1×10 −10 m 2 /s, and precipitation rates, k Nd =1.0×10 −5 and k Dy =5.0×10 −3 m 9 mol −3 s −1 , which enabled a numerical model to be established for predicting the distribution of products in the reaction medium. Our proof‐of‐concept study validates reaction‐diffusion coupling as an effective and versatile approach for critical materials separations, without relying on ligands, membranes, resins, or other specialty chemicals.
Modern cosmic microwave background (CMB) analysis pipelines regularly employ complex time-domain filters, beam models, masking, and other techniques during the production of sky maps and their corresponding angular power spectra. However, these processes can generate couplings between multipoles from the same spectrum and from different spectra, in addition to the typical power attenuation. Within the context of pseudo-C ℓ based, MASTER-style analyses, the net effect of the time-domain filtering is commonly approximated by a multiplicative transfer function, F ℓ , that can fail to capture mode mixing and is dependent on the spectrum of the signal. To address these shortcomings, we have developed a simulation-based spectral correction approach that constructs a two-dimensional transfer matrix, ${J}_{{\ell }{\ell }^{\prime} }$, which contains information about mode mixing in addition to mode attenuation. We demonstrate the application of this approach on data from the first flight of the Spider balloon-borne CMB experiment.
This paper presents a numerical assessment of acoustic installation effects in the tandem cylinder (TC) experiments conducted in the NASA Langley Quiet Flow Facility (QFF), an open-jet, anechoic wind tunnel. Calculations that couple the Computational Fluid Dynamics (CFD) and Computational Aeroacoustics (CAA) of the TC configuration within the QFF are conducted using the CFD simulation results previously obtained at NASA LaRC. The coupled simulations enable the assessment of installation effects associated with several specific features in the QFF facility that may have impacted the measured acoustic signature during the experiment. The CFD-CAA coupling is based on CFD data along a suitably chosen surface, and employs a technique that was recently improved to account for installed configurations involving acoustic backscatter into the CFD domain. First, a CFD-CAA calculation is conducted for an isolated TC configuration to assess the coupling approach, as well as to generate a reference solution for subsequent assessments of QFF installation effects. Direct comparisons between the CFD-CAA calculations associated with the various installed configurations allow the assessment of the effects of each component (nozzle, collector, etc.) or feature (confined vs. free jet flow, etc.) characterizing the NASA LaRC QFF facility.
Tristructural isotropic (TRISO)-coated particle fuel is a robust nuclear fuel and determining its reliability is critical for the success of advanced nuclear technologies. However, TRISO failure probabilities are small and the associated computational models are expensive. Here we used coupled active learning, multifidelity modeling, and subset simulation to estimate the failure probabilities of TRISO fuels using several 1D and 2D models. With multifidelity modeling, we replaced expensive high-fidelity (HF) model evaluations with information fusion from two low-fidelity (LF) models. For the 1D TRISO models, we considered three multifidelity modeling strategies: only Kriging, Kriging LF prediction plus Kriging correction, and deep neural network (DNN) LF prediction plus Kriging correction. While the results across these multifidelity modeling strategies compared satisfactorily, strategies employing information fusion from two LF models called the HF model least often. Next, for the 2D TRISO model, we considered two multifidelity modeling strategies: DNN LF prediction plus Kriging correction (data-driven) and 1D TRISO LF prediction plus Kriging correction (physics-based). The physics-based strategy, as expected, consistently required the fewest calls to the HF model. However, the data-driven strategy had a lower overall simulation time since the DNN predictions are instantaneous, and the 1D TRISO model requires a non-negligible simulation time.
The Large Surface Explosion Coupling Experiment (LSECE) project will include two surface energetic events (chemical explosions) conducted as close as possible to ground zero (GZ) of the Source Physics Experiment (SPE) Phase II Dry Alluvium (DAG) site in Yucca Flat at the Nevada National Security Site (NNSS). Seismic and acoustic data from these events will complement DAG buried source data and improve interpretation of speed-of-sound signals from explosions for national security programs. A satellite image of the region near GZ along with some of the recording stations that we will use for the far- eld simulation is shown in Figure 1. The working design of the HE package is shown in the Figure 2 and is roughly a cube of C-4 bricks.
Energy storage allows flexible use and management of excess electricity and intermittently available renewable energy. Cryogenic energy storage (CES) is a promising storage alternative with a high technology readiness level and maturity, but the round-trip efficiency is often moderate and the Levelized Cost of Storage (LCOS) remains high. The complex flowsheets with intricate thermodynamics at cryogenic temperatures as well as the presence of multiple loops and refrigeration cycles pose considerable challenges for rigorous model-based design and optimization of CES systems. We present an optimization strategy that couples rigorous process simulation and Bayesian optimization with flowsheet decomposition and identification of hidden coupling constraints to optimally design standalone CES systems. Further refinement is done via a local search using the limited-memory Broyden–Fletcher–Goldfarb–Shanno algorithm. Here our results indicate that it is possible to achieve more than 52% round-trip efficiency and an LCOS of $153/MWh for a standalone 100 MW/400 MWh CES system limited to short-term storage with daily charging–discharging. However, a detailed techno-economic assessment reveals that the LCOS considering total capital investment may exceed $267/MWh when all direct and indirect costs of installation and operation are considered.
This report describes the development, modeling, and results of a reduced-order thermal power dispatch pressurized water reactor (RO-TPD-PWR) power plant simulator that incorporates coupled electrical and thermal power dispatch to an industrial process. The simulator is built on models that have been developed in previous work, including a Rancor Microworld model previously developed by INL and the University of Idaho [Ref] and reactor core and secondary system models that have been published previously [Ref]. The simulator also includes a scalable model of a high temperature electrolysis (HTE) system that has been developed at INL. A particular benefit of the Rancor Microworld model is that it is compatible with a human-system interface (HMI) that has already been developed by the Human System Simulation Laboratory (HSSL) at INL that can emulate the control room of a PWR for realistic human-in-the-loop studies with mock nuclear power plant operating procedures. The results of these human-in-the-loop tests will be compared with similar tests using a commercial full-scope generic PWR (GPWR) simulator that has been modified to include thermal power dispatch, referred to as the TPD-GPWR [ , ]. The simplicity of the RO-TPD-PWR helps to better understand and verify the results of the TPD-GPWR, which employs >1,000 variables and has much greater complexity. A second key purpose of the RO-TPD-PWR simulator is that its relative simplicity allows it scale between lab-size equipment (~100 kW) and full-scale nuclear power plants (~1 GW) by adjusting only a few parameters in the models, such as fluid masses and heat exchanger areas. The RO-TPD-PWR simulator will be used to guide the controls of lab tests, so that equipment in those tests, such as the Thermal Energy Distribution System (TEDS) and pilot-scale high temperature electrolysis (HTE) units at scales of 50-250 kW, may operate in a manner that is relevant to full-scale PWR operations.
A number of coupled particle-element and hybrid particle-element methods have been developed for the simulation of hypervelocity impact problems, to avoid certain disadvantages associated with the use of pure continuum based or pure particle based methods. To date these methods have employed spherical particles. In recent work a hybrid formulation has been extended to the ellipsoidal particle case. A model formulation approach based on Lagrange's equations, with particles entropies serving as generalized coordinates, avoids the angular momentum conservation problems which have been reported with ellipsoidal smooth particle hydrodynamics models.
Geothermal energy is sustainable and gaining momentum as a solution to energy crises and environmental issues. However, challenges like production temperature and thermal breakthrough can impact geothermal project efficiency. One innovative solution to alleviate the thermal breakthrough is to inject polymer-based materials that are encapsulated in microcapsules into fractures to modify fracture permeability and prevent preferential flow. In our study, we utilized a coupled computational fluid dynamics and discrete element method to simulate the transport of microcapsules under various scenarios controlled by microcapsule size, microcapsule concentration, and fracture roughness. For a smooth fracture, the results indicate that small microcapsules can travel through a smooth fracture regardless of their concentrations. Large microcapsules can transport through a smooth fracture when present in lower concentrations. However, medium and mixed-size microcapsules tend to cause the sealing of a smooth fracture, irrespective of their concentrations. For a rough fracture, the transport of microcapsules is complicated by their interactions with the rough fracture walls. The presence of two sealing positions in a rough fracture adds further complexity to this transport phenomenon. The size and concentration of microcapsules control one sealing location, while the rough fracture walls determine the other sealing location. The rough walls substantially affect microcapsule transport, rendering the role of microcapsule size and concentration less significant. The simulation results suggest that complex fracture surfaces significantly elevate the occurrence of sealing behavior. To mitigate sealing behavior within more complex fractures, it would be beneficial to use smaller and lower concentrations of microcapsules.
Here, we present a novel mass-momentum consistent coupling between a geometric volume-of-fluid scheme and an incompressible flow solver with differing directional-splitting approaches. The advection of the volume fraction is performed using a direction-split algorithm, whereas the momentum advection algorithm uses a traditional unsplit, fractional-step approach. Both algorithms employ finite-volume discretizations based on Cartesian meshes. In solving the mass-momentum consistency problem, momentum fluxes at the cell faces are weighted by the density fluxes based on the already advected volume fraction. The success of our approach lies on introducing a Favre-averaged velocity interpolation at the liquid/gas interface along with a minmod slope limiter. Mesh-convergence studies show that when the minmod slope limiter is used, the two-phase solver retains an accuracy between first and second order, but when a purely upwind scheme is considered, its accuracy drops to first order. Finally, after considering several validation problems, the solver is shown to agree well with reference numerical and experimental data while retaining its robustness and efficiency.
This paper presents the development of novel models that integrate induction machines with HVAC equipment, such as pumps, heat pumps, and chillers, to analyze the impact of electrical parameters on the operational performance of thermo-fluid systems. The proposed model employs a coupling technique that captures the dynamic interactions between induction machines and HVAC systems. By integrating electrical, thermal, and mechanical dynamics, the models provide a comprehensive framework for simulating real-world scenarios, including interactions with the electrical grid. This achievement was made possible through the development of a Computationally Efficient and Accurate Induction Machine (CEAIM) model. Implemented using the equation-based Modelica language, the CEAIM model has been validated against experimental results, manufacturer data sheets, and various operating conditions. Its performance has been compared with existing induction machine models in the Modelica Standard Library (MSL), demonstrating superior accuracy and computational efficiency. The CEAIM model predicts torque, speed, and power consumption with a coefficient of determination (R 2 ) ranging from 0.98 to 1 and a coefficient of variation of root mean square error (CVRMSE) between 0.27% and 6.67%. Additionally, CEAIM scales more efficiently than conventional MSL models, with a slower computational growth rate in large-scale simulations. After thorough validation of the CEAIM model, it was coupled with HVAC equipment as this approach provides a detailed multi-dimensional view of capturing electrical transients and mechanical performance. To support this, a case study was conducted to showcase its capabilities.
Significant progress has been made in the design, construction and operation of induction coupled devices which simulate the open cycle Gas Core Nuclear Rocket. These devices incorporate solid feed of the plasma forming material (uranium and uranium simulating materials), permeable walls, seeding of the propellant, and transpiration cooled, choked flow nozzles. Operating parameters and performance data of devices employing these design features are discussed. A uranium plasma experiment is included. In addition, operating data of several devices which compare theoretical and actual performance at a variety of powers, pressures, frequencies, and sizes are discussed.