BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Andrew Christlieb

Andrew Christlieb appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NAData Structures and AlgorithmsApplicationscond-mat.softDistributed, Parallel, and Cluster Computingmath.DSNumerical Analysis

Catalog footprint

What is connected

11works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

11 published item(s)

preprint2022arXiv

Competition and Complexity in Amphiphilic Polymer Morphology

We analyze the competitive evolution of codimension one and two morphologies within the $H^{-1}$ gradient flow of the strong Functionalized Cahn-Hilliard equation. On a slow time scale a sharp hypersurface reduction yields a degenerate Mullins-Sekerka evolution for both codimension one and two hypersurfaces, leading to a geometric flow that depends locally on curvatures couples to the dynamic value of the spatially constant far-field chemical potential. Both codimension one and two morphologies admit two classes of bifurcations, one leads to pearling, a short-wavelength in-plane modulation of interfacial width, the other flips motion by curvature to the locally-ill posed motion against curvature, which leads to fingering instabilities. We present a bifurcation diagram for the morphological competition, and compare our results quantitatively to simulations of the full system and qualitatively to simulations of self-consistent mean field models and laboratory experiments; illuminating the role of the pearling bifurcation in the development of complex network morphologies.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Kernel-Based Explicit Unconditionally Stable Scheme for Hamilton-Jacobi Equations on Nonuniform Meshes

In \cite{christlieb2019kernel}, the authors developed a class of high-order numerical schemes for the Hamilton-Jacobi (H-J) equations, which are unconditionally stable, yet take the form of an explicit scheme. This paper extends such schemes, so that they are more effective at capturing sharp gradients, especially on nonuniform meshes. In particular, we modify the weighted essentially non-oscillatory (WENO) methodology in the previously developed schemes by incorporating an exponential basis and adapting the previously developed nonlinear filters used to control oscillations. The main advantages of the proposed schemes are their effectiveness and simplicity, since they can be easily implemented on higher-dimensional nonuniform meshes. We perform numerical experiments on a collection of examples, including H-J equations with linear, nonlinear, convex and non-convex Hamiltonians. To demonstrate the flexibility of the proposed schemes, we also include test problems defined on non-trivial geometry.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A WENO-based Method of Line Transpose Approach for Vlasov Simulations

In this paper, a high order implicit Method of Line Transpose (MOL$^T$ ) method based on a weighted essentially non-oscillatory (WENO) methodology is developed for one-dimensional linear transport equations and further applied to the Vlasov-Poisson (VP) simulations via dimensional splitting. In the MOL$^T$ framework, the time variable is first discretized by a diagonally implicit strong-stability-preserving Runge-Kutta method, resulting in a boundary value problem (BVP) at the discrete time levels. Then an integral formulation coupled with a high order WENO methodology is employed to solve the BVP. As a result, the proposed scheme is high order accurate in both space and time and free of oscillations even though the solution is discontinuous or has sharp gradients. Moreover, the scheme is able to take larger time step evolution compared with an explicit MOL WENO scheme with the same order of accuracy. The desired positivity-preserving (PP) property of the scheme is further attained by incorporating a newly proposed high order PP limiter. We perform numerical experiments on several benchmarks including linear advection, solid body rotation problem; and on the Landau damping, two-stream instabilities, bump-on-tail, and plasma sheath by solving the VP system. The efficacy and efficiency of the proposed scheme is numerically verified.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Method of lines transpose: Energy gradient flows using direct operator inversion for phase field models

In this work, we develop an $\mathcal{O}(N)$ implicit real space method in 1D and 2D for the Cahn Hilliard (CH) and vector Cahn Hilliard (VCH) equations, based on the Method Of Lines Transpose (MOL$^\text{T}$) formulation. This formulation results in a semi-discrete time stepping algorithm, which we prove is gradient stable in the $H^{-1}$ norm. The spatial discretization follows from dimensional splitting, and an $\mathcal{O}(N)$ matrix-free solver, which applies fast convolution to the modified Helmholtz equation. We propose a novel factorization technique, in which fourth order spatial derivatives are incorporated into the solver. The splitting error is included in the nonlinear fixed point iteration, resulting in a high order, logically Cartesian (line-by-line) update. Our method is fast, but not restricted to periodic boundaries like the fast Fourier transform (FFT). The basic solver is implemented using the Backward Euler formulation, and we extend this to both backward difference (BDF) stencils, implicit Runge Kutta (SDIRK) and spectral deferred correction (SDC) frameworks to achieve high orders of temporal accuracy. We demonstrate with numerical results that the CH, and VCH equations maintain gradient stability in one and two spatial dimensions. We also explore time-adaptivity, so that meta-stable states and ripening events can be simulated both quickly and efficiently.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Method of lines transpose: an efficient A-stable solver for wave propagation

Building upon recent results obtained in [7,8,9], we describe an efficient second order, A-stable scheme for solving the wave equation, based on the method of lines transpose (MOL$^T$), and the resulting semi-discrete (i.e. continuous in space) boundary value problem. In [7], A-stable schemes of high order were derived, and in [9] a high order, fast $\mathcal{O}(N)$ spatial solver was derived, which is matrix-free and is based on dimensional-splitting. In this work, are interested in building a wave solver, and our main concern is the development of boundary conditions. We demonstrate all desired boundary conditions for a wave solver, including outflow boundary conditions, in 1D and 2D. The scheme works in a logically Cartesian fashion, and the boundary points are embedded into the regular mesh, without incurring stability restrictions, so that boundary conditions are imposed without any reduction in the order of accuracy. We demonstrate how the embedded boundary approach works in the cases of Dirichlet and Neumann boundary conditions. Further, we develop outflow and periodic boundary conditions for the MOL$^T$ formulation. Our solver is designed to couple with particle codes, and so special attention is also paid to the implementation of point sources, and soft sources which can be used to launch waves into waveguides.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A High Order Time Splitting Method Based on Integral Deferred Correction for Semi-Lagrangian Vlasov Simulations

Semi-Lagrangian schemes with various splitting methods, and with different reconstruction/interpolation strategies have been applied to kinetic simulations. For example, the order of spatial accuracy of the algorithms proposed in {[Qiu and Christlieb, J. Comp. Phys., 2010]} is very high (as high as ninth order). However, the temporal error is dominated by the operator splitting error, which is second order for Strang splitting. It is therefore important to overcome such low order splitting error, in order to have numerical algorithms that achieve higher orders of accuracy in both space and time. In this paper, we propose to use the integral deferred correction (IDC) method to reduce the splitting error. Specifically, the temporal order accuracy is increased by $r$ with each correction loop in the IDC framework, where $r=1,\,2$ for coupling the first order splitting and the Strang splitting, respectively. The proposed algorithm is applied to the Vlasov-Poisson system, the guiding center model, and two dimensional incompressible flow simulations in the vorticity stream-function formulation. We show numerically that the IDC procedure can automatically increase the order of accuracy in time. We also investigate numerical stability of the proposed algorithm via performing Fourier analysis to a linear model problem.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

High Order Maximum Principle Preserving Semi-Lagrangian Finite Difference WENO schemes for the Vlasov Equation

In this paper, we propose the parametrized maximum principle preserving (MPP) flux limiter, originally developed in [Z. Xu, Math. Comp., (2013), in press], to the semi- Lagrangian finite difference weighted essentially non-oscillatory scheme for solving the Vlasov equation. The MPP flux limiter is proved to maintain up to fourth order accuracy for the semi-Lagrangian finite difference scheme without any time step restriction. Numerical studies on the Vlasov-Poisson system demonstrate the performance of the proposed method and its ability in preserving the positivity of the probability distribution function while maintaining the high order accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Multiscale Sub-linear Time Fourier Algorithm for Noisy Data

We extend the recent sparse Fourier transform algorithm of (Lawlor, Christlieb, and Wang, 2013) to the noisy setting, in which a signal of bandwidth N is given as a superposition of k << N frequencies and additive noise. We present two such extensions, the second of which exhibits a novel form of error-correction in its frequency estimation not unlike that of the beta-encoders in analog-to-digital conversion (Daubechies et al, 2006). The algorithm runs in time O(k log(k) log(N/k)) on average, provided the noise is not overwhelming. The error-correction property allows the algorithm to outperform FFTW, a highly optimized software package for computing the full discrete Fourier transform, over a wide range of sparsity and noise values, and is to the best of our knowledge novel in the sparse Fourier transform context.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Adaptive sub-linear Fourier algorithms

We present a new deterministic algorithm for the sparse Fourier transform problem, in which we seek to identify k << N significant Fourier coefficients from a signal of bandwidth N. Previous deterministic algorithms exhibit quadratic runtime scaling, while our algorithm scales linearly with k in the average case. Underlying our algorithm are a few simple observations relating the Fourier coefficients of time-shifted samples to unshifted samples of the input function. This allows us to detect when aliasing between two or more frequencies has occurred, as well as to determine the value of unaliased frequencies. We show that empirically our algorithm is orders of magnitude faster than competing algorithms.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Parallel Semi-Implicit Time Integrators

In this paper, we further develop a family of parallel time integrators known as Revisionist Integral Deferred Correction methods (RIDC) to allow for the semi-implicit solution of time dependent PDEs. Additionally, we show that our semi-implicit RIDC algorithm can harness the computational potential of multiple general purpose graphical processing units (GPUs) in a single node by utilizing existing CUBLAS libraries for matrix linear algebra routines in our implementation. In the numerical experiments, we show that our implementation computes a fourth order solution using four GPUs and four CPUs in approximately the same wall clock time as a first order solution computed using a single GPU and a single CPU.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Scandalously Parallelizable Mesh Generation

We propose a novel approach which employs random sampling to generate an accurate non-uniform mesh for numerically solving Partial Differential Equation Boundary Value Problems (PDE-BVP's). From a uniform probability distribution U over a 1D domain, we sample M discretizations of size N where M>>N. The statistical moments of the solutions to a given BVP on each of the M ultra-sparse meshes provide insight into identifying highly accurate non-uniform meshes. Essentially, we use the pointwise mean and variance of the coarse-grid solutions to construct a mapping Q(x) from uniformly to non-uniformly spaced mesh-points. The error convergence properties of the approximate solution to the PDE-BVP on the non-uniform mesh are superior to a uniform mesh for a certain class of BVP's. In particular, the method works well for BVP's with locally non-smooth solutions. We present a framework for studying the sampled sparse-mesh solutions and provide numerical evidence for the utility of this approach as applied to a set of example BVP's. We conclude with a discussion of how the near-perfect paralellizability of our approach suggests that these strategies have the potential for highly efficient utilization of massively parallel multi-core technologies such as General Purpose Graphics Processing Units (GPGPU's). We believe that the proposed algorithm is beyond embarrassingly parallel; implementing it on anything but a massively multi-core architecture would be scandalous.

0 citations0 reviewsNo code linkedOpen access