首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到10条相似文献,搜索用时 109 毫秒
1.
In this paper, the domain‐free discretization method (DFD) is extended to simulate the three‐dimensional compressible inviscid flows governed by Euler equations. The discretization strategy of DFD is that the discrete form of governing equations at an interior point may involve some points outside the solution domain. The functional values at the exterior‐dependent points are updated at each time step by extrapolation along the wall normal direction in conjunction with the wall boundary conditions and the simplified momentum equation in the vicinity of the wall. Spatial discretization is achieved with the help of the finite element Galerkin approximation. The concept of ‘osculating plane’ is adopted, with which the local DFD can be easily implemented for the three‐dimensional case. Geometry‐adaptive tetrahedral mesh is employed for three‐dimensional calculations. Finally, we validate the DFD method for three‐dimensional compressible inviscid flow simulations by computing transonic flows over the ONERA M6 wing. Comparison with the reference experimental data and numerical results on boundary‐conforming grid was displayed and the results show that the present DFD results compare very well with the reference data. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

2.
This paper presents a local domain‐free discretization (DFD) method for the simulation of unsteady flows over moving bodies governed by the incompressible Navier–Stokes equations. The discretization strategy of DFD is that the discrete form of partial differential equations at an interior point may involve some points outside the solution domain. All the mesh points are classified as interior points, exterior dependent points and exterior independent points. The functional values at the exterior dependent points are updated at each time step by the approximate form of solution near the boundary. When the body is moving, only the status of points is changed and the mesh can stay fixed. The issue of ‘freshly cleared nodes/cells’ encountered in usual sharp interface methods does not pose any particular difficulty in the presented method. The Galerkin finite‐element approximation is used for spatial discretization, and the discrete equations are integrated in time via a dual‐time‐stepping scheme based on artificial compressibility. In order to validate the present method for moving‐boundary flow problems, two groups of flow phenomena have been simulated: (1) flows over a fixed circular cylinder, a harmonic in‐line oscillating cylinder in fluid at rest and a transversely oscillating cylinder in uniform flow; (2) flows over a pure pitching airfoil, a heaving–pitching airfoil and a deforming airfoil. The predictions show good agreement with the published numerical results or experimental data. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

3.
This paper is the first endeavour to present the local domain‐free discretization (DFD) method for the solution of compressible Navier–Stokes/Euler equations in conservative form. The discretization strategy of DFD is that for any complex geometry, there is no need to introduce coordinate transformation and the discrete form of governing equations at an interior point may involve some points outside the solution domain. The functional values at the exterior dependent points are updated at each time step to impose the wall boundary condition by the approximate form of solution near the boundary. Some points inside the solution domain are constructed for the approximate form of solution, and the flow variables at constructed points are evaluated by the linear interpolation on triangles. The numerical schemes used in DFD are the finite element Galerkin method for spatial discretization and the dual‐time scheme for temporal discretization. Some numerical results of compressible flows over fixed and moving bodies are presented to validate the local DFD method. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

4.
This paper details an approach to modelling gas–solid fluidized beds using the two‐fluid granular temperature model. Details concerning the difficulties associated with the boundary conditions, particularly for curved boundaries, are described along with a novel means of obtaining the internal stress of the solid‐phase, in part, by solving an implicit equation. This results in a scheme that is stable even when the solid volume fraction is close to maximum packing. A transient, mixed finite element discretization is used to solve the multi‐phase equations with a discontinuous finite element representation of the granular temperature and continuity equations. A new solution method is proposed to solve the coupled momentum and continuity equations based on Arnoldi iteration. Two fluidized beds are modelled, one in the bubbling regime and the other in the slugging regime. These simulations are compared with experiments. Copyright © 2001 John Wiley & Sons, Ltd.  相似文献   

5.
The paper deals with the use of the discontinuous Galerkin finite element method (DGFEM) for the numerical solution of viscous compressible flows. We start with a scalar convection–diffusion equation and present a discretization with the aid of the non‐symmetric variant of DGFEM with interior and boundary penalty terms. We also mention some theoretical results. Then we extend the scheme to the system of the Navier–Stokes equations and discuss the treatment of stabilization terms. Several numerical examples are presented. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

6.
The deformations of both electrically charged and uncharged incompressible axisymmetric Newtonian viscous liquid droplets acting under the effects of surface tension (but not gravity) are studied using a non‐conforming, discontinuous Galerkin finite element procedure with moving meshes. The full Navier–Stokes equations are discretized and solved in an Eulerian manner with a simple predictor–corrector Lagrangian updating of the free boundary location coincident with the droplet surface at each solution iteration. By using linear Crouzeix–Raviart basis functions for the velocity and piecewise constant pressures, results are presented both for the simple (oscillatory) relaxation of elongated electrically charged/uncharged droplets to a sphere and for the deformation to steady state/Coulomb explosions of initially slightly oblate/prolate spheroidal droplets charged beyond/to the Rayleigh limit. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

7.
Using a non‐conforming C0‐interior penalty method and the Galerkin least‐square approach, we develop a continuous–discontinuous Galerkin finite element method for discretizing fourth‐order incompressible flow problems. The formulation is weakly coercive for spaces that fail to satisfy the inf‐sup condition and consider discontinuous basis functions for the pressure field. We consider the results of a stability analysis through a lemma which indicates that there exists an optimal or quasi‐optimal least‐square stability parameter that depends on the polynomial degree used to interpolate the velocity and pressure fields, and on the geometry of the finite element in the mesh. We provide several numerical experiments illustrating such dependence, as well as the robustness of the method to deal with arbitrary basis functions for velocity and pressure, and the ability to stabilize large pressure gradients. We believe the results provided in this paper contribute for establishing a paradigm for future studies of the parameter of the Galerkin least square method for second‐gradient theory of incompressible flow problems. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

8.
A new finite element method is developed to simulate time‐dependent viscoelastic shear‐thinning flows characterized by the generalized Oldroyd‐B model. The focus of the algorithm is improved stability through a free‐energy dissipative scheme by using low‐order piecewise‐constant finite element approximations for stress. The algorithm is further modified by incorporating a pressure‐projection method, a DG‐upwinding scheme, a symmetric interior penalty DG method to solve the elliptic pressure‐update equation and a geometric multigrid preconditioner. The improved stability and cost to accuracy is compared when using higher order discontinuous bilinear approximation, where in addition, we consider the influence of a slope limiter for these elements. The algorithm is applied to the 2D start‐up‐driven cavity problem, and the stability of the free energy is illustrated and compared between element choices. An application of the model to modelling blood in small arterioles and channels is considered by simulating pulsatile blood flow through a stenotic arteriole. The individual influences of viscoelasticity and shear‐thinning within the generalized Oldroyd‐B model are investigated by comparing results to the Newtonian, generalized Newtonian and Oldroyd‐B models. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

9.
A numerical method for the simulation of compressible two‐phase flows is presented in this paper. The sharp‐interface approach consists of several components: a discontinuous Galerkin solver for compressible fluid flow, a level‐set tracking algorithm to follow the movement of the interface and a coupling of both by a ghost‐fluid approach with use of a local Riemann solver at the interface. There are several novel techniques used: the discontinuous Galerkin scheme allows locally a subcell resolution to enhance the interface resolution and an interior finite volume Total Variation Diminishing (TVD) approximation at the interface. The level‐set equation is solved by the same discontinuous Galerkin scheme. To obtain a very good approximation of the interface curvature, the accuracy of the level‐set field is improved and smoothed by an additional PNPM‐reconstruction. The capabilities of the method for the simulation of compressible two‐phase flow are demonstrated for a droplet at equilibrium, an oscillating ellipsoidal droplet, and a shock‐droplet interaction problem at Mach 3. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

10.
We present a spectral‐element discontinuous Galerkin thermal lattice Boltzmann method for fluid–solid conjugate heat transfer applications. Using the discrete Boltzmann equation, we propose a numerical scheme for conjugate heat transfer applications on unstructured, non‐uniform grids. We employ a double‐distribution thermal lattice Boltzmann model to resolve flows with variable Prandtl (Pr) number. Based upon its finite element heritage, the spectral‐element discontinuous Galerkin discretization provides an effective means to model and investigate thermal transport in applications with complex geometries. Our solutions are represented by the tensor product basis of the one‐dimensional Legendre–Lagrange interpolation polynomials. A high‐order discretization is employed on body‐conforming hexahedral elements with Gauss–Lobatto–Legendre quadrature nodes. Thermal and hydrodynamic bounce‐back boundary conditions are imposed via the numerical flux formulation that arises because of the discontinuous Galerkin approach. As a result, our scheme does not require tedious extrapolation at the boundaries, which may cause loss of mass conservation. We compare solutions of the proposed scheme with an analytical solution for a solid–solid conjugate heat transfer problem in a 2D annulus and illustrate the capture of temperature continuities across interfaces for conductivity ratio γ > 1. We also investigate the effect of Reynolds (Re) and Grashof (Gr) number on the conjugate heat transfer between a heat‐generating solid and a surrounding fluid. Steady‐state results are presented for Re = 5?40 and Gr = 105?106. In each case, we discuss the effect of Re and Gr on the heat flux (i.e. Nusselt number Nu) at the fluid–solid interface. Our results are validated against previous studies that employ finite‐difference and continuous spectral‐element methods to solve the Navier–Stokes equations. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号