Search NASA⌕ Search

Engineering topics

McCallen, David

Publications and source records attributed to McCallen, David.

15 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Seismic response of vertical dry storage casks under three-dimensional earthquake motions

Ensuring the long-term seismic safety of dry storage casks (DSCs) is becoming increasingly critical as these systems evolve from temporary to de facto permanent repositories for spent nuclear fuels. Traditional seismic soil–structure interaction (SSI) assessment methods use one-dimensional deconvolution or simplified boundary conditions to model incident waves. Although computationally appealing, simplifying assumptions may alter the seismic risk by neglecting the full complexity of three-dimensional (3D) wave propagation effects. To address this challenge, this paper introduces a novel high-fidelity computational framework that leverages the Domain Reduction Method (DRM) with perfectly matched layers (PML) to accurately transfer complex, 3D seismic wavefields from regional-scale fault-rupture simulations into local-scale finite element models of DSCs. Using broadband, physics-based ground motions from a generic M w 7.0 strike-slip event, both single-cask and multi-cask configurations were investigated under near- and far-field conditions. Emphasis is placed on capturing complex SSI, spatial variability in the ground motion, and nonlinear phenomena such as cask rocking and sliding. Numerical results demonstrate that near-field conditions, where forward directivity and fling-step effects dominate, lead to significantly higher DSC rocking and sliding. Far-field cases, by contrast, generally exhibit modest responses. Incorporating SSI tends to amplify or alter DSC response spectra and introduce response variability, which underscores the need for site-specific evaluations and robust modeling approaches to ensure the seismic integrity of DSCs in interim spent fuel storage installations.

Das, Tonmoy↗
DOE OSTI2024

Ground-motions site and event specificity: Insights from assessing a suite of simulated ground motions in the San Francisco Bay Area

This article presents the results of a research that is part of a larger collaborative effort between the Lawrence Berkeley National Laboratory and the Pacific Earthquake Engineering Research Center, funded by the US Department of Energy Office of Cybersecurity, Energy Security and Emergency Response. The main objective of this study is to assess a suite of near and far-field simulated ground motions obtained from 20 realizations of an M7 Hayward Fault earthquake in the San Francisco Bay Area, California USA, and inform the selection of rupture simulation parameters leading to strong motions. To this aim, comparisons are conducted with NGA-W2 and directivity ground-motion models and a selected population of records. An archetypal steel moment-resisting frame is utilized to assess infrastructure response distributions. The analyses carried out for each simulated event and subdomain with consistent properties in terms of shallow shear-wave velocity proved to be instrumental for better interpreting the differences between simulated motions and empirical models. The main reasons identified for variances between simulations and empirical relationships included (1) directivity effects fully captured by the simulations across the full breadth of rupture models; (2) site vicinity to ruptures that incorporate large-slip patches, particularly if these are in the forward-directivity direction; and (3) presence of geologic structures that can “trap” seismic waves and produce ground motions with large amplitude and long signal duration. The analyses carried out in this work provide a path for interpreting ground-motion site and event specificity obtained from a suite of physics-based simulations, differing only in the rupture model characterization, to inform the selection of simulation scenarios for site-specific engineering analyses under strong excitations. Evidence from this work points to the possibility that current hazard models may underestimate ground-motion intensities in areas where the combined effect of directivity and site conditions results in large ground-motion amplitudes.

58 GEOSCIENCES↗
DOE OSTI2024

Transformational Regional-Scale Earthquake Simulations with the DOE EarthQuake SIMulation Exascale Framework

Earthquakes present worldwide risk to economic and human safety. The 2023 earthquakes in Turkiye provided a reminder of the potential for catastrophic consequences with 50,700 deaths and 15.7 million people affected. The ability to predict ground motions and infrastructure damage for earthquakes continues to be a challenging problem for scientists and engineers. Until now, estimates of ground motions have been performed empirically by looking at sparse data from past earthquakes. This approach can provide statistical information on intensity amplitudes but cannot inform site-specific ground motions essential to developing the most effective resilience. Interest has grown in large-scale computational models to simulate earthquakes at regional scale. The U.S. Department of Energy EarthQuake SIMulation (EQSIM) framework was developed for regional-scale earthquake simulations at unprecedented fidelity, taking advantage of emerging GPU-accelerated systems. This article describes the EQSIM workflow and demonstrates regional-scale simulations with the new computational capability available to scientists in their quest to mitigate future disasters.

58 GEOSCIENCES↗
DOE OSTI2024

Applicability of 1D site response analysis to shallow sedimentary basins: A critical evaluation through physics-based 3D ground motion simulations

One-dimensional site response analysis (1D SRA) remains the standard practice in considering the effect of local soil deposits and predicting site-specific ground motions, although its range of applicability to realistic seismic wavefields is still in question. In this 1D approach, horizontal and vertical ground shaking are assumed to be induced by vertically propagating shear and compressional waves, respectively. A recent study based on analytical two-dimensional (2D) plane waves and simple point source earthquake simulations has shown two mechanistic limitations in this 1D modelling technique for general inclined seismic waves, that is, systematic over-prediction of the vertical motion and wave trapping in the 1D soil column. In this article, we evaluate in detail the applicability of this 1D modelling approach to realistic three-dimensional (3D) simulated seismic wavefields in shallow sedimentary basins. Linear-viscoelastic 1D SRA predictions using two types of input motions that are commonly used in practice—rock outcrop and in-column motions, are compared with the reference true site response results from 3D earthquake simulations in terms of various measures in the frequency and time domain. It is shown that the horizontal motion in the 3D seismic wavefield exhibits dominant shear wave propagation phenomenon, while the vertical motion is a combined effect of compressional and shear waves and can be over-predicted by the 1D approach when the incident seismic waves are inclined. Direct evidence of the wave refraction process that leads to the vertical motion over-prediction is provided. 1D SRA with in-column inputs can yield motions that have significantly longer duration compared to the true 3D site response solution due to trapped waves, casting in doubt the frequent need for increased soil damping in existing site studies to compensate for wave attenuation due to scattering alone. Sensitivity investigation on the increase of soil profile damping by a multiplier D mul shows D mul values compatible with those found in the literature for both horizontal and vertical motions. It is shown that the level of D mul optimized for a best match of the spectral acceleration is dependent on the characteristic of the input motion and a larger D mul is typically required for the vertical component. In contrast, 1D SRA with outcrop motions predicts motions with shorter significant duration due to its inability to capture the basin-edge generated surface waves. A suite of ground motion simulations was performed to assess the sensitivity of the observations to the basin geologic structure including the velocity gradient, rock-basin impedance contrast and basin depth. The analysis results show that the accuracy of the simplified 1D procedure is dependent on the wavefield composition of both the input motions and the true 3D site response solution. While the horizontal motions in shallow sedimentary basins can, to the first order, be reasonably captured by the simplified 1D approach, 1D SRA for the vertical component is in general not reliable and contributions from inclined shear waves should be accounted for in site-specific evaluation of the vertical design ground motion.

1D assumption↗
DOE OSTI2024

Regional-scale fault-to-structure earthquake simulations with the EQSIM framework: Workflow maturation and computational performance on GPU-accelerated exascale platforms

Continuous advancements in scientific and engineering understanding of earthquake phenomena, combined with the associated development of representative physics-based models, is providing a foundation for high-performance, fault-to-structure earthquake simulations. However, regional-scale applications of high-performance models have been challenged by the computational requirements at the resolutions required for engineering risk assessments. The EarthQuake SIMulation (EQSIM) framework, a software application development under the US Department of Energy (DOE) Exascale Computing Project, is focused on overcoming the existing computational barriers and enabling routine regional-scale simulations at resolutions relevant to a breadth of engineered systems. This multidisciplinary software development—drawing upon expertise in geophysics, engineering, applied math and computer science—is preparing the advanced computational workflow necessary to fully exploit the DOE’s exaflop computer platforms coming online in the 2023 to 2024 timeframe. Achievement of the computational performance required for high-resolution regional models containing upward of hundreds of billions to trillions of model grid points requires numerical efficiency in every phase of a regional simulation. This includes run time start-up and regional model generation, effective distribution of the computational workload across thousands of computer nodes, efficient coupling of regional geophysics and local engineering models, and application-tailored highly efficient transfer, storage, and interrogation of very large volumes of simulation data. This article summarizes the most recent advancements and refinements incorporated in the workflow design for the EQSIM integrated fault-to-structure framework, which are based on extensive numerical testing across multiple graphics processing unit (GPU)-accelerated platforms, and demonstrates the computational performance achieved on the world’s first exaflop computer platform through representative regional-scale earthquake simulations for the San Francisco Bay Area in California, USA.

58 GEOSCIENCES↗
DOE OSTI2023

Simulation-based characterization of the variability of earthquake risk to buildings in the near-field

Recent advancements in high performance computing platforms and computational workflow for regional-scale simulations are enabling unprecedented modeling of fault-to-structure earthquake processes. Regional simulations resolving ground motions at frequencies relevant to engineered systems are becoming computationally viable and provide a new capability to improve understanding of the geographical distribution and intensity of risk to buildings and critical infrastructure. As computational capabilities advance, it is essential to move beyond illustrative single rupture realizations for scenario earthquake events towards the development of a full suite of rupture realizations that appropriately characterize the range of risk to building systems. The work described in this article investigates the application of a suite of fault rupture realizations with the objective of assessing near-fault, site-specific seismic demand variability for building structures. A representative high-performance regional-scale computational model is utilized to execute ground motion and building response simulations based on 18 kinematic rupture realizations of an M7 strike-slip scenario earthquake. The fault rupture models for the scenario earthquake are created by systematically perturbing the hypocenter location and stochastically generating rupture parameters (slip, rise time, rake angle) to represent a breadth of ground motion intensities resulting from the spatial and temporal variabilities of an earthquake rupture process. The resulting seismic demand variability for three-story (short period) and forty-story (long period) steel moment-resisting frame buildings is characterized in terms of the median and distribution of peak inter-story drift ratio for a range of near-fault sites. The full suite of 18 fault rupture realizations and approximately 280,000 nonlinear dynamic building simulations indicate that the three-story building undergoes higher median seismic demand and significantly greater variability of demand at a given site than the forty-story building, which has important implications for the level of certainty in predicting building performance during an earthquake. The simulations performed provide deeper insight into the relationship between fault rupture parameterization and building response, which is essential information for developing a representative suite of rupture realizations for specific earthquake scenarios.

58 GEOSCIENCES↗
DOE OSTI2023

Characteristics and selection of near-fault simulated earthquake ground motions for nonlinear analysis of buildings

Earthquake-induced ground shaking near rupturing faults is highly sensitive to the rupture characteristics, seismic wave propagation patterns and site conditions, and field recordings of near-fault shaking are relatively sparse. These challenges complicate the assessment of the seismic performance of near-fault structures. A common approach to representing near-fault ground motion in engineering analysis is to explicitly consider and select records with strong directivity pulses (pulse records). We use three-dimensional high-resolution physics-based earthquake simulations to test this approach in the context of scenario-based ground motion record selection, and to study the important characteristics of near-fault ground shaking. We highlight the deficiencies associated with classifying near-fault simulated records as “pulse” or “non-pulse,” based on the presence of a single dominating pulse in the velocity time history. We show that this approach is inadequate for characterizing near-fault shaking on soft soils which can be dominated by both forward rupture directivity and basin amplification effects. We conduct ground motion selection experiments for the analysis of near-fault structures with and without explicit classification of the pulse features in the records, and evaluate the bias in the predicted structural demands. We find that the maximum interstory drift demands on building structures imposed by unscaled site-specific simulated ground motion records selected based on relevant spectral shape features are not sensitive to the classification of records as pulse/non-pulse. Therefore, with regard to predicting the maximum interstory drifts in near-fault buildings, we do not find justification for the binary pulse classification of near-fault records.

42 ENGINEERING↗
DOE OSTI2021

Refinements to the Graves–Pitarka Kinematic Rupture Generator, Including a Dynamically Consistent Slip-Rate Function, Applied to the 2019 Mw 7.1 Ridgecrest Earthquake

The main objective of this study is to develop physics–based constraints on the spatiotemporal variation of the slip–rate function using a simplified dynamic rupture model. First, we performed dynamic rupture modeling of the 2019 M w 7.1 Ridgecrest, California, earthquake, to analyze the effects of depth–dependent stress and material friction on slip rate. Then, we used our modeling results to guide refinements to the slip–rate function that were implemented in the Graves–Pitarka kinematic rupture generation technique. The dynamic ruptures were computed on a surface–rupturing, planar strike–slip fault that includes a weak (negative to low–stress–drop) zone in the upper 4 km of the crust. Below the weak zone, we placed high–stress–drop patches designed to mirror the large–slip areas seen in various rupture model inversions of the event. The locations of the high–stress–drop patches and the hypocenter were varied in multiple realizations to investigate how changing the dynamic conditions affected the resulting rupture kinematics, in particular, the slip rate. From these simulations, we observed a systematic change in the shape of the slip–rate function from Kostrov type below the weak zone to a predominantly symmetric shape within the weak zone, along with a depth–dependent reduction of peak slip rate. We generalized these shallow rupture features into a depth–dependent parametric variation of the slip–rate function and implemented it in the Graves–Pitarka kinematic rupture model generator. The performance of the updated kinematic approach was then verified in 0–4 Hz simulations of the M w 7.1 Ridgecrest earthquake, which showed that incorporating the depth–dependent variation in the shape of the slip–rate function improves the fit to the observed near–fault ground motions in the 0.5–3 s period range.

58 GEOSCIENCES↗
DOE OSTI2021

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↗
DOE OSTI2021

Coupling of regional geophysics and local soil-structure models in the EQSIM fault-to-structure earthquake simulation framework

Accurate understanding and quantification of the risk to critical infrastructure posed by future large earthquakes continues to be a very challenging problem. Earthquake phenomena are quite complex and traditional approaches to predicting ground motions for future earthquake events have historically been empirically based whereby measured ground motion data from historical earthquakes are homogenized into a common data set and the ground motions for future postulated earthquakes are probabilistically derived based on the historical observations. This procedure has recognized significant limitations, principally due to the fact that earthquake ground motions tend to be dictated by the particular earthquake fault rupture and geologic conditions at a given site and are thus very site-specific. Historical earthquakes recorded at different locations are often only marginally representative. There has been strong and increasing interest in utilizing large-scale, physics-based regional simulations to advance the ability to accurately predict ground motions and associated infrastructure response. However, the computational requirements for simulations at frequencies of engineering interest have proven a major barrier to employing regional scale simulations. In a U.S. Department of Energy Exascale Computing Initiative project, the EQSIM application development is underway to create a framework for fault-to-structure simulations. This framework is being prepared to exploit emerging exascale platforms in order to overcome computational limitations. This article presents the essential methodology and computational workflow employed in EQSIM to couple regional-scale geophysics models with local soil-structure models to achieve a fully integrated, complete fault-to-structure simulation framework. Here, the computational workflow, accuracy and performance of the coupling methodology are illustrated through example fault-to-structure simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Variability of near-fault seismic risk to reinforced concrete buildings based on high-resolution physics-based ground motion simulations

Broadband physics-based simulated earthquake ground motions are utilized to characterize the regional-scale seismic risk to modern reinforced concrete (RC) structures. A highly dense dataset of ground motions covering a 100-km urn:x-wiley:00988847:media:eqe3413:eqe3413-math-0001 40-km domain was generated using kinematic fault rupture models with varying rupture characteristics to represent shallow crustal earthquakes and resolved up to frequencies of 5 Hz. Over 40,000 nonlinear response history simulations of short- and mid-rise RC special moment frame buildings were conducted using simulation models that are capable of representing nonlinear behavior and component deterioration effects. The spatial variability of structural risk within a single earthquake scenario and between different rupture scenarios is examined, and the regions of strongest directivity effects and highest structural demands are identified. The structural demands may vary by factor of up to 8.0 at very short distances from the fault, and the large dispersion in the demands decreases significantly beyond a distance of 15 km. The interstory drift and member rotation demands are substantially impacted by important features of the geological structure and the characteristics of the rupture scenarios, particularly the presence of localized high-slip regions. The frequency characteristics of the structures are found to play an important role in determining the effects of near-fault ground motions on the structural response and expected damage. The results of this study suggest that the simulated ground motions, particularly those generated using the hybrid rupture approach, may offer reasonable structural risk estimates for low-frequency structures and conservative estimates for high-frequency structures.

58 GEOSCIENCES↗
DOE OSTI2020

EQSIM—A multidisciplinary framework for fault-to-structure earthquake simulations on exascale computers, part II: Regional simulations of building response

The existing observational database of the regional-scale distribution of strong ground motions and measured building response for major earthquakes continues to be quite sparse. As a result, details of the regional variability and spatial distribution of ground motions, and the corresponding distribution of risk to buildings and other infrastructure, are not comprehensively understood. Utilizing high-performance computing platforms, emerging high-resolution, physics-based ground motion simulations can now resolve frequencies of engineering interest and provide detailed synthetic ground motions at high spatial density. This provides an opportunity for new insight into the distribution of infrastructure seismic demands and risk. In the work presented herein, the EQSIM fault-to-structure computational framework described in a companion paper, McCallen et al., is employed to investigate the regional-scale response of buildings to large earthquakes. A representative M = 7.0 strike-slip event is used to explore the distribution and amplitude of building demand, and comparisons are made between building response computed with fault-to-structure simulations and building response computed with existing measured near-fault earthquake records. New information on the distribution and variability of building response from high-performance parallel simulations is described and analyzed, and favorable first comparisons between building response predicted with both fault-to-structure simulations and real ground motions records are presented.

58 GEOSCIENCES↗
DOE OSTI2020

EQSIM—A multidisciplinary framework for fault-to-structure earthquake simulations on exascale computers part I: Computational models and workflow

Computational simulations have become central to the seismic analysis and design of major infrastructure over the past several decades. Most major structures are now “proof tested” virtually through representative simulations of earthquake-induced response. More recently, with the advancement of high-performance computing (HPC) platforms and the associated massively parallel computational ecosystems, simulation is beginning to play a role in increased understanding and prediction of ground motions for earthquake hazard assessments. However, the computational requirements for regional-scale geophysics-based ground motion simulations are extreme, which has restricted the frequency resolution of direct simulations and limited the ability to perform the large number of simulations required to numerically explore the problem parametric space. In this article, recent developments toward an integrated, multidisciplinary earth science-engineering computational framework for the regional-scale simulation of both ground motions and resulting structural response are described with a particular emphasis on advancing simulations to frequencies relevant to engineered systems. This multidisciplinary computational development is being carried out as part of the US Department of Energy (DOE) Exascale Computing Project with the goal of achieving a computational framework poised to exploit emerging DOE exaflop computer platforms scheduled for the 2022–2023 timeframe.

58 GEOSCIENCES↗
DOE OSTI2020

Validation of (not-historical) large-event near-fault ground-motion simulations for use in civil engineering applications

Ground-motion simulations generated from physics-based wave propagation models are gaining increasing interest in the engineering community for their potential to inform the performance-based design and assessment of infrastructure residing in active seismic areas. A key prerequisite before the ground-motion simulations can be used with confidence for application in engineering domains is their comprehensive and rigorous investigation and validation. This article provides a four-step methodology and acceptance criteria to assess the reliability of simulated ground motions of not historical events, which includes (1) the selection of a population of real records consistent with the simulated scenarios, (2) the comparison of the distribution of Intensity Measures (IMs) from the simulated records, real records, and Ground-Motion Prediction Equations (GMPEs), (3) the comparison of the distribution of simple proxies for building response, and (4) the comparison of the distribution of Engineering Demand Parameters (EDPs) for a realistic model of a structure. Specific focus is laid on near-field ground motions (<10km) from large earthquakes (M w 7), for which the database of real records for potential use in engineering applications is severely limited. The methodology is demonstrated through comparison of (2490) near-field synthetic records with 5 Hz resolution generated from the Pitarka et al (2019) kinematic rupture model with a population of (38) pulse-like near-field real records from multiple events and, when applicable, with NGA-W2 GMPEs. Lastly, the proposed procedure provides an effective method for informing and advancing the science needed to generate realistic ground-motion simulations, and for building confidence in their use in engineering domains.

42 ENGINEERING↗