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

Yaman Güçlü

Yaman Güçlü appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NAphysics.comp-phmath-phmath.MPNumerical Analysisphysics.class-phphysics.plasm-ph

Catalog footprint

What is connected

5works
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

5 published item(s)

preprint2022arXiv

A mimetic discretization of the macroscopic Maxwell equations in Hamiltonian form

A mimetic spectral element discretization, utilizing a novel Galerkin projection Hodge star operator, of the macroscopic Maxwell equations in Hamiltonian form is presented. The idea of splitting purely topological and metric dependent quantities is natural in the Hamiltonian modeling framework as the Poisson bracket is metric free with the Hamiltonian containing all metric information. This idea may be incorporated into the mimetic spectral element method by directly discretizing the Poincaré duality structure. This "split exterior calculus mimetic spectral element method" yields spatially discretized Maxwell's equations which are Hamiltonian and exactly and strongly conserve Gauss's laws. Moreover, the new discrete Hodge star operator is itself of interest as a partition of the purely topological and metric dependent portions of the Hodge star operator. As a simple test case, the numerical results of applying this method to a one-dimensional version of Maxwell's equations are given.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Arbitrarily high order Convected Scheme solution of the Vlasov-Poisson system

The Convected Scheme (CS) is a `forward-trajectory' semi-Lagrangian method for solution of transport equations, which has been most often applied to the kinetic description of plasmas and rarefied neutral gases. In its simplest form, the CS propagates the solution by advecting the `moving cells' along their characteristic trajectories, and by remapping them on the mesh at the end of the time step. The CS is conservative, positivity preserving, simple to implement, and not subject to time step restriction to maintain stability. Recently [Y. Güçlü and W.N.G. Hitchon, 2012] a new methodology was introduced for reducing numerical diffusion, based on a modified equation analysis: the remapping error was compensated by applying small corrections to the final position of the moving cells prior to remapping. While the spatial accuracy was increased from 2nd to 4th order, the new scheme retained the important properties of the original method, and was shown to be simple and efficient for constant advection problems. Here the CS is applied to the solution of the Vlasov-Poisson system: the Vlasov equation is split into two constant advection equations, one in configuration space and one in velocity space, and high order time accuracy is achieved by proper composition of the operators. The splitting procedure enables us to use the constant advection solver, which we extend to arbitrarily high order of accuracy: a new improved procedure is given, which makes the calculation of the corrections straightforward. Focusing on periodic domains, we describe a spectrally accurate scheme based on the fast Fourier transform; the proposed implementation is strictly conservative and positivity preserving. The ability to correctly reproduce the system dynamics, as well as resolving small-scale features in the solution, is shown in classical 1D-1V test cases, both in the linear and the non-linear regimes.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Picard integral formulation of weighted essentially non-oscillatory schemes

High-order temporal discretizations for hyperbolic conservation laws have historically been formulated as either a method of lines (MOL) or a Lax-Wendroff method. In the MOL viewpoint, the partial differential equation is treated as a large system of ordinary differential equations (ODEs), where an ODE tailored time-integrator is applied. In contrast, Lax-Wendroff discretizations immediately convert Taylor series in time to discrete spatial derivatives. In this work, we propose the Picard integral formulation (PIF), which is based on the method of modified fluxes, and is used to derive new Taylor and Runge-Kutta (RK) methods. In particular, we construct a new class of conservative finite difference methods by applying WENO reconstructions to the so-called "time-averaged" fluxes. Our schemes are automatically conservative under any modification of the fluxes, which is attributed to the fact that classical WENO reconstructions conserve mass when coupled with forward Euler time steps. The proposed Lax-Wendroff discretization is constructed by taking Taylor series of the flux function as opposed to Taylor series of the conserved variables. The RK discretization differs from classical MOL formulations because we apply WENO reconstructions to time-averaged fluxes rather than taking linear combinations of spatial derivatives of the flux. In both cases, we only need one projection onto the characteristic variables per time step. The PIF is generic, and lends itself to a multitude of options for further investigation. At present, we present two canonical examples: one based on Taylor, and the other based on the classical RK method. Stability analyses are presented for each method. The proposed schemes are applied to hyperbolic conservation laws in one- and two-dimensions and the results are in good agreement with current state of the art methods.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A high order cell-centered semi-Lagrangian scheme for multi-dimensional kinetic simulations of neutral gas flows

The term `Convected Scheme' (CS) refers to a family of algorithms, most usually applied to the solution of Boltzmann's equation, which uses a method of characteristics in an integral form to project an initial cell forward to a group of final cells. As such the CS is a `forward-trajectory' semi-Lagrangian scheme. For multi-dimensional simulations of neutral gas flows, the cell-centered version of this semi-Lagrangian (CCSL) scheme has advantages over other options due to its implementation simplicity, low memory requirements, and easier treatment of boundary conditions. The main drawback of the CCSL-CS to date has been its high numerical diffusion in physical space, because of the 2$^{\text{nd}}$ order remapping that takes place at the end of each time step. By means of a Modified Equation Analysis, it is shown that a high order estimate of the remapping error can be obtained a priori, and a small correction to the final position of the cells can be applied upon remapping, in order to achieve full compensation of this error. The resulting scheme is 4$^{\text{th}}$ order accurate in space while retaining the desirable properties of the CS: it is conservative and positivity-preserving, and the overall algorithm complexity is not appreciably increased. Two monotone (i.e. non-oscillating) versions of the fourth order CCSL-CS are also presented: one uses a common flux-limiter approach; the other uses a non-polynomial reconstruction to evaluate the derivatives of the density function. The method is illustrated in simple one- and two-dimensional examples, and a fully 3D solution of the Boltzmann equation describing expansion of a gas into vacuum through a cylindrical tube.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

High-order multiderivative time integrators for hyperbolic conservation laws

Multiderivative time integrators have a long history of development for ordinary differential equations, and yet to date, only a small subset of these methods have been explored as a tool for solving partial differential equations (PDEs). This large class of time integrators include all popular (multistage) Runge-Kutta as well as single-step (multiderivative) Taylor methods. (The latter are commonly referred to as Lax-Wendroff methods when applied to PDEs.) In this work, we offer explicit multistage multiderivative time integrators for hyperbolic conservation laws. Like Lax-Wendroff methods, multiderivative integrators permit the evaluation of higher derivatives of the unknown in order to decrease the memory footprint and communication overhead. Like traditional Runge-Kutta methods, multiderivative integrators admit the addition of extra stages, which introduce extra degrees of freedom that can be used to increase the order of accuracy or modify the region of absolute stability. We describe a general framework for how these methods can be applied to two separate spatial discretizations: the discontinuous Galerkin (DG) method and the finite difference essentially non-oscillatory (FD-WENO) method. The two proposed implementations are substantially different: for DG we leverage techniques that are closely related to generalized Riemann solvers; for FD-WENO we construct higher spatial derivatives with central differences. Among multiderivative time integrators, we argue that multistage two-derivative methods have the greatest potential for multidimensional applications, because they only require the flux function and its Jacobian, which is readily available. Numerical results indicate that multiderivative methods are indeed competitive with popular strong stability preserving time integrators.

0 citations0 reviewsNo code linkedOpen access