Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical”

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,549 records · Page 86

Application of similitude principle to the numerical simulation of solar atmospheric dynamics

Numerical simulation has become an essential part of solar physics because the complex nonlinear characteristics of solar phenomena makes analytical solutions difficult to obtain. Realistic simulation of the birth and decay of an active region still is not possible because of the wide range of spatial and time scales that must be considered. Therefore, proper scaling rules must be recognized for the development of appropriate models. In this paper, the similitude principle is applied to develop scaling rules. It is found that these rules are highly dependent on the physical nature of the specific problem under consideration. A set of 'similitude critiques' is presented for some specific physical conditions. Numerical examples of coronal dynamic response and active region dynamics are used to demonstrate these ideas.

Wu, S. T.↗

Numerical experiments on stability and transition in plane channel flow

Various secondary and tertiary instabilities in plane channael flow are explored via time-dependent numerical simulations using the incompressible Navier-Stokes equations. Comparisons are made between transitional flows at Reynolds numbers 1500, 5000, and 8000. The lambda vortex, detached shear layer, and inverted vortex regions are identified, and the origin of the latter is explained. The laminar breakdown of the Re = 1500 flow is computed with high resolution, and the nature of its ensuing hairpin eddies is clarified by numerical particle paths. The potential of center-mode, rather than wall-mode transitions is proposed, and the resulting flow structure is described.

Zang, Thomas A.↗

Numerical simulation of flow through a two-strut scramjet inlet

A three-dimensional, Reynolds-averaged Navier-Stokes code has been used to numerically analyze flow through a two-strut, supersonic combustion ramjet (scramjet) inlet configuration. It solves the governing equations in full conservation form using either a fully explicit or explicit-implicit method. An algebraic, two-layer eddy-viscosity model is used for turbulent flow calculations. The analysis allows inclusion of end effects that are caused by the aft placement of the cowl on the underside of the inlet. A special grid has been developed to accommodate the struts embedded in the inlet flowfield. Detailed numerical results are presented here for the two-strut configuration, and a comparison is made with the available experimental results.

Kumar, Ajay↗

On some numerical schemes for transonic flow problems

New second-order-accurate finite-difference approximations for a class of nonlinear PDEs of mixed type, which includes the two-dimensional low-frequency transonic small-disturbance equation (TSD) and full-potential equation (FP), are presented. For the TSD equation, the scheme is implemented via a time-splitting algorithm; the inclusion of flux limiters keeps the total variation nonincreasing and eliminates spurious oscillations near shocks. Global-linear-stability, TVD, and entropy-stability results are proved. Numerical results for the flow over a thin airfoil are presented. Current techniques used to solve the TSD equation may easily be extended to second-order accuracy by this method. For the FP equation, the new scheme requires no subsonic/supersonic switching and no numerical flux biasing. Global linear stability for all values of the Mach number is proved.

Mostrel, Marco Mosche↗

Numerical solutions of forward-flight rotor flow using an upwind method

A finite-volume upwind algorithm for solving the 3-D Euler equations with a moving grid has been developed for computing helicopter forward-flight rotor flows. The computed pressure distributions and shock positions of high-speed rotor flow are compared with various experimental data as well as with other numerical results, and the agreement is encouraging. A comparison of quasi-steady solutions with unsteady solutions reveals that when a shock occurs in the flowfield, the assumption of quasi-steady flow may fail due to the time-lag of the shock motion. Similarly, three-dimensional effects cannot be neglected. Sufficient subiterations for each time step are required to avoid numerical lag effects in using the present method. The redistribution of the residual due to the coordinate transformation is discussed. For high-order MUSCL-type schemes, a coordinate-independent solution can be obtained by interpolating primitive variables.

Chen, C. L.↗

Numerical simulation of hyperbolic heat conduction with convection boundary conditions and pulse heating effects

The paper describes the numerical simulation of hyperbolic heat conduction with convection boundary conditions. The effects of a step heat loading, a sudden pulse heat loading, and an internal heat source are considered in conjunction with convection boundary conditions. Two methods of solution are presened for predicting the transient behavior of the propagating thermal disturbances. In the first method, MacCormack's predictor-corrector method is employed for integrating the hyperbolic system of equations. Next, the transfinite element method, which employs specially tailored elements, is used for accurately representing the transient response of the propagating thermal wave fronts. The agreement between the results of various numerical test cases validate the representative behavior of the thermal wave fronts. Both methods represent hyperbolic heat conduction behavior by effectively modeling the sharp discontinuities of the propagating thermal disturbances.

Glass, David E.↗

Review of FD-TD numerical modeling of electromagnetic wave scattering and radar cross section

Applications of the finite-difference time-domain (FD-TD) method for numerical modeling of electromagnetic wave interactions with structures are reviewed, concentrating on scattering and radar cross section (RCS). A number of two- and three-dimensional examples of FD-TD modeling of scattering and penetration are provided. The objects modeled range in nature from simple geometric shapes to extremely complex aerospace and biological systems. Rigorous analytical or experimental validatons are provided for the canonical shapes, and it is shown that FD-TD predictive data for near fields and RCS are in excellent agreement with the benchmark data. It is concluded that with continuing advances in FD-TD modeling theory for target features relevant to the RCS problems and in vector and concurrent supercomputer technology, it is likely that FD-TD numerical modeling will occupy an important place in RCS technology in the 1990s and beyond.

Taflove, Allen↗

A Newton/upwind method and numerical study of shock wave/boundary layer interactions

The objective of the paper is two-fold. First, an upwind/central differencing method for solving the steady Navier-Stokes equations is described. The symmetric line relation method is used to solve the resulting algebraic system to achieve high computational efficiency. The grid spacings used in the calculations are determined from the triple-deck theory, in terms of Mach and Reynolds numbers and other flow parameters. Thus the accuracy of the numerical solutions is improved by comparing them with experimental, analytical, and other computational results. Secondly, the shock wave/boundary layer interactions are studied numerically, with special attention given to the flow separation. The concept of free interaction is confirmed. Although the separated region varies with Mach and Reynolds numbers, it is found that the transverse velocity component behind the incident shock, which has not been identified heretofore, is also an important parameter. A small change of this quantity is sufficient to eliminate the flow separation entirely.

Liou, Meng-Sing↗

The numerical solution of the Navier-Stokes equations for 3-dimensional, unsteady, incompressible flows by compact schemes

The present numerical method for the solution of unsteady, incompressible three-dimensional flow Navier-Stokes equations using velocity-vorticity variables and irregular Cartesian grids proceeds by solving: (1) equations of Cauchy-Riemann type for the velocity; and (2) transport-diffusion equations for the vorticity, whose solenoidal vorticity components are generated by a Poisson equation for an appropriate scalar potential. Iterations are used to solve the finite difference equations, facilitating the use of vector and parallel-computing methods; numerical experiments have verified the method's second-order spatial and temporal accuracy.

Gatski, T. B.↗

A numerical analysis of Stefan problems for generalized multi-dimensional phase-change structures using the enthalpy transforming model

An enthalpy transforming scheme is proposed to convert the energy equation into a nonlinear equation with the enthalpy, E, being the single dependent variable. The existing control-volume finite-difference approach is modified so it can be applied to the numerical performance of Stefan problems. The model is tested by applying it to a three-dimensional freezing problem. The numerical results are in agreement with those existing in the literature. The model and its algorithm are further applied to a three-dimensional moving heat source problem showing that the methodology is capable of handling complicated phase-change problems with fixed grids.

Cao, Yiding↗

A numerical and experimental study of coaxial jets

An algebraic stress model and the standard k-epsilon model is applied to predict the mean and turbulence quantities for axisymmetric, nonswirling coaxial jets without confinement. To investigate the effects of numerical (false) diffusion on the predicted results, three different discretization schemes, namely, hybrid, power-law, and the flux-spline, are employed. In addition, an experimental study is conducted to provide data of good quality, especially near the inlet, for model assessment. The results show that the use of the algebraic stress model leads to better agreement between the numerical results and experimental data.

Nikjooy, M.↗

On the ability of rough surface scattering approximations to predict hydrodynamic modulation of the ocean radar cross section - A numerical study

Scattering from the ocean surface at microwave frequencies is often described by the Bragg scattering model or the two-scale model. Here, the ability of these models to predict cross section modulation is investigated by comparing Bragg and two-scale cross-sections with cross sections calculated by a numerical electromagnetic approach. The numerical method is used as the standard of comparison. For slightly rough surfaces, the Bragg model is found to be very accurate, but for much rougher surfaces, the two-scale model is superior and should be used for calculating both the absolute cross section and the change in cross section due to spatial or temporal modulation of the ocean wave height spectrum. These results are relevant to the interpretation of SAR images of ocean surface phenomena which perturb the wave height spectrum at ripple wavelengths.

Durden, Stephen L.↗

Development of a numerical solution to the time dependent kinetic equation

A numerical solution was developed for the time dependent Fokker-Planck equation for arbitrary distributions of electrons injected into a magnetized plasma. The code which includes energy loss and pitch angle scattering due to Coulomb collisions and changes in pitch angle due to inhomogeneous magnetic fields was calibrated and tested. The numerical method is versatile so that other scattering or radiation terms can be easily included. Using this code many processes associated with the impulsive phase of solar flares will be investigated.

Hamilton, Russell J.↗

Numerical solutions of the compressible 3-D boundary-layer equations for aerospace configurations with emphasis on LFC

The application of stability theory in Laminar Flow Control (LFC) research requires that density and velocity profiles be specified throughout the viscous flow field of interest. These profile values must be as numerically accurate as possible and free of any numerically induced oscillations. Guidelines for the present research project are presented: develop an efficient and accurate procedure for solving the 3-D boundary layer equation for aerospace configurations; develop an interface program to couple selected 3-D inviscid programs that span the subsonic to hypersonic Mach number range; and document and release software to the LFC community. The interface program was found to be a dependable approach for developing a user friendly procedure for generating the boundary-layer grid and transforming an inviscid solution from a relatively coarse grid to a sufficiently fine boundary-layer grid. The boundary-layer program was shown to be fourth-order accurate in the direction normal to the wall boundary and second-order accurate in planes parallel to the boundary. The fourth-order accuracy allows accurate calculations with as few as one-fifth the number of grid points required for conventional second-order schemes.

Harris, Julius E.↗

Numerical experiments on transition control in wall-bounded shear flows

Results are presented from a numerical simulation of transition control in plane channel and boundary layer flows. The analysis is based on a pseudo-spectral/finite difference semi-implicit solution procedure employed to numerically integrate the time-dependent, three-dimensional, incompressible Navier-Stokes equations in a doubly periodic domain. In the channel flow, it was found that the active periodic suction/blowing method was effective in controlling strongly three-dimensional disturbances. In the boundary layer, the preliminary analysis indicated that in the early stages, passive control by suction is as effective as active control to suppress instabilities. The current work is focused on a detailed comparison of active and passive control by suction/blowing in the boundary layer.

Biringen, S.↗

Numerical investigations in three-dimensional internal flows

In order to develop an understanding of flow fields that exist in sidewall compression inlet models, an explicit, time-accurate full Navier-Stokes code was used to model them. Both laminar and turbulent boundary layers were numerically simulated. These simulations were examined as they were influenced by Mach number, changes in sidewall compression angle, hypersonic viscous interaction effects, and thick entering boundary layers. Cause and effect relationships were established. The numerical simulations were compared with the limited data available to aid in understanding both the simulations and the experiments.

Rose, William C.↗

Numerical solutions of the linearized Euler equations for unsteady vortical flows around lifting airfoils

A linearized unsteady aerodynamic analysis is presented for unsteady, subsonic vortical flows around lifting airfoils. The analysis fully accounts for the distortion effects of the nonuniform mean flow on the imposed vortical disturbances. A frequency domain numerical scheme which implements this linearized approach is described, and numerical results are presented for a large variety of flow configurations. The results demonstrate the effects of airfoil thickness, angle of attack, camber, and Mach number on the unsteady lift and moment of airfoils subjected to periodic vortical gusts. The results show that mean flow distortion can have a very strong effect on the airfoil unsteady response, and that the effect depends strongly upon the reduced frequency, Mach number, and gust wave numbers.

Scott, James R.↗

Materials properties numerical database system established and operational at CINDAS/Purdue University

A computerized comprehensive numerical database system on the mechanical, thermophysical, electronic, electrical, magnetic, optical, and other properties of various types of technologically important materials such as metals, alloys, composites, dielectrics, polymers, and ceramics has been established and operational at the Center for Information and Numerical Data Analysis and Synthesis (CINDAS) of Purdue University. This is an on-line, interactive, menu-driven, user-friendly database system. Users can easily search, retrieve, and manipulate the data from the database system without learning special query language, special commands, standardized names of materials, properties, variables, etc. It enables both the direct mode of search/retrieval of data for specified materials, properties, independent variables, etc., and the inverted mode of search/retrieval of candidate materials that meet a set of specified requirements (which is the computer-aided materials selection). It enables also tabular and graphical displays and on-line data manipulations such as units conversion, variables transformation, statistical analysis, etc., of the retrieved data. The development, content, accessibility, etc., of the database system are presented and discussed.

Ho, C. Y.↗