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

Steve Shkoller

Steve Shkoller appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.APmath-phmath.MPphysics.comp-phphysics.flu-dynmath.NA

Catalog footprint

What is connected

23works
6topics
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

23 published item(s)

preprint2020arXiv

Affine motion of 2d incompressible fluids surrounded by vacuum and flows in ${\rm SL}(2,{\mathbb R})$

The affine motion of two-dimensional (2d) incompressible fluids surrounded by vacuum can be reduced to a completely integrable and globally solvable Hamiltonian system of ordinary differential equations for the deformation gradient in ${\rm SL}(2,{\mathbb R})$. In the case of perfect fluids, the motion is given by geodesic flow in ${\rm SL}(2,{\mathbb R})$ with the Euclidean metric, while for magnetically conducting fluids (MHD), the motion is governed by a harmonic oscillator in ${\rm SL}(2,{\mathbb R})$. A complete classification of the dynamics is given including rigid motions, rotating eddies with stable and unstable manifolds, and solutions with vanishing pressure. For perfect fluids, the displacement generically becomes unbounded, as $t\to\pm\infty$. For MHD, solutions are bounded and generically quasi-periodic and recurrent.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Formation of point shocks for 3D compressible Euler

We consider the 3D isentropic compressible Euler equations with the ideal gas law. We provide a constructive proof of shock formation from smooth initial datum of finite energy, with no vacuum regions, with nontrivial vorticity present at the shock, and under no symmetry assumptions. We prove that for an open set of Sobolev-class initial data which are a small $L^ \infty $ perturbation of a constant state, there exist smooth solutions to the Euler equations which form a generic stable shock in finite time. The blow up time and location can be explicitly computed, and solutions at the blow up time are smooth except for a single point, where they are of cusp-type with Hölder $C^ {\frac{1}{3}}$ regularity. Our proof is based on the use of modulated self-similar variables that are used to enforce a number of constraints on the blow up profile, necessary to establish the stability in self-similar variables of the generic shock profile.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Shock formation and vorticity creation for 3d Euler

We analyze the shock formation process for the 3d non-isentropic Euler equations with the ideal gas law, in which sounds waves interact with entropy waves to produce vorticity. Building on our theory for isentropic flows in [3,4], we give a constructive proof of shock formation from smooth initial data. Specifically, we prove that there exist smooth solutions to the non-isentropic Euler equations which form a generic stable shock with explicitly computable blowup time, location, and direction. This is achieved by establishing the asymptotic stability of a generic shock profile in modulated self-similar variables, controlling the interaction of wave families via: (i) pointwise bounds along Lagrangian trajectories, (ii) geometric vorticity structure, and (iii) high-order energy estimates in Sobolev spaces.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A multiscale model for Rayleigh-Taylor and Richtmyer-Meshkov instabilities

We develop a novel multiscale model of interface motion for the Rayleigh-Taylor instability (RTI) and Richtmyer-Meshkov instability (RMI) for two-dimensional, inviscid, compressible flows with vorticity, which yields a fast-running numerical algorithm that produces both qualitatively and quantitatively similar results to a resolved gas dynamics code, while running approximately two orders of magnitude (in time) faster. Our multiscale model is founded upon a new compressible-incompressible decomposition of the velocity field $u=v+w$. The incompressible component $w$ of the velocity is also irrotational and is solved using a new asymptotic model of the Birkhoff-Rott singular integral formulation of the incompressible Euler equations, which reduces the problem to one spatial dimension. This asymptotic model, called the higher-order $z$-model, is derived using small nonlocality as the asymptotic parameter, allows for interface turn-over and roll-up, and yields a significant simplification for the equation describing the evolution of the amplitude of vorticity. This incompressible component $w$ of the velocity controls the small scale structures of the interface and can be solved efficiently on fine grids. Meanwhile, the compressible component of the velocity $v$ remains continuous near contact discontinuities and can be computed on relatively coarse grids, while receiving subgrid scale information from $w$. We first validate the incompressible higher-order $z$-model by comparison with classical RTI experiments as well as full point vortex simulations. We then consider both the RTI and the RMI problems for our multiscale model of compressible flow with vorticity, and show excellent agreement with our high-resolution gas dynamics solutions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A model for Rayleigh-Taylor mixing and interface turn-over

We first develop a new mathematical model for two-fluid interface motion, subjected to the Rayleigh-Taylor (RT) instability in two-dimensional fluid flow, which in its simplest form, is given by $ h_{tt}(α,t) = A g\, Λh - \fracσ{ρ^++ρ^-} Λ^3 h - A \partial_α(H h_t h_t) $, where $Λ= H \partial_ α$ and $H$ denotes the Hilbert transform. In this so-called $h$-model, $A$ is the Atwood number, $g$ is the acceleration, $ σ$ is surface tension, and $ρ^\pm$ denotes the densities of the two fluids. Under a certain stability condition, we prove that this so-called $h$-model is both locally and globally well-posed. Numerical simulations of the $h$-model show that the interface can quickly grow due to nonlinearity, and then stabilize when the lighter fluid is on top of the heavier fluid and acceleration is directed downward. In the unstable case of a heavier fluid being supported by the lighter fluid, we find good agreement for the growth of the mixing layer with experimental data in the "rocket rig" experiment of Read of Youngs. We then derive an RT interface model with a general parameterization $z(α,t)$ such that $ z_{tt}= Λ\bigg{[}\frac{A}{|\partial_αz|^2}H\left(z_t\cdot (\partial_αz)^\perp H(z_t\cdot (\partial_αz)^\perp)\right) + A g z_2 \bigg{]} \frac{(\partial_αz)^\perp}{|\partial_αz|^2} +z_t\cdot (\partial_αz)^\perp\left(\frac{(\partial_αz_t)^\perp}{|\partial_αz|^2}-\frac{(\partial_αz)^\perp 2(\partial_αz\cdot \partial_αz_t)}{|\partial_αz|^4}\right)$. This more general RT $z$-model allows for interface turn-over. Numerical simulations of the $z$-model show an even better agreement with the predicted mixing layer growth for the "rocket rig" experiment.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Local well-posedness and Global stability of the Two-Phase Stefan problem

The two-phase Stefan problem describes the temperature distribution in a homogeneous medium undergoing a phase transition such as ice melting to water. This is accomplished by solving the heat equation on a time-dependent domain, composed of two regions separated by an a priori unknown moving boundary which is transported by the difference (or jump) of the normal derivatives of the temperature in each phase. We establish local-in-time well-posedness and a global-in-time stability result for arbitrary sufficiently smooth domains and small initial temperatures. To this end, we develop a higher-order energy with natural weights adapted to the problem and combine it with Hopf-type inequalities. This extends the previous work by Hadzic and Shkoller [31,32] on the one-phase Stefan problem to the setting of two-phase problems, and simplifies the proof significantly.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the impossibility of finite-time splash singularities for vortex sheets

In fluid dynamics, an interface splash singularity occurs when a locally smooth interface self-intersects in finite time. By means of elementary arguments, we prove that such a singularity cannot occur in finite time for vortex sheet evolution, i.e. for the two-phase incompressible Euler equations. We prove this by contradiction; we assume that a splash singularity does indeed occur in finite time. Based on this assumption, we find precise blow-up rates for the components of the velocity gradient which, in turn, allow us to characterize the geometry of the evolving interface just prior to self-intersection. The constraints on the geometry then lead to an impossible outcome, showing that our assumption of a finite-time splash singularity was false.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Well-posedness and decay to equilibrium for the Muskat problem with discontinuous permeability

We first prove local-in-time well-posedness for the Muskat problem, modeling fluid flow in a two-dimensional inhomogeneous porous media. The permeability of the porous medium is described by a step function, with a jump discontinuity across the fixed-in-time curve $(x_1,-1+f(x_1))$, while the interface separating the fluid from the vacuum region is given by the time-dependent curve $(x_1,h(x_1,t))$. Our estimates are based on a new methodology that relies upon a careful study of the PDE system, coupling Darcy's law and incompressibility of the fluid, rather than the analysis of the singular integral contour equation for the interface function $h$. We are able to develop an existence theory for any initial interface given by $h_0 \in H^2$ and any permeability curve-of-discontinuity that is given by $f \in H^{2.5}$. In particular, our method allows for both curves to have (pointwise) unbounded curvature. In the case that the permeability discontinuity is the set $f=0$, we prove global existence and decay to equilibrium for small initial data. This decay is obtained using a new energy-energy dissipation inequality that couples tangential derivatives of the velocity in the bulk of the fluid with the curvature of the interface. To the best of our knowledge, this is the first global existence result for the Muskat problem with discontinuous permeability.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Well-posedness for the classical Stefan problem and the zero surface tension limit

We develop a framework for a unified treatment of well-posedness for the Stefan problem with or without surface tension. In the absence of surface tension, we establish well-posedness in Sobolev spaces for the classical Stefan problem. We introduce a new velocity variable which extends the velocity of the moving free-boundary into the interior domain. The equation satisfied by this velocity is used for the analysis in place of the heat equation satisfied by the temperature. Solutions to the classical Stefan problem are then constructed as the limit of solutions to a carefully chosen sequence of approximations to the velocity equation, in which the moving free-boundary is regularized and the boundary condition is modified in a such a way as to preserve the basic nonlinear structure of the original problem. With our methodology, we simultaneously find the required stability condition for well-posedness and obtain new estimates for the regularity of the moving free-boundary. Finally, we prove that solutions of the Stefan problem with positive surface tension $σ$ converge to solutions of the classical Stefan problem as $σ\to 0$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A priori estimates for solutions to the relativistic Euler equations with a moving vacuum boundary

We study the relativistic Euler equations on the Minkowski spacetime background. We make assumptions on the equation of state and the initial data that are relativistic analogs of the well-known physical vacuum boundary condition, which has played an important role in prior work on the non-relativistic compressible Euler equations. Our main result is the derivation, relative to Lagrangian (also known as co-moving) coordinates, of local-in-time a priori estimates for the solution. The solution features a fluid-vacuum boundary, transported by the fluid four-velocity, along which the hyperbolicity of the equations degenerates. In this context, the relativistic Euler equations are equivalent to a degenerate quasilinear hyperbolic wave-map-like system that cannot be treated using standard energy methods.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Global stability of steady states in the classical Stefan problem

The classical one-phase Stefan problem (without surface tension) allows for a continuum of steady state solutions, given by an arbitrary (but sufficiently smooth) domain together with zero temperature. We prove global-in-time stability of such steady states, assuming a sufficient degree of smoothness on the initial domain, but without any a priori restriction on the convexity properties of the initial shape. This is an extension of our previous result [28] in which we studied nearly spherical shapes.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On the splash singularity for the free-surface of a Navier-Stokes fluid

In fluid dynamics, an interface splash singularity occurs when a locally smooth interface self-intersects in finite time. We prove that for $d$-dimensional flows, $d=2$ or $3$, the free-surface of a viscous water wave, modeled by the incompressible Navier-Stokes equations with moving free-boundary, has a finite-time splash singularity. In particular, we prove that given a sufficiently smooth initial boundary and divergence-free velocity field, the interface will self-intersect in finite time.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Regularity of the velocity field for Euler vortex patch evolution

We consider the vortex patch problem for both the 2-D and 3-D incompressible Euler equations. In 2-D, we prove that for vortex patches with $H^{k-0.5}$ Sobolev-class contour regularity, $k \ge 4$, the velocity field on both sides of the vortex patch boundary has $H^k$ regularity for all time. In 3-D, we establish existence of solutions to the vortex patch problem on a finite-time interval $[0,T]$, and we simultaneously establish the $H^{k-0.5}$ regularity of the two-dimensional vortex patch boundary, as well as the $H^k$ regularity of the velocity fields on both sides of vortex patch boundary, for $k \ge 3$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Solvability and regularity for an elliptic system prescribing the curl, divergence, and partial trace of a vector field on Sobolev-class domains

We provide a self-contained proof of the solvability and regularity of a Hodge-type elliptic system, wherein the divergence and curl of a vector field are prescribed in an open, bounded, Sobolev-class domain, and either the normal component or the tangential components of the vector field are prescribed on the boundary. The proof is based on a regularity theory for vector elliptic equations set on Sobolev-class domains and with Sobolev-class coefficients.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Well-posedness of the Muskat problem with $H^2$ initial data

We study the dynamics of the interface between two incompressible fluids in a two-dimensional porous medium whose flow is modeled by the Muskat equations. For the two-phase Muskat problem, we establish global well-posedness and decay to equilibrium for small $H^2$ perturbations of the rest state. For the one-phase Muskat problem, we prove local well-posedness for $H^2$ initial data of arbitrary size. Finally, we show that solutions to the Muskat equations instantaneously become infinitely smooth.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Global stability and decay for the classical Stefan problem

The classical one-phase Stefan problem describes the temperature distribution in a homogeneous medium undergoing a phase transition, such as ice melting to water. This is accomplished by solving the heat equation on a time-dependent domain whose boundary is transported by the normal derivative of the temperature along the evolving and a priori unknown free-boundary. We establish a global-in-time stability result for nearly spherical geometries and small temperatures, using a novel hybrid methodology, which combines energy estimates, decay estimates, and Hopf-type inequalities.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On the finite-time splash and splat singularities for the 3-D free-surface Euler equations

We prove that the 3-D free-surface incompressible Euler equations with regular initial geometries and velocity fields have solutions which can form a finite-time "splash" (or "splat") singularity first introduced in [9], wherein the evolving 2-D hypersurface, the moving boundary of the fluid domain, self-intersects at a point (or on surface). Such singularities can occur when the crest of a breaking wave falls unto its trough, or in the study of drop impact upon liquid surfaces. Our approach is founded upon the Lagrangian description of the free-boundary problem, combined with a novel approximation scheme of a finite collection of local coordinate charts; as such we are able to analyze a rather general set of geometries for the evolving 2-D free-surface of the fluid. We do not assume the fluid is irrotational, and as such, our method can be used for a number of other fluid interface problems, including compressible flows, plasmas, as well as the inclusion of surface tension effects.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Space-time Smooth Artificial Viscosity Method For Nonlinear Conservation Laws

We introduce a new methodology for adding localized, space-time smooth, artificial viscosity to nonlinear systems of conservation laws which propagate shock waves, rarefactions, and contact discontinuities, which we call the $C$-method. We shall focus our attention on the compressible Euler equations in one space dimension. The novel feature of our approach involves the coupling of a linear scalar reaction-diffusion equation to our system of conservation laws, whose solution $C(x,t)$ is the coefficient to an additional (and artificial) term added to the flux, which determines the location, localization, and strength of the artificial viscosity. Near shock discontinuities, $C(x,t)$ is large and localized, and transitions smoothly in space-time to zero away from discontinuities. Our approach is a provably convergent, spacetime-regularized variant of the original idea of Richtmeyer and Von Neumann, and is provided at the level of the PDE, thus allowing a host of numerical discretization schemes to be employed. We demonstrate the effectiveness of the $C$-method with three different numerical implementations and apply these to a collection of classical problems: the Sod shock-tube, the Osher-Shu shock-tube, the Woodward-Colella blast wave and the Leblanc shock-tube. First, we use a classical continuous finite-element implementation using second-order discretization in both space and time, FEM-C. Second, we use a simplified WENO scheme within our $C$-method framework, WENO-C. Third, we use WENO with the Lax-Friedrichs flux together with the $C$-equation, and call this WENO-LF-C. All three schemes yield higher-order discretization strategies, which provide sharp shock resolution with minimal overshoot and noise, and compare well with higher-order WENO schemes that employ approximate Riemann solvers, outperforming them for the difficult Leblanc shock tube experiment.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Global existence and decay for solutions of the Hele-Shaw flow with injection

We study the global existence and decay to spherical equilibrium of Hele-Shaw flows with surface tension. We prove that without injection of fluid, perturbations of the sphere decay to zero exponentially fast. On the other hand, with a time-dependent rate of fluid injection into the Hele-Shaw cell, the distance from the moving boundary to an expanding sphere (with time-dependent radius) also decays to zero but with an algebraic rate, which depends on the injection rate of the fluid.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Well-posedness of the free-boundary compressible 3-D Euler equations with surface tension and the zero surface tension limit

We prove that the 3-D compressible Euler equations with surface tension along the moving free-boundary are well-posed. Specifically, we consider isentropic dynamics and consider an equation of state, modeling a liquid, given by Courant and Friedrichs as $p(ρ) = αρ^ γ- β$ for consants $γ>1$ and $ α, β> 0$. The analysis is made difficult by two competing nonlinearities associated with the potential energy: compression in the bulk, and surface area dynamics on the free-boundary. Unlike the analysis of the incompressible Euler equations, wherein boundary regularity controls regularity in the interior, the compressible Euler equation require the additional analysis of nonlinear wave equations generating sound waves. An existence theory is developed by a specially chosen parabolic regularization together with the vanishing viscosity method. The artificial parabolic term is chosen so as to be asymptotically consistent with the Euler equations in the limit of zero viscosity. Having solutions for the positive surface tension problem, we proceed to obtain a priori estimates which are independent of the surface tension parameter. This requires choosing initial data which satisfy the Taylor sign condition. By passing to the limit of zero surface tension, we prove the well-posedness of the compressible Euler system without surface on the free-boundary, and without derivative loss.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Well-posedness in smooth function spaces for the moving-boundary 3-D compressible Euler equations in physical vacuum

We prove well-posedness for the 3-D compressible Euler equations with moving physical vacuum boundary, with an equation of state given by the so-called gamma gas-law for gamma > 1. The physical vacuum singularity requires the sound speed c to go to zero as the square-root of the distance to the moving boundary, and thus creates a degenerate and characteristic hyperbolic free-boundary system wherein the density vanishes on the free-boundary, the uniform Kreiss--Lopatinskii condition is violated, and manifest derivative loss ensues. Nevertheless, we are able to establish the existence of unique solutions to this system on a short time-interval, which are smooth (in Sobolev spaces) all the way to the moving boundary, and our estimates have no derivative loss with respect to initial data. Our proof is founded on an approximation of the Euler equations by a degenerate parabolic regularization obtained from a specific choice of a degenerate artificial viscosity term, chosen to preserve as much of the geometric structure of the Euler equations as possible. We first construct solutions to this degenerate parabolic regularization using a new higher-order Hardy-type inequality; we then establish estimates for solutions to this degenerate parabolic system which are independent of the artificial viscosity parameter. Solutions to the compressible Euler equations are found in the limit as the artificial viscosity tends to zero. Our regular solutions can be viewed as degenerate viscosity solutions. Out methodology can be applied to many other systems of degenerate and characteristic hyperbolic systems of conservation laws.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A priori estimates for the free-boundary 3-D compressible Euler equations in physical vacuum

We prove a priori estimates for the three-dimensional compressible Euler equations with moving {\it physical} vacuum boundary, with an equation of state given by $p(ρ) = C_γρ^γ$ for $γ>1$. The vacuum condition necessitates the vanishing of the pressure, and hence density, on the dynamic boundary, which creates a degenerate and characteristic hyperbolic {\it free-boundary} system to which standard methods of symmetrizable hyperbolic equations cannot be applied.

0 citations0 reviewsNo code linkedOpen access