Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Earthquakes”

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,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Geomechanical characterization of Rock Valley carbonates

In Rock Valley, Nevada, the USA, past earthquake sequences and shallow faulting in Paleozoic carbonates remain poorly understood. Carbonates are typically more ductile rocks that do not experience large stress drops or fracture coalescence, which contradicts previous observations in the region. As part of the Source Physics Experiment, an effort has been made to experimentally characterize the petrophysical and geomechanical properties of carbonates from Rock Valley. Well core and outcrop samples were used to determine the difference between shallow-buried Tertiary limestones and deeply buried Paleozoic limestones and dolostones. The carbonates possessed porosities between 3 and 9%, with V P and V S in the Tertiary carbonates around 2700 and 1800 m/s, respectively, and 6000 and 3100 m/s, respectively, in the Paleozoic carbonates. Microstructural characterization revealed that the Paleozoic carbonates contained significantly more pre-existing damage and alteration than the Tertiary carbonates, though the deformation varies from localized to diffuse. Unconfined Brazilian tests and triaxial tests with confining pressures between 0 and 50 MPa showed that, although samples all experienced failure, the dolostones typically failed at greater stresses and possessed greater E/ν ratios than the limestones. Furthermore, the Hoek–Brown failure criterion was used to construct failure envelopes with the compressive and tensile failure tests. Velocity, density, and porosity measurements were used to construct a hypothetical velocity-depth profile and compare the results with the theoretical petrophysical measurements. Using brittleness analysis, the likely failure conditions were determined to be low-porosity dolomitic rock for fault nucleation in the Paleozoic basement.

Brittleness↗
DOE OSTI2022

Rapid Evaluation and Response to Impacts on Critical End-Use Loads Following Natural Hazard-Driven Power Outages: A Modular and Responsive Geospatial Technology

The disparate nature of data for electric power utilities complicates the emergency recovery and response process. The reduced efficiency of response to natural hazards and disasters can extend the time that electrical service is not available for critical end-use loads, and in extreme events, leave the public without power for extended periods. This article presents a methodology for the development of a semantic data model for power systems and the integration of electrical grid topology, population, and electric distribution line reliability indices into a unified, cloud-based, serverless framework that supports power system operations in response to extreme events. An iterative and pragmatic approach to working with large and disparate datasets of different formats and types resulted in improved application runtime and efficiency, which is important to consider in real time decision-making processes during hurricanes and similar catastrophic events. This technology was developed initially for Puerto Rico, following extreme hurricane and earthquake events in 2017 and 2020, but is applicable to utilities around the world. Given the highly abstract and modular design approach, this technology is equally applicable to any geographic region and similar natural hazard events. In addition to a review of the requirements, development, and deployment of this framework, technical aspects related to application performance and response time are highlighted.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

A general approach to seismic inversion with automatic differentiation

Imaging Earth structure or seismic sources from seismic data involves minimizing a target misfit function, and is commonly solved through gradient-based optimization. The adjoint-state method has been developed to compute the gradient efficiently; however, its implementation can be time-consuming and difficult. We develop a general seismic inversion framework to calculate gradients using reverse-mode automatic differentiation. The central idea is that adjoint-state methods and reverse-mode automatic differentiation are mathematically equivalent. Here, the mapping between numerical PDE simulation and deep learning allows us to build a seismic inverse modeling library, ADSeismic, based on deep learning frameworks, which supports high performance reverse-mode automatic differentiation on CPUs and GPUs. We demonstrate the performance of ADSeismic on inverse problems related to velocity model estimation, rupture imaging, earthquake location, and source time function retrieval. ADSeismic has the potential to solve a wide variety of inverse modeling applications within a unified framework.

58 GEOSCIENCES↗
DOE OSTI2023

Effect of a micro-scale dislocation pileup on the atomic-scale multi-variant phase transformation and twinning

In this paper, we perform concurrent atomistic–continuum (CAC) simulations to assess the contribution of the internal stress induced by the microscale dislocation pileup at an atomically structured interface to the atomic-scale phase transformations (PTs), reverse PTs, and twinning. The main novelty of this work is to unify the atomistic description of the interface and the coarse-grained (CG) description of the lagging dislocations away from the interface within one single framework. Our major findings are: (a) the interface dynamically responds to a pileup by forming steps/ledges, the height of which is proportional to the number of dislocations arriving at the interface; (b) the pileup-induced internal stress concentration profile follows neither the classical Eshelby model nor the super-dislocation model alone, but a combination of them; (c) when the pre-sheared sample is compressed, a direct square-to-hexagonal PT occurs ahead of the pileup tip and eventually grows into a wedge shape. The two variants of the hexagonal phases form a twin with respect to each other; (d) upon a further increase of the loading, part of the newly formed hexagonal phase transforms back to the square phase. The square product phase resulting from this reverse PT forms a twin with respect to the initial square phase. All phase boundaries (PBs) and twin boundaries (TBs) are stationary and correspond to zero thermodynamic Eshelby driving forces; and (e) the microscale dislocation pileup-induced internal shear stress and the structural change at the atomic-scale interface reduces the stress required for initiating a PT by a factor of 5.5, comparing with that in the sample containing no dislocations. This work is the first characterization of the behavior of PTs/twinning resulting from the reaction between a microscale dislocation slip and an atomically structured interface. In conclusion, the gained knowledge will advance our understanding of how the multi-phase material behaves in many complex physical processes, such as the synthesis of multi-phase high-entropy alloys or superhard ceramics under high-pressure torsion, deep mantle earthquakes in geophysics, and so on, which all involve dislocation slip, PTs, twinning, and their interactions across from the atomistic to the microscale and beyond.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Understanding collective human movement dynamics during large-scale events using big geosocial data analytics

Conventional approaches for modeling human mobility pattern often focus on human activity and movement dynamics in their regular daily lives and cannot capture changes in human movement dynamics in response to large-scale events. With the rapid advancement of information and communication technologies, many researchers have adopted alternative data sources (e.g., cell phone records, GPS trajectory data) from private data vendors to study human movement dynamics in response to large-scale natural or societal events. Big geosocial data such as georeferenced tweets are publicly available and dynamically evolving as real-world events are happening, making it more likely to capture the real-time sentiments and responses of populations. However, precisely-geolocated geosocial data is scarce and biased toward urban population centers. In this research, we developed a big geosocial data analytical framework for extracting human movement dynamics in response to large-scale events from publicly available georeferenced tweets. The framework includes a two-stage data collection module that collects data in a more targeted fashion in order to mitigate the data scarcity issue of georeferenced tweets; in addition, a variable bandwidth kernel density estimation(VB-KDE) approach was adopted to fuse georeference information at different spatial scales, further augmenting the signals of human movement dynamics contained in georeferenced tweets. To correct for the sampling bias of georeferenced tweets, we adjusted the number of tweets for different spatial units (e.g., county, state) by population. To demonstrate the performance of the proposed analytic framework, we chose an astronomical event that occurred nationwide across the United States, i.e., the 2017 Great American Eclipse, as an example event and studied the human movement dynamics in response to this event. Finally, this analytic framework can easily be applied to other types of large-scale events such as hurricanes or earthquakes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Review of analytical studies on seismic fluid-structure interaction of base-supported cylindrical tanks

Analytical solutions for seismic fluid-structure interaction (FSI) of tanks are used by analysts and engineers for a number of purposes, including preliminary sizing, risk assessment, and verification of numerical models. Analytical solutions for seismic FSI of tanks were first developed in the 1930s and have been extended and/or modified in the following nine decades to accommodate different seismic inputs and boundary conditions. This report reviews published analytical solutions for rigid and flexible, base-supported, cylindrical tanks subjected to unidirectional horizontal motion of a small amplitude. The studies parsed the fluid-structure (FSI) response into an impulsive and a convective component, and addressed frequencies of lateral motions of the tank (i.e., impulsive frequencies), hydrodynamic pressures, frequencies of waves (i.e., convective frequencies), wave heights, and reactions at the base. The solutions were categorized into exact and approximate solutions based on the methodologies used. Herein, FSI responses are calculated and compared for a range of tank dimensions using the analytical solutions from the different studies. The responses presented here are normalized to be unitless, and can be used for tanks with different dimensions and mechanical properties, and subjected to different input motions. A sample steel water tank is analyzed using the normalized solutions for two earthquake ground motions with different frequency contents. Ignoring tank flexibility may significantly underestimate seismic responses of the tank and its contained fluid. The approximate solutions, which have been widely used for decades, may not be sufficiently accurate in the modern era. Calculation errors in the reviewed studies are identified and corrected, and differences are quantified for the sample tank.

42 ENGINEERING↗
DOE OSTI2025

Neural operators for stochastic modeling of nonlinear structural system response to natural hazards

Traditionally, neural networks have been employed to learn the mapping between finite-dimensional Euclidean spaces. However, recent research has opened up new horizons, focusing on the utilization of deep neural networks to learn operators capable of mapping infinite-dimensional function spaces. Here, in this work, we employ two state-of-the-art neural operators, the deep operator network (DeepONet) and the Fourier neural operator (FNO) for the prediction of the nonlinear time history response of structural systems exposed to natural hazards, such as earthquakes and windstorms. Specifically, we propose two architectures, a self-adaptive FNO and a fast Fourier transform-based DeepONet (DeepFNOnet), where we employ a FNO beyond the DeepONet to learn the discrepancy between the ground truth and the solution predicted by the DeepONet. To demonstrate the efficiency and applicability of the architectures, two problems are considered. In the first, we use the proposed model to predict the seismic nonlinear dynamic response of a six-story shear building subject to stochastic ground motions. In the second problem, we employ the operators to predict the wind-induced nonlinear dynamic response of a high-rise building while explicitly accounting for the stochastic nature of the wind excitation. In both cases, the trained metamodels achieve high accuracy while being orders of magnitude faster than their corresponding high-fidelity models.

DeepONet↗
DOE OSTI2020

Plate motion in sheared granular fault system

Plate motion near the fault gouge layer, and the elastic interplay between the gouge layer and the plate under stick-slip conditions, is key to understanding the dynamics of sheared granular fault systems. Here, a two-dimensional implementation of the combined finite-discrete element method (FDEM), which merges the finite element method (FEM) and the discrete element method (DEM), is used to explicitly simulate a sheared granular gouge fault system. In this work we focus on investigating the influence of normal load, driving shear velocity and plate stiffness on the velocities and displacements in the direction parallel to the shear direction (x-direction) measured at locations on the upper and lower plates just adjacent to the gouge. The simulations show that during slip phases the magnitudes of the measured velocities on the upper and lower plates are proportional to the normal load and may be inversely proportional to the square root of the plate's shear modulus. Whereas, the driving shear velocity does not show distinct influence on the measured velocities. Additionally, large slip velocities are generally associated with large macroscopic friction coefficient drops. For the models subjected to smaller normal loads, larger shear velocities and with stiffer shear plates, the same magnitude of slip velocity could cause a larger drop of macroscopic friction coefficient. During stick phases, the velocities of the upper and lower plates are respectively slightly greater and slightly smaller than half of the driving shear velocity and are both in the same direction of shear. The shear strain rate of the gouge is calculated from this velocity difference between the upper and lower plate during stick phases and thus the gouge effective shear modulus can be calculated. The results show that the gouge effective shear modulus increases proportionally with normal load, while the influence of shear velocity and plate stiffness on gouge effective shear modulus is minor. The simulations address the dynamics of a laboratory-scale fault gouge system and may aid in revealing the complexities of earthquake frictional dynamics.

58 GEOSCIENCES↗
DOE OSTI2022

Post-hazard condition assessment of nuclear piping-equipment systems: Novel approach to feature extraction and deep learning

Over the past decade, the use of artificial intelligence techniques in the field of health-monitoring has gained significant interest, especially for structures such as building and bridges. However, applications to industrial systems such as equipment-piping systems in nuclear plants have not been explored. Here, in this paper, it is shown that the existing techniques developed for buildings and bridges cannot be extended directly to equipment-piping systems as the response of such systems is governed by multiple localized modes unlike that in buildings and bridges. This paper proposes a new approach that consists of three key aspects: (i) a novel vector of degradation-sensitive features extracted from measured data, (ii) using a deep Artificial Neural Network (ANN) for diagnosis of degradation location and degradation severity, and (iii) consideration of uncertainty in degradation severity when training the ANN. Degradation in piping-equipment systems can occur due to flow-accelerated erosion and corrosion. These locations can potentially exhibit damage such as localized yielding or initiation of cracking due to an external event such as an earthquake. Moreover, such locations can at times go undetected by current inspection techniques. Therefore, a robust framework is needed for detection of degradation after a seismic event. This manuscript proposes a proof-of-concept framework, which utilizes data collected from sensors to generate a deep ANN database for predicting degraded locations and severity in a piping-equipment system. Degradation severity is classified as minor, moderate, and severe. In the suggested methodology, a novel vector of degradation-sensitive features is extracted from the sensor data to train the ANN. A simple piping-equipment system is selected to demonstrate feature extraction as a means to simplify pattern recognition, explore the design and parameters of an ANN, and develop a sensor placement strategy. The effectiveness of the proposed framework is demonstrated on a realistic primary safety system of a two-loop nuclear reactor. It is shown that the proposed post-hazard condition assessment framework is able to detect degraded locations along with the severity levels, including minor degradation, with considerably higher accuracy.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Osmotic control of the spacing of parallel shear cracks in shale growing subcritically in geologic past

The geological genesis of natural cracks in sedimentary rocks such as shale is a problem that needs to be understood to improve the technology of hydraulic fracturing as well as deep sequestration of harmful fluids. Why are the vertical natural cracks roughly parallel and equidistant, and why is the spacing roughly 10 cm rather than 1 cm or 100 cm? Fracture mechanics of critical cracks cannot answer this question. Neither can the material heterogeneity. The growth of critical parallel cracks is impossible because the relative crack face displacements would immediately localize into one crack, leading to an earthquake. The cracks must have formed, on the tectonic time scale, by a slow growth of subcritical shear cracks governed by the Charles-Evans law. The idea advanced here is that what controls the crack spacing is the balance between the reduction, due to shear dilatancy, of the concentration of ions such as Na + and Cl - in each fracture process zone (PFZ), which decelerates the cracks, and the restoration of ion concentration by diffusion of ions from the space between the cracks into the FPZ. This diffusion of water is driven mainly by the osmotic pressure gradient, which offsets the deceleration and depends strongly on the crack spacing. A simple analytical solution of the steady state is rendered possible by approximating the ion concentration profiles between adjacent cracks by parabolic arcs. Applying this theory to Woodford shale yields the approximate crack spacing of 10 cm, which is realistic. Furthermore, the stability of unlimited parallel mode II frictional crack growth is proven by examining the second variation of the free energy. Water concentration drop in the FPZ due to shear dilatancy and its restoration by water diffusion from the inter-crack space have similar effect, although probably much weaker.

42 ENGINEERING↗
DOE OSTI2020

Seismic tomography of compressional wave velocity and attenuation structure for Makushin Volcano, Alaska

Located in the central portion of the Alaska-Aleutian subduction zone, Makushin Volcano is among the most active volcanoes in the United States and has been classified as a very high threat by the National Volcano Early Warning System based on its eruptive history and proximity to the City of Unalaska and international air routes. In 2015, five standalone seismic stations and three mini seismic arrays of 15 stations each were deployed on Unalaska Island to supplement the Alaska Volcano Observatory permanent seismic network. This temporary array was operational for one year. Taking advantage of the increased azimuthal coverage and the array's increased earthquake detection capability, we developed body-wave Vp and Vp/Vs seismic images of the velocity structure beneath the volcano. Our results show a complex structure with the upper 4 km of the crust. The shallow high-Vp features possibly delineate remnant magma pathways or conduits. Low-Vp regions are found east of the caldera at approximately 4–8 km depth. This is in agreement with previous seismic tomographic work and InSAR-based geodetic models, which had identified this region as a possible long-term source of magma. A three-dimensional compressional wave attenuation model complemented the velocity models. In general, we observe regions of corresponding low P-wave velocity and high attenuation east of the summit caldera between 5 and 7 km depth, supporting the presence of a melt-rich zone. The absolute Vp values (5.5–6.0 km/s) in this region are consistent with the presence of a magma body with a mafic signature. Here, we suggest that a triggering mechanism that involves mafic injections into a more evolved shallow magma chamber is representative of the current magmatic system below the volcano, similarly to what petrological studies have inferred for past eruptions.

58 GEOSCIENCES↗
DOE OSTI2021

Recent seismicity at Ceboruco Volcano (Mexico)

Magma and related hydrothermal fluid movement, and their interaction with solid rock, in active volcanic regions, generate a wide variety of seismic waves whose characterization can mitigate the risk of a potential eruption. Located in the western region of the Trans-Mexican Volcanic Belt, Ceboruco Volcano, whose last eruptive period was 1870–1875, is considered to be one of the most hazardous volcanoes in Mexico. We have conducted a detailed study of the seismicity in the surroundings of Ceboruco's volcanic edifice to assess the current state of this volcano. A dense temporary seismic network with 25 seismic stations in an area of 16 km × 16 km was deployed between November 2016 and July 2017, as part of the P-24 project of the CeMIEGeo consortium; this effort has allowed the detection of 81 earthquakes concentrated beneath the crater with depths between 4 and 8 km. In this study, we observe that the recorded seismicity occurs in swarms, and we specifically identify four sequences that we characterize in detail via the first focal mechanisms available for this volcano. Our results suggest a change in the local seismicity distribution compared to earlier observations, which reported seismic activity near the volcano edifice associated with fluid migration along zones of weakness related to the extensional stresses of the Tepic-Zacoalco rift. The changes in seismic patterns and obtained focal mechanisms are consistent with observed fluid effects at many geothermal sites worldwide, but also could suggest resumption of activity at this currently dormant volcano.

58 GEOSCIENCES↗
DOE OSTI2023

Hydrocarbon production reduces natural methane seeps in the Santa Barbara channel

Abundant natural hydrocarbon seepage in the Santa Barbara Channel includes large amounts of gas, mainly methane. Numerous seeps are focused along a faulted anticlinal axis where in 1982 ARCO placed two large (900 m 2 each) steel gas capture tents on the seabed adjacent to Platform Holly. The seep rate into the tent at high tide is about 5 percent less compared with low tide (2m tidal range). An earthquake in 2013 (M = 4.8) did not affect seepage rates into the steel tents, indicating that gas was not released by fault dilation. The CO 2 and H 2 S content of the captured gas are anomalously low compared to typical Holly gas due to bacterial and sea water-gas interaction at the seep tent. Platform Holly production had an overall effect of reducing gas seepage into the tents over time. Wells beneath the tents at 1-km depth had a much greater influence on the seepage rate than more distal wells. The nearby wells had a direct effect on seepage rate when they were put on production (seepage decreased) versus shut in (seepage increased). The response time was within hours, indicating a direct hydraulic connection between the reservoir at 1 km subsea and the seabed. Permeability of the approximate 1-km fault pathway between the reservoir and the seep tent is calculated to be about 3 × 10 -14 m 2 (30 millidarcys). Gradually, the seep rate declined from peak flow in 1990 (4.5×10 4 m 3 dy-1) to less than 7.1×10 3 m 3 dy -1 in 2013. In June of 2013, a horizontal well was completed high on the structural anticline between Platform Holly and directly beneath the seep tents. Within three weeks of turning on this new well, seepage into the tents completely ceased, ending more than 30 years of natural methane seepage captured by the tents.

58 GEOSCIENCES↗
DOE OSTI2020

Structural modeling and dynamic analysis of condensate storage tanks in nuclear power plants

Condensate storage tanks (CSTs) in nuclear power plants (NPPs) are classified as critical equipment capable of surviving strong shaking in a design basis earthquake to assure the ability to subsequently provide cooling water. Structural modeling and dynamic analysis of CSTs are complicated due to fluid-structure interaction (FSI) and coupling issues, which could have significant effects on seismic response. The capability of a number of approximate models to capture seismic response is investigated using dynamic analysis results from detailed finite element models. Modal analyses and time history analyses are carried out using both 2D and 3D models to predict the dynamic behavior of CSTs as a function of ground motion intensity. Although 2D simplified models can be used to quickly evaluate the dynamic response of CSTs when there is substantial margin to failure, a detailed 3D model is required when it is necessary to examine a limit state associated with the failure mode of the tank. The findings show good agreement between the natural frequencies of convective modes determined from the 2D and 3D CST models. The difference between the frequencies of impulsive modes were found to be equal to or less than 25%.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Validation of numerical models for seismic fluid-structure-interaction analysis of nuclear, safety-related equipment

Seismic design, qualification, and risk assessment of nuclear safety-related vessels filled with liquid will have to consider the interaction between the vessel (tank), the contained liquid, and submerged components, if any. Seismic fluid-structure-interaction (FSI) analysis of nuclear vessels will rely on numerical models, which are required to be verified and validated. This study validates previously-verified numerical models using test data generated from earthquake-simulator experiments involving a cylindrical tank. Two solvers in LS-DYNA for FSI simulations are used for the numerical analysis: Arbitrary-Lagrangian-Eulerian (ALE) and Incompressible Computational Fluid Dynamics (ICFD). Numerical and test results are compared for responses critical to seismic design of advanced reactor vessels: hydrodynamic pressures on the tank wall, reactions at the support, and wave heights of the contained liquid. Analysis is performed for one-, two-, and three-directional seismic inputs with a range of intensity, and rocking motions. The accuracy of the numerical results and different methods for outputting wave heights are discussed. Recommendations for validation of seismic FSI numerical models of advanced reactors are provided. As a result, the validation exercise presented in this study is broadly applicable to cylindrical tanks, regardless of industry sector.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2021

Seismic isolation: A pathway to standardized advanced nuclear reactors

Standardizing advanced nuclear reactors is a pathway to substantially reducing their overnight capital cost and achieving parity with other power sources, including renewables and fossil fuels. The seismic load case has thwarted standardization of nuclear power plants because site-specific seismic hazard and local near-surface geology has triggered soil-structure-interaction analysis, design, equipment qualification, regulatory review, and licensing, ensuring that each build is different. To achieve standardized or site-independent certified advanced reactor designs, the impact of the seismic load case on the engineering and construction cost and time must be substantially mitigated. Seismic isolation is a mature technology that has been used for more than 30 years in non-nuclear sectors to substantially reduce earthquake demands in buildings and other infrastructure. In this paper, seismic isolation is used to enable standardization of advanced reactor designs, aimed at the complete re-use of a site-independent, certified design and repeated procurement of safety-class equipment. A pathway to standardized designs using seismic isolation is demonstrated for two fundamentally different advanced reactors: a molten salt reactor and a high temperature gas reactor. Each reactor building is equipped with three specialized pieces of safety-class equipment, namely, a reactor vessel, a steam generator, and a control rod drive mechanism housing that is attached to the reactor head. Analysis is performed per ASCE and ASME standards to design the buildings and the equipment for two base conditions: conventional (fixed base) and base isolated. The impact of the seismic load case is characterized for the reinforced concrete walls in the buildings and for the equipment, measured using vessel wall thickness and horizontal accelerations. Here the analysis results show that the fixed-base buildings, designed for a site of low seismic hazard (peak ground acceleration, PGA = 0.15 g) could be constructed at a site of much greater seismic hazard (PGA = 0.7 g) if seismic base isolation is employed. Importantly, the scope of the site-specific analysis, design, and qualification would be limited to the seismic isolators and the isolated substructure, drastically reducing plant-specific engineering, review, and licensing, and time to construction start. Regulatory challenges and opportunities with standardized reactor designs are identified.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Mid-height seismic isolation of equipment in nuclear power plants: Numerical simulations and design recommendations

Seismic isolation can mitigate earthquake demands on safety-class equipment and is being considered for application to some advanced nuclear reactors. For tall, slender vessels that could represent an advanced reactor, a steam generator or a heat exchanger, mid-height seismic isolation has been shown to be beneficial and practical. Here, this paper focuses on numerical modeling of a mid-height isolated tall, slender vessel with two primary goals 1) investigate and quantify the benefits of mid-height seismic isolation for a range of support structure stiffnesses, isolation systems, and seismic inputs, and 2) provide analysis and design recommendations for isolation of safety-class equipment. Results of response-history analysis were in good agreement with experimental measurements and demonstrated that mid-height isolation can substantially reduce seismic demands on tall, slender vessels for a range of support structure stiffnesses, isolation systems, and seismic inputs. Importantly, the reductions in horizontal spectral accelerations in the mid-height isolated vessel from the non-isolated condition were not affected by the stiffness of the support structure. Recommendations are made for analysis and design of isolated equipment and testing of isolators used for equipment protection in nuclear facilities, which are also applicable to other industries.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Series of semi-Markov processes to model infrastructure resilience under multihazards

Civil infrastructure systems are subjected to multiple hazards, including natural and anthropogenic, that disrupt their function or the level of service offered. Estimating the function recovery of these systems (or how soon normalcy of operations will be restored) when subjected to repeated hazard events by considering the inter-event dependencies is an important problem in multihazard infrastructure resilience. However, this problem has been less addressed in the field. This paper proposes a series of semi-Markov processes model to capture the inter-event dependencies in infrastructure recovery when subjected to successive hazard events. Recovery after each new hazard event is represented by a unique semi-Markov process that models the reduced recovery rates and the increased recovery times caused by the system’s incomplete recovery from the preceding event. Two novel formulations of the inter-event dependency modeling, namely Maximal Effects Dependency (considers the worst impact of two successive hazard events) and Cumulative Effects Dependency (considers the aggregated impacts of two successive hazard events), are proposed and discussed. The model is demonstrated by considering the following applications: Three-state system subjected to deterministic and random occurrences of identical hazard events; and Multihazard resilience of a building in Charleston, SC, considering earthquake and hurricane hazards. Results indicate that considering inter-event dependencies in recovery modeling can lead to lesser-predicted resilience, thereby affecting resilience-based decision-making.

42 ENGINEERING↗