Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Ordering transformation”

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 685 records · Page 38

Wavelets, period-doubling, and time-frequency localization with application to organization of convection over the tropical western Pacific

In this paper, preliminary results in using orthogonal and continuous wavelet transform (WT) to identify period doubling and time-frequency localization in both synthetic and real data are presented. First, the Haar WT is applied to synthetic time series derived from a simple nonlinear dynamical system- a first-order quadratic difference equation. Second, the complex Morlet WT is used to study the time-frequency localization of tropical convection based on a high-resolution Japanese Geostationary Meteorological Satellite infrared (IR) radiance dataset. The Haar WT of the synthetic time series indicates the presence and distinct separation of multiple frequencies in a period-doubling sequence. The period-doubling process generates a multiplicity of intermediate frequencies, which are manifested in the nonuniformity in time with respect to the phase of oscillations in the lower frequencies. Wavelet transform also enables the detection of extremely weak signals in high-order subharmonics resulting from the period-doubling bifurcations. These signals are either undetected or considered statistically insignificant by traditional Fourier analysis. The Morlet WT of the IR radiance dataset indicates the presence of multiple timescales, which are localized in both frequency and time. There are two regimes in the variation of IR radiance, corresponding to the wet and dry periods. Multiple timescales, ranging from semidiurnal, diurnal, synoptic, to intraseasonal with embedding structures, are active in the wet regime. In particular, synoptic variability is more prominent during the wet phase of an intensive intraseasonal cycle. These are not only consistent with, but also show more details than, previous findings by using other techniques. The phase-locking relationships among the oscillations with different time-scales suggest that both synoptic and intraseasonal variations may be mixed oscillations due to the interaction of self-excited oscillations in the tropical atmosphere and external forcings such as annual and diurnal solar radiation variations. Both examples show that WT is a powerful tool for analysis of phenomena involving multiscale interactions that exhibit localization in both frequency and time. A discussion on the caveats in the use of WT in geophysical data analysis is also presented.

Weng, Hengyi↗

The evolution and growth of Al2O3 scales on beta-NiAl

The formation and growth of Al2O3 scales on (beta)-NiAl were studied using electron microscopy and other analytical techniques to gain an understanding of the oxidation properties of (beta)-NiAl and of alumina-forming alloys, in general. The transient and mature stages of oxidation were studied as well as the transformation stage during which the oxide scale transforms from metastable Al2O3 phases to the thermodynamically stable alpha-Al2O3 phase. The transient oxidation stages were studied at 800 deg C and for short times at 1100C. At 800C, the scales consist predominantly of delta-Al2O3 which forms by cation vacancy ordering in the defective spinel lattice of gamma-Al2O3. At 1100C, a fast-growth morphology of theta-Al2O3 forms as a surface layer over delta-Al2O3. For both oxidation temperatures, the scales are often epitaxially oriented with respect to the metal. The transient scales grow by outward cation diffusion as evidenced by surface growth morphologies. The transformation to alpha-Al2O3 occurs within 1 hour at 1100C by a nucleation and radial growth process. The large volume decrease associated with the transformation results in a highly strained alpha-Al2O3 microstructure. A change in scale growth mechanism from outward cation to inward anion diffusion allows transient surface morphologies to be smoothed by surface diffusion. The mature stage of oxidation involves the growth of an alpha-Al2O3 scale having the lacey morphology formed as a result of the gamma yields alpha transformation. Growth of the scale occurs by counterdiffusion along grain boundaries resulting in ridges formed by impingement of alpha-Al2O3 nuclei during the transformation stage. Also scale growth occurs by inward oxygen diffusion through healed cracks; the cracks result from transformation stresses. The measured growth rates of scales having the lacey morphology are an order of magnitude less than fine-grained alpha-Al2O3 scales. Metal orientations were found to have a large effect on oxide morpholoies during all stages of oxidation.

Doychak, J. K.↗

Efficient Fourier transforms for transverse momentum dependent distributions

Hadron production at low transverse momenta in semi-inclusive deep inelastic scattering can be described by transverse momentum dependent (TMD) factorization. This formalism has also been widely used to study the Drell-Yan process and back-to-back hadron pair production in $e^+e^-$ collisions. These processes are the main ones for extractions of TMD parton distribution functions and TMD fragmentation functions, which encode important information about nucleon structure and hadronization. One of the most widely used TMD factorization formalism in phenomenology formulates TMD observables in coordinate $b_\perp$-space, the conjugate space of the transverse momentum. The Fourier transform from $b_\perp$-space back into transverse momentum space is sufficiently complicated due to oscillatory integrands that it requires a careful and computationally intensive numerical treatment in order to avoid potentially large numerical errors. Within the TMD formalism, the azimuthal angular dependence is analytically integrated and the two-dimensional $b_\perp$ integration reduces to a one-dimensional integration over the magnitude $b_\perp$. In this paper we develop a fast numerical Hankel transform algorithm for such a $b_\perp$-integration that improves the numerical accuracy of TMD calculations in all standard processes. Libraries for this algorithm are implemented in Python 2.7 and 3, C++, as well as FORTRAN77. All packages are made available open source.

97 MATHEMATICS AND COMPUTING↗

The role of topography in geodetic gravity field modelling

Masses associated with the topography, bathymetry, and its isostatic compensation are a dominant source of gravity field variations, especially at shorter wavelengths. On global scales the topographic/isostatic effects are also significant, except for the lowest harmonics. In practice, though, global effects need not be taken into account as such effects are included in the coefficients of the geopotential reference fields. On local scales, the short-wavelength gravity variations due to the topography may, in rugged terrain, be an order of magnitude larger than other effects. In such cases, explicit or implicit terrain reduction procedures are mandatory in order to obtain good prediction results. Such effects may be computed by space-domain integration or by fast Fourier transformation (FFT) methods. Numerical examples are given for areas of the Canadian Rockies. In principle, good knowledge of the topographic densities is required to produce the smoothest residual field. Densities may be determined from sample measurements or by gravimetric means, but both are somewhat troublesome methods in practice. The use of a standard density, e.g., 2.67 g/cu cm, may often yield satisfactory results and may be put within a consistent theoretical framework. The independence of density assumptions is the key point of the classical Molodensky approach to the geodetic boundary value problem. The Molodensky solutions take into account that land gravity field observations are done on a non-level surface. Molodensky's problem may be solved by integral expansions or more effective FFT methods, but the solution should not be intermixed with the use of terrain reductions. The methods are actually complimentary and may both be required in order to obtain the smoothest possible signal, least prone to aliasing and other effects coming from sparse data coverage, typical of rugged topography.

Forsberg, R.↗

A new approach to the solution of large, full matrix equations: A two-dimensional potential flow feasibility study

An approach to the solution of matrix problems resulting from integral equations of mathematical physics is presented. Based on the inherent smoothness in such equations, the problem is reformulated using a set of orthogonal basis vectors, leading to an equivalent coefficient problem which can be of lower order without significantly impairing the accuracy of the solution. This approach was evaluated using a two-dimensional Neumann problem describing the inviscid, incompressible flow over an airfoil. Two different kinds of mode functions were investigated, namely eigenfunction series and Fourier series. The method using Fourier series was found preferable. It uses all of the coefficients from a Fast Fourier Transform algorithm in an approximate method which exploits the known structure of the transformed coefficient matrix and very promising results for the flow over a realistic airfoil are obtained. On the basis of the results presented here, an order of magnitude reduction in this computer time can be expected for such problems as compared with the time for a direct matrix solution.

James, R. M.↗

Representations of the language recognition problem for a theorem prover

Two representations of the language recognition problem for a theorem prover in first order logic are presented and contrasted. One of the representations is based on the familiar method of generating sentential forms of the language, and the other is based on the Cocke parsing algorithm. An augmented theorem prover is described which permits recognition of recursive languages. The state-transformation method developed by Cordell Green to construct problem solutions in resolution-based systems can be used to obtain the parse tree. In particular, the end-order traversal of the parse tree is derived in one of the representations. An inference system, termed the cycle inference system, is defined which makes it possible for the theorem prover to model the method on which the representation is based. The general applicability of the cycle inference system to state space problems is discussed. Given an unsatisfiable set S, where each clause has at most one positive literal, it is shown that there exists an input proof. The clauses for the two representations satisfy these conditions, as do many state space problems.

Minker, J.↗

Quantifying the Dynamics of Protein Self-Organization Using Deep Learning Analysis of Atomic Force Microscopy Data

The dynamics of protein self-assembly on the inorganic surface and the resultant geometric patterns are visualized using high-speed atomic force microscopy. The time dynamics of the classical macroscopic descriptors such as 2D fast Fourier transforms, correlation, and pair distribution functions are explored using the unsupervised linear unmixing, demonstrating the presence of static ordered and dynamic disordered phases and establishing their time dynamics. Here, the deep learning (DL)-based workflow is developed to analyze detailed particle dynamics and explore the evolution of local geometries. Finally, we use a combination of DL feature extraction and mixture modeling to define particle neighborhoods free of physics constraints, allowing for a separation of possible classes of particle behavior and identification of the associated transitions. Overall, this work establishes the workflow for the analysis of the self-organization processes in complex systems from observational data and provides insight into the fundamental mechanisms.

36 MATERIALS SCIENCE↗

Investigation of the background in coherent J/ψ production at the EIC

Understanding various fundamental properties of nucleons and nuclei is among the most important scientific goals at the upcoming Electron-Ion Collider (EIC). With the unprecedented opportunity provided by the next-generation machine, the EIC might provide definitive answers to many standing puzzles and open questions in modern nuclear physics. We investigate one of the golden measurements proposed at the EIC, which is to obtain the spatial gluon density distribution within a lead (Pb) nucleus. The proposed experimental process is the exclusive J/ψ vector-meson production off the Pb nucleus: e + Pb → e' + J/ψ + Pb'. The Fourier transformation of the momentum transfer |t| distribution of the coherent diffraction is the transverse gluon spatial distribution. In order to measure it, the experiment has to overcome an overwhelmingly large background arising from the incoherent diffractive production, where the nucleus Pb' mostly breaks up into fragments of particles in the far-forward direction close to the hadron-going beam rapidity. We systematically study the rejection of incoherent J/ψ production by vetoing products from these nuclear breakups—protons, neutrons, and photons—which is based on the BeAGLE event generator and the most up-to-date EIC Far-forward Interaction Region design. The achieved vetoing efficiency, the ratio between the numbers of vetoed events and total incoherent events, ranges from about 80% to 99% depending on |t|. Assuming a 5% smearing applied to the reconstructed |t| resolution in the Sar t re model, this vetoing efficiency can suppress the incoherent background to at least the first minimum of the coherent |t| distribution. Experimental and accelerator machine challenges as well as potential improvements are discussed herein.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A Semi-Analytical Approach for State-Space Electromagnetic Transient Simulation

Here, this paper proposes a semi-analytical approach for efficient and accurate electromagnetic transient (EMT) simulation of a power grid. The approach first derives a high-order semi-analytical solution (SAS) of the grid’s state-space EMT model using the differential transformation (DT), and then evaluates the solution over enlarged, variable time steps to significantly accelerate the simulations while maintaining its high accuracy on detailed fast EMT dynamics. The approach also addresses switches during large time steps by using a limit violation detection algorithm with a binary search-enhanced quadratic interpolation. Case studies are conducted on EMT models of the IEEE 39-bus system and large-scale systems to demonstrate the merits of the new simulation approach against traditional numerical methods.

electromagnetic transient↗

Coupling morphological and magnetic anisotropy for assembling tetragonal colloidal crystals

Morphological and magnetic anisotropy can be combined in colloidal assembly to create unconventional secondary structures. We show here that magnetite nanorods interact along a critical angle, depending on their aspect ratios and assemble into body-centered tetragonal colloidal crystals. Under a magnetic field, size-dependent attractive and repulsive domains develop on the ends and center of the nanorods, respectively. Our joint experiment-computational multiscale study demonstrates the presence of a critical angle in the attractive domain, which defines the equilibrium bonding states of interacting rods and leads to the formation of non–close-packed yet hard-contact tetragonal crystals. Small-angle x-ray scattering measurement attributes the perfect tetragonal phase to the slow assembly kinetics. The crystals exhibit brilliant structural colors, which can be actively tuned by changing the magnetic field direction. These highly ordered frameworks and well-defined three-dimensional nanochannels may offer new opportunities for manipulating nanoscale chemical transformation, mass transportation, and wave propagation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

TPGITF user's manual

The primary purpose of TPG interface (TPGITF) software is to create test patterns, according to a mathematical model, for a real device. In order for the device actually to be tested or to be accurately simulated, a means is required to transform the TPGITF generated test into forms acceptable as input to logic test equipment or to a software logic simulator. The MACRODATA-200 Logic Test System and the LOGSIM logic simulator are the particular systems supported by the TPGITF software. A description is presented of the TPGITF software processor. It is assumed that the reader is familiar with the TPG, LOGSIM, and TOIL systems. Flow charts are given.

Pitts, E. R.↗

Implicit finite-difference procedures for the computation of vortex wakes

Implicit finite-difference procedures for the primitive form of the incompressible Navier-Stokes and the compressible Euler equations are used to compute vortex wake flows. The partial differential equations in strong conservation-law form are transformed to cluster grid points in regions with large changes in vorticity. In addition to clustering, fourth-order accurate, spatial difference operators are used to help resolve the flow-field gradients. The use of implicit time-differencing permits large time steps to be taken since temporal variations are typically small. Computational efficiency is achieved by approximate factorization. Both two-dimensional and preliminary three-dimensional calculations are described and qualitatively compared with existing experimental data.

Steger, J. L.↗

An Eigensystem Realization Algorithm (ERA) for modal parameter identification and model reduction

A method, called the Eigensystem Realization Algorithm (ERA), is developed for modal parameter identification and model reduction of dynamic systems from test data. A new approach is introduced in conjunction with the singular value decomposition technique to derive the basic formulation of minimum order realization which is an extended version of the Ho-Kalman algorithm. The basic formulation is then transformed into modal space for modal parameter identification. Two accuracy indicators are developed to quantitatively identify the system modes and noise modes. For illustration of the algorithm, examples are shown using simulation data and experimental data for a rectangular grid structure.

Juang, J. N.↗

Multigrid applications to three-dimensional elliptic coordinate generation

The multigrid algorithm was applied to solve the coupled set of elliptic quasilinear partial differential equations associated with three-dimensional coordinate generation. The results indicate that the multigrid scheme is more than twice as fast as conventional relaxation schemes on moderate-size grids. Convergence factors of order 0.90 per work unit were achieved on 36,000-point grids. The paper covers the form of transformation, develops the set of generation equations, and gives details on the multigrid approach used. Included are a development of the full-approximation storage scheme, details of the smoothing-rate analysis, and a section devoted to rational programming techniques applicable to the multigrid algorithm.

Thames, F. C.↗

An eigensystem realization algorithm for modal parameter identification and model reduction

A method called the eigensystem realization algorithm is developed for modal parameter identification and model reduction of dynamic systems from test data. A new approach is introduced in conjunction with the singular-value decomposition technique to derive the basic formulation of minimum order realization which is an extended version of the Ho-Kalman algorithm. The basic formulation is then transformed into modal space for modal parameter identification. Two accuracy indicators are developed to quantitatively identify the system and noise modes. For illustration of the algorithm, an example is shown using experimental data from the Galileo spacecraft.

Juang, J.-N.↗

Adiabatic properties of pulsating DA white dwarfs. I - The treatment of the Brunt-Vaisala frequency and the region of period formation

The fundamental issue of the region of period formaton in a degenerate star is examined, with special attention given to the treatment of the Brunt-Vaisala frequency. It is shown that, in order to obtain reliable numerical results in degenerate stellar models, the Brunt-Vaisala frequency must be appropriately transformed, because it is defined in terms of a difference between two numbers which become nearly equal in highly degenerate matter, causing serious numerical problems and systematic errors. An alternative expression is derived, which is valid for multicomponent nonideal partially degenerate and partially ionized plasmas such as those encountered in white dwarf envelopes. This expression is used to compute the period structure of the same white dwarf considered by Pesnell (1987). It is shown that the implicit numerical differencing used in the Lagrangian pulsation code of Pesnell leads to very serious difficulties when used with models of degenerate stars.

Brassard, P.↗

An unsymmetric Lanczos algorithm for damped structural dynamics systems

A one-sided, unsymmetric block Lanczos algorithm is proposed for the model reduction of structural dynamics systems with unsymmetric damping and/or stiffness matrices. The algorithm is a three-term iteration scheme, which transforms the system matrix into an almost skew-symmetric, block-tridiagonal form. The Lanczos reduced-order model is guaranteed to be stable if the full-order system is stable. For unstable systems, a shifting method is available. Also, the algorithm offers flexibility in the choice of starting vectors and thus can yield more accurate reduced-order models. A linear system example and a plane truss structure example are used to show the efficacy of the proposed method.

Su, Tzu-Jeng↗

Structures of aspartic acid-96 in the L and N intermediates of bacteriorhodopsin: analysis by Fourier transform infrared spectroscopy

The light-induced difference Fourier transform infrared spectrum between the L or N intermediate minus light-adapted bacteriorhodopsin (BR) was measured in order to examine the protonated states and the changes in the interactions of carboxylic acids of Asp-96 and Asp-115 in these intermediates. Vibrational bands due to the protonated and unprotonated carboxylic acid were identified by isotope shift and band depletion upon substitution of Asp-96 or -115 by asparagine. While the signal due to the deprotonation of Asp-96 was clearly observed in the N intermediate, this residue remained protonated in L. Asp-115 was partially deprotonated in L. The C = O stretching vibration of protonated Asp-96 of L showed almost no shift upon 2H2O substitution, in contrast to the corresponding band of Asp-96 or Asp-115 of BR, which shifted by 9-12 cm-1 under the same conditions. In the model system of acetic acid in organic solvents, such an absence of the shift of the C = O stretching vibration of the protonated carboxylic acid upon 2H2O substitution was seen only when the O-H of acetic acid is hydrogen-bonded. The non-hydrogen-bonded monomer showed the 2H2O-dependent shift. Thus, the O-H bond of Asp-96 enters into hydrogen bonding upon conversion of BR to L. Its increased hydrogen bonding in L is consistent with the observed downshift of the O-H stretching vibration of the carboxylic acid of Asp-96.

NASA Discipline Exobiology↗