Propagación numérica de ondas de choque acústicas Informe semestral Roberto Velasco Segura 1 Director: Dr. Pablo Luis Rendón Garrido 1 1 Centro de Ciencias Aplicadas y Desarrollo Tecnológico Universidad Nacional Autónoma de México Posgrado en Ciencias Físicas, 2012
28
Embed
Propagación numérica de ondas de choque acústicas ...lya.fciencias.unam.mx/computocientifico/archivos/presentacion.pdf · Ley de Fourier Tiempos de relajación cortos X0 X0 es
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Propagación numérica de ondas de choqueacústicas
Informe semestral
Roberto Velasco Segura1
Director: Dr. Pablo Luis Rendón Garrido1
1Centro de Ciencias Aplicadas y Desarrollo TecnológicoUniversidad Nacional Autónoma de México
Posgrado en Ciencias Físicas, 2012
Marco teórico Esquema numérico Resultados preliminares Administración de información Perspectivas
Objetivo
Modelar numéricamente propagación de ondas de choqueen la atmósfera terrestre.Particularmente: cáusticas de ondas de choque generadaspor un avión.
Propagación numérica de ondas de choque acústicas 2/27
Marco teórico Esquema numérico Resultados preliminares Administración de información Perspectivas
Contenido
1 Marco teórico
2 Esquema numérico
3 Resultados preliminares
4 Administración de información
5 Perspectivas
Propagación numérica de ondas de choque acústicas 3/27
Marco teóricoHamilton, Blackstock (ed.), Nonlinear acoustics, 1998,Academic Press.
Marco teórico Esquema numérico Resultados preliminares Administración de información Perspectivas
¿Qué ecuación vamos a modelar numéricamente?gas heterogéneono linealdisipativados o tres dimensiones(conveniente para tratamiento numérico)
Propagación numérica de ondas de choque acústicas 5/27
Marco teórico Esquema numérico Resultados preliminares Administración de información Perspectivas
Conservación de masa
∂ρ
∂t+∇ ·
(ρ~u)
= 0
(no tiene aproximaciones)
Propagación numérica de ondas de choque acústicas 6/27
Marco teórico Esquema numérico Resultados preliminares Administración de información Perspectivas
Conservación de momento
∂ρ~u∂t
+∇ · (ρ~u ⊗ ~u) = ~B +∇ · σ
Hipótesis de fluido newtoniano: la fuerza de fricción esproporcional a la velocidad relativa.La dependencia de µ y µB, con respecto a lascoordenadas espaciales es baja.~B = 0.
∂ρ~u∂t
+∇ · (ρ~u ⊗ ~u) = −∇p + µ∇2~u +
(µB +
13µ
)∇(∇ · ~u
)
Propagación numérica de ondas de choque acústicas 7/27
Marco teórico Esquema numérico Resultados preliminares Administración de información Perspectivas
Conservación de energía
∂
∂t
(12ρu2 + ρE
)+∇ · ~u
(12ρu2 + ρE
)= ∇ ·
(σ · ~u − ~q
)Heredamos la hipótesis de fluido newtoniano.El cambio de κ con respecto a las coordenadas espacialeses bajo.Ley de Fourier para conducción térmica.
ρTDsDt
= κ∇2T + µB(∇ · ~u)2 +12µ
(∂ui
∂xj+∂uj
∂xi+
23δij∂uk
∂xk
)2
Propagación numérica de ondas de choque acústicas 8/27
Marco teórico Esquema numérico Resultados preliminares Administración de información Perspectivas
Ecuaciones de estado
p = p(ρ, s)
T = T (ρ, s)
Los tiempos de relajación son más cortos que los cambiosinvolucrados en nuestras soluciones.
Propagación numérica de ondas de choque acústicas 9/27
Marco teórico Esquema numérico Resultados preliminares Administración de información Perspectivas
Paso de mecánica de fluidos a acústica
Se asume un valor de equilibrio y una perturbación, para cadauna de las variables.
ρ = ρ0 + ρ′
p = p0 + p′
T = T0 + T ′
s = s0 + s′
∃ ε� 1 tal que ε ∼ ρ′
ρ0∼ p′
p0∼ T ′
T0
Las variables ~u y s se manejan diferente
ε ∼ uc0
Propagación numérica de ondas de choque acústicas 10/27
Marco teórico Esquema numérico Resultados preliminares Administración de información Perspectivas
Primer orden
Ecuación de onda lineal
�2p′ = 0
Hipótesis principales:Se cumplen leyes de conservación y ecuación de estado aorden ε.s′ = 0.
(involucra una sola variable)
Propagación numérica de ondas de choque acústicas 11/27
Marco teórico Esquema numérico Resultados preliminares Administración de información Perspectivas
Segundo orden
Ecuación de Westervelt
�2p′ +δ
c40
∂3p′
∂t3 =−βρ0c4
0
∂2(p′)2
∂t2
Hipótesis:Fluido newtonianoNo hay fuerzasexternasLey de FourierTiempos derelajación cortosX ′
X0es del orden de ε
Se cumplen leyes deconservación y ecuación deestado a orden ε2 y s′
µ, µB y κ son constantesNos encontramos lejos de lasparedes, y esto implica ∇× ~u = 0La no linealidad acumuladadomina a la local
Propagación numérica de ondas de choque acústicas 12/27
Marco teórico Esquema numérico Resultados preliminares Administración de información Perspectivas
Adimensionalización
�2p = −Θ∂3p∂ t3− ∂2p2
∂ t2
Nuevas variables
p =βp′
c20ρ0
x =xλ
y =yλ
z =zλ
t = tc0
λ
Parámetro adimensional
Θ =δ
c0λ
Propagación numérica de ondas de choque acústicas 13/27
Marco teórico Esquema numérico Resultados preliminares Administración de información Perspectivas
Comparación
Una dimensiónlineal
ut + c0ux = 0
no lineal (Burgers)
ut + uux = νuxx
Dos o más dimensioneslineal
�2p = 0
no lineal (Westervelt)
�2p = −Θ∂3p∂t3 −
∂2p2
∂t2
Propagación numérica de ondas de choque acústicas 14/27
Esquema numérico
Marco teórico Esquema numérico Resultados preliminares Administración de información Perspectivas
Aproximaciones de diferencias finitas
∂2p∂x2 '
Pni−1,j − 2Pn
i,j + Pni+1,j
∆x2 + O((∆x)2)
∂2p∂y2 '
Pni,j−1 − 2Pn
i,j + Pni,j+1
∆y2 + O((∆y)2)
∂2p∂t2 '
Pn−1i,j − 2Pn
i,j + Pn+1i,j
∆t2 + O((∆t)2)
∂3p∂t3 '
−Pn−2i,j + 3Pn−1
i,j − 3Pni,j + Pn+1
i,j
∆t3 + O(∆t)
Propagación numérica de ondas de choque acústicas 16/27
Marco teórico Esquema numérico Resultados preliminares Administración de información Perspectivas
Esquema numérico
Pn+1i,j =
−K1 ±√
K 21 − 4K0
2
Explícito2DSegundo orden espacialPrimer orden temporalMalla cartesiana ∆x = ∆yImplementado en paralelo sobre un GPU
Karamalis et al., Fast ultrasound image simulation using thewestervelt equation, 2010, Springer.
Propagación numérica de ondas de choque acústicas 17/27
Resultados preliminares
Marco teórico Esquema numérico Resultados preliminares Administración de información Perspectivas
Observaciones
Para amplitudes bajas 10−6 los resultados numéricosreproducen cualitativamente el caso lineal.Para amplitudes grandes 10−1 se observa formación deondas de choque en distancias muy cortas.Se observa interacción de ondas de choque provenientesde direcciones no paralelas.
Propagación numérica de ondas de choque acústicas 19/27
Administración de información
Marco teórico Esquema numérico Resultados preliminares Administración de información Perspectivas
http://dondecae.net/t2
ReferenciasFichasConceptosIdeas sueltasResultados
Propagación numérica de ondas de choque acústicas 21/27
Perspectivas
Marco teórico Esquema numérico Resultados preliminares Administración de información Perspectivas
Marco teórico
Introducir heterogeneidades buscando modelarcondiciones atmosféricasAlternativas a la ecuación de Westervelt
Partiendo directamente de las ecuaciones de conservación