Search NASA⌕ Search

SEARCH · Search NASA

Results for “step 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 181 records · Page 10

A study of turbulent flow computations in an angled duct with a step

This paper presents a method for effective computation of turbulent flows in an angled duct which has a step at the angled section by using a grid generation technique. An effect of flow curvature was analyzed on the separating and reattaching flows. The computations are based on the Reynolds averaged transport equations which are transformed for use on a generalized curvilinear coordinate system. A special care was taken on the treatment of pressure corrections. Wall function boundary conditions for k-epsilon model were developed.

Maruszewski, J. P.↗

Computation of steady and unsteady vortex dominated flows

The unsteady Euler equations have been written in the conservation form for the relative motion in a rotating frame of reference. The resulting equations are solved by using a central-difference finite-volume scheme with four-stage Runge-Kutta time stepping. For steady flow problems, local time stepping is used and for unsteady flow problems, the minimum global time stepping is used. A three-dimensional fully vectorized computer program has been developed and applied to steady and unsteady maneuvering delta wing. The capability of the three-dimensional program has been demonstrated for a rigid sharp-edged delta wing undergoing uniform rolling in a conical flow and rolling oscillations in a locally conical flow.

Kandil, Osama A.↗

Turbulent atmospheric flow over a backward-facing step

The phenomenon of atmospheric shear layer separation over a man-made structure such as a building (modeled as a backward-facing step) has been analyzed theoretically by (1) solving the two-dimensional equations of motion in the two variables, stream function and vorticity, and by (2) employing an approximate integral technique. Boundary conditions for the undisturbed flow are that of the turbulent atmospheric shear flow over a rough terrain. In the first approach a two-equation model of turbulence was used. In the second approach an approximate technique was utilized in an attempt to describe the details of the flow in the recirculation zone behind the step. The results predict velocity profiles in sufficient detail that the presence of the corner eddy in the region of negative surface pressure gradient is evident. The magnitude of the reversed flow velocity in the recirculation eddy has been found to agree with that found from experiments. Also, a surface eddy viscosity distribution has been an outgrowth of the method which realistically follows the magnitude of the surface pressure gradient distribution as found experimentally.

Kaul, U. K.↗

High Energy Boundary Conditions for a Cartesian Mesh Euler Solver

Inlets and exhaust nozzles are common place in the world of flight. Yet, many aerodynamic simulation packages do not provide a method of modelling such high energy boundaries in the flow field. For the purposes of aerodynamic simulation, inlets and exhausts are often fared over and it is assumed that the flow differences resulting from this assumption are minimal. While this is an adequate assumption for the prediction of lift, the lack of a plume behind the aircraft creates an evacuated base region thus effecting both drag and pitching moment values. In addition, the flow in the base region is often mis-predicted resulting in incorrect base drag. In order to accurately predict these quantities, a method for specifying inlet and exhaust conditions needs to be available in aerodynamic simulation packages. A method for a first approximation of a plume without accounting for chemical reactions is added to the Cartesian mesh based aerodynamic simulation package CART3D. The method consists of 3 steps. In the first step, a components approach where each triangle is assigned a component number is used. Here, a method for marking the inlet or exhaust plane triangles as separate components is discussed. In step two, the flow solver is modified to accept a reference state for the components marked inlet or exhaust. In the third step, the flow solver uses these separated components and the reference state to compute the correct flow condition at that triangle. The present method is implemented in the CART3D package which consists of a set of tools for generating a Cartesian volume mesh from a set of component triangulations. The Euler equations are solved on the resulting unstructured Cartesian mesh. The present methods is implemented in this package and its usefulness is demonstrated with two validation cases. A generic missile body is also presented to show the usefulness of the method on a real world geometry.

Pandya, Shishir↗

Effect of wall suction on the stability of compressible subsonic flows over smooth two-dimensional backward-facing steps

The effect of suction on the stability of compressible flows over backward-facing steps is investigated. Mach numbers up to 0.8 are considered. The results show that continuous suction stabilizes the flow outside the separation bubble, but it destabilizes the flow inside it. Nevertheless, the overall N factor decreases as the suction level increases due to the considerable reduction of the separation bubble. For the same suction flow rate, properly distributed suction strips stabilize the flow more than continuous suction. The size of the separation bubble, and hence its effect on the instability can be considerably reduced by placing strips with high suction velocities in the separation region.

Al-Maaitah, Ayman A.↗

Three-dimensional numerical predictions of the flow behind a rearward-facing-step in a supersonic combustor

A CFD code was developed to compute the mixing and combustion of hydrogen fuel in the turbulent flow fields of scramjets. The code is based on the complete Reynold's time-averaged three-dimensional Navier-Stokes equations, and chemical species equations. The code includes a global finite rate hydrogen-air chemical reaction model and the Baldwin-Lomax (1978) turbulence model. The validity of this CFD code was demonstrated by comparing computed results from the code with existing data for three configurations, each containing an important feature of the scramjet combustor flow. The code was then applied to the reacting flow field in a complex combustor model in which hydrogen fuel was injected transversely into a supersonic hot airstream through several fuel injectors.

Uenishi, K.↗

Development of separation due to interaction between a shock wave and a turbulent boundary layer perturbed by rarefaction waves

The interaction of a turbulent boundary layer with rarefaction and shock waves in flows past inclined steps was investigated experimentally and theoretically. The experiments were carried out in supersonic wind tunnels with test sections of 0.6 x 0.6 m and 0.2 x 0.2 m, with adiabatic conditions on the surface model. Based on detailed measurements of pressure fields, velocity, and surface friction, combined with results of optical visualization and analysis of limiting flow lines, three characteristic flow regimes are identified. These are (1) nonseparated flow at small step angles, (2) formation of a local separation zone with a free separation point at moderate step angles, and (3) formation of a large-scale separated flow with a fixed separation point at sufficiently large step angles. The quantitative gasdynamic schemes and surface pressure distributions for these flow regimes are shown for Mach 2.85.

Zheltovodov, A. A.↗

Stereo Particle Image Velocimetry Measurements of Transition Downstream of a Forward-Facing Step in a Swept-Wing Boundary Layer

Stereo particle image velocimetry measurements were performed downstream of a forward-facing step in a stationary-crossflow dominated flow. Three different step heights were studied with the same leading-edge roughness configuration to determine the effect of the step on the evolution of the stationary-crossflow. Above the critical step height, which is approximately 68% of the boundary-layer thickness at the step, the step caused a significant increase in the growth of the stationary crossflow. For the largest step height studied (68%), premature transition occurred shortly downstream of the step. The stationary crossflow amplitude only reached approximately 7% of U(sub e) in this case, which suggests that transition does not occur via the high-frequency secondary instabilities typically associated with stationary crossflow transition. The next largest step of 60% delta still caused a significant impact on the growth of the stationary crossflow downstream of the step, but the amplitude eventually returned to that of the baseline case, and the transition front remained the same. The smallest step height (56%) only caused a small increase in the stationary crossflow amplitude and no change in the transition front. A final case was studied in which the roughness on the leading edge of the model was enhanced for the lowest step height case to determine the impact of the stationary crossflow amplitude on transition. The stationary crossflow amplitude was increased by approximately four times, which resulted in premature transition for this step height. However, some notable differences were observed in the behavior of the stationary crossflow mode, which indicate that the interaction mechanism which results in the increased growth of the stationary crossflow downstream of the step may be different in this case compared to the larger step heights.

Eppink, Jenna L.↗

Arsenic Removal by Continuous Flow-Through Iron Electrocoagulation and Downstream Steps under Environmentally Relevant Conditions

Arsenic (As) contamination in drinking water is a major health concern. This study investigated As(III) and As(V) removal by continuous flow-through iron electrocoagulation (EC) under varying iron doses and water chemistry conditions. Effective arsenic removal (>90% of 100 µg/L) was achieved with a low iron dose of 5 mg/L and a short EC reactor residence time of 11 s that was followed by flocculation and settling. Removal was primarily through adsorption of As(III) and As(V) to the EC-generated Fe(III) oxyhydroxide solids ferrihydrite and lepidocrocite. For influents containing As(III), oxidation of As(III) to As(V) also occurred. The removal processes occurred in both the EC reactor and downstream treatment. EC-generated reactive solids with high surface area and adsorption capacity provided better arsenic removal than with preformed adsorbents. The removed arsenic remained largely immobile in the EC-generated solids. The extent of arsenic removal was generally similar across pH 6 to 8. Phosphate (1 mg/L as P) inhibited the removal of both arsenic species. Dissolved silica (20 mg/L as SiO2) suppressed As(III) removal but had little effect on the amount of As(V) removed. In addition to thermodynamic effects on adsorption equilibria, pH and co-occurring phosphate and silica influenced removal rates by affecting iron oxidation, nucleation, and aggregation kinetics. The results integrate mechanistic understanding with treatment performance, providing insight relevant to both fundamental research and the application of EC in water treatment.

Arsenic, electrocoagulation, iron oxide, phosphate↗

Heat Transfer in a Liquid Metal Flowing Turbulently through a Channel with a Step Function Boundary Temperature

An analytical heat transfer solution is derived and evaluated for the general case of a turbulently flowing liquid metal which suddenly encounters a step-function boundary temperature in a channel system. Local Nusselt moduli, dimensionless mixed-mean fluid temperatures, and arithmetic-mean Nusselt moduli are given as functions of Reynolds and Prandtl moduli and a dimensionless axial-distance modulus. These solutions are compared with known solutions of more specific systems as well as with a set of experimental liquid-metal heat transfer data for a thermal entrance region.

Poppendiek, H. F.↗

Flight-measured base pressure coefficients for thick boundary-layer flow over an aft-facing step for Mach numbers from 0.4 to 2.5

A 0.56-inch thick aft-facing step was located 52.1 feet from the leading edge of the left wing of an XB-70 airplane. A boundary-layer rake at a mirror location on the right wing was used to obtain local flow properties. Reynolds numbers were near 10 to the 8th power, resulting in a relatively thick boundary-layer. The momentum thickness ranged from slightly thinner to slightly thicker than the step height. Surface static pressures forward of the step were obtained for Mach numbers near 0.9, 1.5, 2.0, and 2.4. The data were compared with thin boundary-layer results from flight and wind-tunnel experiments and semiempirical relationships. Significant differences were found between the thick and the thin boundary-layer data.

Goecke, S. A.↗

Toward particle-resolved accuracy in Euler–Lagrange simulations of multiphase flow using machine learning and pairwise interaction extended point-particle (PIEP) approximation

This study presents two different machine learning approaches for the modeling of hydrodynamic force on particles in a particle-laden multiphase flow. Results from particle-resolved direct numerical simulations (PR-DNS) of flow over a random array of stationary particles for eight combinations of particle Reynolds number (Re) and volume fraction (φ) are used in the development of the models. The first approach follows a two-step process. In the first flow prediction step, the perturbation flow due to a particle is obtained as an axisymmetric superposable wake using linear regression. In the second force prediction step, the force on a particle is evaluated in terms of the perturbation flow induced by all its neighbors using the generalized Faxén form of the force expression. In the second approach, the force data on all the particles from the PR-DNS simulations are used to develop an artificial neural network (ANN) model for direct prediction of force on a particle. Due to the unavoidable limitation on the number of fully resolved particles in the PR-DNS simulations, direct force prediction with the ANN model tends to over-fit the data and performs poorly in the prediction of test data. In contrast, due to the millions of grid points used in the PR-DNS simulations, accurate flow prediction is possible, which then allows accurate prediction of particle force. This hybridization of multiphase physics and machine learning is particularly important, since it blends the strength of each, and the resulting pairwise interaction extended point-particle model cannot be developed by either physics or machine learning alone.

42 ENGINEERING↗

Numerical analysis of laminar and turbulent incompressible flows using the finite element Fluid Dynamics Analysis Package (FIDAP)

The purpose of the study is the evaluation of the numerical accuracy of FIDAP (Fluid Dynamics Analysis Package). Accordingly, four test problems in laminar and turbulent incompressible flows are selected and the computational results of these problems compared with other numerical solutions and/or experimental data. These problems include: (1) 2-D laminar flow inside a wall-driven cavity; (2) 2-D laminar flow over a backward-facing step; (3) 2-D turbulent flow over a backward-facing step; and (4) 2-D turbulent flow through a turn-around duct.

Sohn, Jeong L.↗

Reaction Pathways over ZnZrO 2 -Based Catalysts and Catalytic Sorbents

Reactive capture and conversion (RCC) is a process intensification approach that integrates CO 2 capture and hydrogenation within a single unit, removing the CO 2 purification and storage steps of traditional process flow schemes. This alters the catalytic step from a traditional steady-state (SS) flow process to a transient capture and conversion cycle, which could lead to product distributions distinct from those observed in conventional SS experiments. Such differences are investigated in the combined capture and hydrogenation of carbon dioxide to methanol over a ZnZrO 2 catalyst and a ZnZrO 2 + NaNO 3 /Mg 3 AlO x catalytic sorbent (CS) using fixed-bed kinetic measurements, in situ diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS), and steady-state isotopic transient kinetic analysis-DRIFTS (SSITKA-DRIFTS). Under SS conditions, ZnZrO 2 produced methanol through sequential hydrogenation of HCOO* and CH 3 O* intermediates. On the contrary, CO was attributed primarily to CO 2 dissociation at oxygen vacancies, as supported by isotopic shifts and measured reaction orders. For the CS, isotopic switching experiments suggested that monodentate carbonate species (CO 3 2− , abbreviated as m-CO 3 2− ) act as active intermediates that can be hydrogenated to HCOO* and subsequently to CH 3 O. Under RCC conditions, in situ DRIFTS and isotopic experiments reveal that m-CO 3 2− species formed during the CO 2 capture step follow two competing routes upon H 2 exposure: (i) direct hydrogenation to methane on the sorbent domain or (ii) migration of m-CO 3 2− to the ZnZrO 2 domain, where they are hydrogenated to methanol through the HCOO pathway. Overall, RCC enables carbonate hydrogenation routes not observed under SS cofeed conditions. Thus, the reaction pathways and rates during RCC can be different from operation under conventional SS conditions, and the product distribution is determined here by competition between carbonate hydrogenation on sorbent sites and migration to ZnZrO 2 for methanol synthesis.

CCUS↗

Efficient and Flexible Sensitivity Matrix Computation for Adaptive Electrical Capacitance Volume Tomography

Electrical capacitance tomography is a widely used sensor modality for flow imaging in many industrial settings. Adaptive Electrical Capacitance Volume Tomography (AECVT) extends the capabilities of traditional ECT by enabling direct volumetric imaging and an improved resolution. Construction of the sensitivity matrix is a necessary step to obtain flow images. This step requires computation of the electric field inside the sensing domain, which is done via a typical field solver such as the finite element method. In this work, we present an efficient and flexible method to construct the sensitivity matrix for Adaptive Electrical Capacitance Volume Tomography (AECVT) based on individual electrode segment excitations and their judicious combination to form desired matrix elements. We illustrate how the proposed method yields the same sensitivity matrix as the traditional method but at a much lower computational cost. Once all segment contributions are obtained, we also indicate how the proposed method, unlike the traditional approach, can generate the sensitivity matrix on demand for an arbitrary combination of synthetic electrodes and obviating the need for any additional field computations. Finally, we present image reconstruction results for two different experimental scenarios where the mutual capacitance data and the corresponding sensitivity vectors are obtained through the proposed measurement combination scheme.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Method of Simulating Flow-Through Area of a Pressure Regulator

The flow-through area of a pressure regulator positioned in a branch of a simulated fluid flow network is generated. A target pressure is defined downstream of the pressure regulator. A projected flow-through area is generated as a non-linear function of (i) target pressure, (ii) flow-through area of the pressure regulator for a current time step and a previous time step, and (iii) pressure at the downstream location for the current time step and previous time step. A simulated flow-through area for the next time step is generated as a sum of (i) flow-through area for the current time step, and (ii) a difference between the projected flow-through area and the flow-through area for the current time step multiplied by a user-defined rate control parameter. These steps are repeated for a sequence of time steps until the pressure at the downstream location is approximately equal to the target pressure.

Hass, Neal E.↗

Effect of wall cooling on the stability of compressible subsonic flows over smooth humps and backward-facing steps

The effect of wall cooling on the two-dimensional linear stability of subsonic flows over two-dimensional surface imperfections is investigated. Results are presented for flows over smooth humps and backward-facing steps with Mach numbers up to 0.8. The results show that, whereas cooling decreases the viscous instability, it increases the shear-layer instability and hence it increases the growth rates in the separation region. The coexistence of more than one instability mechanism makes a certain degree of wall cooling most effective. For the Mach numbers 0.5 and 0.8, the optimum wall temperatures are about 80 pct and 60 pct of the adiabatic wall temperature, respectively. Increasing the Mach number decreases the effectiveness of cooling slightly and reduces the optimum wall temperature.

Al-Maaitah, Ayman A.↗

Effect of wall cooling on the stability of compressible subsonic flows over smooth humps and backward-facing steps

The effect of wall cooling on the two-dimensional linear stability of subsonic flows over two-dimensional surface imperfections is investigated. Results are presented for flows over smooth humps and backward-facing steps with Mach numbers up to 0.8. The results show that, whereas cooling decreases the viscous instability, it increases the shear-layer instability and hence it increases the growth rates in the separation region. The coexistence of more than one instability mechanism makes a certain degree of wall cooling most effective. For the Mach numbers 0.5 and 0.8, the optimum wall temperatures are about 80 pct and 60 pct of the adiabatic wall temperature, respectively. Increasing the Mach number decreases the effectiveness of cooling slightly and reduces the optimum wall temperature.

Al-Maaitah, Ayman A.↗