Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical techniques”

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

Phase space simulation of collisionless stellar systems on the massively parallel processor

A numerical technique for solving the collisionless Boltzmann equation describing the time evolution of a self gravitating fluid in phase space was implemented on the Massively Parallel Processor (MPP). The code performs calculations for a two dimensional phase space grid (with one space and one velocity dimension). Some results from calculations are presented. The execution speed of the code is comparable to the speed of a single processor of a Cray-XMP. Advantages and disadvantages of the MPP architecture for this type of problem are discussed. The nearest neighbor connectivity of the MPP array does not pose a significant obstacle. Future MPP-like machines should have much more local memory and easier access to staging memory and disks in order to be effective for this type of problem.

White, Richard L.↗

An introduction to three-dimensional climate modeling

The development and use of three-dimensional computer models of the earth's climate are discussed. The processes and interactions of the atmosphere, oceans, and sea ice are examined. The basic theory of climate simulation which includes the fundamental equations, models, and numerical techniques for simulating the atmosphere, oceans, and sea ice is described. Simulated wind, temperature, precipitation, ocean current, and sea ice distribution data are presented and compared to observational data. The responses of the climate to various environmental changes, such as variations in solar output or increases in atmospheric carbon dioxide, are modeled. Future developments in climate modeling are considered. Information is also provided on the derivation of the energy equation, the finite difference barotropic forecast model, the spectral transform technique, and the finite difference shallow water waved equation model.

Washington, W. M.↗

Construction of explicit and implicit symmetric TVD schemes and their applications

A one-parameter family of explicit and implicit total variation diminishing (TVD) schemes is developed which permits incorporation of an expanded group of slope and flux limiters. The numerical technique is intended for use in calculations which include a time-differencing scheme and an optional Lax-Wendroff scheme. Methods of extending the TVD models to nonlinear scalar equations and systems of hyperbolic conservation equations are described. Sample results are presented from calculations of shocked flows around NACA 0012 and NACA 0018 airfoils.

Yee, H. C.↗

Cutoff wavenumbers and modes for annular-cross-section waveguide with eccentric inner conductor of small radius

Analytical expressions are derived for the cutoff wavenumbers and the corresponding modes in annular-cross-section waveguides having inner conductors of small radius. Waveguides with circular and rectangular outer boundary are considered. In the case of the circular eccentric annular waveguide, comparison is made between the values of cutoff wavenumbers computed from the expressions derived in this paper and data obtained by a more rigorous numerical technique.

Davidovitz, Marat↗

Numerical simulation of a supersonic reacting mixing layer

In order to arrive at physical models that can adequately describe supersonic combustion, and develop accurate and efficient numerical techniques for the solution of such models' governing equations, a computer program has been developed for the study of reacting flows which considers the multicomponent diffusion and convection of important chemical species, as well as their finite state reaction and the interaction of the fluid mechanics and the chemistry that occurred. The code employs a hybrid Chebyshev pseudospectral technique for integration of the models' resulting governing equations; the program is here used to study a spatially developing and reacting mixing layer.

Drummond, J. Philip↗

Fuel-air mixing and combustion in a two-dimensional Wankel engine

A two-equation turbulence model, an algebraic grid generalization method, and an approximate factorization time-linearized numerical technique are used to study the effects of mixture stratification at the intake port and gaseous fuel injection on the flow field and fuel-air mixing in a two-dimensional rotary engine model. The fuel distribution in the combustion chamber is found to be a function of the air-fuel mixture fluctuations at the intake port. It is shown that the fuel is advected by the flow field induced by the rotor and is concentrated near the leading apex during the intake stroke, while during compression, the fuel concentration is highest near the trailing apex and is lowest near the rotor. It is also found that the fuel concentration near the trailing apex and rotor is small except at high injection velocities.

Shih, T. I.-P.↗

Computing Radiation Characteristics Of Phased Arrays

Radiated powers, directivities, and radiation patterns calculated. Computer program devised to generalize results and obtain efficient numerical technique for computing directivity and radiation pattern of generalized array. Characteristics of generalized array used in program include arbitrary element locations; and polarizations and excitations of elements. Written in FORTRAN IV.

Acosta, Roberto J.↗

Dynamic response of laminated composite plates using a three-dimensional hybrid-stress finite-element formulation

A method of analysis of dynamic response of laminated composite plates is presented. The analysis is carried by using a hybrid-stress finite element numerical technique. By using this approach, the response of simply supported laminated plates subjected to sinusoidal loading are investigated. For the solution of the finite element equations of motion of free vibrations and dynamic response problems, two effective methods of solution, the space iteration method and the Newmark direct integration method are used. These two methods are discussed here.

Liou, W. J.↗

Review of fatigue and fracture research at NASA Langley Research Center

Most dynamic components in helicopters are designed with a safe-life constant-amplitude testing approach that has not changed in many years. In contrast, the fatigue methodology in other industries has advanced significantly in the last two decades. Recent research at the NASA Langley Research Center and the U.S. Army Aerostructures Directorate at Langley are reviewed relative to fatigue and fracture design methodology for metallic components. Most of the Langley research was directed towards the damage tolerance design approach, but some work was done that is applicable to the safe-life approach. In the areas of testing, damage tolerance concepts are concentrating on the small-crack effect in crack growth and measurement of crack opening stresses. Tests were conducted to determine the effects of a machining scratch on the fatigue life of a high strength steel. In the area of analysis, work was concentrated on developing a crack closure model that will predict fatigue life under spectrum loading for several different metal alloys including a high strength steel that is often used in the dynamic components of helicopters. Work is also continuing in developing a three-dimensional, finite-element stress analysis for cracked and uncracked isotropic and anisotropic structures. A numerical technique for solving simultaneous equations called the multigrid method is being pursued to enhance the solution schemes in both the finite-element analysis and the boundary element analysis. Finally, a fracture mechanics project involving an elastic-plastic finite element analysis of J-resistance curve is also being pursued.

Everett, Richard A., Jr.↗

Creep rupture analysis of a beam resting on high temperature foundation

A simplified uniaxial strain controlled creep damage law is deduced with the use of experimental observation from a more complex strain dependent law. This creep damage law correlates the creep damage, which is interpreted as the density variation in the material, directly with the accumulated creep strain. Based on the deduced uniaxial strain controlled creep damage law, a continuum mechanical creep rupture analysis is carried out for a beam resting on a high temperature elastic (Winkler) foundation. The analysis includes the determination of the nondimensional time for initial rupture, the propagation of the rupture front with the associated thinning of the beam, and the influence of creep damage on the deflection of the beam. Creep damage starts accumulating in the beam as soon as the load is applied, and a creep rupture front develops at and propagates from the point at which the creep damage first reaches its critical value. By introducing a series of fundamental assumptions within the framework of technical Euler-Bernoulli type beam theory, a governing set of integro-differential equations is derived in terms of the nondimensional bending moment and the deflection. These governing equations are subjected to a set of interface conditions at the propagating rupture front. A numerical technique is developed to solve the governing equations together with the interface equations, and the computed results are presented and discussed in detail.

Gu, Randy J.↗

Studies of HZE particle interactions and transport for space radiation protection purposes

The main emphasis is on developing general methods for accurately predicting high-energy heavy ion (HZE) particle interactions and transport for use by researchers in mission planning studies, in evaluating astronaut self-shielding factors, and in spacecraft shield design and optimization studies. The two research tasks are: (1) to develop computationally fast and accurate solutions to the Boltzmann (transport) equation; and (2) to develop accurate HZE interaction models, from fundamental physical considerations, for use as inputs into these transport codes. Accurate solutions to the HZE transport problem have been formulated through a combination of analytical and numerical techniques. In addition, theoretical models for the input interaction parameters are under development: stopping powers, nuclear absorption cross sections, and fragmentation parameters.

Townsend, Lawrence W.↗

Nonequilibrium flow computations. 1: An analysis of numerical formulations of conservation laws

Modern numerical techniques employing properties of flux Jacobian matrices are extended to general, nonequilibrium flows. Generalizations of the Beam-Warming scheme, Steger-Warming and van Leer Flux-vector splittings, and Roe's approximate Riemann solver are presented for 3-D, time-varying grids. The analysis is based on a thermodynamic model that includes the most general thermal and chemical nonequilibrium flow of an arbitrary gas. Various special cases are also discussed.

Liu, Yen↗

Mixing control in a plane shear layer

An investigation of mixing processes in a plane shear layer by a direct numerical simulation with a high-order numerical technique, the isoparametric spectral element method, is presented. The computational domain includes a region with the splitter plate that allows to start the flow as a laminar (Blasius) boundary layers with different velocities and concentrations, merging in the shear layer and undergoing transition to highly-unsteady structures with intensive mixing. The augmentation of mixing as a function of a downstream location, the intensity and type of vortical structures in the layer is investigated by an introduction of disturbance at the high-velocity inflow. The character of the dynamic mixing is displayed along with time averaged profiles and statistical characteristics.

Korczak, K. Z.↗

Equilibrium gas flow computations. II - An analysis of numerical formulations of conservation laws

Modern numerical techniques employing properties of flux Jacobian matrices are extended to general, equilibrium gas laws. Generalizations of the Beam-Warming scheme, Steger-Warming and van Leer flux-vector splittings, and Roe's approximate Riemann solver are presented for three-dimensional, time-varying grids. The approximations inherent in previous generalizations are discussed.

Vinokur, Marcel↗

Efficient solutions of two-dimensional incompressible steady viscous flows

A simple, efficient, and robust numerical technique is provided for solving two dimensional incompressible steady viscous flows at moderate to high Reynolds numbers. The proposed approach employs an incremental multigrid method and an extrapolation procedure based on minimum residual concepts to accelerate the convergence rate of a robust block-line-Gauss-Seidel solver for the vorticity-stream function Navier-Stokes equations. Results are presented for the driven cavity flow problem using uniform and nonuniform grids and for the flow past a backward facing step in a channel. For this second problem, mesh refinement and Richardson extrapolation are used to obtain useful benchmark solutions in the full range of Reynolds numbers at which steady laminar flow is established.

Morrison, J. H.↗

Implications of truncating semi-infinite physical domains on the accuracy of the solutions to the N-S equations

Numerical techniques for solving the Navier-Stokes equations of laminar incompressible flow over a forward-facing step on a semiinfinite flat plate are evaluated analytically and by means of sample computations. A primitive-variable method with pseudocompressibility and full mapping of the semiinfinite domain is compared with a vorticity/stream-function method with a truncated domain, applying a flux-corrected explicit finite-difference procedure to solve the discretized transport equations. The results are presented in graphs and compared with published experimental data. The solutions obtained with the truncated-domain methods are shown to be inaccurate, even in cases when they are smooth.

Dekruif, Jeffery S.↗

Error indicators and accuracy improvements of finite element solutions

Practical and reliable estimators of the discretization errors in engineering problems are developed. Error indicators for identifying the regions or elements of the solution domain which are likely to have the largest discretization errors are presented, and a simple computational procedure for improving the accuracy of the finite element solutions for shell problems is given. The similarities between the proposed procedure and a preconditioned conjugate gradient (PCG) technique are identified and exploited to generate pointwise error indicators from the PCG technique. Numerical examples in the linear static analysis of shells are presented.

Noor, Ahmed K.↗

Three-dimensional unsteady transonic viscous-inviscid interaction using the Euler and boundary-layer equations

The objective of this study is the development of a numerical technique which can provide three-dimensional, time-accurate, compressible, turbulent flow solutions in a practical and relatively economical way. The approach taken is that of the method of viscous-inviscid interaction. The Euler equations are assumed to govern the outer inviscid portion of the flow, and the viscous layer close to the solid wall is described by a set of integral boundary-layer equations. The viscous solutions are obtained in a direct fashion with a weighted-average phase error scheme. The method of equivalent sources is used for viscous-inviscid coupling. Steady-state and unsteady computations for an AGARD airfoil and a wing show that satisfactory engineering solutions can be obtained for attached, high Reynolds number flows using this method. Quasi-unsteady interactions are shown to produce similar results to those provided by true-unsteady interactions. Considerable computer resources can be saved for some cases by using quasi-unsteady interactions.

Whitfield, David L.↗