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

Daniel Appelö

Daniel Appelö contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysisphysics.acc-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
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

9 published item(s)

preprint2022arXiv

An Energy-Based Discontinuous Galerkin Method with Tame CFL Numbers for the Wave Equation

We extend and analyze the energy-based discontinuous Galerkin method for second order wave equations on staggered and structured meshes. By combining spatial staggering with local time-stepping near boundaries, the method overcomes the typical numerical stiffness associated with high order piecewise polynomial approximations. In one space dimension with periodic boundary conditions and suitably chosen numerical fluxes, we prove bounds on the spatial operators that establish stability for CFL numbers $c \frac {Δt}{h} < C$ independent of order when stability-enhanced explicit time-stepping schemes of matching order are used. For problems on bounded domains and in higher dimensions we demonstrate numerically that one can march explicitly with large time steps at high order temporal and spatial accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

An energy-based summation-by-parts finite difference method for the wave equation in second order form

We develop a new finite difference method for the wave equation in second order form. The finite difference operators satisfy a summation-by-parts (SBP) property. With boundary conditions and material interface conditions imposed weakly by the simultaneous-approximation-term (SAT) method, we derive energy estimates for the semi-discretization. In addition, error estimates are derived by the normal mode analysis. The proposed method is termed as energy-based because of its similarity with the energy-based discontinuous Galerkin method. When imposing the Dirichlet boundary condition and material interface conditions, the traditional SBP-SAT discretization uses a penalty term with a mesh-dependent parameter, which is not needed in our method. Furthermore, numerical dissipation can be added to the discretization through the boundary and interface conditions. We present numerical experiments that verify convergence and robustness of the proposed method.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

EM-WaveHoltz: A flexible frequency-domain method built from time-domain solvers

A novel approach to computing time-harmonic solutions of Maxwell&#39;s equations by time-domain simulations is presented. The method, EM-WaveHoltz, results in a positive definite system of equations which makes it amenable to iterative solution with the conjugate gradient method or with GMRES. Theoretical results guaranteeing the convergence of the method away from resonances is presented. Numerical examples illustrating the properties of EM-WaveHoltz are given.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Extensions and Analysis of an Iterative Solution of the Helmholtz Equation via the Wave Equation

In this paper we extend analysis of the WaveHoltz iteration -- a time-domain iterative method for the solution of the Helmholtz equation. We expand the previous analysis of energy conserving problems and prove convergence of the WaveHoltz iteration for problems with impedance boundary conditions in a single spatial dimension. We then consider interior Dirichlet/Neumann problems with damping in any spatial dimension, and show that for a sufficient level of damping the WaveHoltz iteration converges in a number of iteration independent of the frequency. Finally, we present a discrete analysis of the WaveHoltz iteration for a family of higher order time-stepping schemes. We show that the fixed-point of the discrete WaveHoltz iteration converges to the discrete Helmholtz solution with the order of the time-stepper chosen. We present numerical examples and demonstrate that it is possible to completely remove time discretization error from the WaveHoltz solution through careful analysis of the discrete iteration together with updated quadrature formulas.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Universal AMG Accelerated Embedded Boundary Method Without Small Cell Stiffness

We develop a universally applicable embedded boundary finite difference method, which results in a symmetric positive definite linear system and does not suffer from small cell stiffness. Our discretization is efficient for the wave, heat and Poisson&#39;s equation with Dirichlet boundary conditions. When the system needs to be inverted we can use the conjugate gradient method, accelerated by algebraic multigrid techniques. A series of numerical tests for the wave, heat and Poisson&#39;s equation and applications to shape optimization problems verify the accuracy, stability, and efficiency of our method. Our fast computational techniques can be extended to moving boundary problems (e.g. Stefan problem), to the Navier-Stokes equations, and to the Grad-Shafranov equations for which problems are posed on domains with complex geometry and fast simulations are very important.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A Hermite Method with a Discontinuity Sensor for Hamilton-Jacobi Equations

We present a Hermite interpolation based partial differential equation solver for Hamilton-Jacobi equations. Many Hamilton-Jacobi equations have a nonlinear dependency on the gradient, which gives rise to discontinuities in the derivatives of the solution, resulting in kinks. We built our solver with two goals in mind: 1) high order accuracy in smooth regions and 2) sharp resolution of kinks. To achieve this, we use Hermite interpolation with a smoothness sensor. The degrees-of freedom of Hermite methods are tensor-product Taylor polynomials of degree $m$ in each coordinate direction. The method uses $(m+1)^d$ degrees of freedom per node in $d$-dimensions and achieves an order of accuracy $(2m+1)$ when the solution is smooth. To obtain sharp resolution of kinks, we sense the smoothness of the solution on each cell at each timestep. If the solution is smooth, we march the interpolant forward in time with no modifications. When our method encounters a cell over which the solution is not smooth, it introduces artificial viscosity locally while proceeding normally in smooth regions. We show through numerical experiments that the solver sharply captures kinks once the solution losses continuity in the derivative while achieving $2m+1$ order accuracy in smooth regions.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Re-evaluation of Spin-Orbit Dynamics of Polarized e+ e- Beams in High Energy Circular Accelerators and Storage Rings: an approach based on a Bloch equation

We give an overview of our current/future analytical and numerical work on the spin polarization in high-energy electron storage rings. Our goal is to study the possibility of polarization for the CEPC and FCC-ee. Our work is based on the so-called Bloch equation for the polarization density introduced by Derbenev and Kondratenko in 1975. We also give an outline of the standard approach, the latter being based on the Derbenev-Kondratenko formulas.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Hermite-Discontinuous Galerkin Overset Grid Methods for the Scalar Wave Equation

We present high order accurate numerical methods for the wave equation that combines efficient Hermite methods with eometrically flexible discontinuous Galerkin methods by using overset grids. Near boundaries we use thin boundary fitted curvilinear grids and in the volume we use Cartesian grids so that the computational complexity of the solvers approach a structured Cartesian Hermite method. Unlike many other overset methods we do not need to add artificial dissipation but we find that the built in dissipation of the Hermite and discontinuous Galerkin methods is sufficient to maintain stability. By numerical experiments we demonstrate the stability, accuracy, efficiency and applicability of the methods to forward and inverse problems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Bloch equation for spin dynamics in electron storage rings: computational and theoretical aspects

In this paper we describe our work on spin polarization in high-energy electron storage rings which we base on the Bloch equation for the polarization density and which aims towards the e-/e+ option of the proposed Future Circular Collider (FCC-ee) and the proposed Circular Electron Positron Collider (CEPC). The Bloch equation takes into account non spin-flip and spin-flip effects due to synchrotron radiation including the spin-diffusion effects and the Sokolov-Ternov effect with its Baier-Katkov generalization as well as the kinetic-polarization effect. This mathematical model is an alternative to the standard mathematical model based on the Derbenev-Kondratenko formulas. For our numerical and analytical studies of the Bloch equation we develop an approximation to the latter to obtain an effective Bloch equation. This is accomplished by finding a third mathematical model based on a system of stochastic differential equations underlying the Bloch equation and by approximating that system via the method of averaging from perturbative ODE theory. We also give an overview of our algorithm for numerically integrating the effective Bloch equation. This discretizes the phase space using spectral methods and discretizes time via the additive Runge-Kutta method which is a high-order semi-implicit method. We also discuss the relevance of the third mathematical model for spin tracking.

0 citations0 reviewsNo code linkedOpen access