×

Close

- Finite Element Methods - FEM
- Note
**Jawaharlal Nehru Technological University Anantapur (JNTU) College of Engineering (CEP), Pulivendula, Pulivendula, Andhra Pradesh, India - JNTUACEP**- 11 Topics
**103 Views**- 4 Offline Downloads
- Uploaded 1 year ago

Touch here to read

Page-2

- The A bstract problem - ( 5 - 8 )
- Examples - ( 9 - 28 )
- The finite element method inits simplest form - ( 29 - 35 )
- Examples of finite elements - ( 36 - 52 )
- General properties of finite elements - ( 53 - 59 )
- Interpolation theory in solvable spaces - ( 60 - 66 )
- Applications to second-order problems over polygonal domains - ( 67 - 76 )
- Numerical integration - ( 77 - 94 )
- The obstacle problem - ( 95 - 102 )
- Conforming finite element method for the plate problem - ( 103 - 112 )
- Non- conforming methods for the plate problem - ( 113 - 128 )

Topic:

1. The Abstract Problem 2 Definition 1.1. Let V be a normed linear space. A bilinear form a(·, ·) on V is said to be V-elliptic if there exists a constant α > 0 such that for all v ∈ V. a(v, v) ≥ α||v||2 . (1.4) 2 Theorem 1.1. Let V be a Banach space and K a closed convex subset of V. Let a(·, ·) be V-elliptic. Then there exists a unique solution for the problem (P). Further this solution is characterised by the property: (1.5) a(u, v − u) ≥ f (v − u) for all v ∈ K. Remark 1.1. The inequalities (1.5) are known as variational inequalities. Proof. The V-ellipticity of a(·, ·) clearly implies that if a(v, v) = 0 then v = 0. This together with the symmetry and bilinearity of a(·, ·) shows that a(·, ·) defines an inner-product on V. Further the continuity and the V-ellipticity of a(·, ·) shows that the norm 1 v ∈ V → a(v, v) 2 (1.6) defined by the inner-product is equivalent to the existing norm on V. Thus V acquires the structure of a Hilbert space and we apply the Riesz representation theorem to obtain the following: for all f ∈ V ′ , there exists σ f ∈ V such that (1.7) f (v) = a(σ f, v) for all v ∈ V. The map σ : V ′ → V given by f 7→ σ f is linear. Now, 1 a(v, v) − f (v) 2 1 = a(v, v) − a(σ f, v) 2 J(v) =

1. The Abstract Problem 3 1 1 = a(v − σ f, v − σ f ) − a(σ f, σ f ). 2 2 3 The symmetry of a(·, ·) is essential in obtaining the last equality. For a given f , since σ f is fixed, J is minimised if and only if a(v − σ f, v − σ f ) is minimised. But this being the distance between v and σ f , our knowledge of Hilbert space theory tells us that since K is a closed convex subset, there exists a unique element u ∈ K such that this minimum is obtained. This proves the existence and uniqueness of the solution, which is merely the projection of σ f over K. We know that this projection is characterised by the inequalities: (1.8) a(σ f − u, v − u) ≤ 0 for all v ∈ K. Geometrically, this means that the angle between the vectors (σ f −u) and (v − u) is obtuse. See Fig. 1.1. Figure 1.1: Thus, a(σ f, v − u) ≤ a(u, v − u) which by virtue of (1.7) is precisely the relation (1.5). This completes the proof. We can state the following Corollary 1.1. (a) If K is a non-empty closed convex cone with vertex at origin 0, then the solution of (P) is characterised by: a(u, v) ≥ f (v) for all v ∈ K (1.9) a(u, u) = f (u).

1. The Abstract Problem 4 (b) If K is a subspace of V then the solution is characterised by (1.10) a(u, v) = f (v) for all v ∈ K. Remark 1.2. The relations (1.5), (1.9) and (1.10) are all called variational formulations of the problem (P). Proof. (a) If K is a cone with vertex at 0, then for u, v ∈ K, u + v ∈ K. (cf. Fig. 1.2). If u is the solution to (P), then for all v ∈ K applying (1.5) to (u + v) we get a(u, v) ≥ f (v) for all v ∈ K. In particular this applies to u itself. Setting v = 0 in (1.5) we get −a(u, u) ≥ − f (u) which gives the reverse inequality necessary to complete the proof of (1.9). Conversely, if (1.9) holds, we get (1.5) by just subtracting one inequality from the other. Figure 1.2: (b) Applying (a) to K, since any subspace is a cone with vertex at 0, we get (b) immediately. For if v ∈ K, then −v ∈ K and applying (1.9) both to v and −v we get (1.10). This completes the proof. Remark 1.3. The solution u of (P) corresponding to f ∈ V ′ (for a fixed a(·, ·)) defines a map V ′ → V. Since this solution is the projection of 4

1. The Abstract Problem 5 σ f on K, it follows that the above map is linear if and only if K is a subspace. The problems associated with variational inequalities are therefore nonlinear in general. Exercise 1.1. Let V be as in Theorem 1.1. For f1 , f2 ∈ V ′ , let u1 , u2 5 be the corresponding solutions of (P). If || · ||∗ denotes the norm in V ′ , prove that 1 ||u1 − u2 || ≤ || f1 − f2 ||∗ . α Remark 1.4. The above exercise shows, in particular, the continuous dependence of the solution of f , in the sense described above. This together with the existence and uniqueness establishes that the problem (P) is “well-posed” in the terminology of partial differential equations. Exercise 1.2. If V is a normed linear space, K a given convex subset of V and J : V → R any functional which is once differentiable everywhere, then (i) if u ∈ K is such that J(u) = Min J(v), u satisfies, J ′ (u)(v − u) ≥ 0 v∈K for all v ∈ K. (ii) Conversely, if u ∈ K such that J ′ (u)(v − u) ≥ 0 for all v ∈ K, and J is everywhere twice differentiable with J ′′ satisfying J ′′ (v)(w, w) ≥ α||w||2 , for all v, w ∈ K and some α ≥ 0, then J(u) = Min J(v). v∈K Exercise 1.3 (1 ). Apply the previous exercise to the functional 1 a(v, v) = f (v) 2 with a(·, ·) and f as in Theorem 1.1. If K is a subspace of V, show that J ′ (u)(v) = 0 for all v ∈ K. In particular if K = V, J ′ (u) = 0. J(v) = It was essentially the symmetry of the bilinear form which provided the Hilbert space structure in Theorem 1.1. We now drop the symmetry assumption on a(·, ·) but we assume V to be a Hilbert space. In addition we assume that K = V. Theorem 1.2 (LAX-MILGRAM LEMMA). Let V be a Hilbert space. 6 1 Exercises (1.2)(i) and 1.3 together give relations (1.5)

## Leave your Comments