BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Garth N. Wells

Garth N. Wells contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
Computational Engineering, Finance, and ScienceNumerical Analysismath.NAMathematical Softwarephysics.geo-phphysics.flu-dyncond-mat.mtrl-scimath.OCphysics.comp-phSymbolic ComputationSystems and Control

Trust snapshot

Quick read

Trust 21 - Emerging
14works
0followers
11topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2016arXiv

Automatic calibration of damping layers in finite element time domain simulations

Matched layers are commonly used in numerical simulations of wave propagation to model (semi-)infinite domains. Attenuation functions describe the damping in layers, and provide a matching of the wave impedance at the interface between the domain of interest and the absorbing region. Selecting parameters in the attenuation functions is non-trivial. In this work, an optimisation procedure for automatically calibrating matched layers is presented. The procedure is based on solving optimisation problems constrained by partial differential equations with polynomial and piecewise-constant attenuation functions. We show experimentally that, for finite element time domain simulations, piecewise-constant attenuation function are at least as efficient as quadratic attenuation functions. This observation leads us to introduce consecutive matched layers as an alternative to perfectly matched layers, which can easily be employed for problems with arbitrary geometries. Moreover, the use of consecutive matched layers leads to a reduction in computational cost compared to perfectly matched layers. Examples are presented for acoustic, elastodynamic and electromagnetic problems. Numerical simulations are performed with the libraries FEniCS/DOLFIN and dolfin-adjoint, and the computer code to reproduce all numerical examples is made freely available.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

CO$_2$ dissolution in a background hydrological flow

We investigate the long time steady-state dissolution of CO$_{2}$ in a deep saline aqquifer in the presence of a background hydrological flow. In steady-state, the distribution of CO$_2$ in the groundwater upstream of the aquifer involves a balance between three competing effects: (i) the buoyancy-driven flow of CO$_2$ saturated water; (ii) the diffusion of CO$_2$ from saturated to under-saturated water; and (iii) the advection associated with the oncoming background flow. This leads to three limiting regimes. In the limit of very slow diffusion, a nearly static intrusion of dense fluid may extend a finite distance upstream, balanced by the pressure gradient associated with the oncoming background flow. In the limit of fast diffusion relative to the flow, a gradient zone may become established in which the along aquifer diffusive flux balances the advection associated with the background flow. However, if the buoyancy-driven flow speed exceeds the background hydrological flow speed, then a third, intermediate regime may become established. In this regime, a convective recirculation develops upstream of the anticline involving the vertical diffusion of CO$_2$ from an upstream propagating flow of dense CO$_2$ saturated water into the downstream propagating flow of CO$_2$ unsaturated water. For each limiting case, we find analytical solutions for the distribution of CO$_2$ upstream of the anticline, and test our analysis with full numerical simulations. A key result is that, although there may be very different controls on the distribution and extent of CO$_2$ bearing water upstream of the anticline, in each case the dissolution rate is given by the product of the background volume flux and the difference in concentration between the CO$_2$ saturated water and the original aquifer water upstream.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Three-field block-preconditioners for models of coupled magma/mantle dynamics

For a prescribed porosity, the coupled magma/mantle flow equations can be formulated as a two-field system of equations with velocity and pressure as unknowns. Previous work has shown that while optimal preconditioners for the two-field formulation can be obtained, the construction of preconditioners that are uniform with respect to model parameters is difficult. This limits the applicability of two-field preconditioners in certain regimes of practical interest. We address this issue by reformulating the governing equations as a three-field problem, which removes a term that was problematic in the two-field formulation in favour of an additional equation for a pressure-like field. For the three-field problem, we develop and analyse new preconditioners and we show numerically that they are optimal in terms of problem size and less sensitive to model parameters, compared to the two-field preconditioner. This extends the applicability of optimal preconditioners for coupled mantle/magma dynamics into parameter regimes of physical interest.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Torsion of a cylinder of partially molten rock with a spherical inclusion: theory and simulation

The processes that are involved in migration and extraction of melt from the mantle are not yet fully understood. Gaining a better understanding of material properties of partially molten rock could help shed light on the behavior of melt on larger scales in the mantle. In this study, we simulate three-dimensional torsional deformation of a partially molten rock that contains a rigid, spherical inclusion. We compare the computed porosity patterns to those found in recent laboratory experiments. The laboratory experiments show emergence of melt-rich bands throughout the rock sample, and pressure shadows around the inclusion. The numerical model displays similar melt-rich bands only for a small bulk-to-shear-viscosity ratio (five or less). The results are consistent with earlier two-dimensional numerical simulations; however, we show that it is easier to form melt-rich bands in three dimensions compared to two. The addition of strain-rate dependence of the viscosity causes a distinct change in the shape of pressure shadows around the inclusion. This change in shape presents an opportunity for experimentalists to identify the strain-rate dependence and therefore the dominant deformation mechanism in torsion experiments with inclusions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Analysis of block-preconditioners for models of coupled magma/mantle dynamics

This article considers the iterative solution of a finite element discretisation of the magma dynamics equations. In simplified form, the magma dynamics equations share some features of the Stokes equations. We therefore formulate, analyse and numerically test a Elman, Silvester and Wathen-type block preconditioner for magma dynamics. We prove analytically and demonstrate numerically the optimality of the preconditioner. The presented analysis highlights the dependence of the preconditioner on parameters in the magma dynamics equations that can affect convergence of iterative linear solvers. The analysis is verified through a range of two- and three-dimensional numerical examples on unstructured grids, from simple illustrative problems through to large problems on subduction zone-like geometries. The computer code to reproduce all numerical examples is freely available as supporting material.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Unified Form Language: A domain-specific language for weak formulations of partial differential equations

We present the Unified Form Language (UFL), which is a domain-specific language for representing weak formulations of partial differential equations with a view to numerical approximation. Features of UFL include support for variational forms and functionals, automatic differentiation of forms and expressions, arbitrary function space hierarchies for multi-field problems, general differential operators and flexible tensor algebra. With these features, UFL has been used to effortlessly express finite element methods for complex systems of partial differential equations in near-mathematical notation, resulting in compact, intuitive and readable programs. We present in this work the language and its construction. An implementation of UFL is freely available as an open-source software library. The library generates abstract syntax tree representations of variational problems, which are used by other software libraries to generate concrete low-level implementations. Some application examples are presented and libraries that support UFL are highlighted.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Energy stable and momentum conserving hybrid finite element method for the incompressible Navier-Stokes equations

A hybrid method for the incompressible Navier--Stokes equations is presented. The method inherits the attractive stabilizing mechanism of upwinded discontinuous Galerkin methods when momentum advection becomes significant, equal-order interpolations can be used for the velocity and pressure fields, and mass can be conserved locally. Using continuous Lagrange multiplier spaces to enforce flux continuity across cell facets, the number of global degrees of freedom is the same as for a continuous Galerkin method on the same mesh. Different from our earlier investigations on the approach for the Navier--Stokes equations, the pressure field in this work is discontinuous across cell boundaries. It is shown that this leads to very good local mass conservation and, for an appropriate choice of finite element spaces, momentum conservation. Also, a new form of the momentum transport terms for the method is constructed such that global energy stability is guaranteed, even in the absence of a point-wise solenoidal velocity field. Mass conservation, momentum conservation and global energy stability are proved for the time-continuous case, and for a fully discrete scheme. The presented analysis results are supported by a range of numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A FEniCS-Based Programming Framework for Modeling Turbulent Flow by the Reynolds-Averaged Navier-Stokes Equations

Finding an appropriate turbulence model for a given flow case usually calls for extensive experimentation with both models and numerical solution methods. This work presents the design and implementation of a flexible, programmable software framework for assisting with numerical experiments in computational turbulence. The framework targets Reynolds-averaged Navier-Stokes models, discretized by finite element methods. The novel implementation makes use of Python and the FEniCS package, the combination of which leads to compact and reusable code, where model- and solver-specific code resemble closely the mathematical formulation of equations and algorithms. The presented ideas and programming techniques are also applicable to other fields that involve systems of nonlinear partial differential equations. We demonstrate the framework in two applications and investigate the impact of various linearizations on the convergence properties of nonlinear solvers for a Reynolds-averaged Navier-Stokes model.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Analysis of an interface stabilised finite element method: The advection-diffusion-reaction equation

Analysis of an interface stabilised finite element method for the scalar advection-diffusion-reaction equation is presented. The method inherits attractive properties of both continuous and discontinuous Galerkin methods, namely the same number of global degrees of freedom as a continuous Galerkin method on a given mesh and the stability properties of discontinuous Galerkin methods for advection dominated problems. Simulations using the approach in other works demonstrated good stability properties with minimal numerical dissipation, and standard convergence rates for the lowest order elements were observed. In this work, stability of the formulation, in the form of an inf-sup condition for the hyperbolic limit and coercivity for the elliptic case, is proved, as is order $k+1/2$ order convergence for the advection-dominated case and order $k+1$ convergence for the diffusive limit in the $L^{2}$ norm. The analysis results are supported by a number of numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Automated code generation for discontinuous Galerkin methods

A compiler approach for generating low-level computer code from high-level input for discontinuous Galerkin finite element forms is presented. The input language mirrors conventional mathematical notation, and the compiler generates efficient code in a standard programming language. This facilitates the rapid generation of efficient code for general equations in varying spatial dimensions. Key concepts underlying the compiler approach and the automated generation of computer code are elaborated. The approach is demonstrated for a range of common problems, including the Poisson, biharmonic, advection--diffusion and Stokes equations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

DOLFIN: Automated Finite Element Computing

We describe here a library aimed at automating the solution of partial differential equations using the finite element method. By employing novel techniques for automated code generation, the library combines a high level of expressiveness with efficient computation. Finite element variational forms may be expressed in near mathematical notation, from which low-level code is automatically generated, compiled and seamlessly integrated with efficient implementations of computational meshes and high-performance linear algebra. Easy-to-use object-oriented interfaces to the library are provided in the form of a C++ library and a Python module. This paper discusses the mathematical abstractions and methods used in the design of the library and its implementation. A number of examples are presented to demonstrate the use of the library in application code.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Optimal control with stochastic PDE constraints and uncertain controls

The optimal control of problems that are constrained by partial differential equations with uncertainties and with uncertain controls is addressed. The Lagrangian that defines the problem is postulated in terms of stochastic functions, with the control function possibly decomposed into an unknown deterministic component and a known zero-mean stochastic component. The extra freedom provided by the stochastic dimension in defining cost functionals is explored, demonstrating the scope for controlling statistical aspects of the system response. One-shot stochastic finite element methods are used to find approximate solutions to control problems. It is shown that applying the stochastic collocation finite element to the formulated problem leads to a coupling between stochastic collocation points when a deterministic optimal control is considered or when moments are included in the cost functional, thereby obviating the primary advantage of the collocation method over the stochastic Galerkin method for the considered problem. The application of the presented methods is demonstrated through a number of numerical examples. The presented framework is sufficiently general to also consider a class of inverse problems, and numerical examples of this type are also presented.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Optimisations for quadrature representations of finite element tensors through automated code generation

We examine aspects of the computation of finite element matrices and vectors which are made possible by automated code generation. Given a variational form in a syntax which resembles standard mathematical notation, the low-level computer code for building finite element tensors, typically matrices, vectors and scalars, can be generated automatically via a form compiler. In particular, the generation of code for computing finite element matrices using a quadrature approach is addressed. For quadrature representations, a number of optimisation strategies which are made possible by automated code generation are presented. The relative performance of two different automatically generated representations of finite element matrices is examined, with a particular emphasis on complicated variational forms. It is shown that approaches which perform best for simple forms are not tractable for more complicated problems in terms of run time performance, the time required to generate the code or the size of the generated code. The approach and optimisations elaborated here are effective for a range of variational forms.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Phase field model for coupled displacive and diffusive microstructural processes under thermal loading

A non-isothermal phase field model that captures both displacive and diffusive phase transformations in a unified framework is presented. The model is developed in a formal thermodynamic setting, which provides guidance on admissible constitutive relationships and on the coupling of the numerous physical processes that are active. Phase changes are driven by temperature-dependent free-energy functions that become non-convex below a transition temperature. Higher-order spatial gradients are present in the model to account for phase boundary energy, and these terms necessitate the introduction of non-standard terms in the energy balance equation in order to satisfy the classical entropy inequality point-wise. To solve the resulting balance equations, a Galerkin finite element scheme is elaborated. To deal rigorously with the presence of high-order spatial derivatives associated with surface energies at phase boundaries in both the momentum and mass balance equations, some novel numerical approaches are used. Numerical examples are presented that consider boundary cooling of a domain at different rates, and these results demonstrate that the model can qualitatively reproduce the evolution of microstructural features that are observed in some alloys, especially steels. The proposed model opens a number of interesting possibilities for simulating and controlling microstructure pattern development under combinations of thermal and mechanical loading.

0 citations0 reviewsNo code linkedOpen access