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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,585 records · Page 88

All sources NASA NTRS DOE OSTI

NASA NTRS1983

An analysis of the two-dimensional acoustic field in a nonuniform duct carrying compressible flow

An analytical/numerical linear acoustic solution in a nonuniform two-dimensional duct carrying a compressible mean flow is developed. A quasi-one dimensional mean flow model is employed together with a consistent expression for the cross-flow velocity. The acoustic solution is obtained using the Wave Envelope Method. Numerical results are compared with those of an existing Wave Envelope solution which includes a more general mean flow model. In addition, the singular behavior of the solutions derived in the current work for nearly sonic flow conditions is examined explicitly.

Uenishi, K.↗

NASA NTRS1983

Modeling of steady, rotational, transonic winds from rotating stars and galaxies

The theory of steady transonic winds from condensed bodies is extended to general, two-dimensional, axisymmetric systems. A stream function is used to reduce the gasdynamics equations to a single, second-order differential equation plus an algebraic equation for the density. The approach extends Parker's (1958) quasi-one-dimensional theory, which uses Bernoulli's theorem, to a complete two-dimensional calculation which includes vorticity and rotation. The conceptual basis of the stream function approach is described in detail, and a numerical method for solving the resulting equations is presented. The applications illustrate solutions for stellar (spherical source surface) and galactic (oblate spheroidal source surface) models. Among other things, it is found that for rapidly rotating stellar models the Coriolis force dominates the centrifugal terms and the streamlines bend toward the rotation axis as a consequence of the axisymmetry. For both the rotating and nonrotating galaxy models, the streamlines bend toward the equator.

Kopriva, D. A.↗

NASA NTRS1983

Propagation of solar disturbances - Theories and models

Recent theoretical developments and construction of several models for the propagation of solar disturbances from the sun and their continuation throughout heliospheric space are discussed. Emphasis centers on physical mechanisms as well as mathematical techniques (i.e., analytical and numerical methods). This outline will lead to a discussion of the state-of-the-art of theoretically based modeling efforts in this area. It is shown that the fundamental theory for the study of propagation of disturbances in heliospheric space is centered around the self-consistent analysis of wave and mass motion within the context of magnetohydrodynamics in which the small scale structures will be modified by kinetic effects. Finally, brief mention is made of some interesting problems for which attention is needed for advancement of the understanding of the physics of large scale propagation of solar disturbances in heliospheric space.

Wu, S. T.↗

NASA NTRS1982

The use of solution adaptive grids in solving partial differential equations

The grid point distribution used in solving a partial differential equation using a numerical method has a substantial influence on the quality of the solution. An adaptive grid which adjusts as the solution changes provides the best results when the number of grid points available for use during the calculation is fixed. Basic concepts used in generating and applying adaptive grids are reviewed in this paper, and examples illustrating applications of these concepts are presented.

Anderson, D. A.↗

NASA NTRS1983

An interacting boundary layer model for cascades

A laminar, incompressible interacting boundary layer model is developed for two-dimensional cascades. In the limit of large cascade spacing these equations reduce to the interacting boundary layer equations for a single body immersed in an infinite stream. A fully implicit numerical method is used to solve the governing equations, and is found to be at least as efficient as the same technique applied to the single body problem. Solutions are then presented for a cascade of finite flat plates and a cascade of finite sine-waves, with cusped leading and trailing edges.

Davis, R. T.↗

NASA NTRS1983

A least squares finite element scheme for transonic flow around harmonically oscillating airfoils

The present investigation shows that a finite element scheme with a weighted least squares variational principle is applicable to the problem of transonic flow around a harmonically oscillating airfoil. For the flat plate case, numerical results compare favorably with the exact solution. The obtained numerical results for the transonic problem, for which an exact solution is not known, have the characteristics of known experimental results. It is demonstrated that the performance of the employed numerical method is independent of equation type (elliptic or hyperbolic) and frequency. The weighted least squares principle allows the appropriate modeling of singularities, which such a modeling of singularities is not possible with normal least squares.

Cox, C. L.↗

NASA NTRS1983

Computation of turbulent flows over backward-facing step

A numerical method for computing incompressible turbulent flows is presented. The method is tested by calculating laminar recirculating flows and is applied in conjunction with a modified Kappa-epsilon model to compute the flow over a backward-facing step. In the laminar regime, the computational results are in good agreement with the experimental data. The turbulent flow study shows that the reattachment length is underpredicted by the standard Kappa-epsilon model. The addition of a term to the standard model that accounts for the effects of rotation on turbulent flow improves the results in the recirculation region and increases the computed reattachment length.

Mansour, N. N.↗

NASA NTRS1983

The mechanics of delamination in fiber-reinforced composite materials. Part 2: Delamination behavior and fracture mechanics parameters

Based on theories of laminate anisotropic elasticity and interlaminar fracture, the complete solution structure associated with a composite delamination is determined. Fracture mechanics parameters characterizing the interlaminar crack behavior are defined from asymptotic stress solutions for delaminations with different crack-tip deformation configurations. A numerical method employing singular finite elements is developed to study delaminations in fiber composites with any arbitrary combinations of lamination, material, geometric, and crack variables. The special finite elements include the exact delamination stress singularity in its formulation. The method is shown to be computationally accurate and efficient, and operationally simple. To illustrate the basic nature of composite delamination, solutions are shown for edge-delaminated (0/-0/-0/0) and (+ or - 0/+ or - 0/90/90 deg) graphite-epoxy systems under uniform axial extenstion. Three-dimensional crack-tip stress intensity factors, associated energy release rates, and delamination crack-closure are determined for each individual case. The basic mechanics and mechanisms of composite delamination are studied, and fundamental characteristics unique to recently proposed tests for interlaminar fracture toughness of fiber composite laminates are examined.

Wang, S. S.↗

NASA NTRS1983

Acoustic shocks in a variable area duct containing near sonic flows

Acoustic shock waves in a variable area duct which contains near sonic flows are considered. The problem is modeled after an aeroengine inlet. Area variation of a duct and high Mach number mean the flow reduces acoustical energy yielding substantial noise reduction. One possible reason for this is acoustic shock. The use of an explicit accurate numerical method which captures shocks is described. Comparison of the results are made with an existing asymptotic theory for Mach numbers close to unity. When shock occurs reduction of sound pressure levels are shown by example.

Hariharan, S. I.↗

NASA NTRS1984

Delta wings with shock-free cross flow

In order to have a high level of maneuverability, supersonic delta wings should have a cross flow that is free of embedded shock waves. The conical cross flow sonic surface differs from that of plane transonic flow in many aspects. Well-known properties such as the monotone law are not true for conical cross flow sonic surfaces. By using a local analysis of the cross flow sonic line, relevant conditions for smooth cross flow are obtained. A technique to artificially construct a smooth sonic surface and an efficient numerical method to calculate the flow field are used to obtain cones with smooth cross flow.

Sritharan, S. S.↗

NASA NTRS1983

Cumulative reports and publications through 31 December 1983

All reports for the calendar years 1975 through December 1983 are listed by author. Since ICASE reports are intended to be preprints of articles for journals and conference proceedings, the published reference is included when available. Thirteen older journal and conference proceedings references are included as well as five additional reports by ICASE personnel. Major categories of research covered include: (1) numerical methods, with particular emphasis on the development and analysis of basic algorithms; (2) computational problems in engineering and the physical sciences, particularly fluid dynamics, acoustics, structural analysis, and chemistry; and (3) computer systems and software, especially vector and parallel computers, microcomputers, and data management.

Source record↗

NASA NTRS1984

Analysis of inviscid and viscous flows in cascades with an explicit multiple-grid algorithm

A rapid technique is used for calculating inviscid and viscous flows in turbomachinery cascades. The Euler and thin-layer Navier-Stokes equations are solved using the original explicit MacCormack algorithm. The Baldwin-Lomax eddy viscosity model is used for turbulent flows. Convergence to a steady state is accelerated by use of a variable time-step and a multiple-grid scheme. Computer time is reduced through vectorization. Details of the numerical method are presented along with computed results for two low-speed wind tunnel turning vanes, a space shuttle fuel pump turbine rotor, and a supersonic inflow compressor rotor. The method can predict subtle viscous flow phenomena in cascades and is fast enough to be used as a design tool.

Chima, R. V.↗

NASA NTRS1984

Finite-difference fluid dynamics computer mathematical models for the design and interpretation of experiments for space flight

Numerical methods are used to design a spherical baroclinic flow model experiment of the large scale atmosphere flow for Spacelab. The dielectric simulation of radial gravity is only dominant in a low gravity environment. Computer codes are developed to study the processes at work in crystal growing systems which are also candidates for space flight. Crystalline materials rarely achieve their potential properties because of imperfections and component concentration variations. Thermosolutal convection in the liquid melt can be the cause of these imperfections. Such convection is suppressed in a low gravity environment. Two and three dimensional finite difference codes are being used for this work. Nonuniform meshes and implicit iterative methods are used. The iterative method for steady solutions is based on time stepping but has the options of different time steps for velocity and temperature and of a time step varying smoothly with position according to specified powers of the mesh spacings. This allows for more rapid convergence. The code being developed for the crystal growth studies allows for growth of the crystal as the solid-liquid interface. The moving interface is followed using finite differences; shape variations are permitted. For convenience in applying finite differences in the solid and liquid, a time dependent coordinate transformation is used to make this interface a coordinate surface.

Roberts, G. O.↗

NASA NTRS1984

Space structure vibration modes: How many exist? Which ones are important?

This report attempts to shed some light on the two issues raised in the title, namely, how many vibration modes does a real structure have, and which of these modes are important? The surprise-free answers to these two questions are, respectively, an infinite number and the first several modes. The author argues that the absurd subspace (all but the first billion modes) is not a strength of continuum modeling, but, in fact, a weakness. Partial differential equations are not real structures, only mathematical models. This note also explains (1) that the PDE model and the finite element model are, in fact, the same model, the latter being a numerical method for dealing with the former, (2) that modes may be selected on dynamical grounds other than frequency alone, and (3) that long slender rods are useful as primitive cases but dangerous to extrapolate from.

Hughes, P. C.↗

NASA NTRS1984

Numerical simulation of boundary-layer transition

The transition to turbulence in boundary layers was investigated by direct numerical solution of the nonlinear, three-dimensional, incompressible Navier-Stokes equations in the half-infinite domain over a flat plate. Periodicity was imposed in the streamwise and spanwise directions. A body force was applied to approximate the effect of a nonparallel mean flow. The numerical method was spectra, based on Fourier series and Jacobi polynomials, and used divergence-free basis functions. Extremely rapid convergence was obtained when solving the linear Orr-Sommerfeld equation. The early nonlinear and three-dimensional stages of transition, in a boundary layer disturbed by a vibrating ribbon, were successfully simulated. Excellent qualitative agreement was observed with either experiments or weakly nonlinear theories. In particular, the breakdown pattern was staggered or nonstaggered depending on the disturbance amplitude.

Spalart, P. R.↗

NASA NTRS1984

An experimental and theoretical analysis of the aerodynamic characteristics of a biplane-winglet configuration

Improving the aerodynamic characteristics of an airplane with respect to maximizing lift and minimizing induced and parasite drag are of primary importance in designing lighter, faster, and more efficient aircraft. Previous research has shown that a properly designed biplane wing system can perform superiorly to an equivalent monoplane system with regard to maximizing the lift-to-drag ratio and efficiency factor. Biplanes offer several potential advantages over equivalent monoplanes, such as a 60-percent reduction in weight, greater structural integrity, and increased roll response. The purpose of this research is to examine, both theoretically and experimentally, the possibility of further improving the aerodynamic characteristics of the biplanes configuration by adding winglets. Theoretical predictions were carried out utilizing vortex-lattice theory, which is a numerical method based on potential flow theory. Experimental data were obtained by testing a model in the Pennsylvania State University's subsonic wind tunnel at a Reynolds number of 510,000. The results showed that the addition of winglets improved the performance of the biplane with respect to increasing the lift-curve slope, increasing the maximum lift coefficient, increasing the efficiency factor, and decreasing the induced drag. A listing of the program is included in the Appendix.

Gall, P. D.↗

NASA NTRS1984

Predictive thermal-vacuum modeling by variable coefficients

In computer-assisted thermal-vacuum testing, an accurate mathematical thermal model is necessary. This model can indicate what effect a given change in control temperature will produce on the test article. A three-node thermal model is used to illustrate some numerical methods for increasing the accuracy by generating more accurate coefficients for the model.

Chuvala, J. T.↗

NASA NTRS1984

Nonlinear computation of wing-body-vertical tail-wake flows at low supersonic speeds

A numerical method based on the conservation form of the full potential equation has been applied to the problem of three-dimensional supersonic flows with embedded subsonic regions. The governing equation is cast in a nonorthogonal coordinate system, and the theory of characteristics is used to accurately monitor the type-dependent flow field. A conservative switching scheme is employed to transition from the supersonic marching procedure to a subsonic relaxation algorithm and vice versa. The newly developed computer program can handle arbitrary geometries with fuselage, wing, vertical tail and wake components at combined angle of attack and sideslip. Results are presented for a low supersonic Mach numbers flow over the Shuttle orbiter (including the OMS pods and vertical tail), and for flows over a realistic fighter type configuration. Comparisons with experimental data are shown to be in good agreement for various cases.

Szema, K.-Y.↗