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

Jon Wilkening

Jon Wilkening contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math-phmath.APmath.MPnlin.SIphysics.flu-dynmath.DSmath.NAnlin.PSNumerical Analysis

Trust snapshot

Quick read

Trust 19 - UnverifiedVerification L1Unclaimed author
5works
0followers
9topics
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

5 published item(s)

preprint2022arXiv

Numerical Algorithms for Water Waves with Background Flow over Obstacles and Topography

We present two accurate and efficient algorithms for solving the incompressible, irrotational Euler equations with a free surface in two dimensions with background flow over a periodic, multiply-connected fluid domain that includes stationary obstacles and variable bottom topography. One approach is formulated in terms of the surface velocity potential while the other evolves the vortex sheet strength. Both methods employ layer potentials in the form of periodized Cauchy integrals to compute the normal velocity of the free surface, are compatible with arbitrary parameterizations of the free surface and boundaries, and allow for circulation around each obstacle, which leads to multiple-valued velocity potentials but single-valued stream functions. We prove that the resulting second-kind Fredholm integral equations are invertible, possibly after a physically motivated finite-rank correction. In an angle-arclength setting, we show how to avoid curve reconstruction errors that are incompatible with spatial periodicity. We use the proposed methods to study gravity-capillary waves generated by flow around several elliptical obstacles above a flat or variable bottom boundary. In each case, the free surface eventually self-intersects in a splash singularity or collides with a boundary. We also show how to evaluate the velocity and pressure with spectral accuracy throughout the fluid, including near the free surface and solid boundaries. To assess the accuracy of the time evolution, we monitor energy conservation and the decay of Fourier modes and compare the numerical results of the two methods to each other. We implement several solvers for the discretized linear systems and compare their performance. The fastest approach employs a graphics processing unit (GPU) to construct the matrices and carry out iterations of the generalized minimal residual method (GMRES).

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Harmonic Stability of Standing Water Waves

A numerical method is developed to study the stability of standing water waves and other time-periodic solutions of the free-surface Euler equations using Floquet theory. A Fourier truncation of the monodromy operator is computed by solving the linearized Euler equations about the standing wave with initial conditions ranging over all Fourier modes up to a given wave number. The eigenvalues of the truncated monodromy operator are computed and ordered by the mean wave number of the corresponding eigenfunctions, which we introduce as a method of retaining only accurately computed Floquet multipliers. The mean wave number matches up with analytical results for the zero-amplitude standing wave and is helpful in identifying which Floquet multipliers collide and leave the unit circle to form unstable eigenmodes or rejoin the unit circle to regain stability. For standing waves in deep water, most waves with crest acceleration below $A_c=0.889$ are found to be linearly stable to harmonic perturbations; however, we find several bubbles of instability at lower values of $A_c$ that have not been reported previously in the literature. We also study the stability of several new or recently discovered time-periodic gravity-capillary or gravity waves in deep or shallow water, finding several examples of large-amplitude waves that are stable to harmonic perturbations and others that are not. A new method of matching the Floquet multipliers of two nearby standing waves by solving a linear assignment problem is also proposed to track individual eigenvalues via homotopy from the zero-amplitude state to large-amplitude standing waves.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Computation of Time-Periodic Solutions of the Benjamin-Ono Equation

We present a spectrally accurate numerical method for finding non-trivial time-periodic solutions of non-linear partial differential equations. The method is based on minimizing a functional (of the initial condition and the period) that is positive unless the solution is periodic, in which case it is zero. We solve an adjoint PDE to compute the gradient of this functional with respect to the initial condition. We include additional terms in the functional to specify the free parameters, which, in the case of the Benjamin-Ono equation, are the mean, a spatial phase, a temporal phase and the real part of one of the Fourier modes at $t=0$. We use our method to study global paths of non-trivial time-periodic solutions connecting stationary and traveling waves of the Benjamin-Ono equation. As a starting guess for each path, we compute periodic solutions of the linearized problem by solving an infinite dimensional eigenvalue problem in closed form. We then use our numerical method to continue these solutions beyond the realm of linear theory until another traveling wave is reached. By experimentation with data fitting, we identify the analytical form of the solutions on the path connecting the one-hump stationary solution to the two-hump traveling wave. We then derive exact formulas for these solutions by explicitly solving the system of ODE's governing the evolution of solitons using the ansatz suggested by the numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Global paths of time-periodic solutions of the Benjamin-Ono equation connecting pairs of traveling waves

We classify all bifurcations from traveling waves to non-trivial time-periodic solutions of the Benjamin-Ono equation that are predicted by linearization. We use a spectrally accurate numerical continuation method to study several paths of non-trivial solutions beyond the realm of linear theory. These paths are found to either re-connect with a different traveling wave or to blow up. In the latter case, as the bifurcation parameter approaches a critical value, the amplitude of the initial condition grows without bound and the period approaches zero. We then prove a theorem that gives the mapping from one bifurcation to its counterpart on the other side of the path and exhibits exact formulas for the time-periodic solutions on this path. The Fourier coefficients of these solutions are power sums of a finite number of particle positions whose elementary symmetric functions execute simple orbits (circles or epicycles) in the unit disk of the complex plane. We also find examples of interior bifurcations from these paths of already non-trivial solutions, but we do not attempt to analyze their analytic structure.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Practical Error Estimates for Reynolds' Lubrication Approximation and its Higher Order Corrections

Reynolds' lubrication approximation is used extensively to study flows between moving machine parts, in narrow channels, and in thin films. The solution of Reynolds' equation may be thought of as the zeroth order term in an expansion of the solution of the Stokes equations in powers of the aspect ratio $ε$ of the domain. In this paper, we show how to compute the terms in this expansion to arbitrary order on a two-dimensional, $x$-periodic domain and derive rigorous, a-priori error bounds for the difference between the exact solution and the truncated expansion solution. Unlike previous studies of this sort, the constants in our error bounds are either independent of the function $h(x)$ describing the geometry, or depend on $h$ and its derivatives in an explicit, intuitive way. Specifically, if the expansion is truncated at order $2k$, the error is $O(ε^{2k+2})$ and $h$ enters into the error bound only through its first and third inverse moments $\int_0^1 h(x)^{-m} dx$, $m=1,3$ and via the max norms $\big\|\frac{1}{\ell!} h^{\ell-1} \partial_x^\ell h\big\|_\infty$, $1\le\ell\le2k+2$. We validate our estimates by comparing with finite element solutions and present numerical evidence that suggests that even when $h$ is real analytic and periodic, the expansion solution forms an asymptotic series rather than a convergent series.

0 citations0 reviewsNo code linkedOpen access