Κατά την εκβολή ιξωδοελαστικών ρευστών μπορεί να δημιουργηθούν αστάθειες οι οποίες είναι πολύ σημαντικές σε διάφορες βιομηχανικές εφαρμογές κατά τις οποίες λαμβάνει χώρα μορφοποίηση πολυμερών. Για την προσομοίωση αυτής της διεργασίας χρησιμοποιήσαμε τη μέθοδο των μικτών πεπερασμένων στοιχείων μαζί με μια ελλειπτική μέθοδο δημιουργίας του πλέγματος. Για την μοντελοποίηση της ιξωδοελαστικής συμπεριφοράς του ρευστού χρησιμοποιήσαμε το μοντέλο Phan-Thien-Tanner και για την επίλυση του χρησιμοποιήθηκε η τεχνική EVSS-G σε συνδυασμό με την μέθοδο SUPG για την διακριτοποίηση της καταστατικής εξίσωσης. Πραγματοποιήσαμε γραμμική ανάλυση ευστάθειας στην διδιάστατη αυτή ροή σύμφωνα με την οποία το σύστημα γίνεται ασταθές πάνω από μια κρίσιμη τιμή του αριθμού Wi όπου τα πεδία ταχύτητας, πίεσης και τάσεων αρχίζουν να αυξάνουν εκθετικά με το χρόνο με ή χωρίς σύγχρονη περιοδική μεταβολή. Η αστάθεια γεννάται σε μια περιοχή γύρω από το σημείο εκβολής που είναι και ιδιάζον σημείο για τη ροή και η διαταραχή εκτείνεται 2-5 ακτίνες μακριά από αυτό, ανάλογα το λόγο των τριχοειδών προς τις ιξώδεις δυνάμεις, αριθμός Ca, ενώ αφήνει το υπόλοιπο ρευστό σχεδόν αδιατάρακτο. Πραγματοποιούμε παραμετρική ανάλυση μεταβάλλοντας τον αριθμό Ca καθώς και τις σχετικές ρεολογικές παραμέτρους του καταστατικού μοντέλου των τάσεων προκειμένου να μελετήσουμε την επίδραση των παραμέτρων αυτών στην ευστάθεια της ροής κατά την εκβολή.