Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model inversion”

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,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Configuration mixing in 28 Mg and the 26 Mg (t,p) 28 Mg reaction

We have studied the 26 Mg(t,p) 28 Mg reaction. The nucleus 28 Mg lies between 24 Mg and 32 Mg, the latter of which is central to the “island of inversion” characterized by low-lying neutron fp-shell configurations. The present results show that the second-excited 0 + state in 28 Mg contains large fp-shell occupation for the neutrons. The experiment was performed in inverse kinematics using the HELIcal Orbit Spectrometer (HELIOS) at Argonne National Laboratory. Shell-model calculations using the SDPF-MU interaction yielded level energies and two-neutron transfer amplitudes that were used in one-step distorted-wave Furthermore, born approximation calculations to provide theoretical predictions of the proton angular distributions. In many, but not all cases, the data were in good agreement with theoretical predictions.

20 ≤ A ≤ 38↗
NASA NTRS1995

Tests of two convection theories for red giant and red supergiant envelopes

Two theories of stellar envelope convection are considered here in the context of red giants and red supergiants of intermediate to high mass: Boehm-Vitense's standard mixing-length theory (MLT) and Canuto & Mazzitelli's new theory incorporating the full spectrum of turbulence (FST). Both theories assume incompressible convection. Two formulations of the convective mixing length are also evaluated: l proportional to the local pressure scale height (H(sub P)) and l proportional to the distance from the upper boundary of the convection zone (z). Applications to test both theories are made by calculating stellar evolutionary sequences into the red zone (z). Applications to test both theories are made by calculating stellar evolutionary sequences into the red phase of core helium burning. Since the theoretically predicted effective temperatures for cool stars are known to be sensitive to the assigned value of the mixing length, this quantity has been individually calibrated for each evolutionary sequence. The calibration is done in a composite Hertzsprung-Russell diagram for the red giant and red supergiant members of well-observed Galactic open clusters. The MLT model requires the constant of proportionality for the convective mixing length to vary by a small but statistically significant amount with stellar mass, whereas the FST model succeeds in all cases with the mixing lenghth simply set equal to z. The structure of the deep stellar interior, however, remains very nearly unaffected by the choices of convection theory and mixing lenghth. Inside the convective envelope itself, a density inversion always occurs, but is somewhat smaller for the convectively more efficient MLT model. On physical grounds the FST model is preferable, and seems to alleviate the problem of finding the proper mixing length.

Stothers, Richard B.↗
DOE OSTI2023

Machine learning-based real-time kinetic profile reconstruction in DIII-D

Abstract Kinetic equilibrium reconstruction plays a vital role in the physical analysis of plasma stability and control in fusion tokamaks. However, the traditional approach is subjective and prone to human biases. To address this, the consistent automatic kinetic equilibrium reconstruction (CAKE) method was introduced, providing objective results. Nonetheless, its offline nature limits its application in real-time plasma control systems (PCSs). To address this limitation, we present RTCAKENN, a machine learning model that approximates 7 CAKE-level output profiles, namely pressure, inverse q , toroidal current density, electron temperature and density, carbon ion impurity temperature and rotation profiles, using real-time available inputs. The deep neural network consists of an encoder layer, where the scalars and interdependent inputs such as plasma boundary coordinates and motional Stark effect data are encoded using multi-layer perceptrons (MLPs), while profile inputs are encoded by 1D convolutional layers. The encoded data is passed through a MLP for latent feature extraction, before being decoded in the decoding layers, which consist of upsampling and convolutional layers. RTCAKENN has been implemented in the DIII-D PCS and our model achieves accuracy comparable to CAKE and surpasses existing real-time alternatives. Through clever dropout training, RTCAKENN exhibits robustness and can operate even in the absence of Thomson scattering data or charge exchange recombination data. It executes in under 8 ms in the real-time environment, enabling future application in real-time control and analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRSDate not supplied

Accuracy of Center of Pressure Determination via Motion Capture

BACKGROUND: This study was conducted to support the stability assessment for tasks in lunar gravity and exercises on a Vibration Isolation and Stabilization (VIS) system in microgravity based on the dynamic feasibility criterion of whether the calculated position of the center of pressure (COP) falls within the base of support (BOS) which outlines the subject’s feet. Motion capture data combined with biomechanical modeling and simulation allows the forces and moments between the human and the VIS platform to be computed and the position of the COP as well as the location and shape of the BOS to be determined. The goal of this study was to assess the accuracy of the COP trajectory calculated using motion capture-based data. METHODS AND RESULTS: To obtain the dynamic quantities from which COP is calculated, motion capture data is first collected in the 1g lab environment by recording the trajectories of passive retroreflective markers placed on a subject during exercise or performance of a given task. The OpenSim [1] inverse kinematics (IK) tool is used to fit a scaled subject model to recorded marker trajectories while minimizing marker error to obtain joint angles. Then, a custom OpenSim plugin [2] is used to determine the subject’s time-varying moment of inertia and its time derivative, center of mass (CM) position, velocity, and acceleration, as well as the angular momentum and its time derivative relative to the subject’s CM. Some of these quantities are not needed for modeling tasks performed on a stationary lunar surface but, due to the moving exercise platform, are needed to model VIS response to the subject’s motion. Hand positions, used in calculating a cable force if present, are recorded as well. These quantities are used to calculate the total force (F ⃗^((plate) )) and moment (M ⃗^((plate) )) exerted by the lunar surface or the VIS plate on the subject’s shoe soles. COP is then calculated from the following equations: r_x^((cop) )= M_z^((plate) )/F_y^((plate) ) and r_z^((cop) )= 〖-M〗_x^((plate) )/F_y^((plate) ), where the y axis is normal to the surface. COP accuracy for feasibility assessments is then determined by whether it falls within the BOS, which is also computed by the plugin. To study the accuracy of COP calculated from motion capture, we first investigated whether COP remained within the BOS, as it must, for exercises performed in the 1g lab environment. Standard exercises such as back squat and deadlift were analyzed, as well as more explosive exercises including hang clean and press. Cases in which the COP exited the BOS indicated that COP accuracy required further investigation. In this study, an exercise device with cables was used, so cable force modeling accuracy should also be considered. In a separate study, we collected motion capture and force plate data for twenty-seven motions not involving an exercise device. About a third were genuine countermeasures exercises (e.g., hang clean and press), some were relevant for lunar tasks (e.g., object pick up), and the rest were of a “unit test” nature (e.g., swaying back and forth or side to side). Motion capture-based COP positions were compared with force plate measured COP. We found that while force plate measured COP remained within the BOS, motion capture-based COP was observed to briefly exit the BOS on occasion. Techniques to mitigate IK artifacts and filtering of calculated data could be used to improve the agreement of calculated and measured results, resolving excursions from the BOS within this dataset. The mean error between calculated and measured COP was found to be less than 6 mm. Additionally, we derived and investigated equations for the COP in terms of the cable force, cable location, as well as the human CM position, acceleration, and angular momentum with respect to the CM, and analyzed them for sensitivity to errors in individual quantities. Several were found, but the most significant one was that when the vertical force on the feet approaches zero, indicating a near-detachment or ‘jump off’ condition, errors are amplified. This is consistent with the observation that in the absence of pressure, the concept of the center of pressure would become meaningless.

C A Bell↗
DOE OSTI2022

Adversarial autoencoder ensemble for fast and probabilistic reconstructions of few-shot photon correlation functions for solid-state quantum emitters

Second-order photon correlation measurements [g (2) (τ) functions] are widely used to classify single-photon emission purity in quantum emitters or to measure the multiexciton quantum yield of emitters that can simultaneously host multiple excitations – such as quantum dots – by evaluating the value of g (2) (τ = 0). Accumulating enough photons to accurately calculate this value is time consuming and could be accelerated by fitting of few-shot photon correlations. Here, we develop an uncertainty-aware, deep adversarial autoencoder ensemble (AAE) that reconstructs noise-free g (2) (τ) functions from noise-dominated, few-shot inputs. The model is trained with simulated g (2) (τ) functions that are facilely generated by Poisson sampling time bins. The AAE reconstructions are performed orders-of-magnitude faster, with reconstruction errors and estimates of g (2) (τ = 0) that are lower in variance and similar in accuracy compared to Maximum likelihood estimation and Levenberg-Marquardt least-squares fitting approaches, for simulated and experimentally measured few-shot g (2) (τ) functions (~100 two-photon events) of InP/ZnS/ZnSe and CdS/CdSe/CdS quantum dots. The deep-ensemble model comprises eight individual autoencoders, allowing for probabilistic reconstructions of noise-free g (2) (τ) functions, and we show that the predicted variance scales inversely with number of shots, with comparable uncertainties to computationally intensive Markov chain Monte Carlo sampling. Furthermore, this work demonstrates the advantage of machine learning models to perform uncertainty-aware, fast, and accurate reconstructions of simple Poisson-distributed photon correlation functions, allowing for on-the-fly reconstructions and accelerated materials characterization of solid-state quantum emitters.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Deep Image Prior Enabled Full Waveform Inversion (Final Technical Report)

MS Student Naveen Gupta worked on the problem of full waveform inversion (FWI) using neural networks as shown in Figure 1. Our goal was to learn a neural network to represent the subsurface velocity model, which when fed into the FWI module (implemented using a numerical forward model of wave equations) produces amplitude estimates that match with ground-truth observations of amplitude. We used neural networks to solve the inverse problem of estimating velocity distributions for a given seismic amplitude data such that, once trained, our neural network model can generate a distribution of velocity profiles for different random vectors fed as inputs to the neural network model.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Morphogenic Growth 3D Printing

Inspired by nature's morphogenesis, a new 3D printing process –growth printing (GP)– takes advantage of a self‐propagating curing front to produce 3D polymeric parts following a growth‐like development plan. The propagation of the curing front is driven by the exothermic polymerization of dicyclopentadiene (DCPD), which transforms the liquid resin into a stiff polymer as it propagates at 1 mm s −1 . GP is triggered when a heated initiator contacts the uncured liquid resin in an open container. The initiator nucleates the frontal polymerization reaction and the isotropic radial propagation of the growth front. Simultaneously, the initiator is moved up across the free surface of the resin, pulling the cured object out of the uncured resin. The motion trajectory of the initiator with respect to the free resin surface controls the growth morphology of the 3D part. An inverse design algorithm is developed to produce 3D parts by modeling the reaction‐diffusion‐driven solidification process. This process has substantial energy savings and high printing speeds.

3D printing↗
DOE OSTI2024

Evolution of the Electronic and Excitonic Properties in 2D Ruddlesden–Popper Perovskites Induced by Bifunctional Ligands

2D Ruddlesden–Popper metal‐halide perovskites exhibit structural diversity due to a variety of choices of organic ligands. Incorporating bifunctional ligands in such materials is particularly intriguing since it can result in novel electronic properties and functions. However, an in‐depth understanding of the effects of bifunctional ligands on perovskite structures and, consequently, their electronic and excitonic properties, is still lacking. Here, in this work, n = 1 2D perovskites built with organic ligands containing ─CN, ─OH, ─COOH, ─phenyl (Ph), and ─CH 3 functional groups are investigated using ultraviolet and inverse photoemission spectroscopies, density functional theory calculations, and tight‐binding model analyses. The experimentally determined electronic gaps of the ─CN, ─COOH, ─Ph, and ─CH 3 based perovskites exhibit a strong correlation with the in‐plane Pb─I─Pb bond angle, while the ─OH based perovskite deviates from the linear trend. Based on the band structure calculations, this anomaly is attributed to the out‐of‐plane dispersion, caused predominantly by significant interlayer electronic coupling that is present in ─OH based perovskites. These results highlight the complex and diverse impacts of organic ligands on electronic properties, especially in terms of the involvement of strong interlayer electronic coupling. The impact of the bifunctional ligands on the evolution of the exciton binding energy is also addressed.

2D perovskites↗
DOE OSTI2021

An adaptive Hessian approximated stochastic gradient MCMC method

Bayesian approaches have been successfully integrated into training deep neural networks. One popular family is stochastic gradient Markov chain Monte Carlo methods (SG-MCMC), which have gained increasing interest due to their ability to handle large datasets and the potential to avoid overfitting. Although standard SG-MCMC methods have shown great performance in a variety of problems, they may be inefficient when the random variables in the target posterior densities have scale differences or are highly correlated. Here, we present an adaptive Hessian approximated stochastic gradient MCMC method to incorporate local geometric information while sampling from the posterior. The idea is to apply stochastic approximation (SA) to sequentially update a preconditioning matrix at each iteration. The preconditioner possesses second-order information and can guide the random walk of a sampler efficiently. Instead of computing and saving the full Hessian of the log posterior, we use limited memory of the samples and their stochastic gradients to approximate the inverse Hessian-vector multiplication in the updating formula. Moreover, by smoothly optimizing the preconditioning matrix via SA, our proposed algorithm can asymptotically converge to the target distribution with a controllable bias under mild conditions. To reduce the training and testing computational burden, we adopt a magnitude-based weight pruning method to enforce the sparsity of the network. Our method is user-friendly and demonstrates better learning results compared to standard SG-MCMC updating rules. The approximation of inverse Hessian alleviates storage and computational complexities for large dimensional models. Numerical experiments are performed on several problems, including sampling from 2D correlated distribution, synthetic regression problems, and learning the numerical solutions of heterogeneous elliptic PDE. The numerical results demonstrate great improvement in both the convergence rate and accuracy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Multi-energy reconstructions, central electron temperature measurements, and early detection of the birth and growth of runaway electrons using a versatile soft x-ray pinhole camera at MST

A multi-energy soft x-ray pinhole camera has been designed, built, and deployed at the Madison Symmetric Torus to aid the study of particle and thermal transport, as well as MHD stability physics. This novel imaging diagnostic technique employs a pixelated x-ray detector in which the lower energy threshold for photon detection can be adjusted independently on each pixel. The detector of choice is a PILATUS3 100 K with a 450 μm thick silicon sensor and nearly 100 000 pixels sensitive to photon energies between 1.6 and 30 keV. An ensemble of cubic spline smoothing functions has been applied to the line-integrated data for each time-frame and energy-range, obtaining a reduced standard-deviation when compared to that dominated by photon-noise. The multi-energy local emissivity profiles are obtained from a 1D matrix-based Abel-inversion procedure. Central values of Te can be obtained by modeling the slope of the continuum radiation from ratios of the inverted radial emissivity profiles over multiple energy ranges with no a priori assumptions of plasma profiles, magnetic field reconstruction constraints, high-density limitations, or need of shot-to-shot reproducibility. In tokamak plasmas, a novel application has recently been tested for early detection, 1D imaging, and study of the birth, exponential growth, and saturation of runaway electrons at energies comparable to 100 × Te,0; thus, early results are also presented.

Delgado-Aparicio, L. F. (ORCID:0000000167394380)↗
DOE OSTI2025

Symmetry instability induced by topological phase transitions

The symmetry-topology interplay dictates how to define order parameters and classify material ordered phases. However, current understanding of this interplay has been predominately approached from a one-sided perspective, with topological states being classified within the constraints imposed by specific fixed symmetries. Here, in this study, we complete this full circle by demonstrating spontaneous symmetry breaking that results from a periodic alteration of topological phases induced by light in a centrosymmetric Dirac material ZrTe 5 . The distinguishing feature is the observation of robust correlation and striking anomalies in the fluence and temperature dependence of key transport parameters. First, both shift current 𝐽 𝑠 and displacement current 𝐽 𝑑 , arising from interband transition and infrared phonon driving, respectively, along with charge carrier pumping, exhibit similar behaviors. Second, they all peak at similar low pump fluence, followed by a subsequent reduction as the fluence further increases. This behavior cannot be explained by conventional energetically allowed, direct excitations. Third, all the three observables exhibit anomalies when they approach the topological phase transition temperature. These results highlight the unique low-energy pumping behaviors in ZrTe 5 , characterized by reversible fluence dependence and a “hinge-like” interaction that connects various electronic and lattice observables, including phonons, charge carriers, and currents. Our findings, supported by model analysis, provide key insights into the fragility of crystalline (inversion) and time-reversal symmetries during the dynamics of topological phase transitions. This fragility drives spontaneous symmetry breaking, evidenced by the synchronized emergence of off-resonant infrared phonons and broken-symmetry photocurrents.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

OASIS: Offsetting Active Reconstruction Attacks in Federated Learning

Federated Learning (FL) has garnered significant attention for its potential to protect user privacy while enhancing model training efficiency. For that reason, FL has found its use in various domains, from health care to industrial engineering, especially where data cannot be easily exchanged due to sensitive information or privacy laws. However, recent research has demonstrated that FL protocols can be easily compromised by active reconstruction attacks executed by dishonest servers. These attacks involve the malicious modification of global model parameters, allowing the server to obtain a verbatim copy of users' private data by inverting their gradient updates. Tackling this class of attack remains a crucial challenge due to the strong threat model. In this paper, we propose a defense mechanism, namely OASIS, based on image augmentation that effectively counteracts active reconstruction attacks while preserving model performance. We first uncover the core principle of gradient inversion that enables these attacks and theoretically identify the main conditions by which the defense can be robust regardless of the attack strategies. We then construct our defense with image augmentation showing that it can undermine the attack principle. Comprehensive evaluations demonstrate the efficacy of the defense mechanism highlighting its feasibility as a solution.

deep neural networks↗
NASA NTRS1980

Fundamental studies in geodynamics

Progress in modeling instantaneous plate kinematics is reviewed, with emphasis on recently developed models of present day plate motions derived by the systematic inversion of globally distributed data sets. Rivera plate motions, the Caribbean South American boundary, Indian plate deformation, Pacific-North America, seismicity and subduction processes, and the study of slow earthquakes and free oscillations are discussed.

Anderson, D. L.↗
NASA NTRS1982

Emission mechanisms of the March 5th 1979 gamma ray transient

This paper reviews the spectral and temporal constraints on the emission mechanisms of the initial burst and subsequent pulses of the March 5, 1979 gamma-transient. Interpretation of the hard tail above 300 keV as due to inverse comptonization, when combined with the synchrotron-pair annihilation model of Ramaty et al. (1981), makes it possible to derive from first principles the intrinsic luminosity of the source. This is in general agreement with that expected from N49 in LMC. The pulsed emissions have characteristics difficult to explain in terms of conventional models.

Liang, E. P. T.↗
NASA NTRS1986

An analysis of the role of drift waves in equatorial spread F

An account is given of results of rocket measurements of the wave number spectrum of equatorial spread F irregularities, with emphasis on wavelengths less than 100 m. The measurements were made from two sounding rockets launched from Peru as part of Project Condor. The Condor density fluctuation spectra display a break at a wavelength near 100 m, identical to that found in the PLUMEX experiment (Kelley et al., 1982). The Condor data also confirm a subrange in which the density and the wave potential obey the Boltzmann relation - a strong indication of the presence of low-frequency electrostatic waves with finite wavelength parallel to the magnetic field, perhaps low-frequency drift waves as proposed by Kelley et al. The Condor data are also consistent with the previous conjecture that drift waves only exist above 300 km altitude. To investigate the difference in spectra observed over two altitude ranges, the data must be fitted to a form for the power spectrum taken from Keskinen and Ossakow (1981). The fitted spectrum, along with empirically determined growth and dissipation rates, is used to calculate the energy pumped into the spectrum at long wavelengths as well as the energy dissipated at shorter wavelengths. It is found that the energy is balanced by classical collisional effects in the low-altitude case, but energy balance in the high-altitude case requires an enhanced dissipation of about 500 times that due to classical diffusion. The model is consistent with, but does not uniquely imply, an inverse cascade of drift wave turbulence in equatorial spread F.

Labelle, J.↗
NASA NTRS1988

Estimation of canopy parameters for inhomogeneous vegetation canopies from reflectance data. III - TRIM: A model for radiative transfer in heterogeneous three-dimensional canopies

A model for radiative transfer in heterogeneous three-dimensional canopies such as those found in forests is proposed. Its use in estimating important biophysical variables such as leaf area index and canopy architecture from bidirectional canopy reflectance data is discussed. The model and its use in estimating canopy parameters through its inversion are validated with measured canopy reflectance data for corn canopies.

Goel, Narendra S.↗
NASA NTRS1989

An observational study of stratocumulus entrainment and thermodynamics

The stratocumulus (SC) entrainment and thermodynamics are studied using data from the Dynamics and Chemistry of Marine Stratocumulus (DYCOMS) experiment. The rate of entrainment of air from the free troposphere into the cloud-topped PBL is estimated using a technique based on the measurement of ozone flux and mean distribution. The average measured value of entrainment rate was found to be 3.0 mm/s with a range of 1.0 to 5.0 mm/s for cloudy cases. Thermodynamic budgets are constructed for eight DYCOMS cases. It was found that the divergence of the solar radiative flux is an important component of the boundary-layer energetics during midday, the factor which must be accounted for in mixed-layer models. The net longwave radiative flux profiles show good agreement with theoretical models, but an unambiguous partitioning of the flux divergence between inversion and mixed layers could not be established.

Kawa, S. R.↗