Study of matter behavior dynamics governed by the interaction with laser pulses and external strong currents

Abstract

This thesis studies the physical phenomena that occur during the deposition of energy in matter and consists of two parts. In the first part, the nanosecond pulsed laser irradiation of thin metal films-substrate systems is investigated for thermoelastic, melting and ablation regimes, using the finite element method (FEM). The numerical simulations are compared and validated with experimental results, obtained by dynamic imaging interferometry and white-light tomographic interferometry. In the second part, the sequential stages of explosion of a Z-pinch copper wire from solid to plasma formation and its plasma expansion are investigated using multiphysics coupled numerical simulations, validated by experimental results obtained by interferometry, shadowgraphy, schlieren and diffraction imaging techniques. In more details, for the first part coupled thermal-structural, transient models based on FEM are developed to give a comprehensive spatiotemporal numerical solution of the physical phenomena occurring in laser matter interactions. Temperature dependent material properties of metal coatings deposited on glass substrates are used for the simulations. Initially, a 2D axisymmetric model is developed to study the generation and propagation of laser generated ultrasounds, for laser fluences below the melting threshold of metallic film-substrates. Afterwards, a 3D quarter-symmetric FEM model is developed and used for all regimes of interest. The developed FEM modeling provides a simultaneous analysis of the thermal and structural parameters, as defined by the solution of the heat conduction and wave propagation equations. The wave equation determines the displacements of the target imposed by the laser energy deposition, while the heat conduction equation predicts the temperature distribution induced by absorption of the laser pulse in the target. The 3D model computes the phase changes of matter by taking into account the latent heats of melting and vaporization, depending on the laser fluence. Ablation is simulated by the ‘killing’ of the elements that exceed the boiling temperature and are subsequently deactivated. The attenuation of the laser irradiance at the target surface due to plasma absorption is also taken into account. With regard to loading conditions, a heat source term is used, representing the laser energy absorbed by the sample and described by a Gaussian distribution in time and space. The Lagrangian mesh is locally adaptive, depending on the simulation needs. Moreover, the elastic and plastic behaviour of thin films is being investigated by taking into account stress-strains temperature-dependent curves of different materials until fracture into the 3D model. Furthermore, the finite element model was further extended to simulate the presence of surface and solid volume defects (gaps) and their subsequent influence on the generation and propagation of ultrasonic waves (surface acoustic waves). For the second part that concerns the matter behaviour dynamics, governed by the interaction with external strong currents, a 3D coupled mechanical/thermal FEM model, simulating the initial stages of explosion of a Z-pinch metallic copper wire, was initially developed. The model conveyed a simultaneous analysis of the thermal and structural parameters, as defined by the solution of the heat conduction and mechanical motion equations. The mechanical equation determines the displacements of the wire imposed by an alternative nanosecond pulsed current while the heat equation predicts the temperature distribution. The source term of the heat generation rate is the Joule heating term and is used as loading condition, while the pulsed imposed current is provided from experiments.In order to investigate the crucial physical phenomena that take place in the initial stages of wire explosion, a 3D electromagnetic-thermal-structural hydrodynamic FEM model is developed. Large volumetric deformation is considered by taking into account the material’s hydrodynamic behavior via equations of state (Gruneisen and multiphase tabular) and the elastoplastic behavior is also considered by a flow stress constitutive model (Johnson-Cook). A Lagrangian transient analysis has been carried out. The loading condition of the simulation is a nanosecond pulsed alternative current, which is provided from experiments. Τhe sequential stages of explosion of a Z-pinch copper wire from solid to plasma formation and its plasma expansion are investigated. The expansion of the exploded material is being investigated with a magnetohydrodynamic code that uses as initial condition the density distribution and radius instabilities from the above aforementioned 3D FEM model.

Η παρούσα διατριβή μελετάει τα φυσικά φαινόμενα που λαμβάνουν χώρα κατά την εναπόθεση ενέργειας στην ύλη και αποτελείται από δύο μέρη στα οποία παρουσιάζονται δύο διαφορετικοί μηχανισμοί εναπόθεσης. Αναφορικά με τη δυναμική συμπεριφορά της ύλης όταν αλληλεπιδρά με παλμούς λέιζερ, η ακτινοβόληση λεπτών μεταλλικών φιλμ-υποστρωμάτων από παλμικό λέιζερ διάρκειας νανοδευτερολέπτων ερευνάται για την θερμοελαστική περιοχή, την περιοχή τήξης και την περιοχή φωτοαποδόμησης χρησιμοποιώντας τη μέθοδο των πεπερασμένων στοιχείων (ΜΠΣ). Οι αριθμητικές προσομοιώσεις συγκρίνονται και επικυρώνονται με τα πειραματικά αποτελέσματα που προκύπτουν από την τεχνική δυναμικής συμβολομετρίας απεικόνισης. Αναφορικά με τη δυναμική συμπεριφορά της ύλης όταν αλληλεπιδρά με εξωτερικά ισχυρά ρεύματα, τα διαδοχικά στάδια της έκρηξης ενός χάλκινου σύρματος Ζ-pinch από στερεά κατάσταση ως τη δημιουργία πλάσματος καθώς και η εξάπλωση του πλάσματος διερευνώνται αναπτύσσοντας multiphysics αριθμητικές προσομοιώσεις συζευγμένων πεδιών οι οποίες επικυρώνονται από πειραματικά αποτελέσματα που λαμβάνονται από οπτικές τεχνικές απεικόνισης (συμβολομετρία, σκιαγράφηση, Schlieren, απεικόνιση περίθλασης Fraunhofer).Αναλυτικότερα για το πρώτο μέρος, αναπτύσσεται συζευγμένο θερμικό-μηχανικό μοντέλο δυναμικά μεταβαλλόμενο στο χρόνο που βασίζεται στη ΜΠΣ με σκοπό να δώσει μια ολοκληρωμένη χωροχρονική αριθμητική επίλυση των φυσικών φαινομένων που συμβαίνουν κατά τη διάρκεια της αλληλεπίδραση ύλης με λέιζερ. Για τις προσομοιώσεις λαμβάνονται υπόψη θερμοκρασιακά εξαρτώμενες ιδιότητες υλικών, των μεταλλικών επικαλύψεων που εναποτίθενται σε γυάλινα υποστρώματα. Αρχικά, αναπτύχθηκε δισδιάστατο 2D αξονοσυμμετρικό μοντέλο με σκοπό να μελετηθεί η δημιουργία και διάδοση υπερήχων που παράγονται για πυκνότητες ροής ακτινοβολίας λέιζερ κάτω από το όριο τήξης του μεταλλικού φιλμ-υποστρώματος. Στη συνέχεια τρισδιάστατο 3D τέταρτο-συμμετρικό μοντέλο πεπερασμένων στοιχείων αναπτύχθηκε για όλες τις περιοχές ενδιαφέροντος (θερμοελαστική, τήξης, φωτοαποδόμησης/πλάσματος). Η αναπτυχθείσα μεθοδολογία παρέχει μια ταυτόχρονη ανάλυση των θερμικών και μηχανικών παραμέτρων όπως ορίζεται από την επίλυση των εξισώσεων της θερμικής αγωγιμότητας και κυματικής διάδοσης. Η κυματική εξίσωση καθορίζει τις μετατοπίσεις του στόχου που επιβάλλονται από την εναπόθεση ενέργειας λέιζερ, ενώ η εξίσωση αγωγής θερμότητας προβλέπει την κατανομή της θερμοκρασίας που προκαλείται από την απορρόφηση της ενέργειας του παλμικού λέιζερ στον στόχο. Το 3D μοντέλο υπολογίζει τις αλλαγές φάσης της ύλης, λαμβάνοντας υπόψη τις λανθάνουσες θερμότητες τήξης και εξάτμισης, ανάλογα με την πυκνότητα ροής της ακτινοβολίας λέιζερ. Η φωτοαποδόμηση προσομοιώνεται από την τεχνική «killing» των στοιχείων που υπερβαίνουν τη θερμοκρασία βρασμού και στη συνέχεια απενεργοποιούνται. Η εξασθένηση της ακτινοβολίας λέιζερ στην επιφάνεια-στόχο λόγω της απορρόφησης του πλάσματος λαμβάνεται επίσης υπόψη. Όσον αφορά τις συνθήκες φόρτισης, χρησιμοποιείται ως πηγή θερμότητας, η ενέργεια του λέιζερ που απορροφάται από το δείγμα και περιγράφεται από μια Γκαουσιανή κατανομή σε χρόνο και χώρο. Το Λαγκρανζιανό πλέγμα είναι τοπικά προσαρμοζόμενο ανάλογα με τις ανάγκες της προσομοίωσης. Επιπλέον, η ελαστική και πλαστική συμπεριφορά των λεπτών υμενίων ερευνάται λαμβάνοντας υπόψη, θερμοκρασιακά εξαρτώμενες καμπύλες τάσεων-παραμορφώσεων από διαφορετικά υλικά, μέχρι τη θραύση τους, στο 3D μοντέλο. Επιπλέον, το μοντέλο πεπερασμένων στοιχείων επεκτείνεται περαιτέρω, προκειμένου να προσομοιωθεί η παρουσία των επιφανειακών και των ελαττωμάτων του όγκου στο στερεό (κενά) και να μελετηθεί η μετέπειτα επιρροή τους στην παραγωγή και διάδοση των υπερηχητικών κυμάτων (επιφανειακών ακουστικών κυμάτων). Αναφορικά με τη δυναμική συμπεριφορά της ύλης όταν αλληλεπιδρά με εξωτερικά ισχυρά ρεύματα, αναπτύχθηκε αρχικά ένα 3-D συζευγμένο μηχανικό-θερμικό μοντέλο πεπερασμένων στοιχείων, που προσομοιώνει τα αρχικά στάδια της έκρηξης ενός Ζ-pinch μεταλλικού χάλκινου σύρματος. Το μοντέλο λαμβάνει υπόψη τις θερμοκρασιακά εξαρτώμενες θερμικές και μηχανικές παραμέτρους μέσω της επίλυσης των εξισώσεων της θερμικής αγωγιμότητας και της κλασσικής μηχανικής. Η μηχανική εξίσωση προσδιορίζει τις μετατοπίσεις του σύρματος που επιβάλλεται από εναλλασόμενο παλμικό ρεύμα διάρκειας νανοδευτερολέπτων, ενώ η εξίσωση αγωγής θερμότητας προβλέπει την κατανομή της θερμοκρασίας. Ο όρος που λειτουργεί ως πηγή του ρυθμού παραγωγής θερμότητας είναι ο όρος της θέρμανσης Joule και χρησιμοποιείται ως συνθήκη φόρτισης, ενώ το παλμικό εναλλασσόμενο ρεύμα παρέχεται από πειράματα. Για να αντιμετωπιστούν τα φυσικά φαινόμενα που λαμβάνουν χώρα στα αρχικά στάδια της έκρηξης του σύρματος πιο ρεαλιστικά αναπτύσσεται ένα 3D ηλεκτρομαγνητικό-θερμικό-μηχανικό υδροδυναμικό μοντέλο βασιζόμενο στη ΜΠΣ. Μεγάλες ογκομετρικές παραμορφώσεις εξετάζονται λαμβάνοντας υπόψη την υδροδυναμική συμπεριφορά του υλικού μέσω καταστατικών εξισώσεων (Gruneisen και καταστατική σε μορφή πίνακα λαμβάνοντας υπόψη τις αλλαγές φάσης), ενώ η ελαστοπλαστική συμπεριφορά του υλικού μελετάται θεωρώντας επίσης καταστατικό μοντέλο τάσεων (Johnson-Cook). Πραγματοποιείται Λαγκρανζιανή δυναμική ανάλυση στο χρόνο, ενώ η πηγή φόρτισης της προσομοίωσης είναι παλμικό εναλλασσόμενο ρεύμα, διάρκειας νανοδευτερολέπτων το οποίο παρέχεται από πειράματα. Μελετώνται τα διαδοχικά στάδια της έκρηξης χάλκινου σύρματος Ζ-pinch από στερεά κατάσταση ως τη δημιουργία πλάσματος. Στη συνέχεια η διαστολή και διάδοση του πλάσματος στο χωροχρόνο διερευνάται με μαγνητοϋδροδυναμικό κώδικα ο οποίος χρησιμοποιεί ως αρχικές συνθήκες την κατανομή πυκνότητας και τις αστάθειες στην τιμή της ακτίνας από το προαναφερθέν 3D μοντέλο με τη ΜΠΣ.