Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulation modelling”

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

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

At least 1,117 records · Page 62

Comparison of Deterministic and Statistical Models for Water Quality Compliance Forecasting in the San Joaquin River Basin, California

Model selection for water quality forecasting depends on many factors including analyst expertise and cost, stakeholder involvement and expected performance. Water quality forecasting in arid river basins is especially challenging given the importance of protecting beneficial uses in these environments and the livelihood of agricultural communities. In the agriculture-dominated San Joaquin River Basin of California, real-time salinity management (RTSM) is a state-sanctioned program that helps to maximize allowable salt export while protecting existing basin beneficial uses of water supply. The RTSM strategy supplants the federal total maximum daily load (TMDL) approach that could impose fines associated with exceedances of monthly and annual salt load allocations of up to $1 million per year based on average year hydrology and salt load export limits. The essential components of the current program include the establishment of telemetered sensor networks, a web-based information system for sharing data, a basin-scale salt load assimilative capacity forecasting model and institutional entities tasked with performing weekly forecasts of river salt assimilative capacity and scheduling west-side drainage export of salt loads. Web-based information portals have been developed to share model input data and salt assimilative capacity forecasts together with increasing stakeholder awareness and involvement in water quality resource management activities in the river basin. Two modeling approaches have been developed simultaneously. The first relies on a statistical analysis of the relationship between flow and salt concentration at three compliance monitoring sites and the use of these regression relationships for forecasting. The second salt load forecasting approach is a customized application of the Watershed Analysis Risk Management Framework (WARMF), a watershed water quality simulation model that has been configured to estimate daily river salt assimilative capacity and to provide decision support for real-time salinity management at the watershed level. Analysis of the results from both model-based forecasting approaches over a period of five years shows that the regression-based forecasting model, run daily Monday to Friday each week, provided marginally better performance. However, the regression-based forecasting model assumes the same general relationship between flow and salinity which breaks down during extreme weather events such as droughts when water allocation cutbacks among stakeholders are not evenly distributed across the basin. A recent test case shows the utility of both models in dealing with an exceedance event at one compliance monitoring site recently introduced in 2020.

54 ENVIRONMENTAL SCIENCES↗

An evolving Coupled Model Intercomparison Project phase 7 (CMIP7) and Fast Track in support of future climate assessment

The Coupled Model Intercomparison Project (CMIP) coordinates community-based efforts to answer key and timely climate science questions, facilitate delivery of relevant multi-model simulations through shared infrastructure, and support national and international climate assessments. Generations of CMIP have evolved through extensive community engagement from punctuated phasing into more continuous support for the design of experimental protocols, infrastructure for data publication and access, and public delivery of climate information. We identify four fundamental research questions motivating a seventh phase of coupled model intercomparison relating to patterns of sea surface temperature change, changing weather, the water–carbon–climate nexus, and tipping points. Key CMIP7 advances include an expansion of baseline experiments, a focus on CO 2 -emissions-driven experiments, sustained support for community MIPs, periodic updating of historical forcings and diagnostics requests, and a collection of prioritized experiments, or the “Assessment Fast Track”, drawn from community MIPs to support climate research, assessment, and service goals across prediction and projection, characterization, attribution, and process understanding.

Environmental sciences↗

Mathematical analysis of laboratory microbial experiments demonstrating deterministic chaotic dynamics

Presented is a system of four ordinary differential equations and a mathematical analysis of microbiological experiments in a four-component chemostat—nutrient n, rods r, cocci c, and predators p. The analysis is consistent with the conclusion that previous experiments produced features of deterministic chaotic and classical dynamics depending on dilution rate. The surrogate model incorporates as much experimental detail as possible, but necessarily contains unmeasured parameters. The objective is to understand better the differences between model simulations and experimental results in complex microbial populations. The key methodology for simulation of chaotic dynamics, consistent with the measured dilution rate and microbial volume averages, was to cause the preference of p for r vs. c to vary with the r and c concentrations, to make r more competitive for nutrient than c, and to recycle some dying p biomass, leading to a modified version of the Monod kinetics model. Our mathematical model demonstrated that the occurrence of chaotic dynamics requires a predator, p, preference for r versus c to increase significantly with increases in r and c populations. Also included is a discussion of several generalizations of the existing model and a possible involvement of the minimum energy dissipation principle. This principle appears fundamental to thermodynamic systems including living systems. Several new experiments are suggested.

59 BASIC BIOLOGICAL SCIENCES↗

Quantifying UAS Observation Error Variance Used in Data Assimilation Systems and Its Impact on Predictive Skill

Observation error determines the weights of the observations and background state used in data assimilation to generate analyses. Quantifying observation error is critical for the optimal assimilation of observational data sets. Uncrewed Aircraft System (UAS) observations have shown potential benefits in filling observational gaps in the lower atmosphere; however, characterization of their error characteristics has been limited. To optimize the use of UAS observations in numerical weather prediction, UAS observation error is estimated based on the 3‐cornered hat diagnostic approach which uses three independent estimates of the atmospheric state. This approach is applied to data from the 2018 Lower Atmospheric Profiling Studies at Elevation‐a Remotely‐piloted Aircraft Team Experiment field campaign using collocated UAS and rawinsonde observations along with output from a set of convection‐permitting model simulations. The estimated observation error values for UAS temperature, wind, and relative humidity measurements were found to be only weakly dependent on height AGL with mean values equal to 0.5°C, 0.8 m s −1 , and 3%, respectively. Only the newly estimated observation error for temperature differed from that previously used to assimilate commercial aircraft observations into global models (1.0°C). However, using this reduced temperature observation error produced more accurate mesoscale analyses and forecasts of both terrain‐driven flows and convection initiation generated by colliding outflow boundaries within the San Luis Valley of Colorado.

54 ENVIRONMENTAL SCIENCES↗

Sparse Linear Solvers for Large-scale Electromagnetic Transient Simulations

Linear solvers form the basis for electromagnetic transient (EMT) simulations. There is a need to speed up EMT simulations as larger regions are analyzed using EMT simulations. For the same, the performance of linear solvers plays an important role. Exploiting the sparsity of the matrices generated in EMT simulations could assist with speed-up. Scalability is also crucial as power grids expand, demanding solutions capable of accommodating the increasing system size. Recent studies from the North American Electric Reliability Corporation (NERC) increasingly emphasize that EMT simulation models of the power grid will grow larger with the inclusion of power electronics components. Parallelisms in sparsity patterns exploit modern central processing units (CPUs), multi-core CPUs, and graphics processing units (GPUs) architectures in sparse solver designs. Therefore, this paper explores publicly available existing linear solvers and investigates their efficiency in large-scale power grid simulations. A large-scale power grid is developed by increasing the size of the IEEE 39 bus test system to up to 39000 bus systems.

Hsu, Kuan-Chieh↗

Improving Fission Products at CARIBU: Near Field Detection (Q4/FY21 Quarterly Progress Report)

Since the last quarterly, we have been working on improving the GEANT4 simulation model, which we currently use for the 4-pi gas counter efficiency calculations, to streamline the work when calculating multi-level decay properties. We have also been preparing for the next branching-ratio measurement, which will be for 161 Tb. Because travel is still limited and restricted, we have worked with a summer student on performing simulations to estimate our detector’s efficiency for different thresholds.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Analysis of the HAE activity in the TJ-II stellarator using a Landau closure model

The aim of this study is to analyze the stability of helical Alfvén eigenmodes (HAEs) in TJ-II discharges and the stabilizing effect of the energetic particles generated by the neutral beam injector (NBI) on pressure gradient-driven modes (PGDMs). HAE and PGDM stability is studied using the linear version of the gyro-fluid code FAR3d and the continuous structure by the STELLGAP code. First, Alfvén eigenmode (AE) and PGDM activity observed in the experiments is reproduced by the simulations, identifying unstable m/n = 4/7 − 2/3 and 7/12 − 5/8 HAEs triggered around ρ = 0.66 showing a frequency of 209 and 204 kHz, respectively, as well as 5/3 PGDM. Next, a parametric study is performed with respect to the thermal ion density and iota profile in the middle-outer plasma region to verify the robustness of the simulation results with respect to the uncertainty of experimental profiles. The analysis confirms that experimental uncertainty does not cause large deviations in the simulation results, showing the destabilization of the same HAEs for all the configurations tested. The simulations also indicate the decay of the 5/3 PGDM growth rate as the energetic particle (EP) population in the plasma increases, consistent with the experiment. Stability analysis of the n = 3, 7, 11, n = 5, 9, 13, n = 6, 10, 14, and n = 8, 12 helical families is performed with respect to the NBI operational regime for different EP energies, β as well as deposition profiles. The most unstable configuration is the radially localized on-axis NBI operation (stiff EP density profile gradients nearby the magnetic axis). Using the simulation model that reproduces the observed Alfvén activity, we extend the study to analyze NBI performance within a theoretical framework. It shows that increasing NBI voltage (which raises EP energy) leads to a degradation in NBI performance for a given power (related to EP β and their density). To achieve better NBI operation, higher voltage must be balanced with lower injection power, ensuring stable AEs while keeping the same EP β.

Plasma confinement↗

Identification of the molecular-weight growth reaction network in counterflow flames of the C 3 H 4 isomers allene and propyne

The reaction networks responsible for aromatics formation in counterflow flames of the C 3 H 4 isomers allene and propyne are identified through a combined experimental and modeling study. Mole fraction profiles of near-atmospheric pressure (933 mbar) diffusion flames fueled by the C 3 H 4 isomers are analyzed by means of a newly assembled, chemically detailed kinetic mechanism. The experiment consists of a counterflow burner system that is coupled to a high-resolution time-of-flight molecular-beam mass spectrometer with single-photon ionization via synchrotron-generated vacuum-ultraviolet photons. Flame-sampled, mass-specific photoionization efficiency curves are used to identify the presence of aliphatically substituted aromatic species in addition to the commonly considered pericondensed ring structures. The new mechanism describes the formation and growth of aromatics through repetitive sequences of radical–radical and radical–molecule reactions that include C 1 6 intermediates. Higher concentrations of aromatic species are observed in the allene flame and the new mechanism captures the observed experimental trends very accurately. The results indicate the importance of the aliphatically substituted aromatics and of ring-enlargement reactions for the growth reactions. According to the model simulations, radical+radical recombination and PAH-radical+molecule reactions play an important role in PAH growth.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A new metrics framework for quantifying and intercomparing atmospheric rivers in observations, reanalyses, and climate models

We present a new atmospheric river (AR) analysis and benchmarking tool, namely Atmospheric River Metrics Package (ARMP). It includes a suite of new AR metrics that are designed for quick analysis of AR characteristics via statistics in gridded climate datasets such as model output and reanalysis. This package can be used for climate model evaluation in comparison with reanalysis and observational products. Integrated metrics such as mean bias and spatial pattern correlation are efficient for diagnosing systematic AR biases in climate models. For example, the package identifies the fact that, in CMIP5 and CMIP6 (Coupled Model Intercomparison Project Phases 5 and 6) models, AR tracks in the South Atlantic are positioned farther poleward compared to ERA5 reanalysis, while in the South Pacific, tracks are generally biased towards the Equator. For the landfalling AR peak season, we find that most climate models simulate a completely opposite seasonal cycle over western Africa. This tool can also be used for identifying and characterizing structural differences among different AR detectors (ARDTs). For example, ARs detected with the Mundhenk algorithm exhibit systematically larger size, width, and length compared to the TempestExtremes (TE) method. The AR metrics developed from this work can be routinely applied for model benchmarking and during the development cycle to trace performance evolution across model versions or generations and set objective targets for the improvement of models. They can also be used by operational centers to perform near-real-time climate and extreme event impact assessments as part of their forecast cycle.

58 GEOSCIENCES↗

Emulation With Uncertainty Quantification of Regional Sea‐Level Change Caused by the Antarctic Ice Sheet

Abstract Projecting regional sea‐level change under various climate‐change scenarios typically involves running forward simulations of the Earth's gravitational, rotational and deformational (GRD) response to ice‐mass change, which requires substantial computational cost if applied to probabilistic frameworks requiring thousands to millions of samples. Here we build emulators of regional sea‐level change at 27 coastal locations, due to the GRD effects associated with future Antarctic Ice Sheet mass change over the 21st century. The emulators are evaluated against a numerical sea‐level model applied to an ensemble of ice‐sheet model simulations of the Antarctic Ice Sheet through 2100. We build a physics‐based emulator using a recent sensitivity kernel approach and compare it to machine learning based emulators (neural network and conditional variational autoencoder methods). In order to quantify uncertainty, we derive well‐calibrated prediction intervals for regional sea‐level change via split‐conformal inference and linear regression, and show that Monte Carlo dropout does not yield well‐calibrated uncertainties in this instance. We also demonstrate substantial gains in computational efficiency using both the physics‐based emulator and neural networks in comparison to the numerical model for the complete regional sea‐level solution. Overall, we find the physics‐based emulator modestly outperforms the machine learning emulators for this problem.

58 GEOSCIENCES↗

Control Design for a Marine Hydrokinetic Cycloturbine Vehicle

Abstract Marine Hydrokinetic (MHK) cycloturbines generate sustainable power by exploiting tidal currents. By powering the turbines and using pitching foils for control, a vehicle comprised of MHK cycloturbines also has the ability to station keep and maneuver. The vehicle consists of four counter-rotating cycloturbines, with hydrofoils oriented perpendicular to the flow in a paddlewheel configuration. Lift and drag generated from these foils sum together to produce thrust. An experimentally tuned simulation model that solves the six-degrees-of-freedom rigid body equations of motion for the MHK vehicle subject to hydrodynamic, hydrostatic, and propulsive forces is used to aid the design of vehicle controllers. Global feedback controllers are initially designed by applying classical control methods to an approximate linear model of the system dynamics. A higher performing nonlinear controller is designed using the nonlinear dynamic inversion (NDI) method. NDI accounts for the nonlinearities of the MHK system and therefore is suitable for a wide range of operating conditions. The response of the classical and NDI controllers to speed, depth, roll, pitch, and yaw commands are evaluated and compared in simulation. The classical controller outperforms the NDI controller for small amplitude maneuvers, although the degradation with NDI is minor. However, in the nonlinear operating regime the NDI controller outperforms the classical controller and the classical controller exhibits instability.

Automation & Control Systems↗

A global study of hygroscopicity-driven light-scattering enhancement in the context of other in situ aerosol optical properties

The scattering and backscattering enhancement factors (ƒ(RH) and ƒ b (RH)) describe how aerosol particle light scattering and backscattering, respectively, change with relative humidity (RH). They are important parameters in estimating direct aerosol radiative forcing (DARF). In this study we use the dataset presented in Burgos et al. (2019) that compiles ƒ(RH) and ƒ b (RH) measurements at three wavelengths (i.e., 450, 550 and 700 nm) performed with tandem nephelometer systems at multiple sites around the world. We present an overview of ƒ(RH) and ƒ b (RH) based on both long-term and campaign observations from 23 sites representing a range of aerosol types. The scattering enhancement shows a strong variability from site to site, with no clear pattern with respect to the total scattering coefficient. In general, higher ƒ(RH) is observed at Arctic and marine sites, while lower values are found at urban and desert sites, although a consistent pattern as a function of site type is not observed. The backscattering enhancement ƒ b (RH) is consistently lower than ƒ(RH) at all sites, with the difference between ƒ(RH) and ƒ b (RH) increasing for aerosol with higher ƒ(RH). This is consistent with Mie theory, which predicts higher enhancement of the light scattering in the forward than in the backward direction as the particle takes up water. Our results show that the scattering enhancement is higher for PM 1 than PM 10 at most sites, which is also supported by theory due to the change in scattering efficiency with the size parameter that relates particle size and the wavelength of incident light. At marine-influenced sites this difference is enhanced when coarse particles (likely sea salt) predominate. For most sites, ƒ(RH) is observed to increase with increasing wavelength, except at sites with a known dust influence where the spectral dependence of ƒ(RH) is found to be low or even exhibit the opposite pattern. The impact of RH on aerosol properties used to calculate radiative forcing (e.g., single-scattering albedo, ω 0 , and backscattered fraction, b) is evaluated. The single-scattering albedo generally increases with RH, while b decreases. The net effect of aerosol hygroscopicity on radiative forcing efficiency (RFE) is an increase in the absolute forcing effect (negative sign) by a factor of up to 4 at RH = 90 % compared to dry conditions (RH < 40 %). Because of the scarcity of scattering enhancement measurements, an attempt was made to use other more commonly available aerosol parameters (i.e., ω 0 and scattering Ångström exponent, α sp ) to parameterize ƒ(RH). The majority of sites (75 %) showed a consistent trend with ω 0 (higher ƒ(RH = 85 %) for higher ω 0 ), while no clear pattern was observed between ƒ(RH = 85 %) and αsp. This suggests that aerosol ω 0 is more promising than α sp as a surrogate for the scattering enhancement factor, although neither parameter is ideal. Nonetheless, the qualitative relationship observed between ω 0 and ƒ(RH) could serve as a constraint on global model simulations.

54 ENVIRONMENTAL SCIENCES↗

Interpretable Machine Learning for Characterizing Electric Vehicle Charging Behavior: Insights from Real-World Data

As electric vehicle (EV) adoption rises globally, concerns about the impact on aging electrical grids grow, particularly regarding the charging behavior of EV drivers. This study analyzes real-world driving and charging data from Ford battery electric vehicles (BEVs) collected between 2018 and 2019 to develop interpretable models that characterize charging behavior and quantify influencing factors. Prior research has relied on assumptions regarding driver behavior, often overlooking actual charging patterns. By employing generalized linear mixed models (GLMMs), this work offers insights into how various elements, such as next trip distance and state of charge (SOC), influence charging decisions. The dataset comprises over three million park-trip pairs from 1,997 vehicles, revealing that features related to driving behavior significantly dictate charging behavior, while infrastructure and regional factors have lesser impacts. The findings suggest that existing simulation models may oversimplify EV charging behavior assumptions. This work utilizes real-world EV driving and charging data to train interpretable models that describe charging behavior and quantify the factors most associated with how drivers use charging infrastructure. This research underscores the need for interpretable, data-driven methodologies to inform future EV infrastructure planning and grid management.

29 - ENERGY PLANNING, POLICY AND ECONOMY↗

Acquire and Test IDAES Framework Components

The cooperation under the Grid Modernization Initiative and the Applied Energy Tri-Laboratory umbrellas requires mutual understanding of capabilities and a common effort to identify gaps existing in the available computational modeling and analysis tools. In this work, we report acquisition and testing of the National Energy Technology Laboratory’s Institute for the Design of Advanced Energy Systems (IDAES) framework, particularly including both electric power systems production cost model simulator Prescient as well as system optimization IDAES itself. Further, we consider the methods in which developments in Prescient and IDAES for grid-energy simulation and optimization may be included in the Framework for Optimization of ResourCes and Economics (FORCE) ecosystem that has been under collaborative development by Idaho National Laboratory, Argonne National Laboratory, and Oak Ridge National Laboratory. We propose two integration methods to introduce Prescient and IDAES into existing FORCE automated workflows. One integration involves using Prescient as a market reduced-order model and IDAES models similar to how existing Modelica models are used. The other integration involves wrapping the Prescient-IDAES dispatch optimization process with the stochastic capacity optimization workflow in FORCE.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Optimizing radar scan strategies for tracking isolated deep convection using observing system simulation experiments

Abstract. Optimizing radar observation strategies is one of the most important considerations in pre-field campaign periods. This is especially true for isolated convective clouds that typically evolve faster than the observations captured by operational radar networks. This study investigates uncertainties in radar observations of the evolution of the microphysical and dynamical properties of isolated deep convective clouds developing in clean and polluted environments. It aims to optimize the radar observation strategy for deep convection through the use of high-spatiotemporal cloud-resolving model simulations, which resolve the evolution of individual convective cells every 1 min, coupled with a radar simulator and a cell tracking algorithm. The radar simulation settings are based on the Tracking Aerosol Convection Interactions ExpeRiment (TRACER) and Experiment of Sea Breeze Convection, Aerosols, Precipitation and Environment (ESCAPE) field campaigns held in the Houston, TX, area but are generalizable to other field campaigns focusing on isolated deep convection. Our analysis produces the following four outcomes. First, a 5–7 m s−1 median difference in maximum updrafts of tracked cells is shown between the clean and polluted simulations in the early stages of the cloud lifetimes. This demonstrates the importance of obtaining accurate estimates of vertical velocity from observations if aerosol impacts are to be properly resolved. Second, tracking of individual cells and using vertical cross section scanning every minute capture the evolution of precipitation particle number concentration and size represented by polarimetric observables better than the operational radar observations that update the volume scan every 5 min. This approach also improves multi-Doppler radar updraft retrievals above 5 km above ground level for regions with updraft velocities greater than 10 m s−1. Third, we propose an optimized strategy composed of cell tracking by quick (1–2 min) vertical cross section scans from more than one radar in addition to the operational volume scans. We also propose the use of a single-RHI (range height indicator) updraft retrieval technique for cells close to the radars, for which multi-Doppler radar retrievals are still challenging. Finally, increasing the number of deep convective cells sampled by such observations better represents the median maximum updraft evolution with sample sizes of more than 10 deep cells, which decreases the error associated with sampling the true population to less than 3 m s−1.

54 ENVIRONMENTAL SCIENCES↗

Radial–axial transport coordination enhances sugar translocation in the phloem vasculature of plants

Abstract Understanding mass transport of photosynthates in the phloem of plants is necessary for predicting plant carbon allocation, productivity, and responses to water and thermal stress. Several hypotheses about optimization of phloem structure and function and limitations of phloem transport under drought have been proposed and tested with models and anatomical data. However, the true impact of radial water exchange of phloem conduits with their surroundings on mass transport of photosynthates has not been addressed. Here, the physics of the Munch mechanism of sugar transport is re-evaluated to include local variations in viscosity resulting from the radial water exchange in two dimensions (axial and radial) using transient flow simulations. Model results show an increase in radial water exchange due to a decrease in sap viscosity leading to increased sugar front speed and axial mass transport across a wide range of phloem conduit lengths. This increase is around 40% for active loaders (e.g. crops) and around 20% for passive loaders (e.g. trees). Thus, sugar transport operates more efficiently than predicted by previous models that ignore these two effects. A faster front speed leads to higher phloem resiliency under drought because more sugar can be transported with a smaller pressure gradient.

58 GEOSCIENCES↗

Optimizing Price-Informed Operation of a Battery Storage System in an Office Building

New prescriptive efficiency requirements amended to US model energy codes historically have been evaluated using an average, blended electricity rate, which obscures demand charges. Post 2019, new measures can be evaluated using a representative time-of-use (TOU) tariff yet more sophisticated analysis methods are needed to consider a variety of TOU tariffs and assess price-informed control. To address these needs, this study couples building prototype simulation model with varying-in-sophistication battery storage operating strategies for different electricity TOU tariffs. The analysis compares the impact of operating a battery storage system following simpler rule-of-thumb methods versus a semi-optimized priced-informed heuristic approach. The investigation demonstrates that the heuristic approach results in greater electricity cost savings and is practical to implement.

Lei, Xuechen↗