Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Code Analysis”

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 505 records · Page 28

Sierra/Solid Mechanics 5.16 User's Guide

Sierra/SolidMechanics (Sierra/SM) is a Lagrangian, three-dimensional code for finite element analysis of solids and structures. It provides capabilities for explicit dynamic, implicit quasistatic and dynamic analyses. The explicit dynamics capabilities allow for the efficient and robust solution of models with extensive contact subjected to large, suddenly applied loads. For implicit problems, Sierra/SM uses a multi-level iterative solver, which enables it to effectively solve problems with large deformations, nonlinear material behavior, and contact. Sierra/SM has a versatile library of continuum and structural elements, and a large library of material models. The code is written for parallel computing environments enabling scalable solutions of extremely large problems for both implicit and explicit analyses. It is built on the SIERRA Framework, which facilitates coupling with other SIERRA mechanics codes. This document describes the functionality and input syntax for Sierra/SM.

36 MATERIALS SCIENCE↗

Sierra/SolidMechanics 4.56 User's Guide

Siera/SolidMechanics (Sierra / SM) is a Lagrangian, three-dimensional code for finite element analysis of solids and structures. It provides capabilities for explicit dynamic, implicit quasistatic and dynamic analyses. The explicit dynamics capabilities allow for the efficient and robust solution of models with extensive contact subjected to large, suddenly applied loads. For implicit problems, Sierra / SM uses a multi-level iterative solver, which enables it to effectively solve problems with large deformations, nonlinear material behavior, and contact. Sierra / SM has a versatile library of continuum and structural elements, and a large library of material models. The code is written for parallel computing environments enabling scalable solutions of extremely large problems for both implicit and explicit analyses. It is built on the SIERRA Framework, which facilitates coupling with other SIERRA mechanics codes . This document describes the functionality and input syntax for Sierra/SM.

97 MATHEMATICS AND COMPUTING↗

Sierra/SolidMechanics 4.58 User's Guide

Sierra / SolidMechanics (Sierra / SM) is a Lagrangian, three-dimensional code for finite element analysis of solids and structures. It provides capabilities for explicit dynamic, implicit quasistatic and dynamic analyses. The explicit dynamics capabilities allow for the efficient and robust solution of models with extensive contact subjected to large, suddenly applied loads. For implicit problems, Sierra / SM uses a multi-level iterative solver, which enables it to effectively solve problems with large deformations, nonlinear material behavior, and contact. Sierra / SM has a versatile library of continuum and structural elements, an d a large library of material models. The code is written for parallel computing environments enabling scalable solutions of extremely large problems for both implicit and explicit analyses. It is built on the SIERRA Framework, which facilitates coupling with other SIERRA mechanics codes . This document describes the functionality and input syntax for Sierra / SM.

36 MATERIALS SCIENCE↗

Development of phenomena identification and ranking table for Westinghouse lead fast reactor’s safety

The Westinghouse Lead-cooled Fast Reactor (LFR) is a medium-size, passively safe, economic, Gen-IV nuclear reactor. An important effort within the Westinghouse LFR program is the development of the safety analysis methodology, which comprises computer code development, model development, and experimental testing. A key initial task associated with the development of the safety analysis methodology is the identification of processes and phenomena that affect the plant's capability to meet selected safety performance indicators. This is accomplished through the development of a Phenomena Identification and Ranking Table (PIRT) for selected accident scenarios which, for this specific PIRT effort, included selected postulated design basis accidents and hypothetical beyond design basis accidents in LFRs. This paper describes the role of PIRT in the development of the Westinghouse LFR safety analysis methodology and the process used in the PIRT development. Specifically, the Westinghouse LFR PIRT assessed importance of pertinent phenomena and identified gaps in their knowledge-base by evaluating current modeling capabilities and data available for validation. The ultimate goal was to provide guidance on computer code development and validation efforts and to prioritize testing to support LFR design and licensing. The key phenomena and processes that are deemed highly important for the safety performance indicators, but for which the state of knowledge is low, are presented. The testing program and analyses development are planned to address significant phenomena in the PIRT.

Lead fast reactor↗

System code evaluation of near-term accident tolerant claddings during pressurized water reactor station blackout accidents

Following the Fukushima Daiichi nuclear accident in 2011, researches on Accident-tolerant fuels (ATFs) are currently of high interest in not only the nuclear industry but also governmental and international organizations. In this work, a quantitative evaluation of the performance of monolithic FeCrAl cladding and Cr-coated Zircaloy cladding has been performed for Pressurized Water Reactor (PWR) Station Blackout (SBO) accidents. A generic PWR model has been built in system thermal-hydraulics code TRACE based on the Surry Nuclear Power Station with counter-current natural circulation modelling capability for hotleg and steam generator U-tube components during the accidents. The base model results are then compared to MELCOR and RELAP simulations to verify the system component implementation in TRACE. Two PWR SBO scenarios were investigated, including: short-term SBO and long-term SBO with early reactor coolant pump (RCP) seal failure. These scenarios are defined to be very similar to the accidents studied in the State-of-the-Art Reactor Consequence Analysis (SOARCA) project. TRACE code is modified to reflect the oxidation kinetics of FeCrAl and Cr-coating. Larson-Miller creep rupture model is also implemented in TARCE using its built-in control systems to simulate the creep rupture of hotlegs. Additionally, the comparison between the TRACE models with and without the counter-current flow modeling resulted in significant difference when comparing ATF cladding to Zircaloy for short term SBO, while it marginal impacted the performance during long term SBO with RCP seal failure. For short term SBO, both ATF cladding underwent hot leg creep rupture ~20 min after Zircaloy cladding. While Zircaloy and Cr-coated cladding had generated significant amount of hydrogen gas (>10 kg) before the creep rupture event, FeCrAl cladding had only generated <0.5 kg of hydrogen gas. For long term SBO with RCP seal failure, significant hydrogen generation and fuel melting was predicted before hot leg creep rupture for the ATF cladding while providing only 10–20 min additional coping time compared to Zircaloy.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A Scoping Review of Mixed Initiative Visual Analytics in the Automation Renaissance

Artificial agents are increasingly integrated into data analysis workflows, carrying out tasks that were primarily done by humans. Our research explores how the introduction of automation recalibrates the dynamic between humans and automating technology. To explore this question, we conducted a scoping review encompassing twenty years of mixed-initiative visual analytic systems. To describe and contrast the relationship between humans and automation, we developed an integrated taxonomy to delineate the objectives of these mixed-initiative visual analytics tools, how much automation they support, and the assumed roles of humans. Here, we describe our qualitative approach of integrating existing theoretical frameworks with new codes we developed. Our analysis shows that the visualization research literature lacks consensus on the definition of mixed-initiative systems and explores a limited potential of the collaborative interaction landscape between people and automation. Our research provides a scaffold to advance the discussion of human-AI collaboration during visual data analysis. Our integrated taxonomy is available in the form of a web application on https://smonadjemi.github.io/miva.

Monadjemi, Shayan [ORNL] (ORCID:0000000293855969)↗

Code Coverage Status of ARC Code-DIF3D

The Argonne Reactor Code (ARC) software system supports users in their fast reactor design goals by providing neutronic, thermal-hydraulic, and structural analysis capabilities. DIF3D plays a pivotal role in the ARC system as the primary homogenized assembly neutronic calculation methodology for fast reactor problems. Over its 40 years history, ARC software usage with DIF3D has been applied to numerous fast and thermal spectrum reactor analysis projects with good to excellent comparison against experiments. With continued improvement of computation resources, many of the geometry modeling capabilities in DIF3D that were primarily used in low order schemes are not really needed anymore. Today, the diffusion and transport capabilities of DIF3D-VARIANT are primarily used in the reactor design process with some scattered usage of DIF3D-FD and DIF3D-Nodal. In recent work, the DIF3D software verification was completed for DIF3D-FD and DIF3D-VARIANT on the geometry options used in the Versatile Test Reactor project. While we can be confident that these capabilities of DIF3D are well used and thus trusted, it does not demonstrate that all possible input options of DIF3D are actually working, but just those that were tested as part of VTR are and that they are correct. Thus, the purpose of the present work is to identify a set of test problems for DIF3D and assess the code coverage of DIF3D for those test problems. The goal is to document what parts of the existing DIF3D code are touched by the set of test problems and which are not. Because the verification work done on DIF3D-VARIANT and DIF3D-FD was focused on the most common uses of DIF3D for fast reactor analysis, the code coverage assessment of those capabilities is the highest priority. This will ensure that nothing is being missed by the existing verification test problems that DIF3D relies upon. The DIF3D-Nodal capability will also be inspected for code coverage as part of this work to further ensure that regular regression testing of DIF3D will trap any likely errors the end user might experience with the DIF3D software. The code coverage analysis of DIF3D was performed with the Code Coverage Tool of the Intel Fortran compiler which requires modifications to the compilation of DIF3D. The detailed coverage tables are given for each submodule of DIF3D separately, and for the submodules which are primarily developed for DIF3D, most of the source files could be at least partially touched. Most of the uncovered parts/files could be easily ignored, because they are either for error message and debugging output or obviously not needed by DIF3D. Out of the entire source codes of DIF3D, only a few uncovered modules deserve further investigation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Code Coverage Status of the ARC Code DIF3D

The Argonne Reactor Code (ARC) software system supports users in their fast reactor design goals by providing neutronic, thermal-hydraulic, and structural analysis capabilities. DIF3D plays a pivotal role in the ARC system as the primary homogenized assembly neutronic calculation methodology for fast reactor problems. Over its 40 years history, ARC software usage with DIF3D has been applied to numerous fast and thermal spectrum reactor analysis projects with good to excellent comparison against experiments. With continued improvement of computation resources, many of the geometry modeling capabilities in DIF3D that were primarily used in low order schemes are not really needed anymore. Today, the diffusion and transport capabilities of DIF3D-VARIANT are primarily used in the reactor design process with some scattered usage of DIF3D-FD and DIF3D-Nodal. In recent work, the DIF3D software verification was completed for DIF3D-FD and DIF3D-VARIANT on the geometry options used in the Versatile Test Reactor project. While we can be confident that these capabilities of DIF3D are well used and thus trusted, it does not demonstrate that all possible input options of DIF3D are actually working, but just those that were tested as part of VTR are and that they are correct. Thus, the purpose of the present work is to identify a set of test problems for DIF3D and assess the code coverage of DIF3D for those test problems. The goal is to document what parts of the existing DIF3D code are touched by the set of test problems and which are not. Because the verification work done on DIF3D-VARIANT and DIF3D-FD was focused on the most common uses of DIF3D for fast reactor analysis, the code coverage assessment of those capabilities is the highest priority. This will ensure that nothing is being missed by the existing verification test problems that DIF3D relies upon. The DIF3D-Nodal capability will also be inspected for code coverage as part of this work to further ensure that regular regression testing of DIF3D will trap any likely errors the end user might experience with the DIF3D software. The code coverage analysis of DIF3D was performed with the Code Coverage Tool of the Intel Fortran compiler which requires modifications to the compilation of DIF3D. The detailed coverage tables are given for each submodule of DIF3D separately, and for the submodules which are primarily developed for DIF3D, most of the source files could be at least partially touched. Most of the uncovered parts/files could be easily ignored, because they are either for error message and debugging output or obviously not needed by DIF3D. Out of the entire source codes of DIF3D, only a few uncovered modules deserve further investigation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Coupling of Laminar-Turbulent Transition with RANS Computational Fluid Dynamics

This project combines several new concepts to create a boundary layer transition prediction capability that is suitable for analyzing modern hypersonic flight vehicles. The first new concept is the use of ''optimization'' methods to detect the hydrodynamic instabilities that cause boundary layer transition; the use of this method removes the need for many limiting assumptions of other methods and enables quantification of the interactions between boundary layer instabilities and the flow field imperfections that generate them. The second new concept is the execution of transition analysis within a conventional hypersonics CFD code, using the same mesh and numerical schemes for the transition analysis and the laminar flow simulation. This feature enables rapid execution of transition analysis with less user oversight required and no interpolation steps needed.

42 ENGINEERING↗

Leveraging NREL's ResStock & ComStock Dataset to Evaluate Building Stock Electrification: Preprint

Residential and commercial buildings accounted for 40% of U.S. energy consumption in 2022 and represent a significant opportunity for decarbonization through energy efficiency and electrification, and for grid planning. Building stock energy modeling is a powerful tool that can evaluate what-if scenarios as utilities, municipalities, policymakers, building owners and others work towards equitable building decarbonization and climate goals. This presentation will highlight several high-impact use cases of the National Renewable Energy Laboratory (NREL)'s highly granular, bottom-up building stock energy modeling tools, ResStock and ComStock. These use cases cover a wide range of project scale, from neighborhood electrification analysis and municipality long-term energy planning, to state energy code development and national policy evaluation. This presentation will showcase specific real-world applications for which ResStock and ComStock have been utilized across the country, including California codes and standards cost-effectiveness analysis, New York City affordable housing electrification cost gap analysis, and California targeted electrification and gas decommissioning analysis. For each use case, this presentation will illustrate how ResStock and ComStock played a crucial role in accurately characterizing regional building stocks, providing discrete and aggregated end-use load shapes, and calculating lifecycle consumption, emissions, and costs for a variety of building electrification strategies and scenarios. Finally, this presentation will demonstrate how the data provided by ResStock and ComStock can help unlock significant outcomes for these use cases, including but not limited to, customer bill impact, incentive and program design, and energy equity analyses.

building stock modeling↗

Upgrade of the Materials Analysis Particle Probe (MAPP-U) to decipher the impact of lithium-based surfaces on NSTX-U plasma behavior

Understanding the plasma wall interaction (PWI) remains a critical issue for the feasibility of thermonuclear magnetic fusion energy solutions. Key issues with PWI mechanisms in fusion tokamak reactors include: evolution of surface chemistry and its role on hydrogen retention. In particular how low-Z coatings such as lithium can impact the behavior of plasma at the edge and in the core. PMI (plasma-material interactions) are particularly important for strategies that involve low-recycling regimes and the use of lithium PFS (plasma-facing surfaces) to attain them, as in the case of NSTX-U. Recent reports have indicated the importance of access to the evolving plasma-facing surface during and in-between plasma discharges. Changes in surface chemistry and morphology due to ion bombardment and the difficulty of diagnosing plasma-facing surfaces, especially reactive surfaces, complicate the development of a predictive understanding of the wall and its interaction with the plasma. Consequently, this impairs the ability to design advanced PFC materials for future plasma-burning fusion reactors and appropriate PMI code validation. The Materials Analysis Particle Probe (MAPP) is an established and on-going PMI probe diagnostic system compatible with the highly chemically reactive system of lithium and boron coatings adopted by the NSTX-U research program. MAPP is the first PMI diagnostic to capture the surface physics and chemistry in-vacuo in a fusion tokamak system and correlate this data to controlled plasma shots. Currently MAPP captures this information at a fixed radial location at the NSTX-U outboard divertor region. The MAPP diagnostic has enabled understanding of the near-surface and surface chemistry of complex evolving lithiated and borated carbon-based PFC surfaces retention and transport of hydrogen. Coupled to atomistic simulations in collaboration with P. Krstic of Stony Brook U. MAPP has been very successful in achieving high-impact scientific research in its current grant period evidenced by two invited review articles and over 20 peer-reviewed manuscripts and over 40 contributed and invited presentations at both national and international conferences.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Efficient analysis routines for single and double peaked Type 2 AGN spectra

Driven by the imminent need to rapidly process and classify millions of AGN spectra drawn from next generation astronomical facilities, we present a spectral fitting routine for Type 2 AGN spectra optimized for high volume processing, using the quasar spectral fitting library (qsfit). We analyse an optically selected sample of 813 luminous Type 2 AGN spectra at z < 0.83 from the Sloan Digital Sky Survey (SDSS) to qualify its performance. We report a median narrow line H α/H β Balmer decrement of 4.5 ± 0.8, alluding to the presence of dust in the narrow line region (NLR). We publish a specialized QSFIT fitting routine for high signal-to-noise ratio spectra and general fitting routine for double peaked Type 2 AGN spectra applied on a subsample of 45 spectra from our parent sample. We report a median red and blue peak velocity separation of 390 ± 60 kms −1 . No trend is found for red or blue peaks to exhibit systematically different luminosity or ionization properties. Emission line diagnostics show that the double peaks in all sources are illuminated by an AGN-powered ionizing continuum. Finally, we examine the morphology of host galaxies of our double peaked sample. We find double peaked Type 2 AGN reside in merging systems at a comparable frequency to single peaked AGN. This suggests that the double peaked AGN phenomenon is likely to have a bi-conical outflow origin in the majority of cases. We publicly release the code used for spectral analysis and produced catalogues used in this work.

79 ASTRONOMY AND ASTROPHYSICS↗

Liquid hydrogen storage system for heavy duty trucks: Configuration, performance, cost, and safety

In this work, we investigate the potential of liquid hydrogen storage (LH 2 ) on-board Class-8 heavy duty trucks to resolve many of the range, weight, volume, refueling time and cost issues associated with 350 or 700-bar compressed H 2 storage in Type-3 or Type-4 composite tanks. We present and discuss conceptual storage system configurations capable of supplying H 2 to fuel cells at 5-bar with or without on-board LH 2 pumps. Structural aspects of storing LH 2 in double walled, vacuum insulated, and low-pressure Type-1 tanks are investigated. Structural materials and insulation methods are discussed for service at cryogenic temperatures and mitigation of heat leak to prevent LH 2 boiloff. Failure modes of the liner and shell are identified and analyzed using the regulatory codes and detailed finite element (FE) methods. The conceptual systems are subjected to a Failure modes and effects analysis (FMEA) and a safety, codes, and standards (SCS) review to rank failures and identify safety gaps. The results indicate that the conceptual systems can reach 19.6% usable gravimetric capacity, 40.9 g-H 2 /L usable volumetric capacity and $174-183/kg-H 2 cost (2016 USD) when manufactured 100,000 systems annually.

08 HYDROGEN↗

Analysis of Hydrology and Interim Measure Performance in the Chromium Plume at Los Alamos National Laboratory - 20498

Hexavalent chromium, Cr(VI), is present in the regional aquifer beneath Los Alamos National Laboratory (LANL) at concentrations greater than the New Mexico groundwater standard (50 ppb). The source of the Cr is blowdown from a power plant where potassium dichromate was used as a corrosion inhibitor from 1956 to 1972. Blowdown comingled with effluent from another source and migrated approximately 3 miles down a canyon, primarily as surface flow, until reaching an alluvial infiltration zone. The infiltrating water moved vertically through a thick (approximately 900 ft), complex stratigraphic sequence of unsaturated tuffs, basalt, and unconsolidated alluvial deposits. Horizontal transport also took place within spatially limited perched groundwater horizons in the vadose zone (VZ). Breakthrough at the regional aquifer water table is thought to have occurred in multiple hydraulic windows or 'drip points' located approximately 3 miles from the initial release site at the power plant. The regional aquifer is highly heterogeneous and anisotropic, and displays unconfined behavior in an upper zone and confined behavior at depth. The Cr(VI) contamination has so far remained confined to the upper approximately 25 m of the aquifer. The ambient hydraulic gradient in the regional aquifer is toward the east/southeast, but the water table is extremely flat in the Cr plume area compared to surrounding areas. There is some local variability in flow direction, potentially related to effects of seasonal pumping from nearby Los Alamos county water-supply wells, which draw from the regional aquifer several hundred feet below the water table. An interim measure (IM) is currently operating at the site to achieve hydraulic control of plume migration using a series of extraction wells, ion-exchange treatment, and injection along the downgradient portion of the plume. The migration of Cr contamination in the aquifer has been previously evaluated using a calibrated probabilistic numerical model, which has also been used to design the spatial configuration of extraction and injection wells for the IM, and to identify locations for monitoring wells. The model is built using the Finite Element Heat and Mass (FEHM) transfer code, and the Model Analysis and Decision Support (MADS) software is used for analyses. To date, the portion of the IM system that has operated in a continuous manner has been very effective at hydraulic control. Modifications to the flow field associated with the influence of IM pumping and injection are evaluated here, with modeling constraints to be provided by recent observations of tracers injected into both monitoring wells and IM injection wells, as well as opportunistic observations of geochemical signatures of treated ion-exchange effluent injected into IM injection wells. The tracer and geochemical data provide insights into modified flow directions within the Cr plume area that cannot be otherwise obtained. Results of the modeling work are being used to inform adaptive management of IM operations and will be a valuable tool for development and refinement of long-term remediation strategies. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

An independent analysis of bias sources and variability in wind plant pre-construction energy yield estimation methods

The wind resource assessment community has long had the goal of reducing the bias between wind plant pre-construction energy yield assessment (EYA) and the observed annual energy production (AEP). This comparison is typically made between the 50% probability of exceedance (P50) value of the EYA and the long-term corrected operational AEP (hereafter OA P50), and is known as the P50 bias. The industry has critically lacked an independent analysis of bias reduction investigated across multiple consultants to identify the greatest sources of uncertainty and variance in the EYA process and the best opportunities for uncertainty reduction. The present study addresses this gap by benchmarking consultant methodologies against each other and against operational data at a scale not seen before in industry collaborations. We consider data from 10 wind plants and evaluate discrepancies between eight consultancies in the steps taken from estimates of gross to net energy. Consultants tend to overestimate the gross energy produced at the turbines and then compensate by further overestimating downstream losses, leading to a mean P50 bias near zero, still with significant variability among the individual wind plants. Within our data sample, we find that consultant estimates of all loss categories, except environmental losses, tend to reduce the project-to-project variability of the P50 bias. The disagreement between consultants, however, remains flat throughout the addition of losses. Finally, we find that differences in consultants’ estimates of project performance can lead to differences up to $10/MWh in the levelized cost of energy for a wind plant.

Todd, Austin C.↗

Dataset for 'Ombadi, M. & Varadharajan, C. (2022). Urbanization and aridity mediate distinct salinity response to floods in rivers and streams across the Contiguous United States, Water Research'

This package contains data sets and code used to obtain the results in Ombadi, M., & Varadharajan, C. (2022). Urbanization and aridity mediate distinct salinity response to floods in rivers and streams across the Contiguous United States. Water Research, 118664. The folder "data" contains 259 .csv files, each of which has daily time series of concurrent streamflow (Q) and specific conductance (SC) for each of the sites used in this study originally downloaded from the USGS National Water Information System (NWIS; USGS, 2016). The number of data points in each of the files is at least 3650 (i.e. 10 years of daily measurements). The folder "RF_single_data" contains 259 .csv files, each of which include data used to train and test the Random Forest models at individual sites for predicting SC during days of floods. The folder "RF_regional_data" contains 3 .csv files, each of which include scaled data compiled from all sites within each climate zone (arid, temperate and wet). "metadata.csv" contains the physical properties of the 259 catchments corresponding to the sites used in this study; this data was extracted from GAGES-II dataset (Falcone et al., 2010). "RF_implementation.ipynb" is a Jupyter notebook with the code needed to implement the analysis using Random Forest models either for individual sites or for the regional models (for each climate zone). The code utilizes the data in the two folders: "RF_single_data" and "RF_regional_data" and the metadata.csv file.

54 ENVIRONMENTAL SCIENCES↗

Data for “Tree root nutrient uptake kinetics vary with nutrient availability, environmental conditions, and root traits: A global analysis”

This data package contains data and code used in the paper “Tree root nutrient uptake kinetics vary with nutrient availability, environmental conditions, and root traits: A global analysis”. The central product is a global dataset of root inorganic nutrient uptake rates and kinetics parameters covering temperate, boreal, and sub/tropical tree species, representing a collection of nutrient uptake data from published studies. This dataset enables tree investigation of root nutrient uptake rates across species, space, and experimental conditions. The data can also be combined with supplementary data on root and soil traits or with external datasets (e.g. R scripts contained within use data from FRED 3.0; (Iversen et al., 2021)). Contained within is the main nutrient data “uptake_data.csv” as well as 4 additional .csv files that link uptake data to supplementary measurements, source references, taxonomic information, and additional nutrient uptake measurements across nutrient gradients, and 1 .csv file that records meta-analysis results for plotting with the R scripts. There are seven R scripts that support data analysis and creation of the figures in the related publication.

54 ENVIRONMENTAL SCIENCES↗

Experimental and numerical studies of seismic fluid-structure interaction in a base-supported cylindrical vessel

Seismic design and qualification of advanced reactors will rely heavily on the use of verified and validated numerical models capable of capturing the interaction of the vessel, its contained fluid, and the internal equipment: fluid-structure interaction (FSI) analysis. Analytical solutions can be used for preliminary sizing and design of such vessels but their application is limited to simple geometries and boundary conditions, and small amplitude, translational (and rotational) inputs. To validate numerical models for seismic FSI analysis in finite element codes, a comprehensive set of experiments was performed on a liquid-filled cylindrical vessel, using a 6 degree-of-freedom earthquake simulator. Results in terms of sloshing frequency, damping ratio in sloshing modes, and hydrodynamic responses (wave height, hydrodynamic pressure, base shear, and base moment) for multi-directional earthquake simulator inputs are reported and compared with analytical solutions for liquid-filled vessels. The impact of seismic (base) isolation on hydrodynamic responses was studied using earthquake simulator inputs generated using a virtual isolation system. Data from the experiments are used to validate a numerical model of the fluid-filled vessel using the Arbitrary Lagrangian Eulerian (ALE) solver in the commercial finite element program LS-DYNA. Validation studies are presented for multi-directional seismic inputs, including rocking motions. Lagrangian modeling approaches using an elastic material formulation for the fluid are also investigated and their limitations and possible applications are identified. Here, the results are broadly applicable to the seismic response of base supported, liquid-filled vessels.

58 GEOSCIENCES↗