Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical computation”

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,567 records · Page 87

Rapid calculation of radiative heating rates and photodissociation rates in inhomogeneous multiple scattering atmospheres

The solution of the generalized two-stream approximation for radiative transfer in homogeneous multiple scattering atmospheres is extended to vertically inhomogeneous atmospheres in a manner which is numerically stable and computationally efficient. It is shown that solar energy deposition rates, photolysis rates, and infrared cooling rates all may be calculated with the simple modifications of a single algorithm. The accuracy of the algorithm is generally better than 10 percent, so that other uncertainties, such as in absorption coefficients, may often dominate the error in calculation of the quantities of interest to atmospheric studies.

Toon, Owen B.↗

A complete spherical harmonic approach to luni-solar tides

In this work a spherical harmonic theory of ocean tides is presented. The theory is based on Laplace tide equations modified to include turbulence with constant eddy viscosity, linearized bottom friction, and oceanic loading and self-gravitation. Variable bathymetry is also treated in harmonic terms, and no-flow boundary conditions are applied at continental coastlines. The tide and boundary constraint equations are reduced to matrix form and solved by a weighted least-squares procedure. Five zonal luni-solar tides, ranging in period from 14 days to 18.6 yr, are investigated using the theory; such tides have typically been difficult to compute using traditional numerical approaches. The polar motion and changes in the length of day induced by these long-period tides are calculated. Tidal solutions are compared extensively with results from other tidal theories and from recent satellite and sea-level observations. The greatest limitation to accurate prediction of zonal tides - for any theory - appears to be the marginal failure of all tide theories to conserve mass globally; the use of additional mass constraints may be warranted.

Dickman, S. R.↗

Numerical study of vortical flow over a sideslipping delta wing

The three-dimensional vortical viscous flow past a sideslipping delta wing at a high angle of attack is numerically investigated. A computational method for calculating low-speed viscous flowfields is developed. The time-iterative method uses an implicit upwind-relaxation finite-difference algorithm with a nonsingular eigensystem to solve the preconditioned, three-dimensional, incompressible Navier-Stokes equations in curvilinear coordinates. An algebraic turbulence model is implemented to account for the eddy viscosity. The technique of local time stepping is incorporated to accelerate the rate of convergence to a steady-state solution. Computed results are compared with experimental data.

Hsu, C.-H.↗

Probabilistic structural analysis of a truss typical for space station

A three-bay, space, cantilever truss is probabilistically evaluated using the computer code NESSUS (Numerical Evaluation of Stochastic Structures Under Stress) to identify and quantify the uncertainties and respective sensitivities associated with corresponding uncertainties in the primitive variables (structural, material, and loads parameters) that defines the truss. The distribution of each of these primitive variables is described in terms of one of several available distributions such as the Weibull, exponential, normal, log-normal, etc. The cumulative distribution function (CDF's) for the response functions considered and sensitivities associated with the primitive variables for given response are investigated. These sensitivities help in determining the dominating primitive variables for that response.

Pai, Shantaram S.↗

Emission from small dust particles in diffuse and molecular cloud medium

Infrared Astronomy Satellite (IRAS) observations of the whole galaxy has shown that long wavelength emission (100 and 60 micron bands) can be explained by thermal emission from big grains (approx 0.1 micron) radiating at their equilibrium temperature when heated by the InterStellar Radiation Field (ISRF). This conclusion has been confirmed by continuum sub-millimeter observations of the galactic plane made by the EMILIE experiment at 870 microns (Pajot et al. 1986). Nevertheless, shorter wavelength observations like 12 and 25 micron IRAS bands, show an emission from the galactic plane in excess with the long wavelength measurements which can only be explained by a much hotter particles population. Because dust at equilibrium cannot easily reach high temperatures required to explain this excess, this component is thought to be composed of very small dust grains or big molecules encompassing thermal fluctuations. Researchers present here a numerical model that computes emission, from Near Infrared Radiation (NIR) to Sub-mm wavelengths, from a non-homogeneous spherical cloud heated by the ISRF. This model fully takes into account the heating of dust by multi-photon processes and back-heating of dust in the Visual/Infrared Radiation (VIS-IR) so that it is likely to describe correctly emission from molecular clouds up to large A sub v and emission from dust experiencing temperature fluctuations. The dust is a three component mixture of polycyclic aromatic hydrocarbons, very small grains, and classical big grains with independent size distributions (cut-off and power law index) and abundances.

Bernard, J. P.↗

Low Reynolds number two-equation modeling of turbulent flows

A k-epsilon model that accounts for viscous and wall effects is presented. The proposed formulation does not contain the local wall distance thereby making very simple the application to complex geometries. The formulation is based on an existing k-epsilon model that proved to fit very well with the results of direct numerical simulation. The new form is compared with nine different two-equation models and with direct numerical simulation for a fully developed channel flow at Re = 3300. The simple flow configuration allows a comparison free from numerical inaccuracies. The computed results prove that few of the considered forms exhibit a satisfactory agreement with the channel flow data. The model shows an improvement with respect to the existing formulations.

Michelassi, V.↗

Development of UMARC (University of Maryland Advanced Rotorcraft Code)

The University of Maryland Advanced Rotorcraft Code (UMARC) is a user-friendly, FEM-based comprehensive helicopter rotor simulation code of high numerical robustness and computational efficiency. UMARC formulates the rotor-fuselage equations using Hamilton's principle, and are discretized using finite elements in space and time. The FEM formulation allows the code to analyze a wide variety of rotor designs. Dynamic inflow modeling is used for unsteady wake inflow computations. Predicted stability, response, and blade-load data are validated with experimental data for several configurations, including representative articulated, hingeless, and bearingless rotors.

Bir, Gunjit↗

Turbulence modeling in aerodynamic shear flows - Status and problems

This paper briefly summarizes the status and problems of turbulence modeling for aerodynamical applications. For complex flows the 'approach of choice' is (increasingly) full second-order (Reynolds stress equation) closure. These closures have not yet developed to anywhere near their full potential, significant further research is required especially regarding length-scale equations, representation of pressure-strain correlations, and wall region treatments. Recent developments in computer capability, algorithms, numerical simulations, theory and quantitative flow visualization should assist in and hasten this research. Several problem areas such as shock interaction and discrete dynamic instabilities of turbulent flows may require mega-to-large eddy simulation or theoretical adjuncts.

Bushnell, D. M.↗

Periodic response of multi-disk rotors with bearing clearances

The forced steady state response of a multi-disk rotor system involving a clearance at one of the bearings is determined by using a harmonic balance approach. The impedance method is applied to each of the harmonic steady state components in order to reduce the system to its displacement at the nonlinear bearing support. The results reveal the interrelated roles of the bearing clearance, mass eccentricity and side force in producing dangerous subharmonics. The significant effects of the strong nonlinearity of a bearing clearance are studied as related to the various system parameters. The results show that the approach developed in this study is computationally superior to numerical integration methods in analyzing multi-disk rotor systems with strong nonlinearity.

Kim, Y. B.↗

Effects of measurement errors on microwave antenna holography

The effects of measurement errors appearing during the implementation of the microwave holographic technique are investigated in detail, and many representative results are presented based on computer simulations. The numerical results are tailored for cases applicable to the utilization of the holographic technique for the NASA's Deep Space Network antennas, although the methodology of analysis is applicable to any antenna. Many system measurement topics are presented and summarized.

Rochblatt, David J.↗

A numerical and experimental study of tangential jet blowing applied to bodies at high angles of attack

The effects of tangential jet blowing on the vortical structures developing about the fuselage of a missile or aircraft at high angle of attack are investigated numerically and experimentally. Computations are carried out for fully laminar and fully turbulent flows, and the results are compared with flow visualization data. It is shown that turbulent computations are more accurate in describing the flow behavior. The actuator plane method is shown to be sufficiently accurate in modeling the essential flow physics of a tangential jet.

Font, G. I.↗

Fractional-Step, Finite-Volume Computation Of Flow

Method of solving Navier-Stokes equations of incompressible flow in general nonorthogonal curvilinear coordinates incorporates fractional-step and finite-volume approaches. Developed in continuing effort to achieve accuracy without need for excessive computation time in numerical simulation of time-dependent, three-dimensional flows bounded by surfaces of complicated shape.

Kwak, Dochan↗

Rotor-to-stator Partial Rubbing and Its Effects on Rotor Dynamic Response

Results from experimental and analytical studies on rotor to stationary element partial rubbings at several locations and their effects on rotor dynamic responses are presented. The mathematical model of a rubbing rotor is given. The computer program provides numerical results which agree with experimentally obtained rotor responses.

Muszynska, Agnes↗

Linear and nonlinear dynamic analysis by boundary element method

An advanced implementation of the direct boundary element method (BEM) applicable to free-vibration, periodic (steady-state) vibration and linear and nonlinear transient dynamic problems involving two and three-dimensional isotropic solids of arbitrary shape is presented. Interior, exterior, and half-space problems can all be solved by the present formulation. For the free-vibration analysis, a new real variable BEM formulation is presented which solves the free-vibration problem in the form of algebraic equations (formed from the static kernels) and needs only surface discretization. In the area of time-domain transient analysis, the BEM is well suited because it gives an implicit formulation. Although the integral formulations are elegant, because of the complexity of the formulation it has never been implemented in exact form. In the present work, linear and nonlinear time domain transient analysis for three-dimensional solids has been implemented in a general and complete manner. The formulation and implementation of the nonlinear, transient, dynamic analysis presented here is the first ever in the field of boundary element analysis. Almost all the existing formulation of BEM in dynamics use the constant variation of the variables in space and time which is very unrealistic for engineering problems and, in some cases, it leads to unacceptably inaccurate results. In the present work, linear and quadratic isoparametric boundary elements are used for discretization of geometry and functional variations in space. In addition, higher order variations in time are used. These methods of analysis are applicable to piecewise-homogeneous materials, such that not only problems of the layered media and the soil-structure interaction can be analyzed but also a large problem can be solved by the usual sub-structuring technique. The analyses have been incorporated in a versatile, general-purpose computer program. Some numerical problems are solved and, through comparisons with available analytical and numerical results, the stability and high accuracy of these dynamic analysis techniques are established.

Ahmad, Shahid↗

Modeling of compressible turbulent shear flows

Despite all the recent developments in computer technologies and numerical algorithms, full numerical simulations of turbulent flows are feasible only at moderate Reynolds numbers and for flows with relatively simple geometries. The main goal of this research is to develop new second order moment closures for compressible turbulence. It has been shown that the models based on the extension of those developed originally for incompressible flows fail to adequately predict turbulent flows at high Mach numbers. In this attempt, the compressibility effects are explicitly considered. A successful development of these models that directly takes into account the compressibility effects may have a range of technological implications in the design of supersonic and hypersonic vehicles.

Liou, William W.↗

Development of new flux splitting schemes

Maximizing both accuracy and efficiency has been the primary objective in designing a numerical algorithm for computational fluid dynamics (CFD). This is especially important for solutions of complex three dimensional systems of Navier-Stokes equations which often include turbulence modeling and chemistry effects. Recently, upwind schemes have been well received for their capability in resolving discontinuities. With this in mind, presented are two new flux splitting techniques for upwind differencing. The first method is based on High-Order Polynomial Expansions (HOPE) of the mass flux vector. The second new flux splitting is based on the Advection Upwind Splitting Method (AUSM). The calculation of the hypersonic conical flow demonstrates the accuracy of the splitting in resolving the flow in the presence of strong gradients. A second series of tests involving the two dimensional inviscid flow over a NACA 0012 airfoil demonstrates the ability of the AUSM to resolve the shock discontinuity at transonic speed. A third case calculates a series of supersonic flows over a circular cylinder. Finally, the fourth case deals with tests of a two dimensional shock wave/boundary layer interaction.

Liou, Meng-Sing↗

A parallel-vector equation solver for unsymmetric matrices on supercomputers

A parallel-vector unsymmetric equation solver is presented. The solver exploits both vector and parallel capabilities provided by modern, high-performance supercomputers. A special storage scheme and loop-unrolling technique are used to optimize the vector performance. A parallel FORTRAN language is used to develop the solver on the CRAY 2 and CRAY Y-MP multiple processing computer environment. Three numerical examples are presented which demonstrate the efficiency and accuracy of this equation solver. The first two examples demonstrate the improved performance, and the third example utilizes the proposed solver to solve a highly nonlinear, unsymmetric finite element formulation for panel flutter.

Qin, J.↗

Minimizing distortion in truss structures - A Hopfield network solution

Distortions in truss structures can result from random errors in element lengths that are typical of a manufacturing process. These distortions may be minimized by an optimal selection of elements from those available for placement between the prescribed nodes - a combinatorial optimization problem requiring significant investment of computational resource for all but the smallest problems. The present paper describes a formulation in which near-optimal element assignments are obtained as minimum-energy stable states, of an analogous Hopfield neural network. This requires mapping of the optimization problem into an energy function of the appropriate Liapunov form. The computational architecture is ideally suited to a parallel processor implementation and offers significant savings in computational effort. A numerical implementation of the approach is discussed with reference to planar truss problems.

Fu, B.↗