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

Gideon Simpson

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

ResearcherUnclaimed source record
math.APmath.NAnlin.PScond-mat.mtrl-sciNumerical Analysisphysics.geo-phmath-phmath.MPmath.PRmath.SPBiological Physicscond-mat.stat-mechmath.DSphysics.chem-phphysics.comp-phphysics.flu-dynphysics.opticsphysics.plasm-ph

Catalog footprint

What is connected

30works
18topics
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

30 published item(s)

preprint2022arXiv

A Numerical Method for a Nonlocal Diffusion Equation with Additive Noise

We consider a nonlocal evolution equation representing the continuum limit of a large ensemble of interacting particles on graphs forced by noise. The two principle ingredients of the continuum model are a nonlocal term and Q-Wiener process describing the interactions among the particles in the network and stochastic forcing respectively. The network connectivity is given by a square integrable function called a graphon. We prove that the initial value problem for the continuum model is well-posed. Further, we construct a semidiscrete (discrete in space and continuous in time) and a fully discrete schemes for the nonlocal model. The former is obtained by a discontinuous Galerkin method and the latter is based on further discretizing time using the Euler-Maruyama method. We prove convergence and estimate the rate of convergence in each case. For the semidiscrete scheme, the rate of convergence estimate is expressed in terms of the regularity of the graphon, Q-Wiener process, and the initial data. We work in generalized Lipschitz spaces, which allows to treat models with data of lower regularity. This is important for applications as many interesting types of connectivity including small-world and power-law are expressed by graphons that are not smooth. The error analysis of the fully discrete scheme, on the other hand, reveals that for some models common in applied science, one has a higher speed of convergence than that predicted by the standard estimates for the Euler-Maruyama method. The rate of convergence analysis is supplemented with detailed numerical experiments, which are consistent with our analytical results. As a by-product, this work presents a rigorous justification for taking continuum limit for a large class of interacting dynamical systems on graphs subject to noise.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Iterate Averaging, the Kalman Filter, and 3DVAR for Linear Inverse Problem

It has been proposed that classical filtering methods, like the Kalman filter and 3DVAR, can be used to solve linear statistical inverse problems. In the work of Iglesias, Lin, Lu, & Stuart (2017), error estimates were obtained for this approach. By optimally tuning a regularization parameter in the filters, the authors were able to show that the mean squared error could be systematically reduced. Building on the aforementioned work of Iglesias, Lin, Lu, & Stuart, we prove that by (i) considering the problem in a weaker norm and (ii) applying simple iterate averaging of the filter output, 3DVAR will converge in mean square, unconditionally on the choice of parameter. Without iterate averaging, 3DVAR cannot converge by running additional iterations with a fixed choice of parameter. We also establish that the Kalman filter's performance in this setting cannot be improved through iterate averaging. We illustrate our results with numerical experiments that suggest our convergence rates are sharp.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Sampling from Rough Energy Landscapes

We examine challenges to sampling from Boltzmann distributions associated with multiscale energy landscapes. The multiscale features, or "roughness," corresponds to highly oscillatory, but bounded, perturbations of a smooth landscape. Through a combination of numerical experiments and analysis we demonstrate that the performance of Metropolis Adjusted Langevin Algorithm can be severely attenuated as the roughness increases. In contrast, we prove that Random Walk Metropolis is insensitive to such roughness. We also formulate two alternative sampling strategies that incorporate large scale features of the energy landscape, while resisting the impact of fine scale roughness; these also outperform Random Walk Metropolis. Numerical experiments on these landscapes are presented that confirm our predictions. Open questions and numerical challenges are also highlighted.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Transient probability currents provide upper and lower bounds on non-equilibrium steady-state currents in the Smoluchowski picture

Probability currents are fundamental in characterizing the kinetics of non-equilibrium processes. Notably, the steady-state current $J_{ss}$ for a source-sink system can provide the exact mean-first-passage time (MFPT) for the transition from source to sink. Because transient non-equilibrium behavior is quantified in some modern path sampling approaches, such as the "weighted ensemble" strategy, there is strong motivation to determine bounds on $J_{ss}$ -- and hence on the MFPT -- as the system evolves in time. Here we show that $J_{ss}$ is bounded from above and below by the maximum and minimum, respectively, of the current as a function of the spatial coordinate at any time $t$ for one-dimensional systems undergoing over-damped Langevin (i.e., Smoluchowski) dynamics and for higher-dimensional Smoluchowski systems satisfying certain assumptions when projected onto a single dimension. These bounds become tighter with time, making them of potential practical utility in a scheme for estimating $J_{ss}$ and the long-timescale kinetics of complex systems. Conceptually, the bounds result from the fact that extrema of the transient currents relax toward the steady-state current.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Existence and Stability Properties of Radial Bound States for Schrödinger-Poisson with an External Coulomb Potential in Three Space Dimensions

We consider radial solutions to the Schrödinger-Poisson system in three dimensions with an external smooth potential with Coulomb-like decay. Such a system can be viewed as a model for the interaction of dark matter with a bright matter background in the non-relativistic limit. We find that there are infinitely many critical points of the Hamiltonian, subject to fixed mass, and that these bifurcate from solutions to the associated linear problem at zero mass. As a result, each branch has a different topological character defined by the number of zeros of the radial states. We construct numerical approximations to these nonlinear states along the first several branches. The solution branches can be continued, numerically, to large mass values, where they become asymptotic, under a rescaling, to those of the Schrödinger-Poisson problem with no external potential. Our numerical computations indicate that the ground state is orbitally stable, while the excited states are linearly unstable for sufficiently large mass.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Local structure of singular profiles for a Derivative Nonlinear Schrödinger Equation

The Derivative Nonlinear Schrödinger equation is an $L^2$-critical nonlinear dispersive equation model for Alfvén waves in space plasmas. Recent numerical studies on an $L^2$-supercritical extension of this equation provide evidence of finite time singularities. Near the singular point, the solution is described by a universal profile that solves a nonlinear elliptic eigenvalue problem depending only on the strength of the nonlinearity. In the present work, we describe the deformation of the profile and its parameters near criticality, combining asymptotic analysis and numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Theoretical Examination of Diffusive Molecular Dynamics

Diffusive molecular dynamics is a novel model for materials with atomistic resolution that can reach diffusive time scales. The main ideas of diffusive molecular dynamics are to first minimize an approximate variational Gaussian free energy of the system with respect to the mean atomic coordinates (averaging over many vibrational periods), and to then to perform a diffusive step where atoms and vacancies (or two species in a binary alloy) flow on a diffusive time scale via a master equation. We present a mathematical framework for studying this algorithm based upon relative entropy, or Kullback-Leibler divergence. This adds flexibility in how the algorithm is implemented and interpreted. We then compare our formulation, relying on relative entropy and absolute continuity of measures, to existing formulations. The main difference amongst the equations appears in a model for vacancy diffusion, where additional entropic terms appear in our development.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Asymptotic Stability of high-dimensional Zakharov-Kuznetsov solitons

We prove that solitons (or solitary waves) of the Zakharov-Kuznetsov (ZK) equation, a physically relevant high dimensional generalization of the Korteweg-de Vries (KdV) equation appearing in Plasma Physics, and having mixed KdV and nonlinear Schrödinger (NLS) dynamics, are strongly asymptotically stable in the energy space in the physical region. We also prove that the sum of well-arranged solitons is stable in the same space. Orbital stability of ZK solitons is well-known since the work of de Bouard. Our proofs follow the ideas by Martel and Martel and Merle, applied for generalized KdV equations in one dimension. In particular, we extend to the high dimensional case several monotonicity properties for suitable half-portions of mass and energy; we also prove a new Liouville type property that characterizes ZK solitons, and a key Virial identity for the linear and nonlinear part of the ZK dynamics, obtained independently of the mixed KdV-NLS dynamics. This last Virial identity relies on a simple sign condition, which is numerically tested for the two and three dimensional cases, with no additional spectral assumptions required. Possible extensions to higher dimensions and different nonlinearities could be obtained after a suitable local well-posedness theory in the energy space, and the verification of a corresponding sign condition.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Generalized Parallel Replica Dynamics

Metastability is a common obstacle to performing long molecular dynamics simulations. Many numerical methods have been proposed to overcome it. One method is parallel replica dynamics, which relies on the rapid convergence of the underlying stochastic process to a quasi-stationary distribution. Two requirements for applying parallel replica dynamics are knowledge of the time scale on which the process converges to the quasi-stationary distribution and a mechanism for generating samples from this distribution. By combining a Fleming-Viot particle system with convergence diagnostics to simultaneously identify when the process converges while also generating samples, we can address both points. This variation on the algorithm is illustrated with various numerical examples, including those with entropic barriers and the 2D Lennard-Jones cluster of seven atoms.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Algorithms for Kullback-Leibler Approximation of Probability Measures in Infinite Dimensions

In this paper we study algorithms to find a Gaussian approximation to a target measure defined on a Hilbert space of functions; the target measure itself is defined via its density with respect to a reference Gaussian measure. We employ the Kullback-Leibler divergence as a distance and find the best Gaussian approximation by minimizing this distance. It then follows that the approximate Gaussian must be equivalent to the Gaussian reference measure, defining a natural function space setting for the underlying calculus of variations problem. We introduce a computational algorithm which is well-adapted to the required minimization, seeking to find the mean as a function, and parameterizing the covariance in two different ways: through low rank perturbations of the reference covariance; and through Schrödinger potential perturbations of the inverse reference covariance. Two applications are shown: to a nonlinear inverse problem in elliptic PDEs, and to a conditioned diffusion process. We also show how the Gaussian approximations we obtain may be used to produce improved pCN-MCMC methods which are not only well-adapted to the high-dimensional setting, but also behave well with respect to small observational noise (resp. small temperatures) in the inverse problem (resp. conditioned diffusion).

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Kullback-Leibler Approximation for Probability Measures on Infinite Dimensional Spaces

In a variety of applications it is important to extract information from a probability measure $μ$ on an infinite dimensional space. Examples include the Bayesian approach to inverse problems and possibly conditioned) continuous time Markov processes. It may then be of interest to find a measure $ν$, from within a simple class of measures, which approximates $μ$. This problem is studied in the case where the Kullback-Leibler divergence is employed to measure the quality of the approximation. A calculus of variations viewpoint is adopted and the particular case where $ν$ is chosen from the set of Gaussian measures is studied in detail. Basic existence and uniqueness theorems are established, together with properties of minimising sequences. Furthermore, parameterisation of the class of Gaussians through the mean and inverse covariance is introduced, the need for regularisation is explained, and a regularised minimisation is studied in detail. The calculus of variations framework resulting from this work provides the appropriate underpinning for computational algorithms.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Local Existence Theory for Derivative Nonlinear Schrödinger Equations with Non-Integer Power Nonlinearities

We study a derivative nonlinear Schrödinger equation, allowing non-integer powers in the nonlinearity, $|u|^{2σ} u_x$. Making careful use of the energy method, we are able to establish short-time existence of solutions with initial data in the energy space, $H^1$. For more regular initial data, we establish not just existence of solutions, but also well-posedness of the initial value problem. These results hold for real-valued $σ\geq 1,$ while prior existence results in the literature require integer-valued $σ$ or $σ$ sufficiently large ($σ\geq 5/2$), or use higher-regularity function spaces.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Petviashvilli's Method for the Dirichlet Problem

We examine the Petviashvilli method for solving the equation $ ϕ- Δϕ= |ϕ|^{p-1} ϕ$ on a bounded domain $Ω\subset \mathbb{R}^d$ with Dirichlet boundary conditions. We prove a local convergence result, using spectral analysis, akin to the result for the problem on $\mathbb{R}$ by Pelinovsky & Stepanyants, 2004. We also prove a global convergence result by generating a suite of nonlinear inequalities for the iteration sequence, and we show that the sequence has a natural energy that decreases along the sequence.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The parallel replica method for simulating long trajectories of Markov chains

The parallel replica dynamics, originally developed by A.F. Voter, efficiently simulates very long trajectories of metastable Langevin dynamics. We present an analogous algorithm for discrete time Markov processes. Such Markov processes naturally arise, for example, from the time discretization of a continuous time stochastic dynamics. Appealing to properties of quasistationary distributions, we show that our algorithm reproduces exactly (in some limiting regime) the law of the original trajectory, coarsened over the metastable states.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Dynamics near a minimal-mass soliton for a Korteweg-de Vries equation

We study soliton solutions to a generalized Korteweg - de Vries (KdV) equation with a saturated nonlinearity, following the line of inquiry of the authors for the nonlinear Schrödinger equation (NLS). KdV with such a nonlinearity is known to possess a minimal-mass soliton. We consider a small perturbation of a minimal-mass soliton and identify a system of ODEs, which models the behavior of the perturbation for short times. This connects nicely to a work of Comech, Cuccagna & Pelinovsky (2007). These ODEs form a simple dynamical system with a single unstable hyperbolic fixed point with two possible dynamical outcomes. A particular feature of the dynamics are that they are non-oscillatory. This distinguishes the KdV problem from the analogous NLS one.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Focusing Singularity in a Derivative Nonlinear Schrödinger Equation

We present a numerical study of a derivative nonlinear Schrödinger equation with a general power nonlinearity, $|ψ|^{2σ}ψ_x$. In the $L^2$-supercritical regime, $σ>1$, our simulations indicate that there is a finite time singularity. We obtain a precise description of the local structure of the solution in terms of blowup rate and asymptotic profile, in a form similar to that of the nonlinear Schrödinger equation with supercritical power law nonlinearity.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Behavior of a Model Dynamical System with Applications to Weak Turbulence

We experimentally explore solutions to a model Hamiltonian dynamical system derived in Colliander et al., 2012, to study frequency cascades in the cubic defocusing nonlinear Schrödinger equation on the torus. Our results include a statistical analysis of the evolution of data with localized amplitudes and random phases, which supports the conjecture that energy cascades are a generic phenomenon. We also identify stationary solutions, periodic solutions in an associated problem and find experimental evidence of hyperbolic behavior. Many of our results rely upon reframing the dynamical system using a hydrodynamic formulation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Numerical Analysis of Parallel Replica Dynamics

Parallel replica dynamics is a method for accelerating the computation of processes characterized by a sequence of infrequent events. In this work, the processes are governed by the overdamped Langevin equation. Such processes spend much of their time about the minima of the underlying potential, occasionally transitioning into different basins of attraction. The essential idea of parallel replica dynamics is that the exit time distribution from a given well for a single process can be approximated by the minimum of the exit time distributions of $N$ independent identical processes, each run for only 1/N-th the amount of time. While promising, this leads to a series of numerical analysis questions about the accuracy of the exit distributions. Building upon the recent work in Le Bris et al., we prove a unified error estimate on the exit distributions of the algorithm against an unaccelerated process. Furthermore, we study a dephasing mechanism, and prove that it will successfully complete.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Analytical and Numerical Results on the Positivity of Steady State Solutions of a Thin Film Equation

We consider an equation for a thin-film of fluid on a rotating cylinder and present several new analytical and numerical results on steady state solutions. First, we provide an elementary proof that both weak and classical steady states must be strictly positive so long as the speed of rotation is nonzero. Next, we formulate an iterative spectral algorithm for computing these steady states. Finally, we explore a non-existence inequality for steady state solutions from the recent work of Chugunova, Pugh, & Taranets.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Broad Band Solitons in a Periodic and Nonlinear Maxwell System

We consider the one-dimensional Maxwell equations with low contrast periodic linear refractive index and weak Kerr nonlinearity. In this context, wave packet initial conditions with a single carrier frequency excite infinitely many resonances. On large but finite time-scales, the coupled evolution of backward and forward waves is governed by nonlocal equations of resonant nonlinear geometrical optics. For the special class of solutions which are periodic in the fast phase, these equations are equivalent to an infinite system of nonlinear coupled mode equations, the so called it extended nonlinear coupled mode equations, xNLCME. Numerical studies support the existence of long-lived spatially localized coherent structures, featuring a slowly varying envelope and a train of carrier shocks. In this paper we explore, by analytical, asymptotic and numerical methods, the existence and properties of spatially localized structures of the xNLCME system, which arises for a refractive index profile consisting of periodic array of Dirac delta functions. We consider the limit of small amplitude solutions with frequencies near a band-edge. In this case, stationary xNLCME is well-approximated by an infinite system of coupled, stationary, nonlinear Schrödinger equations, the extended nonlinear Schrödinger system, xNLS. We embed xNLS in a one-parameter family of equations, xNLS$^ε$, which interpolates between infinitely many decoupled NLS equations ($ε=0$) and xNLS ($ε=1$). Using bifurcation methods we show existence of solutions for a range of $ε\in(-ε_0,ε_0)$ and, by a numerical continuation method, establish the continuation of certain branches all the way to $ε=1$. Finally, we perform time-dependent simulations of truncated xNLCME and find the small-amplitude solitons to be robust to both numerical errors and the NLS approximation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Embedded Eigenvalues and the Nonlinear Schrodinger Equation

A common challenge to proving asymptotic stability of solitary waves is understanding the spectrum of the operator associated with the linearized flow. The existence of eigenvalues can inhibit the dispersive estimates key to proving stability. Following the work of Marzuola & Simpson, we prove the absence of embedded eigenvalues for a collection of nonlinear Schrodinger equations, including some one and three dimensional supercritical equations, and the three dimensional cubic-quintic equation. Our results also rule out nonzero eigenvalues within the spectral gap and, in 3D, endpoint resonances. The proof is computer assisted as it depends on the sign of certain inner products which do not readily admit analytic representations. Our source code is available for verification at http://www.math.toronto.edu/simpson/files/spec_prop_asad_simpson_code.zip.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Ill-posedness of degenerate dispersive equations

In this article we provide numerical and analytical evidence that some degenerate dispersive partial differential equations are ill-posed. Specifically we study the K(2,2) equation $u_t = (u^2)_{xxx} + (u^2)_{x}$ and the "degenerate Airy" equation $u_t = 2 u u_{xxx}$. For K(2,2) our results are computational in nature: we conduct a series of numerical simulations which demonstrate that data which is very small in $H^2$ can be of unit size at a fixed time which is independent of the data's size. For the degenerate Airy equation, our results are fully rigorous: we prove the existence of a compactly supported self-similar solution which, when combined with certain scaling invariances, implies ill-posedness (also in $H^2$).

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Coherent Structures and Carrier Shocks in the Nonlinear Periodic Maxwell Equations

We consider the one-dimensional propagation of electromagnetic waves in a weakly nonlinear and low-contrast spatially inhomogeneous medium with no energy dissipation. We focus on the case of a periodic medium, in which dispersion enters only through the (Floquet-Bloch) spectral band dispersion associated with the periodic structure; chromatic dispersion (time-nonlocality of the polarization) is neglected. Numerical simulations show that for initial conditions of wave-packet type (a plane wave of fixed carrier frequency multiplied by a slow varying, spatially localized function) very long-lived spatially localized coherent soliton-like structures emerge, whose character is that of a slowly varying envelope of a train of shocks. We call this structure an envelope carrier-shock train. The structure of the solution violates the oft-assumed nearly monochromatic wave packet structure, whose envelope is governed by the nonlinear coupled mode equations (NLCME). The inconsistency and inaccuracy of NLCME lies in the neglect of all (infinitely many) resonances except for the principle resonance induced by the initial carrier frequency. We derive, via a nonlinear geometrical optics expansion, a system of nonlocal integro-differential equations governing the coupled evolution of backward and forward propagating waves. These equations incorporate effects of all resonances. In a periodic medium, these equations may be expressed as a system of infinitely many coupled mode equations, which we call the extended nonlinear coupled mode system (xNLCME). Truncating xNLCME to include only the principle resonances leads to the classical NLCME. Numerical simulations of xNLCME demonstrate that it captures both large scale features, related to third harmonic generation, and fine scale carrier shocks features of the nonlinear periodic Maxwell equations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Solitary Wave Benchmarks in Magma Dynamics

We present a model problem for benchmarking codes that investigate magma migration in the Earth's interior. This system retains the essential features of more sophisticated models, yet has the advantage of possessing solitary wave solutions. The existence of such exact solutions to the nonlinear problem make it an excellent benchmark problem for combinations of solver algorithms. In this work, we explore a novel algorithm for computing high quality approximations of the solitary waves and use them to benchmark a semi-Lagrangian Crank-Nicholson scheme for a finite element discretization of the time dependent problem.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Spectral Analysis for Matrix Hamiltonian Operators

In this work, we study the spectral properties of matrix Hamiltonians generated by linearizing the nonlinear Schrödinger equation about soliton solutions. By a numerically assisted proof, we show that there are no embedded eigenvalues for the three dimensional cubic equation. Though we focus on a proof of the 3d cubic problem, this work presents a new algorithm for verifying certain spectral properties needed to study soliton stability. Source code for verification of our comptuations, and for further experimentation, are available at http://www.math.toronto.edu/simpson/files/spec_prop_code.tgz.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Vortex Collapse for the L2-Critical Nonlinear Schrödinger Equation

The focusing cubic nonlinear Schrödinger equation in two dimensions admits vortex solitons, standing wave solutions with spatial structure, Qm(r,theta) = e^{i m theta} Rm(r). In the case of spin m = 1, we prove there exists a class of data that collapse with the vortex soliton profile at the log-log rate. This extends the work of Merle and Raphaël, (the case m = 0,) and suggests that the L2 mass that may be concentrated at a point during generic collapse may be unbounded. Difficulties with m >= 2 or when breaking the spin symmetry are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A Multiscale Model of Partial Melts 1: Effective Equations

In this paper a model for partial melts is constructed using two-scale homogenization theory. While this technique is well known to the mathematics and materials communities, it is relatively novel to problems in the solid Earth. This approach begins with a grain scale model of the medium, coarsening it into a macroscopic one. The emergent model is in good agreement with previous work, including D. McKenzie's, and serves as verification. This methodology also yields a series of Stokes problems whose solutions provide constitutive relations for permeability and viscosity. A numerical investigation of these relations appears in a companion paper.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A Multiscale Model of Partial Melts 2: Numerical Results

In a companion paper, equations for partially molten media were derived using two-scale homogenization theory. One advantage of homogenization is that material properties, such as permeability and viscosity, readily emerge. A caveat is that the dependence of these parameters upon the microstructure is not self-evident. In particular, one seeks to relate them to the porosity. In this paper, we numerically solve ensembles of the cell problems from which these quantities emerge. Using this data, we estimate relationships between the parameters and the porosity. In particular, the bulk viscosity appears to be inversely proportional to the porosity. Finally, we synthesize these numerical estimates with the models. Our hybrid numerical--analytical model predicts that the compaction length vanishes with porosity.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Arrest of Langmuir wave collapse by quantum effects

The arrest of Langmuir-wave collapse by quantum effects, first addressed by Haas and Shukla [Phys. Rev. E 79, 066402 (2009)] using a Rayleigh-Ritz trial-function method is revisited, using rigorous estimates and systematic asymptotic expansions. The absence of blow up for the so-called quantum Zakharov equations is proved in two and three dimensions, whatever the strength of the quantum effects. The time-periodic behavior of the solution for initial conditions slightly in excess of the singularity threshold for the classical problem is established for various settings in two space dimensions. The difficulty of developing a consistent perturbative approach in three dimensions is also discussed, and a semi-phenomenological model is suggested for this case.

0 citations0 reviewsNo code linkedOpen access