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 289 records · Page 16

Prior-to-failure extension of flaws under monotonic and pulsating loadings: Inelastic fatigue

An equation governing the prior to failure crack propagation is proposed. For a rate-sensitive solid containing two-dimensional crack and subject to the tensile mode of fracture, the differential equations are integrated numerically for the loads increasing monotonically in time. The resulting integral curves sigma = sigma(l) and l = l(t), i.e. 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 analogue computer. The fourth power law proposed by Paris is confirmed only within a certain range of high-cycle fatigue propagation and for a rate-insensitive solid. For a more pronounced rate dependency, induced by the viscosity of a solid and/or in the proximity of the final instability point, the growth is markedly enhanced.

Wnuk, M. P.↗

Some procedures for computerized electronic data processing of absorption measurements from artificial earth satellites

The processing of data obtained from solar absorption radiation measurements is discussed. The position of the satellite was obtained by numerical integration of the differential equations of motion using initial conditions. The position of the sun was calculated as a function of time, and the tangential elevation was determined approximately from the positions of the satellite and the sun. The coefficients of an approximation formula and of a data smoothing process were determined, and the inversion of an Abel integral equation is solved analytically.

Trinkkeller, B.↗

Historical evolution of vortex-lattice methods

A review of the beginning and some orientation of the vortex-lattice method were given. The historical course of this method was followed in conjunction with its field of computational fluid dynamics, spanning the period from L.F. Richardson's paper in 1910 to 1975. The following landmarks were pointed out: numerical analysis of partial differential equations, lifting-line theory, finite-difference method, 1/4-3/4 rule, block relaxation technique, application of electronic computers, and advanced panel methods.

Deyoung, J.↗

A formulation consideration for orthogonal collocation procedures

Numerical instabilities often arise in the use of high-ordered collocation approximations for numerically solving parabolic partial differential equations. These problems may be reduced by formulations involving evaluation of collocation polynomials rather than combination of the polynomials into a power series. As an illustration, two formulations using shifted Legendre polynomials of order 26 and less are compared.

Lashmet, P. K.↗

A three-dimensional, time-dependent model of Mobile Bay

A three-dimensional, time-variant mathematical model for momentum and mass transport in estuaries was developed and its solution implemented on a digital computer. The mathematical model is based on state and conservation equations applied to turbulent flow of a two-component, incompressible fluid having a free surface. Thus, bouyancy effects caused by density differences between the fresh and salt water, inertia from thare river and tidal currents, and differences in hydrostatic head are taken into account. The conservation equations, which are partial differential equations, are solved numerically by an explicit, one-step finite difference scheme and the solutions displayed numerically and graphically. To test the validity of the model, a specific estuary for which scaled model and experimental field data are available, Mobile Bay, was simulated. Comparisons of velocity, salinity and water level data show that the model is valid and a viable means of simulating the hydrodynamics and mass transport in non-idealized estuaries.

Pitts, F. H.↗

Application of the method of integral relations to laminar boundary layers in three dimensions

The method of integral relations is extended to general three-dimensional compressible laminar boundary layer flows. The transformation employed to transform the basic three-dimensional compressible boundary layer equations into quasi-incompressible form is an extension of the Howarth transformation. The resulting system of differential equations is integrated numerically by the method of integral relations as proposed by Dorodnitsyn. To demonstrate the accuracy of the method, it is applied to calculation of the parabolic flow over a flat plate and the boundary flow over an infinite yawed cylinder, for which solutions are known. It is then applied to the flow over a flat plate disturbed by a cylinder normal to the plate, for which a finite-difference solution is available for comparison. It is finally applied to calculating the crossflow velocity variation for supersonic flow over a swept wedge.

Holt, M.↗

Body-fitted coordinates systems transformations

Two computer programs generate two-dimensional body-fitted coordinate systems and coordinate transformation. Programs are useful in fields requiring accurate numerical representation of boundary conditions and accurate numerical solutions of partial differential equations.

Mastin, C. W.↗

Axially grooved heat pipe study

A technology evaluation study on axially grooved heat pipes is presented. The state-of-the-art is reviewed and present and future requirements are identified. Analytical models, the Groove Analysis Program (GAP) and a closed form solution, were developed to facilitate parametric performance evaluations. GAP provides a numerical solution of the differential equations which govern the hydrodynamic flow. The model accounts for liquid recession, liquid/vapor shear interaction, puddle flow as well as laminar and turbulent vapor flow conditions. The closed form solution was developed to reduce computation time and complexity in parametric evaluations. It is applicable to laminar and ideal charge conditions, liquid/vapor shear interaction, and an empirical liquid flow factor which accounts for groove geometry and liquid recession effects. The validity of the closed form solution is verified by comparison with GAP predictions and measured data.

Source record↗

An element formulation for perturbed motion about the center of mass

The perturbed motion of a rigid body about its center of mass, is formulated in terms of the six elements: l, the magnitude of the angular momentum vector; h, the total energy; delta and epsilon, two linear functions of the independent variable; and psi(1) and theta (1), two Euler angles that orientate the inertial frame with respect to the unperturbed solution. Solutions from the element formulation and the original Euler equations are numerically compared using shuttle-type data. For applied torques smaller than a given magnitude, the element formulation produced the following results: (1) larger step sizes in the numerical integration of the differential equations, resulting in an overall computational time-saving, and (2) more significant figures of accuracy in the computation of the variables describing the state of the rigid body.

Donaldson, J. D.↗

Calculation of CO concentration for liquid fueled gas turbine combustor

The extensive computation time required for the numerical integration of the differential equations obtained in studies of the concentrations of pollutants emitted by gas turbine combustors, can be reduced significantly by assuming the heat releasing hydrocarbon reactions to be in local equilibrium. In determining the CO and NO concentrations produced in spray combustion, it is, therefore, tempting to use the same local equilibrium assumption in order to simplify the kinetic calculations. An investigation of the validity of the local equilibrium assumption is conducted by using a simple analytical model, and then by actually carrying out the kinetic and local equilibrium calculations for typical case.

Patil, P. B.↗

A real-time subsynoptic, meso, and microscale severe storm forecast system

A system of differential equations is integrated numerically in space and time over several different matrices in an effort to simulate the atmospheric wave structures which organize severe local storms. Preliminary results with case studies indicate that dynamical fields produced by the numerical simulations can be translated into very fine scale space and time zones where severe storm forecast indices can be developed. These fine scale indices are now available in real-time when run on Langley's STAR 100 computer system.

Kaplan, M. L.↗

Numerical solution of potential flow about arbitrary 2-dimensional multiple bodies

A procedure for the finite-difference numerical solution of the lifting potential flow about any number of arbitrarily shaped bodies is given. The solution is based on a technique of automatic numerical generation of a curvilinear coordinate system having coordinate lines coincident with the contours of all bodies in the field, regardless of their shapes and number. The effects of all numerical parameters involved are analyzed and appropriate values are recommended. Comparisons with analytic solutions for single Karman-Trefftz airfoils and a circular cylinder pair show excellent agreement. The technique of application of the boundary-fitted coordinate systems to the numerical solution of partial differential equations is illustrated.

Thompson, J. F.↗

The turbulent recirculating flow field in a coreless induction furnace. A comparison of theoretical predictions with measurements

A mathematical representation for the electromagnetic force field and the fluid flow field in a coreless induction furnace is presented. The fluid flow field was represented by writing the axisymmetric turbulent Navier-Stokes equation, containing the electromagnetic body force term. The electromagnetic body force field was calculated by using a technique of mutual inductances. The kappa-epsilon model was employed for evaluating the turbulent viscosity and the resultant differential equations were solved numerically. Theoretically predicted velocity fields are in reasonably good agreement with the experimental measurements reported by Hunt and Moore; furthermore, the agreement regarding the turbulent intensities are essentially quantitative. These results indicate that the kappa-epsilon model provides a good engineering representation of the turbulent recirculating flows occurring in induction furnaces. At this stage it is not clear whether the discrepancies between measurements and the predictions, which were not very great in any case, are attributable either to the model or to the measurement techniques employed.

El-Kaddah, N.↗

Survey of new vector computers: The CRAY 1S from CRAY research; the CYBER 205 from CDC and the parallel computer from ICL - architecture and programming

Problems which can arise with vector and parallel computers are discussed in a user oriented context. Emphasis is placed on the algorithms used and the programming techniques adopted. Three recently developed supercomputers are examined and typical application examples are given in CRAY FORTRAN, CYBER 205 FORTRAN and DAP (distributed array processor) FORTRAN. The systems performance is compared. The addition of parts of two N x N arrays is considered. The influence of the architecture on the algorithms and programming language is demonstrated. Numerical analysis of magnetohydrodynamic differential equations by an explicit difference method is illustrated, showing very good results for all three systems. The prognosis for supercomputer development is assessed.

Gentzsch, W.↗

Grid spacing control with variation diminishing splines

Methods used to specify and control two and three dimensional grids on which numerical solutions of partial differential equations may be obtained were studied. While initially focusing on grid generation, the research will evolve into a consideration of the interaction of grid generation with the solution of a partial differential equation. The multisurface method of grid generation was used to continuously patch a grid onto an existing grid. In the resulting grid the elements of the Jacobian matrix must be continuous across the boundary between the original grid and the patched grid. Programs were written which accept as input the coordinates of the original grid and the desired new boundary, and then use the three-surface or four-surface version of the multisurface method to extend the original grid out to a new boundary.

Smith, P. W.↗

General curvilinear coordinate systems

The basic ideas of the construction and use of numerically-generated boundary-fitted coordinate systems for the numerical solution of partial differential equations are discussed. With such coordinate systems, all computation can be done on a fixed square grid in the rectangular transformed region regardless of the shape or movement of the physical boundaries. A number of different types of configurations for the transformed region and the basic transformation relations from a cartesian system to a general curvilinear system are given. The material of this paper is applicable to all types of coordinate system generation.

Thompson, J. P.↗

Algebraic grid generation

The numerical solution of partial differential equations about irregular geometries and with varying characteristic scales has created the need for coordinate systems and associated transformations which reflect both geometric and physical requirements. The process of finding coordinate transformations in discrete representations is called 'grid generation'. The present investigation is concerned with three algebraic grid generation methods. The methods include transfinite interpolation, the multisurface method, and the two-boundary technique. Interpolation formulas, in terms of homotopic mappings and constraints in terms of point positions and/or derivatives are the essential elements of the techniques. The methods are relatively simple to understand, they are explicit and do not require extensive computational effort, and they have a high degree of generality.

Smith, R. E.↗

The use of solution adaptive grids in solving partial differential equations

The grid point distribution used in solving a partial differential equation using a numerical method has a substantial influence on the quality of the solution. An adaptive grid which adjusts as the solution changes provides the best results when the number of grid points available for use during the calculation is fixed. Basic concepts used in generating and applying adaptive grids are reviewed in this paper, and examples illustrating applications of these concepts are presented.

Anderson, D. A.↗