Search NASA⌕ Search

SEARCH · Search NASA

Results for “Modeling & Simulation (M&S)”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Measurement-based quantum simulation of Abelian lattice gauge theories

Numerical simulation of lattice gauge theories is an indispensable tool in high energy physics, and their quantum simulation is expected to become a major application of quantum computers in the future. In this work, for an Abelian lattice gauge theory in d d spacetime dimensions, we define an entangled resource state (generalized cluster state) that reflects the spacetime structure of the gauge theory. We show that sequential single-qubit measurements with the bases adapted according to the former measurement outcomes induce a deterministic Hamiltonian quantum simulation of the gauge theory on the boundary. Our construction includes the (2+1) ( 2 + 1 ) -dimensional Abelian lattice gauge theory simulated on three-dimensional cluster state as an example, and generalizes to the simulation of Wegner’s lattice models M_{(d,n)} M ( d , n ) that involve higher-form Abelian gauge fields. We demonstrate that the generalized cluster state has a symmetry-protected topological order with respect to generalized global symmetries that are related to the symmetries of the simulated gauge theories on the boundary.Our procedure can be generalized to the simulation of Kitaev’s Majorana chain on a fermionic resource state. We also study the imaginary-time quantum simulation with two-qubit measurements and post-selections, and a classical-quantum correspondence, where the statistical partition function of the model M_{(d,n)} M ( d , n ) is written as the overlap between the product of two-qubit measurement bases and the wave function of the generalized cluster state.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

ARM shortwave spectrometers to study the clear-cloud transition zone and mixing processes

The proposed research is a collaborative effort between NASA/Goddard Space Flight Center and the Hebrew University of Jerusalem. While the NASA team focused on analyzing ground-based hyper-spectral radiance observations to understand cloud edge properties and their connection to mixing processes, the Hebrew University team tackled the problem through cloud modeling activities. By approximating the shortwave spectra in the cloud-clear transition zone as a linear combination of purely clear and purely cloudy spectra we can characterize the variations of cloud optical thickness and cloud droplet effective radius in the transition zone. When applying this method to the measurements of a ground-based shortwave spectroradiometer at the ARM’s SGP site, representing continental conditions, and MAGIC field campaign between Log Angeles, California and Honolulu, Hawaii, representing maritime scenarios, we found that cloud optical depth consistently decreases in both cases, but droplet size decreases much more substantially for the continental regime, suggesting different mixing processes for the continental and maritime conditions. The investigation and measurements of radiation clouds were coupled with a unique cloud modeling. A novel spectral bin microphysics was developed and implemented to the System of Atmospheric Modeling (SAM). In order to resolve cloud transition zones with high spatial gradients of microphysical variables a unique high resolution (10 m) was used in simulations. The model calculates droplet size distributions in each grid point. The model output was transferred to the NASA/GSFC team for utilization in radiative calculations and testing of both radiative algorithm and model representation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Validation of the EDGES Low-band Antenna Beam Model

The response of the antenna is a source of uncertainty in measurements with the Experiment to Detect the Global Epoch of Reionization Signature (EDGES). We aim to validate the electromagnetic beam model of the low-band (50–100 MHz) dipole antenna with comparisons between models and against data. We find that simulations of a simplified model of the antenna over an infinite perfectly conducting ground plane are, with one exception, robust to changes in the numerical electromagnetic solver code or algorithm. For simulations of the antenna with the actual finite ground plane and realistic soil properties, we find that two out of three numerical solvers agree well. Applying our analysis pipeline to a simulated drift-scan observation from an early EDGES low-band instrument that had a 10 m × 10 m ground plane, we find residual levels after fitting and removing a five-term foreground model from the simulated data binned in local sidereal time (LST) average about 250 mK with ±40 mK variation between numerical solvers. A similar analysis of the primary 30 m × 30 m sawtooth ground plane reduced the LST-averaged residuals to about 90 mK with ±10 mK between the two viable solvers. More broadly we show that larger ground planes generally perform better than smaller ground planes. Simulated data have a power that is within 4% of real observations, a limitation of net accuracy of the sky and beam models. We observe that residual spectral structures after foreground model fits match qualitatively between simulated data and observations, suggesting that the frequency dependence of the beam is reasonably represented by the models. We find that a soil conductivity of 0.02 S m{sup -1} and relative permittivity of 3.5 yield good agreement between simulated spectra and observations. This is consistent with the soil properties reported by Sutinjo et al. for the Murchison Radio-astronomy Observatory, where EDGES is located.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Simulating High-Frequency Seismograms in Realistic Earth Models to Better Understand Source Discrimination Based on Differential Magnitudes ( M L− M c)

Discriminating low-yield underground nuclear explosions from small earthquakes is a key task in monitoring nuclear test ban treaties. P/S amplitude ratios have been an effective discriminant for moderate-sized events recorded at regional distances, but it is unclear if they are as effective in discriminating small seismic events recorded at local distances (<150 km). The difference between local magnitude (M L ) and coda duration magnitude (M c ) has been proposed as a new discriminant that may complement P/S amplitude ratios at local distances. Here, in this work, we calculate high-frequency (up to ∼4 Hz) synthetic seismograms at epicentral distances of 0–30 km in realistic models of the Salt Lake basin (Utah, United States) to better understand how variations in source type and depth affect M L −M c values. The Earth models incorporate simplified 1D and deterministic 3D structures, small-wavelength stochastic velocity perturbations, and surface topography. Coda waves are enhanced for the more complicated models compared to the base 1D model, but still underpredict observed durations by about a factor of two, which results in overprediction of amplitude to duration ratios (i.e., M L −M c values) for a near-surface explosion and a 7 km deep earthquake. For both source types, the predicted M L and M c values decrease as source depth increases, and M L −M c shows only minor variation with depth; however, M L −M c is on average ∼0.5 units smaller for explosions than earthquakes. This finding may imply that M L −M c has sensitivity to source type, in addition to being a depth discriminant, but more modeling is needed given the limitations of the current study. Future modeling should incorporate higher-frequency (≳5 Hz) simulations over a larger distance range (0–150 km), where M L and M c are commonly measured, while honoring low shear velocities (<300 m/s) near the surface and sampling a wider range of earthquake and explosion source mechanisms.

Hutchings, Sean J. [Univ. of Utah, Salt Lake City,↗
NASA NTRS1991

A nested grid, nonhydrostatic, elastic model using a terrain-following coordinate transformation - The radiative-nesting boundary conditions

This study presents a nested-grid nonhydrostatic and elastic model using a terrain-following coordinate transformation as well as a unique application of grid-nesting techniques to the time-splitting elastic model. A simulation of the 10-m-high Witch of Agnesi Mountain provides the control to test this new model. The results show that the model produces the same solution as that derived from a simple linear analytic model. It is demonstrated that the new nested-grid model improves model resolution without resorting to the costly method of placing a fine-resolution grid over the entire domain. Since the wave reflection from the boundaries of the fine-grid model is well controlled, the boundary of the nested fine-grid model can be located even at the wave-active region. The model can be used to simulate various weather systems in which scale interactions are important.

Chen, Chaing↗
DOE OSTI2024

Dynamic data-driven multiscale modeling for predicting the degradation of a 316L stainless steel nuclear cladding material

Here, we have developed a long short-term memory stacked ensemble (LSTM-SE) surrogate modeling approach that can provide rapid predictions of microstructural evolution and the resultant mechanical properties of American Iron and Steel Institute (AISI) 316L series stainless steel (316LSS) fuel cladding under conditions of varying temperature and radiation dose rate. To acquire training data, we developed and implemented a kinetic Monte Carlo (KMC) model to simulate precipitation kinetics of M 23 C 6 , γ', and G phases within SS316L cladding. Experimentally reported precipitation kinetics of SS316L in literature were linked to the kinetic parameters of the simulated precipitation in our KMC model. The model was then used to simulate microstructure evolution under synthetically generated treatments of varying temperature and radiation dose rate, for periods of up to 3000 hours. Changes in volume fraction, number density, and particle size of precipitates were recorded, and particle area fractions were correlated using statistical methods to develop the surrogate model. Simultaneously, the mechanical properties of the simulated microstructures were evaluated using microstructure-based finite element method (FEM) analysis to determine the elastic modulus, yield stress, ultimate tensile strength, and elongation to failure of the aged microstructures. Using this approach, our surrogate model can predict precipitation behavior within 0.25% volume fraction and mechanical properties within 6% relative error from the values predicted by the KMC and FEM models using 50 training simulations as input. The trained recurrent neural network-based model can return estimations of precipitation kinetics and mechanical properties ~1000 times faster than the physics-based codes. This work demonstrates, as a proof of concept, that reactor material service lifetimes under variable service conditions can be predicted for a statistics-based model from a practicably obtainable dataset.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Dynamic Data-Driven Multiscale Modeling for Predicting Microstructure and Mechanical Properties of 316L Stainless Steel

We have developed a long short-term memory stacked ensemble (LSTM-SE) surrogate modeling approach that can provide rapid predictions of microstructural evolution and the resultant mechanical properties of American Iron and Steel Institute (AISI) 316L series stainless steel (316LSS) fuel cladding under conditions of varying temperature and radiation dose rate. To acquire training data, we developed and implemented a kinetic Monte Carlo (KMC) model to simulate precipitation kinetics of M 23 C 6 , γ', and G phases within SS316L cladding. Experimentally reported precipitation kinetics of SS316L in literature were linked to the kinetic parameters of the simulated precipitation in our KMC model. The model was then used to simulate microstructure evolution under synthetically generated treatments of varying temperature and radiation dose rate, for periods of up to 3000 hours. Changes in volume fraction, number density, and particle size of precipitates were recorded, and particle area fractions were correlated using statistical methods to develop the surrogate model. Simultaneously, the mechanical properties of the simulated microstructures were evaluated using microstructure-based finite element method (FEM) analysis to determine the elastic modulus, yield stress, ultimate tensile strength, and elongation to failure of the aged microstructures. Using this approach, our surrogate model can predict precipitation behavior within 0.25% volume fraction and mechanical properties within 6% relative error from the values predicted by the KMC and FEM models using 50 training simulations as input. The trained recurrent neural network-based model can return estimations of precipitation kinetics and mechanical properties ~1000 times faster than the physics-based codes. This work demonstrates, as a proof of concept, that reactor material service lifetimes under variable service conditions can be predicted for a statistics-based model from a practicably obtainable dataset.

36 MATERIALS SCIENCE↗
NASA NTRS2012

Sensitivity Studies of Dust Ice Nuclei Effect on Cirrus Clouds with the Community Atmosphere Model CAM5

In this study the effect of dust aerosol on upper tropospheric cirrus clouds through heterogeneous ice nucleation is investigated in the Community Atmospheric Model version 5 (CAM5) with two ice nucleation parameterizations. Both parameterizations consider homogeneous and heterogeneous nucleation and the competition between the two mechanisms in cirrus clouds, but differ significantly in the number concentration of heterogeneous ice nuclei (IN) from dust. Heterogeneous nucleation on dust aerosol reduces the occurrence frequency of homogeneous nucleation and thus the ice crystal number concentration in the Northern Hemisphere (NH) cirrus clouds compared to simulations with pure homogeneous nucleation. Global and annual mean shortwave and longwave cloud forcing are reduced by up to 2.0+/-0.1Wm (sup−2) (1 uncertainty) and 2.4+/-0.1Wm (sup−2), respectively due to the presence of dust IN, with the net cloud forcing change of −0.40+/-0.20W m(sup−2). Comparison of model simulations with in situ aircraft data obtained in NH mid-latitudes suggests that homogeneous ice nucleation may play an important role in the ice nucleation at these regions with temperatures of 205-230 K. However, simulations overestimate observed ice crystal number concentrations in the tropical tropopause regions with temperatures of 190- 205 K, and overestimate the frequency of occurrence of high ice crystal number concentration (greater than 200 L(sup-1) and underestimate the frequency of low ice crystal number concentration (less than 30 L(sup-1) at NH mid-latitudes. These results highlight the importance of quantifying the number concentrations and properties of heterogeneous IN (including dust aerosol) in the upper troposphere from the global perspective.

Liu, Xiaohong↗
NASA NTRS2013

Multi-sensor Cloud Retrieval Simulator and Remote Sensing from Model Parameters : Synthetic Sensor Radiance Formulation - Pt. 1

In this paper we describe a general procedure for calculating synthetic sensor radiances from variable output from a global atmospheric forecast model. In order to take proper account of the discrepancies between model resolution and sensor footprint, the algorithm takes explicit account of the model subgrid variability, in particular its description of the probability density function of total water (vapor and cloud condensate.) The simulated sensor radiances are then substituted into an operational remote sensing algorithm processing chain to produce a variety of remote sensing products that would normally be produced from actual sensor output. This output can then be used for a wide variety of purposes such as model parameter verification, remote sensing algorithm validation, testing of new retrieval methods and future sensor studies.We show a specific implementation using the GEOS-5 model, the MODIS instrument and the MODIS Adaptive Processing System (MODAPS) Data Collection 5.1 operational remote sensing cloud algorithm processing chain (including the cloud mask, cloud top properties and cloud optical and microphysical properties products). We focus on clouds because they are very important to model development and improvement.

Simulations↗
NASA NTRSDate not supplied

Efficient Variable-Pitch Propeller Aerodynamic Model Development for Vectored-Thrust eVTOL Aircraft

This paper describes a variable-pitch propeller testing and aerodynamic modeling methodology for electric vertical takeoff and landing (eVTOL) aircraft using wind tunnel testing. Propellers used for eVTOL aircraft propulsion systems experience a wide range of operating conditions resulting in significant variation of axial thrust and torque, as well as off-axis forces and moments, across a typical flight envelope. An experimental design is developed to collect informative wind tunnel data in an efficient manner using response surface methods. System identification methods are then applied to isolated propeller wind tunnel data to develop a mathematical model of the propeller valid throughout the transition envelope of a vectored-thrust eVTOL aircraft. Multiple explanatory variable definitions are postulated and compared using modeling and prediction performance metrics. Modeling results validated against data withheld from the modeling process indicate good predictive capability and agree with theoretical expectations. All identified propeller models are provided, which allows the models to be used in future eVTOL aircraft flight dynamics simulations.

aerodynamic modeling↗
NASA NTRS1993

Analysis of the GOES 6.7 micrometer channel observations during FIRE 2

Clouds form in moist environments. FIRE Phase II Cirrus Implementation Plan (August, 1990) noted the need for mesoscale measurements of upper tropospheric water vapor content. These measurements are needed for initializing and verifying numerical weather prediction models and for describing the environment in which cirrus clouds develop and dissipate. Various instruments where deployed to measure the water vapor amounts of the upper troposphere during FIRE II (e.g. Raman lidar, CLASS sonds and new cryogenic frost hygrometer on-board aircraft). The formation, maintenance and dissipation of cirrus clouds involve the time variation of the water budget of the upper troposphere. The GOES 6.7 mu m radiance observations are sensitive to the upper tropospheric relative humidity, and therefore proved extremely valuable in planning aircraft missions during the field phase of FIRE II. Warm 6.7 mu m equivalent black body temperatures indicate a relatively dry upper troposphere and were associated with regions generally free of cirrus clouds. Regions that were colder, implying more moisture was available may or may not have had cirrus clouds present. Animation of a time sequence of 6.7 mu m images was particularly useful in planning various FIRE missions. The 6.7 mu m observations can also be very valuable in the verification of model simulations and describing the upper tropospheric synoptic conditions. A quantitative analysis of the 6.7 mu m measurement is required to successfully incorporate these satellite observations into describing the upper tropospheric water vapor budget. Recently, Soden and Bretherton (1993) have proposed a method of deriving an upper tropospheric humidity based on observations from the GOES 6.7 mu m observations. The method is summarized in the next section. In their paper they compare their retrieval method to radiance simulations. Observations were also compared to ECMWF model output to assess the model performance. The FIRE experiment provides a unique opportunity to further verify the GOES upper tropospheric relative humidity retrieval scheme by providing (1) aircraft observations to cross-validate the calibration of the GOES 6.7 mu m channel, (2) accurate upper tropospheric water vapor concentrations for verification, and (3) veritical variability of upper tropospheric water vapor.

Soden, B. J.↗
DOE OSTI2021

Stellar Mass Black Hole Formation and Multimessenger Signals from Three-dimensional Rotating Core-collapse Supernova Simulations

We present self-consistent 3D core-collapse supernova simulations of a 40 ${M}_{\odot }$ progenitor model using the isotropic diffusion source approximation for neutrino transport and an effective general relativistic potential up to `0.9 s postbounce. We consider three different rotational speeds with initial angular velocities of ${{\rm{\Omega }}}_{0}=0$, 0.5, and 1 rad s -1 and investigate the impact of rotation on shock dynamics, black hole (BH) formation, and gravitational wave (GW) signals. The rapidly rotating model undergoes an early explosion at ~250 ms postbounce and shows signs of the low $T/| W| $ instability. We do not find BH formation in this model within ~460 ms postbounce. In contrast, we find BH formation at 776 ms postbounce and 936 ms postbounce for the nonrotating and slowly rotating models, respectively. The slowly rotating model explodes at ~650 ms postbounce, and the subsequent fallback accretion onto the proto–neutron star (PNS) results in BH formation. In addition, the standing accretion shock instability induces rotation of the PNS in the model that started with a nonrotating progenitor. Assuming conservation of specific angular momentum during BH formation, this corresponds to a BH spin parameter of a = J/M = 0.046. However, if no explosion sets in, all the angular momentum will eventually be accreted by the BH, resulting in a nonspinning BH. The successful explosion of the slowly rotating model drastically slows down the accretion onto the PNS, allowing continued cooling and contraction that results in an extremely high GW frequency (f ~ 3000 Hz) at BH formation, while the nonrotating model generates GW signals similar to our corresponding 2D simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Uncertainty quantification in elastic constants of SiC f /SiC m tubular composites using global sensitivity analysis

Silicon carbide fiber and silicon carbide matrix (SiC f /SiC m ) tubes produced through the chemical vapor infiltration process have become a candidate cladding material in nuclear applications. The performance of this composite is influenced by many variables such as braiding angle, porosity, material properties, etc., which vary over a range of values due to the inherent fluctuations in the manufacturing process. In this study, the variability in elastic constants of SiC f /SiC m composite has been quantified through multiscale finite element (FE) simulations, variable screening, and high-fidelity surrogate modeling. The key variables dominantly affecting the elastic constants of SiC f /SiC m tubes were identified using global sensitivity analysis. A surrogate to the high-fidelity FE-based model was used in Monte Carlo simulations to generate a hundred thousand samples from which the uncertainty in elastic constants was assessed. It turned out that the coefficient of variation was less than 10%.

Materials Science↗
DOE OSTI2022

More Realistic Intermediate Depth Dry Firn Densification in the Energy Exascale Earth System Model (E3SM)

Abstract Earth system models account for seasonal snow cover, but many do not accommodate the deeper snowpack on ice sheets (aka firn) that slowly transforms to ice under accumulating snowfall. To accommodate and resolve firn depths of up to 60 m in the Energy Exascale Earth System Model's land surface model (ELM), we add 11 layers to its snowpack and evaluate three dry snow compaction equations in multi‐century simulations. After comparing results from ELM simulations (forced with atmospheric reanalysis) with empirical data, we find that implementing into ELM a two‐stage firn densification model produces more accurate dry firn densities at intermediate depths of 20–60 m. Compared to modeling firn using the equations in the (12 layer) Community Land Model (version 5), switching to the two‐stage firn densification model (with 16 layers) significantly decreases root‐mean‐square errors in upper 60 m dry firn densities by an average of 41 kg m −3 (31%). Simulations with three different firn density parameterizations show that the two‐stage firn densification model should be used for applications that prioritize accurate upper 60 m firn air content (FAC) in regions where the mean annual surface temperature is greater than roughly −31°C. Because snow metamorphism, firn density, and FAC are major components in modeling ice sheet surface albedo, melt water retention, and climatic mass balance, these developments advance broader efforts to simulate the response of land ice to atmospheric forcing in Earth system models.

58 GEOSCIENCES↗
DOE OSTI2021

Effect of the Addition of a Low Equivalent Stress Mechanism to the Analysis of Geomechanical Behavior of the SPR West Hackberry Site

Sandia National Laboratories has long used the Munson-Dawson (M-D) model to predict the geomechanical behavior of salt caverns used to store oil at the Strategic Petroleum Reserve (SPR). Salt creep causes storage caverns to deform inward, thus losing volume. This loss of volume affects the salt above and around the caverns, puts stresses and strains on borehole casings, and creates surface subsidence which affects surface infrastructure. Therefore, accurate evaluation of salt creep behavior drives decisions about cavern operations. Parameters for the M-D model are typically fit against laboratory creep tests, but nearly all historic creep tests have been performed at equivalent stresses of 8 MPa or higher. Creep rates at lower equivalent stresses are very slow, such that tests take months or years to run, and the tests are sensitive to small temperature perturbations (<0.1°C). A recent collaboration between US and German researchers, recently characterized the creep behavior at low equivalent (deviatoric) stresses (<8 MPa) of salt from the Waste Isolation Pilot Plant. In addition, the M-D model was recently extended to include a low stress creep “mechanism”. This paper details new simulations of SPR caverns that use this extended M-D model, called the M-D Viscoplastic model. The results show that the inclusion of low stress creep significantly alters the prediction of steady-state cavern closure behavior and indicate that low stress creep is the dominant displacement mechanism at the dome scale. The implications for evaluating cavern and well integrity are demonstrated by investigating three phenomena: the extent of stress changes around the cavern; the predicted vertical strains applied to wellbore casings; and the evaluation of oscillating stress changes around the cavern due to oil sale cycles and their potential effect on salt fatigue.

58 GEOSCIENCES↗
DOE OSTI2023

Realistic large eddy and dispersion simulation experiments during project sagebrush phase 1

Typical dispersion simulations use mesoscale models at grid spacing greater than 1 km to simulate pollutant transport. Advances in computing power have allowed Large Eddy Simulations (LES), (grid spacing less than 100 m) to be run with boundary conditions in the grayscale regime (grid size of 1 km–100 m). We used the WRF mesoscale/LES model and the HYSPLIT particle dispersion model to simulate SF 6 transport from Project Sagebrush Phase 1, using a grayscale-aware scheme with 4 nests of 3 km, 1 km, 200 m, 67 m. The WRF/LES simulation was in satisfactory agreement with meteorological observations but the HYSPLIT SF6 simulations indicated insufficient sub-grid turbulence kinetic energy, especially in the lowest 100 m. This resulted in over predictions of the surface SF6 concentration and insufficient vertical transport. Best agreement with the Sagebrush observations was obtained with the WRF/LES large-scale and sub-grid turbulence, enhanced with a HYSPLIT parameterized turbulence profile.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A reinvestigation into Munson’s model for room closure in bedded rock salt

Accurate predictions of room closure are important for hazardous waste repositories in rock salt formations, such as the Waste Isolation Pilot Plant (WIPP). When Munson and co-workers simulated several room closure experiments conducted at the WIPP during the 1980’s and 1990’s, their simulated closure curves closely agreed with the closure measurements. A careful review of their work, however, raised concerns and prompted the reinvestigation in this paper. Here, to begin the reinvestigation, Munson’s legacy Room D closure simulation was reasonably recreated in a current-day finite element code. Next, special care was taken to obtain numerically converged results, re-introduce the anhydrite strata intermittently ignored by Munson, and calibrate the Munson–Dawson (M–D) constitutive model for salt as much as possible from laboratory test measurements. When this new model was used to simulate Room D’s closure, it under-predicted the horizontal and vertical closure rates by 2.34x and 3.10x, respectively, at 5.7 years after room excavation. As a result, the M–D model was extended to capture the newly established creep behavior at low equivalent stresses (<8 MPa) and replace the Tresca with the Hosford equivalent stress. Simulations using the new M–D model over-predicted the horizontal closure rate by 1.15x and under-predicted the vertical closure rate by 1.08x at 5.7 years, averaged over three room closure experiments. Although further improvements could be made, the new model has a stronger scientific foundation than Munson’s legacy model and appears ready for careful engineering use.

58 GEOSCIENCES↗