• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    A Study on the Far Wake of Elliptic Cylinders

    2017-03-13 05:01:01SubhankarSenandSanjayMittal
    關(guān)鍵詞:包辦婚姻幼子秋瑾

    Subhankar Sen and Sanjay Mittal

    1 Introduction

    Most studies concerning external flow past two-dimensional bluff bodies concentrate on regions immediately downstream of the body while the far wake is of no less physical and practical significance. The practical application of far wake analysis includes the non-acoustic detection of submarines near thermoclines (Vorobieff et al., 2002). As an outcome of diffusion and viscous dissipation of Karman vortices, the vortex street might be expected to disappear in far downstream (Zdravkovich, 1997). For flow past a circular cylinder for Reynolds number, Re<200, (Tritton, 1977) commented that `the vortex street continues to all distances downstream'. As pointed out by Cimbala et al. (1988) earlier,this statement deems invalid as the flow characteristics in the near (small multiple of characteristic dimension) and far (large multiple of characteristic dimension) wake vary drastically. The pioneering experimental work of Taneda (1959) demonstrated the formation of `secondary vortex street' in the far wake, where the scales of vortical structures are multiples of those in the primary or Karman vortex street. The wake transition in general, can occur in any staggered vortex street (Vorobieff et al., 2002).With downstream distance, the spatial transition of the wake obeys the following sequence: primary or Karman vortex street → nearly parallel shear layers → secondary vortex street.

    Developments in the wake of two-dimensional symmetric bluff bodies (eg. circular,elliptic cylinders and vertical flat plate) with Re have been studied earlier by several researchers. These studies, in general, do not consider the effect of blockage, B on the wake state. Blockage is the ratio of cross-stream projection of the body and width of the experimental apparatus or computational domain. In the present work, we study the evolution of far wake of elliptic cylinders with angle of attack, α while Re is fixed. The effect of channel confinement on the development of secondary vortex street is also studied. Geometry of an ellipse is characterized by its aspect ratio, AR defined as the ratio of lengths of its minor and major axes denoted by b and a, respectively. Two extreme configurations of an ellipse are the flat plate (AR = 0) and circular cylinder (AR = 1). A review on the earlier studies on unsteady flow past stationary elliptic cylinders and evolution of far wake with Re for various symmetric bluff bodies is provided below.

    Lugt and Haussling (1974) presented results at Re = 15, 30 and 200 for the flow past elliptic cylinders of AR = 0.1 and 0.2 at 450incidence. They employed finite-difference to discretize the unsteady streamfunction-vorticity (ψ-ω) equations. Patel (1981) followed a semi-analytical method for the unsteady ψ-ω equations. For impulsively started elliptic cylinders of various eccentricities and α between 00and 900, he presented detailed results at Re = 60, 100 and 200. Based on numerical study, Park et al. (1989) identified for 25≤Re≤600, various flow regimes for an elliptic cylinder of AR≈0.15. For elliptic cylinders with major axes oriented normal to the incoming stream, Johnson et al. (2001)studied the evolution of wake vortex structures with Re and AR. The parameter space for the numerical study was 30≤Re≤200 and 0.01≤AR≤1. At constant Re, a monotonic increase of the drag coefficient with decreasing AR was noted. For symmetric elliptic cylinders of AR = 0.6, 0.8, 1.0 and 1.2, Kim and Sengupta (2005) presented numerical results highlighting the effects of Re and AR on various integral parameters, such as the aerodynamic coefficients and shedding frequency. The Reynolds numbers studied were 200, 400 and 1000. Experiments at high Re were conducted by Taneda (1972), Modi and Dixit (1975), Ota and Nishiyama (1986) and Ota et al. (1987). Numerical investigations in the high Re regime include those of Panikker and Lavan (1975), Mittal and Balachandar (1996), Badr et al. (2001) and Li et al. (2006). The initial separation of laminar boundary layer for elliptic cylinders of various thicknesses and orientations have been numerically investigated by Sen et al. (2012). They proposed the wake topology for inclined elliptic cylinders and also explored the relationship between separation Reynolds number and angle of incidence.

    Following flow visualization in water tank, Taneda (1959) first observed for a circular cylinder and flat plate, downstream decay of the primary or Karman vortex street and further downstream, eventual formation of the secondary vortex street. Taneda (1959)conjectured that the formation of secondary vortex street is an outcome of changing hydrodynamic stability. By artificially decelerating the Karman street of a circular cylinder, Durgin and Karlsson (1971) demonstrated that the formation of a new(secondary) vortex street is preceded by the deformation and streamwise alignment of primary vortices to form shear layers (referred to as `calm region' by them). An alternate mechanism leading to the formation of secondary vortex street was proposed by Matsui and Okude (1981) who suggested that pairing of Karman vortices rather than hydrodynamic instability of the wake profile leads to the formation of large scale vortices.Based on experimental investigation (flow visualization and hot-wire anemometry) of near and far wakes of two-dimensional bluff bodies, Cimbala et al. (1988) concluded that the far wake structures originate from hydrodynamic instability of the developing mean wake profile. The scale or frequency of the Karman vortices does not affect these large structures. The experimental work of Williamson and Prasad (1993a,b) shows that the far wake structures are sensitive to free-stream disturbances, they originate from hydrodynamic instability and in contrast to the observations made by Cimbala et al.(1988), the far wake frequency is closely linked to the near wake frequency through the frequency of free-stream disturbances. For Re between 70 and 154, Karasudani and Funakoshi (1994) experimentally investigated the breakdown and rearrangement of the primary vortex street of a circular cylinder.

    For a zero-thickness flat plate oriented normal to the flow, Najjar and Balachandar (1998)numerically investigated at Re = 250, the origin of low frequency (roughly one-tenth of Karman frequency) unsteadiness present in the wake and its effects on the Karman frequency, drag and lift amplitudes. Based on numerical flow visualization, they conjectured that a lack of synchronization between the streamwise and spanwise vortices in the wake leads to the low frequency unsteadiness. Inoue and Yamazaki (1999) studied numerically both the unforced and sub-harmonic (of shedding frequency) forced far wakes of a circular cylinder for 140≤Re≤1000. While the mechanism involving merger of Karman vortices played a key role in forced wake transition, no clear evidence was found that this mechanism also leads to unforced wake transition. For an elliptic cylinder of AR= 0.5 with α = 900, Johnson et al. (2001) identified six different wake patterns for 30≤Re≤250. The study indicates that the transition Re for the onset of shedding in far wake is associated with a peak in the Strouhal number and a minimum in the drag coefficient. In a subsequent numerical investigation, Johnson et al. (2004) demonstrated that the frequency of the secondary structures is not a multiple of the primary frequency and hence argued that the secondary structures arise out of an instability in the mean velocity profile. The two-dimensional numerical simulations and experiments conducted by Vorobieff et al. (2002) for Re<1000 show that the streamwise distance for the onset of second wake measured from a circular cylinder decays as Re-1/2. Saha (2007) investigated the developments with Re in the far wake of a vertical flat plate of AR = 0.125. The finitedifference based two-dimensional analysis of Saha (2007) for Re = 30-175, identified various frequencies and also demonstrated that the steady far wake undergoes unsteady transition at Re = 145 marked by the presence of large scale vortical structures. Kumar and Mittal (2012) employed direct numerical simulations as well as global linear stability analysis to investigate the origin of secondary vortex street in the far wake of a circular cylinder at Re = 150. They identified the streamwise location of the dominant modes in the wake that lead to the formation of secondary vortex street. They concluded that the far wake structures arise out of a convective instability. The study demonstrates that suppression of primary vortex street via a slip splitter plate placed suitably, can also suppress the large scale structures in far wake.

    For elliptic cylinders of AR = 0.2, 0.5 and 0.8, we report in the paper, results demonstrating the transitions in the far wake as function of α. For a fixed Re of 200, α varies between 00and 900, in steps of 150. The present numerical investigation also explores the dependence of integral parameters on cylinder thickness and α. The effect of blockage (B = 0.20) on the development of secondary vortex street is studied for the AR =0.2, α = 900configuration. A stabilized Petrov Galerkin finite-element method using equal order bilinear interpolation for velocity and pressure is employed.

    2 The governing equations

    Here u = (u, v), f and σ denote the fluid velocity, body force per unit volume and the Cauchy stress tensor, respectively. The stress is the sum of its isotropic and deviatoric parts:

    where p, I, μ and σ are the pressure, identity tensor, dynamic viscosity of the fluid and strain rate tensor, respectively. Both, the Dirichlet and Neumann-type boundary conditions are accounted for and are represented as

    respectively, where Γgand Γhare complementary subsets of the boundary Γ, n is its unit normal vector and h is the surface traction vector. In the present simulations concerning low blockage, we have employed the towing tank boundary condition (see Figure 1a) on the lateral walls of the domain. This involves free-stream speed condition on the upstream as well as lateral boundaries. No-slip boundary condition is applied on the surface of the cylinder. At the downstream boundary, a Neumann condition for velocity is specified that corresponds to stress-free condition. For high blockage, the inlet flow is fully developed and sidewalls represent no-slip boundaries.

    The initial condition on the velocity is specified on Ω at t = 0: u(x,0)= u0on Ω,

    where u0is divergence-free, i.e. u0satisfies Equation (2).

    3 The finite-element formulation

    The spatial domain Ω is discretized into non-overlapping subdomains Ωe, e =1,2, ..., nelwhere nelis the number of elements. Letandbe the finite dimensional trial function spaces for velocity and pressure, respectively and the corresponding weighting function spaces are denoted byand. The stabilized finite-element formulation of the conservation Equations (1) and (2) is written as follows:

    而秋瑾與王廷均是包辦婚姻。這個(gè)婚姻不可謂不門當(dāng)戶對,一個(gè)是厘金局總辦千金,一個(gè)是湘潭首富之子?;橐龀跗?,詩詞唱和,生活安逸。但此后由于思想的差距,二人漸行漸遠(yuǎn)。王廷均是不是深愛著秋瑾因沒有情書書信流傳不得而知,但是對秋瑾一直不離不棄,赴京捐官,也是因?yàn)榍镨纱私ㄗh。直至后期秋瑾只身留學(xué),王卻帶著兩個(gè)幼子回到鄉(xiāng)中。而且堅(jiān)持不與秋瑾離婚,并在秋瑾死后兩個(gè)月,積郁而死。這兩顆越來越遠(yuǎn)的心,昭示著有一個(gè)包辦婚姻的悲劇。

    In the variational formulation given by Equation 5, the first three terms and the right hand side constitute the Galerkin formulation of the problem. The first series of element level integrals are the SUPG (streamline-upwind/Petrov-Galerkin) and the PSPG (pressurestabilizing/Petrov-Galerkin) stabilization terms added to the variational formulations of the momentum and the continuity equations, respectively. At high Re, in an advection dominated flow, Galerkin formulation of the flow equations lead to oscillations in the velocity field. This numerical instability is overcome by adding SUPG stabilization terms.The SUPG formulation for convection dominated flows was introduced by Hughes and Brooks (1979) and Brooks and Hughes (1982). PSPG stabilization terms are added to the formulation to enable the use of equal order interpolation for velocity and pressure.Hughes et al. (1986) introduced the pressure stabilization methods in the context of Stokes flow and Tezduyar et al. (1992) generalized the method to flows at finite Re. More details of the finite-element formulation can be found in Tezduyar et al. (1992).

    4 Problem set-up

    Figure 1: Re = 200 unsteady flow past elliptic cylinders: (a) problem definition for the unbounded case. Positive sign convention for the aerodynamic forces and moment is shown in (b). The moment is calculated about the center of the cylinder.

    The problem statement for unbounded flow past an elliptic cylinder is shown in Figure 1a.The center of the cylinder coincides with the origin of the Cartesian coordinate system.The angle of attack is measured clockwise with respect to the direction of incoming flow.To simulate unbounded flow, the domain width, H is fixed at 100a (thus, B = 0.01 for α =900and B<0.01 for other values of α) and the towing tank condition is used at the sidewalls of the domain. The Reynolds number and Strouhal number, St are based on a and free-stream speed, U. For confined flow (B = 0.20), the velocity is set to zero at the sidewalls and a parabolic velocity profile is specified at the inlet: u = 1-(2y/H)2, v = 0. Re and St for this flow are based on the centerline speed. For all the computations, distances of the upstream and downstream boundaries measured from the center of the cylinder are Lu= 45a and Ld= 80a, respectively. Figure 1b illustrates the positive sign convention for the aerodynamic forces and moment. The instantaneous force (Cd, Cl) and moment (Cm)coefficients are defined as

    Here nxand nyare the x and y components, respectively of the unit vector n normal to the cylinder boundary Γcyland r is the radius vector of any arbitrary point located on the cylinder surface measured from the center of cylinder.

    5 The finite-element mesh

    For a representative AR of 0.5 and α of 450, Figure 2 shows the non-uniform finiteelement mesh consisting of bilinear quadrilateral elements and its close-up near the cylinder. The number of nodes and elements are 107976 and 107168, respectively and remains constant for all the meshes used. The mesh has been constructed by combining five blocks; a central square block containing the cylinder and four neighbouring rectangular blocks located to the left, right, top and bottom of the central block. To efficiently resolve the secondary structures in the far wake, a sufficiently fine mesh is designed for the right side rectangular block. This block contains 18954 nodes and 18792 elements. For every mesh, the number of nodes on the cylinder, Ntis kept constant at 464 and radial thickness of the first layer of elements located on the cylinder,0.0005a.

    Figure 2: Unsteady flow past elliptic cylinders: (a) the non-uniform finite-element mesh for a 450 inclined elliptic cylinder of AR = 0.5 and (b) close-up of the central square block containing the cylinder. The mesh contains 107976 nodes and 107168 bilinear quadrilateral elements.

    6 Validation of method and convergence of results

    6.1 Comparison with the earlier studies

    In order to ascertain validity of the computed results, predicted time-averaged drag coefficient,and vortex shedding frequency of elliptic cylinders of AR = 0.25 and 0.5 for α = 900are compared with those obtained by Johnson et al. (2001) at Re = 175 and 200, respectively (see Table 1). The comparison in general reveals satisfactory agreement(the maximum difference is about 5% forat Re = 200) between the two sets of results.

    Table 1: Unsteady flow past elliptic cylinders of various AR with major axis oriented normal to the flow: comparison of the predicted and St with those obtained by Johnson et al. (2001).

    Table 1: Unsteady flow past elliptic cylinders of various AR with major axis oriented normal to the flow: comparison of the predicted and St with those obtained by Johnson et al. (2001).

    ?

    Detailed validation of computed results for unsteady flow past stationary elliptic cylinders has been recently discussed in Sourav and Sen (2017).

    6.2 Effect of spatial resolution

    To establish mesh independence of the computed results, the flow past a 450inclined elliptic cylinder of AR= 0.2 is computed at Re = 200 on several meshes. Out of the many meshes, information on meshes M1 and M2 are listed in Table 2. The number of nodes and elements of M2 is roughly twice the ones used for M1. The non-dimensional time step size, ΔtU/a for the computations is 0.005. The results summarized in Table 2 indicate that the spatial resolution of M1 is sufficient to accurately predict the unsteady flow past elliptic cylinders.

    Table 2: Flow past an elliptic cylinder of AR= 0.2 at α = 450: effect of mesh resolution on the flow characteristics at Re = 200. The other parameters for this study are: Lu = 45a, Ld= 80a, = 0.0005a and ΔtU/a = 0.005. For M1, Nt = 464 and for M2, Nt = 672. The subscript r.m.s. implies root mean square of a quantity.

    Table 2: Flow past an elliptic cylinder of AR= 0.2 at α = 450: effect of mesh resolution on the flow characteristics at Re = 200. The other parameters for this study are: Lu = 45a, Ld= 80a, = 0.0005a and ΔtU/a = 0.005. For M1, Nt = 464 and for M2, Nt = 672. The subscript r.m.s. implies root mean square of a quantity.

    ?

    6.3 Effect of time step size

    To examine the effect of time step size on the flow characteristics, the Re = 200 flow past a 450inclined elliptic cylinder of AR = 0.2 is computed on mesh M1 using ΔtU/a = 0.005 and 0.0025. Table 3 summarizes the results and establishes the adequacy of ΔtU/a =0.005 for the present computations. Based on the spatial and temporal resolution tests, all computations in the present work are performed on mesh M1 using ΔtU/a = 0.005.

    Table 3: Flow past an elliptic cylinder of AR = 0.2 at α = 450: effect of time step size on the flow characteristics at Re = 200. The computations are carried out on Mesh M1.

    7 Results

    Results are presented for unsteady flow at Re = 200 past elliptic cylinders of aspect ratios 0.2, 0.5 and 0.8 as well as angles of attack between 00and 900. Each numerical computation uses the corresponding steady flow solution at Re = 40 as the initial condition. The linearized algebraic equation system has been solved by GMRES or Generalized Minimal RESidual method (Saad and Schultz, 1986) in conjunction with a diagonal preconditioner.

    7.1 Flow characteristics: the integral parameters

    Figure 3: Re = 200 unsteady flow past elliptic cylinders of AR = 0.2-0.8: variation of the time-averaged (first row) and r.m.s. (second row) values of aerodynamic coefficients with α.

    For any orientation of an elliptic cylinder, amplitude of lift oscillation is the largest while moment oscillates with the smallest amplitude (Sen, 2010). Thus, relative amplitudes of the aerodynamic coefficients suggest Clrms> Cdrms, Cmrms. Figures 3(a,d) indicate that the time-averaged and r.m.s. drag of elliptic cylinders continue to increase with increasing α and attain the maxima when α = 900. For α≥600,and Cdrmsincrease with decreasing AR. As seen from Figure 3b,of each cylinder initially increases with increasing α and then decreases. The average lift and moment vanish for symmetric configurations of α = 00and 900. When α is fixed, the lift generated by a cylinder increases with decreasing thickness. Figure 3e demonstrates that Clrmsfor AR = 0.2 displays a non-monotonic variation with α. Interestingly, a monotonic rise is seen for AR≥0.5. With α,of each cylinder displays non-monotonic variation (Figure 3c) and r.m.s. moment coefficient, Cmrmsdisplays monotonic rise (Figure 3f) for α≥300.

    To analyze the primary vortex shedding characteristics, Fast Fourier Transform (FFT) is performed on the signals of aerodynamic coefficients. For each AR and α, the lift and moment coefficients are found to oscillate at the primary vortex shedding or Strouhal frequency. The drag coefficient in general, oscillates at twice the primary frequency.Exceptions to this are found for AR≤0.5 at low α. When α≤450, FFTs of Cdand Clfor AR= 0.2 lead to the same value of St. The power spectra of Cdand Clfor AR = 0.2, α = 450reveal three and two distinct peaks, respectively (Figures 4a and 4b). The dominant peak of each spectrum corresponds to the same frequency (St = 0.2431) and the drag-lift phase diagram consists of a single loop (Figure 4c). For AR = 0.5, values of St obtained from FFTs of Cdand Clconverge when α ≤300. Figures 4(d,e) illustrate for AR = 0.5 and α =900, the power spectra of Cdand Cl. For this case, the drag oscillates at twice the shedding frequency. Phase diagram for this configuration (Figure 4f) is of the shape of a distorted figure eight characterized by the presence of a cross-over point. Phase diagrams of Figures 4(c,f) are closed and ascertain periodic nature of flow. Dominance of the amplitude of lift over drag is also apparent from these plots.

    Figure 4: Unsteady flow past elliptic cylinders at Re = 200: non-dimensional power spectra of the (a) drag, (b) lift coefficients and (c) the drag-lift phase diagram for AR =0.2, α = 450. Figures (d), (e) and (f) respectively plot the same quantities for AR = 0.5 and α = 900.

    Figure 5a indicates that St of an elliptic cylinder decays monotonically with increasing α.St of a thick cylinder (AR = 0.8) is relatively less sensitive to α. The decaying trend of St with α is opposite to the one observed for a square cylinder. Slow monotonic growth of St with increasing α was reported by Sohankar et al. (1998) for a square cylinder with 00≤α≤450. For α≤450, St decreases monotonically with decreasing AR. The trend reverses for α≥750. The predicted vortex shedding frequencies are in close agreement (maximum deviation being 2.26%) with those obtained by Sheard (2007) for AR = 0.5, α = 150and 300. The intervortex spacing along each row of the vortex street or primary wavelength,a1is inversely proportional to the local St. In particular, a1/a = 1/St (Vorobieff et al.,2002). Figure 5b suggests monotonic growth of a1/a with increasing α. a1/a of a cylinder attains the maximum when α = 900.

    Figure 5: Re = 200 unsteady flow past elliptic cylinders of AR = 0.2-0.8 in the unbounded medium: variation of (a) Strouhal number and (b) primary wavelength, a1/a with α.

    7.2 The evolution of wake with α

    We study the wake transition for a thin cylinder (AR = 0.2). To identify various frequencies, time histories of cross-flow component of velocity are recorded at sixteen locations in the wake along the x axis. The velocity probes are located at approximate downstream distances of 2a, 6a, 10a, 15a, 20a, 25a, 30a, 35a, 40a,45a, 50a, 55a, 60a, 65a, 70a and 75a, respectively measured from the cylinder center.To resolve the frequencies, about 218data points are used for performing FFT of the cross-flow component of velocity.

    7.2.1 The instantaneous vorticity field

    Figure 6 shows that the flow is steady and vortex shedding is absent at low α, such as 00and 150. The flow becomes unsteady at some critical α between 150and 300and Karman vortex shedding from the cylinder ensues. When α = 300, the wake is comprised of two parts. The first part close to the cylinder is the Karman vortex street where the alternately shed vortices appear as concentric circles. In the second part, the vortices are deformed to elliptical cross-section and orient themselves with major axes parallel to the free-stream as observed by Durgin and Karlsson (1971). The mechanisms leading to the deformation of Karman vortices were explained by Durgin and Karlsson (1971) and later, Karasudani and Funakoshi (1994). These deformed vortices form two parallel and opposite signed shear layers on either side of the wake centerline that extend downstream and appear not to interact mutually. A markedly different wake pattern is observed as α approaches 450.The shear layers become unstable, begin mutual interaction and terminate well ahead of the domain exit. Also, width of individual shear layers increase relative to the width at α= 300. The secondary vortex street characterized by large scale vortical structures, width and wavelength than the primary, is seen to develop at x ≈ 50a.

    Figure 6: Re = 200 unsteady flow past an elliptic cylinder of AR = 0.2 in an unbounded medium: instantaneous vorticity field for various values of α between 00 and 900.

    Rapid decay of the streamwise extent of the shear layers is observed for α = 600. The far wake vorticity distribution for this configuration is characterized by the one-sided convection of the secondary vortex pairs (comprising of vortices of opposite signs) away from the x axis. With increasing α, breakdown of the primary vortex street occurs closer to the cylinder. The shear layers gradually shrink in streamwise extent and attain the minimum when α = 900. Consequently, the location of formation of secondary vortex street moves upstream. For α ≥750, the secondary vortices interact strongly with each other and group of three vortices form with two of them having the same sign. It is shown in Section 7.2.2 that a tertiary frequency of low magnitude appears in the far downstream for α≥750. When α = 900, the secondary structures are visible sufficiently close to the cylinder, the distance being of the order of 20a. The developments in the wake with α for a fixed Re appear analogous to the cases of transition of wake with Re for a cylinder of fixed orientation (Johnson et al., 2001; Johnson et al., 2004; Saha, 2007). Though not shown here, it is found that the location of formation of secondary vortices moves far downstream with increasing AR.

    7.2.2 The frequency distribution

    For α = 00and 150, the flow is steady and vortex-shedding frequency is absent. The transverse velocity component for α = 300is time periodic in the near wake and the fluctuations diminish at large distance downstream implying steady nature of the far wake (Figure 7a). The coexistence of unsteady near wake and steady far wake was earlier pointed out by Saha (2007) for the flow past a vertical flat plate at Re = 100. The wake state for α = 300is characterized by the primary frequency alone. This is evident from Figure 7b which also demonstrates the weakening and eventual disappearance of the shedding frequency in the far wake.

    Figure 7: Re = 200 unsteady flow past a 300 inclined elliptic cylinder of AR = 0.2 in an unbounded medium: (a) time series of v velocity signals at x/a = 2 and 75, (b) normalized power spectra of v velocity at selected downstream locations along the x axis.

    Figure 8 shows for α = 450and 600, frequency characteristics at selected downstream locations. The secondary frequency first appears at α = 450and when α is increased to 600,this frequency is realized relatively close to the cylinder. At the location of probes, Figure 9 shows the power spectra of cross-flow component of velocity for α = 750. Associated with α≥750are several additional frequencies. The most prominent among them is the low magnitude tertiary frequency that appears in the far wake. The tertiary frequency becomes the most dominant for x≥60a.

    Figure 8: Re = 200 unsteady flow past an inclined elliptic cylinder of AR = 0.2 in an unbounded medium: normalized power spectra of v velocity at various downstream locations for α = (a) 450 and (b) 600.

    Figure 9: Re = 200 unsteady flow past a 750 inclined elliptic cylinder of AR = 0.2 in an unbounded medium: normalized power spectra of v velocity at various downstream locations.

    It is found that the values of primary shedding frequency obtained by performing FFTs of the time series of Cland v are in excellent agreement. Values of various frequencies present in the wake are listed in Table 4. The primary frequency of a cylinder has the largest magnitude and the magnitude of tertiary frequency, if it exists, is the smallest. The secondary frequency has a magnitude intermediate to the ones for primary and tertiary. A non-monotonic variation of the secondary frequency with α is apparent from the table. In contrast, the primary frequency exhibits monotonic decrease with increasing α (also see Figure 5a). It is also evident from the table that the primary frequency is not a multiple of the secondary for α = 450, 750and 900. Interestingly, for α = 600, the secondary frequency is a sub-harmonic of the primary. The primary frequency in most cases is not a multiple of the secondary frequency suggesting that the secondary vortices do not form out of interaction of the shed primary vortices (also see Karasudani and Funakoshi (1994) and Johnson et al. (2004)). The distribution of dominant frequencies is shown in Figure 10 for 300≤α≤900. Among the primary, secondary and tertiary frequencies, secondary frequency exhibits the least sensitivity to α.

    Table 4: Flow past an elliptic cylinder of AR = 0.2 for α = 300-900: summary of various frequencies present in the wake at Re = 200.

    Figure 10: Re = 200 unsteady flow past an elliptic cylinder of AR = 0.2 in unbounded medium: distribution of dominant frequencies along the wake centerline for 300≤α≤900.

    7.3 Effect of blockage on the far wake

    We next investigate the effect of channel confinement on the wake structures of the AR =0.2 cylinder at Re = 200. To meet this objective, a blockage of 0.20 is used (see Section 4). FFT of cross-flow component of velocity measured at all the sixteen probes leads to the same value of primary frequency (St = 0.2232). This value is significantly higher than the one for unbounded flow (St = 0.1569 from Table 4).As apparent from Figure 11, the primary frequency prevails over all locations along the wake centerline and the wake state is characterized by the absence of any other significant frequency. The immediate effect of blockage is therefore to suppress the formation of secondary vortex street. Absence of secondary vortices in the far wake is also evident from Figure 12 showing the instantaneous vorticity field.

    Figure 11: Confined flow (B = 0.20) past a perpendicularly oriented elliptic cylinder of AR = 0.2 at Re = 200: frequency distribution at various locations along the x axis.

    Figure 12: Instantaneous vorticity field at Re = 200 for the confined flow past an elliptic cylinder of AR = 0.2 and α = 900. For this figure, Ld = 80a and B = 0.20.

    In order to investigate if the imposition of blockage pushes the onset of secondary wake further downstream, computations are performed by doubling Ld, i.e. Ld= 160a while other parameters remain unaltered. Apart from the previously employed sixteen velocity probes, eleven more are placed along the wake centerline at approximate downstream distances of 82.5a, 90a, 97.5a, 105a, 112.5a, 120a, 127.5a, 135a, 142.5a, 150a and 157.5a measured from the cylinder center. The frequency spectra presented in Figure 13 for certain streamwise locations along the x axis, demonstrate monofrequency (St =0.2232) nature of the wake and establish that imposition of blockage can suppress the formation of secondary vortices.

    Figure 13: Confined flow past a perpendicularly oriented elliptic cylinder of AR = 0.2 at Re = 200: normalized power spectra of cross-flow component of velocity component at various locations along the x axis. For this figure, Ld = 160a and B = 0.20.

    8 Conclusions

    The transition of wake of elliptic cylinders with varying angle of attack (00≤α≤900) is investigated using a stabilized finite-element method. Results are presented for twodimensional simulations at Re = 200. Stationary elliptic cylinders of aspect ratios 0.2, 0.5 and 0.8 are considered. For each AR and α, the lift coefficient oscillates with the largest amplitude. The smallest amplitude oscillations are in general associated with moment coefficient. The r.m.s. value of lift exceeds the ones for drag and moment. Irrespective of AR, the time-averaged and r.m.s. drag increase monotonically with increasing α. In contrast, the time-averaged lift and moment initially increase and then decrease. With increasing α, the primary shedding frequency of each cylinder continues to decay and the primary wavelength exhibits a reverse trend. In general, the shedding frequency obtained via FFT of lift is half the value obtained from FFT of drag. Interestingly, both FFTs lead to the same frequency for AR≤0.5 at low α. For AR = 0.2, the flow is steady for α = 00and 150. An unsteady wake characterized by Karman vortices and associated primary frequency, is seen as α increases to 300. The large scale secondary structures appear first when α = 450. The secondary frequency is also captured first at this incidence angle. For α≥750, analysis of power spectra shows the existence of a low magnitude tertiary frequency in the far wake. Compared to the primary and tertiary counterparts, secondary frequency displays lesser sensitivity to α. For most of the angles of attack studied, the secondary frequency is not a sub-harmonic of primary, signifying that the secondary structures do not form due to coalescence of primary vortices. For each AR, the location of formation of secondary structures moves upstream with increasing angle of attack.When α is fixed, similar behaviour is observed with decreasing cylinder thickness. For confined flow, the instantaneous vorticity field and power spectra of transverse velocity suggest that blockage suppresses the formation of secondary vortices.

    Badr, H. M., Dennis, S. C. R. and Kocabiyik, S., (2001). Numerical simulation of the unsteady flow over an elliptic cylinder at different orientations. Int. J. Num. Meth. Fluids,vol. 37, pp. 905-931.

    Brooks, A. N. and Hughes, T. J. R., (1982). Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations. Comput. Methods Appl. Mech. Eng., vol. 32, pp.199-259.

    Cimbala, J. M., Nagib, H. M. and Roshko, A, (1988). Large structure in the far wakes of two-dimensional bluff bodies. J. Fluid Mech., vol. 190, pp. 265-298.

    Durgin, W. W. and Karlsson, S. K. F., (1971). On the phenomenon of vortex street breakdown. J. Fluid Mech., vol. 48, pp. 507-527.

    Hughes, T. J. R. and Brooks, A. N., (1979). A multi-dimensional upwind scheme with no crosswind diffusion. Finite Element Methods for Convection Dominated Flows (ed. T.J. R. Hughes), vol. 47, pp. 19-35.

    Hughes, T. J. R., Franca, L. P. and Balestra, M., (1986). A new finite element formulation for computational fluid dynamics: V. Circumventing the Babuska-Brezzi condition: A stable Petrov-Galerkin formulation of the Stokes problem accommodating equal-order interpolations. Comput. Methods Appl. Mech. Eng., vol. 59, pp. 85-99.

    Inoue, O. and Yamazaki, T., (1999). Secondary vortex streets in two-dimensional cylinder wakes. Fluid Dyn. Res., vol. 25, pp. 1-18.

    Johnson, S. A., Thompson, M. C. and Hourigan, K., (2001). Flow past elliptical cylinders at low Reynolds numbers. In Proc. 14thAustralasian Fluid Mechanics Conference, Adelaide, pp. 343-346.

    Johnson, S. A., Thompson, M. C. and Hourigan, K., (2004). Predicted low frequency structures in the wake of elliptical cylinders. European J. Mech. B Fluids, vol. 23, pp.229-239.

    Karasudani, T. and Funakoshi, M., (1994). Evolution of a vortex street in the far wake of a cylinder. Fluid Dyn. Res., vol. 14, pp. 331-352.

    Kim, M. S. and Sengupta, A., (2005). Unsteady viscous flow over elliptic cylinders at various thickness with different Reynolds numbers. J. Mech. Sci. Tech. (KSME Int. J.),vol. 19, pp. 877-886.

    Kumar, B. and Mittal, S., (2012). On the origin of the secondary vortex street. J. Fluid Mech., vol. 711, pp. 641-666.

    Li, Z., Davidson, J. H. and Mantell, S., (2006). Numerical simulation of flow field and heat transfer of streamlined cylinders in cross flow. Trans. ASME: J. Fluids Eng., vol. 128,pp. 564-570.

    Lugt, H. J. and Haussling, H. J.,(1974). Laminar flow past an abruptly accelerated elliptic cylinder at 450incidence. J. Fluid Mech., vol. 65, pp. 711-734.

    Matsui, T. and Okude, M., (1981). Vortex pairing in a Karman vortex street. In Proc.Seventh Biennial Symp. Turbulence.

    Mittal, R. and Balachandar, S., (1996). Direct numerical simulation of flow past elliptic cylinders. J. Comput. Phys., vol. 124, pp. 351-367.

    Modi, V. J. and Dikshit, A. K., (1975). Near-wakes of elliptic cylinders in subcritical flow. AIAA J., vol. 13, pp. 490-497.

    Najjar, F. M. and Balachandar, S., (1998). Low-frequency unsteadiness in the wake of a normal flat plate. J. Fluid Mech., vol. 370, pp. 101-147.

    Ota, T. and Nishiyama, H., (1986). Flow around two elliptic cylinders in tandem arrangement. Trans. ASME: J. Fluids Eng., vol. 108, pp. 98-103.

    Ota, T., Nishiyama, H. and Taoka, Y., (1987). Flow around an elliptic cylinder in the critical Reynolds number regime. Trans. ASME: J. Fluids Eng., vol. 109, pp. 149-155.

    Panikker, P. K. G. and Lavan, Z., (1975). Flow past impulsively started bodies using Green's functions. J. Comput. Phys., vol. 18, pp. 46-65.

    Park, J. K., Park, S. O. and Hyun, J. M., (1989). Flow regimes of unsteady laminar flow past a slender elliptic at incidence. Int. J. Heat Fluid Flow, vol. 10, pp. 311-317.

    Patel, V. A., 1981. Flow around the impulsively started elliptic cylinder at various angles of attack. Comput. Fluids, vol. 9, pp. 435-462.

    Saad, Y.; Schultz, M. H., (1986). GMRES: A generalised minimal residual algorithm for solving nonsymmetric linear systems. SIAM J. Sci. Stat. Comput., vol. 7, pp. 856-869.

    Saha, A. K., (2007). Far-wake characteristics of two-dimensional flow past a normal flat plate. Phys. Fluids, vol. 19, pp. 128110-1-128110-4.

    Sen, S., (2010). Flow past stationary and vibrating cylinders of various cross-sections at low Reynolds numbers. PhD Thesis, Indian Institute of Technology Kanpur.

    Sen, S., Mittal, S. and Biswas, G., (2012). Steady separated flow past elliptic cylinders using a stabilized finite-element method. Comput. Model. Eng. Sci., vol. 86, pp 1-27.

    Sheard, G. J., (2007). Cylinders with elliptic cross-section: wake stability with variation in angle of incidence. In Proc. IUTAM Symp. Unsteady separated flow and their control.

    Sohankar, A., Norberg, C. and Davidson, L., (1998). Low-Reynolds-number flow around a square cylinder at incidence: study of blockage, onset of vortex shedding and outlet boundary condition. Int. J. Num. Meth. Fluids, vol. 26, pp. 39-56.

    Taneda, S., (1959). Downstream development of the wakes behind cylinders. J. Phys.Soc. Japan, vol. 14, pp. 843-848.

    Taneda, S., (1972). The development of the lift of an impulsively started elliptic cylinder at incidence. J. Phys. Soc. Japan, vol. 33, pp. 1706-1711.

    Tezduyar, T. E., Mittal, S., Ray, S. E. and Shih, R., (1992). Incompressible flow computations with stabilized bilinear and linear equal-order-interpolation velocitypressure elements. Comput. Methods Appl. Mech. Eng., vol. 95, pp. 221-242.

    Tritton, D. J., (1977). Physical~Fluid~Dynamics. Van Nostrand Reinhold (UK).

    Vorobieff, P., Goergiev, D. and Ingber, M. S., (2002). Onset of the second wake:dependence on the Reynolds number. Phys. Fluids, vol. 14, pp. L53-L56.

    Williamson, C. H. K. and Prasad, A., (1993a). Wave interactions in the far wake of a body. Phys. Fluids, vol. 5, pp. 1854-1856.

    Williamson, C. H. K. and Prasad, A., (1993b). A new mechanism for oblique wave resonance in the `natural' far wake. J. Fluid Mech., vol. 256, pp. 269-313.

    Zdravkovich, M. M., (1997). Flow Around Circular Cylinders, vol. 1. Oxford University Press.

    猜你喜歡
    包辦婚姻幼子秋瑾
    信仰的真與執(zhí)著——讀《秋瑾集》
    兩岸中秋
    作品集
    青年生活(2019年14期)2019-10-21 02:04:56
    “包辦婚姻”
    幸福的包辦婚姻
    巧分牛
    秋瑾年譜的幾個(gè)問題
    不滿包辦婚姻 埃男子揮刀自宮
    生命的祭奠
    意林(2008年7期)2008-05-14 22:02:40
    亚洲情色 制服丝袜| 免费搜索国产男女视频| 在线观看日韩欧美| 可以免费在线观看a视频的电影网站| 久久人人爽av亚洲精品天堂| 国产亚洲精品一区二区www| 亚洲电影在线观看av| 大型黄色视频在线免费观看| 变态另类丝袜制服| 69精品国产乱码久久久| 色尼玛亚洲综合影院| 国产精品亚洲美女久久久| 欧美亚洲日本最大视频资源| 美女高潮喷水抽搐中文字幕| 国产色视频综合| 午夜福利成人在线免费观看| 亚洲国产日韩欧美精品在线观看 | 91字幕亚洲| 久久中文看片网| 国产精品 国内视频| 久久久国产成人精品二区| 精品国产乱子伦一区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 欧美日本亚洲视频在线播放| 琪琪午夜伦伦电影理论片6080| 真人一进一出gif抽搐免费| 欧美中文综合在线视频| 一本大道久久a久久精品| 99久久久亚洲精品蜜臀av| 国产一卡二卡三卡精品| 久久久久久国产a免费观看| 99国产精品一区二区蜜桃av| 亚洲国产中文字幕在线视频| 亚洲va日本ⅴa欧美va伊人久久| 黑人巨大精品欧美一区二区蜜桃| 女生性感内裤真人,穿戴方法视频| a级毛片在线看网站| 香蕉国产在线看| av超薄肉色丝袜交足视频| 久久性视频一级片| 非洲黑人性xxxx精品又粗又长| 91麻豆av在线| 国产日韩一区二区三区精品不卡| 国产亚洲精品av在线| 亚洲国产精品合色在线| 国产精品1区2区在线观看.| 色哟哟哟哟哟哟| 国产野战对白在线观看| 久久精品国产亚洲av香蕉五月| 久久精品影院6| 精品免费久久久久久久清纯| 高清毛片免费观看视频网站| 国产麻豆69| 亚洲av电影不卡..在线观看| 成人三级做爰电影| 国产av精品麻豆| 制服丝袜大香蕉在线| 一本久久中文字幕| 欧美一级毛片孕妇| 亚洲国产精品久久男人天堂| 色婷婷久久久亚洲欧美| 免费看a级黄色片| 岛国视频午夜一区免费看| 丰满的人妻完整版| 久久久国产精品麻豆| 中国美女看黄片| 久久草成人影院| 久久青草综合色| 日韩中文字幕欧美一区二区| 日韩精品免费视频一区二区三区| 国产精品自产拍在线观看55亚洲| 一边摸一边抽搐一进一小说| 中文字幕人妻熟女乱码| 色av中文字幕| 法律面前人人平等表现在哪些方面| 精品一区二区三区av网在线观看| 亚洲一区二区三区不卡视频| 在线观看免费视频日本深夜| 成人国产一区最新在线观看| 日本 av在线| 十八禁网站免费在线| 叶爱在线成人免费视频播放| 午夜免费激情av| 在线av久久热| 亚洲情色 制服丝袜| 超碰成人久久| 久久久久久人人人人人| 男人舔女人下体高潮全视频| 国产精品 国内视频| 国产片内射在线| 在线观看66精品国产| 黄色女人牲交| 日韩有码中文字幕| 国产成+人综合+亚洲专区| 亚洲成人精品中文字幕电影| 老鸭窝网址在线观看| 亚洲男人的天堂狠狠| 成人国产综合亚洲| 99国产精品免费福利视频| tocl精华| 精品国产国语对白av| 精品人妻1区二区| 黄频高清免费视频| 久久影院123| 亚洲精品国产区一区二| 亚洲人成77777在线视频| 国产国语露脸激情在线看| 男男h啪啪无遮挡| 亚洲av成人不卡在线观看播放网| 亚洲欧美一区二区三区黑人| 正在播放国产对白刺激| 久久人人97超碰香蕉20202| 91成年电影在线观看| 黄色片一级片一级黄色片| 丝袜美足系列| 日韩 欧美 亚洲 中文字幕| 亚洲黑人精品在线| 人人澡人人妻人| 悠悠久久av| 女性生殖器流出的白浆| 亚洲成国产人片在线观看| 久久精品亚洲熟妇少妇任你| 国产一卡二卡三卡精品| 很黄的视频免费| 久久精品亚洲熟妇少妇任你| 18禁国产床啪视频网站| bbb黄色大片| 国产精品亚洲美女久久久| 男人舔女人下体高潮全视频| 午夜精品在线福利| 久久国产亚洲av麻豆专区| 日本五十路高清| 亚洲七黄色美女视频| www.精华液| 久久国产精品影院| 女性生殖器流出的白浆| 日韩欧美在线二视频| 亚洲人成网站在线播放欧美日韩| av有码第一页| 老司机福利观看| 一级a爱视频在线免费观看| 深夜精品福利| aaaaa片日本免费| 欧美 亚洲 国产 日韩一| 午夜福利一区二区在线看| 亚洲人成电影观看| av电影中文网址| 久久人人精品亚洲av| 性色av乱码一区二区三区2| 久久伊人香网站| 在线观看www视频免费| 亚洲精品粉嫩美女一区| 成人三级做爰电影| 一区二区三区高清视频在线| 日日干狠狠操夜夜爽| 欧美日韩亚洲国产一区二区在线观看| 不卡一级毛片| 久久人妻熟女aⅴ| 91精品三级在线观看| 免费高清在线观看日韩| 19禁男女啪啪无遮挡网站| 国产99白浆流出| av视频免费观看在线观看| 国产av一区二区精品久久| 99国产综合亚洲精品| 99精品久久久久人妻精品| av天堂久久9| 丝袜美足系列| 日韩中文字幕欧美一区二区| 精品人妻在线不人妻| 国产免费av片在线观看野外av| 看免费av毛片| 9热在线视频观看99| 国产高清激情床上av| 国产亚洲精品久久久久5区| 可以在线观看的亚洲视频| 久久人人爽av亚洲精品天堂| 女人爽到高潮嗷嗷叫在线视频| 欧美一区二区精品小视频在线| 一进一出好大好爽视频| 国产一区二区三区在线臀色熟女| 一本久久中文字幕| 欧洲精品卡2卡3卡4卡5卡区| 欧美乱码精品一区二区三区| 国产一区二区三区视频了| 欧美午夜高清在线| 国产成人精品无人区| 两人在一起打扑克的视频| 国产精品电影一区二区三区| 动漫黄色视频在线观看| 色综合亚洲欧美另类图片| 久久久久久久久久久久大奶| 欧美日韩亚洲国产一区二区在线观看| 欧美乱色亚洲激情| 午夜视频精品福利| 香蕉丝袜av| 成熟少妇高潮喷水视频| 大型黄色视频在线免费观看| 桃红色精品国产亚洲av| 亚洲最大成人中文| 琪琪午夜伦伦电影理论片6080| 久久久久国产一级毛片高清牌| 亚洲人成网站在线播放欧美日韩| 久久久久久大精品| 亚洲人成77777在线视频| 99国产极品粉嫩在线观看| 国产在线观看jvid| 一夜夜www| 99久久综合精品五月天人人| 久久精品亚洲精品国产色婷小说| 亚洲少妇的诱惑av| 国产精品久久久久久精品电影 | 欧美色欧美亚洲另类二区 | 精品电影一区二区在线| 午夜福利一区二区在线看| 母亲3免费完整高清在线观看| 一卡2卡三卡四卡精品乱码亚洲| 天天一区二区日本电影三级 | 女警被强在线播放| 婷婷丁香在线五月| av电影中文网址| 神马国产精品三级电影在线观看 | 在线观看午夜福利视频| 成人国语在线视频| 韩国av一区二区三区四区| 伊人久久大香线蕉亚洲五| 久久久国产精品麻豆| 日本撒尿小便嘘嘘汇集6| 久久 成人 亚洲| 亚洲精品国产区一区二| 欧美激情极品国产一区二区三区| 成人手机av| 久久久久久久久免费视频了| 日本黄色视频三级网站网址| 精品乱码久久久久久99久播| 丝袜人妻中文字幕| 国产成人av激情在线播放| 黄网站色视频无遮挡免费观看| 欧美老熟妇乱子伦牲交| 久久精品91无色码中文字幕| 亚洲在线自拍视频| 色婷婷久久久亚洲欧美| 天堂√8在线中文| 免费久久久久久久精品成人欧美视频| 日韩精品青青久久久久久| 久久久久久久久久久久大奶| 亚洲av成人一区二区三| 亚洲,欧美精品.| 又大又爽又粗| 97人妻天天添夜夜摸| 色综合站精品国产| 久久人妻av系列| 99国产极品粉嫩在线观看| 色播亚洲综合网| 亚洲国产欧美网| 91麻豆精品激情在线观看国产| 亚洲成国产人片在线观看| 国产高清视频在线播放一区| av网站免费在线观看视频| 国产精品电影一区二区三区| 婷婷六月久久综合丁香| 日韩一卡2卡3卡4卡2021年| 国产精品98久久久久久宅男小说| 中文字幕另类日韩欧美亚洲嫩草| 一级毛片高清免费大全| 91精品三级在线观看| 午夜福利成人在线免费观看| 黑人巨大精品欧美一区二区蜜桃| 男人操女人黄网站| 精品少妇一区二区三区视频日本电影| 国产伦人伦偷精品视频| 亚洲色图综合在线观看| 亚洲成人国产一区在线观看| 欧美一级毛片孕妇| 亚洲 国产 在线| 欧美成人一区二区免费高清观看 | 一区二区三区激情视频| 久久精品国产亚洲av香蕉五月| 在线观看免费日韩欧美大片| 精品国产国语对白av| 久久精品国产亚洲av香蕉五月| 国产av又大| 亚洲中文日韩欧美视频| 亚洲国产精品sss在线观看| 又黄又爽又免费观看的视频| 亚洲专区中文字幕在线| 欧美黑人精品巨大| av福利片在线| 亚洲国产日韩欧美精品在线观看 | 久久久久久免费高清国产稀缺| 免费看美女性在线毛片视频| 亚洲第一电影网av| 侵犯人妻中文字幕一二三四区| 久久香蕉国产精品| 亚洲第一欧美日韩一区二区三区| 午夜激情av网站| 国产精品久久久av美女十八| 啪啪无遮挡十八禁网站| 91大片在线观看| www.自偷自拍.com| 国产伦人伦偷精品视频| 国产精品一区二区在线不卡| 大型黄色视频在线免费观看| 天天添夜夜摸| 国产精品美女特级片免费视频播放器 | 亚洲成人久久性| 国产精品自产拍在线观看55亚洲| 多毛熟女@视频| 精品电影一区二区在线| 亚洲一码二码三码区别大吗| 日本三级黄在线观看| 日日爽夜夜爽网站| 国产成人系列免费观看| 一级作爱视频免费观看| 又黄又粗又硬又大视频| 两个人免费观看高清视频| 午夜福利在线观看吧| 午夜久久久久精精品| 欧美日本亚洲视频在线播放| 亚洲成av人片免费观看| 不卡av一区二区三区| 天堂动漫精品| 精品久久蜜臀av无| 老司机靠b影院| 日本a在线网址| 国产片内射在线| www.自偷自拍.com| 国产成人影院久久av| 在线观看舔阴道视频| 宅男免费午夜| 国产亚洲欧美精品永久| 12—13女人毛片做爰片一| 国产成人精品无人区| 精品国内亚洲2022精品成人| 在线观看一区二区三区| 无人区码免费观看不卡| 一区二区三区国产精品乱码| 此物有八面人人有两片| 亚洲aⅴ乱码一区二区在线播放 | 亚洲精品粉嫩美女一区| 国产成人啪精品午夜网站| 嫩草影视91久久| 午夜福利视频1000在线观看 | 精品国产乱码久久久久久男人| 欧美大码av| 国产三级在线视频| 男人操女人黄网站| 色老头精品视频在线观看| 亚洲成av片中文字幕在线观看| 欧美一级毛片孕妇| 国产成人欧美| 亚洲无线在线观看| www.熟女人妻精品国产| 亚洲一区二区三区色噜噜| av在线天堂中文字幕| 久久伊人香网站| 成人国产一区最新在线观看| 亚洲七黄色美女视频| 国内精品久久久久精免费| 免费少妇av软件| 一本久久中文字幕| 好看av亚洲va欧美ⅴa在| 丰满的人妻完整版| 亚洲人成伊人成综合网2020| 久热爱精品视频在线9| 露出奶头的视频| 亚洲欧洲精品一区二区精品久久久| 亚洲精品国产区一区二| 大型av网站在线播放| 久久 成人 亚洲| 欧美黄色片欧美黄色片| 一级毛片女人18水好多| 久久久久久久午夜电影| 可以在线观看毛片的网站| www.999成人在线观看| 久久国产精品人妻蜜桃| АⅤ资源中文在线天堂| 精品高清国产在线一区| 天天添夜夜摸| 亚洲最大成人中文| 亚洲aⅴ乱码一区二区在线播放 | av电影中文网址| 黑丝袜美女国产一区| 亚洲午夜理论影院| 成在线人永久免费视频| 91九色精品人成在线观看| 久久久精品欧美日韩精品| 国产亚洲精品av在线| 99国产综合亚洲精品| 亚洲专区字幕在线| 一个人免费在线观看的高清视频| 在线国产一区二区在线| 成人免费观看视频高清| 一区二区三区国产精品乱码| 91av网站免费观看| 老熟妇乱子伦视频在线观看| a在线观看视频网站| 精品国产亚洲在线| 可以在线观看的亚洲视频| 午夜福利高清视频| 免费不卡黄色视频| 国产一区在线观看成人免费| 欧美中文综合在线视频| 91麻豆精品激情在线观看国产| 精品第一国产精品| 美女高潮到喷水免费观看| 中文字幕av电影在线播放| 久久久久久大精品| 亚洲va日本ⅴa欧美va伊人久久| 搡老妇女老女人老熟妇| 亚洲精品一卡2卡三卡4卡5卡| 又大又爽又粗| 女人被躁到高潮嗷嗷叫费观| x7x7x7水蜜桃| 啦啦啦免费观看视频1| 国产成人啪精品午夜网站| 国产欧美日韩一区二区三| 亚洲天堂国产精品一区在线| 99热只有精品国产| 99国产精品99久久久久| 亚洲人成77777在线视频| 咕卡用的链子| 狠狠狠狠99中文字幕| 亚洲一区中文字幕在线| 桃红色精品国产亚洲av| 每晚都被弄得嗷嗷叫到高潮| 午夜老司机福利片| 亚洲少妇的诱惑av| 免费久久久久久久精品成人欧美视频| 亚洲第一青青草原| 91精品国产国语对白视频| 国产av精品麻豆| 韩国av一区二区三区四区| 亚洲精品中文字幕一二三四区| 国语自产精品视频在线第100页| 麻豆av在线久日| 国产一级毛片七仙女欲春2 | 亚洲av成人一区二区三| 久久人妻福利社区极品人妻图片| 国产乱人伦免费视频| 在线观看午夜福利视频| 十分钟在线观看高清视频www| 在线国产一区二区在线| 欧美激情极品国产一区二区三区| 最近最新中文字幕大全电影3 | 免费看十八禁软件| 多毛熟女@视频| 午夜福利欧美成人| 国产国语露脸激情在线看| 18禁观看日本| 真人一进一出gif抽搐免费| 亚洲电影在线观看av| 亚洲黑人精品在线| 一本大道久久a久久精品| 国产一区二区在线av高清观看| 久久久久国产一级毛片高清牌| 韩国精品一区二区三区| 后天国语完整版免费观看| 一区二区三区激情视频| 日本 欧美在线| 欧美乱码精品一区二区三区| 免费高清视频大片| 高清在线国产一区| 国产又色又爽无遮挡免费看| 精品一区二区三区四区五区乱码| 一进一出抽搐gif免费好疼| 香蕉国产在线看| 亚洲专区国产一区二区| 69av精品久久久久久| 亚洲最大成人中文| 制服丝袜大香蕉在线| 亚洲午夜理论影院| 成人免费观看视频高清| 久久青草综合色| 91精品三级在线观看| 九色国产91popny在线| 日韩欧美一区二区三区在线观看| av福利片在线| 最好的美女福利视频网| 亚洲精品国产色婷婷电影| 亚洲久久久国产精品| 韩国av一区二区三区四区| 久久亚洲真实| 999精品在线视频| 国产麻豆69| 日韩精品免费视频一区二区三区| 黄片播放在线免费| 精品久久久久久久毛片微露脸| 一边摸一边抽搐一进一出视频| 日韩大码丰满熟妇| 亚洲 国产 在线| 91麻豆av在线| www.999成人在线观看| 91九色精品人成在线观看| 日本 欧美在线| 老汉色∧v一级毛片| 久久久精品国产亚洲av高清涩受| 久久国产精品男人的天堂亚洲| 法律面前人人平等表现在哪些方面| 人人妻人人爽人人添夜夜欢视频| 亚洲午夜理论影院| 18美女黄网站色大片免费观看| 妹子高潮喷水视频| 久久香蕉激情| 国产亚洲精品第一综合不卡| 制服人妻中文乱码| videosex国产| 欧美精品亚洲一区二区| 国产国语露脸激情在线看| 久久狼人影院| 窝窝影院91人妻| 成人亚洲精品一区在线观看| 国产精品免费视频内射| 国产成人一区二区三区免费视频网站| 一区二区日韩欧美中文字幕| 精品久久久久久久人妻蜜臀av | 国内精品久久久久精免费| 黑丝袜美女国产一区| 变态另类成人亚洲欧美熟女 | 男女午夜视频在线观看| www.熟女人妻精品国产| 色综合站精品国产| 亚洲欧美精品综合一区二区三区| 一进一出抽搐动态| 国产xxxxx性猛交| 中文字幕久久专区| 黄频高清免费视频| 国产午夜福利久久久久久| 国产区一区二久久| 日韩欧美一区二区三区在线观看| 色综合婷婷激情| 18美女黄网站色大片免费观看| 男人舔女人的私密视频| 久久伊人香网站| 久久久久久大精品| 国产精品99久久99久久久不卡| 国产麻豆成人av免费视频| 搡老岳熟女国产| 亚洲情色 制服丝袜| 日韩欧美免费精品| 国产亚洲av高清不卡| av片东京热男人的天堂| 91麻豆精品激情在线观看国产| 亚洲国产日韩欧美精品在线观看 | 久久午夜综合久久蜜桃| 久久精品人人爽人人爽视色| 亚洲国产高清在线一区二区三 | 久久国产亚洲av麻豆专区| 亚洲男人天堂网一区| 亚洲色图综合在线观看| videosex国产| 搞女人的毛片| www.www免费av| 每晚都被弄得嗷嗷叫到高潮| 99在线视频只有这里精品首页| ponron亚洲| 夜夜夜夜夜久久久久| 国产xxxxx性猛交| 欧洲精品卡2卡3卡4卡5卡区| 午夜免费观看网址| 精品一区二区三区av网在线观看| 长腿黑丝高跟| 亚洲av成人一区二区三| 色播亚洲综合网| 亚洲精品在线美女| 在线观看午夜福利视频| 午夜精品国产一区二区电影| 精品一区二区三区四区五区乱码| 中文字幕人成人乱码亚洲影| 自拍欧美九色日韩亚洲蝌蚪91| 国产不卡一卡二| 欧美另类亚洲清纯唯美| avwww免费| 法律面前人人平等表现在哪些方面| 少妇被粗大的猛进出69影院| 91大片在线观看| 亚洲精品粉嫩美女一区| 色综合婷婷激情| 久久精品国产亚洲av高清一级| 9热在线视频观看99| 桃红色精品国产亚洲av| 999久久久精品免费观看国产| 在线观看日韩欧美| 久久久久久人人人人人| 97超级碰碰碰精品色视频在线观看| 亚洲国产毛片av蜜桃av| 亚洲欧美精品综合久久99| 精品人妻1区二区| 精品午夜福利视频在线观看一区| 午夜福利欧美成人| 国产日韩一区二区三区精品不卡| 亚洲精品美女久久av网站| 国产成人免费无遮挡视频| 精品久久久久久久毛片微露脸| 老司机深夜福利视频在线观看| 老鸭窝网址在线观看| 中文亚洲av片在线观看爽| 国产欧美日韩一区二区精品| 人人妻人人澡欧美一区二区 | 天天躁狠狠躁夜夜躁狠狠躁| 亚洲欧美精品综合久久99| 欧美日韩中文字幕国产精品一区二区三区 | 欧美日本亚洲视频在线播放| 亚洲在线自拍视频| 男女下面插进去视频免费观看| 国产精品久久久久久亚洲av鲁大| 久久久久精品国产欧美久久久| 国产精品av久久久久免费| 亚洲中文字幕一区二区三区有码在线看 | 黑人巨大精品欧美一区二区mp4| 最新在线观看一区二区三区| 日韩三级视频一区二区三区| 黄色成人免费大全| 亚洲免费av在线视频| 亚洲熟妇中文字幕五十中出|