Extension of a discontinuous Galerkin finite element method to viscous rotor flow simulations

Heavy vibratory loading of rotorcraft is relevant for many operational aspects of helicopters, such as the structural life span of (rotating) components, operational availability, the pilot's comfort, and the effectiveness of weapon targeting systems. A precise understanding of the source of these vibrational loads has important consequences in these application areas. Moreover, in order to exploit the full potential offered by new vibration reduction technologies, current analysis tools need to be improved with respect to the level of physical modeling of flow phenomena which contribute to the vibratory loads. In this paper, a computational fluid dynamics tool for rotorcraft simulations based on first-principles flow physics is extended to enable the simulation of viscous flows. Viscous effects play a significant role in the aerodynamics of helicopter rotors in high-speed flight. The new model is applied to three-dimensional vortex flow and laminar dynamic stall. The applications clearly demonstrate the capability of the new model to perform on deforming and adaptive meshes. This capability is essential for rotor simulations to accomodate the blade motions and to enhance vortex resolution

Discontinuous Galerkin approximations in computational mechanics: hybridization, exact geometry and degree adaptivity

Discontinuous Galerkin (DG) discretizations with exact representation of the geometry and local polynomial degree adaptivity are revisited. Hybridization techniques are employed to reduce the computational cost of DG approximations and devise the hybridizable discontinuous Galerkin (HDG) method. Exact geometry described by non-uniform rational B-splines (NURBS) is integrated into HDG using the framework of the NURBS-enhanced finite element method (NEFEM). Moreover, optimal convergence and superconvergence properties of HDG-Voigt formulation in presence of symmetric second-order tensors are exploited to construct inexpensive error indicators and drive degree adaptive procedures. Applications involving the numerical simulation of problems in electrostatics, linear elasticity and incompressible viscous flows are presented. Moreover, this is done for both high-order HDG approximations and the lowest-order framework of face-centered finite volumes (FCFV).Peer ReviewedPostprint (author's final draft

Spectral/hp element methods: recent developments, applications, and perspectives

The spectral/hp element method combines the geometric flexibility of the classical h-type finite element technique with the desirable numerical properties of spectral methods, employing high-degree piecewise polynomial basis functions on coarse finite element-type meshes. The spatial approximation is based upon orthogonal polynomials, such as Legendre or Chebychev polynomials, modified to accommodate C0-continuous expansions. Computationally and theoretically, by increasing the polynomial order p, high-precision solutions and fast convergence can be obtained and, in particular, under certain regularity assumptions an exponential reduction in approximation error between numerical and exact solutions can be achieved. This method has now been applied in many simulation studies of both fundamental and practical engineering flows. This paper briefly describes the formulation of the spectral/hp element method and provides an overview of its application to computational fluid dynamics. In particular, it focuses on the use the spectral/hp element method in transitional flows and ocean engineering. Finally, some of the major challenges to be overcome in order to use the spectral/hp element method in more complex science and engineering applications are discussed

Theoretical and numerical comparison of hyperelastic and hypoelastic formulations for Eulerian non-linear elastoplasticity

The aim of this paper is to compare a hyperelastic with a hypoelastic model describing the Eulerian dynamics of solids in the context of non-linear elastoplastic deformations. Specifically, we consider the well-known hypoelastic Wilkins model, which is compared against a hyperelastic model based on the work of Godunov and Romenski. First, we discuss some general conceptual differences between the two approaches. Second, a detailed study of both models is proposed, where differences are made evident at the aid of deriving a hypoelastic-type model corresponding to the hyperelastic model and a particular equation of state used in this paper. Third, using the same high order ADER Finite Volume and Discontinuous Galerkin methods on fixed and moving unstructured meshes for both models, a wide range of numerical benchmark test problems has been solved. The numerical solutions obtained for the two different models are directly compared with each other. For small elastic deformations, the two models produce very similar solutions that are close to each other. However, if large elastic or elastoplastic deformations occur, the solutions present larger differences.Comment: 14 figure

Extension of the discontinuous Galerkin finite element method to viscous rotor flow simulations

Heavy vibratory loading of rotorcraft is relevant for many operational aspects of helicopters, such as the structural life span of (rotating) components, op- erational availability, the pilot’s comfort, and the ef- fectiveness of weapon targeting systems. A precise understanding of the source of these vibrational loads has important consequences in these application ar- eas. Moreover, in order to exploit the full poten- tial offered by new vibration reduction technologies, current analysis tools need to be improved with re- spect to the level of physical modeling of flow phe- nomena which contribute to the vibratory loads. In this paper, a computational fluid dynamics tool for rotorcraft simulations based on first-principles flow physics is extended to enable the simulation of vis- cous flows. Viscous effects play a significant role in the aerodynamics of helicopter rotors in high-speed flight. The new model is applied to three-dimensional vortex flow and laminar dynamic stall. The applica- tions clearly demonstrate the capability of the new model to perform on deforming and adaptive meshes. This capability is essential for rotor simulations to accomodate the blade motions and to enhance vor- tex resolution

SOLID-SHELL FINITE ELEMENT MODELS FOR EXPLICIT SIMULATIONS OF CRACK PROPAGATION IN THIN STRUCTURES

Crack propagation in thin shell structures due to cutting is conveniently simulated using explicit finite element approaches, in view of the high nonlinearity of the problem. Solidshell elements are usually preferred for the discretization in the presence of complex material behavior and degradation phenomena such as delamination, since they allow for a correct representation of the thickness geometry. However, in solid-shell elements the small thickness leads to a very high maximum eigenfrequency, which imply very small stable time-steps. A new selective mass scaling technique is proposed to increase the time-step size without affecting accuracy. New ”directional” cohesive interface elements are used in conjunction with selective mass scaling to account for the interaction with a sharp blade in cutting processes of thin ductile shells