The Squirmer model and the Boundary Element Method
|
|
- Magnus Sutton
- 5 years ago
- Views:
Transcription
1 Hauptseminar : Active Matter The Squirmer model and the Boundary Element Method Miru Lee April 6, 017 Physics epartment, University of Stuttgart 1 Introduction A micro scale swimmer exhibits directional motion in the fluid without any external force and torque. Our main goal is to understand and explain mathematically how it moves at very low Reynolds number in a Newtonian fluid. Micro-organisms such as Opalina moving via mechanical motions of densely packed cilia covering the whole surface can be modeled as if they swim via a motion of an "envelope" covering the surface as shown in figures 1a and 1b. The motion of the envelope, which gives rise to local fluid motion on the surface (i.e., the boundary conditions of the flow field), is called squirming motion [1,, and the swimmers that can be understood by means of the squirming motion are called squirmers. We will show an analytic solution to the squirmer model problem under the assumptions that the flow field due to the squirmer is axisymmetric and the shape of the body is a sphere. Later, we will introduce the rigid body and the far-field approximations to gain intuition about the flow field. This will allow us to identify four distinct types of squirmer. On the other hand, a chemically active particle can also be a micro-swimmer. Basic mechanism of which a chemically active particle may take advantage is the so-called phoretic mechanism; a local flow field in the vicinity of the particle s surface, which is created by the short-range interaction between the surface and local field gradients around the particle (e.g., solute density or temperature field), drives the motion of the particle. If the particle is able to induce the local chemical gradients around its surface itself, it can move autonomously [3, 4. We will focus on self-diffusiophoresis to explain the directional motion of a chemically active particle, and we will show that the mathematical structure of a simplified model for self-diffusiophoresis is very similar with that of the squirmer even though the basic mechanisms seem very different from each other. Therefore, this model can be regarded as a chemical version of the squirmer. If a particle, however, does not possess a particularly simple shape (i.e., spherical shape, axial symmetry, etc.), analytic solutions cannot be found, and thus numerical approaches have to be considered. There are several numerical methods to solve partial differential equations, and what we will present here is the boundary element method (BEM). The main idea of the method is that; 1) find a fundamental solution of the relevant equation, describing a field due to a point source, ) express the general solution in terms of boundary distributions of the fundamental solutions, and 3) compute the densities of the distributions to satisfy the given boundary conditions [5. We will use Laplace s equation as an example to explain the methodology of BEM because Laplace s equation is not only relevant to the chemically active particle but also less cumbersome than the Stokes equation, though the letter can also be solved using the same methodology. Squirmer model The squirmer at low Reynolds number in the incompressible flow of the unbounded Newtonian fluid is assumed to induce the axisymmetric flow and to have a sphere shape in order to simplify the 1
2 (a) Opalina: a micro-organism having cilia around the J. Fluid Mech. (1971) body. (b) The squirming motion: the cilia motions are assumed to be a periodic wave motion [1 Figure 1: Opalina (a) and the squirming motion of Opalina (b). mathematical structure. Then, the motion of the envelope is replaced by the flow field on the surface, that is, a boundary condition [. Before dealing with the problem, we have to choose a frame. Figure shows the two possible options: the fixed laboratory frame and the so-moving frame, attached to the moving swimmer. We will work in the co-moving frame, and employ spherical coordinates..1 Stream function At low Reynolds number, the Stokes equation for the incompressible flow of the Newtonian fluid is given by: u 1 p = 0 (1) µ where u and p are the flow field and the hydrostatic pressure, and µ is the viscosity of the fluid. For an axisymmetric flow field with the axis of symmetry being the z axis, one can introduce the stream function ψ(r, θ) in terms of which the flow fields is [6: u r (r, θ) = 1 r ψ (cos θ) The flow field does not depend on the azimuthal angle φ. Taking the curl on Eq.(1) and using p = 0 leads to Substituting Eq.() into Eq.(3) gives where, in spherical coordinates, the operator E is given by The general solution of Eqs.(4) and (5) is [6: Figure : The laboratory frame and the comoving frame. u θ (r, θ) = 1 ψ r r. () u = 0. (3) E (E ψ) = 0, (4) E = r + 1 cos θ r (cos θ). (5) ψ(r, θ) = (α n r n + β n r n+1 + γ n r n+ + δ n r n+3 )G n (cos θ), n=0 + (α nr n + β nr n+1 + γ nr n+ + δ nr n+3 )H n (cos θ) n= (6)
3 where G n and H n are Gegenbauer functions of the first kind and the second kind, respectively. 1 The two functions can also be re-presented in terms of Legendre polynomials of the first and the second kinds, P n and Q n : G n (cos θ) = P n (cos θ) P n (cos θ), H n (cos θ) = Q n (cos θ) Q n (cos θ) n 1 n 1 for n. (7) For n, G n and H n are defined as. Flow field Using G 0 (cos θ) = H 1 (cos θ) = 1, G 1 (cos θ) = H 0 (cos θ) = cos θ. (8) P n 1 (cos θ) = dg n(cos θ) d(cos θ), Q n 1(cos θ) = dh n(cos θ) d(cos θ) and substituting Eq.(6) in to Eq.() leads to u r (r, θ) = (α n r n + β n r n 1 + γ n r n + δ n r n+1 )P n 1 (cos θ) n=1 (α nr n + β nr n 1 + γ nr n + δ nr n+1 )Q n 1 (cos θ), n= u θ (r, θ) = [nα n r n (n 1)β n r n 1 + (n + )γ n r n (n 3)δ n r n+1 G n(cos θ) n=0 + [nα nr n (n 1)β nr n 1 + (n + )γ nr n (n 3)δ nr n+1 H n(cos θ). n= (9) For our squirmer, the flow field should be a bounded function at every point. Legendre functions of the second kind and Gegenbauer functions of the second kind diverge at cos θ = ±1. In addition, G 0 (cos θ) and G 1(cos θ) become infinite at = 0. Hence, in order for u(r, θ) to be finite for all θ, the related coefficients must vanish [6, i.e., α n = β n = γ n = δ n = 0 for all n, β 0 = γ 0 = δ 0 = 0, α 1 = γ 1 = δ 1 = 0. (10).3 Boundary conditions The system is a particle in an unbounded fluid. Thus, we have boundary conditions at infinity and on the surface of the particle. We begin with the boundary conditions at infinity. The fluid at infinity in the co-moving frame is moving in the negative z-direction with the velocity U. Consequently, the corresponding boundary condition is [6 which implies u(r ) Uê z, (11) { U for n =, α n = 0 otherwise, and γ n = 0 for all n. (1) 1 Note that the summation of the second line on Eq.(6) starts from n =. This is because of the definitions of G n and H n at n = 0 and n = 1. Gn(cos θ) For n, G n(cos θ) always contains term, so is well behaved. 3
4 As a result, the flow field becomes ( U + β u r (r, θ) = β 1 r P 0(cos θ) r 3 + δ r ( u θ (r, θ) = U β r 3 + δ ) G (cos θ) r n= ) P 1 (cos θ) ( βn+1 n= ( nβn+1 r n+ + (n )δ n+1 r n ) P n (cos θ), r n+ + δ n+1 r n ) Gn+1 (cos θ). The squirmer model is defined by the velocities on the surface as functions of θ, i.e., u r (θ) = f(θ) and u θ (θ) = g(θ), with f and g known functions. Since it is possible to expand any well behaving function of θ in terms of P n (cos θ) and G n (cos θ)/, the flow on the surface can be written as [ u r (a, θ) = f(θ) = A n P n (cos θ), u θ (a, θ) = g(θ) = B n V n (cos θ) (14) n=0 n=1 where V n (cos θ) = G n+1(cos θ), and {A n 0 } and {B n 1 } are known quantities. Evaluating Eq.(13) at r = a and matching with Eq.(14) provides the rest of the coefficients {β n } and {δ n }: [( ) n β 1 = a A 0, β n+1 = a n+ 1 A n B n, δ n+1 = an ( β = a3 3 A 1 + B 1 + δ a ), U = 1 3 A B 1 δ 3a. (na n B n ), However, β and δ are still unknown because U is undetermined yet. Thus, we need an additional condition..4 Force-free condition Let us consider the net force and the net torque acting on the swimmer. Given that the swimmer moves with the constant velocity, the net force and the net torque on it must vanish [7, 8. There are two forces and torques acting on the body: the external force and torque (F ext and T ext ) and the hydrodynamic force and torque (F hydro and T hydro ), hence F net = F ext + F hydro = 0, T net = T ext + T hydro = 0. Since the external force and torque are zero, the hydrodynamic force and torque, which are given by [6 F hydro = Π ˆndS, T hydro = r (Π ˆn)dS, (17) have to be zero. Using (13), one can find that the first term of Eq.(17) reduces to At the same time, due to the axial symmetry, (13) (15) (16) F z = 4πµδ. (18) T hydro 0. The force-free condition thus implies δ = 0. Putting everything together, one arrives at: [ u r (r, θ) = a r A a 3 0P 0 (cos θ) + 3 r 3 (A 1 + B 1 ) U P 1 (cos θ) [( n a n + r n n a n+ ) ( a n+ an r n+ A n + )B rn+ r n n P n (cos θ), [ 1 a 3 u θ (r, θ) = 3 r 3 (A 1 + B 1 ) + U V 1 (cos θ) + 1 ) ( [(n an+ n(n ) a n (n )an rn+ r n B n + )A r n an+ r n+ n V n (cos θ), U = 1 3 (B 1 A 1 ). (19) 4
5 Note that the force and torque free condition results in vanishing 1/r term, making 1/r term be the leading order..5 The laboratory frame The flow field in the lab frame is given by: v(r) = u(r) + U + R Ω (0) where R r + r c (see figure 3). In Eq.(0), Ω is zero due to the axisymmetric flow field, leaving v(r) = u(r) + U. (1) Accordingly, the flow field in the lab frame is given by: Figure 3: The relation between the lab frame and the co-moving frame. v r (R) = a r A 0P 0 (cos θ) + 3 r 3 (A 1 + B 1 )P 1 (cos θ) [( n a n + r n n a n+ ) ( a n+ r n+ A n + a 3 a 3 v θ (R) = 1 3 r 3 (A 1 + B 1 )V 1 (cos θ) + 1 [(n an+ U = 1 3 (B 1 A 1 ). (n )an rn+ r n ) B n + an rn+ r n )B n P n (cos θ) ( n(n ) a n )A r n an+ r n+ n V n (cos θ) () Note that R r, unless one chooses the origin of the lab frame to coincide with the current center of the squirmer..6 Reduced models We will consider examples of which there is only a tangential velocity at the surface of the particle, which would cover the case of an impermeable, rigid sphere with just tangential motion, that is, slip. This leads to u r (a, θ) = n=0 A n P n (cos θ) = 0 A n = 0 for all n. (3) Additionally, since each mode B n contributes to the flow field with the order of 1/r n and 1/r n+, the first and the second term (i.e., B 1, B modes) dominates the flow field at a point far away from the body; B 1 mode corresponds to a uniform flow, while B mode accounts for a disturbance. Thus, for the sake of gaining physical intuition, we will only examine B 1 and B modes. Accordingly, the flow field in the co-moving frame is approximately given by ( a 3 u r (r, θ) = 3 r 3 P 1(cos θ) ) ( a 4 ) 3 cos θ B 1 + r 4 a r B P (cos θ), ( 1 a 3 u θ (r, θ) = 3 r 3 V 1(cos θ) + ) 3 B 1 + a4 r 4 B V (cos θ), U = 3 B 1. Since B 1 is directly connected to U, it is clear that B determines the appearance of the flow field. Hence, by introducing the dimensionless parameter ξ := B /B 1, one can evaluate a configuration of the flow field; given ξ, the appearance of the flow field will be the same regardless of the value of B 1. (4) 5
6 (a) Puller (b) Pusher (c) Neutral (d) Shaker Figure 4: Flow field due to a squirmer for various values of ξ. The top row: in the co-moving frame. The bottom row : in the lab frame. [( a 3 u r (r, θ) = B 1 3 r 3 P 1(cos θ) ) ( a 4 ) 3 cos θ + ξ r 4 a r P (cos θ), [( 1 a 3 u θ (r, θ) = B 1 3 r 3 V 1(cos θ) + ) 3 B 1 + ξ a4 r 4 V (cos θ). By varying ξ, four distinct flow configurations are established (see figure 4) [7: (5) ξ > 0 puller, ξ < 0 pusher, ξ = 0 neutral, ξ shaker. puller: it pulls the fluid to move forward (panel (a)). pusher: it pushes the fluid to move forward (panel (b)). neutral: it neither pushes nor pulls. It moves by sucking and releasing the fluid (panel (c)). shaker: it actuates the fluid without directional motion (panel (d))..7 Self-diffusiophoretic particle Self-diffusiophoresis of a Janus particle involves local fluid motion around the surface of the particle (i.e., the boundary condition of the flow field on the surface), and the fluid motion is generated by the chemical interaction between an anisotropically distributed solute and the surface, which the so-called phoretic mechanism. The gradients of the solute density are induced by the particle, e.g., releasing a solute from a certain fraction of its surface (as shown in figure 5) [4. For this model, we assume that the solute density c(r) reaches the steady-state in a very short time scale compared to the motion of the particle. Therefore, the solute density is governed by Laplace s equation [3: c(r) = 0. (6) Figure 5: Self-diffusiophoretic Janus particle releases a solute from parts of the surface to induce chemical composition gradients. [4 6
7 The boundary condition for the solute density c(r) on the surface (r = a) of the particle is given: { k on the active surface, ˆn c(a, θ) = 0 on the inert surface. where and k is the diffusion coefficient of the solute and the solute flux constant. normal vector of the surface. At infinity, the boundary condition is (7) ˆn denotes a The slip velocity at the surface of the particle is given by c(r ) = c. (8) v s = b c(r) r=a (9) where is the projection of the gradient onto the particle surface, leaving only the tangential component. Eq.(9) involves the effect of the chemical interaction between the particle and the anisotropically diffused solute in the vicinity of it. The surface mobility b, which accounts for the strength and the type (i.e., repulsive or attractive) of the interaction, is assumed to be constant on the surface [3, 4. We further assume that the self-diffusiophoresis is spherical, axisymmetric, and thus the corresponding flow field is also axisymmetric. Consequently, the particle moves along the z direction without rotation, that is, U = Uê z and Ω = 0. The Stokes equation is the same as that of the squirmer model, that is, Eq.(13) with the slip velocity Eq.(9) [4. Therefore, one already has the solution, from section.4. 3 Boundary Element Method If the particle has an arbitrary shape, one has to approach the problem numerically. The Boundary Elements Method (BEM) is a tool to solve partial differential equations. The working mechanism of BEM will be illustrated for the case of Laplace s equation. That of the Stokes equation is similar with Laplace s equation, but the details are more complex due to the vectorial character of the field. Nevertheless, one can use the same mathematical structure to solve the Stokes equation. 3.1 Green s function for Laplace s equation Since the Laplacian is a linear operator, the solution for Laplace s equation can be expressed in terms of a superposition of fundamental solutions; the fundamental solutions are Green s functions and Green s function dipoles. Therefore, understanding Green s function is of importance. A Green s function G(r, r 0 ) is defined as satisfying Laplace s equation except at a single point r 0, namely a singular point or a pole [5: G(r, r 0 ) + δ(r r 0 ) = 0. (30) Eq.(30) is Poisson s equation, and the Green s function describes a field at a point r due to a point source located at r 0. In three dimensional free space, it has the explicit form: 1 G(r, r 0 ) = 4π r r 0. (31) Figure 6: A field due to a pair of a source and a sink of identical strength; the white circle and the black circle denote the source and the sink respectively. The colormap indicates the strength of the field. 7
8 3. Green s function dipole Figure 6 shows a field due to a pair of two sources of identical strength but opposite sign (i.e., a source and a sink). When the distance between the source and the sink is negligible compared to the distance between the origin and at a point r, the source and the sink form the so-called "source dipole". The field due to a source dipole is shown in figure 6, and can be expressed as [5: V dip (r, r 0 ) = d G (r, r 0 ), (3) where d is the vector from the sink to the source. G (r, r 0 ) is called Green s function dipole, and mathematically can be derived by taking the gradient of the Green s function G(r, r 0 ) with respect to the coordinates the pole r 0 : G (r, r 0 ) := 0 G(r, r 0 ). (33) Note that the Green s function dipole is a vector. In free space, it is given by (form Eq.(31)): 3.3 Green s identities To understand how the general solution of Laplace s equation can be recasted into a form of a integral representation, it is helpful to understand Green s identities beforehand. Any two functions f(r) and g(r), which are at least twice differentiable, satisfy Green s first identity [5. g f = (g f) g f. (35) Changing the order of f and g leads to f g = (f g) f g. (36) By subtracting the two equations, one can obtain Green s second identity: G (r, r 0 ) = r r 0 4π r r 0 3. (34) Figure 7: An arbitrarily chosen volume V c. denotes the boundary surface. ˆn is unit normal vector pointing into V c [5. g f f g = (g f f g). (37) If the two functions satisfy Laplace s equation, the second identity becomes: (g f f g) = 0. (38) This is called the reciprocal relation. 3 Integrating the second identity over an arbitrary volume V c, as depicted in figure 7, leads to (g f f g)dv = (g f f g)dv = (g f f g) ˆndS. (39) V c V c Note that the Gauss theorem is applied in the second equality and the negative sign comes from the choice of ˆn (as shown in figure 7). 3 In the case of Stokes equation, the analogous relation is the Lorentz reciprocal theorem. 8
9 3.4 Integral representation Assuming that f(r) satisfies Laplace s equation and replacing g(r) with the Green s function G(r, r 0 ) in Eq.(37) leads to f(r) G(r, r 0 ) = (G(r, r 0 ) f(r) f(r) G(r, r 0 )). (40) Substituting Eq.(30) into Eq.(40) and integrating over the control volume results in f(r 0 )δ(r r 0 )dv = G(r, r 0 )[ˆn f(r)ds + f(r)[ˆn G(r, r 0 )ds. V c (41) If r 0 / V c, the left-hand side of Eq.(41) becomes zero, leaving 0 = G(r, r 0 )[ˆn f(r)ds + f(r)[ˆn G(r, r 0 )ds. (4) On the other hand, if r 0 V c, Eq.(41) gives f(r 0 ) = G(r, r 0 )[ˆn f(r)ds + f(r)[ˆn G(r, r 0 )ds. (43) Since the two arguments of the Green s function is interchangeable (as one can check from Eq.(31)), one can rewrite Eq.(43) as: f(r 0 ) = G(r 0, r)[ˆn f(r)ds + f(r)[ˆn G(r 0, r)ds. (44) Eq.(44) is a boundary-integral representation in terms of the boundary distributions of the Green s function and the Green s function dipole. The first term of Eq.(44) is called the single-layer potential, and the second term is called the double-layer potential [ Boundary integral equation To use Eq.(44), we need information about both sets of the boundary values, that is, f(r) and ˆn f(r) on the boundary. In practice, only one of two sets is known. For the case of the self-diffusiophoretic particle, ˆn c(r) on the boundary of the particle is given. However, obtaining the other set of the boundary values is possible by taking the limit of Eq.(44) such that r 0 approaches the boundary [5: lim f(r 0) = lim G(r 0, r)[ˆn f(r)ds + lim f(r)[ˆn G(r 0, r)ds. (45) r 0 r 0 r 0 The single-layer potential integral of Eq.(45) is continuous as the singular point approaches and crosses the boundary, but the double-layer potential integral is not [5. By separating the integral region of the double-layer potential as shown in figure 8, one arrives at: lim r 0 f(r)[ˆn G(r 0, r)ds = P V f(r)[ˆn G(r 0, r)ds + Pole f(r)[ˆn G(r 0, r)ds, (46) where P V denotes the so-called principal value integral, which integrates over the boundary surface except for the region where the pole is located as depicted in figure 8. The second term in Eq.(46), therefore, integrates the small area around the pole r 0, and it can be set to give the result [5: Pole f(r)[ˆn G(r 0, r)ds 1 f(r 0). (47) 9
10 Consequently, Eq.(45) turns to f(r 0 ) = + P V G(r 0, r)[ˆn f(r)ds f(r)[ˆn G(r 0, r)ds. (48) Since ˆn f(r) is known, Eq.(48) can be expressed compactly as f(r 0 ) = where P V f(r)[ˆn G(r 0, r)ds + Φ(r 0 ), Φ(r 0 ) G(r 0, r)[ˆn f(r)ds. (49) Figure 8: A schematic picture of Eq.(46). The principal value integral corresponds to integrating the boundary following the orange line (i.e., the boundary surface). Thus, it integrates over the whole boundary surface except for the area around the singular point. This equation is called a Fredholm integral equation of the second kind for the boundary values of f(r) [ How to use Eq.49 in a numerical scheme? To apply Eq.(49) in a numerical scheme, one has to implement a so-called "discretization". Figure 9 represents a discretized boundary. Each triangle is called a boundary element. The boundary elements can be realized such that the element is small enough so that the value of f(r) over that region can be approximated to be constant f(r 0 Σ n ) f n. (50) Therefore, taking Eq.(49) with r 0 Σ n for n = 1,,...N leads to Figure 9: A discretized boundary. f 1 = f =..., f N = N i=1 N i=1 N i=n f i a (1) i + Φ (1) i, f i a () i + Φ () i, f i a (N) i + Φ (N) i. (51) where a (k) i and Φ (k) i are defined as P V a (k) i = ˆn G(r 0 Σ k, r)ds, Φ (k) i Σ i = ˆn f(r) r Σi Σ i G(r 0 Σ k, r)ds. (5) Eq.(51) consists linear system of N equations for the N unknowns for f i. Accordingly, the discretized form of f(r) is determined. Knowing f(r), one has all the information needed for using Eq.(44), and therefore, f(r) at any point is computable by interpolation. 10
11 4 Summary We have explored the squirmer model, and BEM. The squirmer model is a model to explain the ciliated self-propelling micro-organisms at low Reynolds number. In order to understand such kind of swimmers, we approximated the flow produced by the cilia motion with a "squirming" fluid velocity, providing the boundary conditions on the surface. We solved the Stokes equation analytically, in the axisymmetric flow field with a spherical shaped body. Then, introducing the rigid body and the far-field approximations allows us to identify the four distinct types of squirmer. We also studied self-diffusiophoretic particles, which is driven by the phoretic mechanism, and showed that the simple basic model of self-diffusiophoresis reduces to a squirmer model. Lastly, for the case of which an analytic approach is not available, we introduced and explained the boundary element method. The first step was to re-expressing the function of our interest into boundary integral representation in terms of two sets of the boundary values of the function. In practice, only one set is given, so the other set has to be obtained. In order to find the unknown set, we established an integral equation. It can be solved by discretizing the boundary into the boundary elements. After determining the two sets of the boundary values, one could compute the function at any point via the integral representation. Even though we illustrated the methodology of Laplace s equation, the same structure can be applied to solving the Stokes equation. References [1 J. Blake, A spherical envelope approach to ciliary propulsion, Journal of Fluid Mechanics, vol. 46, pp , [ M. Lighthill, On the squirming motion of nearly spherical deformable bodies through liquids at very small Reynolds numbers, Communications on Pure and Applied Mathematics, vol. 5, pp , 195. [3 W. Uspal, M. N. Popescu, S. ietrich, and M. Tasinkevych, Self-propulsion of a catalytically active particle near a planar wall: from reflection to sliding and hovering, Soft Matter, vol. 11, pp , 015. [4 S. Michelin and E. Lauga, Phoretic self-propulsion at finite Péclet numbers, Journal of Fluid Mechanics, vol. 747, pp , 014. [5 C. Pozrikidis, A practical guide to boundary element methods with the software library BEMLIB. CRC Press, 00. [6 J. Happel and H. Brenner, Low Reynolds number hydrodynamics: with special applications to particulate media. Springer Science & Business Media, 01, vol. 1. [7 M. T. ownton and H. Stark, Simulation of a model microswimmer, Journal of Physics: Condensed Matter, vol. 1, p , 009. [8 O. S. Pak and E. Lauga, Generalized squirming motion of a sphere, Journal of Engineering Mathematics, vol. 88, pp. 1 8,
arxiv: v1 [cond-mat.soft] 6 Apr 2018
Effective squirmer models for self-phoretic chemically active spherical colloids M.N. Popescu, 1,2, W.E. Uspal, 1,2 Z. Eskandari, 1,2 M. Tasinkevych, 3 and S. Dietrich 1,2 1 Max-Planck-Institut für Intelligente
More informationSupplemental Information for Self-propulsion of a catalytically active particle near a planar wall: from reflection to sliding and hovering
Electronic Supplementary Material (ESI) for Soft Matter. This journal is The Royal Society of Chemistry 1 Supplemental Information for Self-propulsion of a catalytically active particle near a planar wall:
More informationParticle self-diffusiophoresis near solid walls and interfaces
Particle self-diffusiophoresis near solid walls and interfaces Darren G. CROWDY Tel.: +44(0)2075948587; Fax: +44(0)2075948517; Email: d.crowdy@imperial.ac.uk Department of Mathematics, Imperial College
More informationFluid dynamics is the study of flows of fluid. However, before we begin discussing the dynamics of
Hydrodynamics of Newtonian Fluids Cameron Stewart Universität Stuttgart cameron.n.stewart@gmail.com April 5, 2017 Abstract In this review we provide an overview of the hydrodynamics of viscous, Newtonian
More informationarxiv: v1 [physics.flu-dyn] 14 Mar 2014
Phoretic self-propulsion at finite Péclet numbers arxiv:143.361v1 [physics.flu-dyn] 14 Mar 214 Sébastien Michelin 1, and Eric Lauga 2, 1 LadHyX Département de Mécanique, Ecole Polytechnique CNRS, 91128
More informationarxiv: v1 [cond-mat.soft] 27 May 2015
Self-Diffusiophoretic Colloidal Propulsion Near a Solid Boundary arxiv:55.772v cond-mat.soft] 27 May 25 Ali Mozaffari,, Nima Sharifi-Mood, 2, Joel Koplik, 3 and Charles Maldarelli, Benjamin Levich Institute
More informationarxiv: v1 [cond-mat.soft] 9 Aug 2016
arxiv:1608.02973v1 [cond-mat.soft] 9 Aug 2016 Chemically active colloids near osmotic-responsive walls with surface-chemistry gradients M N Popescu, W E Uspal and S Dietrich Max-Planck-Institut für Intelligente
More informationGeneral Solution of the Incompressible, Potential Flow Equations
CHAPTER 3 General Solution of the Incompressible, Potential Flow Equations Developing the basic methodology for obtaining the elementary solutions to potential flow problem. Linear nature of the potential
More informationSquirming Sphere in an Ideal Fluid
Squirming Sphere in an Ideal Fluid D Rufat June 8, 9 Introduction The purpose of this wor is to derive abalytically and with an explicit formula the motion due to the variation of the surface of a squirming
More informationMesoscale fluid simulation of colloidal systems
Mesoscale fluid simulation of colloidal systems Mingcheng Yang Institute of Physics, CAS Outline (I) Background (II) Simulation method (III) Applications and examples (IV) Summary Background Soft matter
More informationMotion of a model swimmer near a weakly deforming interface
Downloaded from https://www.cambridge.org/core. Purdue University Libraries, on 14 Jul 17 at 19:39:36, subject to the Cambridge Core terms of use, available at https://www.cambridge.org/core/terms. https://doi.org/1.117/jfm.17.85
More informationEmergence of collective dynamics in active biological systems -- Swimming micro-organisms --
12/08/2015, YITP, Kyoto Emergence of collective dynamics in active biological systems -- Swimming micro-organisms -- Norihiro Oyama John J. Molina Ryoichi Yamamoto* Department of Chemical Engineering,
More informationNumber of pages in the question paper : 05 Number of questions in the question paper : 48 Modeling Transport Phenomena of Micro-particles Note: Follow the notations used in the lectures. Symbols have their
More informationClassical Field Theory: Electrostatics-Magnetostatics
Classical Field Theory: Electrostatics-Magnetostatics April 27, 2010 1 1 J.D.Jackson, Classical Electrodynamics, 2nd Edition, Section 1-5 Electrostatics The behavior of an electrostatic field can be described
More informationBasic concepts in viscous flow
Élisabeth Guazzelli and Jeffrey F. Morris with illustrations by Sylvie Pic Adapted from Chapter 1 of Cambridge Texts in Applied Mathematics 1 The fluid dynamic equations Navier-Stokes equations Dimensionless
More informationSELF-DIFFUSIOPHORESIS AND BIOLOGICAL MOTILITY. Department of Physics & Astronomy University of Pennsylvania
SELF-DIFFUSIOPHORESIS AND BIOLOGICAL MOTILITY Department of Physics & Astronomy University of Pennsylvania Informal Seminar, University of Oxford, 8th June 2011 ACTIN BASED PROPULSION Listeria monocytogenes
More informationFluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition
Fluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition C. Pozrikidis m Springer Contents Preface v 1 Introduction to Kinematics 1 1.1 Fluids and solids 1 1.2 Fluid parcels and flow
More information2.29 Numerical Fluid Mechanics Fall 2011 Lecture 7
Numerical Fluid Mechanics Fall 2011 Lecture 7 REVIEW of Lecture 6 Material covered in class: Differential forms of conservation laws Material Derivative (substantial/total derivative) Conservation of Mass
More information6.2 Governing Equations for Natural Convection
6. Governing Equations for Natural Convection 6..1 Generalized Governing Equations The governing equations for natural convection are special cases of the generalized governing equations that were discussed
More information3 Chapter. Gauss s Law
3 Chapter Gauss s Law 3.1 Electric Flux... 3-2 3.2 Gauss s Law (see also Gauss s Law Simulation in Section 3.10)... 3-4 Example 3.1: Infinitely Long Rod of Uniform Charge Density... 3-9 Example 3.2: Infinite
More informationF11AE1 1. C = ρν r r. r u z r
F11AE1 1 Question 1 20 Marks) Consider an infinite horizontal pipe with circular cross-section of radius a, whose centre line is aligned along the z-axis; see Figure 1. Assume no-slip boundary conditions
More information1 Fundamentals. 1.1 Overview. 1.2 Units: Physics 704 Spring 2018
Physics 704 Spring 2018 1 Fundamentals 1.1 Overview The objective of this course is: to determine and fields in various physical systems and the forces and/or torques resulting from them. The domain of
More information1. Poisson-Boltzmann 1.1. Poisson equation. We consider the Laplacian. which is given in spherical coordinates by (2)
1. Poisson-Boltzmann 1.1. Poisson equation. We consider the Laplacian operator (1) 2 = 2 x + 2 2 y + 2 2 z 2 which is given in spherical coordinates by (2) 2 = 1 ( r 2 ) + 1 r 2 r r r 2 sin θ θ and in
More informationREE Internal Fluid Flow Sheet 2 - Solution Fundamentals of Fluid Mechanics
REE 307 - Internal Fluid Flow Sheet 2 - Solution Fundamentals of Fluid Mechanics 1. Is the following flows physically possible, that is, satisfy the continuity equation? Substitute the expressions for
More informationDetailed Outline, M E 521: Foundations of Fluid Mechanics I
Detailed Outline, M E 521: Foundations of Fluid Mechanics I I. Introduction and Review A. Notation 1. Vectors 2. Second-order tensors 3. Volume vs. velocity 4. Del operator B. Chapter 1: Review of Basic
More informationarxiv: v1 [physics.flu-dyn] 1 Feb 2015
Stokesian swimmers and active particles B. U. Felderhof Institut für Theorie der Statistischen Physik RWTH Aachen University Templergraben 55 52056 Aachen Germany arxiv:1502.00242v1 [physics.flu-dyn] 1
More informationChapter 6: Incompressible Inviscid Flow
Chapter 6: Incompressible Inviscid Flow 6-1 Introduction 6-2 Nondimensionalization of the NSE 6-3 Creeping Flow 6-4 Inviscid Regions of Flow 6-5 Irrotational Flow Approximation 6-6 Elementary Planar Irrotational
More informationp + µ 2 u =0= σ (1) and
Master SdI mention M2FA, Fluid Mechanics program Hydrodynamics. P.-Y. Lagrée and S. Zaleski. Test December 4, 2009 All documentation is authorized except for reference [1]. Internet access is not allowed.
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationarxiv: v1 [physics.flu-dyn] 29 Apr 2013
Accepted. 1 Low-Reynolds number swimming in a capillary tube arxiv:1304.7671v1 [physics.flu-dyn] 29 Apr 2013 L. Zhu 1, E. Lauga 2 and L. Brandt 1 1 Linné Flow Centre, KTH Mechanics, S-100 44 Stockholm,
More informationBäcklund Transformations: a Link Between Diffusion Models and Hydrodynamic Equations
Copyright 2014 Tech Science Press CMES, vol.103, no.4, pp.215-227, 2014 Bäcklund Transformations: a Link Between Diffusion Models and Hydrodynamic Equations J.R. Zabadal 1, B. Bodmann 1, V. G. Ribeiro
More informationDNS of colloidal dispersions using the smoothed profile method: formulation and applications
Hokusai, 1831 Workshop III: High Performance and Parallel Computing Methods and Algorithms for Multiphase/Complex Fluids Institute for Mathematical Sciences, NUS, Singapore 2 6 March 2015 DNS of colloidal
More informationJun-ichi Takada
Wave Theory II Numerical Simulation of Waves 2) Representation of Wave Function by Using Green Function I) Green s Theorem and Free Space Green Function Jun-ichi Takada takada@ide.titech.ac.jp)) This lecture
More informationElectrodynamics Exam Solutions
Electrodynamics Exam Solutions Name: FS 215 Prof. C. Anastasiou Student number: Exercise 1 2 3 4 Total Max. points 15 15 15 15 6 Points Visum 1 Visum 2 The exam lasts 18 minutes. Start every new exercise
More informationd 1 µ 2 Θ = 0. (4.1) consider first the case of m = 0 where there is no azimuthal dependence on the angle φ.
4 Legendre Functions In order to investigate the solutions of Legendre s differential equation d ( µ ) dθ ] ] + l(l + ) m dµ dµ µ Θ = 0. (4.) consider first the case of m = 0 where there is no azimuthal
More informationNina Aguillon 1, Astrid Decoene 1, Benoît Fabrèges 1, Bertrand Maury 1 and Benoît Semin 2
ESAIM: PROCEEDINGS, December 2012, Vol. 38, p. 36-53 F. Coquel, M. Gutnic, P. Helluy, F. Lagoutière, C. Rohde, N. Seguin, Editors MODELLING AND SIMULATION OF 2D STOKESIAN SQUIRMERS Nina Aguillon 1, Astrid
More informationChapter 5. Magnetostatics
Chapter 5. Magnetostatics 5.4 Magnetic Vector Potential 5.1.1 The Vector Potential In electrostatics, E Scalar potential (V) In magnetostatics, B E B V A Vector potential (A) (Note) The name is potential,
More informationProblem Set Number 01, MIT (Winter-Spring 2018)
Problem Set Number 01, 18.377 MIT (Winter-Spring 2018) Rodolfo R. Rosales (MIT, Math. Dept., room 2-337, Cambridge, MA 02139) February 28, 2018 Due Thursday, March 8, 2018. Turn it in (by 3PM) at the Math.
More informationarxiv: v2 [cond-mat.soft] 17 Dec 2015
Large-scale simulation of steady and time-dependent active suspensions with the force-coupling method Blaise Delmotte a,b,, Eric E. Keaveny c,, Franck Plouraboué a,b, Eric Climent a,b arxiv:5.9v [cond-mat.soft]
More informationThe purpose of this lecture is to present a few applications of conformal mappings in problems which arise in physics and engineering.
Lecture 16 Applications of Conformal Mapping MATH-GA 451.001 Complex Variables The purpose of this lecture is to present a few applications of conformal mappings in problems which arise in physics and
More informationReview of Vector Analysis in Cartesian Coordinates
Review of Vector Analysis in Cartesian Coordinates 1 Scalar: A quantity that has magnitude, but no direction. Examples are mass, temperature, pressure, time, distance, and real numbers. Scalars are usually
More informationFluid Dynamics Exercises and questions for the course
Fluid Dynamics Exercises and questions for the course January 15, 2014 A two dimensional flow field characterised by the following velocity components in polar coordinates is called a free vortex: u r
More informationSolutions to Laplace s Equations- II
Solutions to Laplace s Equations- II Lecture 15: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay Laplace s Equation in Spherical Coordinates : In spherical coordinates
More informationUNIVERSITY OF EAST ANGLIA
UNIVERSITY OF EAST ANGLIA School of Mathematics May/June UG Examination 2007 2008 FLUIDS DYNAMICS WITH ADVANCED TOPICS Time allowed: 3 hours Attempt question ONE and FOUR other questions. Candidates must
More informationFluid mechanics of nodal flow due to embryonic primary cilia: Electronic Supplementary Material
Fluid mechanics of nodal flow due to embryonic primary cilia: Electronic Supplementary Material D. J. Smith, J. R. Blake and E. A. Gaffney December 2, 27 The volume flow rate due to a conically rotating
More informationMagnetostatics. Lecture 23: Electromagnetic Theory. Professor D. K. Ghosh, Physics Department, I.I.T., Bombay
Magnetostatics Lecture 23: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay Magnetostatics Up until now, we have been discussing electrostatics, which deals with physics
More informationSTEADY VISCOUS FLOW THROUGH A VENTURI TUBE
CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 12, Number 2, Summer 2004 STEADY VISCOUS FLOW THROUGH A VENTURI TUBE K. B. RANGER ABSTRACT. Steady viscous flow through an axisymmetric convergent-divergent
More informationNumber of pages in the question paper : 06 Number of questions in the question paper : 48 Modeling Transport Phenomena of Micro-particles Note: Follow the notations used in the lectures. Symbols have their
More informationChapter 4. Electrostatic Fields in Matter
Chapter 4. Electrostatic Fields in Matter 4.1. Polarization 4.2. The Field of a Polarized Object 4.3. The Electric Displacement 4.4. Linear Dielectrics 4.5. Energy in dielectric systems 4.6. Forces on
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationEntropy generation and transport
Chapter 7 Entropy generation and transport 7.1 Convective form of the Gibbs equation In this chapter we will address two questions. 1) How is Gibbs equation related to the energy conservation equation?
More informationThe Kelvin transformation in potential theory and Stokes flow
IMA Journal of Applied Mathematics (2009 74, 427 438 doi:10.1093/imamat/hxn027 Advance Access publication on September 9, 2008 The Kelvin transformation in potential theory and Stokes flow GEORGE DASSIOS
More informationSteady and unsteady diffusion
Chapter 5 Steady and unsteady diffusion In this chapter, we solve the diffusion and forced convection equations, in which it is necessary to evaluate the temperature or concentration fields when the velocity
More informationChapter 5. Magnetostatics
Chapter 5. Magnetostatics 5.1 The Lorentz Force Law 5.1.1 Magnetic Fields Consider the forces between charges in motion Attraction of parallel currents and Repulsion of antiparallel ones: How do you explain
More informationCHAPTER 8 ENTROPY GENERATION AND TRANSPORT
CHAPTER 8 ENTROPY GENERATION AND TRANSPORT 8.1 CONVECTIVE FORM OF THE GIBBS EQUATION In this chapter we will address two questions. 1) How is Gibbs equation related to the energy conservation equation?
More informationThermocapillary Migration of a Drop
Thermocapillary Migration of a Drop An Exact Solution with Newtonian Interfacial Rheology and Stretching/Shrinkage of Interfacial Area Elements for Small Marangoni Numbers R. BALASUBRAMANIAM a AND R. SHANKAR
More informationHomework Two. Abstract: Liu. Solutions for Homework Problems Two: (50 points total). Collected by Junyu
Homework Two Abstract: Liu. Solutions for Homework Problems Two: (50 points total). Collected by Junyu Contents 1 BT Problem 13.15 (8 points) (by Nick Hunter-Jones) 1 2 BT Problem 14.2 (12 points: 3+3+3+3)
More informationContinuum Mechanics Lecture 5 Ideal fluids
Continuum Mechanics Lecture 5 Ideal fluids Prof. http://www.itp.uzh.ch/~teyssier Outline - Helmholtz decomposition - Divergence and curl theorem - Kelvin s circulation theorem - The vorticity equation
More informationMaxwell s equations for electrostatics
Maxwell s equations for electrostatics October 6, 5 The differential form of Gauss s law Starting from the integral form of Gauss s law, we treat the charge as a continuous distribution, ρ x. Then, letting
More informationThe 3 dimensional Schrödinger Equation
Chapter 6 The 3 dimensional Schrödinger Equation 6.1 Angular Momentum To study how angular momentum is represented in quantum mechanics we start by reviewing the classical vector of orbital angular momentum
More information2. FLUID-FLOW EQUATIONS SPRING 2019
2. FLUID-FLOW EQUATIONS SPRING 2019 2.1 Introduction 2.2 Conservative differential equations 2.3 Non-conservative differential equations 2.4 Non-dimensionalisation Summary Examples 2.1 Introduction Fluid
More informationWeek 2 Notes, Math 865, Tanveer
Week 2 Notes, Math 865, Tanveer 1. Incompressible constant density equations in different forms Recall we derived the Navier-Stokes equation for incompressible constant density, i.e. homogeneous flows:
More informationLegendre s Equation. PHYS Southern Illinois University. October 13, 2016
PHYS 500 - Southern Illinois University October 13, 2016 PHYS 500 - Southern Illinois University Legendre s Equation October 13, 2016 1 / 10 The Laplacian in Spherical Coordinates The Laplacian is given
More informationQuick Recapitulation of Fluid Mechanics
Quick Recapitulation of Fluid Mechanics Amey Joshi 07-Feb-018 1 Equations of ideal fluids onsider a volume element of a fluid of density ρ. If there are no sources or sinks in, the mass in it will change
More informationChapter 5. The Differential Forms of the Fundamental Laws
Chapter 5 The Differential Forms of the Fundamental Laws 1 5.1 Introduction Two primary methods in deriving the differential forms of fundamental laws: Gauss s Theorem: Allows area integrals of the equations
More informationCOLLECTIVE HYDRODYNAMICS OF SWIMMING MICROORGANISMS
COLLECTIVE HYDRODYNAMICS OF SWIMMING MICROORGANISMS Timothy J Pedley DAMTP, University of Cambridge APS-DFD, Salt Lake City, 2007 OUTLINE 1 Bioconvection Mechanisms Model Instability of a uniform suspension
More informationSpherical squirmers: models for swimming micro-organisms
IMA Journal of Applied Mathematics Advance Access published July 11, 2016 IMA Journal of Applied Mathematics (2016) 81, 488 521 doi:10.1093/imamat/hxw030 Spherical squirmers: models for swimming micro-organisms
More informationSupplementary Figures
Supplementary Figures Supplementary Fig. 1. Fabrication process of the micro-scallop. (a) The negative mold of the micro-scallop is made by 3-D printing. The mold for the hinge is much shallower and narrower
More informationHaus, Hermann A., and James R. Melcher. Electromagnetic Fields and Energy. Englewood Cliffs, NJ: Prentice-Hall, ISBN:
MIT OpenCourseWare http://ocw.mit.edu Haus, Hermann A., and James R. Melcher. Electromagnetic Fields and Energy. Englewood Cliffs, NJ: Prentice-Hall, 1989. ISBN: 9780132490207. Please use the following
More information2.29 Numerical Fluid Mechanics Fall 2011 Lecture 5
.9 Numerical Fluid Mechanics Fall 011 Lecture 5 REVIEW Lecture 4 Roots of nonlinear equations: Open Methods Fixed-point Iteration (General method or Picard Iteration), with examples Iteration rule: x g(
More informationlim = F F = F x x + F y y + F z
Physics 361 Summary of Results from Lecture Physics 361 Derivatives of Scalar and Vector Fields The gradient of a scalar field f( r) is given by g = f. coordinates f g = ê x x + ê f y y + ê f z z Expressed
More informationVortex motion. Wasilij Barsukow, July 1, 2016
The concept of vorticity We call Vortex motion Wasilij Barsukow, mail@sturzhang.de July, 206 ω = v vorticity. It is a measure of the swirlyness of the flow, but is also present in shear flows where the
More informationEquations of linear stellar oscillations
Chapter 4 Equations of linear stellar oscillations In the present chapter the equations governing small oscillations around a spherical equilibrium state are derived. The general equations were presented
More informationComputation of Unsteady Flows With Moving Grids
Computation of Unsteady Flows With Moving Grids Milovan Perić CoMeT Continuum Mechanics Technologies GmbH milovan@continuummechanicstechnologies.de Unsteady Flows With Moving Boundaries, I Unsteady flows
More informationFluid Mechanics II Viscosity and shear stresses
Fluid Mechanics II Viscosity and shear stresses Shear stresses in a Newtonian fluid A fluid at rest can not resist shearing forces. Under the action of such forces it deforms continuously, however small
More informationMAE 101A. Homework 7 - Solutions 3/12/2018
MAE 101A Homework 7 - Solutions 3/12/2018 Munson 6.31: The stream function for a two-dimensional, nonviscous, incompressible flow field is given by the expression ψ = 2(x y) where the stream function has
More informationHorizontal buoyancy-driven flow along a differentially cooled underlying surface
Horizontal buoyancy-driven flow along a differentially cooled underlying surface By Alan Shapiro and Evgeni Fedorovich School of Meteorology, University of Oklahoma, Norman, OK, USA 6th Baltic Heat Transfer
More information. (70.1) r r. / r. Substituting, we have the following equation for f:
7 Spherical waves Let us consider a sound wave in which the distribution of densit velocit etc, depends only on the distance from some point, ie, is spherically symmetrical Such a wave is called a spherical
More informationChapter 3: Newtonian Fluids
Chapter 3: Newtonian Fluids CM4650 Michigan Tech Navier-Stokes Equation v vv p t 2 v g 1 Chapter 3: Newtonian Fluid TWO GOALS Derive governing equations (mass and momentum balances Solve governing equations
More informationarxiv: v1 [physics.flu-dyn] 12 Sep 2016
Stresslets induced by active swimmers arxiv:609.0375v [physics.flu-dyn] Sep 06 Eric Lauga, and Sébastien Michelin, Department of Applied Mathematics and Theoretical Physics, University of Cambridge, CB3
More informationPEMP ACD2505. M.S. Ramaiah School of Advanced Studies, Bengaluru
Two-Dimensional Potential Flow Session delivered by: Prof. M. D. Deshpande 1 Session Objectives -- At the end of this session the delegate would have understood PEMP The potential theory and its application
More informationarxiv: v1 [physics.flu-dyn] 4 May 2010
Jet propulsion without inertia Saverio E. Spagnolie and Eric Lauga Department of Mechanical and Aerospace Engineering, University of California San Diego, 95 Gilman Drive, La Jolla CA 993-4. (Dated: May
More informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
More informationThe Navier-Stokes equations for an incompressible fluid in indicial notation are. i u i = 0 (9.1)
Chapter 9 Viscous flow The Navier-Stokes equations for an incompressible fluid in indicial notation are i u i = 0 9.1 ρ t + u j j u i = i p + µ 2 j u i 9.2 When these are applied to practical problems,
More informationASTR 320: Solutions to Problem Set 2
ASTR 320: Solutions to Problem Set 2 Problem 1: Streamlines A streamline is defined as a curve that is instantaneously tangent to the velocity vector of a flow. Streamlines show the direction a massless
More informationDielectrics - III. Lecture 22: Electromagnetic Theory. Professor D. K. Ghosh, Physics Department, I.I.T., Bombay
Dielectrics - III Lecture 22: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay We continue with our discussion of dielectric medium. Example : Dielectric Sphere in a uniform
More informationENGR Heat Transfer II
ENGR 7901 - Heat Transfer II External Flows 1 Introduction In this chapter we will consider several fundamental flows, namely: the flat plate, the cylinder, the sphere, several other body shapes, and banks
More informationThis is an author-deposited version published in : Eprints ID : 14290
Open Archive TOULOUSE Archive Ouverte (OATAO) OATAO is an open access repository that collects the work of Toulouse researchers and makes it freely available over the web where possible. This is an author-deposited
More informationStatistical Mechanics of Active Matter
Statistical Mechanics of Active Matter Umberto Marini Bettolo Marconi University of Camerino, Italy Naples, 24 May,2017 Umberto Marini Bettolo Marconi (2017) Statistical Mechanics of Active Matter 2017
More informationChapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
More informationReview of Electrostatics
Review of Electrostatics 1 Gradient Define the gradient operation on a field F = F(x, y, z) by; F = ˆx F x + ŷ F y + ẑ F z This operation forms a vector as may be shown by its transformation properties
More informationCONVECTIVE HEAT TRANSFER
CONVECTIVE HEAT TRANSFER Mohammad Goharkhah Department of Mechanical Engineering, Sahand Unversity of Technology, Tabriz, Iran CHAPTER 3 LAMINAR BOUNDARY LAYER FLOW LAMINAR BOUNDARY LAYER FLOW Boundary
More informationfree space (vacuum) permittivity [ F/m]
Electrostatic Fields Electrostatic fields are static (time-invariant) electric fields produced by static (stationary) charge distributions. The mathematical definition of the electrostatic field is derived
More informationBoundary Conditions in Fluid Mechanics
Boundary Conditions in Fluid Mechanics R. Shankar Subramanian Department of Chemical and Biomolecular Engineering Clarkson University The governing equations for the velocity and pressure fields are partial
More informationSupplementary Information. Directed flow of micromotors through alignment interactions with micropatterned ratchets
Supplementary Information Directed flow of micromotors through alignment interactions with micropatterned ratchets Jaideep Katuri 1,2, David Caballero 1,3,4, Raphael Voituriez 5,6, Josep Samitier 1,3,4,
More informationSome Basic Plane Potential Flows
Some Basic Plane Potential Flows Uniform Stream in the x Direction A uniform stream V = iu, as in the Fig. (Solid lines are streamlines and dashed lines are potential lines), possesses both a stream function
More informationIII.3 Momentum balance: Euler and Navier Stokes equations
32 Fundamental equations of non-relativistic fluid dynamics.3 Momentum balance: Euler and Navier tokes equations For a closed system with total linear momentum P ~ with respect to a given reference frame
More informationBioconvection revisited
Bioconvection revisited A progress report... Jean-Luc Thiffeault Peter Mueller Natalie Cook Department of Mathematics University of Wisconsin Madison AIMS Conference on Dynamical Systems, Orlando, FL 1
More informationComplex Analysis MATH 6300 Fall 2013 Homework 4
Complex Analysis MATH 6300 Fall 2013 Homework 4 Due Wednesday, December 11 at 5 PM Note that to get full credit on any problem in this class, you must solve the problems in an efficient and elegant manner,
More informationFORMULA SHEET FOR QUIZ 2 Exam Date: November 8, 2017
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.07: Electromagnetism II November 5, 207 Prof. Alan Guth FORMULA SHEET FOR QUIZ 2 Exam Date: November 8, 207 A few items below are marked
More information