Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “inverse prediction”

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 379 records · Page 21

Shear reflectivity compensation in full-waveform inversion using least-squares reverse-time migration

SUMMARY The computational cost of elastic-waveform inversion is too high for inverting PP reflections, while using acoustic full-waveform inversion (FWI) is inaccurate because it does not depend on the shear modulus/velocity/impedance that affects elastic PP wavefield amplitudes. To solve this problem, we develop a waveform inversion method that uses acoustic least-squares reverse-time migration (LSRTM) to compensate the shear reflectivity for acoustic FWI. Our method is based on the quasi-elastic-wave equation developed by Chapman et al. (2014). The quasi-elastic-wave equation uses a linearized acoustic-wave equation with shear modulus μ as a virtual source to correct the acoustic PP wavefield amplitudes toward elastic ones. Our waveform inversion method inverts for elastic parameters by minimizing the L2 norm of the difference between recorded and predicted PP reflections modelled using the quasi-elastic-wave equation. Numerical tests on synthetic and field data show that our method can properly handle the amplitudes of elastic PP reflections and provides an accurate estimate of the P- and S-wave velocities/impedances and, in some cases, the density. The method does not need the computationally expensive numerical solution to the elastic-wave equation. It also gives a better estimate of elastic parameters than a pure LSRTM method for elastic PP reflections.

Feng, Zongcai↗

Terahertz emission spectroscopy of ultrafast exciton shift current in the noncentrosymmetric semiconductor CdS

The charge-neutral exciton has been predicted to carry genuine photocurrent due to the geometric Berry phase of the electronic bands, if the inversion symmetry in a crystal is broken. Here we detect such exciton shift current in a prototypical polar semiconductor CdS by using terahertz emission spectroscopy. A distinct peak emerges in the photocurrent spectra at the energy of the exciton resonance, which is demonstrated to result from the distinct displacements of electrons and holes in real space within the excitons to produce a finite transient charge current at subpicosecond time scale. Our findings elucidate the Berry phase physics of the charge-neutral photoexcitations and also shed light on the novel energy harvesting mechanism by exciton generation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

High-Efficiency Solar-To-Fuel Photoelectrochemistry in Disordered Photonic Glass Electrodes (Final Technical Report)

This project investigated how photonic glass (PG) photoelectrodes—disordered arrangements of dielectric scatterers—can serve as scalable, tunable platforms for light trapping in photoelectrochemical (PEC) solar-to-fuel systems. By leveraging disorder-driven optical phenomena such as multiple scattering resonances and light localization, PG structures offer an alternative to conventional photonic crystals and inverse opals that require high structural precision. The scientific goals were to twofold: (1) develop approaches to predictive models for high performance PG electrodes based on light absorption simulations, and (2) fabricate, characterize, and optimize PG-based photoelectrodes for solar-to-hydrogen and solar-to-fuel photoelectrochemical applications. To overcome the complexity of ensemble optical simulations for disordered materials, the researchers developed a machine-learning-accelerated emulation of all configurations in the design space. With this approach, PG photoelectrodes based on a TiO2 semiconductor were designed to enhance PEC currents of up to one hundred times higher than the equivalent ultra-thin film photoanodes and several times higher than the equivalent photonic crystal. The research also explored integrated systems for electrochemical hydrogen production based on replacing water oxidation with the specific glycerol oxidation electrocatalysis. Overall, the project outlined an approach to a simple-to-fabricate photoelectrode system to drive photoelectrochemical reactions relevant to solar photochemical energy conversion.

14 SOLAR ENERGY↗

Range of earth structure nonuniqueness implied by body wave observations.

The Herglotz-Wiechert integral for the direct inversion of ray parameter versus distance curves can be manipulated to find the envelope of all possible models consistent with geometrical body wave observations (travel time and ray parameter versus distance). Such an extremal inversion approach has been used to find the uncertainty bounds for the velocity structure in the mantle and core. It is found, for example, that there is an uncertainty of plus or minus 40 km in the radius of the inner core boundary, plus or minus 18 km at the core-mantle boundary, and plus or minus 35 km at the 435-km transition zone. The velocity uncertainty is about plus or minus 0.08 km/sec for P and S waves in the lower mantle and about plus or minus 0.20 km/sec in the core. Experiments with various combinations of ray types in the core indicate that rather crude observations of SKKS-SKS travel times confine the range of possible models far more dramatically than do the most precise estimates of PmKP travel times. Comparisons of results from extremal inversion and linearized perturbation inversions indicate that body wave behavior is too strongly nonlinear for linearized schemes to be effective for predicting uncertainty.

Wiggins, R. A.↗

Electron transport and the magnetic field in the Crab Nebula

Observational data obtained at radio, optical, and X-ray frequencies on the size of the emission region in the Crab Nebula are used to examine various energy-transport models. It is shown that the data cannot be fitted by diffusion models, even when the energy-injection region is extended. A more favorable comparison is made with an extended source model in which bulk motion is the transport mechanism. Assuming a point source of electrons and a uniform nebular magnetic field, a perpendicular magnetic-field component of less than 0.0008 gauss is determined. Combining this with the lower limit obtained from gamma-ray observations, a perpendicular field component of between 0.0005 and 0.0008 gauss is obtained. A break in the size dependence is predicted at about 12 microns and 45 keV. Above 45 keV, the size should decrease as the inverse square root of proton energy.

Wolff, R. S.↗

Photolysis of methane and the ionosphere of Uranus

Photochemical calculations for Uranus predict an extensive region of condensation of acetylene, ethane and methane in the vicinity of the temperature inversion layer. This could explain why ethane was not detected on Uranus, unlike Neptune which has a much warmer inversion layer. Subsequent snow-out of the condensibles is expected to result in reduced visibility in the troposphere. Ionospheric calculations for the equatorial region to be probed by Voyager, indicate peak electron concentrations on the order of 5,000 per cu cm, if dynamical effects are important. Upper limit to the electron peak is 30,000 per cu cm. Exospheric temperatures as high as 200-250K are conceivable.

Atreya, S. K.↗

Near-infrared observations of IRAS sources in and near dense cores

We report 0.4 to 20 micron photometry of 27 Infrared Astronomy Satellite (IRAS) point sources associated with dense cores in nearby dark clouds. It is found that these objects have a bimodal distribution of spectral slopes. The objects in the group with steep spectral slopes are typically within one half power radius of the core peak and are not visible on the POSS. The remainder of the objects are typically further from their cores, are optically visible and have shallow spectral slopes (s is approximately 0.6). Most of the sources in this group which have been previously identified are T Tauri stars. Both the groups of objects have essentially the same luminosity function with median luminosity of 1 to 2 L(solar), which is similiar to the luminosity function for T Tauri stars. The near-infrared (J,H,K,L) colors of the objects have been used to estimate the visual extinction to the stars. This indicates that the typical steep spectrum source has extinction A(upsilon) of approximately 30 magnitudes, which is larger by a factor of 3 to 4 than a uniform core can provide, thus the density must rise steeply in the vicinity of the star. For those objects where optical estimates of the extinction are available the optical estimates are typically smaller than our IR estimates by a factor of 1.6. For multiple scattering by grains that scatter primarily in the forward direction, this implies a grain albedo of 0.4 at V, consistent with theoretical predictions for bare grains. Using this derived extinction, the FIR emission from a spherical shell with an inverse power law density which contributes the correct amount of extinction has been modelled. It is shown that to be consistent with both the observations of the density of the core at a scale of 1.5 x 10(17) cm and the extinction derived from the near-IR colors there must be a circumstellar hole of radius 10-100 A.U. for all density laws with exponents greater than 1.2. It is also shown that a component in addition to the star and circumstellar shell is necessary to match the observed spectrum of steep spectrum sources.

Fuller, G. A.↗

Rotational and Translational Components of Motion Parallax: Observers' Sensitivity and Implications for Three-Dimensional Computer Graphics

The motion of objects during motion parallax can be decomposed into 2 observer-relative components: translation and rotation. The depth ratio of objects in the visual field is specified by the inverse ratio of their angular displacement (from translation) or equivalently by the inverse ratio of their rotations. Despite the equal mathematical status of these 2 information sources, it was predicted that observers would be far more sensitive to the translational than rotational component. Such a differential sensitivity is implicitly assumed by the computer graphics technique billboarding, in which 3-dimensional (3-D) objects are drawn as planar forms (i.e., billboards) maintained normal to the line of sight. In 3 experiments, observers were found to be consistently less sensitive to rotational anomalies. The implications of these findings for kinetic depth effect displays and billboarding techniques are discussed.

Kaiser, Mary K.↗

A Hysteresis Model for Piezoceramic Materials

This paper addresses the modeling of nonlinear constitutive relations and hysteresis inherent to piezoceramic materials at moderate to high drive levels. Such models are, necessary to realize the, full potential of the materials in high performance control applications, and a necessary prerequisite is the development of techniques which permit control implementation. The approach employed here is based on the qualification of reversible and irreversible domain wall motion in response to applied electric fields. A comparison with experimental data illustrates that because the resulting ODE model is physics-based, it can be employed for both characterization and prediction of polarization levels throughout the range of actuator operation. Finally, the ODE formulation is amenable to inversion which facilitates the development of an inverse compensator for linear control design.

Smith, Ralph C.↗

Development of Solid State Laser Materials for Application in Lasers for Atmospheric Ozone and Water Vapor Sensing

We have grown neodymium doped mixed apatite crystals, (Sr(x)Ba(l-x)5(PO4)3F, Sr5(P(1-x)V(x)O4)3F, and Ba5(P(1-x)V(x)O4)3F, and spectroscopically studied them as potential gain media for a laser source for atmospheric water sensing operating at 944.11 nm0. We conclude that an appropriate apatite host material for a 944.11 nm laser should be a mixture of Sr5(PO4)3F with a small fraction of Ba5(PO4)3F. The precise wavelength tuning around 944.11 nm can be accomplished by varying the host composition, temperature, and threshold population inversion. In apatite crystals of mixed composition, the Amplified Spontaneous Emission (ASE) loss at 1.06 microns is predicted to be significantly smaller than that in the end members.

Noginov, Makhail A.↗

Gamma-Ray Observations of the Supernova Remnant RX J0852.0-4622 with the Fermi Large Area Telescope

We report on gamma-ray observations of the supernova remnant (SNR) RX J0852.04622 with the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope. In the Fermi-LAT data, we find a spatially extended source at the location of the SNR. The extension is consistent with the SNR size seen in other wavelengths such as X-rays and TeV gamma rays, leading to the identification of the gamma-ray source with the SNR. The spectrum is well described as a power law with a photon index of = 1.85 0.06 (stat)+0.18 0.19 (sys), which smoothly connects to the H.E.S.S. spectrum in the TeV energy band. We discuss the gamma-ray emission mechanism based on multiwavelength data. The broadband data can be fit well by a model in which the gamma rays are of hadronic origin. We also consider a scenario with inverse Compton scattering of electrons as the emission mechanism of the gamma rays. Although the leptonic model predicts a harder spectrum in the Fermi-LAT energy range, the model can fit the data considering the statistical and systematic errors.

Tanaka, T.↗

Comparison of the Tangent Linear Properties of Tracer Transport Schemes Applied to Geophysical Problems.

A number of geophysical applications require the use of the linearized version of the full model. One such example is in numerical weather prediction, where the tangent linear and adjoint versions of the atmospheric model are required for the 4DVAR inverse problem. The part of the model that represents the resolved scale processes of the atmosphere is known as the dynamical core. Advection, or transport, is performed by the dynamical core. It is a central process in many geophysical applications and is a process that often has a quasi-linear underlying behavior. However, over the decades since the advent of numerical modelling, significant effort has gone into developing many flavors of high-order, shape preserving, nonoscillatory, positive definite advection schemes. These schemes are excellent in terms of transporting the quantities of interest in the dynamical core, but they introduce nonlinearity through the use of nonlinear limiters. The linearity of the transport schemes used in Goddard Earth Observing System version 5 (GEOS-5), as well as a number of other schemes, is analyzed using a simple 1D setup. The linearized version of GEOS-5 is then tested using a linear third order scheme in the tangent linear version.

4DVAR↗

Probabilistic Fatigue Damage Prognosis Using a Surrogate Model Trained Via 3D Finite Element Analysis

Utilizing inverse uncertainty quantification techniques, structural health monitoring can be integrated with damage progression models to form probabilistic predictions of a structure's remaining useful life. However, damage evolution in realistic structures is physically complex. Accurately representing this behavior requires high-fidelity models which are typically computationally prohibitive. In the present work, a high-fidelity finite element model is represented by a surrogate model, reducing computation times. The new approach is used with damage diagnosis data to form a probabilistic prediction of remaining useful life for a test specimen under mixed-mode conditions.

Leser, Patrick E.↗

Transfer Function Models and Sensitivity Analysis

In many situations, real or induced flaws such as tight cracks with known morphology cannot be manufactured in part configuration specimens or in real parts. Typically, fatigue cracks are manufactured in simple geometry specimens such as flat plates, dog-bone shaped flat or round specimens. If a nondestructive evaluation (NDE) technique is required to provide a reliably detectable target flaw size denoted as α_(90/95) for induced flaws in real part, then the direct method for qualifying the NDE procedure is to use the appropriate induced flaw specimens and perform probability of detection analysis using these flaws. This can be described as direct POD demonstration testing, which may follow guidelines of MIL-HDBK-1823. This paper considers a case, where induced flaws are not available in part configuration specimens. Therefore, a direct POD demonstration study cannot be undertaken. In such situation, general practice for NDE procedure qualification is to use artificial flaws in simple geometry and part configuration specimens, and induced flaws in the chosen simple geometry specimens. Signal response data is taken on all sets of artificial and induced flaws. NDE procedure testing on induced flaw in simple geometry specimen is called NDE demonstration testing here. A transfer function NDE procedure qualification method for forward case calculates predicted induced flaw size for demonstration using a chosen target flaw size. Another transfer function method for inverse case, calculates the target flaw size using a given demonstration flaw size. The transfer function analysis assumes relationship of artificial flaw signal responses in real parts and simple geometry specimens; and induced flaw responses in simple geometry specimens to induced flaws in real parts. The signal response transfer relationships model should be defined before transfer function models can be devised. Assuming that the signal response transfer relationships model is valid, forward and inverse case transfer function methods have been devised. Because of lack of signal response data from induced flaws in real part, 90/95% POD/confidence (P/C) cannot be demonstrated directly. However, the transfer function method may be assessed using simulation to evaluate whether the resulting target flaw size or demonstration flaw size provides adequate confidence to the assumed signal response transfer relationships model. Therefore, the transfer function approach is a risk assessment approach. Both the signal response transfer relationships model and the transfer function model are important in managing risk in results provided by the transfer function NDE technique qualification or assessment. The signal response transfer relationships model needs to be validated with empirical data and then transfer function model needs to be validated for desired P/C on case-by-case basis.

Nondestructive evaluation↗

Magnetic spectral signatures in the Earth's magnetosheath and plasma depletion layer

Correlations between plasma properties and magnetic fluctuations in the sub-solar magnetosheath downstream of a quasi-perpendicular shock have been found and indicate that mirror and ion cyclotronlike fluctuations correlate with the magnetosheath proper and plasma depletion layer, respectively (Anderson and Fueselier, 1993). We explore the entire range of magnetic spectral signatures observed from the Active Magnetospheric Particle Tracer Explorers/Charge Composition Explorer (AMPTE/CCE)spacecraft in the magnetosheath downstream of a quasi-perpendicular shock. The magnetic spectral signatures typically progress from predominantly compressional fluctuations,delta B(sub parallel)/delta B perpendicular to approximately 3, with F/F (sub p) less than 0.2 (F and F (sub p) are the wave frequency and proton gyrofrequency, respectively) to predominantly transverse fluctuations, delta B(sub parallel)/delta B perpendicular to approximately 0.3, extending up to F(sub p). The compressional fluctuations are characterized by anticorrelation between the field magnitude and electron density, n(sub e), and by a small compressibility, C(sub e) identically equal to (delta n(sub e)/n(sub e)) (exp 2) (B/delta B(sub parallel)) (exp 2) approximately 0.13, indicative of mirror waves. The spectral characteristics of the transverse fluctuations are in agreement with predictions of linear Vlasov theory for the H(+) and He(2+) cyclotron modes. The power spectra and local plasma parameters are found to vary in concert: mirror waves occur for beta(s ub parallel p) (beta (sub parallel p) identically = 2 mu(sub zero) n(sub p) kT (sub parallel p) / B(exp 2) approximately = 2, A(sub p) indentically = T(sub perpendicular to p)/T(sub parallel p) - 1 approximately = 0.4, whereas cyclotron waves occur for beta (sub parallel p) approximately = 0.2 and A(sub p) approximately = 2. The transition from mirror to cyclotron modes is predicted by linear theory. The spectral characteristics overlap for intermediate plasma parameters. The plasma observations are described by A(sub p) = 0.85 beta(sub parallel P) (exp - 0.48) with a log regression coefficient of -0.74. This inverse A(sub p) - beta(sub parallel p) correlation corresponds closely to the isocontours of maximum ion anisotropy instability growth, gamma (sub m)/omega(sub p) = 0.01, for the mirror and cyclotron modes. The agreement of observed properties and predictions of local theory suggests that the spectral signatures reflect the local plasma environment and that the anisotropy instabilities regulate A(sub p). We suggest that the spectral characteristics may provide a useful basis for ordering observations in the magnetosheath and that the A(sub p) - beta(sub parallel p) inverse correlation may be used as a beta-dependent upper limit on the proton anisotropy to represent kinetic effects.

Anderson, Brian J.↗

Updates to the Regional Seismic Travel Time (RSTT) Model: 2. Path-dependent Travel-time Uncertainty

Abstract The regional seismic travel time (RSTT) model and software were developed to improve travel-time prediction accuracy by accounting for three-dimensional crust and upper mantle structure. Travel-time uncertainty estimates are used in the process of associating seismic phases to events and to accurately calculate location uncertainty bounds (i.e. event location error ellipses). We improve on the current distance-dependent uncertainty parameterization for RSTT using a random effects model to estimate slowness (inverse velocity) uncertainty as a mean squared error for each model parameter. The random effects model separates the error between observed slowness and model predicted slowness into bias and random components. The path-specific travel-time uncertainty is calculated by integrating these mean squared errors along a seismic-phase ray path. We demonstrate that event location error ellipses computed for a 90% coverage ellipse metric (used by the Comprehensive Nuclear-Test-Ban Treaty Organization International Data Centre (IDC)), and using the path-specific travel-time uncertainty approach, are more representative (median 82.5% ellipse percentage) of true location error than error ellipses computed using distance-dependent travel-time uncertainties (median 70.1%). We also demonstrate measurable improvement in location uncertainties using the RSTT method compared to the current station correction approach used at the IDC (median 74.3% coverage ellipse).

58 GEOSCIENCES↗

Stimulated Compton conversion of Langmuir waves by relativistic electron beams

The scattering of a Langmuir wave by a relativistic electron to produce a transversely polarized high-frequency EM wave by an inverse Compton process is investigated theoretically. The results of numerical computations are presented in diagrams and graphs, and it is shown that although net growth rates are predicted for several classes of Langmuir spectra, they are at least several orders of magnitude too small to account for the intense high-frequency emission seen in recent beam-plasma experiments (Benford et al., 1980; Kato et al., 1983) unless some unknown coherence mechanism is focusing the inverse-Compton emission. The implications of the mechanism for astrophysical beam-plasma systems such as quasars and radio galaxies are indicated.

Newman, D. L.↗

Multi-fidelity Bayesian neural networks: Algorithms and applications

Here we propose a new class of Bayesian neural networks (BNNs) that can be trained using noisy data of variable fidelity, and we apply them to learn function approximations as well as to solve inverse problems based on partial differential equations (PDEs). These multi-fidelity BNNs consist of three neural networks: The first is a fully connected neural network, which is trained following the maximum a posteriori probability (MAP) method to fit the low-fidelity data; the second is a Bayesian neural network employed to capture the cross-correlation with uncertainty quantification between the low- and high-fidelity data; and the last one is the physics-informed neural network, which encodes the physical laws described by PDEs. For the training of the last two neural networks, we first employ the mean-field variational inference (VI) to maximize the evidence lower bound (ELBO) to obtain informative prior distributions for the hyperparameters in the BNNs, and subsequently we use the Hamiltonian Monte Carlo (HMC) method to estimate accurately the posterior distributions for the corresponding hyperparameters. We demonstrate the accuracy of the present method using synthetic data as well as real measurements. Specifically, we first approximate a one- and four-dimensional function, and then infer the reaction rates in one- and two-dimensional diffusion-reaction systems. Moreover, we infer the sea surface temperature (SST) in the Massachusetts and Cape Cod Bays using satellite images and in-situ measurements. Taken together, our results demonstrate that the present method can capture both linear and nonlinear correlation between the low- and high-fidelity data adaptively, identify unknown parameters in PDEs, and quantify uncertainties in predictions, given a few scattered noisy high-fidelity data. Finally, we demonstrate that we can effectively and efficiently reduce the uncertainties and hence enhance the prediction accuracy with an active learning approach, using as examples a specific one-dimensional function approximation and an inverse PDE problem.

97 MATHEMATICS AND COMPUTING↗