Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical integration”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Effects of anisotropic conduction and heat pipe interaction on minimum mass space radiators

Equations are formulated for the two dimensional, anisotropic conduction of heat in space radiator fins. The transverse temperature field was obtained by the integral method, and the axial field by numerical integration. A shape factor, defined for the axial boundary condition, simplifies the analysis and renders the results applicable to general heat pipe/conduction fin interface designs. The thermal results are summarized in terms of the fin efficiency, a radiation/axial conductance number, and a transverse conductance surface Biot number. These relations, together with those for mass distribution between fins and heat pipes, were used in predicting the minimum radiator mass for fixed thermal properties and fin efficiency. This mass is found to decrease monotonically with increasing fin conductivity. Sensitivities of the minimum mass designs to the problem parameters are determined.

Baker, Karl W.↗
NASA NTRS1994

Algorithms For Integrating Nonlinear Differential Equations

Improved algorithms developed for use in numerical integration of systems of nonhomogenous, nonlinear, first-order, ordinary differential equations. In comparison with integration algorithms, these algorithms offer greater stability and accuracy. Several asymptotically correct, thereby enabling retention of stability and accuracy when large increments of independent variable used. Accuracies attainable demonstrated by applying them to systems of nonlinear, first-order, differential equations that arise in study of viscoplastic behavior, spread of acquired immune-deficiency syndrome (AIDS) virus and predator/prey populations.

Freed, A. D.↗
NASA NTRS1982

A numerical method for integrating the kinetic equations of droplet spectra evolution by condensation/evaporation and by coalescence/breakup processes

An extension of the method of moments is developed for the numerical integration of the kinetic equations of droplet spectra evolution by condensation/evaporation and by coalescence/breakup processes. The number density function n sub k (x,t) in each separate droplet packet between droplet mass grid points (x sub k, x sub k+1) is represented by an expansion in orthogonal polynomials with a given weighting function. In this way droplet number concentrations, liquid water contents and other moments in each droplet packet are conserved and the problem of solving the kinetic equations is replaced by one of solving a set of coupled differential equations for the number density function moments. The method is tested against analytic solutions of the corresponding kinetic equations. Numerical results are obtained for different coalescence/breakup and condensation/evaporation kernels and for different initial droplet spectra. Also droplet mass grid intervals, weighting functions, and time steps are varied.

Emukashvily, I. M.↗
NASA NTRS1966

Numerical Studies of the Non-linear Vlasov Equation

The subject of this dissertation is the numerical integration of the initial-value problem for the non-linear Vlasov equation. The Vlasov equation is used to describe the dynamics of a "collisionless", one-dimensional, classical electron gas confined between two perfectly reflecting boundaries. Only the long-range Coulomb interactions of the electrons are considered; effects associated with the discrete structure are neglected. The numerical results obtained for non-linear Landau damping compare well with similar results obtained by Knorr. A general statement of the results on stable initial conditions is: As the degree of non-linearity of the initial conditions is increased, the deviation from linear Landau damping appears sooner and is more severe. In some cases damping was observed to cease. Curves showing the time dependence of the damping decrement are derived and compared with predictions of non-linear theories. New results obtained in this study include the observation that for strongly non-linear cases, the damping of the electric field causes an initially Maxwellian fo (v, O) to develop a peak in the neighborhood of the phase velocity; strong growth of the second harmonic is seen after fo (v, t) develops such a peak. Also new in this study is the interpretation of the development of a certain class of strongly unstable initial conditions as approaching an inhomogeneous equilibrium.

NUMERICAL INTEGRATION↗
NASA NTRS1975

Integrated potential formulation of unsteady supersonic aerodynamics for interacting wings

A numerical integrated velocity potential method for the determination of unsteady aerodynamic forces on arbitrary interacting wings and tails in supersonic flow has been developed. Constant Mach number has been assumed throughout the flow field. Normalwash and sidewash integrals have been derived. The upwash integral remains to be derived. Singular integrals in the expressions for the velocity components have been evaluated in closed form. Lifting surfaces are represented by triangular elements defined by arbitrarily spaced characteristic lines and the true surface edges. The wake field is represented by rectangular strip elements. Velocity potential distributions and generalized aerodynamic coefficients have been compared.

Appa, K.↗
NASA NTRS1979

Shuttle program. Onorbit navigation integrator results for typical shuttle orbits

Three types of navigation onorbit numerical integrators were evaluated: (1) power integrators with no delta-V incorporation, just coasting (using Taylor series expansion integrators); (2) coasting integrators using the Cowell method of special perturbations; and (3) coasting integrator using the Pines variation of parameter perturbation method. Results show that the super G integrator is a very simple and effective for 2 and 4 second time steps. Since IMU delta-V data can be easily incorporated in the integration scheme, its use as the standard onorbit navigation propagator for the maintenance of the current state was implemented in the onboard navigation software. The Pines formulation method with a Runge-Kutta-Gill fourth-order integrator method produces excellent results up to 300 second time steps. On orbit prediction with this method was implemented in the onboard onorbit navigation scheme. The Runge-Kutta third order, using Cowell's method, is an excellent general purpose determination integrator for time steps up to a 60 second duration.

Olszewski, O. W.↗
NASA NTRS1974

Nonlinearities in analyses of unsteady flow around oscillating wings

An analysis of nonlinearities involved in unsteady subsonic viscous flow around wings at high Reynolds numbers is presented. The finite element method was applied for the numerical integration of the vorticity transport equations. Physical and numerical aspects of the nonlinearities are discussed. Magnification of nonlinearities in the convection terms near wing surfaces due to changing geometry at leading and trailing edges of the wing chord is determined. As a bubble is generated near the leading edge, the timewise variation of vorticity and pressure distributions over the surface is shown. The different procedures for the integration of the vorticity transport equations in three dimensions along spanwise, chordwise, and normal directions are discussed.

Bratanow, T.↗
NASA NTRS1975

A three-dimensional viscoelastic model of a strike slip fault

An analytic approximation to the Green's function for the displacements due to a strike slip point source in an elastic layer over a viscoelastic half-space is developed. This approximate Green's function is useful because it can be analytically integrated over the fault surface. Comparison with a numerical integration of the exact solution integral indicates that the approximation is quite good. The approximate Green's function is integrated analytically to obtain the displacements due to a finite rectangular strike slip fault in an elastic layer over a viscoelastic half-space. Ground displacements and angle changes from a model survey net are computed to illustrate the viscoelastic relaxation which follows a fracture in the elastic region.

Rundle, J. B.↗
NASA NTRS1977

Unsteady supersonic aerodynamic theory by the method of potential gradient

A generalized solution of the hyperbolic wave equation has been derived. The method used has been further developed to relate the velocity components at a field point to the potential gradient distribution in the dependence domain. Singular integrals have been evaluated in closed form, while numerical integration methods are suggested for treating more complex but analytic functions. Idealization of the lifting surfaces by trapezoidal elements with two sides parallel to the streamlines is computationally efficient because line integrations along the other two sides need only be considered. Furthermore, all the integrands vanish on the Mach cone and the need for determining the hyperbolic curves of intersection of the cone with the lifting surface is avoided. Generalized aerodynamic coefficients for three AGARD planforms have been calculated and compared with the available results.

Jones, W. P.↗
NASA NTRS1972

The influence of laser ranging on selenodetic control.

In this paper, a mathematical model is presented which can be used to perform an adjustment using laser distances to improve coordinates of points of the earth and on the moon. The observation equations are given and the orientation of the earth and the moon is defined in terms of three Eulerian angles. Parameters related to the orientation of the two bodies are the six initial conditions in each case, for the numerical integration of the three angles and their time derivatives. Numerical experiments were performed for the purpose of investigating the expected accuracies of the various parameters in the adjustment model under differing conditions. The results of the experiments indicate that the statistics for both the coordinates of points on earth and on the moon, and for the orientation parameters of the two bodies, are favorable.

Mueller, I. I.↗
NASA NTRS1973

Prior-to-failure extension of flaws under monotonic and pulsating loadings.

An equation governing the prior-to-failure crack propagation is proposed. For a rate-sensitive solid containing a two-dimensional crack and subject to the tensile mode of fracture the differential equations are integrated numerically for loads increasing monotonically in time. The resulting integral curves of load vs crack length and length vs time indicate that the growth of cracks in the subcritical range is strongly rate-dependent. The fatigue growth, viewed as a sequence of slow growth periods, is simulated on an EAI 380 analog computer. The fourth-power law proposed by Paris is confirmed only within certain range of high-cycle fatigue propagation and for a rate-insensitive solid. Otherwise, i.e., for a more pronounced rate dependence induced by viscosity of a solid and/or in the proximity of the final instability point, the growth is markedly enhanced.

Wnuk, M. P.↗
NASA NTRS1981

Sound diffraction around screens and wedges for arbitrary point source locations

A development is presented for the Green's function for a point source in the vicinity of a rigid wedge. The diffraction contributions to the Green's function for arbitrary source and listener location is expressed in a form which can be readily evaluated using the Laguerre technique for numerical integration. The present approach offers the advantages of efficient numerical evaluation and of relatively straightforward reduction to well-known analytical approximations in limiting cases. Comparisons with previously obtained experimental and numerical results obtained by Ambaud and Bergassoli (1972) are presented. The comparison with the experimental results is excellent; the advantages of the present numerical technique, vis a vis that of Ambaud and Bergassoli, are pointed out.

Hadden, W. J., Jr.↗
NASA NTRS1990

Downdating a time-varying square root information filter

A new method to efficiently downdate an estimate and covariance generated by a discrete time Square Root Information Filter (SRIF) is presented. The method combines the QR factor downdating algorithm of Gill and the decentralized SRIF algorithm of Bierman. Efficient removal of either measurements or a priori information is possible without loss of numerical integrity. Moreover, the method includes features for detecting potential numerical degradation. Performance on a 300 parameter system with 5800 data points shows that the method can be used in real time and hence is a promising tool for interactive data analysis. Additionally, updating a time-varying SRIF filter with either additional measurements or a priori information proceeds analogously.

Muellerschoen, Ronald J.↗
NASA NTRS1991

Solution of the Navier-Stokes equations for a driven cavity

The flow field in a lid driven cavity is determined by integration of the incompressible Navier-Stokes equations. The numerical integration is accomplished via an operator splitting method known as the theta-scheme. This splitting separates the problem into the solution of a quasi-stokes problem and a nonlinear convection problem. Some details of solution methods used for the two subproblems and results obtained for the driven cavity are described. The schemes developed for the quasi-Stokes problem are more advanced at this stage than those for the nonlinear problem. However, the approaches used for both parts are outlined. As a model problem, a two dimensional square cavity with sides of unit length and a lid moving with unit velocity from left to right is considered. The Navier-Stokes equations are discretized in space on a uniform staggered or MAC mesh. The time discretization is accomplished via the theta-scheme.

Semeraro, B. D.↗
NASA NTRS2002

Designing Adaptive Low Dissipative High Order Schemes

Proper control of the numerical dissipation/filter to accurately resolve all relevant multiscales of complex flow problems while still maintaining nonlinear stability and efficiency for long-time numerical integrations poses a great challenge to the design of numerical methods. The required type and amount of numerical dissipation/filter are not only physical problem dependent, but also vary from one flow region to another. This is particularly true for unsteady high-speed shock/shear/boundary-layer/turbulence/acoustics interactions and/or combustion problems since the dynamics of the nonlinear effect of these flows are not well-understood. Even with extensive grid refinement, it is of paramount importance to have proper control on the type and amount of numerical dissipation/filter in regions where it is needed.

Yee, H. C.↗