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.
This is a presentation for the May 31 - June 3, 2023, SIAM conference with the goal to present advances in mathematical modeling and a case study as part of the Carbon Capture Simulation Initiative.
The University of Miami (Subcontractor)’s Research Scientist James Benedict (with oversight by PI Amy Clement) was tasked with producing global climate model simulations as part of the Cloud-Feedback Model Intercomparison Project (CFMIP) Tier 2 protocol, managing the output of these simulations, providing assistance (as requested) to staff at Lawrence Livermore National Lab (LLNL) regarding model output and setup, and meeting with LLNL scientists and team members to coordinate work on the project. The simulations represent an important contribution to the CFMIP data repository and will advance understanding of a wide range of critical cloud, circulation, and precipitation responses to climate change. The Subcontractor completed all proposed simulations, formatted the model output to be compliant with CFMIP protocols, and published the model data to the CFMIP repository. When requested, the Subcontractor provided assistance to LLNL for model output and configuration queries. Benedict and/or Clement also met in person or virtually 1-2 times per year with project scientists to coordinate work and review results.
Following the launch of multiple global navigation satel-lite system (GNSS) reflectometry (GNSS-R) missions, theSignals of Opportunity (SoOp) method has proven to be apowerful tool for geophysical parameter retrieval for land ap-plications such as soil moisture. Having demonstrated thefeasibility of the SoOp techniques at P- and S-band, the devel-opment of SoOp measurements beyond the GNSS frequencyregime is highly anticipated. The SoOp Coherent Bistatic(SCoBi) model and simulator, developed in 2017 and open-sourced in 2018, has been made available to provide multi-frequency, fully polarimetric SoOp simulations for ground-based applications through the joint use of analytical wavetheory and distorted Borne approximation to evaluate landcontributions from multilayer dielectric profiles composed ofsoil moisture, vegetation, and surface roughness effects. Thispaper describes the advancement of SCoBi from a ground-and airborne-based model to a spaceborne model. This ex-tension allows for fully polarimetric, complex delay-Dopplermap (DDM) simulations through evaluation of the coherentsuperposition of electric fields emerging from a grid of ori-ented facets. The model generates a grid of facets by de-termining the geometry of contributing elements from digitalelevation models, with each element providing its contribu-tion under a flat-earth assumption. This module will enablethe analysis of fully polarimetric scattering from frequenciesavailable across the ultra-high frequency (UHF) regime.
Advanced control algorithms for building systems are known to have significant potential in reducing energy consumption while optimizing thermal comfort. The success of such algorithms is critically contingent on several different types of sensor systems, which are in turn, used for continuous monitoring, identification and estimation of several important building states, such as temperatures, humidity, air quality, power consumption and occupancy status. Nonidealities in any of these sensors can lead to significant performance degradation of the control functionalities, and may lead to unwanted sub-optimal building operation. In this paper, we provide a simulation example with a high-fidelity building model, for a particular use-case of advanced optimization-based control in buildings, i.e., occupancy-based controls. We show how imperfections in occupancy sensing can offset performance. Subsequently, we discuss a novel learning-based architecture to efficiently evaluate the impact of sensor nonidealities for building systems, in context of advanced control algorithms.
This report summarizes the design and application of an outdoor power line sensor testbed (OPLST) with a real-time simulator and power meter to compare a potential transformer (PT) and current transformer (CT) versus an advanced line-post sensor. The OPLST was created to validate advanced medium-voltage (20/34.5 kV) outdoor power line sensors used in electrical distribution systems. Electrical utilities have installed metering and relay-protection transformers like PTs and CTs for several decades. The PTs/CTs are iron core measurement transformers based on the electromagnetic induction principle and provide reliable data in normal grid operation. However, new outdoor power line sensors (OPLS) using other technologies like voltage divider, Rogowski coil and optical principles have become available and may have favorable performance and costs compared to PTs/CTs. Therefore, the importance of testing these technologies with the PT/CT, to compare the measured phase voltage/current at different power grid scenarios is crucial to understand the performance of these new OPLS. For this study a G&W Electric Model CVS-36-O power line sensor was chosen as the OPLS to test with voltage/current signals. An OPAL-RT Technologies Model OP4510 real-time simulator and Schweitzer Engineering Laboratories Model SEL-735 power meter were installed with the 20/34.5 kV OPLST to compare the measured transient events collected from an advanced OPLS and the PT/CT. This system is installed at the Distributed Energy Communications and Control (DECC) lab, Oak Ridge National Laboratory (ORNL). The simulator generated different power grid scenarios (electrical faults, capacitor bank operation, service restoration, etc.), and its analog-output signals were connected to the voltage/current amplifiers that feed the 20/34.5 kV aerial cable loop through the PT/CT devices. Additional PT/CT devices were also wired with the medium voltage aerial cable loop to measure the phase current/voltage signals and server as references. After each test, common format for transient data exchange (COMTRADE) files were collected from the SEL-735 power meter and used to compare the performance of the OPLS with the PT/CT. The behavior of analog signals, harmonic components, total harmonic distortion and crest factors were assessed and found favorable for the G&W sensor as compared to the reference PT/CT.
We report on the development progress of a hybrid fluid-kinetic code for simulating fluids and plasmas in a wide range of environments, such as laser–matter interactions, inertial confinement fusion, magnetic confinement fusion, and pulsed power. The suite of numerical tools under development utilizes heterogeneous computer architectures and leverages the benefits of particle–based simulation techniques. By working to combine the kinetic particle-in-cell (PIC) model with a particle-based fluid simulation technique, such as smoothed particle hydrodynamics, we are developing a flexible framework capable of accurately modeling complex flows within and between kinetic and fluid regimes. The TriForce code is under development as a C++ framework for parallel, 3D, particle-based, hybrid fluid-kinetic plasma simulations. The fluid half of TriForce will be based upon the meshless smoothed-particle-hydrodynamics (SPH) approach, well-suited for shear, mixing, and turbulence, whereas the kinetic half resembles a traditional particle-in-cell (PIC) code; other particle-based approaches to fluid modeling that do use a mesh are also possible to use and are under investigation. Maxwell’s electromagnetic field equations are solved either via explicit or implicit algorithms or approximated via resistive magnetohydrodynamics (MHD) using an Ohm’s law and resulting induction equation (extended MHD is under development). A primary goal of enabling direct comparisons, from the same code, between results from the variants of MHD and implicit electromagnetic solutions is to improve our fundamental understanding of systems with magnetic fields. The code is under development to recover results from both radiation-MHD and fully kinetic codes in those limits, and is continuing to be developed from other follow-on grants to operate in between where both descriptions may co-exist and interact. For certain applications, it is desired for a simulation to contain fluid ions and electrons as well as kinetic ions and electrons. Typically, it is too computationally intensive to model a full-scale ICF or HEDP experiment fully kinetically since many cycles are expended with very small time steps on modeling the fluid part of a material that is well treated by the fluid approximation. In this case, many traditional PIC particles can be replaced with a single fluid particle representing the thermal part of the distribution function, and there are fewer needed kinetic particles, which describe the non-thermal part and can be sub-cycled relative to the fluid particle advance. Furthermore, a pure fluid code may, depending on the problem, simply lack many physically important details that are beyond the scope of its reduced approximations and assumptions. In this report, we summarize the objectives achieved in the development of the collisional and kinetic half of the code, and the physics problems to which the code has been applied in the areas of advanced and innovative fusion concepts, pulsed power, and magneto-inertial fusion.
Score based generative models are a new class of generative models that have been shown to accurately generate high dimensional calorimeter datasets. Recent advances in generative models have used images with 3D voxels to represent and model complex calorimeter showers. Point clouds, however, are likely a more natural representation of calorimeter showers, particularly in calorimeters with high granularity. Furthermore, point clouds preserve all of the information of the original simulation, more naturally deal with sparse datasets, and can be implemented with more compact models and data files. In this work, two state-of-the-art score based models are trained on the same set of calorimeter simulation and directly compared.
There is a critical need for the development and verification of practically useful multiscale modeling strategies for simulating the mechanical response of multiphase metallic materials with heterogeneous microstructures. In this contribution, we present data-driven reduced order models for effective yield strength and strain partitioning in such microstructures. These models are built employing the recently developed framework of Materials Knowledge Systems that employ 2-point spatial correlations (or 2-point statistics) for the quantification of the heterostructures and principal component analyses for their low-dimensional representation. The models are calibrated to a large collection of finite element (FE) results obtained for a diverse range of microstructures with various sizes, shapes, and volume fractions of the phases. The performance of the models is evaluated by comparing the predictions of yield strength and strain partitioning in two-phase materials with the corresponding predictions from a classical self-consistent model as well as results of full-field FE simulations. The reduced-order models developed in this work show an excellent combination of accuracy and computational efficiency, and therefore present an important advance towards computationally efficient microstructure-sensitive multiscale modeling frameworks. - Graphical abstract:.
The first measurements of fast-ion D-alpha (FIDA) radiation have been acquired on C-2U, Tri Alpha Energy’s advanced, beam-driven field-reversed configuration (FRC). These measurements are also forward modeled by FIDASIM. This is the first measurement and simulation of FIDA carried out on an FRC topology. FIDA measurements are made of Doppler-shifted Balmer-alpha light from neutralized fast ions using a bandpass filter and photomultiplier tube. One adjustable line-of-sight measured signals at eight locations and eight times during the FRC lifetime over 26 discharges. Filtered signals include only the highest energy ions (>6 keV) and share some salient features with the FIDASIM result. Highly Doppler-shifted beam radiation is also measured with a high-speed camera and is spatially well-correlated with FIDASIM.
Scientific Machine Learning (SciML) and Artificial Intelligence (AI) will have broad use and transformative effects across the Department of Energy. Accordingly, the January 2018 Basic Research Needs workshop identified six Priority Research Directions (PRDs). The first three PRDs describe foundational research themes that correspond to the need for domain-awareness (PRD #1), interpretability (PRD #2), and robustness (PRD #3). The other three PRDs describe capability research themes and correspond to the three major use cases for massive scientific data analysis (PRD #4), machine learning-enhanced modeling and simulation (PRD #5), and intelligent automation and decision-support for complex systems (PRD #6). The Priority Research Directions provide a sound basis for a coherent, long-term research and development strategy in SciML and AI. Over the last decade, DOE investments in applied mathematics have laid the groundwork for the type of basic research that will underpin key advances in the six PRDs. Such advances will build on the work from leading researchers in optimization, linear algebra, high-performance solvers and algorithms, multiscale modeling and simulation, complex systems research, uncertainty quantification, and the new basic research areas that will emerge from the pursuit of transformative technologies.
Scientific Machine Learning (SciML) and Artificial Intelligence (AI) will have broad use and transformative effects across the Department of Energy. Accordingly, the January 2018 Basic Research Needs workshop identified six Priority Research Directions (PRDs). The first three PRDs describe foundational research themes that correspond to the need for domain-awareness (PRD #1), interpretability (PRD #2), and robustness (PRD #3). The other three PRDs describe capability research themes and correspond to the three major use cases for massive scientific data analysis (PRD #4), machine learning-enhanced modeling and simulation (PRD #5), and intelligent automation and decision-support for complex systems (PRD #6). The Priority Research Directions provide a sound basis for a coherent, long-term research and development strategy in SciML and AI. Over the last decade, DOE investments in applied mathematics have laid the groundwork for the type of basic research that will underpin key advances in the six PRDs. Such advances will build on the work from leading researchers in optimization, linear algebra, high-performance solvers and algorithms, multiscale modeling and simulation, complex systems research, uncertainty quantification, and the new basic research areas that will emerge from the pursuit of transformative technologies.
Recent renewed interest in Nuclear Electric Propulsion (NEP) at NASA has driven the need for design space exploration and analysis of alternative architectures to evaluate performance and feasibility of NEP concepts. It has been over a decade since NASA last invested in Nuclear Electric Propulsion development efforts, and the landscape has shifted in both mission objectives and technology readiness levels. This also means that many of the engineering analyses and design frameworks tailored to NEP date back to the Jupiter Icy Moons Orbiter (JIMO) days of the early 2000s. Thus, a need was identified for updating and modernizing the modeling and simulation toolsets for NEP. This paper will detail the development efforts in the Marshall Space Flight Center’s Advanced Concepts Office, with a particular focus on the development of the NEP stage of the NEP-Chem vehicle. The paper will cover the various models used, and their current development status, as well as provide an overview of the integrated vehicle model and plans for further development.
Recent renewed interest in Nuclear Electric Propulsion (NEP) at NASA has driven the need for design space exploration and analysis of alternative architectures to evaluate performance and feasibility of NEP concepts. It has been over a decade since NASA last invested in Nuclear Electric Propulsion development efforts, and the landscape has shifted in both mission objectives and technology readiness levels. This also means that many of the engineering analyses and design frameworks tailored to NEP date back to the Jupiter Icy Moons Orbiter (JIMO) days of the early 2000s. Thus, a need was identified for updating and modernizing the modeling and simulation toolsets for NEP. This paper will detail the development efforts in the Marshall Space Flight Center’s Advanced Concepts Office, with a particular focus on the development of the NEP stage of the NEP-Chem vehicle. The paper will cover the various models used, and their current development status, as well as provide an overview of the integrated vehicle model and plans for further development.
Dust aerosols have a wide size distribution from less than 0.1 to over 100 μm and dominate Earth's atmospheric aerosol mass. However, most Earth system models (ESMs) inadequately represent dust aerosols larger than 10 μm in diameter, limiting the accuracy of the simulated dust cycle and climate impacts. Here, we introduce a new modeling framework that captures the full observed size distribution of dust aerosols, incorporating recent advances into a mineral-resolved version of the Community ESM, while addressing known issues in previous versions. Comprehensive evaluation against diverse observations of bulk dust and component minerals demonstrates that the model reproduces the observed dust cycle across particle sizes. Incorporating the previously unrepresented large-dust fractions substantially alters dust budget estimates, highlighting potential changes in simulated climate impacts and underscoring the importance of comprehensive size-resolved dust modeling. Despite these advancements, uncertainties persist. Our results indicate that a size-dependent reduction in settling velocity is required to reproduce the observed dust size distribution downwind of source regions. Specifically, in the new model, the gravitational settling velocity of dust particles larger than 10 μm in diameter must be reduced by as much as 85% to achieve agreement with observations. This empirical reduction serves as a constraint on physics-based models of dust settling. Future developments should address misrepresented physical processes that hinder accurate modeling of the large dust aerosol transport. Expanding observational data sets covering the full-size distribution is also essential to better constrain the dust cycle and improve the representation of dust optical properties and climate effects.
Dust aerosols have a wide size distribution from less than 0.1 to over 100 μm and dominate Earth's atmospheric aerosol mass. However, most Earth system models (ESMs) inadequately represent dust aerosols larger than 10 μm in diameter, limiting the accuracy of the simulated dust cycle and climate impacts. Here, we introduce a new modeling framework that captures the full observed size distribution of dust aerosols, incorporating recent advances into a mineral-resolved version of the Community ESM, while addressing known issues in previous versions. Comprehensive evaluation against diverse observations of bulk dust and component minerals demonstrates that the model reproduces the observed dust cycle across particle sizes. Incorporating the previously unrepresented large-dust fractions substantially alters dust budget estimates, highlighting potential changes in simulated climate impacts and underscoring the importance of comprehensive size-resolved dust modeling. Despite these advancements, uncertainties persist. Our results indicate that a size-dependent reduction in settling velocity is required to reproduce the observed dust size distribution downwind of source regions. Specifically, in the new model, the gravitational settling velocity of dust particles larger than 10 μm in diameter must be reduced by as much as 85% to achieve agreement with observations. This empirical reduction serves as a constraint on physics-based models of dust settling. Future developments should address misrepresented physical processes that hinder accurate modeling of the large dust aerosol transport. Expanding observational data sets covering the full-size distribution is also essential to better constrain the dust cycle and improve the representation of dust optical properties and climate effects.
Electronic transistors are fundamental building blocks of large scale integrated circuits in modern advanced electronic equipments, and their sizes have been down-scaled to nanometers. Modeling and simulations in the framework of quantum dynamics have emerged as important tools to study functional characteristics of these nano-devices. This work explores the effects of geometric shapes of semiconductor–insulator interfaces, phonon–electron interactions, and quantum tunneling of three-dimensional (3D) nano-transistors. First, we propose a two-scale energy functional to describe the electron dynamics in a dielectric continuum of device material. Coupled governing equations, i.e., Poisson–Kohn–Sham (PKS) equations, are derived by the variational principle. Additionally, it is found that at a given channel cross section area and gate voltage, the geometry that has the smallest perimeter of the channel cross section offers the largest channel current, which indicates that ultra-thin nanotransistors may not be very efficient in practical applications. Moreover, we introduce a new method to evaluate quantum tunneling effects in nanotransistors without invoking the comparison of classical and quantum predictions. It is found that at a given channel cross section area and gate voltage, the geometry that has the smallest perimeter of the channel cross section has the smallest quantum tunneling ratio, which indicates that geometric defects can lead to higher geometric confinement and larger quantum tunneling effect. Furthermore, although an increase in the phonon–electron interaction strength reduces channel current, it does not have much impact to the quantum tunneling ratio. Finally, advanced numerical techniques, including second order elliptic interface methods, have been applied to ensure computational accuracy and reliability of the present PKS simulation.
The ParFlow Sand Tank model is an open source application designed to allow users to interactively simulate and visualize groundwater movement through the subsurface. The app is designed for both research and education; teaching hydrogeology concepts and making it easy explore and run sophisticated groundwater simulations. Our goal is to support increased accessibility and usability of research grade hydrology tools for research and teaching. The Sand Tank application simulates groundwater and surface water fluxes as well as contaminant transport in real time using the integrated physical hydrology model ParFlow (Kollet & Maxwell, 2006; Maxwell & Miller, 2005; Osei-Kuffuor et al., 2014) and the particle tracking code EcoSlim (Maxwell et al., 2019). ParFlow is a numerical hydrology model that simulates spatially distributed groundwater and surface water flow. It is a well established research tool with more than 90 publications documenting its development use to advance our understanding of groundwater dynamics and groundwater surface water interactions from the hillslope to the continental scale e.g. (Condon et al., 2020; Condon & Maxwell, 2019; Maxwell & Condon, 2016). It is designed for efficient parallel computation and has been run on many platforms spanning from laptops to supercomputers. However, one of the challenges of ParFlow is that it requires significant training and hydrologic expertise to develop simulations. The Sand Tank application makes this model accessible to anyone for education and exploration. Our application uses ParFlow for its simulation backend and ParaView for the data loading and processing. The communication infrastructure relies on the ParaViewWeb framework. We use model templates deployed in Docker images to setup the Sand Tank framework. Users can build the application locally or interact with it through our web deployment. When interacting with a template users can interactively change model parameters like subsurface processes or pump/inject water into the subsurface and watch the system respond to their changes in real time as the simulation runs. Additionally, our template setup will allow more advanced users to build custom templates of increasing complexity for both research and educational purposes.
ABSTRACT Urban areas experience the impact of natural disasters, such as heatwaves and flash floods, disparately in different neighbourhoods across a city. The demand for precise urban hydrometeorological and hydroclimatological modelling to examine this disparity, and the interacting challenges posed by climate change and urbanisation, has thus surged. The Weather Research and Forecasting (WRF) model has served such operational and research purposes for decades. Recent advancements in WRF, including enhanced numerical schemes and sophisticated urban atmospheric‐hydrological parameterizations, have empowered the simulation of urban geophysical processes at high resolution (~1 km), but even this resolution misses significant urban microclimate variability. This study applies the large‐eddy simulations (LES) mode within WRF, coupled with single‐layer urban canopy models (SLUCM), to enable even finer‐scale modelling (150 m) of the Urban Heat Island (UHI) effect in the Baltimore metropolitan area. We run nine scenarios to evaluate various methods of initializing soil moisture and various spinup lead times, and to assess the impact of WRF's Mosaic approach in depicting subgrid‐scale processes. We evaluate the scenarios by comparing the WRF simulated land surface temperature (LST) against Landsat LST and the WRF simulated hourly 2‐m air temperatures (AT) with observations from eight weather stations across the domain. Results underscore the paramount influence of the lead spinup time on the spatiotemporal distribution of simulated soil moisture, consequently shaping WRF's efficacy in predicting the UHI. Furthermore, interpolating soil moisture‐related parameters from the parent for child domain initialization yields a notable reduction in mean and root‐mean‐squared errors. This improvement was particularly evident in simulations with the longest spinup time, affirming the importance of carefully designing the initialization of soil moisture for improved urban temperature predictions.