Benchmark test cases for the National Combustion Code

Author(s):  
K.-H. Chen ◽  
A. Norris ◽  
A. Quealy ◽  
N.-S. Liu
Keyword(s):  
Author(s):  
Klaus-Peter Helbig ◽  
Pedro J. Coelho ◽  
Joaquim Chamico ◽  
John Marakis ◽  
Aristide Mbiock ◽  
...  

1998 ◽  
Vol 06 (01n02) ◽  
pp. 29-43 ◽  
Author(s):  
Mark R. Fallat ◽  
Stan E. Dosso

This paper applies simulated annealing matched-field inversion to the Workshop '97 geoacoustic inversion benchmark test cases. The unknown parameters in these test cases include various combinations of sediment and basement geoacoustic properties (layer thickness, compressional and shear speeds, attenuations, and densities) and geometric parameters (water depth, and source range and depth). The sensitivity of the acoustic field data to the different parameters varies substantially. The simulated annealing inversion developed here minimizes a measure of mismatch based on the average Bartlett processor for acoustic data at a number of ranges. The replica pressure fields are computed by summing the modal components provided by an elastic normal mode model, which allows replica fields at multiple ranges to be calculated efficiently. A total of 15 different geoacoustic inversions are considered (three realizations each for five test cases). Excellent results are obtained in a reasonable amount of computation time for all cases. In particular, the final mismatch values are exceedingly small and relatively sensitive parameters (geometric parameters and sediment properties) are generally estimated to high precision. Relatively insensitive parameters (shear speeds and, in some cases, basement properties) are somewhat less well determined, but generally good estimates are obtained.


2015 ◽  
Author(s):  
Κωνσταντίνος Στόκος

Σκοπός της διδακτορικής διατριβής ήταν η ανάπτυξη μιας αριθμητικής-υπολογιστικής μεθοδολογίας για την προσομοίωση γενικά μη μόνιμων ροών με ταυτόχρονη μεταφορά θερμότητας.Ο αναπτυχθείς επιλύτης είναι κατάλληλος για την προσομοίωση δισδιάστατων ή τρισδιάστατων, ασυμπίεστων, στρωτών ή τυρβωδών ροών συνεκτικού και Νευτώνειου ρευστού. Για την απαραίτητη σύζευξη των πεδίων πίεσης και ταχύτητας εφαρμόζεται η προσέγγιση της τεχνητής συμπιεστότητας. Για την πρόβλεψη ανωστικών ροών του ασυμπίεστου ρευστού, η μεταβολή της πυκνότητας λόγω θερμοκρασιακών διαφορών προσομοιώνεται σύμφωνα με την προσέγγιση Boussinesq. Για τις τυρβώδεις ροές χρησιμοποιούνται δύο ευρέως διαδεδομένα RANS μοντέλα τύρβης, το μοντέλο τύρβης k-ω SST στη χαμηλών και υψηλών αριθμών Re έκδοσή του και το υψηλών αριθμών Re k-ε μοντέλο τύρβης.Ο επιλύτης εφαρμόζει μια κεντρο-κομβική (node-centered) πεπερασμένων όγκων τεχνική διακριτοποίησης, χρησιμοποιώντας ακμο-βασική (edge-based) και διαφανή (transparent) προσέγγιση σε υβριδικά αριθμητικά πλέγματα. Για τον υπολογισμό των μη-συνεκτικών όρων αναπτύχθηκε εκ του μηδενός ο προσεγγιστικός Riemann επιλύτης του Roe. Οι συνεκτικοί όροι διακριτοποιούνται χρησιμοποιώντας ένα κεντρικό σχήμα. Η χρονική διακριτοποίηση επιτυγχάνεται μέσω ενός πλήρως πεπλεγμένου σχήματος για τη χρονοπροέλαση στον ψευδο-χρόνο και το φυσικό χρόνο.Όλες οι εξισώσεις μέσης ροής (εξισώσεις συνέχειας, ορμής και ενέργειας) επιλύονται ισχυρά συζευγμένες. Η ισχυρή σύζευξη των εξισώσεων προτιμήθηκε μετά από σύγκριση με την ασθενώς συζευγμένη επίλυση, σύγκριση την οποία δεν έχουμε βρει στη βιβλιογραφία να παρουσιάζεται για το συγκεκριμένο σύστημα εξισώσεων. Η ισχυρή σύζευξη έδειξε γρηγορότερη σύγκλιση για σημαντικά μη-συνεκτικά φαινόμενα και μας επέτρεψε τη χρήση υψηλότερων αριθμών CFL κάνοντας τη σύγκλιση ακόμα γρηγορότερη.Για τη μοντελοποίηση της θερμικής ακτινοβολίας εφαρμόστηκαν δύο μεθοδολογίες. Η πρώτη αποτελεί μια αναλυτική μεθοδολογία που βασίζεται στη χρήση γωνιακών παραγόντων και η δεύτερη αποτελεί μια πεπερασμένων όγκων μεθοδολογία.Μετά από την ανάπτυξη της αριθμητικής μεθοδολογίας σε υπολογιστικό κώδικα επιλύθηκαν αρκετές περιπτώσεις αναφοράς (benchmark test cases) για την επαλήθευση της ορθής λειτουργίας του κώδικα και την αξιολόγηση των δυνατοτήτων του. Τα αποτελέσματα είναι αρκετά ενθαρρυντικά όπως συμβαίνει και με την παραλληλοποίηση του επιλύτη, που παρουσίασε σημαντική μείωση του υπολογιστικού χρόνου.Τέλος, προχωρήσαμε στην προσομοίωση δύο σεναρίων πυρκαγιάς σε αεριζόμενη σήραγγα. Οι περιπτώσεις πυρκαγιάς προσομοιώθηκαν είτε λαμβάνοντας υπόψη την ακτινοβολία και την αγωγή στο τοίχωμα είτε όχι. Η θερμική ακτινοβολία και η αγωγή στο τοίχωμα έδειξαν ότι είναι απαραίτητες για την πρόβλεψη ρεαλιστικών αποτελεσμάτων.


2018 ◽  
Vol 115 (35) ◽  
pp. E8125-E8134 ◽  
Author(s):  
Xiangyang Guan ◽  
Cynthia Chen

A generic modeling framework to infer the failure-spreading process based on failure times of individual nodes is proposed and tested in four simulation studies: one for cascading failures in interdependent power and transportation networks, one for influenza epidemics, one benchmark test case for congestion cascade in a transportation network, and one benchmark test case for cascading power outages. Four general failure-spreading mechanisms—external, temporal, spatial, and functional—are quantified to capture what drives the spreading of failures. With the failure time of each node given, the proposed methodology demonstrates remarkable capability of inferring the underlying general failure-spreading mechanisms and accurately reconstructing the failure-spreading process in all four simulation studies. The analysis of the two benchmark test cases also reveals the robustness of the proposed methodology: It is shown that a failure-spreading process embedded by specific failure-spreading mechanisms such as flow redistribution can be captured with low uncertainty by our model. The proposed methodology thereby presents a promising channel for providing a generally applicable framework for modeling, understanding, and controlling failure spreading in a variety of systems.


1994 ◽  
Vol 144 ◽  
pp. 503-505
Author(s):  
R. Erdélyi ◽  
M. Goossens ◽  
S. Poedts

AbstractThe stationary state of resonant absorption of linear, MHD waves in cylindrical magnetic flux tubes is studied in viscous, compressible MHD with a numerical code using finite element discretization. The full viscosity tensor with the five viscosity coefficients as given by Braginskii is included in the analysis. Our computations reproduce the absorption rates obtained by Lou in scalar viscous MHD and Goossens and Poedts in resistive MHD, which guarantee the numerical accuracy of the tensorial viscous MHD code.


Sign in / Sign up

Export Citation Format

Share Document