A DLM/FD Method for Simulating Particle Motion in Viscoelastic Fluid
|
|
- Rosamond McDonald
- 6 years ago
- Views:
Transcription
1 A Method for Simulating Particle Motion in Viscoelastic Fluid Tsorng-Wha Pan and Roland Glowinski Department of Mathematics Universit of Houston Houston, Teas Ma 3, 6 International Workshop on Fluid-Structure Interaction Problems Institute for Mathematical Sciences, NUS 6 This work is supported b NSF. (UH) Ma 3, 6 / 58
2 Outline.. 3. (UH) Ma 3, 6 / 58
3 Fig. : (Left) Spheres falling in % polo in 98% water draft, kiss and chain. The chain since the forces in a viscoelastic fluid are aggregative. Chained spheres form a long bod which is stable when its long ais is vertical. (Middle) Spheres are falling with their long ais vertical and along the stream, stable and permanentl chained. (Right) Clinders falling in a viscoelastic fluid turn their long ais parallel to the flow b the elasticit, the eact opposite of a Newtonian fluid. (D.D. Jospeh. Interrogations of Direct Simulation of Solid-Liquid Flow. efluids, ). (UH) Ma 3, 6 3 / 58
4 Particle motion in viscoelastic fluid Ω Fig. : An eample of flow region with one ball in IR 3. Ω is a bounded region in IR d (d =, 3) containing a viscoelastic fluid of either Oldrod-B or FENE-CR tpe and a ball B(t). B(t) is the particle region with G(t) as the mass center. Γ = Ω; γ(t) = B(t). B (UH) Ma 3, 6 4 / 58
5 Fluid flow governing equations ( ) u ρ f t +(u )u =ρ f g p+µ D(u)+ σ p in Ω \ B(t), () u = in Ω \ B(t), () u(, ) = u (), Ω \ B(), (with u = ), (3) u = g on Γ, (with g n dγ = ), (4) Γ C t +(u )C ( u)c C( u)t = f (C) (C I) in Ω \ B(t), (5) λ C(, ) = C (), Ω \ B(), C = C L on Γ. (6) u is the flow velocit, p is the pressure, and g is the gravit. D(u) = u + ( u) t, (v )w = { d j= v j w i j } d i=. On Γ, Γ is the upstream part and n is the outer normal unit vector. (UH) Ma 3, 6 5 / 58
6 Fluid flow governing equations For the fluid, ρ f is the densit, η is the viscosit, µ = η λ /λ is the solvent viscosit, η = η µ is the elastic viscosit, λ is the relaation time, and λ is the retardation time. The polmeric stress tensor σ p in () is given b σ p = η λ f (C)(C I), where the conformation tensor C is smmetric and positive definite and I is the identit matri. Setting f equal to unit corresponds to the Oldrod-B model while L f (C) = L tr(c) corresponds to the FENE-CR model, where tr(c) is the trace of the conformation tensor C and L is the maimum etension of the immersed polmer coils and referred to as the etensibilit of the immersed polmer coils. The Oldrod-B model then is a special case associated with infinite etensibilit. (UH) Ma 3, 6 6 / 58
7 Particle motion The motion of particle B satisfies the Euler-Newton s equations: v(, t) = V(t) + ω(t) G, {, t} B(t) (, T ), (7) dg = V, (8) dt dv M p dt = M d(i p ω) p g + F H, = T H, (9) dt G() = G, V() = V, ω() = ω. () M p is the mass, I p is the inertia, G is the mass center. V is the velocit of the mass center, and ω is the particle angular velocit. The hdrodnamical forces and torque are F H = σn dγ, T H = G σn dγ. γ γ (UH) Ma 3, 6 7 / 58
8 formulation For simplicit, we consider one solid ball B(t) in the fluid. To obtain a fictitious domain formulation formulation for above problem () (), we define first the following functional spaces V g (t) = { v v (H (Ω)) d, v = g (t) on Γ}, L (Ω) = { q q L (Ω), q d = }, V CL (t) = { C C (H (Ω)) d d, C = C L (t) on Γ }, V C = { C C (H (Ω)) d d, C = on Γ }, Λ(t) = H (B(t)) d. For <, > B(t), it can be defined as follows: <µ, v> B(t) = (µ v+ µ: v)d, µ H (B(t)) d, v H (Ω) d. B(t) Ω (UH) Ma 3, 6 8 / 58
9 formulation A distributed Lagrange multiplier/fictitious domain () formulation of the equations ()-() reads as follows: For a.e. t >, find u(t) V g (t), p(t) L (Ω), C(t) V C L (t), V(t) IR d, G(t) IR d, ω(t) IR d, λ(t) Λ(t) such that ρ f { u + (u )u} v d + µ D(u) : D(v) d Ω t Ω p v d v ( σ p ) d Ω Ω +( ρ dv f ρ s ){M p dt Y + I dω p dt θ} () =<λ, v Y θ G> B(t) +ρ f g vd+( ρ f )M p g Y, Ω ρ s {v, Y, θ} (H (Ω)) d IR d IR d, a.e. t (, T ), D particulate flow of Oldrod-B fluids: Singh, Joseph, Hesla, Glowinski, Pan. J. Non-Newtonian Fluid Mech. 9 () & Hao, Pan, Glowinski, Joseph. J. Non-Newtonian Fluid Mech. 56 (9). (UH) Ma 3, 6 9 / 58
10 Ω formulation q u(t) d =, q L (Ω), a.e. t (, T ), () <µ, u(, t) V(t) ω G> B(t) =, µ Λ(t), a.e. t (,T ), (3) ( ) C +(u )C ( u)c C( u)t :sd (4) Ω t f (C) = (C I):sd, s V C, and C = I in B(t), Ω λ dg dt = V, (5) C(, ) = C (), Ω, (6) G() = G, V() = V, ω() = ω, B() = B, (7) { u (), Ω \ B, u(, ) = V + ω G (8), B. (UH) Ma 3, 6 / 58
11 formulation In relation () we can replace D(u) : D(v)d b Ω u: vd due to that u is divergence free and satisfies the Ω Dirichlet boundar condition on Γ. For d =, ω(t) = {,, ω(t)} and θ = {,, θ} in (), (3) and (8). The gravit g in () can be absorbed into the pressure term. The Lagrange multiplier λ defined over B in () can be viewed as an etra bod force to maintain the rigid bod motion inside B. The conformation tensor C inside the rigid particle is etended as the identit tensor I as in (4) since the polmeric stress tensor σ p is zero inside the rigid particle. (UH) Ma 3, 6 / 58
12 Lie s scheme Consider the following initial value problem: dφ dt + A(φ) = on (, T ), φ() = φ (9) with < T < +. We suppose that operator A has a decomposition such as A = J j= A j with J. (UH) Ma 3, 6 / 58
13 Lie s scheme Let τ(> ) be a time-discretization step, we denote nτ b t n. With φ n denoting an approimation of φ(t n ), the Lie s scheme reads as follows: For n, assuming that φ n is known (with φ = φ ), compute φ n+ via dφ dt + A j(φ) = on (t n, t n+ ), φ(t n ) = φ n+(j )/J ; φ n+j/j = φ(t n+ ), () for j =,..., J. The Lie s scheme is first order accurate, but its low order of accurac is compensated b its simplicit, making it (relativel) eas to implement, and b its robustness. Chorin, Hughes, McCracken, Marsden. Comm. Pure Appl. Math. 3 (978) (UH) Ma 3, 6 3 / 58
14 Lie s scheme The Lie s operator scheme allows us to decouple the following difficulties: The incompressibilit condition, and the related unknown pressure. The advection terms. The rigid-bod motion in B h (t), and the related DLM λ h. (UH) Ma 3, 6 4 / 58
15 Lie s scheme The constitutive equation satisfied b the conformation tensor C is split with J = 3 for now to show how the factorization approach works out. Suppose that C n and u are known, we compute dc dt + (u )C = on (tn, t n+ ), () C(t n ) = C n ; C n+/3 = C(t n+ ), dc dt ( u)c C( u)t + f (Cn+/3 ) C = on (t n, t n+ ), λ () C(t n ) = C n+/3 ; C n+/3 = C(t n+ ), dc dt = f (Cn+/3 ) I on (t n, t n+ ), λ (3) C(t n ) = C n+/3 ; C n+ = C(t n+ ). (UH) Ma 3, 6 5 / 58
16 Lie s scheme We have derived the following two equivalent equations based on the factorization approach with s = f (C n+/3 ) for equations () and (): Lemma For a matri A and C = AA t, given the velocit u, λ (> ) and a constant s, (a). if A satisfies the equation da + (u )A =, then C satisfies dt the equation dc + (u )C = ; dt (b). if A satisfies the equation da dt + s A ( u)a =, then C λ satisfies the equation dc dt + s λ C ( u)c C( u) t =. Lozinski & Owens. J. Non-Newtonian Fluid Mech. (3). (UH) Ma 3, 6 6 / 58
17 FE approimation Fig. 3: Eamples of uniform finite element mesh for the flow region in IR (left) and IR 3 (right). Ω IR d (d = or 3). h is a space discretization step. T h is a uniform finite element mesh of Ω (See Fig. 3). T h is a uniform finite element mesh twice coarser than T h. (UH) Ma 3, 6 7 / 58
18 FE approimation The finite dimensional spaces for approimating V g (t), (H (Ω))d, L (Ω), L (Ω), V C L (t), and V C, respectivel, are defined b V gh (t) ={v h v h (C (Ω)) d, v h E (P ) d, E T h, v h Γ = g h (t)}, V h ={v h v h (C (Ω)) d, v h E (P ) d, E T h, v h Γ = }, L h ={q h q h C (Ω), q h E P, E T h }, L h ={q h q h L h, q h d = }, Ω V CLh (t) ={s h s h (C (Ω)) d d, s h E (P ) d d, E T h, s h Γ =C Lh (t)}, h V Ch ={s h s h (C (Ω)) d d, s h E (P ) d d, E T h, s h Γ = }, h where P is the space of the polnomials in d variables of degree, g h (t) is an approimation of g satisfing Γ g h(t) ndγ =, and Γ h = { Γ, g h(, t) n() < }. (UH) Ma 3, 6 8 / 58
19 Discrete Lagrange multiplier Choose { i } K i= as a set of points that covers B(t) evenl, and then define the discrete Lagrange multiplier space b K Λ h (t)={µ µ = µ j δ( j ), µ j IR d, j =,..., K}, (4) j= where δ( j ) is the Dirac measure at j. Fig. 4: An eample of collocation points on the particle surface. (UH) Ma 3, 6 9 / 58
20 Discrete Lagrange multiplier Then instead of the scalar product of (H (B h (t)) d, we shall use < µ, v > Bh (t) defined b <µ, v> Bh (t) = K µ j v( j ), µ Λ h (t), v V gh (t) or V h. (5) j= Using the above scalar product implies that the rigid bod motion of B(t) is forced via a collocation. (UH) Ma 3, 6 / 58
21 Appling the Lie s scheme to the discrete analogue of the problem ()-(8) with the C = AA t (and then equations for A) and the backward Euler s to some subproblems, we obtain u = u h, C = C h, G = G, V = V, ω = ω given, (6) for n, u n, C n, G n, V n, ω n being known, we compute u n+ 5, and p n+ 5 via the solution of u n+ 5 u n ρ f vd p n+ 5 vd=, v Vh, Ω t Ω (7) q u n+ 5 d=, q L h ; u n+ 5 V n+ g h, p n+ 5 L h. Ω (UH) Ma 3, 6 / 58
22 Net, we compute u n+ 5 and A n+ 5 via the solution of du(t) ρ f vd+ (u n+ 5 )u(t) vd=, v V n+, Ω dt h ; Ω u(t n ) = u n+ 5, u(t) V h, u(t) = g h (t n+ ) on Γ n+, [t n, t n+ ]; da(t) :sd+ (u n+ 5 )A(t):sd=, s VAh ; Ω dt Ω A(t n ) = A n, where A n (A n ) t = C n A(t) V n+ A Lh, t [t n, t n+ ]; and set u n+ 5 = u(t n+ ) and A n+ 5 = A(t n+ ), where Γ n+, = { Γ, g h (t n+ )() n() < }, V h = {v h v h (C (Ω)) d, v h E (P ) d, E T h, }, and V n+, h = {v V h, v = on Γ n+, }. (8) (9) (UH) Ma 3, 6 / 58
23 Then, compute u n+ 3 5 and A n+ 3 5 via the solution of and set u n+ 3 5 u n+ 5 ρ f vd+αµ u n+ 3 5 : vd=, Ω t Ω v V h ; u n+ 3 5 V n+ Ω ( An+ 3 5 A n+ 5 t g h, ( u n+ 3 5 )A n f (An+ 5 (A n+ 5 ) t ) A n+ 3 5 ) : sd =, λ s V Ah ; A n+ 3 5 V n+ A Lh, (3) (3) C n+ 3 5 = A n+ 3 5 (A n+ 3 5 ) t + t f (An+ 3 5 (A n+ 3 5 ) t ) λ I. (3) (UH) Ma 3, 6 3 / 58
24 Then, predict the position and the translation velocit of the center of mass as follows: Take V n+ 3 5, = V n and G n+ 3 5, = G n ; then predict the new position and translation velocit via the following sub-ccling and predicting-correcting technique For k =,,..., N, compute ˆV n+ 3 5,k = V n+ 3 5,k + ( ρ f /ρ s ) M p F r (G n+ 3 5,k ) t/n, (33) Ĝ n+ 3 5,k = G n+ 3 5,k + ( t/4n)(ˆv n+ 3 5,k + V n+ 3 5,k ), (34) V n+ 3 5,k = V n+ 3 5,k + ( ρ f /ρ s ) M p (F r (Ĝn+ 3 5,k ) (35) +F r (G n+ 3 5,k )) t/4n, G n+ 3 5,k = G n+ 3 5,k + ( t/4n)(v n+ 3 5,k + V n+ 3 5,k ), (36) end do; let V n+ 3 5 = V n+ 3 5,N, G n+ 3 5 = G n+ 3 5,N. (UH) Ma 3, 6 4 / 58
25 Net compute {u n+ 4 5, λ n+ 4 5, V n+ 4 5, ω n+ 4 5 } via the solution of u n+ 4 5 u n+ 3 5 ρ f vd + βµ u n+ 4 5 : vd Ω t [ Ω ] +( ρ f V n+ 4 5 V n+ 3 5 ω n+ 4 5 ω n ) M p Y + I p θ ρ s t t =< λ n+ 4 5, v Y θ G n+ 3 5 >B n+ 3 +( ρ 5 f /ρ s )M p g Y, h v V h, Y IR d, θ IR d, < µ, u n+ 4 5 V n+ 4 5 ω n+ 4 5 G n+ 3 5 > n+ 3 =, µ Λ n+ 3 5 B 5 h ; h u n+ 4 5 Vg n+ h, λ n+ 4 n+ 5 Λ 3 5 h, and set C n+ 4 5 = C n+ 3 5, and then let C n+ 4 5 = I in B n+ 3 5 h. (37) (UH) Ma 3, 6 5 / 58
26 Then take V n+, = V n+ 4 5 and G n+, = G n+ 3 5 ; and predict the final position and translation velocit as follows: For k =,,..., N, compute ˆV n+,k = V n+,k + ( ρ f /ρ s ) M F r (G n+,k ) t/n, (38) Ĝ n+,k = G n+,k + ( t/4n)(ˆv n+,k + V n+,k ), (39) V n+,k = V n+,k + ( ρ f /ρ s ) M (F r (Ĝn+,k ) (4) +F r (G n+,k )) t/4n, G n+,k = G n+,k + ( t/4n)(v n+,k + V n+,k ), (4) end do; let V n+ = V n+,n, G n+ = G n+,n. (UH) Ma 3, 6 6 / 58
27 Finall, compute u n+ via the solution of u n+ u n+ 4 5 ρ f vd + γµ u n+ : vd, Ω t Ω = η v ( f (C n+ 4 5 )(C n+ 4 5 I))d, v Vh ; u n+ V n+ λ Ω g h. (4) We complete the final step b setting C n+ = C n+ 4 5, and ω n+ = ω n In the above, Vg n+ h = V gh (t n+ ), Λ n+s h = Λ h (t n+s ), V n+ A Lh = V ALh (t n+ ), B n+s h = B h (t n+s ) where the spaces V ALh (t) and V Ah for A are defined similar to those C Lh (t) and C h, and α + β + γ =, where α, β, γ. (UH) Ma 3, 6 7 / 58
28 Solution strategies In the algorithm (6)-(4), we have obtained a sequence of simpler subproblems via Lie s scheme: A degenerated quasi-stokes problem (7) is solved b an Uzawa/ preconditioned conjugate gradient algorithm operating in the space L h as in Glowinski: Handbook of Analsis (3), Glowinski, Pan & Periau: Comput. Appl. Mech. Engrg. (998). The advection problems (8) and (9) are solved b a wave-like equation (Dean & Glowinski: C.R. Acad. Sci. Paris, Série (997)) which is an eplicit and does not introduce numerical dissipation. Since the advection problem is decoupled from the other ones, we can choose a proper sub-time step so that the CFL condition is satisfied. Problem (37), concerning the rigid bod motion enforcement, is a saddle point problem and is solved b a conjugate gradient developed in Glowinski, Pan, Hesla & Joseph: Int. J. Multiphase Flow (999). (UH) Ma 3, 6 8 / 58
29 Solution strategies Problems (3) and (4) are classical elliptic problems which can be solved b a matri-free fast solver from FISHPAK b Adams et al.: NCAR (98). Here we have chosen the use of uniform mesh so that a matri-free fast solver can be used to solve classical elliptic problems arising in problems (3), (37) and (4), which is especiall important for three dimensional flow problems we have proposed here since the memor cost can be reduced a lot. In (33)-(36) and (38)-(4), we have used a predicting-correcting scheme to obtain the position of the mass center and the translation velocit of the particle. In order to let the short range repulsive force to be activated effectivel, we have, tpicall, chosen N = sub-time steps to move the particle in one time step in the simulation. (UH) Ma 3, 6 9 / 58
30 Solution strategies Problem (3) gives a simple equation at each grid point which can be solved easil if we use the trapezoidal quadrature rule to compute the integrals as in Hao, Pan, Glowinski & Joseph: J. Non-Newtonian Fluid Mech (9). The discretized conformation tensor C n is clearl smmetric, and positive definite, b construction. We have split the constitutive equation for C b the Lie s scheme and the factorization to obtain (9), (3) and (3). In Lozinski and Owens: J. Non-Newtonian Fluid Mech. (3), the have split the equation with J =. Here we have decoupled the advection term so that we can appl our wave-like equation to solve the advection problem for A in (9). To us, (9) and (3) are simpler subproblems. (UH) Ma 3, 6 3 / 58
31 in an Oldrod-B fluid The computational domain is Ω = (,)(,6) initiall, then it moves verticall with the center of the disk. The disk diameter is d =.5 and its densit is ρ s=.. The initial position of the disk center is (.35,.5). The fluid densit is ρ f = and its viscosit is η =.. The relaation time and retardation time are λ = and λ = λ /4, respectivel. The particle Renolds number is Re = ρ f Ud, the Debra η number is De = λ U, the Mack number is d M = DeRe, and the elasticit number is E = De/Re where U is the averaged terminal speed Fig. 5. Computational domain Ω and a disk at its initial position. (UH) Ma 3, 6 3 / 58
32 in an Oldrod-B fluid V h=/96, dt=.4 h=/96, dt=. h=/44, dt= t V h=/96, dt=.4 h=/96, dt=. h=/44, dt= t Fig. 6. Histories of the -component (left) and -component (right) of the translational velocit: Re =.8, M =.3, De =.57, E = 3.. Hao, Pan, Glowinski, Joseph. J. Non-Newtonian Fluid Mech. 56 (9). (UH) Ma 3, 6 3 / 58
33 in viscoelastic fluids The computational domain is Ω = (,)(,6) initiall, then it moves verticall with the center of the disk. The disk diameters are d =.5 and the densit is ρ s=.. The initial position of the disk centers are (.35,.5) and (.65,.5). The fluid densit is ρ f = and its viscosit is η =.. The relaation time λ is either.5 and ; and the retardation time λ is λ /4. The maimal etension L of the immersed polmer coils is 5,, and Fig. 7. Computational domain Ω and two disks at the initial position. (UH) Ma 3, 6 33 / 58
34 Two disk in viscoelastic fluids Oldrod-B FENE-CR: L=5 FENE-CR: L= FENE-CR: L= Oldrod-B FENE-CR: L=5 FENE-CR: L= FENE-CR: L=5.5 Fig. 8. Trajectories of two disks: (left) λ =.5 & E=.6 and (right) λ = & E=3.. (UH) Ma 3, 6 34 / 58
35 Two disk in viscoelastic fluids Horizontal velocit Horizontal velocit Oldrod-B FENE-CR: L=5 FENE-CR: L= FENE-CR: L= t Oldrod-B FENE-CR: L=5 FENE-CR: L= FENE-CR: L= t Fig. 9. Horizontal velocit of two disks: (top) λ =.5 & E=.6 and (bottom) λ = & E=3.. (UH) Ma 3, 6 35 / 58
36 Vertical velocit Vertical velocit Two disk in viscoelastic fluids Oldrod-B FENE-CR: L=5 FENE-CR: L= FENE-CR: L= t Oldrod-B FENE-CR: L=5 FENE-CR: L= FENE-CR: L= t Fig.. Vertical velocit of two disks: (top) λ =.5 & E=.6 (Oldrod-B: Re=., M=.55, De=.33; FENE-CR of L=5: Re=.97, M=.49, De=.36) and (bottom) λ = & E=3. (Oldrod-B: Re=.8, M=.37, De=.664; FENE-CR of L=5: Re=.95, M=.35, De=.65). (UH) Ma 3, 6 36 / 58
37 Two disk in viscoelastic fluids t= t=5 t=.5 t=5 t= t=5 t=3 t=35 t=4 t=45 t=5 t= t=6 t=65 t=7 t=75 t=8 t=85 t=9 t=95 t= t= t=5 t= t=5.5 t= t=5 t=3 t=35 t=4 t=45 t= t=55 t=6 t=65 t=7 t=75 t=8 t=85 t=9 t=95 t= t= t=5 t= t= t=5 t=5 t=3 t=35 t=4 t=45 t=5.5-5 t= t=6 t=65 t=7 t=75 t=8 t=85 t=9 t=95 t= t= t=5 t= t=5 t= t=5 t=3 t=35 t=4 t=45 t= t=55 t=6 t=65 t=7 t=75 t=8 t=85 t=9 t=95 t=.5 Fig.. Positions of two disks : (left two) FENE-CR of L =5, E =.6, (middle left two) Oldrod-B, E =.6, (middle right two) FENE-CR of L =5, E = 3., and (right two) Oldrod-B, E = 3.. (UH) Ma 3, 6 37 / 58
38 Particle in an Oldrod-B fluid Si disks form a long chain after drafting, kissing and chaining, which is different from the well known phenomenon called drafting, kissing and tumbling for disks settling in Newtonian fluid. The dimension of computational domain is Ω = (,)(,7) initiall, then it moves verticall with the center of the lowest disk. The disk diameters are d =.5 and the densit is ρ s=.. The fluid densit is ρ f = and its viscosit is η =.6. The relaation time and retardation time are λ =.3 and λ = λ /8, respectivel. The mesh sizes for the velocit field, etra stress tensor and pressure are h = /96, /96, and /48, respectivel; and the time step is.4. The averaged terminal speed is U =.535. The particle Renolds number is Re = ρ f Ud η =.476. The Debra number is De = λ U d =.798. The Mack number is M = DeRe =.343. The elasticit number is E = De/Re =5.48. (UH) Ma 3, 6 38 / 58
39 Particle in an Oldrod-B fluid Fig. : The snap shots of particle position and velocit field at t =,,, 4, 6, 8,, 4, 6, 8, and 3 (from left to right). (UH) Ma 3, 6 39 / 58
40 Particle in a FENE-CR fluid Si disks form a long chain after drafting, kissing and chaining, which is different from the well known phenomenon called drafting, kissing and tumbling for disks settling in Newtonian fluid. In the FENE-CR model, the polmer etension limit is set to be L = 5 in the function f (C) = L L/(L L tr(c)) where C is the conformation tensor. The viscoelastic fluid with a shorter polmer etension limit can not hold all s together. The dimension of computational domain is Ω = (,)(,7) initiall, then it moves verticall with the center of the lowest disk. The disk diameters are d =.5 and the densit is ρ s=.. The fluid densit is ρ f = and its viscosit is η =.6. The relaation time and retardation time are λ =.3 and λ = λ /8, respectivel. The mesh sizes for the velocit field, etra stress tensor and pressure are h = /96, /96, and /48, respectivel; and the time step is.4. The averaged terminal speed is U =.37. The particle Renolds number is Re = ρ f Ud η =.66. The Debra number is De = λ U d = The Mack number is M = DeRe =.944. The elasticit number is E = De/Re =5.48. (UH) Ma 3, 6 4 / 58
41 Particle in a FENE-CR fluid Fig. 3: The snap shots of particle position and velocit field at t =, 4, 6, 8,, 8,, 4, 6, and 3 (from left to right) in a FENE-CR fluid with L = 5. (UH) Ma 3, 6 4 / 58
42 Horizontal velocit Vertical velocit Comparison: case Oldrod-B FENE-CR: L=5 FENE-CR: L= FENE-CR: L= t Oldrod-B FENE-CR: L=5 FENE-CR: L= FENE-CR: L= t Oldrod-B FENE-CR: L=5 FENE-CR: L= FENE-CR: L=.5 Fig. 4: Histories of the particle velocit (left) and trajectories of si disks (right). (UH) Ma 3, 6 4 / 58
43 Comparison: case.5 t=..5 t=..5 t=.5 t=.5 t=4.5 t=4.5 t=6.5 t=6.5 t=8.5 t=8.5 t=.5 t= OL FENE OL FENE OL FENE OL FENE OL FENE OL FENE Fig. 5: Comparison of numerical results in Oldrod-B fluid (OL) and FENE-CR fluid (FENE) for L = 5: t =,, 4, 6, 8, and (from left to right). (UH) Ma 3, 6 43 / 58
44 Comparison: case t= t= t= t= t= t= t= t= t= t= t= t= OL FENE OL FENE OL FENE OL FENE OL FENE OL FENE Fig. 6: Comparison of numerical results in Oldrod-B fluid (OL) and FENE-CR fluid (FENE) for L = 5: t =,, 4, 6, 8, and (from left to right). (UH) Ma 3, 6 44 / 58
45 Comparison t= t= t= t= t= t= t= t= t= t= t=3 -.5 t= OL FENE OL FENE OL FENE OL FENE OL FENE OL FENE Fig. 7: Comparison of numerical results in Oldrod-B fluid (OL) and FENE-CR fluid (FENE) for L = 5: t =,, 4, 6, 8, and 3 (from left to right). (UH) Ma 3, 6 45 / 58
46 in viscoelastic fluids The computational domain is Ω = (,)(,5) initiall, then it moves verticall with the center of the disk. The disk diameters are d =.5 and the densit is ρ s=.5,.75, or.. The initial position of the disk centers are (.6, ), (.88, ) and (.6, ). The fluid densit is ρ f = and its viscosit is η =.6. The relaation time λ is either.75 or.5; and the retardation time λ is λ /8. The maimal etension L of the immersed polmer coils is 5,, and Fig. 8. Computational domain Ω and three disks at the initial position. (UH) Ma 3, 6 46 / 58
47 in viscoelastic fluids L=5 L= L=5 Oldrod-B L=5 L= L=5 Oldrod-B -7-8 L=5 L= L=5 Oldrod-B -7-8 Fig. 9: The disk trajectories: ρs =.5 (left),.75 (middle), and. (right) for E=3. (λ =.75). (UH) Ma 3, 6 47 / 58
48 in viscoelastic fluids L=5 L= L=5 Oldrod-B -.5 L=5 L= L=5 Oldrod-B L=5 L= L=5 Oldrod-B Fig. : The disk trajectories: ρs =.5 (left),.75 (middle), and. (right) for E=6.4 (λ =.5). (UH) Ma 3, 6 48 / 58
49 in an Oldrod-B fluid in a vertical clinder The dimension of computational domain is Ω = (,D)(,D)(,4D) initiall, then it moves verticall with the center of the ball. The diameter and height of the clinder inside Ω are D and 4D, respectivel. The ball diameter is d = and its densit is ρ s=3.58. The fluid densit is ρ f =.868 and its viscosit is η =3. The relaation time and retardation time are λ =. and λ = λ /9, respectivel. The ratio of the ball diameter and the clinder diameter is α = d/d=.,., or.5. The particle Renolds number is Re = ρ f Ud, the η Debra number is De = λ U, the Mack number is d M = DeRe, and the elasticit number is E = De/Re where U is the averaged terminal speed. Fig.. Computational domain Ω and a ball for the case α =.. (UH) Ma 3, 6 49 / 58
50 in an Oldrod-B fluid in a Verical velocit α =. vertical clinder Bodart & Crochet [994] Hu, Patankar & Zhu [] h=/8 & dt=.5 h=/8 & dt=.95 h=/6 & dt= t Fig.. The comparison of the vertical speed for the case α =.: U =5.54, Re =.8977, M =.8345, De =.7757, E =.864 for t =.5 and h = /6. Bodart & Crochet. J. non-newtonian Fluid Mech. 54 (994) Hu, Patankar & Zhu. J. Comput. Phs. 69 (). (UH) Ma 3, 6 5 / 58
51 in an Oldrod-B fluid in a Veritcal velocit α =. vertical clinder Bodart & Crochet [994] Hu, Patankar & Zhu [] h=/48 & dt=.45 h=/64 & dt=.5 h=/96 & dt= t Fig. 3. The comparison of the vertical speed for the case α =.: U =.9, Re =.754, M =.6557, De =.695, E =.864 for t =.9 and h = /96. Bodart & Crochet. J. non-newtonian Fluid Mech. 54 (994) Hu, Patankar & Zhu. J. Comput. Phs. 69 (). (UH) Ma 3, 6 5 / 58
52 in an Oldrod-B fluid in a Vertical velocit -5 - α =.5 vertical clinder Bodart & Crochet [994] h=/3 & dt=.5 h=/48 & dt=. h=/64 & dt= t Fig. 4. The comparison of the vertical speed for the case α =.5: U =3.9, Re =.64, M =.4, De =.956, E =.864 for t =.5 and h = /64. Bodart & Crochet. J. non-newtonian Fluid Mech. 54 (994) Hu, Patankar & Zhu. J. Comput. Phs. 69 (). (UH) Ma 3, 6 5 / 58
53 Particle in an Oldrod-B fluid Two balls form a vertical chain after drafting, kissing and chaining, which is different from the well known phenomenon called drafting, kissing and tumbling for disks settling in Newtonian fluid. The dimension of computational domain is Ω = (,)(,)(,4) initiall, then it moves verticall with the center of the lower ball. The ball diameters are all d =. and their densit is ρ s=.. The fluid densit is ρ f = and its viscosit is η =.6. The relaation time and retardation time are λ = and λ = λ /8, respectivel. The mesh sizes for the velocit field, etra stress tensor and pressure are h = /64, /64, and /3, respectivel; and the time step is.5. The averaged terminal speed is U =.797. The particle Renolds number is Re = ρ f Ud η =.546. The Debra number is De = λ U d =.546. The Mack number is M = DeRe =.546. The elasticit number is E = De/Re =. (UH) Ma 3, 6 53 / 58
54 t=.5 t=3 t= t=33 t=35 t=6 Particle in an Oldrod-B fluid Fig. 5: (Left) The particle position at t =.5, 3,, 33, 35, and 6 (from left to right and top to bottom), (middle) eperimental results, and (right) the trajectories of two balls (UH) Ma 3, 6 54 / 58
55 in an Oldrod-B fluid A long bod like ellipsoid turns its broadside to the direction in Oldrod-B fluid. The dimension of computational domain is Ω = (,)(,)(,4) initiall, then it moves verticall with the center of the lower ball. The semi-ais of the ellipsoid are.5,.65, and.65 and its densit is ρ s=.. The fluid densit is ρ f = and its viscosit is η =.6. The relaation time and retardation time are λ = and λ = λ /8, respectivel. The mesh sizes for the velocit field, etra stress tensor and pressure are h = /64, /64, and /3, respectivel; and the time step is.. The averaged terminal speed is U =.7. The particle Renolds number is Re = ρ f Ud η = The Debra number is De = λ U d =.48. The Mack number is M = DeRe =.4. The elasticit number is E = De/Re =4.6. (UH) Ma 3, 6 55 / 58
56 t= t=5 t=5 t=75 t=5 t=3 in an Oldrod-B fluid Fig. 6: (Left) The particle positions at t =, 5, 5, 75, 5, and 3 (from left to right and top to bottom), (middle) eperimental results, and (right) the trajector of the ellipsoid mass center (UH) Ma 3, 6 56 / 58
57 Some future works Chain of man balls. The effect of coil etension L on the cluster of particles. The orientation of long bod in viscoelastic fluid. (UH) Ma 3, 6 57 / 58
58 Thank ou for our attention! (UH) Ma 3, 6 58 / 58
3D DLM/FD METHODS FOR SIMULATING THE MOTION OF SPHERES IN BOUNDED SHEAR FLOWS OF OLDROYD-B FLUIDS
3D DLM/FD METHODS FOR SIMULATING THE MOTION OF SPHERES IN BOUNDED SHEAR FLOWS OF OLDROYD-B FLUIDS A Dissertation Presented to the Faculty of the Department of Mathematics University of Houston In Partial
More informationTsorng-Whay Pan. phone: (713) Web page: pan/
Tsorng-Whay Pan Department of Mathematics University of Houston Houston, TX 77204 e-mail: pan@math.uh.edu phone: (713) 743-3448 Web page: www.math.uh.edu/ pan/ Education: 1990 Ph. D., Mathematics University
More informationA formulation for fast computations of rigid particulate flows
Center for Turbulence Research Annual Research Briefs 2001 185 A formulation for fast computations of rigid particulate flows By N. A. Patankar 1. Introduction A formulation is presented for the direct
More informationJournal of Computational Physics 169, (2001) doi: /jcph , available online at on
Journal of Computational Physics 169, 363 426 (2001) doi:10.1006/cph.2000.6542, available online at http://www.idealibrary.com on A Fictitious Domain Approach to the Direct Numerical Simulation of Incompressible
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 informationDIRECT NUMERICAL SIMULATION OF LIQUID- SOLID FLOW
DIRECT NUMERICAL SIMULATION OF LIQUID- SOLID FLOW http://www.aem.umn.edu/solid-liquid_flows Sponsored by NSF-Grand Challenge Grant Fluid Mechanics & CFD Computer Scientists D.D. Joseph Y. Saad R. Glowinski
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 informationarxiv: v1 [physics.flu-dyn] 6 Jul 2017
arxiv:77.957v [physics.flu-dyn] 6 Jul 27 A 3D DLM/FD method for simulating the motion of spheres in a bounded shear flow of Oldroyd-B fluids Shang-Huan Chiu a, Tsorng-Whay Pan a,, Roland Glowinski a,b
More informationDiscrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction
Discrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction Problem Jörg-M. Sautter Mathematisches Institut, Universität Düsseldorf, Germany, sautter@am.uni-duesseldorf.de
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 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 informationFunctions of Several Variables
Chapter 1 Functions of Several Variables 1.1 Introduction A real valued function of n variables is a function f : R, where the domain is a subset of R n. So: for each ( 1,,..., n ) in, the value of f is
More informationCurved spacetime tells matter how to move
Curved spacetime tells matter how to move Continuous matter, stress energy tensor Perfect fluid: T 1st law of Thermodynamics Relativistic Euler equation Compare with Newton =( c 2 + + p)u u /c 2 + pg j
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 informationAn improved ADI-DQM based on Bernstein polynomial for solving two-dimensional convection-diffusion equations
Mathematical Theor and Modeling ISSN -50 (Paper) ISSN 5-05(Online) Vol 3, No.1, 013 An improved ADI-DQM based on Bernstein polnomial for solving two-dimensional convection-diffusion equations A.S.J. Al-
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 informationSurvey of Wave Types and Characteristics
Seminar: Vibrations and Structure-Borne Sound in Civil Engineering Theor and Applications Surve of Wave Tpes and Characteristics Xiuu Gao April 1 st, 2006 Abstract Mechanical waves are waves which propagate
More informationAn ellipsoidal particle in tube Poiseuille flow
J. Fluid Mech. (217), vol. 822, pp. 664 688. c Cambridge Universit Press 217 doi:1.117/jfm.217.298 An ellipsoidal particle in tube Poiseuille flow 664 Haibo Huang 1, and Xi-Yun Lu 1 1 Department of Modern
More informationUNSTEADY LOW REYNOLDS NUMBER FLOW PAST TWO ROTATING CIRCULAR CYLINDERS BY A VORTEX METHOD
Proceedings of the 3rd ASME/JSME Joint Fluids Engineering Conference Jul 8-23, 999, San Francisco, California FEDSM99-8 UNSTEADY LOW REYNOLDS NUMBER FLOW PAST TWO ROTATING CIRCULAR CYLINDERS BY A VORTEX
More informationMechanics Departmental Exam Last updated November 2013
Mechanics Departmental Eam Last updated November 213 1. Two satellites are moving about each other in circular orbits under the influence of their mutual gravitational attractions. The satellites have
More informationJoule Heating Effects on MHD Natural Convection Flows in Presence of Pressure Stress Work and Viscous Dissipation from a Horizontal Circular Cylinder
Journal of Applied Fluid Mechanics, Vol. 7, No., pp. 7-3, 04. Available online at www.jafmonline.net, ISSN 735-357, EISSN 735-3645. Joule Heating Effects on MHD Natural Convection Flows in Presence of
More informationChapter 6 2D Elements Plate Elements
Institute of Structural Engineering Page 1 Chapter 6 2D Elements Plate Elements Method of Finite Elements I Institute of Structural Engineering Page 2 Continuum Elements Plane Stress Plane Strain Toda
More informationENGI 4430 Advanced Calculus for Engineering Faculty of Engineering and Applied Science Problem Set 3 Solutions [Multiple Integration; Lines of Force]
ENGI 44 Advanced Calculus for Engineering Facult of Engineering and Applied Science Problem Set Solutions [Multiple Integration; Lines of Force]. Evaluate D da over the triangular region D that is bounded
More informationModified DLM method for finite-volume simulation of particle flow
Modified DLM method for finite-volume simulation of particle flow A. M. Ardekani, S. Dabiri, and R. H. Rangel Department of Mechanical and Aerospace Engineering, University of California, Irvine, CA 92697-3975,
More informationPotential flow of a second-order fluid over a tri-axial ellipsoid
Potential flow of a second-order fluid over a tri-axial ellipsoid F. Viana 1, T. Funada, D. D. Joseph 1, N. Tashiro & Y. Sonoda 1 Department of Aerospace Engineering and Mechanics, Universit of Minnesota,
More informationA Fluctuating Immersed Boundary Method for Brownian Suspensions of Rigid Particles
A Fluctuating Immersed Boundary Method for Brownian Suspensions of Rigid Particles Aleksandar Donev Courant Institute, New York University APS DFD Meeting San Francisco, CA Nov 23rd 2014 A. Donev (CIMS)
More informationLecture 7: Weak Duality
EE 227A: Conve Optimization and Applications February 7, 2012 Lecture 7: Weak Duality Lecturer: Laurent El Ghaoui 7.1 Lagrange Dual problem 7.1.1 Primal problem In this section, we consider a possibly
More informationDirect Simulations of 2D Fluid Particle Flows in Biperiodic Domains
Direct Simulations of 2D Fluid Particle Flows in Biperiodic Domains B. Maury Abstract We propose a method to simulate the motion of 2D rigid particles in a viscous, incompressible fluid. Within the Arbitrary
More informationDynamics. Basilio Bona. Semester 1, DAUIN Politecnico di Torino. B. Bona (DAUIN) Dynamics Semester 1, / 18
Dynamics Basilio Bona DAUIN Politecnico di Torino Semester 1, 2016-17 B. Bona (DAUIN) Dynamics Semester 1, 2016-17 1 / 18 Dynamics Dynamics studies the relations between the 3D space generalized forces
More information10.2 The Unit Circle: Cosine and Sine
0. The Unit Circle: Cosine and Sine 77 0. The Unit Circle: Cosine and Sine In Section 0.., we introduced circular motion and derived a formula which describes the linear velocit of an object moving on
More informationKrylov Integration Factor Method on Sparse Grids for High Spatial Dimension Convection Diffusion Equations
DOI.7/s95-6-6-7 Krlov Integration Factor Method on Sparse Grids for High Spatial Dimension Convection Diffusion Equations Dong Lu Yong-Tao Zhang Received: September 5 / Revised: 9 March 6 / Accepted: 8
More informationData Repository Hampel et al., page 1/5
GSA DATA REPOSITORY 2138 Data Repositor Hampel et al., page 1/5 SETUP OF THE FINITE-ELEMENT MODEL The finite-element models were created with the software ABAQUS and consist of a 1-km-thick lithosphere,
More informationDynamics and control of mechanical systems
JU 18/HL Dnamics and control of mechanical sstems Date Da 1 (3/5) 5/5 Da (7/5) Da 3 (9/5) Da 4 (11/5) Da 5 (14/5) Da 6 (16/5) Content Revie of the basics of mechanics. Kinematics of rigid bodies coordinate
More informationApproximation of fluid-structure interaction problems with Lagrange multiplier
Approximation of fluid-structure interaction problems with Lagrange multiplier Daniele Boffi Dipartimento di Matematica F. Casorati, Università di Pavia http://www-dimat.unipv.it/boffi May 30, 2016 Outline
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 informationMath 21a Homework 07 Solutions Spring, 2014
Math a Homework 7 Solutions Spring, 4. valuate the iterated integral. a) Stewart.7 # 6 ) e d d d We perform the iterated integral: e d d d e d d e d [ e [ ] 4 e + 4e. Note that we ve twice done an integral
More informationNONLINEAR DYNAMICS AND CHAOS. Numerical integration. Stability analysis
LECTURE 3: FLOWS NONLINEAR DYNAMICS AND CHAOS Patrick E McSharr Sstems Analsis, Modelling & Prediction Group www.eng.o.ac.uk/samp patrick@mcsharr.net Tel: +44 83 74 Numerical integration Stabilit analsis
More informationDirect Simulation of the Sedimentation of Elliptic Particles in Oldroyd-B Fluids
Direct Simulation of the Sedimentation of Elliptic Particles in Oldroyd-B Fluids by P. Y. Huang 1, H. H. Hu 2 and D. D. Joseph 1 1 Department of Aerospace Engineering and Mechanics and the Minnesota Supercomputer
More informationVector Fields. Field (II) Field (V)
Math 1a Vector Fields 1. Match the following vector fields to the pictures, below. Eplain our reasoning. (Notice that in some of the pictures all of the vectors have been uniforml scaled so that the picture
More informationComputational Fluid Dynamics (CFD, CHD)*
1 / 1 Computational Fluid Dnamics (CFD, CHD)* PDE (Shocks 1st); Part I: Basics, Part II: Vorticit Fields Rubin H Landau Sall Haerer, Producer-Director Based on A Surve of Computational Phsics b Landau,
More informationFluid Mechanics II. Newton s second law applied to a control volume
Fluid Mechanics II Stead flow momentum equation Newton s second law applied to a control volume Fluids, either in a static or dnamic motion state, impose forces on immersed bodies and confining boundaries.
More informationNumerical Analysis of Heat Transfer in the Unsteady Flow of a non- Newtonian Fluid over a Rotating cylinder
Bulletin of Environment, Pharmacolog and Life Sciences Bull. Env.Pharmacol. Life Sci., Vol 4 [Spl issue 1] 215: 318-323 214 Academ for Environment and Life Sciences, India Online ISSN 2277-188 Journal
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 informationAdditional Topics in Differential Equations
6 Additional Topics in Differential Equations 6. Eact First-Order Equations 6. Second-Order Homogeneous Linear Equations 6.3 Second-Order Nonhomogeneous Linear Equations 6.4 Series Solutions of Differential
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 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 informationA distributed Lagrange multiplier/ ctitious domain method for particulate ows
PERGAMN International Journal of Multiphase Flow 25 (1999) 755±794 www.elsevier.com/locate/ijmulflow A distributed Lagrange multiplier/ ctitious domain method for particulate ows R. Glowinski a, *, T.-W.
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 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 information1 Exercise: Linear, incompressible Stokes flow with FE
Figure 1: Pressure and velocity solution for a sinking, fluid slab impinging on viscosity contrast problem. 1 Exercise: Linear, incompressible Stokes flow with FE Reading Hughes (2000), sec. 4.2-4.4 Dabrowski
More informationFEM techniques for nonlinear fluids
FEM techniques for nonlinear fluids From non-isothermal, pressure and shear dependent viscosity models to viscoelastic flow A. Ouazzi, H. Damanik, S. Turek Institute of Applied Mathematics, LS III, TU
More informationUNCORRECTED. To recognise the rules of a number of common algebraic relations: y = x 1 y 2 = x
5A galler of graphs Objectives To recognise the rules of a number of common algebraic relations: = = = (rectangular hperbola) + = (circle). To be able to sketch the graphs of these relations. To be able
More informationDynamic Simulation of Shear-induced Particle Migration in a Two-dimensional Circular Couette Device *
Chin. J. Chem. Eng., 15(3) 333 338 (2007) Dynamic Simulation of Shear-induced Particle Migration in a Two-dimensional Circular Couette Device * YU Zhaosheng( 余钊圣 ) a,b, **, SHAO Xueming( 邵雪明 ) a and Roger
More informationMat 267 Engineering Calculus III Updated on 9/19/2010
Chapter 11 Partial Derivatives Section 11.1 Functions o Several Variables Deinition: A unction o two variables is a rule that assigns to each ordered pair o real numbers (, ) in a set D a unique real number
More informationLECTURE NOTES - VIII. Prof. Dr. Atıl BULU
LECTURE NOTES - VIII «LUID MECHNICS» Istanbul Technical Universit College of Civil Engineering Civil Engineering Department Hdraulics Division CHPTER 8 DIMENSIONL NLYSIS 8. INTRODUCTION Dimensional analsis
More informationMAE4700/5700 Finite Element Analysis for Mechanical and Aerospace Design
MAE4700/5700 Finite Element Analsis for Mechanical and Aerospace Design Cornell Universit, Fall 2009 Nicholas Zabaras Materials Process Design and Control Laborator Sible School of Mechanical and Aerospace
More informationTwo-Dimensional Formulation.
Two-Dimensional Formulation mi@seu.edu.cn Outline Introduction( 概论 ) Two vs. Three-Dimensional Problems Plane Strain( 平面应变 ) Plane Stress( 平面应力 ) Boundar Conditions( 边界条件 ) Correspondence between Plane
More informationExercise solutions: concepts from chapter 5
1) Stud the oöids depicted in Figure 1a and 1b. a) Assume that the thin sections of Figure 1 lie in a principal plane of the deformation. Measure and record the lengths and orientations of the principal
More informationTriple Integrals. y x
Triple Integrals. (a) If is an solid (in space), what does the triple integral dv represent? Wh? (b) Suppose the shape of a solid object is described b the solid, and f(,, ) gives the densit of the object
More informationand ( x, y) in a domain D R a unique real number denoted x y and b) = x y = {(, ) + 36} that is all points inside and on
Mat 7 Calculus III Updated on 10/4/07 Dr. Firoz Chapter 14 Partial Derivatives Section 14.1 Functions o Several Variables Deinition: A unction o two variables is a rule that assigns to each ordered pair
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 informationConsider a slender rod, fixed at one end and stretched, as illustrated in Fig ; the original position of the rod is shown dotted.
4.1 Strain If an object is placed on a table and then the table is moved, each material particle moves in space. The particles undergo a displacement. The particles have moved in space as a rigid bod.
More informationEMA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 4 Pure Bending
EA 3702 echanics & aterials Science (echanics of aterials) Chapter 4 Pure Bending Pure Bending Ch 2 Aial Loading & Parallel Loading: uniform normal stress and shearing stress distribution Ch 3 Torsion:
More information8. THEOREM If the partial derivatives f x. and f y exist near (a, b) and are continuous at (a, b), then f is differentiable at (a, b).
8. THEOREM I the partial derivatives and eist near (a b) and are continuous at (a b) then is dierentiable at (a b). For a dierentiable unction o two variables z= ( ) we deine the dierentials d and d to
More informationCoordinate geometry. + bx + c. Vertical asymptote. Sketch graphs of hyperbolas (including asymptotic behaviour) from the general
A Sketch graphs of = a m b n c where m = or and n = or B Reciprocal graphs C Graphs of circles and ellipses D Graphs of hperbolas E Partial fractions F Sketch graphs using partial fractions Coordinate
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More informationSimultaneous Orthogonal Rotations Angle
ELEKTROTEHNIŠKI VESTNIK 8(1-2): -11, 2011 ENGLISH EDITION Simultaneous Orthogonal Rotations Angle Sašo Tomažič 1, Sara Stančin 2 Facult of Electrical Engineering, Universit of Ljubljana 1 E-mail: saso.tomaic@fe.uni-lj.si
More information7. TURBULENCE SPRING 2019
7. TRBLENCE SPRING 2019 7.1 What is turbulence? 7.2 Momentum transfer in laminar and turbulent flow 7.3 Turbulence notation 7.4 Effect of turbulence on the mean flow 7.5 Turbulence generation and transport
More informationA numerical approximation with IP/SUPG algorithm for P-T-T viscoelastic flows
Available online at www.tjnsa.com J. Nonlinear Sci. Appl. 6 (2016, 152 161 Research Article A numerical approximation with IP/SUPG algorithm for P-T-T viscoelastic flows Lei Hou a, Yunqing Feng a,, Lin
More information2.2 SEPARABLE VARIABLES
44 CHAPTER FIRST-ORDER DIFFERENTIAL EQUATIONS 6 Consider the autonomous DE 6 Use our ideas from Problem 5 to find intervals on the -ais for which solution curves are concave up and intervals for which
More informationModelling Rayleigh Taylor instability of a sedimenting suspension of several thousand circular particles in a direct numerical simulation
J. Fluid Mech. (2001), vol. 434, pp. 23 37. Printed in the United Kingdom c 2001 Cambridge University Press 23 Modelling Rayleigh Taylor instability of a sedimenting suspension of several thousand circular
More informationCHAPTER 4 BOUNDARY LAYER FLOW APPLICATION TO EXTERNAL FLOW
CHAPTER 4 BOUNDARY LAYER FLOW APPLICATION TO EXTERNAL FLOW 4.1 Introduction Boundary layer concept (Prandtl 1904): Eliminate selected terms in the governing equations Two key questions (1) What are the
More informationMathematics 309 Conic sections and their applicationsn. Chapter 2. Quadric figures. ai,j x i x j + b i x i + c =0. 1. Coordinate changes
Mathematics 309 Conic sections and their applicationsn Chapter 2. Quadric figures In this chapter want to outline quickl how to decide what figure associated in 2D and 3D to quadratic equations look like.
More informationWave Phenomena Physics 15c
Wave Phenomena Phsics 15c Lecture 13 Multi-Dimensional Waves (H&L Chapter 7) Term Paper Topics! Have ou found a topic for the paper?! 2/3 of the class have, or have scheduled a meeting with me! If ou haven
More informationFully Discrete Energy Stable High Order Finite Difference Methods for Hyperbolic Problems in Deforming Domains: An Initial Investigation
Full Discrete Energ Stable High Order Finite Difference Methods for Hperbolic Problems in Deforming Domains: An Initial Investigation Samira Nikkar and Jan Nordström Abstract A time-dependent coordinate
More informationVector Penalty-Projection Methods for the Solution of Unsteady Incompressible Flows
Author manuscript, published in "5th International Symposium on Finite Volumes for Complex Applications, Aussois : France 2008)" Vector Penalty-Projection Methods for the Solution of Unsteady Incompressible
More informationStability and Geometric Conservation Laws
and and Geometric Conservation Laws Dipartimento di Matematica, Università di Pavia http://www-dimat.unipv.it/boffi Advanced Computational Methods for FSI May 3-7, 2006. Ibiza, Spain and Introduction to
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 informationVector and scalar penalty-projection methods
Numerical Flow Models for Controlled Fusion - April 2007 Vector and scalar penalty-projection methods for incompressible and variable density flows Philippe Angot Université de Provence, LATP - Marseille
More informationStudy of rotation of ellipsoidal particles in combined simple shear flow and magnetic fields
Study of rotation of ellipsoidal particles in combined simple shear flow and magnetic fields Jie Zhang 1, Cheng Wang 1 1. Department of Mechanical and Aerospace Engineering, Missouri University of Science
More informationSurface Integrals (6A)
Surface Integrals (6A) Surface Integral Stokes' Theorem Copright (c) 2012 Young W. Lim. Permission is granted to cop, distribute and/or modif this document under the terms of the GNU Free Documentation
More informationUNIT 4 HEAT TRANSFER BY CONVECTION
UNIT 4 HEAT TRANSFER BY CONVECTION 4.1 Introduction to convection 4. Convection boundar laers 4..1 Hdrodnamic boundar laer over flat plate 4.. Thermal boundar laer 4..3 Concentration boundar laer 4.3 Dimensional
More informationMATHEMATICS 200 December 2014 Final Exam Solutions
MATHEMATICS 2 December 214 Final Eam Solutions 1. Suppose that f,, z) is a function of three variables and let u 1 6 1, 1, 2 and v 1 3 1, 1, 1 and w 1 3 1, 1, 1. Suppose that at a point a, b, c), Find
More information10. The dimensional formula for c) 6% d) 7%
UNIT. One of the combinations from the fundamental phsical constants is hc G. The unit of this epression is a) kg b) m 3 c) s - d) m. If the error in the measurement of radius is %, then the error in the
More informationWave Phenomena Physics 15c
Wave Phenomena Phsics 15c Lecture 13 Multi-Dimensional Waves (H&L Chapter 7) Term Paper Topics! Have ou found a topic for the paper?! 2/3 of the class have, or have scheduled a meeting with me! If ou haven
More informationRobust Monolithic - Multigrid FEM Solver for Three Fields Formulation Rising from non-newtonian Flow Problems
Robust Monolithic - Multigrid FEM Solver for Three Fields Formulation Rising from non-newtonian Flow Problems M. Aaqib Afaq Institute for Applied Mathematics and Numerics (LSIII) TU Dortmund 13 July 2017
More informationModeling and Analysis of Dynamic Systems
Modeling and Analysis of Dynamic Systems by Dr. Guillaume Ducard Fall 2016 Institute for Dynamic Systems and Control ETH Zurich, Switzerland 1/21 Outline 1 Lecture 4: Modeling Tools for Mechanical Systems
More informationFragment Orbitals for Transition Metal Complexes
Fragment Orbitals for Transition Metal Complees Introduction so far we have concentrated on the MO diagrams for main group elements, and the first 4 lectures will link into our Main Group chemistr course
More informationMultivariable Calculus
Multivariable Calculus In thermodynamics, we will frequently deal with functions of more than one variable e.g., P PT, V, n, U UT, V, n, U UT, P, n U = energy n = # moles etensive variable: depends on
More informationPHYSICS PART II SECTION- I. Straight objective Type
PHYSICS PAT II SECTION- I Straight objective Tpe This section contains 9 multiple choice questions numbered to 1. Each question has choices,, (C) and, out of which ONLY ONE is correct.. A parallel plate
More informationGeneral Physics I. Lecture 3: Newton's Laws. Prof. WAN, Xin ( 万歆 )
General Physics I Lecture 3: Newton's Laws Prof. WAN, Xin ( 万歆 ) xinwan@zju.edu.cn http://zimp.zju.edu.cn/~xinwan/ What Causes Changes of Motion? We define the interaction of a body with its environment
More informationA modular, operator splitting scheme for fluid-structure interaction problems with thick structures
Numerical Analysis and Scientific Computing Preprint Seria A modular, operator splitting scheme for fluid-structure interaction problems with thick structures M. Bukač S. Čanić R. Glowinski B. Muha A.
More informationCitation Key Engineering Materials, ,
NASITE: Nagasaki Universit's Ac Title Author(s) Interference Analsis between Crack Plate b Bod Force Method Ino, Takuichiro; Ueno, Shohei; Saim Citation Ke Engineering Materials, 577-578, Issue Date 2014
More informationTowards Higher-Order Schemes for Compressible Flow
WDS'6 Proceedings of Contributed Papers, Part I, 5 4, 6. ISBN 8-867-84- MATFYZPRESS Towards Higher-Order Schemes for Compressible Flow K. Findejs Charles Universit, Facult of Mathematics and Phsics, Prague,
More informationExact Differential Equations. The general solution of the equation is f x, y C. If f has continuous second partials, then M y 2 f
APPENDIX C Additional Topics in Differential Equations APPENDIX C. Eact First-Order Equations Eact Differential Equations Integrating Factors Eact Differential Equations In Chapter 6, ou studied applications
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 informationMath 21a: Multivariable calculus. List of Worksheets. Harvard University, Spring 2009
Math 2a: Multivariable calculus Harvard Universit, Spring 2009 List of Worksheets Vectors and the Dot Product Cross Product and Triple Product Lines and Planes Functions and Graphs Quadric Surfaces Vector-Valued
More informationInternational Journal of Applied Mathematics and Physics, 3(2), July-December 2011, pp Global Research Publications, India
International Journal of Applied Mathematics and Phsics, 3(), Jul-December 0, pp. 55-67 Global Research Publications, India Effects of Chemical Reaction with Heat and Mass Transfer on Peristaltic Flow
More informationy=1/4 x x=4y y=x 3 x=y 1/3 Example: 3.1 (1/2, 1/8) (1/2, 1/8) Find the area in the positive quadrant bounded by y = 1 x and y = x3
Eample: 3.1 Find the area in the positive quadrant bounded b 1 and 3 4 First find the points of intersection of the two curves: clearl the curves intersect at (, ) and at 1 4 3 1, 1 8 Select a strip at
More informationA Modified Distributed Lagrange Multiplier/Fictitious Domain Method for Particulate Flows with Collisions
A Modified Distributed Lagrange Multiplier/Fictitious Domain Metod for Particulate Flows wit Collisions P. Sing Department of Mecanical Engineering New Jersey Institute of Tecnology University Heigts Newark,
More information