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 1,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

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.↗
NASA NTRS1993

Flow downstream of the heliospheric terminal shock - The magnetic field on the heliopause

Modeling the kinematic magnetic field in the solar wind beyond the terminal shock shows that a ridge of magnetic pressure is produced just inside the heliopause. This ridge is sufficiently large that it will cause the layer immediately inside the heliopause to thicken, pushing the heliopause outward and slightly affecting its position relative to the terminal shock. However, the ridge is far too thin to cause an important change in the distance of the terminal shock from the sun. We show that these conclusions are a simple consequence of geometrical arguments for incompressible, steady, laminar flows. Moreover, the heliopause magnetic field originates on the terminal shock near the substagnation point. Consequently, the heliospheric current sheet field reversals are painted onto the inside surface of the heliopause. Alternate magnetic polarity strips will be oppositely directed relative to the interstellar magnetic field, implying that reconnection inevitably occurs on a fine some near the nose of the heliosphere. This suggests that the heliopause is a leaky, diffuse surface.

Nerney, Steven↗
NASA NTRS1980

WIND: Computer program for calculation of three dimensional potential compressible flow about wind turbine rotor blades

A computer program is presented which numerically solves an exact, full potential equation (FPE) for three dimensional, steady, inviscid flow through an isolated wind turbine rotor. The program automatically generates a three dimensional, boundary conforming grid and iteratively solves the FPE while fully accounting for both the rotating cascade and Coriolis effects. The numerical techniques incorporated involve rotated, type dependent finite differencing, a finite volume method, artificial viscosity in conservative form, and a successive line overrelaxation combined with the sequential grid refinement procedure to accelerate the iterative convergence rate. Consequently, the WIND program is capable of accurately analyzing incompressible and compressible flows, including those that are locally transonic and terminated by weak shocks. The program can also be used to analyze the flow around isolated aircraft propellers and helicopter rotors in hover as long as the total relative Mach number of the oncoming flow is subsonic.

Dulikravich, D. S.↗
NASA NTRS1981

Nonparallel stability of three-dimensional flows

The linear stability of three-dimensional incompressible, isothermal, nonparallel boundary-layer flows has been investigated. The method of multiple scales is used to derive the partial-differential equations that describe the spatial modulations of the amplitude, phase and wavenumber of a disturbance. Group velocities are used to determine the disturbance growth direction. The envelope method is used to calculate the logarithmic amplitude growth rate N. The theory is applied to the flows over a swept-back tapered wing with boundary-layer suction. Results of such analysis for the X-21 wing are discussed. It is found that the nonparallel effects for this wing is substantial.

Padhye, A. R.↗
NASA NTRS1991

Verification of the Proteus two-dimensional Navier-Stokes code for flat plate and pipe flows

The Proteus Navier-Stokes Code is evaluated for two-dimensional/axisymmetric, viscous, incompressible, internal and external flows. The particular cases to be discussed are laminar and turbulent flows over a flat plate, laminar and turbulent dveloping pipe flows and turbulent pipe flow with swirl. Results are compared with exact solutions, empirical correlations and experimental data. A detailed description of the code set-up, including boundary conditions, intitial conditions, grid size and grid packing is given for each case.

Conley, Julianne M.↗
NASA NTRS1993

Modeling bristle lift-off in idealized brush seal configurations

We attempt in this paper to develop a model for the flow through brush seals and determine their elastic behavior in order to predict the dependence of brush/journal clearance on geometry and operating conditions. Several idealizations regarding brush seal configuration, flow conditions, and elastic behavior are made in the analysis in order to determine closed form parametric dependence. This formulation assumes that there is no initial interference between the bristle tip and the rotor. Also, interbristle, bristle-backing plate, and bristle-rotor friction is neglected. The bristle bundle or the brush seal as it is alternately called is assumed homogeneous and isotropic on a macroscopic scale so that a physical property like permeability is uniform. The fluid is assumed to be homogeneous, incompressible, viscous, and flowing under steady conditions. A schematic of a brush seal is shown. If the nominal bristle-shaft interference is absent then under static conditions the bristles may deflect axially due to the imposed pressure differential. This axial deflection may create a clearance permitting leakage flow in excess of that which occurs through the porous matrix formed by the bristle bundles. Under dynamic conditions the Couette flow created by shaft motion could be strong enough to cause bristle deflection and once again a clearance may develop. The paper proposes a means to determine this clearance (or at least describe its parametric dependence on geometry and operating conditions) under static as well as dynamic conditions.

Modi, Vijay↗
NASA NTRS1990

A defect stream function formulation for compressible turbulent boundary layers

Progress to date on the development of a method for turbulent, wall-bounded flow which uses the defect stream function formulation in the outer layer and an analytic law of the wall and wake formulation in the inner region is reviewed. This two-formulation approach avoids the need to computationally resolve the high-gradient inner layer. One of the most appealing recent developments is the transformation of the compressible governing equation for the defect stream function into a linear, second-order differential equation which has analytic solutions for many problems of practical interest. Numerical and analytic results for incompressible and compressible flows are shown to be in excellent agreement with experimental results. In this paper the two-formulation approach is applied to primitive-variable computations. Excellent comparisons with experiment are presented for two compressible flat plate flows.

Barnwell, Richard W.↗
NASA NTRS1974

Steady simple shear flow past a circular cylinder at moderate Reynolds numbers - A numerical solution

The two-dimensional steady flow of an incompressible viscous fluid past a circular cylinder, placed symmetrically in a simple shear field, has been studied for both the stationary and the freely rotating case by solving numerically the Navier-Stokes equations for values of the Reynolds number R in the range from 0.047 to 70. At R = 0.047, the results obtained are in substantial agreement with the analytic small-R perturbation solution given by Robertson and Acrivos (1970). Inertia effects were found, however, to play a significant role even at R = 1, and hence the calculated flow pattern for R greater than or equal to 1 differs significantly from that of the creeping-flow solution. Specifically, for the freely rotating case, the region of closed streamlines decreases rapidly in extent with increasing R, two symmetrically placed wakes are formed on either side of the cylinder, and the dimensionless rotational speed of the freely suspended cylinder decreases as the reciprocal of the square root of R.

Kossack, C. A.↗
NASA NTRS1994

A spectrally accurate boundary-layer code for infinite swept wings

This report documents the development, validation, and application of a spectrally accurate boundary-layer code, WINGBL2, which has been designed specifically for use in stability analyses of swept-wing configurations. Currently, we consider only the quasi-three-dimensional case of an infinitely long wing of constant cross section. The effects of streamwise curvature, streamwise pressure gradient, and wall suction and/or blowing are taken into account in the governing equations and boundary conditions. The boundary-layer equations are formulated both for the attachment-line flow and for the evolving boundary layer. The boundary-layer equations are solved by marching in the direction perpendicular to the leading edge, for which high-order (up to fifth) backward differencing techniques are used. In the wall-normal direction, a spectral collocation method, based upon Chebyshev polynomial approximations, is exploited. The accuracy, efficiency, and user-friendliness of WINGBL2 make it well suited for applications to linear stability theory, parabolized stability equation methodology, direct numerical simulation, and large-eddy simulation. The method is validated against existing schemes for three test cases, including incompressible swept Hiemenz flow and Mach 2.4 flow over an airfoil swept at 70 deg to the free stream.

Pruett, C. David↗
NASA NTRS1972

An analytical investigation of the impingement jets on curved deflectors

Numerical solutions are obtained for the cases of straight circular jets impinging on axisymmetric curved surfaces and plane jets impinging symmetrically on two-dimensional curved surfaces. These geometries are representative of some types of thrust reversers for transport aircraft. The solutions are based on the assumptions of incompressible and potential flow. The velocity field, pressure distribution at the deflector surface and reverser effectiveness are predicted for deflector turning angles of 15 to 75 deg, deflector width to jet diameter ratios of 1.5 to 2.0, and ratios of deflector clearance to jet diameter of 1.0 to 3.0. Reverser effectiveness is found to be a maximum for a ratio of deflector clearance to jet diameter of about 2.0. The effect of back pressuring due to the presence of the deflector is predicted. Experimental verification of the theoretical predictions is obtained. A compressible solution obtained for a limited number of cases indicates that the incompressible solution is satisfactory for jet exit Mach numbers less than 0.8.

Schnurr, N. M.↗
NASA NTRS1990

Experimental investigation of turbulent flow through a circular-to-rectangular transition duct

An incompressible, turbulent, swirl-free flow through a circular-to-rectangular transition duct was studied experimentally. The cross-sectional geometry all along the duct was defined using the equation of a superellipse. The three mean velocity components and the six Reynolds stress components were measured at two axial stations downstream from the transition. It is shown that a secondary flow vortex pair which develops along the duct sidewalls significantly distorts the mean and turbulence fields. At the duct exit, the flow is not in local equilibrium, but recovers to local equilibrium conditions in the rectangular extension duct. Analysis demonstrates that conventional wall functions are not applicable at all streamwise locations in the duct.

Davis, D. O.↗
NASA NTRS1991

Experimental investigation of turbulent flow through a circular-to-rectangular transition duct

Steady, incompressible, turbulent, swirl-free flow through a circular-to-rectangular transition duck was studied experimentally. The cross-sectional area remains the same at the exit as at the inlet, but varies through the transition section to a maximum value approximately 15 percent above the inlet value. The cross-sectional geometry everywhere along the duct is defined by the equation of a superellipse. Mean and turbulence data were accumulated utilizing pressure and hot-wire instrumentation at five stations along the test section. Data are presented for operating bulk Reynolds numbers of 88,000 and 390,000. Measured quantities include total and static pressure, the three components of the mean velocity vector, and the six components of the Reynolds stress tensor. In addition to the transition duct measurements, a hot-wire technique which relies on the sequential use of single rotatable normal and slant-wire probes was proposed. The technique is applicable for measurement of the total mean velocity vector and the complete Reynolds stress tensor when the primary flow is arbitrarily skewed relative to a plane which lies normal to the probe axis of rotation.

Davis, David O.↗
NASA NTRS1981

Stability of flow over plates with porous suction strips

This paper addresses the stability of two-dimensional, incompressible boundary-layer flow over plates with suction through porous strips. The mean flow is calculated using linearized triple-deck, closed-form solutions. The stability results of the triple-deck theory are shown to be in good agreement with those of the interacting boundary layers. Then different configurations of number, spacing, and mass flow rate through such porous strips are analyzed and compared with nonsimilar uniform-suction stability results from the point of view of applicability to laminar flow control.

Reed, H. L.↗
DOE OSTI2020

High-order partitioned spectral deferred correction solvers for multiphysics problems

We present an arbitrarily high-order, conditionally stable, partitioned spectral deferred correction (SDC) method for solving multiphysics problems using a sequence of pre-existing single-physics solvers. This method extends the work in [1], [2], which used implicit-explicit Runge-Kutta methods (IMEX) to build high-order, partitioned multiphysics solvers. We consider a generic multiphysics problem modeled as a system of coupled ordinary differential equations (ODEs), coupled through coupling terms that can depend on the state of each subsystem; therefore the method applies to both a semi-discretized system of partial differential equations (PDEs) or problems naturally modeled as coupled systems of ODEs. The sufficient conditions to build arbitrarily high-order partitioned SDC schemes are derived. Based on these conditions, various of partitioned SDC schemes are designed. The stability of the first-order partitioned SDC scheme is analyzed in detail on a coupled, linear model problem. We show that the scheme is conditionally stable, and under conditions on the coupling strength, the scheme can be unconditionally stable. We demonstrate the performance of the proposed partitioned solvers on several classes of multiphysics problems with moderate coupling strength. They include a stiff linear system of ODEs, advection-diffusion-reaction systems, and fluid-structure interaction problems with both incompressible and compressible flows, where we verify the design order of the SDC schemes and study various stability properties. We also directly compare the accuracy, stability, and cost of the proposed partitioned SDC solver with the partitioned IMEX method in [1], [2] on this suite of test problems. The results suggest that the high-order partitioned SDC solvers are more robust than the partitioned IMEX solvers for the numerical examples considered in this work, while the IMEX methods require fewer implicit solves.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

FreeMHD: Validation and verification of the open-source, multi-domain, multi-phase solver for electrically conductive flows

The extreme heat fluxes in the divertor region of tokamaks may require an alternative to solid plasma-facing components, for the extraction of heat and the protection of the surrounding walls. Flowing liquid metals are proposed as an alternative, but raise additional challenges that require investigation and numerical simulations. Free surface designs are desirable for plasma-facing components, but steady flow profiles and surface stability must be ensured to limit undesirable interactions with the plasma. Previous studies have mainly used steady-state, 2D, or simplified models for internal flows and have not been able to adequately model free-surface liquid metal (LM) experiments. Therefore, FreeMHD has been recently developed as an open-source magnetohydrodynamics (MHD) solver for free-surface electrically conductive flows subject to a strong external magnetic field. The FreeMHD solver computes incompressible free-surface flows with multi-region coupling for the investigation of MHD phenomena involving fluid and solid domains. The model utilizes the finite-volume OpenFOAM framework under the low magnetic Reynolds number approximation. FreeMHD is validated using analytical solutions for the velocity profiles of closed channel flows with various Hartmann numbers and wall conductance ratios. Next, experimental measurements are then used to verify FreeMHD, through a series of cases involving dam breaking, 3D magnetic fields, and free-surface LM flows. These results demonstrate that FreeMHD is a reliable tool for the design of LM systems under free surface conditions at the reactor scale. Furthermore, it is flexible, computationally inexpensive, and can be used to solve fully 3D transient MHD flows.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

SIERRA Multimechanics Module: Aria Thermal Theory Manual - Version 4.58

Aria is a Galerkin finite element based program for solving coupled-physics problems described by systems of PDEs and is capable of solving nonlinear, implicit, transient and direct-to-steady state problems in two and three dimensions on parallel architectures. The suite of physics currently supported by Aria includes thermal energy transport, species transport, and electrostatics as well as generalized scalar, vector and tensor transport equations. Additionally, Aria includes support for manufacturing process flows via the incompressible Navier-Stokes equations specialized to a low Reynolds number (R e < 1) regime. Enhanced modeling support of manufacturing processing is made possible through use of either arbitrary Lagrangian-Eulerian (ALE) and level set based free and moving boundary tracking in conjunction with quasi-static nonlinear elastic solid mechanics for mesh control. Coupled physics problems are solved in several ways including fully-coupled Newtons method with analytic or numerical sensitivities, fully-coupled Newton-Krylov methods and a loosely-coupled nonlinear iteration about subsets of the system that are solved using combinations of the aforementioned methods. Error estimation, uniform and dynamic-adaptivity and dynamic load balancing are some of Arias more advanced capabilities.

42 ENGINEERING↗
DOE OSTI2020

SIERRA Multimechanics Module: Aria User Manual - Version 4.58

Aria is a Galerkin finite element based program for solving coupled-physics problems described by systems of PDEs and is capable of solving nonlinear, implicit, transient and direct-to-steady state problems in two and three dimensions on parallel architectures. The suite of physics currently supported by Aria includes thermal energy transport, species transport, and electrostatics as well as generalized scalar, vector and tensor transport equations. Additionally, Aria includes support for manufacturing process flows via the incompressible Navier-Stokes equations specialized to a low Reynolds number ( %3C 1 ) regime. Enhanced modeling support of manufacturing processing is made possible through use of either arbitrary Lagrangian-Eulerian (ALE) and level set based free and moving boundary tracking in conjunction with quasi-static nonlinear elastic solid mechanics for mesh control. Coupled physics problems are solved in several ways including fully-coupled Newtons method with analytic or numerical sensitivities, fully-coupled Newton-Krylov methods and a loosely-coupled nonlinear iteration about subsets of the system that are solved using combinations of the aforementioned methods. Error estimation, uniform and dynamic -adaptivity and dynamic load balancing are some of Arias more advanced capabilities.

97 MATHEMATICS AND COMPUTING↗