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 973 records · Page 54

Mapped Moments to a Cartesian Grid (MMCG) Value-Added Product Report

Objective analysis (OA) is a method of mapping unstructured data to a structured grid. In the context of scanning radar data, OA is used to interpolate data in antenna coordinates (range, azimuth, and elevation) onto a regularly spaced Cartesian grid (Trapp and Doswell 2000). The U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility’s Mapped Moments to a Cartesian Grid (MMCG) Value-Added Product (VAP) uses the Python ARM Radar Toolkit (Py-ART), a data model-driven interactive architecture for working with weather radar data, to map the data to a Cartesian grid (Helmus and Collis 2016). MMCG, with Py-ART built in, has the ability to take radar data in antenna coordinates and map the gates to a Cartesian grid using inverse distance weight functions such as Cressman (square) and Barnes (exponential), but also can filter the data during the interpolation. MMCG also allows arbitrary formulations for the radius of influence, which are matched to particular radar scanning strategies. This creates a complex parameter space for optimizing the retention of storm structure detail while minimizing artifacts. MMCG takes data processed with ARM’s Corrected Precipitation Radar Moments in Antenna Coordinates (CMAC) VAP and maps it to a Cartesian grid as the output product. A variety of fields that have been mapped to the Cartesian grid are then saved to plots to complement each grid file.

54 ENVIRONMENTAL SCIENCES↗

Applications of generalized inverses to estimation in dynamic systems.

Generalized matrix inverses are used to obtain an estimation procedure for estimation of the state vector of a dynamic system. This sequential procedure is studied analytically with respect to the choice of an arbitrary vector. The covariance matrix of the estimator is determined and compared to the optimal Kalman type procedure. A numerical example illustrates the procedure and compares it to the optimal one.

Prochaska, B. J.↗

A New Normalized Difference Cloud Retrieval Technique Applied to Landsat Radiances Over the Oklahoma ARM Site

We suggest a new approach to cloud retrieval, using a normalized difference of nadir reflectivities (NDNR) constructed from a non-absorbing and absorbing (with respect to liquid water) wavelength. Using Monte Carlo simulations we show that this quantity has the potential of removing first order scattering effects caused by cloud side illumination and shadowing at oblique Sun angles. Application of the technique to TM (Thematic Mapper) radiance observations from Landsat-5 over the Southern Great Plains site of the ARM (Atmospheric Radiation Measurement) program gives very similar regional statistics and histograms, but significant differences at the pixel level. NDNR can be also combined with the inverse NIPA (Nonlocal Independent Pixel Approximation) of Marshak (1998) which is applied for the first time on overcast Landsat scene subscenes. We demonstrate the sensitivity of the NIPA-retrieved cloud fields on the parameters of the method and discuss practical issues related to the optimal choice of these parameters.

Orepoulos, Lazaros↗

Viscous Design of TCA Configuration

The goal in this effort is to redesign the baseline TCA configuration for improved performance at both supersonic and transonic cruise. Viscous analyses are conducted with OVERFLOW, a Navier-Stokes code for overset grids, using PEGSUS to compute the interpolations between overset grids. Viscous designs are conducted with OVERDISC, a script which couples OVERFLOW with the Constrained Direct Iterative Surface Curvature (CDISC) inverse design method. The successful execution of any computational fluid dynamics (CFD) based aerodynamic design method for complex configurations requires an efficient method for regenerating the computational grids to account for modifications to the configuration shape. The first section of this presentation deals with the automated regridding procedure used to generate overset grids for the fuselage/wing/diverter/nacelle configurations analysed in this effort. The second section outlines the procedures utilized to conduct OVERDISC inverse designs. The third section briefly covers the work conducted by Dick Campbell, in which a dual-point design at Mach 2.4 and 0.9 was attempted using OVERDISC; the initial configuration from which this design effort was started is an early version of the optimized shape for the TCA configuration developed by the Boeing Commercial Airplane Group (BCAG), which eventually evolved into the NCV design. The final section presents results from application of the Natural Flow Wing design philosophy to the TCA configuration.

Krist, Steven E.↗

Performance of JAXA’s SGLI standard ocean color products for oceanic to coastal waters: chlorophyll a concentration and light absorption coefficients of colored dissolved organic matter

We present the performance of JAXA’s SGLI standard algorithms for estimating chlorophyll a (chl a) concentration and the light absorption coefcients of colored dissolved organic matter (CDOM) using recently compiled high-quality data obtained from oceanic to coastal waters. Prior to the evaluation of the algorithms, we frst compare two forward models (Gordon et al. in J Geophys Res 93(D9):10909-10924, 1988: G88 and Park and Ruddick in Appl Opt 44(7):1236–1249, 2005: PR05) using a wide range of inherent optical properties (IOPs) to (1) examine if the water refectance is appropriately reproduced and (2) correct measured refectance in terms of its bidirectionality. Based on a good reproducibility of water refectance using the PR05, the optimized IOPs are further used for explaining errors in estimates of chl a concentration and CDOM absorption when using the SGLI, the GSM (Maritorena et al. in Appl Opt 41:2705–2714, 2002), and the QAA (Lee et al. in Appl Opt 41:5755, 2002) inversion algorithms. Results show that the mean error for estimating chl a concentration using the SGLI algorithm is 110% for our dataset. Although this error is lower than that of the GSM and the QAA algorithms, a signifcant underestimate at chl a higher than 2.0 mg m−3 is observed, which is further shown by independent match-up analyses. Another SGLI CDOM product includes the mean error of 50% and shows deviation at high CDOM range (>1.0 per m). A similar trend is observed for the GSM but not for the QAA when a global relationship of CDOM to colored detrital matter is applied. The sources of errors are discussed for potentially improving the retrievals

Second-generation GLobal Imager (SGLI)↗

A general rough-surface inversion algorithm: Theory and application to SAR data

Rough-surface inversion has significant applications in interpretation of SAR data obtained over bare soil surfaces and agricultural lands. Due to the sparsity of data and the large pixel size in SAR applications, it is not feasible to carry out inversions based on numerical scattering models. The alternative is to use parameter estimation techniques based on approximate analytical or empirical models. Hence, there are two issues to be addressed, namely, what model to choose and what estimation algorithm to apply. Here, a small perturbation model (SPM) is used to express the backscattering coefficients of the rough surface in terms of three surface parameters. The algorithm used to estimate these parameters is based on a nonlinear least-squares criterion. The least-squares optimization methods are widely used in estimation theory, but the distinguishing factor for SAR applications is incorporating the stochastic nature of both the unknown parameters and the data into formulation, which will be discussed in detail. The algorithm is tested with synthetic data, and several Newton-type least-squares minimization methods are discussed to compare their convergence characteristics. Finally, the algorithm is applied to multifrequency polarimetric SAR data obtained over some bare soil and agricultural fields. Results will be shown and compared to ground-truth measurements obtained from these areas. The strength of this general approach to inversion of SAR data is that it can be easily modified for use with any scattering model without changing any of the inversion steps. Note also that, for the same reason it is not limited to inversion of rough surfaces, and can be applied to any parameterized scattering process.

Moghaddam, M.↗

Mechanistic modeling study on process optimization and precursor utilization with atmospheric spatial atomic layer deposition

Spatial atomic layer deposition (SALD) is a promising technology with the aim of combining the advantages of excellent uniformity and conformity of temporal atomic layer deposition (ALD), and an industrial scalable and continuous process. In this manuscript, an experimental and numerical combined model of atmospheric SALD system is presented. To establish the connection between the process parameters and the growth efficiency, a quantitative model on reactant isolation, throughput, and precursor utilization is performed based on the separation gas flow rate, carrier gas flow rate, and precursor mass fraction. The simulation results based on this model show an inverse relation between the precursor usage and the carrier gas flow rate. With the constant carrier gas flow, the relationship of precursor usage and precursor mass fraction follows monotonic function. The precursor concentration, regardless of gas velocity, is the determinant factor of the minimal residual time. The narrow gap between precursor injecting heads and the substrate surface in general SALD system leads to a low Péclet number. In this situation, the gas diffusion act as a leading role in the precursor transport in the small gap rather than the convection. Fluid kinetics from the numerical model is independent of the specific structure, which is instructive for the SALD geometry design as well as its process optimization.

42 ENGINEERING↗

An efficient hybrid visible light communication/radio frequency system for vehicular applications

This paper focuses on improving the traffic system for Intelligent Transportation System (ITS) by using an imaging receiver instead of photodiode as a single receiver. Here, two simple traffic models have been proposed and analyzed in order to optimize the design characteristics such as signal to noise ratio, required power, received information, bit error rate, and modulation technique. First, a standalone visible light communication (VLC) with different modulation techniques is considered. The On–OFF Keying, L-Pulse Position Modulation (L-PPM), and Inverse L-Pulse Position Modulation (I-L-PPM) are investigated. Then, a hybrid communication system utilizing VLC and radio frequency (RF) is proposed for position-based services. The VLC/RF system combines the benefits of both systems and offers long distance transmission which is an important concern in ITS and in order to improve the link reliability in infrastructure to vehicle communications. This means that the probability of having a communication link between the transmitter and receiver is maintained without interruption. Then, the performance is finally evaluated. The obtained simulation results show a considerable increase in the received information using the proposed hybrid VLC/RF system compared to a VLC system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Photodiode-based machine learning for optimization of laser powder bed fusion parameters in complex geometries

We report the quality of parts produced through laser powder bed fusion additive manufacturing can be irregular, with complex geometries sometimes exhibiting dimensional inaccuracies and defects. For optimal part quality, laser process parameters should be selected carefully prior to printing and adjusted during the print if necessary. This is challenging since approaches to control and optimize the build parameters need to take into account the part geometry, the material, and the complex physics of laser powder bed fusion. This work describes a data-driven approach using experimental diagnostics for the optimization of laser process parameters prior to printing. A training dataset is generated by collecting high speed photodiode signal data while printing simple parts containing key geometry features with various process parameter strategies. Supervised learning approaches are employed to train both a forward model and an inverse model. The forward model takes as inputs track-wise geometry features and laser parameters and outputs the photodiode signal along the scan path. The inverse model takes as inputs the geometry features and photodiode signal and predicts the laser parameters. Given the part geometry and a desired photodiode signal, the inverse model can thus determine the required laser parameters. Two test parts which contain defect-prone features are used to assess the validity of the inverse model. The use of the model leads to improved part quality (higher dimensional accuracy, reduced dross, reduced distortion) for both test geometries.

36 MATERIALS SCIENCE↗

Divertor heat load estimates on NSTX and DIII-D using new and open-source 2D inversion analysis code

A thermography inversion algorithm has been developed in the open-source Python-based computer code, HYPERION, to calculate the heat flux incident on plasma-facing components (PFCs) in axisymmetric tokamaks. The chosen mesh size at the surface significantly affects the calculated transient heat flux results. The calculated transient heat flux will exceed the real value when the mesh size tends to zero but will underestimate the real value when the mesh size is large. A criterion for determining the appropriate mesh size for the transient heat flux calculation will be discussed. The numerical scheme for HYPERION uses a 2D fully implicit finite-difference approach, allowing temperature-dependent thermal properties of PFC materials. The inversion algorithm is benchmarked against established heat flux calculation codes, TACO and THEODOR, based on thermography data from NSTX and DIII-D respectively. The primary benefits of HYPERION compared to TACO and THEODOR are that it is open-source and it allows for the optimization of mesh thickness along the substrate. The algorithm also accounts for the thermal properties of thin surface layers that characteristically form on PFCs due to plasma-material interactions. The agreement between HYPERION and THEODOR is excellent, as the percent difference between the codes is ~5% on average in the case of the DIII-D data for moderate to high heat flux. Verification tests with TACO show slightly higher average percent differences of 8% and 12%. In using HYPERION to study filaments in heat flux, the initial results indicate that small ELMs filaments significantly broaden the divertor heat flux, and decrease divertor peak flux. Compared to the inter-ELM, the small ELM filaments decrease the divertor peak surface temperature. With intermittent divertor filaments, the divertor heat flux width is comparable with that found in L-mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Application and performance of an ML-EM algorithm in NEXT

The goal of the NEXT experiment is the observation of neutrinoless double beta decay in (136)Xe using a gaseous xenon TPC with electroluminescent amplification and specialized photodetector arrays for calorimetry and tracking. The NEXT Collaboration is exploring a number of reconstruction algorithms to exploit the full potential of the detector. This paper describes one of them: the Maximum Likelihood Expectation Maximization (ML-EM) method, a generic iterative algorithm to find maximum-likelihood estimates of parameters that has been applied to solve many different types of complex inverse problems. In particular, we discuss a bi-dimensional version of the method in which the photosensor signals integrated over time are used to reconstruct a transverse projection of the event. First results show that, when applied to detector simulation data, the algorithm achieves nearly optimal energy resolution (better than 0.5% FWHM at the Q value of (136)Xe) for events distributed over the full active volume of the TPC.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Infrared Observations of the Orion Capsule During EFT-1 Hypersonic Reentry

High-resolution infrared observations of the Orion capsule during its atmospheric reentry on December 5, 2015 were made from a US Navy NP-3D. This aircraft, equipped with a long-range optical sensor system, tracked the capsule from Mach 10 to 7 from a distance of approximately 60 nmi. Global surface temperatures of the capsule's thermal heatshield were derived from near infrared intensity measurements. The global surface temperature measurements complemented onboard instrumentation and were invaluable to the interpretation of the in-depth thermocouple measurements which rely on inverse heat transfer methods and material response codes to infer the desired surface temperature from the sub-surface measurements. The full paper will address the motivations behind the NASA Engineering Safety Center sponsored observation and highlight premission planning processes with an emphasis on aircraft placement, optimal instrument configuration and sensor calibrations. Critical aspects of mission operations coordinated from the NASA Johnson Spaceflight Center and integration with the JSC Flight Test Management Office will be discussed. A summary of the imagery that was obtained and processed to global surface temperature will be presented. At the capsule's point of closest approach relative to the imaging system, the spatial resolution was estimated to be approximately 15-inches per pixel and was sufficient to identify localized temperature increases associated with compression pad support hardware on the heatshield. The full paper will discuss the synergy of the quantitative imagery derived temperature maps with in-situ thermocouple measurements. Comparison of limited onboard surface thermocouple data to the image derived surface temperature will be presented. The two complimentary measurements serve as an example of the effective leveraging of resources to advance the understanding of high Mach number environments associated with an ablated heatshield and provide unique data for the validation of design tools and numerical flight simulation techniques. Collaborative opportunities and technology investments in support of planned observations of NASA's next Orion flight test in 2018 will be explored in the full manuscript.

Horvath, Thomas J.↗

Error analysis of optical correlators

With the growing interest in using binary phase only filters (BPOF) in optical correlators that are implemented on magnetooptic spatial light modulators, an understanding of the effect of errors in system alignment and optical components is critical in obtaining optimal system performance. We present simulations of optical correlator performance degradation in the presence of eight errors. We break these eight errors into three groups: 1) alignment errors, 2) errors due to a combination of component imperfections and alignment errors, and 3) errors which result solely from non-ideal components. Under the first group, we simulate errors in the distance from the object to the first principle plane of the transform lens, the distance from the second principle plane of the transform lens to the filter plane, and rotational misalignment of the input mask with the filter mask. Next we consider errors which result from a combination of alignment and component imperfections. These include errors in the transform lens, the phase compensation lens, and the inverse Fourier transform lens. Lastly we have the component errors resulting from the choice of spatial light modulator. These include contrast error and phase errors caused by the non-uniform flatness of the masks. The effects of each individual error are discussed, and the result of combining all eight errors under assumptions of reasonable tolerances and system parameters is also presented. Conclusions are drawn as to which tolerances are most critical for optimal system performance.

Ma, Paul W.↗

Strategy Developed for Selecting Optimal Sensors for Monitoring Engine Health

Sensor indications during rocket engine operation are the primary means of assessing engine performance and health. Effective selection and location of sensors in the operating engine environment enables accurate real-time condition monitoring and rapid engine controller response to mitigate critical fault conditions. These capabilities are crucial to ensure crew safety and mission success. Effective sensor selection also facilitates postflight condition assessment, which contributes to efficient engine maintenance and reduced operating costs. Under the Next Generation Launch Technology program, the NASA Glenn Research Center, in partnership with Rocketdyne Propulsion and Power, has developed a model-based procedure for systematically selecting an optimal sensor suite for assessing rocket engine system health. This optimization process is termed the systematic sensor selection strategy. Engine health management (EHM) systems generally employ multiple diagnostic procedures including data validation, anomaly detection, fault-isolation, and information fusion. The effectiveness of each diagnostic component is affected by the quality, availability, and compatibility of sensor data. Therefore systematic sensor selection is an enabling technology for EHM. Information in three categories is required by the systematic sensor selection strategy. The first category consists of targeted engine fault information; including the description and estimated risk-reduction factor for each identified fault. Risk-reduction factors are used to define and rank the potential merit of timely fault diagnoses. The second category is composed of candidate sensor information; including type, location, and estimated variance in normal operation. The final category includes the definition of fault scenarios characteristic of each targeted engine fault. These scenarios are defined in terms of engine model hardware parameters. Values of these parameters define engine simulations that generate expected sensor values for targeted fault scenarios. Taken together, this information provides an efficient condensation of the engineering experience and engine flow physics needed for sensor selection. The systematic sensor selection strategy is composed of three primary algorithms. The core of the selection process is a genetic algorithm that iteratively improves a defined quality measure of selected sensor suites. A merit algorithm is employed to compute the quality measure for each test sensor suite presented by the selection process. The quality measure is based on the fidelity of fault detection and the level of fault source discrimination provided by the test sensor suite. An inverse engine model, whose function is to derive hardware performance parameters from sensor data, is an integral part of the merit algorithm. The final component is a statistical evaluation algorithm that characterizes the impact of interference effects, such as control-induced sensor variation and sensor noise, on the probability of fault detection and isolation for optimal and near-optimal sensor suites.

Source record↗

Optimized Response Functions for 2-dimensional Observations of Solar Oscillations

Response functions are calculated for solar oscillation observations made in Doppler velocity which is averaged over a grid of 30" x 30" pixels. In a simulation of an analysis scheme proposed for the HAO/SPO Fourier Tachometer the responses for the individual pixels are combined using Chebychev weighting functions. It is shown how a technique developed for geophysical inversion may be used to find linear combinations concentrated within a fairly narrow range of l and m values of these transformed responses. The extension to bidirectional observations is discussed.

Christensen-Dalsgaard, J.↗

Tokamak elongation – how much is too much? Part 2. Numerical results

We report the analytic theory presented in Paper I is converted into a form convenient for numerical analysis. A fast and accurate code has been written using this numerical formulation. The results are presented by first defining a reference set of physical parameters based on experimental data from high performance discharges. Scaling relations of maximum achievable elongation (κ max ) versus inverse aspect ratio (ε) are obtained numerically for various values of poloidal beta (β p ), wall radius (b/a) and feedback capability parameter (γ τ w ) in ranges near the reference values. It is also shown that each value of κ max occurs at a corresponding value of optimized triangularity (δ), whose scaling is also determined as a function of ε. The results show that the theoretical predictions of κ max are slightly higher than experimental observations for high performance discharges, as measured by high average pressure. The theoretical δ values are noticeably lower. We suggest that the explanation is associated with the observation that high performance involves not only n = 0 MHD stability, but also n ≥ 1 MHD modes described by β N in the Troyon limit and transport as characterized by τ E . Operation away from the n = 0 MHD optimum may still lead to higher performance if there are more than compensatory gains in β N and τ E . Unfortunately, while the empirical scaling of β N and τ E with the elongation (κ) has been determined, the dependence on δ has still not been quantified. This information is needed in order to perform more accurate overall optimizations in future experimental designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗