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.

SourceResearch connectionFrom yearThrough year

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

At least 1,009 records · Page 56

All sources NASA NTRS DOE OSTI

NASA NTRS1987

A numerical study of transient, thermally-conductive solar wind

A numerical analysis of transient solar wind starting at the solar surface and arriving at 1 AU is performed by an implicit numerical method. The model hydrodynamic equations include thermal conduction terms for both steady and unsteady simulations. Simulation results show significant influence of thermal conduction on both steady and time-dependent solar wind. Higher thermal conduction results in higher solar wind speed, higher temperature, but lower plasma density at 1 AU. Higher base temperature at the solar surface gives lower plasma speed, lower temperature, but higher density at 1 AU. Higher base density, on the other hand, gives lower velocity, lower temperature, but higher density at 1 AU.

Han, S. M.↗

NASA NTRS1987

Use of simulated satellite radiances from a mesoscale numerical model to understand kinematic and dynamic processes

Techniques by which mesoscale numerical weather prediction model output and radiative transfer codes are combined to simulate the radiance fields that a given passive temperature/moisture satellite sensor would see if viewing the evolving model atmosphere are introduced. The goals are to diagnose the dynamical atmospheric processes responsible for recurring patterns in observed satellite radiance fields, and to develop techniques to anticipate the ability of satellite sensor systems to depict atmospheric structures and provide information useful for numerical weather prediction (NWP). The concept of linking radiative transfer and dynamical NWP codes is demonstrated with time sequences of simulated radiance imagery in the 24 TIROS vertical sounder channels derived from model integrations for March 6, 1982.

Kalb, Michael↗

NASA NTRS1987

A Numerical Study of the Effects of Curvature and Convergence on Dilution Jet Mixing

An analytical program was conducted to assemble and assess a three-dimensional turbulent viscous flow computer code capable of analyzing the flow field in the transition liners of small gas turbine engines. This code is of the TEACH type with hybrid numerics, and uses the power law and SIMPLER algorithms, an orthogonal curvilinear coordinate system, and an algebraic Reynolds stress turbulence model. The assessments performed in this study, consistent with results in the literature, showed that in its present form this code is capable of predicting trends and qualitative results. The assembled code was used to perform a numerical experiment to investigate the effects of curvature and convergence in the transition liner on the mixing of single and opposed rows of cool dilution jets injected into a hot mainstream flow.

Holdeman, J. D.↗

NASA NTRS1986

Numerical simulation of boundary layers. Part 3: Turbulence and relaminarization in sink flows

Direct numerical simulations of sink-flow boundary layers, with acceleration parameters K between 1.5 and 3.0 x 10 to the -6 power, are presented. The three-dimensional, time-dependent Navier-Stokes equations are solved numerically using a spectral method, with about 10 to the -6 degrees of freedom. The flow is assumed to be statistically steady, and self-similar. A multiple-scale approximation and periodic conditions are applied to the fluctuations. The turbulence is studied using instantaneous and statistical results. Good agreement with the experiments of Jones and Launder is observed. Two effects of the favorable pressure gradient are to extend the logarithmic layer, and to alter the energy balance of the turbulence near the edge of the boundary layer. At very low Reynolds number the logarithmic layer is shortened and slightly displaced, but wall-layer streaks are present even at the lowest values of the Reynolds number for which turbulence can be sustained. Large quiescent patches appear in the flow. Relaminarization occurs at K = 3.0 x 10 to the -6 power, corresponding to a Reynolds number, based on momentum thickness, of about 330.

Spalart, Philippe R.↗

NASA NTRS1987

Numerical modeling of troposphere-induced gravity wave propagation

Sources of internal gravity waves (IGW) in the upper atmosphere are assumed to be meteorological processes in the troposphere. These sources are vertically and horizontally inhomogeneous and time dependent. In order to describe the IGW propagation from such sources, a numerical solution of a system of hydrodynamical equations is required. In addition, it is necessary to take into account the influence of the altitude latitude inhomogeneity of the temperature and wind fields on the IGW propagation as well as the processes of dissipation. An algorithm is proposed for numerical modelling of the IGW propagation over a limited area from tropospheric local sources to the upper atmosphere. The algorithm takes into account all the above features. A spectral grid method is used with the expansion of wave fields into the Fourier series over longitude. The upper limit conditions were obtained from the requirement of a limited energy dissipation rate in an atmospheric column. The no slip (zero velocity) condition was used at the Earth's surface.

Gavrilov, N. M.↗

NASA NTRS1988

Numerical Simulation of Subsonic and Transonic Propeller Flow

The numerical simulation of 3-D transonic flow about a system of propeller blades is investigated. In particular, it is shown that the use of helical coordinates significantly simplifies the form of the governing equation when the propeller system is assumed to be surrounded by an irrotational flow field of an inviscid fluid. The unsteady small disturbance equation, valid for lightly loaded blades and expressed in helical coordinates, is derived from the general blade-fixed potential equation, given for an arbitrary coordinate system. The use of a coordinate system which inherently adapts to the mean flow results in a disturbance equation requiring relatively few terms to accurately model the physics of the flow. Furthermore, the helical coordinate system presented here is novel in that it is periodic in the circumferential direction while, simultaneously, maintaining orthogonal properties at the mean blade locations. The periodic characteristic allows a complete cascade of blades to be treated, and the orthogonality property affords straightforward treatment of blade boundary conditions. An ADI numerical scheme is used to compute the solution of the steady flow as an asymptotic limit of an unsteady flow. As an example of the method, solutions are presented for subsonic and transonic flow about a 5 percent thick bicircular arc blade of an 8-bladed cascade. Both high and low advance ratio cases are computed and include a lifting as well as nonlifting cases. The nonlifting solutions obtained are compared to solutions from a Euler code.

Snyder, Aaron↗

NASA NTRS1987

Estimation of geopotential from satellite-to-satellite range rate data: Numerical results

A technique for high-resolution geopotential field estimation by recovering the harmonic coefficients from satellite-to-satellite range rate data is presented and tested against both a controlled analytical simulation of a one-day satellite mission (maximum degree and order 8) and then against a Cowell method simulation of a 32-day mission (maximum degree and order 180). Innovations include: (1) a new frequency-domain observation equation based on kinetic energy perturbations which avoids much of the complication of the usual Keplerian element perturbation approaches; (2) a new method for computing the normalized inclination functions which unlike previous methods is both efficient and numerically stable even for large harmonic degrees and orders; (3) the application of a mass storage FFT to the entire mission range rate history; (4) the exploitation of newly discovered symmetries in the block diagonal observation matrix which reduce each block to the product of (a) a real diagonal matrix factor, (b) a real trapezoidal factor with half the number of rows as before, and (c) a complex diagonal factor; (5) a block-by-block least-squares solution of the observation equation by means of a custom-designed Givens orthogonal rotation method which is both numerically stable and tailored to the trapezoidal matrix structure for fast execution.

Thobe, Glenn E.↗

NASA NTRS1988

Numerical study of delta wing leading edge blowing

Spanwise and tangential leading edge blowing as a means of controlling the position and strength of the leading edge vortices are studied by numerical solution of the three-dimensional Navier-Stokes equations. The leading edge jet is simulated by defining a permeable boundary, corresponding to the jet slot, where suitable boundary conditions are implemented. Numerical results are shown to compare favorably with experimental measurements. It is found that the use of spanwise leading edge blowing at moderate angle of attack magnifies the size and strength of the leading edge vortices, and moves the vortex cores outboard and upward. The increase in lift primarily comes from the greater nonlinear vortex lift. However, spanwise blowing causes earlier vortex breakdown, thus decreasing the stall angle. The effects of tangential blowing at low to moderate angles of attack tend to reduce the pressure peaks associated with leading edge vortices and to increase the suction peak around the leading edge, so that the integrated value of the surface pressure remains about the same. Tangential leading edge blowing in post-stall conditions is shown to re-establish vortical flow and delay vortex bursting, thus increasing C sub L sub max and stall angle.

Yeh, David↗

NASA NTRS1988

Numerical methods for incompressible viscous flows with engineering applications

A numerical scheme has been developed to solve the incompressible, 3-D Navier-Stokes equations using velocity-vorticity variables. This report summarizes the development of the numerical approximation schemes for the divergence and curl of the velocity vector fields and the development of compact schemes for handling boundary and initial boundary value problems.

Rose, M. E.↗

NASA NTRS1988

An improvement in the numerical integration procedure used in the NASA Marshall engineering thermosphere model

A proposed replacement scheme for the integration of the barometric and diffusion equations in the NASA Marshall Engineering Thermosphere (MET) model is presented. This proposed integration scheme is based on Gaussian Quadrature. Extensive numerical testing reveals it to be faster, more accurate and more reliable than the present integration scheme (a modified form of Simpson's Rule) used in the MET model. Numerous graphical examples are provided, along with a listing of a modified form of the MET model in which subroutine INTEGRATE (using Simpson's Rule) is replaced by subroutine GAUSS (which uses Gaussian Quadrature). It is recommended that the Gaussian Quadrature integration scheme, as used here, be used in the MET model.

Hickey, Michael Philip↗

NASA NTRS1987

Numerical methods for one-dimensional reaction-diffusion equations arising in combustion theory

A review of numerical methods for one-dimensional reaction-diffusion equations arising in combustion theory is presented. The methods reviewed include explicit, implicit, quasi-linearization, time linearization, operator-splitting, random walk and finite-element techniques and methods of lines. Adaptive and nonadaptive procedures are also reviewed. These techniques are applied first to solve two model problems which have exact traveling wave solutions with which the numerical results can be compared. This comparison is performed in terms of both the wave profile and computed wave speed. It is shown that the computed wave speed is not a good indicator of the accuracy of a particular method. A fourth-order time-linearized, Hermitian compact operator technique is found to be the most accurate method for a variety of time and space sizes.

Ramos, J. I.↗

NASA NTRS1987

Numerical simulations of the subsynoptic features associated with the AVE-SESAME I case. I - The preconvective environment

The extensive diagnostic calculations made possible by the AVE-SESAME I database are used in combination with numerical simulations from the Mesoscale Atmospheric Simulation System (MASS) model to examine the dynamics of the meso-alpha-scale features during the preconvective period from 1130 to 2030 UTC on April 10, 1979. The version of the MASS model used in this investigation is presented, and an overview of the general synoptic conditions present at the time of model initialization is presented along with the data used to initialize the model. The dynamical processes present in the numerical simulations are presented and compared with analyses of the observational data from this and other investigations of this case. The relative importance of the adiabatic and diabatic processes in creating and then initiating the release of the convective instability is discussed.

Zack, John W.↗

NASA NTRS1987

A comparison of analytical and numerical models for steadily driven magnetic reconnection

The effect of boundary conditions in both analytical and numerical solutions of steadily driven reconnection on the reconnection of antiparallel magnetic fields is considered, with special attention given to the mathematical problem of properly specifying boundary conditions for the MHD equation. A unified formulation developed by Priest and Forbes (1986) for steady state reconnection, which includes the Petschek solution and the Sonnerup (1970) solution as special cases, is used to reinterpret the previous numerical experiments of steadily driven reconnection. It is shown that many contradictory features of these experiments were caused by the use of boundary conditions which are different from those required by Petschek's (1964) theory.

Forbes, T. G.↗

NASA NTRS1988

K-epsilon turbulence model assessment with reduced numerical diffusion for coaxial jets

A numerical study has been conducted to evaluate the performance of the k-epsilon turbulence model for axisymmetric, unconfined, swirling, and nonswirling coaxial jet flows. Two lower order schemes, hybrid and power-law, and two higher order schemes, flux-spline and bounded skew upwind differencing, were employed in this investigation. The predicted results indicate that the higher order numerical schemes have greater potential for future model improvements and complex flow calculations in terms of storage, accuracy, and execution time. For the nonswirling flow, computations using an algebraic stress model have also been presented.

Nikjooy, M.↗

NASA NTRS1987

A numerical study of the effects of curvature and convergence on dilution jet mixing

Holdeman, J. D.↗

NASA NTRS1988

A numerical method for acoustic oscillations in tubes

A numerical method to obtain the neutral curve for the onset of acoustic oscillations in a helium-filled tube is described. Such oscillations can cause a serious heat loss in the plumbing associated with liquid helium dewars. The problem is modelled by a second-order, ordinary differential eigenvalue problem for the pressure perturbation. The numerical method to find the eigenvalues and track the resulting points along the neutral curve is tailored to this problem. The results show that a tube with a uniform temperature gradient along it is much more stable than one where the temperature suddenly jumps from the cold to the hot value in the middle of the tube.

Gary, John M.↗

NASA NTRS1987

Numerical simulation of single and multiple jet interaction within a combustor

Single and multiple jet interaction in a combustor is simulated numerically to determine the effects of various k-epsilon models and boundary conditions as part of an attempt to understand the physics of the flow field inside the preburner of the Space Shuttle Main Engine. Due to the complexity of the physical and chemical phenomena involved, the model is simplified to two dimensions and a single-component nonreacting fluid. The computer code used is the modified TEACH code, which is a steady state isothermal fluid code using the SIMPLE algorithm. The physical validity and numerical accuracy of the model are established through comparisons with experimentally determined isothermal flow fields in a combustor.

Bai, S. Don↗

NASA NTRS1988

Incompressible and anelastic flow simulations on numerically generated grids

In the numerical simulation of incompressible and anelastic flows, it is necessary to solve an elliptic equation at each time step. When the boundaries of such flows are nonrectangular, it may be advantageous to solve the equations on a new, numerically generated coordinate grid, in which the property of orthogonality has been preserved. Flow equations in general curvilinear coordinates maintaining the conservative form are given for both anelastic models using the momentum equations, and for incompressible models, using the vorticity equation. The general problem of grid-generation in two dimensions is presented, and a quasi-conformal transformation technique is discussed in detail. Some examples of grids generated by this technique are exhibited. Three examples of the flow of a stratified fluid over obstacles are presented, in which the grid-generation permits some new results to be obtained.

Sharman, R. D.↗