Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ORTHOGONAL FUNCTION”

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

Quantifying vorticity in magnetic particle suspensions driven by symmetric and asymmetric multiaxial fields.

We recently reported two methods of inducing vigorous fluid vorticity in magnetic particle suspensions. The first method employs symmetry-breaking rational fields. These fields are comprised of two orthogonal ac components whose frequencies form a rational number and an orthogonal dc field that breaks the symmetry of the biaxial ac field to create the parity required to induce deterministic vorticity. The second method is based on rational triads, which are fields comprised of three orthogonal ac components whose frequency ratios are rational (e.g., 1 : 2 : 3). For each method a symmetry theory has been developed that enables the prediction of the direction and sign of vorticity as functions of the field frequencies and phases. However, this theory has its limitations. It only applies to those particular phase angles that give rise to fields whose Lissajous plots, or principal 2-d projections thereof, have a high degree of symmetry. Nor can symmetry theory provide a measure of the magnitude of the torque density induced by the field. In this paper a functional of the multiaxial magnetic field is proposed that not only is consistent with all of the predictions of the symmetry theories, but also quantifies the torque density. This functional can be applied to fields whose Lissajous plots lack symmetry and can thus be used to predict a variety of effects and trends that cannot be predicted from the symmetry theories. These trends include the dependence of the magnitude of the torque density on the various frequency ratios, the unexpected reversal of flow with increasing dc field amplitude for certain symmetry-breaking fields, and the existence of off-axis vorticity for rational triads, such as 1 : 3 : 5, that do not have the symmetry required to analyze by symmetry theory. As a result, experimental data are given that show the degree to which this functional is successful in predicting observed trends.

36 MATERIALS SCIENCE↗

Optimization of magneto-resistive response of ion-irradiated exchange biased films through zigzag arrangement of magnetization

Exchange coupled ferromagnetic-antiferromagnetic Ni{sub 81}Fe{sub 19}/Ir{sub 23}Mn{sub 77} films with a zigzag alignment of magnetization are prepared by local ion irradiation. The anisotropic magneto-resistive behavior of the magnetic thin film structures is correlated to the magnetic structure and modeled. A unique uniaxial field sensitivity along the net magnetization alignment is obtained through the orthogonally modulated and magnetic domain wall stabilized magnetic ground state. Controlling local thin film magnetization distributions and, thus, the overall magnetization response opens unique ways to tailor the magneto-resistive sensitivity of functional magnetic thin film devices.

36 MATERIALS SCIENCE↗

Interaction of elliptically polarised cross-degenerate cnoidal waves in an isotropic gyrotropic medium with spatial dispersion of cubic nonlinearity

Three unusual classes of particular analytical solutions to a system of four nonlinear equations are found for slowly varying complex amplitudes of circularly polarised components of the electric field. The system describes the self-action and interaction of two elliptically polarised plane waves collinearly propagating in an isotropic medium with second-order frequency dispersion and spatial dispersion of cubic nonlinearity. The solutions correspond to self-consistent combinations of two elliptically polarised cnoidal waves whose mutually orthogonal polarisation components vary in accordance with pairwise identical laws during propagation. At the same time, the amplitudes of the component with the same circular polarisation are proportional to two different elliptic Jacobi functions with the same periods. (nonlinear optical phenomena)

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Efficient optimization method for finding minimum energy paths of magnetic transitions

Here, efficient algorithms for the calculation of minimum energy paths of magnetic transitions are implemented within the geodesic nudged elastic band (GNEB) approach. While an objective function is not available for GNEB and a traditional line search can, therefore, not be performed, the use of limited memory Broyden–Fletcher–Goldfarb–Shanno (LBFGS) and conjugate gradient algorithms in conjunction with orthogonal spin optimization (OSO) approach is shown to greatly outperform the previously used velocity projection and dissipative Landau–Lifschitz dynamics optimization methods. The implementation makes use of energy weighted springs for the distribution of the discretization points along the path and this is found to improve performance significantly. The various methods are applied to several test problems using a Heisenberg-type Hamiltonian, extended in some cases to include Dzyaloshinskii–Moriya and exchange interactions beyond nearest neighbours. Minimum energy paths are found for magnetization reversals in a nano-island, collapse of skyrmions in two-dimensional layers and annihilation of a chiral bobber near the surface of a three-dimensional magnet. The LBFGS-OSO method is found to outperform the dynamics based approaches by up to a factor of 8 in some cases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Fast Generation of body conforming grids for 3-D

A fast algorithm was developed for accurately generating boundary conforming, three dimensional, consecutively refined, computational grids applicable to arbitrary axial turbomachinery geometry. The method is based on using a single analytic function to generate two dimensional grids on a number of coaxial axisymmetric surfaces positioned between the hub and the shroud. These grids are of the "O" type and are characterized by quasi-orthogonality, geometric periodicity, and an adequate resolution throughout the flowfield. Due to the built in additional nonorthogonal coordinate stretching and shearing, the grid lines leaving the trailing of the blade end at downstream infinity, thus simplifying the numerical treatment of the three dimensional trailing vortex sheet.

Dulikravich, O.↗

Boundary-fitted coordinate systems for numerical solution of partial differential equations - A review

A comprehensive review of methods of numerically generating curvilinear coordinate systems with coordinate lines coincident with all boundary segments is given. Some general mathematical framework and error analysis common to such coordinate systems is also included. The general categories of generating systems are those based on conformal mapping, orthogonal systems, nearly orthogonal systems, systems produced as the solution of elliptic and hyperbolic partial differential equations, and systems generated algebraically by interpolation among the boundaries. Also covered are the control of coordinate line spacing by functions embedded in the partial differential operators of the generating system and by subsequent stretching transformation. Dynamically adaptive coordinate systems, coupled with the physical solution, and time-dependent systems that follow moving boundaries are treated. References reporting experience using such coordinate systems are reviewed as well as those covering the system development.

Thompson, J. F.↗

Reduced-Order Aerodynamic Modeling Based on CFD Frequency Responses from Multisine Inputs

A system identification analysis was performed to determine a reduced-order model (ROM) of a computational fluid dynamics (CFD) solver in support of linear aeroservoelastic model development and feedback control design. The approach was applied to the FUN3D code for the half-span wind tunnel test article used in the NASA-Boeing collaboration called the Integrated Adaptive Wing Technology Maturation (IAWTM) project. In a transonic flow condition, multiple inputs (11 structural mode displacements and 3 control surface deflections) were simultaneously excited with orthogonal phase-optimized multisines while multiple outputs (the corresponding 14 generalized aerodynamic forces) were recorded. From these recorded times series, the matrix of frequency responses was computed and subsequently fit using rational function approximations (RFAs). It was found that the entire (14 x 14) matrix of frequency responses could be determined from a single CFD run and that results generally followed trends predicted using other methods. Differences were attributed to the modeling fidelity and nonlinearities from structural mode and control surface interactions at higher reduced frequencies. More accurate fits of the RFAs to the frequency response data were obtained by making two CFD runs, one with only structural mode excitations and one with only control surface excitations, which reduced the degree of nonlinearity in the modeling data.

Aeroservoelasticity↗

Kramers nodal line in the charge density wave state of YTe 3 and the influence of twin domains

Recent studies have focused on the relationship between charge density wave (CDW) collective electronic ground states and nontrivial topological states. YTe 3 , a nonmagnetic quasi-two-dimensional chalcogenide, has been reported to exhibit a CDW state below 334 K. Using angle-resolved photoemission spectroscopy (ARPES) and density functional theory (DFT), we establish that YTe 3 is a CDW-induced Kramers nodal line (KNL) metal, a recently proposed topological state of matter. Scanning tunneling microscopy and low energy electron diffraction reveal two orthogonal domains, each with a unidirectional CDW and a similar wave vector (𝐪 CDW ). When the influence of twin domains is considered, the effective band structure (EBS) computations that utilize DFT-calculated bands using a noncentrosymmetric structure determined by x-ray crystallography, show excellent agreement with ARPES. The noncentrosymmetry of YTe 3 is established by Raman spectroscopy. The Fermi surface and ARPES intensity plots show weak shadow bands displaced by 𝐪 CDW from the main bands. Furthermore, these are linked to CDW modulation, as the EBS calculation confirms. Bilayer split main and shadow bands suggest the existence of crossings, according to theory and experiment. DFT bands, including spin-orbit coupling, indicate existence of a KNL along the Σ direction from multiple crossings of bands dispersing perpendicular to it. Additionally, doubly degenerate bands are only found along the KNL at all energies, with some bands dispersing through the Fermi level.

Angle-resolved photoemission spectroscopy↗

Disconjugacy, regularity of multi-indexed rationally extended potentials, and Laguerre exceptional polynomials

The power of the disconjugacy properties of second-order differential equations of Schrödinger type to check the regularity of rationally extended quantum potentials connected with exceptional orthogonal polynomials is illustrated by re-examining the extensions of the isotonic oscillator (or radial oscillator) potential derived in kth-order supersymmetric quantum mechanics or multistep Darboux-Bäcklund transformation method. The function arising in the potential denominator is proved to be a polynomial with a nonvanishing constant term, whose value is calculated by induction over k. The sign of this term being the same as that of the already known highest degree term, the potential denominator has the same sign at both extremities of the definition interval, a property that is shared by the seed eigenfunction used in the potential construction. By virtue of disconjugacy, such a property implies the nodeless character of both the eigenfunction and the resulting potential.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

GRID2D/3D: A computer program for generating grid systems in complex-shaped two- and three-dimensional spatial domains. Part 2: User's manual and program listing

An efficient computer program, called GRID2D/3D, was developed to generate single and composite grid systems within geometrically complex two- and three-dimensional (2- and 3-D) spatial domains that can deform with time. GRID2D/3D generates single grid systems by using algebraic grid generation methods based on transfinite interpolation in which the distribution of grid points within the spatial domain is controlled by stretching functions. All single grid systems generated by GRID2D/3D can have grid lines that are continuous and differentiable everywhere up to the second-order. Also, grid lines can intersect boundaries of the spatial domain orthogonally. GRID2D/3D generates composite grid systems by patching together two or more single grid systems. The patching can be discontinuous or continuous. For continuous composite grid systems, the grid lines are continuous and differentiable everywhere up to the second-order except at interfaces where different single grid systems meet. At interfaces where different single grid systems meet, the grid lines are only differentiable up to the first-order. For 2-D spatial domains, the boundary curves are described by using either cubic or tension spline interpolation. For 3-D spatial domains, the boundary surfaces are described by using either linear Coon's interpolation, bi-hyperbolic spline interpolation, or a new technique referred to as 3-D bi-directional Hermite interpolation. Since grid systems generated by algebraic methods can have grid lines that overlap one another, GRID2D/3D contains a graphics package for evaluating the grid systems generated. With the graphics package, the user can generate grid systems in an interactive manner with the grid generation part of GRID2D/3D. GRID2D/3D is written in FORTRAN 77 and can be run on any IBM PC, XT, or AT compatible computer. In order to use GRID2D/3D on workstations or mainframe computers, some minor modifications must be made in the graphics part of the program; no modifications are needed in the grid generation part of the program. The theory and method used in GRID2D/3D is described.

Bailey, R. T.↗

GRID2D/3D: A computer program for generating grid systems in complex-shaped two- and three-dimensional spatial domains. Part 1: Theory and method

An efficient computer program, called GRID2D/3D was developed to generate single and composite grid systems within geometrically complex two- and three-dimensional (2- and 3-D) spatial domains that can deform with time. GRID2D/3D generates single grid systems by using algebraic grid generation methods based on transfinite interpolation in which the distribution of grid points within the spatial domain is controlled by stretching functions. All single grid systems generated by GRID2D/3D can have grid lines that are continuous and differentiable everywhere up to the second-order. Also, grid lines can intersect boundaries of the spatial domain orthogonally. GRID2D/3D generates composite grid systems by patching together two or more single grid systems. The patching can be discontinuous or continuous. For continuous composite grid systems, the grid lines are continuous and differentiable everywhere up to the second-order except at interfaces where different single grid systems meet. At interfaces where different single grid systems meet, the grid lines are only differentiable up to the first-order. For 2-D spatial domains, the boundary curves are described by using either cubic or tension spline interpolation. For 3-D spatial domains, the boundary surfaces are described by using either linear Coon's interpolation, bi-hyperbolic spline interpolation, or a new technique referred to as 3-D bi-directional Hermite interpolation. Since grid systems generated by algebraic methods can have grid lines that overlap one another, GRID2D/3D contains a graphics package for evaluating the grid systems generated. With the graphics package, the user can generate grid systems in an interactive manner with the grid generation part of GRID2D/3D. GRID2D/3D is written in FORTRAN 77 and can be run on any IBM PC, XT, or AT compatible computer. In order to use GRID2D/3D on workstations or mainframe computers, some minor modifications must be made in the graphics part of the program; no modifications are needed in the grid generation part of the program. This technical memorandum describes the theory and method used in GRID2D/3D.

Shih, T. I.-P.↗

Calipso Observations of the Dependence of Homo- and Heterogeneous Ice Nucleation in Cirrus Clouds on Latitude, Season and Surface Condition

There are two fundamental mechanisms through which cirrus clouds form; homo- and heterogeneous ice nucleation (henceforth hom and het). The relative contribution of each mechanism to ice crystal production often determines the microphysical and radiative properties of a cirrus cloud. A new satellite remote sensing method is described in this study to estimate cirrus cloud ice particle number concentration and the relative contribution of hom and het to cirrus cloud formation as a function of altitude, latitude, season and surface type (e.g. land vs. ocean). This method uses co-located observations from the Infrared Imaging Radiometer (IIR) and from the CALIOP (Cloud and Aerosol Lidar with Orthogonal Polarization) lidar aboard the CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation) polar orbiting satellite, employing IIR channels at 10.6 μm and 12.05 μm. The method is applied here to single-layered clouds of visible optical depth between about 0.3 and 3. Two years of Version 3 data have been analyzed for the years 2008 and 2013, with each season characterized in terms of 532 nm cirrus cloud centroid altitude and temperature, the cirrus cloud ice particle number concentration, effective diameter, layer-average ice water content and visible optical depth. Using a conservative criterion for hom cirrus, on average, the sampled cirrus clouds formed through hom occur about 43% of the time in the Arctic and 50% of the time in the Antarctic, and during winter at mid-latitudes in the Northern Hemisphere, hom cirrus occur 37% of the time. Elsewhere (and during other seasons in the Northern Hemisphere mid-latitudes), this hom cirrus fraction is lower. Processes that could potentially explain these observations are discussed, as well as the potential relevancy of these results to ice nucleation studies, climate modeling and jet-stream dynamics.

Mitchell, David L.↗

Polymers depolymerizable by metathesis of a cleavable unit

Novel polymers are depolymerizable by metathesis of a cleavable unit. As an example, a series of linear and crosslinked polyurethanes were prepared that can be selectively depolymerized under mild conditions. Two unique polyols were synthesized bearing unsaturated units in a configuration designed to favor ring-closing metathesis to five- and six-membered cycloalkenes. These polyols were co-polymerized with toluene diisocyanate to generate linear polyurethanes and trifunctional hexamethylene- and diphenylmethane-based isocyanates to generate crosslinked polyurethanes. The polyol design is such that the ring-closing metathesis reaction cleaves the backbone of the polymer chain. Upon exposure to dilute solutions of Grubbs' catalyst under ambient conditions, the polyurethanes were rapidly depolymerized to low molecular weight, soluble products bearing vinyl and cycloalkene functionalities. These functionalities enabled further re-polymerization by traditional strategies for polymerization of double bonds. This general approach can be expanded to develop a range of chemically recyclable condensation polymers that are readily depolymerized by orthogonal metathesis chemistry.

Jones, Brad Howard↗

Transitional regimes of natural convection in a differentially heated cubical cavity under the effects of wall and molecular gas radiation

The transition to unsteadiness and the dynamics of weakly turbulent natural convection, coupled to wall or gas radiation in a differentially heated cubical cavity with adiabatic lateral walls, are studied numerically. The working fluid is air with small contents of water vapor and carbon dioxide whose infrared spectral radiative properties are modelled by the absorption distribution function model. A pseudo spectral Chebyshev collocation method is used to solve the flow field equations and is coupled to a direct ray tracing method for radiation transport. Flow structures are identified by means of either the proper orthogonal decomposition or the dynamic mode decomposition methods. We first retrieve the classical mechanism of transition to unsteadiness without radiation, characterized by counter-rotating streamwise-oriented vortices generated at the exit of the vertical boundary layers. Wall radiation through a transparent medium leads to a homogenization of lateral wall temperatures and the resulting transition mechanism is similar to that obtained with perfectly conducting lateral walls. The transition is due to an unstable stratification upstream the vertical boundary layers and is characterized by periodically oscillating transverse rolls of axis perpendicular to the main flow. When molecular gas radiation is accounted for, no periodic solution is found and the transition to unsteadiness displays complex structures with chimneys-like rolls whose axes are again parallel to the main flow. The origin of this instability is probably due to centrifugal forces, as suggested previously for the case without radiation. Above the transition to unsteadiness, at Ra = 3 × 10{sup 8}, it is shown that both wall and gas radiation significantly intensify turbulent fluctuations, decrease the thermal stratification in the core of the cavity, and increase the global circulation.

42 ENGINEERING↗

All-electron Kohn–Sham density functional theory on hierarchic finite element spaces

In this work, a real space formulation of the Kohn–Sham equations is developed, making use of the hierarchy of finite element spaces from different polynomial order. The focus is laid on all-electron calculations, having the highest requirement onto the basis set, which must be able to represent the orthogonal eigenfunctions as well as the electrostatic potential. A careful numerical analysis is performed, which points out the numerical intricacies originating from the singularity of the nuclei and the necessity for approximations in the numerical setting, with the ambition to enable solutions within a predefined accuracy. In this context the influence of counter-charges in the Poisson equation, the requirement of a finite domain size, numerical quadratures and the mesh refinement are examined as well as the representation of the electrostatic potential in a high order finite element space. The performance and accuracy of the method is demonstrated in computations on noble gases. In addition the finite element basis proves its flexibility in the calculation of the bond-length as well as the dipole moment of the carbon monoxide molecule.

97 MATHEMATICS AND COMPUTING↗

Device convolution effects on the collective scattering signal of the E Multiplication-Sign B mode from Hall thruster experiments: 2D dispersion relation

The effect of the collective light scattering diagnostic transfer function is considered in the context of the dispersion relation of the unstable E Multiplication-Sign B mode previously reported. This transfer function is found to have a contribution to the measured frequencies and mode amplitudes which is more or less significant depending on the measurement wavenumbers and angles. After deconvolution, the experimental data are found to be possibly compatible with the idea that the mode frequency in the jet frame (after subtraction of the Doppler effect due to the plasma motion along the thruster axis) is independent of the orientation of the wave vector in the plane orthogonal to the local magnetic field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Localization of quantum Bernoulli noises

The family (∂{sub k},∂{sub k}{sup *}){sub k≥0} of annihilation and creation operators acting on square integrable functionals of a Bernoulli process Z= (Z{sub k}){sub k⩾0} can be interpreted as quantum Bernoulli noises. In this note we consider the operator family (ℓ{sub k},ℓ{sub k}{sup *}){sub k≥0}, where ℓ{sub k}=∂{sub k}E{sub k} with E{sub k} being the conditional expectation (operator) given σ-field σ(Z{sub j}; 0 ⩽j⩽k). We show that ℓ{sub k} (resp. ℓ{sub k}{sup *}) is essentially a kind of localization of the annihilation operator ∂{sub k} (resp. creation operator ∂{sub k}{sup *}). We examine properties of the family (ℓ{sub k},ℓ{sub k}{sup *}){sub k≥0} and prove, among other things, that ℓ{sub k} and ℓ{sub k}{sup *} satisfy a local canonical anti-communication relation and (ℓ{sub k}{sup *}){sub k≥0} forms a mutually orthogonal operator sequence although each ℓ{sub k} is not a projection operator. We find that the operator series Σ{sub k=0}{sup ∞}ℓ{sub k}{sup *}Xℓ{sub k} converges in the strong operator topology for each bounded operator X acting on square integrable functionals of Z. In particular we get an explicit sum of the operator series Σ{sub k=0}{sup ∞}ℓ{sub k}{sup *}ℓ{sub k}. A useful norm estimate on Σ{sub k=0}{sup ∞}ℓ{sub k}{sup *}Xℓ{sub k} is also obtained. Finally we show applications of our main results to quantum dynamical semigroups and quantum probability.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Bayesian estimation of Karhunen–Loève expansions; A random subspace approach

One of the most widely-used statistical procedures for dimensionality reduction of high dimensional random fields is Principal Component Analysis (PCA), which is based on the Karhunen-Lo eve expansion (KLE) of a stochastic process with finite variance. The KLE is analogous to a Fourier series expansion for a random process, where the goal is to find an orthogonal transformation for the data such that the projection of the data onto this orthogonal subspace is optimal in the L 2 sense, i.e, which minimizes the mean square error. In practice, this orthogonal transformation is determined by performing an SVD (Singular Value Decomposition) on the sample covariance matrix or on the data matrix itself. Sampling error is typically ignored when quantifying the principal components, or, equivalently, basis functions of the KLE. Furthermore, it is exacerbated when the sample size is much smaller than the dimension of the random field. In this paper, we introduce a Bayesian KLE procedure, allowing one to obtain a probabilistic model on the principal components, which can account for inaccuracies due to limited sample size. The probabilistic model is built via Bayesian inference, from which the posterior becomes the matrix Bingham density over the space of orthonormal matrices. We use a modified Gibbs sampling procedure to sample on this space and then build a probabilistic Karhunen-Lo eve expansions over random subspaces to obtain a set of low-dimensional surrogates of the stochastic process. We illustrate this probabilistic procedure with a finite dimensional stochastic process inspired by Brownian motion.

97 MATHEMATICS AND COMPUTING↗