Quantitative Derivation of the Gross‐Pitaevskii Equation

Starting from first‐principle many‐body quantum dynamics, we show that the dynamics of Bose‐Einstein condensates can be approximated by the time‐dependent nonlinear Gross‐Pitaevskii equation, giving a bound on the rate of the convergence. Initial data are constructed on the bosonic Fock space applying an appropriate Bogoliubov transformation on a coherent state with expected number of particles N. The Bogoliubov transformation plays a crucial role; it produces the correct microscopic correlations among the particles. Our analysis shows that, on the level of the one‐particle reduced density, the form of the initial data is preserved by the many‐body evolution, up to a small error that vanishes as N^?1/2 in the limit of large N.© 2015 Wiley Periodicals, Inc.     

Computational fluid dynamics simulation of fluid particle fragmentation in turbulent flows

A simulation methodology is presented that allows detailed studies of the breakup mechanism of fluid particles in turbulent flows. The simulations, based on large eddy and volume of fluid simulations, agree very well with high-speed measurements of the breakup dynamics with respect to deformation time and length scales, and also the resulting size of the daughter fragments. The simulations reveal the size of the turbulent vortices that contribute to the breakup and how fast the interaction and energy transfer occurs. It is concluded that the axis of the deformed particle and the vortex core axis are aligned perpendicular to each other, and that breakup sometimes occurs due to interaction with two vortices at the same time. Analysis of the energy transfer from the continuous phase turbulence to the fluid particles reveals that the deformed particle attains it maximum in interfacial energy before the breakup is finalized. Similar to transition state theory in chemistry this implies that an activation barrier exists. Consequently, by considering the dynamics of the phenomenon, more energy than required at the final stage needs to be transferred from the turbulent vortices for breakup to occur. This knowledge helps developing new, more physical sound models for the breakup phenomenon required to solve scale separation problems in computational fluid dynamics simulations.     

Convergence to equilibrium of biased plane Partitions

We study a single‐flip dynamics for the monotone surface in (2 + 1) dimensions obtained from a boxed plane partition. The surface is analyzed as a system of non‐intersecting simple paths. When the flips have a non‐zero bias we prove that there is a positive spectral gap uniformly in the boundary conditions and in the size of the system. Under the same assumptions, for a system of size M, the mixing time is shown to be of order M up to logarithmic corrections. © 2010 Wiley Periodicals, Inc. Random Struct. Alg., 39, 83–114, 2011     

Sharp asymptotics for Kawasaki dynamics on a finite box with open boundary

In this paper we study the metastable behavior of the lattice gas in two and three dimensions subject to Kawasaki dynamics in the limit of low temperature and low density. We consider the local version of the model, where particles live on a finite box and are created, respectively, annihilated at the boundary of the box in a way that reflects an infinite gas reservoir. We are interested in how the system nucleates, i.e., how it reaches a full box when it starts from an empty box. Our approach combines geometric and potential theoretic arguments. In two dimensions, we identify the full geometry of the set of critical droplets for the nucleation, compute the average nucleation time up to a multiplicative factor that tends to one, show that the nucleation time divided by its average converges to an exponential random variable, express the proportionality constant for the average nucleation time in terms of certain capacities associated with simple random walk, and compute the asymptotic behavior of this constant as the system size tends to infinity. In three dimensions, we obtain similar results but with less control over the geometry and the constant. A special feature of Kawasaki dynamics is that in the metastable regime particles move along the border of a droplet more rapidly than they arrive from the boundary of the box. The geometry of the critical droplet and the sharp asymptotics for the average nucleation time are highly sensitive to this motion.     

Turbulent particle dispersion in an electrostatic precipitator

The behaviour of charged particles in turbulent gas flow in electrostatic precipitators (ESPs) is crucial information to optimise precipitator efficiency. This paper describes a strongly coupled calculation procedure for the rigorous computation of particle dynamics during ESP taking into account the statistical particle size distribution. The turbulent gas flow and the particle motion under electrostatic forces are calculated by using the commercial computational fluid dynamics (CFD) package FLUENT linked to a finite volume solver for the electric field and ion charge. Particle charge is determined from both local electrical conditions and the cell residence time which the particle has experienced through its path. Particle charge density and the particle velocity are averaged in a control volume to use Lagrangian information of the particle motion in calculating the gas and electric fields. The turbulent particulate transport and the effects of particulate space charge on the electrical current flow are investigated. The calculated results for poly-dispersed particles are compared with those for mono-dispersed particles, and significant differences are demonstrated.     

ANALYSIS OF SENSITIVITY AND UNCERTAINTY IN AN INDIVIDUAL‐BASED MODEL OF A THREATENED WILDLIFE SPECIES

Sensitivity analysis—determination of how prediction variables affect response variables—of individual‐based models (IBMs) are few but important to the interpretation of model output. We present sensitivity analysis of a spatially explicit IBM (HexSim) of a threatened species, the Northern Spotted Owl (NSO; Strix occidentalis caurina) in Washington, USA. We explored sensitivity to HexSim variables representing habitat quality, movement, dispersal, and model architecture; previous NSO studies have well established sensitivity of model output to vital rate variation. We developed “normative” (expected) model settings from field studies, and then varied the values of ≥ 1 input parameter at a time by ±10% and ±50% of their normative values to determine influence on response variables of population size and trend. We determined time to population equilibration and dynamics of populations above and below carrying capacity. Recovery time from small population size to carrying capacity greatly exceeded decay time from an overpopulated condition, suggesting lag time required to repopulate newly available habitat. Response variables were most sensitive to input parameters of habitat quality which are well‐studied for this species and controllable by management. HexSim thus seems useful for evaluating potential NSO population responses to landscape patterns for which good empirical information is available.     

Mean‐Field Evolution of Fermionic Mixed States

In this paper we study the dynamics of fermionic mixed states in the mean‐field regime. We consider initial states that are close to quasi‐free states and prove that, under suitable assumptions on the initial data and on the many‐body interaction, the quantum evolution of such initial data is well approximated by a suitable quasi‐free state. In particular, we prove that the evolution of the reduced one‐particle density matrix converges, as the number of particles goes to infinity, to the solution of the time‐dependent Hartree‐Fock equation. Our result holds for all times and gives effective estimates on the rate of convergence of the many‐body dynamics towards the Hartree‐Fock evolution.© 2015 Wiley Periodicals, Inc.     

The brazilian nut effect by void filling: an analytic model

We propose an analytic model of the Brazilian nut effect (BNE) that utilizes void filling as the primary mechanism behind the rise to the surface of one large buried intruder particle in an externally driven confined mixture with many small particles. When the intruder rises upward, it creates a void underneath it that is immediately filled by small particles preventing the intruder from sinking back to its previous position. Even though the external driving is only along the vertical direction, the small particles are able to move transversely into the void due to the Janssen effect, which postulates that the magnitudes of the vertically directed and transversely directed forces in a confined granular system are directly proportional to each other. The Janssen effect allows us to calculate the transverse speed distribution in which the small particles fill up the void and the temporal dynamics of the intruder particle in the vertically shaken container. We determine the time‐dependent behavior of the intruder vertical position h, its rise velocity dh/dt, and the phase (dh/dt vs. h) as a function of particle size ratio Φ, container diameter D_c, kinetic friction coefficient μ, and packing fraction c. Finally, we show that the predictions of our BNE model agree well with published experimental and simulation results. © 2010 Wiley Periodicals, Inc. Complexity 16: 9–16, 2011     

平面射流中纳米粒子积聚的矩方法

应用大涡模拟方法求解平面湍射流场,矩方法求解纳米粒子的一般动力学方程．通过对每种情况3 000个时间步的平均,得到了Schmidt数和Damkohler数对纳米粒子动力学特性的影响．结果发现, 当气体参数不变时,Schmidt数的变化只对直径小于1 nm的颗粒数密度的分布产生影响．直径小的颗粒其颗粒数密度沿流动方向下降迅速,而具有大Schmidt数的颗粒,沿横向的分布较窄．较小的颗粒容易发生积聚和扩散,并且体积增长较快,因而颗粒多分散性较为明显．小的颗粒积聚时间尺度能增强颗粒的碰撞和积聚频率,导致颗粒尺寸迅速增大．Damkohler数越大,颗粒的多分散也越明显．     

A comparison between the material point method and the finite element method in view of extrusion problems

An appropriate method for the simulation of continuous forming processes is the material point method (MPM) [1],[2] which combines the viewpoints of fluid dynamics and solid mechanics. The MPM and related methods [3] are derived from the particle‐in‐cell methods [4]. Bodies are discretised by Lagragian particles with pointwise mass distributions. The differential equations in their weak form are solved on temporary meshes built of standard finite elements. At the end of each time step the particle positions are updated and the mesh is replaced by a new mesh with a regular shape. The state variables at the nodes of the new mesh are extracted from the state variables at the particles by a transfer algorithm. When particles pass element boundaries, numerical difficulties might be observed. These are eliminated by a smooth approximation of nodal data from material point data. The modified MPM has been implemented together with the FEM in one programme because the similarities of the methods outbalance the differences. On the basis of numerical examples the results of both methods are compared. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Existence and uniqueness of conservative solutions for nonlocal fragmentation models

The aim of the paper is to investigate the well‐posedness of an integrodifferential equation describing multiple fragmentation processes, where the fragmentation rate is size and position dependent and new particles are randomly distributed in the space according to some probability density. An old method of Reuter and Lederman based on some approximation techniques is modified and applied to analyse the dynamics of the problem. The results of earlier papers on local models are extended to nonlocal models and the conservativeness of the solutions is investigated. In particular we prove the uniqueness of conservative solutions. Copyright © 2010 John Wiley & Sons, Ltd.     

INFERRING A BIOPOLITICAL CONSENSUS VIEW OF STOCHASTIC DYNAMICS FOR MANAGEMENT OF A TRANSBOUNDARY FISHERY

ABSTRACT. Management of trans‐boundary fisheries is a complicated problem with biological, legal, economic and political implications. We propose a simple stochastic differential‐equation model to describe a biopolitical consensus view of fish stock dynamics. Estimates of the drift and diffusion terms of three stochastic differential equations are obtained using data from the southern bluefin tuna (SBT) fishery with a method based on the Kolmogorov‐Smirnov statistic. We refer to these estimated equations as alternative biopolitical consensus views of SBT stock dynamics. Each of these is used to generate a time series of optimal harvest that achieves the objective of maximizing the present value of expected fishery returns. These time series of optimal harvests are then compared to actual harvests for the period 1981 1997.     

On interactive chains with finite populations

An Interactive Markov Chain is a population process in which each individuals's transitions depend on the population's distribution over the various states. We investigate a certain aspect of such process’ dynamics for a fixed population size. Conditions for convergence to steady‐state regardless of population size are provided.     

Weak Dual Pairs and Jetlet Methods for Ideal Incompressible Fluid Models in $$n \ge 2$$ Dimensions

We review the role of dual pairs in mechanics and use them to derive particle-like solutions to regularized incompressible fluid systems. In our case we have a dual pair resulting from the action of diffeomorphisms on point particles (essentially by moving the points). We then augment our dual pair by considering the action of diffeomorphisms on Taylor series, also known as jets. The augmented weak dual pairs induce a hierarchy of particle-like solutions and conservation laws with particles carrying a copy of a jet group. We call these augmented particles jetlets. The jet groups serve as finite-dimensional models of the diffeomorphism group itself, and so the jetlet particles serve as a finite-dimensional model of the self-similarity exhibited by ideal incompressible fluids. The conservation law associated to jetlet solutions is shown to be a shadow of Kelvin’s circulation theorem. Finally, we study the dynamics of infinite time particle mergers. We prove that two merging particles at the zeroth level in the hierarchy yield dynamics which asymptotically approach that of a single particle in the first level in the hierarchy. This merging behavior is then verified numerically as well as the exchange of angular momentum which must occur during a near collision of two particles. The resulting particle-like solutions suggest a new class of meshless methods which work in dimensions \(n \ge 2\) and which exhibit a shadow of Kelvin’s circulation theorem. More broadly, this provides one of the first finite-dimensional models of self-similarity in ideal fluids.     

Bogoliubov equations and functional mechanics

The functional classical mechanics based on the probability approach, where a particle is described not by a trajectory in the phase space but by a probability distribution, was recently proposed for solving the irreversibility problem, i.e., the problem of matching the time reversibility of microscopic dynamics equations and the irreversibility of macrosystem dynamics. In the framework of functional mechanics, we derive Bogoliubov-Boltzmann-type equations for finitely many particles. We show that a closed equation for a one-particle distribution function can be rigorously derived in functional mechanics without any additional assumptions required in the Bogoliubov method. We consider the possibility of using diffusion processes and the Fokker-Planck-Kolmogorov equation to describe isolated particles.     

A group theory derivation of the relativistic path integral and the “history-string” dynamics

We derive the Feynman path integral for relativistic elementary particles using group theory considerations. We apply an approach in which choosing a symmetry group (or semigroup) allows deriving the kinematics and dynamics of a particle including the state space and the propagator from it. The quantum properties of a particle appear from intertwining two representations of the symmetry (semi)group, one of which describes local properties of the particle and the other describes the particle as a whole. The path-integrallike dynamics appears when the symmetry semigroup has a structure similar to that of the relativistic analogue of the Galilei group (in which the Lorentz-invariant “proper time” plays the role of time) with translations replaced with the semigroup of trajectories (parameterized paths). The classical action in the weight functional of the path integral is defined by this semigroup up to couplings to gauge and/or gravitational fields. The obtained formalism is suitable for describing not only pointlike particles but also nonlocal objects of the “history-string” type, which, as previously shown, allow explaining quark confinement.     

An observation of a nascent fractal pattern in MD simulation for a fragmentation of an fcc lattice

To seek for a possible origin of fractal pattern in nature, we perform a molecular dynamics simulation for a fragmentation of an infinite fcc lattice. The fragmentation is induced by the initial condition of the model that the lattice particles have the Hubble-type radial expansion velocities. As time proceeds, the average density decreases and density fluctuation develops. By using the box counting method, it is found that the frequency-size plot of the density follows instantaneously a universal power-law for each Hubble constant up to the size of a cross-over. This cross-over size corresponds to the maximum size of fluctuation and is found to obey a dynamical scaling law as a function of time. This instantaneous generation of a nascent fractal is purely of dynamical origin and it shows us a new formation mechanism of a fractal patterns different from the traditional criticality concept.     

Complex dynamics of a discrete fractional‐order Leslie‐Gower predator‐prey model

A proposed discretized form of fractional‐order prey‐predator model is investigated. A sufficient condition for the solution of the discrete system to exist and to be unique is determined. Jury stability test is applied for studying stability of equilibrium points of the discretized system. Then, the effects of varying fractional order and other parameters of the systems on its dynamics are examined. The system undergoes Neimark‐Sacker and flip bifurcation under certain conditions. We observe that the model exhibits chaotic dynamics following stable states as the memory parameter α decreases and step size h increases. Theoretical results illustrate the rich dynamics and complexity of the model. Numerical simulation validates theoretical results and demonstrates the presence of rich dynamical behaviors include S‐asymptotically bounded periodic orbits, quasi‐periodicity, and chaos. The system exhibits a wide range of dynamical behaviors for fractional‐order α key parameter.     

Nonstationary probabilities and time correlation functions for an asymmetric exclusion process

We give the complete solution of the master equation for a system of interacting particles with finite density. We obtain the solution using a new form of the Bethe ansatz for an asymmetric simple exclusion process on the ring. We first find the one-point time correlation function for the discrete version of the process. __________ Translated from Teoreticheskaya i Matematicheskaya Fizika, Vol. 146, No. 3, pp. 499–508, March, 2006.     

On the global existence of weak solutions for the Cucker-Smale-Navier-Stokes system with shear thickening

We study the large-time dynamics of Cucker-Smale (C-S) flocking particles interacting with non-Newtonian incompressible fluids. Dynamics of particles and fluids were modeled using the kinetic Cucker-Smale equation for particles and non-Newtonian Navier-Stokes system for fluids, respectively and these two systems are coupled via the drag force, which is the main flocking (alignment) mechanism between particles and fluids. We present a global existence theory for weak solutions to the coupled Cucker-Smale-Navier-Stokes system with shear thickening. We also use a Lyapunov functional approach to show that sufficiently regular solutions approach flocking states exponentially fast in time.