SEARCH · Search NASA
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.
Iterative method for systems of two-point boundary value differential equations based on contraction mapping principle
Developments in coercive field engineering and waveguide writing via ion-exchange in the ferroelectric ionic-conductor KTiOPO4 (KTP) have led to breakthroughs within the field of nonlinear optics and quantum photonics. However, the microscopic diffusion dynamics of ion-exchange in KTP are poorly understood and show signs of non-classical diffusion dynamics. The precise control of ion-exchange gradients is critical for developing KTP into a mature, integrated nonlinear optical material but requires a level of control that, in turn, demands a deeper understanding of the ion-exchange process. We employ direct ToF-SIMS mapping to untangle the complex diffusion dynamics of the commonly used ion-exchange dopant ions: Rb, Ba, and K. We map periodic diffusion gradients of two exchanges, one dominated by Rb-ions and the other by Ba-ions. For the first time, we show the interdependent nature of the Ba- and Rb ion in the ion-exchange diffusion gradient. The low incorporation of the Ba ion is mapped and quantified, and dopant–crystal interactions are identified. This work highlights the complex diffusion dynamics of dopants in the KTP structure and lays the foundation for deterministic control of diffusion gradients for the next generation of integrated nonlinear and quantum photonic devices.
A recently developed algorithm for demodulation and decoding of a pulse-position- modulation (PPM) signal is suitable as a basis for designing a single hardware decoding apparatus to be capable of handling any PPM order. Hence, this algorithm offers advantages of greater flexibility and lower cost, in comparison with prior such algorithms, which necessitate the use of a distinct hardware implementation for each PPM order. In addition, in comparison with the prior algorithms, the present algorithm entails less complexity in decoding at large orders. An unavoidably lengthy presentation of background information, including definitions of terms, is prerequisite to a meaningful summary of this development. As an aid to understanding, the figure illustrates the relevant processes of coding, modulation, propagation, demodulation, and decoding. An M-ary PPM signal has M time slots per symbol period. A pulse (signifying 1) is transmitted during one of the time slots; no pulse (signifying 0) is transmitted during the other time slots. The information intended to be conveyed from the transmitting end to the receiving end of a radio or optical communication channel is a K-bit vector u. This vector is encoded by an (N,K) binary error-correcting code, producing an N-bit vector a. In turn, the vector a is subdivided into blocks of m = log2(M) bits and each such block is mapped to an M-ary PPM symbol. The resultant coding/modulation scheme can be regarded as equivalent to a nonlinear binary code. The binary vector of PPM symbols, x is transmitted over a Poisson channel, such that there is obtained, at the receiver, a Poisson-distributed photon count characterized by a mean background count nb during no-pulse time slots and a mean signal-plus-background count of ns+nb during a pulse time slot. In the receiver, demodulation of the signal is effected in an iterative soft decoding process that involves consideration of relationships among photon counts and conditional likelihoods of m-bit vectors of coded bits. Inasmuch as the likelihoods of all the m-bit vectors of coded bits mapping to the same PPM symbol are correlated, the best performance is obtained when the joint mbit conditional likelihoods are utilized. Unfortunately, the complexity of decoding, measured in the number of operations per bit, grows exponentially with m, and can thus become prohibitively expensive for large PPM orders. For a system required to handle multiple PPM orders, the cost is even higher because it is necessary to have separate decoding hardware for each order. This concludes the prerequisite background information. In the present algorithm, the decoding process as described above is modified by, among other things, introduction of an lbit marginalizer sub-algorithm. The term "l-bit marginalizer" signifies that instead of m-bit conditional likelihoods, the decoder computes l-bit conditional likelihoods, where l is fixed. Fixing l, regardless of the value of m, makes it possible to use a single hardware implementation for any PPM order. One could minimize the decoding complexity and obtain an especially simple design by fixing l at 1, but this would entail some loss of performance. An intermediate solution is to fix l at some value, greater than 1, that may be less than or greater than m. This solution makes it possible to obtain the desired flexibility to handle any PPM order while compromising between complexity and loss of performance.
The accuracy of adaptively chosen, mapped polynomial approximations is studied for functions with steep gradients or discontinuities. It is shown that, for steep gradient functions, one can obtain spectral accuracy in the original coordinate system by using polynomial approximations in a transformed coordinate system with substantially fewer collocation points than are necessary using polynomial expansion directly in the original, physical, coordinate system. It is also shown that one can avoid the usual Gibbs oscillation associated with steep gradient solutions of hyperbolic pde's by approximation in suitably chosen coordinate systems. Continuous, high gradient solutions are computed with spectral accuracy (as measured in the physical coordinate system). Discontinuous solutions associated with nonlinear hyperbolic equations can be accurately computed by using an artificial viscosity chosen to smooth out the solution in the mapped, computational domain. Thus, shocks can be effectively resolved on a scale that is subgrid to the resolution available with collocation only in the physical domain. Examples with Fourier and Chebyshev collocation are given.
The accuracy of adaptively chosen, mapped polynomial approximations is studied for functions with steep gradients or discontinuities. It is shown that, for steep gradient functions, one can obtain spectral accuracy in the original coordinate system by using polynomial approximations in a transformed coordinate system with substantially fewer collocation points than are necessary using polynomial expansion directly in the original, physical, coordinate system. It is also shown that one can avoid the usual Gibbs oscillation associated with steep gradient solutions of hyperbolic pde's by approximation in suitably chosen coordinate systems. Continuous, high gradient solutions are computed with spectral accuracy (as measured in the physical coordinate system). Discontinuous solutions associated with nonlinear hyperbolic equations can be accurately computed by using an artificial viscosity chosen to smooth out the solution in the mapped, computational domain. Thus, shocks can be effectively resolved on a scale that is subgrid to the resolution available with collocation only in the physical domain. Examples with Fourier and Chebyshev collocation are given.
We have designed a cubic spline wavelet decomposition for the Sobolev space H(sup 2)(sub 0)(I) where I is a bounded interval. Based on a special 'point-wise orthogonality' of the wavelet basis functions, a fast Discrete Wavelet Transform (DWT) is constructed. This DWT transform will map discrete samples of a function to its wavelet expansion coefficients in O(N log N) operations. Using this transform, we propose a collocation method for the initial value boundary problem of nonlinear PDE's. Then, we test the efficiency of the DWT transform and apply the collocation method to solve linear and nonlinear PDE's.
Imaging laser-induced rotational dynamics is an important and active field due to its applications in capturing reactions in the molecular frame and in molecular imaging. Experimental measurement of the molecular orientation distribution, as a function of the Euler angles, has only been demonstrated for special cases when the detectable signal is generated along the molecular symmetry axis. Here we developed the general theory that maps the probability density distribution of the molecular orientation to the atom-pair angular distributions for nonlinear molecules. With the theory, the molecular orientation distribution can be retrieved from the measured atom-pair angular distribution, which we demonstrate experimentally using ultrafast electron diffractive imaging of impulsively aligned trifluoro-iodomethane molecules. Here, the retrieved molecular orientation distribution is in good agreement with direct numerical simulations of the time-dependent Schrodinger equation using the experimental conditions. Unlike the existing retrieval methods, the retrieval method does not require solving Schrodinger equation, works for any alignment method, and is in principle applicable to asymmetric top molecules.
Doppler tracking data from the Pioneer Venus Orbiter (PVO) have been used to estimate the anomalous gravity field in the region of Venus west of Beta Regio. The analysis invokes a Kalman filter-smoother to solve the nonlinear spacecraft state estimation problem and a linear Bayesian estimator to perform the geophysical inversion. The topographic map for this region, derived from the PVO radar, has been filtered to have the same distortions and degree of smoothing as the gravity map. The undulations of the gravity are about 0.2 times as large as expected from the topography, on the assumption that the latter is uncompensated. A comparison of the gravity and topography by means of the spectral admittance is consistent with Airy compensation at a depth of 50 km, if the surface material has a density of 2.6 g/cu cm. However, this is not a unique interpretation.
Efficient compression of image data requires the understanding of the noise characteristics of sensors as well as the redundancy expected in imagery. Herein, the techniques of Differential Pulse Code Modulation (DPCM) are reviewed and modified for information-preserving data compression. The modifications include: mapping from intensity to an equal variance space; context dependent one and two dimensional predictors; rationale for nonlinear DPCM encoding based upon an image quality model; context dependent variable length encoding of 2x2 data blocks; and feedback control for constant output rate systems. Examples are presented at compression rates between 1.3 and 2.8 bits per pixel. The need for larger block sizes, 2D context dependent predictors, and the hope for sub-bits-per-pixel compression which maintains spacial resolution (information preserving) are discussed.
The increasing electrification of aircraft propulsion systems is leading to new control architectures being developed to address integration between electric machines and gas-based turbine engines. For hybrid-electric propulsion systems, current conceptual architectures often couple electric machines with the shafts of gas turbine engines and introduce energy storage. Leveraging the electrical power system of hybridized engines, Turbine Electrified Energy Management (TEEM) is a recent control approach that improves transient operability in an effort to enable more efficient and lighter weight turbomachinery. This study seeks to expand TEEM’s application beyond traditional proportional-integral (PI) control by presenting linear model predictive control (MPC) schemes to execute the TEEM concept. Through constraint selection and cost function design, transient operability goals for TEEM are considered with no external logic or saturation. Unique to the designs are the use of a washout filter, which simplifies transient detection and motor activation logic. The proposed architectures are implemented with both centralized MPC and distributed MPC approaches, and comparisons are drawn to a benchmark PI controller simulated on a nonlinear turbofan engine model at one ground condition and one cruise condition. Performance is evaluated using compressor maps, stall margin performance, and two novel metrics: transient stack usage and transient excursion integral. Results reveal the linear MPC scheme performs comparably to the baseline controller and can be implemented in at least two distinct configurations with potential for further modifications, thus establishing the groundwork for future investigations.
The increasing electrification of aircraft propulsion systems is leading to new control architectures being developed to address integration between electric machines and gas-based turbine engines. For hybrid-electric propulsion systems, current conceptual architectures often couple electric machines with the shafts of gas turbine engines and introduce energy storage. Leveraging the electrical power system of hybridized engines, Turbine Electrified Energy Management (TEEM) is a recent control approach that improves transient operability in an effort to enable more efficient and lighter weight turbomachinery. This study seeks to expand TEEM’s application beyond traditional proportional integral (PI) control by presenting linear model predictive control (MPC) schemes to execute the TEEM concept. Through constraint selection and cost function design, transient operability goals for TEEM are considered with no external logic or saturation. Unique to the designs are the use of a washout filter, which simplifies transient detection and motor activation logic. The proposed architectures are implemented with both centralized MPC and distributed MPC approaches, and comparisons are drawn to a benchmark PI controller simulated on a nonlinear turbofan engine model at one ground condition and one cruise condition. Performance is evaluated using compressor maps, stall margin performance, and two novel metrics: transient stack usage and transient excursion integral. Results reveal the linear MPC scheme performs comparably to the baseline controller and can be implemented in at least two distinct configurations with potential for further modifications, thus establishing the groundwork for future investigations.
The increasing electrification of aircraft propulsion systems is leading to new control architectures being developed to address integration between electric machines and gas-based turbine engines. For hybrid-electric propulsion systems, current conceptual architectures often couple electric machines with the shafts of gas turbine engines and introduce energy storage. Leveraging the electrical power system of hybridized engines, Turbine Electrified Energy Management (TEEM) is a recent control approach that improves transient operability in an effort to enable more efficient and lighter weight turbomachinery. This study seeks to expand TEEM’s application beyond traditional proportional-integral (PI) control by presenting linear model predictive control (MPC) schemes to execute the TEEM concept. Through constraint selection and cost function design, transient operability goals for TEEM are considered with no external logic or saturation. Unique to the designs are the use of a washout filter, which simplifies transient detection and motor activation logic. The proposed architectures are implemented with both centralized MPC and distributed MPC approaches, and comparisons are drawn to a benchmark PI controller simulated on a nonlinear turbofan engine model at one ground condition and one cruise condition. Performance is evaluated using compressor maps, stall margin performance, and two novel metrics: transient stack usage and transient excursion integral. Results reveal the linear MPC scheme performs comparably to the baseline controller and can be implemented in at least two distinct configurations with potential for further modifications, thus establishing the groundwork for future investigations.
Abstract Predicting complex dynamics in physical applications governed by partial differential equations in real-time is nearly impossible with traditional numerical simulations due to high computational cost. Neural operators offer a solution by approximating mappings between infinite-dimensional Banach spaces, yet their performance degrades with system size and complexity. We propose an approach for learning neural operators in latent spaces, facilitating real-time predictions for highly nonlinear and multiscale systems on high-dimensional domains. Our method utilizes the deep operator network architecture on a low-dimensional latent space to efficiently approximate underlying operators. Demonstrations on material fracture, fluid flow prediction, and climate modeling highlight superior prediction accuracy and computational efficiency compared to existing methods. Notably, our approach enables approximating large-scale atmospheric flows with millions of degrees, enhancing weather and climate forecasts. Here we show that the proposed approach enables real-time predictions that can facilitate decision-making for a wide range of applications in science and engineering.
An analysis of phenomena causing nonlinearities in the transformation from Landsat multispectral scanner coordinates to ground coordinates is presented. Experimental results comparing rms errors at ground control points indicated a slight improvement when a nonlinear (8-parameter) transformation was used instead of an affine (6-parameter) transformation. Using a preliminary ground truth map of a test site in Alabama covering the Mobile Bay area and six Landsat images of the same scene, several classification methods were assessed. A methodology was developed for automatic change detection using classification/cluster maps. A coding scheme was employed for generation of change depiction maps indicating specific types of changes. Inter- and intraseasonal data of the Mobile Bay test area were compared to illustrate the method. A beginning was made in the study of data compression by applying a Karhunen-Loeve transform technique to a small section of the test data set. The second part of the report provides a formal documentation of the several programs developed for the analysis and assessments presented.
The pitch motion of a spacecraft in the planar elliptic restricted three-body system is studied. Previous studies laid the foundation for spacecraft stability analysis with a small perturbation to the zero pitch motion. In this study, a cell mapping approach that combines analytical and numerical techniques is used to study the global behavior of the full nonlinear spacecraft attitude in which coupling between orbital dynamics and attitude occurs. The effect of gravity gradient torques, orbital eccentricity, and the spacecraft configuration at different Lagrangian points is analyzed. Multiple-period periodic solutions and invariant surfaces are presented for different cases. Reference trajectories around the Lagrangian points are also considered to study coupled dynamics.
Constructing surrogate models for uncertainty quantification (UQ) on complex partial differential equations (PDEs) having inherently high-dimensional O(10 n ), n ≥ 2, stochastic inputs (e.g., forcing terms, boundary conditions, initial conditions) poses tremendous challenges. The “curse of dimensionality” can be addressed with suitable unsupervised learning techniques used as a pre-processing tool to encode inputs onto lower-dimensional subspaces while retaining its structural information and meaningful properties. In this work, we review and investigate thirteen dimension reduction methods including linear and nonlinear, spectral, blind source separation, convex and non-convex methods and utilize the resulting embeddings to construct a mapping to quantities of interest via polynomial chaos expansions (PCE). Here, we refer to the general proposed approach as manifold PCE (m-PCE), where manifold corresponds to the latent space resulting from any of the studied dimension reduction methods. To investigate the capabilities and limitations of these methods we conduct numerical tests for three physics-based systems (treated as black-boxes) having high-dimensional stochastic inputs of varying complexity modeled as both Gaussian and non-Gaussian random fields to investigate the effect of the intrinsic dimensionality of input data. We demonstrate both the advantages and limitations of the unsupervised learning methods and we conclude that a suitable m-PCE model provides a cost-effective approach compared to alternative algorithms proposed in the literature, including recently proposed expensive deep neural network-based surrogates and can be readily applied for high-dimensional UQ in stochastic PDEs.
The navigation of imagery from polar orbiting weather satellites includes the correction for geometric distortions due to earth shape/earth rotation, satellite orbit variations, and satellite attitude along with the resampling of the satellite image to a selected geographic map projection. The routine image navigation procedure also compensates for the distortion of the satellite imagery due to the nonlinear scanning of the sensor system. This paper reviews general methods for performing this image navigation, ranging from a method that assumes no orbital information and, thus, relies on nominal orbital parameter values and image corrections computed by matching ground control points (GCPs), to a method that uses high-quality satellite ephemeris data to make the correction with a limited number of GCPs. A procedure to optimize the image navigation by using a spatial remapping, or interpolation, is introduced and outlined. Recommendations are made for people interested in the processing of AVHRR imagery.
The chaos in a one-dimensional system, which would be nonlinear stationary Alfven waves in the absence of an external driver, is characterized. The evolution equations are numerically integrated for the transverse wave magnetic field amplitude and phase using the derivative nonlinear Schroedinger equation (DNLS), including resistive wave damping and a long-wavelength monochromatic, circularly polarized driver. A Poincare map analysis shows that, for the nondissipative (Hamiltonian) case, the solutions near the phase space (soliton) separatrices of this system become chaotic as the driver amplitude increases, and 'strong' chaos appears when the driver amplitude is large. The dissipative system exhibits a wealth of dynamical behavior, including quasiperiodic orbits, period-doubling bifurcations leading to chaos, sudden transitions to chaos, and several types of strange attractors.