A numerical simulation of the transition to turbulence in a two-dimensional flow

A numerical simulation of the transition to turbulence is performed using a finite element method. The unsteady Navier-Stokes equations are discretized using a standard Galerkin approximation and a loading strategy for increasing the Reynolds number. The numerical results are then analysed at different Reynolds numbers showing a transition from a steady-state solution to a weakly chaotic one. Phase space diagrams are presented showing the presence of strange attractors. The dimension and Lyapunov's exponents of these attractors are computed and compared with existing results in the literature.     

Experimental and numerical flow visualization in a transonic turbine

This paper focuses on the visualization of both experimental and numerical results and presents research in a highly-loaded cold-flow transonic turbine under continuous and engine-representative conditions. Special focus was placed on blade row interaction at app. 10600 rpm. While the first step was the investigation of a single stage machine (stator-rotor), the second step extended the test rig to a one-and-a-half stage configuration (stator-rotor-stator). Measurements were carried out in the transonic test turbine at Graz University of Technology using Particle-Image-Velocimetry and Laser-Doppler-Anemometry. The main results of these experiments are discussed and compared to numerical simulations.     

Conservative numerical simulation of multi-component transport in two-dimensional unsteady shallow water flow

An explicit finite volume model to simulate two-dimensional shallow water flow with multi-component transport is presented. The governing system of coupled conservation laws demands numerical techniques to avoid unrealistic values of the transported scalars that cannot be avoided by decreasing the size of the time step. The presence of non conservative products such as bed slope and friction terms, and other source terms like diffusion and reaction, can make necessary the reduction of the time step given by the Courant number. A suitable flux difference redistribution that prevents instability and ensures conservation at all times is used to deal with the non-conservative terms and becomes necessary in cases of transient boundaries over dry bed. The resulting method belongs to the category of well-balanced Roe schemes and is able to handle steady cases with flow in motion. Test cases with exact solution, including transient boundaries, bed slope, friction, and reaction terms are used to validate the numerical scheme. Laboratory experiments are used to validate the techniques when dealing with complex systems as the κ–?$\kappa ''–''?$ model. The results of the proposed numerical schemes are compared with the ones obtained when using uncoupled formulations.     

Experimental and numerical dye washout flow visualization

Flow visualization in realistic models is very important for the study of pathological vessel enlargements (aneurysms). Furthermore, flow visualization may help in treatment decisions. However, the most interesting parameter, the wall shear stress, is difficult to measure in vivo. This parameter can be provided by computational fluid dynamics. However, the numerical methods don’t visualize the results as does of the dye washout method — a method often used in flow studies. This experimental method simulates the cine angiograms acquired during contrast agent injection used in medicine. In this paper we present the dye washout visualization of CFD results and compare these results with the conventional dye washout experiments in the same aneurysm model under steady flow conditions.     

Microscopic-based fluid flow simulation of invasion on a two-dimensional lattice. II. Mobilization and cohesion

An algorithm for modeling secondary invasion processes in porous media is presented. Mobilization of trapped defender fluid is accomplished through interfacial interaction rules. Cohesive forces are also included within the defender phase. A series of simulation runs are performed using two-dimensional lattices and examined to determine optimal conditions for secondary invasions that sweep the trapped defender phase from the porous medium.     

Turbulent flame visualization using direct numerical simulation

Combustion phenomena are of high scientific and technological interest, in particular for energy generation and transportation systems. Direct Numerical Simulations (DNS) have become an essential and well established research tool to investigate the structure of turbulent flames, since they do not rely on any approximate turbulence models. In this work two complementary DNS codes are employed to investigate different types of fuels and flame configurations. The code is π³ is a 3-dimensional DNS code using a low-Mach number approximation. Chemistry is described through a tabulation, using two coordinates to enter a database constructed for example with 29 species and 141 reactions for methane combustion. It is used here to investigate the growth of a turbulent premixed flame in a methane-air mixture (Case 1). The second code,Sider is an explicit three-dimensional DNS code solving the fully compressible reactive Navier-Stokes equations, where the chemical processes are computed using a complete reaction scheme, taking into account accurate diffusion properties. It is used here to compute a hydrogen/air turbulent diffusion flame (Case 2), considering 9 chemical species and 38 chemical reactions.     

Flow visualization of a longitudinal vortex in drag-reducing surfactant flow

Vortical structures behind a wing-type vortex generator (WVG) immersed in a drag-reducing flow were visualized with a fluorescence dye illuminated with a planer laser sheet to investigate the enhancement mechanism of heat transfer on the WVG-installed target wall. Particle image velocimetry (PIV) system was also used to measure the flow velocity profiles. Both the flow visualization and the velocity measurement revealed that small amounts of surfactant, Cetyl-Trimethyl-Ammonium-Chloride/ Sodium Salicylate (CTAC/NaSal), added in water changed drastically the flow patterns together with the wall heat transfer distributions due to the reduction of momentum exchange in the vortical flows behind WVG.     

Transition to chaos in a confined two-dimensional fluid flow

For a two-dimensional fluid in a square domain with no-slip walls, new direct numerical simulations reveal that the transition from steady to chaotic flow occurs through a sequence of various periodic and quasiperiodic flows, similar to the well-known Ruelle-Takens-Newhouse scenario. For all solutions beyond the ground state, the phenomenology is dominated by a domain-filling circulation cell, whereas the associated symmetry is reduced from the full symmetry group of the square to rotational symmetry over an angle pi. The results complement both laboratory experiments in containers with rigid walls and numerical simulations on double-periodic domains.     

3D simulation and visualization of unsteady wake flow behind a cylinder

In this paper we focus on 3D simulation of unsteady wake flow behind a circular cylinder. We show that in addition to accurate formulations and sufficiently-refined meshes, efficient computing methods are essential components of an effective simulation strategy. We use the Multi-Domain Method (MDM) we developed recently in computation of two cases. At Reynolds number 300, we demonstrate how the MDM enables us to use highly-refined meshes to capture wake patterns which we otherwise cannot fully represent. At Reynolds number 140, we show that with the MDM we can extend our computations sufficiently downstream, and with sufficient accuracy, to successfully capture the second phase of the Karman vortex street, which has been observed in laboratory experiments, and which has double the spacing between the vortices compared to the first phase.     

Lattice—BGK simulation of a two—dimensional channel flow around a square cylinder

The confined flow around a square cylinder mounted inside a two-dimensional channel (blockage ratio $\be=1/8$) was investigated in detail by a newly developed incompressible nonuniform lattice-BGK model. It is found that the vortex shedding behind the cylinder induces periodicity in the flow field, and the periodicity of the flow will lose for $Re>$300. A detailed analysis for a range of Reynolds numbers between 1 and 500 was presented. Quantitative comparisons with other methods show that the model gives accurate results for complex flows.     

Three-dimensional direct numerical simulation of turbulent channel flow catalytic combustion of hydrogen over platinum

The turbulent catalytic combustion of a fuel-lean hydrogen/air mixture (equivalence ratio ? = 0.24) was investigated by means of three-dimensional direct numerical simulation (DNS) in a platinum-coated plane channel with a prescribed wall temperature of 960 K and an incoming Reynolds number, based on the channel height, of 5700. Heat transfer from the hot catalytic walls laminarized the flow, as manifested by the progressive suppression of the high vorticity components of the flow aligned parallel to the channel walls at increasing streamwise distances. The impact of turbulence suppression on the mass transfer towards or away from the catalytic wall was subsequently assessed. Far upstream where high turbulence fluctuations persisted, the instantaneous local transverse gradient of the limiting hydrogen reactant (a quantity proportional to the catalytic reaction rate) as well as the instantaneous hydrogen concentration at the wall exhibited strong fluctuations by up to 300%, a result of finite-rate chemistry induced by the high inrush events towards the catalytic walls. Fourier analysis of the reaction rate fluctuations yielded peak frequencies of less than 1 kHz, values comparable to the thermal response frequencies of typical materials in commercial catalytic geometries. This has direct implications on the thermal stress of the reactor walls as well as on the decoupling between flow and solid thermal modeling currently used in practical catalytic reactors. Far downstream, the dampening of turbulence resulted in weaker hydrogen concentration fluctuations with nearly symmetric distributions. Finally, computed transverse turbulent species fluxes indicated inherent weaknesses of near-wall turbulence models in describing turbulent transport of species with disparate molecular diffusivities.     

Particle-based simulation and visualization of fluid flows through porous media

Abstract  

We propose a method of fluid simulation where boundary conditions are designed in such a way that fluid flow through porous media, pipes, and chokes can be realistically simulated. Such flows are known to be low Reynolds number incompressible flows and occur in many real life situations. To obtain a high quality fluid surface, we include a scalar value in isofunction. The scalar value indicates the relative position of each particle with respect to the fluid surface.     

辐射烧蚀CH薄膜非平衡理论计算与实验比对分析

用一维非平衡辐射输运程序和实验提供的黑腔辐射流,对CH薄膜靶的辐射烧蚀进行了数值模拟计算。计算结果表明,非平衡的黑腔辐射场对CH薄膜靶辐射烧蚀的影响十分明显。出靶外界面的辐射流总强度、等效辐射温度随时间变化、400eV光子出界面延迟时间随靶厚度变化规律等方面的数值模拟结果与实验结果基本符合。     

辐射烧蚀CH薄膜非平衡理论计算与实验比对分析

 用一维非平衡辐射输运程序和实验提供的黑腔辐射流，对CH薄膜靶的辐射烧蚀进行了数值模拟计算。计算结果表明，非平衡的黑腔辐射场对CH薄膜靶辐射烧蚀的影响十分明显。出靶外界面的辐射流总强度、等效辐射温度随时间变化、400eV光子出界面延迟时间随靶厚度变化规律等方面的数值模拟结果与实验结果基本符合。     

Vorticity dynamics and sound generation in two-dimensional fluid flow

An approximate solution to the two-dimensional incompressible fluid equations is constructed by expanding the vorticity field in a series of derivatives of a Gaussian vortex. The expansion is used to analyze the motion of a corotating Gaussian vortex pair, and the spatial rotation frequency of the vortex pair is derived directly from the fluid vorticity equation. The resulting rotation frequency includes the effects of finite vortex core size and viscosity and reduces, in the appropriate limit, to the rotation frequency of the Kirchhoff point vortex theory. The expansion is then used in the low Mach number Lighthill equation to derive the far-field acoustic pressure generated by the Gaussian vortex pair. This pressure amplitude is compared with that of a previous fully numerical simulation in which the Reynolds number is large and the vortex core size is significant compared to the vortex separation. The present analytic result for the far-field acoustic pressure is shown to be substantially more accurate than previous theoretical predictions. The given example suggests that the vorticity expansion is a useful tool for the prediction of sound generated by a general distributed vorticity field.