Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Linearity preservation”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1972

Decay of isotropic turbulence generated by a mechanically agitated grid.

Experimental study of weak isotropic turbulence, created by a mechanically agitated grid, has indicated that in the absence of large linear-momentum wakes the energy of turbulence relaxes very quickly into a stable self-preserving structure, which, depending on the initial Reynolds number of turbulence, decays at different constant inverse powers of time. Both the longitudinal correlation coefficients and the corresponding spectral distributions, except for the difference in the parametric constants, are of the same functional type as those found previously for a passive grid.

Ling, S. C.↗
NASA NTRS1976

Feature combinations and the divergence criterion

Classifying large quantities of multidimensional remotely sensed agricultural data requires efficient and effective classification techniques and the construction of certain transformations of a dimension reducing, information preserving nature. The construction of transformations that minimally degrade information (i.e., class separability) is described. Linear dimension reducing transformations for multivariate normal populations are presented. Information content is measured by divergence.

Decell, H. P., Jr.↗
NASA NTRS1978

Linear dimension reduction and Bayes classification

An explicit expression for a compression matrix T of smallest possible left dimension K consistent with preserving the n variate normal Bayes assignment of X to a given one of a finite number of populations and the K variate Bayes assignment of TX to that population was developed. The Bayes population assignment of X and TX were shown to be equivalent for a compression matrix T explicitly calculated as a function of the means and covariances of the given populations.

Decell, H. P., Jr.↗
DOE OSTI2023

Quantum Langevin theory for two coupled phase-conjugated electromagnetic waves

We provide a general macroscopic phenomenological formula of quantum Langevin equations for two coupled phase-conjugated electromagnetic fields with linear loss (gain) and complex nonlinear coupling coefficient. The macroscopic phenomenological formula is obtained from the coupling matrix to preserve the field commutation relations and correlations, which does not require knowing the microscopic details of light-matter interaction and internal atomic structures. To validate this phenomenological formula, we take spontaneous four-wave mixing in a double-Λ four-level atomic system as an example to numerically confirm that our macroscopic phenomenological result is consistent with that obtained from the microscopic Heisenberg-Langevin theory. We find that a complex-valued nonlinear coupling coefficient can lead to noises even without linear gain or loss. Lastly, we apply the quantum Langevin equations to study the effects of linear gain and loss, complex phase mismatching, as well as complex nonlinear coupling coefficient in entangled photon pair (biphoton) generation, particularly to their temporal quantum correlations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

One-Body Properties and Their Perturbative Accuracy with Aufbau Suppressed Coupled Cluster Theory

In this work we derived and implemented the calculation of the one-body reduced density matrix for Aufbau suppressed coupled cluster theory, from which excited state natural orbitals and one-body properties, like atomic populations and dipole moments, are obtained. We utilized the natural orbitals to refine the ASCC solution for simple valence and Rydberg systems, exploring the process of repeatedly solving the ASCC equations in successive natural orbital bases to achieve independence from the starting molecular orbitals. For dipole moments in small molecules where high-level comparison data is available, we find that the accuracy of ASCC essentially matches that of linear response and equation-of-motion coupled cluster as long as care is taken to preserve the response's perturbative completeness.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Extraordinary Magnetic Response of an Anisotropic 2D Antiferromagnet via Site Dilution

A prominent characteristic of 2D magnetic systems is the enhanced spin fluctuations, which reduce the ordering temperature. We report that a magnetic field of only 1000th of the Heisenberg superexchange interaction can induce a crossover, which for practical purposes is the effective ordering transition, at temperatures about 6 times the Néel transition in a site-diluted two-dimensional anisotropic quantum antiferromagnet. Such a strong magnetic response is enabled because the system directly enters the antiferromagnetically ordered state from the isotropic disordered state, skipping the intermediate anisotropic stage. The underlying mechanism is achieved on a pseudospin-half square lattice realized in the [(SrIrO3)1/(SrTiO3)2] superlattice thin film that is designed to linearly couple the staggered magnetization to external magnetic fields by virtue of the rotational symmetry-preserving Dzyaloshinskii–Moriya interaction. Furthermore, our model analysis shows that the skipping of the anisotropic regime despite finite anisotropy is due to the enhanced isotropic fluctuations under moderate dilution.

2D antiferromagnet↗
NASA NTRS2011

Turbulence Modeling for Shock Wave/Turbulent Boundary Layer Interactions

Accurate aerodynamic computational predictions are essential for the safety of space vehicles, but these computations are of limited accuracy when large pressure gradients are present in the flow. The goal of the current project is to improve the state of compressible turbulence modeling for high speed flows with shock wave / turbulent boundary layer interactions (SWTBLI). Emphasis will be placed on models that can accurately predict the separated region caused by the SWTBLI. These flows are classified as nonequilibrium boundary layers because of the very large and variable adverse pressure gradients caused by the shock waves. The lag model was designed to model these nonequilibrium flows by incorporating history effects. Standard one- and two-equation models (Spalart Allmaras and SST) and the lag model will be run and compared to a new lag model. This new model, the Reynolds stress tensor lag model (lagRST), will be assessed against multiple wind tunnel tests and correlations. The basis of the lag and lagRST models are to preserve the accuracy of the standard turbulence models in equilibrium turbulence, when the Reynolds stresses are linearly related to the mean strain rates, but create a lag between mean strain rate effects and turbulence when nonequilibrium effects become important, such as in large pressure gradients. The affect this lag has on the results for SWBLI and massively separated flows will be determined. These computations will be done with a modified version of the OVERFLOW code. This code solves the RANS equations on overset grids. It was used for this study for its ability to input very complex geometries into the flow solver, such as the Space Shuttle in the full stack configuration. The model was successfully implemented within two versions of the OVERFLOW code. Results show a substantial improvement over the baseline models for transonic separated flows. The results are mixed for the SWBLI assessed. Separation predictions are not as good as the baseline models, but the over prediction of the peak heat flux downstream of the reattachment shock that plagues many models is reduced.

Lillard, Randolph P.↗
DOE OSTI2021

Adiabatic matching of particle bunches in a plasma-based accelerator in the presence of ion motion

Witness beam stability and preservation of its ultra-low emittance have been identified as critical challenges toward realizing a TeV-class, plasma-based linear collider. In fact, the witness bunch parameters required by a future TeV-class collider have been expected to trigger hosing instability and background ion motion, leading to emittance degradation and, potentially, to bunch loss. Recently, it has been shown that ion motion suppresses the hosing instability, and proper longitudinal bunch shaping can eliminate the ion-motion-induced emittance growth. Here, we propose and analyze a plasma-based method to generate the shaped bunches that enable emittance preservation in the presence of ion motion. The method is based on an adiabatic matching procedure, where a bunch with an initially untapered profile is injected in the plasma accelerator stage with an energy low enough that ion motion effects are initially small. As the bunch accelerates, it is adiabatically compressed, ion motion is gradually triggered, and the bunch slowly but continuously readjusts itself in the ion motion-perturbed wakefield acquiring the desired taper. The production of tapered witness bunch profiles that minimize energy spread and preserve emittance for collider-relevant parameters could enable the use of plasma-based accelerators for high-energy physics applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A method for bounding high-order finite element functions: Applications to mesh validity and bounds-preserving limiters

We introduce a novel method for bounding high-order multi-dimensional polynomials in finite element approximations. The method involves precomputing optimal piecewise-linear bounding boxes for polynomial basis functions, which can then be used to locally bound any combination of these basis functions. This approach can be applied to any element/basis type at any approximation order, can provide local (i.e., subcell) extremum bounds to a desired level of accuracy, and can be evaluated efficiently on-the-fly in simulations. Furthermore, we show that this approach generally yields more accurate bounds in comparison to traditional methods based on convex hull properties (e.g., Bernstein polynomials). Furthermore, the efficacy of this technique is shown in applications such as mesh validity checks and optimization for high-order curved meshes, where positivity of the element Jacobian determinant can be ensured throughout the entire element, and continuously bounds-preserving limiters for hyperbolic systems, which can enforce maximum principle bounds across the entire solution polynomial.

Bounding box↗
NASA NTRS2011

Hyperspectral Remote Sensing of Atmosphere and Surface Properties

Atmospheric Infrared Sounder (AIRS), Infrared Atmospheric Sounding Interferometer (IASI), and Cross-track Infrared Sounder (CrIS) are all hyper-spectral satellite sensors with thousands of spectral channels. Top of atmospheric radiance spectra measured by these sensors contain high information content on atmospheric, cloud, and surface properties. Exploring high information content contained in these high spectral resolution spectra is a challenging task due to computation e ort involved in modeling thousands of spectral channels. Usually, only very small fractions (4{10 percent) of the available channels are included in physical retrieval systems or numerical weather forecast (NWP) satellite data assimilations. We will describe a method of simultaneously retrieving atmospheric temperature, moisture, cloud, and surface properties using all available spectral channels without sacrificing computational speed. The essence of the method is to convert channel radiance spectra into super-channels by an Empirical Orthogonal Function (EOF) transformation. Because the EOFs are orthogonal to each other, about 100 super-channels are adequate to capture the information content of the radiance spectra. A Principal Component-based Radiative Transfer Model (PCRTM) developed at NASA Langley Research Center is used to calculate both the super-channel magnitudes and derivatives with respect to atmospheric profiles and other properties. There is no need to perform EOF transformations to convert super channels back to spectral space at each iteration step for a one-dimensional variational retrieval or a NWP data assimilation system. The PCRTM forward model is also capable of calculating radiative contributions due to multiple-layer clouds. The multiple scattering effects of the clouds are efficiently parameterized. A physical retrieval algorithm then performs an inversion of atmospheric, cloud, and surface properties in super channel domain directly therefore both reducing the computational need and preserving the information content of the IASI measurements. The inversion algorithm is based on a non-linear Levenberg-Marquardt method with climatology covariance matrices and a priori information as constraints. One advantage of this approach is that it uses all information content from the hyper-spectral data so that the retrieval is less sensitive to instrument noise and eliminates the need for selecting a sub-set of the channels.

Liu, Xu↗
NASA NTRS2006

Sulfur during the Transition from Anoxic to Oxic Atmospheres

The invention of oxygenic photosynthesis was likely accompanied by the introduction of large amounts of O2 and complementary reduced gases (chiefly CH4) into the atmosphere. To first approximation the venting of O2 and CH4 are stochiometrically linked. We therefore present a suite of numerical photochemical models that address the anoxic-oxic transition in an atmosphere driven by large linked inputs of biogenic 02 and CH4. We find in general that, in steady state, there are two solutions, one oxic and the other anoxic. The anoxic solution appears to be linearly stable. If volcanic SO2 fluxes are large, S disproportionates into oxidized (H2S04) and reduced (S8) exit channels. As elemental sulfur is insoluble it provides a means of preserving photochemical mass-independent fractionation (MIF). On the other hand, if the source of volcanic SO2 is smaller than today, all S can leave the atmosphere as S8. Under these conditions there would be no MIF signal. The oxic solution appears to be linearly unstable. In the oxic solutions S is invariably oxidized to sulfate, and the MIF signal would be absent. The transitional atmosphere is relatively unstable and is also the most photochemically active. Consequently it is the transitional atmosphere, not the oxic or anoxic atmospheres, that has the lowest CH4 levels and weakest greenhouse warming. As a practical matter we expect the transitional atmospheres to vary strongly in response to diurnal and seasonal biological forcing.

Zahnle, Kevin↗
DOE OSTI2022

Local bilinear computation of Jacobi sets

Abstract We propose a novel method for the computation of Jacobi sets in 2D domains. The Jacobi set is a topological descriptor based on Morse theory that captures gradient alignments among multiple scalar fields, which is useful for multi-field visualization. Previous Jacobi set computations use piecewise linear approximations on triangulations that result in discretization artifacts like zig-zag patterns. In this paper, we utilize a local bilinear method to obtain a more precise approximation of Jacobi sets by preserving the topology and improving the geometry. Consequently, zig-zag patterns on edges are avoided, resulting in a smoother Jacobi set representation. Our experiments show a better convergence with increasing resolution compared to the piecewise linear method. We utilize this advantage with an efficient local subdivision scheme. Finally, our approach is evaluated qualitatively and quantitatively in comparison with previous methods for different mesh resolutions and across a number of synthetic and real-world examples.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Canonical and noncanonical Hamiltonian operator inference

Here, a method for the nonintrusive and structure-preserving model reduction of canonical and noncanonical Hamiltonian systems is presented. Based on the idea of operator inference, this technique is provably convergent and reduces to a straightforward linear solve given snapshot data and gray-box knowledge of the system Hamiltonian. Examples involving several hyperbolic partial differential equations show that the proposed method yields reduced models which, in addition to being accurate and stable with respect to the addition of basis modes, preserve conserved quantities well outside the range of their training data.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1999

The Reconstruction Problem Revisited

The role of reconstruction in avoiding oscillations in upwind schemes is reexamined, with the aim of providing simple, concise proofs. In one dimension, it is shown that if the reconstruction is any arbitrary function bounded by neighboring cell averages and increasing within a cell for increasing data, the resulting scheme is monotonicity preserving, even though the reconstructed function may have overshoots and undershoots at the cell edges and is in general not a monotone function. In the special case of linear reconstruction, it is shown that merely bounding the reconstruction between neighboring cell averages is sufficient to obtain a monotonicity preservinc,y scheme. In two dimensions, it is shown that some ID TVD limiters applied in each direction result in schemes that are not positivity preserving, i.e. do not give positive updates when the data are positive. A simple proof is given to show that if the reconstruction inside the cell is bounded by the neighboring cell averages (including corner neighbors), then the scheme is positivity preserving. A new limiter that enforces this condition but is not as dissipative as the Minmod limiter is also presented.

Suresh, Ambaby↗
DOE OSTI2020

SympNets: Intrinsic structure-preserving symplectic networks for identifying Hamiltonian systems

Here, we propose new symplectic networks (SympNets) for identifying Hamiltonian systems from data based on a composition of linear, activation and gradient modules. In particular, we define two classes of SympNets: the LA-SympNets composed of linear and activation modules, and the G-SympNets composed of gradient modules. Correspondingly, we prove two new universal approximation theorems that demonstrate that SympNets can approximate arbitrary symplectic maps based on appropriate activation functions. We then perform several experiments including the pendulum, double pendulum and three-body problems to investigate the expressivity and the generalization ability of SympNets. The simulation results show that even very small size SympNets can generalize well, and are able to handle both separable and non-separable Hamiltonian systems with data points resulting from short or long time steps. In all the test cases, SympNets outperform the baseline models, and are much faster in training and prediction. We also develop an extended version of SympNets to learn the dynamics from irregularly sampled data. This extended version of SympNets can be thought of as a universal model representing the solution to an arbitrary Hamiltonian system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Site-controlled telecom-wavelength single-photon emitters in atomically-thin MoTe2

Abstract Quantum emitters (QEs) in two-dimensional transition metal dichalcogenides (2D TMDCs) have advanced to the forefront of quantum communication and transduction research. To date, QEs capable of operating in O-C telecommunication bands have not been demonstrated in TMDCs. Here we report site-controlled creation of telecom QEs emitting over the 1080 to 1550 nm telecommunication wavelength range via coupling of 2D molybdenum ditelluride (MoTe 2 ) to strain inducing nano-pillar arrays. Hanbury Brown and Twiss experiments conducted at 10 K reveal clear photon antibunching with 90% single-photon purity. The photon antibunching can be observed up to liquid nitrogen temperature (77 K). Polarization analysis further reveals that while some QEs display cross-linearly polarized doublets with ~1 meV splitting resulting from the strain induced anisotropic exchange interaction, valley degeneracy is preserved in other QEs. Valley Zeeman splitting as well as restoring of valley symmetry in cross-polarized doublets are observed under 8 T magnetic field.

36 MATERIALS SCIENCE↗
NASA NTRS1991

Model order reduction applied to a hot-bench simulation of an aeroelastic wind-tunnel model

Simulations of an aeroelastically scaled wind-tunnel model were developed for hot-bench testing of a digital controller. The digital controller provided active flutter-suppression, rolling-maneuver-load alleviation, and plant estimation. To achieve an acceptable time scale for the hot-bench application, the mathematical model of the wind-tunnel model was reduced from 220 states to approximately 130 states while assuring that the required accuracy was preserved for all combinations of 10 inputs and 56 outputs. The reduction was achieved by focussing on a linear, aeroelastic submodel of the full mathematical model and by applying a method based on the internally balanced realization of a dynamic system. The error-bound properties of the internally balanced realization significantly contribute to its utility in the model reduction process. The reduction method and the results achieved are described.

Buttrill, Carey S.↗
NASA NTRS1992

Coherent-backscatter effect - A vector formulation accounting for polarization and absorption effects and small or large scatterers

Previous theoretical work on the coherent-backscatter effect in the context of speckle time autocorrelation has gone beyond the diffusion approximation and the assumption of isotropic (point) scatterers. This paper extends the theory to include the effects of polarization and absorption, and to give the angular line shape. The results are expressions for angular variations valid for small and large scatterers and linear and circular polarizations, in lossless or lossy media. Calculations show that multiple anisotropic scattering results in the preservation of incident polarization. Application to a problem in radar astronomy is considered. It is shown that the unusual radar measurements (high reflectivity and polarization ratios) of Jupiter's icy Galilean satellites can be explained by coherent backscatter from anisotropic (forward) scatterers.

Peters, Kenneth J.↗