Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerics”

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 955 records · Page 53

Numerical weather prediction in low latitudes

Based on the results of a number of numerical prediction experiments, the differential heating between land and ocean is an important and critical factor for investigation of phenomenon such as the onset of monsoons over the Indian subcontinent. The pre-onset period during the month of May shows a rather persistent flow field in the monsoon region. At low levels the circulation exhibits anticyclonic excursions over the Arabian Sea, flowing essentially parallel to the west coast of India from the north. Over the Indian subcontinent the major feature is a shallow heat low over northern India. As the heat sources commence a rapid northwestward movement toward the southern edge of the Tibetan Plateau, an interesting configuration of the large-scale divergent circulation occurs. A favorable configuration for a rapid exchange of energy from the divergent to the rotational kinetic energy develops. Strong low level monsoonal circulations evolve, attendant with that the onset of monsoon rains occurs. In order to test this observational sequence, a series of short-range numerical prediction experiments were initiated to define the initial heat sources.

Krishnamurti, T. N.↗

Advances in parameterization of cumulus convection in numerical models

Observational studies over the last decade have revealed new sources of difficulty in the already complex problem of parameterization of precipitating cumulus convection in large scale numerical weather prediction models. The structure of mesoscale updrafts and downdrafts and the difficulties they present for cumulus parameterization is reviewed. An approach is proposed which follows Kuo (1974) but includes cumulus and mesoscale downdrafts in the vertical distributions of heating and moistening. The development of a realistic life cycle behavior in the numerical prediction of a mesoscale convective complex will be shown.

Molinari, J.↗

Advances in numerical and applied mathematics

This collection of papers covers some recent developments in numerical analysis and computational fluid dynamics. Some of these studies are of a fundamental nature. They address basic issues such as intermediate boundary conditions for approximate factorization schemes, existence and uniqueness of steady states for time dependent problems, and pitfalls of implicit time stepping. The other studies deal with modern numerical methods such as total variation diminishing schemes, higher order variants of vortex and particle methods, spectral multidomain techniques, and front tracking techniques. There is also a paper on adaptive grids. The fluid dynamics papers treat the classical problems of imcompressible flows in helically coiled pipes, vortex breakdown, and transonic flows.

South, J. C., Jr.↗

Development of a numerical procedure to map a general 3-d body onto a near-circle

Conformal mapping is a classical technique utilized for solving problems in aerodynamics and hydrodynamics. Conformal mapping is utilized in the construction of grids around airfoils, engine inlets and other aircraft configurations. These shapes are transformed onto a near-circle image for which the equations of fluid motion are discretized on the mapped plane and solved numerically by utilizing the appropriate techniques. In comparison to other grid-generation techniques such as algerbraic or differential type, conformal mapping offers an analytical and accurate form even if the grid deformation is large. One of the most appealing features is that the grid can be constrained to remain orthogonal to the body after the transformation. Hence, the grid is suitable for analyzing the supersonic flow past a blunt object. The associated shock as a coordinate surface adjusts its position in the course of computation until convergence is reached. The present work applied conformal mapping to 3-D bodies with no axis of symmetry such as the Aerobraking Flight Experiment (AFE) vehicle, transforming the AFE shape onto a near-circle image. A numerical procedure and code are used to generate grids around the AFE body.

Hommel, M. J.↗

Numerical simulation of boundary layers. Part 2: Ribbon-induced transition in Blasius flow

The early three-dimensional stages of transition in Blasius boundary layers are studied by numerical solution of the Navier-Stokes equations. A finite-amplitude two-dimensional wave and random low-amplitude three-dimensional disturbances are introduced. Rapid amplification of the three-dimensional components is observed and leads to transition. For intermediate amplitudes of the two-dimensional wave the breakdown is of subharmonic type, and the dominant spanwise wave number increases with the amplitude. For high amplitudes the energy of the fundamental mode is comparable to the energy of the subharmonic mode, but never dominates it; the breakdown is of mixed type. Visualizations, energy histories, and spectra are presented. The sensitivity of the results to various physical and numerical parameters is studied. Agreement with experimental and theoretical results is discussed.

Spalart, P.↗

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.↗

Comparison of uniform perturbation and numerical solutions for some potential flows past slender bodies

Approximate solutions for potential flow past an axisymmetric slender body and past a thin airfoil are calculated using a uniform perturbation method and then compared with either the exact analytical solution or the solution obtained using a purely numerical method. The perturbation method is based upon a representation of the disturbance flow as the superposition of singularities distributed entirely within the body, while the numerical (panel) method is based upon a distribution of singularities on the surface of the body. It is found that the perturbation method provides very good results for small values of the slenderness ratio and for small angles of attack. Moreover, for comparable accuracy, the perturbation method is simpler to implement, requires less computer memory, and generally uses less computation time than the panel method. In particular, the uniform perturbation method yields good resolution near the regions of the leading and trailing edges where other methods fail or require special attention.

Wong, T.-C.↗

The influence of surface roughness of deserts on the July circulation - A numerical study

The effect of the low surface roughness characteristics of deserts on atmospheric circulation in July is examined using numerical simulations with the GCM of the Goddard Laboratory for Atmospheric Science (GLAS). Identical sets of simulations were carried out with the model starting from the initial state of the atmosphere on June 15, for the years 1979 and 1980. The first simulation included a surface roughness factor of 45 cm, and the second set had a surface roughness factor of 0.02 cm for desert regions, and 45 cm for all other land. A comparative analysis of the numerical data was carried out in order to study the variations for the desert regions. It is shown that rainfall in the Sahara desert was reduced significantly in the data set with the nonuniform surface roughness factor in comparison with the other data set. The inter-tropical convergence zone (ITCZ) moved southward to about 15 degrees, which was close to its observed location at about 10 degrees N. In other deserts, the North American Great Plains, Rajputana in India, and the Central Asian desert, no similar changes were observed. Detailed contour maps of the weather conditions in the different desert regions are provided.

Sud, Y. C.↗

Numerical simulation of transonic flow over porous airfoils

A numerical study was made to examine the effect of a porous surface on the aerodynamic performance of a transonic airfoil. The pressure jump across the normal shock wave on the upper surface of the airfoil was reduced by making the surface below the shock porous. The weakened shock is preceded by an oblique shock at the upstream end of the porous surface where air is blown out of the cavity. The lambda shock structure shown in the numerical result qualitatively agrees with that observed in the wind tunnel. According to the present analysis, the porous airfoil has a smaller drag and a higher lift than the solid airfoil.

Chow, C.-Y.↗

Numerical simulation of self-excited oscillations in a ramjet inlet-diffuser flow

This paper describes numerical simulations of self-excited oscillations in a two-dimensional transonic inlet-diffuser flow by solving the Navier-Stokes equations with a two-equation turbulence model. The calculated amplitudes of oscillations for the terminal shock and the velocity fields compare well with experimental measurements; however, the predicted frequency of oscillations is about 50 percent higher. The formation of a pair of downstream-traveling, counter-rotating vortices at each cycle of velocity fluctuations, as reported experimentally, is vividly revealed by the numerical results.

Hsieh, T.↗

Numerical schemes for a model for nonlinear dispersive waves

A description is given of a number of numerical schemes to solve an evolution equation (Korteweg-deVries) that arises when modelling the propagation of water waves in a channel. The discussion also includes the results of numerical experiments made with each of the schemes. It is suggested, on the basis of these experiments, that one of the schemes may have (discrete) solitary-wave solutions.

Bona, J. L.↗

Numerical simulation of unsteady flow in an axisymmetric shear layer

Numerical simulation of unsteady flow in an axisymmetric subsonic shear layer is accomplished by solving the time-dependent compressible Navier-Stokes equations. The objective of the effort is to investigate by numerical means the influence of various flow parameters on the shear layer behavior. The parameters investigated include the velocity ratio of two streams, total temperature, and nozzle lip thickness. The computations are performed on a CRAY-IS computer using McCormack's explicit finite difference scheme. The computed results generally show qualitative agreement with experimental data.

Scott, J. N.↗

A numerical study of three-dimensional incompressible flow around multiple posts

The incompressible laminar flow around single and double rows of an infinite number of posts between two flat plates has been calculated numerically. A C-grid with periodic boundary conditions has been used. The angle of attack, measured from the line normal to the row of posts, was varied between zero and less than 90 deg. The pressure loading was computed for all of the posts in each of these cases. Most of these calculations have been carried out on the Numerical Aerodynamic Simulation Cray 2 at NASA Ames Research Center.

Rogers, S. E.↗

Numerical study of vortex breakdown

The incompressible axisymmetric steady Navier-Stokes equations and the Euler equations are solved numerically to model the breakdown of a vortex. Although the solutions obtained for the Euler equations show a 'vortex breakdown-like' structure, their behavior is very different from that of the Navier-Stokes equations. The details of the numerical algorithmms used are presented, and the results obtained are compared to those in the literature.

Hafez, M.↗

Cirrus clouds. I - A cirrus cloud model. II - Numerical experiments on the formation and maintenance of cirrus

A simplified cirrus cloud model is presented which may be used to investigate the role of various physical processes in the life cycle of a cirrus cloud. The model is a two-dimensional, time-dependent, Eulerian numerical model where the focus is on cloud-scale processes. Parametrizations are developed to account for phase changes of water, radiative processes, and the effects of microphysical structure on the vertical flux of ice water. The results of a simulation of a thin cirrostratus cloud are given. The results of numerical experiments performed with the model are described in order to demonstrate the important role of cloud-scale processes in determining the cloud properties maintained in response to larger scale forcing. The effects of microphysical composition and radiative processes are considered, as well as their interaction with thermodynamic and dynamic processes within the cloud. It is shown that cirrus clouds operate in an entirely different manner than liquid phase stratiform clouds.

Starr, D. OC.↗

Numerical-perturbation technique for stability of flat-plate boundary layers with suction

A numerical-perturbation scheme is proposed for determining the stability of flows over plates with suction through a finite number of porous suction strips. The basic flow is calculated as the sum of the Blasius flow and closed-form linearized triple-deck solutions of the flow due to the strips. A perturbation technique is used to determine the increment a(ij) in the complex wavenumber at a given location x(j) due to the presence of a strip centered at x(i). The end result is a set of influence coefficients that can be used to determine the growth rates and amplification factors for any suction levels without repeating the calculations. The numerical-perturbation results are verified by comparison with interacting boundary layers for the case of six strips and the experimental data of Reynolds and Saric for single- and multiple-strip configurations. The influence coefficient form of the solution suggests a scheme for optimizing the strip configuration. The results show that one should concentrate the suction near branch I of the neutral stability curve, a conclusion verified by the experiments.

Reed, H. L.↗

Numerical simulation of tip vortices of wings in subsonic and transonic flows

A multi block zonal algorithm which solves the thin-layer Navier-Stokes and the Euler equations is used to numerically simulate the formation and roll-up of the tip vortex in both subsonic and transonic flows. Four test cases which used small and large aspect ratio wings have been considered to examine the influence of the tip-cap shape, the tip planform and the free-stream Mach number. It appears that both the tip-planform and the tip-cap shape have some influence on the formation of the tip vortex, but its subsequent roll-up seems to be more influenced by the tip-planform shape. In general, a good definition of the formation and the roll-up of the tip vortex has been observed for all the cases considered here. Comparions of the numerical results with the limited, available experimental data show good agreement with both the surface pressures and the tip-vortex strength.

Srinivasan, G. R.↗

Numerical difficulties associated with using equality constraints to achieve multi-level decomposition in structural optimization

There has been recent interest in multidisciplinary multilevel optimization applied to large engineering systems. The usual approach is to divide the system into a hierarchy of subsystems with ever increasing detail in the analysis focus. Equality constraints are usually placed on various design quantities at every successive level to ensure consistency between levels. In many previous applications these equality constraints were eliminated by reducing the number of design variables. In complex systems this may not be possible and these equality constraints may have to be retained in the optimization process. In this paper the impact of such a retention is examined for a simple portal frame problem. It is shown that the equality constraints introduce numerical difficulties, and that the numerical solution becomes very sensitive to optimization parameters for a wide range of optimization algorithms.

Thareja, R.↗