Visualizza wikitesto di Raggi acustici in mare: curvatura strato isotermo


[[File:265dtc.jpg|thumb|right |300px|Caso tipico di doppia curvatura dei raggi acustici: si generano zone d'ombra non illuminate dal suono ( sezione grigia )]]

Dall'andamento delle traiettorie dei raggi acustici in mare dipendono le capacità di scoperta dei sonar.

Le diverse situazioni che si possono verificare a seguito dei molteplici andamenti delle  traiettorie dei raggi acustici si estendono dall'impossibilità di qualsiasi rilevamento sonar, nelle zone d'ombra, alle capacità di scoperta dei bersagli a distanze eccezionali nei casi di canalizzazione del suono.
  
Per queste ragioni lo studio dell'andamento delle traiettorie dei raggi acustici in mare risulta di primaria importanza nella conduzione del sonar.

Gli esercizi svolti  a titolo d'esempio, nella voce, sono stati estrapolati  dal testo "FONDAMENTALS OF SONAR" di J.W.Horton traducendo le formule di calcolo in routine software. 


== La curvatura dei raggi acustici in mare==


I raggi acustici dovuti all'energia  emessa in mare da una sorgente di suono sono soggetti a curvatura dovuta a due cause diverse:

*la prima, più significativa, dipende dalla variazione della velocità del suono in funzione della temperatura dell'acqua con la profondità e da altre cause minori.

*la seconda,  meno rilevante, dipende dalla variazione della velocità del suono dovuta alla densità del mezzo in funzione della profondità; in questo caso si considera costante la temperatura nell'intervallo di quota interessato ai computi.

Entrambe le situazioni sono governate dalla legge di Snell:

<math> \left ( \frac{\cos \delta _{{1}}}{c_{{1}}} \right) = \left( \frac{\cos \delta _{{2}}}{c_{{2}}} \right)</math>

[[File:242dtc.jpg|thumb|left |400px|figura 01]]

La legge di Snell è illustrata in modo semplice nel grafico di figura 01 dove, invece di una curva, è mostrata una spezzata:

Dalla figura si vede che nel tratto di mare compreso tra quota <math> Z = 0 </math> e quota <math> Z = 100 </math> m la velocità del suono è di <math> {c_{{1}}} = 1480  </math>  m/S. ed  il raggio acustico indicato in celeste si sviluppa secondo l'inclinazione <math> {\delta_{{1}}} </math>, una volta raggiunta la quota <math> Z = 100 </math> m il raggio s'immette nel tratto di mare compreso tra <math> Z = 100 </math>  m e <math> Z = 200 </math>  m dove la velocità del suono e di <math> {c_{{2}}} = 1440  </math>  m/S.  in questo tratto il raggio assume l'inclinazione <math> {\delta_{{2}}} </math> come mostra il tratto viola.

In questo esempio è ipotizzata, ad arte, una variazione della velocità del suono  a scalino, in realtà la variazione della velocità del suono in mare è continua e genera andamenti dei raggi acustici dal percorso curvilineo.

Dalla legge di Snell derivano gli algoritmi per il tracciamento dei raggi acustici in mare illustrati nella presente voce.

==Le grandezze delle variabili utilizzate nelle fasi di calcolo a seguire==

Le grandezze utilizzate per le computazioni sono espresse con unità di misura anglosassoni:

Temperature: (°F ) in gradi Fahrenheit  ( °F = °C x 9/5 + 32 )

Profondità: (ft) in feet   ( ft = mt x 3.281 ) 

Distanze: (yd) in yard      ( yd = m x 1.094 )

Gli angoli, in grandezza universale, sono espressi in radianti

==L'obiettivo dei calcoli==
[[File:241dtc.jpg|thumb|right| 400px|]]

I calcoli hanno come obiettivo la valutazione di n coppie di coordinate in grado di consentire il
tracciamento di una curva, la traiettoria del percorso del suono, che si sviluppa secondo una o più archi di circonferenza così come indicato in figura 1. 

In questo particolare procedimento l'ascissa <math>(S)</math> , distanza dalla sorgente acustica espressa in yard, è la variabile dipendente mentre l'ordinata <math>(D)</math>, variabile indipendente, è la profondità dei punti di calcolo in feet.

L'arco denominato " angolo di radenza" è l'angolo formato tra la tangente del cerchio e un segmento parallelo all'asse delle ascisse passante per il punto di tangenza.

La valutazione del valore di quest'angolo è fondamentale nel definire l'andamento della <math>S=f(D).</math>
 
Per ciascuna coppia di coordinate il valore della variabile <math>D </math> è stabilito dall'operatore
mentre il valore di <math>S</math> ne è una funzione : <math>S = f(D).</math>  

Più elevato è il numero delle coppie più punti definiscono la traiettoria.

I calcoli devono essere eseguiti a 9 decimali.

==  Variazione della temperatura dell’acqua in funzione della profondità==


[[File:2433dtc.jpg|thumb|left| 400px|figura 2]]


Prima d’iniziare il percorso di calcolo accennato è utile mostrare il grafico di figura 2, detto batitermogramma, che mostra, in questo caso specifico, come varia la temperatura dell’acqua in mare,  misurata con [[Batitermografo a perdere|batitermografo]], in funzione della profondità e di conseguenza come varia con la temperatura, la velocità del suono

Lo strato isotermo si estende, con un tratto verticale,  dalla superficie del mare a <math> 65 </math> feet  a seguire, con un tratto inclinato,  la temperatura decresce con l'aumentare della profondità.

Dal grafico assumeremo, nello strato isotermo, la temperatura  T = 70 °F  da utilizzare pe il calcolo della velocità del suono che n'é una variabile.

==Procedura di calcolo nello strato a temperatura costante==

Il calcolo di <math> S = f(D) </math> mira alla determinazione delle coordinate dei punti della traiettoria <math>D \ ; \ S(D)</math>, ed è svolto secondo i seguenti passi:

* p0- si assumono un numero di valori (<math>D</math>) (ordinate dei punti della traiettoria) necessari alla definizione del tracciato

* p1- si assume il gradiente di velocità del suono nello strato isotermo pari a: <math>(g1) = 0.0182 </math> (la velocità del suono. nello strato isotermo, varia soltanto in funzione della pressione con la quota)

* p2- si stabilisce di collocare il trasduttore alla prof. di <math> 15 ft </math> (valore utilizzato per i trasduttori nelle navi di superficie)

* p3- si fissa un valore della temperatura <math>(T) </math> dell'acqua  
 
* p4- si fissa la salinità <math>(sa)</math>, in parti per mille, in base ad ipotesi realistiche

* p5- si stabilisce la profondità <math>(d1)</math>  dello strato nel quale eseguire i calcoli (quota del trasduttore)

* p6- si calcola il valore della velocità del suono <math>C </math> alla quota del trasduttore

* p7- si calcola il raggio <math>(r1)</math> dell'arco di cerchio che forma la traiettoria 

* p8- si calcola <math> \cos \theta</math> ;  coseno dell'angolo di radenza

* p9- si ricava l'angolo <math> \theta</math>

* p10- si calcola infine <math> S = f(D) </math> con l'espressione:  <math> S = [r1 \cdot \sen (z)] / 3 </math>

== Le funzioni di calcolo==

Algoritmi per il calcolo dell’ascissa <math> S(D) </math>:

L'ascissa  si computa con <math> S = [r1 \cdot \sen (z)] / 3 </math>

Dove:

<math>z = \arccos \theta</math>  

<math>\theta </math> = angolo di radenza

<math> \cos \theta = (1 + Dd) / r1 </math> (coseno dell’angolo di radenza)

<math> r1 = (C / g1) </math> (raggio della traiettoria)



<math>C = 4411 + (11.25) \cdot (T) - 0.045 \cdot (T ^ 2) + g1 \cdot d1 + 4.3 \cdot (sa - 34) </math> velocità del suono in corrispondenza del  trasduttore, dove:

<math>T </math> = temperatura nello strato isotermo

<math> g1  </math> = gradiente di velocità nello strato isotermo

<math> d1 </math> = quota del trasduttore

<math>sa </math> = salinità nello strato

==Applicazione delle funzioni di calcolo==
In questo esercizio consideriamo un '''raggio acustico orizzontale''' emesso dal trasduttore avente, ovviamente, l'angolo di radenza <math>\theta = 0 </math>, assumeremo inoltre per la temperatura costante nello strato il valore <math> T = 70 </math>°F indicato in figura 3.

Si voglia calcolare la coppia di coordinate <math>D; S(D)</math>
di un punto della traiettoria per
<math>  D = 10 </math > ft (variabile indipendente)  

<math>d1 = 15 </math> ampiezza dello strato   isotermo (da 0 a 15 ft), dalla superficie al trasduttore d’emissione

<math>g1 = 0.0182  </math> gradiente assunto per la strato isotermo  

<math> Dd = D - d1 </math>  = <math> - 5 </math>

<math>Sa = 35</math> parti per mille

<math>T = 70 ^ \circ F </math>

<math> C = 4411 + (11.25) \cdot (70) - 0.045 \cdot (70 ^ 2) + 0.0182 \cdot 15 + 4.3 \cdot (35 - 34) </math> = <math>4982 </math> ft/s.

<math> r1  = 4982/0.0182 </math>  = <math> 273736 </math> yd


<math> \cos \theta = 1 - 4 / 273736 </math> = <math> 0.99998173</math>


<math> z  = \arccos \ (0.99998173) </math> = <math> 0.006044842 </math><ref> Si osservi che l'angolo di radenza <math>\theta = 0 </math>, all'altezza del trasduttore d'emissione, si modifica con il decrescere del valore <math> D  </math>  </ref>
 

<math> S = (273736  \cdot \sin \ (0.00604442) / 3 </math> = <math> 551 </math> yd


Queste le coordinate del punto: 

<math>  D = 10 </math> ft ; <math> S = 551 </math> yd

==Il tracciamento della traiettoria==
[[File:24555dtc.jpg|thumb|right |400px|figura 3]]

Per una presentazione grafica della traiettoria del raggio acustico si dovrebbero ripetere i calcoli precedentemente fatti per un numero molto elevato di coordinate <math>D ; S</math> pur non ottenendo una soddisfacente presentazione dell'arco di curva.

Il problema si risolve con  l'impiego di un eseguibile di calcolo sviluppato in Visual Basic; si tratta di una routine impostata  secondo gli algoritmi illustrati in precedenza in grado di eseguire le computazione a passi di piccole frazioni della variabile indipendente <math>D</math>.

Il risultato dell’utilizzo dell’eseguibile è mostrato in figura 3 dove a sinistra del tracciato del raggio acustico è disposto il grafico dell'andamento della temperatura già illustrato in figura 2

La variabile indipendente <math>D</math> decresce da <math> D = 15 </math> a <math> D = 0 </math> feet,  la funzione <math>S(D)</math> cresce da <math>S = 0 </math> a <math> S = 880 </math> yard

Sulla curva è evidenziato con un cerchietto rosso il punto calcolato nell'esercizio precedente; <math>D = 10 </math> feet e <math>S = 551 </math> yard.

Nel valutare visivamente il tracciato si deve osservare che l'accentuata curvatura del raggio  acustico dipende dalla differenza di scala tra le coordinate:  <math> D </math> è in feet mentre <math> S </math> è in yard.

== L’eseguibile di calcolo==

[[File:244dtc.jpg|thumb|left |400px|figura 4]]

Come accennato in precedenza  sull'utilizzo di un eseguibile di calcolo
se ne propone l’impiego scaricandolo all'indirizzo riportato nella sezione Collegamenti esterni


Il programma di calcolo consente, oltre il tracciamento della curva di figura 3, anche la possibilità di eseguire calcoli per punti illustrando in tal modo tutto il processo manuale di elaborazione numerica svolto in precedenza presentando un tabulato con funzioni e dati così come riportato in figura 4.

==note==
<references/>

==Bibliografia==

G.Pazienza, ''Fondamenti della localizzazione marina'', La Spezia, Studio grafico Restani, 1970.
	
W. Horton, ''Foundamentals of Sonar'', United States Naval Institute,Annapolis Maryland, 1959

A. De Dominics Rotondi, ''Principi di elettroacustica subacquea'' , Elettronica San Giorgio-Elsag S.p.A., Genova, 1990.

C. Del Turco, '' Sonar- Principi - Tecnologie – Applicazioni '', Tip. Moderna La Spezia, 1992

==Collegamenti esterni==

[http://www.sonar-info.info/wikiraggiisotermo/wikiraggiisotermo.html  Wikiraggiisotermo]

o su 

[https://github.com/2021dtc/prova10/archive/main.zip  raggiisotermo]

[[Categoria: sonar]]