Automated parameters for troubled-cell indicators using outlier detection

In Vuik and Ryan (2014) we studied the use of troubled-cell indicators for discontinuity detection in nonlinear hyperbolic partial differential equations and introduced a new multiwavelet technique to detect troubled cells. We found that these methods perform well as long as a suitable, problem-dependent parameter is chosen. This parameter is used in a threshold which decides whether or not to detect an element as a troubled cell. Until now, these parameters could not be chosen automatically. The choice of the parameter has impact on the approximation: it determines the strictness of the troubled-cell indicator. An inappropriate choice of the parameter will result in detection (and limiting) of too few or too many elements. The optimal parameter is chosen such that the minimal number of troubled cells is detected and the resulting approximation is free of spurious oscillations. In this paper we will see that for each troubled-cell indicator the sudden increase or decrease of the indicator value with respect to the neighboring values is important for detection. Indication basically reduces to detecting the outliers of a vector (one dimension) or matrix (two dimensions). This is done using Tukey's boxplot approach to detect which coefficients in a vector are straying far beyond others (Tukey, 1977). We provide an algorithm that can be applied to various troubled-cell indication variables. Using this technique the problem-dependent parameter that the original indicator requires is no longer necessary as the parameter will be chosen automatically

Historical roots of Agile methods: where did “Agile thinking” come from?

The appearance of Agile methods has been the most noticeable change to software process thinking in the last fifteen years [16], but in fact many of the “Agile ideas” have been around since 70’s or even before. Many studies and reviews have been conducted about Agile methods which ascribe their emergence as a reaction against traditional methods. In this paper, we argue that although Agile methods are new as a whole, they have strong roots in the history of software engineering. In addition to the iterative and incremental approaches that have been in use since 1957 [21], people who criticised the traditional methods suggested alternative approaches which were actually Agile ideas such as the response to change, customer involvement, and working software over documentation. The authors of this paper believe that education about the history of Agile thinking will help to develop better understanding as well as promoting the use of Agile methods. We therefore present and discuss the reasons behind the development and introduction of Agile methods, as a reaction to traditional methods, as a result of people's experience, and in particular focusing on reusing ideas from histor

Discontinuous high-order finite-volume/finite-element method for inviscid compressible flows

The discontinuous, hybrid control-volume/finite-element method merges the desirable conservative properties and intuitive physical formulation of the finite-volume technique, with the capability of local arbitrary high-order accuracy distinctive of the discontinuous finite-element method. This relatively novel scheme has been previously applied to the solution of advection-diffusion problems and the shallow-water equations, and is in the present work extended to the Euler equations. The derivation of the method is presented in the general multi-dimensional case, and selected numerical problems are solved in the one- and two-dimensional case

Tutorial on Hybridizable Discontinous Galerkin (HDG) for second-order elliptic problems

The HDG is a new class of discontinuous Galerkin (DG) methods that shares favorable properties with classical mixed methods such as the well known Raviart-Thomas methods. In particular, HDG provides optimal convergence of both the primal and the dual variables of the mixed formulation. This property enables the construction of superconvergent solutions, contrary to other popular DG methods. In addition, its reduced computational cost, compared to other DG methods, has made HDG an attractive alternative for solving problems governed by partial differential equations. A tutorial on HDG for the numerical solution of second-order elliptic problems is presented. Particular emphasis is placed on providing all the necessary details for the implementation of HDG methods.Peer ReviewedPreprin

Gender, war and militarism: making and questioning the links

The gender dynamics of militarism have traditionally been seen as straightforward, given the cultural mythologies of warfare and the disciplining of ‘masculinity’ that occurs in the training and use of men's capacity for violence in the armed services. However, women's relation to both war and peace has been varied and complex. It is women who have often been most prominent in working for peace, although there are no necessary links between women and opposition to militarism. In addition, more women than ever are serving in many of today's armies, with feminists rather uncertain on how to relate to this phenomenon. In this article, I explore some of the complexities of applying gender analyses to militarism and peace work in sites of conflict today, looking most closely at the Israeli feminist group, New Profile, and their insistence upon the costs of the militarized nature of Israeli society. They expose the very permeable boundaries between the military and civil society, as violence seeps into the fears and practices of everyday life in Israel. I place their work in the context of broader feminist analysis offered by researchers such as Cynthia Enloe and Cynthia Cockburn, who have for decades been writing about the ‘masculinist’ postures and practices of warfare, as well as the situation of women caught up in them. Finally, I suggest that rethinking the gendered nature of warfare must also encompass the costs of war to men, whose fundamental vulnerability to psychological abuse and physical injury is often downplayed, whether in mainstream accounts of warfare or in more specific gender analysis. Feminists need to pay careful attention to masculinity and its fragmentations in addressing the topic of gender, war and militarism

An immersed discontinuous Galerkin method for compressible Navier-Stokes equations on unstructured meshes

We introduce an immersed high-order discontinuous Galerkin method for solving the compressible Navier-Stokes equations on non-boundary-fitted meshes. The flow equations are discretised with a mixed discontinuous Galerkin formulation and are advanced in time with an explicit time marching scheme. The discretisation meshes may contain simplicial (triangular or tetrahedral) elements of different sizes and need not be structured. On the discretisation mesh the fluid domain boundary is represented with an implicit signed distance function. The cut-elements partially covered by the solid domain are integrated after tessellation with the marching triangle or tetrahedra algorithms. Two alternative techniques are introduced to overcome the excessive stable time step restrictions imposed by cut-elements. In the first approach the cut-basis functions are replaced with the extrapolated basis functions from the nearest largest element. In the second approach the cut-basis functions are simply scaled proportionally to the fraction of the cut-element covered by the solid. To achieve high-order accuracy additional nodes are introduced on the element faces abutting the solid boundary. Subsequently, the faces are curved by projecting the introduced nodes to the boundary. The proposed approach is verified and validated with several two- and three-dimensional subsonic and hypersonic low Reynolds number flow applications, including the flow over a cylinder, a space capsule and an aerospace vehicle

Local Error Analysis of Discontinuous Galerkin Methods for Advection-Dominated Elliptic Linear-Quadratic Optimal Control Problems

This paper analyzes the local properties of the symmetric interior penalty upwind discontinuous Galerkin (SIPG) method for the numerical solution of optimal control problems governed by linear reaction-advection-diffusion equations with distributed controls. The theoretical and numerical results presented in this paper show that for advection-dominated problems the convergence properties of the SIPG discretization can be superior to the convergence properties of stabilized finite element discretizations such as the streamline upwind Petrov Galerkin (SUPG) method. For example, we show that for a small diffusion parameter the SIPG method is optimal in the interior of the domain. This is in sharp contrast to SUPG discretizations, for which it is known that the existence of boundary layers can pollute the numerical solution of optimal control problems everywhere even into domains where the solution is smooth and, as a consequence, in general reduces the convergence rates to only first order. In order to prove the nice convergence properties of the SIPG discretization for optimal control problems, we first improve local error estimates of the SIPG discretization for single advection-dominated equations by showing that the size of the numerical boundary layer is controlled not by the mesh size but rather by the size of the diffusion parameter. As a result, for small diffusion, the boundary layers are too “weak” to pollute the SIPG solution into domains of smoothness in optimal control problems. This favorable property of the SIPG method is due to the weak treatment of boundary conditions, which is natural for discontinuous Galerkin methods, while for SUPG methods strong imposition of boundary conditions is more conventional. The importance of the weak treatment of boundary conditions for the solution of advection dominated optimal control problems with distributed controls is also supported by our numerical results

Discontinuous Galerkin method for multifluid Euler equations

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/106454/1/AIAA2013-2595.pd