Simulaci n num rica de reflectores ac sticos
This presentation is the property of its rightful owner.
Sponsored Links
1 / 61

Simulación numérica de reflectores acústicos PowerPoint PPT Presentation


  • 54 Views
  • Uploaded on
  • Presentation posted in: General

II Jornadas Interdisciplinarias de Acústica. Simulación numérica de reflectores acústicos. Federico Miyara Laboratorio de Acústica y Electroacústica UNR [email protected] Laboratorio de Acústica y Electroacústica. Introducción.

Download Presentation

Simulación numérica de reflectores acústicos

An Image/Link below is provided (as is) to download presentation

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.


- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -

Presentation Transcript


II Jornadas Interdisciplinarias de Acústica

Simulación numérica de reflectores acústicos

Federico MiyaraLaboratorio de Acústica y Electroacústica [email protected]

Laboratorio

de Acústica y

Electroacústica


Introducción


Al igual que la luz, las ondas sonoras son capaces de reflejarse en ciertas superficies.

La condición para ello es que la superficie sea rígida y de gran densidad superficial.


Si además la superficie es muy grande con respecto a la longitud de onda, la reflexión es especular.


En el caso óptico, las superficies reflectantes típicas de la vida diaria son, efectivamente, mucho mayores que la longitud de onda y las reflexiones son especulares.

Es el caso de los espejos y otras superficies pulidas.


En el caso acústico, la longitud de onda correspondiente a sonidos audibles está entre 1,7 cm y 17 m.

Por ejemplo, la longitud de onda del tono fundamental de la voz hablada varía entre 1,5 m y 3,5 m.


Para reflectores del tamaño de una pared típica, la reflexión no es puramente especular.

Esto significa que una onda sonora unidireccional no se refleja en una sola dirección.


Los reflectores pequeños con respecto a la longitud de onda son frecuentes en la vida diaria.

Por ejemplo muebles, ventanas, paredes.


En acústica arquitectónica tienen aplicación directa con la finalidad de direccionar las reflexiones del sonido.


El formalismo que permite evaluar el comportamiento real de estos reflectores es la teoría de la difracción.

Es matemáticamente compleja y sólo se conocen soluciones aproximadas para algunos casos.


Una posible alternativa para obtener soluciones prácticas al problema es la simulación numérica.


Para ello se dividirá el panel reflector en pequeños elementos y se considerará cada uno de ellos como un emisor secundario de acuerdo con el principio de Huygens.


R

F


Veremos primero el modelo de radiación acústica de una fuente simple.


Consideremos un pequeño radiador acústico F con simetría esférica

p(r, t)

r

F


Se demuestra que la solución de la ecuación de onda es una cuasi onda


p

p(r, t1)

t1

1/r

r

p

1/r

t2 > t1

p(r, t2)

r

p

1/r

t3 > t2

p(r, t3)

r

p

t4 > t3

1/r

p(r, t4)

r


Un posible modelo físico de un radiador con simetría esférica es la esfera pulsante,cuya superficie se mueve a una velocidad u(t) senoidal


Por una cuestión de simplicidad de cálculos trabajaremos en realidad con

u(t) = U ejt


Ue jt

a


Se demuestra que la presión sonora a una distancia r y en un instante t es

Dondeo: densidad del airec : velocidad del sonido


Se puede reescribir como

R(r)

T(t)


O bien:

R(r)

T(t)


Si a/c << 1, lo cual se cumple sia <<  , se puede aproximar

Pero 4a2U es el caudal de aire G movido por la fuente, llamado poder de la fuente


Es decir,

Una fuente de este tipo se conoce como fuente simple.


p(r)

t1

t2

t3

r


La notación exponencial permite separar el comportamiento temporal del espacial, mediante el concepto de fasor


NOTA: El poder G de una fuente simple puede obtenerse a partir de la presión Po medida a una distancia rodespejándola de la anterior:


Si sólo interesa la magnitud,


Se puede demostrar que una fuente pequeña, no necesariamente esférica, que pulsa en forma tal que toda su superficie oscile en fase moviendo un caudal G,se comporta, para grandes distancias, como fuente simple


Esto permite calcular el campo lejano de una fuente de cualquier forma que oscila en fase en toda su superficie como el de una fuente simple de igual caudal total


Vamos a aplicar esta propiedad a cada elemento en que subdividimos el reflector y luego superpondremos los campos lejanos de todos ellos


y

x

b

R

r2

yR

a

xR

y

x

zF

xF

zR

yF

z

r1

a

b

F


Para poder considerar al elemento dxdy como una fuente simple deberemos asignarle un caudal elemental dG


Podemos obtener el caudal desplazado por F en la posición del elemento dxdy en ausencia del reflectorcalculando la velocidad y multiplicándola por el área proyectada


Se demuestra, por ser una onda de simetría esférica, que


La presencia del reflector equivale a que hubiera una fuente imagen del otro lado


R

x

a

r2

x

z

r1

r1

a

F

F'


La velocidad que genera F' hacia el lado derecho debe multiplicarse por el área proyectada sobre la perpendicular del rayo


R

x

a

r2

x

z

r1

r1

a

F

F'


Pero para que esa superficie pueda asimilarse a una fuente simple, debe suponerse que hay otra superficie que se mueve en forma opuesta.


En efecto, toda la superficie de la fuente debe moverse en fase para que su campo lejano sea asimilable al de una fuente simple.

(Los bordes se desprecian por ser su área nula).


Por lo tanto, el poder de la fuente elemental, es decir el caudal que la abandona, es el doble del caudal que atraviesa el plano del reflector


Entonces


A partir del caudal dG de la fuente elemental se puede obtener su efecto en R


Es decir:


Finalmente, la presión debida a la reflexión se obtiene integrando:


Esta integral no tiene solución explícita excepto si se lleva a coordenadas polares y se integra sobre un plano infinito, en cuyo caso se obtiene la reflexión especular


Si a, b < , queda el recurso de integrar numéricamente.

Para ello transformamos la integral en una sumatoria, donde reemplazamos:

dx =a/N

dy =b/M


Para que los elementos en que se divide el reflector se puedan considerar fuentes simples deben ser mucho menores que la longitud de onda   =  2c/:

a/N  <<  

b/M  <<  


Esto implica que


Un problema es que la cantidad de cálculos requeridos se vuelve muy grande, haciendo que el cálculo sea lento, ya que se requiere evaluar el integrando 4NM veces.


Por ejemplo, para f = 1000 Hz y un panel de 2 m3 m se requiere

N >> 8,7M >> 5,8

Esto lleva a alrededor de 20000 cálculos


AGRADECIMIENTOSMaría Victoria GómezFrancisco y Andrés Miyara Cristina BiassoniAdAACINTRA


  • Login