Dual Functions and Heat Expansions
Radial heat polynomial transformations in series expansions of heat equation
HEAT EQUATION↗
SEARCH · Engineering Papers
Results for “SERIES EXPANSION”
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.
Fast and accurate implementation of Fourier spectral approximations of nonlocal diffusion operators and its applications
This work is concerned with the Fourier spectral approximation of various integral differential equations associated with some linear nonlocal diffusion and peridynamic operators under periodic boundary conditions. For radially symmetric kernels, the nonlocal operators under consideration are diagonalizable in the Fourier space so that the main computational challenge is on the accurate and fast evaluation of their eigenvalues or Fourier symbols consisting of possibly singular and highly oscillatory integrals. For a large class of fractional power-like kernels, we propose a new approach based on reformulating the Fourier symbols both as coefficients of a series expansion and solutions of some simple ODE models. We then propose a hybrid algorithm that utilizes both truncated series expansions and high order Runge–Kutta ODE solvers to provide fast evaluation of Fourier symbols in both one and higher dimensional spaces. It is shown that this hybrid algorithm is robust, efficient and accurate. As applications, we combine this hybrid spectral discretization in the spatial variables and the fourth-order exponential time differencing Runge–Kutta for temporal discretization to offer high order approximations of some nonlocal gradient dynamics including nonlocal Allen–Cahn equations, nonlocal Cahn–Hilliard equations, and nonlocal phase-field crystal models. Numerical results show the accuracy and effectiveness of the fully discrete scheme and illustrate some interesting phenomena associated with the nonlocal models.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
General-altitude transformations between geocentric and geodetic coordinates
Formulas for the general-altitude (height above the ellipsoid) transformation from geocentric to geodetic coordinates and vice versa are derived. The set of four formulas is expressed in each of two useful forms: series expansions in powers of the earth's flattening and series expansions in powers of the earth's eccentricity. The error incurred in these expansions is of the order of one part in 30 million.-
Long, S. A. T.↗
Scalar and tensor spherical harmonics expansion of the velocity correlation in homogeneous anisotropic turbulence
The representation theory of the rotation group is applied to construct a series expansion of the correlation tensor in homogeneous anisotropic turbulence. The resolution of angular dependence is the main analytical difficulty posed by anisotropic turbulence; representation theory parametrises this dependence by a tensor analogue of the standard spherical harmonics expansion of a scalar. As a result, the series expansion is formulated in terms of explicitly constructed tensor bases with scalar coefficients determined by angular moments of the correlation tensor.
97 MATHEMATICS AND COMPUTING↗
Disentangling the f(R)-duality
Motivated by UV realisations of Starobinsky-like inflation models, we study generic exponential plateau-like potentials to understand whether an exact f(R)-formulation may still be obtained when the asymptotic shift-symmetry of the potential is broken for larger field values. Potentials which break the shift symmetry with rising exponentials at large field values only allow for corresponding f(R)-descriptions with a leading order term R{sup n} with 1<n<2, regardless of whether the duality is exact or approximate. The R{sup 2}-term survives as part of a series expansion of the function f(R) and thus cannot maintain a plateau for all field values. We further find a lean and instructive way to obtain a function f(R) describing m{sup 2}ϕ{sup 2}-inflation which breaks the shift symmetry with a monomial, and corresponds to effectively logarithmic corrections to an R+R{sup 2} model. These examples emphasise that higher order terms in f(R)-theory may not be neglected if they are present at all. Additionally, we relate the function f(R) corresponding to chaotic inflation to a more general Jordan frame set-up. In addition, we consider f(R)-duals of two given UV examples, both from supergravity and string theory. Finally, we outline the CMB phenomenology of these models which show effects of power suppression at low-ℓ.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
An Inverse Factorial Series for a General Gamma Ratio and Related Properties of the Nørlund–Bernoulli Polynomials
The inverse factorial series expansion for the ratio of products of gamma functions whose arguments are linear functions of the variable is found. A recurrence relation for the coefficients in terms of the Nørlund–Bernoulli polynomials is provided, and the half-plane of convergence is determined. The results obtained naturally supplement a number of previous investigations of the gamma ratios, which began in the 1930-ies. The expansion obtained in this paper plays a crucial role in the study of the behavior of the delta-neutral Fox’s H-function in the neighborhood of its finite singular point. A particular case of the inverse factorial series expansion is used in deriving a possibly new identity for the Nørlund–Bernoulli polynomials.
97 MATHEMATICS AND COMPUTING↗
Explicit expressions of quantum mechanical rotation operators for spins 1 to 2
Quantum mechanical rotation operators are the subject of quantum mechanics, mathematics and pulsed magnetic resonance spectroscopies, namely NMR, EPR and ENDOR. They are also necessary for spin based quantum information systems. The rotation operators of spin 1/2 are well known and can be found in related textbooks. But rotation operators of other spins greater than 1/2 can be found numerically by evaluating the series expansions of exponential operator obtained from Schrödinger equation, or by evaluating Wigner-d formula or by evaluating recently established expressions in polynomial forms discussed in the text. In this work, explicit symbolic expressions of x, y and z components of rotation operators for spins 1 to 2 are worked out by evaluating series expansion of exponential operator for each element of operators and utilizing linear curve fitting process. The procedures gave out exact expressions of each element of the rotation operators. The operators of spins greater than 2 are under study and will be published in a separate paper.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Transport and optical conductivity in the Hubbard model: A high-temperature expansion perspective
We derive analytical expressions for the spectral moments of the dynamical response functions of the Hubbard model using the high-temperature series expansion. We consider generic dimension d as well as the infinite- d limit, arbitrary electron density n, and both finite and infinite repulsion U. We use moment-reconstruction methods to obtain the one-electron spectral function, the self-energy, and the optical conductivity. They are all smooth functions at high temperature and, at large U, they are featureless with characteristic widths of the order of the lattice hopping parameter t. In the infinite-d limit, we compare the series expansion results with accurate numerical renormalization group and interaction expansion quantum Monte Carlo results. We find excellent agreement down to surprisingly low temperatures, throughout most of the bad-metal regime, which applies for T≳ (1 - n)D, the Brinkman-Rice scale. The resistivity increases linearly in T at high temperature without saturation. This results from the 1/T behavior of the compressibility or kinetic energy, which play the role of the effective carrier number. In contrast, the scattering time (or diffusion constant) saturates at high T. We find that σ( n,T) ≈ (1 - n)σ(n = 0,T) to a very good approximation for all n, with σ(n = 0,T) ∝ t/T at high temperatures. The saturation at small n occurs due to a compensation between the density dependence of the effective number of carriers and that of the scattering time. The T dependence of the resistivity displays a kneelike feature which signals a crossover to the intermediate-temperature regime where the diffusion constant (or scattering time) starts increasing with decreasing T. At high temperatures, the thermopower obeys the Heikes formula, while the Wiedemann-Franz law is violated with the Lorenz number vanishing as 1/T 2 . The relevance of our calculations to experiments probing high-temperature transport in materials with strong electronic correlations or ultracold atomic gases in an optical lattice is briefly discussed.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Radar altimeter mean return waveforms from near-normal-incidence ocean surface scattering
Under assumptions common in radar altimetry, the mean backscattered return power for a short-pulse radar and near-normal-incidence scattering from a rough ocean surface is given by the convolution of several terms. For a nearly Gaussian transmitted pulse shape scattered from a nearly Gaussian distributed sea surface, a small-argument series expansion of one of the terms within the convolution leads to a several-term power series expansion for the mean return waveform. Specific expressions are given for the first four terms. These results, which require much less computer time than would the otherwise necessary numerical convolution, are useful for data analysis from current or past radar altimeters and for design studies of future systems. Several representative results are presented for an idealized Seasat radar altimeter.
Hayne, G. S.↗
Structural Reliability Analysis and Optimization: Use of Approximations
This report is intended for the demonstration of function approximation concepts and their applicability in reliability analysis and design. Particularly, approximations in the calculation of the safety index, failure probability and structural optimization (modification of design variables) are developed. With this scope in mind, extensive details on probability theory are avoided. Definitions relevant to the stated objectives have been taken from standard text books. The idea of function approximations is to minimize the repetitive use of computationally intensive calculations by replacing them with simpler closed-form equations, which could be nonlinear. Typically, the approximations provide good accuracy around the points where they are constructed, and they need to be periodically updated to extend their utility. There are approximations in calculating the failure probability of a limit state function. The first one, which is most commonly discussed, is how the limit state is approximated at the design point. Most of the time this could be a first-order Taylor series expansion, also known as the First Order Reliability Method (FORM), or a second-order Taylor series expansion (paraboloid), also known as the Second Order Reliability Method (SORM). From the computational procedure point of view, this step comes after the design point identification; however, the order of approximation for the probability of failure calculation is discussed first, and it is denoted by either FORM or SORM. The other approximation of interest is how the design point, or the most probable failure point (MPP), is identified. For iteratively finding this point, again the limit state is approximated. The accuracy and efficiency of the approximations make the search process quite practical for analysis intensive approaches such as the finite element methods; therefore, the crux of this research is to develop excellent approximations for MPP identification and also different approximations including the higher-order reliability methods (HORM) for representing the failure surface. This report is divided into several parts to emphasize different segments of the structural reliability analysis and design. Broadly, it consists of mathematical foundations, methods and applications. Chapter I discusses the fundamental definitions of the probability theory, which are mostly available in standard text books. Probability density function descriptions relevant to this work are addressed. In Chapter 2, the concept and utility of function approximation are discussed for a general application in engineering analysis. Various forms of function representations and the latest developments in nonlinear adaptive approximations are presented with comparison studies. Research work accomplished in reliability analysis is presented in Chapter 3. First, the definition of safety index and most probable point of failure are introduced. Efficient ways of computing the safety index with a fewer number of iterations is emphasized. In chapter 4, the probability of failure prediction is presented using first-order, second-order and higher-order methods. System reliability methods are discussed in chapter 5. Chapter 6 presents optimization techniques for the modification and redistribution of structural sizes for improving the structural reliability. The report also contains several appendices on probability parameters.
Grandhi, Ramana V.↗
Time-dependent mixed convection heat transfer from a sphere in a micro-gravity environment
A fundamental problem of interest for crystal growth in micro-gravity applications involves the mixed convection heat transfer from a sphere in a uniform flow of fluid at a differing temperature. Under the combined influence of the imposed free stream as well as an induced buoyancy force due to thermal expansion of the fluid, the heat transfer from the sphere will be different from that of either the pure forced convection flow or the pure free convection flow. For the present study, the method of matched asymptotic expansions is applied to the laminar flow problem of an impulsively heated, impulsively started sphere in an originally quiescent fluid. Time series expansions are developed for the dependent variables by acknowledging the existence of two distinct regions: one, an inner region, near the sphere, in which viscous effects are significant; and two, an outer region in which the fluid may be treated as inviscid. The time series expansions are developed in terms of the Reynolds number and Richardson number (Buoyancy Parameter), and the relevant heat transfer and drag coefficients are calculated and plotted.
Hommel, Mark J.↗
Gravitational waves from rotating and precessing rigid bodies. 2: General solutions and computationally useful formulae
The classical mechanics results for free precession which are needed in order to calculate the weak field, slow-motion, quadrupole-moment gravitational waves are reviewed. Within that formalism, algorithms are given for computing the exact gravitational power radiated and waveforms produced by arbitrary rigid-body freely-precessing sources. The dominant terms are presented in series expansions of the waveforms for the case of an almost spherical object precessing with a small wobble angle. These series expansions, which retain the precise frequency dependence of the waves, may be useful for gravitational astronomers when freely-precessing sources begin to be observed.
Zimmerman, M.↗
Gravitational waves from rotating and precessing rigid bodies. II - General solutions and computationally useful formulas
A rigid, freely precessing Newtonian body emits gravitational radiation. In this paper the classical-mechanics results for free precession which are needed in order to calculate the weak-field slow-motion quadrupole-moment gravitational waves are reviewed. Within that formalism, algorithms for computing the exact gravitational power radiated and waveforms produced by arbitrary rigid-body freely precessing source are given. Also presented are the dominant terms in series expansions of the waveforms for the case of an almost-spherical object precessing with a small wobble angle. These series expansions, which retain the precise frequency dependence of the waves, may be useful for gravitational astronomers when freely precessing sources begin to be observed.
Zimmermann, M.↗
On the spin separation of algebraic two-component relativistic Hamiltonians
The separation of the spin-free and spin-dependent terms of a given relativistic Hamiltonian is usually facilitated by the Dirac identity. However, this is no longer possible for the recently developed exact two-component relativistic Hamiltonians derived from the matrix representation of the Dirac equation in a kinetically balanced basis. This stems from the fact that the decoupling matrix does not have an explicit form. To resolve this formal difficulty, we first define the spin-dependent term as the difference between a two-component Hamiltonian corresponding to the full Dirac equation and its one-component counterpart corresponding to the spin-free Dirac equation. The series expansion of the spin-dependent term is then developed in two different ways. One is in the spirit of the Douglas-Kroll-Hess (DKH) transformation and the other is based on the perturbative expansion of a two-component Hamiltonian of fixed structure, either the two-step Barysz-Sadlej-Snijders (BSS) or the one-step exact two-component (X2C) form. The algorithms for constructing arbitrary order terms are proposed for both schemes and their convergence patterns are assessed numerically. Truncating the expansions to finite orders leads naturally to a sequence of novel spin-dependent Hamiltonians. In particular, the order-by-order distinctions among the DKH, BSS, and X2C approaches can nicely be revealed. The well-known Pauli, zeroth-order regular approximation, and DKH1 spin-dependent Hamiltonians can also be recovered naturally by appropriately approximating the decoupling and renormalization matrices. On the practical side, the sf-X2C+so-DKH3 Hamiltonian, together with appropriately constructed generally contracted basis sets, is most promising for accounting for relativistic effects in two steps, first spin-free and then spin-dependent, with the latter applied either perturbatively or variationally.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Computation of magnetic fields from field components on a plane grid
An algorithm is presented to calculate the field components of the magnetic field [B x (x, y, z), B y (x, y, z), B z (x, y, z)] at a point (x, y, z) in space, from the knowledge of the components [B x (x, y) = 0, z), B y (x, y = 0, z), B z (x, y = 0,z)] on a “reference plane”, which is normal to the y-axis at . The algorithm, which is a general one and is not restricted to fields with mid-plane symmetry is based on the Maclaurin series expansion of the magnetic field components at any point in space in terms of the distance (y) of the point from the reference plane. The coefficients of the Maclaurin series expansion are expressed in terms of the on-plane field components and their partial derivatives with respect to spatial coordinates (). The field components are usually generated from magnetic field measurements on a rectangular grid on the plane. This algorithm was employed in 1986 in the RAYTRACE computer code to help calculate the optical properties of magnets and of the Alternating Gradient Synchrotron (AGS) at the Brookhaven National Laboratory (BNL). Here, a general mathematical formulation of this algorithm based on the differential algebraic method was presented by Makino in 2011. This paper presents the step by step derivation of the algorithm and provides the necessary formulas to be introduced by the reader in any computer code which requires the field components generated by magnetic devices. In addition provides an example of the use of the algorithm and its limitations as applied to a Halbach type magnet with or without median plane symmetry.
43 PARTICLE ACCELERATORS↗
Algorithm to calculate off-plane magnetic field from an on-plane field map
An algorithm is presented to calculate the field components [B x (x,y,z), B y (x,y,z), B z (x,y,z)] of the magnetic field at a point (x,y,z) in space, from the knowledge of the components [B x (x,y = 0,z), B y (x,y = 0,z), B z (x,y = 0,z)] on a ”reference plane”, which is normal to the y-axis at y=0. The algorithm, which is a general one and is not restricted to fields with mid-plane symmetry is based on the Taylor series expansion of the magnetic field components at any point in space in terms of the distance (y) of the point from the reference plane. The coefficients of the Taylor series expansion are expressed in terms of the on-plane field components and their partial derivatives with respect to spatial coordinates (x,z). The field components are usually generated from magnetic field measurements on a rectangular grid on the plane. The required partial derivatives of the field components on the plane can be computed with two methods; numerically, or by fitting a polynomial function to the experimentally measured values of the field components and subsequently take the required partial derivatives of the polynomial function. Each of these two method which calculates the partial derivatives of the field components on the plane, is applied to calculate the off-plane field components of the magnetic field generated by a modified Halbach magnet [1] of inner radius of 4.4 cm and it is found that each method provides accurate results with a relative error <1% deviation from the ”ideal” fields if the y distance of the point from the reference plane is smaller than 2 cm. This algorithm is part of the RAYTRACE computer code [2] and has been applied [3] on a dipole magnet with median plane symmetry.
43 PARTICLE ACCELERATORS↗
The symbolic computation of series solutions to ordinary differential equations using trees (extended abstract)
Algorithms previously developed by the author give formulas which can be used for the efficient symbolic computation of series expansions to solutions of nonlinear systems of ordinary differential equations. As a by product of this analysis, formulas are derived which relate to trees to the coefficients of the series expansions, similar to the work of Leroux and Viennot, and Lamnabhi, Leroux and Viennot.
Grossman, Robert↗