Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “inverse modeling”

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 1,171 records · Page 65

Inverse Bremsstrahlung Absorption

Inverse bremsstrahlung absorption was measured based on transmission through a finite-length plasma that was thoroughly characterized using spatially resolved Thomson scattering. Expected absorption was then calculated using the diagnosed plasma conditions while varying the absorption model components. To match data, it is necessary to account for: (1) the Langdon effect; (2) laser-frequency (rather than plasma-frequency) dependence in the Coulomb logarithm, as is typical of bremsstrahlung theories but not transport theories; and (3) a correction due to ion screening. Radiation-hydrodynamic simulations of inertial confinement fusion implosions have to date used a Coulomb logarithm from the transport literature and no screening correction. Here, we anticipate that updating the model for collisional absorption will substantially revise our understanding of laser-target coupling for such implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Enhanced Planar Antenna Efficiency Through Magnetic Thin-Films

This work proposes to use magnetic material as the substrate of planar antennas to overcome the platform effect caused by the conducting ground plane. The upper bound of the radiation efficiency of an electric-current-driven low-profile antenna is theoretically derived, which is inversely proportional to the Gilbert damping factor of the magnetic material. Meanwhile, the improvement of radiation due to the use of magnetic material is demonstrated by a three-dimensional (3D) multiphysics and multiscale time-domain model. The simulation results match the theoretical derivation, showing 25% radiation efficiency from a planar antenna backed by a FeGaB thin film with 2.56 m thickness. Furthermore, for conductive ferromagnetic materials, it is shown that the eddy current loss can be well suppressed by laminating the thin film into multiple layers. Here, the radiation efficiency of the modeled antenna with a conductive ferromagnetic substrate is improved from 2.2% to 11.8% by dividing the substrate into 10 layers, with a ferromagnetic material fill factor of 93%.

42 ENGINEERING↗

The Effect of Spatially Varying Collision Frequency on the Development of the Rayleigh–Taylor Instability

Abstract The Rayleigh–Taylor (RT) instability is ubiquitously observed, yet has traditionally been studied using ideal fluid models. Collisionality can vary strongly across the fluid interface, and previous work demonstrates the necessity of kinetic models to completely capture dynamics in certain collisional regimes. Where previous kinetic simulations used spatially and temporally constant collision frequency, this work presents five-dimensional (two spatial, three velocity dimensions) continuum-kinetic simulations of the RT instability using a more realistic spatially varying collision frequency. Three cases of collisional variation are explored for two Atwood numbers: low to intermediate, intermediate to high, and low to high. The low-to-intermediate case exhibits no RT instability growth, while the intermediate-to-high case is similar to a fluid-limit kinetic case with interface widening biased toward the lower-collisionality region. A novel contribution of this work is the low-to-high collisionality case that shows significantly altered instability growth through an upward movement of the interface and damped spike growth due to increased free-streaming particle diffusion in the lower region. Contributions to the energy flux from the non-Maxwellian portions of the distribution function are not accessible to fluid models and are greatest in magnitude in the spike and regions of low collisionality. Increasing the Atwood number results in greater RT instability growth and reduced upward interface movement. Deviation of the distribution function from Maxwellian is inversely proportional to collision frequency and concentrated around the fluid interface. The linear phase of RT instability growth is well described by theoretical linear growth rates accounting for viscosity and diffusion.

79 ASTRONOMY AND ASTROPHYSICS↗

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↗

Spectrum and extension of the inverse-Compton emission of the Crab Nebula from a combined Fermi -LAT and H.E.S.S. analysis

The Crab Nebula is a unique laboratory for studying the acceleration of electrons and positrons through their non-thermal radiation. Observations of very-high-energy γ rays from the Crab Nebula have provided important constraints for modelling its broadband emission. We present the first fully self-consistent analysis of the Crab Nebula’s γ-ray emission between 1 GeV and ∼100 TeV, that is, over five orders of magnitude in energy. Using the open-source software package GAMMAPY, we combined 11.4 yr of data from the Fermi Large Area Telescope and 80 h of High Energy Stereoscopic System (H.E.S.S.) data at the event level and provide a measurement of the spatial extension of the nebula and its energy spectrum. We find evidence for a shrinking of the nebula with increasing γ-ray energy. Furthermore, we fitted several phenomenological models to the measured data, finding that none of them can fully describe the spatial extension and the spectral energy distribution at the same time. Especially the extension measured at TeV energies appears too large when compared to the X-ray emission. Our measurements probe the structure of the magnetic field between the pulsar wind termination shock and the dust torus, and we conclude that the magnetic field strength decreases with increasing distance from the pulsar. We complement our study with a careful assessment of systematic uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗

A machine learning approach to determine the elastic properties of printed fiber-reinforced polymers

This work focuses on the simultaneous determination of the elastic constants and the fiber orientation state for a short fiber-reinforced polymer composite by performing a minimum of experimental tests. Here we introduce a methodology that enables the inverse determination of fiber orientation state and the in-situ polymer properties by performing tensile tests at the composite coupon level. We demonstrate the approach for the extrusion deposition additive manufacturing (EDAM) process to illustrate one application of the methodology, but the development is such that it can be applied to short fiber-reinforced polymer (SFRP) systems processed via other methods. Currently, developing composites additive manufacturing digital twins require extensive material characterization. In particular, the mechanical characterization of the orthotropic elastic properties of a composite involves extensive sample preparation and testing, therefore the elasticity tensor is generally populated using a micromechanics model. This, however, requires measuring the fiber orientation state in addition to knowing the constituent material properties. Experimentally measuring the fiber orientation state can be tedious and time consuming. Further, optical methods are limited to resolving the orientation of cylindrical fibers or cluster of non-cylindrical fibers, and computed tomography (CT) methods scan regions of volume that are much smaller than a full printed bead. Therefore, we propose a methodology, accelerated by machine learning, to identify the anisotropic mechanical properties and fiber orientation state at the same time. Early results show that inference of the fiber orientation and composite properties is possible with as few as three tensile tests. Our results show that a combination of the choice of the micromechanics model and reliable set of experiments can yield the nine elastic constants, as well as, the fiber orientation state.

36 MATERIALS SCIENCE↗

A Novel Method for Diagnosing Land–Atmosphere Coupling Sensitivity in a Single-Column Model

The response of boundary layer properties and cloudiness to changes in surface evaporative fraction (EF) is investigated in a single-column model to quantify the locally coupled impact of subgrid surface variations on the atmosphere during summer. Sensitive coupling days are defined when the model atmosphere exhibits large variations across a range of EFs centered on the analyzed value. Coupling sensitivity exists as both positive feedback (cloudiness increases with EF) and negative feedback (clouds increase with decreasing EF) regimes. The positive regime manifests in shallow convection situations, which are capped by a strengthened inversion and subsidence, restricting the vertical extent of convection to just above the boundary layer. Surfaces with larger EF (greater surface latent heat flux) can inject more moisture into the vertically confined system, lowering the cloud base and an increasing cloud liquid water path (LWP). Negative feedback regimes tend to manifest when large-scale deep convection, such as from mesoscale convective systems and fronts, is advected through the domain, where convection strengthens over surfaces with a lower EF (greater surface sensible heat flux). The invigoration of these systems by the land surface leads to an increase in LWP through strengthened updrafts and stronger coupling between the boundary layer and the free atmosphere. These results apply in the absence of heterogeneity-induced mesoscale circulations, providing a one-dimensional dynamical perspective on the effect of surface heterogeneity. Finally, this study provides a framework of intermediate complexity, lying between parcel theory and high-resolution coupled land–atmosphere modeling, and therefore isolates the relevant first-order processes in land–atmosphere interactions.

54 ENVIRONMENTAL SCIENCES↗

B-PINNs: Bayesian Physics-informal Neural Networks for Forward and Inverse PDE Problems with Noisy Data

We propose a Bayesian physics-informed neural network (B-PINN) to solve both forward and inverse nonlinear problems described by partial differential equations (PDEs) and noisy data. In this Bayesian framework, the Bayesian neural network (BNN) combined with a PINN for PDEs serves as the prior while the Hamiltonian Monte Carlo (HMC) or the variational inference (VI) could serve as an estimator of the posterior. B-PINNs make use of both physical laws and scattered noisy measurements to provide predictions and quantify the aleatoric uncertainty arising from the noisy data in the Bayesian framework. Compared with PINNs, in addition to uncertainty quantification, B-PINNs obtain more accurate predictions in scenarios with large noise due to their capability of avoiding overfitting. We conduct a systematic comparison between the two different approaches for the B-PINNs posterior estimation (i.e., HMC or VI), along with dropout used for quantifying uncertainty in deep neural networks. Our experiments show that HMC is more suitable than VI with mean field Gaussian approximation for the B-PINNs posterior estimation, while dropout employed in PINNs can hardly provide accurate predictions with reasonable uncertainty. Finally, we replace the BNN in the prior with a truncated Karhunen-Loève (KL) expansion combined with HMC or a deep normalizing flow (DNF) model as posterior estimators. The KL is as accurate as BNN and much faster but this framework cannot be easily extended to high-dimensional problems unlike the BNN based framework.

Non linear PDEs, Noisy data, Bysian physics inform↗

Physics-informed machine learning

Despite great progress in simulating multiphysics problems using the numerical discretization of partial differential equations (PDEs), one still cannot seamlessly incorporate noisy data into existing algorithms, mesh generation remains complex, and high-dimensional problems governed by parameterized PDEs cannot be tackled. Moreover, solving inverse problems with hidden physics is often prohibitively expensive and requires different formulations and elaborate computer codes. Machine learning has emerged as a promising alternative, but training deep neural networks requires big data, not always available for scientific problems. Instead, such networks can be trained from additional information obtained by enforcing the physical laws (for example, at random points in the continuous space-time domain). Such physics-informed learning integrates (noisy) data and mathematical models, and implements them through neural networks or other kernel-based regression networks. Moreover, it may be possible to design specialized network architectures that automatically satisfy some of the physical invariants for better accuracy, faster training and improved generalization. Furthermore, we review some of the prevailing trends in embedding physics into machine learning, present some of the current capabilities and limitations and discuss diverse applications of physics-informed learning both for forward and inverse problems, including discovering hidden physics and tackling high-dimensional problems.

97 MATHEMATICS AND COMPUTING↗

Advancing electrochemical impedance analysis through innovations in the distribution of relaxation times method

Electrochemical impedance spectroscopy (EIS) is a key tool across various scientific disciplines, including energy sciences, chemistry, and biology, enabling the analysis of electrochemical systems. However, conventional methods for interpreting EIS data are often complex and model dependent. The distribution of relaxation times (DRT) offers a non-parametric approach that simplifies the interpretation process by providing a timescale interpretation of EIS data. This article provides a comprehensive review of current methods for DRT inversion. Additionally, a survey of practitioners highlights key challenges in the field. Here, the findings underscore the need for standardized DRT analysis and benchmarks, as well as the development of automated analysis tools. These advancements would improve the usability and interpretability of EIS data. Ultimately, implementing these improvements could not only propel the field forward but also expand the application of DRT in scientific research by making it accessible to a broader range of researchers, including those without specialized expertise in programming or statistics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molecular-Weight-Dependent Infiltration and Adsorption of Polymers into Nanochannels

The mesoporous silica shell coating hydrogenolysis nano-catalysts alters the molecular weight distributions of cleaved polymer chains compared to catalysts without a shell. The shell, composed of radially aligned narrow cylindrical nanopores, reduces the formation of low-valued gaseous products and increases the median molecular weight of the product, thus enhancing the value of the products for polymer upcycling. Here, to understand the role of the mesoporous shell, we have studied the spatial distribution of polystyrene chains, used as a model polymer, in the nanochannels in both the melt phase and solution phase. In the melt, we observed from small-angle X-ray scattering experiments that the infiltration rate of the polymer into the nanochannels is inversely proportional to the molecular weight, which is consistent with theory. In theta solution experiments using UV–vis spectroscopy, we found that the shell significantly enhances polymer adsorption compared to nanoparticles without pores. In addition, the degree of polymer adsorption is not a monotonic function of molecular weight but initially increases with the molecular weight before eventually decreasing. The molecular weight for the peak adsorption increases with the pore diameter. This adsorption behavior is rationalized as resulting from a balance between the mixing entropy gain by surface adsorption and the conformational entropy penalty incurred by chains confined in the nanochannels. The spatial distribution of polymer chains in the nanochannels is visualized by energy-dispersive X-ray spectroscopy (EDX), and inverse Abel-transformed data reveals a less uniform polymer distribution along the primary pore axis for longer chains.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nernst–Ettingshausen effect in thin Pt and W films at low temperatures

As spin caloritronic measurements become increasingly common techniques for characterizing material properties, it is important to quantify potentially confounding effects. We report measurements of the Nernst–Ettingshausen response from room temperature to 5 K in thin film wires of Pt and W, metals commonly used as inverse spin Hall detectors in spin Seebeck characterization. Johnson–Nyquist noise thermometry is used to assess the temperature change in the metals with heater power at low temperatures, and the thermal path is analyzed via finite-element modeling. The Nernst–Ettingshausen response of W is found to be approximately temperature-independent, while the response of Pt increases at low temperatures. Furthermore, these results are discussed in the context of theoretical expectations and the possible role of magnetic impurities in Pt.

36 MATERIALS SCIENCE↗

Non-perturbative production of fermionic dark matter from fast preheating

We investigate non-perturbative production of fermionic dark matter in the early universe. We study analytically the gravitational production mechanism accompanied by the coupling of fermions to the background inflaton field. The latter leads to the variation of effective fermion mass during preheating and makes the resulting spectrum and abundance sensitive to its parameters. Assuming fast preheating that completes in less than the inflationary Hubble time and no oscillations of the inflaton field after inflation, we find an abundant production of particles with energies ranging from the inflationary Hubble rate to the inverse duration of preheating. The produced fermions can account for all observed dark matter in a broad range of parameters. As an application of our analysis, we study non-perturbative production of fermionic dark matter in the model of Palatini Higgs inflation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Observation of the near-threshold intruder 0 – resonance in 12 Be

A resonant state at $3.21^{+0.12}_{–0.04}$ MeV, located just above the one-neutron separation threshold, was observed for the first time in 12 Be from the 11 Be(d,p) 12 Be one-neutron transfer reaction in inverse kinematics. This state is assigned a spin-parity of 0 – according to the systematics of the level scheme of the N=8 isotones and decay-width analysis. Gamow coupled-channel and Gamow shell-model calculations show the importance of the continuum coupling, which dramatically influences the excitation energy and ordering of low-lying states. Various exotic structures associated with cross-shell intruding configurations in 12 Be and in its isotonic nucleus 11 Li are comparably discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Narrow Slot Gasket Grooves And Step Changes In Width

We provide corrections to the slot capacitance and inverse inductance per unit length for slot gasket groove geometries using an approximate conformal mapping approach. We also provide corrections for abrupt step changes in slot width along with boundary discontinuity conditions for implementation in the various slot models.

42 ENGINEERING↗

Numerical solution of large scale Hartree–Fock–Bogoliubov equations

The Hartree–Fock–Bogoliubov (HFB) theory is the starting point for treating superconducting systems. However, the computational cost for solving large scale HFB equations can be much larger than that of the Hartree–Fock equations, particularly when the Hamiltonian matrix is sparse, and the number of electrons N is relatively small compared to the matrix size N b . We first provide a concise and relatively self-contained review of the HFB theory for general finite sized quantum systems, with special focus on the treatment of spin symmetries from a linear algebra perspective. We then demonstrate that the pole expansion and selected inversion (PEXSI) method can be particularly well suited for solving large scale HFB equations. For a Hubbard-type Hamiltonian, the cost of PEXSI is at most $\mathcal{O}$(N b 2 ) for both gapped and gapless systems, which can be significantly faster than the standard cubic scaling diagonalization methods. We show that PEXSI can solve a two-dimensional Hubbard-Hofstadter model with N b up to 2.88 × 10 6 , and the wall clock time is less than 100 s using 17 280 CPU cores. Finally, this enables the simulation of physical systems under experimentally realizable magnetic fields, which cannot be otherwise simulated with smaller systems.

97 MATHEMATICS AND COMPUTING↗

High Spatial Melt Rate Variability Near the Totten Glacier Grounding Zone Explained by New Bathymetry Inversion

Totten Glacier is a fast-moving East Antarctic outlet with the potential for significant future sea-level contributions. We deployed four autonomous phase-sensitive radars on its ice shelf to monitor ice-ocean interactions near its grounding zone and made active source seismic observations to constrain gravity-derived bathymetry models. We observe an asymmetry in basal melting with mean melt rates along the grounding zone differing by up to 20 m/a. Our new bathymetry model reveals that this melt rate asymmetry coincides with an asymmetry in water column thickness and that the low-melting ice-shelf portion is shielded from the main cavity circulation. A 2-year record yields year-to-year melt rate variability of 7–9 m/a with no seasonal cycle. Our results highlight the key role of bathymetry near grounding lines for accurate modeling of ice-shelf melt, and the importance of sustained multi-year monitoring, especially at ice-shelf cavities where the dominant melt rate drivers vary primarily inter-annually.

54 ENVIRONMENTAL SCIENCES↗

Non-volatile electric-field control of inversion symmetry

Competition between ground states at phase boundaries can lead to significant changes in properties under stimuli, particularly when these ground states have different crystal symmetries. A key challenge is to stabilize and control the coexistence of symmetry-distinct phases. Using BiFeO 3 layers confined between layers of dielectric TbScO 3 as a model system, we stabilize the mixed-phase coexistence of centrosymmetric and non-centrosymmetric BiFeO 3 phases at room temperature with antipolar, insulating and polar semiconducting behaviour, respectively. Application of orthogonal in-plane electric (polar) fields results in reversible non-volatile interconversion between the two phases, hence removing and introducing centrosymmetry. Counterintuitively, we find that an electric field 'erases' polarization, resulting from the anisotropy in octahedral tilts introduced by the interweaving TbScO 3 layers. Consequently, this interconversion between centrosymmetric and non-centrosymmetric phases generates changes in the non-linear optical response of over three orders of magnitude, resistivity of over five orders of magnitude and control of microscopic polar order. Here our work establishes a platform for cross-functional devices that take advantage of changes in optical, electrical and ferroic responses, and demonstrates octahedral tilts as an important order parameter in materials interface design.

36 MATERIALS SCIENCE↗