Diffusion of tagged particles in a crowded medium: supplementary material. Paper submitted to Europhys. Lett.
|
|
- Reynold Sparks
- 5 years ago
- Views:
Transcription
1 Diffusion of tagged particles in a crowded medium: supplementar material. Paper submitted to Europhs. Lett. Marta Galanti 1,2,3, Duccio Fanelli 2, Amos Maritan 4, Francesco Piazza 3 1 Dipartimento di Sistemi e Informatica,Università di Firenze and IF, Via S. Marta 3, IT Florence, Ital 2 Dipartimento di Fisica e Astronomia, Università di Firenze and IF, Via Sansone 1, IT Sesto Fiorentino, Ital 3 Université d Orléans and Centre de Biophsique Moléculaire CBM, CRS-UPR 4301, Rue Charles Sadron, Orléans, France 4 Dipartimento di Fisica e Astronomia G. Galilei Università di Padova, CISM and IF, via Marzolo 8, 35131, Padova, Ital PACS numbers: I. DERIVATIO OF THE MEA FIELD EQUATIOS I TWO DIMESIO As stated in the main tet, the equations governing the evolution of the particles, which we obtained for a 1D setting, can be etended to higher dimensions. In the following we will provide a detailed derivation for the mean field equations in 2D. We will first generalized the procedure described in the main bod of the paper and then turn to consider and alternative approach inspired to the work of Landman and collaborators see e.g. Simpson, Landman and Hughes Phs Rev E To progress in the analsis we assume each site of the the two-dimensional lattice to be labelled with two indices i, j. In the 2D case, the nearest neighbour sites the selected particle can jump to are four. This is at variance with the 1D geometr where each site has just two adjacent neighbors. The binar variables at time k are hence indicated as m i,j k and n i,j k. The stochastic process reads: m i,j k + 1 m i,j k = z + i 1 m i 1,jk[1 m i,j k][1 n i,j k] + z i+1 m i+1,jk[1 m i,j k][1 n i,j k] +z + j 1 m i,j 1k[1 m i,j k][1 n i,j k] + z j+1 m i,j+1k[1 m i,j k][1 n i,j k] z + i m i,jk[1 m i+1,j k][1 n i+1,j k] z i m i,jk[1 m i 1,j k][1 n i 1,j k] z + i m i,jk[1 m i,j+1 k][1 n i,j+1 k] z i m i,jk[1 m i,j 1 k][1 n i,j 1 k] 1 where the stochastic variables z ± are defined in analog to the one dimensional case. The equation governing the evolution of n i,j can be equivalentl modified. After introducing the one-bod occupanc probabilities and assuming a mean-field factorization for the two-bod and three-bod correlations, one eventuall ends up with ρ i,j k = m i,j k 2 φ i,j k = n i,j k 3 ρ i,j k + 1 ρ i,j k = q ρ i 1,j k + ρ i+1,j k + ρ i,j 1 k + ρ i,j+1 k [1 ρ i,j k] [1 φ i,j k] q ρ i,j k[4 ρ i 1,j k + ρ i+1,j k ρ i,j 1 k + ρ i,j+1 k φ i 1,j k + φ i+1,j k 4 φ i,j 1 k + φ i,j+1 k φ i+1,j kρ i+1,j k + φ i 1,j kρ i 1,j k + φ i,j+1 kρ i,j+1 k + φ i,j 1 kρ i,j 1 k] φ i,j k + 1 φ i,j k = w φ i 1,j k + φ i+1,j k + φ i,j 1 k + φ i,j+1 k [1 φ i,j k] [1 ρ i,j k] w φ i,j k[4 φ i 1,j k + φ i+1,j k φ i,j 1 k + φ i,j+1 k ρ i 1,j k + ρ i+1,j k ρ i,j 1 k + ρ i,j+1 k + φ i+1,j kρ i+1,j k + φ i 1,j kρ i 1,j k + φ i,j+1 kρ i,j+1 k + φ i,j 1 kρ i,j 1 k] 5
2 2 Assume the concentration of the tagged particles to be small, namel ρ i,j 1. Then the following approimated relations are found: ρ i,j k + 1 ρ i,j k = qρ i 1,j k + ρ i+1,j k + ρ i,j 1 k + ρ i,j+1 k [1 φ i,j k] qρ i,j k[4 φ i 1,j k + φ i+1,j k + φ i,j 1 k + φ i,j+1 k] φ i,j k + 1 φ i,j k = wφ i 1,j k + φ i+1,j k 2φ i,j k + φ i,j 1 k + φ i,j+1 k 2φ i,j k After introducing the continuous variables ρ,, t = lim ρ i,jk, φ,, t = lim φ i,jk a, t 0 a, t 0 6 where a and t respectivel stand for the linear size of the lattice site and the characteristic time step of the microscopic dnamics. B defining the diffusion coefficients according to the standard practice see main bod of the paper, we get the sought generalized model: ρ = 2 [D ρ 1 φρ] + 2D ρ ρ φ φ = D φ 2 φ 7 As anticipated, we shall now turn to discussing an alternative derivation of the above equation. To this end, we define the variable 1, if the site i, j is occupied b an agent of the species m at time step k γ i,j k = 2, if the site i, j is occupied b an agent of the species n at time step k 0 if the site i, j is empt at the time step k 8 We write then the master equation for the motion of the tagged particles as: P 1 γ i,j k + 1 = 1 P 1 γ i,j k = 1 = α [P 2 γ i,j k = 1, γ i+1,j k = 0 + P 2 γ i,j k = 1, γ i 1,j k = 0 ] + P γ 2 i,j k = 1, γ i,j+1 k = 0 + P 2 γ i,j k = 1, γ i,j 1 k = 0 + α [P 2 γ i,j k = 0, γ i+1,j k = 1 ] + P γ 2 i,j k = 0, γ i 1,j k = 1 + P 2 γ i,j k = 0, γ i,j+1 k = 1 + P 2 γ i,j k = 0, γ i,j 1 k = 1 9 where P 2 stands for a joint probabilit, while P 1 is the probabilit of a single event and α is the rate of success of the selected jump. The master equation for the population n is similar, the variable γ assuming the value 2 instead of 1 and β labeling the associated jump rate in principle the two population can have different jump probabilities. In the mean field limit, we factorize the joint probabilities P 2 in the master equations as in P 2 γ i,j k = 1, γ i,j+1 k = 0 = P 1 γ i,j k = 1 P 1 γ i,j+1 k = 0. To perform the continuum limit we emplo a Talor epansion for the three different probabilit functions P 1 as:
3 3 Pγ i±1,j k = 1 = ρ ± a ρ a2 2 ρ 2 + oa2 Pγ i,j±1 k = 1 = ρ ± a ρ a2 2 ρ 2 + oa2 Pγ i±1,j k = 2 = φ ± a φ a2 2 φ 2 + oa2 Pγ i,j±1 k = 2 = φ ± a φ a2 2 φ 2 + oa2 Pγ i±1,j k = 0 = µ ± a µ a2 2 µ 2 + oa2 Pγ i,j±1 k = 0 = µ ± a µ a2 2 µ 2 + oa2. 10 where again a represents the linear size of each site. Making use of the relation µ = 1 ρ φ and defining as usual we obtain αa 2 lim a, t 0 t = D ρ βa 2 lim a, t 0 t = D φ, ρ = D ρ1 φ 2 ρ + ρ 2 φ for the tagged species. Under the hpotesis of low concentration of the tagged agents, we finall get: φ = D φ 2 φ for the evolution of the bulk densit. II. O THE VALIDITY OF THE MEA FIELD APPROXIMATIO: COMPARIG STOCHASTIC AD MEA FIELD SIMULATIOS I 2D To test the adequac of the proposed mean field model we have performed a campaign of numerical simulations, with reference to the 2D setting. More specificall we have implemented a montecarlo scheme to solve the stochastic process under scrutin and so trace the evolution of the tagged particle in time. At each time iteration, all crowders and the tagged particle can update their position, moving at random, and with equal probabilit, in one of the four allowed directions, provided the selected target site is unoccupied. If the site at destination is occupied, the move is rejected and the particles keep their original positions. The order of selection of the particles is, at each iteration, randomized. B averaging over man independent realizations, one can reconstruct the normalized histogram of the position visited b the crowders at a given time t and compare it with the densit profile ρ obtained upon integration of the mean field sstem 7. To carr out the numerical integration of the above partial differential equations we assumed a forward difference approimation in time and replaced the spatial derivatives b centered approimations. The result of the comparison is reported in figures 1a and 1b, for two choices of the initial conditions. In figure 1a, the tagged particle is initiall positioned in the middle of a two dimensional waterbag, filled with crowders, with average densit equal to φ 0. In figure 1a, the tagged particle is instead positioned, at time t = 0, in the center of an empt region, a square of assigned size. The crowders are instead assumed to occup an adjacent domain with average uniform densit φ 0. In both cases, the agreement between stochastic and mean field simulations is satisfing. III. ALTERATIVE DERIVATIO OF THE MODEL USIG A COARSE GRAIED PICTURE In this final section we discuss an alternative derivation of model 10, which assumes a coarse-grained decription of the scrutinized problem. The derivation follows a different philosph: it is here carried out in one dimension, but readil generalizes to the relevant d = 3 setting. We consider the phsical space to be partitioned in Ω patches, also called urns. Each patch has a maimum carring capacit - it can be filled with particles at most. Labelling m i the number of tagged particles contained in urn i, and with n i the corresponding number of crowders, one can write: n i + m i + v i = where v i stands for the number of vacancies, the empt cases in patch i that can be eventuall filled b incoming particles. The ecluded-volume prescription is here implemented b requiring that particles can move onl into the nearest-neighbor patches that ehibit vacancies, i
4 ρ 0,2 tagged particles, of tpe the crowders or a vacanc belonging to the i-patch. This is stochastic process governed, under the Markov hpothesis, b a Master equation for the probabilit P n, m, t of finding the sstem in a given state specified b the 2Ω dimensional vector n, m = n 1,..., n Ω, m 1..., m Ω at time t. The Master equation reads: , P n,m,t = n n [T n, m n, mp n, m + 12 T n, m n, m P n, m T n, m n, mp n, m T n, m n, mp n, m] where T a b is the rate of transition from a state a to a compatible configuration b. The allowed transitions are those that take place between neighboring patches as dictated b the chemical reactions 11. For eample, the transition probabilit associated with the second of equations 11 reads ρ FIG. 1: Comparison between stochastic thin solid lines and mean field simulations thick solid line. The stocastic simulations are averaged over realizations. Here φ 0 = 0.5. Two dimensional snapshots of the stochastic simulations are also displaed. The crowders are plotted as small black online circles. The tagged particle is represented b the filled red online square. Upper panel: initial condition see first snapshot, top left panel originating a super-diffusive transient. Lower panel initial condition see first snapshot, top left panel, causing a sub-diffusive transient. as eemplified b the following chemical reactions zω M i + V j M j + V i 11 zω i + V j j + V i Here z is the number of nearest-neighbor patches and M i, i, V i are respectivel a particle of tpe M the T n i 1, n j +1 n i, n j = n i v j zω = n i zω 1 n j m j. 13 The transition rates bring into the equation an eplicit dependence on the amount of molecules per patch, the so-called sstem size. To proceed in the analsis, we make use of van Kampen sstem size epansion [? ], which enables one to separate the site-dependent mean concentration φ i t from the corresponding fluctuations ξ i in the epression of the discrete number densit of species. The fluctuations become less influent as the number of the agents is increased, an observation which translates in the following van Kampen ansatz: n i t = φ it + ξ i. 14 In the following we will also assume just one tagged particle, the analsis etending straightforwardl to the case where a bunch of diluted particles is assumed to be dispersed in the background of crowders. Since the tagged particle belongs to one of the patches, it is convenient to look at the evolution of P k n, t = P n, 0, 0,..., 0, 1, 0,..., 0, 0, t }{{}}{{} k 1 Ω k in the master equation 12. P k n, t is the probabilit that the target particle be in the k-patch, for a particular configuration n of species. The Master equation can be hence written in the following compact form
5 5 P k n, t = + j i 1,i+1 j i 1,i+1 ɛ j ɛ+ i 1T n i 1, n j + 1 n i, n j P k n i, n j, t 1 zω 1 n j P k + j i 1,i+1 1 zω 1 n k P j 15 where use has been made of the shift operators: equation 13 takes the form ɛ ± i f..., n i,... = f..., n i ± 1,... Under the van Kampen prescription [? ], one can epand the transition rates in power of 1/. For eample, T n i 1, n j + 1 n i, n j = { φ i 1 φ j + 1 [ ξi 1 φ j ξ j φ i ] + 1 zω [ ξ iξ j m j φ i ] + 1 [ ] } mj ξ 3 i 2 and also epress the shift operators in terms of differential operators: ɛ j ɛ+ i 1 = O 1 ξ j 2 ξ j [ ɛ j ɛ+ i 1T n i 1, n j + 1 n i, n j P n i, n j, t = 1 ξ j zω ξ i1 φ j + ξ j φ i Π k ξ, t [ ξ j ] zω φ i1 φ j Π k ξ, t 2 ] ξ j zω φ i1 φ j Π k ξ, t + O otice that m i / cannot be approimated as a continuum-like densit, the continuum limit being not appropriate for the case of a single tracer. We then define a new probabilit distribution Π k ξ, τ, function of the vector ξ and the scaled time τ = t Ω. In terms of the new probabilit distribution Π k ξ, τ the left hand side of 15 becomes The leading order contribution in 1 gives: 1 Ω Π k φ i = which ields Π k φ i = z and finall: zω j {i 1,i+1} φ i 1 φ j Π k Π k ξ j 16 Π k 2φ i φ i 1 φ i+1 17 P k = 1 Ω Ω Π k φi + 1 Π k Ω. φ i = 2 φ i 18 where is the discrete Laplacian operator defined as φ i = 2 z j i φ j φ i, where j i means a summation
6 6 over the sites, j, which are nearest-neighbors of site i. B taking the size of the patches to zero, one recovers the standard diffusion equation for species φ, in agreement with the result reported in the main bod of the paper. Consider now the following identities: z j {i 1,i+1} = z = z z = 2 ξ j ξ i 1 φ j ξ j φ i Π k ξ i 1 φ i 1 ξ i 1 φ i Π k + 1 z ξ i 1 φ i+1 ξ i+1 φ i Π k +1 ξ i 1 φ i 1 ξ i 1 φ i Π k + ξ i 1 φ i+1 ξ i+1 φ i Π k ξ i+1 1 φ i ξ i φ i+1 Π k z ξ i Π k ξ i 1 1 φ i ξ i φ i 1 Π k and z j {i 1,i+1} ξ j 2 φ i 1 φ j Π k = 2 φ z ξ 2 i 1 φ i 1 Π k φ i ξi 2 i 1 φ i+1 Π k φ ξi+1 2 i 1 φ i+1 Π k + φ 2 ξi 1 2 i 1 φ i 1 Π k φ 22 i 1 φ i 1 Π k φ 22 i 1 φ i+1 Π k 1 +1 = 2 2 z ξ 2 2φ i + φ i 1 + φ i+1 2φ i φ i+1 + φ i 1 Π k + 2φ i 1 φ i 1 Π k i φ i 1 φ i+1 Π k +1 Making use of the above relations, at the net to netto-leading corrections one eventuall gets Π k = ξ i Π k + i+1 B i,j Π k 2 ξ i 2z ξ i=i 1 i ξ j 1 φ k Π k 1 2 φ k+1 φ k 1 Π k + 1 φ k Π k+1 + z Here B represents the diffusion matri, whose entries are B i,i = 2φ i + φ i 1 + φ i+1 2φ i φ i+1 + φ i 1 B i,i 1 = 2φ i 1 φ i 1 B i,i+1 = 2φ i 1 φ i To provide a mean-field description of the inspected problem, we consider the probabilit function of the
7 7 tagged agent integrated over the fluctuations of the - particles. In formulae: ρ k t = Π k dξ whose evolution is governed b ρ k = z 1 φk ρ k 1 2 φ k+1 φ k 1 ρ k + 1 φ k ρ k+1 = 2 ρ k φ k ρ k + ρ k φ k The the last epression involves the discrete laplacian defined above. In the continuum limit, and considering a straightforward generalization to higher dimensions, one gets ρr, t = D ρ 1 φr, t 2 ρr, t + D ρ ρr, t 2 φr, t where D ρ is the diffusion coefficient of the tagged particle. One can finall write the non-linear equation for ρ as a Fokker Plank equation: ρr, t = D1 2 φr, tρr, t +2D ρr, t φr, t Hence, b neglecting the role of fluctuations, which amounts to operating in the mean-field limit, a nonlinear partial differential equation is found for the densit of the tagged species, coupled to a standard diffusion equation for the background densit: φr,t = D 2 φr, t ρr,t = 2 D1 φr, tρr, t + 2D ρr, t φr, t 20 This sstem constitutes the generalization of model 10 to higher dimensions. It is worth emphasising that the second of eqs. 20 can be also cast in the alternative form: ρr, t = D 1 φr, t ρr, t + Dρr, t φr, t
Stochastic models, patterns formation and diffusion
Stochastic models, patterns formation and diffusion Duccio Fanelli Francesca Di Patti, Tommaso Biancalani Dipartimento di Energetica, Università degli Studi di Firenze CSDC Centro Interdipartimentale per
More informationInterspecific Segregation and Phase Transition in a Lattice Ecosystem with Intraspecific Competition
Interspecific Segregation and Phase Transition in a Lattice Ecosstem with Intraspecific Competition K. Tainaka a, M. Kushida a, Y. Ito a and J. Yoshimura a,b,c a Department of Sstems Engineering, Shizuoka
More information1.1 The Equations of Motion
1.1 The Equations of Motion In Book I, balance of forces and moments acting on an component was enforced in order to ensure that the component was in equilibrium. Here, allowance is made for stresses which
More informationMMJ1153 COMPUTATIONAL METHOD IN SOLID MECHANICS PRELIMINARIES TO FEM
B Course Content: A INTRODUCTION AND OVERVIEW Numerical method and Computer-Aided Engineering; Phsical problems; Mathematical models; Finite element method;. B Elements and nodes, natural coordinates,
More informationCONTINUOUS SPATIAL DATA ANALYSIS
CONTINUOUS SPATIAL DATA ANALSIS 1. Overview of Spatial Stochastic Processes The ke difference between continuous spatial data and point patterns is that there is now assumed to be a meaningful value, s
More informationParameterized Joint Densities with Gaussian and Gaussian Mixture Marginals
Parameterized Joint Densities with Gaussian and Gaussian Miture Marginals Feli Sawo, Dietrich Brunn, and Uwe D. Hanebeck Intelligent Sensor-Actuator-Sstems Laborator Institute of Computer Science and Engineering
More informationOn Maximizing the Second Smallest Eigenvalue of a State-dependent Graph Laplacian
25 American Control Conference June 8-, 25. Portland, OR, USA WeA3.6 On Maimizing the Second Smallest Eigenvalue of a State-dependent Graph Laplacian Yoonsoo Kim and Mehran Mesbahi Abstract We consider
More information16.5. Maclaurin and Taylor Series. Introduction. Prerequisites. Learning Outcomes
Maclaurin and Talor Series 6.5 Introduction In this Section we eamine how functions ma be epressed in terms of power series. This is an etremel useful wa of epressing a function since (as we shall see)
More informationInitial Value Problems for. Ordinary Differential Equations
Initial Value Problems for Ordinar Differential Equations INTRODUCTION Equations which are composed of an unnown function and its derivatives are called differential equations. It becomes an initial value
More informationKINEMATIC RELATIONS IN DEFORMATION OF SOLIDS
Chapter 8 KINEMATIC RELATIONS IN DEFORMATION OF SOLIDS Figure 8.1: 195 196 CHAPTER 8. KINEMATIC RELATIONS IN DEFORMATION OF SOLIDS 8.1 Motivation In Chapter 3, the conservation of linear momentum for a
More informationStrain Transformation and Rosette Gage Theory
Strain Transformation and Rosette Gage Theor It is often desired to measure the full state of strain on the surface of a part, that is to measure not onl the two etensional strains, and, but also the shear
More informationRoger Johnson Structure and Dynamics: The 230 space groups Lecture 3
Roger Johnson Structure and Dnamics: The 23 space groups Lecture 3 3.1. Summar In the first two lectures we considered the structure and dnamics of single molecules. In this lecture we turn our attention
More informationUncertainty and Parameter Space Analysis in Visualization -
Uncertaint and Parameter Space Analsis in Visualiation - Session 4: Structural Uncertaint Analing the effect of uncertaint on the appearance of structures in scalar fields Rüdiger Westermann and Tobias
More informationLecture 4.2 Finite Difference Approximation
Lecture 4. Finite Difference Approimation 1 Discretization As stated in Lecture 1.0, there are three steps in numerically solving the differential equations. They are: 1. Discretization of the domain by
More informationA Function of Two Random Variables
akultät Inormatik Institut ür Sstemarchitektur Proessur Rechnernete A unction o Two Random Variables Waltenegus Dargie Slides are based on the book: A. Papoulis and S.U. Pillai "Probabilit random variables
More information5. Nonholonomic constraint Mechanics of Manipulation
5. Nonholonomic constraint Mechanics of Manipulation Matt Mason matt.mason@cs.cmu.edu http://www.cs.cmu.edu/~mason Carnegie Mellon Lecture 5. Mechanics of Manipulation p.1 Lecture 5. Nonholonomic constraint.
More informationINF Introduction to classifiction Anne Solberg Based on Chapter 2 ( ) in Duda and Hart: Pattern Classification
INF 4300 151014 Introduction to classifiction Anne Solberg anne@ifiuiono Based on Chapter 1-6 in Duda and Hart: Pattern Classification 151014 INF 4300 1 Introduction to classification One of the most challenging
More informationSection 1.5 Formal definitions of limits
Section.5 Formal definitions of limits (3/908) Overview: The definitions of the various tpes of limits in previous sections involve phrases such as arbitraril close, sufficientl close, arbitraril large,
More informationStability Analysis of a Geometrically Imperfect Structure using a Random Field Model
Stabilit Analsis of a Geometricall Imperfect Structure using a Random Field Model JAN VALEŠ, ZDENĚK KALA Department of Structural Mechanics Brno Universit of Technolog, Facult of Civil Engineering Veveří
More informationChapter 3. Theory of measurement
Chapter. Introduction An energetic He + -ion beam is incident on thermal sodium atoms. Figure. shows the configuration in which the interaction one is determined b the crossing of the laser-, sodium- and
More informationAnalytic Geometry in Three Dimensions
Analtic Geometr in Three Dimensions. The Three-Dimensional Coordinate Sstem. Vectors in Space. The Cross Product of Two Vectors. Lines and Planes in Space The three-dimensional coordinate sstem is used
More informationMATHEMATICS 200 December 2011 Final Exam Solutions
MATHEMATICS December 11 Final Eam Solutions 1. Consider the function f(, ) e +4. (a) Draw a contour map of f, showing all tpes of level curves that occur. (b) Find the equation of the tangent plane to
More informationarxiv: v1 [cond-mat.stat-mech] 17 Nov 2017
Hopping in the crowd to unveil networ topology Malbor Asllani 1, Timoteo Carletti 1, Francesca Di Patti, Duccio Fanelli and Francesco Piazza 3 1 naxys, Namur Institute for Complex Systems, University of
More information1 Differential Equations for Solid Mechanics
1 Differential Eqations for Solid Mechanics Simple problems involving homogeneos stress states have been considered so far, wherein the stress is the same throghot the component nder std. An eception to
More informationConservation of Linear Momentum for a Differential Control Volume
Conservation of Linear Momentum for a Differential Control Volume When we applied the rate-form of the conservation of mass equation to a differential control volume (open sstem in Cartesian coordinates,
More informationy R T However, the calculations are easier, when carried out using the polar set of co-ordinates ϕ,r. The relations between the co-ordinates are:
Curved beams. Introduction Curved beams also called arches were invented about ears ago. he purpose was to form such a structure that would transfer loads, mainl the dead weight, to the ground b the elements
More informationx y plane is the plane in which the stresses act, yy xy xy Figure 3.5.1: non-zero stress components acting in the x y plane
3.5 Plane Stress This section is concerned with a special two-dimensional state of stress called plane stress. It is important for two reasons: () it arises in real components (particularl in thin components
More informationIntegration. 5.1 Antiderivatives and Indefinite Integration. Suppose that f(x) = 5x 4. Can we find a function F (x) whose derivative is f(x)?
5 Integration 5. Antiderivatives and Indefinite Integration Suppose that f() = 5 4. Can we find a function F () whose derivative is f()? Definition. A function F is an antiderivative of f on an interval
More information6.869 Advances in Computer Vision. Prof. Bill Freeman March 1, 2005
6.869 Advances in Computer Vision Prof. Bill Freeman March 1 2005 1 2 Local Features Matching points across images important for: object identification instance recognition object class recognition pose
More informationCSE 546 Midterm Exam, Fall 2014
CSE 546 Midterm Eam, Fall 2014 1. Personal info: Name: UW NetID: Student ID: 2. There should be 14 numbered pages in this eam (including this cover sheet). 3. You can use an material ou brought: an book,
More informationAnd similarly in the other directions, so the overall result is expressed compactly as,
SQEP Tutorial Session 5: T7S0 Relates to Knowledge & Skills.5,.8 Last Update: //3 Force on an element of area; Definition of principal stresses and strains; Definition of Tresca and Mises equivalent stresses;
More informationPolarization and Related Antenna Parameters
ANTENTOP- 01-007, # 009 Polarization and Related Antenna Parameters Feel Yourself a Student! Dear friends, I would like to give to ou an interesting and reliable antenna theor. Hours searching in the web
More informationAdsorption from a one-dimensional lattice gas and the Brunauer Emmett Teller equation
Proc. Natl. Acad. Sci. USA Vol. 93, pp. 1438 1433, December 1996 Chemistry Adsorption from a one-dimensional lattice gas and the Brunauer Emmett Teller equation (Ising modelreference systemgibbs ecess
More information6. Vector Random Variables
6. Vector Random Variables In the previous chapter we presented methods for dealing with two random variables. In this chapter we etend these methods to the case of n random variables in the following
More information520 Chapter 9. Nonlinear Differential Equations and Stability. dt =
5 Chapter 9. Nonlinear Differential Equations and Stabilit dt L dθ. g cos θ cos α Wh was the negative square root chosen in the last equation? (b) If T is the natural period of oscillation, derive the
More informationComputation fundamentals of discrete GMRF representations of continuous domain spatial models
Computation fundamentals of discrete GMRF representations of continuous domain spatial models Finn Lindgren September 23 2015 v0.2.2 Abstract The fundamental formulas and algorithms for Bayesian spatial
More informationBASE VECTORS FOR SOLVING OF PARTIAL DIFFERENTIAL EQUATIONS
BASE VECTORS FOR SOLVING OF PARTIAL DIFFERENTIAL EQUATIONS J. Roubal, V. Havlena Department of Control Engineering, Facult of Electrical Engineering, Czech Technical Universit in Prague Abstract The distributed
More informationChapter 6. Nonlinear Equations. 6.1 The Problem of Nonlinear Root-finding. 6.2 Rate of Convergence
Chapter 6 Nonlinear Equations 6. The Problem of Nonlinear Root-finding In this module we consider the problem of using numerical techniques to find the roots of nonlinear equations, f () =. Initially we
More informationResearch Article Development of a Particle Interaction Kernel Function in MPS Method for Simulating Incompressible Free Surface Flow
Journal of Applied Mathematics Volume 2, Article ID 793653, 6 pages doi:.55/2/793653 Research Article Development of a Particle Interaction Kernel Function in MPS Method for Simulating Incompressible Free
More informationCopyright, 2008, R.E. Kass, E.N. Brown, and U. Eden REPRODUCTION OR CIRCULATION REQUIRES PERMISSION OF THE AUTHORS
Copright, 8, RE Kass, EN Brown, and U Eden REPRODUCTION OR CIRCULATION REQUIRES PERMISSION OF THE AUTHORS Chapter 6 Random Vectors and Multivariate Distributions 6 Random Vectors In Section?? we etended
More informationSummary of Random Variable Concepts March 17, 2000
Summar of Random Variable Concepts March 17, 2000 This is a list of important concepts we have covered, rather than a review that devives or eplains them. Tpes of random variables discrete A random variable
More informationFACULTY OF MATHEMATICAL STUDIES MATHEMATICS FOR PART I ENGINEERING. Lectures AB = BA = I,
FACULTY OF MATHEMATICAL STUDIES MATHEMATICS FOR PART I ENGINEERING Lectures MODULE 7 MATRICES II Inverse of a matri Sstems of linear equations Solution of sets of linear equations elimination methods 4
More information5.6. Differential equations
5.6. Differential equations The relationship between cause and effect in phsical phenomena can often be formulated using differential equations which describe how a phsical measure () and its derivative
More informationOn the Extension of Goal-Oriented Error Estimation and Hierarchical Modeling to Discrete Lattice Models
On the Etension of Goal-Oriented Error Estimation and Hierarchical Modeling to Discrete Lattice Models J. T. Oden, S. Prudhomme, and P. Bauman Institute for Computational Engineering and Sciences The Universit
More informationLATERAL BUCKLING ANALYSIS OF ANGLED FRAMES WITH THIN-WALLED I-BEAMS
Journal of arine Science and J.-D. Technolog, Yau: ateral Vol. Buckling 17, No. Analsis 1, pp. 9-33 of Angled (009) Frames with Thin-Walled I-Beams 9 ATERA BUCKING ANAYSIS OF ANGED FRAES WITH THIN-WAED
More informationBACKWARD FOKKER-PLANCK EQUATION FOR DETERMINATION OF MODEL PREDICTABILITY WITH UNCERTAIN INITIAL ERRORS
BACKWARD FOKKER-PLANCK EQUATION FOR DETERMINATION OF MODEL PREDICTABILITY WITH UNCERTAIN INITIAL ERRORS. INTRODUCTION It is widel recognized that uncertaint in atmospheric and oceanic models can be traced
More informationMAT389 Fall 2016, Problem Set 2
MAT389 Fall 2016, Problem Set 2 Circles in the Riemann sphere Recall that the Riemann sphere is defined as the set Let P be the plane defined b Σ = { (a, b, c) R 3 a 2 + b 2 + c 2 = 1 } P = { (a, b, c)
More informationApplications of Gauss-Radau and Gauss-Lobatto Numerical Integrations Over a Four Node Quadrilateral Finite Element
Avaiable online at www.banglaol.info angladesh J. Sci. Ind. Res. (), 77-86, 008 ANGLADESH JOURNAL OF SCIENTIFIC AND INDUSTRIAL RESEARCH CSIR E-mail: bsir07gmail.com Abstract Applications of Gauss-Radau
More informationCONSERVATION LAWS AND CONSERVED QUANTITIES FOR LAMINAR RADIAL JETS WITH SWIRL
Mathematical and Computational Applications,Vol. 15, No. 4, pp. 742-761, 21. c Association for Scientific Research CONSERVATION LAWS AND CONSERVED QUANTITIES FOR LAMINAR RADIAL JETS WITH SWIRL R. Naz 1,
More informationThe Plane Stress Problem
. 4 The Plane Stress Problem 4 Chapter 4: THE PLANE STRESS PROBLEM 4 TABLE OF CONTENTS Page 4.. INTRODUCTION 4 3 4... Plate in Plane Stress............... 4 3 4... Mathematical Model.............. 4 4
More informationPerturbation Theory for Variational Inference
Perturbation heor for Variational Inference Manfred Opper U Berlin Marco Fraccaro echnical Universit of Denmark Ulrich Paquet Apple Ale Susemihl U Berlin Ole Winther echnical Universit of Denmark Abstract
More informationINSTRUCTIONS TO CANDIDATES:
NATIONAL NIVERSITY OF SINGAPORE FINAL EXAMINATION FOR THE DEGREE OF B.ENG ME 444 - DYNAMICS AND CONTROL OF ROBOTIC SYSTEMS October/November 994 - Time Allowed: 3 Hours INSTRCTIONS TO CANDIDATES:. This
More information1.6 ELECTRONIC STRUCTURE OF THE HYDROGEN ATOM
1.6 ELECTRONIC STRUCTURE OF THE HYDROGEN ATOM 23 How does this wave-particle dualit require us to alter our thinking about the electron? In our everda lives, we re accustomed to a deterministic world.
More information2: Distributions of Several Variables, Error Propagation
: Distributions of Several Variables, Error Propagation Distribution of several variables. variables The joint probabilit distribution function of two variables and can be genericall written f(, with the
More informationAll parabolas through three non-collinear points
ALL PARABOLAS THROUGH THREE NON-COLLINEAR POINTS 03 All parabolas through three non-collinear points STANLEY R. HUDDY and MICHAEL A. JONES If no two of three non-collinear points share the same -coordinate,
More informationDIGITAL CORRELATION OF FIRST ORDER SPACE TIME IN A FLUCTUATING MEDIUM
DIGITAL CORRELATION OF FIRST ORDER SPACE TIME IN A FLUCTUATING MEDIUM Budi Santoso Center For Partnership in Nuclear Technolog, National Nuclear Energ Agenc (BATAN) Puspiptek, Serpong ABSTRACT DIGITAL
More informationarxiv: v1 [quant-ph] 22 Jan 2019
Article TWO-QUBITS I A LARGE-S EVIROMET arxiv:1901.07416v1 [quant-ph] 22 Jan 2019 Eliana Fiorelli 1,2, ID, Alessandro Cuccoli 3,4 and Paola Verrucchi 5,3,4 * 1 School of Physics and Astronomy, University
More informationroth t dive = 0 (4.2.3) divh = 0 (4.2.4) Chapter 4 Waves in Unbounded Medium Electromagnetic Sources 4.2 Uniform plane waves in free space
Chapter 4 Waves in Unbounded Medium 4. lectromagnetic Sources 4. Uniform plane waves in free space Mawell s equation in free space is given b: H rot = (4..) roth = (4..) div = (4..3) divh = (4..4) which
More information11.1 Double Riemann Sums and Double Integrals over Rectangles
Chapter 11 Multiple Integrals 11.1 ouble Riemann Sums and ouble Integrals over Rectangles Motivating Questions In this section, we strive to understand the ideas generated b the following important questions:
More informationEulerian Angles I. THE EULERIAN ANGLES. FIG. 1: First Rotation in the Euler Scheme. l2h2:euler2.tex. P(x,y,z)=P(x. N y. φ x x
Eulerian Angles I. THE EULERIA AGLES There is no question that once we leave diatomic molecules, life gets difficult. Most chemistr concerns polatomics, so sticking with diatomics is possibl ecessivel
More information4Cubic. polynomials UNCORRECTED PAGE PROOFS
4Cubic polnomials 4.1 Kick off with CAS 4. Polnomials 4.3 The remainder and factor theorems 4.4 Graphs of cubic polnomials 4.5 Equations of cubic polnomials 4.6 Cubic models and applications 4.7 Review
More informationQuasi-Stationary Simulation: the Subcritical Contact Process
Brazilian Journal of Physics, vol. 36, no. 3A, September, 6 685 Quasi-Stationary Simulation: the Subcritical Contact Process Marcelo Martins de Oliveira and Ronald Dickman Departamento de Física, ICEx,
More information1.4 Optimization of an orthogonally stiffened plate considering fatigue constraints
1.4 Optimization of an orthogonall stiffened plate considering fatigue constraints Luis M.C. Simões 1, József Farkas and Karol Jarmai 1 Universit of Coimbra, Portugal,lcsimoes@dec.uc.pt Universit of Miskolc,
More informationThe dynamics of small particles whose size is roughly 1 µmt or. smaller, in a fluid at room temperature, is extremely erratic, and is
1 I. BROWNIAN MOTION The dynamics of small particles whose size is roughly 1 µmt or smaller, in a fluid at room temperature, is extremely erratic, and is called Brownian motion. The velocity of such particles
More informationChapter 2 Overview of the Simulation Methodology
Chapter Overview of the Simulation Methodolog. Introduction This chapter presents an overview of the simulation methodolog that comprises both the art and science involved in simulating phsical phenomena.
More informationMAE 323: Chapter 4. Plane Stress and Plane Strain. The Stress Equilibrium Equation
The Stress Equilibrium Equation As we mentioned in Chapter 2, using the Galerkin formulation and a choice of shape functions, we can derive a discretized form of most differential equations. In Structural
More informationAE/ME 339. K. M. Isaac Professor of Aerospace Engineering. December 21, 2001 topic13_grid_generation 1
AE/ME 339 Professor of Aerospace Engineering December 21, 2001 topic13_grid_generation 1 The basic idea behind grid generation is the creation of the transformation laws between the phsical space and the
More informationCHAPTER 2: Partial Derivatives. 2.2 Increments and Differential
CHAPTER : Partial Derivatives.1 Definition of a Partial Derivative. Increments and Differential.3 Chain Rules.4 Local Etrema.5 Absolute Etrema 1 Chapter : Partial Derivatives.1 Definition of a Partial
More information4452 Mathematical Modeling Lecture 13: Chaos and Fractals
Math Modeling Lecture 13: Chaos and Fractals Page 1 442 Mathematical Modeling Lecture 13: Chaos and Fractals Introduction In our tetbook, the discussion on chaos and fractals covers less than 2 pages.
More informationN coupled oscillators
Waves 1 1 Waves 1 1. N coupled oscillators towards the continuous limit. Stretched string and the wave equation 3. The d Alembert solution 4. Sinusoidal waves, wave characteristics and notation T 1 T N
More information4 Strain true strain engineering strain plane strain strain transformation formulae
4 Strain The concept of strain is introduced in this Chapter. The approimation to the true strain of the engineering strain is discussed. The practical case of two dimensional plane strain is discussed,
More informationDerivatives 2: The Derivative at a Point
Derivatives 2: The Derivative at a Point 69 Derivatives 2: The Derivative at a Point Model 1: Review of Velocit In the previous activit we eplored position functions (distance versus time) and learned
More informationFunctions with orthogonal Hessian
Functions with orthogonal Hessian B. Dacorogna P. Marcellini E. Paolini Abstract A Dirichlet problem for orthogonal Hessians in two dimensions is eplicitly solved, by characterizing all piecewise C 2 functions
More informationSubmarine sand ripples formation in a viscous fluid: 2D and 3D linear stability analysis
Marine Sandwave and River Dune Dnamics 1 & April 4 - Enschede, the Netherlands Submarine sand ripples formation in a viscous fluid: D and 3D linear stabilit analsis V. Langlois (1) and A. Valance (1) Groupe
More informationEffective dynamics of many-particle systems with dynamical constraint
Michael Herrmann Effective dnamics of man-particle sstems with dnamical constraint joint work with Barbara Niethammer and Juan J.L. Velázquez Workshop From particle sstems to differential equations WIAS
More informationEstimation Of Linearised Fluid Film Coefficients In A Rotor Bearing System Subjected To Random Excitation
Estimation Of Linearised Fluid Film Coefficients In A Rotor Bearing Sstem Subjected To Random Ecitation Arshad. Khan and Ahmad A. Khan Department of Mechanical Engineering Z.. College of Engineering &
More informationSolutions to Two Interesting Problems
Solutions to Two Interesting Problems Save the Lemming On each square of an n n chessboard is an arrow pointing to one of its eight neighbors (or off the board, if it s an edge square). However, arrows
More informationA spatially second-order accurate implicit numerical method for the space and time fractional Bloch-Torrey equation
Noname manuscript No. (will be inserted b the editor) A spatiall second-order accurate implicit numerical method for the space and time fractional Bloch-Torre equation J. Song Q. Yu F. Liu I. Turner Received:
More informationInteraction Network Analysis
CSI/BIF 5330 Interaction etwork Analsis Young-Rae Cho Associate Professor Department of Computer Science Balor Universit Biological etworks Definition Maps of biochemical reactions, interactions, regulations
More informationMath 214 Spring problem set (a) Consider these two first order equations. (I) dy dx = x + 1 dy
Math 4 Spring 08 problem set. (a) Consider these two first order equations. (I) d d = + d (II) d = Below are four direction fields. Match the differential equations above to their direction fields. Provide
More informationOn Information and Sufficiency
On Information and Sufficienc Huaiu hu SFI WORKING PAPER: 997-02-04 SFI Working Papers contain accounts of scientific work of the author(s) and do not necessaril represent the views of the Santa Fe Institute.
More informationThe first change comes in how we associate operators with classical observables. In one dimension, we had. p p ˆ
VI. Angular momentum Up to this point, we have been dealing primaril with one dimensional sstems. In practice, of course, most of the sstems we deal with live in three dimensions and 1D quantum mechanics
More informationSchool of Computer and Communication Sciences. Information Theory and Coding Notes on Random Coding December 12, 2003.
ÉCOLE POLYTECHNIQUE FÉDÉRALE DE LAUSANNE School of Computer and Communication Sciences Handout 8 Information Theor and Coding Notes on Random Coding December 2, 2003 Random Coding In this note we prove
More informationChapter 11 Three-Dimensional Stress Analysis. Chapter 11 Three-Dimensional Stress Analysis
CIVL 7/87 Chapter - /39 Chapter Learning Objectives To introduce concepts of three-dimensional stress and strain. To develop the tetrahedral solid-element stiffness matri. To describe how bod and surface
More informationUniform continuity of sinc x
Uniform continuit of sinc Peter Haggstrom www.gotohaggstrom.com mathsatbondibeach@gmail.com August 3, 03 Introduction The function sinc = sin as follows: is well known to those who stud Fourier theor.
More informationThe standard form of the equation of a circle is based on the distance formula. The distance formula, in turn, is based on the Pythagorean Theorem.
Unit, Lesson. Deriving the Equation of a Circle The graph of an equation in and is the set of all points (, ) in a coordinate plane that satisf the equation. Some equations have graphs with precise geometric
More informationLévy stable distribution and [0,2] power law dependence of. acoustic absorption on frequency
Lév stable distribution and [,] power law dependence of acoustic absorption on frequenc W. Chen Institute of Applied Phsics and Computational Mathematics, P.O. Box 89, Division Box 6, Beijing 88, China
More informationMultiple Integration
Contents 7 Multiple Integration 7. Introduction to Surface Integrals 7. Multiple Integrals over Non-rectangular Regions 7. Volume Integrals 4 7.4 Changing Coordinates 66 Learning outcomes In this Workbook
More informationS-I-R Model with Directed Spatial Diffusion
S-I-R Model with Directed Spatial Diffusion Fabio A. Milner Ruijun Zhao Department of Mathematics, Purdue Universit, West Lafaette, IN 4797-267 Abstract A S-I-R epidemic model is described in which susceptible
More informationDirectional derivatives and gradient vectors (Sect. 14.5). Directional derivative of functions of two variables.
Directional derivatives and gradient vectors (Sect. 14.5). Directional derivative of functions of two variables. Partial derivatives and directional derivatives. Directional derivative of functions of
More informationUnit 12 Study Notes 1 Systems of Equations
You should learn to: Unit Stud Notes Sstems of Equations. Solve sstems of equations b substitution.. Solve sstems of equations b graphing (calculator). 3. Solve sstems of equations b elimination. 4. Solve
More information3.0 PROBABILITY, RANDOM VARIABLES AND RANDOM PROCESSES
3.0 PROBABILITY, RANDOM VARIABLES AND RANDOM PROCESSES 3.1 Introduction In this chapter we will review the concepts of probabilit, rom variables rom processes. We begin b reviewing some of the definitions
More informationThe distortion observed in the bottom channel of Figure 1 can be predicted from the full transport equation, C t + u C. y D C. z, (1) x D C.
1 8. Shear Dispersion. The transport models and concentration field solutions developed in previous sections assume that currents are spatially uniform, i.e. u f(,y,). However, spatial gradients of velocity,
More informationDEVELOPMENT AND USE OF OPERATIONAL MODAL ANALYSIS
Università degli Studi di Napoli Federico II DEVELOPMENT AND USE OF OPERATIONAL MODAL ANALYSIS Examples of Civil, Aeronautical and Acoustical Applications Francesco Marulo Tiziano Polito francesco.marulo@unina.it
More informationComputation of Total Capacity for Discrete Memoryless Multiple-Access Channels
IEEE TRANSACTIONS ON INFORATION THEORY, VOL. 50, NO. 11, NOVEBER 2004 2779 Computation of Total Capacit for Discrete emorless ultiple-access Channels ohammad Rezaeian, ember, IEEE, and Ale Grant, Senior
More informationTHE HEATED LAMINAR VERTICAL JET IN A LIQUID WITH POWER-LAW TEMPERATURE DEPENDENCE OF DENSITY. V. A. Sharifulin.
THE HEATED LAMINAR VERTICAL JET IN A LIQUID WITH POWER-LAW TEMPERATURE DEPENDENCE OF DENSITY 1. Introduction V. A. Sharifulin Perm State Technical Universit, Perm, Russia e-mail: sharifulin@perm.ru Water
More informationarxiv: v1 [hep-ph] 19 Jan 2018
version. Eclusivit of p-n interactions within p-a collisions Thomas A. Trainor CENPA 359, Universit of Washington, Seattle, WA 9895 (Dated: October 3, 8) A geometric Glauber model applied to p-p or p-a
More informationLecture 5: Kinetic theory of fluids
Lecture 5: Kinetic theory of fluids September 21, 2015 1 Goal 2 From atoms to probabilities Fluid dynamics descrines fluids as continnum media (fields); however under conditions of strong inhomogeneities
More informationTime-Frequency Analysis: Fourier Transforms and Wavelets
Chapter 4 Time-Frequenc Analsis: Fourier Transforms and Wavelets 4. Basics of Fourier Series 4.. Introduction Joseph Fourier (768-83) who gave his name to Fourier series, was not the first to use Fourier
More informationVelocity Limit in DPD Simulations of Wall-Bounded Flows
Velocit Limit in DPD Simulations of Wall-Bounded Flows Dmitr A. Fedosov, Igor V. Pivkin and George Em Karniadakis Division of Applied Mathematics, Brown Universit, Providence, RI 2912 USA Abstract Dissipative
More information