Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Taylor 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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

STEP and STEPSPL: Computer programs for aerodynamic model structure determination and parameter estimation

The successful parametric modeling of the aerodynamics for an airplane operating at high angles of attack or sideslip is performed in two phases. First the aerodynamic model structure must be determined and second the associated aerodynamic parameters (stability and control derivatives) must be estimated for that model. The purpose of this paper is to document two versions of a stepwise regression computer program which were developed for the determination of airplane aerodynamic model structure and to provide two examples of their use on computer generated data. References are provided for the application of the programs to real flight data. The two computer programs that are the subject of this report, STEP and STEPSPL, are written in FORTRAN IV (ANSI l966) compatible with a CDC FTN4 compiler. Both programs are adaptations of a standard forward stepwise regression algorithm. The purpose of the adaptation is to facilitate the selection of a adequate mathematical model of the aerodynamic force and moment coefficients of an airplane from flight test data. The major difference between STEP and STEPSPL is in the basis for the model. The basis for the model in STEP is the standard polynomial Taylor's series expansion of the aerodynamic function about some steady-state trim condition. Program STEPSPL utilizes a set of spline basis functions.

Batterson, J. G.↗
NASA NTRS1989

The art of spacecraft design: A multidisciplinary challenge

Actual design turn-around time has become shorter due to the use of optimization techniques which have been introduced into the design process. It seems that what, how and when to use these optimization techniques may be the key factor for future aircraft engineering operations. Another important aspect of this technique is that complex physical phenomena can be modeled by a simple mathematical equation. The new powerful multilevel methodology reduces time-consuming analysis significantly while maintaining the coupling effects. This simultaneous analysis method stems from the implicit function theorem and system sensitivity derivatives of input variables. Use of the Taylor's series expansion and finite differencing technique for sensitivity derivatives in each discipline makes this approach unique for screening dominant variables from nondominant variables. In this study, the current Computational Fluid Dynamics (CFD) aerodynamic and sensitivity derivative/optimization techniques are applied for a simple cone-type forebody of a high-speed vehicle configuration to understand basic aerodynamic/structure interaction in a hypersonic flight condition.

Abdi, F.↗
NASA NTRS1993

Annular suspension and pointing system with controlled DC electromagnets

The Annular Suspension and Pointing System (ASPS) developed by the Flight System division of Sperry Corporation is a six-degree of freedom payload pointing system designed for use with the space shuttle. This magnetic suspension and pointing system provides precise controlled pointing in six-degrees of freedom, isolation of payload-carrier disturbances, and end mount controlled pointing. Those are great advantages over the traditional mechanical joints for space applications. In this design, we first analyzed the assumed model of the single degree ASPS bearing actuator and obtained the plant dynamics equations. By linearizing the plant dynamics equations, we designed the cascade and feedback compensators such that a stable and satisfied result was obtained. The specified feedback compensator was computer simulated with the nonlinearized plant dynamics equations. The results indicated that an unstable output occurred. In other words, the designed feedback compensator failed. The failure of the design is due to the Taylor's series expansion not converging.

Vu, Josephine Lynn↗
NASA NTRS1973

Rapid iterative reanalysis for automated design

A method for iterative reanalysis in automated structural design is presented for a finite-element analysis using the direct stiffness approach. A basic feature of the method is that the generalized stiffness and inertia matrices are expressed as functions of structural design parameters, and these generalized matrices are expanded in Taylor series about the initial design. Only the linear terms are retained in the expansions. The method is approximate because it uses static condensation, modal reduction, and the linear Taylor series expansions. The exact linear representation of the expansions of the generalized matrices is also described and a basis for the present method is established. Results of applications of the present method to the recalculation of the natural frequencies of two simple platelike structural models are presented and compared with results obtained by using a commonly applied analysis procedure used as a reference. In general, the results are in good agreement. A comparison of the computer times required for the use of the present method and the reference method indicated that the present method required substantially less time for reanalysis. Although the results presented are for relatively small-order problems, the present method will become more efficient relative to the reference method as the problem size increases. An extension of the present method to static reanalysis is described, ana a basis for unifying the static and dynamic reanalysis procedures is presented.

Bhatia, K. G.↗
DOE OSTI2024

Calculation of machine precision second order derivatives using dual-complex numbers

It is well known that both complex and dual numbers can be employed to obtain machine precision first-order derivatives; however, neither, on their own, can compute machine precision 2nd order derivatives. To address this limitation, it is demonstrated in this paper that combined dual-complex numbers can be used to compute machine precision 1st and 2nd order derivatives. The dual-complex approach is simpler than utilizing multicomplex or hyper-dual numbers as existing dual libraries can be used as is or easily augmented to accept complex numbers, and the complexity of developing, integrating, and deploying multicomplex or hyper-dual libraries is avoided. The efficacy of this approach is demonstrated for both univariate and multivariate functions. Finally, source code examples using the Python, Julia, and Mathematica languages are provided as supplemental material.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Validation of Nuclear Data Sensitivity Calculations by the MCNP PERT Card [Poster]

This work validates a method of estimating sensitivities that utilized the MCNP PERT card and tallies. The agreement of the PERT card and Brute Force results proves the validity of calculating nuclear data sensitivities by use of the MCNP PERT card, especially at higher energies. The overall agreement will be easier to understand once an uncertainty analysis is completed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1986

Carbon/Carbon Pistons for Internal Combustion Engines

Carbon/carbon piston performs same function as aluminum pistons in reciprocating internal combustion engines while reducing weight and increasing mechanical and thermal efficiencies of engine. Carbon/carbon piston concept features low piston-to-cylinder wall clearance - so low piston rings and skirts unnecessary. Advantages possible by negligible coefficient of thermal expansion of carbon/carbon.

Taylor, A. H.↗
NASA NTRS1995

An approach to the development of numerical algorithms for first order linear hyperbolic systems in multiple space dimensions: The constant coefficient case

Two methods for developing high order single step explicit algorithms on symmetric stencils with data on only one time level are presented. Examples are given for the convection and linearized Euler equations with up to the eighth order accuracy in both space and time in one space dimension, and up to the sixth in two space dimensions. The method of characteristics is generalized to nondiagonalizable hyperbolic systems by using exact local polynominal solutions of the system, and the resulting exact propagator methods automatically incorporate the correct multidimensional wave propagation dynamics. Multivariate Taylor or Cauchy-Kowaleskaya expansions are also used to develop algorithms. Both of these methods can be applied to obtain algorithms of arbitrarily high order for hyperbolic systems in multiple space dimensions. Cross derivatives are included in the local approximations used to develop the algorithms in this paper in order to obtain high order accuracy, and improved isotropy and stability. Efficiency in meeting global error bounds is an important criterion for evaluating algorithms, and the higher order algorithms are shown to be up to several orders of magnitude more efficient even though they are more complex. Stable high order boundary conditions for the linearized Euler equations are developed in one space dimension, and demonstrated in two space dimensions.

Goodrich, John W.↗
NASA NTRS2000

Solving ODE Initial Value Problems With Implicit Taylor Series Methods

In this paper we introduce a new class of numerical methods for integrating ODE initial value problems. Specifically, we propose an extension of the Taylor series method which significantly improves its accuracy and stability while also increasing its range of applicability. To advance the solution from t (sub n) to t (sub n+1), we expand a series about the intermediate point t (sub n+mu):=t (sub n) + mu h, where h is the stepsize and mu is an arbitrary parameter called an expansion coefficient. We show that, in general, a Taylor series of degree k has exactly k expansion coefficients which raise its order of accuracy. The accuracy is raised by one order if k is odd, and by two orders if k is even. In addition, if k is three or greater, local extrapolation can be used to raise the accuracy two additional orders. We also examine stability for the problem y'= lambda y, Re (lambda) less than 0, and identify several A-stable schemes. Numerical results are presented for both fixed and variable stepsizes. It is shown that implicit Taylor series methods provide an effective integration tool for most problems, including stiff systems and ODE's with a singular point.

Scott, James R.↗
NASA NTRS1984

Stress-derivative control of keystoning deformation in finite element codes

The stress derivative technique for control of keystoning deformation in under-integrated finte elements is based on expansion of the stress in a Taylor series about the element center and retention of additional terms beyond the constant stress term. It has the advantage over other control techniques that keystoning resistance is provided by actual rather than artificial material properties. Application of this technique to the quadrilateral ring elements used for modelling solids of revolution subjected to axisymmetric loads is described. In a cylindrical coordinate system additional terms appear in the formulation which must be dealt with in arriving at a workable keystoning control scheme.

Schulz, J. C.↗
NASA NTRS1988

Expansion tube test time predictions

The interaction of an interface between two gases and strong expansion is investigated and the effect on flow in an expansion tube is examined. Two mechanisms for the unsteady Pitot-pressure fluctuations found in the test section of an expansion tube are proposed. The first mechanism depends on the Rayleigh-Taylor instability of the driver-test gas interface in the presence of a strong expansion. The second mechanism depends on the reflection of the strong expansion from the interface. Predictions compare favorably with experimental results. The theory is expected to be independent of the absolute values of the initial expansion tube filling pressures.

Gourlay, Christopher M.↗
NASA NTRS1994

Thermonuclear runaways in nova outbursts. 2: Effect of strong, instantaneous, local fluctuations

In an attempt to understand the manner in which nova outbursts are initiated on the surface of a white dwarf, we investigate the effects fluctuations have on the evolution of a thermonuclear runaway. Fluctuations in temperature density, or the composition of material in the burning shell may arise due to the chaotic flow field generated by convection when it occurs, or by the accretion process itself. With the aid of two-dimensional reactive flow calculations, we consider cases where a strong fluctutation in temperature arises during the early, quiescent accretion phase or during the later, more dynamic, explosion phase. In all cases we find that an instantaneous, local temperature fluctuation causes the affected material to become Rayleigh-Taylor unstable. The rapid rise and subsequent expansion of matter immediately cools the hot blob, which prevents the lateral propagation of burning. This suggests that local temperature fluctuations do not play a significant role in directly initiating the runaway, especially during the early stages. However, they may provide an efficient mechanism of mixing core material into the envelope (thereby pre-enriching the fuel for subsequent episodes of explosive hydrogen burning) and of mixing substantial amounts of the radioactive nucleus N-13 into the surface layers, making novae potential gamma-ray sources. This suggests that it is the global not the local, evolution of the core-envelope interface to high temperatures which dominates the development of the runaway. We also present a possible new scenario for the initiation of nova outbursts based on our results.

Shankar, Anurag↗
DOE OSTI2023

Interface models for three-dimensional Rayleigh–Taylor instability

We derive interface models for three-dimensional Rayleigh–Taylor instability (RTI), making use of a novel asymptotic expansion in the non-locality of the fluid flow. These interface models are derived for the purpose of studying universal features associated with RTI such as the Froude number in single-mode RTI, the predicted quadratic growth of the interface amplitude under multi-mode random perturbations, the optimal (viscous) mixing rates induced by the RTI and the self-similarity of horizontally averaged density profiles and the remarkable stabilization of the mixing layer growth rate which arises for the three-fluid two-interface heavy–light–heavy configuration, in which the addition of a third fluid bulk slows the growth of the mixing layer to a linear rate. Our interface models can capture the formation of small-scale structures induced by severe interface roll-up, reproduce experimental data in a number of different regimes and study the effects of multiple interface interactions even as the interface separation distance becomes exceedingly small. Compared with traditional numerical schemes used to study such phenomena, our models provide a computational speed-up of at least two orders of magnitude.

Mechanics↗
NASA NTRS1994

The delayed-detonation model of a type Ia supernovae. 1: The deflagration phase

The nature of the 'delayed detonation' mechanism of Khokhlov for the explosion of Type Ia supernovae is investigated by using two-dimensional numerical hydrodynamics simulations. A new algorithm is used to treat the deflagration front. Assuming that it propagates locally at the laminar flame speed, the deflagration is insufficient to unbind the star. Expansion shuts of the flame; much of this small production of iron group nuclei occurs at lower densities, which reduces the electron-capture problem. The burning front does become wrinkled, but the wavelength of the instability is much larger than the computational grid size and is resolved; this is consistent with previous analysis. Because the degenerate star has an adiabatic exponent only slightly above 4/3, the energy released by deflagration drives a pulsation of large amplitude. During the first expansion phase, adiabatic cooling shuts off the burning, and a Rayleigh-Taylor instability then gives mixing of high-entropy ashes with low-entropy fuel. During the first contraction phase, compressional heating reignites the material. This paper deals with the deflagration phase, from the onset of burning, through expansion and quenching of the flame, to the first contraction.

Arnett, David↗
NASA NTRS1974

Reduction of acoustic disturbances in the test section of supersonic wind tunnels by laminarizing their nozzle and test section wall boundary layers by means of suction

The feasibility of quiet, suction laminarized, high Reynolds number (Re) supersonic wind tunnel nozzles was studied. According to nozzle wall boundary layer development and stability studies, relatively weak area suction can prevent amplified nozzle wall TS (Tollmien-Schlichting) boundary layer oscillations. Stronger suction is needed in and shortly upstream of the supersonic concave curvature nozzle area to avoid transition due to amplified TG (Taylor-Goertler) vortices. To control TG instability, moderately rapid and slow expansion nozzles require smaller total suction rates than rapid expansion nozzles, at the cost of larger nozzle length Re and increased TS disturbances. Test section mean flow irregularities can be minimized with suction through longitudinal or highly swept slots (swept behind local Mach cone) as well as finely perforated surfaces. Longitudinal slot suction is optimized when the suction-induced crossflow velocity increases linearly with surface distance from the slot attachment line toward the slot (through suitable slot geometry). Suction in supersonic blowdown tunnels may be operated by one or several individual vacuum spheres.

Pfenninger, W.↗