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 289 records · Page 16

Simultaneous Retrieval of Aerosol and Surface Optical Properties from Combined Airborne- and Ground-Based Direct and Diffuse Radiometric Measurements

This paper presents a new method for simultaneously retrieving aerosol and surface reflectance properties from combined airborne and ground-based direct and diffuse radiometric measurements. The method is based on the standard Aerosol Robotic Network (AERONET) method for retrieving aerosol size distribution, complex index of refraction, and single scattering albedo, but modified to retrieve aerosol properties in two layers, below and above the aircraft, and parameters on surface optical properties from combined datasets (Cloud Absorption Radiometer (CAR) and AERONET data). A key advantage of this method is the inversion of all available spectral and angular data at the same time, while accounting for the influence of noise in the inversion procedure using statistical optimization. The wide spectral (0.34-2.30 m) and angular range (180 ) of the CAR instrument, combined with observations from an AERONET sunphotometer, provide sufficient measurement constraints for characterizing aerosol and surface properties with minimal assumptions. The robustness of the method was tested on observations made during four different field campaigns: (a) the Southern African Regional Science Initiative 2000 over Mongu, Zambia, (b) the Intercontinental Transport Experiment-Phase B over Mexico City, Mexico (c) Cloud and Land Surface Interaction Campaign over the Atmospheric Radiation Measurement (ARM) Central Facility, Oklahoma, USA, and (d) the Arctic Research of the Composition of the Troposphere from Aircraft and Satellites (ARCTAS) over Elson Lagoon in Barrow, Alaska, USA. The four areas are dominated by different surface characteristics and aerosol types, and therefore provide good test cases for the new inversion method.

Gatebe, C. K.↗

Propulsion-Airframe Integration for Conceptual Redesign of a Low-Boom Supersonic Transport

A low-boom supersonic transport was designed for a cruise Mach of 1.7 and 40 passengers. This low-boom aircraft, referred to as the Mach 1.7 40-PAX concept, was generated using computational fluid dynamics (CFD) based sonic boom analysis at the start of overland cruise (SOC). The engine for the Mach 1.7 40-PAX concept was designed using the Numerical Propulsion System Simulation (NPSS) and modeled as a flow-through nacelle for CFD analysis. To understand how the engine plume affects the undertrack ground signature, the Mach 1.7 40-PAX concept is redesigned using an aeropropulsive CFD simulation. A process is developed for approximation of the NPSS engine at SOC by a CFD engine for aeropropulsive CFD simulation. The generated CFD engine has the identical nozzle boundary conditions and approximately the same mass flow and thrust as those of the NPSS engine at SOC. Then, the configuration with the CFD engines is optimized to approximately restore the low-boom characteristics of the Mach 1.7 40-PAX concept. Finally, the CFD simulation data for the optimized concept with the CFD engines is used to calibrate the low-fidelity aerodynamic analyses for mission analysis of this concept. The cyclic dependency of the involved disciplinary analyses is resolved using an iteration method for a consistent coupling of the mission analysis, NPSS engine analysis, and low-boom redesign using the aeropropulsive CFD simulation. This low-boom redesign study is used as an example to demonstrate how the propulsion-airframe integration could be implemented for conceptual design of low-boom supersonic transports.

Low-boom supersonic transport↗

Propulsion-Airframe Integration for Conceptual Redesign of a Low-Boom Supersonic Transport

A low-boom supersonic transport was designed for a cruise Mach of 1.7 and 40 passengers. This low-boom aircraft, referred to as the Mach 1.7 40-PAX concept, was generated using computational fluid dynamics (CFD) based sonic boom analysis at the start of overland cruise (SOC). The engine for the Mach 1.7 40-PAX concept was designed using the Numerical Propulsion System Simulation (NPSS) and modeled as a flow-through nacelle for CFD analysis. To understand how the engine plume affects the undertrack ground signature, the Mach 1.7 40-PAX concept is redesigned using an aeropropulsive CFD simulation. A process for approximation of the NPSS engine at SOC by a CFD engine is developed. The generated CFD engine has the identical nozzle boundary conditions and approximately the same mass flow and thrust as those of the NPSS engine at SOC. Then, the configuration with the CFD engines is optimized to approximately restore the low-boom characteristics of the Mach 1.7 40-PAX concept. Finally, the CFD simulation data for the optimized concept with the CFD engines is used to calibrate the low-fidelity aerodynamic analyses for mission analysis of this concept. The cyclic dependency of the involved disciplinary analyses is resolved using an iteration method for a consistent coupling of the mission analysis, NPSS engine analysis, and low-boom redesign using the aeropropulsive CFD simulation. This low-boom redesign study is used as an example to demonstrate how the propulsion-airframe integration could be implemented for conceptual design of low-boom supersonic transports.

Low-boom supersonic transport↗

Optimal Mission Profile for Low-Boom Supersonic Aircraft

An optimal mission profile, called the low-boom mission profile, is proposed for low-boom supersonic aircraft. The low-boom mission profile uses a fixed AoA for the cruise segment until reaching the cruise ceiling. The fixed AoA is determined by the cruise condition for the low-boom shaping of the aircraft. The derived CL trajectory is not necessarily the optimal trajectory for the cruise range when unrestricted supersonic flight is allowed. The undertrack sonic boom ground noise level is reduced during the cruise because of the increasing altitude from the fixed AoA trajectory and the approximate invariance of nondimensional overpressure at 3BL below a supersonic aircraft with respect to cruise altitude change when the AoA and cruise Mach are fixed. The approximate invariance of nondimensional overpressure at 3BL below a supersonic aircraft is verified using both Euler and Reynolds-averaged Navier-Stokes CFD solvers. A low-boom supersonic transport concept is used to illustrate that, compared to the fixed cruise altitude trajectory, the low-boom mission profile can increase the cruise range of a low-boom aircraft and reduce the undertrack sonic boom ground noise levels for the entire cruise segment. One significant contribution of this paper is the implication that an aircraft OML satisfying a sonic boom ground noise level constraint at SOC will satisfy the same noise level constraint over the entire cruise segment, as long as a fixed AoA trajectory is flown.

Low-boom characteristics↗

Inference of phase field fracture models

The phase field approach to modeling fracture uses a diffuse damage field to represent cracks. This representation mollifies singularities that arise in computations with sharp interface models and some of the resultant difficulties in the mathematical and numerical treatment of fracture. Phase field fracture models have proven effective at representing crack propagation, branching, and merging. Specific formulations, beginning with brittle fracture, have also been shown to converge to classical solutions. Extensions to cover the range of material failure, including ductile and cohesive fracture, lead to an array of possible models. There exists a large body of literature focusing on this class of models and on the impact of model form on the predicted crack evolution. However, there have not been systematic studies into how optimal models may be chosen. Here, we take a first step in this direction by developing formal methods for identification of the best parsimonious model of phase field fracture given full-field data on the damage and deformation fields. We consider some of the main models that have been used for the degradation of elastic response due to damage and its propagation. Our approach builds upon Variational System Identification (VSI), a weak form variant of the Sparse Identification of Nonlinear Dynamics (SINDy). Furthermore, in this first communication we focus on synthetically generated data but we also consider central issues associated with the use of experimental full-field data, such as data sparsity and noise.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Optimizing hierarchical membrane/catalyst systems for oxidative coupling of methane using additive manufacturing

The advantage of a membrane/catalyst system in oxidative coupling of methane (OCM) compared to conventional reactive systems is that by introducing oxygen to the OCM catalytic sites through a membrane, parasitic gas phase reactions of O 2 (g), responsible for lowering product selectivity, can be avoided. The design and fabrication of membrane/catalyst systems has, however, been hampered by low volumetric chemical conversion rates, high capital cost, and difficulties in codesigning membrane and catalyst properties to optimize the performance. We solve these issues by developing a dual-layer additive manufacturing process, based on phase inversion, to design, fabricate and optimize a hollow fiber membrane/catalyst system for OCM. We demonstrate the approach though a case study using BaCe 0.8 Gd 0.2 O 3-δ (BCG) as the basis of both the catalyst and separation layers. We show that by using the manufacturing approach we can codesign the membrane thickness and catalyst surface area so that the flux of oxygen transport through the membrane and methane activation rates in the catalyst layer match each other. Here, we demonstrate that this “rate matching” is critical for maximizing the performance, with the membrane/catalyst system significantly overperforming conventional reactor designs under identical conditions.

36 MATERIALS SCIENCE↗

Propulsion-Airframe Integration for Conceptual Redesign of a Low-Boom Supersonic Transport

A low-boom supersonic transport was generated in a previous low-boom multidisciplinary optimization (MDO) study that used computational fluid dynamics (CFD) off-body pressure to compute the undertrack sonic boom ground signature and calibrated low-fidelity aerodynamics analyses to compute the mission performance metrics. This low-boom aircraft, referred to as the Mach 1.7 40-PAX concept, can carry 40 passengers for a low-boom overland mission with cruise Mach 1.7 for the airport pairs over the continental US and has the potential to achieve an undertrack sonic boom ground noise level below 70 perceived level of decibels at the start of overland cruise (SOC). The engine for the Mach 1.7 40-PAX concept was designed using the Numerical Propulsion System Simulation (NPSS) for the mission analysis and modeled as a flow-through nacelle for CFD-based sonic boom analysis. To understand how the engine plume affects the undertrack ground signature, the Mach 1.7 40-PAX concept is redesigned in this paper after replacing the flow-through nacelles by CFD engines for sonic boom analysis using aeropropulsive CFD simulation. A process is developed for approximation of the NPSS engine at SOC by a CFD engine for aeropropulsive CFD simulation. The generated CFD engine has the identical nozzle boundary conditions and approximately the same mass flow and thrust as those of the NPSS engine at SOC. Then, the outer mold line of the configuration with the CFD engines is redesigned to approximately restore the low-boom characteristics of the Mach 1.7 40-PAX concept at SOC. Finally, the CFD simulation data for the redesigned concept with the CFD engines is used to calibrate the low-fidelity aerodynamic analyses for mission analysis of this concept. The cyclic dependency of the involved disciplinary analyses is resolved using an iteration method for a consistent coupling of the mission analysis, NPSS engine analysis, and low-boom redesign using the aeropropulsive CFD simulation. This low-boom redesign study is used as an example to demonstrate how the propulsion-airframe integration could be implemented for conceptual design of supersonic transports that satisfy both the low-boom and mission performance requirements.

Low-boom supersonic transport↗

A staged deep learning approach to spatial refinement in 3D temporal atmospheric transport

High-resolution spatiotemporal simulations effectively capture the complexities of atmospheric plume dispersion in complex terrain. However, their high computational cost makes them impractical for applications requiring rapid responses or iterative processes, such as optimization, uncertainty quantification, or inverse modeling. To address this challenge, this work introduces the Dual-Stage Temporal Three-dimensional UNet Super-resolution (DST3D-UNet-SR) model, a highly efficient deep learning model for plume dispersion predictions. DST3D-UNet-SR is composed of two sequential modules: the temporal module (TM), which predicts the transient evolution of a plume in complex terrain from low-resolution temporal data, and the spatial refinement module (SRM), which subsequently enhances the spatial resolution of the TM predictions. We train DST3D-UNet-SR using a comprehensive dataset derived from high-resolution large eddy simulations (LES) of plume transport. We propose the DST3D-UNet-SR model to significantly accelerate LES of three-dimensional (3D) plume dispersion by three orders of magnitude. Additionally, the model demonstrates the ability to dynamically adapt to evolving conditions through the incorporation of new observational data, substantially improving prediction accuracy in high-concentration regions near the source.

3D temporal sequences↗

LaSDI: Parametric Latent Space Dynamics Identification

Enabling fast and accurate physical simulations with data has become an important area of computational physics to aid in inverse problems, design-optimization, uncertainty quantification, and other various decision-making applications. This paper presents a data-driven framework for parametric latent space dynamics identification procedure that enables fast and accurate simulations. The parametric model is achieved by building a set of local latent space model and designing an interaction among them. An individual local latent space dynamics model achieves accurate solution in a trust region. By letting the set of trust region to cover the whole parameter space, our model shows an increase in accuracy with an increase in training data. Herein we introduce two different types of interaction mechanisms, i.e., point-wise and region-based approach. Both linear and nonlinear data compression techniques are used. We illustrate the framework of Latent Space Dynamics Identification (LaSDI) enable a fast and accurate solution process on various partial differential equations, i.e., Burgers’ equations, radial advection problem, and nonlinear heat conduction problem, achieving 0 (100)x speed-up and 0 (1)% relative error with respect to the corresponding full order models.

97 MATHEMATICS AND COMPUTING↗

Evaluation of dual-weighted residual and machine learning error estimation for projection-based reduced-order models of steady partial differential equations

Projection-based reduced-order models (pROMs) show great promise as a means to accelerate many-query applications such as forward error propagation, solving inverse problems, and design optimization. In order to deploy pROMs in the context of high-consequence decision making, accurate error estimates are required to determine the region(s) of applicability in the parameter space. The following paper considers the dual-weighted residual (DWR) error estimate for pROMs and compares it to another promising pROM error estimate, machine learned error models (MLEM). Here, we show how DWR can be applied to ROMs and then evaluate DWR on two partial differential equations (PDEs): a two-dimensional linear convection–reaction–diffusion equation, and a three-dimensional static hyper-elastic beam. It is shown that DWR is able to estimate errors for pROMs extrapolating outside of their training set while MLEM is best suited for pROMs used to interpolate within the pROM training set.

42 ENGINEERING↗

Sensitivity enhancement of surface thermal lens technique with a short-wavelength probe beam: Experiment

Surface thermal lens is a highly sensitive photothermal technique to measure low absorption losses of various solid materials. In such applications, the sensitivity of surface thermal lens is a key parameter for measuring extremely low absorption. In this paper, we experimentally investigated the influence of probe beam wavelength on the sensitivity of surface thermal lens for measuring the low absorptance of optical laser components. Three probe lasers with wavelength 375 nm, 633 nm, and 1570 nm were used, respectively, to detect the surface thermal lens amplitude of a highly reflective coating sample excited by a cw modulated Gaussian beam at 1064 nm. The experimental results showed that the maximum amplitude of surface thermal lens signal obtained at corresponding optimized detection distance was inversely proportional to the wavelength of the probe beam, as predicted by previous theoretical model. The sensitivity of surface thermal lens could, therefore, be improved by detecting surface thermal lens signal with a short-wavelength probe beam.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

TRINIDI: Time-of-Flight Resonance Imaging With Neutrons for Isotopic Density Inference

Accurate reconstruction of 2D and 3D isotope densities is a desired capability with great potential impact in applications such as evaluation and development of next-generation nuclear fuels. Neutron time-of-flight (TOF) resonance imaging offers a potential approach by exploiting the characteristic neutron absorption spectra of each isotope. However, it is a major challenge to compute quantitatively accurate images due to a variety of confounding effects such as severe Poisson noise, background scatter, beam non-uniformity, absorption non-linearity, and extended source pulse duration. We present the TRINIDI algorithm which is based on a two-step process in which we first estimate the neutron flux and background counts, and then reconstruct the areal densities of each isotope and pixel. Both components are based on the inversion of a forward model that accounts for the highly non-linear absorption, energy-dependent emission profile, and Poisson noise, while also modeling the substantial spatio-temporal variation of the background and flux. Further, to do this, we formulate the non-linear inverse problem as two optimization problems that are solved in sequence. We demonstrate on both synthetic and measured data that TRINIDI can reconstruct quantitatively accurate 2D views of isotopic areal density that can then be reconstructed into quantitatively accurate 3D volumes of isotopic volumetric density.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

SU-G-201-14: Is Maximum Skin Dose a Reliable Metric for Accelerated Partial Breast Irradiation with Brachytherapy?

Purpose: To evaluate the reliability of the maximum point dose (Dmax) to the skin surface as a dosimetric constraint, we investigated the correlation between Dmax at the skin surface and dose metrics at various definitions of skin thickness. Methods: 42 patients treated with APBI using a Strut Adjusted Volume Implant (SAVI) applicator between 2010 and 2014 were retrospectively reviewed. Target (PTV-EVAL) and organs at risk (OARs: skin, lung, and ribs) were delineated on a CT following NSABP B-39 guidelines. Six skin structures were contoured: a rind 3cm external to the body surface and 1, 2, 3, 4, and 5mm thick rinds deep to the body surface. Inverse planning simulated annealing optimization was used to deliver 32–34Gy in 8-10 fractions to the target while minimizing OAR doses. Dmax, D0.1cc, D1.0cc, and D2.0cc to the various skin structures were calculated. Linear regressions between the metrics were evaluated using the coefficient of determination (R{sup 2}). Results: The average±SD PTV-EVAL volume and cavity-to-skin distances were 71.1±28.5cc and 6.9±5.0mm. The target V90 and V95 were 97.3±2.3% and 95.1±3.2%. The Dmax to the skin structures were 78.7±10.2% (skin surface), 82.2±10.7% (skin-1mm), 89.4±12.6% (skin-2mm), 97.9±15.4% (skin-3mm), 114.1±32.5% (skin-4mm), and 157.0±85.3% (skin-5mm). Linear regression analysis showed D1.0cc and D2.0cc to the skin 1mm and Dmax to the skin-4mm and 5mm were poorly correlated with other metrics (R{sup 2}=0.413±0.204). Dmax to the skin surface was well correlated (R{sup 2}=0.910±0.047) and D1.0cc to the skin-3mm was strongly correlated with all subsurface skin layers (R{sup 2}=0.935±0.050). Conclusion: Dmax to the skin surface is a relevant metric for breast skin dose. Contouring discontinuities in the skin with a 1mm subsurface rind and the active dwells in the skin 4 and 5mm introduced significant variations in skin DVH. D0.1cc, D1.0cc, and D2.0cc to a 3mm skin rind are more robust metrics in breast brachytherapy.

60 APPLIED LIFE SCIENCES↗

Rapid Optimization Library

ROL provides interfaces to and implementations of algorithms for gradient-based unconstrained and constrained optimization. ROL can be used to optimize the response of any client simulation code that evaluates scalar-valued response functions. If the client code can provide gradient information for the response function, ROL will take advantage of it, resulting in faster runtimes. ROL's interfaces are matrix-free, in other words ROL only uses evaluations of scalar-valued and vector-valued functions. ROL can be used to solve optimal design problems and inverse problems based on a variety of simulation software.

Kouri, DrewKouri↗

High-Dimensional Spectral Sampling

A wide variety of applications in science and engineering depend on effectively and efficiently exploring high- dimensional parameter spaces. Common examples are uncertainty quantification, inverse problems, or various optimization approaches. In such applications, the first step is to create an initial uniform, random sampling of the space to create a baseline of knowledge. The simulations evaluated at these samples are then used for subsequent processing. In its most generic form, the goal of sampling is to produce the maximal amount of information with the minimal number of samples. However, despite being the first step and thus influencing the performance of the entire workflow, the quality of the initial sample designs is rarely analyzed or optimized. This project develops a new spectral sampling theory for analyzing and creating experiment designs in high-dimensional spaces. Further, new algorithms for estimating sample properties and sample synthesis are designed. Finally, applications of the proposed spectral sampling tools in surrogate modeling, scientific machine learning and image analysis are presented.

97 MATHEMATICS AND COMPUTING↗

Characterizing Baselines and Change in Gas Hydrate Systems using EM Methods

The objective of this project was to advance our understanding of gas hydrate systems in nature by characterizing their electrical properties in the field and in the laboratory. In the laboratory measurements, methane hydrate was synthesized from pure water ice and flash frozen seawater, with varying amounts of sand or silt added. Electrical conductivity was determined by impedance spectroscopy, using equivalent circuit modeling to separate the effects of electrodes and to gain insight into conduction mechanisms. Silt and sand increase the conductivity of pure hydrate, inferred to be contaminant NaCl contributing to conduction in hydrate, to a peak conductivity in agreement with peak resistivities observed in well logs through massive hydrate (3,000--10,000 Ωm). The addition of silt and sand lowers the conductivity of hydrate synthesized from seawater, by an amount consistent with Archie's Law. All samples were characterized using cryogenic scanning electron microscopy and energy dispersive spectroscopy, which shows good connectivity of salt and brine phases. Electrical conductivity measurements of pure hydrate and hydrate mixed with silt during pressure-induced dissociation supports previous conclusions that sediment increases dissociation rate. In order to characterize gas hydrate systems in the field, we collected 360 line kilometers of controlled-source electromagnetic data on Walker Ridge 313, Orca Basin (WR100), Mad Dog (GC781), and Green Canyon 955 in the Gulf of Mexico, all areas with known or seismically inferred gas hydrate deposits and which have be drilled or targeted for future drilling. We deep-towed an EM transmitter that generates an alternating electric field which propagates through the seafloor geology. Data were recorded on 6 receivers towed behind the transmitter at distances between 550 and 1550 m. In the presence of conductive geology, the electric fields will be attenuated, and conversely, in resistive geology the fields will be preserved. Our data were inverted using a 2D inversion method that first optimizes the model-data misfit, then finds the smoothest model fitting the data. This ensures that resistivity structures present in the final model are likely necessary. At each of the proposed drilling sites we found increased resistivity, interpreted as increased hydrate concentrations. However, not only were the primary drilling sites not always more resistive than the alternate sites, at WR313 the strongest resistors were not at the locations targeted for drilling.

03 NATURAL GAS↗

Geometric Scale-up Experiments on Fluidization of Geldart B Glass Beads

The objective of this work is to provide a valuable database from controlled experiments for validating computational models. Recently, coarse-grained techniques such as particle-in-cell (PIC) or coarse-grained discrete element modeling (DEM) have gained popularity due to their computational efficiency while modeling large-scale systems; however, the influence of model parameters and their sensitivities at different geometric scales and flow conditions remain to be analyzed. These datasets are critical for the multiphase flow research community to assess predictive capability of modeling techniques as well as elucidate the hydrodynamic behavior in these systems. This study performed fluidization experiments using three different test sections with internal diameters of 2.5, 4, and 6 in. The operating conditions, bed material, and range of flow velocities at the inlet were constant in all the units, which were not hydrodynamically scaled. Glass beads having a Sauter Mean Diameter of 332 μm were used. Superficial velocity was varied from 2.97 to 5.35 times the minimum fluidization velocity and the initial static bed height was 0.1524 m. The order in which the experiments were performed was randomized and replicates were included to provide uncertainty in measurements. Statistics of differential pressure and bed height from these tests were reported. Future plans include validating PIC methodology in the open-source software, MFiX (Multiphase Flow with Interphase Exchanges) using results from this study. This could further be extended to determine optimal model parameters using inverse techniques such as deterministic calibration or Bayesian inference.

20 FOSSIL-FUELED POWER PLANTS↗

Geometric Scale-up Experiments on Fluidization of Geldart B Glass Beads

The objective of this work is to provide a valuable database from controlled experiments for validating computational models. Recently, coarse-grained techniques such as particle-in-cell (PIC) or coarse-grained discrete element modeling (DEM) have gained popularity due to their computational efficiency while modeling large-scale systems; however, the influence of model parameters and their sensitivities at different geometric scales and flow conditions remain to be analyzed. These datasets are critical for the multiphase flow research community to assess predictive capability of modeling techniques as well as elucidate the hydrodynamic behavior in these systems. This study performed fluidization experiments using three different test sections with internal diameters of 2.5, 4, and 6 in. The operating conditions, bed material, and range of flow velocities at the inlet were constant in all the units, which were not hydrodynamically scaled. Glass beads having a Sauter Mean Diameter of 332 μm were used. Superficial velocity was varied from 2.97 to 5.35 times the minimum fluidization velocity and the initial static bed height was 0.1524 m. The order in which the experiments were performed was randomized and replicates were included to provide uncertainty in measurements. Statistics of differential pressure and bed height from these tests were reported. Future plans include validating PIC methodology in the open-source software, MFiX (Multiphase Flow with Interphase Exchanges) using results from this study. This could further be extended to determine optimal model parameters using inverse techniques such as deterministic calibration or Bayesian inference.

42 ENGINEERING↗