Topological Derivatives in Shape Optimization
|
|
- Jocelyn Ward
- 6 years ago
- Views:
Transcription
1 Topological Derivatives in Shape Optimization Jan Sokołowski Institut Élie Cartan 28 mai 2012
2 Shape optimization Well posedness of state equation - nonlinear problems, compressible Navier-Stokes equations Stability of solutions to state equation with respect to the shape - nonlinear problems, compressible Navier-Stokes equations Formula for the shape gradient of shape functional - nonlinear problems, compressible Navier-Stokes equations Asymptotic analysis and the formula for the topological derivative of shape functional - anisotropic elasticity, compressible Navier-Stokes equations Numerical methods for the shape optimization problems - anisotropic elasticity, compressible Navier-Stokes equations
3 Compressible Navier-Stokes Equations weak renormalized global solutions of compressible Navier-Stokes equations : shape optimization problems are well posed strong local solutions of compressible Navier-Stokes equations : work and drag shape functionals are differentiable P.I. Plotnikov and J.S., Compressible Navier-Stokes Equations. Theory and Shape Optimization. Monografie Matematyczne, Vol. 73, Springer, Birkhauser Basel, 470 pages, 2013, to appear in 2012.
4 Anisotropic Elasticity asymptotic analysis of solutions to PDE s in singularly perturbed domains, eg., inclusions and cavities in anisotropic elasticity : topological derivatives of energy functionals combine the descriptions from applied mathematics and mechanics closed formulae in continuous case, easy to implement numerically
5 Asymptotic Analysis vs Shape Sensitivity asymptotic analysis is complementary tool for the shape sensitivity analysis : the way to evaluate the polarization tensors classical asymptotic analysis : the influence of boundary singularities on the energy functionals (Russian school) the regularity required for the asymptotic analysis shape sensitivity analysis leads to the same result on the boundary under weaker regularity assumptions, since the shape gradient is a distribution topological derivatives in the interior are obtained as the singular limits of shape derivatives conclusion : the topological derivative on the boundary of the energy functionals for elliptic system can be recovered from the shape derivative
6 REFERENCES J.S. and A. Zochowski : SICON, T. Lewinski and J.S. : IJSS, 2003 (mechanics) S.A. Nazarov and J.S. : JMPA, 2003 (mathematics) monograph : A.A. Novotny and J.S. to be published by Springer (mathematics and mechanics)
7 Topological derivative singular and regular domain perturbations : asymptotic analysis elementary examples relation between shape gradient and the topological derivative on the boundary energy for scalar equation energy for elasticity system multiscale elasticity
8 Topological derivative Let us consider a bounded domain Ω R d, with d 2, which is subject to a nonsmooth perturbation confined in a small region ω ε ( x) = x + εω of size ε, as shown in fig. 1. Here, x is an arbitrary point of Ω and ω is a fixed domain of R d. We introduce a characteristic function x χ(x), x R d, associated to the unperturbed domain, namely χ = 1 Ω, such that (1) Ω = χ, R d where Ω is the Lebesgue measure of Ω. Then, we define a characteristic function associated to the topologically perturbed domain of the form x χ ε ( x;x), x R d. In the case of a perforation, for instance, χ ε ( x) = 1 Ω 1 ωε ( x) and the perforated domain is obtained as Ω ε = Ω \ ω ε.
9 Topological derivative We assume that a given shape functional ψ(χ ε ( x)), associated to the topologically perturbed domain, admits the following topological asymptotic expansion (2) ψ(χ ε ( x)) = ψ(χ) + f (ε)t ( x) + o(f (ε)), where ψ(χ) is the shape functional associated to the original (unperturbed) domain, f (ε) is a positive function such that f (ε) 0, when ε 0. The function x T ( x) is called the topological derivative of ψ at x. Therefore, this derivative can be seen as a first order correction of ψ(χ) to approximate ψ(χ ε ( x)). In fact, after rearranging (2) we have (3) ψ(χ ε ( x)) ψ(χ) f (ε) = T ( x) + o(f (ε)) f (ε). The limit passage ε 0 in the above expression leads to ψ(χ ε ( x)) ψ(χ) (4) T ( x) = lim. ε 0 f (ε)
10 Topological derivative Since we are dealing with singular perturbations associated to nucleation of holes, the shape functionals ψ(χ ε ( x)) and ψ(χ) are associated to topologically different domains. Therefore, the above limit is not trivial to be calculated. In particular, we need to perform an asymptotic analysis of the shape functional ψ(χ ε ( x)) with respect to the small parameter ε. The notation for the shape functional we are going to introduce is the following : (5) ψ(χ ε ) := J Ωε (u ε ) and ψ(χ ε ) := J χε (u ε ) for singular and regular perturbations, respectively, where u ε is solution to the perturbed boundary values problems. FIGURE: The topological derivative concept.
11 Topological derivative In order to fix these ideas, we present simple examples. The first one concerns singular domain perturbation. The second example is associated to regular perturbation of the integrand.
12 Exemple 1 Example 1 Let us consider a shape functional of the form (6) ψ(χ ε ( x)) := g(x), Ω ε where the function g : R d R is assumed to be continuous, χ ε ( x) = 1 Ω 1 ωε ( x) and Ω ε = Ω \ ω ε. Since ω ε 0 as ε 0, we have (7) ψ(χ) := g(x). We are looking for an asymptotic expansion of the form (2), namely ψ(χ ε ( x)) = g(x) + g(x) g(x) Ω ε ω ε ω ε = g(x) g(x) Ω ω ε (8) = ψ(χ) ω ε g( x) + o( ω ε ). For the simple example, the term g( x) is called the topological derivative of the shape functional ψ, that is, T ( x) = g( x), and we recognize f (ε) = ω ε ε d. Ω
13 Example 2 Example 2 Let us now consider a shape functional of the form (9) ψ(χ ε ( x)) := g ε (x), Ω where χ ε ( x) = 1 Ω (1 γ)1 ωε ( x) and g ε = χ ε g is defined as { g(x) if x Ω \ ωε ( x) (10) g ε (x) := γg(x) if x ω ε ( x), with g standing by a continuous function of the form g : R d R. It means that the set Ω includes the interface ω ε. This is the case of shape functionals depending on the characteristic function of small sets. Since ω ε 0 as ε 0, we have (11) ψ(χ) := g(x). Ω
14 Example 2 We are looking for an asymptotic expansion of the form (2), namely ψ(χ ε ( x)) = g(x) + γg(x) Ω\ω ε ω ε = g(x) + γg(x) + g(x) g(x) Ω\ω ε ω ε ω ε ω ε = g(x) (1 γ) g(x) Ω ω ε (12) = ψ(χ) ω ε (1 γ)g( x) + o( ω ε ). For the simple example, the term (1 γ)g( x) is called the topological derivative of the shape functional ψ, that is, T ( x) = (1 γ)g( x), and we recognize f (ε) = ω ε ε d. In addition, we note that for the limit case γ 0, we have T ( x) = g( x). It means that the first example can be seen as the singular limit of this last one.
15 Relation Between Shape and Topological Changes The topological derivative can be considered as the singular limit of the shape derivative, so it is a generalization of the classical tool in shape optimization. In order to fix this idea, let us present an example concerning the relation between shape and topological derivatives.
16 Relation Between Shape and Topological Changes Example 3 Let as consider a bar under torsion effects. The cross section of the shaft is represented by an open bounded domain Ω R 2, as shown in fig. 2. Following the Prandtl s approach, we consider the complementary dissipation energy shape functional of the form (13) ψ(χ) = 1 u 2 bu, 2 Ω where the scalar field u is solution to the following boundary value problem Find u, such that (14) u = b in Ω, u = 0 on Ω. with the constant b used to denote a rigid twist of the cross-section of the shaft. The shape derivative of ψ(χ) is given by (15) dψ(χ;v) = 1 ( n u) 2 n V. 2 Ω where V is the shape change velocity field Ω
17 Shape and topological derivatives The topological asymptotic expansion for the above problem is given by (16) ψ(χ ε ( x)) = ψ(χ) + πε 2 u( x) 2 + o(ε 2 ) x Ω. where we recognize the topological derivative as T ( x) = u( x) 2. This last result can be extend to the boundary by using the odd extension of the problem, as shown in fig. 2. In addition, since u = 0 on Ω, then its tangential derivative vanishes on the boundary, namely τ u = 0 on Ω. Therefore, the topological asymptotic expansion on the boundary leads to (17) ψ(χ ε ( x)) = ψ(χ) πε2 ( n u( x)) 2 + o(ε 2 ) x Ω.
18 Shape and topological derivatives It means that in this simple example, the topological derivative T ( x) = 1 2 ( nu( x)) 2 coincides with the density of the shape derivative at the point x Ω. FIGURE: Cross section of a shaft under torsion effects : original and singulary perturbed domains.
19 Energy Change-Poisson problem It is sufficient to define V (18) { V = 0 on Ω V = n on B ε, where n = (x x)/ε, with x B ε, is the normal unit vector field pointing toward the center of the circular hole B ε. FIGURE: Topologically perturbed domain by the nucleation of a small circular hole.
20 Energy Change The shape functional in the unperturbed domain Ω is given by (19) ψ(χ) := J Ω (u) = 1 2 Ω u 2 bu + qu, Ω Γ N where the scalar function u is the solution to the variational problem : Find u U, such that (20) u η = bη qη η V, Γ N Ω In the above equation, b is a source-term assumed to be constant everywhere. The set U and the space V are respectively defined as Ω (21) (22) U := {ϕ H 1 (Ω) : ϕ ΓD = u}, V := {ϕ H 1 (Ω) : ϕ ΓD = 0}.
21 Energy Change In addition, Ω = Γ D Γ N with Γ D Γ N =, where Γ D and Γ N are Dirichlet and Neumann boundaries, respectively. Thus u is a Dirichlet data on Γ D and q is a Neumann data on Γ N, both functions are assumed to be smooth enough. FIGURE: The Poisson equation defined in the unperturbed domain Ω.
22 Energy Change The strong formulation associated to the variational problem (20) leads to the so-called Poisson equation, namely : Find u, such that u = b in Ω, (23) u = u on Γ D, n u = q on Γ N,
23 Energy Change Now, let us state the same problem in the perturbed domain Ω ε. In this case, the total potential energy reads (24) ψ(χ ε ) := J Ωε (u ε ) = 1 u ε 2 bu ε + qu ε, 2 Ω ε Ω ε Γ N where the scalar function u ε solves the variational problem : Find u ε U ε, such that (25) u ε η = bη qη η V ε. Ω ε Ω ε Γ N
24 Energy Change The set U ε and the space V ε must be defined according to the boundary condition on B ε. In particular, we have (26) (27) U ε := {ϕ H 1 (Ω ε ) : ϕ ΓD = u, βϕ Bε = 0}, V ε := {ϕ H 1 (Ω ε ) : ϕ ΓD = 0, βϕ Bε = 0}, with β {0,1}. This notation have to be interpreted as follows : When β = 1, we have homogeneous Dirichlet boundary condition on B ε, since u ε = 0 and η = 0 on B ε. When β = 0, u ε and η are free on the boundary of the hole B ε, then we have homogeneous Neumann boundary condition on B ε.
25 Energy Change FIGURE: The Poisson equation defined in the perturbed domain Ω ε.
26 Energy Change The strong formulation associated to the variational problem (25) reads Find u ε, such that u ε = b in Ω ε, (28) u ε = u on Γ D, n u ε = q on Γ N, βu ε + (1 β) n u ε = 0 on B ε.
27 Shape sensitivity analysis. Proposition Let J Ωε (u ε ) be the shape functional defined by (24). Then, the derivative of this functional with respect to the small parameter ε is given by (29) J Ωε (u ε ) = Σ ε V, Ω ε where V is the shape change velocity field defined through (18) and Σ ε is a generalization of the classical Eshelby energy-momentum tensor given by (30) Σ ε = 1 2 ( u ε 2 2bu ε )I u ε u ε.
28 Shape sensitivity analysis. Proof By making use of Reynolds transport theorem, the shape derivative (the material or total derivative with respect to the parameter ε) of the functional (24) is given by (31) J Ωε (u ε ) = 1 (( u ε 2 ) + u ε 2 divv) 2 Ω ε b( u ε + u ε divv) + q u ε, Ω ε Γ N since Γ N remains fixed according to (18). Next, by using the concept of material derivative of spatial fields we find that the first term of the above right-hand side integral can be written as ( u ε 2 ) = ( u ε u ε ) (32) = 2 u ε u ε 2 V u ε u ε = 2 u ε u ε 2( u ε u ε ) V,
29 Proof From these last results we obtain J Ωε (u ε ) = 1 (( u ε 2 2bu ε )I 2 u ε u ε ) V 2 Ω ε (33) + u ε u ε b u ε + q u ε. Ω ε Ω ε Γ N where we have made use of the identity divv = I V. Since u ε is a variation of u ε in the direction of the velocity field V, then u ε V ε. Finally, by taking u ε as test function in the variational problem (25), we have that the last three terms of the above equation vanish. Now, we can prove that the shape derivative of the functional J Ωε (u ε ) can be written in terms of quantities concentrated on the boundary Ω ε. In fact, the following result holds :
30 Boundary shape gradient Let J Ωε (u ε ) be the shape functional defined by (24). Then, the derivative of this functional with respect to the small parameter ε is given by (34) J Ωε (u ε ) = Σ ε n V, Ω ε where V is the shape change velocity field defined through (18) and tensor Σ ε is given by (30). PROOF We have (35) J Ωε (u ε ) = (36) Σ ε n V divσ ε V. Ω ε Ω ε divσ ε = 1 2 ( u ε u ε 2bu ε ) div( u ε u ε ) = ( u ε ) u ε b u ε u ε div u ε ( u ε ) u ε = ( u ε ) u ε ( u ε ) u ε ( u ε + b) u ε. Thus, since u ε = ( u ε ) and u ε = b, we observe that (37) divσ ε V = 0 V, Ω ε
31 Corollary We have (38) J Ωε (u ε ) = = Σ ε n V Ω ε Σ ε n V + B ε Ω Σ ε n V. Since we are dealing with an uniform expansion of the circular hole, then by taking into account the associated velocity field defined through (18), namely, V = n on B ε and V = 0 on Ω, we finally obtain (39) d dε ψ(χ ε) = J Ωε (u ε ) = Σ ε n n. B ε
32 How to evaluate topological derivative? Let ψ(χ ε ( x)) be the shape functional associated to the topologically perturbed domain, which admits, for ε small enough, the topological asymptotic expansion of the form (40) ψ(χ ε ( x)) = ψ(χ) + f (ε)t ( x) + R(f (ε)), where ψ(χ) is the shape functional associated to the original (unperturbed) domain, the positive function f (ε) is such that f (ε) 0, with ε 0, and the function T ( x) is the topological derivative of the shape functional ψ. We assume that the remainder R(f (ε)) = o(f (ε)) has the following additional property R (f (ε)) 0, when ε 0. Then, the topological derivative can be written as 1 (41) T ( x) = lim ε 0 f (ε) d dε ψ(χ ε( x)), where d dε ψ(χ ε( x)) is the derivative of ψ(χ ε ( x)) with respect to the small parameter ε.
33 Proof Let us calculate the total derivative of the expansion (40) with respect to the real parameter ε, that is (42) After division by f (ε) we have d dε ψ(χ ε( x)) = f (ε)t ( x) + R (f (ε))f (ε). (43) T ( x) = 1 d f (ε) dε ψ(χ ε( x)) + R (f (ε)). Finally, by taking the limit ε 0 in the above expression, the result holds since R (f (ε)) 0, when ε 0
34 Asymptotics of solutions. Neumann case For the Neumann case β = 0, let us start with the ansatz for an expansion of u ε in the form u ε (x) = u(x) + εw(ε 1 x) + ũ ε (x) = u( x) + u( x) (x x) + 1 u(y)(x x) (x x) 2 (44) + εw(ε 1 x) + ũ ε (x), where y is an intermediate point between x and x. On the boundary of the hole B ε we have n u ε Bε = 0. Thus, the normal derivative of the above expansion, evaluated on B ε, leads to (45) u( x) n ε u(y)n n + ε n w(ε 1 x) + n ũ ε (x) = 0.
35 Thus, we can choose w such that (46) ε n w(ε 1 x) = u( x) n on B ε. In the new variable ξ = ε 1 x, which implies ξ w(ξ ) = ε w(ε 1 x), the following exterior problem is considered, and formally obtained with ε 0 : Find w, such that (47) ξ w = 0 in R 2 \ B 1, w 0 at, ξ w n = u( x) n on B 1. The above boundary value problem admits an explicit solution, namely (48) w(ε 1 x) = ε u( x) (x x). x x 2
36 Now we can construct ũ ε in such a way that it compensates the discrepancies introduced by the higher-order terms in ε as well as by the boundary-layer w on the exterior boundary Ω. It means that the remainder ũ ε must be the solution to the following boundary value problem : (49) Find ũ ε, such that ũ ε = 0 in Ω ε, ũ ε = εw on Γ D, n ũ ε = ε n w on Γ N, n ũ ε = ε u(y)n n on B ε, where, under appropriate regularity assumption, clearly ũ ε = O(ε). However, this estimate can be improved, namely ũ ε = O(ε 2 ), and the expansion for u ε reads (50) u ε (x) = u(x) + ε2 x x 2 u( x) (x x) + O(ε2 ).
37 FIGURE: Polar coordinate system (r,θ) centered at the point x Ω.
38 Topological derivative evaluation From an orthonormal curvilinear coordinate system n and τ on the boundary B ε, the gradient u ε can be decomposed into its normal and tangential components, that is (51) u ε Bε = ( n u ε )n + ( τ u ε )τ. Therefore, on the boundary of the hole B ε we observe that Σ ε n n Bε = 1 2 ( u ε 2 2bu ε )In n ( u ε u ε )n n = 1 2 ( u ε 2 2bu ε ) ( u ε n) 2 = 1 2 (( nu ε ) 2 + ( τ u ε ) 2 2bu ε ) ( n u ε ) 2 (52) = 1 2 (( τu ε ) 2 ( n u ε ) 2 2bu ε ). Now, we shall analyze the behavior of the total potential energy with respect to the parameter ε by taking into account each type of boundary condition on B ε, namely, Neumann.
39 PROOF The shape derivative of the cost functional (39) reads (53) d dε ψ(χ ε) = J Ωε (u ε ) = 1 (( τ u ε ) 2 2bu ε ), 2 B ε where we have used (52). In addition, we have that the following expansion for u ε holds in the neighborhood of the hole (50) u ε (x) = u(x) + ε2 u( x) (x x) + o(ε) x x 2 (54) = u( x) + u( x) (x x) + ε2 u( x) (x x) + o(ε), x x 2 which can be written in a polar coordinate system (r,θ) centered at the point x (see fig. 6), namely (55) u ε (r,θ) = ϕ 0 + r2 + ε 2 (ϕ 1 cosθ + ϕ 2 sinθ) + o(ε), r where we have introduced the notation u( x) = ϕ 0 and u( x) = (ϕ 1,ϕ 2 ).
40 The tangential derivative of u ε is obtained in the following way τ u ε (x) = 1 r θ u ε (r,θ) (56) = r2 + ε 2 r 2 (ϕ 1 sinθ ϕ 2 cosθ) + o(ε), τ u ε (x) Bε = 2(ϕ 1 sinθ ϕ 2 cosθ) + o(ε), where the last expression was obtained by taking r = ε in the previous one. In addition, considering the above expansion in (53) and after explicitly evaluating the integral on the boundary of the hole B ε, we obtain d dε ψ(χ ε) = 1 2π (4(ϕ 1 sinθ ϕ 2 cosθ) 2 2bϕ 0 )εdθ + o(ε 2 ) 2 0 (57) = 2πε((ϕ1 2 + ϕ2 2 ) bϕ 0) + o(ε 2 ), Therefore, the above result together with the relation between shape and topological derivative given by (41) leads to 1 ( ) (58) T = lim ε 0 f 2πε((ϕ1 2 (ε) + ϕ2 2 ) bϕ 0) + o(ε 2 ).
41 Now, in order to identify the main term of the above expansion, we choose (59) f (ε) = πε 2, which leads to (60) T = (ϕ ϕ2 2 ) + bϕ 0. Recalling that u( x) = ϕ 0 and u( x) = (ϕ 1,ϕ 2 ), the final formula for the topological derivative becomes (61) T ( x) = u( x) 2 + bu( x). This leads to the topological asymptotic expansion of the energy shape functional, namely (62) ψ(χ ε ( x)) = ψ(χ) πε 2 ( u( x) 2 bu( x)) + o(ε 2 ).
42 The Navier Problem The shape functional in the unperturbed domain Ω is given by (63) ψ(χ) := J Ω (u) = 1 2 Ω σ(u) u s b u q u, Ω Γ N where the vector function u is the solution to the variational problem : (64) Find u U, such that σ(u) η s = b η + Ω Ω q η Γ N η V, with σ(u) = C u s.
43 In the above equation, C is the constitutive tensor given by (65) C = E ((1 ν)i + νi I), 1 ν2 where I and I are the second and fourth order identity tensors, respectively, E is the Young modulus and ν the Poisson ratio, both considered constants everywhere. For the sake of simplicity, we also assume that the thickness of the elastic body is constant and equal to one. The set U and the space V are respectively defined as (66) (67) U := {ϕ H 1 (Ω;R 2 ) : ϕ ΓD = u}, V := {ϕ H 1 (Ω;R 2 ) : ϕ ΓD = 0}. In addition, Ω = Γ D Γ N with Γ D Γ N =, where Γ D and Γ N are Dirichlet and Neumann boundaries, respectively. Thus u is a Dirichlet data on Γ D and q is a Neumann data on Γ N, both assumed to be smooth enough.
44 FIGURE: The Navier system defined in the unperturbed domain Ω.
45 The strong system associated to the variational problem (64) reads : Find u, such that divσ(u) = b in Ω, (68) σ(u) = C u s u = u on Γ D, σ(u)n = q on Γ N. Since the Young modulus E and the Poisson ratio ν are considered constants, the above boundary value problem reduces itself to the well-known Navier equations, namely (69) µ u (λ + µ) (divu) = b in Ω, with the Lame s coefficients µ and λ respectively given by (70) µ = E 2(1 + ν) and λ = νe 1 ν 2.
46 Now, let us state the same problem in the perturbed domain Ω ε. In this case, the total potential energy reads (71) ψ(χ ε ) := J Ωε (u ε ) = 1 σ(u ε ) u s ε b u ε q u ε, 2 Ω ε Ω ε Γ N where the vector function u ε solves the variational problem : (72) Find u ε U ε, such that σ(u ε ) η s = b η + q η η V ε, Ω ε Ω ε Γ N with σ(u ε ) = C u s ε.
47 The set U ε and the space V ε are defined as (73) (74) U ε := {ϕ H 1 (Ω ε ;R 2 ) : ϕ ΓD = u}, V ε := {ϕ H 1 (Ω ε ;R 2 ) : ϕ ΓD = 0}. We have homogeneous Neumann condition on the boundary of the hole. FIGURE: The Navier system defined in the perturbed domain Ω ε.
48 The strong system associated to the variational problem (72) reads : Find u ε, such that divσ(u ε ) = b in Ω ε, σ(u (75) ε ) = C u s ε u ε = u on Γ D, σ(u ε )n = q on Γ N, σ(u ε )n = 0 on B ε.
49 Shape Sensitivity Analysis Let J Ωε (u ε ) be the shape functional defined by (71). Then, the derivative of this functional with respect to the small parameter ε is given by (76) J Ωε (u ε ) = Σ ε V, Ω ε where V is the shape change velocity field defined through (18) and Σ ε is a generalization of the classical Eshelby energy-momentum tensor given by (77) Σ ε = 1 2 (σ(u ε) u s ε 2b u ε )I u ε σ(u ε ).
50 PROOF By Reynolds transport theorem the shape derivative (the material or total derivative with respect to the parameter ε) of the functional (71) is given by J Ωε (u ε ) = 1 ((σ(u ε ) u s ε) + (σ(u ε ) u s 2 ε)divv) Ω ε (78) b ( u ε + u ε divv) q u ε. Ω ε Γ N Next, by using the concept of material derivative of spatial fields we find that the first term of the above right-hand side integral can be written as (79) (σ(u ε ) u s ε) = 2σ(u ε ) u s ε 2σ(u ε ) ( u ε V) s = 2σ(u ε ) u s ε 2 u ε σ(u ε ) V, since σ(u ε ) = C u s ε. From this last result we obtain J Ωε (u ε ) = 1 ((σ(u ε ) u s ε 2b u ε )I 2 u ε σ(u ε )) V 2 Ω ε (80) + σ(u ε ) u s ε b u ε q u ε. Ω ε Ω ε Γ N
51 Boundary shape gradient Let J Ωε (u ε ) be the shape functional defined by (71). Then, the derivative of this functional with respect to the small parameter ε is given by (81) J Ωε (u ε ) = Σ ε n V, Ω ε with V standing by the shape change velocity field defined through (18) and tensor Σ ε is given by (77).
52 PROOF By making use of the other version of the Reynolds transport theorem the shape derivative of the functional (71) is (82) J Ωε (u ε ) = 1 (σ(u ε ) u s 2 ε) + 1 (σ(u ε ) u s Ω ε 2 ε)n V Ω ε b u ε (b u ε )n V q u ε. Ω ε Ω ε Γ N Next, by using the concept of shape derivatives of spatial fields, we find that the first term of the above right-hand side integral can be written as (83) (σ(u ε ) u s ε) = 2σ(u ε ) ( u ε) s, since σ(u ε ) = σ (u ε ) and ( ) is the partial derivative of ( ) with respect to ε. Now, let us use the relation between material and spatial derivatives of vector fields, namely ϕ = ϕ ( ϕ)v, which leads to (84) (σ(u ε ) u s ε) = 2σ(u ε ) ( u ε ( u ε )V) s = 2σ(u ε ) u s ε 2σ(u ε ) (( u ε )V) s.
53 From this last result we obtain J Ωε (u ε ) = 1 (σ(u ε ) u s ε 2b u ε )n V 2 Ω ε (σ(u ε ) (( u ε )V) s b ( u ε )V) Ω ε (85) + σ(u ε ) u s ε b u ε q u ε. Ω ε Ω ε Γ N Since u ε is a variation of u ε in the direction of the velocity field V, then u ε V ε. Now, by taking into account that u ε is the solution to the variational problem (72), we have that the last three terms of the above equation vanish. By using a tensor relation and after applying the divergence theorem, we have J Ωε (u ε ) = 1 (σ(u ε ) u s ε 2b u ε )n V ( u ε )V σ (u ε )n 2 Ω ε Ω ε (86) + (divσ(u ε ) + b) ( u ε )V. Ω ε
54 Then, since σ (u ε ) = σ(u ε ), we can rewrite the above equation in the compact form (87) J Ωε (u ε ) = Σ ε n V + Ω ε (divσ(u ε ) + b) ( u ε )V. Ω ε Considering that u ε is also solution to the strong system, namely divσ(u ε ) = b, the last term of the above equation vanishes, which leads to the required result. According to the obtained result we have (88) J Ωε (u ε ) = = Σ ε n V Ω ε Σ ε n V + B ε Ω Σ ε n V. Since we are dealing with an uniform expansion of the circular hole B ε, then in view of the associated velocity field (18), namely, V = n on B ε and V = 0 on Ω, we finally obtain d (89) dε ψ(χ ε) = J Ωε (u ε ) = Σ ε n n. B ε It means that the shape derivative of the cost functional can be represented by a boundary integral or distributions supported on B ε
55 Asymptotic Analysis FIGURE: Polar coordinate system (r,θ) centered at the point x Ω.
56 Let us start with the ansatz for an expansion of u ε in the form (90) u ε (x) = u(x) + w ε (x) + ũ ε (x). After applying the operator σ we have σ(u ε (x)) = σ(u(x)) + σ(w ε (x)) + σ(ũ ε (x)) (91) = σ(u( x)) + σ(u(y))(x x) + σ(w ε (x)) + σ(ũ ε (x)), where y is an intermediate point between x and x. On the boundary of the hole B ε we have σ(u ε )n Bε = 0. Thus, the normal projection of the above expansion, evaluated on B ε, leads to (92) σ(u( x))n ε( σ(u(y))n)n + σ(w ε (x))n + σ(ũ ε (x))n = 0.
57 Thus, we can choose σ(w ε ) such that (93) σ(w ε (x))n = σ(u( x))n on B ε. Now, the following exterior boundary value problem is considered and formally obtained with ε 0 : Find σ(w ε ), such that divσ(w (94) ε ) = 0 in R 2 \ B ε, σ(w ε ) 0 at, σ(w ε )n = σ(u( x))n on B ε.
58 The above boundary value problem admits an explicit solution (see, for instance, the book by Little 1973 ), which can be written in a polar coordinate system (r,θ) centered at the point x (see fig. 9) as σ rr ε 2 ) (95) (w ε (r,θ)) = ϕ 1 r 2 ϕ 2 (4 ε2 r 2 3ε4 r 4 cos2θ, (96) (97) σ θθ (w ε (r,θ)) = ϕ 1 ε 2 ε 4 r 2 3ϕ 2 r 4 cos2θ, ( σ rθ (w ε (r,θ)) = ϕ 2 2 ε2 r 2 3ε4 r 4 The coefficients ϕ 1 and ϕ 2 are given by ) sin2θ. (98) ϕ 1 = 1 2 (σ 1(u( x)) + σ 2 (u( x))), ϕ 2 = 1 2 (σ 1(u( x)) σ 2 (u( x))), where σ 1 (u( x)) and σ 2 (u( x)) are the eigenvalues of tensor σ(u( x)), which can be expressed as (99) σ 1,2 (u( x)) = 1 ( ) trσ(u( x)) ± 2σ 2 D (u( x)) σ D (u( x)), with σ D (u( x)) standing by the deviatoric part of the stress tensor σ(u( x)),
59 namely (100) σ D (u( x)) = σ(u( x)) 1 2 trσ(u( x))i. In addition, σ rr (ϕ), σ θθ (ϕ) and σ rθ (ϕ) are the components of tensor σ(ϕ) in the polar coordinate system, namely, σ rr (ϕ) = e r σ(ϕ)e r, σ θθ (ϕ) = e θ σ(ϕ)e θ and σ rθ (ϕ) = σ θr (ϕ) = e r σ(ϕ)e θ, with e r and e θ used to denote the canonical basis associated to the polar coordinate system (r,θ), such that, e r = e θ = 1 and e r e θ = 0. See fig. 9. Now we can construct ũ ε in such a way that it compensates the discrepancies introduced by the higher-order terms in ε as well as by the boundary-layer w ε on the exterior boundary Ω.
60 It means that the remainder ũ ε must be the solution to the following boundary value problem : Find ũ ε, such that divσ(ũ ε ) = 0 in Ω ε, σ(ũ (101) ε ) = C ũ s ε ũ ε = w ε on Γ D, σ(ũ ε )n = σ(w ε )n on Γ N, σ(ũ ε )n = εh on B ε. where ũ ε = O(ε 2 ). Finally, the expansion for σ(u ε ) in the polar coordinate system (r,θ) reads ( ) ) σ rr (u ε (r,θ)) = ϕ 1 1 ε2 r 2 + ϕ 2 (1 4 ε2 r 2 + 3ε4 r 4 cos2θ + O(ε 2 ), ( ) ) σ θθ (u ε (r,θ)) = ϕ ε2 r 2 ϕ 2 (1 + 3 ε4 r 4 cos2θ + O(ε 2 ), ( σ rθ (u ε (r,θ)) = ϕ ε2 r 2 3ε4 r 4 ) sin2θ + O(ε 2 ), where we have used the fact that σ rr (u( x)) = ϕ 1 + ϕ 2 cos2θ, σ θθ (u( x)) = ϕ 1 ϕ 2 cos2θ and σ rθ (u( x)) = ϕ 2 sin2θ.
61 Topological Derivative Evaluation From an orthonormal curvilinear coordinate system n and τ defined on the boundary B ε (see fig. 9), the tensors σ(u ε ) and u ε can be decomposed as (102) (103) σ(u ε ) Bε = σ nn (u ε )(n n) + σ nτ (u ε )(n τ) + σ τn (u ε )(τ n) + σ ττ (u ε )(τ τ), u ε Bε = n u n ε(n n) + τ u n ε(n τ) + n u τ ε(τ n) + τ u τ ε(τ τ). Therefore, we observe that (104) σ(u ε )n Bε = σ nn (u ε )n + σ τn (u ε )τ = 0, which implies in (105) σ nn (u ε ) = σ τn (u ε ) = σ nτ (u ε ) = 0 on B ε. since σ(u ε ) = σ(u ε ).
62 In addition, Σ ε n n Bε = 1 2 (σ(u ε) u s ε 2b u ε )In n u ε σ(u ε )n n = 1 2 (σ(u ε) u s ε 2b u ε )n n σ(u ε )n ( u ε )n = 1 2 (σ(u ε) u s ε 2b u ε ) (106) = 1 2 (σ ττ (u ε ) τ u τ ε 2b u ε ). On the other hand, the constitutive tensor C is invertible, namely, (107) C 1 = 1 ((1 + ν)i νi I)), E which implies in u s ε = C 1 σ(u ε ) τ u τ ε = 1 E (σ ττ (u ε ) νσ nn (u ε )) = 1 E σ ττ (u ε ) on B ε.
63 Therefore, we have (108) Σ ε n n Bε = 1 2 ( ) 1 E (σ ττ (u ε )) 2 2b u ε. The shape derivative of the cost functional (89) reads d (109) dε ψ(χ ε) = J Ωε (u ε ) = 1 ( ) 1 2 E (σ ττ (u ε )) 2 2b u ε. From formula (102), we have that the following expansion for σ ττ (u ε ) holds in the neighborhood of the hole (110) σ ττ (u ε (x)) = σ θθ (u ε (r,θ)) ( = ϕ ε2 r 2 B ε ) ϕ 2 (1 + 3 ε4 r 4 ) cos2θ + O(ε), σ ττ (u ε (x)) Bε = 2ϕ 1 4ϕ 2 cos2θ + O(ε).
64 The following expansion for u ε holds on the boundary of the hole (111) u ε (x) Bε = ϕ 0 + O(ε), where ϕ 0 = u( x). We obtain d 2π ( ) 1 dε ψ(ε) = E (2ϕ 1 4ϕ 2 cos2θ) 2 2b ϕ 0 εdθ + O(ε 2 ) ( ) 2 (112) = 2πε E (ϕ ϕ2 2 ) b ϕ 0 + O(ε 2 ), This leads to 1 (113) T = lim ε 0 f (ε) We choose f (ε) = πε 2, which results in [ ( ) ] 2 2πε E (ϕ ϕ2 2 ) b ϕ 0 + O(ε 2 ). (114) T = 2 E (ϕ ϕ2 2 ) + b ϕ 0. Recalling that ϕ 1 = (σ 1 (u( x)) + σ 2 (u( x)))/2, ϕ 2 = (σ 1 (u( x)) σ 2 (u( x)))/2 and ϕ 0 = u( x), the final formula for the topological derivative evaluated at x Ω becomes :
65 In terms of the principal stresses σ 1 (u( x)) and σ 2 (u( x)) (115) T ( x) = 1 [ (σ 1 (u( x)) + σ 2 (u( x))) 2 + 2E In terms of the stress tensor σ(u( x)) 2(σ 1 (u( x)) σ 2 (u( x))) 2] + b u( x). (116) T ( x) = 1 [ ] 4σ(u( x)) σ(u( x)) tr 2 σ(u( x)) +b u( x). 2E In terms of the stress tensor σ(u( x)) and the strain tensor u s ( x) (117) T ( x) = Pσ(u( x)) u s ( x) + b u( x), where P is the Pólya-Szegö polarization tensor, given in this particular case by the following isotropic fourth order tensor (118) P = ν I 1 3ν 2(1 ν 2 ) I I.
66 Multiscale Elasticity Problem The homogenization-based multi-scale constitutive framework presented, among others, in Sanchez-Palencia monograph, 1980, is adopted here in the estimation of the macroscopic elastic response from the knowledge of the underlying microstructure. The main idea behind this well-established family of constitutive theories is the assumption that any point x of the macroscopic continuum is associated to a local Representative Volume Element (RVE) whose domain is denoted by Ω µ, with boundary Ω µ, as shown in fig. 10. x FIGURE: Macroscopic continuum with a locally attached microstructure.
67 The Homogenized Elasticity Tensor Crucial to the developments presented here is the closed form of the homogenized elasticity tensor C for the multi-scale model defined in the above. The components of the homogenized elasticity tensor C, in the orthonormal basis {e i }, for i = 1,2,3, of the Euclidean space (refer to fig. 10), can be written as (119) (C) ijkl = 1 (σ µ (u µkl )) ij. V µ Ω µ where V µ denotes the total volume of the RVE, e.g. the volume of the cub in the sketch shown in fig. 10. The canonical microscopic displacement field u µkl is solution to the equilibrium equation of the form Sanchez-Palencia, 1980 (120) σ µ (u µkl ) η s = 0 η V µ, Ω µ
68 The Homogenized Elasticity Tensor We assume that the microscopic stress tensor field σ µ (u µkl ) satisfies (121) σ µ (u µkl ) = C µ u s µ kl, where C µ is the microscopic constitutive tensor given by (122) C µ = E ( µ I + ν ) µ I I, 1 + ν µ 1 2ν µ with I and I the second and fourth order identity tensors, respectively, E µ the Young modulus and ν µ Poisson ratio of the RVE. Without loss of generality, u µkl (y), with y Ω µ, may be decomposed into a sum (123) u µkl (y) := u + (e k e l )y + ũ µkl (y), of a constant (rigid) RVE displacement coinciding with the macroscopic displacement field u at the point x Ω, a linear field (e k e l )y, and a canonical microscopic displacement fluctuation field ũ µkl (y).
69 The Homogenized Elasticity Tensor The microscopic displacement fluctuation field ũ µkl is solution to the following canonical set of variational problems Sanchez-Palencia, 1980 Findũ µkl V µ, such that (124) σ µ (ũ µkl ) η s + C µ (e k s e l ) η s = 0 η V µ, Ω µ Ω µ with σ µ (u µkl ) = C µ u µkl. The complete characterization of the multi-scale constitutive model is obtained by defining the subspace V µ U µ of kinematically admissible displacement fluctuations. In general, different choices produce different macroscopic responses for the same RVE. The analysis can be focussed eg, on media with periodic microstructure. In this case, the geometry of the RVE cannot be arbitrary and must represent a cell whose periodic repetition generates the macroscopic continuum. In addition, the displacement fluctuations must satisfy periodicity on the boundary of the RVE. Accordingly, we have (125) V µ := { ϕ U µ : ϕ(y + ) = ϕ(y ) (y +,y ) P }, where P is the set of pairs of points, defined by a one-to-one periodicity correspondence, lying on opposing sides of the RVE boundary.
70 The Homogenized Elasticity Tensor Finally, the minimally constrained space of kinematically admissible displacements U µ is defined as { } (126) U µ := ϕ H 1 (Ω µ ;R 3 ) : ϕ = 0, ϕ s n = 0. Ω µ Ω µ where n is the outward unit normal to the boundary Ω µ and s denotes the symmetric tensor product between vectors. From the definition of the homogenized elasticity tensor (119), we have (C) ijkl = 1 e i σ µ (u µkl )e j V µ Ω µ = 1 (127) σ µ (u µkl ) (e i e j ). V µ Ω µ On the other hand, the additive decomposition (123) allows us to write (128) e i e j = ((e i e j )y) s since (e i e j )y = u µij (y) ũ µij (y) u. = (u µij (y) ũ µij (y) u) s = (u µij (y) ũ µij (y)) s.
71 The Homogenized Elasticity Tensor Therefore, by combining these two results, we obtain (129) (C) ijkl = 1 σ µ (u µkl ) (u µij ũ µij ) s V µ Ω µ = 1 σ µ (u µkl ) u s µ V ij, µ Ω µ since u µkl satisfies the equilibrium equation (120) and ũ µij V µ.
3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
More informationa x Questions on Classical Solutions 1. Consider an infinite linear elastic plate with a hole as shown. Uniform shear stress
Questions on Classical Solutions. Consider an infinite linear elastic plate with a hole as shown. Uniform shear stress σ xy = T is applied at infinity. Determine the value of the stress σ θθ on the edge
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 informationWeak Convergence Methods for Energy Minimization
Weak Convergence Methods for Energy Minimization Bo Li Department of Mathematics University of California, San Diego E-mail: bli@math.ucsd.edu June 3, 2007 Introduction This compact set of notes present
More information16.20 HANDOUT #2 Fall, 2002 Review of General Elasticity
6.20 HANDOUT #2 Fall, 2002 Review of General Elasticity NOTATION REVIEW (e.g., for strain) Engineering Contracted Engineering Tensor Tensor ε x = ε = ε xx = ε ε y = ε 2 = ε yy = ε 22 ε z = ε 3 = ε zz =
More informationFundamentals of Linear Elasticity
Fundamentals of Linear Elasticity Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy
More informationOn Friedrichs inequality, Helmholtz decomposition, vector potentials, and the div-curl lemma. Ben Schweizer 1
On Friedrichs inequality, Helmholtz decomposition, vector potentials, and the div-curl lemma Ben Schweizer 1 January 16, 2017 Abstract: We study connections between four different types of results that
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationDetermination of thin elastic inclusions from boundary measurements.
Determination of thin elastic inclusions from boundary measurements. Elena Beretta in collaboration with E. Francini, S. Vessella, E. Kim and J. Lee September 7, 2010 E. Beretta (Università di Roma La
More informationPEAT SEISMOLOGY Lecture 2: Continuum mechanics
PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a
More informationA short review of continuum mechanics
A short review of continuum mechanics Professor Anette M. Karlsson, Department of Mechanical ngineering, UD September, 006 This is a short and arbitrary review of continuum mechanics. Most of this material
More informationBasic concepts to start Mechanics of Materials
Basic concepts to start Mechanics of Materials Georges Cailletaud Centre des Matériaux Ecole des Mines de Paris/CNRS Notations Notations (maths) (1/2) A vector v (element of a vectorial space) can be seen
More information16.21 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive Equations
6.2 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive quations Constitutive quations For elastic materials: If the relation is linear: Û σ ij = σ ij (ɛ) = ρ () ɛ ij σ ij =
More information2 Basic Equations in Generalized Plane Strain
Boundary integral equations for plane orthotropic bodies and exterior regions G. Szeidl and J. Dudra University of Miskolc, Department of Mechanics 3515 Miskolc-Egyetemváros, Hungary Abstract Assuming
More informationSEMM Mechanics PhD Preliminary Exam Spring Consider a two-dimensional rigid motion, whose displacement field is given by
SEMM Mechanics PhD Preliminary Exam Spring 2014 1. Consider a two-dimensional rigid motion, whose displacement field is given by u(x) = [cos(β)x 1 + sin(β)x 2 X 1 ]e 1 + [ sin(β)x 1 + cos(β)x 2 X 2 ]e
More informationTopological Derivative
Outline Theory & Outline About the lecturer 2009 Ph.D. in Mathematics at the Faculty of Mathematics, University of Lodz, Poland 2010 Ph.D. in Applied Mathematics at the Elie Cartan Institute, Université
More informationME185 Introduction to Continuum Mechanics
Fall, 0 ME85 Introduction to Continuum Mechanics The attached pages contain four previous midterm exams for this course. Each midterm consists of two pages. As you may notice, many of the problems are
More informationVariational principles in mechanics
CHAPTER Variational principles in mechanics.1 Linear Elasticity n D Figure.1: A domain and its boundary = D [. Consider a domain Ω R 3 with its boundary = D [ of normal n (see Figure.1). The problem of
More informationMechanics of materials Lecture 4 Strain and deformation
Mechanics of materials Lecture 4 Strain and deformation Reijo Kouhia Tampere University of Technology Department of Mechanical Engineering and Industrial Design Wednesday 3 rd February, 206 of a continuum
More informationON LIQUID CRYSTAL FLOWS WITH FREE-SLIP BOUNDARY CONDITIONS. Chun Liu and Jie Shen
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS Volume 7, Number2, April2001 pp. 307 318 ON LIQUID CRYSTAL FLOWS WITH FREE-SLIP BOUNDARY CONDITIONS Chun Liu and Jie Shen Department
More informationEnhanced resolution in structured media
Enhanced resolution in structured media Eric Bonnetier w/ H. Ammari (Ecole Polytechnique), and Yves Capdeboscq (Oxford) Outline : 1. An experiment of super resolution 2. Small volume asymptotics 3. Periodicity
More informationAlternative numerical method in continuum mechanics COMPUTATIONAL MULTISCALE. University of Liège Aerospace & Mechanical Engineering
University of Liège Aerospace & Mechanical Engineering Alternative numerical method in continuum mechanics COMPUTATIONAL MULTISCALE Van Dung NGUYEN Innocent NIYONZIMA Aerospace & Mechanical engineering
More informationUniform estimates for Stokes equations in domains with small holes and applications in homogenization problems
Uniform estimates for Stokes equations in domains with small holes and applications in homogenization problems Yong Lu Abstract We consider the Dirichlet problem for the Stokes equations in a domain with
More informationarxiv: v1 [math.oc] 21 Nov 2018
Toplogical derivative for nonlinear magnetostatic problem Samuel Amstutz Peter Gangl November 22, 218 arxiv:1811.8715v1 [math.oc] 21 Nov 218 Abstract The topological derivative represents the sensitivity
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationIntroduction to the J-integral
Introduction to the J-integral Instructor: Ramsharan Rangarajan February 24, 2016 The purpose of this lecture is to briefly introduce the J-integral, which is widely used in fracture mechanics. To the
More informationMechanics PhD Preliminary Spring 2017
Mechanics PhD Preliminary Spring 2017 1. (10 points) Consider a body Ω that is assembled by gluing together two separate bodies along a flat interface. The normal vector to the interface is given by n
More informationTHE STOKES SYSTEM R.E. SHOWALTER
THE STOKES SYSTEM R.E. SHOWALTER Contents 1. Stokes System 1 Stokes System 2 2. The Weak Solution of the Stokes System 3 3. The Strong Solution 4 4. The Normal Trace 6 5. The Mixed Problem 7 6. The Navier-Stokes
More informationFig. 1. Circular fiber and interphase between the fiber and the matrix.
Finite element unit cell model based on ABAQUS for fiber reinforced composites Tian Tang Composites Manufacturing & Simulation Center, Purdue University West Lafayette, IN 47906 1. Problem Statement In
More informationUniform estimates for Stokes equations in domains with small holes and applications in homogenization problems
Uniform estimates for Stokes equations in domains with small holes and applications in homogenization problems Yong Lu Abstract We consider the Dirichlet problem for the Stokes equations in a domain with
More informationContinuum mechanics V. Constitutive equations. 1. Constitutive equation: definition and basic axioms
Continuum mechanics office Math 0.107 ales.janka@unifr.ch http://perso.unifr.ch/ales.janka/mechanics Mars 16, 2011, Université de Fribourg 1. Constitutive equation: definition and basic axioms Constitutive
More informationContinuum Mechanics. Continuum Mechanics and Constitutive Equations
Continuum Mechanics Continuum Mechanics and Constitutive Equations Continuum mechanics pertains to the description of mechanical behavior of materials under the assumption that the material is a uniform
More information3.2 Hooke s law anisotropic elasticity Robert Hooke ( ) Most general relationship
3.2 Hooke s law anisotropic elasticity Robert Hooke (1635-1703) Most general relationship σ = C ε + C ε + C ε + C γ + C γ + C γ 11 12 yy 13 zz 14 xy 15 xz 16 yz σ = C ε + C ε + C ε + C γ + C γ + C γ yy
More informationA MULTISCALE APPROACH IN TOPOLOGY OPTIMIZATION
1 A MULTISCALE APPROACH IN TOPOLOGY OPTIMIZATION Grégoire ALLAIRE CMAP, Ecole Polytechnique The most recent results were obtained in collaboration with F. de Gournay, F. Jouve, O. Pantz, A.-M. Toader.
More informationContinuum Mechanics Fundamentals
Continuum Mechanics Fundamentals James R. Rice, notes for ES 220, 12 November 2009; corrections 9 December 2009 Conserved Quantities Let a conseved quantity have amount F per unit volume. Examples are
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 informationLecture 2: A Strange Term Coming From Nowhere
Lecture 2: A Strange Term Coming From Nowhere Christophe Prange February 9, 2016 In this lecture, we consider the Poisson equation with homogeneous Dirichlet boundary conditions { u ε = f, x ε, u ε = 0,
More informationEXISTENCE AND UNIQUENESS FOR A THREE DIMENSIONAL MODEL OF FERROMAGNETISM
1 EXISTENCE AND UNIQUENESS FOR A THREE DIMENSIONAL MODEL OF FERROMAGNETISM V. BERTI and M. FABRIZIO Dipartimento di Matematica, Università degli Studi di Bologna, P.zza di Porta S. Donato 5, I-4126, Bologna,
More informationPoint estimates for Green s matrix to boundary value problems for second order elliptic systems in a polyhedral cone
Maz ya, V. G., Roßmann, J.: Estimates for Green s matrix 1 ZAMM Z. angew. Math. Mech. 00 2004 0, 1 30 Maz ya, V. G.; Roßmann, J. Point estimates for Green s matrix to boundary value problems for second
More informationA Locking-Free MHM Method for Elasticity
Trabalho apresentado no CNMAC, Gramado - RS, 2016. Proceeding Series of the Brazilian Society of Computational and Applied Mathematics A Locking-Free MHM Method for Elasticity Weslley S. Pereira 1 Frédéric
More informationExtension and Representation of Divergence-free Vector Fields on Bounded Domains. Tosio Kato, Marius Mitrea, Gustavo Ponce, and Michael Taylor
Extension and Representation of Divergence-free Vector Fields on Bounded Domains Tosio Kato, Marius Mitrea, Gustavo Ponce, and Michael Taylor 1. Introduction Let Ω R n be a bounded, connected domain, with
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 double layer potential
The double layer potential In this project, our goal is to explain how the Dirichlet problem for a linear elliptic partial differential equation can be converted into an integral equation by representing
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationLinear Cosserat elasticity, conformal curvature and bounded stiffness
1 Linear Cosserat elasticity, conformal curvature and bounded stiffness Patrizio Neff, Jena Jeong Chair of Nonlinear Analysis & Modelling, Uni Dui.-Essen Ecole Speciale des Travaux Publics, Cachan, Paris
More informationSOLUTIONS TO SELECTED PROBLEMS FROM ASSIGNMENTS 3, 4
SOLUTIONS TO SELECTED POBLEMS FOM ASSIGNMENTS 3, 4 Problem 5 from Assignment 3 Statement. Let be an n-dimensional bounded domain with smooth boundary. Show that the eigenvalues of the Laplacian on with
More informationStress of a spatially uniform dislocation density field
Stress of a spatially uniform dislocation density field Amit Acharya November 19, 2018 Abstract It can be shown that the stress produced by a spatially uniform dislocation density field in a body comprising
More informationInverse Transport Problems and Applications. II. Optical Tomography and Clear Layers. Guillaume Bal
Inverse Transport Problems and Applications II. Optical Tomography and Clear Layers Guillaume Bal Department of Applied Physics & Applied Mathematics Columbia University http://www.columbia.edu/ gb23 gb23@columbia.edu
More informationUNIVERSITY OF SASKATCHEWAN ME MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich
UNIVERSITY OF SASKATCHEWAN ME 313.3 MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich A CLOSED BOOK EXAMINATION TIME: 3 HOURS For Marker s Use Only LAST NAME (printed): FIRST
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 informationA FULLY COUPLED MULTISCALE SHELL FORMULATION FOR THE MODELLING OF FIBRE REINFORCED LAMINATES
ECCM-6 TH EUROPEAN CONFERENCE ON COMPOSITE MATERIALS, Seville, Spain, 22-26 June 24 A FULLY COUPLED MULTISCALE SHELL FORMULATION FOR THE MODELLING OF FIBRE REINFORCED LAMINATES J. Främby, J. Brouzoulis,
More informationLECTURE 5 APPLICATIONS OF BDIE METHOD: ACOUSTIC SCATTERING BY INHOMOGENEOUS ANISOTROPIC OBSTACLES DAVID NATROSHVILI
LECTURE 5 APPLICATIONS OF BDIE METHOD: ACOUSTIC SCATTERING BY INHOMOGENEOUS ANISOTROPIC OBSTACLES DAVID NATROSHVILI Georgian Technical University Tbilisi, GEORGIA 0-0 1. Formulation of the corresponding
More informationLaplace s Equation. Chapter Mean Value Formulas
Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic
More informationCHAPTER II MATHEMATICAL BACKGROUND OF THE BOUNDARY ELEMENT METHOD
CHAPTER II MATHEMATICAL BACKGROUND OF THE BOUNDARY ELEMENT METHOD For the second chapter in the thesis, we start with surveying the mathematical background which is used directly in the Boundary Element
More informationAn example of nonuniqueness for the continuous static unilateral contact model with Coulomb friction
An example of nonuniqueness for the continuous static unilateral contact model with Coulomb friction Un exemple de non-unicité pour le modèle continu statique de contact unilatéral avec frottement de Coulomb
More informationQuantitative Homogenization of Elliptic Operators with Periodic Coefficients
Quantitative Homogenization of Elliptic Operators with Periodic Coefficients Zhongwei Shen Abstract. These lecture notes introduce the quantitative homogenization theory for elliptic partial differential
More informationInverse Design (and a lightweight introduction to the Finite Element Method) Stelian Coros
Inverse Design (and a lightweight introduction to the Finite Element Method) Stelian Coros Computational Design Forward design: direct manipulation of design parameters Level of abstraction Exploration
More informationLecture 3. Vector fields with given vorticity, divergence and the normal trace, and the divergence and curl estimates
Lecture 3. Vector fields with given vorticity, divergence and the normal trace, and the divergence and curl estimates Ching-hsiao (Arthur) Cheng Department of Mathematics National Central University Taiwan,
More informationREMARKS ON THE VANISHING OBSTACLE LIMIT FOR A 3D VISCOUS INCOMPRESSIBLE FLUID
REMARKS ON THE VANISHING OBSTACLE LIMIT FOR A 3D VISCOUS INCOMPRESSIBLE FLUID DRAGOŞ IFTIMIE AND JAMES P. KELLIHER Abstract. In [Math. Ann. 336 (2006), 449-489] the authors consider the two dimensional
More informationLecture 8: Tissue Mechanics
Computational Biology Group (CoBi), D-BSSE, ETHZ Lecture 8: Tissue Mechanics Prof Dagmar Iber, PhD DPhil MSc Computational Biology 2015/16 7. Mai 2016 2 / 57 Contents 1 Introduction to Elastic Materials
More informationOn the well-posedness of the Prandtl boundary layer equation
On the well-posedness of the Prandtl boundary layer equation Vlad Vicol Department of Mathematics, The University of Chicago Incompressible Fluids, Turbulence and Mixing In honor of Peter Constantin s
More informationarxiv: v2 [math.ap] 25 Aug 2016
Team organization may help swarms of flies to become invisible in closed waveguides Lucas Chesnel 1, Sergei A. Nazarov2, 3, 4 arxiv:1510.04540v2 [math.ap] 25 Aug 2016 1 INRIA/Centre de mathématiques appliquées,
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 informationFurther Linear Elasticity
Torsion of cylindrical bodies Further Linear Elasticity Problem Sheet # 1. Consider a cylindrical body of length L, the ends of which are subjected to distributions of tractions that are statically equivalent
More informationA new closed-form model for isotropic elastic sphere including new solutions for the free vibrations problem
A new closed-form model for isotropic elastic sphere including new solutions for the free vibrations problem E Hanukah Faculty of Mechanical Engineering, Technion Israel Institute of Technology, Haifa
More informationThis introductory chapter presents some basic concepts of continuum mechanics, symbols and notations for future reference.
Chapter 1 Introduction to Elasticity This introductory chapter presents some basic concepts of continuum mechanics, symbols and notations for future reference. 1.1 Kinematics of finite deformations We
More informationConstitutive Equations
Constitutive quations David Roylance Department of Materials Science and ngineering Massachusetts Institute of Technology Cambridge, MA 0239 October 4, 2000 Introduction The modules on kinematics (Module
More informationHydrodynamic Limit with Geometric Correction in Kinetic Equations
Hydrodynamic Limit with Geometric Correction in Kinetic Equations Lei Wu and Yan Guo KI-Net Workshop, CSCAMM University of Maryland, College Park 2015-11-10 1 Simple Model - Neutron Transport Equation
More informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6.1 Mathematik Preprint Nr. 155 A posteriori error estimates for stationary slow flows of power-law fluids Michael Bildhauer,
More information7 The Navier-Stokes Equations
18.354/12.27 Spring 214 7 The Navier-Stokes Equations In the previous section, we have seen how one can deduce the general structure of hydrodynamic equations from purely macroscopic considerations and
More informationTechnische Universität Graz
Technische Universität Graz Toplogical derivative for nonlinear magnetostatic problem S. Amstutz, P. Gangl Berichte aus dem Institut für Angewandte Mathematik Bericht 218/4 Technische Universität Graz
More informationON SINGULAR PERTURBATION OF THE STOKES PROBLEM
NUMERICAL ANALYSIS AND MATHEMATICAL MODELLING BANACH CENTER PUBLICATIONS, VOLUME 9 INSTITUTE OF MATHEMATICS POLISH ACADEMY OF SCIENCES WARSZAWA 994 ON SINGULAR PERTURBATION OF THE STOKES PROBLEM G. M.
More informationFinite Elements for Elastic Shell Models in
Elastic s in Advisor: Matthias Heinkenschloss Computational and Applied Mathematics Rice University 13 April 2007 Outline Elasticity in Differential Geometry of Shell Geometry and Equations The Plate Model
More informationArchimedes Center for Modeling, Analysis & Computation. Singular solutions in elastodynamics
Archimedes Center for Modeling, Analysis & Computation Singular solutions in elastodynamics Jan Giesselmann joint work with A. Tzavaras (University of Crete and FORTH) Supported by the ACMAC project -
More informationLinearized theory of elasticity
Linearized theory of elasticity Arie Verhoeven averhoev@win.tue.nl CASA Seminar, May 24, 2006 Seminar: Continuum mechanics 1 Stress and stress principles Bart Nowak March 8 2 Strain and deformation Mark
More informationApplied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.
Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R
More informationApplied Numerical Mathematics. High-order numerical schemes based on difference potentials for 2D elliptic problems with material interfaces
Applied Numerical Mathematics 111 (2017) 64 91 Contents lists available at ScienceDirect Applied Numerical Mathematics www.elsevier.com/locate/apnum High-order numerical schemes based on difference potentials
More informationEXISTENCE AND REGULARITY OF SOLUTIONS FOR STOKES SYSTEMS WITH NON-SMOOTH BOUNDARY DATA IN A POLYHEDRON
Electronic Journal of Differential Equations, Vol. 2017 (2017), No. 147, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE AND REGULARITY OF SOLUTIONS FOR
More informationSHARP BOUNDARY TRACE INEQUALITIES. 1. Introduction
SHARP BOUNDARY TRACE INEQUALITIES GILES AUCHMUTY Abstract. This paper describes sharp inequalities for the trace of Sobolev functions on the boundary of a bounded region R N. The inequalities bound (semi-)norms
More informationINSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES. Note on the fast decay property of steady Navier-Stokes flows in the whole space
INSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES Note on the fast decay property of stea Navier-Stokes flows in the whole space Tomoyuki Nakatsuka Preprint No. 15-017 PRAHA 017 Note on the fast
More informationWe briefly discuss two examples for solving wave propagation type problems with finite differences, the acoustic and the seismic problem.
Excerpt from GEOL557 Numerical Modeling of Earth Systems by Becker and Kaus 2016 1 Wave propagation Figure 1: Finite difference discretization of the 2D acoustic problem. We briefly discuss two examples
More informationOPTIMAL CONTROL AND STRANGE TERM FOR A STOKES PROBLEM IN PERFORATED DOMAINS
PORTUGALIAE MATHEMATICA Vol. 59 Fasc. 2 2002 Nova Série OPTIMAL CONTROL AND STRANGE TERM FOR A STOKES PROBLEM IN PERFORATED DOMAINS J. Saint Jean Paulin and H. Zoubairi Abstract: We study a problem of
More informationProceedings of Meetings on Acoustics
Proceedings of Meetings on Acoustics Volume 19, 2013 http://acousticalsociety.org/ ICA 2013 Montreal Montreal, Canada 2-7 June 2013 Engineering Acoustics Session 1aEA: Thermoacoustics I 1aEA7. On discontinuity
More informationSome Remarks on the Reissner Mindlin Plate Model
Some Remarks on the Reissner Mindlin Plate Model Alexandre L. Madureira LNCC Brazil Talk at the Workshop of Numerical Analysis of PDEs LNCC February 10, 2003 1 Outline The 3D Problem and its Modeling Full
More informationu xx + u yy = 0. (5.1)
Chapter 5 Laplace Equation The following equation is called Laplace equation in two independent variables x, y: The non-homogeneous problem u xx + u yy =. (5.1) u xx + u yy = F, (5.) where F is a function
More informationEFFECTIVE CHARACTERISTICS OF POROUS MEDIA AS A FUNCTION OF POROSITY LEVEL
AMS Subject Classification Index: 74Q20 EFFECTIVE CHARACTERISTICS OF POROUS MEDIA AS A FUNCTION OF POROSITY LEVEL Irini DJERAN-MAIGRE and Sergey V.KUZNETSOV INSA de Lyon, URGC 34 Avenue des Arts, 69621
More informationGeometry for Physicists
Hung Nguyen-Schafer Jan-Philip Schmidt Tensor Analysis and Elementary Differential Geometry for Physicists and Engineers 4 i Springer Contents 1 General Basis and Bra-Ket Notation 1 1.1 Introduction to
More informationLINEAR FLOW IN POROUS MEDIA WITH DOUBLE PERIODICITY
PORTUGALIAE MATHEMATICA Vol. 56 Fasc. 2 1999 LINEAR FLOW IN POROUS MEDIA WITH DOUBLE PERIODICITY R. Bunoiu and J. Saint Jean Paulin Abstract: We study the classical steady Stokes equations with homogeneous
More informationRegularity for Poisson Equation
Regularity for Poisson Equation OcMountain Daylight Time. 4, 20 Intuitively, the solution u to the Poisson equation u= f () should have better regularity than the right hand side f. In particular one expects
More informationRecovery of anisotropic metrics from travel times
Purdue University The Lens Rigidity and the Boundary Rigidity Problems Let M be a bounded domain with boundary. Let g be a Riemannian metric on M. Define the scattering relation σ and the length (travel
More informationElements of Continuum Elasticity. David M. Parks Mechanics and Materials II February 25, 2004
Elements of Continuum Elasticity David M. Parks Mechanics and Materials II 2.002 February 25, 2004 Solid Mechanics in 3 Dimensions: stress/equilibrium, strain/displacement, and intro to linear elastic
More information2.073 Notes on finite-deformation theory of isotropic elastic-viscoplastic solids
2.073 Notes on finite-deformation theory of isotropic elastic-viscoplastic solids Lallit Anand Department of Mechanical Engineering Massachusetts Institute of Technology Cambridge, MA 02139, USA April
More informationHeterogeneous structures studied by interphase elasto-damaging model.
Heterogeneous structures studied by interphase elasto-damaging model. Giuseppe Fileccia Scimemi 1, Giuseppe Giambanco 1, Antonino Spada 1 1 Department of Civil, Environmental and Aerospace Engineering,
More informationOther state variables include the temperature, θ, and the entropy, S, which are defined below.
Chapter 3 Thermodynamics In order to complete the formulation we need to express the stress tensor T and the heat-flux vector q in terms of other variables. These expressions are known as constitutive
More informationTraces, extensions and co-normal derivatives for elliptic systems on Lipschitz domains
Traces, extensions and co-normal derivatives for elliptic systems on Lipschitz domains Sergey E. Mikhailov Brunel University West London, Department of Mathematics, Uxbridge, UB8 3PH, UK J. Math. Analysis
More information4 Divergence theorem and its consequences
Tel Aviv University, 205/6 Analysis-IV 65 4 Divergence theorem and its consequences 4a Divergence and flux................. 65 4b Piecewise smooth case............... 67 4c Divergence of gradient: Laplacian........
More informationReceiver. Johana Brokešová Charles University in Prague
Propagation of seismic waves - theoretical background Receiver Johana Brokešová Charles University in Prague Seismic waves = waves in elastic continuum a model of the medium through which the waves propagate
More informationIntroduction to fracture mechanics
Introduction to fracture mechanics Prof. Dr. Eleni Chatzi Dr. Giuseppe Abbiati, Dr. Konstantinos Agathos Lecture 6-9 November, 2017 Institute of Structural Engineering, ETH Zu rich November 9, 2017 Institute
More information6.1 Formulation of the basic equations of torsion of prismatic bars (St. Venant) Figure 6.1: Torsion of a prismatic bar
Module 6 Torsion Learning Objectives 6.1 Formulation of the basic equations of torsion of prismatic bars (St. Venant) Readings: Sadd 9.3, Timoshenko Chapter 11 e e 1 e 3 Figure 6.1: Torsion of a prismatic
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 information