C. Bozkaya  M. Tezer‐Sezgin 《国际流体数值方法杂志》2009,59(2):215-234

The two‐dimensional time‐dependent Navier–Stokes equations in terms of the vorticity and the stream function are solved numerically by using the coupling of the dual reciprocity boundary element method (DRBEM) in space with the differential quadrature method (DQM) in time. In DRBEM application, the convective and the time derivative terms in the vorticity transport equation are considered as the nonhomogeneity in the equation and are approximated by radial basis functions. The solution to the Poisson equation, which links stream function and vorticity with an initial vorticity guess, produces velocity components in turn for the solution to vorticity transport equation. The DRBEM formulation of the vorticity transport equation results in an initial value problem represented by a system of first‐order ordinary differential equations in time. When the DQM discretizes this system in time direction, we obtain a system of linear algebraic equations, which gives the solution vector for vorticity at any required time level. The procedure outlined here is also applied to solve the problem of two‐dimensional natural convection in a cavity by utilizing an iteration among the stream function, the vorticity transport and the energy equations as well. The test problems include two‐dimensional flow in a cavity when a force is present, the lid‐driven cavity and the natural convection in a square cavity. The numerical results are visualized in terms of stream function, vorticity and temperature contours for several values of Reynolds (Re) and Rayleigh (Ra) numbers. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

    

S. Turek  A. Ouazzi  J. Hron 《国际流体数值方法杂志》2009,59(4):387-403

We investigate a special technique called ‘pressure separation algorithm’ (PSepA) (see Applied Mathematics and Computation 2005; 165 :275–290 for an introduction) that is able to significantly improve the accuracy of incompressible flow simulations for problems with large pressure gradients. In our numerical studies with the computational fluid dynamics package FEATFLOW ( www.featflow.de ), we mainly focus on low‐order Stokes elements with nonconforming finite element approximations for the velocity and piecewise constant pressure functions. However, preliminary numerical tests show that this advantageous behavior can also be obtained for higher‐order discretizations, for instance, with Q₂/P₁ finite elements. We analyze the application of this simple, but very efficient, algorithm to several stationary and nonstationary benchmark configurations in 2D and 3D (driven cavity and flow around obstacles), and we also demonstrate its effect to spurious velocities in multiphase flow simulations (‘static bubble’ configuration) if combined with edge‐oriented, resp., interior penalty finite element method stabilization techniques. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

    

J. McKernan  G. Papadakis  J. F. Whidborne 《国际流体数值方法杂志》2009,59(8):907-925

This paper examines the performance of optimal linear quadratic state and output feedback controllers in stabilizing two‐dimensional perturbations in a plane Poiseuille flow. The synthesis of the controllers is based on a linearized model of the flow using a new set of interpolating polynomials in the wall‐normal direction, which automatically satisfy the homogeneous Dirichlet and Neumann boundary conditions at the walls and eliminate spurious eigenvalues. The controllers are implemented into a non‐linear Navier–Stokes solver, which is modified to compute the evolution of the flow perturbations. Two cases are examined, one with small initial disturbances that do not violate the linearity assumptions and the other with much larger disturbances that trigger the non‐linear convection terms. For the smallest disturbances, the solver accurately reproduced the results of the linear simulations of open‐ and closed‐loop systems. The simulations for the larger disturbances without control showed a rapid initial growth but the flow soon reached a saturated state in agreement with previous findings in the literature. The large initial growth is a consequence of the non‐normal nature of the system dynamics. The state feedback and output feedback controllers were able to reduce significantly the perturbation energy. For the larger disturbances, the energy calculated from the state variables is well below the energy evaluated by direct integration of the velocity field. This is probably due to the non‐linear terms transferring energy to harmonics of the considered wavenumber, which are not sensed by the linear controller. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

    

Quentin Araud  Pascal Finaud‐Guyot  Vincent Guinot  Robert Mosé  José Vazquez 《国际流体数值方法杂志》2012,70(12):1590-1604

Discontinuous Galerkin (DG) methods have shown promising results for solving the two‐dimensional shallow water equations. In this paper, the classical Runge–Kutta (RK) time discretisation is replaced by the eigenvector‐based reconstruction (EVR) that allows the second‐order time accuracy to be achieved within a single time‐stepping procedure. Moreover, the EVRDG approach yields stable solutions near drying and wetting fronts, whereas the classical RKDG approach yields instabilities. The proposed EVRDG technique is compared with the original RKDG approach on various test cases with analytical solutions. The EVRDG solutions are shown to be as accurate as those obtained with the RKDG scheme. Besides, the EVRDG scheme is 1.6 times faster than the RKDG method. Simulating dambreaks involving dry beds confirms that EVRDG scheme gives correct solutions, whereas the RKDG method yields instabilities. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

    

A. Montlaur  S. Fernandez‐Mendez  A. Huerta 《国际流体数值方法杂志》2012,70(5):603-626

The spatial discretization of unsteady incompressible Navier–Stokes equations is stated as a system of differential algebraic equations, corresponding to the conservation of momentum equation plus the constraint due to the incompressibility condition. Asymptotic stability of Runge–Kutta and Rosenbrock methods applied to the solution of the resulting index‐2 differential algebraic equations system is analyzed. A critical comparison of Rosenbrock, semi‐implicit, and fully implicit Runge–Kutta methods is performed in terms of order of convergence and stability. Numerical examples, considering a discontinuous Galerkin formulation with piecewise solenoidal approximation, demonstrate the applicability of the approaches and compare their performance with classical methods for incompressible flows. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

    

Bashar Albaalbaki  Roger E. Khayat 《国际流体数值方法杂志》2012,69(11):1762-1785

A comparison among three weakly nonlinear approaches for thermo‐gravitational instability in a Newtonian fluid layer heated from below is presented. First, the dynamical systems describing the time evolution of the problem from different weakly nonlinear approaches, namely, the Lorenz model, the amplitude equations and the perturbation expansion approaches are obtained. Next, the steady states and their stability, as well as the transient behaviour are obtained from each dynamical system. The similarity and difference among the three models are emphasized. The role of each of the nondimensional groups, the Rayleigh number and the Prandtl number is compared for the three models. The different approaches lead to similar behaviours when the Rayleigh number is just above its critical value and Prandtl number is high. However, only the dynamical system obtained from the amplitude equations is able to reflect the role of the Prandtl number. On the other hand, the amplitude equations and perturbation expansion techniques are not suitable for predicting the uniform oscillatory behaviour observed frequently in Rayleigh–Bénard convection. The novelty of the current work lies in studying the critical differences in the findings of the three popular approaches to investigate weakly nonlinear thermal convection for the first time. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

    

Z. Skalk 《ZAMM - Journal of Applied Mathematics and Mechanics / Zeitschrift für Angewandte Mathematik und Mechanik》2011,91(9):724-732

Let Ω ⊆ ℝ³ be a uniformly regular domain of the class C³ or Ω = ℝ³. Let A denote the Stokes operator and {E_λ; λ > 0} be the resolution of identity of A. We show as the main result of the paper that if w is a nonzero global weak solution to the Navier‐Stokes equations in Ω satisfying the strong energy inequality, then there exists a nonnegative finite number a = a(w) such that for every ε > 0 [lim_{t rightarrow infty} frac {||(E_{a+varepsilon}‐E_{a‐varepsilon}) w(t)||} {||w(t)||} = 1, ] where we put E_a‐ε = 0 if a‐ε < 0. Thus, every nonzero global weak solution satisfying the strong energy inequality exhibits large‐time energy concentration in a particular frequency. Moreover, the solutions with the exponentially decreasing energy are characterized by the positivity of a. In Appendix, we present some further results describing in detail the large‐time behavior of w.  相似文献   

    

S. J. Baxter  H. Power  K. A. Cliffe  S. Hibberd 《国际流体数值方法杂志》2010,62(5):530-564

Gravity‐driven Stokes flow down an inclined plane over and around multiple obstacles is considered. The flow problem is formulated in terms of a boundary integral equation and solved using the boundary element method. A Hermitian radial basis function (RBF) is used for the interpolation of the free surface, generation of the unit normal and curvature, and to prescribe the far‐field conditions. For flow over an obstacle, hemispheres are taken. For flow around an obstacle, circular cylinders are modelled and the contact angle condition on the obstacle/free surface intersection specified using the RBF formulation. Explicit profiles are produced for flow over and around two obstacles placed in various locations relative to one another. Interaction due to two obstacles is given by comparisons made with the profiles for flow over and around individual obstacles. In general, when the obstacles are separated by a sufficiently large distance the flow profiles are identical to a single obstacle analysis. For flow over and around two obstacles in‐line with the incident flow, effects of the governing parameters are examined, with variations in plane inclination angle, Bond number, obstacle size, and in the case of obstacles intersecting the free surface, static contact angle is considered. Finally flows over and around three obstacles are modelled. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

    

Jihn‐Sung Lai  Wen‐Dar Guo  Gwo‐Fong Lin  Yih‐Chi Tan 《国际流体数值方法杂志》2010,62(8):927-944

This study extends the upstream flux‐splitting finite‐volume (UFF) scheme to shallow water equations with source terms. Coupling the hydrostatic reconstruction method (HRM) with the UFF scheme achieves a resultant numerical scheme that adequately balances flux gradients and source terms. The proposed scheme is validated in three benchmark problems and applied to flood flows in the natural/irregular river with bridge pier obstructions. The results of the simulations are in satisfactory agreement with the available analytical solutions, experimental data and field measurements. Comparisons of the present results with those obtained by the surface gradient method (SGM) demonstrate the superior stability and higher accuracy of the HRM. The stability test results also show that the HRM requires less CPU time (up to 60%) than the SGM. The proposed well‐balanced UFF scheme is accurate, stable and efficient to solve flow problems involving irregular bed topography. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

    

Wenjing Yan  Yaling He  Yichen Ma 《国际流体数值方法杂志》2010,62(6):632-646

This paper is concerned with the problem of the shape reconstruction of two‐dimensional flows governed by the Navier–Stokes equations. Our objective is to derive a regularized Gauss–Newton method using the corresponding operator equation in which the unknown is the geometric domain. The theoretical foundation for the Gauss–Newton method is given by establishing the differentiability of the initial boundary value problem with respect to the boundary curve in the sense of a domain derivative. The numerical examples show that our theory is useful for practical purpose and the proposed algorithm is feasible. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献