7_geos_31_2_elementos
TRANSCRIPT
-
7/21/2019 7_GEOS_31_2_Elementos
1/18
Elementos Finitos de Orden Superior para Propagacin de OndasSsmicas
Jons D. De Basabe
Divisin de Ciencias de la Tierra
CICESE
RESUMEN
En las lmas dcadas se han realizado grandes avances en los mtodos de elementos
nitos, tanto prccos como analcos. Estos avances son relevantes para sus aplicaciones
en sismologa. Dos mtodos en parcular que han atraido la atencin de invesgadores en
el rea de sismologa numrica y computacional son: el de elementos espectrales y el de
Galerkin disconnuo, por las ventajas que ofrecen para discrezar un medio heterogneo
y obtener sismogramas sintcos con un orden elevado de precisin. En este arculo sedescriben estos mtodos, se presentan los princiales resultados del anlisis de precisin
y estabilidad numrica y se muestran sismogramas sintcos para el problema de Lamb.
Con base en los resultados analcos y numricos, concluimos que los mtodos de orden
4 superior enen mnima dispersin y anisotropa numrica, por lo cual son ideales para
propagaciones numricas de gran escala. La razn de muestreo puede ser de hasta 4 5
nodos por longitud de onda.
INTRODUCCION
En las lmas cinco dcadas se ha
invesgado intensamente la aplicacin
de disntos mtodos numricos para
la propagacin de ondas ssmicas. La
principal movacin de esto es el hecho
de que no existen soluciones analcas
para modelos geosicos de inters, por
ejemplo, modelos con una distribucin
arbitraria de heterogeneidades o con
topograa. La mayora de los mtodos
numricos para propagacin de ondas que
se han ulizado hasta ahora se puedenclasicar en alguna de las siguientes
categoras:
1. Los que ulizan alguna aproximacin
asintca de la ecuacin de onda,
como el mtodo de trazado de rayos
(Cerven, 2001; Madrid, 2008),
2. Los que aproximan las derivadas
espaciales en la ecuacin de onda a
travs de la transformada de Fourier,
como el mtodo pseudoespectral
(Koslo y Baysal, 1982; Fornberg,
1987),
3. Los que evalan las derivadas
espaciales y temporales de la
ecuacin de onda a travs de series
de Taylor truncadas, como el mtodo
de diferencias nitas con mallado
simple (Alterman y Karal, 1968;
Alford et al., 1974; Kelly et al., 1976)
o el de diferencias nitas con mallado
intercalado (Madariaga, 1976; Virieux,
1984, 1986), y
GEOS, Vol. 31, No. 2
228
-
7/21/2019 7_GEOS_31_2_Elementos
2/18
4. Los que discrezan alguna forma
dbil de la ecuacin de onda, como el
mtodo de elementos nitos (FEM, por
sus siglas en ingls) (Lysmer y Drake,
1972; Smith, 1975; Marfurt, 1984) o el
de volumen nito (Dormy y Tarantola,
1995; Dumbser et al., 2007a).
Describir los mritos y limitaciones
de cada uno de estos mtodos queda
fuera del alcance de los objevos de
este arculo. Para una descripcin de
todos estos mtodos, ver Carcione et
al. (2002). Basta decir que los mtodos
basados en la ecuacin de onda, como elpseudoespectral, el de diferencias nitas
y el de elementos nitos, enen la ventaja
de que incluyen en las simulaciones las
ondas directas, reejadas (incluyendo
reexiones mlples), transmidas,
converdas, refractadas, difractadas y
superciales. Asimismo, los mtodos que
ulizan directamente la ecuacin de onda,
como los de diferencias nitas y elementos
nitos, ofrecen mayor exibilidad para
incluir heterogeneidades en el modeloy para incorporar cualquier po de
condiciones de frontera. En parcular,
el FEM ofrece la mayor exibilidad para
ulizar discrezaciones que se apegan a la
topograa y a las disconnuidades de los
parmetros del modelo sico.
En este arculo enfocaremos nuestra
atencin en los FEM, y en parcular en los
de orden superior. Aunque los FEM han
sido ulizados en sismologa desde nalesde los aos sesenta (Chopra et al., 1969),
la popularidad de este po de mtodos fue
en un principio muy limitada debido a que
las versiones propuestas (Lysmer y Drake,
1972; Smith, 1975) tenan una precisin
limitada y consuman considerablemente
ms recursos computacionales que las
implementaciones de diferencias nitas.
La precisin de estos mtodos ha sido
estudiada en Mullen y Belytschko (1982)
y en Marfurt (1984) en funcin de su
dispersin numrica. Mullen y Belytschko
(1982) estudiaron la dispersin numrica
del FEM con diferentes pos de elementos
y concluyeron que los rectangulares
introducen un mnimo de dispersin. Por
su parte, Marfurt (1984) compar el FEM
con el de diferencias nitas y concluy
que el error por dispersin numrica
de ambos mtodos es comparable. Sin
embargo, en el FEM es necesario obtener
la inversa de una matriz, conocida comola matriz de masa, en cada paso de las
iteraciones, para avanzar el sistema en la
dimensin temporal. Este autor propuso
diferentes estrategias para superar estas
limitaciones, por ejemplo resolver el
sistema en el domino de la frecuencia,
ulizar una matriz de masa diagonalizada
o una combinacin de las matrices de
masa consistente y diagonalizada, pero los
resultados no fueron alentadores. En una
connuacin de este anlisis (Marfurt,1990) tambin consider el efecto en
la dispersin numrica que resulta de
ulizar polinomios de orden superior
para las aproximaciones y concluy que
esto introducira ondas espurias que
contaminaran aun ms los resultados. Por
otro lado, la popularidad de los mtodos
de diferencias nitas se incrementaba
por la llegada del mtodo de mallado
intercalado (Virieux, 1984, 1986; Levander,
1988), el cual reduce drscamente ladispersin numrica (Moczo et al., 2000).
A pesar de todo esto, la exibilidad del FEM
para discrezar geometras complicadas
sigui atrayendo a invesgadores en el
rea de sismologa de terremotos. Por
ejemplo, Bao et al. (1998) realizaron una
De Basabe: Elementos nitos de orden superior para propagacin de ondas ssmicas
229
-
7/21/2019 7_GEOS_31_2_Elementos
3/18
implementacin a gran escala de este
mtodo para simular eventos ssmicos
en el sur de California y evaluar el riesgo
en reas urbanas. En su implementacin
usaron una matriz de masa diagonalizada
para evadir el costo computacional de
calcular su inversa, y controlaron el
error asociado a la dispersin numrica
ulizando entre 10 y 20 elementos por
longitud de onda.
En la dcada de los noventa hubo avances
importantes en el desarrollo del FEM
que permieron que el mtodo fuera de
nuevo compevo con el ubicuo mtodo
de diferencias nitas, superandolo
en precisin y con una eciencia
computacional comparable. Dichos
avances se presentaron en dos frentes;
por un lado, se desarrollaron tcnicas para
diagonalizar la matriz de masa, conocidas
como tcnicas de acumulacin de masa,
que no comprometen la precisin (Cohen,
2002). Por otro lado, se empezaron a
ulizar polinomios de orden superior,
inspirados en el mtodo pseudoespectral,
para reducir la dispersin y anisotropa
numrica (Cohen, 2002). Las tcnicas de
acumulacin de masa estn basadas en el
uso de los mismos nodos para la cuadratura
numrica y para los nodos que denen los
polinomios de aproximacin dentro de
los elementos; los ms comunes son los
de Gauss-Lobao-Legendre (GLL), pero
existen otras opciones dependiendo del
mtodo parcular de elementos nitos.
Como se mencion anteriormente, el uso
de polinomios de orden superior goz de
escasa popularidad, debido en parte a que
se crea que introducan errores numricos,
al basarse en el anlisis preliminar de
Marfurt (1990). Adicionalmente, hay
un costo computacional considerable
siempre que la proporcin de elementos
diferentes de cero en la matriz de masa es
proporcional al orden de los polinomios,
de tal manera que el costo computacional
para almacenar e inverr esta matriz
aumenta rpidamente. Estas limitaciones
fueron superadas ulizando los nodos
y tcnicas de cuadratura mencionadas
arriba; asimismo, se ha demostrado
que, al contrario de lo que sealaban los
anlisis preliminares, los polinomios de
orden superior reducen la dispersin y
anisotropa numrica (ver la seccin 4
para ms detalles).
En este arculo enfocaremos nuestra
atencin en dos mtodos que se han
vuelto muy populares en la literatura
y que combinan las caracterscas
mencionadas anteriormente, es decir,
que ulizan polinomios de orden superior
y una matriz de masa diagonal. Estos
mtodos son el de elementos espectrales
y el de Galerkin disconnuo. En las
siguientes dos secciones se describen
cada uno de estos mtodos. En la seccin
4 se describe su dispersin y estabilidad
numrica, y posteriormente se comparan
los resultados de aplicar estos mtodos al
problema de Lamb (Lamb, 1904).
ELEMENTOS ESPECTRALES
El Mtodo de Elementos Espectrales
(SEM, por sus siglas en Ingls),
originalmente desarrollado para dinmica
de uidos (Patera, 1984), ha sido aplicado
exitosamente para la ecuacin de onda.
Ha logrado superar las limitaciones deeciencia, ofreciendo una mayor precisin
que el mtodo de diferencias nitas y
una mayor exibilidad geomtrica. Este
mtodo fue ulizado originalmente para
la ecuacin de onda acsca por Seriani y
Priolo (1994), quienes ulizaron los nodos
GEOS, Vol. 31, No. 2
230
-
7/21/2019 7_GEOS_31_2_Elementos
4/18
de Gauss-Lobao-Chebyshev (GLC), y por
Tordjman et al. (1994) y Tordjman(1995),
quienes ulizaron los nodos de Gauss-
Lobao-Legendre (GLL). El mtodo fue
ulizado para discrezar la ecuacin de
onda elsca por Komatsch y Viloe
(1998), con lo cual se increment su
popularidad rpidamente dentro de la
comunidad sismolgica (Komatsch y
Tromp, 1999, 2002a, b; Priolo, 2001;
Komatsch et al., 2002, 2004, 2005;
Madec et al., 2009; Chaljub et al., 2007).
Este mtodo tambin se ha ulizado para
problemas en sismologa global y regional,
arrojando resultados notablemente
cercanos a los datos empricos (Komatschet al., 2002, 2004).
La principal caractersca de SEM es que
uliza polinomios de orden superior
para representar el campo de onda, con
la ventaja de que se puede ulizar una
razn de muestreo muy baja de entre 4 y
5 nodos por longitud de onda. La matriz
de masa se puede diagonalizar usando los
nodos de GLL y las reglas de cuadratura
correspondientes. Esta estrategia es
equivalente a la tcnica de concentracin
de masa, pero con la ventaja de que no se
sacrica la precisin.
Este mtodo, al igual que la versin clsica
del FEM, est basado en la forma dbil de
la ecuacin de onda, la cual se desarrolla
a connuacin. El modelo sico para
la propagacin de ondas elscas est
basado en la ecuacin del movimiento,
dada por (usando la notacin de Einstein)
Aki y Richards (2002)
en
(1)
donde d son las dimensiones (2 3),
es el dominio en el espacio,
(0,T] es el dominio en el empo,, es el
vector de desplazamiento, es el tensor
de esfuerzo, y la fuente est dada por el
tensor de momento
y el vector de
fuerzas.En la ecuacin anterior se uliza
la notacin abreviada .El
tensor de esfuerzo se puede escribir en
funcin del desplazamiento por medio de
la ley de Hooke, la cual est dada en un
medio elsco e isotrpico por
(2)
donde l y m son los parmetros de
Lam, los cuales pueden ser funciones
de las variables espaciales en un medio
heterogneo, y
es la delta de Kronecker.Sustuyendo la ley de Hooke en la ecuacin
de movimiento se obene la ecuacin de
onda elsca:
(3)
donde , , ,
es la frontera de ,
y
es unvector unitario ortogonal a . La forma
dbil de la ecuacin de onda elsca se
obene al mulplicar la ecuacin (3) por
un vector arbitrario v, integrar sobre Wy
aplicar el teorema de Gauss, para obtener
(4)
donde
Observar que el espacio vectorial es un
espacio de Sobolev de primer orden en el
cual se encuentra la solucin dbil. Para
De Basabe: Elementos nitos de orden superior para propagacin de ondas ssmicas
231
-
7/21/2019 7_GEOS_31_2_Elementos
5/18
discrezar esta ecuacin en 2D es necesario
introducir un subespacio de de
dimensiones nitas
donde
. Las funciones
, conocidas como funciones
base, se denen ulizando la parcin
del dominio W en subdominios,
llamados elementos, y los nodos dentro
de los elementos. Sustuyendo u por
la aproximacin
dada por la
siguiente combinacin lineal
,, , ,
,
(5)
donde
y
son los coecientes
de las aproximaciones de las componenteshorizontal y vercal del vector de
desplazamientos, y sustuyendo
se obene el siguiente
sistema de ecuaciones diferenciales
ordinarias
. (6)
De manera similar, sustuyendo
se obene
(7)
donde las matrices en las ecuaciones (6) y
(7) estn dadas por
(8)
(9)
(10)
(11)
(12)
(13)
(14)
y y son las componentes en x y
z de f. Los coecientes
y
y se
pueden obtener de las ecuaciones (6) y
(7) ulizando un mtodo para resolver
ecuaciones diferenciales ordinarias,
como los de diferencias nitas, Newmark,
Runge-Kua o Lax-Wendro.
Galerkin Discontinuo
El Mtodo de Galerkin Disconnuo (DGM,
por sus siglas en Ingls) fue originalmente
propuesto por Reed y Hill (1973) para
solucionar la ecuacin hiperblica de
transporte de neutrones. La estrategiageneral para resolver ecuaciones
diferenciales hiperblicas ulizando la
forma de primer orden con el DGM fue
desarrollada en una serie de arculos
por Cockburn y sus colaboradores, en
Cockburn y Shu (1989); Cockburn et al.
(1989, 1990) y Cockburn y (1998). En
estos arculos, ellos proponen discrezar
la ecuacin diferencial en el espacio con
una formulacin basada en ujos y en el
empo usando mtodos de Runge-Kua.Como regla general, el orden del mtodo
de Runge-Kua debe de ser el mismo que
el de la discrezacin en el espacio.
En el contexto de propagacin de ondas
ssmicas, recientemente se ha propuesto
sustuir el mtodo de Runge-Kua por la
estretegia ADER para avanzar el sistema
en el empo en: Dumbser y Munz (2006);
Kser y Dumbser (2006); Dumbser y Kser
(2006); Kaser et al. (2007a,b); de la Puenteet al. (2007); Dumbser et al. (2007a, b).
Esta estrategia est basada en series de
Taylor en la variable temporal, las cuales
estn truncadas para obtener la misma
precisin en empo que en el espacio. Las
derivadas de orden superior se sustuyen
GEOS, Vol. 31, No. 2
232
-
7/21/2019 7_GEOS_31_2_Elementos
6/18
por derivadas espaciales ulizando el
teorema de Cauchy-Kovalewski. A este
mtodo tambin se le conoce como el
mtodo de Lax-Wendro (Lax y Wendro,
1964). Los mtodos mencionados arriba
ulizan una forma de primer orden de la
ecuacin de onda en velocidad y esfuerzo.
Otros mtodos de Galerkin disconnuo
que se han propuesto para propagacin
de ondas son el simtrico de penalizacin
interior (Grote et al., 2006) y el no-
simtrico de penalizacin interior (Riviere
y Wheeler, 2001).
En este arculo enfocaremos nuestra
atencin en el DGM de penalizacin
interior (IP-DGM), porque ene la ventaja
de que uliza menos grados de libertad
en cada nodo, pues se puede ulizar
con la ecuacin de onda en funcin del
desplazamiento (ecuacin (3)). Este
mtodo est basado en la siguiente forma
dbil de la ecuacin de onda:
Encontrar
{ }
tal que, para todo
,
(15)
donde
(16)
(17)
() *()+,-
(18) *()+,- * +,-,-
y es el vector de traccin, el cual est
dado en un medio isotrpico por
(19)
El parmetro R en la ecuacin (18) es la
penalizacin y Ses un parmetro que tomalos valores de 0, 1 -1 dependiendo de la
forma parcular del IP-DGM: S = 0 para
el Galerkin incompleto con penalizacin
interior (IIPG) (Dawson et al., 2004), S=-1
para el Galerkin simtrico con penalizacin
interior (SIPG) (Darlow, 1980; Grote et al.,
2006) y S =1 para el Galerkin no simtrico
con penalizacin interior (NIPG) (Riviere
et al., 1999, 2001).
De nuevo, como en el caso de SEM,introducimos un espacio nito
, y sustuimos en
la ecuacin (15) las funciones u y v por
combinaciones lineales de las funciones
base del espacio nito para obtener
un sistema de ecuaciones diferenciales
ordinarias como el de las ecuaciones (6) y
(7), pero con las siguientes matrices:
(20)
(21)
(22)
(23)
(24)
(25)
(26)
De Basabe: Elementos nitos de orden superior para propagacin de ondas ssmicas
233
-
7/21/2019 7_GEOS_31_2_Elementos
7/18
Dispersin y Estabilidad Numrica
Los anlisis de dispersin y estabilidad
numrica son herramientas fundamentales
que evalan la aplicabilidad de los mtodosnumricos para propagacin de ondas y
que ayudan a determinar los parmetros
necesarios para realizar simulaciones,
especcamente, incrementos en espacio
y empo. Al analizar la dispersin
numrica se busca obtener un criterio
para determinar el tamao de los
incrementos en el espacio de tal manera
que el error sea mnimo, y al analizar la
estabilidad se busca obtener un rango
para el incremento en el empo dentrodel cual el mtodo sea estable, es decir,
dentro del cual la magnitud del campo de
onda est acotada.
La dispersin numrica es un po de error
que se comporta de una manera similar
a la dispersin sica. Este fenmeno se
presenta debido a que las ondas con altas
frecuencias se transportan a una velocidad
diferente que las de bajas frecuencias,
por lo que el frente de onda se dispersa.En la prcca, los resultados numricos
pueden parecer cualitavamente
correctos, pero los empos de arribo
pueden estar adelantados o atrasados, y
las amplitudes y formas de onda tambin
pueden estar distorsionadas. Para evitar
que los resultados se vean afectados por
este po de error es necesario determinar
cuidadosamente el nmero de nodos por
longitud de onda, de tal manera que la
dispersin sea mnima para las frecuenciasde inters.
En cuanto a la estabilidad numrica, es
importante resaltar que la gran mayora de
los mtodos que se ulizan en Sismologa
para discrezar las derivadas en el empo
son condicionalmente estables. Esto
signica que existe un rango restringido
dentro del cual se puede seleccionar el
incremento en el empo. De otra manera
las amplitudes del campo de onda se
incrementaran indenidamente. Dado
que por lo general las cotas impuestas
por los criterios de estabilidad son ms
restricvas que las impuestas por criterios
de precisin, normalmente se escogen los
incrementos en el empo ms grandes
posibles que permite el criterio de
estabilidad. Las condiciones de estabilidad
se expresan por lo general como una
condicin de Courant-Friedrichs-Lewy
(CFL), la cual es una desigualdad que
acota el incremento en el empo conuna constante por el incremento en el
espacio dividido por la velocidad, donde la
constante depende del mtodo numrico
en cuesn (Courant et al., 1967).
La precisin de SEM para el caso acsco
ha sido estudiada por varios autores, por
ejemplo Seriani y Priolo (1994), Mulder
(1999), Cohen (2002), Ainsworth (2004a)
y De Basabe y Sen (2007). La principal
conclusin a la que llegaron estos autores
es que el mtodo arroja una precisin
aceptable ulizando nicamente 4.5
nodos por longitud de onda en promedio
y polinomios de orden cuatro o superior.
Esta es una gran ventaja si lo comparamos,
por ejemplo, con el mtodo de orden 1
con los mtodos de diferencias nitas con
mallado intercalado, que necesitan por lo
menos 10 nodos por longitud de onda, o
con el de diferencias nitas con mallado
estndar, que requiere ms de 15 nodos
por longitud de onda. Adems, en todos
estos arculos se enfaza la ventaja de
ulizar polinomios de orden superior, ya
que enen menos dispersin y anisotropa
numrica.
GEOS, Vol. 31, No. 2
234
-
7/21/2019 7_GEOS_31_2_Elementos
8/18
Por otro lado, el caso elsco ha sido
estudiado por De Basabe y Sen (2007).
Las conclusiones de estos autores fueron
consistentes con las mencionadas para
el caso acsco. Adems, observaron
que este mtodo es superconvergente
al incrementar el orden polinomial de
las funciones base y que la anisotropa
numrica es nula al usar polinomios
de orden tres o superior. En la gura
1 se encuentran gracados los errores
atribuibles a la dispersin numrica de
SEM contra el logaritmo de la razn de
muestreo d ulizando polinomios de
orden 1 hasta 4 (k=1,...,4). La pendiente de
las rectas indica el orden de convergenciadel mtodo. De esta grca es posible
concluir que el orden de convergencia de
SEM es aproximadamente 2k. Dado que el
orden de convergencia esperada es k, se
dice que este mtodo ene la propiedad
de superconvergencia. La gura 2 muestra
la anisotropa numrica de SEM de orden
2 a 4, donde se encuentra gracada la
dispersin en funcin del ngulo entre
la malla y el ngulo de incidencia. Para
poder visualizar la anisotropa numrica,en la grca se uliz una razn de las
velocidades de la onda P sobre la onda S
de a/b= 10 (donde ay b son las velocidad
de las ondas P y S respecvamente), una
razn de muestreo de d=0.2 (5 nodos por
longitud de onda), y una propagacin de
200 longitudes de onda. En esta grca,
la dispersin est normalizada, por lo
tanto el crculo con radio 1 representa la
ausencia de dispersin; tal es el caso de
los mtodos de orden 3 y 4. El mtodo desegundo orden muestra una signicava
anisotropa numrica asociada a la
dispersin, con un error de hasta 50% a los
30 y 60, mientras que para las ondas que
se propagan paralelas a los ejes (0 y 90)
la dispersin es mnima. Para ms detalles
sobre las guras 1 y 2 ver De Basabe et
al.(2008).
La estabilidad de SEM en los casos
acsco y elsco ulizando el mtodo
clsico de diferencias nitas y los deLax-Wendro de orden superior ha sido
estudiada por De Basabe y Sen (2010).
En ese arculo se uliz el mtodo de
Euler para obtener condiciones de CFL,
y se muestran condiciones necesarias y
sucientes de estabilidad hasta orden 10
en espacio y empo. Como se mencion
anteriormente, las condiciones de CFL
se expresan como una desigualdad de la
forma
(27)
donde Dtes el incremento en el empo,
a es la velocidad de la onda P, h es el
tamao de los elementos, y C es una
constante que depende nicamente de
los mtodos numricos en cuesn. Los
valores de C para SEM de orden 1 a 10
ulizando el mtodo clsico de diferencias
nitas para discrezar las derivadas
temporales se muestran en la tabla 1 parael caso acsco y en la tabla 3 para el caso
elsco. Muchas veces las condiciones
de estabilidad se muestran ulizando el
incremento mnimo en el espacio Dx en
lugar del tamao de los elementosh, de la
siguiente manera:
(28)
Las constantes C de la ecuacin (28) se
muestran en las tablas 2 y 4 para los casos
acsco y elsco respecvamente. Al
ulizar el incremento mnimo en lugar del
tamao de los elementos se observa que
las constantes se acercan asintcamente
a los valores de 0.607 y 0.7 para los casos
acsco y elsco respecvamente.
De Basabe: Elementos nitos de orden superior para propagacin de ondas ssmicas
235
-
7/21/2019 7_GEOS_31_2_Elementos
9/18
Figura 1. Error por dispersin nmerica de SEM
de orden 1 a 4 en el caso elsco (De Basabe et
al., 2008)
Figura 2. Anisotropa numrica de SEM de orden 2
a 4 en el caso elsco para una propagacin de 200
longitudes de onda (De Basabe et al., 2008)
Figura 3. Error por dispersin numrica de la ver-sin SIPG de DGM de orden 1 a 4 en el caso els-
co (De Basabe et al., 2008).
Figure 4. Anosotropa numrica de la versinSIPG de DGM de orden 2 a 4 en el caso elsco
para una propagacin de 200 longitudes de onda
(De Basabe et al., 2008)
GEOS, Vol. 31, No. 2
236
-
7/21/2019 7_GEOS_31_2_Elementos
10/18
Orden SEM SIPG IIPG NIPG
1 0.709 0.408 0.408 0.353
2 0.288 0.182 0.166 0.152
3 0.164 0.108 0.0934 0.0831
4 0.104 0.0725 0.0600 0.0522
5 0.0714 0.0520 0.0419 0.0359
6 0.0516 0.0391 0.0309 0.0262
7 0.0390 0.0305 0.0238 0.0200
8 0.0304 0.0244 0.0188 0.0157
9 0.0244 0.0200 0.0153 0.0127
10 0.0200 0.0167 0.0126 0.0105
Orden SEM SIPG IIPG NIPG
1 0.709 0.408 0.408 0.353
2 0.577 0.365 0.333 0.305
3 0.593 0.393 0.338 0.301
4 0.604 0.420 0.347 0.302
5 0.608 0.443 0.357 0.306
6 0.608 0.461 0.364 0.309
7 0.608 0.476 0.371 0.312
8 0.607 0.488 0.376 0.314
9 0.607 0.498 0.381 0.316
10 0.607 0.506 0.384 0.318
Orden SEM SIPG IIPG NIPG
1 0.816 0.288 0.287 0.272
2 0.333 0.121 0.117 0.113
3 0.189 0.0683 0.0650 0.0622
4 0.120 0.0439 0.0413 0.0392
5 0.0823 0.0306 0.0286 0.0270
6 0.0595 0.0226 0.0210 0.0197
7 0.0449 0.0173 0.0161 0.0150
8 0.0350 0.0137 0.0127 0.0118
9 0.0281 0.0111 0.0103 0.00960
10 0.0230 0.00926 0.00852 0.00792
Orden SEM SIPG IIPG NIPG
1 0.816 0.288 0.287 0.272
2 0.666 0.243 0.235 0.226
3 0.684 0.247 0.235 0.225
4 0.697 0.254 0.239 0.227
5 0.700 0.261 0.244 0.230
6 0.700 0.266 0.247 0.232
7 0.700 0.271 0.251 0.235
8 0.699 0.274 0.253 0.236
9 0.698 0.278 0.256 0.238
10 0.698 0.280 0.258 0.240
Tabla 3. Constantes para las condiciones de
estabilidad en el caso elsco para SEM y las tres
versiones de DGM en funcin del tamao de los
elementos h ulizando a/b =1.4. Las constantes
para SIPG, IIPG y NIPG son para las bases po GLL
(De Basabe y Sen, 2010)
Tabla 4. Constantes para las condiciones deestabilidad en el caso elsco para SEM y las tres
versiones de DGM en funcin del tamao del
incremento mnimo Dx ulizando a/b =1.4. Las
constantes para SIPG, IIPG y NIPG son para las
bases po GLL. (De Basabe y Sen, 2010)
Tabla 1. Constantes para las condiciones de
estabilidad en el caso acsco para SEM y las tres
versiones de DGM en funcin del tamao de los
elementosh. Las constantes para SIPG, IIPG y NIPG
son para las bases po GLL (De Basabe y Sen, 2010)
Tabla 2. Constantes para las condiciones deestabilidad en el caso acsco para SEM y las
tres versiones de DGM en funcin del incremento
mnimo Dx. Las constantes para SIPG, IIPG y NIPG
son para las bases po GLL (De Basabe y Sen, 2010)
De Basabe: Elementos nitos de orden superior para propagacin de ondas ssmicas
237
-
7/21/2019 7_GEOS_31_2_Elementos
11/18
Por su parte, la dispersin numrica de
DGM en el caso acsco ha sido estudiada
por Hu et al. (1999), Stanescu et al.
(2000), Ainsworth (2004b) y Ainsworth
et al., (2006), y el caso elsco ha sido
estudiado por De Basabe et al. (2008).
Las principales conclusiones son que los
mtodos de orden 1 y 2 enen menor
dispersin numrica que los de SEM y que
los de orden superior enen las mismas
ventajas que los de SEM en cuanto
a la superconvergencia y anisotropa
numrica si se escogen cuidadosamente
los parmetros y pos de bases. La gura
3 muestra la convergencia del error de la
dispersin para la versin SIPG de DGM ylos nodos de GLL. Comparando esta gura
con la gura 1 se observa que el orden de
convergencia de este mtodo es igual al de
SEM. Las otras versiones de DGM enen
una convergencia ms lenta con respecto
a la dispersin numrica. La gura 4
muestran la anisotropa numrica de SIPG
ulizando los mismos parmetros que en
la gura 2 y los nodos de GLL. En esta gura
se observa que los mtodos de orden 3 y
4 no enen anisotropa numrica y que laanisotropa del mtodo de segundo orden
es menor que la de SEM de segundo orden
(comparar con la gura 2). La estabilidad
de DGM en los casos acsco y elsco
ha sido estudiada por De Basabe y Sen
(2010), y los resultados se muestran en las
tablas 1 al 4. Al comparar las condiciones
de estabilidad de las tres versiones de
DGM, se observa que la versin simtrica
(SIPG) ene un rango de estabilidad ms
amplio que NIPG y IIPG. Sin embargo, alcontrastar con las condiciones de SEM, se
observa que son casi 20% ms restricvas
en el caso acsco y ms de 60% ms
restricvas en el caso elsco.
Sismogramas Sintticos
En esta seccin, comparamos los
resultados numricos de SEM y DGM
para el problema de Lamb (Lamb, 1904),
ulizando un semiespacio isotrpico
homogneo con velocidades de ondas P y
S, y densidad dadas por 3 km/s, 1.73 km/s
y 2.5 g/cm (ver gura 5). La supercie
libre est ubicada en z=0. La fuente
ssmica ene un componente vercal
ubicada en el punto (0, 0.05 km), con una
ondcula de Ricker centrada en 0.08 s con
frecuencia central de 17.3 Hz. La solucin
exacta de este problema se puede obtener
ulizando el mtodo de Cagniard-De
Hoop (Aki y Richards, 2006, captulo 6).
Para calcular la solucin analca de este
problema, ulizamos el cdigo de Fortran
disponible en la pgina de internet del
proyecto SPICE.1
En las guras 6 y 7 se muestran los
sismogramas sintcos para una estacin
ubicada en la supercie, en el punto (1
km, 0), para SEM y DGM de orden 1, 2, 4
y 8. En las discrezaciones de diferentes
rdenes se ulizaron consistentemente10 nodos por longitud de onda para
la frecuencia central de la onda S. En
estas guras, se puede observar que
los mtodos de orden 1 introducen un
poco de dispersin numrica, la cual
se maniesta principalmente como
oscilaciones espurias despus del arribo
de la onda S. Al observar los errores en la
amplitud y la fase es posible apreciar como
estos errores disminuyen rpidamente
al incrementar el orden polinomial; enparcular, las aproximaciones de orden
8 son casi indisnguibles de la solucin
exacta. Para comparar los mtodos de
orden superior, en la gura 8 se muestran1 Seismic wave Propagaon and Imaging in Complex
media: an European network (SPICE)
hp://www.spice-rtn.org/library/soware/EXDDIR
GEOS, Vol. 31, No. 2
238
-
7/21/2019 7_GEOS_31_2_Elementos
12/18
los sismogramas obtenidos con SEM y
DGM de orden 4 y 8. En esta gura, se
muestra nicamente el arribo de la onda
S, que es donde se observan los errores
ms claramente. En esta escala es posible
apreciar que hay pequeos errores en
la amplitud y la fase, en parcular para
las aproximaciones de SEM. Aunque
aparentemente los resultados de DGM son
ms exactos, esto se debe en parte a que
se ulizaron incrementos en empo ms
pequeos debido a que las condiciones
de estabilidad son ms restricvas. Dicho
de otra manera, los errores observados se
atribuyen principalmente a la discrezacin
en el empo.
Superficie libre
Semiespacio elstico
Exacto
sem1
sem2
sem4
sem8
0.4 0.6 0.8
- 1.0
- 0.5
0.0
0.5
1.0
0.4
- 1.0
- 0.5
0.0
0.5
1.0
0.6 0.8
Exacto
dg1dg2
dg4
dg8
Exactosem4
sem8
dg4
dg8
- 1.0
- 0.5
0.0
0.5
1.0
0.65 0.70 0.75
Figura 5. Geometra del problema de Lamb. Este
problema consiste en un semiespacio elsco
con una supercie libre y una fuenta puntual.
Las lneas disconnuas representan los bordes
arciales introducidos al truncar el dominio para
las simulaciones.
Figura 6. Sismogramas sintcos para el problema
de Lamb ulizando SEM de orden 1, 2, 4 y 8. Eltamao de los elementos para la discrezacin en
el espacio fue seleccionado para obtener 10 nodos
para la longitud de la onda S de la frecuencia
central, y los incrementos en el empo ulizados
son los mximos permidos por las condiciones de
estabilidad.
Figura 7. Sismogramas sintcos para el problema
de Lamb ulizando DGM de orden 1, 2, 4 y 8. El
tamao de los elementos para la discrezacin en
el espacio fue seleccionado para obtener 10 nodos
para la longitud de la onda S de la frecuencia
central, y los incrementos en el empo ulizados
son los mximos permidos por las condiciones de
estabilidad.
Figura 8. Comparacin de los simogramas de
orden 4 y 8 ulizando SEM y DGM.
De Basabe: Elementos nitos de orden superior para propagacin de ondas ssmicas
239
-
7/21/2019 7_GEOS_31_2_Elementos
13/18
CONCLUSIONES
En este arculo hemos descrito
brevemente dos mtodos que se derivan
del de elementos nitos y que ofrecen
ventajas para simulaciones en sismologa.
Estos mtodos son el de elementos
espectrales y el de Galerkin disconnuo.
Para enfazar sus ventajas, hemos
presentado los principales resultados
del anlisis de dispersin y estabilidad
numrica, as como su aplicacin para
el problema de Lamb. A connuacin se
presentan algunas conclusiones:
Los mtodos de primer orden
requieren un nmero elevado de nodos porlongitud de onda debido a su dispersin y
anisotropa numrica, y son ms sensibles
a la razn entre las velocidades de las
ondas P y S.
La dispersin y anisotropa numrica
de los mtodos de orden 4 superior es
mnima, ulizando al menos 4 nodos por
longitud de onda.
Al comparar las tres versiones deIP-DGM, la simtrica (SIPG) ene mejores
propiedades en cuanto a dispersin y
estabilidad, en parcular si se ulizan los
nodos de Gauss o de GLL para construir las
funciones base.
Comparando SEM con IP-DGM, los
dos mtodos ofrecen ventajas similares en
cuanto a dispersin y anisotropa numrica,
pero IP-DGM ene la desventaja de que
sus condiciones de estabilidad son ms
restricvas. Sin embargo, IP-DGM ofrece
ventajas sobre SEM que no se discueron
en este arculo, como la posibilidad de
ulizar otros pos de elementos y bases, y
de ulizar mallas ms generales (e.g. Sun
y Wheeler, 2005; Kser y Dumsler, 2006).
AGRADECIMIENTOS
Agradezco al CONACYT por la beca
otorgada durante mi estancia en la
Universidad de Texas, Ausn, as como por
el apoyo por repatriacin. As mismo, a
Antonio Gonzlez por su revisin experta
que mejor este manuscrito, y a un revisor
annimo por sus observaciones.
REFERENCIAS
Ainsworth, M., 2004a. Discrete dispersion
relaon for hp-version nite element
approximaon at high wave number,SIAM Journal on Numerical Analysis,
42(2), 553-575.
Ainsworth, M., 2004b. Dispersive and
dissipave behaviour of high
order disconnuous Galerkin
nite element methods, Journal of
Computaonal Physics, 198(1), 106-
130.
Ainsworth, M., Monk, P., and Muniz, W.,
2006. Dispersive and dissipaveproperes of disconnuous Galerkin
nite element methods for the
second-order wave equaon,
Journal of Scienc Compung,
27(1), 5-40.
Aki, K. and Richards, P., 2002. Quantave
Seismology, University Science
Books, Sausalito, California, 2nd ed.
Alford, R., Kelly, K., and Boore, D., 1974.Accuracy of nite-dierence
modeling of the acousc wave
equaon, Geophysics, 39(6), 834-
842.
GEOS, Vol. 31, No. 2
240
-
7/21/2019 7_GEOS_31_2_Elementos
14/18
Alterman, Z. and Karal, F., 1968. Propagaon
of elasc waves in layered media by
nite dierence methods, Bullen
of the Seismological Society of
America, 58(1), 367-398.
Bao, H., Bielak, J., Ghaas, O., Kallivokas,
L., OHallaron, D., Shewchuk, J., and
Xu, J., 1998. Large-scale simulaon
of elasc wave propagaon in
heterogeneous media on parallel
computers, Computer Methods in
Applied Mechanics and Engineering,
152(1-2), 85-102.
Carcione, J., Herman, G., and Ten Kroode,
A., 2002. Seismic modeling,Geophysics, 67(4), 1304-1325.
Cerven, V., 2001. Seismic ray theory,
Cambridge Univ Press.
Chaljub, E., Komatsch, D., Viloe, J.-P.,
Capdeville, Y., Valee, B., and Festa,
G., 2007. Spectral element analysis
in seismology, in Advances in Wave
Propagaon in Heterogeneous
Media, Advances in Geophysics, pp.
365-419, eds Wu, R.-S. & Maupin, V.,
Elsevier.
Chopra, A. K., Dibaj, M., Clough, R. W.,
Penzien, J., and Seed, H. B., 1969.
Earthquake analysis of earth
dams, 4th World Conf. Earthquake
Engineering, pp. A5-55.
Cockburn, B. and Shu, C., 1989. TVB Runge-
Kua Local Projecon Disconnuous
Galerkin Finite Element Methodfor Conservaon Laws II: General
Framework, Mathemacs of
Computaon, 52(186), 411-435.
Cockburn, B. and Shu, C.-W., 1998. The
Runge-Kua disconnuous Galerkin
method for conservaon laws V:
muldimensional systems, Journal
of Computaonal Physics, 141(2),
199-224.
Cockburn, B., Lin, S., and Shu, C., 1989.
TVB Runge-Kua local projecon
disconnuous Galerkin nite
element method for conservaon
laws III: one-dimensional systems,
Journal of Computaonal Physics,
84(1), 90-113.
Cockburn, B., Hou, S., and Shu, C.-W., 1990.
The Runge-Kua Local ProjeconDisconnuous Galerkin Finite
Element Method for Conservaon
Laws. IV: The Muldimensional
Case, Mathemacs of Computaon,
54(190), 545-581.
Cohen, G., 2002. Higher-Order Numerical
Methods for Transient Wave
Equaons, Scienc Computaon,
Springer-Verlag.
Courant, R., Friedrichs, K., and Lewy, H.,
1967. On the paral dierence
equaons of mathemacal physics,
IBM Journal of Research and
Development, 11(2), 215-234.
Darlow, B., 1980. A Penalty-Galerkin
Method for Solving the Miscible
Displacement Problem, Ph.D. thesis,
Rice University, Houston, Texas.
Dawson, C., Sun, S., and Wheeler, M.F., 2004. Compable algorithms
for coupled ow and transport,
Computer Methods in Applied
Mechanics and Engineering, 193(23-
26), 2565-2580.
De Basabe: Elementos nitos de orden superior para propagacin de ondas ssmicas
241
-
7/21/2019 7_GEOS_31_2_Elementos
15/18
De Basabe, J. D. and Sen, M. K., 2007.
Grid dispersion and stability criteria
of some common nite-element
methods for acousc and elasc
wave equaons, Geophysics, 72(6),
T81-T95.
De Basabe, J. D. and Sen, M. K., 2010.
Stability of the high-order nite
elements for acousc or elasc
wave propagaon with high-order
me stepping, Geophysical Journal
Internaonal, 181(1), 577-590.
De Basabe, J. D., Sen, M. K., and Wheeler,
M. F., 2008. The interior penalty
disconnuous Galerkin method forelasc wave propagaon: Grid dis-
persion, Geophysical Journal Inter-
naonal, 175(1), 83-93.
de la Puente, J., Kser, M., Dumbser, M.,
and Igel, H., 2007. An arbitrary high-
order disconnuous Galerkin meth-
od for elasc waves on unstructured
meshes -- IV. Anisotropy, Geophysi-
cal Journal Internaonal, 169(3),
1210-1228.
Dormy, E. and Tarantola, A., 1995. Nu-
merical simulaon of elasc wave
propagaon using a nite volume
method, Journal of Geophysical
Research, 100(B2), 2123-2134.
Dumbser, M. and Kser, M., 2006. An
arbitrary high-order disconnuous
Galerkin method for elasc waves
on unstructured meshes -- II. Thethree-dimensional isotropic case,
Geophysical Journal Internaonal,
167(1), 319-336.
Dumbser, M. and Munz, C.-D., 2006.
Building Blocks for Arbitrary High
Order Disconnuous Galerkin
Schemes, Journal of Scienc Com-
pung, 27(1-3), 215-230.
Dumbser, M., Kser, M., and de la Pu-
ente, J., 2007a. Arbitrary high-order
nite volume schemes for seismic
wave propagaon on unstructured
meshes in 2D and 3D, Geophysical
Journal Internaonal, 171(2), 665-
694.
Dumbser, M., Kser, M., and Toro, E.,
2007b. An arbitrary high-order
Disconnuous Galerkin method
for elasc waves on unstructured
meshes -- V. Local me stepping and
p-adapvity, Geophysical Journal
Internaonal, 171(2), 695-717.
Fornberg, B., 1987. The pseudospectral
method: Comparisons with nite
dierences for the elasc wave
equaon, Geophysics, 52(4), 483-
501.
Grote, M., Schneebeli, A., and Schotzau,
D., 2006. Disconnuous Galerkin
Finite Element Method for the
Wave Equaon, SIAM Journal on
Numerical Analysis, 44(6), 2408-
2431.
Hu, F., Hussaini, M., and Rasitarinera,
P., 1999. An analysis of the
disconnuous Galerkin method
for wave propagaon problems,Journal of Computaonal Physics,
151(2), 921-946.
GEOS, Vol. 31, No. 2
242
-
7/21/2019 7_GEOS_31_2_Elementos
16/18
Kser, M. and Dumbser, M., 2006. An
arbitrary high-order disconnuous
Galerkin method for elasc waves
on unstructured meshes -- I. The
two-dimensional isotropic case with
external source terms, Geophysical
Journal Internaonal, 166(2), 855-
877.
Kser, M., Dumbser, M., de la Puente,
J., and Igel, H., 2007a. An arbitrary
high-order disconnuous Galerkin
method for elasc waves on
unstructured meshes - III. Viscoelasc
aenuaon, Geophysical Journal
Internaonal, 168(1), 224-242.
Kser, M., Mai, P., and Dumbser, M.,
2007b. Accurate Calculaon of
Fault-Rupture Models Using the
High-Order Disconnuous Galerkin
Method on Tetrahedral Meshes,
Bullen of the Seismological Society
of America, 97(5), 1570-1586.
Kelly, K., Ward, R., Treitel, S., and Alford,
R., 1976. Synthec seismograms
-A nite-dierence approach,Geophysics, 41(1), 2-27.
Komatsch, D. and Tromp, J., 1999.
Introducon to the spectral element
method for three-dimensional
seismic wave propagaon,
Geophysical Journal Internaonal,
139(3), 806-822.
Komatsch, D. and Tromp, J., 2002a.
Spectral-element simulaons of
global seismic wave propagaon-I.
Validaon, Geophysical Journal
Internaonal, 149(2), 390-412.
Komatsch, D. and Tromp, J., 2002b.
Spectral-element simulaons of
global seismic wave propagaon-
II. Three-dimensional models,
oceans, rotaon and self-gravitaon,
Geophysical Journal Internaonal,
150(1), 303-318.
Komatsch, D. and Viloe, J., 1998.
The spectral-element method:
an ecient tool to simulate the
seismic response of 2D and 3D
geological structures, Bullen of the
Seismological Society of America,
88(2), 368-392.
Komatsch, D., Ritsema, J., and Tromp, J.,2002. The spectral-element method,
Beowulf compung, and global
seismology, Science, 298 (5599),
1737-1742.
Komatsch, D., Liu, Q., Tromp, J., Sss,
P., Sdham, C., and Shaw, J., 2004.
Simulaons of Ground Moon in the
Los Angeles Basin Based upon the
Spectral-Element Method, Bullen
of the Seismological Society ofAmerica, 94(1), 187-206.
Komatsch, D., Tsuboi, S., and Tromp,
J., 2005. The spectral-element
method in seismology, Geophysical
Monograph, 157, 205-227.
Koslo, D. and Baysal, E., 1982. Forward
modeling by a Fourier method,
Geophysics, 47(10), 1402-1412.
Lamb, H., 1904. On the Propagaon ofTremors over the Surface of an Elasc
Solid, Philosophical Transacons of
the Royal Society of London, 203,
1-42.
De Basabe: Elementos nitos de orden superior para propagacin de ondas ssmicas
243
-
7/21/2019 7_GEOS_31_2_Elementos
17/18
Lax, P. D. and Wendro, B., 1964.
Dierence schemes for hyperbolicequaons with high order ofaccuracy, Communicaons on Pure
and Applied Mathemacs, 17(3),381-398.
Levander, A., 1988. Fourth-order nite-dierence P-SV seismograms,
Geophysics, 53, 1425.
Lysmer, J. and Drake, L., 1972. A niteelement method for seismology, in
Methods in Computaonal Physics,vol. 11, pp. 181-216, Academic Press.
Madariaga, R., 1976. Dynamics of an
expanding circular fault, Bullen ofthe Seismological Society of America,
66(3), 639-666.
Madec, R., Komatsch, D., and Diaz,
J., 2009. Energy-conserving localme stepping based on high-order
nite elements for seismic wavepropagaon across a uid-solidinterface, Computer Modeling in
Engineering and Sciences, 49(2),
163-189.
Madrid, J., 2008. A geometrical approachto me evolving wave fronts,
Geophysical Journal Internaonal,172(3), 1117-1122.
Marfurt, K. J., 1984. Accuracy of nite-
dierence and nite-elementmodeling of the scalar and elasc
wave equaons, Geophysics, 49(5),533--549.
Marfurt, K. J., 1990. Analysys of higher
order nite element methods,
in Numerical modeling of seismic
wave propagaon, vol. 13, pp. 516-
520, ed. Marfurt, K. J., Soc. of Expl.
Geophys.
Moczo, P., Kristek, J., and Halada, L., 2000.
3D fourth-order staggered-grid
nite-dierence schemes: stability
and grid dispersion, Bullen of the
Seismological Society of America,
90(3), 587-603.
Mulder, W. A., 1999. Spurious modes
in nite-element discrezaons
of the wave equaon may not be
all that bad, Applied Numerical
Mathemacs, 30(4), 425-445.
Mullen, R. and Belytschko, T., 1982.
Dispersion analysis of nite
element semidiscrezaons of the
two-dimensional wave equaon,Internaonal Journal for Numerical
Methods in Engineering, 18(1), 11-
29.
Patera, A. T., 1984. A spectral element
method for uid dynamics: Laminar
ow in a channel expansion, Journal
of Computaonal Physics, 54(3),
468-488.
Priolo, E., 2001. Earthquake Ground
Moon Simulaon Through the 2-D
Spectral Element Method, Journal
of Computaonal Acouscs, 9(4),
1561-1581.
Reed, W. and Hill, T., 1973. Triangular mesh
methods for the neutron transport
equaon, Los Alamos Scienc
Laboratory Technical Report, LA-
UR-73-479.
Rivire, B. and Wheeler, M. F., 2001.Disconnuous nite element
methods for acousc and elasc
wave problems. part I: Semidiscrete
error esmates, TICAM report, 0(01-
02), 1-12.
GEOS, Vol. 31, No. 2
244
-
7/21/2019 7_GEOS_31_2_Elementos
18/18
Rivire, B., Wheeler, M., and Girault, V.,
1999. Improved energy esmates
for interior penalty, constrained
and disconnuous Galerkin
methods for ellipc problems. Part
I, Computaonal Geosciences, 3(3),
337-360.
Rivire, B., Wheeler, M. F., and Girault,
V., 2001. A Priori Error Esmates
for Finite Element Methods Based
on Disconnuous Approximaon
Spaces for Ellipc Problems, SIAM
Journal on Numerical Analysis,
39(3), 902-931.
Seriani, G. and Priolo, E., 1994. Spectralelement method for acousc wave
simulaon in heterogeneous media,
Finite Elements in Analysis and
Design, 16(3-4), 337-348.
Smith, W. D., 1975. The applicaon of
nite element analysis of body wave
propagaon problems, Geophysical
Journal of the Royal Astronomical
Society, 42(2), 747 -768.
Stanescu, D., Kopriva, D., and Hussaini,
M., 2000. Dispersion Analysis for
Disconnuous Spectral Element
Methods, Journal of Scienc
Compung, 15(2), 149-171.
Sun, S. and Wheeler, M., 2005.
Disconnuous Galerkin methods for
coupled ow and reacve transport
problems, Applied Numerical
Mathemacs, 52(2-3), 273-298.
Tordjman, N., 1995. Elments nis dordre
lev avec condensaon de masse
pour lquaon des ondes, Ph.D.
thesis, Universit de Paris, IX-
Dauphine.
Tordjman, N., Cohen, G., and Joly, P.,
1994. Elments nis dordre lev
avec condensaon de masse pour
lquaon des ondes en dimension
1, Rapport de Recherche 2323,INRIA, Le Chesnay, France.
Virieux, J., 1984. SH wave propagaon
in heterogeneous media: velocity-
stress nite-dierence method,
Geophysics, 49(11), 1933-1957.
Virieux, J., 1986. P-SV wave propagaon
in heterogeneous media: velocity-
stress nite-dierence method,
Geophysics, 51(4), 889-901.
Manuscrito recibido: 31 de mayo, 2011.
Recepcin del manuscrito corregido: 29 de
sepembre, 2011.
Manuscrito aceptado: 10 de diciembre, 2011.
De Basabe: Elementos nitos de orden superior para propagacin de ondas ssmicas
245