共查询到20条相似文献,搜索用时 0 毫秒
1.
The time-dependent Navier–Stokes equations and the energy balance equation for an incompressible, constant property fluid in the Boussinesq approximation are solved by a least-squares finite element method based on a velocity–pressure–vorticity–temperature–heat-flux ( u –P–ω–T– q ) formulation discretized by backward finite differencing in time. The discretization scheme leads to the minimization of the residual in the l2-norm for each time step. Isoparametric bilinear quadrilateral elements and reduced integration are employed. Three examples, thermally driven cavity flow at Rayleigh numbers up to 106, lid-driven cavity flow at Reynolds numbers up to 104 and flow over a square obstacle at Reynolds number 200, are presented to validate the method. 相似文献
2.
A class of stable least-square finite element methods for non-linear hyperbolic problems is developed and some exploratory studies made. The methods are based on modifying the L2-norm of the. residual and a related approximation to the H1-norm of the residual. The effect of the additional terms in these residual functionals is to introduce a dissipative effect proportional to the solution gradient. This acts to stabilize the solution for non-linear hyperbolic problems which generate shocks. Numerical results for a one-dimensional nozzle and shock tube problem demonstrate the accuracy and stability of the method. Results are for an implicit scheme and calculations for linear, quadratic and cubic elements are given. 相似文献
3.
A velocity–vorticity formulation of the Navier–Stokes equations is presented as an alternative to the primitive variables approach. The velocity components and the vorticity are solved for in a fully coupled manner using a Newton method. No artificial viscosity is required in this formulation. The pressure is updated by a method allowing natural imposition of boundary conditions. Incompressible and subsonic results are presented for two-dimensional laminar internal flows up to high Reynolds numbers. 相似文献
4.
Stabilized finite element methods have been shown to yield robust, accurate numerical solutions to both the compressible and incompressible Navier–Stokes equations for laminar and turbulent flows. The present work focuses on the application of higher‐order, hierarchical basis functions to the incompressible Navier–Stokes equations using a stabilized finite element method. It is shown on a variety of problems that the most cost‐effective simulations (in terms of CPU time, memory, and disk storage) can be obtained using higher‐order basis functions when compared with the traditional linear basis. In addition, algorithms will be presented for the efficient implementation of these methods within the traditional finite element data structures. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
5.
Viscoelastic flows remain a demanding class of problems for approximate analysis, particularly at increasing Weissenberg numbers. Part of the difficulty stems from the convective behavior and in the treatment of the stress field as a primary unknown. This latter aspect has led to the use of higher-order piecewise approximations for the stress approximation spaces in recent finite element research. The computational complexity of the discretized problem is increased significantly by this approach but at present it appears the most viable technique for solving these problems. Motivated by recent success in treating mixed systems and convective problems, we formulate here a least squares finite element method for the viscoelastic flow problem. Numerical experiments are conducted to test the method and examine its strengths and limitations. Some difficulties and open issues are identified through the numerical experiments. We consider the use of high degree elements (p refinement) to improve performance and accuracy. 相似文献
6.
P. Wilders 《国际流体数值方法杂志》1988,8(8):977-987
The standard least-squares finite element method for the linearized Euler equations turns out to be inaccurate. This method is studied in detail for a system of composite type, obtained by transformation of the linearized Euler equations. The shortcomings of the method are clarified and an embedding method is constructed. It is shown numerically that this new method is O(h2)-accurate. 相似文献
7.
A finite element solution of the two-dimensional incompressible Navier–Stokes equations has been developed. The present method is a modified velocity correction approach. First an intermediate velocity is calculated, and then this is corrected by the pressure gradient which is the solution of a Poisson equation derived from the continuity equation. The novelty, in this paper, is that a second-order Runge–Kutta method for time integration has been used. Discretization in space is carried out by the Galerkin weighted residual method. The solution is in terms of primitive variables, which are approximated by polynomial basis functions defined on three-noded, isoparametric triangular elements. To demonstrate the present method, two examples are provided. Results from the first example, the driven cavity flow problem, are compared with previous works. Results from the second example, uniform flow past a cylinder, are compared with experimental data. 相似文献
8.
A numerical procedure for solving the time-dependent, incompressible Navier-Stokes equations is presented. The present method is based on a set of finite element equations of the primitive variable formulation, and a direct time integration method which has unique features in its formulation as well as in its evaluation of the contribution of external functions. Particular processes regarding the continuity conditions and the boundary conditions lead to a set of non-linear recurrence equations which represent evolution of the velocities and the pressures under the incompressibility constraint. An iteration process as to the non-linear convective terms is performed until the convergence is achieved in every integration step. Excessively artificial techniques are not introduced into the present solution procedure. Numerical examples with vortex shedding behind a rectangular cylinder are presented to illustrate the features of the proposed method. The calculated results are compared with experimental data and visualized flow fields in literature. 相似文献
9.
The paper deals with the numerical solution of fluid dynamics using the boundary‐domain integral method (BDIM). A velocity–vorticity formulation of the Navier–Stokes equations is adopted, where the kinematic equation is written in its parabolic form. Computational aspects of the numerical simulation of two‐dimensional flows is described in detail. In order to lower the computational cost, the subdomain technique is applied. A preconditioned Krylov subspace method (PKSM) is used for the solution of systems of linear equations. Level‐based fill‐in incomplete lower upper decomposition (ILU) preconditioners are developed and their performance is examined. Scaling of stopping criteria is applied to minimize the number of iterations for the PKSM. The effectiveness of the proposed method is tested on several benchmark test problems. Copyright © 1999 John Wiley & Sons, Ltd. 相似文献
10.
C. T. Shaw 《国际流体数值方法杂志》1991,12(1):81-92
In this paper, a segregated finite element scheme for the solution of the incompressible Navier-Stokes equations is proposed which is simpler in form than previously reported formulations. A pressure correction equation is derived from the momentum and continuity equations, and equal-order interpolation is used for both the velocity components and pressure. Algorithms such as this have been known to lead to checkerboard pressure oscillations; however, the pressure correction equation of this scheme should not produce these oscillations. The method is applied to several laminar flow situations, and details of the methods used to achieve converged solutions are given. 相似文献
11.
We note in this study that the Navier-Stokes equations, when expressed in streamfunction-vorticity form, can be approximated to fourth-order accuracy with stencils extending only over a 3 x 3 square of points. The key advantage of the new compact fourth-order scheme is that it allows direct iteration for low-to-medium Reynolds numbers. Numerical solutions are obtained for the model problem of the driven cavity and compared with solutions available in the literature. For Re ? 7500 point-SOR iteration is used and the convergence is fast. 相似文献
12.
In this paper a total linearization method is derived for solving steady viscous free boundary flow problems (including capillary effects) by the finite element method. It is shown that the influence of the geometrical unknown in the totally linearized weak formulation can be expressed in terms of boundary integrals. This means that the implementation of the method is simple. Numerical experiments show that the iterative method gives accurate results and converges very fast. 相似文献
13.
The steady Navier–Stokes equations in primitive variables are discretized in conservative form by a vertex-centred finite volume method Flux difference splitting is applied to the convective part to obtain an upwind discretization. The diffusive part is discretized in the central way. In its first-order formulation, flux difference splitting leads to a discretization of so-called vector positive type. This allows the use of classical relaxation methods in collective form. An alternating line Gauss–Seidel relaxation method is chosen here. This relaxation method is used as a smoother in a multigrid method. The components of this multigrid method are: full approximation scheme with F-cycles, bilinear prolongation, full weighting for residual restriction and injection of grid functions. Higher-order accuracy is achieved by the flux extrapolation method. In this approach the first-order convective fluxes are modified by adding second-order corrections involving flux limiting. Here the simple MinMod limiter is chosen. In the multigrid formulation the second-order discrete system is solved by defect correction. Computational results are shown for the well known GAMM backward-facing step problem and for a channel with a half-circular obstruction. 相似文献
14.
This paper is concerned with the implementation of Lagrange-Galerkin finite element methods for the Navier-Stokes equations. A scheme is developed to efficiently handle unstructed meshes with local refinement, using a quad-tree-based algorithm for the geometric search. Several difficulties that arise in the construction of the right-hand side are discussed in detail and some useful tricks are proposed. The resulting method is tested on the lid-driven square cavity and the vortex shedding behind a rectangular cylinder and is found to give satisfactory agreement with previous works. A detailed analysis of the effect of time discretization is included. 相似文献
15.
A methodology for the decomposition of the Crouzeix-Raviart finite element into six linear subelements is described. The resulting element is shown to satisfy the Brezzi-Babu?ka compatibility condition. The error bounds are also established. A comparison in accuracy between this and the standard Crouzeix-Raviart element is presented for driven cavity flows. Other results include the execution time for the DCR element and the Crouzeix-Raviart element along with both analytical and numerical integration. It is shown that the decomposed element results in shorter execution times with only marginal changes in accuracy. 相似文献
16.
Finite element schemes for hyperbolic systems are applied to the St. Venant equations for one-dimensional, unsteady, open channel flow. The comparative performances of the characteristic-dissipative-Galerkin, Taylor-Galerkin and least squares finite element schemes are assessed by means of linear Fourier analysis and solution of idealized non-linear wave propagation problems. Of particular interest is the behaviour of these schemes for the regressive wave component in both subcritical and supercritical flows. To assess the quality of the basic solution, the methods are compared without any additional artificial diffusion or shock-capturing formulations. The balanced treatment of both wave components in the characteristic-dissipative-Galerkin method is illustrated. Also, the method displays little sensitivity to parameters variations. The Taylor-Galerkin scheme provides good solutions, although oscillations due to wave dispersion and minimal diffusion of the regressive wave are displayed. Also, this method is somewhat sensitive to the time step increment. The least squares method is considered unsuitable for unsteady, open channel flow problems owing to its inability to propagate a regressive wave in a supercritical flow. 相似文献
17.
Rolf Stenberg 《国际流体数值方法杂志》1990,11(6):935-948
We give a self-contained presentation of our macroelement technique for verifying the stability of finite element discretizations of the Navier–Stokes equations in the velocity–pressure formulation. 相似文献
18.
The numerical simulation of complex flows demands efficient algorithms and fast computer platforms. The use of adaptive techniques permits adjusting the discretisation according to the analysis requirements, but creates variable computational loads that are difficult to manage in a parallel/vector program. This paper describes the approach adopted to implement an adaptive finite element incompressible Navier–Stokes solver on the Cray J90 machine. Performance measurements for the simulation of free and forced convection incompressible flows indicate that the techniques employed result in a fast parallel/vector code. Copyright © 1999 John Wiley & Sons, Ltd. 相似文献
19.
This paper presents a hybrid finite volume/finite element method for the incompressible generalized Newtonian fluid flow (Power-Law model). The collocated (i.e. non-staggered) arrangement of variables is used on the unstructured triangular grids, and a fractional step projection method is applied for the velocity-pressure coupling. The cell-centered finite volume method is employed to discretize the momentum equation and the vertex-based finite element for the pressure Poisson equation. The momentum interpolation method is used to suppress unphysical pressure wiggles. Numerical experiments demonstrate that the current hybrid scheme has second order accuracy in both space and time. Results on flows in the lid-driven cavity and between parallel walls for Newtonian and Power-Law models are also in good agreement with the published solutions. 相似文献
20.
Yang Zuosheng 《国际流体数值方法杂志》1998,28(3):565-568
A complete boundary integral formulation for incompressible Navier–Stokes equations with time discretization by operator splitting is developed using the fundamental solutions of the Helmholtz operator equation with different order. The numerical results for the lift and the drag hysteresis associated with a NACA0012 aerofoil oscillating in pitch show good agreement with available experimental data. © 1998 John Wiley & Sons, Ltd. 相似文献