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 1,135 records · Page 63

Optimal inputs for system identification

Identification criteria are presented for linear dynamic systems with and without process noise. With process noise, the state equations are replaced by the Kalman filter equations. If the identification performance index is expanded in a Taylor's series with respect to the parameters to be identified, then maximizing the weighting factor of the quadratic term with respect to the inputs will insure that an identification algorithm will converge more rapidly and to a more accurate result than with nonoptimal inputs. The expectation of this weighting factor is the Fisher information matrix, and its inverse is a lower bound for the covariance of the parameters. Direct and indirect methods of calculating the information matrix are presented for systems with and without process noise.

Reid, D. B.↗

Strong Upper-Plate Heterogeneity at the Hikurangi Subduction Margin (North Island, New Zealand) Imaged by Adjoint Tomography

We use earthquake-based adjoint tomography to invert for three-dimensional structure of the North Island, New Zealand, and the adjacent Hikurangi subduction zone. The study area, having a shallow depth to the plate interface below the North Island, offers a rare opportunity for imaging material properties at an active subduction zone using land-based measurements. Starting from an initial model derived using ray tomography, we perform iterative model updates using spectral element and adjoint simulations to fit waveforms with periods ranging from 4–30 s. We perform 28 model updates using an L-BFGS optimization algorithm, improving data fit and introducing P- and S-wave velocity changes of up to ±30%. Resolution analysis using point spread functions show that our measurements are most sensitive to heterogeneities in the upper 30 km. The most striking velocity changes coincide with areas related to the active Hikurangi subduction zone. Lateral velocity structures in the upper 5 km correlate well with New Zealand geology. Here, the inversion reveals increased along-strike heterogeneity on the margin. In Cook Strait we observe a low-velocity zone interpreted as deep sedimentary basins. In the central North Island, low-velocity anomalies are linked to surface geology, and we relate velocity structures at depth to crustal magmatic activity below the Taupō Volcanic Zone.

58 GEOSCIENCES↗

Preparation and characterization of nanocomposite polymer electrolytes poly(vinylidone fluoride)/nanoclay

Polymer electrolytes are defined as semi solid electrolytes used as separator in lithium ion battery. Separator used as medium for transfer ions and to prevent electrical short circuits in battery cells. To obtain the optimal battery performance, separator with high porosity and electrolyte uptake is required. This can reduce the resistance in the transfer of ions between cathode and anode. The main objective of this work is to investigate the impact of different solvent (Dimethyl acetamide (DMAc), N-methyl-2-pyrrolidone (NMP) and dimethyl formamide (DMF)), pore forming agent poly(vinylpyrolidone) (PVP) and nanoclay as filler in addition of membrane using phase inversion method on the morphology, porosity, electrolyte uptake and degree of crystallinity. The membrane was prepared by the phase inversion method by adding PVP and Nanoclay using different solvents. The phase inversion method was prepared by dissolving Nanoclay and PVP in solvent for 1-2 hours, and then add the PVDF with stirring for 4 hours at 60°C. The membranes were characterized by porosity test, electrolyte uptake test, scanning electron microscope (SEM), and X-ray diffraction (XRD). The results showed that DMAc as solvent gives the highest value of porosity and electrolyte uptake. The addition of nanoclay and PVP enlarge the size of the pores and reduce the degree of crystallinity. So, the usage of DMAc as solvent is better than NMP or DMF.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Adjustable collimation for dark-field proton radiography contrast enhancement

A dark-field imaging condition established through pre- and post-object beam collimation has been tested for a variety of different settings at the proton radiography facility (pRad) at Los Alamos National Laboratory. This technique facilitates imaging of thin samples and materials with small areal density changes with 800 MeV protons that would otherwise appear nearly transparent for high energy protons. In the configuration that was examined, the pre-object (upstream) collimator is utilized to remove protons with large scattering angles prior to the target object. The post-object collimation entails a combination of a conventional collimator, which eliminates protons with large scattering angles after target interaction, and an inverse collimator. The inverse collimator removes unscattered protons traveling on the beam axis and those with smaller scattering angles up to a desired angle. This establishes a dark-field condition by permitting solely protons that exhibit scattering angles from target interaction within a designated range to contribute to the image. This study explores the potential of combining different sized collimators for pre- and post-object collimation. The objective is to optimize contrast for a range of samples with minimal areal density variations. A series of dark-field conditions are explored employing a remote-controlled collimator wheel and the contrast-to-noise ratio, spatial resolution, and proton transmission are compared for those combinations. It is found that the contrast-to-noise ratio can be substantially enhanced for thin samples and low areal density objects, thereby facilitating capturing high-contrast images of objects that were previously invisible to high energy protons.

43 PARTICLE ACCELERATORS↗

Thin ice clouds in the Arctic: cloud optical depth and particle size retrieved from ground-based thermal infrared radiometry

Multiband downwelling thermal measurements of zenith sky radiance, along with cloud boundary heights, were used in a retrieval algorithm to estimate cloud optical depth and effective particle diameter of thin ice clouds in the Canadian High Arctic. Ground-based thermal infrared (IR) radiances for 150 semitransparent ice clouds cases were acquired at the Polar Environment Atmospheric Research Laboratory (PEARL) in Eureka, Nunavut, Canada (80° N, 86° W). We analyzed and quantified the sensitivity of downwelling thermal radiance to several cloud parameters including optical depth, effective particle diameter and shape, water vapor content, cloud geometric thickness and cloud base altitude. A lookup table retrieval method was used to successfully extract, through an optimal estimation method, cloud optical depth up to a maximum value of 2.6 and to separate thin ice clouds into two classes: (1) TIC1 clouds characterized by small crystals (effective particle diameter ≤ 30 µm), and (2) TIC2 clouds characterized by large ice crystals (effective particle diameter > 30 µm). The retrieval technique was validated using data from the Arctic High Spectral Resolution Lidar (AHSRL) and Millimeter Wave Cloud Radar (MMCR). Inversions were performed over three polar winters and results showed a significant correlation ( R 2 = 0.95) for cloud optical depth retrievals and an overall accuracy of 83 % for the classification of TIC1 and TIC2 clouds. A partial validation relative to an algorithm based on high spectral resolution downwelling IR radiance measurements between 8 and 21µm was also performed. It confirms the robustness of the optical depth retrieval and the fact that the broadband thermal radiometer retrieval was sensitive to small particle (TIC1) sizes.

54 ENVIRONMENTAL SCIENCES↗

Method and apparatus for configuration control of redundant robots

A method and apparatus to control a robot or manipulator configuration over the entire motion based on augmentation of the manipulator forward kinematics is disclosed. A set of kinematic functions is defined in Cartesian or joint space to reflect the desirable configuration that will be achieved in addition to the specified end-effector motion. The user-defined kinematic functions and the end-effector Cartesian coordinates are combined to form a set of task-related configuration variables as generalized coordinates for the manipulator. A task-based adaptive scheme is then utilized to directly control the configuration variables so as to achieve tracking of some desired reference trajectories throughout the robot motion. This accomplishes the basic task of desired end-effector motion, while utilizing the redundancy to achieve any additional task through the desired time variation of the kinematic functions. The present invention can also be used for optimization of any kinematic objective function, or for satisfaction of a set of kinematic inequality constraints, as in an obstacle avoidance problem. In contrast to pseudoinverse-based methods, the configuration control scheme ensures cyclic motion of the manipulator, which is an essential requirement for repetitive operations. The control law is simple and computationally very fast, and does not require either the complex manipulator dynamic model or the complicated inverse kinematic transformation. The configuration control scheme can alternatively be implemented in joint space.

Seraji, Homayoun↗

Highly-scalable, physics-informed GANs for learning solutions of stochastic PDEs

Uncertainty quantification for forward and inverse problems is a central challenge across physical and biomedical disciplines. We address this challenge for the problem of modeling subsurface flow at the Hanford Site by combining stochastic computational models with observational data using physics-informed GAN models. The geographic extent, spatial heterogeneity, and multiple correlation length scales of the Hanford Site require training a computationally intensive GAN model to thousands of dimensions. We develop a highly optimized implementation that scales to {\herogpucount} NVIDIA Volta GPUs. We develop a hierarchical scheme based on a multi-player game-theoretic approach for exploiting domain parallelism, map discriminators and generators to multiple GPUs, and employ efficient communication schemes to ensure training stability and convergence. Our implementation scales to {\heronodecount} nodes on the Summit supercomputer with a \weakscalingeff\% scaling efficiency, achieving peak and sustained half-precision rates of {\heropeakperf} PF/s and {\herosustainedperf} PF/s.

physics informed ML, machine learning, generative ↗

Projected Multi-Agent Consensus Equilibrium (PMACE) With Application to Ptychography

Multi-Agent Consensus Equilibrium (MACE) formulates an inverse imaging problem as a balance among multiple update agents such as data-fitting terms and denoisers. However, each such agent operates on a separate copy of the full image, leading to redundant memory use and slow convergence when each agent affects only a small subset of the full image. In this article, we extend MACE to Projected Multi-Agent Consensus Equilibrium (PMACE), in which each agent updates only a projected component of the full image, thus greatly reducing memory use for some applications. We describe PMACE in terms of an equilibrium problem and an equivalent fixed-point problem and show that in most cases the PMACE equilibrium is not the solution of an optimization problem. To demonstrate the value of PMACE, we apply it to the problem of ptychography, in which a sample is reconstructed from the diffraction patterns resulting from coherent X-ray illumination at multiple overlapping spots. In our PMACE formulation, each spot corresponds to a separate data-fitting agent, with the final solution found as an equilibrium among all the agents. In conclusion, our results demonstrate that the PMACE reconstruction algorithm generates more accurate reconstructions at a lower computational cost than existing ptychography algorithms when the spots are sparsely sampled.

97 MATHEMATICS AND COMPUTING↗

TPSAS-NF1676L-27509-DND

The information about composition of aerosols has a great importance for monitoring and understanding of climate and environment dynamics. Such information can be obtained using situ measurements or chemical transport models. However, in situ sampling has limited spatial and temporal coverage, while estimations have large uncertainties. The present work enables the monitoring of aerosol chemical species from space-borne observations, providing observationally-based results with spatial and temporal coverage. Following the ideas of Schuster et al.(2005, 2009, 2016), we retrieve chemical composition directly from satellite observations without intermediate retrieval of the refractive index (in contrast with Schuster's approach). This approach is expected to reduce the influence of modeling uncertainties, and to provide additional constraints in situations where remote sensing observations do not have enough spectral sensitivity to refractive index. One of principal difficulties is the identification of an adequate model for linking refractive index to chemical composition. Therefore, the initial effort of this work has focused on identifying an optimal physico-chemistry to refractive index conversion model. With that purpose, we first tested the retrieval approach using a simplified volume-weighting model and updated by the Maxwell Garnett mixing model. This concept was incorporated into the GRASP algorithm designed to retrieve an extended set of atmospheric parameters from remote sensing observations. Then a series of sensitivity tests using synthetic data of POLDER/PARASOL were conducted, and followed by inversion of real PARASOL observations over Africa. The sensitivity tests showed that these two models allow the retrieval to distinguish amongst the assumed chemical species. Results obtained from real PARASOL data demonstrated good agreement with the optical characteristics provided by AERONET. In addition, the obtained patterns of chemical component distribution agreed with known physical expectations. We also compared the chemical composition retrievals to GOCART model simulations and obtained good consistencies.

Lei Li↗

SU-F-T-192: Study of Robustness Analysis Method of Multiple Field Optimized IMPT Plans for Head & Neck Patients

Purpose: Proton therapy is more sensitive to uncertainties than photon treatments due to protons’ finite range depending on the tissue density. Worst case scenario (WCS) method originally proposed by Lomax has been adopted in our institute for robustness analysis of IMPT plans. This work demonstrates that WCS method is sufficient enough to take into account of the uncertainties which could be encountered during daily clinical treatment. Methods: A fast and approximate dose calculation method is developed to calculate the dose for the IMPT plan under different setup and range uncertainties. Effects of two factors, inversed square factor and range uncertainty, are explored. WCS robustness analysis method was evaluated using this fast dose calculation method. The worst-case dose distribution was generated by shifting isocenter by 3 mm along x,y and z directions and modifying stopping power ratios by ±3.5%. 1000 randomly perturbed cases in proton range and x, yz directions were created and the corresponding dose distributions were calculated using this approximated method. DVH and dosimetric indexes of all 1000 perturbed cases were calculated and compared with the result using worst case scenario method. Results: The distributions of dosimetric indexes of 1000 perturbed cases were generated and compared with the results using worst case scenario. For D95 of CTVs, at least 97% of 1000 perturbed cases show higher values than the one of worst case scenario. For D5 of CTVs, at least 98% of perturbed cases have lower values than worst case scenario. Conclusion: By extensively calculating the dose distributions under random uncertainties, WCS method was verified to be reliable in evaluating the robustness level of MFO IMPT plans of H&N patients. The extensively sampling approach using fast approximated method could be used in evaluating the effects of different factors on the robustness level of IMPT plans in the future.

60 APPLIED LIFE SCIENCES↗

Evaluation of adsorption and mechanical strength of 13X zeolite mixtures with phyllosilicate binders using molecular dynamics simulation and positron annihilation spectroscopy

There is growing interest in developing zeolites with novel internal structures that have optimal adsorptive capacity and high mechanical strength, while offering advantages, such as being light weight. We integrate computational and experimental methods to explore the effect of binder/zeolite types, and weight percentages on the mechanical strength of 13X zeolite and adsorption capacities of N 2 , H 2 O, and CO 2 for additive manufacturing (AM) applications with the goal of maximizing both adsorption and strength. Zeolite 13X mixtures and phyllosilicate binders (either bentonite or kaolin) are combined using molecular dynamics (MD) simulations to create structures with various binder/zeolite weight percentages. Further, adsorption capabilities and mechanical strength are assessed using the grand canonical Monte Carlo (GCMC) and ReaxFF modules, respectively. Our modeling shows that an optimized zeolite/binder ratio for N 2 adsorption is around 15 wt% for kaolin and roughly 10 wt% for bentonite. The resulting parameters can be applied to facilitate macro-scale computational fluid dynamics (CFD) and finite element method (FEM) simulations of an AM zeolite structure. We also performed Positron Annihilation Lifetime Spectroscopy (PALS) measurements on zeolite samples to explore the effect of changes in the internal volume. The results show an inverse relationship between the free volume and the solid loading and adding a binder changes the morphology of the zeolite-binder compound and decreases open-volume area significantly.

36 MATERIALS SCIENCE↗

Design of Materials with Alchemite

Machine learning models that establish the relationships between materials processing and properties can enable inverse design of materials through active learning. Alchemite is a commercial software that can perform inverse materials design on sparse data. Here we evaluate Alchemite’s performance on a dataset of shape memory alloys and a dataset of heat exchangers compared to baseline random forest models. Alchemite had higher accuracy when making predictions on sparse data and was more accurate or nearly as accurate as random forests on complete datasets while also quantifying uncertainty. The software was also used to suggest processing steps and design parameters to optimize properties and performance; however, physical validation of the suggested design parameters was beyond the scope of this work. Several useful design insights were gained about the impact of the design parameters on properties and performance including the importance of dopant choice and amount for shape memory alloys and the importance of height and weight on the thermal resistance of heat exchangers.

Machine learning↗

Evaluating the feasibility of using downwind methods to quantify point source oil and gas emissions using continuously monitoring fence-line sensors

The dependable reporting of methane (CH 4 ) emissions from point sources, such as fugitive leaks from oil and gas infrastructure, is important for profit maximization (retaining more hydrocarbons), evaluating climate impacts, assessing CH 4 fees for regulatory programs, and validating CH 4 intensity in differentiated gas programs. Currently, there are disagreements between emissions reported by different quantification techniques for the same sources. It has been suggested that downwind CH 4 quantification methods using CH 4 measurements on the fence line of production facilities could be used to generate emission estimates from oil and gas operations at the site level, but it is currently unclear how accurate the quantified emissions are. To investigate the accuracy of downwind methods, this study uses fence-line simulated data collected during controlled-release experiments as input for a non-standard closed-path eddy covariance (EC), the Gaussian plume inverse model (GPIM), and the backward Lagrangian stochastic (bLs) model in a range of atmospheric conditions. This study's EC attempt was unsuccessful due to data collection and instrumentation issues, resulting in invalid results characterized by underestimated emissions, large negative fluxes, and cospectra/ogives that deviated from their ideal shapes. Consequently, the EC results could not be compared with the GPIM and bLS model. The bLs model demonstrated the highest accuracy for single-release single-point emissions, though it exhibited greater uncertainty than GPIM under multi-release conditions. Across the GPIM and bLs model, the most reliable quantification was achieved with 15 min averaging and a narrow 5° wind sector range. Although EC was limited in this context, future studies should consider employing a standard EC system and further optimizing GPIM and bLs approaches – particularly for complex multi-source scenarios – to enhance quantification accuracy and reduce uncertainty.

03 NATURAL GAS↗

Near-optimal ground state preparation

Preparing the ground state of a given Hamiltonian and estimating its ground energy are important but computationally hard tasks. However, given some additional information, these problems can be solved efficiently on a quantum computer. We assume that an initial state with non-trivial overlap with the ground state can be efficiently prepared, and the spectral gap between the ground energy and the first excited energy is bounded from below. With these assumptions we design an algorithm that prepares the ground state when an upper bound of the ground energy is known, whose runtime has a logarithmic dependence on the inverse error. When such an upper bound is not known, we propose a hybrid quantum-classical algorithm to estimate the ground energy, where the dependence of the number of queries to the initial state on the desired precision is exponentially improved compared to the current state-of-the-art algorithm proposed in [Ge et al. 2019]. These two algorithms can then be combined to prepare a ground state without knowing an upper bound of the ground energy. We also prove that our algorithms reach the complexity lower bounds by applying it to the unstructured search problem and the quantum approximate counting problem.

97 MATHEMATICS AND COMPUTING↗

Optimization of Ocean Color Algorithms: Application to Satellite Data Merging

The objective of the program is to develop and validate a procedure for ocean color data merging, which is one of the major goals of the SIMBIOS project. As part of the SIMBIOS Program, we have developed a merging method for ocean color data. Conversely to other methods our approach does not combine end-products like the subsurface chlorophyll concentration (chl) from different sensors to generate a unified product. Instead, our procedure uses the normalized water-leaving radiances L((sub wN)(lambda)) from single or multiple sensors and uses them in the inversion of a semi-analytical ocean color model that allows the retrieval of several ocean color variables simultaneously. Beside ensuring simultaneity and consistency of the retrievals (all products are derived from a single algorithm), this model-based approach has various benefits over techniques that blend end-products (e.g. chlorophyll): 1) It works with single or multiple data sources regardless of their specific bands; 2) It exploits band redundancies and band differences; 3) It accounts for uncertainties in the L((sub wN)(lambda)) data; 4) It provides uncertainty estimates for the retrieved variables.

Ritorena, Stephane↗

Management Trends and Outcomes for Stage I to II Mantle Cell Lymphoma Using the National Cancer Data Base: Ascertaining the Ideal Treatment Paradigm

Purpose: Mantle cell lymphoma (MCL) is a rare, albeit aggressive subset of non-Hodgkin lymphoma, resulting in varied treatment approaches. Given the paucity of data defining the optimal management for early-stage MCL, we conducted an analysis using the National Cancer Data Base (NCDB) to identify practice patterns and outcomes. Methods and Materials: The NCDB was queried for patients with stage I to II MCL diagnosed from 1998 to 2012 receiving chemotherapy (CT) or radiation therapy (RT), or both (CT+RT). Univariate and multivariable analyses for factors associated with treatment selection were completed using logistic regression. Propensity scores with inverse probability treatment weighting (IPTW) were calculated based on the conditional probability of receiving CT+RT. The log-rank test and Cox proportional hazards modeling with IPTW adjustment were conducted for the survival analyses. Results: In total, 2539 patients were identified. The key characteristics were as follows: 69% were male, 71% were aged ≥60 years, 28% had extranodal involvement, and 51% had stage I disease. Of the 2539 patients, 70% underwent CT, 11% underwent RT, and 19% underwent CT+RT. The use of CT+RT decreased from 23.1% to 14.1% in 1998 to 2002 and 2010 to 2012 (P<.001). CT+RT usage was lower for patients with the following characteristics: age ≥60 years, female sex, stage II disease, and the presence of B symptoms. With a median follow-up period of 42.8 months, the unadjusted 3-year overall survival estimates for patients receiving CT, RT, or CT+RT were 67.8%, 72.4%, and 79.8%, respectively (P<.001). After correcting for indication bias through IPTW-adjusted modeling, CT+RT reduced the risk of overall mortality compared with monotherapy (hazard ratio 0.65, P=.029). Conclusions: Although uncommon, patients with stage I-II MCL can have favorable outcomes. Despite a continued decline in the usage of consolidative RT, combined modality therapy improves survival in this cohort compared with monotherapy.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Performance of a Quintuple-GEM Based RICHDetector Prototype

Cerenkov technology is often the optimal choice for particle identification in high energy particle collision applications. Typically, the most challenging regime is at high pseudorapidity (forward) where particle identification must perform well at high laboratory momenta. For the upcoming electron ion collider (EIC), the physics goals require hadron (π, K, p) identification up to ~50 GeV/c. In this region Cerenkov ring-imaging (RICH) is the most viable solution. The speed of light in a radiator medium is inversely proportional to the refractive index. Hence, for particle identification (PID) reaching out to high momenta a small index of refraction is required. Unfortunately, the lowest indices of refraction also result in the lowest light yield ([(dNγ)/dx] ∝ sin2(θC)) driving up the radiator length and thereby the overall detector cost. In this paper we report on a successful test of a compact RICH detector (1 meter radiator) capable of delivering in excess of 10 photoelectrons per ring with a low index radiator gas (CF4). The detector concept is a natural extension of the PHENIX hadron-blind detector (HBD) achieved by adding focusing capability at low wavelength and adequate gain for high efficiency detection of single-electron induced avalanches. Our results indicate that this technology is indeed a viable choice in the forward direction of the EIC. The setup and results are described within.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Optical amplifiers for coherent lidar

We examine application of optical amplification to coherent lidar for the case of a weak return signal (a number of quanta of the return optical field close to unity). We consider the option that has been explored to date, namely, incorporation of an optical amplifier operated in a linear manner located after reception of the signal and immediately prior to heterodyning and photodetection. We also consider alternative strategies where the coherent interaction, the nonlinear processes, and the amplification are not necessarily constrained to occur in the manner investigated to date. We include the complications that occur because of mechanisms that occur at the level of a few, or one, quantum excitation. Two factors combine in the work to date that limit the value of the approach. These are: (1) the weak signal tends to require operation of the amplifier in the linear regime where the important advantages of nonlinear optical processing are not accessed, (2) the linear optical amplifier has a -3dB noise figure (SN(out)/SN(in)) that necessarily degrades the signal. Some improvement is gained because the gain provided by the optical amplifier can be used to overcome losses in the heterodyned process and photodetection. The result, however, is that introduction of an optical amplifier in a well optimized coherent lidar system results in, at best, a modest improvement in signal to noise. Some improvement may also be realized on incorporating more optical components in a coherent lidar system for purely practical reasons. For example, more compact, lighter weight, components, more robust alignment, or more rapid processing may be gained. We further find that there remain a number of potentially valuable, but unexplored options offered both by the rapidly expanding base of optical technology and the recent investigation of novel nonlinear coherent interference phenomena occurring at the single quantum excitation level. Key findings are: (1) insertion of linear optical amplifiers in well optimized conventional lidar systems offers modest improvements, at best, (2) the practical advantages of optical amplifiers, especially fiber amplifiers, such as ease of alignment, compactness, efficiency, lightweight, etc., warrant further investigation for coherent lidar, (3) the possibility of more fully optical lidar systems should be explored, (4) advantages gained by use of coherent interference of optical fields at the level of one, or a few, signal quanta should be explored, (5) amplification without inversion, population trapping, and use of electromagnetic induced transparency warrant investigation in connection with coherent lidar, (6) these new findings are probably more applicable to earth related NASA work, although applications to deep space should not be excluded, and (7) our own work in the Ultrafast Laboratory at UAH along some of the above lines of investigation, may be useful.

Fork, Richard↗