Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “incompressible flow”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Computation of saddle point of attachment

Low-speed flows over a cylinder mounted on a flat plate are studied numerically in order to confirm the existence of a saddle point of attachment in the flow before an obstacle, to analyze the flow characteristics near the saddle point theoretically, and to address the significance of the saddle point of attachment to the construction of external flow structures, the interpretation of experimental surface oil-flow patterns, and the theoretical definition of three-dimensional flow separation. Two numerical codes, one for an incompressible flow and another for a compressible flow, are used for various Mach numbers, Reynolds numbers, grid sizes, and numbers of grid points. It is pointed out that the potential presence of a saddle point of attachment means that a line of 'oil accumulation' from both sides of a skin-friction line emanating outward from a saddle point can be either a line of separation or a line of attachment.

Hung, Ching-Mao↗
NASA NTRS1994

The least-squares finite element method for low-mach-number compressible viscous flows

The present paper reports the development of the Least-Squares Finite Element Method (LSFEM) for simulating compressible viscous flows at low Mach numbers in which the incompressible flows pose as an extreme. Conventional approach requires special treatments for low-speed flows calculations: finite difference and finite volume methods are based on the use of the staggered grid or the preconditioning technique; and, finite element methods rely on the mixed method and the operator-splitting method. In this paper, however, we show that such difficulty does not exist for the LSFEM and no special treatment is needed. The LSFEM always leads to a symmetric, positive-definite matrix through which the compressible flow equations can be effectively solved. Two numerical examples are included to demonstrate the method: first, driven cavity flows at various Reynolds numbers; and, buoyancy-driven flows with significant density variation. Both examples are calculated by using full compressible flow equations.

Yu, Sheng-Tao↗
NASA NTRS2000

Boby-Vortex Interaction, Sound Generation and Destructive Interference

It is generally recognized that interaction of vortices with downstream blades is a major source of noise production. To analyze this problem numerically, a two-dimensional model of inviscid flow together with the method of matched asymptotic expansions is proposed. The method of matched asymptotic expansions is used to match the inner region of incompressible flow to the outer region of compressible flow. Because of incompressibility, relatively simple numerical methods are available to treat multiple vortices and multiple bodies of arbitrary shape. Disturbances from vortices and bodies propagate outward as sound waves. Due to their interactions, either constructive or destructive interference may result. When it is destructive, the combined sound intensity can be reduced, sometimes substantially. In addition, an analytical solution to sound generation by the cascade-vonex interaction is given.

Kao, Hsiao C.↗
NASA NTRS1979

Flow in nonrotating passages of radial inflow turbines

The analysis of irrotational incompressible flow field in the stator unit of a radial inflow turbine is presented. The solution in the combined scroll-nozzle assembly is complicated by the domain geometry and by its multiconnectivity. This model is necessary, however, in order to provide a better understanding of the mutual interaction effects of these two components on the flow field. The finite element method is used in the solution which is limited to the two dimensional case. A substructuring technique is adopted in the computational procedure and results in considerable savings in both computer time and core storage requirements. The results are presented for the flow velocity magnitude and direction in the scroll and through the various nozzles, for two nozzle blade geometries. In addition, the mass flow rates in the different nozzles are computed and their deviations from the mean value determined.

Baskharone, E.↗
NASA NTRS1992

Unsteady Shear Disturbances Within a Two Dimensional Stratified Flow

The origin and evolution of shear disturbances within a stratified, inviscid, incompressible flow are investigated numerically by a Clebsch/Weber decomposition based scheme. In contrast to homogeneous flows, within which vorticity can be redistributed but not generated, the presence of a density stratification can render an otherwise irrotational flow vortical. In this work, a kinematic decomposition of the unsteady Euler equations separates the unsteady velocity field into rotational and irrotational components. The subsequent evolution of these components is used to study the influence various velocity disturbances have on both stratified and homogeneous flows. In particular, the flow within a two-dimensional channel is used to investigate the evolution of rotational disturbances, generated or convected, downstream from an unsteady inflow condition. Contrasting simulations of both stratified and homogeneous flows are used to distinguish between redistributed inflow vorticity and that which is generated by a density stratification.

Yokota, Jeffrey W.↗
NASA NTRS1989

A spectral multi-domain technique applied to stability and transition of incompressible and compressible flows

Techniques for the numerical simulation of compressible and incompressible viscous flows are described and demonstrated. The focus is on a multidomain spectral method based on the flux-balance domain-interface scheme of Macaraeg and Streett (1986). The fundamental principles of the method are reviewed, and it is shown that the multidomain spectral approach facilitates the treatment of instability and the transition to turbulence. Typical numerical results for the Taylor-Couette problem and for bounded compressible high-speed mixing layers are presented in graphs and briefly characterized.

Macaraeg, Michele↗
NASA NTRS1979

Application of a velocity-split Navier-Stokes solution technique to external flow problems

A method for solving the Navier-Stokes equations based on splitting the velocity vector into its rotational and irrotational parts has recently been applied successfully to internal flow computations. In this paper, the applicability of the method to external flows is examined by studying two model problems; laminar incompressible flow over semiinfinite and finite flat plates. For these problems the procedure accurately reproduces the known solutions and is computationally very efficient even at high Reynolds numbers. Computational aspects of the method are discussed along with the possibility of using the procedure to retrofit a viscous capability into existing transonic potential flow codes.

Dwoyer, D. L.↗
NASA NTRS1998

Aerodynamic Design Optimization on Unstructured Meshes Using the Navier-Stokes Equations

A discrete adjoint method is developed and demonstrated for aerodynamic design optimization on unstructured grids. The governing equations are the three-dimensional Reynolds-averaged Navier-Stokes equations coupled with a one-equation turbulence model. A discussion of the numerical implementation of the flow and adjoint equations is presented. Both compressible and incompressible solvers are differentiated and the accuracy of the sensitivity derivatives is verified by comparing with gradients obtained using finite differences. Several simplifying approximations to the complete linearization of the residual are also presented, and the resulting accuracy of the derivatives is examined. Demonstration optimizations for both compressible and incompressible flows are given.

Nielsen, Eric J.↗
DOE OSTI2025

Assessment and Improvement of the SST-Gamma Transition Model in Nalu-Wind

We conduct laminar–turbulent boundary-layer transition simulations using a local correlation-based transition model for two-dimensional incompressible flow and present enhancements to improve the accuracy of transition predictions. Menter’s Galilean-invariant 𝛾 transition model is implemented in the incompressible, unstructured-grid flow solver Nalu-Wind and is validated against experimental data and results from NASA’s flow solvers. The test cases of the AIAA Transition Prediction and Modeling Workshop are investigated, namely, the T3A/T3B flat plates and the NLF(1)-0416 and S809 airfoils. Based on the results, best practices for transition simulations, particularly for an unstructured-grid flow solver, are identified. Additional airfoil simulations are conducted for two wind turbine airfoils, S822 at Reynolds numbers of 𝑂⁡(10 5 ) and DU00-W-212 at Reynolds numbers of 𝑂⁡(10 7 ), to assess the model at low and high Reynolds numbers. Furthermore, through this work, we propose several approaches to enhance transition simulations, including 1) enforcing positivity of the implicit operator for the source terms of the transition model, 2) employing a constant turbulence intensity in stationary external flow simulations, and 3) recommending meshing for unstructured-grid flow solvers. Finally, we provide detailed documentation of the validation and data for the canonical cases to the transition modeling community.

17 WIND ENERGY↗
NASA NTRS1989

Numerical Simulation of Flow Through an Artificial Heart

A solution procedure was developed that solves the unsteady, incompressible Navier-Stokes equations, and was used to numerically simulate viscous incompressible flow through a model of the Pennsylvania State artificial heart. The solution algorithm is based on the artificial compressibility method, and uses flux-difference splitting to upwind the convective terms; a line-relaxation scheme is used to solve the equations. The time-accuracy of the method is obtained by iteratively solving the equations at each physical time step. The artificial heart geometry involves a piston-type action with a moving solid wall. A single H-grid is fit inside the heart chamber. The grid is continuously compressed and expanded with a constant number of grid points to accommodate the moving piston. The computational domain ends at the valve openings where nonreflective boundary conditions based on the method of characteristics are applied. Although a number of simplifing assumptions were made regarding the geometry, the computational results agreed reasonably well with an experimental picture. The computer time requirements for this flow simulation, however, are quite extensive. Computational study of this type of geometry would benefit greatly from improvements in computer hardware speed and algorithm efficiency enhancements.

Rogers, Stuart E.↗
NASA NTRS1990

Experimental investigation of flow about a strut-endwall configuration

An experimental study was conducted to investigate incompressible flow about a strut-endwall configuration positioned within a constant area duct. The endwal boundary layer was tripped, but natural transition occurred on the strut surface. The results indicate that spanwise varying transition on this surface leads to the formation of a secondary vortex which coexists with the horseshoe vortex generated by endwall flow separation upstream of the strut leading edge. Both vortices are similar in strength downstream of the strut trailing edge, and both distort the primary flow and local turbulence structure in the wake-endwall region. The level of distortion is demonstrated by means of axial mean velocity contours, turbulence kinetic energy contours, and Reynolds shear stress contours measured in the cross plane at two streamwise locations. Analysis of the results shows that conventional eddy viscosity and k-epsilon transport equation models are not wholly adequate for predicting this flow situation.

Chang, P. S.↗
NASA NTRS1994

Low-dimensional description of the dynamics in separated flow past thick airfoils

Results are presented for the numerical simulation of unsteady viscous incompressible flow past thick airfoils. Specifically, flow past a NACA 4424 at an angle of attack of 2.5 deg and Reynolds numbers in the range of 1700-4000 has been simulated using the spectral element method. At these conditions the flow is separatedd and an unsteady wake is formed. Application of the method of empirical eigenfunction reveals the structure of the most energetic components of the flow. These are found to occur in pairs that, through phase exchange, are responsible for the vortex shedding. A set of ordinary differential equations is obtained for the amplitudes of these eigenfunctions by a Galerkin projection of the Navier-Stokes equations. The solutions of the model system are compared with the full simulation. The work is of relevance to the transition process and observed routes to chaos in airfoil wakes.

Deane, Anil E.↗
DOE OSTI2021

Ocean Energy Systems Wave Energy Modeling Task 10.4: Numerical Modeling of a Fixed Oscillating Water Column

This paper reports on an ongoing international effort to establish guidelines for numerical modeling of wave energy converters, initiated by the International Energy Agency Technology Collaboration Program for Ocean Energy Systems. Initial results for point absorbers were presented in previous work, and here we present results for a breakwater-mounted Oscillating Water Column (OWC) device. The experimental model is at scale 1:4 relative to a full-scale installation in a water depth of 12.8 m. The power-extracting air turbine is modeled by an orifice plate of 1–2% of the internal chamber surface area. Measurements of chamber surface elevation, air flow through the orifice, and pressure difference across the orifice are compared with numerical calculations using both weakly-nonlinear potential flow theory and computational fluid dynamics. Both compressible- and incompressible-flow models are considered, and the effects of air compressibility are found to have a significant influence on the motion of the internal chamber surface. Recommendations are made for reducing uncertainties in future experimental campaigns, which are critical to enable firm conclusions to be drawn about the relative accuracy of the numerical models. It is well-known that boundary element method solutions of the linear potential flow problem (e.g., WAMIT) are singular at infinite frequency when panels are placed directly on the free surface. This is problematic for time-domain solutions where the value of the added mass matrix at infinite frequency is critical, especially for OWC chambers, which are modeled by zero-mass elements on the free surface. A straightforward rational procedure is described to replace ad-hoc solutions to this problem that have been proposed in the literature.

16 TIDAL AND WAVE POWER↗
NASA NTRS1993

Prediction of stall and post-stall behavior of airfoils at low and high Reynolds numbers

An interactive boundary-layer method, together with the e(super n)-approach to the calculation of transition, has been used to predict the stall and post-stall behavior of airfoils at low and high Reynolds numbers. The turbulence model is based on the Cebeci-Smith algebraic eddy-viscosity formulation with improvements for strong pressure gradient effects and transitional flows at low Reynolds numbers. Comparison of calculated results for incompressible flows indicate good agreement with experiment for a wide range of Reynolds numbers. Preliminary calculations for low Mach number flows with this interactive method with compressibility corrections to the panel method indicate that, at a Mach number of 0.3, the compressibility effect on (C sub Q)max is not negligible.

Cebeci, T.↗
NASA NTRS1992

A coupled Eulerian/Lagrangian method for the solution of three-dimensional vortical flows

A coupled Eulerian/Lagrangian method is presented for the reduction of numerical diffusion observed in solutions of three-dimensional rotational flows using standard Eulerian finite-volume time-marching procedures. A Lagrangian particle tracking method using particle markers is added to the Eulerian time-marching procedure and provides a correction of the Eulerian solution. In turn, the Eulerian solutions is used to integrate the Lagrangian state-vector along the particles trajectories. The Lagrangian correction technique does not require any a-priori information on the structure or position of the vortical regions. While the Eulerian solution ensures the conservation of mass and sets the pressure field, the particle markers, used as 'accuracy boosters,' take advantage of the accurate convection description of the Lagrangian solution and enhance the vorticity and entropy capturing capabilities of standard Eulerian finite-volume methods. The combined solution procedures is tested in several applications. The convection of a Lamb vortex in a straight channel is used as an unsteady compressible flow preservation test case. The other test cases concern steady incompressible flow calculations and include the preservation of turbulent inlet velocity profile, the swirling flow in a pipe, and the constant stagnation pressure flow and secondary flow calculations in bends. The last application deals with the external flow past a wing with emphasis on the trailing vortex solution. The improvement due to the addition of the Lagrangian correction technique is measured by comparison with analytical solutions when available or with Eulerian solutions on finer grids. The use of the combined Eulerian/Lagrangian scheme results in substantially lower grid resolution requirements than the standard Eulerian scheme for a given solution accuracy.

Felici, Helene Marie↗
NASA NTRS2003

CFD-based Design of LOX Pump Inlet Duct for Reduced Dynamic Loads

Numerical simulations have been completed for a variety of designs for a 90 deg elbow duct. The objective is to identify a design that minimizes the dynamic load entering a LOX turbopump located at the elbow exit. Designs simulated to date indicate that simpler duct geometries result in lower losses. Benchmark simulations have verified that the compressible flow code used in this study is applicable to these incompressible flow simulations.

Rothermel, Jeffry↗