Engineering PapersSearch

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.

At least 19 records

Two-point Taylor series expansions

Coefficients calculated for Taylor series expansion about two points - application of Taylor expansion to two-body problem

TWO-BODY PROBLEM

A series expansion of the acoustic power radiated from planar sources

A series expansion in ascending powers of the wavenumber k is derived for the acoustic power delivered by baffled or unbaffled planar sources. This series provides a relatively simple means of derving expressions for the power radiated by a baffled source with a known velocity distribution and can be used for unbaffled plates when the velocity field outside the plate is also known. The terms in the series are calculated from the moments of this velocity distribution in the plane containing the source. If these moments are written as derivaties in wavenumber space, it is shown that a MacLaurin expansion of the Fourier transformed velocity provides an easy technique for computing the first few terms of the acoustic power. Examples are provided for baffled, rectangular plates with various boundary conditions. The arbirarily shaped plate with free boundaries is particularly interesting. It is proven that the volume flow across it surface must be zero and as a result corner and edge mode radiation cannot exist for this kind of source.

Willams, E. G.

Loop series expansions for tensor networks

Belief propagation (BP) can be a useful tool to approximately contract a tensor network, provided that the contributions from any closed loops in the network are sufficiently weak. In this article, we describe how a loop series expansion can be applied to systematically improve the accuracy of a BP approximation to a tensor network contraction, in principle converging arbitrarily close to the exact result. More generally, our result provides a framework for expanding a tensor network as a sum of component networks in a hierarchy of increasing complexity. We benchmark this proposal for the contraction of infinite projected entangled pair states, either representing the ground state of an Affleck-Kennedy-Lieb-Tasaki (AKLT) model or with randomly defined tensors, where it is shown to improve in accuracy over standard BP by several orders of magnitude while incurring only a minor increase in computational cost. These results indicate that the proposed series expansions could be a useful tool to accurately evaluate tensor networks in cases that otherwise exceed the limits of established contraction routines.

Evenbly, Glen [AWS Center for Quantum Computing, P

A note on velocity-related series expansions in the two-body problem

The present note describes a few important series expansions in the two-body problem. They are related to the magnitude v of the velocity vector and are important for the treatment of atmospheric drag by the method of general perturbations. These series have been obtained with computerized Poisson series manipulations. The results are given to order seven in the eccentricity, for both the mean anomaly and the true anomaly.

Broucke, R.

Power series expansions for the frequency and period of the limit cycle of the van der Pol equation

An equation reported by van der Pol (1926) in connection with relaxation-oscillations studies is considered. The equation contains the factor epsilon which can assume values in the range from zero to infinity. The period T(epsilon), or equivalently the frequency nu(epsilon) of the limit cycle has been studied. However, to date there has been little success in discovering the analytical structure of T(epsilon) as a function of epsilon. The present investigation has the objectives to present the Taylor series expansion of nu(epsilon), to locate the singularities which determine the radius of convergence of that expansion, to introduce a new damping variable in terms of which the expansion converges for all epsilon, to form a new expansion for the period T(epsilon) which improves the rate of convergence of the series, to attempt to 'complete' the series, and to compare the obtained results with the numerically determined values of T(epsilon) and with the asymptotic approximation valid for large epsilon.

Andersen, C. M.

A Highly Accurate Technique for the Treatment of Flow Equations at the Polar Axis in Cylindrical Coordinates using Series Expansions

Numerical methods for solving the flow equations in cylindrical or spherical coordinates should be able to capture the behavior of the exact solution near the regions where the particular form of the governing equations is singular. In this work we focus on the treatment of these numerical singularities for finite-differences methods by reinterpreting the regularity conditions developed in the context of pseudo-spectral methods. A generally applicable numerical method for treating the singularities present at the polar axis, when nonaxisymmetric flows are solved in cylindrical, coordinates using highly accurate finite differences schemes (e.g., Pade schemes) on non-staggered grids, is presented. Governing equations for the flow at the polar axis are derived using series expansions near r=0. The only information needed to calculate the coefficients in these equations are the values of the flow variables and their radial derivatives at the previous iteration (or time) level. These derivatives, which are multi-valued at the polar axis, are calculated without dropping the accuracy of the numerical method using a mapping of the flow domain from (0,R)*(0,2pi) to (-R,R)*(0,pi), where R is the radius of the computational domain. This allows the radial derivatives to be evaluated using high-order differencing schemes (e.g., compact schemes) at points located on the polar axis. The proposed technique is illustrated by results from simulations of laminar-forced jets and turbulent compressible jets using large eddy simulation (LES) methods. In term of the general robustness of the numerical method and smoothness of the solution close to the polar axis, the present results compare very favorably to similar calculations in which the equations are solved in Cartesian coordinates at the polar axis, or in which the singularity is removed by employing a staggered mesh in the radial direction without a mesh point at r=0, following the method proposed recently by Mohseni and Colonius (1). Extension of the method described here for incompressible flows or for any other set of equations that are solved on a non-staggered mesh in cylindrical or spherical coordinates with finite-differences schemes of various level of accuracy is immediate.

Constantinescu, George S.

A New Method for Accurate Treatment of Flow Equations in Cylindrical Coordinates Using Series Expansions

The motivation of this work is the ongoing effort at the Center for Turbulence Research (CTR) to use large eddy simulation (LES) techniques to calculate the noise radiated by jet engines. The focus on engine exhaust noise reduction is motivated by the fact that a significant reduction has been achieved over the last decade on the other main sources of acoustic emissions of jet engines, such as the fan and turbomachinery noise, which gives increased priority to jet noise. To be able to propose methods to reduce the jet noise based on results of numerical simulations, one first has to be able to accurately predict the spatio-temporal distribution of the noise sources in the jet. Though a great deal of understanding of the fundamental turbulence mechanisms in high-speed jets was obtained from direct numerical simulations (DNS) at low Reynolds numbers, LES seems to be the only realistic available tool to obtain the necessary near-field information that is required to estimate the acoustic radiation of the turbulent compressible engine exhaust jets. The quality of jet-noise predictions is determined by the accuracy of the numerical method that has to capture the wide range of pressure fluctuations associated with the turbulence in the jet and with the resulting radiated noise, and by the boundary condition treatment and the quality of the mesh. Higher Reynolds numbers and coarser grids put in turn a higher burden on the robustness and accuracy of the numerical method used in this kind of jet LES simulations. As these calculations are often done in cylindrical coordinates, one of the most important requirements for the numerical method is to provide a flow solution that is not contaminated by numerical artifacts. The coordinate singularity is known to be a source of such artifacts. In the present work we use 6th order Pade schemes in the non-periodic directions to discretize the full compressible flow equations. It turns out that the quality of jet-noise predictions using these schemes is especially sensitive to the type of equation treatment at the singularity axis. The objective of this work is to develop a generally applicable numerical method for treating the singularities present at the polar axis, which is particularly suitable for highly accurate finite-differences schemes (e.g., Pade schemes) on non-staggered grids. The main idea is to reinterpret the regularity conditions developed in the context of pseudo-spectral methods. A set of exact equations at the singularity axis is derived using the appropriate series expansions for the variables in the original set of equations. The present treatment of the equations preserves the same level of accuracy as for the interior scheme. We also want to point out the wider utility of the method, proposed here in the context of compressible flow equations, as its extension for incompressible flows or for any other set of equations that are solved on a non-staggered mesh in cylindrical coordinates with finite-differences schemes of various level of accuracy is straightforward. The robustness and accuracy of the proposed technique is assessed by comparing results from simulations of laminar forced-jets and turbulent compressible jets using LES with similar calculations in which the equations are solved in Cartesian coordinates at the polar axis, or in which the singularity is removed by employing a staggered mesh in the radial direction without a mesh point at r = 0.

Constantinescu, G.S.

Determination of Vortex Paths by Series Expansion Technique with Application to Cruciform Wings

A series method of determining two-dimensional vortex paths is considered and applied to the computation of vortex positions behind a slender equal-span cruciform wing at any angle of bank as a function of the distance behind the trailing edge. Calculated paths are shown for four bank angles. For a bank angle of 45 degrees comparison is made with the results of a closed expression given in NACA-TN-2605. For other bank angles water-tank experiments provide qualitative comparison. Satisfactory agreement is found for a sufficient distance downstream to include most practical missile-tail positions. The interference forces on an equal-span cruciform wing are calculated for five angles of bank (including the trivial case of zero bank) from the vortex positions found by use of the series.

Alksne, Alberta Y

A convergent series expansion for hyperbolic systems of conservation laws

The discontinuities piecewise analytic initial value problem for a wide class of conservation laws is considered which includes the full three-dimensional Euler equations. The initial interaction at an arbitrary curved surface is resolved in time by a convergent series. Among other features the solution exhibits shock, contact, and expansion waves as well as sound waves propagating on characteristic surfaces. The expansion waves correspond to he one-dimensional rarefactions but have a more complicated structure. The sound waves are generated in place of zero strength shocks, and they are caused by mismatches in derivatives.

Harabetian, E.

Expansion series of integral functions occurring in unsteady aerodynamics

Attention is given to two real integral functions which occur in the kernel of singular integral equations for subsonic unsteady lifting surfaces. The arguments k, r, and X of the functions correspond to the reduced frequency, spanwise distance, and modified coordinate in the flow direction, respectively. The value of the parameter nu in the functions depends on geometrical conditions. The considered investigation has the objective to present a series for general values of the nonnegative integer nu in order to compute efficiently the integral functions. The approach makes it possible to avoid any approximation or numerical quadrature.

Ueda, T.