Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SIMULATOR”

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,513 records · Page 84

Implementation of an Open-Scenario, Long-Term Space Debris Simulation Approach

This paper provides a status update on the implementation of a flexible, long-term space debris simulation approach. The motivation is to build a tool that can assess the long-term impact of various options for debris-remediation, including the LightForce space debris collision avoidance concept that diverts objects using photon pressure [9]. State-of-the-art simulation approaches that assess the long-term development of the debris environment use either completely statistical approaches, or they rely on large time steps on the order of several days if they simulate the positions of single objects over time. They cannot be easily adapted to investigate the impact of specific collision avoidance schemes or de-orbit schemes, because the efficiency of a collision avoidance maneuver can depend on various input parameters, including ground station positions and orbital and physical parameters of the objects involved in close encounters (conjunctions). Furthermore, maneuvers take place on timescales much smaller than days. For example, LightForce only changes the orbit of a certain object (aiming to reduce the probability of collision), but it does not remove entire objects or groups of objects. In the same sense, it is also not straightforward to compare specific de-orbit methods in regard to potential collision risks during a de-orbit maneuver. To gain flexibility in assessing interactions with objects, we implement a simulation that includes every tracked space object in Low Earth Orbit (LEO) and propagates all objects with high precision and variable time-steps as small as one second. It allows the assessment of the (potential) impact of physical or orbital changes to any object. The final goal is to employ a Monte Carlo approach to assess the debris evolution during the simulation time-frame of 100 years and to compare a baseline scenario to debris remediation scenarios or other scenarios of interest. To populate the initial simulation, we use the entire space-track object catalog in LEO. We then use a high precision propagator to propagate all objects over the entire simulation duration. If collisions are detected, the appropriate number of debris objects are created and inserted into the simulation framework. Depending on the scenario, further objects, e.g. due to new launches, can be added. At the end of the simulation, the total number of objects above a cut-off size and the number of detected collisions provide benchmark parameters for the comparison between scenarios. The simulation approach is computationally intensive as it involves tens of thousands of objects; hence we use a highly parallel approach employing up to a thousand cores on the NASA Pleiades supercomputer for a single run. This paper describes our simulation approach, the status of its implementation, the approach to developing scenarios and examples of first test runs.

laser↗

NASA's GMAO Atmospheric Motion Vectors Simulator: Description and Application to the MISTiC Winds Concept

An atmospheric wind vectors (AMVs) simulator was developed by NASA's GMAO to simulate observations from future satellite constellation concepts. The synthetic AMVs can then be used in OSSEs to estimate and quantify the potential added value of new observations to the present Earth observing system and, ultimately, the expected impact on the current weather forecasting skill. The GMAO AMV simulator is a tunable and flexible computer code that is able to simulate AMVs expected to be derived from different instruments and satellite orbit configurations. As a case study and example of the usefulness of this tool, the GMAO AMV simulator was used to simulate AMVs envisioned to be provided by the MISTiC Winds, a NASA mission concept consisting of a constellation of satellites equipped with infrared spectral midwave spectrometers, expected to provide high spatial and temporal resolution temperature and humidity soundings of the troposphere that can be used to derive AMVs from the tracking of clouds and water vapor features. The GMAO AMV simulator identifies trackable clouds and water vapor features in the G5NR and employs a probabilistic function to draw a subset of the identified trackable features. Before the simulator is applied to the MISTiC Winds concept, the simulator was calibrated to yield realistic observations counts and spatial distributions and validated considering as a proxy instrument to the MISTiC Winds the Himawari-8 Advanced Imager (AHI). The simulated AHI AMVs showed a close match with the real AHI AMVs in terms of observation counts and spatial distributions, showing that the GMAO AMVs simulator synthesizes AMVs observations with enough quality and realism to produce a response from the DAS equivalent to the one produced with real observations. When applied to the MISTiC Winds scanning points, it can be expected that the MISTiC Winds will be able to collect approximately 60,000 wind observations every 6 hours, if considering a constellation composed of 12 satellites (4 orbital planes). In addition, one of the main expected impacts of the MISTiC Winds concept is the ability to derive water vapor feature tracking AMVs below 500-400 hPa, an unique feature among the water vapor AMVs derived from the current Earth observing system.

Carvalho, David↗

NASA's GMAO Atmospheric Motion Vectors Simulator: Description and Application to the MISTiC Winds Concept

An atmospheric wind vectors (AMVs) simulator was developed by NASA's GMAO to simulate observations from future satellite constellation concepts. The synthetic AMVs can then be used in OSSEs to estimate and quantify the potential added value of new observations to the present Earth observing system and, ultimately, the expected impact on the current weather forecasting skill. The GMAO AMV simulator is a tunable and flexible computer code that is able to simulate AMVs expected to be derived from different instruments and satellite orbit configurations. As a case study and example of the usefulness of this tool, the GMAO AMV simulator was used to simulate AMVs envisioned to be provided by the MISTiC Winds, a NASA mission concept consisting of a constellation of satellites equipped with infrared spectral midwave spectrometers, expected to provide high spatial and temporal resolution temperature and humidity soundings of the troposphere that can be used to derive AMVs from the tracking of clouds and water vapor features. The GMAO AMV simulator identifies trackable clouds and water vapor features in the G5NR and employs a probabilistic function to draw a subset of the identified trackable features. Before the simulator is applied to the MISTiC Winds concept, the simulator was calibrated to yield realistic observations counts and spatial distributions and validated considering as a proxy instrument to the MISTiC Winds the Himawari-8 Advanced Imager (AHI). The simulated AHI AMVs showed a close match with the real AHI AMVs in terms of observation counts and spatial distributions, showing that the GMAO AMVs simulator synthesizes AMVs observations with enough quality and realism to produce a response from the DAS equivalent to the one produced with real observations. When applied to the MISTiC Winds scanning points, it can be expected that the MISTiC Winds will be able to collect approximately 60,000 wind observations every 6 hours, if considering a constellation composed of 12 satellites (4 orbital planes). In addition, one of the main expected impacts of the MISTiC Winds concept is the ability to derive water vapor feature tracking AMVs below 500-400 hPa, an unique feature among the water vapor AMVs derived from the current Earth observing system.

Carvalho, David↗

Planetary Rover Simulation for Lunar Exploration Missions

When planning planetary rover missions it is useful to develop intuition and skills driving in, quite literally, alien environments before incurring the cost of reaching said locales. Simulators make it possible to operate in environments that have the physical characteristics of target locations without the expense and overhead of extensive physical tests. To that end, NASA Ames and Open Robotics collaborated on a Lunar rover driving simulator based on the open source Gazebo simulation platform and leveraging ROS (Robotic Operating System) components. The simulator was integrated with research and mission software for rover driving, system monitoring, and science instrument simulation to constitute an end-to-end Lunar mission simulation capability. Although we expect our simulator to be applicable to arbitrary Lunar regions, we designed to a reference mission of prospecting in polar regions. The harsh lighting and low illumination angles at the Lunar poles combine with the unique reflectance properties of Lunar regolith to present a challenging visual environment for both human and computer perception. Our simulator placed an emphasis on high fidelity visual simulation in order to produce synthetic imagery suitable for evaluating human rover drivers with navigation tasks, as well as providing test data for computer vision software development.In this paper, we describe the software used to construct the simulated Lunar environment and the components of the driving simulation. Our synthetic terrain generation software artificially increases the resolution of Lunar digital elevation maps by fractal synthesis and inserts craters and rocks based on Lunar size-frequency distribution models. We describe the necessary enhancements to import large scale, high resolution terrains into Gazebo, as well as our approach to modeling the visual environment of the Lunar surface. An overview of the mission software system is provided, along with how ROS was used to emulate flight software components that had not been developed yet. Finally, we discuss the effect of using the high-fidelity synthetic Lunar images for visual odometry. We also characterize the wheel slip model, and find some inconsistencies in the produced wheel slip behaviour.

Allan, Mark↗

LongRunMIP - Motivation and Design for a Large Collection of Millennial-Length AO-GCM Simulations

LongRunMIP is the first collection of millennial-length simulations of complex coupled climate models and enables investigations of how these models equilibrate in response to radiative perturbations.We present a model intercomparison project, LongRunMIP, the first collection of millennial-length (1000+ year) simulations of complex coupled climate models with a representation of ocean, atmosphere, sea ice, and land surface, and their interactions. Standard model simulations are generally only a few hundred years long. However, modeling the long-term equilibration in response to radiative forcing perturbation is important for understanding many climate phenomena, such as the evolution of ocean circulation, time-and temperature-dependent feedbacks, and the differentiation of forced signal and internal variability. The aim of LongRunMIP is to facilitate research into these questions by serving as an archive for simulations that capture as much of this equilibration as possible. The only requirement to participate in LongRunMIP is to contribute a simulation with elevated, constant CO2 forcing that lasts at least 1000 years. LongRunMIP is a MIP of opportunity in that the simulations were mostly performed prior to the conception of the archive without an agreed-upon set of experiments. For most models, the archive contains a preindustrial control simulation and simulations with an idealized (typically abrupt) CO2 forcing. We collect 2D surface and top-of-atmosphere fields, and 3D ocean temperature and salinity fields. Here, we document the collection of simulations and discuss initial results, including the evolution of surface and deep ocean temperature and cloud radiative effects. As of summer 2019, the collection includes 50 simulations of 15 models by 10 modeling centers. The data of LongRunMIP are publicly available. We encourage submission of more simulations in the future.

Maria Rugenstein↗

Limitations in Representation of Physical Processes Prevents Successful Simulation of PM2.5 During KORUS-AQ

High levels of fine particulate matter (PM2.5) pollution in East Asia often exceed local air quality standards. Observations from the Korea United States-Air Quality (KORUS-AQ) field campaign in May and June 2016 showed that development of extreme pollution (haze) occurred through a combination of long-range transport and favorable meteorological conditions that enhanced local production of PM2.5. Atmospheric models often have difficulty simulating PM2.5 chemical composition during haze, which is of concern for the development of successful control measures. We use observations from KORUS-AQ to examine the ability of the GEOS-Chem chemical transport model to simulate PM2.5 composition throughout the campaign and identify the mechanisms driving the pollution event. At the surface, the model underestimates sulfate by -64% but overestimates nitrate by +36%. The largest underestimate in sulfate occurs during the pollution event, where models typically struggle to generate elevated sulfate concentrations due to missing heterogeneous chemistry in aerosol liquid water in the polluted boundary layer. Hourly surface observations show that the model nitrate bias is driven by an overestimation of the nighttime peak. In the model, nitrate formation is limited by the supply of nitric acid, which is biased by +100% against aircraft observations. We hypothesize that this is due to a large missing sink, which we implement here as a factor of five increase in dry deposition. We show that the resulting increased deposition velocity is consistent with observations of total nitrate as a function of photochemical age. The model does not account for factors such as the urban heat island effect or the heterogeneity of the built-up urban landscape resulting in insufficient model turbulence and surface area over the study area that likely results in insufficient dry deposition. Other species such as NH3 could be similarly affected but were not measured during the campaign. Nighttime production of nitrate is driven by NO2 hydrolysis in the model, while observations show that unexpectedly elevated nighttime ozone (not present in the model) should result in N2O5 hydrolysis as the primary pathway. The model is unable to represent nighttime ozone due to an overly rapid collapse of the afternoon mixed layer and excessive titration by NO. We attribute this to missing nighttime heating driving deeper nocturnal mixing that would be expected to occur in a city like Seoul. This urban heating is not considered in air quality models run at large enough scales to treat both local chemistry and long-range transport. Key model failures in simulating nitrate, mainly overestimated daytime nitric acid, incorrect representation of nighttime chemistry, and an overly shallow and insufficiently turbulent nighttime mixed layer, exacerbate the model’s inability to simulate the buildup of PM2.5 during haze pollution. To address the underestimate in sulfate most evident during the haze event, heterogeneous aerosol uptake of SO2 is added to the model which previously only considered aqueous production of sulfate from SO2 in cloud water. Implementing a simple parameterization of this chemistry improves the model abundance of sulfate but degrades the SO2 simulation implying that emissions are underestimated. We find that improving model simulations of sulfate has direct relevance to determining local vs. transboundary contributions to PM2.5. During the haze pollution event, the inclusion of heterogeneous aerosol uptake of SO2 decreases the fraction of PM2.5 attributable to long-range transport from 66% to 54%. Locally-produced sulfate increased from 1% to 25% of locally-produced PM2.5, implying that local emissions controls would have a larger effect than previously thought. However, this additional uptake of SO2 is coupled to the model nitrate prediction which affects the aerosol liquid water abundance and chemistry driving sulfate-nitrate-ammonium partitioning. An additional simulation of the haze pollution with heterogeneous uptake of SO2 to aerosol and simple improvements to the model nitrate simulation results in 30% less sulfate due to 40% less nitrate and aerosol water, and results in an underestimate of sulfate during the haze event. Future studies need to better consider the impact of model physical processes such as dry deposition and nighttime boundary layer mixing on the simulation of nitrate and the effect of improved nitrate simulations on the overall simulation of secondary inorganic aerosol (sulfate+nitrate+ammonium) in East Asia. Foreign emissions are rapidly changing, increasing the need to understand the impact of local emissions on PM2.5 in South Korea to ensure continued air quality improvements.

Katherine R. Travis↗

Wind-US Simulations of the NASA 1507 Inlet Test Case

Steady-state, Reynolds-Averaged Navier-Stokes (RANS) computational fluid dynamics (CFD) simulations were performed of the NASA 1507 Inlet Test Case from the 6th AIAA Propulsion Aerodynamics Workshop (PAW) held in January 2023 using the Wind-US CFD solver. The simulations were performed with bleed regions modeled as surface patches and the Slater bleed boundary condition models were applied within those regions. Both the constant-plenum-pressure and the fixed-plenum-exit bleed models where applied. The vortex generators were modeled using the Bender-Anderson-Yagle (BAY) vortex generator model. The CFD simulations were performed for inlet engine-face flow ratios varying from supercritical to subcritical inlet operation until inlet unstart. The results of the CFD simulations were compared to wind-tunnel data which included boundary layer profiles, axial static pressure distributions, bleed rates, total pressure recovery performance curves, and engine-face total pressure rake profiles. The CFD simulations were able to match most of the boundary layer rake data near the geometric throat, but the boundary layers on the centerbody were thinner than indicated by the wind-tunnel data. For the supercritical CFD simulations, the engine-face flow ratio was at least 2% to 2.5% lower than indicated by the wind-tunnel data, but still consistent with the level of supercritical bleed. The engine-face flow rates of the wind-tunnel tests had an uncertainty of ±2% which could explain the differences. The supercritical engine-face flow rates of the wind-tunnel data indicated a 5% supercritical bleed rate rather than the reported 7.18% bleed rate. The CFD simulations were able to calculate the maximum total pressure recovery of pt2/pt = 0.907 as observed for the wind-tunnel data. The CFD simulations with the constant-pressure bleed model were able to obtain a much lower engine-face total pressure ratio than indicated by both the wind-tunnel data and the CFD simulations using a fixed-plenum-exit bleed model. The CFD simulations indicated a sharp static pressure rise for the terminal shock rather than a more gradual rise indicated by the wind-tunnel data. The CFD simulations were not able to accurately compute the total pressure variation across the engine-face total pressure rake.

Supersonics↗

Simulating weather regimes: impact of model resolution and stochastic parameterization

The simulation of quasi-persistent regime structures in an atmospheric model with horizontal resolution typical of the Intergovernmental Panel on Climate Change fifth assessment report simulations, is shown to be unrealistic. A higher resolution configuration of the same model, with horizontal resolution typical of that used in operational numerical weather prediction, is able to simulate these regime structures realistically. The spatial patterns of the simulated regimes are remarkably accurate at high resolution. A model configuration at intermediate resolution shows a marked improvement over the low-resolution configuration, particularly in terms of the temporal characteristics of the regimes, but does not produce a simulation as accurate as the very-high-resolution configuration. It is demonstrated that the simulation of regimes can be significantly improved, even at low resolution, by the introduction of a stochastic physics scheme. At low resolution the stochastic physics scheme drastically improves both the spatial and temporal aspects of the regimes simulation. These results highlight the importance of small-scale processes on large-scale climate variability, and indicate that although simulating variability at small scales is a necessity, it may not be necessary to represent the small-scales accurately, or even explicitly, in order to improve the simulation of large-scale climate. It is argued that these results could have important implications for improving both global climate simulations, and the ability of high-resolution limited-area models, forced by low-resolution global models, to reliably simulate regional climate change signals.

Meteorology & Atmospheric Sciences↗

Emulation of dynamic simulators with application to hydrology

Many simulation-intensive tasks in the applied sciences, such as sensitivity analysis, parameter inference or real time control, are hampered by slow simulators. Emulators provide the opportunity of speeding up simulations at the cost of introducing some inaccuracy. An emulator is a fast approximation to a simulator that interpolates between design input–output pairs of the simulator. Increasing the number of design data sets is a computationally demanding way of improving the accuracy of emulation. We investigate the complementary approach of increasing emulation accuracy by including knowledge about the mechanisms of the simulator into the formulation of the emulator. To approximately reproduce the output of dynamic simulators, we consider emulators that are based on a system of linear, ordinary or partial stochastic differential equations with a noise term formulated as a Gaussian process of the parameters to be emulated. This stochastic model is then conditioned to the design data so that it mimics the behavior of the nonlinear simulator as a function of the parameters. The drift terms of the linear model are designed to provide a simplified description of the simulator as a function of its key parameters so that the required corrections by the conditioned Gaussian process noise are as small as possible. The goal of this paper is to compare the gain in accuracy of these emulators by enlarging the design data set and by varying the degree of simplification of the linear model. We apply this framework to a simulator for the shallow water equations in a channel and compare emulation accuracy for emulators based on different spatial discretization levels of the channel and for a standard non-mechanistic emulator. Our results indicate that we have a large gain in accuracy already when using the simplest mechanistic description by a single linear reservoir to formulate the drift term of the linear model. Adding some more reservoirs does not lead to a significant improvement in accuracy. However, the transition to a spatially continuous linear model leads again to a similarly large gain in accuracy as the transition from the non-mechanistic emulator to that based on one reservoir.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Evaluation of the Impact of Horizontal Grid Spacing in Terra Incognita on Coupled Mesoscale–Microscale Simulations Using the WRF Framework

Coupled mesoscale-microscale simulations are required to provide time-varying weather-dependent inflow and forcing for large-eddy simulations under general flow conditions. Such coupling necessarily spans a wide range of spatial scales (i.e., ~10 m to ~10 km). In this paper, we use simulations that involve multiple nested domains with horizontal grid spacings in the terra incognita (i.e., ≲km) that may affect simulated conditions in both the outer and inner domains. We examine the impact on simulated wind speed and turbulence associated with forcing provided by a terrain with grid spacing in the terra incognita. We perform a suite of simulations that use combinations of varying horizontal grid spacings and turbulence parameterization/modeling using the Weather Research and Forecasting (WRF) Model using a combination of planetary boundary layer (PBL) and large-eddy simulation subgrid-scale (LES-SGS) models. The results are analyzed in terms of spectral energy, turbulence kinetic energy, and proper orthogonal decomposition (POD) energy. The results show that the output from the microscale domain depends on the type of turbulence model (e.g., PBL or LES-SGS model) used for a given horizontal grid spacing but is independent of the horizontal grid spacing and turbulence modeling of the parent domain. Simulation using a single domain produced less POD energy in the first few modes compared to a coupled simulation (one-way nesting) for similar horizontal grid spacing, which highlights that coupled simulations are required to accurately pass the mesoscale features into the microscale domain.

97 MATHEMATICS AND COMPUTING↗

Whistler turbulence heating of electrons and ions: Three-dimensional particle-in-cell simuations

In this study, the decay of whistler turbulence in a collisionless, homogeneous, magnetized plasma is studied using three-dimensional particle-in-cell simulations. The simulations are initialized with a narrowband, relatively isotropic distribution of long wavelength whistler modes. A first ensemble of simulations at electron beta ${\beta }_{{\rm{e}}}$ = 0.25 and ion-to-electron mass ratio ${m}_{{\rm{i}}}$/${m}_{{\rm{e}}}$ = 400 is carried out on a domain cube of dimension $L{\omega }_{\mathrm{pi}}$/c = 5.12 where ${\omega }_{\mathrm{pi}}$ is the ion plasma frequency. The simulations begin with a range of dimensionless fluctuating field energy densities, ${\epsilon }_{{\rm{o}}}$, and follow the fluctuations as they cascade to broadband, anisotropic turbulence which dissipates at shorter wavelengths, heating both electrons and ions. The electron heating is stronger and preferentially parallel/antiparallel to the background magnetic field ${{\boldsymbol{B}}}_{{\rm{o}}};$ the ion energy gain is weaker and is preferentially in directions perpendicular to ${{\boldsymbol{B}}}_{{\rm{o}}}$. The important new results here are that, over 0.01 < ${\epsilon }_{{\rm{o}}}$ < 0.25, the maximum rate of electron heating scales approximately as ${\epsilon }_{{\rm{o}}}$, and the maximum rate of ion heating scales approximately as ${\epsilon }_{{\rm{o}}}^{1.5}$. A second ensemble of simulations at ${\epsilon }_{{\rm{o}}}$ = 0.10 and ${\beta }_{{\rm{e}}}$ = 0.25 shows that, over 25 < ${m}_{{\rm{i}}}$/${m}_{{\rm{e}}}\;$< 1836, the ratio of the maximum ion heating rate to the maximum electron heating rate scales approximately as ${m}_{{\rm{e}}}$/${m}_{{\rm{i}}}$.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

SpaceFOM: An Interoperability Standard for Space Systems Simulations

There is a long history of simulation supporting space systems development. This includes relatively simple parametric simulations to more complex trajectory simulations to large scale integrated vehicle simulation. One area of relatively recent development is in the area of distributed or interoperable simulation. Distributed simulation has been in wide use by the US military for years but is being used more widely in the aerospace community. To support large scale distributed simulation, the military community has developed a number of standards to support a priori interoperability between large collections of disparate simulation. For example, the IEEE 1516 High Level Architecture (HLA) and the Real-time Platform Reference Federation Object Model (RPR FOM). While HLA is suitable for space systems, there are a number of design decisions made in the development of the RPR FOM that prevent it from working well for space applications. In order to address these deficiencies, the Simulation Interoperability Standards Organization (SISO) developed a new HLS-compatible interoperability standard to support the needs of complex space systems. This standard is the Space Reference Federation Object Model (SpaceFOM). This paper presents on overview of the SpaceFOM including the fundamentals of the SpaceFOM, the key features of the SpaceFOM, and how the SpaceFOM supports large scale distributed simulation of complex space systems.

SISO↗

Suitability of Synthetic Driving Profiles from Traffic Micro-Simulation for Real-World Energy Analysis: Preprint

A shift towards increased levels of driving automation is generally expected to result in improved safety and traffic congestion outcomes. However, little empirical data exists to estimate the impact that automated driving could have on energy consumption and greenhouse gas emissions. In the absence of empirical data on differences between drive cycles from present day vehicles (primarily operated by humans) and future vehicles (partially or fully operated by computers) one approach is to model both situations over identical traffic conditions. Such an exercise requires traffic micro-simulation to not only accurately model vehicle operation under high levels of automation, but also (and potentially more challenging) vehicle operation under present day human drivers. This work seeks to quantify the ability of a commercial traffic micro-simulation program to accurately model real-world drive cycles in vehicles operated primarily by humans in terms of driving speed, acceleration, and simulated fuel economy. Synthetic profiles from models of freeway and arterial facilities near Atlanta, Georgia, are compared to empirical data collected from real-world drivers on the same facilities. Empirical and synthetic drive cycles are then simulated in a powertrain efficiency model to enable comparison on the basis of fuel economy. Synthetic profiles from traffic micro-simulation were found to exhibit low levels of transient behavior relative to the empirical data. Even with these differences, the synthetic and empirical data in this study agree well in terms of driving speed and simulated fuel economy. The differences in transient behavior between simulated and empirical data suggest that larger stochastic contributions in traffic micro-simulation (relative to those present in the traffic micro-simulation tool used in this study) are required to fully capture the arbitrary elements of human driving. Interestingly, the lack of stochastic contributions from models of human drivers in this study did not result in a significant discrepancy between fuel economy simulations based on synthetic and empirical data; a finding with implications on the potential energy efficiency gains of automated vehicle technology.

33 ADVANCED PROPULSION SYSTEMS↗

Engineering evaluation of the EQSIM simulated ground-motion database: The San Francisco Bay Area region

Ground-motion simulations for infrequent earthquake scenarios are gaining increasing interest in the engineering community for their potential to inform performance-based structural design and assessment, particularly in regions where major earthquakes are expected, but recordings from consistent historical events are not available. However, the absence of empirical data makes the evaluation of such scenarios and the interpretation of the simulation results very challenging. In this context, this paper presents the evaluation of the first EQSIM (v.1.1.0) ground-motion simulations generated for the San Francisco Bay Area (SFBA) region. The current database, which is at the first stages of development, includes eight realizations of a Hayward Fault M w 7 event. The evaluation of the simulated motions is first conducted on the average ground-motion amplitudes through comparison against the NGA-W2 GMPEs and a selected population of real records. A methodology for selecting a population of real records consistent with the simulated event is proposed. The objective of the proposed evaluation is twofold: (1) build confidence in the realistic character of the simulated motions for their use in engineering applications and (2) offer a critical physics-based interpretation of the simulation results that can help improve key features of the simulation models. To further investigate the implications of using the simulated ground motions for site-specific structural assessments, the ASCE7-16 approach is employed for the analysis of two hazard-consistent near-field sites in the SFBA. Results are discussed. Furthermore, this study offers a critical review of the EQSIM (v.1.1.0) SFBA ground-motion simulations and suggestions for improving the earthquake rupture models.

58 GEOSCIENCES↗

Simulation and visualization of energy-related occupant behavior in office buildings

In current building performance simulation programs, occupant presence and interactions with building systems are over-simplified and less indicative of real world scenarios, contributing to the discrepancies between simulated and actual energy use in buildings. Simulation results are normally presented using various types of charts. However, using those charts, it is difficult to visualize and communicate the importance of occupants’ behavior to building energy performance. This study introduced a new approach to simulating and visualizing energy-related occupant behavior in office buildings. First, the Occupancy Simulator was used to simulate the occupant presence and movement and generate occupant schedules for each space as well as for each occupant. Then an occupant behavior functional mockup unit (obFMU) was used to model occupant behavior and analyze their impact on building energy use through co-simulation with EnergyPlus. Finally, an agent-based model built upon AnyLogic was applied to visualize the simulation results of the occupant movement and interactions with building systems, as well as the related energy performance. A case study using a small office building in Miami, FL was presented to demonstrate the process and application of the Occupancy Simulator, the obFMU and EnergyPlus, and the AnyLogic module in simulation and visualization of energy-related occupant behaviors in office buildings. Furthermore, the presented approach provides a new detailed and visual way for policy makers, architects, engineers and building operators to better understand occupant energy behavior and their impact on energy use in buildings, which can improve the design and operation of low energy buildings.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

A Fast Algorithm for Scanning Transmission Electron Microscopy Imaging and 4D-STEM Diffraction Simulations

Scanning transmission electron microscopy (STEM) is an extremely versatile method for studying materials on the atomic scale. Many STEM experiments are supported or validated with electron scattering simulations. However, using the conventional multislice algorithm to perform these simulations can require extremely large calculation times, particularly for experiments with millions of probe positions as each probe position must be simulated independently. Recently, the plane-wave reciprocal-space interpolated scattering matrix (PRISM) algorithm was developed to reduce calculation times for large STEM simulations. Here, we introduce a new method for STEM simulation: partitioning of the STEM probe into “beamlets,” given by a natural neighbor interpolation of the parent beams. This idea is compatible with PRISM simulations and can lead to even larger improvements in simulation time, as well requiring significantly less computer random access memory (RAM). We have performed various simulations to demonstrate the advantages and disadvantages of partitioned PRISM STEM simulations. We find that this new algorithm is particularly useful for 4D-STEM simulations of large fields of view. We also provide a reference implementation of the multislice, PRISM, and partitioned PRISM algorithms.

97 MATHEMATICS AND COMPUTING↗

Gyrokinetic GENE simulations of DIII-D near-edge L-mode plasmas

We present gyrokinetic simulations with the GENE code addressing the near-edge region of an L-mode plasma in the DIII-D tokamak. At radial position ρ = 0.80, simulations with the ion temperature gradient (ITG) increased by 40% above the nominal value give electron and ion heat fluxes that are in simultaneous agreement with the experiment. This gradient increase is consistent with the combined statistical and systematic uncertainty σ of the charge exchange recombination spectroscopy measurements at the 1.6σ level. Multiscale simulations are carried out with a realistic mass ratio and geometry for the first time in the near-edge. These multiscale simulations suggest that the highly unstable ion temperature gradient (ITG) modes of the flux-matched ion-scale simulations suppress electron-scale transport, such that ion-scale simulations are sufficient at this location. At radial position ρ = 0.90, nonlinear simulations show a hybrid state of ITG and trapped electron modes, which was not expected from linear simulations. The nonlinear simulations reproduce the total experimental heat flux with the inclusion of E × B shear effects and an increase in the electron temperature gradient by ∼23%. This gradient increase is compatible with the combined statistical and systematic uncertainty of the Thomson scattering data at the 1.3σ level. These results are consistent with previous findings that gyrokinetic simulations are able to reproduce the experimental heat fluxes by varying input parameters close to their experimental uncertainties, pushing the validation frontier closer to the edge region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

COMPARISONS OF COSMOLOGICAL MAGNETOHYDRODYNAMIC GALAXY CLUSTER SIMULATIONS TO RADIO OBSERVATIONS

Radio observations of galaxy clusters show that there are {mu}G magnetic fields permeating the intracluster medium (ICM), but it is hard to accurately constrain the strength and structure of the magnetic fields without the help of advanced computer simulations. We present qualitative comparisons of synthetic Very Large Array observations of simulated galaxy clusters to radio observations of Faraday rotation measure (RM) and radio halos. The cluster formation is modeled using adaptive mesh refinement magnetohydrodynamic simulations with the assumption that the initial magnetic fields are injected into the ICM by active galactic nuclei (AGNs) at high redshift. In addition to simulated clusters in Xu et al., we present a new simulation with magnetic field injections from multiple AGNs. We find that the cluster with multiple injection sources is magnetized to a similar level as in previous simulations with a single AGN. The RM profiles from simulated clusters, both |RM| and the dispersion of RM ({sigma}{sub RM}), are consistent at a first order with the radial distribution from observations. The correlations between the {sigma}{sub RM} and X-ray surface brightness from simulations are in a broad agreement with the observations, although there is an indication that the simulated clusters could be slightly overdense and less magnetized with respect to those in the observed sample. In addition, the simulated radio halos agree with the observed correlations between the radio power versus the cluster X-ray luminosity and between the radio power versus the radio halo size. These studies show that the cluster-wide magnetic fields that originate from AGNs and are then amplified by the ICM turbulence match observations of magnetic fields in galaxy clusters.

79 ASTRONOMY AND ASTROPHYSICS↗