Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Analytical sensitivity 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 271 records · Page 15

Plant-wide modeling and techno-economic analysis of a direct non-oxidative methane dehydroaromatization process via conventional and microwave-assisted catalysis

Direct non-oxidative methane dehydroaromatization (DHA) process via conventional and microwave (MW)-assisted thermo-catalytic catalysis is studied. Rate models for methane DHA reactions, including the effect of catalyst deactivation, are developed by using the in-house experimental data. Model results for gas concentration profile and catalyst deactivation are in good agreement with the experimental data. This rate model is then used for the development of dynamic multi-scale, multi-physics commercial-scale reactor models. Total number of fixed bed reactors desired for a cyclic steady state process is estimated. Plant-wide models are then developed for conventional and MW-assisted processes for producing products of desired specifications. Techno-economic analysis of the methane DHA process is undertaken. Economics of these methane DHA processes are compared with the typical multi-step natural gas to aromatics production process via methanol synthesis. Sensitivity of internal rate of return (IRR) and net present value (NPV) to various economic and process parameters such as plant scale, desired rate of return, reactor cost, feedstock and utility cost, catalyst variable cost, and MW reactor cost is studied. Here, electric equivalent efficiency of the conventional methane DHA process is found to be 69.2 % and 67.3 % at 750 °C and 800 °C, respectively, while the MW-assisted methane DHA process has the electric equivalent efficiency of 48.9 % at 800 °C. IRRs of the conventional methane DHA process at 750 °C and 800 °C, and MW-assisted process are 15.2 %, 17.5 %, and 18.8 %, respectively for a methane feed flowrate of 19,782 kg/h, while the IRR of the multi-step natural gas to aromatics production process is estimated to be 0 % for the same plant scale. Impact of change in the methane price, electricity price, and catalyst cost is found to be considerable on the process economics, while the cost of the MW reactor is found to have negligible impact.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Delineating the Effects of Counterions on the Structural and Vibrational Properties of U(IV) Lindqvist Polyoxometalate Complexes

Herein we conducted a full investigation into the fundamental structural and vibrational properties of uranium(IV) Peacock−Weakley-type lacunary Lindqvist (W 10 ) polyoxometalate (POM) complexes. We recently demonstrated the importance of the secondary lattice elements in tuning the distortion of the D 4d symmetry in W 10 POM complexes, and here, we synthesized eight UW 10 complexes with different alkali metal counterions and evaluated how the composition and packing of counterion species affected complex structural and vibrational properties. Single-crystal X-ray diffraction analysis on complexes 1−8 revealed changes in structural distortion parameters as a function of differences in counterion configurations, while far-infrared and Raman spectra for 1−8 also demonstrated that vibrational mode frequencies were sensitive to changes in counterion composition and packing. To more effectively compare different counterion configurations, we developed counterion effective ionic radius (eIR) as a new structural parameter, and comparisons between structural distortion parameters and eIR values strongly suggested that modulation by the secondary lattice elements can affect structural and vibrational manifolds within POM complexes. Partial least squares (PLS) analysis was used to quantitatively evaluate correlations observed within this investigation, and PLS statistical models showed a strong correlation between counterion eIR and both structural distortion parameters and vibrational mode frequencies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Development of Direct Injection/Ionization Mass Spectrometry Methods for Whole Molecule Characterization

The objective of this work is to adapt ambient ionization mass spectrometry (AMS) techniques for the rapid analysis of intact uranium complexes, stable strontium, cerium, and explosive compounds. The methods used were “soft ionization” techniques, which facilitate the detection of whole molecule complexes. The soft ionization mass spectrometry (MS) techniques that were investigated include paper spray ionization (PSI), matrix-assisted ionization (MAI), electrospray ionization (ESI) and direct analysis in real time (DART). For the first time, PSI-MS methods were successfully developed for whole molecule uranium-containing analytes (uranyl acetate, uranyl nitrate, and uranyl-tributylphosphate complexes). This was also the first demonstration of uranium complex detection and characterization and one of the few examples of inorganic analysis using MAI techniques. Proof of concept experiments also putatively identified matrix-derived ions and ion complexes that have not previously been described in the literature. Additionally, PSI-MS on cotton swipe samples doped with a multi-element standard containing μg levels of U, Bi, Pb, Cd, Fe, and Zn were directly analyzed without purification, representing a major improvement over existing methods. Both PSI and MAI methods demonstrated limits of detection (LODs) in the 10 - 100’s ng for various uranyl species within a range of 10’s ppm - 100’s ppb, dependent on analytical method and analyte species. AMS methods were also developed for other inorganics, including Ce and Sr, and organic explosive residues to address specific challenges in environmental monitoring and forensics. Further refinement and qualification of the AMS techniques developed within this effort would lead to significant cost reduction and timeliness by facilitating the triage and queueing of samples for subsequent more sensitive and time-consuming analyses.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Advancing Enhanced Raman Spectroscopy to Detect and Measure NOW Pheromones (CPRB 2022 proposal)

Insect sex pheromones are chemical compounds that insects release to attract their partners over distances of hundreds of meters or even kilometers, in complete darkness and without any audible signals. Use of synthetic forms of key compounds have in some cases become an essential component of monitoring and/or managing key pests of agricultural crops, including navel orangeworm (Amyelois transitella) (NOW) in California tree nuts. These pheromone-based strategies can include monitoring, mating disruption, mass trapping, attract-and-kill and push-pull. There are currently multiple commercially available mating disruption products available for NOW and recent studies have demonstrated that they can be effectively used to reduce crop damage by this pest. Just how mating disruption works is not fully established and likely varies across products and target species. For instance, the extent to which synthetic pheromones compete with natural pheromone is not well understood, or in the case of monitoring, how efficaciously the insect follows the diffusing plumes, especially across large blocks and at plot borders. Furthermore, there may be specific conditions under which poor or impeded diffusion of synthetic pheromone diminishes the disruption effect, which could result in some males effectively locating females for reproduction. At the same time, many commercially available synthetic lures are also available, and growers have been effectively using them to track population development – although the attractive range of these lures remains unclear. That is, while pheromone lures can attract many moths, the relationship between trap capture and local populations, much less crop damage, remains unclear. Regardless of the emission source (aerosol puffers, meso emitters, lures etc.), we currently lack the ability to fully understand how these synthetic pheromone compounds diffuse away from their point-source of emission – and subsequently how this might affect the efficacy of mating disruption and/or the accuracy of monitoring efforts. Raman spectroscopy (RS) may offer a means to detect synthetic pheromone because of its analytical properties and recent technological advancements. In combination with nanostructured probes, RS has a demonstrated capability to detect volatiles at extremely low concentrations. While other methods for detecting pheromones have proven to be successful, such analytical methods like mass spectrometry (MS) gas chromatography (GC) require sampling in the fields and are not real-time. Other methods such as EAG/ EAD are also invasive (use antennas) but are of course valuable. Opto-electronic, micro and nano sensors just recently enticed interest. We fall in this latter category, proposing the use of optical spectroscopy (Raman) in a fiber aided by nanostructures to provide a chemical fingerprint without transductions and in real-time (i.e. no sampling). To date, we have been able to generate Raman signatures for the main chemicals in synthetic pheromones (a first for science, to the best of our knowledge), observe their trends at different concentrations, and also provide a clear distinction of such signatures in mixed samples via Principal Component Analysis (PCA) – all of which strengthens our conclusions on the ability of RS to detect emissions in more realistic conditions. We also established processes to enhance the detection sensitivity via preconcentration with Solid Phase Micro Extraction (SPME) fibers and Surface Enhanced Raman Spectroscopy (SERS) to increase the overall optical signals. Furthermore, we explored limitations and requirements for detecting from dispensers, lures, or live moths, including Dept. of Energy and CA State regulatory approvals. Finally, we setup some of our first “orchard in a box” experiments to measure pheromone diffusion along with developing appropriate models.

42 ENGINEERING↗

Study of the Segmental Dynamics and Ion Transport of Solid Polymer Electrolytes in the Semi-crystalline State

Solid polymer electrolytes are promising in fulfilling the requirements for a stable lithium metal anode toward higher energy and power densities. In this work, we investigate the segmental dynamics, ionic conductivity, and crystallinity of a polymer electrolyte consisting of poly(ethylene oxide) (PEO) and lithium triflate salt, in the semi-crystalline state. Using quasi-elastic neutron scattering, the segmental dynamics of PEO chains confined between the crystalline lamellae is quantified, using Cole–Cole analysis. We show that the structural relaxation time, τ 0 , of PEO equilibrated near room temperature is six-fold longer than the same sample that had just cooled down to room temperature. This corresponds to a three-fold smaller ionic conductivity in the equilibrated condition. This work reveals that the segmental dynamics of semi-crystalline polymer electrolytes is very sensitive to thermal history. We demonstrate that quasi-elastic neutron scattering can be used to characterize the ion transport and segmental dynamics in the semi-crystalline state.

36 MATERIALS SCIENCE↗

Thin Film Electrodes for Anodic Stripping Voltammetry: A Mini-Review

Anodic stripping voltammetry (ASV) is a powerful electrochemical analytical technique that allows for the detection and quantification of a variety of metal ion species at very low concentrations in aqueous media. While early, traditional ASV measurements relied on macroscopic electrodes like Hg drop electrodes to provide surfaces suitable for plating/stripping, more recent work on the technique has replaced these electrodes with thin film metal electrodes generated in situ . Such electrodes are plated alongside the analyte species onto the surface of a primary electrode, producing a composite metal electrode from which the analyte(s) can then be stripped, identified, and quantified. In this minireview, we will explore the development and use of these unique electrodes in a variety of different applications. A number of metals (e.g., Hg, Bi, Sn, etc.) have shown promise as thin film ASV electrodes in both acidic and alkaline media, and frequently multiple metals in addition to the analyte of interest are deposited together to optimize the plating/stripping behavior, improving sensitivity. Due to the relatively simple nature of the measurement and its suitability for a wide range of pH, it has been used broadly: To measure toxic metals in the environment, characterize battery materials, and enable biological assays, among other applications. We will discuss these applications in greater detail, as well as provide perspective on future development and uses of these thin film electrodes in ASV measurements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Tutorial on Bayesian analysis of linear shock compression data

Gas gun and other shock compression experiments often produce shock wave velocity measurements that are linearly associated with particle velocity. Traditionally, this empirical relationship is quantified with a single Hugoniot curve that is estimated using least squares regression. However, for downstream modeling and simulation tasks, it is often more useful to have multiple Hugoniot curves in the pressure–volume plane that are consistent with the data. We employ Bayesian uncertainty quantification methods as a framework for propagating measurement uncertainty through to model parameters and predictions. Specifically, this Tutorial shows how to sample multiple Hugoniot curves in the pressure–volume plane that are consistent with the shock wave-particle velocity measurements in a two-step Bayesian approach. First, we obtain an analytical expression for the posterior distribution of the linear model parameters using Bayesian linear regression. Second, we propagate samples from the posterior distribution through the Rankine–Hugoniot equations to yield Hugoniot curves in the pressure–volume plane. The procedure is demonstrated with publicly available data on argon, copper, and nickel, and compared against bootstrapping and linear regression. The Bayesian procedure is shown to be interpretable, computationally inexpensive, and less sensitive than an alternative bootstrapping approach to the removal of the point in the copper dataset that has the largest particle velocity. As a Tutorial on Bayesian methodology for the shock compression community, we provide several derivations and explanations that make this paper self-contained, and make all code and data available at github.com/llnl/BALSCD.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The NanoSIMS-HR: The Next Generation of High Spatial Resolution Dynamic SIMS

The high lateral resolution and sensitivity of the NanoSIMS 50 and 50L series of dynamic SIMS instruments have enabled numerous scientific advances over the past 25 years. Here, in this study, we report on the NanoSIMS-HR, the first major upgrade to the series, and analytical tests in a suite of sample types, including an aluminum sample containing silicon crystals, microalgae, and plant roots colonized with a symbiotic fungus. Significant improvements have been made in the Cs + ion source, high voltage (HV) control, stage reproducibility, and other aspects of the instrument that affect performance. The modified design of the NanoSIMS-HR thermal-ionization Cs + source enables a 5 pA primary ion beam to be focused into a 100 nm spot, a ~2.5-fold increase compared to Cs + sources on previous instruments (~2 pA at 100 nm). The brightness of the new Cs + source enables an ultimate lateral resolution as high as 30 nm and improved detection limits for a given analysis area. Sample stage movement accuracy is higher than 500 nm, enabling many-fold higher throughput automated analyses. With the new HV control, the primary ion beam impact energy can be reduced from 16 to 2 keV, which enables higher depth resolution during depth profiling (a 2-fold improvement), albeit with a 5-fold decrease in lateral resolution. In the NanoSIMS-HR, the secondary ion column and detection system are identical to those used in the previous series, and the isotopic analysis performance is as precise as in previous NanoSIMS instruments.

54 ENVIRONMENTAL SCIENCES↗

In‐situ Analysis of Paste Properties in Resonant Acoustic Mixers for Quality Monitoring

Formulation control is key to achieving consistent target properties of energetic materials, as feedstock variations and slight deviations in the ratios of different ingredients can have major effects on final product properties, particularly in dense pastes with high particle loading >65 vol.%. In large‐scale operations, it is imperative to either correct or remove batches of material that perform outside baseline property specifications as early as possible to avoid unnecessary processing of suboptimal material. Quality monitoring is the practice of measuring material properties during processing using process analytical technologies as opposed to only testing the properties of the final product; it is a key principle in the quality‐by‐design frameworks used for designing formulations and manufacturing processes. Herein, a process analytical technology method for correlating material properties of dense pastes directly after mixing in a Resonant Acoustic Mixer to motor data is developed and used to detect differences in the particle content of dense paste formulations. This method was also capable of detecting variations in powder feedstock properties, such as particle packing efficiency, and is sensitive enough to detect changes of 2 wt.% in the total solids content of the formulation. The techniques presented herein show excellent promise for use as a process analytical technology capable of quantifying formulation effects on material movement modes during resonant acoustic mixing.

Materials science↗

Techno-Economic, Feasibility, and Life Cycle Analysis of Renewable Propane: 2025 Update

To clarify the current and future landscape for renewable propane (RP) production, this work evaluates the value proposition of recovering RP from existing and planned hydroprocessed esters and fatty acids (HEFA) biorefineries and surveys emerging technologies under development or deployment. HEFA biorefineries co-produce a propane-rich fuel gas stream, normally used to meet HEFA process heat requirements, from which propane can be recovered and sold to create an additional revenue stream alongside liquid transportation fuels such as renewable diesel (RD) and sustainable aviation fuel (SAF). This report updates and extends a 2022 analysis of RP recovery from HEFA facilities by escalating capital and operating costs to 2024 prices, incorporating recent policy developments (including the Section 45Z Clean Fuel Production Credit), evaluating RP recovery for both RD- and SAF-focused HEFA facilities at two scales (3,000 and 75,000 barrels per day of feedstock), and quantifying the impact of RP recovery on HEFA liquid-fuel carbon intensity (CI) and associated tax credits using the 45ZCF-GREET model. For a 3,000 BPD RD-focused HEFA facility, approximately 3.5 million gallons per year (MGPY) of RP can be recovered; in this base case, the estimated payback period is 18 months based on the total installed cost of the RP recovery equipment and 36 months based on the total capital investment for the entire RP recovery project. The payback period is slightly shorter for the analogous SAF-focused configuration (approximately 4.3 MGPY RP). Sensitivity analysis shows that CAPEX magnitude, RP recovery plant scale, and CI-driven tax credit valuations are the dominant determinants of project viability. RP recovery may increase the CI of HEFA liquid fuels, which can reduce liquid-fuel tax credits (a key revenue stream for the HEFA biorefinery) and lengthen payback periods. However, RP recovery generally remains economically favorable across a wide range of plausible scenarios and market conditions. The report also summarizes emerging pathways that could expand future RP supply.

09 BIOMASS FUELS↗

Time-resolved luminescence detection of peroxynitrite using a reactivity-based lanthanide probe

Peroxynitrite (ONOO - ) is a powerful and short-lived oxidant formed in vivo, which can react with most biomolecules directly. To fully understand the roles of ONOO - in cell biology, improved methods for the selective detection and real-time analysis of ONOO - are needed. We present a water-soluble, luminescent europium(III) probe for the rapid and sensitive detection of peroxynitrite in human serum, living cells and biological matrices. We have utilised the long luminescence lifetime of the probe to measure ONOO - in a time-resolved manner, effectively avoiding the influence of autofluorescence in biological samples. To demonstrate the utility of the Eu(III) probe, we monitored the production of ONOO - in different cell lines, following treatment with a cold atmospheric plasma device commonly used in the clinic for skin wound treatment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hydrogenolysis of Poly(Ethylene–co–Vinyl Alcohol) and Related Polymer Blends over Ruthenium Heterogeneous Catalysts

The hydrogenolysis of polymers is emerging as a promising approach to deconstruct plastic waste into valuable chemicals. Yet, the complexity of plastic waste, including multilayer packaging, is a significant barrier to handling realistic waste streams. Herein, we reveal fundamental insights into a new chemical route for transforming a previously unaddressed fraction of plastic waste – poly(ethylene–co­­–vinyl alcohol) (EVOH) and related polymer blends – into alkane products. Additionally, we report that Ru/ZrO2 is active for the concurrent hydrogenolysis, hydrogenation, and hydrodeoxygenation of EVOH and its thermal degradation products into alkanes (C1–C35) and water. Detailed reaction data, product analysis, and catalyst characterization reveal that the in–situ thermal degradation of EVOH forms aromatic intermediates that are detrimental to catalytic activity. Increased hydrogen pressure promotes hydrogenation of these aromatics, preventing catalyst deactivation and improving alkane product yields. Calculated apparent rates of C–C scission reveal that the hydrogenolysis of EVOH is slower than low–density polyethylene. We apply these findings to achieve hydrogenolysis of EVOH/polyethylene blends and elucidate the sensitivity of hydrogenolysis catalysts to such blends. Overall, we demonstrate progress towards efficient catalytic processes for the hydroconversion of waste multilayer film plastic packaging into valuable products.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A simple numerical method for evaluating heat dissipation from curved wires with periodic applied heating

In many situations, the dual-purpose heater/thermometer wires used in the three-omega method – one of the most precise and sensitive techniques for measuring the thermal conductivity of thin films and interfaces – must include bends and curves to avoid obstructions on the surface of a sample. Although the three-omega analysis assumes that the heating wire is infinitely long and straight, recent experimental work has demonstrated that in some cases curved-wire geometries can be used without introducing detectable systematic error. We describe a general numerical method that can be used to calculate the temperature of three-omega heating wires with arbitrary wire geometries. Furthermore, this method provides experimentalists with a simple quantitative procedure for calculating how large the systematic error caused by a particular wire asymmetry will be. We show calculations of two useful cases: a straight wire with a single bend of arbitrary angle and a wire that forms a circle. We find that the amplitude of the in-phase temperature oscillations near a wire that forms a circle differs from the prediction using the analytic straight-line source solution by <12%, provided that the thermal penetration depth is less than ten times the radius of curvature of the wire path. The in-phase temperature amplitude 1.5 wire widths away from a 90° bend in a wire is within 11% of the straight-line source prediction for all penetration depths greater than the wire width. Our calculations indicate that the straight-line source solution breaks down significantly when the wire bend angle is less than 45°.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Circuit-level design principles for transmission-mode microwave impedance microscopy

A recently developed technique of transmission-mode microwave impedance microscopy (T-MIM) has greatly extended the capabilities of standard reflection-mode MIM to novel applications, such as the in operando study of nanoscale electro-acoustic devices. As is common for new techniques, systematic design principles for boosting sensitivity and balancing bandwidth are lacking. Here, we show numerically and analytically that the T-MIM signal is proportional to the reflection-mode voltage enhancement factor η of the circuit, as long as the output impedance of the local voltage source is properly treated. We show that this proportionality holds in the currently achievable “weak sampling” regime and beyond, for which we demonstrate a realistic path with commercially available superconducting components and critically coupled impedance matching networks. Furthermore, we demonstrate that for these next-generation designs, the sensitivity is generally maximized at a slightly different frequency from the unloaded S 11 resonance, which can be explained by the maximum power transfer theorem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Nanoscale Compositional Mapping of Commercial LiNi 0.8 Co 0.15 Al 0.05 O 2 Cathodes Using Atom Probe Tomography

Nickel-rich cathodes provide improved specific capacity that lead to higher gravimetric energy density, critical for electric vehicles. However, poor long-term capacity retention at elevated temperatures/high C-rates (the rate of charge and discharge on a battery) stem from material issues: surface phase changes, corrosive side reactions with the electrolyte, ion dissolution and propagation of cracks. Introducing dopants, developing nanoscale surface coatings and graded core-shell structures have all improved the electrochemical performance of nickel-rich cathodes. However, a material level understanding on the effect of Li composition and distribution in Ni-rich cathodes is limited, due to a lack of characterization methods available that can directly image Li at the nanoscale. Hence, it is critical to establish methods such as atom probe tomography (APT), having both nanometer scale spatial resolution and high compositional sensitivity to quantitatively analyze battery cathodes. To fully realize its potential as a method for quantitative compositional analysis of commercial Li-ion batteries, we provide a comprehensive description of the challenges in sample preparation and analyze the dependency of the analysis parameters, specifically laser pulse energy on the measured stoichiometry of elements in a high Ni content cathode material; LiNi 0.8 Co 0.15 Al 0.05 O 2 (NCA). In this work, our findings show the stoichiometry variations cannot be explained by charge state ratios or Ga implantation damage alone during FIB preparation, indicating that additional factors such as crystallographic orientation may need to be considered to achieve quantitative nanoscale compositional analysis of such battery cathodes using APT.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Process Interactions Can Change Process Ranking in a Coupled Complex System Under Process Model and Parametric Uncertainty

For a complex hydrologic system with multiple processes and process interactions, global sensitivity analysis is often used to identify important or influential parameters for model development and improvement. The identification is complicated by process model uncertainty, when a system process can be represented by multiple process models. This study develops a new total-effect process sensitivity index to identify influential processes under model uncertainty. This is done by extending Sobol's total-effect parameter sensitivity index for one system model to total-effect process sensitivity index for multiple system models to account for uncertainty in process models and model parameters. The total-effect process sensitivity index includes not only the first-order process sensitivity index for measuring the importance of individual processes but also higher-order indices that account for process interactions. The total-effect process sensitivity index can identify an influential process that itself and its interactions with other processes influence a model output. Here, the total-effect process sensitivity index is applied to two numerical examples: (a) Sobol's G*-functions with analytical solutions of first-order and total-effect process sensitivity indices, and (b) groundwater flow models with interactions between recharge, geology, and snowmelt processes. The second evaluation shows that, due to second-order and higher-order process interactions, the first-order and total-effect process sensitivity indices give different process ranking. It is thus necessary to estimate both first-order and total-effect process sensitivity indices to appreciate the difference between the first-order impact of a process alone and the overall total-effect impact of the process itself and its interactions with other processes on a model output.

54 ENVIRONMENTAL SCIENCES↗

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↗

Machine Learning-enabled Scalable Performance Prediction of Scientific Codes

Hardware architectures become increasingly complex as the compute capabilities grow to exascale. Here, we present the Analytical Memory Model with Pipelines (AMMP) of the Performance Prediction Toolkit (PPT). PPT-AMMP takes high-level source code and hardware architecture parameters as input and predicts runtime of that code on the target hardware platform, which is defined in the input parameters. PPT-AMMP transforms the code to an (architecture-independent) intermediate representation, then (i) analyzes the basic block structure of the code, (ii) processes architecture-independent virtual memory access patterns that it uses to build memory reuse distance distribution models for each basic block, and (iii) runs detailed basic-block level simulations to determine hardware pipeline usage. PPT-AMMP uses machine learning and regression techniques to build the prediction models based on small instances of the input code, then integrates into a higher-order discrete-event simulation model of PPT running on Simian PDES engine. We validate PPT-AMMP on four standard computational physics benchmarks and present a use case of hardware parameter sensitivity analysis to identify bottleneck hardware resources on different code inputs. We further extend PPT-AMMP to predict the performance of a scientific application code, namely, the radiation transport mini-app SNAP. To this end, we analyze multi-variate regression models that accurately predict the reuse profiles and the basic block counts. We validate predicted SNAP runtimes against actual measured times.

97 MATHEMATICS AND COMPUTING↗