2.1 quản phương trìnhPhương trình quản lý dòng chảy chất lỏng là phương trình toán học mà được gọi là bảo tồn các định luật vật lý. Các phương trình mô tả hành vi dòng chảy ở cấp độ vi, tức là vận tốc, nhiệt độ, mật độ và áp lực. Bằng cách sử dụng các phương trình, ba loại phương trình được xác định; phương trình khối lượng bảo tồn, động lực và năng lượng. Trong ba chiều, các yếu tố nhỏ nhất của các chất lỏng được giới thiệu như là một điểm hoặc hạt ở cấp độ vi mô. Hình 7 minh họa các yếu tố chất lỏng có sáu điều kiện biên. Các phương trình toán học được đại diện theo khối lượng, động lượng và năng lượng thay đổi qua các ranh giới (khuôn mặt). Cũng một điểm () được gọi là trung tâm của nguyên tố. Hình 7: Các yếu tố chất lỏng trong bảo tồn các định luật vật lý2.1.1 khối lượng bảo tồnSự cân bằng khối lượng của nguyên tố sẽ trở thành như các tỷ lệ ngày càng tăng của các khối lượng trong các yếu tố chất lỏng là tương đương với net lệ khối lượng dòng chảy vào các yếu tố. Do đó bảo tồn hàng loạt phương trình toán học trở thành(2,1)Trong phương trình trên, khối lượng đến để các yếu tố cho các dấu hiệu tích cực và dòng chảyvới dấu âm ra đến khối lượng từ ranh giới (xem hình 8). Hình 8: Khối lượng dòng chảy của mỗi khuôn mặt của nguyên tố chất lỏng, (hình ảnh nguồn: Versteeg và Malalasekera 2007)Phương trình nhỏ gọn hơn cũng có thể được viết dưới dạng (2,2)Vì vậy, phương trình (2,2) là phương trình ba chiều và không ổn định khối lượng bảo tồn (liên tục phương trình) tại một thời điểm () trong một chất lỏng nén. Đối với chất lỏng không nén, thay đổi mật độ trong thời gian là zero (mật độ là hằng số), do đó khối lượng bảo tồn phương trình sản lượng (2,3)2.1.2 Đà phương trìnhTheo định luật 2 Newton, vận động lượng thay đổi của một hạt chất lỏng là tương đương với tổng của các lực lượng trên khuôn mặt hạt.Thay đổi thuộc tính của các hạt chất lỏng được đại diện bởi luật pháp bảo tồn động lượng và năng lượng. Trong từ khác, bởi tên gọi của là một giá trị của một tài sản cho một đơn vị khối lượng, phát sinh từ tất cả là (2,4) Do đó, bằng cách thay thế tên gọi, và là một tài sản, mức độ ngày càng tăng của-, - và - đà cho mỗi đơn vị khối lượng của một chất lỏng được viết bằng Hai loại lực lượng tác động về hạt chất lỏng; lực lượng bề mặt và cơ thể lực lượng. Các lực lượng bề mặt là áp lực và độ nhớt cao lực lượng trong khi lực lượng cơ thể bao gồm trọng lực ly tâm,Coriolis và lực lượng điện từ. Ở đây lực lượng bề mặt được coi là một thuật ngữ riêng biệt trong phương trình động lực và cơ thể lực lượng là một thuật ngữ nguồn. Vì vậy lực lượng bề mặt làmột căng thẳng bình thường áp lực và chín thành phần nhớt căng thẳng trên các yếu tố chất lỏng.Bởi tổng kết của tất cả các lực lượng net tại và hướng, các lực lượng tất cả cho một đơn vị khối lượng được cho bởi (2,5) Bằng cách sử dụng phương trình (2,4) và (2,5), x-thành phần của phương trình động lực sẽ trở thành (2.6a) The y-component and z-component of the momentum equations also yield respectively(2.6b) (2.6c) where are the nine viscous stress components of the fluid element.2.1.3 Energy equationFirst law of thermodynamics states that rate of energy change of a fluid particle is equal to net rate of heat added to fluid particle plus net rate of work which is done on fluid particle. According to equation (2.4), the rate of energy increasing of a fluid particle is given by On the other side, the total rate of work that is done on the fluid particle by surface stresses after collecting the pressure terms together becomes Also the rate of heat added to the fluid particle from all boundary condition faces yields (2.8)By defining a source term of as potential energy changes and according to energybalance of fluid particle, the energy equation is given by (2.9)Where energy of a fluid ( ) is the sum of internal (thermal) energy, kinetic energy and gravitational potential energy.2.1.4 Navier-stokes equationThe fluid motion is described by solving the five governing equations; mass conservation, and momentum equations and energy equation. Also four thermodynamic variables of pressure ( ), Density ( ), internal energy ( ) and temperature ( ) can be solved by using the equation of states. For example if we use and as state variables, then the relation between these variables becomes and (2.10)But governing equations include unknown’s viscous stress components according to equations (2.6a-c). These viscous stress components are proportional to deformation rate in a Newtonian fluid. Therefore, by substitution of above viscous stresses into the momentum equations, theNavier-stokes equations are given by Which by re-arrange the shear stress terms, the most useful form of Navier-stokesequations are written,(2.13a)(2.13b)(2.13c)2.2 TurbulenceIn our surrounding, fluid flows are divided in two different regimes; Laminar and turbulent flows. As a good measure, Reynolds number states the critical boundary which at values below that, the fluid flow is laminar with smooth and stable pattern. By increasing the Reynolds number, the irregularity and instability is taken place and more complicated events is occurred and laminar flow is changed to turbulent flow.Most of the flows in our surrounding are turbulent flows which has irregular, random and unstable behaviour with unpredictable eddies (Figure 9). Some example of turbulence flows are outgoing plume of solid rocket jets, smoke from a chimney, river flow, etc. For the first time, von Karman (1937) presented basic definition of turbulence;“Turbulence is an irregular motion which in general makes its appearance in fluids, gaseous or liquid, when they flow past solid surfaces or even when neighbouring streams of the same fluid flow past or over one another” Figure 9: Image of a turbulent jet, (Picture source: Dahm and Dimotakis 1990)As the Figure 10 shows, the flow velocity varies irregularly both in time and position which is denoted by . In this Figure, the horizontal line is the mean axial velocity of centre line of a turbulent single jet and the flow velocity fluctuates in the time variation compare to mean velocity. Therefore the flow velocity equation is given by (2.14)where is a fluctuation velocity component with properties of .35343332313029282726250 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1time (s)Figure 10: Axial velocity component ( on centre line of a turbulent jet, Re=17232 Particle of fluid follows the eddy motions in turbulent flow; mass, momentum and heat are changed effectively. These eddy motions which can be seen in flow visualization, made by rotational flow structures with a wide range of length scales. The largest turbulent eddies interact with the energy that is extracted from mean flow by process of vortex stretching. Smaller eddies are stretched by larger eddies strongly rather than mean flow and it leads to transfer kinetic energy from larger eddies to smaller eddies which this process called energy cascade. The Figure 11 illustrates the kinetic energy handing down from large eddies to small eddies. The characteristic length and velocity of eddies are of the same order of mean flow length scale and velocity . Therefore Reynolds number of large eddy is defined as
(2.15)
Figure 11: Energy cascade illustration, (Picture source: Rundstrom and Moshfegh 2004)
In a turbulent flow, all variable properties contain energy in a wide range of wavenumbers. The Reynolds number of smallest eddies is equal to 1 where length scale is on the order of 0.01-0.1 mm and eddies motion is due to viscosity in this level of eddy’s size (Versteeg and Malalasekera 2007).
2.2.1 Turbulence models
In a century, it was shown experimentally that turbulence problem is always hard to be analysed. But there is a hope to analyse and simulate the turbulent flows more accurate by power of digital computer’s improvement. Therefore turbulence models are improved day by day, to predict turbulent flows in wide range of flows.
Turbulence models solve mean flow equations (2.2), (2.6a-c) and (2.9) in a computational procedure. The velocity field is in three-dimensional as well as dependents to time. Therefore for a CFD code to be useful, it should be able to challenge in some criteria: level of
description (mean quantities, instantaneous quantities), applicability, accuracy (solving wide range of flow conditions), cost and completeness.
Totally turbulence models are classified as
• Classical models: based on time averaged Reynolds equations
1. Zero equation model
2. Two-equation models
3. Reynolds stress equation model
4. Algebraic stress model
• Large eddy simulation (LES)
• Direct numerical simulation (DNS)
Classical models
đang được dịch, vui lòng đợi..
