Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “INCOMPRESSIBLE FLUID”

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 199 records · Page 11

Pseudo-compressibility methods for the incompressible flow equations

Preconditioning methods to accelerate convergence to a steady state for the incompressible fluid dynamics equations are considered. The analysis relies on the inviscid equations. The preconditioning consists of a matrix multiplying the time derivatives. Thus the steady state of the preconditioned system is the same as the steady state of the original system. The method is compared to other types of pseudo-compressibility. For finite difference methods preconditioning can change and improve the steady state solutions. An application to viscous flow around a cascade with a non-periodic mesh is presented.

Turkel, Eli↗

Effect of an extendable slat on the stall behavior of a VR-12 airfoil

Experimental and computational tests were performed on a VR-12 airfoil to determine if the dynamic-stall behavior that normally accompanies high-angle pitch oscillations could be modified by segmenting the forward portion of the airfoil and extending it ahead of the main element. In the extended position the configuration would appear as an airfoil with a leading-edge slat, and in the retracted position it would appear as a conventional VR-12 airfoil. The calculations were obtained from a numerical code that models the vorticity transport equation for an incompressible fluid. These results were compared with test data from the water tunnel facility of the Aeroflightdynamics Directorate at Ames Research Center. Steady and unsteady flows around both airfoils were examined at angles of attack between 0 and 30 deg. The Reynolds number was fixed at 200,000 and the unsteady pitch oscillations followed a sinusoidal motion described by alpha = alpha(sub m) + 10 deg sin(omega t). The mean angle (alpha(sub m)) was varied from 10 to 20 deg and the reduced frequency from 0.05 to 0.20. The results from the experiment and the calculations show that the extended-slat VR-12 airfoil experiences a delay in both static and dynamic stall not experienced by the basic VR-12 airfoil.

Dehugues, P. Plantin↗

Towards enhancing and delaying disturbances in free shear flows

The family of shear flows comprising the jet, wake, and the mixing layer are subjected to perturbations in an inviscid incompressible fluid. By modeling the basic mean flows as parallel with piecewise linear variations for the velocities, complete and general solutions to the linearized equations of motion can be obtained in closed form as functions of all space variables and time when posed as an initial value problem. The results show that there is a continuous as well as the discrete spectrum that is more familiar in stability theory and therefore there can be both algebraic and exponential growth of disturbances in time. These bases make it feasible to consider control of such flows. To this end, the possibility of enhancing the disturbances in the mixing layer and delaying the onset in the jet and wake is investigated. It is found that growth of perturbations can be delayed to a considerable degree for the jet and the wake but, by comparison, cannot be enhanced in the mixing layer. By using moving coordinates, a method for demonstrating the predominant early and long time behavior of disturbances in these flows is given for continuous velocity profiles. It is shown that the early time transients are always algebraic whereas the asymptotic limit is that of an exponential normal mode. Numerical treatment of the new governing equations confirm the conclusions reached by use of the piecewise linear basic models. Although not pursued here, feedback mechanisms designed for control of the flow could be devised using the results of this work.

Criminale, W. O.↗

Particle-mesh techniques

This is an introduction to numerical Particle-Mesh techniques, which are commonly used to model plasmas, gravitational N-body systems, and both compressible and incompressible fluids. The theory behind this approach is presented, and its practical implementation, both for serial and parallel machines, is discussed. This document is based on a four-hour lecture course presented by the author at the NASA Summer School for High Performance Computational Physics, held at Goddard Space Flight Center.

Macneice, Peter↗

Toward a theory of interstellar turbulence. 2: Strong alfvenic turbulence

We continue to investigate the possibility that interstellar turbulence is caused by nonlinear interactions among shear Alfven waves. Here, we restrict attention to the symmetric case where the oppositely directed waves carry equal energy fluxes. This precludes application to the solar wind in which the outward flux significantly exceeds the ingoing one. All our detailed calculations are carried out for an incompressible magnetized fluid. In incompressible magnetohydrodynamics (MHD), nonlinear interactions only occur between oppositely direct waves. We present a theory for the strong turbulence of shear Alfven waves. It has the following main characteristics. (1) The inertial-stage energy spectrum exhibits a critical balance between linear wave periods and nonlinear turnover timescales. (2) The 'eddies' are elongated in the direction of the field on small spatial scales; the parallel and perpendicular components of the wave vector, k(sub z) and k(perpendicular) are related by k(sub z) approximately equals k(sub perpendicular to)(exp 2/3) L(exp -1/3), where L is the outer scale of the turbulence. (3) The 'one-dimensional' energy spectrum is proportional to k(sub perpendicular)(exp -5/3)-an anisotropic Kolmogorov energy spectrum. Shear Alfvenic turbulence mixes specific entropy as a passive contaminant. This gives rise to an electron density power spectrum whose form mimics the energy spectrum of the turbulence. Radio wave scattering by these electron density fluctuations produces anisotropic scatter-broadened images. Damping by ion-neutral collisions restricts Alfvenic turbulence to highly ionized regions of the interstellar medium.

Goldreich, P.↗

Upwind-Differencing Computations Of Incompressible Flows

Report discusses application of upwind-differencing numerical-integration scheme in conjunction with method of pseudocompressibility to computations of flows of incompressible fluids, with emphasis on steady-state flows. Three test cases presented: driven flow in square cavity, flow over backward-facing step, and flow around circular cylinder.

Rogers, Stuart↗

On the Nonlinear Evolution of a Stationary Cross-Flow Vortex in a Fully Three-Dimensional Boundary Layer Flow

We consider the nonlinear stability of a fully three-dimensional boundary layer flow in an incompressible fluid and derive an equation governing the nonlinear development of a stationary cross-flow vortex. The amplitude equation is a novel integro-differential equation which has spatial derivatives of the amplitude occurring in the kernal function. It is shown that the evolution of the cross-flow vortex is strongly coupled to the properties of an unsteady wall layer which is in fact driven by an unknown slip velocity, proportional to the amplitude of the cross-flow vortex. The work is extended to obtain the corresponding equation for rotating disk flow. A number of special cases are examined and the numerical solution for one of cases, and further analysis, demonstrates the existence of finite-distance as well as focussing type singularities. The numerical solutions also indicate the presence of a new type of nonlinear wave solution for a certain set of parameter values.

Gajjar, J. S. B.↗

The Evolution of Finite Amplitude Wavetrains in Plane Channel Flow

We consider a viscous incompressible fluid flow driven between two parallel plates by a constant pressure gradient. The flow is at a finite Reynolds number, with an 0(l) disturbance in the form of a traveling wave. A phase equation approach is used to discuss the evolution of slowly varying fully nonlinear two dimensional wavetrains. We consider uniform wavetrains in detail, showing that the development of a wavenumber perturbation is governed by Burgers equation in most cases. The wavenumber perturbation theory, constructed using the phase equation approach for a uniform wavetrain, is shown to be distinct from an amplitude perturbation expansion about the periodic flow. In fact we show that the amplitude equation contains only linear terms and is simply the heat equation. We review, briefly, the well known dynamics of Burgers equation, which imply that both shock structures and finite time singularities of the wavenumber perturbation can occur with respect to the slow scales. Numerical computations have been performed to identify areas of the (wavenumber, Reynolds number, energy) neutral surface for which each of these possibilities can occur. We note that the evolution equations will breakdown under certain circumstances, in particular for a weakly nonlinear secondary flow. Finally we extend the theory to three dimensions and discuss the limit of a weak spanwise dependence for uniform wavetrains, showing that two functions are required to describe the evolution. These unknowns are a phase and a pressure function which satisfy a pair of linearly coupled partial differential equations. The results obtained from applying the same analysis to the fully three dimensional problem are included as an appendix.

Hewitt, R. E.↗

Multigrid, Fractional-Step Computation Of Flow

Speed of computer code solving three-dimensional Navier-Stokes equations of flow of incompressible fluid by fractional-step method increased significantly by use of multigrid procedures. In method, equations solved on general nonorthogonal curvilinear coordinate grid, using volume fluxes. At each time step, computations performed in two fractional steps. In first step, equations of conservation of momentum solved by use of gradient of pressure from previous time step via explicit approximate-factorization method, yielding approximate flow field that does not satisfy equation of conservation of mass. In second step, discrete Poisson-like equation with Neumann-type boundary conditions, formed by combining equations of conservation of momentum and mass, solved iteratively.

Kwak, Dochan↗

Convergence Acceleration of the Navier-Stokes Equations Through Time-Derivative Preconditioning

Chorin's method of artificial compressibility is extended to both compressible and incompressible fluids by using physical arguments to define artificial fluid properties that make up a local preconditioning matrix. In particular, perturbation expansions are used to provide appropriate temporal derivatives for the equations of motion at both low speeds and low Reynolds numbers. These limiting forms are then combined into a single function that smoothly merges into the physical time derivatives at high speeds so that the equations are left unchanged at transonic, high Reynolds number conditions. The effectiveness of the resulting preconditioning procedures for the Navier-Stokes equations is demonstrated for a wide speed and Reynolds number ranges by means of stability results and computational solutions. Nevertheless, the preconditioned equations sometimes fail to provide a solution for applications for which the non-preconditioned equations converge. Often this is because the reduced dissipation in the preconditioned equations results in an unsteady solution while the more dissipative non-preconditioned equations result in a steady state. Problems of this type represent a computational challenge; it is important to distinguish between non-convergence of algorithms, and the non-existence of steady state solutions.

Merkle, Charles L.↗

Structural Stability and gamma-Transitions

In this work, we review some new discoveries and physical interpretations that concern the secular and dynamical instabilities and the associated structural changes in rotating, self-gravitating, incompressible fluids. For such fluids, we have been able to examine in depth and to understand physically a variety of classic results. Here we discuss the structural stability and the lambda-transitions of astrophysical fluids. We conclude with a summary of all the different types of lambda-transitions found in the course of our investigation and with a discussion of the implications of our results for superfluids, catastrophes, thermodynamical phase transitions, and the breaking of symmetry and topology.

Christodoulou, Dimitris M.↗

A Parallel Multilevel Spectral Element Scheme

A parallel multilevel strategy is developed using spectral (p) finite elements. Hierarchic bases are particularly well suited since the element matrices and vectors are nested and the multilevel projections easiliy performed. Since the basis degree is used to specify the multigrid level, an EBE strategy is natural br the multilevel technique. Results are presented for two candidate nonlinear elliptic transport problems: the augmented drift-diffusion equations of semiconductor device modeling and the stream function-vorticity equations of incompressible fluid dynamics.

Davis, M. B.↗

The Role of Multiphysics Simulation in Multidisciplinary Analysis

This article describes the applications of the Spectrum(Tm) Solver in Multidisciplinary Analysis (MDA). Spectrum, a multiphysics simulation software based on the finite element method, addresses compressible and incompressible fluid flow, structural, and thermal modeling as well as the interaction between these disciplines. Multiphysics simulation is based on a single computational framework for the modeling of multiple interacting physical phenomena. Interaction constraints are enforced in a fully-coupled manner using the augmented-Lagrangian method. Within the multiphysics framework, the finite element treatment of fluids is based on Galerkin-Least-Squares (GLS) method with discontinuity capturing operators. The arbitrary-Lagrangian-Eulerian method is utilized to account for deformable fluid domains. The finite element treatment of solids and structures is based on the Hu-Washizu variational principle. The multiphysics architecture lends itself naturally to high-performance parallel computing. Aeroelastic, propulsion, thermal management and manufacturing applications are presented.

Rifai, Steven M.↗

Multiphysics Simulation of Active Hypersonic Lip Cooling

This article describes the application of the Multidisciplinary Analysis (MDA) solver, Spectrum, in analyzing a hydrogen-cooled hypersonic cowl leading-edge structure. Spectrum, a multiphysics simulation code based on the finite element method, addresses compressible and incompressible fluid flow, structural, and thermal modeling, as well as the interactions between these disciplines. Fluid-solid-thermal interactions in a hydrogen impingement-cooled leading edge are predicted using Spectrum. Two- and semi-three-dimensional models are considered for a leading edge impingement coolant, concept under either specified external heat flux or aerothermodynamic heating from a Mach 5 external flow interaction. The solution accuracy is demonstrated from mesh refinement analysis. With active cooling, the leading edge surface temperature is drastically reduced from 1807 K of the adiabatic condition to 418 K. The internal coolant temperature profile exhibits a sharp gradient near channel/solid interface. Results from two different cooling channel configurations are also presented to illustrate the different behavior of alternative active cooling schemes.

Melis, Matthew E.↗

Towards a Modern Theory of Multiphase Filtration Flow

An alternative theoretical model of joint filtration flow of immiscible incompressible fluids is presented. The model takes into account relaxation processes due to the interchange of the fluids between pores of difference sizes which is driven by capillary forces. The fluids occupy connected regions in a four-dimensional space formed by three coordinates and the pore length scale. When the fluid exchange between pores of given sizes is effected by way of successive flow through pores of all the intermediate sizes, the pressure within each region is governed by a hyperbolic equation, the role of time being played by the pore linear scale. Pressure jumps across hypersurfaces separating the regions equal corresponding values of the capillary pressure. A supplementary condition at any such hypersurface requires the speed of its displacement in the four-dimensional space to coincide with the normal velocity components of both the adjoining fluids. As a result, a principally new statement of multiphase filtration flow problems is gained with allowance for capillary relaxation in the porous space.

Buyevich, Yu A.↗

Transient Growth Theory Prediction of Optimal Placing of Passive and Active Flow Control Devices for Separation Delay in LPT Airfoils

An analysis of the non-modal growth of perturbations in a boundary layer in the presence of a streamwise pressure gradient is presented. The analysis is based on PSE equations for an incompressible fluid. Examples with Falkner-Skan profiles indicate that a favorable pressure gradient decreases the non-modal growth while an unfavorable pressure gradient leads to an increase of the amplification. It is suggested that the transient growth mechanism be utilized to choose optimal parameters of tripping elements on a low-pressure turbine (LPT) airfoil. As an example, a boundary layer flow with a streamwise pressure gradient corresponding to the pressure distribution over a LPT airfoil is considered. It is shown that there is an optimal spacing of the tripping elements and that the transient growth effect depends on the starting point. At very low Reynolds numbers, there is a possibility to enhance the transient energy growth by means of wall cooling.

Tumin, Anatoli↗

Optimal Disturbances in Boundary Layers Subject to Streamwise Pressure Gradient

An analysis of the non-modal growth of perturbations in a boundary layer in the presence of a streamwise pressure gradient is presented. The analysis is based on PSE equations for an incompressible fluid. Examples with Falkner- Skan profiles indicate that a favorable pressure gradient decreases the non-modal growth while an unfavorable pressure gradient leads to an increase of the amplification. It is suggested that the transient growth mechanism be utilized to choose optimal parameters of tripping elements on a low-pressure turbine (LPT) airfoil. As an example, a boundary-layer flow with a streamwise pressure gradient corresponding to the pressure distribution over a LPT airfoil is considered. It is shown that there is an optimal spacing of the tripping elements and that the transient growth effect depends on the starting point. The amplification is found to be small at the LPT s very low Reynolds numbers, but there is a possibility to enhance the transient energy growth by means of wall cooling.

Ashpis, David E.↗

Effect of the Acceleration of Elongated Bodies of Revolution Upon the Resistance in a Compressible Flow

The problem of the motion of an elongated body of revolution in an incompressible fluid may, as is known, be solved approximately with the aid of the distribution of sources along the axis of the body. In determining the velocity field, the question of whether the body moves uniformly or with an acceleration is no factor in the problem. The presence of acceleration must be taken into account in determining the pressures acting on the body. The resistance of the body arising from the accelerated motion may be computed either directly on the basis of these pressures or with the aid of the so-called associated masses (inertia coefficients). A different condition holds in the case of the motion of bodies in a compressible gas. In this case the finite velocity of sound must be taken into account.

Frankl, F. I.↗