Numerical Solutions to Three-Dimensional Flow Problems

Author(s):  
Bastian E. Rapp
2001 ◽  
Vol 432 ◽  
pp. 127-166 ◽  
Author(s):  
K. W. BRINCKMAN ◽  
J. D. A. WALKER

Unsteady separation processes at large finite, Reynolds number, Re, are considered, as well as the possible relation to existing descriptions of boundary-layer separation in the limit Re → ∞. The model problem is a fundamental vortex-driven three-dimensional flow, believed to be relevant to bursting near the wall in a turbulent boundary layer. Bursting is known to be associated with streamwise vortex motion, but the vortex/wall interactions that drive the near-wall flow toward breakdown have not yet been fully identified. Here, a simulation of symmetric counter-rotating vortices is used to assess the influence of sustained pumping action on the development of a viscous wall layer. The calculated solutions describe a three-dimensional flow at finite Re that is independent of the streamwise coordinate and consists of a crossflow plane motion, with a developing streamwise flow. The unsteady problem is constructed to mimic a typical cycle in turbulent wall layers and numerical solutions are obtained over a range of Re. Recirculating eddies develop rapidly in the near-wall flow, but these eddies are eventually bisected by alleyways which open up from the external flow region to the wall. At sufficiently high Re, an oscillation was found to develop in the streamwise vorticity field near the alleyways with a concurrent evolution of a local spiky behaviour in the wall shear. Above a critical value of Re, the oscillation grows rapidly in amplitude and eventually penetrates the external flow field, suggesting the onset of an unstable wall-layer breakdown. Local zones of severely retarded streamwise velocity are computed which are reminiscent of the low-speed streaks commonly observed in turbulent boundary layers. A number of other features also bear a resemblance to observed coherent structure in the turbulent wall layer.


2014 ◽  
Vol 748 ◽  
pp. 433-456 ◽  
Author(s):  
Giuliano De Stefano ◽  
Oleg V. Vasilyev

AbstractThe wavelet-based eddy capturing approach is extended to three-dimensional bluff body flows, where the flow geometry is enforced through Brinkman volume penalization. The wavelet-collocation/volume-penalization combined method is applied to the simulation of vortex shedding flow behind an isolated stationary prism with square cross-section. Wavelet-based direct numerical simulation is conducted at low supercritical Reynolds number, where the wake develops fundamental three-dimensional flow structures, while wavelet-based adaptive large-eddy simulation supplied with the one-equation localized dynamic kinetic-energy-based model is performed at moderately high Reynolds number. The present results are in general agreement with experimental findings and numerical solutions provided by classical non-adaptive methods. This study demonstrates that the proposed hybrid methodology for modelling bluff body flows is feasible, accurate and efficient.


1968 ◽  
Vol 8 (04) ◽  
pp. 331-340 ◽  
Author(s):  
K.H. Coats

Coats, K.H., Member AIME, The U. of Texas, Austin, Texas Abstract This paper describes a generalized analysis for calculating three-phase, three-dimensional flow in reservoirs. The analysis handles pressure maintenance type problems where fluid compressibility effects are negligible. A separate analysis for depletion type problems is described in another paper. The calculations consist of numerical, simultaneous solution of the three-flow equations using the iterative alternating direction technique of Douglas and Rachford. The mathematical details are fully described in the Appendix. The analysis is a computerized mathematical model that accounts for gravity, and capillary and viscous forces, and allows arbitrary reservoir heterogeneity, geometry, well locations and rates. A unique aspect of the analysis is the simultaneous solution of only as many difference equations in each grid block of the reservoir as there are mobile phases present. Thus, while the analysis handles phases present. Thus, while the analysis handles three-phase flow, the efficiency of the calculations (in a typical problem where three phases actually coexist only in a minor portion of the reservoir) is four to eight times greater than that of an analysis solving three equations in every block. The program may be applied to two-phase flow problems and to one-, two- or three-dimensional flow problems with negligible loss in efficiency, compared to programs specifically written for these sub cases. This paper also describes several applications of the analysis which illustrate some effects of gravitational and capillary forces in waterflooding of a heterogeneous reservoir. Another application indicates the utility of the program in simulating the fillup stage of water injection into a reservoir containing an initial free gas phase. Computer times and costs for the applications performed are given to indicate the current expense performed are given to indicate the current expense of three-dimensional, three-phase reservoir simulation. Introduction Under pressure maintenance by water and/or gas injection, fluid compressibility effects are generally negligible in producing operations. Although gas compressibility may be appreciable, the maintenance of pressure results in negligible time variation of gas density. In addition, the spatial variation of gas density is usually small in relation to the gas density itself. Producing schemes of pattern or flank waterflood and/or crestal gas injection, therefore, may be simulated with an analysis which presumes fluid incompressibility. presumes fluid incompressibility. The computing efficiency of a numerical model for simulating incompressible fluid flow is as much as 50 percent greater than that of a compressible flow model. Therefore, an analysis for numerically simulating three-dimensional flow of three incompressible, immiscible phases was developed and programmed. programmed. THE MODEL The equations describing three-phase, incompressible flow are the continuity equation and Darcy's law for each flowing phase. Combining these equations and introducing capillary pressures gives the three flow equations: ................ (1a) ............... (1b) SPEJ p. 331


1981 ◽  
Vol 109 ◽  
pp. 259-275 ◽  
Author(s):  
G. D. Mallinson ◽  
A. D. Graham ◽  
G. De Vahl Davis

A numerical and experimental study has been made of the three-dimensional flow and heat transfer by natural convection in a closed, rectangular thermosyphon. At low Rayleigh numbers, the flows in the two halves of the cavity remain separate, with heat transfer across the mid-height plane occurring only by conduction. At increasing Rayleigh numbers, an exchange process of increasing complexity occurs. The numerical solutions were used to explore this process and to predict flow patterns which were found to resemble closely those observed during previous investigations of cylindrical thermosyphons. The results were verified by a flow visualization study. Heat-transfer rates are presented and augment previous data for higher values of the governing parameters.


1995 ◽  
Vol 117 (2) ◽  
pp. 208-218 ◽  
Author(s):  
C. J. Freitas

This paper summarizes the results of a series of five benchmark simulations which were completed using commercial Computational Fluid Dynamics (CFD) codes. These simulations were performed by the vendors themselves, and then reported by them in ASME’s CFD Triathlon Forum and CFD Biathlon Forum. The first group of benchmarks consisted of three laminar flow problems. These were the steady, two-dimensional flow over a backward-facing step, the low Reynolds number flow around a circular cylinder, and the unsteady three-dimensional flow in a shear-driven cubical cavity. The second group of benchmarks consisted of two turbulent flow problems. These were the two-dimensional flow around a square cylinder with periodic separated flow phenomena, and the steady, three-dimensional flow in a 180-degree square bend. All simulation results were evaluated against existing experimental data and thereby satisfied item 10 of the Journal’s policy statement for numerical accuracy. The objective of this exercise was to provide the engineering and scientific community with a common reference point for the evaluation of commercial CFD codes.


1985 ◽  
Vol 107 (2) ◽  
pp. 277-284 ◽  
Author(s):  
Qinghuan Wang ◽  
Genxing Zhu ◽  
Chung-Hua Wu

Progress in the development of quasi-three-dimensional and full three-dimensional numerical solutions for steady subsonic rotational flow through turbomachines is presented. An iterative calculation between the flow on a mean hub-to-tip S2 stream surface and a number of blade-to-blade S1 stream surfaces gives the quasi-three-dimensional solution, which is very easily extended to give full three-dimensional solution by merely calculating a few more S2 surface flows and relaxing the restriction that S1 surfaces are surfaces of revolution. A new S2–S1 iteration scheme has been developed and employed in the present code. The governing equations on the S1 and S2 surfaces are expressed in terms of general nonorthogonal curvilinear coordinates so that they are body-fitted without any coordinate transformation and are solved by either matrix method or line-relaxation method. An automatic computing system is used, which first computes the quasi-three-dimensional flow for blade design and then computes the full three-dimensional flow for the blade row just designed. The results obtained by applying this computing system to the design and determination of full three-dimensional flow field of a two-stage axial compressor and a high subsonic compressor stator are obtained and shows clearly the amount of the twist of the general S1 surfaces and the difference in the flow field between the quasi-three-dimensional and full three-dimensional solutions.


Sign in / Sign up

Export Citation Format

Share Document