Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time integrator”

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 685 records · Page 38

Goldstone radar contributions to Mars Pathfinder landing safety

Goldstar radar can provide topography 'profiles', statistical surface roughness, and radar images within a few degrees of the sub-Earth point. Goldstone/Very Large Array (VLA) bistatic radar observations can image the whole disk of Mars with integration times on the order of ten min before pixel smearing occurs. Data from all these radar techniques can be useful for observing the local surface conditions relating to landing safety issues for Mars Pathfinder.

Slade, Martin A.↗

A Lagrangian dynamic subgrid-scale model turbulence

A new formulation of the dynamic subgrid-scale model is tested in which the error associated with the Germano identity is minimized over flow pathlines rather than over directions of statistical homogeneity. This procedure allows the application of the dynamic model with averaging to flows in complex geometries that do not possess homogeneous directions. The characteristic Lagrangian time scale over which the averaging is performed is chosen such that the model is purely dissipative, guaranteeing numerical stability when coupled with the Smagorinsky model. The formulation is tested successfully in forced and decaying isotropic turbulence and in fully developed and transitional channel flow. In homogeneous flows, the results are similar to those of the volume-averaged dynamic model, while in channel flow, the predictions are superior to those of the plane-averaged dynamic model. The relationship between the averaged terms in the model and vortical structures (worms) that appear in the LES is investigated. Computational overhead is kept small (about 10 percent above the CPU requirements of the volume or plane-averaged dynamic model) by using an approximate scheme to advance the Lagrangian tracking through first-order Euler time integration and linear interpolation in space.

Meneveau, C.↗

Navier-Stokes simulations of the effects of suction holes on a flat plate boundary layer

Direct numerical simulations (DNS) of the Navier-Stokes equations are employed to explore the effects of suction holes on transition in a laminar flat plate boundary layer. The Navier-Stokes equations are cast in vorticity-velocity form. Periodicity is imposed in spanwise direction; all other spatial derivatives are discretized with fourth order compact differences. An explicit fourth order Runge-Kutta scheme is employed for the time-integration of the vorticity transport equations. Suction is applied through a row of holes aligned in spanwise direction. For low suction strengths, each hole generates a pair of stable streamwise vortices. When the suction strength exceeds a critical value, the vortices become unstable. For high suction strengths, vortex shedding occurs right at the suction holes. Our numerical findings agree well with experimental observations.

Meitz, Hubert L.↗

User's Guide for ECAP2D: an Euler Unsteady Aerodynamic and Aeroelastic Analysis Program for Two Dimensional Oscillating Cascades, Version 1.0

This guide describes the input data required for using ECAP2D (Euler Cascade Aeroelastic Program-Two Dimensional). ECAP2D can be used for steady or unsteady aerodynamic and aeroelastic analysis of two dimensional cascades. Euler equations are used to obtain aerodynamic forces. The structural dynamic equations are written for a rigid typical section undergoing pitching (torsion) and plunging (bending) motion. The solution methods include harmonic oscillation method, influence coefficient method, pulse response method, and time integration method. For harmonic oscillation method, example inputs and outputs are provided for pitching motion and plunging motion. For the rest of the methods, input and output for pitching motion only are given.

Reddy, T. S. R.↗

ICASE/LaRC Workshop on Benchmark Problems in Computational Aeroacoustics (CAA)

The proceedings of the Benchmark Problems in Computational Aeroacoustics Workshop held at NASA Langley Research Center are the subject of this report. The purpose of the Workshop was to assess the utility of a number of numerical schemes in the context of the unusual requirements of aeroacoustical calculations. The schemes were assessed from the viewpoint of dispersion and dissipation -- issues important to long time integration and long distance propagation in aeroacoustics. Also investigated were the effect of implementation of different boundary conditions. The Workshop included a forum in which practical engineering problems related to computational aeroacoustics were discussed. This discussion took the form of a dialogue between an industrial panel and the workshop participants and was an effort to suggest the direction of evolution of this field in the context of current engineering needs.

Hardin, Jay C.↗

Application of traditional CFD methods to nonlinear computational aeroacoustics problems

This paper describes an implementation of a high order finite difference technique and its application to the category 2 problems of the ICASE/LaRC Workshop on Computational Aeroacoustics (CAA). Essentially, a popular Computational Fluid Dynamics (CFD) approach (central differencing, Runge-Kutta time integration and artificial dissipation) is modified to handle aeroacoustic problems. The changes include increasing the order of the spatial differencing to sixth order and modifying the artificial dissipation so that it does not significantly contaminate the wave solution. All of the results were obtained from the CM5 located at the Numerical Aerodynamic Simulation Laboratory. lt was coded in CMFortran (very similar to HPF), using programming techniques developed for communication intensive large stencils, and ran very efficiently.

Chyczewski, Thomas S.↗

The use of staggered scheme and an absorbing buffer zone for computational aeroacoustics

Various problems from those proposed for the Computational Aeroacoustics (CAA) workshop were studied using second and fourth order staggered spatial discretizations in conjunction with fourth order Runge-Kutta time integration. In addition, an absorbing buffer zone was used at the outflow boundaries. Promising results were obtained and provide a basis for application of these techniques to a wider variety of problems.

Nark, Douglas M.↗

Computational aeroacoustics on massively parallel computers

A computational fluid dynamics (CFD) technique is employed to solve aeroacoustics problems on massively parallel computers. The algorithm is based on a 4th order accurate central finite difference and 4th order accurate 4-stage Runge-Kutta time integration method that solves the 3-D full Navier-Stokes/Euler equations on a curvilinear coordinate system. The code has been developed to predict noise radiation from ducted fans. A specific application of this code is made here to the oscillating circular piston problem (Workshop Category 4, Problem 2) by solving the full Euler equations with nonreflecting boundary conditions in the axisymmetric mode. The features of the code that apply to this particular problem are described in this paper.

Oezyoeruek, Yusuf↗

Comparison of spatial numerical operators for duct-nozzle acoustics

A production Navier-Stokes/Euler CFD code, NASTD, developed for aircraft flowfield analysis has been modified to analyze acoustic fields associated with propulsion exhaust systems. The modified code has been applied to the Category 5 nozzle problem using six different spatial discretization schemes combined with a third-order, compact storage Runge-Kutta time integration. NASTD was found capable of tracking pressure disturbances normalized by the freestream value of order 10(exp -6), even with lower-order schemes, for the benchmark problem.

Cain, A. B.↗

A computational procedure for multibody systems including flexible beam dynamics

A computational procedure suitable for the solution of equations of motions for flexible multibody systems has been developed. The flexible beams are modeled using a fully nonlinear theory which accounts for both finite rotations and large deformations. The present formulation incorporates physical measures of conjugate Cauchy stress and covariant strain increments. As a consequence, the beam model can easily be interfaced with real-time strain measurements and feedback control systems. A distinct feature of the present work is the computational preservation of total energy for undamped systems; this is obtained via an objective strain increment/stress update procedure combined with an energy-conserving time integration algorithm which contains an accurate update of angular orientations. The procedure is demonstrated via several example problems.

Downer, J. D.↗

Initialization of soil-water content in regional-scale atmospheric prediction models

The purpose of this study is to demonstrate the feasibility of determining the soil-water content fields required as initial conditions for land surface components within atmospheric prediction models. This is done using a model of the hydrologic balance and conventional meteorological observations, land cover, and soils information. A discussion is presented of the subgrid-scale effects, the integration time, and the choice of vegetation type on the soil-water content patterns. Finally, comparisons are made between two The Pennsylvania State University/National Center for Atmospheric Research mesoscale model simulations, one using climatological fields and the other one using the soil-moisture fields produced by this new method.

Smith, Christopher B.↗

Use of networked workstations for parallel nonlinear structural dynamic simulations of rotating bladed-disk assemblies

The principal objective of this research is to investigate, develop and demonstrate coarse-grained, parallel-processing strategies for nonlinear dynamic simulations for rotating bladed-disk assemblies. The parallel -processing strategies addressed include numerical algorithms for parallel nonlinear solutions and techniques to effect load balancing among processors. The parallel environment employed is a distributed-memory, coarse-grained one consisting of networked workstations. A parallel explicit time integration method has been implemented for transient nonlinear solutions of rotationg bladed-disk assemblies. Automatic domain partitioning techniques have been investigated for load balancing among processors. Advanced computing environments, data structures and interactive computer graphics all contribute to an integrated parallel finite element analysis system to facilitate more efficient and powerful dynamic simulations.

Hsieh, Shang-Hsien↗

An Occultation by Saturn's Rings on 1991 October 2-3 October 2-3 Observed with the Hubble Space Telescope

An occultation of the star GSC 6323-01396 (V = 11.9) by Saturn's rings was observed with the High-Speed Photometer on the Hubble Space Telescope (HST) on 1991 October 2-3. This occultation occurred when Saturn was near a stationary point, so the apparent motion of Saturn relative to the star was dominated by the HST orbital motion (8 km/s). Data were recorded simultaneously at effective wavelengths of 3200 and 7500 A, with an integration time of 0.15 s. Fifteen segments of occultation data, totaling 6.8 h, were recorded in 13 successive orbits during the 20.0 h interval from UTC 1991 October 2, 19:35 until UTC 1991 October 3, 15:35. Occultations by 43 different features throughout the classical rings were unambiguously identified in the light curve, with a second occultation by 24 of them occurring due to spacecraft orbital parallax during this extremely slow event. Occultation times for features currently presumed circular were measured and employed in a geometrical model for the rings. This model, relating the observed occultation times to feature radii and longitudes, is presented here and is used in a least-squares fit for the pole direction and radius scale of Saturn's ring system.

Elliot, J. L.↗

Airborne measurements of sulfur dioxide, dimethyl sulfide, carbon disulfide, and carbonyl sulfide by isotope dilution gas chromatography/mass spectrometry

A gas chromatograph/mass spectrometer is described for determining atmospheric sulfur dioxide, carbon disulfide, dimethyl sulfide, and carbonyl sulfide from aircraft and ship platforms. Isotopically labelled variants of each analyte were used as internal standards to achieve high precision. The lower limit of detection for each species for an integration time of 3 min was 1 pptv for sulfur dioxide and dimethyl sulfide and 0.2 pptv for carbon disulfide and carbonyl sulfide. All four species were simultaneously determined with a sample frequency of one sample per 6 min or greater. When only one or two species were determined, a frequency of one sample per 4 min was achieved. Because a calibration is included in each sample, no separate calibration sequence was needed. Instrument warmup was only a few minutes. The instrument was very robust in field deployments, requiring little maintenance.

Bandy, Alan R.↗

Measurement of tropospheric OH and HO2 by laser-induced fluorescence at low pressure

The hydroxyl radical (OH) is the primary oxidant in the atmosphere, responsible for many photochemical reactions that affect both regional air quality and global climate change. Because of its high reactivity, abundances of OH in the troposphere are less than 1 part per trillion by volume (pptv) and thus difficult to measure accurately. This paper describes an instrument for the sensitive detection of OH in the troposphere using low-pressure laser-induced fluorescence. Ambient air is expanded into a low pressure detection chamber, and OH is both excited and detected using the A(sup 2) Epsilon(+)(v prime = 0) yields X(sup 2)Pi(v double prime = 0) transition near 308 nm. An injector upstream of the detection axis allows for the addition of reagent NO to convert ambient HO2 to OH using the fast reaction HO2 + NO yields OH + NO2. Using recent advances in laser and detector technologies, this prototype instrument is able to detect less than 1 x 10(exp 5) molecules/cu cm (0.004 pptv) of OH with an integration time of 30 s with negligible interferences.

Stevens, P. S.↗

Plasmasphere dynamics in the duskside bulge region: A new look at old topic

Data acquired during several multiday periods in 1982 at ground stations Siple, Halley, and Kerguelen and on satellites Dynamics Explorer 1, International Sun Earth Explorer 1, and GEOS 2 have been used to investigate thermal plasma structure and dynamics in the duskside plasmasphere bulge region of the Earth. The distribution of thermal plasma in the dusk bulge sector is difficult to describe realistically, in part because of the time integral manner in which the thermal plasma distribution depends upon on the effects of bulk cross-B flow and interchange plasma flows along B. While relatively simple MHD models can be useful for qualitatively predicting certain effects of enhanced convection on a quiet plasmasphere, such as an initial sunward entrainment of the outer regions, they are of limited value in predicting the duskside thermal plasma structures that are observed. Furthermore, use of such models can be misleading if one fails to realize that they do not address the question of the formation of the steep plasmapause profile or provide for a possible role of instabilities or other irreversible processes in plasmapause formation. Our specific findings, which are based both upon the present case studies and upon earlier work, include the following: (1) during active periods the plasmasphere appears to become divided into two entities, a main plasmasphere and a duskside bulge region. (2) in the aftermath of an increase in convection activity, the main plasmasphere tends (from a statistical point of view) to become roughly circular in equatorial cross section, with only a slight bulge at dusk; (3) the abrupt westward edge of the duskside bulge observed from whistlers represents a state in the evolution of sunward extending streamers; (4) in the aftermath of a weak magnetic storm, 10 to 30% of the plasma 'removed' from the outer plasmasphere appears to remain in the afternoon-dusk sector beyond the main plasmasphere. (5) outlying dense plasma structures may circulate in the outer duskside magetosphere for many days following an increase in convection, unless there is extremely deep quieting; (6) a day-night plasmatrough boundary may be identified in equatorial satellite data; (7) factor-of-2-to-10 density irregularities appear near the plasmatrough from the ionosphere at L = 4.6, predominantly bidirectional field aligned and equatorially trapped light ion pitch angle distributions give away to a predominantly isotropic distribution (as seen by DE 1) when the plasma density reaches a level a factor of about 3 below the satured plasmasphere level; (9) some outlying dense plasma structures are effectively detached from the main plasmasphere, while others appear to be connected to that body.

Carpenter, D. L.↗

Flow simulations about steady-complex and unsteady moving configurations using structured-overlapped and unstructured grids

The limiting factor in simulating flows past realistic configurations of interest has been the discretization of the physical domain on which the governing equations of fluid flow may be solved. In an attempt to circumvent this problem, many Computational Fluid Dynamic (CFD) methodologies that are based on different grid generation and domain decomposition techniques have been developed. However, due to the costs involved and expertise required, very few comparative studies between these methods have been performed. In the present work, the two CFD methodologies which show the most promise for treating complex three-dimensional configurations as well as unsteady moving boundary problems are evaluated. These are namely the structured-overlapped and the unstructured grid schemes. Both methods use a cell centered, finite volume, upwind approach. The structured-overlapped algorithm uses an approximately factored, alternating direction implicit scheme to perform the time integration, whereas, the unstructured algorithm uses an explicit Runge-Kutta method. To examine the accuracy, efficiency, and limitations of each scheme, they are applied to the same steady complex multicomponent configurations and unsteady moving boundary problems. The steady complex cases consist of computing the subsonic flow about a two-dimensional high-lift multielement airfoil and the transonic flow about a three-dimensional wing/pylon/finned store assembly. The unsteady moving boundary problems are a forced pitching oscillation of an airfoil in a transonic freestream and a two-dimensional, subsonic airfoil/store separation sequence. Accuracy was accessed through the comparison of computed and experimentally measured pressure coefficient data on several of the wing/pylon/finned store assembly's components and at numerous angles-of-attack for the pitching airfoil. From this study, it was found that both the structured-overlapped and the unstructured grid schemes yielded flow solutions of comparable accuracy for these simulations. This study also indicated that, overall, the structured-overlapped scheme was slightly more CPU efficient than the unstructured approach.

Newman, James C., III↗

Analysis of an astrometric Fizeau interferometer for GAIA

The concepts related to the operation and design of the global astrometric interferometer for astrophysics (GAIA) bring together solutions chosen for the astrometry satellite and interferometric techniques. Like the Hipparcos satellite, GAIA is a continuously scanning instrument for which the integration time on any observed object is limited by the field of view of the detector. If a final astrometric accuracy of 10 microarcsec is aimed at, a field of 1 deg in diameter is needed. A design is presented for the proposed 2.6 m baseline Fizeau interferometer with two 40 cm apertures and overall dimensions compatible with the size of the Ariane 5 payload shroud. It has a 0.9 deg diffraction limited field of view. The response of the optical system to small perturbations on each optical element is given in terms of the fringe visibility, which is shown to be dependent on the sub-aperture spot separation. The robustness of the design to thermal, mechanical and manufacturing errors is discussed. The unavoidable distortion present in wide field optical systems is analyzed in terms of displacement of the interference fringes.

Loiseau, Sacha↗