SEARCH · Search NASA

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 253 records · Page 14

All sources NASA NTRS DOE OSTI

NASA NTRS1986

Numerical analysis of some supersonic viscous flows related to inlet and nozzle systems

A numerical method originally developed for three-dimensional supersonic inlet flow calculations is extended and applied to the study of two-dimensional and three-dimensional flows associated with arbitrary propulsion systems. The method is based on the forward spatial marching solution of a reduced form of the three-dimensional steady Navier-Stokes equations in which streamwise pressure gradients are retained in both the subsonic and supersonic regions. The present paper briefly describes the analysis and then shows three applications. In the first application, a wall transpiration study has been performed for the two-dimensional shock wave/turbulent boundary layer interaction flow field with application to an inlet configuration. The second application treats a rectangular high speed inlet with a swept sideplate including the effects of sideplate spillage. Finally, the method is utilized to analyze the interaction of an under-expanded supersonic jet with an ambient flow. Computed results are examined and compared with available experimental measurements. It is demonstrated that the present numerical method is capable of numerically simulating complex two- and three-dimensional flows relevant to hypersonic propulsion systems in a manner which both shows good agreement with data when such data is available, and which shows the complex flow features in the absence of data.

Kim, Y.-N.↗

NASA NTRS1984

Overview: Applications of numerical optimization methods to helicopter design problems

There are a number of helicopter design problems that are well suited to applications of numerical design optimization techniques. Adequate implementation of this technology will provide high pay-offs. There are a number of numerical optimization programs available, and there are many excellent response/performance analysis programs developed or being developed. But integration of these programs in a form that is usable in the design phase should be recognized as important. It is also necessary to attract the attention of engineers engaged in the development of analysis capabilities and to make them aware that analysis capabilities are much more powerful if integrated into design oriented codes. Frequently, the shortcoming of analysis capabilities are revealed by coupling them with an optimization code. Most of the published work has addressed problems in preliminary system design, rotor system/blade design or airframe design. Very few published results were found in acoustics, aerodynamics and control system design. Currently major efforts are focused on vibration reduction, and aerodynamics/acoustics applications appear to be growing fast. The development of a computer program system to integrate the multiple disciplines required in helicopter design with numerical optimization technique is needed. Activities in Britain, Germany and Poland are identified, but no published results from France, Italy, the USSR or Japan were found.

Miura, H.↗

NASA NTRS1986

Optimization methods and silicon solar cell numerical models

An optimization algorithm for use with numerical silicon solar cell models was developed. By coupling an optimization algorithm with a solar cell model, it is possible to simultaneously vary design variables such as impurity concentrations, front junction depth, back junction depth, and cell thickness to maximize the predicted cell efficiency. An optimization algorithm was developed and interfaced with the Solar Cell Analysis Program in 1 Dimension (SCAP1D). SCAP1D uses finite difference methods to solve the differential equations which, along with several relations from the physics of semiconductors, describe mathematically the performance of a solar cell. A major obstacle is that the numerical methods used in SCAP1D require a significant amount of computer time, and during an optimization the model is called iteratively until the design variables converge to the values associated with the maximum efficiency. This problem was alleviated by designing an optimization code specifically for use with numerically intensive simulations, to reduce the number of times the efficiency has to be calculated to achieve convergence to the optimal solution.

Girardini, K.↗

NASA NTRS1987

Numerical simulation of shock wave diffraction by TVD schemes

An upwind total variation diminishing (TVD) scheme and a predictor-corrector symmetric TVD scheme were used to numerically simulate the blast wave diffraction on a stationary object. The objective is to help design an optimum configuration so that lateral motion is minimized and at the same time vortex shedding and flow separation are reduced during a blast wave encounter. Results are presented for a generic configuration for both a coarse grid and a fine grid to illustrate the global and local diffraction flow fields. Numerical experiments for the shock wave reflection on a wedge are also included to validate the current approach. Numerical study indicated that these TVD schemes are more stable and produced higher shock resolution than classical shock capturing methods such as the explicit MacCormack scheme.

Young, Victor Y. C.↗

NASA NTRS1987

Vortical flow aerodynamics - Physical aspects and numerical simulation

Progress in the numerical simulation of vortical flow due to three-dimensional flow separation about flight vehicles at high angles of attack and quasi-steady flight conditions is surveyed. Primary emphasis is placed on Euler and Reynolds-averaged Navier-Stokes methods where the vortices are 'captured' as a solution to the governing equations. A discussion of the relevant flow physics provides a perspective from which to assess numerical solutions. Current numerical prediction capabilities and their evolutionary development are surveyed. Future trends and challenges are identified and discussed.

Newsome, Richard W.↗

NASA NTRS1987

The use of a numerical filter to correct airborne temperature measurements for the effects of sensor lag

A numerical filter for transforming measured temperature signals into a close approximation of the actual temperature signal is described. The filter is derived by minimizing the mean-square error of the system, and assuming a knowledge of the characteristics of the sensing element and its housing. The equation representing the frequency-response function of the numerical filter is given. Input and output spectra for a filter applied to a case with negligible noise and a noise level of 1.5 percent of the total power in the input spectrum are analyzed, and the numerical weights for these two cases are calculated. Phase angle and gain for the entire system are examined. It is noted that the filter can enhance spectral components as high as 8 Hz with little phase and gain degradation over the bandwidth.

Ritter, John A.↗

NASA NTRS1985

Three-dimensional convection in horizontal cylinders - Numerical solutions and comparison with experimental and analytical results

The results of a three-dimensional numerical simulation of Boussinesq free convection in a horizontal differentially heated cylinder are presented. The computation was based on a Samarskii-Andreyev scheme (described by Leong, 1981) and a false-transient advancement in time, with vorticity, velocity, and temperature as dependent variables. Solutions for velocity and temperature distributions were obtained for Rayleigh numbers (based on the radius) Ra = 74-18,700, thus covering the core- and boundary-layer-driven regimes. Numerical solutions are compared with asymptotic analytical solutions and experimental data. The numerical results well represent the complex three-dimensional flows found experimentally.

Smutek, C.↗

NASA NTRS1987

Mesoscale lake-effect snowstorms in the vicinity of Lake Michigan - Linear theory and numerical simulations

Mesoscale lake-effect snowstorms in the vicinity of Lake Michigan are studied by a linear steady-state analytic model and a nonlinear time-dependent numerical model with parameterized subgrid-scale physics. The solutions of the linear model show that the orientation of the mean wind field to the surface heating pattern is crucial to the shapes of the disturbances. The results indicate that the relative warmth of the lake surface can induce three updraft centers under a westerly wind, two updraft centers/bands under a northwesterly wind, and a convergence band under a northerly wind. Such convergences are caused by the interaction between the mean wind and the local circulations forced by the curved thermal contrasts. The numerical results from the nonlinear model not only produce convergence patterns very similar to those found in the linear theory and in other numerical studies, but also capture the transient property in some of the lake-induced disturbances. All of these results are qualitatively confirmed by satellite images.

Hsu, Hsiao-Ming↗

NASA NTRS1987

The spatially nonuniform convergence of the numerical solution of flows

The spatial distribution of the numerical disturbances that are generated during the numerical solution of a flow is examined. It is shown that the distribution of the disturbances is not uniform. In regions where the structure of a flow is simple, the magnitude of the generated disturbances is small and their decay is fast. However, in complex flow regions, as in separation and vortical areas, large magnitude disturbances appear and their decay may be very slow. The observed nonuniformity of the numerical disturbances makes possible the reduction of the calculation time by application of what may be called the partial-grid calculation technique, in which a major part of the calculation procedure is applied in selective subregions, where the velocity disturbances are large, and not within the whole grid. This technique is expected to prove beneficial in large-scale calculations such as the flow about complete aircraft configurations at high angle of attack. Also, it has been shown that if the Navier-Stokes equations are written in a generalized coordinate system, then in regions in which the grid is fine, such as near solid boundaries, the norms become infinitesimally small, because in these regions the Jacobian has very large values. Thus, the norms, unless they are unscaled by the Jacobians, reflect only the changes that happen at the outer boundaries of the computation domain, where the value of the Jacobian approaches unity, and not in the whole flow field.

Panaras, Argyris G.↗

NASA NTRS1986

Numerical experiments with a symmetric high-resolution shock-capturing scheme

Characteristic-based explicit and implicit total variation diminishing (TVD) schemes for the two-dimensional compressible Euler equations have recently been developed. This is a generalization of recent work of Roe and Davis to a wider class of symmetric (non-upwind) TVD schemes other than Lax-Wendroff. The Roe and Davis schemes can be viewed as a subset of the class of explicit methods. The main properties of the present class of schemes are that they can be implicit, and, when steady-state calculations are sought, the numerical solution is independent of the time step. In a recent paper, a comparison of a linearized form of the present implicit symmetric TVD scheme with an implicit upwind TVD scheme originally developed by Harten and modified by Yee was given. Results favored the symmetric method. It was found that the latter is just as accurate as the upwind method while requiring less computational effort. Currently, more numerical experiments are being conducted on time-accurate calculations and on the effect of grid topology, numerical boundary condition procedures, and different flow conditions on the behavior of the method for steady-state applications. The purpose here is to report experiences with this type of scheme and give guidelines for its use.

Yee, H. C.↗

NASA NTRS1987

A numerical study of transient heat and mass transfer in crystal growth

A numerical analysis of transient heat and solute transport across a rectangular cavity is performed. Five nonlinear partial differential equations which govern the conservation of mass, momentum, energy and solute concentration related to crystal growth in solution, are simultaneously integrated by a numerical method based on the SIMPLE algorithm. Numerical results showed that the flow, temperature and solute fields are dependent on thermal and solutal Grashoff number, Prandtl number, Schmidt number and aspect ratio. The average Nusselt and Sherwood numbers evaluated at the center of the cavity decrease markedly when the solutal buoyancy force acts in the opposite direction to the thermal buoyancy force. When the solutal and thermal buoyancy forces act in the same direction, however, Sherwood number increases significantly and yet Nusselt number decreases. Overall effects of convection on the crystal growth are seen to be an enhancement of growth rate as expected but with highly nonuniform spatial growth variations.

Han, Samuel Bang-Moo↗

NASA NTRS1987

Numerical and analytical study of fluid dynamic forces in seals and bearings

A numerical model based on a transformed, conservative form of the three dimensional Navier-Stokes equation and an analytical model based on lumped fluid parameters are presented and compared with studies of modeled rotor bearing seal systems. The rotor destabilizing factors are related to the rotative character of the flow field. It is shown that these destabilizing factors can be reduced through a decrease in the fluid average circumferential velocity. However, the rotative character of the flow field is a complex three dimensional system with bifurcated secondary flow patterns that significantly alter the fluid circumferential velocity. By transforming the Navier-Stokes equations to those for a rotating observer and using the numerical code PHOENICS-84 with a nonorthogonal body fitted grid, several numerical experiments were carried out to demonstrate the character of this complex flow field. In general, fluid injection and/or preswirl of the flow field opposing the shaft rotation significantly intensified these secondary recirculation zones and thus reduced the average circumferential velocity; injection or preswirl in the direction of rotation significantly weakened these zones.

Tam, L. T.↗

NASA NTRS1987

Numerical simulations of aerodynamic contribution of flows about a space-plane-type configuration

The slightly supersonic viscous flow about the space-plane under development at the National Aerospace Laboratory (NAL) in Japan was simulated numerically using the LU-ADI algorithm. The wind-tunnel testing for the same plane also was conducted with the computations in parallel. The main purpose of the simulation is to capture the phenomena which have a great deal of influence to the aerodynamic force and efficiency but is difficult to capture by experiments. It includes more accurate representation of vortical flows with high angles of attack of an aircraft. The space-plane shape geometry simulated is the simplified model of the real space-plane, which is a combination of a flat and slender body and a double-delta type wing. The comparison between experimental results and numerical ones will be done in the near future. It could be said that numerical results show the qualitatively reliable phenomena.

Matsushima, Kisa↗

NASA NTRS1987

A numerical simulation of vortex breakdown

A numerical simulation of vortex breakdown using the time-dependent Navier-Stokes equations is performed. Unlike previous studies, the numerical algorithm, formulated in terms of the velocity and vorticity, is not restricted by axisymmetry conditions. The vortex is parameterized in terms of the Reynolds number and Rossby number. The resulting breakdown structure is analyzed using contour plots of velocity, vorticity and pressure as well as axial, radial, and swirl velocity profiles at various streamwise locations. The relationship of these results to experimentally observed structures and previous numerical results is discussed.

Spall, R. E.↗

NASA NTRS1988

A Three-Dimensional, Time-Dependent Numerical Modeling of Super-Sonic, Super-Alfvenic MHD Flow

A three-dimensional, time-dependent, MHD disturbance propagating through a steady inhomogeneous MHD flow is analyzed by a numerical model based on the Lax-Wendroff finite difference method. The initial steady MHD flow is assumed to be supersonic and super-Alfvenic. The background steady MHD flow is obtained by dropping the time-dependent terms in the unsteady MHD equations. A transient MHD disturbance corresponding to a MHD shock is generated at the lower radial boundary surface. The resulting development and propagation of 3-D MHD disturbances interacting with the background steady MHD flow are numerically simulated in this first application to the field of interplanetary physics. Detailed descriptions of the mathematical formulation and numerical methodology are presented. Example computations are performed to demonstrate the validity of the formulation and the capabilities of the model.

Han, S. M.↗

NASA NTRS1987

Numerical analysis of discrete fractional integrodifferential structural dampers

This paper develops solution algorithms enabling the handling of the dynamic response of nonlinear structures contained discretely attached dampers modeled by fractional integrodifferential operators of the Grunwald-Liouville-Riemann type. The development consists of two levels of formulation, namely: (1) numerical approximations of fractional operators and, (2) the establishment of global level implicit schemes enabling the solution to nonlinear structural formulations. To generalize the overall results, error estimates are derived for the fractional operator approximation algorithm. These enable an ongoing optimization of solution efficiency for a given error tolerance. To benchmark the scheme, the results of several numerical experiments are presented. These illustrate the numerical characteristics of the overall formulation.

Padovan, Joseph↗

NASA NTRS1988

Supersonic combusting flow - A comparison of numerical and experimental results

The flow field of an unconfined coaxial supersonic burner (SSB) was obtained numerically, using a CFD code. The explicit MacCormack algorithm, with point implicit treatment of chemistry source terms, was employed in the CFD code. Combustion of hydrogen and air was simulated by a two-step finite-rate combustion model whereas turbulence was accounted for by a Prandtl mixing length scheme. Boundary conditions and some important features of the numerical scheme were discussed. The qualitative features of the numerical solution were compared with those obtained by regular as well as schlieren photographs. General agreements of the CFD solution and the experimental photographs were observed.

Chitsomboon, T.↗

NASA NTRS1987

Numerical and analytical study of fluid dynamic forces in seals and bearings

Study results for rotor/bearing/seal system models are compared with both a numerical model based on a transformed, conservative version of the three-dimensional Navier-Stokes equation, and an analytical model based on 'lumped' fluid parameters. Rotor-destabilizing factors are related to flow field rotation; it is shown that these factors can be reduced through a decrease in the fluid average circumferential velocity. A transformation of the Navier-Stokes equations to those for a rotating observer and a numerical code with a nonorthogonal body-fitted grid are used to carry out numerical experiments demonstrating the character of the complex flow field.

Tam, L. T.↗