Coupled fluid–solid interaction under shock wave loading

This paper considers the treatment of fluid–solid interaction problems under shock wave loading, where the solid experiences large bulk Lagrangian displacements. This work addresses the issues associated with using a level set as a generalized interface for fluid–solid coupling where the fluid–solid interface is embedded in an unstructured fluid grid. We outline the formulation used for the edge‐based unstructured‐grid Euler solver. The identification of the fluid–solid interface on the unstructured fluid mesh uses a super‐sampled ??² projection technique, which in conjunction with a Lagrangian interface position, permits fast identification of the interface and the concomitant imposition of boundary conditions. The use of a narrow‐band approach for the identification of the wetted interface is presented with the details of the construction of interface conditions. A series of two and three‐dimensional shock‐body computations are presented to demonstrate the validity of the current approach on problems with static and dynamic interfaces, including projectile/shock interaction simulations. Copyright © 2010 John Wiley & Sons, Ltd.     

Discontinuous Galerkin methods for the Navier–Stokes equations using solenoidal approximations

An interior penalty method and a compact discontinuous Galerkin method are proposed and compared for the solution of the steady incompressible Navier–Stokes equations. Both compact formulations can be easily applied using high‐order piecewise divergence‐free approximations, leading to two uncoupled problems: one associated with velocity and hybrid pressure, and the other one only concerned with the computation of pressures in the elements interior. Numerical examples compare the efficiency and the accuracy of both proposed methods. Copyright © 2009 John Wiley & Sons, Ltd.     

Nonlinear Normal Modes of a Parametrically Excited Cantilever Beam

We investigate theoretically thenonlinear normal modes of a vertical cantilever beam excited by aprincipal parametric resonance. We apply directly the method ofmultiple scales to the governing nonlinear nonautonomousintegral-partial-differential equation and associated boundary conditions.In the absence of damping, it is shown that the system has nonlinear normal modes, as defined by Rosenberg, even in the presence of the parametric excitation.We calculate the spatial correction to the linear mode shapedue to the effects of the inertia and curvature nonlinearities andthe parametric excitation. We compare the result obtained withthe direct approach with that obtained using a single-mode Galerkindiscretization.The deviation between the two predictions increases as the oscillationamplitude increases.     

Are upwind techniques in multi-phase flow models necessary?

Two alternatives of primary variables are compared for two-phase flow in heterogeneous media by solving fully established benchmarks. The first combination utilizes pressure of the wetting fluid and saturation of the non-wetting fluid as primary variables, while the second employs capillary pressure of the wetting fluid and pressure of the non-wetting fluid. While the standard Galerkin finite element method (SGFEM) is known to fail in the physical reproduction of two-phase flow in heterogeneous media (unless employing a fully upwind correction), the second scheme with capillary pressure as a primary variable without applying an upwind technique produces correct physical fluid behaviour in heterogeneous media, as observed from experiments.     

An eigenvector‐based linear reconstruction approach for time stepping in discontinuous Galerkin scheme used to solve shallow water equations

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.     

On the stability and dissipation of wall boundary conditions for compressible flows

Characteristic formulations for boundary conditions have demonstrated their effectiveness to handle inlets and outlets, especially to avoid acoustic wave reflections. At walls, however, most authors use simple Dirichlet or Neumann boundary conditions, where the normal velocity (or pressure gradient) is set to zero. This paper demonstrates that there are significant differences between characteristic and Dirichlet methods at a wall and that simulations are more stable when using walls modelled with a characteristic wave decomposition. The derivation of characteristic methods yields an additional boundary term in the continuity equation, which explains their increased stability. This term also allows to handle the two acoustic waves going towards and away from the wall in a consistent manner. Those observations are confirmed by stability matrix analysis and one‐ and two‐dimensional simulations of acoustic modes in cavities. Copyright © 2009 John Wiley & Sons, Ltd.     

Bifurcations of a Thin Plate-Strip Excited Transversally and Axially

The complex vibrations and bifurcations of plates modeled as systemswith infinite degrees-of-freedom are considered. Both theBubnov–Galerkin with high-order approximations and finite differencemethods with approximation O(h ⁴)are applied. In addition, the calculation ofthe Lyapunov exponents of the system is performed, and the results arecompared to those derived by Bennetin's method. Some examples of newnonlinear phenomena exhibited by the considered systems are reported.     

Efficient computation of mean drag for the subcritical flow past a circular cylinder using general Galerkin G2

General Galerkin (G2) is a new computational method for turbulent flow, where a stabilized Galerkin finite element method is used to compute approximate weak solutions to the Navier–Stokes equations directly, without any filtering of the equations as in a standard approach to turbulence simulation, such as large eddy simulation, and thus no Reynolds stresses are introduced, which need modelling. In this paper, G2 is used to compute the drag coefficient c_D for the flow past a circular cylinder at Reynolds number Re=3900, for which the flow is turbulent. It is found that it is possible to approximate c_D to an accuracy of a few percent, corresponding to the accuracy in experimental results for this problem, using less than 10⁵ mesh points, which makes the simulations possible using a standard PC. The mesh is adaptively refined until a stopping criterion is reached with respect to the error in a chosen output of interest, which in this paper is c_D. Both the stopping criterion and the mesh‐refinement strategy are based on a posteriori error estimates, in the form of a space–time integral of residuals times derivatives of the solution of a dual problem, linearized at the approximate solution, and with data coupling to the output of interest. Copyright © 2008 John Wiley & Sons, Ltd.