首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到10条相似文献,搜索用时 109 毫秒
1.
The unsteady compressible flow equations are solved using a stabilized finite‐element formulation with C0 elements. In 2D, the performance of three‐noded linear and six‐noded quadratic triangular elements is compared. In 3D, the relative performance is evaluated for 6‐noded linear and 18‐noded quadratic wedge elements. Results are compared for the solutions to Euler, laminar, and turbulent flows at different Mach numbers for several flow problems. The finite‐element meshes considered for comparison have same location of nodes for the linear and quadratic interpolations. For the turbulent flow, the Spalart–Allmaras model is used for closure. It is found that the quadratic elements yield better performance than the linear elements. This is attributed to accurate representation of the stabilization terms that involve second‐order derivatives in the formulation. When these terms are dropped from the formulation with quadratic interpolation, the numerical results are similar to those obtained with linear interpolation. The absence of these terms result in added numerical diffusion that seems to give the effect of a relatively reduced Reynolds number. For the same location of nodes, the computations with the linear triangular and wedge elements are approximately 20% and 100% faster than those with quadratic triangular and wedge elements, respectively. However, if the same quadrature rule for numerical integration is used for both interpolations, the computations with quadratic elements are approximately 20% and 45% faster in 2D and 3D, respectively. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

2.
Least-squares stabilization stands out among the numerous approaches that have been proposed for relaxing resolution requirements of Galerkin computations for acoustics, by combining substantial improvement in performance with extremely simple implementation. The Galerkin/least-squares and Galerkin-gradient/least-squares methods are quite similar for structured meshes of linear finite elements. A series of numerical tests compares the two methods for several configurations with different kinds of boundary conditions employing structured and unstructured meshes. Various definitions of the resolution-dependent stability parameters are considered, along with different definitions of the mesh size upon which they depend.  相似文献   

3.
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.  相似文献   

4.
CFD modelling of ‘real‐life’ thermo‐fluid processes often requires solutions in complex three‐dimensional geometries, which can result in meshes containing aspects that are badly distorted. Cell‐centred finite volume methods (CC‐FV), typical of most commercial CFD tools, are computationally efficient, but can lead to convergence problems on meshes that feature cells with highly non‐orthogonal shapes. The control volume‐finite element method (CVFE) uses a vertex‐based approach and handles distorted meshes with relative ease, but is computationally expensive. A combined vertex‐based—cell‐centre technique (CFVM), detailed in this paper, allows solutions on distorted meshes where purely cell‐centred solutions procedures fail. The method utilizes the ability of the vertex‐based approach to resolve the flow field on a distorted mesh, enabling well established cell‐centred physical models to be employed in the solution of other transported quantities. The vertex‐based flow code is verified against a manufactured 3D solution and error norms are compared on meshes with various degrees of distortion. The CFVM method is validated with benchmark solutions for thermally driven flow and turbulent flow. Finally, the method is illustrated on three‐dimensional turbulent flow over an aircraft wing on a distorted mesh where purely cell‐centred techniques fail. The CFVM is relatively straightforward to embed within generic CC based CFD tools allowing it to be employed in a wide variety of processing applications. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

5.
A new mixed‐interpolation finite element method is presented for the two‐dimensional numerical simulation of incompressible magnetohydrodynamic (MHD) flows which involve convective heat transfer. The proposed method applies the nodal shape functions, which are locally defined in nine‐node elements, for the discretization of the Navier–Stokes and energy equations, and the vector shape functions, which are locally defined in four‐node elements, for the discretization of the electromagnetic field equations. The use of the vector shape functions allows the solenoidal condition on the magnetic field to be automatically satisfied in each four‐node element. In addition, efficient approximation procedures for the calculation of the integrals in the discretized equations are adopted to achieve high‐speed computation. With the use of the proposed numerical scheme, MHD channel flow and MHD natural convection under a constant applied magnetic field are simulated at different Hartmann numbers. The accuracy and robustness of the method are verified through these numerical tests in which both undistorted and distorted meshes are employed for comparison of numerical solutions. Furthermore, it is shown that the calculation speed for the proposed scheme is much higher compared with that for a conventional numerical integration scheme under the condition of almost the same memory consumption. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

6.
Curved geometries and the corresponding near-surface fields typically require a large number of linear computational elements. High-order numerical solvers have been primarily used with low-order meshes. There is a need for curved, high-order computational elements. Typical near-surface meshes consist of hexahedral and/or prismatic elements. The present work studies the employment of quadratic meshes that are relatively coarse for field simulations. Directionally quadratic high-order elements are proposed for the near-surface field regions. The quadratic meshes are compared with the conventional low-order ones in terms of accuracy and efficiency. The cases considered include closed surface volume calculations, as well as computation of gradients of several analytic fields. A special method of adaptive local quadratic meshes is proposed and evaluated. Truncation error analysis for quadratic grids yields comparison with the conventional linear hexahedral/prismatic meshes, which are subject to typical distortions such as stretching, skewness, and torsion.  相似文献   

7.
Discontinuous Galerkin (DG) methods allow high‐order flow solutions on unstructured or locally refined meshes by increasing the polynomial degree and using curved instead of straight‐sided elements. DG discretizations with higher polynomial degrees must, however, be stabilized in the vicinity of discontinuities of flow solutions such as shocks. In this article, we device a consistent shock‐capturing method for the Reynolds‐averaged Navier–Stokes and kω turbulence model equations based on an artificial viscosity term that depends on element residual terms. Furthermore, the DG method is combined with a residual‐based adaptation algorithm that targets at resolving all flow features. The higher‐order and adaptive DG method is applied to a fully turbulent transonic flow around the second Vortex Flow Experiment (VFE‐2) configuration with a good resolution of the vortex system.Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

8.
A discontinuous Galerkin method for the solution of the immiscible and incompressible two‐phase flow problem based on the nonsymmetric interior penalty method is presented. Therefore, the incompressible Navier–Stokes equation is solved for a domain decomposed into two subdomains with different values of viscosity and density as well as a singular surface tension force. On the basis of a piecewise linear approximation of the interface, meshes for both phases are cut out of a structured mesh. The discontinuous finite elements are defined on the resulting Cartesian cut‐cell mesh and may therefore approximate the discontinuities of the pressure and the velocity derivatives across the interface with high accuracy. As the mesh resolves the interface, regularization of the density and viscosity jumps across the interface is not required. This preserves the local conservation property of the velocity field even in the vicinity of the interface and constitutes a significant advantage compared with standard methods that require regularization of these discontinuities and cannot represent the jumps and kinks in pressure and velocity. A powerful subtessellation algorithm is incorporated to allow the usage of standard time integrators (such as Crank–Nicholson) on the time‐dependent mesh. The presented discretization is applicable to both the two‐dimensional and three‐dimensional cases. The performance of our approach is demonstrated by application to a two‐dimensional benchmark problem, allowing for a thorough comparison with other numerical methods. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

9.
The acoustic perturbation equations (APE) are suitable to predict aerodynamic noise in the presence of a non‐uniform mean flow. As for any hybrid computational aeroacoustics approach, a first computational fluid dynamics simulation is carried out from which the mean flow characteristics and acoustic sources are obtained. In a second step, the APE are solved to get the acoustic pressure and particle velocity fields. However, resorting to the finite element method (FEM) for that purpose is not straightforward. Whereas mixed finite elements satisfying an appropriate inf–sup compatibility condition can be built in the case of no mean flow, that is, for the standard wave equation in mixed form, these are difficult to implement and their good performance is yet to be checked for more complex wave operators. As a consequence, strong simplifying assumptions are usually considered when solving the APE with FEM. It is possible to avoid them by resorting to stabilized formulations. In this work, a residual‐based stabilized FEM is presented for the APE at low Mach numbers, which allows one to deal with the APE convective and reaction terms in its full extent. The key of the approach resides in the design of the matrix of stabilization parameters. The performance of the formulation and the contributions of the different terms in the equations are tested for an acoustic pulse propagating in sheared‐solenoidal mean flow, and for the aeolian tone generated by flow past a two‐dimensional cylinder. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

10.
In this paper we study the stability and performance of the quadrilateral finite element ??1–??0 (bili‐ near/constant) for the Stokes equations. We set up a framework to show the stability of the element for a wide range of meshes with macroelement patches. We apply the new theory to show the stability of ??1–??0 elements on some previously studied meshes and on some newly suggested meshes. Nevertheless such earlier and newly suggested meshes are not effective in practice, compared to the traditional unstable meshes for the ??1–??0 element. The new theory leads naturally to a general idea in treating instability of square ??1–??0 elements by the local stabilization on macroelement patches of larger, but fixed sizes. The good performance of the traditional ??1–??0 square elements with filtering can be kept in some cases after the local stabilization. Some numerical tests are provided to support the theory and to show the performance of stabilized ??1–??0 elements. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

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

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