Engineering Papers⌕ Search

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.

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

At least 1,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

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.↗
NASA NTRS1989

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.↗
NASA NTRS1989

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↗
NASA NTRS1989

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↗
NASA NTRS1989

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.↗
NASA NTRS1989

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↗
NASA NTRS1989

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.↗
NASA NTRS1989

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.↗
NASA NTRS1989

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.↗
NASA NTRS1987

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.↗
NASA NTRS1987

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.↗
NASA NTRS1988

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.↗
NASA NTRS1990

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.↗
NASA NTRS1989

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.↗
NASA NTRS1989

A PBL-radiation model for application to regional numerical weather prediction

Often in the short-range limited-area numerical weather prediction (NWP) of extratropical weather systems the effects of planetary boundary layer (PBL) processes are considered secondarily important. However, it may not be the case for the regional NWP of mesoscale convective systems over the arid and semi-arid highlands of the southwestern and south-central United States in late spring and summer. Over these dry regions, the PBL can grow quite high up into the lower middle troposphere (600 mb) due to very effective solar heating and hence a vigorous air-land thermal interaction can occur. The interaction representing a major heat source for regional dynamical systems can not be ignored. A one-dimensional PBL-radiation model was developed. The model PBL consists of a constant-flux surface layer superposed with a well-mixed (Ekman) layer. The vertical eddy mixing coefficients for heat and momentum in the surface layer are determined according to the surface similarity theory, while their vertical profiles in the Ekman layer are specified with a cubic polynomial. Prognostic equations are used for predicting the height of the nonneutral PBL. The atmospheric radiation is parameterized to define the surface heat source/sink for the growth and decay of the PBL. A series of real-data numerical experiments has been carried out to obtain a physical understanding how the model performs under various atmospheric and surface conditions. This one-dimensional model will eventually be incorporated into a mesoscale prediction system. The ultimate goal of this research is to improve the NWP of mesoscale convective storms over land.

Chang, Chia-Bo↗
NASA NTRS1990

Thermal radiation characteristics of nonisothermal cylindrical enclosures using a numerical ray tracing technique

Analysis of energy emitted from simple or complex cavity designs can lead to intricate solutions due to nonuniform radiosity and irradiation within a cavity. A numerical ray tracing technique was applied to simulate radiation propagating within and from various cavity designs. To obtain the energy balance relationships between isothermal and nonisothermal cavity surfaces and space, the computer code NEVADA was utilized for its statistical technique applied to numerical ray tracing. The analysis method was validated by comparing results with known theoretical and limiting solutions, and the electrical resistance network method. In general, for nonisothermal cavities the performance (apparent emissivity) is a function of cylinder length-to-diameter ratio, surface emissivity, and cylinder surface temperatures. The extent of nonisothermal conditions in a cylindrical cavity significantly affects the overall cavity performance. Results are presented over a wide range of parametric variables for use as a possible design reference.

Baumeister, Joseph F.↗
NASA NTRS1989

Numerical algorithms for computations of feedback laws arising in control of flexible systems

Several continuous models will be examined, which describe flexible structures with boundary or point control/observation. Issues related to the computation of feedback laws are examined (particularly stabilizing feedbacks) with sensors and actuators located either on the boundary or at specific point locations of the structure. One of the main difficulties is due to the great sensitivity of the system (hyperbolic systems with unbounded control actions), with respect to perturbations caused either by uncertainty of the model or by the errors introduced in implementing numerical algorithms. Thus, special care must be taken in the choice of the appropriate numerical schemes which eventually lead to implementable finite dimensional solutions. Finite dimensional algorithms are constructed on a basis of a priority analysis of the properties of the original, continuous (infinite diversional) systems with the following criteria in mind: (1) convergence and stability of the algorithms and (2) robustness (reasonable insensitivity with respect to the unknown parameters of the systems). Examples with mixed finite element methods and spectral methods are provided.

Lasiecka, Irena↗
NASA NTRS1989

Numerical solution of large Lyapunov equations

A few methods are proposed for solving large Lyapunov equations that arise in control problems. The common case where the right hand side is a small rank matrix is considered. For the single input case, i.e., when the equation considered is of the form AX + XA(sup T) + bb(sup T) = 0, where b is a column vector, the existence of approximate solutions of the form X = VGV(sup T) where V is N x m and G is m x m, with m small is established. The first class of methods proposed is based on the use of numerical quadrature formulas, such as Gauss-Laguerre formulas, applied to the controllability Grammian. The second is based on a projection process of Galerkin type. Numerical experiments are presented to test the effectiveness of these methods for large problems.

Saad, Youcef↗