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

David Shirokoff

David Shirokoff contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NAphysics.comp-phphysics.flu-dynnlin.CDastro-ph.COhep-phmath.DS

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
6works
0followers
7topics
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)

preprint2015arXiv

A Fourier penalty method for solving the time-dependent Maxwell's equations in domains with curved boundaries

We present a high order, Fourier penalty method for the Maxwell's equations in the vicinity of perfect electric conductor boundary conditions. The approach relies on extending the smooth non-periodic domain of the equations to a periodic domain by removing the exact boundary conditions and introducing an analytic forcing term in the extended domain. The forcing, or penalty term is chosen to systematically enforce the boundary conditions to high order in the penalty parameter, which then allows for higher order numerical methods. We present an efficient numerical method for constructing the penalty term, and discretize the resulting equations using a Fourier spectral method. We demonstrate convergence orders of up to 3.5 for the one-dimensional Maxwell's equations, and show that the numerical method does not suffer from dispersion (or pollution) errors. We also illustrate the approach in two dimensions and demonstrate convergence orders of 2.5 for transverse magnetic modes and 1.5 for the transverse electric modes. We conclude the paper with numerous test cases in dimensions two and three including waves traveling in an irregular waveguide, and scattering off of a windmill-like geometry.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Sharp-Interface Active Penalty Method for the Incompressible Navier-Stokes Equations

The volume penalty method provides a simple, efficient approach for solving the incompressible Navier-Stokes equations in domains with boundaries or in the presence of moving objects. Despite the simplicity, the method is typically limited to first order spatial accuracy. We demonstrate that one may achieve high order accuracy by introducing an active penalty term. One key difference from other works is that we use a sharp, unregularized mask function. We discuss how to construct the active penalty term, and provide numerical examples, in dimensions one and two. We demonstrate second and third order convergence for the heat equation, and second order convergence for the Navier-Stokes equations. In addition, we show that modifying the penalty term does not significantly alter the time step restriction from that of the conventional penalty method.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

An efficient method for the incompressible Navier-Stokes equations on irregular domains with no-slip boundary conditions, high order up to the boundary

Common efficient schemes for the incompressible Navier-Stokes equations, such as projection or fractional step methods, have limited temporal accuracy as a result of matrix splitting errors, or introduce errors near the domain boundaries (which destroy uniform convergence to the solution). In this paper we recast the incompressible (constant density) Navier-Stokes equations (with the velocity prescribed at the boundary) as an equivalent system, for the primary variables velocity and pressure. We do this in the usual way away from the boundaries, by replacing the incompressibility condition on the velocity by a Poisson equation for the pressure. The key difference from the usual approaches occurs at the boundaries, where we use boundary conditions that unequivocally allow the pressure to be recovered from knowledge of the velocity at any fixed time. This avoids the common difficulty of an, apparently, over-determined Poisson problem. Since in this alternative formulation the pressure can be accurately and efficiently recovered from the velocity, the recast equations are ideal for numerical marching methods. The new system can be discretized using a variety of methods, in principle to any desired order of accuracy. In this work we illustrate the approach with a 2-D second order finite difference scheme on a Cartesian grid, and devise an algorithm to solve the equations on domains with curved (non-conforming) boundaries, including a case with a non-trivial topology (a circular obstruction inside the domain). This algorithm achieves second order accuracy (in L-infinity), for both the velocity and the pressure. The scheme has a natural extension to 3-D.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Bouncing droplets on a billiard table

In a set of experiments, Couder et. al. demonstrate that an oscillating fluid bed may propagate a bouncing droplet through the guidance of the surface waves. We present a dynamical systems model, in the form of an iterative map, for a droplet on an oscillating bath. We examine the droplet bifurcation from bouncing to walking, and prescribe general requirements for the surface wave to support stable walking states. We show that in addition to walking, there is a region of large forcing that may support the chaotic bouncing of the droplet. Using the map, we then investigate the droplet trajectories for two different wave responses in a square (billiard ball) domain. We show that for waves which are quickly damped in space, the long time trajectories in a square domain are either non-periodic dense curves, or approach a quasiperiodic orbit. In contrast, for waves which extend over many wavelengths, at low forcing, trajectories tend to approach an array of circular attracting sets. As the forcing increases, the attracting sets break down and the droplet travels throughout space.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Renormalized waves and thermalization of the Klein-Gordon equation: What sound does a nonlinear string make?

We study the thermalization of the classical Klein-Gordon equation under a u^4 interaction. We numerically show that even in the presence of strong nonlinearities, the local thermodynamic equilibrium state exhibits a weakly nonlinear behavior in a renormalized wave basis. The renormalized basis is defined locally in time by a linear transformation and the requirement of vanishing wave-wave correlations. We show that the renormalized waves oscillate around one frequency, and that the frequency dispersion relation undergoes a nonlinear shift proportional to the mean square field. In addition, the renormalized waves exhibit a Planck like spectrum. Namely, there is equipartition of energy in the low frequency modes described by a Boltzmann distribution, followed by a linear exponential decay in the high frequency modes.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Flat-top oscillons in an expanding universe

Oscillons are extremely long lived, oscillatory, spatially localized field configurations that arise from generic initial conditions in a large number of non-linear field theories. With an eye towards their cosmological implications, we investigate their properties in an expanding universe. We (1) provide an analytic solution for one dimensional oscillons (for the models under consideration) and discuss their generalization to 3 dimensions, (2) discuss their stability against long wavelength perturbations and (3) estimate the effects of expansion on their shapes and life-times. In particular, we discuss a new, extended class of oscillons with surprisingly flat tops. We show that these flat topped oscillons are more robust against collapse instabilities in (3+1) dimensions than their usual counterparts. Unlike the solutions found in the small amplitude analysis, the width of these configurations is a non-monotonic function of their amplitudes.

0 citations0 reviewsNo code linkedOpen access