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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Galerkin approximation for inverse problems for nonautonomous nonlinear distributed systems

An abstract framework and convergence theory is developed for Galerkin approximation for inverse problems involving the identification of nonautonomous nonlinear distributed parameter systems. A set of relatively easily verified conditions is provided which are sufficient to guarantee the existence of optimal solutions and their approximation by a sequence of solutions to a sequence of approximating finite dimensional identification problems. The approach is based on the theory of monotone operators in Banach spaces and is applicable to a reasonably broad class of nonlinear distributed systems. Operator theoretic and variational techniques are used to establish a fundamental convergence result. An example involving evolution systems with dynamics described by nonstationary quasilinear elliptic operators along with some applications are presented and discussed.

Banks, H. T.↗
NASA NTRS1988

Inversion of helioseismology data

The paper presents an algorithm for performing an inversion of helioseismology data. The method is characterized by high speed, stability, and versatility. The existence of a method for producing an analog of the Backus-Gilbert optimal kernals is demonstrated. The method performs well, even when a moderate amount of observational noise is included.

Jeffrey, W.↗
NASA NTRS1987

Kinematically redundant robot manipulators

Research on control, design and programming of kinematically redundant robot manipulators (KRRM) is discussed. These are devices in which there are more joint space degrees of freedom than are required to achieve every position and orientation of the end-effector necessary for a given task in a given workspace. The technological developments described here deal with: kinematic programming techniques for automatically generating joint-space trajectories to execute prescribed tasks; control of redundant manipulators to optimize dynamic criteria (e.g., applications of forces and moments at the end-effector that optimally distribute the loading of actuators); and design of KRRMs to optimize functionality in congested work environments or to achieve other goals unattainable with non-redundant manipulators. Kinematic programming techniques are discussed, which show that some pseudo-inverse techniques that have been proposed for redundant manipulator control fail to achieve the goals of avoiding kinematic singularities and also generating closed joint-space paths corresponding to close paths of the end effector in the workspace. The extended Jacobian is proposed as an alternative to pseudo-inverse techniques.

Baillieul, J.↗
NASA NTRS2009

Automation of a Versatile Crane (the LSMS) for Lunar Outpost Construction, Maintenance and Inspection

Devices for lifting, translating and precisely placing payloads are critical for efficient Earth-based construction operations. Both recent and past studies have demonstrated that devices with similar functionality will be needed to support lunar outpost operations. Although several designs have been developed for Earth based applications, these devices lack unique design characteristics necessary for transport to and use on the harsh lunar surface. These design characteristics include: a) lightweight components, b) compact packaging for launch, c) automated deployment, d) simple in-field reconfiguration and repair, and e) support for tele-operated or automated operations. Also, because the cost to transport mass to the lunar surface is very high, the number of devices that can be dedicated to surface operations will be limited. Thus, in contrast to Earth-based construction, where many single-purpose devices dominate a construction site, a lunar outpost will require a limited number of versatile devices that provide operational benefit from initial construction through sustained operations. The first generation test-bed of a new high performance device, the Lunar Surface Manipulation System (LSMS) has been designed, built and field tested. The LSMS has many unique features resulting in a mass efficient solution to payload handling on the lunar surface. Typically, the LSMS device mass is estimated at approximately 3% of the mass of the heaviest payload lifted at the tip, or 1.8 % of the mass of the heaviest mass lifted at the elbow or mid-span of the boom for a high performance variant incorporating advanced structural components. Initial operational capabilities of the LSMS were successfully demonstrated during field tests at Moses Lake, Washington using a tele-operated approach. Joint angle sensors have been developed for the LSMS to improve operator situational awareness. These same sensors provide the necessary information to support fully automated operations, greatly expanding the operational versatility of the LSMS. This paper develops the equations describing the forward and inverse relation between LSMS joint angles and Cartesian coordinates of the LSMS tip. These equations allow a variety of schemes to be used to maneuver the LSMS to optimize the maneuver. One such scheme will be described in detail that eliminates undesirable swinging of the payload at the conclusion of a maneuver, even when the payload is suspended from a passive rigid link. The swinging is undesirable when performing precision maneuvers, such as aligning an object for mating or positioning a camera. Use of the equations described here enables automated control of the LSMS greatly improving its operational versatility.

Doggett, William R.↗
DOE OSTI2017

Selective Shielding Bone Marrow: An Approach to Create Personal Protective Equipment for Gamma Radiation

Many casualties of the Hiroshima and Nagasaki atomic bombs, and many of the firefighters who first responded to the Chernobyl nuclear power plant accident, became ill with acute radiation syndrome (ARS). Most of the people who do not recover from ARS will die within a few weeks to a few months after exposure, with the primary cause of death being the destruction of the body's bone marrow (BM). This is exemplified by the Chernobyl firefighters, the majority of whom received whole body doses between 0.8 and 10 Gy, and the United Nations Scientific Committee on the Effects of Atomic Radiation (UNSCEAR) final report on the disaster concluded that the underlying cause of death of all 28 firefighters who succumbed to ARS was from bone marrow failure. BM is comprised of hematopoietic stem cells (HSCs) which are responsible for the constant renewal of blood cells. Due to their high rate of proliferation, HSCs are especially vulnerable to ionizing radiation, but are endowed with remarkable regenerative potential. Owing to their central role in blood production, lethal irradiation of HSCs leads to death from severe anemia, infection and internal bleeding. This relationship between high doses of radiation and HSC apoptosis has led to the use of HSC transplantation as a life-saving intervention in cases of acute exposure. This is demonstrable in medical practice, where thousands of individuals have undergone supra-lethal Total Body Irradiation (TBI) for purposes of cancer therapy, and were rescued by well-matched BM transplantation. Indeed, life-threatening damage may be reversed by BM transplantation in individuals receiving doses of radiation of up to 10 Gy. At doses > 10 Gy, damage to gastrointestinal (GI) tissue may become a limiting factor in survival. Remarkably, in the common procedure of BM transplantation, the number of HSCs extracted from a single active BM site containing less than 5% of the donor's BM tissue is sufficient to support the complete reconstitution of the HSC compartment in a lethally irradiated recipient. This capacity of BM to expand and replenish is due to the high regenerative potential of the HSCs it harbors. In contrast to mature blood cells, which are dispersed throughout the body, HSCs are confined to the bones, allowing for effective targeted shielding. This has been confirmed in several animal models, where sparing the HSC-rich area of a subject receiving otherwise full-body irradiation is sufficient to support hematopoietic functions and allow survival. Approximately 48% of the human body's active BM is contained within the lumbar vertebrae and the pelvic region. This high concentration of BM makes the pelvic region an ideal area to shield for preventing the serious hematopoietic effects of radiation exposure. Stem-cell rich tissues like active BM are also more sensitive to the stochastic effects (i.e. cancer) of radiation exposure. The high amount of BM in the pelvic region, combined with the presence of the sensitive gastrointestinal system and female gonads in this same area of the body allows for the possibility of significantly reducing stochastic effects of radiation by limiting radiation exposure to this vital area. Thus, shielding this region holds great promise for both acute and protracted exposures. Thus, a specialized radiation shielding device for the protection of active BM concentrations is presented here. This belt-like selective shielding device focuses on the protection of BM that is present in the pelvis and hip bones, protecting the medullary volume from which BM is commonly extracted for transplantation (i.e. the iliac bones) while allowing relatively unhindered movement of the wearer. In order to optimize the use of shielding materials towards the protection of active BM, the shielding uniquely brings into account the natural shielding properties of human tissue, with thickness being inversely related to the thickness and radiation attenuation of the underlying tissue at each point surrounding the area being protected. The findings resulting from both experimental and simulation testing of this selective shielding device are presented here.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2022

Transient Storage Model Parameter Optimization Using the Simulated Annealing Method

Hyporheic exchange in streams is critical to ecosystem functions such as nutrient cycling along river corridors, especially for slowly moving or small stream systems. The transient storage model (TSM) has been widely used for modeling of hyporheic exchange. TSM calibration, for hyporheic exchange, is typically used to estimate four parameters, including the mass exchange rate coefficient, the dispersion coefficient, stream cross-sectional area, and hyporheic zone cross-sectional area. Prior studies have raised concerns regarding the non-uniqueness of the inverse problem for the TSM, that is, the occurrence of different parameter vectors resulting in TSM solution that reproduces the observed in-stream tracer break through curve (BTC) with the same error. This leads to practical non-identifiability in determining the unknown parameter vector values even when global-optimal values exist, and the parameter optimization becomes practically non-unique. To address this problem, we applied the simulated annealing method to calibrate the TSM to BTCs, because it is less susceptible to local minima-induced non-identifiability. A hypothetical (or synthetic) tracer test data set with known parameters was developed to demonstrate the capability of the simulated annealing method to find the global minimum parameter vector, and it identified the “hypothetically-true” global minimum parameter vector even with input data that were modified with up to 10% noise without increasing the number of iterations required for convergence. The simulated annealing TSM was then calibrated using two in-stream tracer tests conducted in East Fork Poplar Creek, Tennessee. In conclusion, simulated annealing was determined to be appropriate for quantifying the TSM parameter vector because of its search capability for the global minimum parameter vector.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1994

Preconditioned domain decomposition scheme for three-dimensional aerodynamic sensitivity analysis

A discrete sensitivity analysis algorithm had previously been developed and applied to two-dimensional aerodynamic optimization problems, where the computational domains were discretized by using single grids. The sparse, unsymmetric systems of linear equations resulting from this algorithm were solved by a direct matrix inversion matrix. However, for large two-dimensional problems and, practically, all three-dimensional problems, direct inversion methods become inapplicable, primarily due to the prohibitive computer storage needed. In an attempt to alleviate such hindrances, the sensitivity analysis with domain decomposition (SADD) scheme was developed. This scheme divides the computational domain into smaller and nonoverlapping subdomains (multiblock grids) that are solved separately. Then, the final solution is constructed from the subdomain solutions. As the number of grid points in the interface boundaries of the subdomains becomes large, the computer memory required to store the effective coefficient matrix of these interface points starts to increase. Presented in this Technical Note is a preconditioned iterative procedure to overcome this particular problem.

Eleshaky, Mohamed E.↗
NASA NTRS2013

Real-Time Minimization of Tracking Error for Aircraft Systems

This technology presents a novel, stable, discrete-time adaptive law for flight control in a Direct adaptive control (DAC) framework. Where errors are not present, the original control design has been tuned for optimal performance. Adaptive control works towards achieving nominal performance whenever the design has modeling uncertainties/errors or when the vehicle suffers substantial flight configuration change. The baseline controller uses dynamic inversion with proportional-integral augmentation. On-line adaptation of this control law is achieved by providing a parameterized augmentation signal to a dynamic inversion block. The parameters of this augmentation signal are updated to achieve the nominal desired error dynamics. If the system senses that at least one aircraft component is experiencing an excursion and the return of this component value toward its reference value is not proceeding according to the expected controller characteristics, then the neural network (NN) modeling of aircraft operation may be changed.

Garud, Sumedha↗
DOE OSTI2015

Optimization of leaf margins for lung stereotactic body radiotherapy using a flattening filter-free beam

Purpose: The authors sought to determine the optimal collimator leaf margins which minimize normal tissue dose while achieving high conformity and to evaluate differences between the use of a flattening filter-free (FFF) beam and a flattening-filtered (FF) beam. Methods: Sixteen lung cancer patients scheduled for stereotactic body radiotherapy underwent treatment planning for a 7 MV FFF and a 6 MV FF beams to the planning target volume (PTV) with a range of leaf margins (−3 to 3 mm). Forty grays per four fractions were prescribed as a PTV D95. For PTV, the heterogeneity index (HI), conformity index, modified gradient index (GI), defined as the 50% isodose volume divided by target volume, maximum dose (Dmax), and mean dose (Dmean) were calculated. Mean lung dose (MLD), V20 Gy, and V5 Gy for the lung (defined as the volumes of lung receiving at least 20 and 5 Gy), mean heart dose, and Dmax to the spinal cord were measured as doses to organs at risk (OARs). Paired t-tests were used for statistical analysis. Results: HI was inversely related to changes in leaf margin. Conformity index and modified GI initially decreased as leaf margin width increased. After reaching a minimum, the two values then increased as leaf margin increased (“V” shape). The optimal leaf margins for conformity index and modified GI were −1.1 ± 0.3 mm (mean ± 1 SD) and −0.2 ± 0.9 mm, respectively, for 7 MV FFF compared to −1.0 ± 0.4 and −0.3 ± 0.9 mm, respectively, for 6 MV FF. Dmax and Dmean for 7 MV FFF were higher than those for 6 MV FF by 3.6% and 1.7%, respectively. There was a positive correlation between the ratios of HI, Dmax, and Dmean for 7 MV FFF to those for 6 MV FF and PTV size (R = 0.767, 0.809, and 0.643, respectively). The differences in MLD, V20 Gy, and V5 Gy for lung between FFF and FF beams were negligible. The optimal leaf margins for MLD, V20 Gy, and V5 Gy for lung were −0.9 ± 0.6, −1.1 ± 0.8, and −2.1 ± 1.2 mm, respectively, for 7 MV FFF compared to −0.9 ± 0.6, −1.1 ± 0.8, and −2.2 ± 1.3 mm, respectively, for 6 MV FF. With the heart inside the radiation field, the mean heart dose showed a V-shaped relationship with leaf margins. The optimal leaf margins were −1.0 ± 0.6 mm for both beams. Dmax to the spinal cord showed no clear trend for changes in leaf margin. Conclusions: The differences in doses to OARs between FFF and FF beams were negligible. Conformity index, modified GI, MLD, lung V20 Gy, lung V5 Gy, and mean heart dose showed a V-shaped relationship with leaf margins. There were no significant differences in optimal leaf margins to minimize these parameters between both FFF and FF beams. The authors’ results suggest that a leaf margin of −1 mm achieves high conformity and minimizes doses to OARs for both FFF and FF beams.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Double sweep preconditioner for optimized Schwarz methods applied to the Helmholtz problem

This paper presents a preconditioner for non-overlapping Schwarz methods applied to the Helmholtz problem. Starting from a simple analytic example, we show how such a preconditioner can be designed by approximating the inverse of the iteration operator for a layered partitioning of the domain. The preconditioner works by propagating information globally by concurrently sweeping in both directions over the subdomains, and can be interpreted as a coarse grid for the domain decomposition method. The resulting algorithm is shown to converge very fast, independently of the number of subdomains and frequency. The preconditioner has the advantage that, like the original Schwarz algorithm, it can be implemented as a matrix-free routine, with no additional preprocessing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1976

Analytical and experimental design and analysis of an optimal processor for image registration

The author has identified the following significant results. A quantitative measure of the registration processor accuracy in terms of the variance of the registration error was derived. With the appropriate assumptions, the variance was shown to be inversely proportional to the square of the effective bandwidth times the signal to noise ratio. The final expressions were presented to emphasize both the form and simplicity of their representation. In the situation where relative spatial distortions exist between images to be registered, expressions were derived for estimating the loss in output signal to noise ratio due to these spatial distortions. These results are in terms of a reduction factor.

Mcgillem, C. D.↗
DOE OSTI2025

Inverse design of photonic surfaces via multi fidelity ensemble framework and femtosecond laser processing

We demonstrate a multi-fidelity (MF) machine learning ensemble framework for the inverse design of photonic surfaces, trained on a dataset of 11,759 samples that we fabricate using high throughput femtosecond laser processing. The MF ensemble combines an initial low fidelity model for generating design solutions, with a high fidelity model that refines these solutions through local optimization. The combined MF ensemble can generate multiple disparate sets of laser-processing parameters that can each produce the same target input spectral emissivity with high accuracy (root mean squared errors < 2%). SHapley Additive exPlanations analysis shows transparent model interpretability of the complex relationship between laser parameters and spectral emissivity. Finally, the MF ensemble is experimentally validated by fabricating and evaluating photonic surface designs that it generates for improved efficiency energy harvesting devices. Our approach provides a powerful tool for advancing the inverse design of photonic surfaces in energy harvesting applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Constraints on methane emissions in North America from future geostationary remote-sensing measurements

The success of future geostationary (GEO) satellite observation missions depends on our ability to design instruments that address their key scientific objectives. Here, an Observation System Simulation Experiment (OSSE) is performed to quantify the constraints on methane (CH 4 ) emissions in North America obtained from shortwave infrared (SWIR), thermal infrared (TIR), and multi-spectral (SWIR+TIR) measurements in geostationary orbit and from future SWIR low-Earth orbit (LEO) measurements. Furthermore, we used an efficient stochastic algorithm to compute the information content of the inverted emissions at high spatial resolution (0.5° × 0.7°) in a variational framework using the GEOS-Chem chemistry-transport model and its adjoint. Our results show that at sub-weekly timescales, SWIR measurements in GEO orbit can constrain about twice as many independent flux patterns than in LEO orbit, with a degree of freedom for signal (DOF) for the inversion of 266 and 115, respectively. Comparisons between TIR GEO and SWIR LEO configurations reveal that poor boundary layer sensitivities for the TIR measurements cannot be compensated for by the high spatiotemporal sampling of a GEO orbit. The benefit of a multi-spectral instrument compared to current SWIR products in a GEO context is shown for sub-weekly timescale constraints, with an increase in the DOF of about 50 % for a 3-day inversion. Our results further suggest that both the SWIR and multi-spectral measurements on GEO orbits could almost fully resolve CH 4 fluxes at a spatial resolution of at least 100 km × 100 km over source hotspots (emissions > 4 × 10 5 kg day -1 ). The sensitivity of the optimized emission scaling factors to typical errors in boundary and initial conditions can reach 30 and 50 % for the SWIR GEO or SWIR LEO configurations, respectively, while it is smaller than 5 % in the case of a multi-spectral GEO system. Our results demonstrate that multi-spectral measurements from a geostationary satellite platform would address the need for higher spatiotemporal constraints on CH 4 emissions while greatly mitigating the impact of inherent uncertainties in source inversion methods on the inferred fluxes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

RLMolLM: Reinforcement Learning-Enhanced Language Model Framework for Inverse Molecular Design

Inverse molecular design faces significant challenges due to vast chemical space and complex property requirements. While language models show promise for molecular generation, they struggle with validity, multi-property optimization, and structural constraints. This work presents RLMolLM, a reinforcement learning framework combining Proximal Policy Optimization (PPO) with genetic algorithms to address these limitations. Our approach optimizes multiple user-specified properties including quantitative estimates of drug-likeness (QED), synthetic accessibility (SA), and ADMET (absorption, distribution, metabolism, excretion, and toxicity) endpoints without requiring complete model retraining, while maintaining capability for scaffold-constrained generation where specific substructures must be preserved. We outperform state-of-the-art methods for molecular optimization, achieving best QED scores across GDB13, Moses, and Zinc datasets with up to 31% improvement over previous methods while maintaining excellent validity, uniqueness, and novelty metrics. For simultaneous multi-property optimization, our framework achieves substantial improvements in ADMET properties including 4.5-fold reduction in hERG toxicity and enhanced Caco-2 permeability compared to Moses dataset. Under structural constraints, the framework significantly improves molecular validity while preserving scaffolds and effectively optimizing properties. In conclusion, this versatile solution advances pharmaceutical and materials molecular design through effective integration of reinforcement learning and genetic algorithms with multi-property optimization and scaffold preservation.

Genetic algorithms↗
DOE OSTI2025

Polynomial Scaling Localized Active Space Unitary Selective Coupled Cluster Singles and Doubles

We present a polynomial-scaling algorithm for the localized active space unitary selective coupled cluster singles and doubles (LAS-USCCSD) method. In this approach, cluster excitations are selected based on a threshold ϵ determined by the absolute gradients of the LAS-UCCSD energy with respect to cluster amplitudes. Using the generalized Wick’s theorem for multireference wave functions, we derive the gradient expression as a polynomial function of one-, two-, and three-body reduced density matrices and 1- and 2-electron integrals, valid for any multireference wave function. The resulting gradient implementation exhibits a memory scaling of 𝒪(N 6 ), with N spin orbitals in the combined active space of all fragments. The variational quantum eigensolver is used to optimize the selected cluster excitations on a quantum simulator. Furthermore, by plotting the energy error, defined as the difference between the LAS-USCCSD and corresponding CASCI energies, against the inverse cluster amplitude selection threshold (ϵ –1 ) for polyene chains containing 2 to 5 π-bond units, we establish a relationship between the energy error and the threshold. To further validate the accuracy of LAS-USCCSD, we computed the cis–trans isomerization energy of stilbene (a 20-qubit system) and the magnetic coupling constant of the tris-hydroxo-bridged chromium dimer [Cr 2 (OH) 3 (NH 3 ) 6 ] 3+ (evaluated as both 12- and 20-qubit systems) using the Qiskit-Qulacs simulator. Assessing such examples is important to determine the practical feasibility of quantum simulations for chemically realistic systems. Toward this goal, with the LAS-USCCSD algorithm we estimated the quantum resources required for simulating an active space of (30e,22o) in [Cr 2 (OH) 3 (NH 3 ) 6 ] 3+ , a size that remains beyond the reach of current quantum simulators for accurate treatment.

Algorithms↗
DOE OSTI2015

A general framework of noise suppression in material decomposition for dual-energy CT

Purpose: As a general problem of dual-energy CT (DECT), noise amplification in material decomposition severely reduces the signal-to-noise ratio on the decomposed images compared to that on the original CT images. In this work, the authors propose a general framework of noise suppression in material decomposition for DECT. The method is based on an iterative algorithm recently developed in their group for image-domain decomposition of DECT, with an extension to include nonlinear decomposition models. The generalized framework of iterative DECT decomposition enables beam-hardening correction with simultaneous noise suppression, which improves the clinical benefits of DECT. Methods: The authors propose to suppress noise on the decomposed images of DECT using convex optimization, which is formulated in the form of least-squares estimation with smoothness regularization. Based on the design principles of a best linear unbiased estimator, the authors include the inverse of the estimated variance–covariance matrix of the decomposed images as the penalty weight in the least-squares term. Analytical formulas are derived to compute the variance–covariance matrix for decomposed images with general-form numerical or analytical decomposition. As a demonstration, the authors implement the proposed algorithm on phantom data using an empirical polynomial function of decomposition measured on a calibration scan. The polynomial coefficients are determined from the projection data acquired on a wedge phantom, and the signal decomposition is performed in the projection domain. Results: On the Catphan{sup ®}600 phantom, the proposed noise suppression method reduces the average noise standard deviation of basis material images by one to two orders of magnitude, with a superior performance on spatial resolution as shown in comparisons of line-pair images and modulation transfer function measurements. On the synthesized monoenergetic CT images, the noise standard deviation is reduced by a factor of 2–3. By using nonlinear decomposition on projections, the authors’ method effectively suppresses the streaking artifacts of beam hardening and obtains more uniform images than their previous approach based on a linear model. Similar performance of noise suppression is observed in the results of an anthropomorphic head phantom and a pediatric chest phantom generated by the proposed method. With beam-hardening correction enabled by their approach, the image spatial nonuniformity on the head phantom is reduced from around 10% on the original CT images to 4.9% on the synthesized monoenergetic CT image. On the pediatric chest phantom, their method suppresses image noise standard deviation by a factor of around 7.5, and compared with linear decomposition, it reduces the estimation error of electron densities from 33.3% to 8.6%. Conclusions: The authors propose a general framework of noise suppression in material decomposition for DECT. Phantom studies have shown the proposed method improves the image uniformity and the accuracy of electron density measurements by effective beam-hardening correction and reduces noise level without noticeable resolution loss.

60 APPLIED LIFE SCIENCES↗