SEARCH · Search NASA

Results for “Pade approximations”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Programs Help Create And Evaluate Markov Models

Pade Approximation With Scaling (PAWS) and Scaled Taylor Exponential Matrix (STEM) computer programs provide flexible, user-friendly, language-based interface for creation and evaluation of Markov models describing behaviors of fault-tolerant reconfigurable computer systems. Produce exact solution for probabilities of system failures and provide conservative estimates of numbers of significant digits in solutions. Also offer as part of bundled package with SURE and ASSIST, two other reliable analysis programs developed by Systems Validation Methods group at Langley Research Center.

Butler, Ricky W.↗

NASA NTRS1992

Renormalization group methods for the Reynolds stress transport equations

The Yakhot-Orszag renormalization group is used to analyze the pressure gradient-velocity correlation and return to isotropy terms in the Reynolds stress transport equations. The perturbation series for the relevant correlations, evaluated to lowest order in the epsilon-expansion of the Yakhot-Orszag theory, are infinite series in tensor product powers of the mean velocity gradient and its transpose. Formal lowest order Pade approximations to the sums of these series produce a rapid pressure strain model of the form proposed by Launder, Reece, and Rodi, and a return to isotropy model of the form proposed by Rotta. In both cases, the model constants are computed theoretically. The predicted Reynolds stress ratios in simple shear flows are evaluated and compared with experimental data. The possibility is discussed of deriving higher order nonlinear models by approximating the sums more accurately. The Yakhot-Orszag renormalization group provides a systematic procedure for deriving turbulence models. Typical applications have included theoretical derivation of the universal constants of isotropic turbulence theory, such as the Kolmogorov constant, and derivation of two equation models, again with theoretically computed constants and low Reynolds number forms of the equations. Recent work has applied this formalism to Reynolds stress modeling, previously in the form of a nonlinear eddy viscosity representation of the Reynolds stresses, which can be used to model the simplest normal stress effects. The present work attempts to apply the Yakhot-Orszag formalism to Reynolds stress transport modeling.

Rubinstein, R.↗

NASA NTRS1991

Optimal control study for the Space Station Solar Dynamic power module

The authors present the design of an optimal control system for the Space Station Freedom's Solar Dynamic Fine Pointing and Tracking (SDFPT) module. A very large state model of six rigid body modes and 272 flexible modes is used in conjunction with classical LQG optimal control to produce a full-order controller which satisfies the requirements. The results obtained are compared with those of a classically designed PID (proportional plus integral plus derivative) controller that was implemented for a six-rigid-body-mode forty-flexible-mode model. A major difficulty with designing LQG controllers for large models is solving the Riccati equation that arises from the optimal formulation. A Riccati solver based on a Pade approximation to the matrix sign function is used. A symmetric version of this algorithm is derived for the special class of Hamiltonion matrices, thereby yielding, for large problems, a nearly twofold speed increase over a previous algorithm.

Papadopoulos, P. M.↗

NASA NTRS1994

Accurate finite difference methods for time-harmonic wave propagation

Finite difference methods for solving problems of time-harmonic acoustics are developed and analyzed. Multidimensional inhomogeneous problems with variable, possibly discontinuous, coefficients are considered, accounting for the effects of employing nonuniform grids. A weighted-average representation is less sensitive to transition in wave resolution (due to variable wave numbers or nonuniform grids) than the standard pointwise representation. Further enhancement in method performance is obtained by basing the stencils on generalizations of Pade approximation, or generalized definitions of the derivative, reducing spurious dispersion, anisotropy and reflection, and by improving the representation of source terms. The resulting schemes have fourth-order accurate local truncation error on uniform grids and third order in the nonuniform case. Guidelines for discretization pertaining to grid orientation and resolution are presented.

Harari, Isaac↗

NASA NTRS1994

Second-order reconstruction of the inflationary potential

To first order in the deviation from scale invariance the inflationary potential and its first two derivatives can be expressed in terms of the spectral indices of the scalar and tensor perturbations, n and n(sub T), and their contributions to the variance of the quadrupole CBR temperature anisotropy, S and T. In addition, there is a 'consistency relation' between these quantities: n(sub T) = (-1/ 7)(T/S). We derive the second-order expressions for the inflationary potential and its first two derivatives and the first-order expression for its third derivative, in terms, of n, n(sub T), S, T, and dn/d ln gamma. We also obtain the second-order consistency relation, n(sub T) = (-1/7)(T/S)(1 + 0.11(T/S) + 0.15(n-1)). As an example we consider the exponential potential, the only known case where exact analytic solutions for the perturbation spectra exist. We reconstruct the potential via Taylor expansion (with coefficients calculated at both first and second order), and introduce the Pade approximate as a greatly improved alternative.

Liddle, Andrew R.↗

NASA NTRS1997

Application of AWE Along with a Combined FEM/MoM Technique to Compute RCS of a Cavity-Backed Aperture in an Infinite Ground Plane Over a Frequency Range

A hybrid Finite Element Method (FEM)/Method of Moments (MoM) technique in conjunction with the Asymptotic Waveform Evaluation (AWE) technique is applied to obtain radar cross section (RCS) of a cavity-backed aperture in an infinite ground plane over a frequency range. The hybrid FEM/MoM technique when applied to the cavity-backed aperture results in an integro-differential equation with electric field as the unknown variable, the electric field obtained from the solution of the integro-differential equation is expanded in Taylor series. The coefficients of the Taylor series are obtained using the frequency derivatives of the integro-differential equation formed by the hybrid FEM/MoM technique. The series is then matched via the Pade approximation to a rational polynomial, which can be used to extrapolate the electric field over a frequency range. The RCS of the cavity-backed aperture is calculated using the electric field at different frequencies. Numerical results for a rectangular cavity, a circular cavity, and a material filled cavity are presented over a frequency range. Good agreement between AWE and the exact solution over the frequency range is obtained.

Reddy, C.J.↗

NASA NTRS1999

Phenomenological Theory of the Translational Relaxation Times in Gases

The exact solution to the classical equations governing the translational dispersion and absorption of sound in a gas obscures its relaxational character because of its mathematical complexity. The approach taken here is to solve the secular equation by the method of Pade approximants, which even to the relatively low order R(sub 11) yields a remarkably close approximation to the exact solution over a wide range of frequency/pressure (f/P) ratios. As a result, translational relaxation can be formulated in terms of a conventional relaxation process with well-defined relaxation times, relaxation strength, collision numbers, additivity relations, etc. To extend the theory to high values of f/P ratio, a model is proposed to account for the noncontinuum behavior of the transport coefficients (viscosity and thermal conductivity) as the molecular mean free path approaches the acoustical enclosure dimensions. The theoretical dispersion and absorption show good agreement with measurements in argon over the classical and transition regions of f/P, but a discrepancy appears at higher values of f/P, where collective propagating modes, assumed in the theory, give way to single-particle modes, prevailing in the experiments.

Zuckerwar, Allan J.↗

NASA NTRS2000

TAILSIM Users Guide

The TAILSIM program uses a 4th order Runge-Kutta method to integrate the standard aircraft equations-of-motion (EOM). The EOM determine three translational and three rotational accelerations about the aircraft's body axis reference system. The forces and moments that drive the EOM are determined from aerodynamic coefficients, dynamic derivatives, and control inputs. Values for these terms are determined from linear interpolation of tables that are a function of parameters such as angle-of-attack and surface deflections. Buildup equations combine these terms and dimensionalize them to generate the driving total forces and moments. Features that make TAILSIM applicable to studies of tailplane stall include modeling of the reversible control System, modeling of the pilot performing a load factor and/or airspeed command task, and modeling of vertical gusts. The reversible control system dynamics can be described as two hinged masses connected by a spring. resulting in a fifth order system. The pilot model is a standard form of lead-lag with a time delay applied to an integrated pitch rate and/or airspeed error feedback. The time delay is implemented by a Pade approximation, while the commanded pitch rate is determined by a commanded load factor. Vertical gust inputs include a single 1-cosine gust and a continuous NASA Dryden gust model. These dynamic models. coupled with the use of a nonlinear database, allow the tailplane stall characteristics, elevator response, and resulting aircraft response, to be modeled. A useful output capability of the TAILSIM program is the ability to display multiple post-run plot pages to allow a quick assessment of the time history response. There are 16 plot pages currently available to the user. Each plot page displays 9 parameters. Each parameter can also be displayed individually. on a one plot-per-page format. For a more refined display of the results the program can also create files of tabulated data. which can then be used by other plotting programs. The TAILSIM program was written straightforwardly assuming the user would want to change the database tables, the buildup equations, the output parameters. and the pilot model parameters. A separate database file and input file are automatically read in by the program. The use of an include file to set up all common blocks facilitates easy changing of parameter names and array sizes.

Hiltner, Dale W.↗

NASA NTRS2002

Diffusion NMR methods applied to xenon gas for materials study

We report initial NMR studies of (i) xenon gas diffusion in model heterogeneous porous media and (ii) continuous flow laser-polarized xenon gas. Both areas utilize the pulsed gradient spin-echo (PGSE) techniques in the gas phase, with the aim of obtaining more sophisticated information than just translational self-diffusion coefficients--a brief overview of this area is provided in the Introduction. The heterogeneous or multiple-length scale model porous media consisted of random packs of mixed glass beads of two different sizes. We focus on observing the approach of the time-dependent gas diffusion coefficient, D(t) (an indicator of mean squared displacement), to the long-time asymptote, with the aim of understanding the long-length scale structural information that may be derived from a heterogeneous porous system. We find that D(t) of imbibed xenon gas at short diffusion times is similar for the mixed bead pack and a pack of the smaller sized beads alone, hence reflecting the pore surface area to volume ratio of the smaller bead sample. The approach of D(t) to the long-time limit follows that of a pack of the larger sized beads alone, although the limiting D(t) for the mixed bead pack is lower, reflecting the lower porosity of the sample compared to that of a pack of mono-sized glass beads. The Pade approximation is used to interpolate D(t) data between the short- and long-time limits. Initial studies of continuous flow laser-polarized xenon gas demonstrate velocity-sensitive imaging of much higher flows than can generally be obtained with liquids (20-200 mm s-1). Gas velocity imaging is, however, found to be limited to a resolution of about 1 mm s-1 owing to the high diffusivity of gases compared with liquids. We also present the first gas-phase NMR scattering, or diffusive-diffraction, data, namely flow-enhanced structural features in the echo attenuation data from laser-polarized xenon flowing through a 2 mm glass bead pack. c2002 John Wiley & Sons, Ltd.

NASA Discipline Life Sciences Technologies↗

NASA NTRS1995

A hybrid neurocomputing/numerical strategy for nonlinear structural analysis

A hybrid neurocomputing/numerical strategy is presented for geometrically nonlinear analysis of structures. The strategy combines model-free data processing capabilities of computational neural networks with a Pade approximants-based perturbation technique to predict partial information about the nonlinear response of structures. In the hybrid strategy, multilayer feedforward neural networks are used to extend the validity of solutions by using training samples produced by Pade approximations to the Taylor series expansion of the response function. The range of validity of the training samples is taken to be the radius of convergence of Pade approximants and is estimated by setting a tolerance on the diverging approximants. The norm of residual vector of unbalanced forces in a given element is used as a measure to assess the quality of network predictions. To further increase the accuracy and the range of network predictions, additional training data are generated by either applying linear regression to weight matrices or expanding the training data by using predicted coefficients in a Taylor series. The effectiveness of the hybrid strategy is assessed by performing large-deflection analysis of a doubly-curved composite panel with a circular cutout, and postbuckling analyses of stiffened composite panels subjected to an in-plane edge shear load. In all the problems considered, the hybrid strategy is used to predict selective information about the structural response, namely the total strain energy and the maximum displacement components only.

Szewczyk, Z. Peter↗

NASA NTRS1993

A hybrid Pade-Galerkin technique for differential equations

A three-step hybrid analysis technique, which successively uses the regular perturbation expansion method, the Pade expansion method, and then a Galerkin approximation, is presented and applied to some model boundary value problems. In the first step of the method, the regular perturbation method is used to construct an approximation to the solution in the form of a finite power series in a small parameter epsilon associated with the problem. In the second step of the method, the series approximation obtained in step one is used to construct a Pade approximation in the form of a rational function in the parameter epsilon. In the third step, the various powers of epsilon which appear in the Pade approximation are replaced by new (unknown) parameters (delta(sub j)). These new parameters are determined by requiring that the residual formed by substituting the new approximation into the governing differential equation is orthogonal to each of the perturbation coordinate functions used in step one. The technique is applied to model problems involving ordinary or partial differential equations. In general, the technique appears to provide good approximations to the solution even when the perturbation and Pade approximations fail to do so. The method is discussed and topics for future investigations are indicated.

Geer, James F.↗

NASA NTRS1984

Exponential-fitted methods for integrating stiff systems of ordinary differential equations: Applications to homogeneous gas-phase chemical kinetics

Conventional algorithms for the numerical integration of ordinary differential equations (ODEs) are based on the use of polynomial functions as interpolants. However, the exact solutions of stiff ODEs behave like decaying exponential functions, which are poorly approximated by polynomials. An obvious choice of interpolant are the exponential functions themselves, or their low-order diagonal Pade (rational function) approximants. A number of explicit, A-stable, integration algorithms were derived from the use of a three-parameter exponential function as interpolant, and their relationship to low-order, polynomial-based and rational-function-based implicit and explicit methods were shown by examining their low-order diagonal Pade approximants. A robust implicit formula was derived by exponential fitting the trapezoidal rule. Application of these algorithms to integration of the ODEs governing homogenous, gas-phase chemical kinetics was demonstrated in a developmental code CREK1D, which compares favorably with the Gear-Hindmarsh code LSODE in spite of the use of a primitive stepsize control strategy.

Pratt, D. T.↗

DOE OSTI2024

Tunable noninteracting free-energy density functionals for high-energy-density physics applications

In this work, we introduce the concept of a tunable noninteracting free-energy density functional and present two examples realized: (i) via a simple one-parameter convex combination of two existing functionals and (ii) via the construction of a generalized gradient approximation (GGA) enhancement factor that contains one free parameter and is designed to satisfy a set of incorporated constraints. Functional (i), constructed as a combination of the local Thomas–Fermi and a pseudopotential-adapted GGA for the noninteracting free-energy, has already demonstrated its practical usability for establishing the high temperature end of the equation of state of deuterium [Phys. Rev. B 104, 144104 (2021)] and CHON resin [Phys. Rev. E 106, 045207 (2022)] for inertial confinement fusion applications. Hugoniot calculations for liquid deuterium are given as another example of how the application of computationally efficient orbital-free density functional theory (OF-DFT) can be utilized with the employment of the developed functionals. Once the functionals have been tuned such that the OF-DFT Hugoniot calculation matches the Kohn–Sham solution at some low-temperature point, agreement with the reference Kohn–Sham results for the rest of the high temperature Hugoniot path is very good with relative errors for compression and pressure on the order of 2% or less.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1993

Rational trigonometric approximations using Fourier series partial sums

A class of approximations (S(sub N,M)) to a periodic function f which uses the ideas of Pade, or rational function, approximations based on the Fourier series representation of f, rather than on the Taylor series representation of f, is introduced and studied. Each approximation S(sub N,M) is the quotient of a trigonometric polynomial of degree N and a trigonometric polynomial of degree M. The coefficients in these polynomials are determined by requiring that an appropriate number of the Fourier coefficients of S(sub N,M) agree with those of f. Explicit expressions are derived for these coefficients in terms of the Fourier coefficients of f. It is proven that these 'Fourier-Pade' approximations converge point-wise to (f(x(exp +))+f(x(exp -)))/2 more rapidly (in some cases by a factor of 1/k(exp 2M)) than the Fourier series partial sums on which they are based. The approximations are illustrated by several examples and an application to the solution of an initial, boundary value problem for the simple heat equation is presented.

Geer, James F.↗

DOE OSTI2026

Comparison between explicit and implicit discretization strategies for a dissipative thermal environment

We investigate strategies for simulating open quantum systems coupled to dissipative baths by comparing explicit wave function-based discretization [via multi-layer multi-configuration time-dependent Hartree (ML-MCTDH)] and the implicit density matrix-based master equation method [via tree tensor network hierarchical equations of motion (TTN-HEOM)]. For dissipative baths characterized by exponentially decaying bath correlation functions, the implicit discretization approach of HEOM—rooted in bath correlation function decompositions—proves significantly more efficient than explicit discretization of the bath into discrete harmonic modes. Explicit methods, like ML-MCTDH, require extensive mode discretization to approximate continuum baths, leading to computational bottlenecks. Case studies for two-level systems and a Fenna–Matthews–Olson complex model highlight TTN-HEOM’s superiority in capturing dissipative dynamics with relaxations with a minimal number of auxiliary modes, while the explicit methods are as exact as the HEOM in pure dephasing regimes. This comparison is enabled by the TENSO package, which has both ML-MCTDH and TTN-HEOM implemented using the same computational structure and propagation strategy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Jets from shocked metal surfaces with grooves: Missing experiments

Many studies have investigated the mass outflows generated when a planar shock transits an imperfect (“defected”) metal surface, where the defects are symmetric triangular or sinusoidal grooves. Yet a fundamental question remains unanswered: how does the quantity of outflow mass and its maximum velocity vary as a function of the groove cross-sectional aspect ratio? We identify two sets of missing experiments that must be addressed to answer the question. The aspect ratio (groove depth over width) is equivalently represented by θ, the cross-sectional half angle, or by η 0 k, the amplitude multiplied by an effective wavenumber. Low θ (high η 0 k) grooves comprise the first set of missing experiments, which are necessary to determine the validity of theoretical predictions of the nonlinear regime (η 0 k ≥ 1, θ < 57.5°). The second set of missing experiments are those in which the volume of the groove (or equivalently, the axial cross-sectional area) has been held constant as θ or η 0 k are varied. Such experiments are necessary to independently measure the effects of variations in groove volume and groove aspect ratio on the resulting jets.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1973

Control law synthesis and sensor design for active flutter suppression.

Methods are presented for representing unsteady aerodynamic loadings, valid for arbitrary motion, for the two-dimensional typical section with a trailing-edge control surface in incompressible flow and a three-dimensional lifting surface with leading- or trailing-edge control surfaces in subsonic compressible flow. Loads for the two-dimensional incompressible case are obtained via analytic continuation of the Theodorsen function into the complex plane, with guidance from a time-domain approximation of the Wagner indicial function. For the three-dimensional case, it is proposed that oscillatory generalized aerodynamic forces be approximated by Pade fractions, thereby permitting an obvious continuation into the complex plane. A theoretical justification for this procedure is briefly outlined, and some practical aspects of its implementation are discussed.

Lyons, M. G.↗

NASA NTRS1979

Spectropolarimetry of Venus and Jupiter clouds - Information content of equivalent widths

The theory of the curve-of-growth for absorption lines formed in an anisotropically scattering medium, particularly the clouds of Venus and Jupiter, is generalized to the case of polarization. When the photon path-length distributions associated with each of Stokes' parameters (Il, Ir, and Q) are evaluated using Pade's approximations, analytical expressions can be derived for the corresponding equivalent widths. It is found that spectra of Q, derivable from the rotational temperature of the Q-lines, are formed at a level closer to the cloud top than the photopolarimetric level, rendering a more accurate location of the cloud top than possible by photopolarization. It is also demonstrated that simultaneous measurements of the equivalent widths of strong I- and Q-lines can provide analytical relations for retrieving the scale-height ratios of ambient gas-to-cloud particles and minor constituent-to-ambient gas molecules without resorting to any other independent measurement.

Buriez, J. C.↗