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 1,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

A simplified Reynolds stress model for unsteady turbulent boundary layers

A simplified Reynolds stress model has been developed for the prediction of unsteady turbulent boundary layers. By assuming that the net transport of Reynolds stresses is locally proportional to the net transport of the turbulent kinetic energy, the time dependent full Reynolds stress model is reduced to a set of ordinary differential equations. These equations contain only time derivatives and can be readily integrated in a time dependent boundary layer or Navier-Stokes code. The turbulent kinetic energy and dissipation rate needed for the model are obtained by solving the k-epsilon equations. This simplified Reynolds stress turbulence model (SRSM) does not use the eddy viscosity assumption, which may not be valid for unsteady turbulent flows. The anisotropy of both the steady and the unsteady turbulent normal stresses can be captured by the SRSM model. Through proper damping of the shear stresses, the present model can be used in the near wall region of turbulent boundary layers. This model has been validated against data for steady and unsteady turbulent boundary layers, including periodic turbulent boundary layers subjected to a mean adverse pressure gradient. For the cases tested, the predicted unsteady velocity and turbulent stress components agree well with the experimental data. Comparison between the predictions from the SRSM model and a k-epsilon model is also presented.

Fan, Sixin↗
NASA NTRS1993

Three-dimensional supersonic vortex breakdown

Three-dimensional supersonic vortex-breakdown problems in bound and unbound domains are solved. The solutions are obtained using the time-accurate integration of the unsteady, compressible, full Navier-Stokes (NS) equations. The computational scheme is an implicit, upwind, flux-difference splitting, finite-volume scheme. Two vortex-breakdown applications are considered in the present paper. The first is for a supersonic swirling jet which is issued from a nozzle into a supersonic uniform flow at a lower Mach number than that of the swirling jet. The second is for a supersonic swirling flow in a configured circular duct. In the first application, an extensive study of the effects of grid fineness, shape and grid-point distribution on the vortex breakdown is presented. Four grids are used in this study and they show a substantial dependence of the breakdown bubble and shock wave on the grid used. In the second application, the bubble-type and helix-type vortex breakdown have been captured.

Kandil, Osama A.↗
NASA NTRS1993

Three-dimensional unstructured grid Euler method applied to turbine blades

Flow through a turbine annular cascade is calculated using a three-dimensional Euler method based on unstructured tetrahedral meshes. The equations are integrated in time using an explicit Runge-Kutta time-stepping scheme. The inviscid flux terms are discretized using a cell-centered finite-volume formulation with upwind flux-difference splitting. The tetrahedral meshes around the turbine blade are generated using an advancing-front technique with forced geometric periodicity between the blades. Good agreement is obtained between the present calculation and the experiment for both surface pressure distribution and flow behavior in the passage between the blades, demonstrating the capability of the present methodology for turbomachinery flow applications.

Kwon, Oh J.↗
NASA NTRS1993

Hypersonic flutter of a curved shallow panel with aerodynamic heating

The general equations describing the nonlinear fluttering oscillations of shallow, curved, heated orthotropic panels have been derived. The formulation takes into account the location of the panel on the surface of a generic hypersonic vehicle, when calculating the aerodynamic loads. It is also shown that third order piston theory produces unsteady aerodynamic loading which is in close agreement with that based upon direct solution of the Euler equations. Results, for simply supported panels, are obtained using Galerkin's method combined with direct numerical integration in time to compute stable limit cycle amplitudes. These results illustrate the sensitivity of the aeroelastic behavior to the unsteady aerodynamic assumptions, temperature, orthotropicity and flow orientation.

Bein, T.↗
NASA NTRS1993

A Navier-Stokes solver for turbomachinery applications

A computer code for solving the Reynolds-averaged full Navier-Stokes equations has been developed and applied using H- and C-type grids. The Baldwin-Lomax eddy-viscosity model is used for turbulence closure. The integration in time is based on an explicit four-stage Runge-Kutta scheme. Local time stepping, variable coefficient implicit residual smoothing, and a full multigrid method have been implemented to accelerate steady-state calculations. A grid independence analysis is presented for a transonic rotor blade. Comparisons with experimental data show that the code is an accurate viscous solver and can give very good blade-to-blade predictions for engineering applications.

Arnone, A.↗
NASA NTRS1993

Two-dimensional CFD modeling of wave rotor flow dynamics

A two-dimensional Navier-Stokes solver developed for detailed study of wave rotor flow dynamics is described. The CFD model is helping characterize important loss mechanisms within the wave rotor. The wave rotor stationary ports and the moving rotor passages are resolved on multiple computational grid blocks. The finite-volume form of the thin-layer Navier-Stokes equations with laminar viscosity are integrated in time using a four-stage Runge-Kutta scheme. The Roe approximate Riemann solution scheme or the computationally less expensive Advection Upstream Splitting Method (AUSM) flux-splitting scheme are used to effect upwind-differencing of the inviscid flux terms, using cell interface primitive variables set by MUSCL-type interpolation. The diffusion terms are central-differenced. The solver is validated using a steady shock/laminar boundary layer interaction problem and an unsteady, inviscid wave rotor passage gradual opening problem. A model inlet port/passage charging problem is simulated and key features of the unsteady wave rotor flow field are identified. Lastly, the medium pressure inlet port and high pressure outlet port portion of the NASA Lewis Research Center experimental divider cycle is simulated and computed results are compared with experimental measurements. The model accurately predicts the wave timing within the rotor passage and the distribution of flow variables in the stationary inlet port region.

Welch, Gerard E.↗
NASA NTRS1993

Study and modeling of finite rate chemistry effects in turbulent non-premixed flames

The development of numerical models that reflect some of the most important features of turbulent reacting flows requires information about the behavior of key quantities in well defined combustion regimes. In turbulent flames, the coupling between turbulent and chemical processes is so strong that it is extremely difficult to isolate the role played by one individual physical phenomenon. Direct numerical simulation (hereafter DNS) allows us to study in detail the turbulence-chemistry interaction in some restricted but completely defined situations. Globally, non-premixed flames are controlled by two limiting regimes: the fast chemistry case, where the turbulent flame can be pictured as a random distribution of local chemical equilibrium problems; and the slow chemistry case, where the chemistry integrates in time the turbulent fluctuations. The Damkoehler number, ratio of a mechanical time scale to chemical time scale, is used to distinguish between these regimes. Today most of the industrial computer codes are able to perform predictions in the hypothesis of local equilibrium chemistry using a presumed shape for the probability density function (pdt) of the conserved scalar. However, the finite rate chemistry situation is of great interest because industrial burners usually generate regimes in which, at some points, the flame is undergoing local extinction or at least non-equilibrium situations. Moreover, this variety of situations strongly influences the production of pollutants. To quantify finite rate chemistry effect, the interaction between a non-premixed flame and a free decaying turbulence is studied using DNS. The attention is focused on the dynamic of extinction, and an attempt is made to quantify the effect of the reaction on the small scale mixing process. The unequal diffusivity effect is also addressed. Finally, a simple turbulent combustion model based on the DNS observations and tractable in real flow configurations is proposed.

Vervisch, Luc↗
NASA NTRS1993

Investigate wave-mean flow interaction and transport in the extratropical winter stratosphere

The grant supported studies using several models along with observations in order to investigate some questions of wave-mean flow interaction and transport in the extratropical winter stratosphere. A quasi-geostrophic wave model was used to investigate the possibility that resonant growth of planetary wave 2 may have played a role in the sudden stratospheric warming of February 1979. The results of the time-dependent integration support the interpretation of resonance during February, 1979. Because of the possibility that the model treatment of critical line interactions exerted a controlling influence on the atmospheric dynamics, a more accurate model was needed for wave-mean flow interaction studies. A new model was adapted from the 3-dimensional primitive equation model developed by K. Rose and G. Brasseur. In its present form the model is global, rather than hemispheric; it contains an infrared cooling algorithm and a parameterized solar heating; it has parameterized gravity wave drag; and the chemistry has been entirely revised.

Smith, Anne K.↗
NASA NTRS1994

Two-dimensional CFD modeling of wave rotor flow dynamics

A two-dimensional Navier-Stokes solver developed for detailed study of wave rotor flow dynamics is described. The CFD model is helping characterize important loss mechanisms within the wave rotor. The wave rotor stationary ports and the moving rotor passages are resolved on multiple computational grid blocks. The finite-volume form of the thin-layer Navier-Stokes equations with laminar viscosity are integrated in time using a four-stage Runge-Kutta scheme. Roe's approximate Riemann solution scheme or the computationally less expensive advection upstream splitting method (AUSM) flux-splitting scheme is used to effect upwind-differencing of the inviscid flux terms, using cell interface primitive variables set by MUSCL-type interpolation. The diffusion terms are central-differenced. The solver is validated using a steady shock/laminar boundary layer interaction problem and an unsteady, inviscid wave rotor passage gradual opening problem. A model inlet port/passage charging problem is simulated and key features of the unsteady wave rotor flow field are identified. Lastly, the medium pressure inlet port and high pressure outlet port portion of the NASA Lewis Research Center experimental divider cycle is simulated and computed results are compared with experimental measurements. The model accurately predicts the wave timing within the rotor passages and the distribution of flow variables in the stationary inlet port region.

Welch, Gerard E.↗
NASA NTRS1994

Dynamic analysis of a magnetic bearing system with flux control

Using measured values of two-dimensional forces in a magnetic actuator, equations of motion for an active magnetic bearing are presented. The presence of geometric coupling between coordinate directions causes the equations of motion to be nonlinear. Two methods are used to examine the unbalance response of the system: simulation by direct integration in time; and determination of approximate steady state solutions by harmonic balance. For relatively large values of the derivative control coefficient, the system behaves in an essentially linear manner, but for lower values of this parameter, or for higher values of the coupling coefficient, the response shows a split of amplitudes in the two principal directions. This bifurcation is sensitive to initial conditions. The harmonic balance solution shows that the separation of amplitudes actually corresponds to a change in stability of multiple coexisting solutions.

Knight, Josiah↗
NASA NTRS1994

Development of advanced Navier-Stokes solver

The objective of research was to develop and validate new computational algorithms for solving the steady and unsteady Euler and Navier-Stokes equations. The end-products are new three-dimensional Euler and Navier-Stokes codes that are faster, more reliable, more accurate, and easier to use. The three-dimensional Euler and full/thin-layer Reynolds-averaged Navier-Stokes equations for compressible/incompressible flows are solved on structured hexahedral grids. The Baldwin-Lomax algebraic turbulence model is used for closure. The space discretization is based on a cell-centered finite-volume method augmented by a variety of numerical dissipation models with optional total variation diminishing limiters. The governing equations are integrated in time by an implicit method based on lower-upper factorization and symmetric Gauss-Seidel relaxation. The algorithm is vectorized on diagonal planes of sweep using two-dimensional indices in three dimensions. Convergence rates and the robustness of the codes are enhanced by the use of an implicit full approximation storage multigrid method.

Yoon, Seokkwan↗
NASA NTRS1996

Two-Dimensional Computational Model for Wave Rotor Flow Dynamics

A two-dimensional (theta,z) Navier-Stokes solver for multi-port wave rotor flow simulation is described. The finite-volume form of the unsteady thin-layer Navier-Stokes equations are integrated in time on multi-block grids that represent the stationary inlet and outlet ports and the moving rotor passages of the wave rotor. Computed results are compared with three-port wave rotor experimental data. The model is applied to predict the performance of a planned four-port wave rotor experiment. Two-dimensional flow features that reduce machine performance and influence rotor blade and duct wall thermal loads are identified. The performance impact of rounding the inlet port wall, to inhibit separation during passage gradual opening, is assessed.

Welch, Gerard E.↗
NASA NTRS1994

Static aeroelastic analysis of wings using Euler/Navier-Stokes equations coupled with improved wing-box finite element structures

Strong interactions between flow about an aircraft wing and the wing structure can result in aeroelastic phenomena which significantly impact aircraft performance. Time-accurate methods for solving the unsteady Navier-Stokes equations have matured to the point where reliable results can be obtained with reasonable computational costs for complex non-linear flows with shock waves, vortices and separations. The ability to combine such a flow solver with a general finite element structural model is key to an aeroelastic analysis in these flows. Earlier work involved time-accurate integration of modal structural models based on plate elements. A finite element model was developed to handle three-dimensional wing boxes, and incorporated into the flow solver without the need for modal analysis. Static condensation is performed on the structural model to reduce the structural degrees of freedom for the aeroelastic analysis. Direct incorporation of the finite element wing-box structural model with the flow solver requires finding adequate methods for transferring aerodynamic pressures to the structural grid and returning deflections to the aerodynamic grid. Several schemes were explored for handling the grid-to-grid transfer of information. The complex, built-up nature of the wing-box complicated this transfer. Aeroelastic calculations for a sample wing in transonic flow comparing various simple transfer schemes are presented and discussed.

Guruswamy, Guru P.↗
NASA NTRS1995

Fluid flow volume measurements using a capacitance/conductance probe system

A probe system has been developed to measure the flow volume of a single fluid passing through an orifice or flow line. The system employs both capacitance and a conductance probe at the orifice, together with phase detection and data acquisition circuitry to measure flow volume and salinity under low or zero gravity conditions. A wide variety of frequencies can be used for the radio frequency (RF) signal source which is chosen primarily by the capacitance of the orifice probe and the fluid passing through the orifice. Rapid measurements are made using the reflected signal from the orifice probe to determine the 'instantaneous' permittivity of the fluid/gas mixture passing through. The 'instantaneous' measurements are integrated over time to determine flow volume. Analysis reveals that a narrow orifice helps to reduce non-linearities caused by differing flow rates. The geometry of 'deflectors' and 'directors' for the flowing fluid are important in obtaining linearity. Measured data shows that a volume measurement accuracy of approximately four percent can be consistently achieved. The prototype hardware system and associated software have been optimized and are available for further applications. The system has immediate application in low or zero gravity environments for measurements of urine or other liquid volumes.

Nguyen, T. X.↗
NASA NTRS1998

Long-Term file activity patterns in a UNIX workstation environment

As mass storage technology becomes more affordable for sites smaller than supercomputer centers, understanding their file access patterns becomes crucial for developing systems to store rarely used data on tertiary storage devices such as tapes and optical disks. This paper presents a new way to collect and analyze file system statistics for UNIX-based file systems. The collection system runs in user-space and requires no modification of the operating system kernel. The statistics package provides details about file system operations at the file level: creations, deletions, modifications, etc. The paper analyzes four months of file system activity on a university file system. The results confirm previously published results gathered from supercomputer file systems, but differ in several important areas. Files in this study were considerably smaller than those at supercomputer centers, and they were accessed less frequently. Additionally, the long-term creation rate on workstation file systems is sufficiently low so that all data more than a day old could be cheaply saved on a mass storage device, allowing the integration of time travel into every file system.

Gibson, Timothy J.↗
NASA NTRS1997

Instability Mechanisms of Thermally-Driven Interfacial Flows in Liquid-Encapsulated Crystal Growth

During the past year, a great deal of effort was focused on the enhancement and refinement of the computational tools developed as part of our previous NASA grant. In particular, the interface mollification algorithm developed earlier was extended to incorporate the effects of surface-rheological properties in order to allow the study of thermocapillary flows in the presence of surface contamination. These tools will be used in the computational component of the proposed research in the remaining years of this grant. A detailed description of the progress made in this area is provided elsewhere. Briefly, the method developed allows for the convection and diffusion of bulk-insoluble surfactants on a moving and deforming interface. The novelty of the method is its grid independence: there is no need for front tracking, surface reconstruction, body-fitted grid generation, or metric evaluations; these are all very expensive computational tasks in three dimensions. For small local radii of curvature there is a need for local grid adaption so that the smearing thickness remains a small fraction of the radius of curvature. A special Neumann boundary condition was devised and applied so that the calculated surfactant concentration has no variations normal to the interface, and it is hence truly a surface-defined quantity. The discretized governing equations are solved subsequently using a time-split integration scheme which updates the concentration and the shape successively. Results demonstrate excellent agreement between the computed and exact solutions.

Haj-Hariri, Hossein↗
NASA NTRS2000

Stellar Calibration of the WIC and SI Imagers and the GEO Photometers on IMAGE/FUV

The FUV instrument on the IMAGE spacecraft comprises three wide-field imagers, the Wide-Band Imaging Camera (WIC) of observing N2 Lyman-Birge-Hopfield (LBH) (140-190 nm) emissions and the Spectrographic Imager (SI), which has a 121.8 nm channel for observing red-shifted HI Lya photons and a 135.6 run channel for observing 01 135.6 nm emissions. In addition, three HI Lya photometers (GEO) are used to monitor the geocorona. The fields of view are 17 degrees x 17 degrees for the WIC imagers, 15 degrees x 15 degrees for the two SI imagers, and 10 diameter for the three GEO photometers. As the IMAGE spacecraft spins every 120 seconds, the GEO photometers sweep out circles on the sky (at 0 degrees and plus or minus 30 degrees with respect to the spacecraft spin plane), and the WIC and SI imagers use the Time Delay Integration (TDI) method to construct images centered on the Earth. Many FUV-bright stars are seen in the WIC, SI and even the GEO data. WE have used archived International Ultraviolet Explorer (IUE) far-ultraviolet flux spectra for 22 of the brightest of these stars to help refine the FUV instrumental sensitivities. The stars chosen range in spectral type form B0V to A11V, with magnitudes ranging from V- 1.3 (a Cru) to V=4.7 (G Cen) (although many more fainter stars are also seen). The initial results of this stellar calibration will be presented and compared with the pre-flight and dayglow modeling results.

Gladstone, R.↗
NASA NTRS1995

Absorbing Boundary Conditions For Optical Pulses In Dispersive, Nonlinear Materials

This paper will present results in computational nonlinear optics. An algorithm will be described that provides absorbing boundary conditions for optical pulses in dispersive, nonlinear materials. A new numerical absorber at the boundaries has been developed that is responsive to the spectral content of the pulse. Also, results will be shown of calculations of 2-D electromagnetic nonlinear waves computed by directly integrating in time the nonlinear vector Maxwell's equations. The results will include simulations of "light bullet" like pulses. Here diffraction and dispersion will be counteracted by nonlinear effects. Comparisons will be shown of calculations that use the standard boundary conditions and the new ones.

Goorjian, Peter M.↗