1 / 77

A U niversal Field Equation for Dispersive Processes

A U niversal Field Equation for Dispersive Processes. J. H. Cushman, Purdue University Former Students M. Park , University of Alabama N. Kleinfelter-Domelle , Brown University M . Moroni , University of Rome, I B. Stroud , Private Industry M. Schoen , Berlin Technical.

shilah
Download Presentation

A U niversal Field Equation for Dispersive Processes

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. A Universal Field Equation for Dispersive Processes J. H. Cushman, Purdue University Former Students M. Park, University of Alabama N. Kleinfelter-Domelle, Brown University M. Moroni, University of Rome, I B. Stroud, Private Industry M. Schoen, Berlin Technical

  2. Anomalous Diffusion in Slit Nanopores Computational Statistical Mechanics: Quasi-static experiments

  3. Slit Pore Expansion Each frame of the animation shows a single configuration generated in a Monte Carlo simulation in the grand-canonicalisostrain ensemble (μ, T, h*, α) for a fixed pore width. Played in succession, the frames show configurations for the quasistatic pore expansion.  It is important to note that the animation does not portray the time evolution of an expanding pore. Periodic boundary conditions are imposed in the plane of the pore. The configurations shown represent small sub-regions of the fluid in an infinite slit pore. The pore walls are five atoms thick, but only a single layer of atoms is rendered for each wall with the remainder of the walls being indicated in outline. As the pore expands you should observe that the configurations are highly ordered for particular wall separations. This order is the result of freezing of the pore fluid.

  4. Simulation Cell Snapshot of a single configuration in the simulation cell.

  5. Thick Nano-Wire Expansion Maps of two dimensional slices in the plane y=0 through the three dimensional, ensemble averaged particle density for a Lennard-Jones 10-wire • reduced chemical potential, μ = -11.7 • reduced temperature, T= 1.0 • wall registration, alpha = 0.0 and reduced pore serpation, h*, as indicated. Points on the rendered density surfaces are displaced towards the viewer in proportion to the density.

  6. Thin Nano-Wire Expansion Maps of two dimensional slices in the plane y=0 through the three dimensional, ensemble averaged particle density for a Lennard-Jones 5-wire • reduced chemical potential, mu = -11.7 • reduced temperature, T= 1.0 • wall registration, alpha = 0.0 and reduced pore serpation, h*, as indicated. Points on the rendered density surfaces are displaced towards the viewer in proportion to the density.

  7. Nano-Wire and Slit Pore Expansion Maps of two dimensional slices through the three dimensional, ensemble averaged particle density in the plane y=0 for reduced chemical potential, μ = -11.7, reduced temperature, T= 1.0,wall registration, α = 0.0 and reduced pore serpation, h*, for (from left to right) a Lennard Jones 5-wire, 10-wire, and an infinite-wire (slit pore). Points on the rendered density surfaces are displaced towards the viewer in proportion to the density.

  8. Anomalous Diffusion When the mean-free-path of a molecule is on the same order as a characteristic dimension of the pore, the classical Fickian expression relating the diffusive flux to the activity gradient breaks down. That is, both diffusivity and the related viscosity become wave vector and frequency dependent.

  9. A Gedankin Experiment

  10. Panel b is especially interesting as it shows sub-diffusion. Here the mean-square displacement scales with d =0.64 as opposed to the Brownian limit of d =1.

  11. General Dispersion theory without scale constraints G is the conditional probability of finding a particle at x at time t given itwas initially at X(0). If all particles are identical and are released from the origin, then G represents the concentration. The Fourier transform of G is Ĝ

  12. Letting G´ in the last slide be g above, then upon differentiating and assuming dP/dt=0 where P is underlying probability and G'(k,0)=0 we obtain the following equation in Fourier space General Dispersion theory without scale constraints For suitably well behaved g(t) we have where Here we’ve introduced the Exponential Differential Displacement:

  13. General Dispersion theory without scale constraints With a little work can be rewritten

  14. General Dispersion theory without scale constraints Now, invert into real space and obtain where Di' is the inverse transform of (Di'Δ)^. This is an equation for dispersion without the need for scale constraints.

  15. At equilibrium the first and second terms are zero

  16. Small wave vector limit-retrieving the classical ADE Assume constant mean velocity so that depends only on k and τ making a convolution. Using a Taylor series expansion about |k|=0, on G^ it can be shown that If the system is stationary, then the mean square displacement can be replaced by twice velocity covariance.

  17. Anomalous Dispersion in Porous Media 3-D Porous Media PTV Experiments

  18. Porous media experimental set-up Particle Tracking Velocimetry (2D-PTV)

  19. Three-Dimensional Particle Tracking Velocimetry

  20. Porous Media Flow

  21. Baricenters from the coupled cameras (Ex. 1) Baricenters recognized on the image acquired by camera  1 (21741 baricenters) Baricenters recognized on the image acquired by camera  2 (16740 baricenters)

  22. 3D-PTV3-D Trajectory reconstruction from two 2-D projected trajectories

  23. 3 dimensional view of trajectories longer than 6 seconds

  24. General Dispersion theory without scale constraints This equation for dispersion is just as applicable in porous media as in the nano-film.

  25. Generalized dispersion coefficient Generalized dispersion coefficient in the transverse directions, k=2p/d, compared with the velocity covariance .

  26. CTRW Master Equation for Porous Media

  27. Anomalous Dispersion of Motile Bacteria

  28. Bacterial Swimming Videos courtesy of Howard Berg, Molecular and Cellular Biology, Harvard University

  29. Particle tracking of microbes

  30. Bacterial mortor and drive train • courtesy of David DeRosier, Brandeis university

  31. Microbe vs. Levy Motion E coli swimming Berg, Phys Today, Jan 2000. / Levy motion with α=1.2

  32. 25 20 15 10 5 0 -5 -10 -15 -20 0 10 20 30 40 LévyversusBrownianmotion a = 2 fractal dimension of the trace of Lévy motion (including Brownian motion) = a DISTANCE FROM ORIGIN a = 1.7 DISTANCE FROM ORIGIN

  33. Levy trajectory dispersion A particle which follows an α-stable Levy motion has a transition density given by the following equation (for a divergence free velocity field) where the general (asymmetric) fractional derivative is defined in Fourier space by the following with M the mixing measure, S the d-dimensional unit sphere

  34. Rewrite the fractional derivative as: Now take the Inverse Fourier transform : Inverse FT Levy trajectory dispersion

  35. Levy trajectory dispersion The equation can be written in the divergence form with convolution-Fickian flux

  36. The Finite-Size Lyapunov Exponent:Determining the exponent for Levy motion Definethe Finite-Size Lyapunov Exponent (FSLE) The FSLE describes the exponential divergence of two trajectories that start a specified distance, r, apart and depends on that initial starting scale as well as the threshold ratio, a. or la: is the expected exponential rate two particles separated by a distance r initially, separate to ar at time t. Td is the expected doubling time, i.e. the time it takes two particles to separate from r to ar.

  37. A Statistical Mechanical Formulation of the Finite-Size Lyapunov Exponent Assume for simplicity that the system is completely ‘expansive’. Let <<>> indicate integration over the product space with respect to a joint probability density ,f, defined on (τ) x  (τ+t). Here is the joint probability particles i and j have separation in at time τ and separation in at time

  38. A Statistical Mechanical Formulation of the Finite-Size Lyapunov Exponent where G(x,y,t) is the probability of a particle going from separation x to separation y in time t. Therefore the a-time is:

  39. FSLE for α-stable Lévy Motion Varying α Slope-1.2 Slope -1.5 Slope -2

  40. For the case when (1+α)/2 is an integer (i.e., α=1), the summation P(2,1)=2/16

  41. Turbulent Flows • Fluid Flow: Eulerian approach Lagrangian approach (1) Laminar Flow: Re = small (2) Turbulent Flow: Re = very large

  42. Turbulent mixing layer growth and internal waves formation: laboratory simulations

  43. The flux through the interface between the mixing layer and the stable layer plays a fundamental role in characterizing and forecasting the quality of water in stratified lakes and in the upper oceans, and the quality of air in the atmosphere.

  44. Experimental set-up

  45. Convective Layer Growth and Internal waves

  46. Only internal waves

  47. Vortex Generating Chromatograph

More Related