Researcher profile

Thomas Lewis

Thomas Lewis contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
6works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

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

6 published item(s)

preprint2023arXiv

Isolating the extreme debris disc signature -- explorations of eccentric extreme debris discs formed by giant impacts

In this work, we used N-body simulations and a radiative transfer package to model the evolution of eccentric debris discs produced by giant impacts between planetary embryos. This included how the morphology and infrared emission of these discs varied with embryo eccentricity and collision true anomaly. We found that eccentric discs inherit the eccentric properties of the centre of mass orbit of the two colliding embryos. However, the orientation of the collision with the respect to this orbit plays a key role in determining how closely the disc material resembles the centre of mass orbit. Additionally, we found that increased eccentricity acted to suppress the formation of certain short-term variations in the disc emission depending on the collision position. These short-term variations have been associated with an observational phenomenon called extreme debris discs. Short-term variability has been suggested as a potential signature for giant impacts.

preprint2022arXiv

A narrow-stencil framework for convergent numerical approximations of fully nonlinear second order PDEs

This paper develops a unified general framework for designing convergent finite difference and discontinuous Galerkin methods for approximating viscosity and regular solutions of fully nonlinear second order PDEs. Unlike the well-known monotone (finite difference) framework, the proposed new framework allows for the use of narrow stencils and unstructured grids which makes it possible to construct high order methods. The general framework is based on the concepts of consistency and g-monotonicity which are both defined in terms of various numerical derivative operators. Specific methods that satisfy the framework are constructed using numerical moments. Admissibility, stability, and convergence properties are proved, and numerical experiments are provided along with some computer implementation details.

preprint2013arXiv

Convergent finite difference methods for one-dimensional fully nonlinear second order partial differential equations

This paper develops a new framework for designing and analyzing convergent finite difference methods for approximating both classical and viscosity solutions of second order fully nonlinear partial differential equations (PDEs) in 1-D. The goal of the paper is to extend the successful framework of monotone, consistent, and stable finite difference methods for first order fully nonlinear Hamilton-Jacobi equations to second order fully nonlinear PDEs such as Monge-Ampère and Bellman type equations. New concepts of consistency, generalized monotonicity, and stability are introduced; among them, the generalized monotonicity and consistency, which are easier to verify in practice, are natural extensions of the corresponding notions of finite difference methods for first order fully nonlinear Hamilton-Jacobi equations. The main component of the proposed framework is the concept of "numerical operator", and the main idea used to design consistent, monotone and stable finite difference methods is the concept of "numerical moment". These two new concepts play the same roles as the "numerical Hamiltonian" and the "numerical viscosity" play in the finite difference framework for first order fully nonlinear Hamilton-Jacobi equations. In the paper, two classes of consistent and monotone finite difference methods are proposed for second order fully nonlinear PDEs. The first class contains Lax-Friedrichs-like methods which also are proved to be stable and the second class contains Godunov-like methods. Numerical results are also presented to gauge the performance of the proposed finite difference methods and to validate the theoretical results of the paper.

preprint2013arXiv

Discontinuous Galerkin finite element differential calculus and applications to numerical solutions of linear and nonlinear partial differential equations

This paper develops a discontinuous Galerkin (DG) finite element differential calculus theory for approximating weak derivatives of Sobolev functions and piecewise Sobolev functions. By introducing numerical one-sided derivatives as building blocks, various first and second order numericaloperators such as the gradient, divergence, Hessian, and Laplacian operator are defined, and their corresponding calculus rules are established. Among the calculus rules are product and chain rules, integration by parts formulas and the divergence theorem. Approximation properties and the relationship between the proposed DG finite element numerical derivatives and some well-known finite difference numerical derivative formulas on Cartesian grids are also established. Efficient implementation of the DG finite element numerical differential operators is also proposed. Besides independent interest in numerical differentiation, the primary motivation and goal of developing the DG finite element differential calculus is to solve partial differential equations. It is shown that several existing finite element, finite difference and DG methods can be rewritten compactly using the proposed DG finite element differential calculus framework. Moreover, new DG methods for linear and nonlinear PDEs are also obtained from the framework.

preprint2012arXiv

Local discontinuous Galerkin methods for one-dimensional second order fully nonlinear elliptic and parabolic equations

This paper is concerned with developing accurate and efficient discontinuous Galerkin methods for fully nonlinear second order elliptic and parabolic partial differential equations (PDEs) in the case of one spatial dimension. The primary goal of the paper to develop a general framework for constructing high order local discontinuous Galerkin (LDG) methods for approximating viscosity solutions of these fully nonlinear PDEs which are merely continuous functions by definition. In order to capture discontinuities of the first order derivative $u_x$ of the solution $u$, two independent functions $q_1$ and $q_2$ are introduced to approximate one-sided derivatives of $u$. Similarly, to capture the discontinuities of the second order derivative $u_{xx}$, four independent functions $p_{1}$, $p_{2}$, $p_{3}$, and $p_{4}$ are used to approximate one-sided derivatives of $q_1$ and $q_2$. The proposed LDG framework, which is based on a nonstandard mixed formulation of the underlying PDE, embeds a given fully nonlinear problem into a mostly linear system of equations where the given nonlinear differential operator must be replaced by a numerical operator which allows multiple value inputs of the first and second order derivatives $u_x$ and $u_{xx}$. An easy to verify criterion for constructing "good" numerical operators is also proposed. It consists of a consistency and a generalized monotonicity. To ensure such a generalized monotonicity, the crux of the construction is to introduce the numerical moment in the numerical operator. The proposed framework extends a companion finite difference framework developed by the authors in [9] and allows for the approximation of fully nonlinear PDEs using high order polynomials and non-uniform meshes.

preprint2012arXiv

Mixed Interior Penalty Discontinuous Galerkin Methods for One-dimensional Fully Nonlinear Second Order Elliptic and Parabolic Equations

This paper is concerned with developing accurate and efficient numerical methods for one-dimensional fully nonlinear second order elliptic and parabolic partial differential equations (PDEs). In the paper we present a general framework for constructing high order interior penalty discontinuous Galerkin (IP-DG) methods for approximating viscosity solutions of these fully nonlinear PDEs. In order to capture discontinuities of the second order derivative $u_{xx}$ of the solution $u$, three independent functions $p_1, p_2$ and $p_3$ are introduced to represent numerical derivatives using various one-sided limits. The proposed DG framework, which is based on a nonstandard mixed formulation of the underlying PDE, embeds a nonlinear problem into a mostly linear system of equations where the nonlinearity has been modified to include multiple values of the second order derivative $u_{xx}$. The proposed framework extends a companion finite difference framework developed by the authors in [9] and allows for the approximation of fully nonlinear PDEs using high order polynomials and non-uniform meshes.In addition to the nonstandard mixed formulation setting, another main idea is to replace the fully nonlinear differential operator by a numerical operator, which is consistent with the differential operator and satisfies certain monotonicity (called g-monotonicity) properties. To ensure such a g-monotonicity, the crux of the construction is to introduce the numerical moment, which plays a critical role in the proposed DG framework. This paper also presents and analyzes numerical results for several numerical test problems which are used to guage the accuracy and efficiency of the proposed DG methods.