Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model inversion”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,765 records · Page 98

Making Invisible Visible: Data-Driven Seismic Inversion With Spatio-Temporally Constrained Data Augmentation

Deep learning and data-driven approaches have shown great potential in scientific domains. The promise of data-driven techniques relies on the availability of a large volume of high-quality training datasets. Due to the high cost of obtaining data through expensive physical experiments, instruments, and simulations, data augmentation techniques for scientific applications have emerged as a new direction for obtaining scientific data recently. However, existing data augmentation techniques originating from computer vision yield physically unacceptable data samples that are not helpful for the domain problems that we are interested in. In this article, we develop new data augmentation techniques based on convolutional neural networks. Specifically, our generative models leverage different physics knowledge (such as governing equations, observable perception, and physics phenomena) to improve the quality of the synthetic data. To validate the effectiveness of our data augmentation techniques, we apply them to solve a subsurface seismic full-waveform inversion using simulated CO 2 leakage data. Our interest is to invert for subsurface velocity models associated with very small CO 2 leakage. We validate the performance of our methods using comprehensive numerical tests. Here via comparison and analysis, we show that data-driven seismic imaging can be significantly enhanced by using our data augmentation techniques. Particularly, the imaging quality has been improved by 15% in test scenarios of general-sized leakage and 17% in small-sized leakage when using an augmented training set obtained with our techniques.

58 GEOSCIENCES↗

The Determination of the Large-Scale Circulation of the Pacific Ocean from Satellite Altimetry using Model Green's Functions

A Green's function method for obtaining an estimate of the ocean circulation using both a general circulation model and altimetric data is demonstrated. The fundamental assumption is that the model is so accurate that the differences between the observations and the model-estimated fields obey a linear dynamics. In the present case, the calculations are demonstrated for model/data differences occurring on very a large scale, where the linearization hypothesis appears to be a good one. A semi-automatic linearization of the Bryan/Cox general circulation model is effected by calculating the model response to a series of isolated (in both space and time) geostrophically balanced vortices. These resulting impulse responses or 'Green's functions' then provide the kernels for a linear inverse problem. The method is first demonstrated with a set of 'twin experiments' and then with real data spanning the entire model domain and a year of TOPEX/POSEIDON observations. Our present focus is on the estimate of the time-mean and annual cycle of the model. Residuals of the inversion/assimilation are largest in the western tropical Pacific, and are believed to reflect primarily geoid error. Vertical resolution diminishes with depth with 1 year of data. The model mean is modified such that the subtropical gyre is weakened by about 1 cm/s and the center of the gyre shifted southward by about 10 deg. Corrections to the flow field at the annual cycle suggest that the dynamical response is weak except in the tropics, where the estimated seasonal cycle of the low-latitude current system is of the order of 2 cm/s. The underestimation of observed fluctuations can be related to the inversion on the coarse spatial grid, which does not permit full resolution of the tropical physics. The methodology is easily extended to higher resolution, to use of spatially correlated errors, and to other data types.

Stammer, Detlef↗

Radiative Transfer and Retrievals in EOF Domain

The Infrared Atmospheric Sounding Interferometer (IASI) is a hyperspectral sensor with 8461 spectral channels and a nominal spectral resolution of 0.25 cm(sup -1). It is computationally intensive to perform radiative transfer calculations and inversions using all these channels. We will present a Principal Component-based Radiative Transfer Model (PCRTM) and a retrieval algorithm which perform all the necessary calculations in EOF domain. Since the EOFs are orthogonal to each other, only about 100 principal components are needed to represent the information content of the 8461 channels. The PCRTM provides the EOF coefficients and associated derivatives with respect to atmospheric and surface parameters needed by the inversion algorithm. The inversion algorithm is based on a non-linear Levenberg-Marquardt method with climatology covariance and a priori information as constraints. The retrieved parameters include atmospheric temperature, moisture and ozone profiles, cloud parameters, surface skin temperature, and surface emissivities. To make the retrieval system even more compact and stable. The atmospheric vertical profiles are compressed into the EOF space as well. The surface emissivities are also compressed into EOF space.

Liu, Xu↗

Ab Initio-Based Bond Order Potential for Arsenene Polymorphs Developed via Hierarchical Reinforcement Learning

Arsenene, a less-explored two-dimensional material, holds the potential for applications in wearable electronics, memory devices, and quantum systems. This study introduces a bond-order potential model with Tersoff formalism, the ML-Tersoff, which leverages multireward hierarchical reinforcement learning (RL), trained on an ab initio data set. This data set covers a spectrum of properties for arsenene polymorphs, enhancing our understanding of its mechanical and thermal behaviors without the complexities of traditional models requiring multiple parameter sets. Our RL strategy utilizes decision trees coupled with a hierarchical reward strategy to accelerate convergence in high-dimensional continuous search spaces. Unlike the Stillinger-Weber approach, which demands separate formalisms for buckled and puckered forms, the ML-Tersoff model concurrently captures multiple properties of the two polymorphs by effectively representing the local environment, thereby avoiding the need for different atomic types. Here, we apply the ML model to understand the mechanical and thermal properties of the arsenene polymorphs and nanostructures. We observe an inverse relationship between the critical strain and temperature in arsenene. Thermal conductivity calculations in nanosheets show good agreement with ab initio data, reflecting a decrease in thermal conductivity attributable to increased anharmonic effects at higher temperatures. We also apply the model to predict the thermal behavior of arsenene nanotubes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Automated Processing of ISIS Topside Ionograms into Electron Density Profiles

Modeling of the topside ionosphere has for the most part relied on just a few years of data from topside sounder satellites. The widely used Bent et al. (1972) model, for example, is based on only 50,000 Alouette 1 profiles. The International Reference Ionosphere (IRI) (Bilitza, 1990, 2001) uses an analytical description of the graphs and tables provided by Bent et al. (1972). The Alouette 1, 2 and ISIS 1, 2 topside sounder satellites of the sixties and seventies were ahead of their times in terms of the sheer volume of data obtained and in terms of the computer and software requirements for data analysis. As a result, only a small percentage of the collected topside ionograms was converted into electron density profiles. Recently, a NASA-funded data restoration project has undertaken and is continuing the process of digitizing the Alouette/ISIS ionograms from the analog 7-track tapes. Our project involves the automated processing of these digital ionograms into electron density profiles. The project accomplished a set of important goals that will have a major impact on understanding and modeling of the topside ionosphere: (1) The TOPside Ionogram Scaling and True height inversion (TOPIST) software was developed for the automated scaling and inversion of topside ionograms. (2) The TOPIST software was applied to the over 300,000 ISIS-2 topside ionograms that had been digitized in the fkamework of a separate AISRP project (PI: R.F. Benson). (3) The new TOPIST-produced database of global electron density profiles for the topside ionosphere were made publicly available through NASA s National Space Science Data Center (NSSDC) ftp archive at . (4) Earlier Alouette 1,2 and ISIS 1, 2 data sets of electron density profiles from manual scaling of selected sets of ionograms were converted fiom a highly-compressed binary format into a user-friendly ASCII format and made publicly available through nssdcftp.gsfc.nasa.gov. The new database for the topside ionosphere established as a result of this project, has stimulated a multitude of new studies directed towards a better description and prediction of the topside ionosphere. Marinov et al. (2004) developed a new model for the upper ion transition height (Oxygen to Hydrogen and Helium) and Bilitza (2004) deduced a correction term for the I N topside electron density model. Kutiev et al. (2005) used this data to develop a new model for the topside ionosphere scale height (TISH) as a function of month, local time, latitude, longitude and solar flux F10.7. Comparisons by Belehaki et al. (2005) show that TISH is in general agreement with scale heights deduced from ground ionosondes but the model predicts post-midnight and afternoon maxima whereas the ionosonde data show a noon maximum. Webb and Benson (2005) reported on their effort to deduce changes in the plasma temperature and ion composition from changes in the topside electron density profile as recorded by topside sounders. Limitations and possible improvements of the IRI topside model were discussed by Coisson et al. (2005) including also the possible use of the NeQuick model, Our project progressed in close collaboration and coordination with the GSFC team involved in the ISIS digitization effort. The digitization project was highly successful producing a large amount of digital topside ionograms. Several no-cost extensions of the TOPIST project were necessary to keep up with the pace and volume of the digitization effort.

Reinisch, bodo W.↗

The Connection between Model Performance on the CCMVal Transport Diagnostics and Simulated Sensitivity of Ozone to Chlorine Change

The chemistry climate models that contributed simulations for past and future ozone evolution to the 2010 Scientific Assessment of Ozone Depletion were subject to extensive evaluation by the SPARC (Stratospheric Processes and their Role in Climate) CCMVal (Chemistry-Climate Model Validation) activity. The sensitivity of ozone to changes in composition and climate varies among the models, but the relationship between these variations and the model evaluations of CCMVal is not obvious. We have learned that the transport evaluation can be used to interpret the comparisons between observed and simulated columns of chlorine reservoirs, hydrochloric acid (HCl) and chlorine nitrate (ClONO2); these comparisons were part of the CCMVal evaluation of chemistry. The simulations with best performance on the transport diagnostics most faithfully reproduce the evolution and seasonal variation of the chlorine reservoirs as observed at NDACC (Network for Detection of Atmospheric Composition Change) stations (NyAlesund 78.9N, Kiruna 67.8N, Harestua 60.2N, Jungfraujoch 46.6N, Toronto 43.6N, Kitt Peak 31.9N, Izana 28.3N, Mauna Loa 19.5N, Lauder 45S and Arrival Heights 77.8S). In the simulations, the HCl in the lower stratosphere depends on total inorganic chlorine (Cly) and partitioning between HCl and ClON02. Total inorganic chlorine depends on the fractional release of chlorine from source gases, and ratio of ClON02 to HCl is inversely dependent on methane and varies quadratically with ozone. Simulated HCl from various models may agree with observations even though Cly is in error, partitioning is in error, or both. Simulated ozone sensitivity to chlorine is shown to be greater for models that produce larger values of chlorine nitrate for background chlorine levels, and vice versa. Comparisons with the NDACC data show why the models with 'best' transport have similar sensitivity to chlorine change. The realistic evolution of the simulated HCl and ClONO2 columns suggests realistic levels of Cly in the lower atmosphere. In addition, the wide range values for the sensitivity of ozone to chlorine obtained from the CCMVal simulations is explained by the wide range in lower atmospheric columns of ClONO2 and the concomitant wide range of levels for chlorine monoxide.

Douglass, Anne↗

Constraints on a massive double-copy and applications to massive gravity

We propose and study a BCJ double-copy of massive particles, showing that it is equivalent to a KLT formula with a kernel given by the inverse of a matrix of massive bi-adjoint scalar amplitudes. For models with a uniform non-zero mass spectrum we demonstrate that the resulting double-copy factors on physical poles and that up to at least 5-particle scattering, color-kinematics duality satisfying numerators always exist. For the scattering of 5 or more particles, the procedure generically introduces spurious singularities that must be cancelled by imposing additional constraints. When massive particles are present, color-kinematics duality is not enough to guarantee a physical double-copy. As an example, we apply the formalism to massive Yang-Mills and show that up to 4-particle scattering the double-copy construction generates physical amplitudes of a model of dRGT massive gravity coupled to a dilaton and a two-form with dilaton parity violating couplings. We show that the spurious singularities in the 5-particle double-copy do not cancel in this example, and the construction fails to generate physically sensible amplitudes. We conjecture sufficient constraints on the mass spectrum, which in addition to massive BCJ relations, guarantee the absence of spurious singularities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Lepto-axiogenesis and the scale of supersymmetry

If the Peccei-Quinn field containing the QCD axion undergoes rotations in the early universe, the dimension-five operator responsible for neutrino masses can generate a lepton asymmetry that ultimately gives rise to the observed baryon asymmetry of the Universe. This lepto-axiogenesis scenario requires a flat potential for the radial direction of the Peccei-Quinn field, naturally realized in supersymmetric models. We carefully compute the efficiency of this mechanism for the Dine-Fischler-Srednicki-Zhitnitsky (DFSZ) and Kim-Shifman-Vainshtein-Zakharov (KSVZ) axion models and place lower bounds on the masses of scalar superpartners required to reproduce the observed baryon asymmetry. For the KSVZ model, we find an efficiency for generation of the asymmetry six times larger than the previously extant computation after including scattering channels involving superpartners. In this case, the superpartner scale should be above 30 TeV for a domain wall number of one; the lower bound weakens for larger domain wall numbers. We find that the superpartner mass scale may also be as low as ~ 30 TeV for the DFSZ model. In all cases, the lower bound on the superpartner masses is inversely proportional to the sum of the squares of the neutrino masses and so can strengthen as the upper bound on the neutrino mass improves. We identify the parameter space where the axion rotation can simultaneously produce axion dark matter via kinetic misalignment; in this case it is possible to put an upper bound of order PeV on the masses of scalar superpartners.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

An asynchronous parallel high-throughput model calibration framework for crystal plasticity finite element constitutive models

Crystal plasticity finite element model (CPFEM) is a powerful numerical simulation in the integrated computational materials engineering toolboxes that relates microstructures to homogenized materials properties and establishes the structure–property linkages in computational materials science. However, to establish the predictive capability, one needs to calibrate the underlying constitutive model, verify the solution and validate the model prediction against experimental data. Bayesian optimization (BO) has stood out as a gradient-free efficient global optimization algorithm that is capable of calibrating constitutive models for CPFEM. Here in this paper, we apply a recently developed asynchronous parallel constrained BO algorithm to calibrate phenomenological constitutive models for stainless steel 304 L, Tantalum, and Cantor high-entropy alloy.

304L stainless steel↗

Drift reduced Landau fluid model for magnetized plasma turbulence simulations in BOUT++ framework

Recently the drift-reduced Landau fluid six-field turbulence model within the BOUT++ framework has been upgraded. In particular, this new model employs a new normalization, adds a volumetric flux-driven source option, the Landau fluid closure for parallel heat flux and a Laplacian inversion solver which is able to capture $n$ = 0 axisymmetric mode evolution in realistic tokamak configurations. As we report here, these improvements substantially extended model's capability to study a wider range of tokamak edge phenomena, and are essential to build a fully self-consistent edge turbulence model capable of both transient (e.g., ELM, disruption) and transport time-scale simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Physics-constrained deep neural network method for estimating parameters in a redox flow battery

Here, in this paper, we present a physics-constrained deep neural network (PCDNN) method for parameter estimation in the zero-dimensional (0D) model of the vanadium redox flow battery (VRFB). In this approach, we use deep neural networks to approximate the model parameters as functions of the operating conditions. This method allows the integration of VRFB computational models as the physical constraints in the parameter learning process, leading to enhanced accuracy of parameter estimation and cell voltage prediction. Using an experimental dataset, we demonstrate that the PCDNN method can estimate model parameters for a range of operating conditions and improve the 0D model prediction of voltage compared to the 0D model prediction with constant operation-condition-independent parameters estimated with traditional inverse methods. We also demonstrate that the PCDNN approach has an improved generalization ability for estimating parameter values for operating conditions not used in the training process.

25 ENERGY STORAGE↗

Summary report of the 4th IAEA Technical Meeting on Fusion Data Processing, Validation and Analysis (FDPVA)

The objective of the Fourth Technical Meeting on Fusion Data Processing, Validation and Analysis was to provide a platform during which a set of topics relevant to fusion data processing, validation and analysis are discussed with the view of extrapolating needs to next step fusion devices such as ITER. The validation and analysis of experimental data obtained from diagnostics used to characterize fusion plasmas are crucial for a knowledge-based understanding of the physical processes governing the dynamics of these plasmas. This paper presents the recent progress and achievements in the domain of plasma diagnostics and synthetic diagnostics data analysis (including image processing, regression analysis, inverse problems, deep learning, machine learning, big data and physics-based models for control) reported at the meeting. The progress in these areas highlight trends observed in current major fusion confinement devices. A special focus is dedicated on data analysis requirements for ITER and DEMO with a particular attention paid to artificial intelligence for automatization and improving reliability of control processes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quantum material topology via defect engineering

Current approaches of discovering new topological quantum materials, such as (magnetic) Weyl materials, by screening material databases or by rational design suffer from experimental challenges in growing the material in the desired geometry, hampering their potential applications. We take an approach of defect-engineering a topology starting from systems with experimentally well-controlled growth . Using SnTe, an experimentally viable topological crystalline insulator, as our model system, we elucidate the different effects of breaking time-reversal, inversion, and mirror symmetries using magnetic dopant defects, to controllably induce a topological phase transition. To this effect, we compute the thermodynamics, magnetism, and band-structure topology of magnetically doped SnTe using a full ab initio approach. In this process, we have discovered that Cr-doped SnTe is a Weyl semimetal. We have computed its classic Fermi-arc experimental signature on the (001) surface and also showed that it has a large anamolous Hall conductivity (AHC) of ~ 250 Ω - 1 cm - 1 . We predict a large Curie temperature ( T c ) of 62 K even at 3.3 % doping, with both T c and AHC being tunable with Cr concentration, suggesting potential for room-temperature applications. Our study in general paves the way for defect-designing room-temperature topological quantum phases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Topological surface states in strained Dirac semimetal thin films

We computationally study the Fermi arc states in a Dirac semimetal, both in a semi-infinite slab and in the thin-film limit. We use Cd 3 As 2 as a model system, and include perturbations that break the C-4 symmetry and inversion symmetry while mirror symmetries remain. The surface states are protected by the mirror symmetries present in the bulk states and thus survive these perturbations. Additionally, the Fermi arc states persist down to very thin films, thinner than presently measured experimentally, but are affected by breaking the symmetry of the Hamiltonian. Recent experimental observations of transport in Cd 3 As 2 films strongly suggest transport by surface states. Our findings are consistent with those observations in that we show that Fermi arc states can survive to rather thin films, and our results also suggest that symmetry-breaking terms that preserve the Fermi arc states nevertheless can have a profound effect in the thin film limit.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Structure retrieval from four-dimensional scanning transmission electron microscopy: Statistical analysis of potential pitfalls in high-dimensional data

Four-dimensional (4D) scanning transmission electron microscopy is one of the most rapidly growing modes of electron microscopy imaging. The advent of fast pixelated cameras and the associated data infrastructure have greatly accelerated this process. Yet conversion of the 4D datasets into physically meaningful structure images in real space remains an open issue. In this work, we demonstrate that it is possible to systematically create filters that will affect the apparent resolution or even qualitative features of the real-space structure image, reconstructing artificially generated patterns. As initial efforts, we explore statistical model selection algorithms, aiming for robustness and reliability of estimated filters. This statistical model selection analysis demonstrates the need for regularization and cross validation of inversion methods to robustly recover structure from high-dimensional diffraction datasets.

36 MATERIALS SCIENCE↗

Almost strong zero modes at finite temperature

Interacting fermionic chains exhibit extended regions of topological degeneracy of their ground states as a result of the presence of Majorana or parafermionic zero modes localized at the edges. In the opposite limit of infinite temperature, the corresponding nonintegrable spin chains, obtained via generalized Jordan-Wigner mapping, are known to host so-called almost strong zero modes, which are long-lived with respect to any bulk excitations. Here we study the fairly unexplored territory that bridges these two extreme cases of zero and infinite temperature. We blend two established techniques for states, the Lanczos series expansion and a tensor network ansatz, uplifting them to the level of operator algebra. This allows us to efficiently simulate large system sizes for arbitrarily long timescales and to extract the temperature-dependent decay rates. We observe that for the Kitaev-Hubbard model, the decay rate of the edge mode depends exponentially on the inverse temperature 𝛽, and on an effective energy scale Δ eff that is greater than the thermodynamic gap of the system Δ.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗