Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Inverse Optimization”

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 865 records · Page 48

Real-time neuromorphic algorithms for inverse kinematics of redundant manipulators

The paper presents an efficient neuromorphic formulation to accurately solve the inverse kinematics problem for redundant manipulators. The approach involves a dynamical learning procedure based on a novel formalism in neural network theory: the concept of 'terminal' attractors. Topographically mapped terminal attractors are used to define a neural network whose synaptic elements can rapidly encapture the inverse kinematics transformations, and, subsequently generalize to compute joint-space coordinates required to achieve arbitrary end-effector configurations. Unlike prior neuromorphic implementations, this technique can also systematically exploit redundancy to optimize kinematic criteria, e.g., torque optimization. Simulations on 3-DOF and 7-DOF redundant manipulators, are used to validate the theoretical framework and illustrate its computational efficacy.

Barhen, Jacob↗

Plasma puff initiation of high Coulomb transfer switches

The plasma-puff triggering mechanism based on a hypocycloidal pinch geometry was investigated to determine the optimal operating conditions for the azimuthally uniform surface flashover which initiates plasma-puff under wide range of fill gas pressure of Ar, He and N2. The optimal fill gas pressure for the azimuthally uniform plasma-puff was about 120 mTorr and 450 Torr for He and N2, and between 120 mTorr and 5 Torr for Ar. The inverse pinch switch was triggered with the plasma-puff and the switching capability under various electrical parameters and working gas pressures of Ar, He and N2 was determined. It was also shown that the azimuthally uniform switching discharges were dependent on the type of fill gas and its fill pressure. A new concept of plasma-focus driven plasma-puff was also discussed in comparison with the hypocycloidal pinch plasma-puff triggering. The main discharge of inverse pinch switch with plasma-focus driven plasma-puff trigger is found to be more azimuthally uniform than that with hypocycloidal pinch plasma-puff trigger in a gas pressure region between 80 mTorr and 1 Torr.

Venable, D. D.↗

Improved Bayesian regularization of inverse problems in vibrations and acoustics using noise-only measurements

Here, this paper studies Tikhonov regularization (ridge regression) parameter selection for problems in vibrations and acoustics. The selection method is based on a popular Bayesian method, but it incorporates measurements of sensor noise. The regularization parameter is closely related to the ratio of system input energy to noise energy, so noise measurements inform the inference procedure and improve parameter identification. In cases where standard Bayesian regularization identifies zero as the optimal regularization parameter, noise measurements guarantee a unique nonzero optimum. Sufficient theoretical criteria are developed for this guarantee. The method is verified in even-determined and under-determined configurations in an acoustic source localization simulation and a vibration load identification experiment. It is shown to yield significant improvements over existing empirical Bayesian regularization. Improvements are larger in the even-determined case and smaller in the under-determined case, wherein the inverse solution is less sensitive to the regularization parameter.

42 ENGINEERING↗

Machine learning inversion of interatomic force constants from single-crystal inelastic neutron scattering

Atomic vibrations govern many macroscopic properties of materials, but experiments to comprehensively probe them remain challenging. Inelastic neutron scattering (INS) is a powerful technique to map phonon dispersions in crystals, especially when leveraging modern time-of-flight (ToF) spectrometers with large detectors. However, efficiently and robustly extracting interatomic force constants (FCs) parameterizing phonon dynamics from experimental spectra remains a bottleneck due to the complexity and high dimensionality of ToF INS datasets. Here, we present a machine learning approach for the direct inversion of FCs from single-crystal INS measurements. The framework leverages synthetic training data generated using universal machine-learned force fields and an efficient physics-based forward model. We benchmark two neural architectures–one emphasizing structured latent representation learning and the other direct, supervised spectral regression–across simulated datasets for two materials under idealized and noisy conditions. The latent-representation model is subsequently applied to experimental single-crystal INS data on germanium. The model is shown to reproduce FCs derived from both first-principles simulations and from iterative optimization, and furthermore achieves reliable inference even from sparse, single-orientation measurements representing short data acquisitions. Analysis of the learned latent space reveals semantically continuous and physically interpretable encodings that support strong cross-domain generalization. By bridging theoretical and experimental domains, we establish a path toward rapid inversion of experimental spectra and data-driven interpretation of temperature-dependent lattice dynamics.

42 ENGINEERING↗

A Spitzer five-band analysis of the Jupiter-sized planet TrES-1

With an equilibrium temperature of 1200 K, TrES-1 is one of the coolest hot Jupiters observed by Spitzer. It was also the first planet discovered by any transit survey and one of the first exoplanets from which thermal emission was directly observed. We analyzed all Spitzer eclipse and transit data for TrES-1 and obtained its eclipse depths and brightness temperatures in the 3.6 μm (0.083% ± 0.024%, 1270 ± 110 K), 4.5 μm (0.094% ± 0.024%, 1126 ± 90 K), 5.8 μm (0.162% ± 0.042%, 1205 ± 130 K), 8.0 μm (0.213% ± 0.042%, 1190 ± 130 K), and 16 μm (0.33% ± 0.12%, 1270 ± 310 K) bands. The eclipse depths can be explained, within 1σ errors, by a standard atmospheric model with solar abundance composition in chemical equilibrium, with or without a thermal inversion. The combined analysis of the transit, eclipse, and radial-velocity ephemerides gives an eccentricity of e=0.033{sub −0.031}{sup +0.015}, consistent with a circular orbit. Since TrES-1's eclipses have low signal-to-noise ratios, we implemented optimal photometry and differential-evolution Markov Chain Monte Carlo (MCMC) algorithms in our Photometry for Orbits, Eclipses, and Transits pipeline. Benefits include higher photometric precision and ∼10 times faster MCMC convergence, with better exploration of the phase space and no manual parameter tuning.

79 ASTRONOMY AND ASTROPHYSICS↗

Force-Free Identification of Minimum-Energy Pathways and Transition States for Stochastic Electronic Structure Theories

Here, the accurate mapping of potential energy surfaces (PESs) is crucial to our understanding of the numerous physical and chemical processes mediated by atomic rearrangements, such as conformational changes and chemical reactions, and the thermodynamic and kinetic feasibility of these processes. Stochastic electronic structure theories, e.g., Quantum Monte Carlo (QMC) methods, enable highly accurate total energy calculations that in principle can be used to construct the PES. However, their stochastic nature poses a challenge to the computation and use of forces and Hessians, which are typically required in algorithms for minimum-energy pathway (MEP) and transition state (TS) identification, such as the nudged elastic band (NEB) algorithm and its climbing image formulation. Here, we present strategies that utilize the surrogate Hessian line-search method, previously developed for QMC structural optimization, to efficiently identify MEP and TS structures without requiring force calculations at the level of the stochastic electronic structure theory. By modifying the surrogate Hessian algorithm to operate in path-orthogonal subspaces and at saddle points, we show that it is possible to identify MEPs and TSs by using a force-free QMC approach. We demonstrate these strategies via two examples, the inversion of the ammonia (NH 3 ) molecule and the nucleophilic substitution (S N 2) reaction F – + CH 3 F → FCH 3 + F – . We validate our results using Density Functional Theory (DFT)- and Coupled Cluster (CCSD, CCSD(T))-based NEB calculations. We then introduce a hybrid DFT-QMC approach to compute thermodynamic and kinetic quantities, free energy differences, rate constants, and equilibrium constants that incorporates stochastically optimized structures and their energies, and show that this scheme improves upon DFT accuracy. Our methods generalize straightforwardly to other systems and other high-accuracy theories that similarly face challenges computing energy gradients, paving the way for highly accurate PES mapping, transition state determination, and thermodynamic and kinetic calculations at significantly reduced computational expense.

Iyer, Gopal R.↗

Hydrothermal synthesis of chiral carbon dots

Abstract Nanocolloids that are cumulatively referred to as nanocarbons, attracted significant attention during the last decade because of facile synthesis methods, water solubility, tunable photoluminescence, easy surface modification, and high biocompatibility. Among the latest development in this reserach area are chiral nanocarbons exemplified by chiral carbon dots (CDots). They are expected to have applications in sensing, catalysis, imaging, and nanomedicine. However, the current methods of CDots synthesis show often contradictory chemical/optical properties and structural information that required a systematic study with careful structural evaluation. Here, we investigate and optimize chiroptical activity and photoluminescence of L‐ and D‐ CDots obtained by hydrothermal carbonization of L‐ and D‐ cysteine, respectively. Nuclear magnetic resonance spectroscopy demonstrates that they are formed via gradual dehydrogenation and condensation reactions of the starting amino acid leading to particles with a wide spectrum of functional groups including aromatic cycles. We found that the chiroptical activity of CDots has an inverse correlation with the synthesis duration and temperature, whereas the photoluminescence intensity has a direct one, which is associated with degree of carbonization. Also, our studies show that the hydrothermal synthesis of cysteine in the presence of boric acid leads to the formation of CDots rather than boron nitride nanoparticles as was previously proposed in several reports. These results can be used to design chiral carbon‐based nanoparticles with optimal chemical, chiroptical, and photoluminescent properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Seismic Characterization of the Blue Mountain Geothermal Field

Subsurface characterization is crucial for geothermal energy exploration and production. Yet hydrothermal reservoirs usually reside in highly fractured and faulted zones where accurate characterization is very challenging because of low signal-to-noise ratios of land seismic data and lack of coherent reflection signals. We perform an active-source seismic characterization for the Blue Mountain geothermal field in Nevada using active seismic data to reveal the elastic medium property complexity and fault distribution at this field. We first employ an unsupervised machine learning method to attenuate groundroll and near-surface guided-wave noise and enhance coherent reflection and scattering signals from noisy seismic data. We then build a smooth initial P-wave velocity model based on an existing magnetotellurics survey result, and use 3D first-arrival traveltime tomography to refine the initial velocity model. We then derive a set of elastic wave velocities and anisotropic parameters using elastic full-waveform inversion, and obtain PP and PS images using elastic reverse-time migration. We identify major faults by analyzing the variations of seismic velocities and anisotropy parameters, and reveal mid- to small-scale faults by applying a supervised machine learning method to the seismic migration images. Our characterization reveals complex velocity heterogeneities and anisotropies, as well as faults, with a high spatial resolution. These results can provide valuable information for optimal placement of future injection and production wells to increase geothermal energy production at the Blue Mountain geothermal power plant.

58 GEOSCIENCES↗

Properties of Equivalent Source Distributions Useful for Inverse Thermal Analysis

A general approach for inverse thermal analysis of heat deposition processes is parametric modeling of such processes using effective or equivalent source distributions (ESDs). ESDs can be adjusted conveniently for calculation of temperature fields, which are constrained according to boundary and temperature constraint conditions. This study applies properties of ESDs useful for parametric modeling of heat deposition processes. Inverse thermal analysis of a friction stir weld is presented, which is prototypical of processes tending to be complex with respect to energy deposition morphology. This analysis demonstrates properties of ESDs providing relatively optimal parametric modeling of energy deposition processes.

36 MATERIALS SCIENCE↗

Co-optimization of nuclear reactor flexible power operation and maintenance scheduling

As flexible power operation of nuclear power plants becomes more attractive due to the reduction in fossil-fueled dispatchable generation on energy grids, finding optimal power production strategies that balance revenue generation with operational concerns becomes more complex. This article presents a general framework to aid operators in designing economically optimal long term dispatch strategies for nuclear power plants. The principal novelty is the linking of estimated system remaining useable life (RUL) to strategic operational decisions. It is shown that, depending on the relationship between the fixed costs from maintenance and the associated lost revenue from an outage, it can be economically optimal in the long term to delay a maintenance outage and not perform this alongside refueling. For a given relationship between power ramping and degradation, optimal strategies were found that discouraged load following in some situations while minimizing unnecessary maintenance. It is shown that heavy load following can cause maintenance and refueling outages to diverge due to their inverse relationships with respect to load following, potentially leading to a significant loss in capacity factor. As a result, this general framework can be applied to specific reactor dispatch allowing operators to adapt operational strategies as future grid conditions change.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

An analysis of the spatio-temporal resolution of the immersed boundary method with direct forcing

The immersed boundary method (IBM) with direct forcing is very popular in the simulation of rigid particulate flows. In the IBM, an interaction force is introduced at the interface between fluid and particle in order to approximate the no-slip boundary condition. The interaction force is calculated through dividing the velocity difference (or error) between fluid and particle at the interface by the time step. Here, a dynamic equation for the velocity difference is derived. Additionally, analyses on the dynamic equation provide a few new findings: (i) The interaction force is the solution of a least-squares error problem, with the direct implication that the Lagrangian marker distribution has no effect on the large scale flow structure once the distribution of Lagrangian markers become saturated along the interface (i.e., each marker remains properly correlated with all its neighbors); (ii) The Lagrangian volume-weight is a relaxation factor to control how fast the velocity error decays to the ideal value of zero; (iii) The optimal choice of the Lagrangian volume-weight is the largest value permissible by a stability condition. A comprehensive convergence analysis with regard to the spatial and temporal resolution is presented for the velocity error and also for the shear-stress and surface pressure. In three simple canonical problems, it is analytically and numerically shown that the IBM results converge to the theoretical solutions obtained with precise imposition of no-slip and no-penetration boundary conditions. It is observed that it is not necessary to match the Lagrangian marker volume-weight to that of the local Eulerian cell volume and in fact this matching leads to lower than optimal computational efficiency. However, it is found that extremely high Eulerian grid resolution and small time step have to be used to obtain high precision simulation results. Especially, the time step should be inversely proportional to the particle Reynolds number for low Reynolds number flows. For high frequency oscillation problems, the grid size needs to be reduced by a factor of the square root of the frequency, and the time step to be reduced by a factor of the frequency. The theoretical findings here can be used to alleviate the technical difficulties in simulating non-spherical particles by not requiring the Lagrangian marker distribution to match the Eulerian grids and also in the implementation of IBM on non-uniform Eulerian grids. The present work also provides simple practical guidance on the choice of temporal and spatial resolution so as to control the simulation error a priori.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Diffusion Probabilistic Modeling for Video Generation

Denoising diffusion probabilistic models are a promising new class of generative models that mark a milestone in high-quality image generation. This paper showcases their ability to sequentially generate video, surpassing prior methods in perceptual and probabilistic forecasting metrics. We propose an autoregressive, end-to-end optimized video diffusion model inspired by recent advances in neural video compression. The model successively generates future frames by correcting a deterministic next-frame prediction using a stochastic residual generated by an inverse diffusion process. We compare this approach against six baselines on four datasets involving natural and simulation-based videos. We find significant improvements in terms of perceptual quality and probabilistic frame forecasting ability for all datasets.

97 MATHEMATICS AND COMPUTING↗

Three-dimensional fracture continuum characterization aided by surface time-domain electromagnetics and hydrogeophysical joint inversion—proof-of-concept

Efficient and safe production of hydraulically fractured reservoirs benefits from the prediction of their geometrical attributes. Geophysical methods have the potential to provide data that are sensitive to fracture geometries, alleviating the typically sparse nature of in situ reservoir observations. Moreover, surface-based methods can be logistically and economically attractive since they avoid operational interference with the injection well infrastructure. Here, we investigate the potential of the surface-based time-domain electromagnetic (EM) method. EM methods can play an important role owing to their sensitivity to injection-induced fluid property changes. Two other advantageous factors are the EM signal-enhancing effect of vertical steel-cased wells and the fact that injected proppants can be enhanced to produce a stronger electrical conductivity contrast with the reservoir’s connate fluid. Nevertheless, an optimal fracture characterization will no doubt require the integration of EM and reservoir injection and production data. We hence carry out our investigations within a hydrogeophysical parameter estimation framework where EM data and injection flow rates are combined in a fully coupled way. Given the interdisciplinary nature of coupled hydrogeophysical inverse modeling, we dedicate one section to laying out key aspects in a didactic manner.

58 GEOSCIENCES↗

Engineering Antifreeze Proteins to Optimally Resist Engulfment by Ice

Antifreeze proteins (AFPs) facilitate the survival of organisms in cold climates by inhibiting the growth and/or recrystallization of ice. To function, AFPs must first bind to ice crystals; bound AFPs must then resist engulfment by using their nonbinding side (NBS) to pin the ice–water interface. Here, we seek to understand how the molecular characteristics of an NBS, such as its ice-phobicity or shape, influence its ability to resist engulfment. By characterizing the free energy barriers that impede the engulfment of model AFPs, we find that the critical supercooling ΔT*, above which an AFP is engulfed, is dictated by an optimal pinning site on the NBS. We further find that the optimal pinning site is determined by an interplay between the contact line perimeter P and a pinning efficiency η, with ΔT* ∝ ηP at the optimal pinning site. For a hemispherical AFP, which displays progressively inward tapering, we find that η increases during engulfment, whereas P decreases; conversely, an NBS with outward tapering can achieve high P, but it suffers from low η. Because the product of η and P determines ΔT*, the inverse correlation between them limits ΔT. To circumvent such limiting behavior, we propose an NBS shape with an outward bulge; by initially tapering outward, a bulged NBS permits higher P, and by subsequently tapering inward, it promotes high η as well. Importantly, we find that ΔT* is enhanced by more than a factor of 2 with an outward bulge of only 1 nm. We also find that the more ice-phobic an NBS is, the more efficiently it pins the ice–water interface, resulting in a higher ΔT. Furthermore, our findings shed light on how the NBS molecular characteristics influence ΔT* and suggest strategies for engineering the NBS to optimally resist engulfment by ice.

Antifreeze↗

Performance analysis of a finite radon transform in OFDM system under different channel models

In this paper, a class of discrete Radon transforms namely Finite Radon Transform (FRAT) was proposed as a modulation technique in the realization of Orthogonal Frequency Division Multiplexing (OFDM). The proposed FRAT operates as a data mapper in the OFDM transceiver instead of the conventional phase shift mapping and quadrature amplitude mapping that are usually used with the standard OFDM based on Fast Fourier Transform (FFT), by the way that ensure increasing the orthogonality of the system. The Fourier domain approach was found here to be the more suitable way for obtaining the forward and inverse FRAT. This structure resulted in a more suitable realization of conventional FFT- OFDM. It was shown that this application increases the orthogonality significantly in this case due to the use of Inverse Fast Fourier Transform (IFFT) twice, namely, in the data mapping and in the sub-carrier modulation also due to the use of an efficient algorithm in determining the FRAT coefficients called the optimal ordering method. The proposed approach was tested and compared with conventional OFDM, for additive white Gaussian noise (AWGN) channel, flat fading channel, and multi-path frequency selective fading channel. The obtained results showed that the proposed system has improved the bit error rate (BER) performance by reducing inter-symbol interference (ISI) and inter-carrier interference (ICI), comparing with conventional OFDM system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Atmospheric retrieval analysis of the directly imaged exoplanet HR 8799b

Directly imaged exoplanets are unexplored laboratories for the application of the spectral and temperature retrieval method, where the chemistry and composition of their atmospheres are inferred from inverse modeling of the available data. As a pilot study, we focus on the extrasolar gas giant HR 8799b, for which more than 50 data points are available. We upgrade our non-linear optimal estimation retrieval method to include a phenomenological model of clouds that requires the cloud optical depth and monodisperse particle size to be specified. Previous studies have focused on forward models with assumed values of the exoplanetary properties; there is no consensus on the best-fit values of the radius, mass, surface gravity, and effective temperature of HR 8799b. We show that cloud-free models produce reasonable fits to the data if the atmosphere is of super-solar metallicity and non-solar elemental abundances. Intermediate cloudy models with moderate values of the cloud optical depth and micron-sized particles provide an equally reasonable fit to the data and require a lower mean molecular weight. We report our best-fit values for the radius, mass, surface gravity, and effective temperature of HR 8799b. The mean molecular weight is about 3.8, while the carbon-to-oxygen ratio is about unity due to the prevalence of carbon monoxide. Our study emphasizes the need for robust claims about the nature of an exoplanetary atmosphere to be based on analyses involving both photometry and spectroscopy and inferred from beyond a few photometric data points, such as are typically reported for hot Jupiters.

79 ASTRONOMY AND ASTROPHYSICS↗

Tau-REx I: A next generation retrieval code for exoplanetary atmospheres

Spectroscopy of exoplanetary atmospheres has become a well established method for the characterization of extrasolar planets. We here present a novel inverse retrieval code for exoplanetary atmospheres. T-REx(Tau Retrieval for Exoplanets) is a line-by-line radiative transfer fully Bayesian retrieval framework. T-REx includes the following features:(1) the optimized use of molecular line lists from the ExoMol project; (2) an unbiased atmospheric composition prior selection, through custom built pattern recognition software; (3) the use of two independent algorithms to fully sample the Bayesian likelihood space: nested sampling as well as a more classical Markov Chain Monte Carlo approach; (4) iterative Bayesian parameter and model selection using the full Bayesian Evidence as well as the Savage–Dickey Ratio for nested models; and (5) the ability to fully map very large parameter spaces through optimal code parallelization and scalability to cluster computing. In this publication we outline the T-REx framework and demonstrate, using a theoretical hot-Jupiter transmission spectrum, the parameter retrieval and model selection. We investigate the impact of signal-to-noise ratio and spectral resolution on the retrievability of individual model parameters, both in terms of error bars on the temperature and molecular mixing ratios as well as its effect on the model’s global Bayesian evidence.

79 ASTRONOMY AND ASTROPHYSICS↗

Thermal diffusivity above the Mott-Ioffe-Regel limit

We present high-resolution thermal diffusivity measurements on several near optimally doped electron- and hole-doped cuprate systems in a temperature range that passes through the Mott-Ioffe-Regel limit, above which the quasiparticle picture fails. Our primary observations are that the inverse thermal diffusivity is linear in temperature and can be fitted to $D^{-1}_Q = aT + b$. The slope a is interpreted through the Planckian relaxation time $τ ≈ ℏ/k_BT$ and a thermal diffusion velocity v B, which is close, but larger than the sound velocity. The intercept b represents a crossover diffusion constant that separates coherent from incoherent quasiparticles. These observations suggest that both phonons and electrons participate in the thermal transport, while reaching the Planckian limit for relaxation time.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗