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.

SourceResearch connectionFrom yearThrough year

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

At least 1,009 records · Page 56

All sources NASA NTRS DOE OSTI

NASA NTRS1992

Calculation of two-dimensional shear-driven cavity flows at high Reynolds numbers

The time-dependent Navier-Stokes equations are numerically integrated for two-dimensional incompressible viscous flow in a shear-driven square cavity. Using a time-splitting method and finite differences on a staggered mesh, the momentum and pressure equations are directly solved by a tensor product method where one finite difference direction is diagonalized by eigenvalue decomposition. The effects of increasing Reynolds number are studied and the developing boundary layer is captured by using a finely clustered mesh. At Re = 30,000 the flow is in a continuously developing unsteady regime. Power spectrum plots indicate that the unsteady flow oscillates with one fundamental frequency and exhibits some characteristics of transition between laminar and turbulent states.

Huser, A.↗

NASA NTRS1992

Chaotic scattering of pitch angles in the current sheet of the magnetotail

The modified Harris field model is used to investigate the process of pitch angle scattering by a current sheet. The relationship between the incoming asymptotic pitch angle alpha(in) and the outgoing asymptotic pitch angle alpha(out) is studied from first principles by numerically integrating the equation of motion. Evidence that charged particles undergo chaotic scattering by the current sheet is found. For fixed alpha(in), it is shown that alpha(out) exhibits sensitive dependence on the energy parameter in certain energy ranges. For a fixed energy parameter value in the same energy ranges, alpha(out) sensitively depends on alpha(in). For other energy values, alpha(out) does not show sensitive dependence on alpha(in) for most phase angles. A distribution of alpha(in) is mapped from the asymptotic region to the midplane, and it is found that the resulting particle distribution should have beam structures with well-collimated pitch angles near each resonance energy value. Implications for the particle distribution functions in the earth's magnetotail are discussed.

Burkhart, G. R.↗

NASA NTRS1992

Libration point staging options for SEI lunar missions - Station keeping implications

Libration point staging is one alternative to lunar orbit rendezvous which is being examined for the Space Exploration Initiative lunar missions. The libration point staging strategy being considered would eliminate the lunar mid-latitude accessibility constraint for long duration missions. Staging options for such a strategy are discussed with particular emphasis on the impact of the libration points' stability. The equations of motion for a satellite at the libration points are numerically integrated to investigate the stability and to determine the station-keeping requirements. These requirements are presented, and overall requirements for libration point staging briefly summarized. The results indicate that the slight eccentricity of the Moon's orbit would result in excessive station-keeping requirements for long duration missions.

Haynes, Davy A.↗

NASA NTRS1992

An analytical evaluation of mascon effects in a semi-analytic theory

A semianalytic method to predict the orbit of a spacecraft orbiting a planet and perturbed by oblateness and a gravity anomaly (mascon) is presented. The Hamiltonian is first expressed in closed form in terms of the nonsingular canonical set of Poincare elements where the mascon is represented by a point-mass. The planet global potential is limited to the second and third harmonics. The mascon potential in harmonics is limited to second degree and first order. The original Lie-Hori method is then applied to transform the canonical equations of motion from the osculating element space to the mean element space. Numerical integration of the new equations is used.

Salama, Ahmed H.↗

NASA NTRS1992

Boundary Heat Fluxes for Spectral Radiation from a Uniform Temperature Rectangular Medium

The effect of spectral behavior is analytically shown for radiation in a 2D rectangular geometry. The solution provides exact boundary heat flux values that can be used for comparison with values obtained from general computer programs. The spectral solution presented can be easily evaluated by numerical integration for complex variations of the spectral absorption coefficient with wavelength.

Siegel, Robert↗

NASA NTRS1991

Hypersonic merged layer blunt body flows with wakes

An attempt is made here to understand the basic physics of the flowfield with wake on a blunt body of revolution under hypersonic rarefied conditions. A merged layer model of flow is envisioned. Full steady-state Navier-Stokes equations in spherical polar coordinate system are computed from the surface with slip and temperature jump conditions to the free stream by the Accelerated Successive Replacement method of numerical integration. Analysis is developed for bodies of arbitrary shape, but actual computations have been carried out for a sphere and sphere-cone body. Particular attention is paid to set the limit of the onset of separation, wake closure, shear-layer impingement, formation and dissipation of the shocks in the flowfield. Validity of the results is established by comparing the present results for sphere with the corresponding results of the SOFIA code in the common region of their validity and with the experimental data.

Jain, Amolak C.↗

NASA NTRS1992

The Space Transportation System summer environment on launch pad

This paper describes a 2D flow and thermalanalysis to determine the solar effect on the Space Shuttle launch components subsequent to the external tank (ET) loading operation in extremely hot conditions. An existing CFD code Parabolic Hyperbolic or Elliptical Numerical Integration Code Series was used in the study. The analysis was done for a 2D slice between planes perpendicular to the longitudinal axis of the STS and passing through the lower portions of the Redesigned Solid Rocket Motors (RSRMs), the ET, and the wing of the Orbiter. The results are presented as local and average values of the heat transfer coefficient, and the Nusselt number, and the surface temperature around the RSRMs and the ET. Solar heating effects increased the surface temperatures of the RSRMs by 9-11 F. Higher prelaunch surface temperatures measured on the east and west RSRMs (in the inboard region between the RSRMs and the ET) during 19 most recent launches of the STS are correlated as a function of the ambient temperature.

Ahmad, R. A.↗

NASA NTRS1993

Development Of A Navier-Stokes Computer Code

Report discusses aspects of development of CENS3D computer code, solving three-dimensional Navier-Stokes equations of compressible, viscous, unsteady flow. Implements implicit finite-difference or finite-volume numerical-integration scheme, called "lower-upper symmetric-Gauss-Seidel" (LU-SGS), offering potential for very low computer time per iteration and for fast convergence.

Yoon, Seokkwan↗

NASA NTRS1993

General-Purpose Time-Accurate Adaptive-Grid Method

Adaptive-grid method provides for generation of time-accurate grids used to solve wide variety of problems numerically. Intended use in computation of compressible flows. Part of solution method including algebraic grid-movement scheme; algorithm solving, by use of explicit numerical-integration scheme, Euler equations of unsteady flow; and temporal coupling routine utilizing grid-prediction/correction method to link dynamic grid with Euler-equation-solving algorithm. Improves ability to adapt grid to multiple features of solution and involves use of coarse grid in prediction stage of temporal coupling reducing central-processing-unit cost of obtaining new grid.

Bockelie, Michael J.↗

NASA NTRS1993

LSENS, a general chemical kinetics and sensitivity analysis code for gas-phase reactions: User's guide

A general chemical kinetics and sensitivity analysis code for complex, homogeneous, gas-phase reactions is described. The main features of the code, LSENS, are its flexibility, efficiency and convenience in treating many different chemical reaction models. The models include static system, steady, one-dimensional, inviscid flow, shock initiated reaction, and a perfectly stirred reactor. In addition, equilibrium computations can be performed for several assigned states. An implicit numerical integration method, which works efficiently for the extremes of very fast and very slow reaction, is used for solving the 'stiff' differential equation systems that arise in chemical kinetics. For static reactions, sensitivity coefficients of all dependent variables and their temporal derivatives with respect to the initial values of dependent variables and/or the rate coefficient parameters can be computed. This paper presents descriptions of the code and its usage, and includes several illustrative example problems.

Radhakrishnan, Krishnan↗

NASA NTRS1992

Antenna analysis using neural networks

Conventional computing schemes have long been used to analyze problems in electromagnetics (EM). The vast majority of EM applications require computationally intensive algorithms involving numerical integration and solutions to large systems of equations. The feasibility of using neural network computing algorithms for antenna analysis is investigated. The ultimate goal is to use a trained neural network algorithm to reduce the computational demands of existing reflector surface error compensation techniques. Neural networks are computational algorithms based on neurobiological systems. Neural nets consist of massively parallel interconnected nonlinear computational elements. They are often employed in pattern recognition and image processing problems. Recently, neural network analysis has been applied in the electromagnetics area for the design of frequency selective surfaces and beam forming networks. The backpropagation training algorithm was employed to simulate classical antenna array synthesis techniques. The Woodward-Lawson (W-L) and Dolph-Chebyshev (D-C) array pattern synthesis techniques were used to train the neural network. The inputs to the network were samples of the desired synthesis pattern. The outputs are the array element excitations required to synthesize the desired pattern. Once trained, the network is used to simulate the W-L or D-C techniques. Various sector patterns and cosecant-type patterns (27 total) generated using W-L synthesis were used to train the network. Desired pattern samples were then fed to the neural network. The outputs of the network were the simulated W-L excitations. A 20 element linear array was used. There were 41 input pattern samples with 40 output excitations (20 real parts, 20 imaginary). A comparison between the simulated and actual W-L techniques is shown for a triangular-shaped pattern. Dolph-Chebyshev is a different class of synthesis technique in that D-C is used for side lobe control as opposed to pattern shaping. The interesting thing about D-C synthesis is that the side lobes have the same amplitude. Five-element arrays were used. Again, 41 pattern samples were used for the input. Nine actual D-C patterns ranging from -10 dB to -30 dB side lobe levels were used to train the network. A comparison between simulated and actual D-C techniques for a pattern with -22 dB side lobe level is shown. The goal for this research was to evaluate the performance of neural network computing with antennas. Future applications will employ the backpropagation training algorithm to drastically reduce the computational complexity involved in performing EM compensation for surface errors in large space reflector antennas.

Smith, William T.↗

NASA NTRS1992

Stiffnesses by (TtN) ensemble molecular dynamics

Calculation of elastic constants requires initializing a lattice of atoms, an h matrix defining the computational cell boundaries, and the parameters defining the interatomic potentials, then numerically integrating the equations of motion for the atoms and the computational cell boundaries, calculating strains and stresses at each timestep, and updating at each timestep the average quantities required for calculation of fluctuations. Enforcements of the DYNAMO FORTRAN code included addition of documentation, improved memory management and data flow, and enabling use of several interatomic potentials. Documentation included headers for each subroutine, which document the parameters input to the routine, algorithms employed, and output.

Tibbits, Patrick↗

NASA NTRS1992

Parabolized Navier-Stokes methods for hypersonic flows

The development and numerical solution of the Parabolized Navier-Stokes (PNS) equations is discussed. Special emphasis is placed on numerical integration algorithms for PNS solvers, most notably recent applications of upwind schemes. Application of this technology to high speed, high Reynolds number flows was the focus, and modeling techniques for real gas effects and turbulence were discussed. Experiences in validating the UPS upwind PNS flow solver for relatively simple body shapes are presented. The test cases range in complexity from two dimensional perfect gas cases to three dimensional finite rate chemistry cases.

Lawrence, S. L.↗

NASA NTRS1992

Precessional quantities for the Earth over 10 Myr

The insolation parameters of the Earth depend on its orbital parameters and on the precession and obliquity. Until 1988, the usually adopted solution for paleoclimate computation consisted in (Bretagnon, 1974) for the orbital elements of the Earth, which was completed by (Berger, 1976) for the computation of the precession and obliquity of the Earth. In 1988, I issued a solution for the orbital elements of the Earth, which was obtained in a new manner, gathering huge analytical computations and numerical integration (Laskar, 1988). In this solution, which will be denoted La88, the precession and obliquity quantities necessary for paleoclimate computations were integrated at the same time, which insure good consistency of the solutions. Unfortunately, due to various factors, this latter solution for the precession and obliquity was not widely distributed (Berger, Loutre, Laskar, 1988). On the other side, the orbital part of the solution La88 for the Earth, was used in (Berger and Loutre, 1991) to derive another solution for precession and obliquity, aimed to climate computations. I also issued a new solution (La90) which presents some slight improvements with respect to the previous one (Laskar, 1990). As previously, this solution contains orbital, precessional, and obliquity variables. The main features of this new solution are discussed.

Laskar, Jacques↗

NASA NTRS1992

The orbital evolution of real asteroids near the 4:1 mean-motion resonance with Jupiter

Numerical integrations of the orbits of ten asteroids with osculating elements near the 4:1 mean-motion resonance with Jupiter have been performed over 200,000 years into the future. A variety of orbital evolutions was found, depending on the start values of the semi-major axis. The orbit of asteroid 1983 RJ4, which lies almost exactly at the resonance center, experiences large variations in eccentricity, evolving into an Earth-crosser on a time-scale of a few 10(exp 4) years. This makes this region a potential source for Apollo objects and meteoritic material, although the width of the resonance region in semi-major axis seems to be very narrow.

Dahlgren, Mats↗

NASA NTRS1992

The effect of secular resonances in the asteroid region between 2.1 and 2.4 AU

The asteroid region between 2.1 and 2.4 AU appears to be depopulated at inclinations i greater than 12 deg. This region is surrounded by the three main secular resonances nu(sub 5), nu(sub 6), and nu(sub 16) and is crossed by higher order secular resonances. Secular resonances appear to overlap in this region. Numerical integrations of the orbits of seventeen fictituous asteroids with initial inclinations 12 deg less than or equal to i less than or equal to 20 deg show the following: (1) this particular asteroid region is not depopulated in our computer experiment on timescales of 2.7 Myrs; (2) inclinations are pumped up by successive crossings through higher order secular resonances while eccentricities are not increased sufficiently to produce planet-crossers; (3) bodies located in the bordering nu(sub 6) resonance with semi-major axes a less than or equal to 2.4 AU become Earth-crossers on a time scale of 1 Myr; and (4) we confirm the result that modes due to higher order secular resonances must be eliminated when proper elements are computed.

Froeschle, Christiane↗

NASA NTRS1992

Asteroid-type orbit evolution near the 5:2 resonance

In this case of the 5:2 commensurability with the motion of Jupiter, an asteroid can reach the orbits of Mars, Earth, and Venus when eccentricity e is greater than 0.41, 0.65, and 0.74, respectively. For individual fictitious asteroids, Ipatov and Yoshikawa obtained a growth in e from 0.15 to 074-0.76. Rates of changes in orbital orientations are different for Mars, Earth, Venus, and the asteroid. Therefore, for corresponding values of e, the asteroid could encounter these planets and leave the gap at those encounters. In order to investigate this hypothesis of the 5:2 Kirkwood gap formation, Ipatov studied the regions of initial data for which the eccentricities of asteroids located near the 5:2 commensurability exceeded 0.41 during evolution. The orbit evolution for 500 fictitious asteroids was investigated by numerical integration of the complete (unaveraged) equations of motion for the three-body problem (Sun-Jupiter-asteroid). The equations of motion were integrated in the time intervals T is greater than or equal to 5(10)(exp 3)t(sub J) (t(sub J) is the heliocentric orbital period of Jupiter) in the planar model, T is greater than or equal to 10(exp 4)t(sub J) at initial inclination 5 deg is less than or equal to i(sub 0) is less than or equal to 20 deg and T = 10(exp 5)t(sub J) at i(sub 0) = 40 deg. The larger interval T was taken at i(sub 0) = 40 deg because in this case for the majority of runs maximum values of e and i were reached in the time delta(t) is greater than 2(10)(exp 4)t(sub J).

Ipatov, S. I.↗

NASA NTRS1992

Dynamical timescales in the Jupiter family

Numerically integrated fictitious comets starting in orbits perihelion tangent to Jupiter have been used to estimate the duration of a typical visit to the observable Jupiter family of comets. The results show values of 3 to 6 thousand years, narrowing the previously estimated interval of 10(exp 3) to 10(exp 4) years.

Lindgren, Mats↗