UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS

(Universidad del Perú Decana de América)

FACULTAD DE CIENCIAS MATEMÁTICAS

ESCUELA PROFESIONAL DE MATEMÁTICA

“Esquema de diferencias finitas sin malla aplicado a la

solución numérica de la ecuación de onda en regiones
espaciales altamente irregulares”

AUTORES

Garay Astudillo, Gian Carlo

Echevarria Ventocilla, Alexander
Rojas Cabrera, Jhon Percy

DOCENTE

Mg. Cruz Huallpa, Alex Armando

LIMA - PERÚ
2024
Agradecimientos

Ustedes han sido siempre el motor que impulsa mis sueños y esperanzas, quienes
estuvieron siempre a mi lado en los dı́as y noches más difı́ciles durante mis horas de
estudio. Siempre han sido mis mejores guı́as de vida. Hoy cuando concluya mis
estudios, les dedico a ustedes este logro amados padres, como una meta más
conquistada. Orgulloso de haberlos elegido como mis padres y que estén a mi lado
en este momento tan importante.
”Gracias por ser quienes son y por creer en mı́”
Índice general

Agradecimientos I

Índice general II

Resumen III

Abstract IV

I Introducción 1

Introducción 1

II Esquema propuesto 4
2.1 Discretización temporal . . . . . . . . . . . . . . . . . . . . . . . . . . 4
2.2 Discretización espacial . . . . . . . . . . . . . . . . . . . . . . . . . . 6

II
Resumen

Si bien la ecuación de onda es un problema clásico para resolver en cursos sobre

ecuaciones diferenciales parciales hiperbólicas, la solución numérica de esta ecuación
en regiones irregulares ha sido poco estudiada. Esto se debe a todos los problemas que
podrı́an surgir al resolver en el tiempo la derivada de segundo orden. Este artı́culo
presenta un esquema generalizado en diferencias finitas aplicado a la solución numéri-
ca de la ecuación de onda, resuelta sobre nodos de puntos para dominios altamente
irregulares. La implementación sobre nodos de puntos permite que este esquema
calcule resultados satisfactorios, con errores cuadráticos menores, en regiones alta-
mente irregulares, como se analiza en el presente manuscrito. Además, los resultados
numéricos muestran la precisión obtenida con un método sencillo y implementación
de bajo costo.
Palabras clave: Método sin malla, Generalización de diferencias finitas, ecua-
ción de la onda, regiones irregulares, soluciones numéricas

III
Abstract

Even when the wave equation is a classical problem to solve in courses about
hyperbolic partial differential equations, the numerical solution for this equation in
irregular regions has been little studied. This is because of all the problems that
could arise when solving the second-order derivative in time. This paper presents
a generalized finite difference scheme applied to the numerical solution of the wave
equation, solved on clouds of points for highly irregular domains. The implementation
over clouds of points allows this scheme to compute satisfactory results, with minor
quadratic errors, over highly irregular regions, as discussed in the present manuscript.
Also, the numerical results show the accuracy obtained with a straightforward and
low-cost implementation
Keywords: Meshless method, Generalize finite difference, Wave equation, Irre-
gular regions, Numerical solution

IV
Capı́tulo I

Introducción

Una amplia gama de procesos fı́sicos están relacionados con la ecuación de onda.
Debido a esto, obtener soluciones para esta ecuación puede ser una cuestión crı́tica
al modelar problemas de la vida real en muchas áreas. Algunos investigadores han
propuesto varios métodos para obtener soluciones numéricas de la ecuación de la
onda. Sin embargo, aun cuando muchos de estos métodos obtienen resultados satis-
factorios, se basan en dominios unidimensionales o rectangulares o tienen un alto
coste computacional. Algunos ejemplo pueden ser encontrados en [2], [6] y [15].

En los últimos años, el diseño de esquemas de diferencias finitas ha sido objeto

de estudio de varios autores. Se han producido avances significativos realizado al
aplicar estos esquemas a soluciones numéricas de diferentes ecuaciones diferenciales
parciales [3], [4], [18] y [19]. Los métodos de diferencias finitas son la facilidad con
la que pueden implementarse y, en sus formas clásicas, el bajo coste computacional
que suponen. requieren para lograr resultados satisfactorios.

En [10], se establece un método basado en diferencias finitas para resolver la

ecuación de onda acústica. El trabajo realizado en este trabajo demuestra que es
posible crear modelos con alta precisión que puedan representar la comportamiento
de situaciones de la vida real, al menos en regiones rectangulares. Otro excelente
trabajo sobre el uso de diferencias finitas para obtener soluciones numéricas a la
ecuación de onda se presentan en [9].

Sobre él, los autores dan una nuevo esquema en diferencias finitas para la ecuación
de onda unidimensional, obtener resultados satisfactorios bajo retroalimentación de

1
Capı́tulo I. Introducción 2

lı́mites; sin embargo, el estudio se limita a regiones unidimensionales. Otro gran

ejemplo es la diferencia finita generalizada sin malla es el método presentado en [8].

Para este trabajo, los autores introdujeron un sistema estable del método sin
malla para resolver la propagación de ondas de agua no lineales y los resultados
numéricos muestran que los métodos de diferencias finitas generalizadas sin malla
pueden resolver numéricamente los problemas propuestos. La distribución de nodos
está confinada a un rectángulo. se presenta en [13] un método efectivo sin malla
para resolver la ecuación de onda elástica sobre dominios irregulares; a pesar de la
complejidad del plan y de la implementación involucrada, los resultados numéricos
muestran que este esquema puede producir resultados satisfactorios para diferentes
problemas de ondas sobre nubes de puntos irregulares.

En sentido general, la ecuación de onda modela la variación de una función escalar

u(→
−r , t), donde →
−
r representa las coordenadas espaciales, y t representa el tiempo;
que podrı́a describir la presión en lı́quidos o el desplazamiento de las partı́culas lejos
de sus puntos de reposo. Esta función escalar debe satisfacer

utt = c2 △u (1.1)

donde △ = ∇2 es el operador espacial laplaciano y c es un coeficiente real no negativo

que representa la velocidad de propagación de la onda.

Figura 1.1: Región 2D irregular

Capı́tulo I. Introducción 3

Para el caso de interés de este artı́culo, los dominios espaciales 2D (como se

muestra en la (I)), por lo tanto u(→
−
r , t) = u(x, y, t) ; esto implica que

∂ 2 u(x, y, t) ∂ 2 u(x, y, t) ∂ 2 u(x, y, t)

 
2
= c2 + (1.2)
∂t ∂x2 ∂y 2

Siguiendo la idea presentada en [5], el método presentado en este artı́culo aplica

un esquema generalizado de diferencias finitas para obtener soluciones de valores
numéricos de problemas de ecuaciones de onda en este tipo de regiones.

Figura 1.2: Distribución arbitraria de nodos

Figura 1.3: Selección de nodos vecinos

Capı́tulo II

Esquema propuesto

Este artı́culo tiene como objetivo obtener un esquema en diferencias finitas para
la ecuación de la onda.

∂ 2u
 2
∂ u ∂ 2u

2
=c + , Ω × [0, T ], c ∈ R (2.1)
∂t2 ∂x2 ∂y 2

con las condiciones iniciales:

u(x, y, 0) = h(x, y, 0), (x, y) ∈ Ω.

u(x, y, t)|∂Ω = h(x, y, t), (x, y) ∈ Ω, t ∈ [0, T ].

∂(x, y, 0)
= g(x, y), (x, y) ∈ Ω. (2.2)
∂t
donde Ω es un dominio plano simplemente conexo, su frontera ∂Ω es un polı́gono de
Jordan, orientado positivamente. Dado que esta ecuación presenta algunas compli-
caciones en la discretización, una buena idea es discretizar el tiempo y el espacio por
separado. En las siguientes secciones se presentan ambas discretizaciones.

2.1. Discretización temporal

La discretización del tiempo es crucial en esta ecuación ya que se tiene una deri-
vada del tiempo de segundo orden. Para empezar, puede ser una buena idea expresar
derivada parcial con respecto al tiempo como una aproximación en diferencias finitas.

4
Capı́tulo II. Esquema propuesto 5

Esto es,
∂ 2 u(x, y, t) uk+1
i − 2uki + uk−1
i
=
∂t2 (△t)2
Entonces, para un nodo arbitrario p0 = (x0 , y0 ), la ecuación de onda se puede escribir
como:
uk+1 − 2uk0 + uk−1
 2
∂ u ∂ 2u

0 0 2
=c +
(△t)2 ∂x2 ∂y 2
Luego, resolviendo para uk+1
0 , obtenemos:

∂ 2u ∂ 2u
 
uk+1
0 = 2uk0 − uk−1
0 + (c△t) 2
+ (2.3)
∂x2 ∂y 2

Se puede notar que esta aproximación requiere pasos de 2 veces para ser calculado.
Entonces, es esencial aproximar el segundo paso de tiempo. adecuadamente. En [11],
se afirma que es posible considerar la condición para la primera derivada parcial en
el tiempo.
∂(x, y, 0)
= g(x, y)
∂t
Donde, si aplicamos el método centrado de diferencias finitas, toma la forma:

u10 − u−1
0
= g(x, y)
2△t

Luego, resolviendo para uk−1

0 , obtenemos:

u−1 1
0 = u0 − 2△tg(x, y) (2.4)

Suponiendo que (2.3) se cumple en k = 0, es decir:

∂ 2u ∂ 2u
 
u10 = 2u00 − u−1
0 + (c△t) 2
+ (2.5)
∂x2 ∂y 2
t=0

Si sustituimos (2.4) en (2.5), obtenemos:

(c△t)2 ∂ 2u ∂ 2u
 
u10 = u00 + △tg(x, y) + + (2.6)
2 ∂x2 ∂y 2
t=0
Capı́tulo II. Esquema propuesto 6

Ahora es posible utilizar la discretización temporal en (2.6) para calcular el segundo

paso de tiempo, y el presentado en (2.3) para calcular todos los siguientes pasos de
tiempo.

2.2. Discretización espacial

Una vez que la ecuación ha sido parcialmente discretizada en el tiempo, realiza-

remos los siguientes pasos:

El primer paso es obtener el esquema de diferencias finitas que se va a aplicar.

Una forma sencilla de lograr esto es considerar la operador lineal de segundo orden.

Lu = Auxx + Buxy + Cuyy + Dux + Euy + F u

considerando A, B, C, D, E y F como funciones dadas.

Dentro de una distribución arbitraria de nodos, como la que se presenta en la I,

es posible aproximar su valor en un nodo central po = (x0 , y0 ) usando valores de u
en algunos nodos vecinos pi = (xi , yi ), i = 1, 2, ..., q. Se puede aplicar un esquema en
diferencias finitas en p0 como una combinación lineal:
q
X
L0 u = Γ0 u(p0 ) + Γ1 u(p1 ) + Γ2 u(p2 ) + . . . + Γq u(pq ) = Γi u(pi )
i=o

donde Γ0 , Γ1 , Γ2 , . . . , Γn son pesos adecuados. Según [14] y [16], un esquema en di-

ferencias finitas p0 es consistente con el operador lineal L0 si el error de truncamiento
local τ satisface:
τ = [Lu]p0 − [L0 u]p0 → 0 (2.7)

donde p1 , p2 , . . . , pq → 0 Expandiendo en series de Taylor, hasta segundo orden,

Capı́tulo II. Esquema propuesto 7

esto se puede reescribir como:

q q
! !
X X
[Lu]p0 − [L0 u]p0 = F (p0 ) − Γi u(p0 ) + D(p0 ) − Γi △xi ux (p0 )
i=0 i=0
q
!
X
+ E(p0 ) − Γi △yi uy (p0 )
i=0
q
!
X Γi (△xi )2
+ A(p0 ) − uxx (p0 )
i=0
2
q
!
X
+ B(p0 ) − Γi △xi △yi uxy (p0 )
i=0
q
!
X Γi (△yi )2
+ C(p0 ) − uyy (p0 ) + ϑ(máx{△xi , △yi })3
i=0
2

donde △xi = xi − x0 y △yi = yi − y0 . Para lograr la consistencia de la condición

(2.7), cada uno de los términos entre paréntesis debe desaparecer simultáneamente;
esto es q q
X Γi (△xi )2 X
A(p0 ) − =0 D(p0 ) − Γi △xi = 0
i=0
2 i=0
q q
X X
B(p0 ) − Γi △xi △yi = 0 E(p0 ) − Γi △yi = 0
i=0 i=0
q q
X Γi (△yi )2 X
C(p0 ) − =0 F (p0 ) − Γi = 0
i=0
2 i=0

Estas condiciones definen el sistema lineal:

 
  Γ0  
1 1 ... 1 F (p0 )
 Γ1  
 
 
0
 △x1 ... △xq   Γ2   D(p0 ) 
   
0 △x1 ... △yq     E(p0 ) 
   
 .  =  (2.8)
 
0 (△xi )2 2 
 
... (△xq )    2A(p0 )

0 △x △y . 
  
i i ... △xq △yq     B(p0 ) 
 
.
 
2 2
0 (△yi ) ... (△yq ) 2C(p0 )
Γq

Para resolver este sistema lineal, es posible separar la primera ecuación del sistema
Capı́tulo II. Esquema propuesto 8

(2.8).
q
X
Γi − F (p0 ) = 0 (2.9)
i=0

Luego, el problema esta definido como

 
  Γ1  
△x1 ... △xq   D(p0 )
 Γ2 
 △x1 ... △yq    E(p0 ) 
   
. 
   
 (△xi )2 ... (△xq )2    = 2A(p )
0 
 (2.10)
  .  
△xi △yi ... △xq △yq  
 .   B(p0 ) 
    

(△yi )2 ... (△yq )2 2C(p0 )

 
Γq

Se puede resolver utilizando la factorización reducida de Cholesky de su normal

ecuaciones para evitar la eficiencia de rango de fila, como en [17], siendo

MT MΓ = MT β

donde  
  Γ1  
△x1 ... △xq   D(p0 )
Γ2 
 △x1 ... △yq   E(p0 ) 
     
  .  
2
M =
 (△xi ) ... (△xq )2 
, Γ =  , β= 2A(p )
  
0 
. 
△xi △yi ... △xq △yq   B(p0 ) 
     
.
(△yi )2 ... (△yq )2 2C(p0 )
 
Γq
El valor de Γ0 se obtiene usando (2.9) suponiendo para el caso de la ecuación de onda,
F (p0 ) = 0. Ahora, el esquema definido por (2.10) se puede utilizar para aproximar
el operador lineal  2
∂ 2u

2 ∂ u
Lu = (c△t) +
∂x2 ∂y 2
en cada nivel de tiempo, tomando A = C = (c△t)2 , y B = D = E = F = 0. Los
coeficientes Γi resultantes, con las discretizaciones temporales (2.3) y (2.6) definen
Capı́tulo II. Esquema propuesto 9

los esquemas generalizados de diferencias finitas

q
1X
u10 = u00 + △tg(x, y) + Γi u0i (2.11)
2 i=0

q
X
uk+1
0 = 2uk0 − u0k−1 + △tg(x, y) + Γi uki (2.12)
i=0

para la ecuación de onda, donde k representa el nivel de tiempo y uki es el aproxi-

mación a la solución en el punto pi = (xi , yi ) en el instante k△t.
Una cuestión crı́tica a considerar es el número de vecinos, q, a usar en el esquema.
En este trabajo, el número de vecinos variará de nodo a nodo. Para cada uno, los
vecinos se seleccionan usando la distancia con el nodo central considerando todos los
nodos en un espacio menor o igual a δ donde δ es la distancia promedio entre los
nodos a lo largo del lı́mite (I).
Amplias pruebas han demostrado que, con este criterio para seleccionar el nodos
vecinos, se considera un promedio de 6 vecinos por cada nube nodo, con un mı́nimo
de 4 y un máximo de 9, en todos los casos. Esta es una reducción significativa en el
número de nodos en comparación con otros esquemas en [8] y [13].
Bibliografı́a

[1] M.B. Abbott, A.D. McCowan, I.R. Warren, Accuracy of short-wave numerical
models, J. Hydraul. Eng. 110 (10) (1984) 1287–1301.

[2] H. Ahmad, A.R. Seadawy, T.A. Khan, P. Thounthong, Analytic approximate

solutions for some nonlinear parabolic dynamical wave equations, J. Taibah
Univ. Sci. 14 (1) (2020) 346–358.

[3] P. Barrera-Sánchez, L. Castellanos-Noda, F.J. Domı́nguez-Mota, G. González-

Flores, A. Pérez-Domı́nguez, Adaptive discrete harmonic grid generation, Math.
Comput. Simul. 79 (6) (2009) 1792–1809.

[4] F.J. Domı́nguez-Mota, S. Mendoza-Armenta, G. Tinoco-Guerrero, J.G. Tinoco-

Ruiz, Numerical solution of Poisson-like equations with Robin boundary con-
ditions using a finite difference scheme defined by an optimality condition, in:
IMACS Series in Computational and Applied Mathematics: MASCOT 2011,
2011.

[5] F.J. Domı́nguez-Mota, J.S. Lucas-Martı́nez, G. Tinoco-Guerrero, A generalized

finite differences scheme used in modeling of a direct and an inverse problem of
advection-diffusion, Int. J. Appl. Math. 33 (4) (2020) 599–608.

[6] L. Feng, F. Liu, I. Turner, Finite difference/finite element method for a novel
2d multi-term time-fractional mixed sub-diffusion and diffusion-wave equation
on convex domains, Commun. Nonlinear Sci. Numer. Simul. 70 (2019) 354–371.

[7] S. Gottlieb, C.-W. Shu, Total variation dimishing Runge-Kutta schemes, Math.
Comput. 67 (221) (1998) 73–85.

10
Bibliografı́a 11

[8] J. Huang, C.-M. Fan, J.-H. Chen, J. Yan, Meshless generalized finite differen-
ce method for the propagation of nonlinear water waves under complex wave
conditions, Mathematics 10 (6) (2022).

[9] J. Liu, B.-Z. Guo, A new semidiscretized order reduction finite difference scheme
for uniform approximation of one-dimensional wave equation, SIAM J. Control
Optim. 58 (4) (2020) 2256–2287.

[10] Y. Liu, M.K. Sen, A new time-space domain high-order finite-difference method
for acoustic wave equation, J. Comput. Phys. 228 (23) (2009) 8779–8806.

[11] A. Peirce, in: Solving the Heat, Laplace and Wave Equations Using Finite Dif-
ference Methods, 2018, Chap. 8.

[12] P.-O. Persson, G. Strang, A simple mesh generator in MATLAB, SIAM Rev. 46
(2)(2004) 329–345.

[13] E. Salete, A. Vargas, A. Garcı́a, J. Benito, F. Ureña, M. Ureña, An effective nu-

meric method for different formulation of the elastic wave propagation problem
in isotropic medium, Appl. Math. Model. 96 (2021) 480–496.

[14] J.C. Strikwerda, Finite Difference Schemes and Partial Differential Equations,
Society for Industrial and Applied Mathematics, 2004.

[15] N. Sweilam, M. Khader, A. Nagy, Numerical solution of two-sided space-

fractional wave equation using finite difference method, J. Comput. Appl. Math.
235 (8) (2011) 2832–2841.

[16] J.W. Thomas, Numerical Partial Differential Equations: Finite Difference

Methods, Texts in Applied Mathematics, vol. 2, Springer, 1995.

[17] G. Tinoco-Guerrero, F.J. Domı́nguez-Mota, J.G. Tinoco-Ruiz, A study of the

stability for a generalized finite-difference scheme applied to the advection-
diffusion equation, Math. Comput. Simul. 176 (2020) 301–311.

[18] J.G. Tinoco-Ruiz, F.J. Domı́nguez-Mota, S. Mendoza-Armenta, G. Tinoco-

Guerrero, Solving Stokes equation in plane irregular region using an optimal
Bibliografı́a 12

consistent finite difference scheme, in: 20th International Meshing Roundtable,

2011.

[19] J.G. Tinoco-Ruiz, F.J. Domı́nguez-Mota, G. Tinoco-Guerrero, P.M. Fernández

Valdés, S. Mendoza-Armenta, Numerical solution of differential equations in
irregular plane regions using quality structured convex grids, Int. J. Model.
Simul. Sci. Comput. 4 (2) (2013).