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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

An efficient numerical algorithm for transverse impact problems

Transverse impact problems in which the elastic and plastic indentation effects are considered, involve a nonlinear integral equation for the contact force, which, in practice, is usually solved by an iterative scheme with small increments in time. In this paper, a numerical method is proposed wherein the iterations of the nonlinear problem are separated from the structural response computations. This makes the numerical procedures much simpler and also efficient. The proposed method is applied to some impact problems for which solutions are available, and they are found to be in good agreement. The effect of the magnitude of time increment on the results is also discussed.

Sankar, B. V.↗
NASA NTRS1982

Reliability of the finite element method for calculating free edge stresses in composite laminates

The interlaminar normal stress distributions along the interface between the +45 deg and -45 deg plies of a graphite/epoxy laminate, obtained by various investigators, were found to disagree in both magnitude and sign. The reliability of the displacement-formulated finite element method in analyzing the edge-stress problem of a composite laminate is investigated. The history of the edge-stress problem is reviewed, and two well-known elasticity problems, one involving a stress discontinuity and one a singularity, are analyzed. The finite element analysis in these problems yields accurate stress distributions everywhere except in two elements closest to the stress discontinuity or singularity. Stress distributions for a + or -45 deg ply laminate near the singularity were similar to those of the two elasticity problems, demonstrating the methods, accuracy for calculating interlaminar stresses in composite laminates. The disagreement between the numerical methods was attributed to the unsymmetric stress tensor at singularity.

Whitcomb, J. D.↗
NASA NTRS1987

A method of calculating the total flow from a given sea surface topography

Using a simple dynamical model of a wind-driven ocean circulation of the Stommel type, and an analytical basis developed to objectively analyze the sea surface height residuals from an altimeter and, in the process, to determine the total flow instead of just the near surface geostrophic component associated with the given sea surface topography. The method is based on first deriving the solution to the forced problem for a given wind stress required to develop a hypothetical true or perfect data field and to establishing the basis for the objective analysis. The stream function and the surface height field for the forced problem are developed in terms of certain characteristic functions with the same expansion coefficients for both fields. These characteristic functions are simply the solutions for a homogeneous elliptic equation for the stream function and the solutions of an inhomogeneous balance equation for the height field. For the objective analysis, using a sample of randomly selected height values from the true data field, the height field characteristic functions are used to fit the given topography in a least squares sense. The resulting expansion coefficients then permit the synthesis of the total flow field via the stream function characteristic modes and the solution is perfectly well behaved even along the equator. The method of solution is easily adaptable to realistic ocean basis by straight forward numerical methods. The analytical basis of the theory and the results for an ideal rectangular basin on a beta plane are described.

Rao, Desiraju B.↗
NASA NTRS1990

Application of the probabilistic approximate analysis method to a turbopump blade analysis

An eigenvalue analysis of a typical space propulsion system turbopump blade is presented using an approximate probabilistic analysis methodology. The methodology was developed originally to investigate the feasibility of computing probabilistic structural response using closed-form approximate models. This paper extends the methodology to structures for which simple closed-form solutions do not exist. The finite element method will be used for this demonstration, but the concepts apply to any numerical method. The results agree with detailed analysis results and indicate the usefulness of using a probabilistic approximate analysis in determining efficient solution strategies.

Thacker, B. H.↗
NASA NTRS1993

Virtual-pulse time integral methodology: A new explicit approach for computational dynamics - Theoretical developments for general nonlinear structural dynamics

The present paper describes a new explicit virtual-pulse time integral methodology for nonlinear structural dynamics problems. The purpose of the paper is to provide the theoretical basis of the methodology and to demonstrate applicability of the proposed formulations to nonlinear dynamic structures. Different from the existing numerical methods such as direct time integrations or mode superposition techniques, the proposed methodology offers new perspectives and methodology of development, and possesses several unique and attractive computational characteristics. The methodology is tested and compared with the implicit Newmark method (trapezoidal rule) through a nonlinear softening and hardening spring dynamic models. The numerical results indicate that the proposed explicit virtual-pulse time integral methodology is an excellent alternative for solving general nonlinear dynamic problems.

Chen, Xiaoqin↗
NASA NTRS1992

Optimal matrix approximants in structural identification

Problems of model correlation and system identification are central in the design, analysis, and control of large space structures. Of the numerous methods that have been proposed, many are based on finding minimal adjustments to a model matrix sufficient to introduce some desirable quality into that matrix. In this work, several of these methods are reviewed, placed in a modern framework, and linked to other previously known ideas in computational linear algebra and optimization. This new framework provides a point of departure for a number of new methods which are introduced here. Significant among these is a method for stiffness matrix adjustment which preserves the sparsity pattern of an original matrix, requires comparatively modest computational resources, and allows robust handling of noisy modal data. Numerical examples are included to illustrate the methods presented herein.

Beattie, C. A.↗
NASA NTRS1998

Spurious Numerical Oscillations in Numerical Simulation of Supersonic Flows Using Shock Capturing Schemes

The numerical simulation of transitional and turbulent processes in hypersonic boundary layers often involves a physical process of a shock-disturbance wave interaction in complex two-dimensional and three-dimensional flow fields. For such simulations, it is required that there be a high order of accuracy in capturing both the shock and the small disturbances. The purpose of this paper is to evaluate the viability of using high order shock capturing schemes to track small disturbances in a multi-dimensional steady hypersonic flow. The numerical methods that are to be studied are the Total Variation Diminishing (TVD) scheme, and Essentially Non-Oscillatory (ENO) scheme. This paper shows that the presence of numerical oscillations in the flow field solution may drastically hinder any attempt at tracking the propagation of any physical disturbances. It has been found that the numerical oscillations that exist for shock capturing methods may be significant enough to pollute a flow field containing small physical disturbances. The effects of the refinement of the grid do not reduce the oscillations, but rather they decrease the wavelength of the oscillations. It is shown that by aligning the shock with the grid, the amplitude of these spurious oscillations may be greatly reduced.

Lee, Theodore K.↗
NASA NTRS1993

Three-dimensional magnetotail equilibria by numerical relaxation techniques

Consideration is given to a numerical method to iteratively approach 3D magnetostatic force equilibria, with modeling emphasis on the development of a suitable model of the earth's magnetotail, including a portion of the inner magnetosphere, i.e., on models which violate the 'tail approximation' commonly employed in analytical models. The numerical approach is compared to methods developed for laboratory plasma physics. The method is applied to the magnetotail outside of 10 R(E), using Tsyganenko's (1987) model for the quiet magnetosphere as an initial condition. The changes of the magnetic field necessary to yield an equilibrium configuration and the resulting distribution of the self-consistently derived pressure are discussed. It is shown that a self-consistent magnetotail equilibrium based on a close approximation to Tsyganenko's requires a region 1 type current system, which is not present in the initial configuration.

Hesse, Michael↗
NASA NTRS1987

Parallel Algorithm Solves Coupled Differential Equations

Numerical methods adapted to concurrent processing. Algorithm solves set of coupled partial differential equations by numerical integration. Adapted to run on hypercube computer, algorithm separates problem into smaller problems solved concurrently. Increase in computing speed with concurrent processing over that achievable with conventional sequential processing appreciable, especially for large problems.

Hayashi, A.↗
NASA NTRS1990

A multiblock multigrid method for the solution of the three-dimensional Euler equations

A general multiblock, multigrid method for the solution of the Euler equations has been developed. Two types of numerical methods were investigated, van Leer's flux-vector-splitting and Roe's flux-difference-splitting, with MUSCL type differencing used in both methods. An explicit two-step method and a multi-stage Runge-Kutta method have been tested. Results are presented for test cases of a channel flow, nozzle exhaust flow, and a transonic wing.

Cannizzaro, Frank E.↗
NASA NTRS1989

An Approximate Model for the Performance and Acoustic Predictions of Counterrotating Propeller Configurations

An approximate method was developed to analyze and predict the acoustics of a counterrotating propeller configuration. The method employs the analytical techniques of Lock and Theodorsen as described by Davidson to predict the steady performance of a counterrotating configuration. Then, a modification of the method of Lesieutre is used to predict the unsteady forces on the blades. Finally, the steady and unsteady loads are used in the numerical method of Succi to predict the unsteady acoustics of the propeller. The numerical results are compared with experimental acoustic measurements of a counterrotating propeller configuration by Gazzaniga operating under several combinations of advance ratio, blade pitch, and number of blades. In addition, a constant-speed commuter-class propeller configuration was designed with the Davidson method and the acoustics analyzed at three advance ratios. Noise levels and frequency spectra were calculated at a number of locations around the configuration. The directivity patterns of the harmonics in both the horizontal and vertical planes were examined, with the conclusion that the noise levels of the even harmonics are relatively independent of direction whereas the noise levels of the odd harmonics are extremely dependent on azimuthal direction in the horizontal plane. The equations of Succi are examined to explain this behavior.

Denner, Brett William↗
NASA NTRS1992

Recovery of the Navier-Stokes equations using a lattice-gas Boltzmann method

A lattice Boltzmann model is presented which gives the complete Navier-Stokes equation and may provide an efficient parallel numerical method for solving various fluid problems. The model uses the single-time relaxation approximation and a particular Maxwell-type distribution. The model eliminates exactly (1) the non-Galilean invariance caused by a density-dependent coefficient in the convection term and (2) a velocity-dependent equation of state.

Chen, Hudong↗
DOE OSTI2025

Evaluation of Flow Routing on the Unstructured Voronoi Meshes in Earth System Modeling

Flow routing is a fundamental process of Earth System Models' (ESMs) river component. Traditional flow routing models rely on Cartesian rectangular meshes, which exhibit limitations, particularly when coupled with unstructured mesh-based ocean components. They also lack the support for regionally refined models. While previous studies have highlighted the potential benefits of unstructured meshes for flow routing, their widespread application and comprehensive evaluation within ESMs remain limited. This study extends the river component of the Energy Exascale Earth System Model to unstructured Voronoi meshes. We evaluated the model's performance in simulating river discharge and water depth across three watersheds spanning the Arctic, temperate, and tropical regions. The results show that while providing several benefits, unstructured mesh-based flow routing can achieve comparable performance to structured mesh-based routing, and their difference is often less than 10%. Although the unstructured mesh-based method could address several existing limitations, this research also shows that additional improvements in the numerical method are needed to fully exploit the advantages of unstructured mesh for hydrologic and ESMs.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1968

Concepts for a theoretical and experimental study of lifting rotor random loads and vibrations, Phase 2

A comparison with NASA conducted simulator studies has shown that the approximate digital method for computing rotor blade flapping responses to random inputs, tentatively suggested in Phase I Report, gives with increasing rotor advance ratio the wrong trend. Consequently, three alternative methods of solution have been considered and are described: (1) an approximate method based on the functional relation between input and output double frequency spectra, (2) a numerical method based on the system responses to deterministic inputs and (3) a perturbation approach. Among these the perturbation method requires the least amount of computation and has been developed in two forms - the first form to obtain the response correlation function and the second for the time averaged spectra of flapping oscillations.

Hohenemser, K. H.↗
NASA NTRS1978

Numerical spatial marching techniques for estimating duct attenuation and source pressure profiles

A numerical method was developed that could predict the pressure distribution of a ducted source from far field pressure inputs. Using an initial value formulation, the two-dimensional homogeneous Helmholtz wave equation (no steady flow) was solved using explicit marching techniques. The Von Neumann method was used to develop relationships which describe how sound frequency and grid spacing effect numerical stability. At the present time, stability considerations limit the approach to high frequency sound. Sample calculations for both hard and soft wall ducts compare favorably to known boundary value solutions. In addition, assuming that reflections in the duct are small, this initial value approach was successfully used to determine the attenuation of a straight soft wall duct. Compared to conventional finite difference or finite element boundary value approaches, the numerical marching technique is orders of magnitude shorter in computation time and required computer storage and can be easily employed in problems involving high frequency sound.

Baumeister, K. J.↗
NASA NTRS1978

Numerical spatial marching techniques for estimating duct attenuation and source pressure profiles

A numerical method is developed that could predict the pressure distribution of a ducted source from far-field pressure inputs. Using an initial value formulation, the two-dimensional homogeneous Helmholtz wave equation (no steady flow) is solved using explicit marching techniques. The Von Neumann method is used to develop relationships which describe how sound frequency and grid spacing effect numerical stability. At the present time, stability considerations limit the approach to high frequency sound. Sample calculations for both hard and soft wall ducts compare favorably to known boundary value solutions. In addition, assuming that reflections in the duct are small, this initial value approach is successfully used to determine the attenuation of a straight soft wall duct. Compared to conventional finite difference or finite element boundary value approaches, the numerical marching technique is orders of magnitude shorter in computation time and required computer storage and can be easily employed in problems involving high frequency sound.

Baumeister, K. J.↗