A method for numerically calculating the area and distribution of water impingement on the leading edge of an airfoil in a cloud
Explore the source record for details and available documents.
Bergrun, N. R.↗
SEARCH · Search NASA
Results for “Numerical Methods”
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.
Recent advances in methods for numerical solution of O.D.E. initial value problems
In the mathematical modeling of physical systems, it is often necessary to solve an initial value problem (IVP), consisting of a system of ordinary differential equations (ODE). A typical program produces approximate solutions at certain mesh points. Almost all existing codes try to control the local truncation error, while the user is really interested in controlling the true or global error. The present investigation provides a review of recent advances regarding the solution of the IVP, giving particular attention to stiff systems. Stiff phenomena are customarily defined in terms of the eigenvalues of the Jacobian. There are, however, some difficulties connected with this approach. It is pointed out that an estimate of the Lipschitz constant proves to be a very practical way to determine the stiffness of a problem.
Bui, T. D.↗
Kinetic theory based new upwind methods for inviscid compressible flows
Two new upwind methods called the Kinetic Numerical Method (KNM) and the Kinetic Flux Vector Splitting (KFVS) method for the solution of the Euler equations have been presented. Both of these methods can be regarded as some suitable moments of an upwind scheme for the solution of the Boltzmann equation provided the distribution function is Maxwellian. This moment-method strategy leads to a unification of the Riemann approach and the pseudo-particle approach used earlier in the development of upwind methods for the Euler equations. A very important aspect of the moment-method strategy is that the new upwind methods satisfy the entropy condition because of the Boltzmann H-Theorem and suggest a possible way of extending the Total Variation Diminishing (TVD) principle within the framework of the H-Theorem. The ability of these methods in obtaining accurate wiggle-free solution is demonstrated by applying them to two test problems.
Deshpande, S. M.↗
A spectral collocation method for compressible, non-similar boundary layers
An efficient and highly accurate algorithm based on a spectral collocation method is developed for numerical solution of the compressible, two-dimensional and axisymmetric boundary layer equations. The numerical method incorporates a fifth-order, fully implicit marching scheme in the streamwise (timelike) dimension and a spectral collocation method based on Chebyshev polynomial expansions in the wall-normal (spacelike) dimension. The spectral collocation algorithm is used to derive the nonsimilar mean velocity and temperature profiles in the boundary layer of a 'fuselage' (cylinder) in a high-speed (Mach 5) flow parallel to its axis. The stability of the flow is shown to be sensitive to the gradual streamwise evolution of the mean flow and it is concluded that the effects of transverse curvature on stability should not be ignored routinely.
Pruett, C. D.↗
Gravitational Biology Facility on Space Station: Meeting the needs of space biology
The Gravitational Biology Facility (GBF) is a set of generic laboratory equipment needed to conduct research on Space Station Freedom (SSF), focusing on Space Biology Program science (Cell and Developmental Biology and Plant Biology). The GBF will be functional from the earliest utilization flights through the permanent manned phase. Gravitational biology research will also make use of other Life Sciences equipment on the space station as well as existing equipment developed for the space shuttle. The facility equipment will be developed based on requirements derived from experiments proposed by the scientific community to address critical questions in the Space Biology Program. This requires that the facility have the ability to house a wide variety of species, various methods of observation, and numerous methods of sample collection, preservation, and storage. The selection of the equipment will be done by the members of a scientific working group (5 members representing cell biology, 6 developmental biology, and 6 plant biology) who also provide requirements to design engineers to ensure that the equipment will meet scientific needs. All equipment will undergo extensive ground based experimental validation studies by various investigators addressing a variety of experimental questions. Equipment will be designed to be adaptable to other space platforms. The theme of the Gravitational Biology Facility effort is to provide optimal and reliable equipment to answer the critical questions in Space Biology as to the effects of gravity on living systems.
Allen, Katherine↗
A fast method of numerical quadrature for p-version finite element matrices
A new technique of numerical quadrature especially suited for p-version finite element matrices is presented. This new technique separates the integrand into two parts, and numerically operates on each part separately. The objective of this scheme is to minimize the computational cost of integrating the entire element matrix as opposed to minimizing the cost of integrating a single function. The efficiency of the new technique is compared with Gaussian quadrature and found to take a small fraction of the computational effort.
Hinnant, Howard E.↗
A method for numerically calculating the area and distribution of water impingement on the leading edge of an airfoil in a cloud
Explore the source record for details and available documents.
VAPOR, WATER↗
Galaxy Flows within 8000 km s −1 from Numerical Action Methods
Explore the source record for details and available documents.
Edward J. Shaya↗
Iterative methods for design sensitivity analysis
A numerical method is presented for design sensitivity analysis, using an iterative-method reanalysis of the structure generated by a small perturbation in the design variable; a forward-difference scheme is then employed to obtain the approximate sensitivity. Algorithms are developed for displacement and stress sensitivity, as well as for eignevalues and eigenvector sensitivity, and the iterative schemes are modified so that the coefficient matrices are constant and therefore decomposed only once.
Belegundu, A. D.↗
Some notes on the numerical solution of shear-lag and mathematically related problems
The analysis of box beams with shear deformation of the flanges can be reduced to the solution of a differential equation. The same equation is met in other problems of stress analysis. No analytical solutions of this equation can be given for practical cases, and numerical methods of evaluation must be used. Available methods are briefly discussed. Two numerical examples show the application of the step-by-step method of integration to shear-lag problems.
Kuhn, Paul↗
The Thin Oil Film Equation
A thin film of oil on a surface responds primarily to the wall shear stress generated on that surface by a three-dimensional flow. The oil film is also subject to wall pressure gradients, surface tension effects and gravity. The partial differential equation governing the oil film flow is shown to be related to Burgers' equation. Analytical and numerical methods for solving the thin oil film equation are presented. A direct numerical solver is developed where the wall shear stress variation on the surface is known and which solves for the oil film thickness spatial and time variation on the surface. An inverse numerical solver is also developed where the oil film thickness spatial variation over the surface at two discrete times is known and which solves for the wall shear stress variation over the test surface. A One-Time-Level inverse solver is also demonstrated. The inverse numerical solver provides a mathematically rigorous basis for an improved form of a wall shear stress instrument suitable for application to complex three-dimensional flows. To demonstrate the complexity of flows for which these oil film methods are now suitable, extensive examination is accomplished for these analytical and numerical methods as applied to a thin oil film in the vicinity of a three-dimensional saddle of separation.
Brown, James L.↗
Comparison of linear inversion methods by examination of the duality between iterative and inverse matrix methods
Linear numerical inversion methods applied to atmospheric remote sounding generally can be categorized in two ways: (1) iterative, and (2) inverse matrix methods. However, these two categories are not unrelated; a duality exists between them. In other words, given an iterative scheme, a corresponding inverse matrix method exists, and conversely. This duality concept is developed for the more familiar linear methods. The iterative duals are compared with the classical linear iterative approaches and their differences analyzed. The importance of the initial profile in all methods is stressed. Calculations using simulated data are made to compare accuracies and to examine the dependence of the solution on the initial profile.
Fleming, H. E.↗
The method of lines in analyzing solids containing cracks
A semi-numerical method is reviewed for solving a set of coupled partial differential equations subject to mixed and possibly coupled boundary conditions. The line method of analysis is applied to the Navier-Cauchy equations of elastic and elastoplastic equilibrium to calculate the displacement distributions in various, simple geometry bodies containing cracks. The application of this method to the appropriate field equations leads to coupled sets of simultaneous ordinary differential equations whose solutions are obtained along sets of lines in a discretized region. When decoupling of the equations and their boundary conditions is not possible, the use of a successive approximation procedure permits the analytical solution of the resulting ordinary differential equations. The use of this method is illustrated by reviewing and presenting selected solutions of mixed boundary value problems in three dimensional fracture mechanics. These solutions are of great importance in fracture toughness testing, where accurate stress and displacement distributions are required for the calculation of certain fracture parameters. Computations obtained for typical flawed specimens include that for elastic as well as elastoplastic response. Problems in both Cartesian and cylindrical coordinate systems are included. Results are summarized for a finite geometry rectangular bar with a central through-the-thickness or rectangular surface crack under remote uniaxial tension. In addition, stress and displacement distributions are reviewed for finite circular bars with embedded penny-shaped cracks, and rods with external annular or ring cracks under opening mode tension. The results obtained show that the method of lines presents a systematic approach to the solution of some three-dimensional mechanics problems with arbitrary boundary conditions. The advantage of this method over other numerical solutions is that good results are obtained even from the use of a relatively coarse grid.
Gyekenyesi, John P.↗
Comparison of characteristic schemes for three-dimensional, steady, isentropic flow.
The accuracy and order of three bicharacteristic schemes based on the numerical method of characteristics for three-dimensional, steady, isentropic, supersonic flow were compared based on numerical calculations. Comparisons were made for source flows and for flows in contoured axisymmetric nozzles. The first two methods were tetrahedral line networks, differing only in the manner in which the cross derivatives were evaluated. The third method was the pentahedral line network proposed by Butler. All three methods were numerically stable and produced efficient and accurate solutions. A discussion of the relative merits and disadvantages of each method is presented.
Cline, M. C.↗
The Finite Analytic Method for steady and unsteady heat transfer problems
A new numerical method called the Finite Analytical Method for solving partial differential equations is introduced. The basic idea of the finite analytic method is the incorporation of the local analytic solution in obtaining the numerical solution of the problem. The finite analytical method first divides the total region of the problem into small subregions in which local analytic solutions are obtained. Then an algebraic equation is derived from the local analytic solution for each subregion relating an interior nodal value at a point P in the subregion to its neighboring nodal values. The assembly of all the local analytic solutions thus provides the finite-analytic numerical solution of the problem. In this paper the finite analytic method is illustrated in solving steady and unsteady heat transfer problems.
Chen, C.-J.↗
Motion and interaction of decaying trailing vortices in spanwise shear wind
A simulation is presented of the drift of trailing vortices in a cross-wind near the ground by an unsteady, two-dimensional, rotational flow field with a concentration of large vorticity in vortical spots (having a finite but small effective size and finite total strength). The problem is analyzed by a combination of the method of matched asymptotic analyses for the decay of the vortical spots and the Euler solution for the unsteady rotational flow. Using the method of averaging, a special numerical method is developed in which the grid size and time step depend only on the length and velocity scales of the background flow and are independent of the effective core size of a vortical spot. The core size can be much smaller than the grid size, whereas the peak velocity in the core is inversely propertional to the spot size. Numerical results are presented to demonstrate the strong interaction between the trajectories of the vortical spots and the change of the vorticity distribution in the background flow field.
Liu, C. H.↗
Improved vortex methods for three-dimensional flows
Robust numerical methods are developed for three-dimensional incompressible vortical flows, using Lagrangian vortex elements. A successful scheme must be able to handle regions of intense vortex stretching and vortex reconnection with reasonable accuracy (without diverging). Here, consideration is given to vortex particles, also commonly called vortons or vortex sticks. The following issues are discussed: (1) use of delta-function elements and weak solutions of the vorticity equation; (2) use of smoothed elements and the choice of the smoothing function; (3) representation of viscous effects and the redistribution of element strength; and (4) conservation laws (are they satisfied?). The various proposed schemes have been tested on flows involving a strong interaction between two vortex rings.
Winckelmans, G.↗
A Numerical Analysis of the Transient Response of an Ablation System Including Effects of Thermal Nonequilibrium, Mass Transfer and Chemical Kinetics
The differential equations governing the transient response of a one-dimensional ablative thermal protection system undergoing stagnation ablation are derived. These equations are for thermal nonequilibrium effects between the pyrolysis gases and the char layer and kinetically controlled chemical reactions and mass transfer between the pyrolysis gases and the char layer. The boundary conditions are written for the particular case of stagnation heating with surface removal by oxidation or sublimation and pyrolysis of the uncharred layer occurring in a plane. The governing equations and boundary conditions are solved numerically using the modified implicit method (Crank-Nicolson method). Numerical results are compared with exact solutions for a number of simplified cases. The comparison is favorable in each instance.
Clark, R. K.↗