Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “inversion methods”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Interpretable inverse-designed cavity for on-chip nonlinear photon pair generation

Inverse design is a powerful tool in wave physics for compact, high-performance devices. To date, applications in photonics have mostly been limited to linear systems and it has rarely been investigated or demonstrated in the nonlinear regime. In addition, the “black box” nature of inverse design techniques has hindered the understanding of optimized inverse-designed structures. We propose an inverse design method with interpretable results to enhance the efficiency of on-chip photon generation rate through nonlinear processes by controlling the effective phase-matching conditions. We fabricate and characterize a compact, inverse-designed device using a silicon-on-insulator platform that allows a spontaneous four-wave mixing process to generate photon pairs at a rate of 1.1 MHz with a coincidence to accidental ratio of 162. Our design method accounts for fabrication constraints and can be used for scalable quantum light sources in large-scale communication and computing applications.

42 ENGINEERING↗
NASA NTRS1982

Backscattered UV radiation - Effects of multiple scattering and the lower boundary of the atmosphere

A method is proposed for the calculation of a multiple-scattering correction to the single-scattering calculation of the radiance of the terrestrial atmosphere resulting from backscattered ultraviolet solar radiation in the spectral region used in the ozone profile inversion. This method uses jointly the usual analytical and Monte Carlo methods. Effects of the lower boundary of the atmosphere, cloud tops, and ground surface are investigated both qualitatively and quantitatively. The ratio of multiple to single scattering is determined, and its importance in ozone profile inversion of backscattered UV solar radiation from the terrestrial atmosphere is evaluated. The polarization of the atmospheric radiance is treated briefly.

Aruga, T.↗
NASA NTRS1976

Computer assisted expansion of 1/Delta series with Newton iteration

The Newton iteration method has been widely applied to the solution of various equations such as Kepler's equation. In this study it is used in planetary and satellite theory as a general procedure for Fourier series inversion. The method is used for the construction of the 1/Delta series either in literal form or in numerical form with small eccentricities and inclinations substituted in advance. This usually results in very compact series. With the Newton iteration procedure and a computerized series manipulation technique, the Fourier series of 1/Delta of the mutual perturbations among most natural satellites can be easily constructed.

Chao, C. C.↗
NASA NTRS1983

Computation of viscous flows over airfoils, including separation, with a coupling approach

Viscous incompressible flows over single or multiple airfoils, with or without separation, were computed using an inviscid flow calculation, with modified boundary conditions, and by a method providing calculation and coupling for boundary layers and wakes, within conditions of strong viscous interaction. The inviscid flow is calculated with a method of singularities, the numerics of which were improved by using both source and vortex distributions over profiles, associated with regularity conditions for the fictitious flows inside of the airfoils. The viscous calculation estimates the difference between viscous flow and inviscid interacting flow, with a direct or inverse integral method, laminar or turbulent, with or without reverse flow. The numerical method for coupling determines iteratively the boundary conditions for the inviscid flow. For attached viscous layers regions, an underrelaxation is locally calculated to insure stability. For separated or separating regions, a special semi-inverse algorithm is used. Comparisons with experiments are presented.

Leballeur, J. C.↗
NASA NTRS1987

Inverse kinematic-based robot control

A fundamental problem which must be resolved in virtually all non-trivial robotic operations is the well-known inverse kinematic question. More specifically, most of the tasks which robots are called upon to perform are specified in Cartesian (x,y,z) space, such as simple tracking along one or more straight line paths or following a specified surfacer with compliant force sensors and/or visual feedback. In all cases, control is actually implemented through coordinated motion of the various links which comprise the manipulator; i.e., in link space. As a consequence, the control computer of every sophisticated anthropomorphic robot must contain provisions for solving the inverse kinematic problem which, in the case of simple, non-redundant position control, involves the determination of the first three link angles, theta sub 1, theta sub 2, and theta sub 3, which produce a desired wrist origin position P sub xw, P sub yw, and P sub zw at the end of link 3 relative to some fixed base frame. Researchers outline a new inverse kinematic solution and demonstrate its potential via some recent computer simulations. They also compare it to current inverse kinematic methods and outline some of the remaining problems which will be addressed in order to render it fully operational. Also discussed are a number of practical consequences of this technique beyond its obvious use in solving the inverse kinematic question.

Wolovich, W. A.↗
DOE OSTI2021

Solving Inverse Stochastic Problems from Discrete Particle Observations Using the Fokker--Planck Equation and Physics-Informed Neural Networks

The Fokker--Planck (FP) equation governing the evolution of the probability density function (PDF) is applicable to many disciplines, but it requires specification of the coefficients for each case, which can be functions of space-time and not just constants and hence require the development of a data-driven modeling approach. When the data available is directly on the PDF, there exist methods for inverse problems that can be employed to infer the coefficients and thus determine the FP equation and subsequently obtain its solution. Herein, we address a more realistic scenario, where only sparse data are given on the particles' positions at a few time instants, which are not sufficient to accurately construct directly the PDF even at those times from existing methods, e.g., kernel estimation algorithms. To this end, we develop a general framework based on physics-informed neural networks (PINNs) that introduces a new loss function using the Kullback--Leibler divergence to connect the stochastic samples with the FP equation to simultaneously learn the equation and infer the multidimensional PDF at all times. In particular, we consider two types of inverse problems, type I, where the FP equation is known but the initial PDF is unknown, and type II, in which, in addition to the unknown initial PDF, the drift and diffusion terms are also unknown. In both cases, we investigate problems with either Brownian or Lévy noise or a combination of both. Here, we demonstrate the new PINN framework in detail in the one-dimensional (1D) case, but we also provide results for up to five dimensions demonstrating that we can infer both the FP equation and dynamics simultaneously at all times with high accuracy using only very few discrete observations of the particles.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2010

Method and apparatus for the guided ablative therapy of fast ventricular arrhythmia

Method and apparatus for guiding ablative therapy of abnormal biological electrical excitation. The excitation from the previous excitatory wave is significant at the beginning of the next excitation. In particular, it is designed for treatment of fast cardiac arrhythmias. Electrical signals are acquired from recording electrodes, and an inverse dipole method is used to identify the site of origin of an arrhythmia. The location of the tip of an ablation catheter is similarly localized from signals acquired from the recording electrodes while electrical pacing energy is delivered to the tip of the catheter close to or in contact with the cardiac tissue. The catheter tip is then guided to the site of origin of the arrhythmia, and ablative radio frequency energy is delivered to its tip to ablate the site.

Cohen, Richard J.↗
NASA NTRS1999

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.↗
NASA NTRS1989

Finite element analysis in fluids; Proceedings of the Seventh International Conference on Finite Element Methods in Flow Problems, University of Alabama, Huntsville, Apr. 3-7, 1989

Recent advances in computational fluid dynamics are examined in reviews and reports, with an emphasis on finite-element methods. Sections are devoted to adaptive meshes, atmospheric dynamics, combustion, compressible flows, control-volume finite elements, crystal growth, domain decomposition, EM-field problems, FDM/FEM, and fluid-structure interactions. Consideration is given to free-boundary problems with heat transfer, free surface flow, geophysical flow problems, heat and mass transfer, high-speed flow, incompressible flow, inverse design methods, MHD problems, the mathematics of finite elements, and mesh generation. Also discussed are mixed finite elements, multigrid methods, non-Newtonian fluids, numerical dissipation, parallel vector processing, reservoir simulation, seepage, shallow-water problems, spectral methods, supercomputer architectures, three-dimensional problems, and turbulent flows.

Chung, T. J.↗
DOE OSTI2023

Estimating source-sink distributions and fluxes of reactive nitrogen and sulfur within a mixed forest canopy

The vertical source-sink distribution of air pollutants within and above forested canopies is necessary for describing the biological, physical, and chemical processes influencing the soil-vegetation-atmosphere exchange. Here, this study implemented inverse modeling methods to estimate the source-sink and flux profiles of reactive nitrogen (N) and sulfur (S) compounds from measurements of the mean concentration profiles of ammonia (NH 3 ), nitric acid (HNO 3 ), sulfur dioxide (SO 2 ), and particulate ammonium (NH 4 + ), nitrate (NO 3 − ), and sulfate (SO 4 2− ) at a forest site in the southern Appalachian Mountains. Three inverse approaches utilizing different approximations to scalar transport within the canopy were developed and evaluated against sensible heat flux measurements. The Eulerian model (EUL), which incorporates vertical velocity skewness, performed well in reproducing the turbulent heat fluxes and was subsequently used to calculate the chemical source-sink and flux profiles. Above-canopy fluxes of NH 3 were downward, indicating that the forest was a net sink of NH 3 . The soil/litter layer was both a source and a sink for NH 3 but the exchange rate at the forest floor was small. Fluxes of HNO 3 , SO 2 , NO 3 − , NH 4 + , and SO 4 2- were uni-directional (deposition only) between the air and the canopy/ground and increased monotonically from the forest floor to the canopy top. Crown foliage dominated the uptake of reactive N and S during the growing season, accounting for 80–90% of the total canopy-scale flux. Fluxes and canopy-ground partitioning estimated using the resistance-based Surface Tiled Aerosol and Gas Exchange (STAGE) model were generally comparable to EUL. The comparison highlights the need for improved parameterizations of litter exchange and NH 3 compensation points in resistance models for forest ecosystems. The findings here benefit the application of critical loads in forest ecosystems and guide further development of resistance-based exchange models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Prospects for Surrogate Neutron Capture Measurements with Radioactive Ion Beams and GODDESS

Neutron capture reactions are responsible for the synthesis of almost all of the elements heavier than iron through the slow s-process, that proceeds close to the line of stability, and the rapid r-process, with very neutron-rich waiting points. Uncertainties in (n,γ) rates in neutron rich nuclei, especially near closed neu- tron shells, can have significant impact [1] on the predictions of final abundances for different astrophysical scenarios for the r process. Understanding (n,γ) rates on neutron-rich fission fragments is also important for nuclear forensics and stockpile stewardship science. Ratkiewicz et al. [2 and references therein] has recently demonstrated that the (d,pγ) reaction is a valid surrogate for (n,γ), where the formation of the compound nucleus from the breakup of the deuteron has been calculated in a reaction model and the subsequent measured gamma-decay probabilities are reproduced with standard level density and strength functions in a Bayesian approach. In parallel to the surrogate validation efforts, we have demonstrated that the (d,pγ) reaction can be measured in inverse kinematics with Gammasphere ORRUBA: Dual Detectors for Experimental Structure Studies (GODDESS) [3] where the Gammasphere array of Compton-suppressed HPGe detectors is coupled to the Oak Ridge Rutgers University Barrel Array of position-sensitive silicon strip detectors. During the commissioning campaign we measured the (d,pγ) reaction with 134 Xe and 95 Mo beams, the latter to demonstrate the surrogate method in inverse kinematics. The present talk will present preliminary results from this campaign including γ-decay probabilities and prospects for surrogate (n,γ) measurements with 143 Ba fission-fragment beams.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Statistical Treatment of Inverse Problems Constrained by Differential Equations-based Models with Stochastic Terms

This paper introduces a statistical treatment of inverse problems constrained by models with stochastic terms. The solution of the forward problem is given by a distribution represented numerically by an ensemble of simulations. The goal is to formulate the inverse problem, in particular the objective function, to find the closest forward distribution (i.e., the output of the stochastic forward problem) that best explains the distribution of the observations in a certain metric. We use proper scoring rules, a concept employed in statistical forecast verification, namely energy, variogram, and hybrid (i.e., combination of the two) scores. We study the performance of the proposed formulation in the context of two applications: a coefficient field inversion for subsurface flow governed by an elliptic partial differential equation with a stochastic source and a parameter inversion for power grid governed by differential-algebraic equations. In both cases we show that the variogram and the hybrid scores produce better parameter inversion results than does the energy score, whereas the energy score leads to better probabilistic predictions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

sMF-BO-2CoGP: A sequential multi-fidelity constrained Bayesian optimization framework for design applications

Bayesian optimization (BO) is an effective surrogate-based method that has been widely used to optimize simulation-based applications. While the traditional Bayesian optimization approach only applies to single-fidelity models, many realistic applications provide multiple levels of fidelity with various levels of computational complexity and predictive capability. In this work, we propose a multi-fidelity Bayesian optimization method for design applications with both known and unknown constraints. The proposed framework, called sMF-BO-2CoGP, is built on a multi-level CoKriging method to predict the objective function. An external binary classifier, which we approximate using a separate CoKriging model, is used to distinguish between feasible and infeasible regions. Finally, the sMF-BO-2CoGP method is demonstrated using a series of analytical examples and a flip-chip application for design optimization to minimize the deformation due to warping under thermal loading conditions.

42 ENGINEERING↗
NASA NTRS1985

The Nagoya cosmic-ray muon spectrometer 3, part 4: Track reconstruction method

One of the greatest problems in measuring particle trajectories with an optical or visual detector system, is the reconstruction of trajectories in real space from their recorded images. In the Nagoya cosmic-ray muon spectrometer, muon tracks are detected by wide gap spark chambers and their images are recorded on the photographic film through an optical system of 10 mirrors and two cameras. For the spatial reconstruction, 42 parameters of the optical system should be known to determine the configuration of this system. It is almost impossible to measure this many parameters directly with usual techniques. In order to solve this problem, the inverse transformation method was applied. In this method, all the optical parameters are determined from the locations of fiducial marks in real space and the locations of their images on the photographic film by the non-linear least square fitting.

Shibata, S.↗
DOE OSTI2023

More Than 30 000-fold Field Enhancement of Terahertz Nanoresonators Enabled by Rapid Inverse Design

The rapid development of 6G communications using terahertz (THz) electromagnetic waves has created a demand for highly sensitive THz nanoresonators capable of detecting these waves. Among the potential candidates, THz nanogap loop arrays show promising characteristics but require significant computational resources for accurate simulation. This requirement arises because their unit cells are 10 times smaller than millimeter wavelengths, with nanogap regions that are 1 000 000 times smaller. To address this challenge, we propose a rapid inverse design method using physics-informed machine learning, employing double deep Q-learning with an analytical model of the THz nanogap loop array. In ~39 h on a middle-level personal computer, our approach identifies the optimal structure through 200 000 iterations, achieving an experimental electric field enhancement of 32 000 at 0.2 THz, 300% stronger than prior results. In conclusion, our analytical model-based approach significantly reduces the amount of computational resources required, offering a practical alternative to numerical simulation-based inverse design for THz nanodevices.

21 terahertz time-domain spectroscopy↗
DOE OSTI2022

Broadband Asymmetric Transmission of Linearly Polarized Mid‐Infrared Light Based on Quasi‐3D Metamaterials

Abstract Metamaterials consisting of subwavelength resonators offer an exciting opportunity for realizing asymmetric transmission (AT) of linearly polarized light. However, to date, only moderate/narrow‐band AT responses have been obtained in metadevices based on stacked planar nanostructures. Here, leveraging a combination of a genetic algorithm (GA) based optimization method and a membrane projection lithography (MPL) fabrication approach, a quasi‐3D metamaterial for broadband AT of linearly polarized mid‐infrared light is demonstrated. Facilitated by the customized GA, an efficient exploration of 3D plasmonic meta‐atoms with broken mirror symmetry in the light propagation direction allows the satisfaction of the rigorous conditions for AT of linearly polarized waves over a broad wavelength range. Confirmed by surface current analysis, the observed AT behavior is attributed to the resonant coupling between the plasmonic nanostructures located on the two orthogonal walls of the MPL cavities. Incorporating an advanced inverse‐design method and a state‐of‐art fabrication technique, the methodology used in the present study provides a promising route for exploiting 3D metamaterials with sophisticated functionalities via effectively exploring the high‐dimensional parametric space offered by true 3D meta‐atoms.

Whiting, Eric B.↗
DOE OSTI2020

Detector Tilt Considerations in Bragg Coherent Diffraction Imaging: A Simulation Study

This paper addresses the three-dimensional signal distortion and image reconstruction issues in X-ray Bragg coherent diffraction imaging (BCDI) in the event of a general non-orthogonal orientation of the area detector with respect to the diffracted beam. Growing interest in novel BCDI adaptations at fourth-generation synchrotron light sources has necessitated improvisations in the experimental configuration and the subsequent data analysis. One such possibly unavoidable improvisation that is envisioned in this paper is a photon-counting area detector whose face is tilted away from the perpendicular to the Bragg-diffracted beam during the acquisition of the coherent diffraction signal. We describe a likely circumstance in which one would require such a detector configuration, along with the experimental precedent at third-generation synchrotrons. Using physically accurate diffraction simulations from synthetic scatterers in the presence of such tilted detectors, we analyze the general nature of the observed signal distortion qualitatively and quantitatively and provide a prescription to correct for it during image reconstruction. Our simulations and reconstructions are based on an adaptation of the known theory of BCDI sampling geometry, as well as the recently developed projection-based methods of wavefield propagation. Such configurational modifications and their numerical remedies are potentially valuable in realizing unconventional coherent diffraction measurement geometries, eventually paving the way for the integration of BCDI into new material characterization experiments at next-generation light sources.

36 MATERIALS SCIENCE↗