Il modello di interazione tra un getto supersonico protostellare ed il mezzo ambiente è stato sviluppato usando il codice numerico FLASH (Fryxell et al. 2000; vedi Sez. 12.1), sviluppato dall'Alliances Center for Astrophysical Thermonuclear Flashes (ASCI) presso l'università di Chicago. Le simulazioni numeriche sono state condotte presso il Sistema di Calcolo per l'Astrofisica Numerica (SCAN) dell'Osservatorio Astronomico di Palermo e presso il CINECA (Bologna), sulla base di proposte di assegnazione di tempo di calcolo accettate.
L'evoluzione di un getto supersonico che viaggia nel mezzo ambiente inizialmente imperturbato è descritta dalle equazioni idrodinamiche che esprimono la conservazione della massa, del momento e dell'energia. Il nostro modello tiene conto delle perdite radiative da plasma otticamente sottile e della conduzione termica.
Per le simulazioni numeriche, abbiamo adottato un sistema di coordinate bidimensionale cilindrico con l'asse del getto coincidente con l'asse .
Assumendo equilibrio di pressione iniziale tra getto e mezzo ambiente, i parametri liberi che determinano l'evoluzione del sistema si riducono a due parametri adimensionali: il numero di Mach, definito come il rapporto tra la velocità del flusso e la velocità del suono nel mezzo ambiente, , ed il contrasto di densità tra mezzo ambiente e getto, .
Per determinare quali condizioni fisiche di un getto protostellare possano riprodurre i dati sperimentali, relativi per esempio al moto proprio e all'emissione X osservati, abbiamo condotto una dettagliata e completa esplorazione dello spazio dei parametri che principalmente influenzano l'evoluzione del sistema getto - ambiente mostrata in Fig. 25.
L'eplorazione della spazio dei parametri permette di vincolare esattamente l'intervallo delle configurazioni fisiche che un getto protostellare può avere per poter riprodurre l'emissione X ed il moto proprio del fronte di shock osservati.
Il modello di getto protostellare che meglio riproduce le osservazioni è quello di un getto con numero di Mach iniziale , con densità iniziale minore rispetto a quella dell'ambiente in cui si propaga, cioè con rapporto di densità ambiente - getto e con densità iniziale cm e temperatura K.
I risultati ottenuti riguardo a questo modello sono stati presentati in Bonito et al. (2004) e (2007).
La Fig. 26 mostra come evolve nel tempo la temperatura, la densità e l'emissione X del getto protostellare.
Alla testa del getto è evidente la presenza di una regione calda e densa immediatamente alle spalle del fronte di shock che corrisponde alla regione dalla quale viene originata l'emissione X prevista dal modello come osservata da Chandra/ACIS-I.
Confrontando la posizione della bolla calda e densa che corrisponde alla sorgente X al variare del tempo, si può notare la presenza di un moto proprio misurabile della regione emettente. In particolare abbiamo ricavato un valore medio della velocità associata al fronte di shock pari a km/s.
Il valore di velocità ricavato è in ottimo accordo con le osservazioni ( km/s), tenedo conto del fatto che HH è inclinato di gradi rispetto alla linea di vista.
Dalle simulazioni numeriche bidimensionali, abbiamo derivato la misura di emissione e la temperatura, , per ogni elemento di fluido; da queste distribuzioni 2-D, considerando la simmetria assiale, abbiamo ricavato le distribuzioni spaziali 3-D della e della . Integrando lungo la linea di vista e tenendo condo dell'assorbimento interstellare, determinato dalle osservazioni, e della risposta strumentale, abbiamo ottenuto la mappa dell'emissione X al piano focale dello strumento scelto da confrontare direttamente con le osservazioni.
Dalla , usando il codice spettrale MEKAL e tenendo in considerazione la risposta spettrale degli strumenti a bordo di XMM-Newton e Chandra, abbiamo sintetizzato gli spettri di piano focale come verrebbero rivelati da tali satelliti.
Dall'analisi spettrale, si può concludere che l'emissione X simulata risulta consistente con l'emissione da un plasma otticamente sottile ad una sola temperatura, sia nel caso di bassa statistica (cioè considerando al più fotoni raggruppati in modo da avere almeno conteggi per canale) che nel caso di alta statistica (con fotoni complessivi raggruppati in modo da averne per canale, caso meno realistico, ma che consente di avere più dettagli).
A partire dalla componente isoterma derivata dall'analisi spettrale, abbiamo ricavato la luminosità X della regione emettente. Abbiamo ottenuto (Bonito et al. 2004 e 2007) che la luminosità X, , varia tra erg/s, valori del tutto consistenti con quello osservato in HH : erg/s.