Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical techniques”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Lubrication of nonconformal contacts

Minimum film thickness results for piezoviscous-rigid regime of lubrication are developed for a compressible Newtonian fluid with Roelands viscosity. The results provide a basis for the analysis and design of a wide range of machine elements operating in the piezoviscous-rigid regime of lubrication. A new numerical method of calculating elastic deformation in contact stresses is developed using a biquadratic polynomial to approximate the pressure distribution on the whole domain analyzed. The deformation of every node is expressed as a linear combination of the nodal pressures whose coefficients can be combined into an influence coefficient matrix. This approach has the advantages of improved numerical accuracy, less computing time and smaller storage size required for influence matrix. The ideal elastohydrodynamic lubrication is extended to real bearing systems in order to gain an understanding of failure mechanisms in machine elements. The improved elastic deformation calculation is successfully incorporated into the EHL numerical scheme. Using this revised numerical technique and the flow factor model developed by Patir and Cheng (1978) the surface roughness effects on the elastohydrodynamic lubrication of point contact is considered. Conditions typical of an EHL contact in the piezoviscous-elastic regime entrained in pure rolling are investigated. Results are compared with the smooth surface solutions. Experiments are conducted to study the transient EHL effects in instrument ball bearings.

Jeng, Y. R.↗
NASA NTRS1986

Computation of the radiation characteristics of a generalized phased array

With the advent of monolithic microwave integrated circuit (MMIC) technology, the phased array has become a key component in the design of advanced antenna systems. Array-fed antennas are used extensively in today's multiple beam satellite antennas. A computer program based on a very efficient numerical technique for calculating the radiated power (Romberg integration), directivity, and radiation pattern of a phased array is described. The formulation developed is very general, and takes into account arbitrary element polarization, E- and H-plane element pattern, element location, and complex element excitation. For comparison purposes sample cases have been presented. Excellent agreement has been obtained for all cases. Also included are a user guide and a copy of the computer program.

Acosta, R. J.↗
NASA NTRS1985

Aerodynamic properties of turbulent combustion fields

Flow fields involving turbulent flames in premixed gases under a variety of conditions are modeled by the use of a numerical technique based on the random vortex method to solve the Navier-Stokes equations and a flame propagation algorithm to trace the motion of the front and implement the Huygens principle, both due to Chorin. A successive over-relaxation hybrid method is applied to solve the Euler equation for flows in an arbitrarily shaped domain. The method of images, conformal transformation, and the integral-equation technique are also used to treat flows in special cases, according to their particular requirements. Salient features of turbulent flame propagation in premixed gases are interpreted by relating them to the aerodynamic properties of the flow field. Included among them is the well-known cellular structure of flames stabilized by bluff bodies, as well as the formation of the characteristic tulip shape of flames propagating in ducts. In its rudimentary form, the mechanism of propagation of a turbulent flame is shown to consist of: (1) rotary motion of eddies at the flame front, (2) self-advancement of the front at an appropriate normal burning speed, and (3) dynamic effects of expansion due to exothermicity of the combustion reaction. An idealized model is used to illustrate these fundamental mechanisms and to investigate basic aerodynamic features of flames in premixed gases. The case of a confined flame stabilized behind a rearward-facing step is given particular care and attention. Solutions are shown to be in satisfactory agreement with experimental results, especially with respect to global properties such as the average velocity profiles and reattachment length.

Hsiao, C. C.↗
NASA NTRS1986

On the solution of integral equations with a generalized cauchy kernel

In this paper a certain class of singular integral equations that may arise from the mixed boundary value problems in nonhomogeneous materials is considered. The distinguishing feature of these equations is that in addition to the Cauchy singularity, the kernels contain terms that are singular only at the end points. In the form of the singular integral equations adopted, the density function is a potential or a displacement and consequently the kernel has strong singularities of the form (t-x) sup-2, x sup n-2 (t+x) sup n, (n or = 2, 0x,tb). The complex function theory is used to determine the fundamental function of the problem for the general case and a simple numerical technique is described to solve the integral equation. Two examples from the theory of elasticity are then considered to show the application of the technique.

Kaya, A. C.↗
NASA NTRS1986

On the solution of integral equations with a generalized cauchy kernal

A certain class of singular integral equations that may arise from the mixed boundary value problems in nonhonogeneous materials is considered. The distinguishing feature of these equations is that in addition to the Cauchy singularity, the kernels contain terms that are singular only at the end points. In the form of the singular integral equations adopted, the density function is a potential or a displacement and consequently the kernal has strong singularities of the form (t-x)(-2), x(n-2) (t+x)(n), (n is = or 2, 0 x, t b). The complex function theory is used to determine the fundamental function of the problem for the general case and a simple numerical technique is described to solve the integral equation. Two examples from the theory of elasticity are then considered to show the application of the technique.

Kaya, A. C.↗
NASA NTRS1985

Measurements of distributed polarized radio sources from spinning spacecraft - Effect of a tilted axial antenna ISEE-3 application and results

An analysis is presented of the system response of a satellite receiver-antenna system to locate a radio source when the satellite is tilted on its axis. The satellite is spin stabilized but experiences a tilt due to either a mechanical misalignment or a shift in the electrical axis caused by parasitic currents in other spacecraft structures. The shorter the antenna, the more significant the effects. Numerical techniques are developed for obtaining the Stokes parameters and the angular parameters of a uniform conical source sensed by a linear antenna in order to derive the average power response of a synthesized dipole to a point on a distributed polarized source. Relative gains are calculated along the antenna at different angles to the source. The techniques are applied to sample ISEE-3 satellite data for Type III solar radio bursts which were sensed by an axial and an equatorial antenna. The two antennas permit localization of the source and quantification of the polarization and angular extent of the source. The resulting high precision in calculations of all three source parameters commends use of the model in analyses of data from the planned ULYSSES mission.

Fainberg, J.↗
NASA NTRS1985

Numerical simulation of a turbulent flame stabilized behind a rearward-facing step

Flow of combustible mixtures in a plane channel past a smooth contraction followed by an abrupt expansion, in a typical dump combustor configuration, is modeled by a two-dimensional numerical technique based on the random vortex method. Both the inert and the reacting case are considered. In the latter, the flame is treated as an interface, self-advancing at a prescribed normal burning speed, while the dynamic effects of expansion due to the exothermicity of combustion are expressed by volumetric source lines delineated by its front. Solutions are shown to be in satisfactory agreement with experimental results, especially with respect to global properties such as the average velocity profiles and the reattachment length. The stochastic turbulent velocity components manifest interesting differences, especially near the walls where three-dimensional effects of turbulence are expected to be of importance.

Hsiao, C. C.↗
NASA NTRS1986

Spline-based distributed system identification with application to large space antennas

A parameter and state estimation technique for distributed models is demonstrated through the solution of a problem generic to large space antenna system identification. Assuming the position of the reflective surface of the maypole (hoop/column) antenna to be approximated by the static two-dimensional, stretched-membrane partial differential equation with variable-stiffness coefficient functions, a spline-based approximation procedure is described that estimates the shape and stiffness functions from data set observations. For given stiffness functions, the Galerkin projection with linear spline-based functions is applied to project the distributed problem onto a finite-dimensional subspace wherein algebraic equations exist for determining a static shape (state) prediction. The stiffness functions are then parameterized by cubic splines and the parameters estimated by an output error technique. Numerical results are presented for data descriptive of a 100-m-diameter maypole antenna.

Banks, H. T.↗
NASA NTRS1986

The origin of the moon and the single-impact hypothesis. I

One of the newer ideas regarding the origin of the moon is concerned with a single-impact hypothesis. It is pointed out that this theory has the advantage of overcoming most of the difficulties with the classical theories. The angular momentum of the earth-moon system can easily be obtained by varying the initial conditions of the impact. A series of three-dimensional numerical simulations of the collision between the earth and an object of about 1/10 its mass is presented. Different impact velocities, impact parameters, and initial internal energies are considered. Attention is given to assumptions, the equation of state, numerical techniques utilizing the momentum equation and the energy conservation equation, tests, and initial conditions and units.

Benz, W.↗
NASA NTRS1987

Approximate Simulation of Turbulence

Computed spectra resemble von Karman spectra of frequencies of interest. Numerical technique yields simulated atmospheric-turbulence spectra that closely approximates von Karman spectra within frequency ranges of interest for aircraft response. Suitable for computation in that they represent stable systems and roll off as f-2 at high frequencies outside range of interest.

Campbell, C. W.↗
NASA NTRS1986

Efficient solutions of two-dimensional incompressible steady viscous flows

A simple, efficient, and robust numerical technique is provided for solving two dimensional incompressible steady viscous flows at moderate to high Reynolds numbers. The proposed approach employs an incremental multigrid method and an extrapolation procedure based on minimum residual concepts to accelerate the convergence rate of a robust block-line-Gauss-Seidel solver for the vorticity-stream function Navier-Stokes equations. Results are presented for the driven cavity flow problem using uniform and nonuniform grids and for the flow past a backward facing step in a channel. For this second problem, mesh refinement and Richardson extrapolation are used to obtain useful benchmark solutions in the full range of Reynolds numbers at which steady laminar flow is established.

Morrison, J. H.↗
NASA NTRS1986

Theoretical kinetic computations in complex reacting systems

Nasa Lewis' studies of complex reacting systems at high temperature are discussed. The changes which occur are the result of many different chemical reactions occurring at the same time. Both an experimental and a theoretical approach are needed to fully understand what happens in these systems. The latter approach is discussed. The differential equations which describe the chemical and thermodynamic changes are given. Their solution by numerical techniques using a detailed chemical mechanism is described. Several different comparisons of computed results with experimental measurements are also given. These include the computation of (1) species concentration profiles in batch and flow reactions, (2) rocket performance in nozzle expansions, and (3) pressure versus time profiles in hydrocarbon ignition processes. The examples illustrate the use of detailed kinetic computations to elucidate a chemical mechanism and to compute practical quantities such as rocket performance, ignition delay times, and ignition lengths in flow processes.

Bittker, David A.↗
NASA NTRS1986

Magnetic field re-arrangement after prominence eruption

It has long been known that magnetic reconnection plays a fundamental role in a variety of solar events. Although mainly invoked in flare problems, large scale loops interconnecting active regions, evolving coronal hole boundaries, the solar magnetic cycle itself, provide different evidence of phenomena which involve magnetic reconnection. A further example might be given by the magnetic field rearrangement which occurs after the eruption of a prominence. Since most often a prominence reforms after its disappearance and may be observed at about the same position it occupied before erupting, the magnetic field has to undergo a temporary disruption of relax back, via reconnection, to a configuration similar to the previous one. The above sequence of events is best observable in the case of two ribbon (2-R) flares but most probably is associated with all filament eruptions. Even if the explanation of the magnetic field rearrangement after 2-R flares in terms of reconnection is generally accepted, the lack of a 3-dimensional model capable of describing the field reconfiguration, has prevented, up to now, a thorough analysis of its topology as traced by H alpha/x ray loops. A numerical technique is presented which enables oneto predict and visualize the reconnected configuration, at any time, and therefore allows one to make a significant comparison of observations and model predictions throughout the whole process.

Kopp, R. A.↗
NASA NTRS1987

Phase space simulation of collisionless stellar systems on the massively parallel processor

A numerical technique for solving the collisionless Boltzmann equation describing the time evolution of a self gravitating fluid in phase space was implemented on the Massively Parallel Processor (MPP). The code performs calculations for a two dimensional phase space grid (with one space and one velocity dimension). Some results from calculations are presented. The execution speed of the code is comparable to the speed of a single processor of a Cray-XMP. Advantages and disadvantages of the MPP architecture for this type of problem are discussed. The nearest neighbor connectivity of the MPP array does not pose a significant obstacle. Future MPP-like machines should have much more local memory and easier access to staging memory and disks in order to be effective for this type of problem.

White, Richard L.↗
NASA NTRS1986

An introduction to three-dimensional climate modeling

The development and use of three-dimensional computer models of the earth's climate are discussed. The processes and interactions of the atmosphere, oceans, and sea ice are examined. The basic theory of climate simulation which includes the fundamental equations, models, and numerical techniques for simulating the atmosphere, oceans, and sea ice is described. Simulated wind, temperature, precipitation, ocean current, and sea ice distribution data are presented and compared to observational data. The responses of the climate to various environmental changes, such as variations in solar output or increases in atmospheric carbon dioxide, are modeled. Future developments in climate modeling are considered. Information is also provided on the derivation of the energy equation, the finite difference barotropic forecast model, the spectral transform technique, and the finite difference shallow water waved equation model.

Washington, W. M.↗
NASA NTRS1987

Construction of explicit and implicit symmetric TVD schemes and their applications

A one-parameter family of explicit and implicit total variation diminishing (TVD) schemes is developed which permits incorporation of an expanded group of slope and flux limiters. The numerical technique is intended for use in calculations which include a time-differencing scheme and an optional Lax-Wendroff scheme. Methods of extending the TVD models to nonlinear scalar equations and systems of hyperbolic conservation equations are described. Sample results are presented from calculations of shocked flows around NACA 0012 and NACA 0018 airfoils.

Yee, H. C.↗
NASA NTRS1987

Cutoff wavenumbers and modes for annular-cross-section waveguide with eccentric inner conductor of small radius

Analytical expressions are derived for the cutoff wavenumbers and the corresponding modes in annular-cross-section waveguides having inner conductors of small radius. Waveguides with circular and rectangular outer boundary are considered. In the case of the circular eccentric annular waveguide, comparison is made between the values of cutoff wavenumbers computed from the expressions derived in this paper and data obtained by a more rigorous numerical technique.

Davidovitz, Marat↗