Outline Intro EDPs Método de Galerkin Ejemplos

Métodos isogeométricos para resolver

ecuaciones en derivadas parciales
Eduardo M. Garau

Facultad de Ingenierı́a Quı́mica - UNL

Consejo Nacional de Investigaciones
Cientı́ficas y Técnicas

Seminario para la Maestrı́a en Matemática de la FaCET (UNT)

San Miguel de Tucumán, 27 de noviembre de 2018

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Outline

1 Splines

2 Forma variacional de EDPs

3 Aproximaciones de Galerkin

4 Algunos ejemplos

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Outline

1 Splines

2 Forma variacional de EDPs

3 Aproximaciones de Galerkin

4 Algunos ejemplos

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Espacio de splines unidimensionales

d ∈ N: grado polinomial
∆ := {ξi }N
i=1 : breakpoints,

a = ξ1 < ξ2 < · · · < ξN−1 < ξN = b.

r := {ri }N−1
i=2 : regularidad en los breakpoints

ri ∈ {−1, 0, 1, . . . , d − 1}.

Sdr (∆): espacio de splines

s : [a, b] → R tal que
s|[ξ ,ξ ∈ Pd , para i = 1, . . . , N − 1.
i i+1 ]

s es ri veces continuamente diferenciable en ξi , para

i = 2, . . . , N − 1.

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Casos particulares

Polinomios a trozos discontinuos de grado d

(rj = −1, ∀ j)

Polinomios a trozos continuos de grado d

(rj = 0, ∀ j)

Splines de máxima regularidad de grado d

(rj = d − 1, ∀ j)

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Algunas definiciones

Supongamos que r1 = rN = −1.

M := (m1 , . . . , mN ) vector de multiplicidad en los breakpoints,

mj := d − rj .

N−1
X
dim Sdr (∆) = mi =: n.
i=1

t := {ti }n+d+1
i=1 : d-open knot vector (sucesión), dada por

t = {ξ1 , . . . , ξ1 , ξ2 , . . . , ξ2 , . . . ξN , . . . , ξN }.
| {z } | {z } | {z }
m1 veces m2 veces mN veces

Notemos que

a = t1 = · · · = td+1 < td+2 ≤ · · · ≤ tn < tn+1 = · · · = tn+d+1 = b.

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Una base para Sdr (∆): B-splines

Propiedades de la base de B-splines B := {B1 , . . . , Bn }

Bi depende sólo del knot vector local {ti , ti+1 , . . . , ti+d+1 }.
supp Bi = [ti , ti+d+1 ]
n
X
Partición de la unidad: Bi (x) = 1, ∀ x ∈ [a, b].
i=1

Ejemplo 1: d = 2 con nodos internos simples

1

0.9

0.8

0.7

0.6

0.5

0.4

0.3

0.2

0.1

0
xi_1=xi_2=xi_3 xi_4 xi_5 xi_6 xi_7 xi_8 xi_9=xi_10=xi_11

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Una base para Sdr (∆): B-splines

Propiedades de la base de B-splines B := {B1 , . . . , Bn }

Bi depende sólo del knot vector local {ti , ti+1 , . . . , ti+d+1 }.
supp Bi = [ti , ti+d+1 ]
n
X
Partición de la unidad: Bi (x) = 1, ∀ x ∈ [a, b].
i=1

Ejemplo 1: d = 2 con nodos internos simples

1

0.9

0.8

0.7

0.6

0.5

0.4

0.3

0.2

0.1

0
xi_1=xi_2=xi_3 xi_4 xi_5 xi_6 xi_7 xi_8 xi_9=xi_10=xi_11

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Una base para Sdr (∆): B-splines

Propiedades de la base de B-splines B := {B1 , . . . , Bn }

Bi depende sólo del knot vector local {ti , ti+1 , . . . , ti+d+1 }.
supp Bi = [ti , ti+d+1 ]
n
X
Partición de la unidad: Bi (x) = 1, ∀ x ∈ [a, b].
i=1

Ejemplo 1: d = 2 con nodos internos simples

1

0.9

0.8

0.7

0.6

0.5

0.4

0.3

0.2

0.1

0
xi_1=xi_2=xi_3 xi_4 xi_5 xi_6 xi_7 xi_8 xi_9=xi_10=xi_11

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Una base para Sdr (∆): B-splines

Propiedades de la base de B-splines B := {B1 , . . . , Bn }

Bi depende sólo del knot vector local {ti , ti+1 , . . . , ti+d+1 }.
supp Bi = [ti , ti+d+1 ]
n
X
Partición de la unidad: Bi (x) = 1, ∀ x ∈ [a, b].
i=1

Ejemplo 1: d = 2 con nodos internos simples

1

0.9

0.8

0.7

0.6

0.5

0.4

0.3

0.2

0.1

0
xi_1=xi_2=xi_3 xi_4 xi_5 xi_6 xi_7 xi_8 xi_9=xi_10=xi_11

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Una base para Sdr (∆): B-splines

Propiedades de la base de B-splines B := {B1 , . . . , Bn }

Bi depende sólo del knot vector local {ti , ti+1 , . . . , ti+d+1 }.
supp Bi = [ti , ti+d+1 ]
n
X
Partición de la unidad: Bi (x) = 1, ∀ x ∈ [a, b].
i=1

Ejemplo 1: d = 2 con nodos internos simples

1

0.9

0.8

0.7

0.6

0.5

0.4

0.3

0.2

0.1

0
xi_1=xi_2=xi_3 xi_4 xi_5 xi_6 xi_7 xi_8 xi_9=xi_10=xi_11

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Una base para Sdr (∆): B-splines

Propiedades de la base de B-splines B := {B1 , . . . , Bn }

Bi depende sólo del knot vector local {ti , ti+1 , . . . , ti+d+1 }.
supp Bi = [ti , ti+d+1 ]
n
X
Partición de la unidad: Bi (x) = 1, ∀ x ∈ [a, b].
i=1

Ejemplo 1: d = 2 con nodos internos simples

1

0.9

0.8

0.7

0.6

0.5

0.4

0.3

0.2

0.1

0
xi_1=xi_2=xi_3 xi_4 xi_5 xi_6 xi_7 xi_8 xi_9=xi_10=xi_11

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Una base para Sdr (∆): B-splines

Propiedades de la base de B-splines B := {B1 , . . . , Bn }

Bi depende sólo del knot vector local {ti , ti+1 , . . . , ti+d+1 }.
supp Bi = [ti , ti+d+1 ]
n
X
Partición de la unidad: Bi (x) = 1, ∀ x ∈ [a, b].
i=1

Ejemplo 1: d = 2 con nodos internos simples

1

0.9

0.8

0.7

0.6

0.5

0.4

0.3

0.2

0.1

0
xi_1=xi_2=xi_3 xi_4 xi_5 xi_6 xi_7 xi_8 xi_9=xi_10=xi_11

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Una base para Sdr (∆): B-splines

Propiedades de la base de B-splines B := {B1 , . . . , Bn }

Bi depende sólo del knot vector local {ti , ti+1 , . . . , ti+d+1 }.
supp Bi = [ti , ti+d+1 ]
n
X
Partición de la unidad: Bi (x) = 1, ∀ x ∈ [a, b].
i=1

Ejemplo 1: d = 2 con nodos internos simples

1

0.9

0.8

0.7

0.6

0.5

0.4

0.3

0.2

0.1

0
xi_1=xi_2=xi_3 xi_4 xi_5 xi_6 xi_7 xi_8 xi_9=xi_10=xi_11

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Una base para Sdr (∆): B-splines

Propiedades de la base de B-splines B := {B1 , . . . , Bn }

Bi depende sólo del knot vector local {ti , ti+1 , . . . , ti+d+1 }.
supp Bi = [ti , ti+d+1 ]
n
X
Partición de la unidad: Bi (x) = 1, ∀ x ∈ [a, b].
i=1

Ejemplo 1: d = 2 con nodos internos simples

1

0.9

0.8

0.7

0.6

0.5

0.4

0.3

0.2

0.1

0
xi_1=xi_2=xi_3 xi_4 xi_5 xi_6 xi_7 xi_8 xi_9=xi_10=xi_11

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Una base para Sdr (∆): B-splines

Propiedades de la base de B-splines B := {B1 , . . . , Bn }

Bi depende sólo del knot vector local {ti , ti+1 , . . . , ti+d+1 }.
supp Bi = [ti , ti+d+1 ]
n
X
Partición de la unidad: Bi (x) = 1, ∀ x ∈ [a, b].
i=1

Ejemplo 2: d = 2 con un nodo interno doble

1

0.9

0.8

0.7

0.6

0.5

0.4

0.3

0.2

0.1

0
xi_1=xi_2=xi_3 xi_4 xi_5=xi_6 xi_7 xi_8 xi_9 xi_10=xi_11=xi_12

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Una base para Sdr (∆): B-splines

Propiedades de la base de B-splines B := {B1 , . . . , Bn }

Bi depende sólo del knot vector local {ti , ti+1 , . . . , ti+d+1 }.
supp Bi = [ti , ti+d+1 ]
n
X
Partición de la unidad: Bi (x) = 1, ∀ x ∈ [a, b].
i=1

Ejemplo 2: d = 2 con un nodo interno doble

1

0.9

0.8

0.7

0.6

0.5

0.4

0.3

0.2

0.1

0
xi_1=xi_2=xi_3 xi_4 xi_5=xi_6 xi_7 xi_8 xi_9 xi_10=xi_11=xi_12

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Una base para Sdr (∆): B-splines

Propiedades de la base de B-splines B := {B1 , . . . , Bn }

Bi depende sólo del knot vector local {ti , ti+1 , . . . , ti+d+1 }.
supp Bi = [ti , ti+d+1 ]
n
X
Partición de la unidad: Bi (x) = 1, ∀ x ∈ [a, b].
i=1

Ejemplo 2: d = 2 con un nodo interno doble

1

0.9

0.8

0.7

0.6

0.5

0.4

0.3

0.2

0.1

0
xi_1=xi_2=xi_3 xi_4 xi_5=xi_6 xi_7 xi_8 xi_9 xi_10=xi_11=xi_12

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Una base para Sdr (∆): B-splines

Propiedades de la base de B-splines B := {B1 , . . . , Bn }

Bi depende sólo del knot vector local {ti , ti+1 , . . . , ti+d+1 }.
supp Bi = [ti , ti+d+1 ]
n
X
Partición de la unidad: Bi (x) = 1, ∀ x ∈ [a, b].
i=1

Ejemplo 2: d = 2 con un nodo interno doble

1

0.9

0.8

0.7

0.6

0.5

0.4

0.3

0.2

0.1

0
xi_1=xi_2=xi_3 xi_4 xi_5=xi_6 xi_7 xi_8 xi_9 xi_10=xi_11=xi_12

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Una base para Sdr (∆): B-splines

Propiedades de la base de B-splines B := {B1 , . . . , Bn }

Bi depende sólo del knot vector local {ti , ti+1 , . . . , ti+d+1 }.
supp Bi = [ti , ti+d+1 ]
n
X
Partición de la unidad: Bi (x) = 1, ∀ x ∈ [a, b].
i=1

Ejemplo 2: d = 2 con un nodo interno doble

1

0.9

0.8

0.7

0.6

0.5

0.4

0.3

0.2

0.1

0
xi_1=xi_2=xi_3 xi_4 xi_5=xi_6 xi_7 xi_8 xi_9 xi_10=xi_11=xi_12

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Una base para Sdr (∆): B-splines

Propiedades de la base de B-splines B := {B1 , . . . , Bn }

Bi depende sólo del knot vector local {ti , ti+1 , . . . , ti+d+1 }.
supp Bi = [ti , ti+d+1 ]
n
X
Partición de la unidad: Bi (x) = 1, ∀ x ∈ [a, b].
i=1

Ejemplo 2: d = 2 con un nodo interno doble

1

0.9

0.8

0.7

0.6

0.5

0.4

0.3

0.2

0.1

0
xi_1=xi_2=xi_3 xi_4 xi_5=xi_6 xi_7 xi_8 xi_9 xi_10=xi_11=xi_12

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Una base para Sdr (∆): B-splines

Propiedades de la base de B-splines B := {B1 , . . . , Bn }

Bi depende sólo del knot vector local {ti , ti+1 , . . . , ti+d+1 }.
supp Bi = [ti , ti+d+1 ]
n
X
Partición de la unidad: Bi (x) = 1, ∀ x ∈ [a, b].
i=1

Ejemplo 2: d = 2 con un nodo interno doble

1

0.9

0.8

0.7

0.6

0.5

0.4

0.3

0.2

0.1

0
xi_1=xi_2=xi_3 xi_4 xi_5=xi_6 xi_7 xi_8 xi_9 xi_10=xi_11=xi_12

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Una base para Sdr (∆): B-splines

Propiedades de la base de B-splines B := {B1 , . . . , Bn }

Bi depende sólo del knot vector local {ti , ti+1 , . . . , ti+d+1 }.
supp Bi = [ti , ti+d+1 ]
n
X
Partición de la unidad: Bi (x) = 1, ∀ x ∈ [a, b].
i=1

Ejemplo 2: d = 2 con un nodo interno doble

1

0.9

0.8

0.7

0.6

0.5

0.4

0.3

0.2

0.1

0
xi_1=xi_2=xi_3 xi_4 xi_5=xi_6 xi_7 xi_8 xi_9 xi_10=xi_11=xi_12

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Una base para Sdr (∆): B-splines

Propiedades de la base de B-splines B := {B1 , . . . , Bn }

Bi depende sólo del knot vector local {ti , ti+1 , . . . , ti+d+1 }.
supp Bi = [ti , ti+d+1 ]
n
X
Partición de la unidad: Bi (x) = 1, ∀ x ∈ [a, b].
i=1

Ejemplo 2: d = 2 con un nodo interno doble

1

0.9

0.8

0.7

0.6

0.5

0.4

0.3

0.2

0.1

0
xi_1=xi_2=xi_3 xi_4 xi_5=xi_6 xi_7 xi_8 xi_9 xi_10=xi_11=xi_12

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Una base para Sdr (∆): B-splines

Propiedades de la base de B-splines B := {B1 , . . . , Bn }

Bi depende sólo del knot vector local {ti , ti+1 , . . . , ti+d+1 }.
supp Bi = [ti , ti+d+1 ]
n
X
Partición de la unidad: Bi (x) = 1, ∀ x ∈ [a, b].
i=1

Ejemplo 2: d = 2 con un nodo interno doble

1

0.9

0.8

0.7

0.6

0.5

0.4

0.3

0.2

0.1

0
xi_1=xi_2=xi_3 xi_4 xi_5=xi_6 xi_7 xi_8 xi_9 xi_10=xi_11=xi_12

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Knot insertion

0.9

0.8

0.7

0.6

0.5
1
0.4
0.9
0.3
0.8
0.2
0.7
0.1
0.6
0
xi_1=xi_2=xi_3 xi_4 xi_5=xi_6 xi_7 xi_8 xi_9 xi_10=xi_11=xi_12
0.5
1
0.4
0.9
0.3
0.8
0.2
0.7
0.1
0.6
0
xi_1=xi_2=xi_3 xi_4 xi_5 xi_6 xi_7 xi_8 xi_9=xi_10=xi_11 0.5

0.4

0.3

0.2

0.1

0
xi_1=xi_2=xi_3 xi_4 xi_5 xi_6 xi_7 xi_8 xi_9 xi10=xi_11=xi_12

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Knot insertion y relación de dos escalas

1 1

0.8 0.8

0.6 0.6

0.4 0.4

0.2 0.2

0 0
0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1

1 1

0.8 0.8

0.6 0.6

0.4 0.4

0.2 0.2

0 0
0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Espacios de splines producto tensor

Sdr11 (∆1 ) and Sdr22 (∆2 ): Dos espacios de splines unidimensionales

B 1 and B 2 : Correspondientes bases de B-splines

Sdr (∆1 , ∆2 ): Espacio de splines producto tensor

s : [0, 1] × [0, 1] → R tal que s ∈ span B, donde

B := {B : [0, 1]×[0, 1] → R | B(x, y) = B 1 (x)B 2 (y ), B 1 ∈ B 1 , B 2 ∈ B 2 }.

(B se llama base de B-splines)

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Funciones de base (B-splines)

d1 = d2 = 2 (Bicuadráticas, globalmente C 1 )
Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Funciones de base (B-splines)

d1 = d2 = 2 (Bicuadráticas, globalmente C 1 )
Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Funciones de base (B-splines)

d1 = d2 = 2 (Bicuadráticas, globalmente C 1 )
Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Funciones de base (B-splines)

d1 = d2 = 2 (Bicuadráticas, globalmente C 1 )
Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Funciones de base (B-splines)

d1 = d2 = 2 (Bicuadráticas, globalmente C 1 )
Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Funciones de base (B-splines)

d1 = d2 = 2 (Bicuadráticas, globalmente C 1 )
Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Funciones de base (B-splines)

d1 = d2 = 2 (Bicuadráticas, globalmente C 1 )
Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Funciones de base (B-splines)

d1 = d2 = 2 (Bicuadráticas, globalmente C 1 )
Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Funciones de base (B-splines)

d1 = d2 = 2 (Bicuadráticas, globalmente C 1 )
Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Funciones de base (B-splines)

d1 = d2 = 2 (Bicuadráticas, globalmente C 1 )
Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Funciones de base (B-splines)

d1 = d2 = 2 (Bicuadráticas, globalmente C 1 )
Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Funciones de base (B-splines)

d1 = d2 = 2 (Bicuadráticas, globalmente C 1 )
Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Funciones de base (B-splines)

d1 = d2 = 2 (Bicuadráticas, globalmente C 1 )
Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Funciones de base (B-splines)

d1 = d2 = 2 (Bicuadráticas, globalmente C 1 )
Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Funciones de base (B-splines)

d1 = d2 = 2 (Bicuadráticas, globalmente C 1 )
Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Funciones de base (B-splines)

d1 = d2 = 2 (Bicuadráticas, globalmente C 1 )
Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Funciones de base (B-splines)

d1 = d2 = 2 (Bicuadráticas, globalmente C 1 )
Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Funciones de base (B-splines)

d1 = d2 = 2 (Bicuadráticas, globalmente C 1 )
Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Funciones de base (B-splines)

d1 = d2 = 2 (Bicuadráticas, globalmente C 1 )
Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Funciones de base (B-splines)

d1 = d2 = 2 (Bicuadráticas, globalmente C 1 )
Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Funciones de base (B-splines)

d1 = d2 = 2 (Bicuadráticas, globalmente C 1 )
Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Funciones de base (B-splines)

d1 = d2 = 2 (Bicuadráticas, globalmente C 1 )
Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Funciones de base (B-splines)

d1 = d2 = 2 (Bicuadráticas, globalmente C 1 )
Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Funciones de base (B-splines)

d1 = d2 = 2 (Bicuadráticas, globalmente C 1 )
Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Funciones de base (B-splines)

d1 = d2 = 2 (Bicuadráticas, globalmente C 1 )
Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Funciones de base (B-splines)

d1 = d2 = 2 (Bicuadráticas, globalmente C 1 )
Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Funciones de base (B-splines)

d1 = d2 = 2 (Bicuadráticas, globalmente C 1 )
Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Funciones de base (B-splines)

d1 = d2 = 2 (Bicuadráticas, globalmente C 1 )
Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Funciones de base (B-splines)

d1 = d2 = 2 (Bicuadráticas, globalmente C 1 )
Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Funciones de base (B-splines)

d1 = d2 = 2 (Bicuadráticas, globalmente C 1 )
Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Funciones de base (B-splines)

d1 = d2 = 2 (Bicuadráticas, globalmente C 1 )
Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Funciones de base (B-splines)

d1 = d2 = 2 (Bicuadráticas, globalmente C 1 )
Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Funciones de base (B-splines)

d1 = d2 = 2 (Bicuadráticas, globalmente C 1 )
Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Funciones de base (B-splines)

d1 = d2 = 2 (Bicuadráticas, globalmente C 1 )
Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Funciones de base (B-splines)

d1 = d2 = 2 (Bicuadráticas, globalmente C 1 )
Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Funciones de base (B-splines)

d1 = d2 = 2 (Bicuadráticas, globalmente C 1 )
Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Outline

1 Splines

2 Forma variacional de EDPs

3 Aproximaciones de Galerkin

4 Algunos ejemplos

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Caso estacionario de la ecuación de difusión con reacción

Consideramos
−∆u + u = f ,

 en Ω

u = 0,

(FF ) sobre Γ0
 ∂u
= ∇u · n = 0,

 sobre Γ1
∂n
∂2u ∂2u
donde ∆u := ∂x12
+ ∂x22
es el laplaciano de u.

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Formulación variacional - Forma débil del problema

Si u es solución clásica de (FF) entonces

− ∆u(x) + u(x) = f (x)

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Formulación variacional - Forma débil del problema

Si u es solución clásica de (FF) entonces

− ∆u(x) v (x) + u(x) v (x) = f (x) v (x)

para toda v ∈ C 1 (Ω), con v |Γ0 ≡ 0.

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Formulación variacional - Forma débil del problema

Si u es solución clásica de (FF) entonces

Z Z Z
− ∆u(x) v (x)dx + u(x) v (x)dx = f (x) v (x)dx
Ω Ω Ω

para toda v ∈ C 1 (Ω), con v |Γ0 ≡ 0.

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Formulación variacional - Forma débil del problema

Si u es solución clásica de (FF) entonces

Z Z Z
− ∆u(x) v (x)dx + u(x) v (x)dx = f (x) v (x)dx
Ω Ω Ω

para toda v ∈ C 1 (Ω), con v |Γ0 ≡ 0.

Usando la Identidad de Green en el primer término:

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Formulación variacional - Forma débil del problema

Si u es solución clásica de (FF) entonces

Z Z Z
− ∆u(x) v (x)dx + u(x) v (x)dx = f (x) v (x)dx
Ω Ω Ω

para toda v ∈ C 1 (Ω), con v |Γ0 ≡ 0.

Usando la Identidad de Green en el primer término:
Z Z Z Z Z
∂u ∂u
− v
|{z} ∂n ds− v ds+ ∇u(x)∇v (x)dx+ u(x)v (x)dx = f (x)v (x)dx
Γ0 ∂n
Γ1 |{z} Ω Ω Ω
=0
=0

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Formulación variacional - Forma débil del problema

Si u es solución clásica de (FF) entonces

Z Z Z
− ∆u(x) v (x)dx + u(x) v (x)dx = f (x) v (x)dx
Ω Ω Ω

para toda v ∈ C 1 (Ω), con v |Γ0 ≡ 0.

Usando la Identidad de Green en el primer término:
Z Z Z Z Z
∂u ∂u
− v
|{z} ∂n ds− v ds+ ∇u(x)∇v (x)dx+ u(x)v (x)dx = f (x)v (x)dx
Γ0 ∂n
Γ1 |{z} Ω Ω Ω
=0
=0

Luego,
Z Z
∀ v ∈ C 1 (Ω) : v |Γ0 ≡ 0.


∇u(x)∇v (x)+u(x)v (x) dx = f (x)v (x)dx,
Ω Ω

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Espacio de Sobolev

El espacio de Sobolev de orden 1 está dado por

H 1 (Ω) := {v ∈ L2 (Ω) | existe ∇v en sentido débil y ∇v ∈ L2 (Ω)}.

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Espacio de Sobolev

El espacio de Sobolev de orden 1 está dado por

H 1 (Ω) := {v ∈ L2 (Ω) | existe ∇v en sentido débil y ∇v ∈ L2 (Ω)}.

Es un espacio de Hilbert, con producto interno dado por

Z


(u, v )H 1 (Ω) := ∇u(x)∇v (x) + u(x)v (x) dx,
Ω

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Espacio de Sobolev

El espacio de Sobolev de orden 1 está dado por

H 1 (Ω) := {v ∈ L2 (Ω) | existe ∇v en sentido débil y ∇v ∈ L2 (Ω)}.

Es un espacio de Hilbert, con producto interno dado por

Z


(u, v )H 1 (Ω) := ∇u(x)∇v (x) + u(x)v (x) dx,
Ω

y ası́, la norma está dada por

Z
2
|∇v (x)|2 + [v (x)]2 dx = k∇v k2L2 + kv k2L2 .


kv kH 1 Ω =
Ω

Nótese que H 1 (Ω) ,→ L2 (Ω).

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Formulación variacional - Forma débil del problema

Consideramos V := {v ∈ H 1 (Ω) | v |Γ0 ≡ 0}, subespacio cerrado de

H 1 (Ω), y ası́, un espacio de Hilbert con el producto interno heredado
de H 1 (Ω):
Z


(u, v )V = (u, v )H 1 (Ω) = ∇u(x)∇v (x) + u(x)v (x) dx.
Ω

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Formulación variacional - Forma débil del problema

Consideramos V := {v ∈ H 1 (Ω) | v |Γ0 ≡ 0}, subespacio cerrado de

H 1 (Ω), y ası́, un espacio de Hilbert con el producto interno heredado
de H 1 (Ω):
Z


(u, v )V = (u, v )H 1 (Ω) = ∇u(x)∇v (x) + u(x)v (x) dx.
Ω

Si f ∈ L2 (Ω), Z
F (v ) := f (x)v (x)dx
Ω

es un funcional lineal y acotado sobre V.

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Formulación variacional - Forma débil del problema

Consideramos V := {v ∈ H 1 (Ω) | v |Γ0 ≡ 0}, subespacio cerrado de

H 1 (Ω), y ası́, un espacio de Hilbert con el producto interno heredado
de H 1 (Ω):
Z


(u, v )V = (u, v )H 1 (Ω) = ∇u(x)∇v (x) + u(x)v (x) dx.
Ω

Si f ∈ L2 (Ω), Z
F (v ) := f (x)v (x)dx
Ω

es un funcional lineal y acotado sobre V.

Forma débil (FD)

u∈V : (u, v )V = F (v ), ∀ v ∈ V.

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Obeservaciones sobre (FD)

Relación entre las soluciones de (FF) y (FD)

Si u es solución de (FF) entonces u es solución de (FV), es
decir, si u es solución clásica entonces es solución débil.
Se puede demostrar que si u es solución débil y u ∈ C 2 (Ω)
entonces u es solución clásica.

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Obeservaciones sobre (FD)

Relación entre las soluciones de (FF) y (FD)

Si u es solución de (FF) entonces u es solución de (FV), es
decir, si u es solución clásica entonces es solución débil.
Se puede demostrar que si u es solución débil y u ∈ C 2 (Ω)
entonces u es solución clásica.

Caracterı́sticas de (FD)
Por el teorema de Representación de Riesz, (FD) tiene una
única solución.
El problema (FD) es estable, es decir, kukV ≤ kf kL2 .

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Outline

1 Splines

2 Forma variacional de EDPs

3 Aproximaciones de Galerkin

4 Algunos ejemplos

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Aproximaciones de Galerkin

Si S es un subespacio de dimensión finita de V, consideramos

Aproximación de Galerkin (G)

U∈S : (U, V )V = F (V ), ∀ V ∈ S.

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Aproximaciones de Galerkin
PN
Si {B1 , B2 , . . . , BN } una base de S, entonces U = j=1 Uj Bj .

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Aproximaciones de Galerkin
PN
Si {B1 , B2 , . . . , BN } una base de S, entonces U = j=1 Uj Bj .

Ues solución de (G) ⇐⇒ (U, Bi )V = F (Bi ), i = 1, . . . , N

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Aproximaciones de Galerkin
PN
Si {B1 , B2 , . . . , BN } una base de S, entonces U = j=1 Uj Bj .

Ues solución de (G) ⇐⇒ (U, Bi )V = F (Bi ), i = 1, . . . , N

N
X
⇐⇒ ( Uj Bj , Bi )V = F (Bi ), i = 1, . . . , N
j=1

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Aproximaciones de Galerkin
PN
Si {B1 , B2 , . . . , BN } una base de S, entonces U = j=1 Uj Bj .

Ues solución de (G) ⇐⇒ (U, Bi )V = F (Bi ), i = 1, . . . , N

N
X
⇐⇒ ( Uj Bj , Bi )V = F (Bi ), i = 1, . . . , N
j=1
N
X
⇐⇒ Uj (Bj , Bi )V = F (Bi ), i = 1, . . . , N
j=1
| {z } | {z }
=:Kij =:Fi

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Aproximaciones de Galerkin
PN
Si {B1 , B2 , . . . , BN } una base de S, entonces U = j=1 Uj Bj .

Ues solución de (G) ⇐⇒ (U, Bi )V = F (Bi ), i = 1, . . . , N

N
X
⇐⇒ ( Uj Bj , Bi )V = F (Bi ), i = 1, . . . , N
j=1
N
X
⇐⇒ Uj (Bj , Bi )V = F (Bi ), i = 1, . . . , N
j=1
| {z } | {z }
=:Kij =:Fi

K11 U1 + K12 U2 + . . . + K1N UN = F1




 K21 U1 + K22 U2 + . . . + K2N UN = F2


⇐⇒ ..



 .

KN1 U1 + KN2 U2 + . . . + KNN UN = FN

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Aproximaciones de Galerkin
PN
Si {B1 , B2 , . . . , BN } una base de S, entonces U = j=1 Uj Bj .

Ues solución de (G) ⇐⇒ (U, Bi )V = F (Bi ), i = 1, . . . , N

N
X
⇐⇒ ( Uj Bj , Bi )V = F (Bi ), i = 1, . . . , N
j=1
N
X
⇐⇒ Uj (Bj , Bi )V = F (Bi ), i = 1, . . . , N
j=1
| {z } | {z }
=:Kij =:Fi

K11 U1 + K12 U2 + . . . + K1N UN = F1




 K21 U1 + K22 U2 + . . . + K2N UN = F2


⇐⇒ ..



 .

KN1 U1 + KN2 U2 + . . . + KNN UN = FN
⇐⇒ K U = F, donde K = (Kij ), U = (Ui ), F = (Fi ).

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

K: Matriz del sistema

K es simétrica, pues Kij = (Bj , Bi )V = (Bi , Bj )V = Kji .

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

K: Matriz del sistema

K es simétrica, pues Kij = (Bj , Bi )V = (Bi , Bj )V = Kji .

Si V = (Vi ) y W = (Wi ) entonces
N N N N N
VT · K W
P P P P P
= Vi (K W)i = Vi Kij Wj = Vi (Bj , Bi )V Wj
i=1 i=1 j=1 i=1 j=1
N
P PN PN PN
= Vi ( j=1 Wj Bj , Bi )V = ( j=1 Wj Bj , i=1 Vi Bi )V .
i=1

PN PN
Si V = i=1 Vi Bi y W = i=1 Wi Bi ,

VT · K W = (V , W )V .

Además,
VT · K V = kV k2V ≥ 0,
y ası́, K es definida positiva, y por lo tanto, invertible.

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

K: Matriz del sistema

K es simétrica, pues Kij = (Bj , Bi )V = (Bi , Bj )V = Kji .

Si V = (Vi ) y W = (Wi ) entonces
N N N N N
VT · K W
P P P P P
= Vi (K W)i = Vi Kij Wj = Vi (Bj , Bi )V Wj
i=1 i=1 j=1 i=1 j=1
N
P PN PN PN
= Vi ( j=1 Wj Bj , Bi )V = ( j=1 Wj Bj , i=1 Vi Bi )V .
i=1

PN PN
Si V = i=1 Vi Bi y W = i=1 Wi Bi ,

VT · K W = (V , W )V .

Además,
VT · K V = kV k2V ≥ 0,
y ası́, K es definida positiva, y por lo tanto, invertible.

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Aproximaciones de Galerkin: Resumen

Si S es un subespacio de dimensión finita de V, consideramos

Aproximación de Galerkin (G)

U∈S : (U, V )V = F (V ), ∀ V ∈ S.
Si {B1 , B2 , . . . , BN } es una base de S, (G) no es mas que el sistema
de ecuaciones lineales

K U = F, donde K = ((Bj , Bi )V ), U = coordenadas(U), F = (F (Bi )).

La matriz K del sistema es simétrica y definida positiva.

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Optimalidad de la aproximación de Galerkin

U es la mejor aproximación a u en el espacio S, es decir,

ku − UkV = mı́n ku − V kV .
V ∈S

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Optimalidad de la aproximación de Galerkin

U es la mejor aproximación a u en el espacio S, es decir,
ku − UkV = mı́n ku − V kV .
V ∈S

Demostración.
(FD) (u, v )V = F (v ), v ∈ V.
(G) (U, V )V = F (V ), V ∈ S ⊂ V.
Restando (FD) - (G),
S

(u − U, V )V = 0, para toda V ∈ S.

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Optimalidad de la aproximación de Galerkin

U es la mejor aproximación a u en el espacio S, es decir,

ku − UkV = mı́n ku − V kV .
V ∈S

Demostración.
Si V ∈ S,

ku − Uk2V = (u − U, u − U)V
= (u − U, u − V )V + (u − U, V − U )V
| {z }
∈S
| {z }
=0
≤ ku − UkV ku − V kV .

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Outline

1 Splines

2 Forma variacional de EDPs

3 Aproximaciones de Galerkin

4 Algunos ejemplos

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Problemas de minimización de funcionales sobre superficies

A familia de superficies admisibles.

J : A → R.

Problema de minimización
Hallar Γ∗ = arg mı́n J(Γ).
Γ∈A

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Superficies mı́nimas

Ejemplo 1: Funcional de área

Z
J(Γ) = dΓ,
Γ

Con Γ ∈ A = A(γ) : superficies con un borde fijo γ.

Análisis Isogeométrico
Outline Intro EDPs Método de Galerkin Ejemplos

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Un algoritmo de quasi-Newton

1 Tomar Γ ∈ A (superficie inical).

2 Hallar u ∈ VΓ tal que

J 00 (Γ)[u][v ] = −J 0 (Γ)[v ], para todo v ∈ VΓ .

3 Actualizar superficie Γ → Γ + Pk (un).
4 Repetir los dos pasos anteriores hasta obtener una solución
satisfactoria.

Análisis Isogeométrico
Outline Intro EDPs Método de Galerkin Ejemplos

Análisis Isogeométrico
Outline Intro EDPs Método de Galerkin Ejemplos

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Discretización espacial

S = span B, donde B = {Φ̂1 , ..., Φ̂N }, (Φ̂i : [0, 1]2 → R) es una

base de splines producto tensor de regularidad global k, (grado
mayor o igual a k + 1), con k ≥ 3.

Superficies admisibles
nX o
S = span ci Φ̂i | ci ∈ R3

Funciones tests: Espacio sobre la superficie Γ := X ([0, 1]2 )

SΓ ⊂ span{Φi = Φ̂i ◦ X −1 | Φ̂i ∈ B}

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

The Willmore energy

Γ superficie admisible
h curvatura media
κ curvatura Gaussiana
b > 0 constante
Z
1 2
W (Γ) = h para ∂Γ = ∅.
Γ 2
Z
1 2
Wb (Γ) = h − bκ para ∂Γ 6= ∅.
Γ 2

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Minimization problem

γ curva en R3 (Borde dado)

ns campo vectorial sobre γ
A = A(γ, ns ) familia de superficies admisibles con borde fijo γ y
normal ns prescripta en el borde.

Problema de minimización sin restricciones

Hallar Γ∗ = arg mı́n W (Γ).
Γ∈A

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Minimization problem with constraints

γ curva en R3 (Borde dado)

ns campo vectorial sobre γ
A = A(γ, ns ) familia de superficies admisibles con borde fijo γ y
normal ns prescripta en el borde.
A0 > 0 dado
V0 > 0 dado

Problema de minimización con restricciones

Hallar Γ∗ = arg mı́n{W (Γ) | Γ ∈ A, A(Γ) = A0 , V (Γ) = V0 }.

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Minimization problem with constraints

γ curva en R3 (Borde dado)

ns campo vectorial sobre γ
A = A(γ, ns ) familia de superficies admisibles con borde fijo γ y
normal ns prescripta en el borde.
A0 > 0 dado
V0 > 0 dado

Problema de minimización con restricciones

Hallar Γ∗ = arg mı́n{W (Γ) | Γ ∈ A, A(Γ) = A0 , V (Γ) = V0 }.

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Ejemplo 1: Un plano perturbado

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Ejemplo 2: Pico 1

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Ejemplo 3: Pico 2

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Ejemplo 4: Globo

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Ejemplo 4: Globo

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Ejemplo 4: Globo

Análisis Isogeométrico
 Outline Intro EDPs Método de Galerkin Ejemplos

Ejemplo 4: Globo

Análisis Isogeométrico
Outline Intro EDPs Método de Galerkin Ejemplos

¡Gracias por su atención!

Análisis Isogeométrico