Universidad Nacional de Trujillo

Facultad de Ciencias Fı́sicas y Matemáticas

Escuela Profesional de Matemáticas

EXISTENCIA Y UNICIDAD DE SOLUCIONES EN PROBLEMAS

VARIACIONALES DE SEGUNDO ORDEN Y SU SOLUCIÓN POR EL
MÉTODO DE RITZ

Proyecto de Tesis para optar el Tı́tulo de Licenciado en Matemática

Presentado por:
Leyli del Pilar Salas Plasencia

Asesor:
Dr. Obidio Elisban Rubio Mercedes

Trujillo - Perú
2022
escudo-eps-converted-to.pdf

Universidad Nacional de Trujillo

Facultad de Ciencias Fı́sicas y Matemáticas
Escuela Profesional de Matemáticas

I. GENERALIDADES
1. Titulo : Existencia y Unicidad de soluciones en Problemas Variacionales de Segundo
Orden y su solución por el Método de Ritz

2. Personal del Proyecto

Investigador Cargo Institución
Leyli del Pilar Salas Plasencia Alumno Escuela Profesional de Matemáticas
Luis Alberto Lara Romero Asesor Departamento Académico de Matemáticas

3. Tipo de investigación:

a) Por el fin que persigue : Básica

b) Por el diseño de contrastación : Explicativa

4. Departamento al que pertence el Proyecto : Matemáticas

5. Lugar e Institución donde se desarrollará el Proyecto

a) Universidad: Universidad Nacional de Trujillo

b) Facultad : Facultad de Ciencias Fı́sicas y Matemáticas

c) Departamento : Matemáticas

6. Duración del proyecto: Un Semestre Académico

a) Cronograma de Trabajo

F-M01.03.03-DRT/PG-01-(Rev.1) 2
escudo-eps-converted-to.pdf

Cronograma de Ejecución
Descripción MESES/ HORAS SEMANALES
Actividades/ Tareas S1 S2 S3 S4 S5 S6 S7 S8 S9 S10 S11 S12 S13 S14 S15 S16
Revisión de la Bibliografı́a X
Espacios funcionales X X
Formulación Variacional X
Lema de Lax-Milgram X X
Aplicación del Método de Ritz X X
Ejemplos de aplicación de Ritz X X
MEF X X
Aplicación del MEF X X
Implementación del MEF X X
Redacción preliminar X X
Redacción del informe X X
Sustentación del informe X

b) Fecha de Inicio- Fecha de Término: 01/11/2022 - 30/04/2023

Horas Semanales: 2 horas con el Asesor y 16 horas en forma independiente.

7. Recursos Disponibles

a) Personal

1. Responsable(s) del Proyecto

2. Personal de apoyo al Proyecto

b) Materiales y Equipos

(1). De Consumo

• Disponibles

(a). Libros especializados de turbulencia.

(b). Acceso a revistas electrónicas

• No disponibles

(a). Cuadernillos espiralados tamaño A4 de 80 hojas cada uno.

(b). Resaltadores Faber - Castell de colores.

(c). Bibliografı́a especializados.

(d). CDs, DVD, PenDrive

(2). De Inversión

• PC Pentium IV 800 de 2.4 Ghz HD 80 Gb 1000 RAM.

• Impresora HP Deskjet 858C (800x300 DPI 220V PUSB)

• Software Cientı́fico Fortran90 y MatLab.

(3). Servicios

F-M01.03.03-DRT/PG-01-(Rev.1) 3
escudo-eps-converted-to.pdf

• Fotocopiado

• Encuadernado y empastado

• Tarifa de servicios de telefonı́a

• Congresos, Cursos, Simposios, Conferencias.

8. Presupuesto

a) Recursos Disponibles

Clasificador Denominación Cantidad Costo

5.3.11.30 Bibliografı́a especializada 4 000.00
5.3.11.36 Internet 700.00
SUB TOTAL 4 700.00

a) Recursos no Disponibles
Clasificador Denominación Cantidad Costo
5.3.11.30 Bienes de Consumo
Bibiliografı́a especializada 06 unidades 2 500
Sofware Cientı́fico 01 unidad 2 000
Materiales de escritorio
Cuadernos espiralados A4 10 unidades 60.00
Transparencias 01 millar 112.00
Papel Bond A4 80 grs. 05 millar 220.00
Folder de plástico A4 40 unidades 402.00
Lápiceros de tinta seca 03 docena 108.00
Lápices HB 01 docena 10.00
Resaltadores Faber-Castell 01 docena 40.00
Correctores Pentel 06 unidades 50.00
Borradores Faber-Castell 06 unidades 20.00
Tajadores de lápiz 06 unidad 12.00
Plumones para pizarra acrĺica 01 docena 96.00
Tinta para plumones 06 unidades 90.00
Grapas estándar 01 caja 15.00
CDs, DVDs 04 cajas 40.00
Cinta de Impresora 04 unidades 84.00
SUB TOTAL 5 829.00

F-M01.03.03-DRT/PG-01-(Rev.1) 4
escudo-eps-converted-to.pdf

Clasificador Denominación Cantidad Costo

5.3.11.36 Tarifa de servicios generales
Internet 1000 horas 2 200.00
SUB TOTAL 2 200.00

Clasificador Denominación Cantidad Costo

5.3.11.39 Otros servicios de terceros
Servicio de transporte 20 pasajes. 2 200.00
Tipeado 200 páginas. 900.00
Impresión 200 páginas 400.00
Fotocopias 900 páginas 250.00
Encuadernación y empastado 10 juegos 300.00
Gastos menores
de rápida cancelación 500.00
SUB TOTAL 4 550.00

Clasificador Denominación Cantidad Costo

5.3.11.51 Equipamiento y bienes
duraderos
PC Multimdia Pentium IV 01 unidad 4 000.00
Impresora HP Deskjet 858C 01 unidad 1 200.00
SUB TOTAL 5 200.00

Tabla de Resumen de Presupuesto Subtotal

Recursos disponibles 4 700.00
Recursos no disponibles 17 779.00
TOTAL 23 608.00
9. Financiamiento

a) Con recursos propios : autofinaciado por el investigador hasta un 25 % por ciento

del total.

b) Con recursos universitarios : 20 % por ciento del total.

c) Sin financiamiento: 55 %.

F-M01.03.03-DRT/PG-01-(Rev.1) 5
escudo-eps-converted-to.pdf

II. PLAN DE INVESTIGACION

1. Situación problemática

En la actualiada se viene desarrollando una teorı́a matemática llamada Teorı́a moderna

de las ecuaciones diferenciales, la cual esta siendo utilizada, no por su simplicidad sino
por su efectiva resolución a problemas reales.

Para ilustrar mejor las cosas, consideremos el siguiente problema de valor de contorno
!
d du
− κ(x) = f (x) x∈Ω (1)
dx dx
u(0) = u(1) = 0 (2)

donde κ ∈ C(Ω), seria natural preguntarnos: dónde se encuentra la solución del problema
(1)- (2)?. pensariamos primero, como serı́a lógico que dicha solución se encuentra en un
espacio de funciones, es decir, en un espacio funcional, ahora, si la función f es continua,
buscariamos la solución en un espacio de funciones continuas, es decir, u ∈ C 2 (Ω) ∩ C0 (Ω)
la cual es llamada solución clásica del problema (1)- (2).

Ahora bien, cabe preguntarse que pasaria si la función f no es continua, cómo superar este
inconveniente, para ello consideremos un espacio de funciones que no son necesariamente
continuas, para las cuales generaliamos el concepto de derivada usual a una derivada
generalizada. Los espacios a los que nos referimos son los espacios L2 (Ω), H 2 (Ω), H0 (Ω)
que tienen derivadas en el sentido generalizado, entonces la solución se busca en el espacio
H 2 (Ω) ∩ H0 (Ω).

Luego de la construcción de los espacios apropiados, el siguiente paso es formular variacio-

nalmente el problema (1)- (2), para ello se considera el espacio de funciones C0∞ (Ω), que
denota el espacio de funciones infinitamente diferenciables y que poseen soporte compacto
sobre Ω. Una vez formulado variacionalmente el problema (1)- (2), estamos listos para
utiliar un método de aproximación, tal como el Método de Ritz, que nos permitira encon-
trar una solución aproximada a la solución analı́tica del problema utilizando funciones
seccionalmente polinomiales lo que ha traido consigo enormes ventajas computacionales y
perfectamente aplicables como es el Método de los Elementos Finitos. Ntouyas y Tsamatos
[10] presentan la solución de problemas de valor inicial de la forma
 ′
′ ′ ′
ρ(t)x (t) + f (t, x(t), x(σ(t), x (t), x (g(t)))) = 0, t ∈ [a, b] (3)

Vojin [1] presenta una derivación elemental de una extensión del método de Ritz para
funciones prueba que no satisfacen las condiciones de contorno. Como en el enfoque de
Babuska-Brezzi, las condiciones de contorno se tratan como restricciones variacionales
donde los multiplicadores de Lagrange son utilizados para eliminarlos. Müller [2] estima
el error del método de Deep Ritz para ecuaciones elı́pticas. Para las condicionees de

F-M01.03.03-DRT/PG-01-(Rev.1) 6
escudo-eps-converted-to.pdf

Dirichlet, se estima el error cuando los valores en la frontera son impuestas a través del
método de penalización de frontera. Basak [3] utiliza principios variacionales, el método de
Rayleigh-Ritz convierte problemas de valor de contorno en un problema de minimizar una
funcional. Yakovlev [4] presenta una cota de error del método de Ritz para el problema
de minimización de la funcional
Z 1h i2 Z 1 Z 1
′
J(u) = u (t) dt + q(t)u (t)dt − 2
2
f (t)u(t)dt (4)
0 0 0

en el espacio W12 (0, 1) en el caso donde la hipótesis de continuida de q(t) se reemplaza por
la condición q 2 (t)t(1−t) ∈ L(0, 1) Akl et.al [5] introduce el método de Rayleigh-Ritz (RR)
con diferentes funciones base y compara este método con otros métodos numéricos para
resolver problemas de valor de contorno de segundo orden para describir como este método
está logrando una alta precisión, usando funciones base lineal, cuadrática, Hermite cúbica,
B-Spline cúbica y funciones polinómicas de diferentes tamaños de paso mostrando como la
elección de las funciones prueba afecta o reduce el error. Deshmukh y Shukla [6] presenta
una técnica para determinar la solución de un problema de valor inicial de segundo orden
(BVP) que frecuentemente se presenta en ingenierı́a, por múltiples métodos variacionales.
Además las soluciones obtenidas por estos métodos se comparan con la solución exacta
de la ecuación diferencial de segundo orden. Presenta el Método de subdominio, Método
de colocación, Método de Mı́nimo Cuadrados, Petrov Galerkin, Método de Galerkin,
Método de Rayleigh-Ritz y el Método de Diferencias Finitas. Las soluciones obtenidas
por los métodos numéricos estan bastante de acuerdo con la solución exacta. Phys [7]
proporciona la teorı́a general de razones de convergencia para el método variacional de
Rayleigh-Ritz para los estados fundamentales del átomo y moléculas. La teorı́a muestra
qué funciones deben agregarse al conjunto básico para mejorar la tasa de convergencia y
brinda fórmulas explı́citas para estimar las correcciones de las energı́as variacionales y las
funciones de onda. Presenta una aplicación al cálculo de CI sobre el estado fundamental
del átomo de helio. Bougoffa [8] considera ecuaciones diferenciales de segundo orden no
lineales con condiciones de contorno no locales. Primero reformula el problema de valor de
contorno como un problema de punto fijo para el operador integral de Fredholm, y luego
presenta un resultado sobre la existencia y unicidad de la solución usando el teorema de
mapeo de contracción. Kaufmann [9] considera la existencia y unicidad de las soluciones
de problemas de valor de contorno de segundo orden
′′
x (t) = f (t, x(t), x(2)(t) ), a≤t≤b (5)

donde la principal herramienta empleada es el teorema del punto fijo de Shauder. Tisdell
[11] investiga la nolinealidad de ecuaciones diferenciales de segundo orden sujeta a varias
condiciones de contorno. Introduce una condición de contorno, a priori en las derivadas
de la solución al problema. En particular, las condiciones de crecimiento cuadrático en
el lado derecho de la ecuación diferencial que no son utilizadas. Esta ideas son aplicadas
para garantizar la existencia de al menos una solución.

F-M01.03.03-DRT/PG-01-(Rev.1) 7
escudo-eps-converted-to.pdf

En el marco general del Método de Ritz iniciamos dando algunas definiciones. Sea V un
espacio de Hilbert, a(., .) : V × V → R una forma bilineal ( procedente de la formulación
variacional de un problema de contorno) y l ∈ V ′ (lado derecho del problema variacional).

Entonces se quiere encontrar u ∈ V tal que

a(u, v) = l(v) ∀v ∈V (6)

Se supone que la forma bilineal a(., .) es acotada y V elı́ptica, es decir, existen constantes
k y α tal que:
|a(u, v)| ≤ k∥u∥V ∥v∥V ∀v ∈V (7)
y
a(v, v) ≥ α∥v∥2V ∀v ∈V (8)
El Método de Ritz busca una aproximación un en un espacio de Hilbert Vn , que es solución
de la ecuación (6), minimizando una funcional J sobre un espacio finito dimensional
Vn . Esto conduce a un sistema de ecuaciones lineales, cuya matriz de coefecientes es
definida positiva, para cualquier elección de la base para Vn . La forma bilineal a(., .) es
una integral que genera a una matriz K esparcida si el soporte de las funciones base es
tal que a(φi , φj ) = 0 cuando (sop φi ) ∩ (sop φj ) = ϕ.

Tenemos en este punto un aspecto importante: escoger bases atractivas desde el punto de
vista computacional, esto es, bases que conduzcan a una matriz de coeficientes esparcida
y con cierta estructura que facilite la resolución del sistema de ecuaciones lineales.

El Método de Ritz consiste en tomar un subespacio cerrado Vn del espacio de Hilbert V ,

entonces Vn es por sı́ mismo un espacio de Hilbert y por lo tanto las hipótesis del Lema
de Lax-Milgram se cumplen para a(..) y l(.) en Vn siendo a(., .) simétrica.

Sea {Vn } una sucesión de subespacios cerrados del espacio de Hilbert V . Entones existe
un único elemento ũn ∈ Vn tal que

q(ũn ) = mı́n q(u) (9)

u∈Vn

donde la funcional q esta dada por

1
q(u) = a(u, u) − l(u) u ∈ Vn (10)
2
Si bien es cierto no tenemos un procedimiento para encontrar ũn ∈ Vn . Sin embargo,
podemos utilizar facilmente la dimensionalidad de Vn para obtener un procedimiento
general para encontrar ũn ∈ Vn . Ası́, asumimos que a(., .) es simétrica con {φi (x)}ni=1
una base para el espacio Vn , es decir, cualquier conjunto de n funciones linealmente
independientes en Vn . Entonces para cada un ∈ Vn se tiene un único desarrollo en la
forma n
un = αi φi (x) (11)
X

i=1

F-M01.03.03-DRT/PG-01-(Rev.1) 8
escudo-eps-converted-to.pdf

donde αi ∈ R. Sustituyendo (16) en (15) y usando la propiedad de bilinealidad de a(., .)

y la propiedad de linealidad de l(.), se tiene la funcional cuadrática en R⋉ :
1
Q(α) = αT Kα − αT B α ∈ R⋉ (12)
2
Observar que la funcional cuadrática (12) está definida en R⋉ , luego el problema se puede
escribir como un problema de minimización en R⋉ en la forma:

Q(α̃) = mı́n⋉ Q(α) (13)

α∈R

Para minimizar la ecuación (13) es necesario hallar el punto estacionario que satisfaga
∇Q(α) = 0 que es equivalente a resolver un sistema de ecuaciones de la forma Kα = B
donde los coeficientes α son tal que corresponden a la función de aproximación
n
ũn = αi φi (x) (14)
X

i=1

Justificación

Muchos modelos matemáticos en ingenieria vienen siendo obtenidos de manera empı́rica,

asumiendo hipótesis generadas en laboratorios. Dichos modelos debido a su importancia
requieren un tratamiento y estudio matemático previo, para garantizar la existencia y
unicidad de soluciones. Existen muchas técnicas para ello, ya sean modelos lineales o no
lineales, uno de ellos son las formulaciones débiles usando el procedimiento de Galerkin,
ello debido a su facilidad al momento de implementar numéricamente lo ha llevado al
desarrollo de uno de los métodos numéricos más utilizados para resolver problemas de
matemáticas e ingenierı́a como es el Método de los Elementos Finitos que viene imple-
mentado en muchos programas libre y comerciales.

En nuestra investigación vamos a utilizar el Método de Ritz que consiste en buscar el

minimizador de una funcional lineal que proviene de la formulación variacional de un
problema eliptico de segundo orden.

Propósito

Se pretende demostrar la existencia y unicidad de problemas variacionales de segundo

orden y buscar su solución aproximada utilizando el Método de Ritz buscando el minimi-
zador de una funcional cuadrática.

Presentaremos teoremas básicos para la formulación variacional de un problema de valor

en la frontera, como también estudiaremos los teoremas fundamentales del método de
Ritz.

El trabajo termina con la aplicación a problemas elı́pticos de segundo orden en una

variable en el plano.

F-M01.03.03-DRT/PG-01-(Rev.1) 9
escudo-eps-converted-to.pdf

2. Enunciado del Problema

Nos planteamos el siguiente problema: Es posible demostrar la existencia y unicidad de

soluciones en problemas variacionales de segundo orden y encontrar su solución por el
Método de Ritz?

3. Objetivos

General Demostrar la existencia y unicidad de soluciones en problemas variacionales de

segundo orden y encontrar su solución aproximada por el método de Ritz

Especı́fico

Escribir la forma variacional de problemas de segundo orden utilizando el método

de residuos ponderados

Aplicar el Lema de Lax-Milgram para demostrar la existencia y unicidad en los

problemas variacionales de segundo orden

Encontrar soluciones aproximadas al problema variacional utiliando el método de

Ritz

Implementar un algoritmo para encontrar soluciones aproximadas al problema va-

riacional

4. Hipótesis

Garantizando las hipótesis del Lema de Lax-Milgram afirmamos que si es posible demos-
trar la existencia y unicidad de soluciones débiles a problemas variacionales de segundo
orden y también podemos aplicar el método de Ritz para encontrar soluciones aproxima-
das al problema variacional.

5. Marco Teórico Conceptual

Presentamos algunas definiciones para dar fundamento al Proyecto de Investigación.

Lema 0.0.1 (Lax-Milgram)

Sea a(., .) una forma bilineal acotada y coerciva y l(.) una forma lineal acotada, entonces
existe un único elemento ũ ∈ V tal que

a(ũ, v) = l(v) ∀v ∈ V (15)

Si además a(., .) es simétrica, entonces ũ es el único minimizador global de la funcional

1
q(u) = a(u, u) − l(u) ∀v ∈ V (16)
2

F-M01.03.03-DRT/PG-01-(Rev.1) 10
escudo-eps-converted-to.pdf

Métodos de Aproximación
Dados los fundamentos matemáticos y la teorı́a básica sobre la formulación variacional
de problemas de contorno, en este capı́tulo se introduce un método para encontrar una
solución aproximada del problema variacional (el método de Galerkin) para posterior-
mente dar un caso especial del mismo, el Método de los Elementos Finitos. El método de
Galerkin resulta ser un caso particular del método de los residuos ponderados (MRP).

Como marco general, sea V un espacio de Hilbert, a(., .) : V × V → R una forma bilineal
′
(procedente de la formulación débil de un problema de contorno) y l ∈ V (lado derecho
del problema débil).

Entonces se quiere encontrar u ∈ V tal que

a(u, v) = l(v) ∀v ∈V (17)

Se supone que la forma bilineal a(., .) es acotada y V elı́ptica, es decir, existen constantes
k y α tal que
|a(u, v)| ≤ k∥u∥V ∥v∥V ∀v ∈V (18)

a(v, v) ≥ α∥v∥2V ∀v ∈V (19)

0.0.1. Método de Ritz (1908)

Walter Ritz (1878 − 1909), fı́sico suizo fue el primero en formalizar el método de los
elementos finitos. Él propuso que las frecuencias de las lı́neas espectrales de los átomos
podı́an ser expresadas por diferencias entre un relativamente pequeño número de ele-
mentos. Ritz desarrolló la formulación matemática del Método de los Elementos Finitos
con base en el cálculo variacional.

6. Metodologı́a de Trabajo

a) Objeto de Estudio

En la investigación se estudiará la forma variacional de los problemas de segundo

orden tipo elı́ptico y como se encuentra soluciones numéricas utilizando el método
de Ritz.

Se estudiará la implementación computacional del método de Ritz y su aplicación

en la solución de problemas de ingenierı́a.

b) Metodologı́a de la Investigación

F-M01.03.03-DRT/PG-01-(Rev.1) 11
escudo-eps-converted-to.pdf

Para el desarrollo de la investigación se seguirá el cronograma de trabajo propuesto

donde iniciamos con una revisión exhaustiva de la bibliografı́a existente del tema
propuesto. La información recolectada será organizada, analizada e interpretada con
miras a que forme parte de los antecedentes para la investigación.

7. Definición de variables de estudio

Variable independiente: Método de Ritz

Variable dependiente: Existencia y Unicidad de Problemas Variacionales de Se-

gundo Orden

F-M01.03.03-DRT/PG-01-(Rev.1) 12
Referencias

[1] Vojin J, Sergiy K. The Ritz Method for Boundary Problems with Essential Condition as
Constraints, Hindawi Publishing Corporation, Advances in Mathematical Physics, 2016,
Article ID 7058017, 12 pages.

[2] Muller J, Zeinhofer M. Error Estimates for the Deep Ritz Method with Boundary Pe-
nalty, Proceeding of Machine Learning Research, 2022; (145): 1-20

[3] Basak P.S. Calculus of Variations and Application of Rayleigh-Ritz Method in Heat and
Mass Transfer Problems, International Journal of Applied Engineering Reseach ISSN 0973-
4562 2022, Volume (17), Number 4, pp. 422-426.

[4] Yakovlev M.N. An error bound of the Ritz method for a singular second-order differential
equation, Journal of Mathematical, 2009, Sciences, Vol. 157, N o 5

[5] Akl M.Y, El-Salam F.A, Holail S. Comparation between the Rayleigh Ritz Method
(RR) and other Numerical Methods for Solving Second Order Boundary Value Problems,
IOSR Journal of Mathematics (IOSR-JM), 2019, Vol. 15, Issue 3 Ser. II, pp 47-60.

[6] Deshmukh S.V, Shukla V.V. Comparing Numerical Solution to a Second Order Boun-
dary Value Problem by Variational Techniques, International Journal of Engineering Re-
search Technology (IJERT), (2021), Vol. 10 (09).

[7] Phys J CH. Rates of convergence and error estimation formulas for the Rayleigh–Ritz
variational method, The Journal of Chemical Physics), (1985), Volume 83, Issue 3.

[8] Bougoffa L, Khanfer A. Existence and Uniqueness Theorems of Second-Order Equations

with Integral Boundary Conditions, Bull. Korean Math. Soc (55) (2018, Vol. 3, pp. 899-911.

[9] Kaufmann ER. Existence and Uniqueness of Solutions for a Second-Order Iterative
Boundary-Value Problem, Electronic Journal of Differential Equations, Vol. 2018 (2018),
Vol. 150, pp. 1-6.

[10] Ntouyas SK, Tsamatos PCh. Existence and Uniqueness of Second Order Boundary
Value Problems, Funkcialaj Ekvacioj, 38 (1995) 59-69.

13
escudo-eps-converted-to.pdf REFERENCIAS

[11] Tisdell ChC. Existence of Solutions to Second-Order Boundary Value Problems without
growth restrictions, Electronic Journal of Qualitative Theory of Differential Equations.
(2016) 92, 1-14.

UNT, 2022

F-M01.03.03-DRT/PG-01-(Rev.1) 14