Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical differentiation”

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 1,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
NASA NTRS1972

A discrete variable approximation to minimum weight panel designs subject to a supersonic flutter-speed constraint.

A numerical technique is presented that determines the minimum-weight thickness distribution for a one-dimensional, simply supported panel subjected to a parallel supersonic flow on one side. An aeroelastic eigenvalue, which characterizes flutter speed, is held constant. The governing differential equations are approximated by sets of difference equations adjoined to the weight function via a penalty function. A conjugate gradient method then solves the resulting sequence of unconstrained minimization problems. Numerical results are obtained without imposing a minimum thickness constraint; the effect of constant inplane tensile and compressive stresses is also investigated.

Pierson, B. L.↗
NASA NTRS2001

Candle Flames in Microgravity

The goal of this work is to study both experimentally and numerically the behavior of a candle flame burning in a microgravity environment. Two space experiments (Shuttle and Mir) have shown the candle flame in microgravity to be small (approximately 1.5 cm diameter), dim blue, and hemispherical. Near steady flames with very long flame lifetimes (up to 45 minutes in some tests) existed for many of the tests. Most of the flames spontaneously oscillated with a period of approximately 1 Hz just prior to extinction). In a previous model of candle flame in microgravity, a porous sphere wetted with liquid fuel simulated the evaporating wick. The sphere, with a temperature equal to the boiling temperature of the fuel, was at the end of an inert cone that had a prescribed temperature. This inert cone produces the quenching effect of the candle wax in the real configuration. Although the computed flame shape resembled that observed in the microgravity experiment, the model was not able to differentiate the effect of wick geometry, e.g., a long vs. a short wick. This paper presents recent developments in the numerical model of the candle flame. The primary focus has been to more realistically account for the actual shape of the candle.

Dietrich, D. L.↗
NASA NTRS1972

Earth zonal harmonics from rapid numerical analysis of long satellite arcs

A zonal geopotential is presented to degree 21 from evaluation of mean elements for 21 satellites including 2 of low inclination. Each satellite is represented by an arc of at least one apsidal rotation. The lengths range from 200 to 800 days. Differential correction of the initial elements in all of the arcs, together with radiation pressure and atmospheric drag coefficients, was accomplished simultaneously with the correction for the zonal harmonics. The satellite orbits and their variations are generated by numerical integration of the Lagrange equations for mean elements. Disturbances due to precession and nutation of the earth's pole, atmospheric drag, radiation pressure and luni-solar gravity are added at from 1- to 8-day intervals in the integrated orbits. The results agree well with recent solutions from other authors using different methods and different satellite sets.

Wagner, C. A.↗
NASA NTRS1973

Aircraft maneuver optimization of reduced-order approximation

Review of recent work in recasting the energy type of approximation to aircraft flight in terms of singular perturbation theory and in extending this theory to three-dimensional maneuvers. Singular perturbations for differential equations arising in optimal control are first examined for a system of fairly general form but low order. The attitude dynamics for optimal flight of a rocket in vacuum is then studied as an introductory, and fairly transparent, example. The question of the choice of variables is then discussed. Finally, optimal aircraft flight in various reduced-order approximations is investigated. In particular, the problem of three-dimensional aircraft flight is formulated for singular perturbation treatment, and possibilities for decoupling into several lower-order problems are illustrated. The use of reduced-order approximation facilitates numerical computations by reducing the number of multiplier initial values that must be determined simultaneously and by improving the conditioning of the differential equations.

Kelley, H. J.↗
DOE OSTI2019

Recursion relation for general 3d blocks

We derive closed-form expressions for all ingredients of the Zamolodchikov-like recursion relation for general spinning conformal blocks in 3-dimensional conformal field theory. This result opens a path to efficient automatic generation of conformal block tables, which has immediate applications in numerical conformal bootstrap program. Our derivation is based on an understanding of null states and conformally-invariant differential operators in momentum space, combined with a careful choice of the relevant tensor structures bases. This derivation generalizes straightforwardly to higher spacetime dimensions d, provided the relevant Clebsch-Gordan coefficients of Spin (d) are known.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Temporally quasiperiodic data, propagating in the laboratory frame, can be rendered periodic by Galilean transformation

For a broad class of distributions of temperature, concentration, or another quantity propagating rectilinearly, we show that temporally quasiperiodic behavior in the laboratory frame can be rendered periodic by Galilean transformation. The approach is illustrated analytically and numerically using as an example a closed-form model distribution generated from a one-dimensional partial differential equation, and a detailed process is developed to determine frame speed from more general quasiperiodic, one-dimensional, temporally- and spatially-discretized data. Furthermore, the approach is extended to two- and three-dimensional rectilinear propagation, and its application to nonrectilinear propagation, along with implications for interpreting noise-corrupted data, are also discussed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

A robust and efficient generalized real cascade model

This paper describes a new separation cascade model that allows for arbitrary cascade connectivity and includes the effects of mixing, material losses and generalized separation element performance. The model also accommodates mixtures with an arbitrary number of component species. The mathematical and numerical framework of this model includes both time-dependent and steady state cascade simulation. Robust numerical approximations of the continuous mathematical model are accomplished using the Sundials suite of nonlinear algebraic and differential-algebraic equation solvers. Verification of the model is accomplished using manufactured solutions while validation is demonstrated by comparison to published experimental results and benchmarks.

36 MATERIALS SCIENCE↗
NASA NTRS1963

Application of Dorodnitsyn's Integral Method to Nonequilibrium Flows over Pointed Bodies

Dorodnitsyn's integral method is used to obtain an approximate solution to the supersonic nonequilibrium flow over pointed bodies with attached shock waves. The partial differential equations governing the flow are converted to an approximate set of ordinary equations, which are solved by numerical integration starting at the body tip. Detailed analytical and numerical results for the first approximation are presented, considering the vibrational relaxation of a diatomic gas over a wedge or cone. It is shown that the first approximation yields: (1) The exact flow-variable gradients at the wedge tip (2) Expressions for the flow-variable gradients at the cone tip which are in agreement with extrapolations of characteristics calculations (3) A good approximate algebraic solution for frozen or equilibrium conical flow (4) An approximate expression for the nonequilibrium-flow stream function which affords a means of obtaining variations across the shock layer of the temperature and vibrational energy. Numerical results for both the wedge and cone compare favorably with identical cases computed by the method of characteristics.

NONEQUILIBRIUM FLOW↗
NASA NTRS1974

Two-demensional analysis of heat and mass transfer in porous media using the strongly implicit procedure

Numerical results of the heat and mass transfer in a porous matrix are presented. The coupled, nonlinear partial differential equations describing this physical phenomenon are solved in finite difference form for two dimensions, using a new iterative technique (the strongly implicit procedure). The influence of the external environment conditions (heating and pressure) is shown to produce two-dimensional flow in the porous matrix. Typical fluid and solid temperature distributions in the porous matrix and internal pressure distributions are presented.

Curry, D. M.↗
NASA NTRS1977

Variants and extensions of a fast direct numerical cauchy-riemann solver, with illustrative applications

Revised and extended versions of a fast, direct (noniterative) numerical Cauchy-Riemann solver are presented for solving finite difference approximations of first order systems of partial differential equations. Although the difference operators treated are linear and elliptic, one significant application of these extended direct Cauchy-Riemann solvers is in the fast, semidirect (iterative) solution of fluid dynamic problems governed by the nonlinear mixed elliptic-hyperbolic equations of transonic flow. Different versions of the algorithms are derived and the corresponding FORTRAN computer programs for a simple example problem are described and listed. The algorithms are demonstrated to be efficient and accurate.

Martin, E. D.↗
NASA NTRS1986

Physical and numerical sources of computational inefficiency in integration of chemical kinetic rate equations: Etiology, treatment and prognosis

The design of a very fast, automatic black-box code for homogeneous, gas-phase chemical kinetics problems requires an understanding of the physical and numerical sources of computational inefficiency. Some major sources reviewed in this report are stiffness of the governing ordinary differential equations (ODE's) and its detection, choice of appropriate method (i.e., integration algorithm plus step-size control strategy), nonphysical initial conditions, and too frequent evaluation of thermochemical and kinetic properties. Specific techniques are recommended (and some advised against) for improving or overcoming the identified problem areas. It is argued that, because reactive species increase exponentially with time during induction, and all species exhibit asymptotic, exponential decay with time during equilibration, exponential-fitted integration algorithms are inherently more accurate for kinetics modeling than classical, polynomial-interpolant methods for the same computational work. But current codes using the exponential-fitted method lack the sophisticated stepsize-control logic of existing black-box ODE solver codes, such as EPISODE and LSODE. The ultimate chemical kinetics code does not exist yet, but the general characteristics of such a code are becoming apparent.

Pratt, D. T.↗
NASA NTRS1988

Evolution of inclined galactic gas disks. I - A cloud-fluid approach

In this paper, numerical models are developed to follow the time-dependent behavior of inclined gas disks in nonspherical galaxies. Differential equations are developed which describe the time rate of change of the orbits of mass elements composing a thin, smoothly warping and twisting gas disk. These equations include multiple moments due to the nonsphericity of the galactic mass distribution, Coriolis effects arising from determination of orbits in a reference frame which tumbles with the surface figure of the galaxy, and viscous forces arising from cloud-cloud collisions. A cloud-fluid approach is used in determining the viscous forces. Physical arguments are used to determine the maximum permissible local value of the coefficient of kinematic viscosity. For small initial inclinations, analytic expressions are obtained for the inclination as a function of time and for the settling time, and the features of the solution are discussed.

Steiman-Cameron, Thomas Y.↗
NASA NTRS1993

Conservative high-order-accurate finite-difference methods for curvilinear grids

Two fourth-order-accurate finite-difference methods for numerically solving hyperbolic systems of conservation equations on smooth curvilinear grids are presented. The first method uses the differential form of the conservation equations; the second method uses the integral form of the conservation equations. Modifications to these schemes, which are required near boundaries to maintain overall high-order accuracy, are discussed. An analysis that demonstrates the stability of the modified schemes is also provided. Modifications to one of the schemes to make it total variation diminishing (TVD) are also discussed. Results that demonstrate the high-order accuracy of both schemes are included in the paper. In particular, a Ringleb-flow computation demonstrates the high-order accuracy and the stability of the boundary and near-boundary procedures. A second computation of supersonic flow over a cylinder demonstrates the shock-capturing capability of the TVD methodology. An important contribution of this paper is the dear demonstration that higher order accuracy leads to increased computational efficiency.

Rai, Man M.↗
NASA NTRS1995

Transport phenomena in the micropores of plug-type phase separators

This study numerically investigates the transport phenomena within and across a porous-plug phase separator. The effect of temperature differential across a single pore and of the sidewall boundary conditions, i.e., isothermal or linear thermal gradient, are presented and discussed. The effects are quantified in terms of the evaporation mass flux across the boundary and the mean surface temperature. A two-dimensional finite element model is used to solve the continuity, momentum, and energy equations for the liquid. Temperature differentials across the pore interface of 1.0, and 1.5 K are examined and their effect on evaporation flux and mean surface temperature is shown. For isothermal side boundary conditions, the evaporation flux across the pore is directly proportional and linear with Delta T. For the case of an imposed linear thermal gradient on the side boundaries, Biot numbers of 0.0, 0.15, and 0.5 are examined. The most significant effect of Biot number is to lower the overall surface temperature and evaporation flux.

Fazah, M. M.↗
NASA NTRS1996

Time-Dependent Thermally-Driven Interfacial Flows in Multilayered Fluid Structures

A computational study of thermally-driven convection in multilayered fluid structures will be performed to examine the effect of interactions among deformable fluid-fluid interfaces on the structure of time-dependent flow in these systems. Multilayered fluid structures in two models configurations will be considered: the differentially heated rectangular cavity with a free surface, and the encapsulated cylindrical liquid bridge. An extension of a numerical method developed as part of our recent NASA Fluid Physics grant will be used to account for finite deformations of fluid-fluid interfaces.

Haj-Hariri, Hossein↗
NASA NTRS1997

Chemical and Thermal Analysis

Work during the past three years has included significant research in several areas aimed at further clarification of the aging and chemical failure mechanism of thermoplastics (PVDF or Tefzel) for pipes. Among the areas investigated were the crystallinity changes associated with both the Coflon and Tefzel after various simulated environmental exposures using X-Ray diffraction analysis. We have found that significant changes in polymer crystallinity levels occur as a function of the exposures. These crystallinity changes may have important consequences on the fracture, fatigue, tensile, and chemical resistance of the materials. We have also noted changes in the molecular weight distribution and the increased crosslinking of the Coflon material using Gel Permeation Chromatographic Analysis. Again these changes may result in variations in the mechanical and chemical properties in the material. We conducted numerous analytical studies with methods including X-ray Diffraction, Gel Permeation Chromatography, Fourier Transform Infrared Spectroscopy, and Differential Scanning Calorimetry. We investigated a plethora of aged samples of both Tefzel and Coflon that were forwarded from MERL. Pressurized tests were performed on powdered PVDF in a modified Fluid A, which we will call A-2. In this case the ethylene diamine concentration was increased to 3 percent in methanol. Coflon pipe sections and powdered Coflon were exposed in pressure cells at 1700 psi at three separate test temperatures.

Bulluck, J. W.↗