Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “inverse test”

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 37 records · Page 2

EGS-Collab Experiment 1: Time-lapse ERT Data and E4D Inversion Files for 10-24-2018 through 11-07-2018 Flow Test

This data submission includes the raw time-lapse ERT (electrical resistivity tomography) monitoring data, flow system data, operator logs, E4D (https://e4d.pnnl.gov) inversion files, and meta-data necessary to reproduce the 4D ERT inversion for the Oct. 24 through Nov. 7 2018 post-stimulation flow test in test bed 1. The tests were done at the Sanford Underground Research Facility at Homestake Mine in South Dakota.

15 GEOTHERMAL ENERGY↗

Dynamics of dark matter misalignment through the Higgs portal

A light singlet scalar field feebly coupled through the super-renormalizable Higgs portal provides a minimal and well-motivated realization of ultra-light bosonic dark matter. We study the cosmological production of dark matter in this model by elucidating the dynamics of two sources of scalar field misalignment generated during the radiation era. For large scalar masses (above $\mathcal{O}$ (10 −3 eV)), dark matter is produced through thermal misalignment, by which the scalar field is driven towards large field values as a result of the finite-temperature effective potential. The dominance of thermal misalignment in this mass range leads to a sharp relic abundance prediction which is, to a significant extent, insensitive to the initial conditions of the scalar field. On the other hand, for low mass scalars (below $\mathcal{O}$ (10 −5 eV)), dark matter is produced via VEV misalignment, which is caused by the induced scalar field vacuum expectation value triggered by the electroweak phase transition. We show that the relic abundance in this low mass range is sensitive to the scalar field initial conditions. In the intermediate mass range, the relic abundance is a consequence of a competition between thermal misalignment and VEV misalignment, which can potentially lead to novel forced resonance effects which cause a recurring enhancement and suppression in the late time oscillation amplitude as a function of the scalar mass. We compare our relic abundance predictions with constraints and projections from equivalence principle and inverse square law tests, stellar cooling, resonant molecular absorption, and observations of extra-galactic background light and diffuse X-ray backgrounds. New experimental ideas are needed to probe most of the cosmologically motivated regions of parameter space.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Resolving the fast ion distribution from imaging neutral particle analyzer measurements

A recently developed imaging neutral particle analyzer (INPA) on the DIII-D tokamak (Du 2018 Nucl. Fusion 58 082006) enables fast ion velocity-space tomography of high fidelity at the interrogated phase space. To accomplish this, the spatial and energy depending fast ( E < 80 keV) neutral flux towards the INPA stripping foils is calculated with FIDASIM and a newly developed code INPASIM simulates the INPA instrumental response to this neutral flux. Included in INPASIM is the neutral-foil interaction, the Larmor orbit tracing between the foil and the phosphor, the phosphor response to the incident ion flux as well as camera focusing. Benefiting from heavy, localized velocity-space weights and excellent signal to noise, computed tomography using the Ridge regression method is able to successfully reconstruct fine-scale velocity-space structures produced by multiple neutral beams separated by as small as ~3 keV in tests. Applying the inversion method to a sawtooth crash event reveals a significant profile flattening of confined passing particles across q = 1 flux surface, as well as a redistribution of fast ions into the trapped orbits at the plasma edge close to the last closed flux surface.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Sensitivity of modelled passive margin stratigraphy to variations in sea level, sediment supply and subsidence

Abstract We produced a 10 Myr synthetic stratigraphic section using a forward stratigraphic model that generates marine deltaic stratigraphy over geological timescales. We recursively fit the model using a Bayesian inversion algorithm to test: (1) if it could be accurately reconstructed; (2) if the parameters used to create it could be recovered; and (3) the sensitivity of the model output to given model parameters and the attendant physical processes. The original synthetic stratigraphic section was produced with cyclical sea‐level variations of 40 and 30 m with 2.4 and 10 Myr periods respectively. Sediment was also supplied cyclically, in 2.4 and 10 Myr cycles with amplitudes of 30 and 80 tons/100 kyr, respectively, varying from a mean of 232 tons/100 kyr. Parameter values were sampled to fit the model using a Markov chain Monte Carlo algorithm, resulting in a ±5 m (1σ) variation between the experimental output and the original. Sea level varied by ±7 m (1σ) within the posterior distribution of parameters. As a result, both the 10 Myr and 2.4 Myr sea‐level cycles could be extracted from the original output. The variation in sediment supply was approximately ±38 tons/100 kyr (1σ) and, as a result, only the larger long‐term supply variations could be accurately recovered in refitting the model. The variation in thermal, flexural and total subsidence across those parameter sets is less than ±10 m (1σ). The original section experienced 150 m of total subsidence at the depocentre. Our results demonstrate the distinct and interpretable imprint of sea level and subsidence on continental margin stratigraphy can be quantified. Moreover, we conclude that sea‐level change produces a defined effect on the geometries of stratigraphic architecture, and that techniques applied for the purpose of delineating sea‐level variation from continental margin strata have a well‐founded conceptual basis.

Schmelz, William J.↗

Field and Model Data Associated with the Manuscript “Drivers of Streamflow Intermittency in Humid Regions: 2. Evaluating Controls on Flow Persistence in an Urbanized Catchment”

This package contains field data, modeling files, and scripts supporting the investigation of the drivers of streamflow intermittency in an urbanized catchment. It includes the field data collected from electrical resistivity tomography (ERT) surveys, distributed temperature sensing (DTS), continuous self-potential (SP) monitoring, groundwater and stilling well. In addition, it contains the data and results of the coupled water- and electrical-flow model developed using the COMSOL Multiphysics and Advanced Terrestrial Simulator (ATS), as well as software files and Jupyter notebooks used to process the data and generate figures in the manuscript submitted for peer review. The data archive is organized in the following directories: 1) Climate Includes hourly precipitation and daily evapotranspiration time series (2024 – 2025) provided as CSV files, alongside a text file detailing dataset units. 2) Coupled_model Field_Application subfolder contains the ATS XML input scripts, data files, output data for the SP site. It also contains the Jupyter notebook (Plot_final_calib.ipynb) to visualize the results of the modeled SP, stream-groundwater exchange and moisture content. The flow model simulation is executed using the ATS XML scripts and the included Python script (generate_data_set.py) to convert ATS output to COMSOL-ready input. COMSOL Multiphysics template (.m can only be used with COMSOL with MATLAB) is executed using the ATS output data to simulate the potential field. 3) Discharge Includes the electrical conductivity (EC) time series (provided as CSV files) from salt slug injections. It also includes the Jupyter notebook (Discharge_process.ipynyb) used to estimate discharge. All discharge measurements collated into rating_curve_processed.csv 4) DTS Contains collated DTS data including raw Stokes and anti-Stokes measurement (provided as .h5 file). It also includes DTS processing.ipynb, a Jupyter notebook for calibrating the DTS data using dts_calibration Python package. cooler_calibration.csv is the DTS calibration CSV used in the calibration sequence. 5) ERT Contains raw resistivity data (provided as CSV files), spatial location of each of the electrodes (provided as CSV files), and files used for the resistivity inversion. 6) Slug_test Includes the slug test data at all the groundwater wells provided as CSV files, as well as the Jupyter notebook (Slug_test.ipynb) for calculating hydraulic conductivity. 7) SP Contains the SP data collected in field at the SP sites (provided as CSV files). 8) Well_data Contains two subfolders: 1) Raw, which provides unprocessed pressure, electrical conductivity and temperature timeseries downloaded from the loggers in all the groundwater and stilling wells, and 2) Processed, which contains sorted, QA/QC timeseries data for each well. The data archive also contains data_process.ipynb, a Jupyter notebook used for field data analysis and generating figures (plotting well, SP, climate, and discharge data, as well as calculating head gradient at sites with nested groundwater wells). Note: Code files (.ipynb, .py, .xml) can be opened in any standard code editor, .exo file can be viewed using Paraview, .h5 files can be opened using HDFView software and h5py Python package, and .resipy file can be opened with the open-source ResIPy software.

ATS↗

Hydrogen Permeation in FeCrAl APMT Alloy for Accident Tolerant Fuel Cladding

Iron-chromium-aluminum (FeCrAl) alloys such as APMT (advanced powder metallurgy tubing) are candidate materials to replace zirconium alloys for the light water reactor fuel cladding. This alloy meets the requirements to be a material more tolerant of high-temperature accidents than the current zirconium alloys. One concern is that the use of FeCrAl may result in an increase in tritium presence in the coolant compared to the current design. The aim of the current research was to obtain effective diffusion coefficients (D eff ) for hydrogen through APMT using the Devanathan-Stachurski cell. Results showed that at 30°C the D eff value was 2.8 × 10 -8 cm 2 /s. Results also showed a linear relationship between the permeated hydrogen flux and the inverse of the test specimen thickness, which demonstrated the validity of the permeation measurements.

36 MATERIALS SCIENCE↗

Fast and sensitive measurements of sub-3 nm particles using Condensation Particle Counters For Atmospheric Rapid Measurements (CPC FARM)

New particle formation (NPF) is the atmospheric process whereby gas molecules react and nucleate to form detectable particles. NPF has a strong impact on Earth's radiative balance as it produces roughly half of global cloud condensation nuclei. However, the time resolution and sensitivity of current instrumentation are inadequate in measuring the size distribution of sub-3 nm particles, the particles critical for understanding NPF. Here we present the Condensation Particle Counters For Atmospheric Rapid Measurements (CPC FARM), a method to measure the concentrations of freshly nucleated particles. The CPC FARM consists of five CPCs operating in parallel, each configured to operate at different detectable particle sizes between 1–3 nm. This study explores two methods to calculate the size distribution from the differential measurements across the CPC channels. The performance of both inversion methods was tested against the size distribution measured by a pair of stepping particle mobility sizers (SMPSs) during an ambient air sampling study in Pittsburgh, PA. Observational results indicate that the CPC FARM is more accurate with higher time resolution and sensitivity in the sub-3 nm range compared to the SMPS.

Cheng, Darren [Carnegie Mellon Univ., Pittsburgh, ↗

InverseBench: Inverse design benchmark suite that contains inverse problems from science and engineering (InverseBench) v0.0.1

A software package that contains three inverse design blackbox problems to investigate the efficiency and accuracy of inverse design machine learning models. The software contains highly accurate forward machine learning models that can be used to assess the inverse predictions. The package also contains separate test data for each problem. The inverse design problems that are in the package are: airfoil inverse design, scalar boundary reconstruction and photonic surfaces inverse design.

Grbcic, Luka [Lawrence Berkeley National Laborator↗

Inversion of downhole electrical measurements for proppant mapping using very fast simulated annealing

Proppant mapping is critical for optimizing fracture treatment design and improving wells’ productivity. An electrode-based resistivity tool concept was developed earlier for proppant mapping in cased-hole wells. An array of insulating gaps is installed and cemented in place as a permanent part of the casing string. The electrical measurements are performed by imposing a voltage across each insulating gap, one at a time, before and after hydraulic fracture operations. The voltages across other insulating gaps near the transmitter gap are recorded. The method relies on direct excitation of the casing, which is expected to overcome the severe limitations of induction tools in cased-hole wells. A forward model based on a finite volume method has been developed to simulate the tool’s response to one or multiple fractures. To enable the implementation of such a practical system in multistage fractured horizontal wells, a fast and robust inversion approach is required. To that end, we have developed a divide-and-conquer approach based on a global optimization algorithm very fast simulated annealing (VFSA). Specifically, the original inverse problem is divided into subproblems and each subproblem can be solved separately using VFSA. The results indicate that our approach can invert the data and output widths and radii of multiple fractures without requiring a large number of forward simulations. The robustness of the inverse solver is also tested by adding Gaussian noise to the synthetic data. We tested example cases that demonstrate that when up to 5% noise is introduced, VFSA still provides very accurate inversion results with moderate uncertainties. Inversion results with some more realistic conditions, e.g., tilted fractures, complex fractures, and so on, are also presented.

Geochemistry & Geophysics↗

Hydraulic Testing of the Lost Lake Aquifer Near Recovery Well RWM001R

An aquifer pumping test was conducted on the Lost Lake Aquifer Zone (LLAZ) at the recently installed recovery well RWM001R in accordance with the approved test plan (Dixon, 2024). The objective of the testing was to determine baseline well performance parameters and aquifer hydraulic conductivity. This testing consisted of two step-drawdown tests to determine well performance properties and two constant pumping rate aquifer tests with recovery periods to determine aquifer hydraulic properties. Well performance parameters determined included specific capacity, well efficiency, and head loss coefficients. The results of the step-drawdown testing are presented in Tables ES1 and ES2 whereas the results of the constant pumping rate testing are presented in Table ES3. The average specific capacity of RWM001R operating at pump capacity (~54 gpm) was 2.2 gpm/ft using data from the final step of both tests. The average well efficiency determined from both step-drawdown tests was 56% using data from the final step of both tests. For both step-drawdown tests, well efficiency was inversely related to the pumping rate with the average efficiency decreasing from approximately 79% to 56%. Aquifer head loss coefficients from both tests were comparable (Table ES2).

54 ENVIRONMENTAL SCIENCES↗

Deep learning-enhanced design for functionally graded auxetic lattices

Auxetic materials with counterintuitive negative Poisson’s ratio have been of significant interest due to potential applications across diverse engineering fields. Functionally grading such auxetics further enables customization of the structural response and harnesses the potential for multi-functional applications. However, a critical challenge in designing functionally graded lattices is to efficiently determine the spatial variation of the functional gradient and the corresponding geometric designs to achieve the desired response. In this paper, a highly efficient deep learning-based inverse design framework for functionally graded tetra-petal auxetics with spatially tailored properties is presented. This framework significantly improves the efficiency of tailoring functionally graded auxetics where many unit cells need to be tailor-designed. The graded tetra-petal auxetics obtained from the inverse design framework are additively manufactured and subjected to impact tests. The results show superior impact performance compared with uniform designs, demonstrating the effectiveness of the proposed inverse design framework, which can be inspirable to promote advanced structures/materials with enhanced impact resistance.

42 ENGINEERING↗

Inversion-based correction of Double-Torsion (DT) subcritical crack growth tests for crack profile geometry

Because of its simplicity and the ability to produce a stable, slow-propagating crack, the Double-Torsion (DT) method has been used widely for investigating the critical and subcritical propagation of a slow-propagating tensile (mode-I) crack. However, to determine the complex relationship between the crack velocity $\mathcal{v_c}$ vs. the strain energy release rate $\mathscr{G}$ (or the stress intensity factor K) from laboratory measurements, several corrections must be made to account for the impact of sample and crack geometry. Particularly, DT test typically produces a crack with a curved edge profile instead of a straight line, causing the local $\mathcal{v_c}$ and $\mathscr{G}$ vary along the crack front. The experimentally measured $\mathcal{v_c}$ and $\mathscr{G}$ data merely reflect collective, averaged behavior of the crack. This makes inversion for the intrinsic, “true” crack growth kinetics necessary, based upon the knowledge of the crack geometry. Simple and effective correction methods have been proposed and validated for the slow, chemical-reaction-controlled part (Region I) of the $\mathcal{v_c}$-$\mathscr{G}$ curve. However, reliable methods for the highly nonlinear, transport-dominated part (Region II) and its sudden transition to the dynamic propagation part (Region III) are still lacking. Here we propose a method for determining the intrinsic $\mathcal{v_c}$-$\mathscr{G}$ relationship across all three Regions based upon DT test data, using a simple model function and its numerical inversion. The performance of this approach is examined and demonstrated using both synthetic and laboratory data for subcritical crack growth in soda lime glass.

42 ENGINEERING↗

Invertible Neural Networks for Airfoil Design

We report the airfoil design problem, in which an engineer seeks a shape with desired performance characteristics, is fundamental to aerodynamics. Design workflows traditionally rely on iterative optimization methods using low-fidelity integral boundary-layer methods as higher-fidelity adjoint-based computational fluid dynamics methods are computationally expensive. Surrogate-based approaches can accelerate the design process but still rely on some iterative inverse design procedure. In this work, we leverage emerging invertible neural network (INN) tools to enable the rapid inverse design of airfoil shapes for wind turbines. INNs are specialized deep-learning models with well-defined inverse mappings. When trained appropriately, INN surrogate models are capable of forward prediction of aerodynamic and structural quantities for a given airfoil shape as well as inverse recovery of airfoil shapes with specified aerodynamic and structural characteristics. The INN approach offers a roughly 100 times speed-up compared to adjoint-based methods for inverse design. We demonstrate the INN tool for inverse design on three test cases of 100 airfoils each that satisfy the performance characteristics close to those of airfoils used in wind-turbine blades. All generated shapes satisfy the desired aerodynamic characteristics, demonstrating the success of the INN approach for inverse design of airfoils.

17 WIND ENERGY↗

Ferroelectricity of wurtzite Al 1−x Hf x N heterovalent alloys

Thin films of aluminum hafnium nitride (Al1−xHfxN) were synthesized via reactive magnetron sputtering for Hf contents up to x = 0.13. X-ray diffraction showed a single c-axis oriented wurtzite phase for all films. Hard x-ray photoelectron spectroscopy demonstrated homogeneous Al:Hf distribution through the thin films and confirmed their insulating character. A collection of complementary tests showed unambiguous polarization inversion, and thus ferroelectricity in multiple samples. Current density vs electric field hysteresis measurements showed distinct ferroelectric switching current peaks, the piezoelectric coefficient d 33,f,meas measured using a double beam laser interferometer (DBLI) showed a reversal in sign with similar magnitude, and anisotropic wet etching confirmed field-induced polarization inversion. This demonstrates the possibility of using tetravalent–and not just trivalent–alloying elements to enable ferroelectricity in AlN-based thin films, highlighting the compositional flexibility of ferroelectricity in wurtzites and greatly expanding the chemistries that can be considered for future devices.

36 MATERIALS SCIENCE↗

Seismic Event Characterization Using Full Moment Tensors on the Hypersphere

Moment tensor solutions provide insights into the deformation that has occurred in the source region of a seismic event and are therefore of great value in identifying different types of seismic sources, such as when monitoring for underground nuclear tests. Despite this utility, inversion of waveforms recorded by seismometers for their full seismic moment tensor is not yet routine, and development of robust methods to classify events based on this information is in its infancy. Here, we assemble an inventory of 1405 full moment tensor solutions that include explosive, earthquake, and collapse events, and investigate the use of anisotropic probability distribution functions on the 5D hypersphere to discriminate between these sources. Using a Bayesian classifier, we obtain optimal success rates of 98.4% across all events and demonstrate that modification of the prior probabilities provides a natural way to alter the balance between not missing desirable events (such as explosions) versus misclassifying large numbers of undesired events (such as earthquakes). The approach is specifically designed to progress from traditional, bipolar event screening metrics to more generalized event identification across multiple types of seismic sources. Despite current databases containing insufficient numbers of events to definitively demonstrate at present, we also find intriguing evidence of subgroupings within individual source populations on the hypersphere, for example, between chemical and nuclear explosions, raising the potential possibility of discriminating between these event types in the future.

Geosciences↗

Hydraulic Testing of the Lost Lake Aquifer Near Recovery Well RWM019

An aquifer pumping test was conducted on the Lost Lake Aquifer Zone (LLAZ) at the recently installed recovery well RWM019 in accordance with the approved test plan (Dixon, 2020). The objective of the testing was to determine baseline well performance parameters and aquifer hydraulic conductivity. This testing consisted of a step-drawdown test to determine well performance properties, a constant pumping rate aquifer test to determine aquifer hydraulic properties, and a post-test aquifer recovery monitoring period also used to estimate aquifer hydraulic properties. Well performance parameters determined included specific capacity, well efficiency, and head loss coefficients. The specific capacity of RWM019 was determined to be approximately 2.6 gpm/ft of drawdown based on the final step of the step-drawdown test. Well efficiency was inversely related to pumping rate and decreased from 87% to 81% over a pumping range of approximately 49 to 79 gpm. The aquifer head loss coefficient was determined to be 2.3 ft/ft 3 /min and the well loss coefficient was determined to be 0.05 min 2 /ft 5 .

42 ENGINEERING↗

Dataset for SIAM MPI23 Project "Model inversion for complex physical systems using low-dimensional surrogates"

Training and testing data for the SIAM MPI23 workshop "Model inversion for complex physical systems using low-dimensional surrogates". The data consists of ensembles of input and output pairs corresponding to queries of a 2D saturated groundwater flow model. The inputs consist of random vectors which map to discretized model parameter fields via a Kosambi-Karhunen-Loève expansion. Output corresponds to discretized pressure fields.

basis adaptation, conditional Karhunen-Loève expan↗

Massively parallel modeling and inversion of electrical resistivity tomography data using PFLOTRAN

Abstract. Electrical resistivity tomography (ERT) is a broadly accepted geophysical method for subsurface investigations. Interpretation of field ERT data usually requires the application of computationally intensive forward modeling and inversion algorithms. For large-scale ERT data, the efficiency of these algorithms depends on the robustness, accuracy, and scalability on high-performance computing resources. In this regard, we present a robust and highly scalable implementation of forward modeling and inversion algorithms for ERT data. The implementation is publicly available and developed within the framework of PFLOTRAN, an open-source, state-of-the-art massively parallel subsurface flow and transport simulation code. The forward modeling is based on a finite-volume discretization of the governing differential equations, and the inversion uses a Gauss–Newton optimization scheme. To evaluate the accuracy of the forward modeling, two examples are first presented by considering layered (1D) and 3D earth conductivity models. The computed numerical results show good agreement with the analytical solutions for the layered earth model and results from a well-established code for the 3D model. Inversion of ERT data, simulated for a 3D model, is then performed to demonstrate the inversion capability by recovering the conductivity of the model. To demonstrate the parallel performance of PFLOTRAN's ERT process model and inversion capabilities, large-scale scalability tests are performed by using up to 131 072 processes on a leadership class supercomputer. These tests are performed for the two most computationally intensive steps of the ERT inversion: forward modeling and Jacobian computation. For the forward modeling, we consider models with up to 122 ×106 degrees of freedom (DOFs) in the resulting system of linear equations and demonstrate that the code exhibits almost linear scalability on up to 10 000 DOFs per process. On the other hand, the code shows superlinear scalability for the Jacobian computation, mainly because all computations are fairly evenly distributed over each process with no parallel communication.

58 GEOSCIENCES↗