Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical method”

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 1,513 records · Page 84

Singularity computations

Direct descriptions of the structure of a singularity would describe the radial and angular distributions of the field quantities as explicitly as practicable along with some measure of the intensity of the singularity. This paper discusses such an approach based on recent development of numerical methods for elastoplastic flow. Attention is restricted to problems where one variable or set of variables is finite at the origin of the singularity but a second set is not.

Swedlow, J. L.↗

A computational fluid dynamic technique valid at the centerline for non-axisymmetric problems in cylindrical coordinates

A technique is described for the numerical solution of non-axisymmetric flow problems posed in cylindrical coordinates when the z-axis is included in the flowfield. The highlight of the technique is the manner in which the singularities at the centerline are handled. Specifically, the governing flowfield equations at r = 0 are put in a special form by applying L'Hospital's Rule. The required radial derivatives are evaluated using a one-sided, second-order accurate, first-difference. This leads to a smooth, convergent calculation of the flowfield at the centerline. This appears to be the first generally applicable numerical method for avoiding coordinate system singularities in the context of a finite-difference scheme, and could have application to many nonaxisymmetric flows. The technique is illustrated by specific results for the time-dependent flowfield inside an internal combustion engine.

Griffin, M. D.↗

Calculated scan characteristics of a large spherical reflector antenna

A previously published numerical method to calculate the radiation properties of parabolic reflectors has been modified to also include very large spherical reflectors. The method has been verified by comparing the calculated and the measured results for a 120-wavelength spherical reflector.

Agrawal, P. K.↗

A theoretical study of three-dimensional barotropic instability with applications to the upper stratosphere

A numerical method for tropospheric baroclinic instability is used to establish the structures, phase speeds and growth rates of the normal modes of a number of idealized barotropically unstable mean zonal wind fields in a quasi-geostrophic Boussinesq framework. The objective was to gain insight into the effects of vertical mean flow variation on barotropically unstable stratosphere. Both horizontally symmetric and asymmetric flows are considered, but in all cases the flows are vertically symmetric. Major findings show that (1) the wave growth rates associated with the strongest expected horizontal shears decrease with increasing vertical shear and are at least 28% smaller for the weakest expected vertical shear than with no shear; (2) differences in the form of either the vertical or the horizontal shear flow can affect this growth rate reduction by about 25%; and (3) the vertical scale of the barotropic unstable waves is of the order of the geometric mean of the vertical mean flow scale and the vertical penetrat

Pfister, L.↗

Large eddy simulation - A predictive approach to turbulent flow computation

The filtering approach (Leonard, 1974) to large eddy simulation is described. Large eddy simulation is then surveyed with reference to subgrid scale models, numerical methods, channel flow, and free shear flows. It is concluded that the premise of LES (i.e., that the small eddies are much easier to model than the large ones) seems to be borne out by computational experience to date. The greatest obstacle in LES at present is the lack of computer power which effectively denies the ability to simulate all the significant eddies in flows of interest.

Ferziger, J. H.↗

Effects of plasma sheath on solar power satellite array

The structure of the plasma sheath and equilibrium voltage distribution of a high-power solar array governs various kinds of plasma-interaction phenomena and array losses. Sheath effects of a linearly-connected array are investigated for GEO. Although the array may be large, the thin-sheath-limit analysis may be invalid, necessitating numerical methods. Three-dimensional computer calculations show that potential barriers and over-lapping sheaths can occur, i.e., structures not predictable under the thin-sheath-limit analysis, but nevertheless controlling the distribution of plasma currents impacting on the array.

Parker, L. W.↗

A viscous flow analysis of the tip vortex generation process

A viscous primary-secondary flow analysis designed for numerical solution by forward marching integration is applied to the tip vortex generation problem. Equations governing streamwise momentum, streamwise vorticity, and a gas law are solved in conjunction with equations governing scalar and vector surface potentials for the secondary velocities. The numerical method used combines a consistently split linearized block implicit (LBI) scheme for parabolic equations and a scalar iterative ADI scheme for elliptic equations. Computed results are presented for a wing of rectangular planform immersed in a high Reynolds number stream at 6 degrees incidence. The physical mechanism for generation of the tip vortex, as revealed by the analysis and computations, is discussed. Detailed flow results for the case of 6 degree incidence are included.

Shamroth, S. J.↗

Optimum design of laminated plates under axial compression

The paper presents a method for designing laminated plates with orthotropic layers under uniaxial and biaxial compression. The objective is to find the fiber direction of all the layers that give the highest buckling stress; therefore, it is necessary to solve an unconstrained maximization problem. The optimization technique used is Powell's conjugate direction method. Numerical calculations were made for plates with three, four, and six layers.

Hirano, Y.↗

An extended soft-cube model for the thermal accommodation of gas atoms on solid surfaces

A numerical soft cube model was developed for calculating thermal accommodation coefficients alpha and trapping fractions f sub t for the interaction of gases incident upon solid surfaces. A semiempirical correction factor c which allows the calculation of alpha and f sub t when the collision times are long compared to the surface oscillator period were introduced. The processes of trapping, evaporation, and detailed balancing were discussed. The numerical method was designed to treat economically and with moderate (+ or - 20 percent) accuracy the dependence of alpha and f sub t on finite and different surface and gas temperatures for a large number of gas/surface combinations. Comparison was made with experiments of rare gases on tungsten and on alkalis, as well as one astrophysical case of H2 on graphite. The dependence of alpha on the soft cube dimensionless parameters is presented graphically.

Burke, J. R.↗

Symbolic-numeric interface: A review

A survey of the use of a combination of symbolic and numerical calculations is presented. Symbolic calculations primarily refer to the computer processing of procedures from classical algebra, analysis, and calculus. Numerical calculations refer to both numerical mathematics research and scientific computation. This survey is intended to point out a large number of problem areas where a cooperation of symbolic and numerical methods is likely to bear many fruits. These areas include such classical operations as differentiation and integration, such diverse activities as function approximations and qualitative analysis, and such contemporary topics as finite element calculations and computation complexity. It is contended that other less obvious topics such as the fast Fourier transform, linear algebra, nonlinear analysis and error analysis would also benefit from a synergistic approach.

Ng, E. W.↗

Three dimensional thermal pollution models. Volume 1: Review of mathematical formulations

A mathematical model package for thermal pollution analyses and prediction is presented. These models, intended as user's manuals, are three dimensional and time dependent using the primitive equation approach. Although they have sufficient generality for application at sites with diverse topographical features; they also present specific instructions regarding data preparation for program execution and sample problems. The mathematical formulation of these models is presented including assumptions, approximations, governing equations, boundary and initial conditions, numerical method of solution, and same results.

Lee, S. S.↗

Mountain winds (revisited)

The prediction of extremely high wind speeds, at ground level on the downstream side of a mountain range, is possible by solving the initial value problem for a two-layered nonlinear shallow water model of the atmosphere. Three different numerical methods are described to find the solutions which may involve shocks: (1) the vonNeumann-Richtmyer artificial viscosity method, (2) a filtering scheme, and (3) a hybrid method.

Isaacson, E.↗

An experimental and analytical investigation of a hypersonic cruise vehicle at Mach 6

Flow field properties over a hypersonic vehicle are measured and compared with results from the numerical analysis of inviscid supersonic flow theory. Good agreement between theory and experiment at zero angle of attack, and significant discrepancy at angles of attack greater than or equal to 5 deg are found. The experiments indicate that turbulent boundary layer without flow separation exists over the model surface at zero angle of attack. Vortex lift-off, flow separation, and shock boundary layer interaction occur over the leeward surface, and near the wing-fuselage junction with angles of attack greater than or equal to 5 deg. Improvements of the existing numerical method to compute hypersonic cruise vehicle flow field at large angle of attack are presented.

Wang, C. R.↗

A definition of the degree of controllability - A criterion for actuator placement

The unsolved problem of how to control the attitude and shape of future very large flexible satellite structures represents a challenging problem for modern control theory. One aspect of this problem is the question of how to choose the number and locations throughout the spacecraft of the control system actuators. Starting from basic physical considerations, this paper develops a concept of the degree of controllability of a control system, and then develops numerical methods to generate approximate values of the degree of controllability for any spacecraft. These results offer the control system designer a tool which allows him to rank the effectiveness of alternative actuator distributions, and hence to choose the actuator locations on a rational basis. The degree of controllability is shown to take a particularly simple form when the satellite dynamics equations are in modal form. Examples are provided to illustrate the use of the concept on a simple flexible spacecraft.

C N Viswanathan↗

Stress analysis of bonded plates and joints

The paper presents a general numerical method of solution to the transverse normal stress and longitudinal shear stress distributions in the adhesive layers of bonded stiffener plates, double lap joints, and strap joints, all of dissimilar orthotropic adherends. The bending deformations and the transverse shear deformations are included in the analytical model. The present influence of bending deformations of adherends on both the normal and shear stresses in adhesive layers is discussed in detail for practical cases.

Yuceoglu, U.↗

A two-dimensional finite-difference solution for the transient thermal behavior of a tubular solar collector

A numerical procedure was established using the finite-difference technique in the determination of the time-varying temperature distribution of a tubular solar collector under changing solar radiancy and ambient temperature. Three types of spatial discretization processes were considered and compared for their accuracy of computations and for selection of the shortest computer time and cost. The stability criteria of this technique were analyzed in detail to give the critical time increment to ensure stable computations. The results of the numerical analysis were in good agreement with the analytical solution previously reported. The numerical method proved to be a powerful tool in the investigation of the collector sensitivity to two different flow patterns and several flow control mechanisms.

Lansing, F. L.↗

A comprehensive comparison between experiment and prediction for a transonic turbulent separated flow

Attempts to predict surface pressure distributions on lifting surfaces have been relatively unsuccessful in the transonic regime when the shock wave is of sufficient strength to produce an extensive region of turbulent separated flow. For these conditions, the viscous flow behavior must be accurately described even to obtain reasonable predictions of surface pressure. The present paper addresses this problem. Detailed comparisons between prediction and experiment are made for a transonic, turbulent boundary-layer separation (freestream Mach number = 0.875) for which the turbulent flow properties (including the turbulent Reynolds stress) had been measured by the laser velocimeter technique from upstream of the separated region through reattachment. The flow was generated on an axisymmetric 'bump' model designed to simulate the flow on an airfoil at transonic conditions. The numerical methods used in the comparisons include the solution of the time-dependent, mass-averaged Navier-Stokes equations, and the solution of the compressible boundary-layer equations by the inverse method. Solutions were obtained for the well established Cebeci-Smith algebraic turbulence model and the more recently developed Wilcox-Rubesin two-equation turbulence model.

Johnson, D. A.↗