Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical integration”

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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Diffusion Of Mass In Evaporating Multicomponent Drops

Report summarizes study of diffusion of mass and related phenomena occurring in evaporation of dense and dilute clusters of drops of multicomponent liquids intended to represent fuels as oil, kerosene, and gasoline. Cluster represented by simplified mathematical model, including global conservation equations for entire cluster and conditions on boundary between cluster and ambient gas. Differential equations of model integrated numerically. One of series of reports by same authors discussing evaporation and combustion of sprayed liquid fuels.

Bellan, Josette↗
NASA NTRS1992

Developing A Navier-Stokes Algorithm For Supercomputers

Report discusses development of algorithm for solution of Navier-Stokes equations of flow on parallel-processing supercomputers. Involves combination of prior techniques to form algorithm to compute flows in complicated three-dimensional configurations. Includes explicit finite-difference numerical-integration scheme applicable to flows represented by hierarchy of mathematical models ranging from Euler to full Navier-Stokes. Of interest to researchers looking for ways to structure problems for greater computational efficiency.

Swisshelm, Julie M.↗
NASA NTRS1992

Examples Of Synthesis Of Dual-Shaped Reflectors

Report presents examples to demonstrate validity and utility of method of synthesis of offset dual-shape reflectors. Method of synthesis described by the authors in previous journal article. Current report reviews derivation of partial differential equations and iterative method of numerical solution. Discusses significance of starting point of numerical integration on each reflector surface; this point could be at center, on outer rim, or at interior point. Emphasizes that one of notable attributes of partial differential equations is speed with which they can be solved.

Galindo, Victor↗
NASA NTRS1991

Asteroid flux and impact cratering rate on Venus

By the end of 1990, 65 Venus-crossing asteroids were recognized; these represent 59 percent of the known Earth-crossing asteroids. Further studies, chiefly numerical integrations of orbit evolution, may reveal one or two more Venus crossers among the set of discovered asteroids. A Venus crosser was defined as an asteroid whose orbit can intersect the orbit of Venus as a result of secular (long range) perturbations. Venus crossers revolving on orbits that currently overlap the orbit of Venus are called Venapol asteroids, and those on orbit that don't overlap are called Venamor asteroids; 42 Venapols and 23 Venamors were recognized. Collision probabilities with Venus for 60 of the known Venus crossers were determined.

Shoemaker, E. M.↗
NASA NTRS1991

Mapping the stability field of Jupiter Trojans

Jupiter Trojans are a remnant of outer solar system planetesimals captured into stable or quasistable libration about the 1:1 resonance with the mean motion of Jupiter. The observed swarms of Trojans may provide insight into the original mass of condensed solids in the zone from which the Jovian planets accumulated, provided that the mechanisms of capture can be understood. As the first step toward this understanding, the stability field of Trojans were mapped in the coordinate proper eccentricity, e(sub p), and libration amplitude, D. To accomplish this mapping, the orbits of 100 particles with e(sub p) in the range of 0 to 0.8 and D in the range 0 to 140 deg were numerically integrated. Orbits of the Sun, the four Jovian planets, and the massless particles were integrated as a full N-body system, in a barycentric frame using fourth order symplectic scheme.

Levison, H. F.↗
NASA NTRS1991

Long-term evolution of a planetesimal swarm in the vicinity of a protoplanet

Many models of planet formation involve scenarios in which one or a few large protoplanets interact with a swarm of much smaller planetesimals. In such scenarios, three-body perturbations by the protoplanet as well as mutual collisions and gravitational interactions between the swarm bodies are important in determining the velocity distribution of the swarm. We are developing a model to examine the effects of these processes on the evolution of a planetesimal swarm. The model consists of a combination of numerical integrations of the gravitational influence of one (or a few) massive protoplanets on swarm bodies together with a statistical treatment of the interactions between the planetesimals. Integrating the planetesimal orbits allows us to take into account effects that are difficult to model analytically or statistically, such as three-body collision cross-sections and resonant perturbations by the protoplanet, while using a statistical treatment for the particle-particle interactions allows us to use a large enough sample to obtain meaningful results.

Kary, David M.↗
NASA NTRS1991

Collision probabilities in the presence of nebular gas drag

We are developing a model to determine what fraction of the planetesimals would have hit a protoplanet on their sunward journey as opposed to having a close approach and passing into an inferior orbit. The model involves direct numerical integration of restricted-three-body orbits using a predictor-corrector integrator. A simple gas drag law with a v(exp 2) dependence was also included in the equations of motion. Runs of 100 to 500 particles were already performed, while some future runs may require several times this number in order to get good impact statistics. All planetesimals start in superior orbits with semi-major axes 5 to 10 R(sub H) from the protoplanets, where R(sub H) is the protoplanet's Hill Sphere radius. The orbit is followed until the planetesimal passed into an inferior orbit at least 10 R(sub H) from the protoplanet. This process typically requires 10(exp 4) to 10(exp 5) orbits.

Kary, David M.↗
NASA NTRS1991

Prediction of noise field of a propfan at angle of attack

A method for predicting the noise field of a propfan operating at an angle of attack to the oncoming flow is presented. The method takes advantage of the high-blade-count of the advanced propeller designs to provide an accurate and efficient formula for predicting their noise field. The formula, which is written in terms of the Airy function and its derivative, provides a very attractive alternative to the use of numerical integration. A preliminary comparison shows rather favorable agreement between the predictions from the present method and the experimental data.

Envia, Edmane↗
NASA NTRS1991

Evaluation of Rotordynamic Coefficients of Look-through Labyrinths by Means of a Three Volume Bulk Flow Model

To describe the compressible, turbulent flow in a labyrinth seal, a three volume bulk flow model is presented. The conservation equations for mass, momentum, and energy are established in every control volume. A perturbation analysis is performed, yielding zeroth order equations for centric rotor position and first order equations describing the flow field for small rotor motions around the seal center. The equations are integrated numerically. From perturbation pressure, the forces on the shaft and the dynamic coefficients are calculated.

Nordmann, R.↗
NASA NTRS1992

High-frequency techniques for RCS prediction of plate geometries

The principal-plane scattering from perfectly conducting and coated strips and rectangular plates is examined. Previous reports have detailed Geometrical Theory of Diffraction/Uniform Theory of Diffraction (GTD/UTD) solutions for these geometries. The GTD/UTD solution for the perfectly conducting plate yields monostatic radar cross section (RCS) results that are nearly identical to measurements and results obtained using the Moment Method (MM) and the Extended Physical Theory of Diffraction (EPTD). This was demonstrated in previous reports. The previous analysis is extended to bistatic cases. GTD/UTD results for the principal-plane scattering from a perfectly conducting, infinite strip are compared to MM and EPTD data. A comprehensive overview of the advantages and disadvantages of the GTD/UTD and of the EPTD and a detailed analysis of the results from both methods are provided. Several previous reports also presented preliminary discussions and results for a GTD/UTD model of the RCS of a coated, rectangular plate. Several approximations for accounting for the finite coating thickness, plane-wave incidence, and far-field observation were discussed. Here, these approximations are replaced by a revised wedge diffraction coefficient that implicitly accounts for a coating on a perfect conductor, plane-wave incidence, and far-field observation. This coefficient is computationally more efficient than the previous diffraction coefficient because the number of Maliuzhinets functions that must be calculated using numerical integration is reduced by a factor of 2. The derivation and the revised coefficient are presented in detail for the hard polarization case. Computations and experimental data are also included. The soft polarization case is currently under investigation.

Balanis, Constantine A.↗
NASA NTRS1992

Numerical investigation of the effects of icing on fixed and rotary wing aircraft

A 2-D multi-element airfoil code was modified to study the effects of icing on the aerodynamic characteristics of high lift systems. In each zone of the flow field, the solver numerically integrates the 2-D compressible Navier-Stokes equations using a time marching scheme. The surface pressure distribution is generated over a GAW 130 airfoil/flap combination for a flap setting of 25 degrees, and an angle of attack equal to 5 degrees, at a freestream Mach number equal to 0.3. A series of calculations were performed to determine the effects of small scale ice build up on the high lift characteristics of this arifoil/flap combination. The appendix summarizes this progress. Joint studies on correlation of a 3-D iced wing code with experimental data reviewed new measured laser Doppler velocimeter data in the separated region behind the leading edge ice shape. A version of the iced wing analysis using the Roe scheme was developed to evaluate the poor correlation between the computed and measured velocities in the separated region. Work on the extension of the wing-alone analysis to wing body configuration began with modifications to the 3-D iced wing analyses to accept externally generated grids and multi-block grids.

Sankar, L. N.↗
NASA NTRS1991

Transient cooling of a square region of radiating medium

A time-accurate numerical solution was carried out for transient radiative cooling of a gray emitting and absorbing medium in a square two-dimensional region. The integro-differential energy equation for transient temperature distributions was solved in two stages. At each time increment, the local radiative source term was obtained by numerical integration of the temperature field using two-dimensional Gaussian integration over rectangular subregions. Then the differential portion of the equation was integrated forward in time by use of the local first and second time derivatives. The results were compared with available limiting case, and excellent agreement was obtained. Transient results are given for a wide range of optical thicknesses of the region. Optimum transient cooling is obtained when the optical side length is about 4.

Siegel, Robert↗
NASA NTRS1991

Two-dimensional hydrodynamic calculations of the nonlinear development of the Goldreich-Schubert-Fricke instability in a rotating annulus

Two-dimensional nonlinear hydrodynamic calculations are presented which may help assess the effectiveness of the instability in transporting angular momentum in the equatorial zones of stars and planets which are stably stratified with respect to convection. The calculations were made by numerically integrating the 2D axisymmetric Navier-Stokes equations, including viscosity and heat conduction. The instability was followed into the nonlinear regime. The maximum rms velocity amplitude was found to correlate well with the product of the linear growth rate and radial length scale of the instability, consistent with the idea that the instability grows to an amplitude such that an eddy turnover time becomes equal to the growth time defined by the inverse of the growth rate. The time scale for angular momentum to be redistributed to a state of marginal stability was consistent with this picture. The results suggest that in physical situations a state of marginal stability will be maintained, since departures from such a state will be rapidly corrected.

Korycansky, D. G.↗
NASA NTRS1992

Tidal evolution of globular clusters. I - Method

Tidal evolution of globular clusters is regulated by both Galactic tidal effects and internal relaxation processes. In order to investigate the tidal evolution of globular clusters, a numerical scheme which utilizes a Fokker-Planck approach as well as direct numerical integration of the restricted three-body problem is developed. In the inner regions of the cluster, stellar orbits are mapped with the cluster's gravitational potential and orbit-averaged diffusion coefficients. In the outer regions, the Galactic tidal field is explicitly included in the direct orbital integration. This method is presented here with some tests on King-Michie models.

Oh, K. S.↗
NASA NTRS1992

Convective instability in differentially rotating disks

A normal mode analysis for nonaxisymmetric perturbations in a thin, differentially rotating disk with a vertical structure that is isothermal and convectively unstable is performed. The vertical gravity is assumed to be external and constant. The perturbation scale is assumed to be much shorter than the radius of the disk but comparable to or less than the thickness. The initial value problem is formulated in shearing coordinates. Dispersion relations are obtained for the three limiting cases of zero shear, axisymmetric perturbations, and small radial wavelengths. The full effects of shear are studied by integrating numerically the initial value problem. Nonaxisymmetric local Fourier modes are found to have a radial wavenumber that increases linearly with time in proportion to the shear times the azimuthal wavenumber. While Coriolis forces exert stabilizing effects on the convective modes, reducing their growth rate and the range of unstable wavelengths, shear has destablizing effects inasmuch as it reduces the epicyclic frequency at a given angular velocity. In a Keplerian disk, perturbations with azimuthal wavelengths about 2 times smaller than vertical wavelengths grow exponentially.

Ryu, Dongsu↗
NASA NTRS1992

Discrete shaped strain sensors for intelligent structures

Design of discrete, highly distributed sensor systems for intelligent structures has been studied. Data obtained indicate that discrete strain-averaging sensors satisfy the functional requirements for distributed sensing of intelligent structures. Bartlett and Gauss-Hanning sensors, in particular, provide good wavenumber characteristics while meeting the functional requirements. They are characterized by good rolloff rates and positive Fourier transforms for all wavenumbers. For the numerical integration schemes, Simpson's rule is considered to be very simple to implement and consistently provides accurate results for five sensors or more. It is shown that a sensor system that satisfies the functional requirements can be applied to a structure that supports mode shapes with purely sinusoidal curvature.

Andersson, Mark S.↗
NASA NTRS1992

Structural tailoring/analysis for hypersonic components - Executive system development

No direct analytical or integrated numerical tool exists today for the optimal design of a generic class of built-up actively cooled composite structure for applications in hypersonic propulsion ducts. The need exists for a numerical tool to perform the comprehensive design/analysis of a panel on the inlet wall under hypersonic flight conditions. Such a tool requires relatively complex multi-disciplinary analysis. One such numerical tool controlled by an executive system has been developed and is named as STAHYC (Structural Tailoring/Analysis for HYpersonic Components). A detailed account of the executive system development of STAHYC along with the results of one example inlet panel design problem is given in this paper.

Narayanan, G. V.↗