Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation Cases”

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 775 records · Page 43

Predicting Frigid Mixed-Phase Clouds for Pristine Coastal Antarctica

Supercooled water is common in the clouds near coastal Antarctica and occasionally occurs at temperatures at or below -30°C. Yet the ice physics in most regional and global numerical models will glaciate out these clouds. This presents a challenge for the simulation of highly supercooled clouds that were observed at McMurdo, Antarctica during the Atmospheric Radiation Measurement (ARM) West Antarctic Radiation Experiment (AWARE) project during 2015–2017. The polar optimized version of the Weather Research and Forecasting model (Polar WRF) with the recently developed two-moment P3 microphysics scheme was used to simulate observed supercooled liquid water cases during March and November 2016. Nudging of the simulations to observed rawinsonde profiles and Antarctic automatic weather station observations provided increased realism and much greater cloud water amounts. Furthermore, sensitivity tests that adjust the ice physics for extremely low ice nucleating particle (INP) concentrations decrease cloud ice and increases the cloud liquid water closer to observed amounts. In these tests, a liquid layer near cloud top is simulated, in agreement with observations. Accurate representation of INP concentrations appears to be critical for the simulation of coastal Antarctic clouds.

54 ENVIRONMENTAL SCIENCES↗

Assessing Disparities of Dengue Virus Transmission Risk across the US-Mexican Border Using a Climate Driven Vector-Epidemiological Model

Dengue fever is a mosquito‐borne viral disease reemerging throughout much of the tropical Americas. Dengue virus transmission is explicitly influenced by climate and the environment through its primary vector, Aedes aegypti. Temperature regulates Ae. aegypti development, survival, and replication rates as well as the incubation period of the virus within the mosquito. Precipitation provides water for many of the preferred breeding habitats of the mosquito, including buckets, old tires, and other places water can collect. Although transmission regularly occurs along the border region in Mexico, dengue virus transmission in bordering Arizona has not occurred. Using NASA's TRMM (Tropical Rainfall Measuring Mission) satellite for precipitation input and Daymet for temperature and supplemental precipitation input, we modeled dengue transmission along a US‐Mexico transect using a dynamic dengue transmission model that includes interacting vector ecology and epidemiological components. Model runs were performed for 5 cities in Sonora, Mexico and southern Arizona. Employing a Monte Carlo approach, we performed ensembles of several thousands of model simulations in order to resolve the model uncertainty arising from using different combinations of parameter values that are not well known. For cities with reported dengue case data, the top model simulations that best reproduced dengue case numbers were retained and their parameter values were extracted for comparison. These parameter values were used to run simulations in areas where dengue virus transmission does not occur or where dengue fever case data was unavailable. Additional model runs were performed to reveal how changes in climate or parameter values could alter transmission risk along the transect. The relative influence of climate variability and model parameters on dengue virus transmission is assessed to help public health workers prepare location specific infection prevention strategies.

Morin, Cory↗

Climate effects of future aerosol reductions for achieving carbon neutrality in China

To limit the global warming to 1.5 °C above pre-industrial levels, beyond which the most dangerous impacts of climate change will occur, achieving carbon neutrality by the mid-21st century is essential. As the largest developing country and a significant contributor to carbon dioxide (CO 2 ) emissions, China has announced its ambitious climate commitment to pass carbon peak before 2030 and to achieve carbon neutrality by 2060. Both climate policies and regional clean air actions have been implemented for reductions in fossil fuel emissions, including the emissions of short-lived aerosols and precursors. Furthermore, in the context of pursuing carbon neutrality in China, aerosol reductions due to clean air actions and pollution control policies are very likely to have a great impact on climate. In this study, climate effects of aerosol reductions due to China's clean air actions under localized future emission scenarios are investigated using the Community Earth System Model version 1 (CESM1). Fully coupled and atmosphere-only experiments in years of carbon peak (2030) and carbon neutrality (2060) are performed with anthropogenic emissions of aerosols and precursors under “Current-goals” (Current) and “Carbon-Neutral” (Neutral) scenarios from the Dynamic Projection for Emissions in China (DPEC) model that consider socioeconomic development, climate policy, and pollution control actions (Figs. S1–S4 online). In addition, a present-day emissions simulation (PD) is conducted as the reference case (Supplementary materials online). Another sensitivity simulation is also conducted, with black carbon (BC) emissions set to follow the Neutral experiment but other emissions kept at the present-day levels (Neutral_BC), to quantify the relative roles of reducing strongly absorbing BC and other aerosols on future climate towards carbon neutrality. Greenhouse gases concentrations are kept at the 2015 levels in all simulations (Text S1 online).

Yang, Yang↗

EMT-TS Hybrid Simulation for Large Power Grids Considering IBR-Driven Dynamics

The escalating integration of inverter-based resources (IBRs) poses new challenges to power systems by introducing fast dynamics with higher frequencies, which may need to be simulated by an electromagnetic transient (EMT) program. As an alternative to conducting EMT simulations for the entire system, which is typically time consuming, hybrid simulation between EMT and phasor-domain transient stability (TS) can greatly reduce the computational burden while preserving the detailed fast dynamics in the EMT zone. This paper establishes an EMT-TS hybrid simulation platform using open-source tools, specifically ParaEMT, GridPACK, and HELICS, which are the EMT simulator, TS simulator, and interface framework, respectively. Case studies on the 240-bus Western Electricity Coordinating Council (WECC) system demonstrate that the developed ParaEMT-HELICS-GridPACK hybrid simulator can accurately capture both slow electromechanical and fast IBR-driven dynamics with a 2.4x speedup.

electromagnetic transient↗

Medium-Scale Methanol Pool Fire Model Validation

In this work, medium scale (30 cm diameter) methanol pool fires were simulated using the latest fire modeling suite implemented in Sierra/Fuego, a low Mach number multiphysics reacting flow code. The sensitivity of model outputs to various model parameters was studied with the objective of providing model validation. This work also assesses model performance relative to other recently published large eddy simulations (LES) of the same validation case. Two pool surface boundary conditions were simulated. The first was a prescribed fuel mass flux and the second used an algorithm to predict mass flux based on a mass and energy balance at the fuel surface. Gray gas radiation model parameters (absorption coefficients and gas radiation sources) were varied to assess radiant heat losses to the surroundings and pool surface. The radiation model was calibrated by comparing the simulated radiant fraction of the plume to experimental data. The effects of mesh resolution were also quantified starting with a grid resolution representative of engineering type fire calculations and then uniformly refining that mesh in the plume region. Simulation data were compared to experimental data collected at the University of Waterloo and the National Institute of Standards and Technology (NIST). Validation data included plume temperature, radial and axial velocities, velocity temperature turbulent correlations, velocity velocity turbulent correlations, radiant and convective heat fluxes to the pool surface, and plume radiant fraction. Additional analyses were performed in the pool boundary layer to assess simulated flame anchoring and the effect on convective heat fluxes. This work assesses the capability of the latest Fuego physics and chemistry model suite and provides additional insight into pool fire modeling for nonluminous, non-sooting flames.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Impacts of Estimated Plume Rise on PM 2.5 Exceedance Prediction During Extreme Wildfire Events: A Comparison of Three Schemes (Briggs, Freitas, and Sofiev)

Plume height plays a vital role in wildfire smoke dispersion and the subsequent effects on air quality and human health. In this study, we assess the impact of different plume rise schemes on predicting the dispersion of wildfire air pollution and the exceedances of the National Ambient Air Quality Standards (NAAQS) for fine particulate matter (PM 2.5 ) during the 2020 western United States wildfire season. Three widely used plume rise schemes (Briggs, 1969; Freitas et al., 2007; Sofiev et al., 2012) are compared within the Community Multiscale Air Quality (CMAQ) modeling framework. The plume heights simulated by these schemes are comparable to the aerosol height observed by the Multi-angle Imaging SpectroRadiometer (MISR) and Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO). The performance of the simulations with these schemes varies by fire case and weather conditions. On average, simulations with higher plume injection heights predict lower aerosol optical depth (AOD) and surface PM 2.5 concentrations near the source region but higher AOD and PM 2.5 in downwind regions due to the faster spread of the smoke plume once ejected. The 2-month mean AOD difference caused by different plume rise schemes is approximately 20 %–30 % near the source regions and 5 %–10 % in the downwind regions. Thick smoke blocks sunlight and suppresses photochemical reactions in areas with high AOD. The surface PM 2.5 difference reaches 70 % on the West Coast of the USA, and the difference is lower than 15 % in the downwind regions. Moreover, the plume injection height affects pollution exceedance (>35 µg m−3) predictions. Higher plume heights generally produce larger downwind PM 2.5 exceedance areas. The PM 2.5 exceedance areas predicted by the three schemes largely overlap, suggesting that all schemes perform similarly during large wildfire events when the predicted concentrations are well above the exceedance threshold. At the edges of the smoke plumes, however, there are noticeable differences in the PM 2.5 concentration and predicted PM 2.5 exceedance region. For the whole period of study, the difference in the total number of exceedance days could be as large as 20 d in northern California and 4 d in the downwind regions. This disagreement among the PM 2.5 exceedance forecasts may affect key decision-making regarding early warning of extreme air pollution episodes at local levels during large wildfire events.

Yunyao Li↗

URANIUM HOLDUP MONITORING WITH COMPTON IMAGING AS FUNCTION OF DEPTH AND MASS

Mass quantification using standoff gamma detection is a logical progression for MC&A applications in uranium processing facilities. While obtaining a ground truth from physical collection and weighing of masses is ideal and routine, it poses a logistical problem from both a standpoint of ALARA programs as well as interruptions to operations. The standard approach is the Generalized Geometry in Holdup (GGH) method that can apply for many cases where geometries of the holdup match simulated models; however, there are many cases where a model may not exist or a feature of the holdup geometry is such that it does not match a model particularly well.1 Gamma imaging at standoff distances and outside process equipment has potential to address these shortcomings and improve mass estimations while reducing personnel exposure; however, algorithms still require a significant degree of benchmarking to build sufficient confidence in algorithms before displacing GGH in some procedures. This paper presents the creation of a set of laboratory standards of depleted U3O8 in fixed size containers covering a range of masses and effective thicknesses up to the infinite thickness for the 185-keV line. These blocks are identically sized aluminum boxes loaded such that a thin, flat face is placed towards the imager, while the back side contains closed cell foam to apply pressure to keep the powder in a fixed position at the front face as the sources are manipulated for measurements. Included is an analysis of the mass estimation using a commercially available Compton imaging instrument where multiple geometries are created showing limitations in current methodologies.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Applications of artificial neural nets in structural mechanics

A brief introduction to the fundamental of Neural Nets is given, followed by two applications in structural optimization. In the first case, the feasibility of simulating with neural nets the many structural analyses performed during optimization iterations was studied. In the second case, the concept of using neural nets to capture design expertise was studied.

Berke, Laszlo↗

Applications of artificial neural nets in structural mechanics

A brief introduction to the fundamental of Neural Nets is given, followed by two applications in structural optimization. In the first case, the feasibility of simulating with neural nets the many structural analyses performed during optimization iterations was studied. In the second case, the concept of using neural nets to capture design expertise was studied.

Berke, L.↗

MOSCATO Solver Development and Integration Plan

During FY21, we conducted ongoing development work for the MOSCATO (Molten Salt Chemistry and Transport) solver. The code development work primarily consisted of transitioning capabilities from the original version of the solver, which was written in OpenFOAM, into Nek5000. In doing so, a fast, highly parallelizable solver was created that is capable of complex chemistry and corrosion simulations for engineering-scale molten salt systems. The Nek5000 version of MOSCATO is now fully featured and capable of higher-fidelity simulations than were previously possible. Demonstration cases including a thermal convection loop have been simulated to test these new capabilities. Although capable of large-scale simulations, MOSCATO is not well-suited to parametric studies of complete reactor geometries. These types of simulations are instead better handled by reduced-order modeling codes such as ORNL’s Mole code. Reduced-order simulation tools like Mole, however, are dependent on high fidelity correlations to account for complex, coupled three-dimensional phenomena that they do not directly simulate. Tools such as MOSCATO must therefore be used to create these correlations, as suitable empirical relationships are not available for most molten salt systems. Toward that end, we used the Nek-derived version of MOSCATO to create new mass transfer correlations for three relevant cases including tubular, tube bank, and subchannel geometries. These new correlations are more accurate than any existing ones and can be readily integrated into any reduced-order modeling tools that are targeting full-scale MSR simulations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A Study of Cost-Saving Potential of Load Flexibility Measures in Grid-Interactive Multifamily Buildings

With recent advances in smart technologies, more and more smart devices are penetrating the residential and commercial buildings market. The introduction of these smart devices is also helping IoT companies emerge with load aggregator roles in the sector. With more utility companies on the track of supporting OpenADR protocols, the aggregators could play a significant role in providing load flexibilities by automatically responding to demand response (DR) events and coordinating load flexibility measures between customers. This would benefit utility companies by reducing stress on the grid during critical peak demand hours as well as customers by allowing them to utilize utility rate structures advantageous to those able to reduce electric usage during high-demand hours. This study evaluates cost and energy savings from adopting multiple load flexibility measures in multifamily buildings. Combinations of different load flexibility measures, including space temperature floating, light dimming, automatic window shading, and water heater temperature floating, are considered. The simulations are performed using OpenStudio®, an open-source U.S. Department of Energy (DOE) simulation platform. For the case study, we used a midrise apartment building with weather conditions from Denver, Colorado. To compare climate zone differences, simulations were also performed for Los Angeles, California, and Chicago, Illinois. Initial results indicate that the application of automated load flexibility measures without careful consideration of dispatching strategies and DR program enrollments could significantly affect the savings. To get meaningful cost savings, aggregators need to encourage tenant awareness to curtail energy usage through occupant behavior in addition to dispatching automatic load flexibility measures. The outcomes from this study are believed to help load aggregators understand the risks and benefits of load flexibility opportunities.

building energy modeling↗

Numerical study of magnetic island coalescence using magnetohydrodynamics with adaptively embedded particle-in-cell model

Collisionless magnetic reconnection typically requires kinetic treatment that is, in general, computationally expensive compared to fluid-based models. In this study, we use the magnetohydrodynamics with an adaptively embedded particle-in-cell (MHD-AEPIC) model to study the interaction of two magnetic flux ropes. This innovative model embeds one or more adaptive PIC regions into a global MHD simulation domain such that the kinetic treatment is only applied in regions where the kinetic physics is prominent. We compare the simulation results among three cases: (1) MHD with adaptively embedded PIC regions, (2) MHD with statically (or fixed) embedded PIC regions, and (3) a full PIC simulation. The comparison yields good agreement when analyzing their reconnection rates and magnetic island separations as well as the ion pressure tensor elements and ion agyrotropy. In order to reach good agreement among the three cases, large adaptive PIC regions are needed within the MHD domain, which indicates that the magnetic island coalescence problem is highly kinetic in nature, where the coupling between the macro-scale MHD and micro-scale kinetic physics is important.

79 ASTRONOMY AND ASTROPHYSICS↗

Check-Cases for Verification of 6-Degree-of-Freedom Flight Vehicle Simulations

The rise of innovative unmanned aeronautical systems and the emergence of commercial space activities have resulted in a number of relatively new aerospace organizations that are designing innovative systems and solutions. These organizations use a variety of commercial off-the-shelf and in-house-developed simulation and analysis tools including 6-degree-of-freedom (6-DOF) flight simulation tools. The increased affordability of computing capability has made highfidelity flight simulation practical for all participants. Verification of the tools' equations-of-motion and environment models (e.g., atmosphere, gravitation, and geodesy) is desirable to assure accuracy of results. However, aside from simple textbook examples, minimal verification data exists in open literature for 6-DOF flight simulation problems. This assessment compared multiple solution trajectories to a set of verification check-cases that covered atmospheric and exo-atmospheric (i.e., orbital) flight. Each scenario consisted of predefined flight vehicles, initial conditions, and maneuvers. These scenarios were implemented and executed in a variety of analytical and real-time simulation tools. This tool-set included simulation tools in a variety of programming languages based on modified flat-Earth, round- Earth, and rotating oblate spheroidal Earth geodesy and gravitation models, and independently derived equations-of-motion and propagation techniques. The resulting simulated parameter trajectories were compared by over-plotting and difference-plotting to yield a family of solutions. In total, seven simulation tools were exercised.

Murri, Daniel G.↗

Coordinated Modeling of Electric Grid and Natural Gas Network Operations

Presentation based upon full report of a Colorado case study and coordination framework, which is available at https://www.nrel.gov/docs/fy20osti/77096.pdf. JISEA analysts Brian Sergi, Omar Guerra, and Bri-Mathias Hodge will present on recent JISEA work on coordination between the natural gas and electricity sectors. Power and gas are becoming increasingly interdependent but weren't designed to function together. How does greater coordination impact system operations with different levels of solar and wind penetrations? This free JISEA presentation will take place in a Webinar on Tuesday, May 11, 2021 at 12 p.m.: Learn how greater coordination between the natural gas and electricity sectors impacts system operations with different levels of solar and wind penetrations. Presentation information includes: Historical and projected data for natural gas consumption and power generation in the United States as demonstrated by 2014 East Coast Polar Vortex, 2021 Texas Winter Storm Uri. Showing coupling points, FERC identified need for better coordination, different levels of coordination (decision making and optimization control), coordination framework. Case study on Colorado Front Range, ramping requirements and gas nominations, real-time dispatch (June and December scenarios), results of analysis of impacts of coordination on unserved load. Total real-time gas offtakes by node, impacts on unserved gas, effect on CO2 emissions. Discussions on clusions drawn from Colorado Case study expanded co-simulation via HELICS, objectives of the HELICS+ natural gas use case, modeling of hydrogen blending, techno economic assessment of blending, and blending impacts on energy content and pressure.

analysis↗

Teleoperation of Construction Robots Using Human Interaction: Case Study Examples: Preprint

The integration of advanced teleoperation methods in construction robotics has the potential to enhance productivity and safety in demanding environments. Among different methods, hand gestures have emerged as an intuitive approach, enabling seamless communication between human operators and robotic systems. This study demonstrates the use of hand gestures for controlling a model robotic excavator and a virtual robotic arm in construction applications. Two case study examples are presented. In the first case study, a series of hand gestures are performed to control a model robotic excavator's bucket and arm movements. The second case study focuses on using hand gestures to control a virtual robotic arm to perform a material handling task in a simulated environment. In both cases, the performed hand gestures are all successfully captured and interpreted by the gesture-based teleoperation method to conduct the corresponding tasks, which demonstrates the application of the gesture-based control method for various types of construction robots. Also, the limitations of the gesture-based teleoperation method (e.g., recognition delay, sensitivity to motion outliers) are discussed. Future work will focus on developing a digital twin system to support gesture-based teleoperation for a wider range of construction robots.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Molecular dynamics and experimental study on the adhesion mechanism of polyvinyl alcohol (PVA) fiber in alkali-activated slag/fly ash

This paper aims to study the adhesion mechanism of polyvinyl alcohol (PVA) fiber within alkali-activated slag/fly ash (AASF) matrix using molecular dynamics (MD) simulation in combination with systematic experimental characterization. The adhesion of PVA to C-(N-)A-S-H gel with different Ca/(Si+Al) and Al/Si ratios was modeled using MD simulation, with the related adsorption enthalpy calculated and the adhesion mechanism explored. The experimentally attained chemical bonding energy of PVA fiber in AASF coincides well with the simulation results. In both cases, the adhesion enhances primarily with increasing Ca/(Si+Al) ratio of C-(N-)A-S-H gel. Additionally, MD simulation indicates preferential element distributions of Ca around PVA molecule, which was confirmed experimentally by the detection of the Ca-rich C-(N-)A-S-H gel in the interfacial transition zone (ITZ). }This study provides further insights into the adhesion mechanism of PVA fiber to C-(N-)A-S-H gel formed in AASF, which is particularly valuable for the future development of PVA-based high-performance alkali-activated composites.

36 MATERIALS SCIENCE↗

Computational modeling of graphite degradation in molten salt reactors: Role of infiltration

Molten salt reactors (MSRs) often employ graphite as a moderator and reflector. An important challenge for deploying graphite in these reactors is that, due to limited experimental data, our understanding of graphite’s structural integrity in molten salt environments remains incomplete. Here, this study addresses heat generation from fuel-bearing salt that has infiltrated open pores in the graphite, driven primarily by pressure differentials. This is one of multiple identified physical and chemical mechanisms through which molten salt could potentially degrade graphite. Thermally driven stresses are quantified using the Molten-Salt Reactor Experiment (MSRE) graphite moderator elements as a case study. Finite element simulations predict stress distributions at varying infiltration levels, indicating that thermal stresses increase with higher infiltration. Rare-event simulations using the parallel subset simulation framework identify the combinations and corresponding ranges of input parameters that lead to stresses above a specified threshold. In particular, combinations involving high infiltration amounts, high power density, and low thermal conductivity tend to induce the highest stresses. Under the inputs and assumptions considered in this work, the magnitudes of the thermally driven stresses are quite low, with a very low likelihood of causing failure due to exceeding the graphite’s tensile strength. Additionally, rare-event simulations were performed for two more scenarios: a scaled-up moderator geometry and a localized hotspot in the original geometry. Both cases resulted in increased susceptibility to failure, though not to a detrimental extent. Furthermore, the combined effects of irradiation and infiltration-induced thermal stresses were evaluated. The results showed that thermal stresses from infiltration were negligible compared to those caused by irradiation. The findings of such a study are inherently component-specific, but the methodology presented here could be used for similar assessments of salt-infiltration effects in other graphite components.

36 - MATERIALS SCIENCE↗

Large-Eddy Simulations of a Single-Injector Cooling Flow Using the High-Order Flux Reconstruction Method

A single-injector cooling flow into a heated crossflow was used as a validation case for large eddy simulations (LES) from two high-order CFD codes with comparisons to a state-of-the-art Reynolds averaged Navier-Stokes (RANS) turbulence model. The focus of this paper is on the validation of LES for one of the high-order CFD codes, namely the GFR (Glenn Flux Reconstruction) code that is being developed at NASA Glenn Research Center. Fourth-order LES were performed for the single-injector cooling flow configuration at blowing ratios of both 1.0 and 2.0, and a fifth-order LES was also performed with a blowing ratio of 1.0. The GFR simulations agree very well with the experiment mean quantities and reasonably well with the experiment turbulence quantities. The LES solutions from GFR and the other high-order CFD code, FDL3DI, agreed very closely for nearly every flow quantity examined, despite the fact that these two codes use completely different numerical methods, different grid geometries and domains, and different inflow boundary conditions. Finally, both LES show a significant accuracy improvement over RANS turbulence models for flows with large temperature gradients, where the standard gradient diffusion approximation is insufficient for temperature predictions.

Large Eddy Simulations↗