Search NASA⌕ Search

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.

At least 469 records · Page 26

Two-dimensional compressible convection extending over multiple scale heights

The theoretical description of the dynamics of a stellar convection zone is considered, taking into account one of the most basic issues by studying compressible convection extending over multiple scale heights. A revised version of a code reported by Graham (1975) is employed. Two-dimensional simulations show that nonlinear compressible convection possesses cellular structures with strong localized downward-directed plumes and broader upflows. The horizontal flows which close the circulation within the cell satisfy an approximate Bernoulli integral along a considerable portion of the horizontal trajectory. Attention is given to details regarding the numerical methods, the properties of the numerical solutions, the overall effects of compressibility on nonlinear convection, and a comparison with anelastic modal solutions.

Hurlburt, N. E.↗

Transient natural convection heat and mass transfer in a rectangular enclosure - A numerical analysis

A numerical analysis of transient heat and solute transport across a rectangular cavity with combined horizontal temperature and concentration gradients is performed by a numerical method based on the SIMPLE. Numerical results show that the average Nusselt and Sherwood numbers both decrease markedly when the solutal and thermal buoyancy forces act in the opposite directions. When the solutal and thermal buoyancy forces act in the same directions, however, the average Sherwood number increases significantly and yet the average Nusselt number decreases slightly.

Han, Samuel S.↗

Dynamic force response of spherical hydrostatic journal bearing for cryogenic applications

Hydrostatic Journal Bearings (HJB's) are reliable and resilient fluid film rotor support elements ideal to replace roller bearings in cryogenic turbomachinery. HJB' will be used for primary space-power applications due to their long lifetime, low friction and wear, large load capacity, large direct stiffness, and damping force coefficients. An analysis for the performance characteristics of turbulent flow, orifice compensated, spherical hydrostatic journal bearings (HJB's) is presented. Spherical bearings allow tolerance for shaft misalignment without force performance degradation and have also the ability to support axial loads. The spherical HJB combines these advantages to provide a bearing design which could be used efficiently on high performance turbomachinery. The motion of a barotropic liquid on the thin film bearing lands is described by bulk-flow mass and momentum equations. These equations are solved numerically using an efficient CFD method. Numerical predictions of load capacity and force coefficients for a 6 recess, spherical HJB in a LO2 environment are presented. Fluid film axial forces and force coefficients of a magnitude about 20% of the radial load capacity are predicted for the case analyzed. Fluid inertia effects, advective and centrifugal, are found to affect greatly the static and dynamic force performance of the bearing studied.

Sanandres, Luis↗

A geometric deformable model for echocardiographic image segmentation

Gradient vector flow (GVF), an elegant external force for parametric deformable models, can capture object boundaries from both sides. A new geometric deformable model is proposed that combines GVF and the geodesic active contour model. The level set method is used as the numerical method of this model. The model is applied for echocardiographic image segmentation.

NASA Program Biomedical Research and Countermeasur↗

UCLA parallel PIC framework

The UCLA Parallel PIC Framework (UPIC) has been developed to provide trusted components for the rapid construction of new, parallel Particle-in-Cell (PIC) codes. The Framework uses object-based ideas in Fortran95, and is designed to provide support for various kinds of PIC codes on various kinds of hardware. The focus is on student programmers. The Framework supports multiple numerical methods, different physics approximations, different numerical optimizations and implementations for different hardware. It is designed with "defensive" programming in mind, meaning that it contains many error checks and debugging helps. Above all, it is designed to hide the complexity of parallel processing. It is currently being used in a number of new Parallel PIC codes.

Norton, Charles D.↗

A study of multiple-shaker modal survey testing

The principal objective was to examine and to assess the practical value of a method of multiple-shaker sinusoidal modal vibration testing known as Asher's method. Numerical studies which simulate the application of Asher's method and a unique experimental implementation of the method were completed. Another objective of the research was to develop and to demonstrate with numerical simulation a quantitative method for determining from transfer function data the number of dominant modes of vibration in sinusoidal structural response.

Hallauer, W. L., Jr.↗

SIAM Conference on Parallel Processing for Scientific Computing, 4th, Chicago, IL, Dec. 11-13, 1989, Proceedings

Attention is given to such topics as an evaluation of block algorithm variants in LAPACK and presents a large-grain parallel sparse system solver, a multiprocessor method for the solution of the generalized Eigenvalue problem on an interval, and a parallel QR algorithm for iterative subspace methods on the CM2. A discussion of numerical methods includes the topics of asynchronous numerical solutions of PDEs on parallel computers, parallel homotopy curve tracking on a hypercube, and solving Navier-Stokes equations on the Cedar Multi-Cluster system. A section on differential equations includes a discussion of a six-color procedure for the parallel solution of elliptic systems using the finite quadtree structure, data parallel algorithms for the finite element method, and domain decomposition methods in aerodynamics. Topics dealing with massively parallel computing include hypercube vs. 2-dimensional meshes and massively parallel computation of conservation laws. Performance and tools are also discussed.

Dongarra, Jack↗

Application of a fractional-step method to incompressible Navier-Stokes equation

A numerical method for computing three dimensional, time dependent incompressible flows is presented. The method is based on a fractional step, or time-splitting, scheme in conjunction with the approximate-factorization technique. The use of velocity boundary conditions for the intermediate velocity field leads to inconsistent numerical solutions. Appropriate boundary conditions for the intermediate velocity field are derived and tested. Numerical solutions for flow inside a driven cavity and over a backward-facing step are presented and compared with experimenal data and other numerical results.

Kim, J.↗

Application of a fractional-step method to incompressible Navier-Stokes equations

A numerical method for computing three dimensional, time dependent incompressible flows is presented. The method is based on a fractional step, or time-splitting, scheme in conjunction with the approximate-factorization technique. The use of velocity boundary conditions for the intermediate velocity field leads to inconsistent numerical solutions. Appropriate boundary conditions for the intermediate velocity field are derived and tested. Numerical solutions for flow inside a driven cavity and over a backward-facing step are presented and compared with experimental data and other numerical results.

Kim, J.↗

Development of finite analytic method for unsteady three-dimensional Navier-Stokes equation

A finite analytic (FA) numerical method is described for the numerical study of three-dimensional fluid flow problems. The FA method incorporates a local linearization scheme and a classical analytic solution technique to formulate the algebraic representation of unsteady three-dimensional Navier Stokes Equations. The twenty-eight point FA formula thus obtained gives all positive coefficients and desired upwind shift. Cubic cavity flows are then solved by the FA method for Reynolds numbers of 100 and 400. A simple variant of SIMPLER algorithm is used to resolve the pressure-velocity coupling problem. The numerical solutions show that the presence of side walls reduce the strength of the primary vortex when compared with the two-dimensional square cavity flows.

Chen, H.-C.↗

Spectral methods on arbitrary grids

Stable and spectrally accurate numerical methods are constructed on arbitrary grids for partial differential equations. These new methods are equivalent to conventional spectral methods but do not rely on specific grid distributions. Specifically, we show how to implement Legendre Galerkin, Legendre collocation, and Laguerre Galerkin methodology on arbitrary grids.

Carpenter, Mark H.↗

Combined Numerical/Analytical Perturbation Solutions of the Navier-Stokes Equations for Aerodynamic Ejector/Mixer Nozzle Flows

In spite of rapid advances in both scalar and parallel computational tools, the large number of variables involved in both design and inverse problems make the use of sophisticated fluid flow models impractical, With this restriction, it is concluded that an important family of methods for mathematical/computational development are reduced or approximate fluid flow models. In this study a combined perturbation/numerical modeling methodology is developed which provides a rigorously derived family of solutions. The mathematical model is computationally more efficient than classical boundary layer but provides important two-dimensional information not available using quasi-1-d approaches. An additional strength of the current methodology is its ability to locally predict static pressure fields in a manner analogous to more sophisticated parabolized Navier Stokes (PNS) formulations. To resolve singular behavior, the model utilizes classical analytical solution techniques. Hence, analytical methods have been combined with efficient numerical methods to yield an efficient hybrid fluid flow model. In particular, the main objective of this research has been to develop a system of analytical and numerical ejector/mixer nozzle models, which require minimal empirical input. A computer code, DREA Differential Reduced Ejector/mixer Analysis has been developed with the ability to run sufficiently fast so that it may be used either as a subroutine or called by an design optimization routine. Models are of direct use to the High Speed Civil Transport Program (a joint government/industry project seeking to develop an economically.viable U.S. commercial supersonic transport vehicle) and are currently being adopted by both NASA and industry. Experimental validation of these models is provided by comparison to results obtained from open literature and Limited Exclusive Right Distribution (LERD) sources, as well as dedicated experiments performed at Texas A&M. These experiments have been performed using a hydraulic/gas flow analog. Results of comparisons of DREA computations with experimental data, which include entrainment, thrust, and local profile information, are overall good. Computational time studies indicate that DREA provides considerably more information at a lower computational cost than contemporary ejector nozzle design models. Finally. physical limitations of the method, deviations from experimental data, potential improvements and alternative formulations are described. This report represents closure to the NASA Graduate Researchers Program. Versions of the DREA code and a user's guide may be obtained from the NASA Lewis Research Center.

DeChant, Lawrence Justin↗

Rendezvous Guidance Technology

The literature on rendezvous guidance technology is reviewed and a bibliography of about 160 documents is presented. The importance of the proper selection of the ascent trajectory and of the orbit for a manned space station or an orbiting launch complex is emphasized. The way the launch delay problem affects the selection of the ascent trajectory and the target orbit is discussed. A specific ascent trajectory and an orbit for a manned space station or an orbiting launch complex is recommended on the premise that launch delay problems will continue to outweigh most of the other orbit selection criteria. A new method for determining the velocity corrections required for mid-course guidance is proposed . The proposed methods consist of an "exact-numerical" solution of the relative equations of motion in the "Shell-coordinate system." The proposed "exact-numerical" method extends the applicability of the mid- course equations to a much greater range than that of the currently available linear methods .

Trajectory↗

An algorithm for calculation of the Jordan canonical form of a matrix

Jordan canonical forms are used extensively in the literature on control systems. However, very few methods are available to compute them numerically. Most numerical methods compute a set of basis vectors in terms of which the given matrix is diagonalized when such a change of basis is possible. Here, a simple and efficient method is suggested for computing the Jordan canonical form and the corresponding transformation matrix. The method is based on the definition of a generalized eigenvector, and a natural extension of Gauss elimination techniques.

Sridhar, B.↗

Aeroelastic stability of coupled flap-lag motion of hingeless helicopter blades at arbitrary advance ratios

Equations for large amplitude coupled flap-lag motion of a hingeless elastic helicopter blade in forward flight are derived. Only a torsionally rigid blade excited by quasi-steady aerodynamic loads is considered. The effects of reversed flow together with some new terms due to radial flow are included. Using Galerkin's method the spatial dependence is eliminated and the equations are linearized about a suitable equilibrium position. The resulting system of homogeneous periodic equations is solved using multivariable Floquet-Liapunov theory, and the transition matrix at the end of the period is evaluated by two separate methods. Computational efficiency of the two numerical methods is compared. Results illustrating the effects of forward flight and various important blade parameters on the stability boundaries are presented.

Friedmann, P.↗

A numerical study of transient heat and mass transfer in crystal growth

A numerical analysis of transient heat and solute transport across a rectangular cavity is performed. Five nonlinear partial differential equations which govern the conservation of mass, momentum, energy and solute concentration related to crystal growth in solution, are simultaneously integrated by a numerical method based on the SIMPLE algorithm. Numerical results showed that the flow, temperature and solute fields are dependent on thermal and solutal Grashoff number, Prandtl number, Schmidt number and aspect ratio. The average Nusselt and Sherwood numbers evaluated at the center of the cavity decrease markedly when the solutal buoyancy force acts in the opposite direction to the thermal buoyancy force. When the solutal and thermal buoyancy forces act in the same direction, however, Sherwood number increases significantly and yet Nusselt number decreases. Overall effects of convection on the crystal growth are seen to be an enhancement of growth rate as expected but with highly nonuniform spatial growth variations.

Han, Samuel Bang-Moo↗