Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “inversion methods”

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 973 records · Page 54

Physics-assisted generative adversarial network for X-ray tomography

X-ray tomography is capable of imaging the interior of objects in three dimensions non-invasively, with applications in biomedical imaging, materials science, electronic inspection, and other fields. The reconstruction process can be an ill-conditioned inverse problem, requiring regularization to obtain satisfactory results. Recently, deep learning has been adopted for tomographic reconstruction. Unlike iterative algorithms which require a distribution that is known a priori , deep reconstruction networks can learn a prior distribution through sampling the training distributions. In this work, we develop a Physics-assisted Generative Adversarial Network (PGAN), a two-step algorithm for tomographic reconstruction. In contrast to previous efforts, our PGAN utilizes maximum-likelihood estimates derived from the measurements to regularize the reconstruction with both known physics and the learned prior. Compared with methods with less physics assisting in training, PGAN can reduce the photon requirement with limited projection angles to achieve a given error rate. The advantages of using a physics-assisted learned prior in X-ray tomography may further enable low-photon nanoscale imaging.

47 OTHER INSTRUMENTATION↗

Rice husk and charred husk amendments increase porewater and plant Si but water management determines grain As and Cd concentration

Abstract Purpose Rice is a staple crop worldwide and a silicon (Si) hyperaccumulator with Si levels reaching 5–10% of its mass; this can result in desilication and Si-deficiency if plant residues are not managed correctly. Rice is also uniquely subject to arsenic (As) and cadmium (Cd) contamination depending on soil conditions. Our goal is to quantify the effects of rice husk (a Si-rich milling byproduct) amendments and different water management strategies on rice uptake of Si, As, and Cd. Methods We employed 4 husk amendment treatments: Control (no husk), Husk (untreated husk), Biochar (husk pyrolyzed at 450 °C), and CharSil (husk combusted at > 1000 °C). Each of these amendments was studied under nonflooded, alternate wetting and drying (AWD), and flooded water management in a pot study. Porewater chemistry and mature plant elemental composition were measured. Results Husk and Biochar treatments, along with flooding, increased porewater and plant Si. Vegetative tissue As decreased with increasing porewater Si, but grain As and plant Cd were primarily controlled by water management. Grain As and Cd were inversely correlated and are simultaneously minimized in a redox potential (Eh) range of 225–275 mV in the studied soil. Ferrihydrite in root iron plaque decreased As translocation from porewater to grain, but amendments were not able to increase plaque ferrihydrite content. Conclusion We conclude moderate husk amendment rates (i.e., 4 years’ worth) with minimal pretreatment strongly increases rice Si content but may not be sufficient to decrease grain As in low Si and As soil.

36 MATERIALS SCIENCE↗

Calibration of elastoplastic constitutive model parameters from full-field data with automatic differentiation-based sensitivities

Here, we present a framework for calibration of parameters in elastoplastic constitutive models that is based on the use of automatic differentiation (AD). The model calibration problem is posed as a partial differential equation-constrained optimization problem where a finite element (FE) model of the coupled equilibrium equation and constitutive model evolution equations serves as the constraint. The objective function quantifies the mismatch between the displacement predicted by the FE model and full-field digital image correlation data, and the optimization problem is solved using gradient-based optimization algorithms. Forward and adjoint sensitivities are used to compute the gradient at considerably less cost than its calculation from finite difference approximations. Through the use of AD, we need only to write the constraints in terms of AD objects, where all of the derivatives required for the forward and inverse problems are obtained by appropriately seeding and evaluating these quantities. We present three numerical examples that verify the correctness of the gradient, demonstrate the AD approach's parallel computation capabilities via application to a large-scale FE model, and highlight the formulation's ease of extensibility to other classes of constitutive models.

42 ENGINEERING↗

Filtering adaptive output feedback control for multivariable nonlinear systems with mismatched uncertainties and unmodeled dynamics

This article synthesizes a filtering adaptive output feedback controller for multivariable nonlinear systems with mismatched uncertainties and unmodeled dynamics. The multivariable nonlinear systems under consideration have both matched and mismatched uncertainties, which satisfy the semiglobal Lipschitz condition. The unmodeled dynamics are bounded-input bounded-output stable. By adopting an estimation/cancellation strategy, a piecewise constant adaptive law drives the estimation error to zero at every time instant, which yields the adaptive parameters; a disturbance rejection control law is designed to compensate the nonlinear uncertainties within the bandwidth of low-pass filters. The matched uncertainties are cancelled directly by adopting their opposite in the control signal, while a dynamic inversion of the system is required to eliminate the effect of the mismatched uncertainties on the output. A feedforward control law is designed to track a given command. Since the virtual reference system defines the best performance that can be achieved by the closed-loop system, the uniform performance bounds are derived for the states and control signals via comparison. Both numerical and practical examples are provided to illustrate the effectiveness of the proposed filtering adaptive output feedback control architecture, comparisons with the model reference adaptive control demonstrates the superiority of the proposed control method.

filtering adaptive control, mismatched uncertainit↗

Electron-beam–controlled deflection of near-infrared laser in semiconductor plasma

A timing method for experiments on the interaction of a near-infrared laser and an ultra-relativistic electron beam via a semiconductor plasma switch is experimentally validated. As an intermediate medium, a thin Si plate is excited by the energetic, intense electron beam to produce a semiconductor plasma, which in turn deflects counter-colliding laser light having 1 μm wavelength. An electron beam of sub-nC charge sufficiently induces the needed electron number density gradient of 1×10 20 cm -3 per tens of μm length at the interaction point. Demonstration during an inverse Compton scattering experiment by a counter-colliding electron beam of 300 pC and 70 MeV with an Nd: YAG laser at a wavelength of 1 μm is reported here.

43 PARTICLE ACCELERATORS↗

Grad–Shafranov equilibria via data-free physics informed neural networks

A large number of magnetohydrodynamic (MHD) equilibrium calculations are often required for uncertainty quantification, optimization, and real-time diagnostic information, making MHD equilibrium codes vital to the field of plasma physics. In this paper, we explore a method for solving the Grad–Shafranov equation by using physics-informed neural networks (PINNs). For PINNs, we optimize neural networks by directly minimizing the residual of the partial differential equation as a loss function. We show that PINNs can accurately and effectively solve the Grad–Shafranov equation with several different boundary conditions, making it more flexible than traditional solvers. This method is flexible as it does not require any mesh and basis choice, thereby streamlining the computational process. We also explore the parameter space by varying the size of the model, the learning rate, and boundary conditions to map various tradeoffs such as between reconstruction error and computational speed. Additionally, we introduce a parameterized PINN framework, expanding the input space to include variables such as pressure, aspect ratio, elongation, and triangularity in order to handle a broader range of plasma scenarios within a single network. Parameterized PINNs could be used in future work to solve inverse problems such as shape optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Neural posterior unfolding

Differential cross section measurements are the currency of scientific exchange in particle and nuclear physics. A key challenge for these analyses is the correction for detector distortions, known as deconvolution or unfolding. Binned unfolding of cross section measurements traditionally rely on the regularized inversion of the response matrix that represents the detector response, mapping pre-detector (`particle level') observables to post-detector (`detector level') observables. In this paper we introduce Neural Posterior Unfolding, a modern, Bayesian approach that leverages normalizing flows for unfolding. By using normalizing flows for neural posterior estimation, NPU offers several key advantages including implicit regularization through the neural network architecture, fast amortized inference that eliminates the need for repeated retraining, and direct access to the full uncertainty in the unfolded result. In addition to introducing NPU, we implement a classical Bayesian unfolding method called Fully Bayesian Unfolding (FBU) in modern Python so it can also be studied. These tools are validated on simple Gaussian examples and then tested on simulated jet substructure examples from the Large Hadron Collider (LHC). We find that the Bayesian methods are effective and worth additional development to be analysis ready for cross section measurements at the LHC and beyond.

Analysis and statistical methods↗

MrHyDE v.1.0

SAND2024-01324O MrHyDE, which stands for Multi-resolution Hybridized Differential Equations, is a general-purpose C++ package for the solution of coupled multiphysics and multiscale systems on massively parallel computing systems. MrHyDE is designed to enable moving beyond forward simulation for multiscale applications which includes optimization, control, uncertainty quantification, and stochastic inversion. The framework provides interfaces to several packages within the Trilinos framework and leverages automatic differentiation to enable adjoint capabilities for large-scale, gradient-based optimization. MrHyDE provides automated multiscale capabilities through a subgrid model interface and multiscale Dirichlet-to-Neumann maps. For extreme-scale applications, MrHyDE provides in situ data-compression algorithms to reduce memory requirements while maintaining performance. MrHyDE is a general-purpose, computational framework for the solution of multiscale and multiphysics applications. It uses a combination of structure-preserving, physics-compatible discretizations, fully implicit methods, multi-resolution schemes, or fully explicit methods. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

SciDAC↗

Extraction and selective promotion of zero-group velocity and cutoff frequency resonances in bi-dimensional waveguides using the electromechanical impedance method

This study showcases the electromechanical impedance (EMI) technique for extracting and promoting zero-group velocity (ZGV) and cutoff frequency resonances in a waveguide structure. Here, we identify the mechanisms of multiple resonances in the EMI spectra via a wave propagation perspective. Both simulation and experiments reveal the fact that sharp resonances in the conductance spectra are associated with either ZGV or cutoff frequency points. Consequently, we design four test configurations to enhance local resonances by aligning induced motions with considered mode shapes. Reasonable agreement between simulation and experiment results is observed. We evaluate the performance of considered configurations in terms of mode enhancement, and configurations that can selectively promote certain mode families are summarized. This study also shines the light on the EMI technique for quantitative non-destructive evaluation (NDE) by potentially supporting the inverse characterization of mechanical properties of host structures.

42 ENGINEERING↗

Hierarchical, Self-Assembled Metasurfaces via Exposure-Controlled Reflow of Block Copolymer-Derived Nanopatterns

Here, nanopatterning for the fabrication of optical metasurfaces entails a need for high-resolution approaches like electron beam lithography that cannot be readily scaled beyond prototyping demonstrations. Block copolymer thin film self-assembly offers an attractive alternative for producing periodic nanopatterns across large areas, yet the pattern feature sizes are fixed by the polymer molecular weight and composition. Here, a general strategy is reported that overcomes the limitation of fixed feature size by treating the copolymer thin film as a hierarchical resist, in which the nanoscale pattern motif is defined by self-assembly. Feature sizes can then be tuned by thermal reflow controlled locally by irradiative crosslinking or chemical alteration using lithographic ultraviolet light or electron beam exposure. Using blends of polystyrene-block-poly(methylmethacrylate) (PS-b-PMMA) with PS and PMMA homopolymer, we demonstrate both self-assembled PS grating and hexagonal hole patterns; exposure-controlled reflow is then used to reduce hole diameter by as much as 50% or increase PS grating linewidth by more than 180%. Transferring these nanopatterns, or their inverse obtained by a lift-off approach, into silicon yields structural colors that may be prescriptively controlled based on nanoscale feature size. Furthermore, patterned exposure enables area-selective feature size control, yielding uniform structural color patterns across centimeter square areas. Electron beam lithography is also used to show that the lithographic resolution of this selective-area control can be extended to the nanoscale dimensions of the self-assembled features. The exposure-controlled reflow approach demonstrated here takes a pivotal step towards fabricating complex, hierarchical optical metasurfaces using scalable self-assembly methods.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

A deep learning-based workflow for fast prediction of 3D state variables in geological carbon storage: A dimension reduction approach

Deep learning (DL) models are extensively used as surrogate models for high-fidelity simulations of multiphase fluid flow in porous media at large scales, enabling fast forecasts of the spatial–temporal evolution of three-dimensional (3D) state variables in geological carbon storage (GCS). However, training these models in high-dimensional space remains computationally demanding and prone to overfitting because of limited training data. This paper presents a novel workflow to address these challenges by integrating dimension reduction (DR) methods. Here, the proposed workflow employed pre-trained DR models to extract the latent variables of geological models and state variables and utilized the multi-layer perceptron (MLP) for constructing mapping functions between the input and output variables in latent spaces. Subsequently, the pre-trained reconstruction models converted the MLP-predicted latent state variables to their original high-dimensional form. Furthermore, we proposed a novel strategy for the DR and reconstruction of 3D saturation fields to account for the unique data characteristics of sparsity, nonuniformity, and discontinuity. The proposed strategy applied PCA and inverse PCA for 2D average saturation fields and developed a DL-based 3D reconstruction model, leveraging three 2D average saturation fields as input to produce a 3D saturation field as output. The pre-training of DR and reconstruction models and training of MLP models were conducted on 84 Gulf of Mexico (GoM) simulations and evaluated on 12 testing simulations. Each simulation contained 720 monthly time steps, with the first 360 months as the injection period and the rest as the post-injection period. The proposed workflow, incorporating DR and DL models, accurately predicts the normalized 3D pressure fields, achieving mean square error (MSE) of 2.92 × 10 -7 compared to the ground truth obtained from a full-physics simulator. Furthermore, the proposed strategy outperformed PCA and convolutional autoencoder (CAE) models on 3D saturation fields, resulting in minor workflow prediction errors with an MSE of 2.93 × 10 -5 . The results suggest the proposed workflow provides sufficient predictive fidelity across temporal and spatial scales, and enables a speedup of 160 times compared to the full-physics simulator, facilitating improved decision-making and risk assessment for large-scale GCS management in real-time scenarios.

3D reconstruction model↗

Evaluation of Aliphatic Alcohols for CO 2 Capture Using the Characteristic Carbonate Frequency

Control over CO 2 capture and utilization are important scientific and technological challenges. Although a variety of amine absorbents are used for capture, releasing the captured CO 2 is often difficult and limits their recyclability. Therefore, it is crucial to control the strength of the CO 2 bond with the absorbent. Furthermore, it is desirable to use a method that can conveniently report the strength of this bond. This motivates exploring adducts of CO 2 with alcohols in the presence of a base, using vibrational spectroscopy to report on the bond strength. Although reactions of alcohols with CO 2 to form alkyl carbonates are known, a systematic study of these adducts has not been conducted. Here we show formation of alkyl carbonates by a series of alcohols spanning the pK a range 9.5 to 16.8. We show experimental and computational results for the frequency of the characteristic asymmetric stretch of the carbonate and demonstrate that it correlates inversely with the pK a of the alcohol. Based on computations of the bond lengths and previous work, we propose that this frequency also correlates inversely with the adduct strength. As a result, this work extends the scope of CO 2 capture reagents and inspires further research in tuning alcohols as reversible absorbents.

Alcohols↗

Drive-pressure optimization in ramp-wave compression experiments through differential evolution

Ramp-wave dynamic-compression experiments are used to examine quasi-isentropic loading paths in materials. The gradual and continuous increase in pressure created by ramp waves make these types of experiments ideal for studying nonequilibrium material behavior, such as solidification kinetics. In ramp-wave compression experiments, the input drive pressure to the experimental setup may be exerted through one of a number of different mechanisms (e.g., magnetic fields, gas-gun-driven impactors, or high-energy lasers) and is generally required for simulating such experiments. Yet, regardless of the specific mechanism, this drive pressure cannot be measured directly (measurements are generally taken at a location near the back of the experimental setup through a transparent window), leading to an inverse problem where one must determine the drive pressure at the front of the experimental setup (i.e., the input) that corresponds to the particle velocity (the output) measured near the back of the experimental setup. Furthermore, we solve this inverse problem using a heuristic optimization algorithm, known as differential evolution, coupled with a multiphysics, hydrodynamics code that simulates the compression of the experimental setup. By running many rounds of forward simulations of the experimental setup, our optimization process iteratively searches for a drive pressure that is optimized to closely reproduce the experimentally measured particle velocity near the back of the experimental setup. While our optimization methodology requires a significant number of hydrodynamics simulations to be conducted, many of these can be performed in parallel, which greatly reduces the time cost of our methodology. One novel aspect of our method for determining the drive pressure is that it does not require physical modeling of the drive mechanism and can thus be broadly applied to many types of ramp-compression experiments, regardless of the drive mechanism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Laser ultrasonic imaging of subsurface defects with the linear sampling method

Laser ultrasonics is a remote nondestructive evaluation technique suitable for real-time monitoring of fabrication processes in semiconductor metrology, advanced manufacturing, and other applications where non-contact, high fidelity measurements are required. Here we investigate laser ultrasonic data processing approaches to reconstruct images of subsurface side drilled holes in aluminum alloy specimens. We demonstrate through simulation that the model-based linear sampling method (LSM) can perform accurate shape reconstruction of single and multiple holes and produce images with well-defined boundaries. We experimentally confirm that LSM produces images that represent the internal geometric features of an object, some of which may be missed by conventional imaging.

42 ENGINEERING↗

Improved Earthquake Source Parameters with 3D Wavespeed Models in California and Nevada

Seismic tomography harnesses earthquake data to explore the inaccessible structure of the Earth. Adjoint waveform tomography (AWT), a method of seismic tomography, updates the tomographic model by optimizing the fit between observed earthquake data and synthetic waveforms. The synthetic data are calculated by solving the wave equation through a given 3D model. An important requirement to calculating synthetics is the source information (location, centroid time, depth, and moment tensor). Errors in source information affect the quality of the synthetics produced, which in turn can limit how structure can be inferred in the AWT workflow. Here, to test the effect of updating source information, we used MTTime (Chiang, 2020), a time-domain full-waveform moment tensor inversion code, to calculate the moment tensors and depths of 118 earthquakes that occurred in California and Nevada over a 20-yr period. We calculated 3D Green’s functions using a 3D seismic wavespeed model of California and Nevada (Doody et al., 2023b). We show that the inverted solutions provide better waveform fits than the Global Centroid Moment Tensor catalog and increase usable, well-correlated data by up to 7%. Therefore, we argue that recalculating source parameters should be considered in AWT workflows, particularly for smaller magnitude events (⁠M w > 5.0).

58 GEOSCIENCES↗

Comparative Assessment of U-Net-Based Deep Learning Models for Segmenting Microfractures and Pore Spaces in Digital Rocks

Segmentation of high-resolution X-ray microcomputed tomography (µCT) images is crucial in digital rock physics (DRP), affecting the characterization and analysis of microscale phenomena in the porous media. The complexity of geological structures and nonideal scanning conditions pose significant challenges to conventional image segmentation approaches. Motivated by the recent increasing popularity of deep learning (DL) techniques in image processing, this work undertakes a comparative study of DL models, specifically U-Net and its variants, for segmenting multiple targets with distinguished features in digital rocks, including discrete fracture networks (DFNs), pore spaces, and solid rock. Particularly, DFNs have a smaller volumetric fraction over others, bringing in a substantial challenge of imbalanced segmentation. The primary focus is to evaluate the architecture and feature enhancement strategies of various DL models, including U-Net, attention U-Net, residual U-Net, U-Net++, and residual U-Net++. The models were designed as 2.5D, utilizing a central 2D image and its two adjacent upper and lower 2D images as input to provide a pseudo-3D context. In addition, because the ground truth of segmentation was unknown for real-world digital rocks, we created a benchmark data set following the inverse operations of segmentation. The data synthesis started from the label images (i.e., solid rock, pore spaces, and DFNs), followed by simulating partial volume blurring, adding random background noise, and introducing ring artifacts to mimic real raw X-ray µCT images. The data set, which included various rock types (i.e., sandstone and artificial data), scanning resolution, and magnitudes of noise and artifacts, was divided into training and testing data sets with a 90% and 10% ratio, respectively. Moreover, in addition to the conventional pixel-wise evaluation metrics, the physics-based metric of the lattice-Boltzmann method (LBM) simulated permeability provided more comprehensive assessments. The results demonstrated that the residual connections, nested architectures, and redesigned skip connections contribute to the model performance and give the residual U-Net++ the highest accuracy. The improvements were mainly on the boundaries and small targets, especially the DFNs, which dominate the interconnectivity and therefore affect the permeability greatly. This study also rigorously evaluated the efficiency and generalization of each model, demonstrating that the sophisticated architectures achieved excellent practicability and maintained robust performance on completely unseen data, ensuring their suitability for diverse and challenging DRP applications.

58 GEOSCIENCES↗

Investigation of acoustic waves under subsurface conditions to improve the predictions of rock mechanical properties and natural fracture characteristics

Mechanical properties and natural fracture characteristics are critical to investigate for subsurface engineering applications, including carbon storage, well drilling, and stimulation, as they govern rock stability, fluid flow, and mechanical behavior under stress. This dissertation integrates experimental and machine learning approaches to enhance the prediction and understanding of these properties by analyzing acoustic wave behavior under varied subsurface conditions. First, the influence of temperature, pore pressure, and supercritical CO2 (scCO2) saturation on poroelastic properties is examined using Gray Berea sandstone samples. The results show that temperature and pore pressure significantly affect the bulk modulus and Biot’s coefficient, while scCO2 saturation impacts rock compressibility, informing strategies for effective geological carbon storage. The study extends this understanding by experimentally evaluating the impact of reservoir depletion on the dynamic mechanical properties of the emerging Caney shale in South Oklahoma with the employment of unsupervised machine learning to predict static mechanical properties across the Caney shale. Integrating petrophysical data and chemostratigraphy, the workflow—featuring K-means clustering, principal component analysis (PCA), and inverse distance weighting (IDW)—improves stratigraphic characterization and the estimation of static-to-dynamic modulus ratios, which is vital for optimizing drilling and stimulation strategies. Finally, the work explores how natural fracture characteristics in shale influence acoustic waveforms and shear wave splitting (SWS) analysis. Experimental data on fractured samples under different stress and temperature conditions, combined with machine learning models such as K-nearest neighbors (KNN) and extreme gradient boosting (XGBoost), reveal key fracture properties impacting SWS and wave propagation. Together, these studies provide a comprehensive framework for linking acoustic wave behavior with rock properties, advancing the methods for monitoring and predicting geomechanical changes. The insights offered valuable implications for safer, more efficient CO2 injection, hydrocarbon extraction, and subsurface management.

Elkholy, Sherif↗

Single-Crystal Growth and Room-Temperature Magnetocaloric Effect of X-Type Hexaferrite Sr 2 Co 2 Fe 28 O 46

X-type hexaferrites have been receiving considerable attention due to their promising applications in many magnetic-electronic fields. However, the growth of single-crystal X-type hexaferrite is still a challenge. Herein we reported, for the first time, the preparation of single crystal X-type hexaferrite Sr 2 Co 2 Fe 28 O 46 (Sr 2 Co 2 X) with high-quality and large size using floating-zone method with laser as the heating source. The crystals show rhombohedral symmetry with space group of R-3m (No. 166, a = 5.8935(1) Å and c = 83.7438(17) Å). Co 2+ and Fe 3+ oxidation states were confirmed by the X-ray absorption near-edge spectroscopy. The prepared Sr 2 Co 2 X exhibits a spin reorientation transition from easy-cone to easy-axis at T 2 of 343 K and a ferrimagnetism–paramagnetism transition at Curie temperature (TC) of ~743 K. The spin reorientation transition was accompanied by magnetocaloric effect (MCE). Both conventional and inverse MCEs were observed near T 2 with a magnetic field applied along the c-axis. The maximum value of the magnetic entropy change along the c-axis was evaluated to be 1.1 J/kg·K for a magnetic field change of 5 T.

36 MATERIALS SCIENCE↗