Researcher profile

John Zweck

John Zweck contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 15 - UnverifiedVerification L1Unclaimed author
3works
0followers
4topics
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

3 published item(s)

preprint2022arXiv

Spatiotemporal dynamics in a twisted, circular waveguide array

We consider the existence and spectral stability of nonlinear discrete localized solutions representing light pulses propagating in a twisted multi-core optical fiber. By considering an even number, $N$, of waveguides, we derive asymptotic expressions for solutions in which the bulk of the light intensity is concentrated as a soliton-like pulses confined to a single waveguide. The leading order terms obtained are in very good agreement with results of numerical computations. Furthermore, as in the model without temporal dispersion, when the twist parameter, $ϕ$, is given by $ϕ= π/N$, these standing waves exhibit optical suppression, in which a single waveguide remains unexcited, to leading order. Spectral computations and numerical evolution experiments suggest that these standing wave solutions are stable for values of the coupling parameter less than a critical value, at which point a spectral instability results from the collision of an internal eigenvalue with the eigenvalues at the origin. This critical value has a maximum when $ϕ= π/N$.

preprint2021arXiv

A higher order moment preserving reduction scheme for the Stochastic Weighted Particle Method

The Stochastic Weighted Particle Method (SWPM) is a Monte Carlo technique developed by Rjasanow and Wagner that generalizes Bird's Direct Simulation Monte Carlo (DSMC) method for solving the Boltzmann equation. To reduce computational cost due to the gradual increase in the number of stochastic particles in the SWPM, Rjasanow and Wagner proposed several particle reduction schemes designed to preserve specified moments of the velocity distribution. Here, we introduce an improved particle reduction scheme that preserves all moments of the velocity distribution up to the second order, as well as the raw and central heat flux both within each group of particles to be reduced and for the entire system. Furthermore, we demonstrate that with the new reduction scheme the scalar fourth-order moment can be computed more accurately at a reduced computational cost.

preprint2021arXiv

Spectral computation of low probability tails for the homogeneous Boltzmann equation

We apply the spectral-Lagrangian method of Gamba and Tharkabhushanam for solving the homogeneous Boltzmann equation to compute the low probability tails of the velocity distribution function, $f$, of a particle species. This method is based on a truncation, $Q^{\operatorname{tr}}(f,f)$, of the Boltzmann collision operator, $Q(f,f)$, whose Fourier transform is given by a weighted convolution. The truncated collision operator models the situation in which two colliding particles ignore each other if their relative speed exceeds a threshold, $g_{\text{tr}}$. We demonstrate that the choice of truncation parameter plays a critical role in the accuracy of the numerical computation of $Q$. Significantly, if $g_{\text{tr}}$ is too large, then accurate numerical computation of the weighted convolution integral is not feasible, since the decay rate and degree of oscillation of the convolution weighting function both increase as $g_{\text{tr}}$ increases. We derive an upper bound on the pointwise error between $Q$ and $Q^{\text{tr}}$, assuming that both operators are computed exactly. This bound provides some additional theoretical justification for the spectral-Lagrangian method, and can be used to guide the choice of $g_{\text{tr}}$ in numerical computations. We then demonstrate how to choose $g_{\text{tr}}$ and the numerical discretization parameters so that the computation of the truncated collision operator is a good approximation to $Q$ in the low probability tails. Finally, for several different initial conditions, we demonstrate the feasibility of accurately computing the time evolution of the velocity pdf down to probability density levels ranging from $10^{-5}$ to $10^{-9}$.