Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “velocity steps”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

Reduction of Gun Erosion and Correlation of Gun Erosion Measurements

Gun barrel erosion is serious problem with two-stage light gas guns. Excessive barrel erosion can lead to poor or failed launches and frequent barrel changes, with the corresponding down time. Also, excessive barrel erosion can limit the maximum velocity obtainable by loading down the hydrogen working gas with eroded barrel material. Guided by a CFD code, the operating conditions of the Ames 0.5-inch gun were modified to reduce barrel erosion. The changes implemented included: (1) reduction in the piston mass, powder mass and hydrogen fill pressure; and (2) reduction in pump tube volume, while maintaining hydrogen mass. The latter change was found, in particular, to greatly reduce barrel erosion. For muzzle velocity ranges of 6.1 - 6.9 km/sec, the barrel erosion was reduced by a factor of 10. Even for the higher muzzle velocity range of 7.0 - 8.2 km/sec, the barrel erosion was reduced by a factor of 4. Gun erosion data from the Ames 0.5-inch, 1.0-inch, and 1.5-inch guns operated over a wide variety of launch conditions was examined and it was found that this data could be correlated using four different parameters: normalized powder charge energy, normalized hydrogen energy density, normalized pump tube volume and barrel diameter. The development of the correlation and the steps used to collapse the experimental data are presented. Over a certain parameter range in the correlation developed, the barrel erosion per shot is found to increase very rapidly. The correlation should prove useful in the selection of gun operating conditions and the design of new guns. Representative shapes of eroded gun barrels are also presented.

Bogdanoff, Dave↗
NASA NTRS2004

Microgravity Spray Cooling Research for High Powered Laser Applications

An extremely powerful laser is being developed at Goddard Space Flight Center for use on a satellite. This laser has several potential applications. One application is to use it for upper atmosphere weather research. In this case, the laser would reflect off aerosols in the upper atmosphere and bounce back to the satellite, where the aerosol velocities could be calculated and thus the upper atmosphere weather patterns could be monitored. A second application would be for the US. Air Force, which wants to use the laser strategically as a weapon for satellite defense. The Air Force fears that in the coming years as more and more nations gain limited space capabilities that American satellites may become targets, and the laser could protect the satellites. Regardless of the ultimate application, however, a critical step along the way to putting the laser in space is finding a way to efficiently cool it. While operating the laser becomes very hot and must be cooled to prevent overheating. On earth, this is accomplished by simply running cool tap water over the laser to keep it cool. But on a satellite, this is too inefficient. This would require too much water mass to be practical. Instead, we are investigating spray cooling as a means to cool the laser in microgravity. Spray cooling requires much less volume of fluid, and thus could be suitable for use on a satellite. We have inherited a 2.2 second Drop Tower rig to conduct our research with. In our experiments, water is pressurized with a compressed air tank and sprayed through a nozzle onto our test plate. We can vary the pressure applied to the water and the temperature of the plate before an experiment trial. The whole process takes place in simulated microgravity in the 2.2 second Drop Tower, and a high speed video camera records the spray as it hits the plate. We have made much progress in the past few weeks on these experiments. The rig originally did not have the capability to heat the test plate, but I did some heat transfer calculations and picked out a heater to order for the rig. I learned QBasic programming language to change the operating code for our drops, allowing us to rapidly cycle the spray nozzle open and closed to study the effects. We have derived an equation for flow rate vs. pressure for our experiment. We have recorded several videos of drops at different pressures, some with heated test plate and some without, and have noticed substantial differences in the liquid behavior. I have also changed the computer program to write a file with temperature vs. time profiles for the test plate, and once the necessary thermocouple comes in (it was ordered last week), we will have temperature profiles to accompany the videos. Once we have these temperature profiles to go with the videos, we will be able to see how the temperature is affected by the spray at different pressures, and how the spray changes its behavior once as the plate changes from hot to cool. With quantitative temperature data, we can then mathematically model the heat transfer from the plate to the cooling spray. Finally, we can look at the differences between trials in microgravity and those in normal earth gravity.

Zivich, Chad P.↗
DOE OSTI2026

Kinetic model development for single step ethanol to butene rich olefin process over Cu-Y/Beta catalysts

Here, this study presents the first intrinsic kinetic model for the single-step conversion of ethanol to butene-rich olefins over bifunctional Cu-Y/Beta catalysts, addressing a critical gap in the design and scale-up of Sustainable Aviation Fuel (SAF) processes. The reaction network comprises ten global steps involving dehydrogenation, aldol condensation, hydrogenation, and dehydration reactions, distributed across Cu and rare-earth (Y) active sites. The model incorporates dual-site functionality (Cu and Y site) and explicitly accounts for key intermediates such as crotonaldehyde and butanal. Reaction rates are formulated using Langmuir–Hinshelwood–Hougen–Watson (LHHW) kinetics. Kinetic parameters are extracted by fitting the model to lab-scale packed-bed reactor data across a wide range of temperatures and space velocities, demonstrating strong agreement in ethanol conversion and product selectivity. The reaction kinetics developed in this work provide a foundational basis for constructing reactor models that enable process optimization and scale-up of ethanol-to-jet fuel technologies.

Cu-Y/Beta catalyst↗
NASA NTRS1985

Numerical simulation of a turbulent flame stabilized behind a rearward-facing step

Flow of combustible mixtures in a plane channel past a smooth contraction followed by an abrupt expansion, in a typical dump combustor configuration, is modeled by a two-dimensional numerical technique based on the random vortex method. Both the inert and the reacting case are considered. In the latter, the flame is treated as an interface, self-advancing at a prescribed normal burning speed, while the dynamic effects of expansion due to the exothermicity of combustion are expressed by volumetric source lines delineated by its front. Solutions are shown to be in satisfactory agreement with experimental results, especially with respect to global properties such as the average velocity profiles and the reattachment length. The stochastic turbulent velocity components manifest interesting differences, especially near the walls where three-dimensional effects of turbulence are expected to be of importance.

Hsiao, C. C.↗
DOE OSTI2024

Towards an understanding of particle-scale flaws and microstructure evolution in cold spray via accumulation of single particle impacts

Cold spray coatings are the sum of countless individual bonding events between single particles impacting on top of one another at high velocities. Thus, the collective behavior of microparticles must be considered to elucidate the origins of coating flaws at the scale of the particles and larger, or the dynamic evolution of the overall coating microstructure. Laser-induced particle impact testing (LIPIT) has been extensively used to study single-particle impacts, and in this work is adapted to study the accumulation of numerous particles with knowledge of each individual particle's impact parameters (particle size, velocity). The method reproducibly deposits stacks of gold particles (>20 particles) with different characteristic spectra of impact velocity. The observation of impact-induced erosion lets us define a critical velocity for material-build-up that is higher than that for single-particle bonding. The quantitative single-particle data are analyzed in a correlative manner to the structure and flaws in the resulting stacks, providing some first statistical connections between, e.g., strain and recrystallization, or aberrant particle characteristics and defects. The results highlight opportunities for the study of many-particle phenomena in microparticle impact—from interaction of particles in cold spray to multi-step erosion processes—with a quantitative view of the behavior of single particles.

Erosion↗
NASA NTRS1981

Development of an engineering analysis of progressive damage in composites during low velocity impact

A computerized, analytical methodology was developed to study damage accumulation during low velocity lateral impact of layered composite plates. The impact event was modeled as perfectly plastic with complete momentum transfer to the plate structure. A transient dynamic finite element approach was selected to predict the displacement time response of the plate structure. Composite ply and interlaminar stresses were computed at selected time intervals and subsequently evaluated to predict layer and interlaminar damage. The effects of damage on elemental stiffness were then incorporated back into the analysis for subsequent time steps. Damage predicted included fiber failure, matrix ply failure and interlaminar delamination.

Humphreys, E. A.↗
NASA NTRS2002

Understanding the Solar Sources of In Situ Observations

The solar wind can, to a good approximation be described as a two component flow with fast, tenuous, quiescent flow emanating from coronal holes, and slow, dense and variable flow associated with the boundary between open and closed magnetic fields. In spite of its simplicity, this picture naturally produces a range of complex heliospheric phenomena, including the presence, location, and orientation of corotating interaction regions and their associated shocks. In this study, we apply a two-step mapping technique, incorporating a magnetohydrodynamic model of the solar corona, to bring in situ observations h m Ulysses, WIND, and ACE back to the solar surface in an effort to determine some intrinsic properties of the quasi-steady solar wind. In particular, we find that a "layer" of approx. 35,000 h n exists between the Coronal Hole Boundary (CHB) and the fast solar wind, where the wind is slow and variable. We also- derive a velocity gradient within large polar coronal boles (that were present during Ulysses rapid latitude scan) as a function of distance from the CHB. We find that nu = 713 km/s + 3.2 d, where d is the angular distance from the CHB boundary in degrees.

Riley, Pete↗
DOE OSTI2023

Convolutional Neural Networks for Problems in Transport Phenomena: A Theoretical Minimum

Convolutional neural network (CNN), a deep learning algorithm, has gained popularity in technological applications that rely on interpreting images (typically, an image is a 2D field of pixels). Transport phenomena is the science of studying different fields representing mass, momentum, or heat transfer. Some of the common fields are species concentration, fluid velocity, pressure, and temperature. Each of these fields can be expressed as an image(s). Consequently, CNNs can be leveraged to solve specific scientific problems in transport phenomena. Herein, we show that such problems can be grouped into three basic categories: (a) mapping a field to a descriptor (b) mapping a field to another field, and (c) mapping a descriptor to a field. After reviewing the representative transport phenomena literature for each of these categories, we illustrate the necessary steps for constructing appropriate CNN solutions using sessile liquid drops as an exemplar problem. If sufficient training data is available, CNNs can considerably speed up the solution of the corresponding problems. Finally, the present discussion is meant to be minimalistic such that readers can easily identify the transport phenomena problems where CNNs can be useful as well as construct and/or assess such solutions.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1992

Investigation of forced unsteady separated flows using velocity-vorticity form of Navier-Stokes equations

The phenomenon of forced unsteady separation and eruption of boundary-layer vorticity is a highly-complex, high-Reynolds number flow phenomenon, which abruptly leads to the formation of a dynamic stall vortex as demonstrated earlier by the authors for a NACA 0015 airfoil undergoing constant rate pitch-up motion. This, as well as the results of other researchers, have convincingly demonstrated a complex vortical structure within the state of unsteady separation prior to the evolution of dynamic stall. This phenomenon of vortex eruption, although observed in studying dynamic stall phenomena, is also associated with transition from laminar to turbulence flow and its generic nature has been stressed by many researchers including the present investigators. An unsteady Navier-Stokes (NS) analysis is developed for arbitrarily maneuvering bodies using velocity-vorticity variables; this formulation is nearly form-invariant under a generalized non-inertial coordinate transformation. A fully-implicit uniformly second-order accurate method is used, with the nonlinear convective terms approximated using a biased third-order upwind differencing scheme to be able to simulate higher-Re flows. No explicit artificial dissipation is added. The numerical method is fully vectorized and currently achieves a computational index of 7 micro-seconds per time step per mesh point, using a single processor on a CRAY Y-MP. The simulation results show that the energetic free shear from the leading edge is responsible for the wall viscous layer to abruptly erupt near the center of the counterclockwise rotating eddy in the unsteady boundary layer. Primary, secondary, tertiary and quaternary vortices have been observed before the dynamic stall vortex evolves and gathers its maximum strength. This study will discuss the simulation results of Reynolds number up to Re = 45,000 and will also discuss the efforts of initial acceleration in a specific maneuver, on the evolution of the stall vortex.

Ghia, K. N.↗
DOE OSTI2022

Polyethylene upcycling to fuels: Narrowing the carbon number distribution in n-alkanes by tandem hydropyrolysis/hydrocracking

The extensive use of plastics in modern life has resulted in a global waste crisis to the environment. Polyethylene (PE) is one of the most popular and hardest plastics to recycle because of its strong C(sp 3 )-C(sp 3 ) bonds. In this study, a tandem conversion process, i.e., hydropyrolysis and subsequent vapor-phase hydrocracking of primary intermediates (C n>5 alkenes and long-chain alkanes), was conducted via a two-step pressurized flow-through fixed-bed reactor over a CoAl 2 O 4 spinel-derived catalyst. The product distribution could be flexibly tuned by regulating operating parameters in the cascade fixed-bed reactor and Co/Al molar ratio in CoAl 2 O 4 spinel catalysts. Under optimal reaction conditions (0.2 MPa H 2 , 550 °C for hydropyrolysis in the 1st reactor, 300–325 °C for hydrocracking in the 2nd reactor, 20 s -1 of gas hourly space velocity), the maximum single-pass yields of gasoline (C 5 -C 12 ) and jet-fuel (C 8 -C 16 ) range n-alkanes reached 86.0 wt% and 68.1 wt%, respectively. The CoAl 2 O 4 spinel catalysts also gained high activity in degrading realistic post-consumer plastics such as linear low-, low-, and high-density PE, and a ~ 73.1 wt% gasoline yield and a ~ 54.7 wt% C 8 -C 16 yield were retained even after 3 cycles in-situ regeneration of deactivated CoAl 2 O 4 catalysts. This work provides an efficient and tunable approach to upcycle PE wastes into liquid fuels with an ideal carbon length.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Development of New Experimental Methods for Correlated Operando Surface/Gas Characterization

The predictive understanding of catalytic surface reactions requires accurate microkinetic models, and while decades of work has been devoted to the elucidation of the reaction steps in these models, many open questions remain. One key issue is a lack of approaches enabling the local spatially resolved assessment of catalytic activity over a surface. In this report, we detail efforts to develop a new diagnostic approach to solve this problem. The approach is based upon laser resonance enhanced multiphoton ionization of reaction products emitted into the gas phase followed by spatially resolved imaging of the resultant ions or electrons. Ion imaging is pursued with a velocity-selected spatially resolved ion imaging microscope, while electron imaging was attempted in a low energy electron microscope. Successful demonstration of the ion imaging microscope coupled with the development of transport simulations shows promise for a revolutionary new tool to assess local catalytic activity

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2010

Validation of Supersonic Film Cooling Modeling for Liquid Rocket Engine Applications

Topics include: upper stage engine key requirements and design drivers; Calspan "stage 1" results, He slot injection into hypersonic flow (air); test articles for shock generator diagram, slot injector details, and instrumentation positions; test conditions; modeling approach; 2-d grid used for film cooling simulations of test article; heat flux profiles from 2-d flat plate simulations (run #4); heat flux profiles from 2-d backward facing step simulations (run #43); isometric sketch of single coolant nozzle, and x-z grid of half-nozzle domain; comparison of 2-d and 3-d simulations of coolant nozzles (run #45); flowfield properties along coolant nozzle centerline (run #45); comparison of 3-d CFD nozzle flow calculations with experimental data; nozzle exit plane reduced to linear profile for use in 2-d film-cooling simulations (run #45); synthetic Schlieren image of coolant injection region (run #45); axial velocity profiles from 2-d film-cooling simulation (run #45); coolant mass fraction profiles from 2-d film-cooling simulation (run #45); heat flux profiles from 2-d film cooling simulations (run #45); heat flux profiles from 2-d film cooling simulations (runs #47, #45, and #47); 3-d grid used for film cooling simulations of test article; heat flux contours from 3-d film-cooling simulation (run #45); and heat flux profiles from 3-d and 2-d film cooling simulations (runs #44, #46, and #47).

Morris, Christopher I.↗
NASA NTRS2018

A Formally Verified Floating-Point Implementation of the Compact Position Reporting Algorithm

The Automatic Dependent Surveillance-Broadcast (ADS-B) system allows aircraft to communicate their current state, including position and velocity information, to other aircraft in their vicinity and to ground stations. The Compact Position Reporting (CPR) algorithm is the ADS-B module responsible for the encoding and decoding of aircraft positions. CPR is highly sensitive to computer arithmetic since it heavily relies on functions that are intrinsically unstable such as floor and modulo. In this paper, a formally-verified double-precision floating-point implementation of the CPR algorithm is presented. The verification proceeds in three steps. First, an alternative version of CPR, which reduces the floating-point rounding error is proposed. Then, the Prototype Verification System (PVS) is used to formally prove that the ideal real-number counterpart of the improved algorithm is mathematically equivalent to the standard CPR definition. Finally, the static analyzer Frama-C is used to verify that the double-precision implementation of the improved algorithm is correct with respect to its operational requirement. The alternative algorithm is currently being considered for inclusion in the revised version of the ADS-B standards document as the reference implementation of the CPR algorithm.

Laura Titolo↗
DOE OSTI2024

HERO WEC V1.0 2024 - WEC-Sim Detailed Simulation Runs and Summary Data

This dataset includes results from simulations of NREL's hydraulic and electric reverse osmosis wave energy converter (HEREO WEC). Simulation runs include 135 wave cases that were based on the updated WEC-Sim model, which is linked below. The data represented in this repository is based on an updated WEC-Sim model using laboratory data to tune and refine the original WEC-Sim model for the V1.0 HERO WEC. The 135 wave cases represent waves with the following wave height and wave period ranges: - Significant Wave Height: 0.25 - 3.75m in 0.25m increments - Wave Period: 5 - 13 sec in 1 sec increments Each run was simulated using a Pierson-Moskowitz irregular wave spectrum with a 100 second ramp time, a total simulation time of 3,100 seconds, and a simulation time-step of 0.005s. A reference table has been included to map each multi condition run (MCR) case with each wave condition. Summary data set includes a spreadsheet and image files with matrices that are associated with data from simulation runs. All matrices cover the same significant wave height and wave periods from the simulation runs, in the same increments. The following matrices are included: - Power Abs: The average absorbed power from the WEC (calculated from anchor reaction force and heave velocity) - Power Hyd: The average hydraulic power output at pump (calculated from pump output flow and pressure) - Power - Hyd ROi: The average hydraulic power measured at the RO system inlet (calculated from RO system pressure and flow (pre-accumulator)) - Flow - Pump out: The average flowrate measured at the pump outlet - Flow - Perm: The average permeate (clean water) production - Flow - RO (pre): The average flowrate measured at the inlet of the RO system before the accumulators - Flow - RO (post): The average flowrate measured after the accumulator bank in the RO system - Pressure - RO: The average pressure measured at the inlet of the RO system This data set has been developed by the National Renewable Energy Laboratory, operated by the Alliance for Sustainable Energy, LLC, for the U.S. Department of Energy (DOE) under Contract No. DE-AC36-08GO28308. Funding provided by the U.S. Department of Energy Office of Energy Efficiency and Renewable Energy Water Power Technologies Office.

16 TIDAL AND WAVE POWER↗
DOE OSTI2022

Impact of particle size and particle-flow-wall coupling on pressurized oxy-combustion in the down-fired burner

Concerns over climate change have led to numerous efforts in developing low-carbon energy technologies. Pressurized oxy-combustion (POC) is a promising candidate to reduce carbon emissions in power generation. Designing an effective burner plays a vital role in developing new coal combustion technologies. Because of the high pressure, the volume fraction of coal particles at the fuel inlet of the burner of a pressurized oxy-combustor is close or even higher than the maximum limit that commercial CFD codes (e.g., ANSYS FLUENT) can handle. At this high particle volume fraction, the interactions among particles, fluid flow, and wall need to be re-evaluated. The present computational work is the first step of a systematic analysis of the particle influence, like releasing method, releasing locating, and particle size, in a pilot-scale POC combustor, developed at Washington University in St. Louis (WUSTL). In a POC process, pulverized coal is burned under elevated pressure and O2-CO2 environment. Specifically, a 15-bar, 100 kWth, POC combustor is modeled employing ANSYS FLUENT, using Reynolds-averaged Navier-Stokes (RANS) modeling. It is revealed that for this pilot-scale, pressurized burner, velocity profiles of the near-wall region in this POC facility exhibit some discrepancy against the near-wall turbulent flow velocity profile. In order to investigate the particle influence in the near-wall region, particle releasing location will be investigated. Numerical simulation results exhibit the coupling effect of turbulence flow and particles in this case. The particle size also demonstrates a great effect on particle trajectory, then further has an impact on flame stability and temperature.

Li, Lei↗
DOE OSTI2024

symfind : Addressing the Fragility of Subhalo Finders and Revealing the Durability of Subhalos

Abstract A major question in ΛCDM is what this theory actually predicts for the properties of subhalo populations. Subhalos are difficult to accurately simulate and to find within simulations, and this propagates into uncertainty in theoretical predictions for satellite galaxies. We present Symfind , a new particle-tracking-based subhalo finder, and demonstrate that it can track subhalos to orders-of-magnitude lower masses than commonly used halo-finding tools, with a focus on Rockstar and consistent-trees . These longer survival times mean that at a fixed peak subhalo mass, we find ≈ 15%–40% more subhalos within the virial radius, R vir , and ≈35%–120% more subhalos within R vir /4 in the Symphony dark-matter-only simulation suite. More subhalos are found as the resolution is increased, in contrast to the Rockstar halo finder, which appears to be converged at smaller subhalo counts. We perform extensive numerical testing. In agreement with idealized simulations, we show that the v max , the maximum circular velocity, is systematically biased low until high resolutions ( n peak ≳ 3 × 10 4 ) are achieved, but that mass loss itself can be resolved at much more modest resolutions ( n peak ≳ 4 × 10 3 ). We show that Rockstar converges to false solutions for the mass function, radial distribution, and disruption masses of subhalos. We argue that our new method can trace resolved subhalos until the point of typical galaxy disruption without invoking post hoc orphan modeling. We outline a concrete set of steps for determining whether other subhalo finders meet the same criteria. We publicly release Symfind catalogs and particle data for the Symphony simulation suite at http://web.stanford.edu/group/gfc/symphony .

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2020

Dose-Dependent Sensorimotor Impairment in Human Ocular Tracking After Acute Low-Dose Alcohol Administration

Changes in oculomotor behaviours are often used as metrics of sensorimotor disruption due to ethanol (EtOH); however, previous studies have focused on deficits at blood-alcohol concentrations (BACs) above about 0.04%.We investigated the dose dependence of the impairment in oculomotor and ocular behaviours caused by EtOH administration across a range of ultra-low BACs (≤0.035%).We took repeated measures of oculomotor and ocular performance from sixteen participants, both pre- and post-EtOH administration. To assess the neurological impacts across a wide range of brain areas and pathways, our protocol measured 21 largely independent performance metrics extracted from a range of behavioural responses ranging from ocular tracking of radial step-ramp stimuli, to eccentric gaze holding, to pupillary responses evoked by light flashes. Our results show significant impairment of pursuit and visual motion processing at 0.015% BAC, reflecting degraded neural processing within extrastriate cortical pathways. However, catch-up saccades largely compensate for the tracking displacement shortfall caused by low pursuit gain, although there still is significant residual retinal slip and thus degraded dynamic acuity. Furthermore, although saccades are more frequent, their dynamics are more sluggish (i.e. show lower peak velocities) starting at BAC levels as low as 0.035%. Small effects in eccentric gaze holding and no effect in pupillary response dynamics were observed at levels below 0.07%, showing the higher sensitivity of the pursuit response to very low levels of blood alcohol, under the conditions of our study.

visual motion processing↗
DOE OSTI2026

Regional Earthquake Ground Motion Simulations for Southern California With EQSIM: Insights From the 2008 Chino Hills, 2024 Highland Park, and 2021 Carson Earthquakes

This study presents physics-based, 3D simulations using the EQSIM framework for several earthquakes in the Los Angeles region. The primary objective was to assess the ability of deterministic physics-based ground motion simulations to reproduce the observed motions from historical events. The selected events included the mathematical equation M w 5.4 2008 Chino Hills, the mathematical equation M w 4.4 2024 Highland Park, and the mathematical equation M w 4.3 2021 Carson events. The simulated motions were evaluated by comparing the recorded and simulated seismograms, as well as the Fourier amplitude spectra, across multiple seismic stations. The SCEC 3D velocity model, CVM-S4.26.M01, was used to represent the regional geology, and ground motion simulations were carried out with a resolution of up to 5 Hz. The results indicate that the simulated motions captured the recorded motions up to approximately 4 Hz. While careful iterations regarding source parameters and corner frequencies were required, and, for the case of the Highland Park event, some of the near-source stations had relatively low accuracy, the present study established a positive step toward the utilization of physics-based simulations in practical applications. The computational efficiencies exhibited by EQSIM, especially on GPU clusters, further supported this assertion, as wall-clock times of simulations involving more than 10 billion grid points were as low as mathematical equation minutes. This permits ensemble simulations for a considered scenario event so that modeling uncertainties (e.g., source and geology) can be bracketed.

EQSIM↗