• Korea-Australia Rheology Journal
• /
• v.20 no.3
• /
• pp.117-125
• /
• 2008

Viscoelastic instabilities are of fundamental importance to understanding the physics of complex fluids and of practical importance to materials processing and fluid characterization. Significant progress has been made over the past 15 years in understanding instabilities in viscoelastic flows with curved streamlines and is reviewed here. Taylor-Couette flow, torsional flow between a cone and plate, and torsional flow between parallel plates have received special attention due to both the basic significance of these flows and their critical role in rheometry. First, we review the criteria for determining when these flows become unstable due to elasticity in the absence of inertia, and discuss the generalization of these criteria to more complex flows with curved streamlines. Then, focusing on experiments and simulations in the Taylor-Couette problem, we review how thermal sensitivity (i.e., the dependence of fluid viscosity and elasticity on temperature) and inertia affect the stability of viscoelastic flows. Finally, we conclude with some general thoughts on unresolved issues and remaining challenges related to viscoelastic instabilities.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.86-87
• /
• 2003

Unstructured grid system is suitable for flows of complex geometries. For problems with moving boundary walls, the grid system must be changed and deformed with time if we use a body fitted grid system. In this paper, a new moving-grid finite-volume method on unstructured grid system is proposed and developed for unsteady compressible flows with shock waves. To assure geometric conservation laws on moving grid system, a control volume on the space-time unified domain is adopted for estimating numerical flux. The method is described and applied for two-dimensional flows.

• Proceedings of the Korean Society of Computational and Applied Mathematics Conference
• /
• 2003.09a
• /
• pp.1-1
• /
• 2003

In this talk, a reduced-order modeling methodology based on centroidal Voronoi tessellations (CVT's)is introduced. CVT's are special Voronoi tessellations for which the generators of the Voronoi diagram are also the centers of mass (means) of the corresponding Voronoi cells. The discrete data sets, CVT's are closely related to the h-means clustering techniques. Even with the use of good mesh generators, discretization schemes, and solution algorithms, the computational simulation of complex, turbulent, or chaotic systems still remains a formidable endeavor. For example, typical finite element codes may require many thousands of degrees of freedom for the accurate simulation of fluid flows. The situation is even worse for optimization problems for which multiple solutions of the complex state system are usually required or in feedback control problems for which real-time solutions of the complex state system are needed. There hava been many studies devoted to the development, testing, and use of reduced-order models for complex systems such as unsteady fluid flows. The types of reduced-ordered models that we study are those attempt to determine accurate approximate solutions of a complex system using very few degrees of freedom. To do so, such models have to use basis functions that are in some way intimately connected to the problem being approximated. Once a very low-dimensional reduced basis has been determined, one can employ it to solve the complex system by applying, e.g., a Galerkin method. In general, reduced bases are globally supported so that the discrete systems are dense; however, if the reduced basis is of very low dimension, one does not care about the lack of sparsity in the discrete system. A discussion of reduced-ordering modeling for complex systems such as fluid flows is given to provide a context for the application of reduced-order bases. Then, detailed descriptions of CVT-based reduced-order bases and how they can be constructed of complex systems are given. Subsequently, some concrete incompressible flow examples are used to illustrate the construction and use of CVT-based reduced-order bases. The CVT-based reduced-order modeling methodology is shown to be effective for these examples and is also shown to be inexpensive to apply compared to other reduced-order methods.

• Journal of Mechanical Science and Technology
• /
• v.15 no.12
• /
• pp.1750-1756
• /
• 2001

Three-dimensional and unsteady flow analysis is a practical target of high performance computation. As recently advances of computers, a numerical prediction by the large eddy simulation (LES) are introduced and evaluated for various engineering problems. Its advanced methods for the complex turbulent flows are discussed by several examples applied for aerodynamic designs, analysis of fluid flow mechanisms and their interaction to complex phenomena. These results of time-dependent and three-dimensional phenomena are visualized by interactive graphics and animations.

• Korea-Australia Rheology Journal
• /
• v.21 no.4
• /
• pp.225-233
• /
• 2009

In various attempts to relate the behaviour of highly-elastic liquids in complex flows to their rheometrical behaviour, obvious candidates for study have been the variation of shear viscosity with shear rate, the two normal stress differences $N_1$ and $N_2$, especially $N_1$, and the extensional viscosity $\eta_E$. In this paper, we shall be mainly interested in 'constant-viscosity' Boger fluids, and, accordingly, we shall limit attention to $N_1$ and $\eta_E$. We shall concentrate on two important flows - axisymmetric contraction flow and "splashing" (particularly that which arises when a liquid drop falls onto the tree surface of the same liquid). Modern numerical techniques are employed to provide the theoretical predictions. It is shown that the two obvious manifestations of viscoelastic rheometrical behaviour can sometimes be opposing influences in determining flow characteristics. Specifically, in an axisymmetric contraction flow, high $\eta_E$ can retard the flow, whereas high $N_1$ can have the opposite effect. In the splashing experiment, high $\eta_E$ can certainly reduce the height of the so-called Worthington jet, thus confirming some early suggestions, but, again, other rheometrical influences can also have a role to play and the overall picture may not be as clear as it was once envisaged.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.2
• /
• pp.668-676
• /
• 1991

Turbulent flows around tube banks and in the diffuser were studied using a non-orthogonal boundary fitted coordinate system and the modified K-.epsilon. turbulence model. In these cases, many problems emerge which stem from the geometrical complexity of the flow domain and the physical complexity of turbulent flow itself. To treat the complex geometry, governing equations were reformulated in a non-orthogonal coordinate system with Cartesian velocity components and discretised by the finite volume method with a non-staggered variable arrangement. The modified K-.epsilon. model of Hanjalic and Launer was applied to solve above two cases under the condition of strong and mild pressure gradient. The results using the modified K-.epsilon. model results in both test cases.

• Journal of Advanced Marine Engineering and Technology
• /
• v.23 no.4
• /
• pp.533-542
• /
• 1999

In this study the flow characteristics of developing turbulent pulsating flows in a square-sec-tional 180。 curved duct are investigated experimentally. The experimental study of air flow in a square-sectional curved duct is carried out to measure axial velocity distribution secondary flow velocity profiles and wall shear stress distributions by using a Laser Doppler Velocimetry system with the data acquisition and processing system of Rotating Machinery Resolver (RMR) and PHASE software at the entrance region of the duct which is divided into 7 sections from the inlet(${{\o}}=0_{\circ}$) to the outlet (${{\o}}=180_{\circ}$) in $30_{\circ}$ intervals. The results obtained from the study are summarized as follows: (1) The time-averaged critical Dean number of turbulent pulsating flow(De ta, cr) is greater than $75{\omega}+$ It is understood that the critical Dean number and the critical Reynolds number are related to the dimensionless angular frequency in a curved duct. (2) Axial velocity profiles of turbulent pulsating flows are of an annular type similar to those of turbulent stead flows. (3) Secondary flows of trubulent pulsating flows are strong and complex at the entrance region. As velocity amplitudes(A1) become larger secondary flows become stronger. (4) Wall shear stress distributions of turbulent pulsating flows in a square-sectional $180_{\circ}$ curved duct are exposed variously in the outer wall and are stabilized in the inner wall without regard to the phase angle.

• Journal of Power System Engineering
• /
• v.2 no.1
• /
• pp.31-38
• /
• 1998

A diffuser, an important equipment to change kinetic energy into pressure energy, has been studied for a long time. Though experimental and theoretical researches have been done, the understanding of energy transfer and detailed mechanism of energy dissipation is unclear. As far as numerical prediction of diffuser flows are concerned, various numerical studies have also been done. On the contrary, many turbulence models have constraint to the applicability of diffuser-like complex flows, because of anisotropy of turbulence near the wall and of local nonequilibrium induced by an adverse pressure gradient. The existing $k-{\epsilon}$ turbulence models have some problems in the case of being applied to complex turbulent flows. The purpose of this paper is to test the applicability of the nonlinear $k-{\epsilon}$ model concerning diffuser-like flows with expansion and streamline curvature. The results show that the nonlinear $k-{\epsilon}$ turbulence model predicted well the coefficient of pressure, velocity profiles and turbulent kinetic energy distributions, however the shear stress prediction was failed.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2008.04a
• /
• pp.222-227
• /
• 2008

Complex terrain like hill, mountain, and escarpment etc. makes complex air flow. This topographic condition will affect not only speed but also turbulence of wind over the complex terrain. In this paper, turbulence intensities are considered to investigate characteristics of wind over cone-type hills. There are five simple hill models with different slope 0.1${\sim}$0.5(tan${\theta}$) for wind tunnel test. It was observed through wind tunnel tests that turbulence intensities of down-slope wind at the end of the 3-Dimensional hills remarkably increased but ones of windward slope wind at the front side of the hills slightly increased. Also, turbulence intensities proportionally increased with slope of the cone-type hills.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1996.05b
• /
• pp.173-178
• /
• 1996

A three-dimensional thermo-hydraulic computer code is developed for simulation of incompressible flows in complex geometries. The computer code employs a body-fitted, nonorthogonal grid system in order to efficiently handle the complex geometries encountered in many engineering applications. The finite volume method is used to discretize the governing equations and the convection term is treated by higher-order bounded schemes. The cell-centered, nonstaggered grid arrangement is adopted and the resulting checkerboard pressure oscillation is avoided by use of momentum interpolation practice. The computer code employs the SIMPLE algorithm for pressure and velocity coupling and the k-$\varepsilon$ turbulence for turbulent calculation. The computer code has been tested through application to a variety of test problems and some results are presented in this paper