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

Calculation of far-field scattering from nonspherical particles using a geometrical optics approach

A numerical method was developed using geometrical optics to predict far-field optical scattering from particles that are symmetric about the optic axis. The diffractive component of scattering is calculated and combined with the reflective and refractive components to give the total scattering pattern. The phase terms of the scattered light are calculated as well. Verification of the method was achieved by assuming a spherical particle and comparing the results to Mie scattering theory. Agreement with the Mie theory was excellent in the forward-scattering direction. However, small-amplitude oscillations near the rainbow regions were not observed using the numerical method. Numerical data from spheroidal particles and hemispherical particles are also presented. The use of hemispherical particles as a calibration standard for intensity-type optical particle-sizing instruments is discussed.

Hovenac, Edward A.↗
NASA NTRS1974

A numerical-perturbation method for the nonlinear analysis of structural vibrations

A numerical-perturbation method is proposed for the determination of the nonlinear forced response of structural elements. Purely analytical techniques are capable of determining the response of structural elements having simple geometries and simple variations in thickness and properties, but they are not applicable to elements with complicated structure and boundaries. Numerical techniques are effective in determining the linear response of complicated structures, but they are not optimal for determining the nonlinear response of even simple elements when modal interactions take place due to the complicated nature of the response. Therefore, the optimum is a combined numerical and perturbation technique. The present technique is applied to beams with varying cross sections.

Nayfeh, A. H.↗
NASA NTRS1974

Multiple zeros of polynomials

For polynomials of higher degree, iterative numerical methods must be used. Four iterative methods are presented for approximating the zeros of a polynomial using a digital computer. Newton's method and Muller's method are two well known iterative methods which are presented. They extract the zeros of a polynomial by generating a sequence of approximations converging to each zero. However, both of these methods are very unstable when used on a polynomial which has multiple zeros. That is, either they fail to converge to some or all of the zeros, or they converge to very bad approximations of the polynomial's zeros. This material introduces two new methods, the greatest common divisor (G.C.D.) method and the repeated greatest common divisor (repeated G.C.D.) method, which are superior methods for numerically approximating the zeros of a polynomial having multiple zeros. These methods were programmed in FORTRAN 4 and comparisons in time and accuracy are given.

Wood, C. A.↗
NASA NTRS2002

Libration Orbit Mission Design: Applications of Numerical & Dynamical Methods

Sun-Earth libration point orbits serve as excellent locations for scientific investigations. These orbits are often selected to minimize environmental disturbances and maximize observing efficiency. Trajectory design in support of libration orbits is ever more challenging as more complex missions are envisioned in the next decade. Trajectory design software must be further enabled to incorporate better understanding of the libration orbit solution space and thus improve the efficiency and expand the capabilities of current approaches. The Goddard Space Flight Center (GSFC) is currently supporting multiple libration missions. This end-to-end support consists of mission operations, trajectory design, and control. It also includes algorithm and software development. The recently launched Microwave Anisotropy Probe (MAP) and upcoming James Webb Space Telescope (JWST) and Constellation-X missions are examples of the use of improved numerical methods for attaining constrained orbital parameters and controlling their dynamical evolution at the collinear libration points. This paper presents a history of libration point missions, a brief description of the numerical and dynamical design techniques including software used, and a sample of future GSFC mission designs.

Bauer, Frank↗
NASA NTRS1982

Random element method for numerical modeling of diffusional processes

The random element method is a generalization of the random vortex method that was developed for the numerical modeling of momentum transport processes as expressed in terms of the Navier-Stokes equations. The method is based on the concept that random walk, as exemplified by Brownian motion, is the stochastic manifestation of diffusional processes. The algorithm based on this method is grid-free and does not require the diffusion equation to be discritized over a mesh, it is thus devoid of numerical diffusion associated with finite difference methods. Moreover, the algorithm is self-adaptive in space and explicit in time, resulting in an improved numerical resolution of gradients as well as a simple and efficient computational procedure. The method is applied here to an assortment of problems of diffusion of momentum and energy in one-dimension as well as heat conduction in two-dimensions in order to assess its validity and accuracy. The numerical solutions obtained are found to be in good agreement with exact solution except for a statistical error introduced by using a finite number of elements, the error can be reduced by increasing the number of elements or by using ensemble averaging over a number of solutions.

Ghoniem, A. F.↗
NASA NTRS1993

A time-accurate finite volume method valid at all flow velocities

A finite volume method to solve the Navier-Stokes equations at all flow velocities (e.g., incompressible, subsonic, transonic, supersonic and hypersonic flows) is presented. The numerical method is based on a finite volume method that incorporates a pressure-staggered mesh and an incremental pressure equation for the conservation of mass. Comparison of three generally accepted time-advancing schemes, i.e., Simplified Marker-and-Cell (SMAC), Pressure-Implicit-Splitting of Operators (PISO), and Iterative-Time-Advancing (ITA) scheme, are made by solving a lid-driven polar cavity flow and self-sustained oscillatory flows over circular and square cylinders. Calculated results show that the ITA is the most stable numerically and yields the most accurate results. The SMAC is the most efficient computationally and is as stable as the ITA. It is shown that the PISO is the most weakly convergent and it exhibits an undesirable strong dependence on the time-step size. The degenerated numerical results obtained using the PISO are attributed to its second corrector step that cause the numerical results to deviate further from a divergence free velocity field. The accurate numerical results obtained using the ITA is attributed to its capability to resolve the nonlinearity of the Navier-Stokes equations. The present numerical method that incorporates the ITA is used to solve an unsteady transitional flow over an oscillating airfoil and a chemically reacting flow of hydrogen in a vitiated supersonic airstream. The turbulence fields in these flow cases are described using multiple-time-scale turbulence equations. For the unsteady transitional over an oscillating airfoil, the fluid flow is described using ensemble-averaged Navier-Stokes equations defined on the Lagrangian-Eulerian coordinates. It is shown that the numerical method successfully predicts the large dynamic stall vortex (DSV) and the trailing edge vortex (TEV) that are periodically generated by the oscillating airfoil. The calculated streaklines are in very good comparison with the experimentally obtained smoke picture. The calculated turbulent viscosity contours show that the transition from laminar to turbulent state and the relaminarization occur widely in space as well as in time. The ensemble-averaged velocity profiles are also in good agreement with the measured data and the good comparison indicates that the numerical method as well as the multipletime-scale turbulence equations successfully predict the unsteady transitional turbulence field. The chemical reactions for the hydrogen in the vitiated supersonic airstream are described using 9 chemical species and 48 reaction-steps. Consider that a fast chemistry can not be used to describe the fine details (such as the instability) of chemically reacting flows while a reduced chemical kinetics can not be used confidently due to the uncertainty contained in the reaction mechanisms. However, the use of a detailed finite rate chemistry may make it difficult to obtain a fully converged solution due to the coupling between the large number of flow, turbulence, and chemical equations. The numerical results obtained in the present study are in good agreement with the measured data. The good comparison is attributed to the numerical method that can yield strongly converged results for the reacting flow and to the use of the multiple-time-scale turbulence equations that can accurately describe the mixing of the fuel and the oxidant.

Kim, S.-W.↗
NASA NTRS1986

Numerical quadrature methods for integrals of singular periodic functions and their application to singular and weakly singular integral equations

High accuracy numerical quadrature methods for integrals of singular periodic functions are proposed. These methods are based on the appropriate Euler-Maclaurin expansions of trapezoidal rule approximations and their extrapolations. They are used to obtain accurate quadrature methods for the solution of singular and weakly singular Fredholm integral equations. Such periodic equations are used in the solution of planar elliptic boundary value problems, elasticity, potential theory, conformal mapping, boundary element methods, free surface flows, etc. The use of the quadrature methods is demonstrated with numerical examples.

Sidi, A.↗
NASA NTRS1973

Computational considerations in application of the finite element method for analysis of unsteady flow around airfoils.

Difficulties related to the application of the finite element method for numerical solution of the Navier-Stokes equation are discussed. The unsteady flow around a stationary and oscillating NACA 0012 airfoil was analyzed. In an effort to establish the range of applicability of the developed numerical method, the accuracy and stability problems in representing unsteady flow patterns and determining pressure distribution around the airfoil were investigated. Error involved in determining the velocity field and in the numerical integration of the discretized equations were analyzed in terms of geometry of the finite element gridwork, boundary conditions for stationary and oscillating airfoils, angles of attack and flow conditions.

Bratanow, T.↗
NASA NTRS1988

A review of various flow-solid interaction analysis methods with emphasis on recent advances in turbulence models and flow analysis methods

A review on various experimental, analytical, and numerical methods for the flow-solid interaction of a nest of cylinders subjected to cross flow is presented in this paper. In order to assess the possibility of initiating a development of major numerical analysis method for flow-solid interaction of a nest of cylinders in cross flows, a review on various turbulence models, numerical methods such as the finite element method and the finite difference method, and the available numerical methods for flow-solid interactions have been included in this paper. A multiple-time-scale turbulence model for complex turbulent flows and a finite element method for high Reynolds number flows have been introduced herein for the same reason.

Kim, S.-W.↗
NASA NTRS1979

Calculation of three-dimensional turbulent subsonic flows in transition ducts

A method for computing three-dimensional turbulent subsonic flow in curved ducts is being developed. A set of tube-like surface oriented coordinates is employed for a general class of geometries applicable to subsonic diffusers with offset bends. The geometric formulation is complex and no previous treatment of this class of viscous flow problems is known to the authors. The duct centerline is a space curve specified by piecewise polynomials. A Frenet frame is located on the centerline at each axial location. The cross sections are described by superellipses imbedded in the Frenet frame. Duct surfaces are also coordinate surfaces, which greatly simplifies the boundary conditions. The resulting coordinates are nonorthogonal. An approximate set of governing equations is employed for viscous flows having strong flow in a primary flow direction. The derivation is coordinate invariant and the resulting equations are expressed in tensor form. These equations are solved by an efficient alternating direction implicit (ADI) method. This numerical method is generally stable and permits solution in difficult geometries using the general tensor formulation.

Levy, R.↗
NASA NTRS1988

Analysis Of Flow Across Cylinders

Report reviews various experimental, theoretical, and numerical methods of analysis of interactions between fluid and one or more cylinders across which it flows. Begins with discussion of mechanisms inducing vibrations in nests of cylinders. Major portion of review devoted to methods for numerical analysis of flow/cylinder interactions.

Kim, Sang-Wook↗
NASA NTRS1999

Integrated Aeroservoelastic Optimization: Status and Direction

The interactions of lightweight flexible airframe structures, steady and unsteady aerodynamics, and wide-bandwidth active controls on modern airplanes lead to considerable multidisciplinary design challenges. More than 25 years of mathematical and numerical methods' development, numerous basic research studies, simulations and wind-tunnel tests of simple models, wind-tunnel tests of complex models of real airplanes, as well as flight tests of actively controlled airplanes, have all contributed to the accumulation of a substantial body of knowledge in the area of aeroservoelasticity. A number of analysis codes, with the capabilities to model real airplane systems under the assumptions of linearity, have been developed. Many tests have been conducted, and results were correlated with analytical predictions. A selective sample of references covering aeroservoelastic testing programs from the 1960s to the early 1980s, as well as more recent wind-tunnel test programs of real or realistic configurations, are included in the References section of this paper. An examination of references 20-29 will reveal that in the course of development (or later modification), of almost every modern airplane with a high authority active control system, there arose a need to face aeroservoelastic problems and aeroservoelastic design challenges.

Livne, Eli↗
NASA NTRS2005

Numerical Modeling of Flow Distribution in Micro-Fluidics Systems

This paper describes an application of a general purpose computer program, GFSSP (Generalized Fluid System Simulation Program) for calculating flow distribution in a network of micro-channels. GFSSP employs a finite volume formulation of mass and momentum conservation equations in a network consisting of nodes and branches. Mass conservation equation is solved for pressures at the nodes while the momentum conservation equation is solved at the branches to calculate flowrate. The system of equations describing the fluid network is solved by a numerical method that is a combination of the Newton-Raphson and successive substitution methods. The numerical results have been compared with test data and detailed CFD (computational Fluid Dynamics) calculations. The agreement between test data and predictions is satisfactory. The discrepancies between the predictions and test data can be attributed to the frictional correlation which does not include the effect of surface tension or electro-kinetic effect.

Majumdar, Alok↗
NASA NTRS1993

Analysis of Radiation-natural Convection Interactions in 1-g and low-g Environments using the Discrete Exchange Factor Method

A new numerical method is presented for the analysis of combined natural convection and radiation heat transfer with applications in many engineering situations such as materials processing, combustion and fire research. Because of the recent interest in the low gravity environment of space, attention is devoted to both 1-g and low-g applications. The two-dimensional mathematical model is represented by a set of coupled nonlinear integro-partial differential equations. Radiative exchange is formulated using the Discrete Exchange Factor method (DEF). This method considers point to point exchange and provides accurate results over a wide range of radiation parameters. Numerical results show that radiation significantly influences the flow and heat transfer in both low-g and 1-g applications. In the low-g environment, convection is weak, and radiation can easily become the dominant heat transfer mode. It is also shown that volumetric heating by radiation gives rise to an intricate cell pattern in the top heated enclosure.

Kassemi, M.↗
NASA NTRS1995

Implementation of a multiblock sensitivity analysis method in numerical aerodynamic shape optimization

A multiblock sensitivity analysis method is applied in a numerical aerodynamic shape optimization technique. The Sensitivity Analysis Domain Decomposition (SADD) scheme which is implemented in this study was developed to reduce the computer memory requirements resulting from the aerodynamic sensitivity analysis equations. Discrete sensitivity analysis offers the ability to compute quasi-analytical derivatives in a more efficient manner than traditional finite-difference methods, which tend to be computationally expensive and prone to inaccuracies. The direct optimization procedure couples CFD analysis based on the two-dimensional thin-layer Navier-Stokes equations with a gradient-based numerical optimization technique. The linking mechanism is the sensitivity equation derived from the CFD discretized flow equations, recast in adjoint form, and solved using direct matrix inversion techniques. This investigation is performed to demonstrate an aerodynamic shape optimization technique on a multiblock domain and its applicability to complex geometries. The objectives are accomplished by shape optimizing two aerodynamic configurations. First, the shape optimization of a transonic airfoil is performed to investigate the behavior of the method in highly nonlinear flows and the effect of different grid blocking strategies on the procedure. Secondly, shape optimization of a two-element configuration in subsonic flow is completed. Cases are presented for this configuration to demonstrate the effect of simultaneously reshaping interfering elements. The aerodynamic shape optimization is shown to produce supercritical type airfoils in the transonic flow from an initially symmetric airfoil. Multiblocking effects the path of optimization while providing similar results at the conclusion. Simultaneous reshaping of elements is shown to be more effective than individual element reshaping due to the inclusion of mutual interference effects.

Lacasse, James M.↗
NASA NTRS1974

Computation of steady axisymmetric flow using a one-dimensional time-dependent method

An iterative numerical method for computing steady, three dimensional, viscous, compressible flow fields, about aerodynamic bodies was studied. In order to develop the iterative method economically, the primary emphasis was directed towards supersonic, axisymmetric flow. However, the technique readily extends to three spatial dimensions. The viscous flow field about a cone-cylinder-flare body was calculated and compared to existing experimental data. Numerical predictions of the cone boundary layer and the flow field shock wave structure agreed with corresponding measurements. A separation was calculated at the cylinder-flare junction in six iterations; however, the size of the vortex did not correspond to the measured size. It was estimated that fifty iterations would be required to properly define the vortex. It was concluded that the iteration technique is of limited value for plane two dimensional and axisymmetrix flows, but of great value for three-dimensional flows.

Walitt, L.↗