A continuous adjoint-based approach for the optimization of wing flapping

Author(s):  
Mingjun Wei ◽  
Min Xu
Keyword(s):  
2017 ◽  
Vol 89 (3) ◽  
pp. 444-456
Author(s):  
Lei Chen ◽  
Jiang Chen

Purpose This paper aims to conduct the optimization of the multi-stage gas turbine with the effect of the cooling air injection based on the adjoint method. Design/methodology/approach Continuous adjoint method is combined with the S2 surface code. Findings The optimization of the stagger angles, stacking lines and the passage can improve the attack angles and restrain the development of the boundary, reducing the secondary flow loss caused by the cooling air injection. Practical implications The aerodynamic performance of the gas turbine can be improved via the optimization of blade and passage based on the adjoint method. Originality/value The results of the first study on the adjoint method applied to the S2 surface through flow calculation including the cooling air effect are presented.


2017 ◽  
Vol 10 (3) ◽  
pp. 1157-1174 ◽  
Author(s):  
Yosuke Niwa ◽  
Hirofumi Tomita ◽  
Masaki Satoh ◽  
Ryoichi Imasu ◽  
Yousuke Sawa ◽  
...  

Abstract. A four-dimensional variational (4D-Var) method is a popular algorithm for inverting atmospheric greenhouse gas (GHG) measurements. In order to meet the computationally intense 4D-Var iterative calculation, offline forward and adjoint transport models are developed based on the Nonhydrostatic ICosahedral Atmospheric Model (NICAM). By introducing flexibility into the temporal resolution of the input meteorological data, the forward model developed in this study is not only computationally efficient, it is also found to nearly match the transport performance of the online model. In a transport simulation of atmospheric carbon dioxide (CO2), the data-thinning error (error resulting from reduction in the time resolution of the meteorological data used to drive the offline transport model) is minimized by employing high temporal resolution data of the vertical diffusion coefficient; with a low 6-hourly temporal resolution, significant concentration biases near the surface are introduced. The new adjoint model can be run in discrete or continuous adjoint mode for the advection process. The discrete adjoint is characterized by perfect adjoint relationship with the forward model that switches off the flux limiter, while the continuous adjoint is characterized by an imperfect but reasonable adjoint relationship with its corresponding forward model. In the latter case, both the forward and adjoint models use the flux limiter to ensure the monotonicity of tracer concentrations and sensitivities. Trajectory analysis for high CO2 concentration events are performed to test adjoint sensitivities. We also demonstrate the potential usefulness of our adjoint model for diagnosing tracer transport. Both the offline forward and adjoint models have computational efficiency about 10 times higher than the online model. A description of our new 4D-Var system that includes an optimization method, along with its application in an atmospheric CO2 inversion and the effects of using either the discrete or continuous adjoint method, is presented in an accompanying paper Niwa et al.(2016).


2021 ◽  
Vol 9 ◽  
Author(s):  
Andrea Maggi ◽  
Dominik Garmatter ◽  
Sebastian Sager ◽  
Martin Stoll ◽  
Kai Sundmacher

A chemical plant layout for the production of syngas from renewable power, H2O and biogas, is presented to ensure a steady productivity of syngas with a constant H2-to-CO ratio under time-dependent electricity provision. An electrolyzer supplies H2 to the reverse water-gas shift reactor. The system compensates for a drop in electricity supply by gradually operating a tri-reforming reactor, fed with pure O2 directly from the electrolyzer or from an intermediate generic buffering device. After the introduction of modeling assumptions and governing equations, suitable reactor parameters are identified. Finally, two optimal control problems are investigated, where computationally expensive model evaluations are lifted viaparareal and necessary objective derivatives are calculated via the continuous adjoint method. For the first time, modeling, simulation, and optimal control are applied to a combination of the reverse water-gas shift and tri-reforming reactor, exploring a promising pathway in the conversion of renewable power into chemicals.


2016 ◽  
Author(s):  
Ιωάννης Καββαδίας

Η παρούσα διδακτορική διατριβή ασχολείται με τη μαθηματική ανάπτυξη, τον προγραμματισμό και την πιστοποίηση των συνεχών συζυγών (continuous adjoint) μεθόδων για χρονικά μόνιμες και μή-μόνιμες τυρβώδεις ροές με έμφαση στην ακρίβεια των υπολογιζομένων παραγώγων ευαισθησίας για συναρτήσεις-στόχους που συναντώνται στην αεροδυναμική. Εξετάζονται προβλήματα βελτιστοποίησης μορφής καθώς και ελέγχου της ροής με χρήση δεσμών ρευστού, τόσο σε ακαδημαϊκά προβλήματα όσο και σε εφαρμογές της βιομηχανίας. Σχετικά με τη διαφόριση των μοντέλων τύρβης, η συνεχής συζυγής μέθοδος επεκτείνεται για να καλύπτει ροές που μοντελοποιούνται με το μοντέλο k−ω SST για πρώτη ϕορά στη σχετική βιβλιογραφία. Η αναλυτική διαφόριση του μοντέλου k−ω SST παρουσιάζει σημαντικές δυσκολίες καθώς περιέχει μη-διαφορίσιμες συναρτήσεις και προτείνεται η κατάλληλη αντιμετώπισή τους. Το συζυγές πρόβλημα διατυπώνεται τόσο για τη χαμηλών (LowRe) όσο και την υψηλών αριθμών (HighRe) (συναρτήσεις τοίχου) Reynolds της τύρβης εκδοχή του μοντέλου k−ω SST . Ταυτοχρόνως, εξετάζονται οι επιπτώσεις του να αμελείται η παραγώγιση του μοντέλου τύρβης κατά την ανάπτυξη της συζυγούς μεθόδου (frozen turbulence assumption). Η συζυγής μέθοδος αναπτύσσεται δύο φορές, στη βάση δύο διαφορετικών μαθηματικών προσεγγίσεων. Και με τις δύο παράγονται οι ίδιες συζυγείς εξισώσεις και οριακές συνθήκες. Διαφορές εντοπίζονται στον τύπο υπολογισμού των παραγώγων ευαισθησίας. Σύμφωνα με την πρώτη, τη μέθοδο των Επιφανειακών Ολοκληρωμάτων (Surface Integral -SI- adjoint formulation), αυτές εκφράζονται αποκλειστικά με επιφανειακά ολοκληρώματα. Αντίθετα, στη δεύτερη, τη μέθοδο των Χωρικών Ολοκληρωμάτων (Field Integral -FI- adjoint formulation), οι παράγωγοι εκφράζονται με συνδυασμό ολοκληρωμάτων τόσο στην επιφάνεια όσο και στο χώρο του υπολογιστικού χωρίου. Ανάμεσα στις δύο συζυγείς διατυπώσεις, τις SI και FI, παρατηρούνται αριθμητικές διαφορές στις υπολογιζόμενες παραγώγους οι οποίες και γίνονται ιδιαίτερα εμφανείς σε μη-επαρκώς πυκνά πλέγματα. Για αυτόν τον λόγο, επανεξετάζεται η ισοδυναμία τους, τόσο αναλυτικά όσο και αριθμητικά. Όπως αναμένονταν, αναλυτικά οι δύο προσεγγίσεις προκύπτουν ισοδύναμες. Η αριθμητική διαφορά τους οφείλεται στον λανθασμένο χειρισμό ενός όρου της SI διατύπωσης ο οποίος εκφράζει τη συμβολή των παραγώγων ευαισθησίας πλέγματος (grid sensitivities) στις παραγώγους της συνάρτησης-στόχου. Για την αντιμετώπισή του, προτείνεται ένας διαφορετικός χειρισμός του όρου αυτού, με τον οποίο είναι δυνατός ο σωστός υπολογισμός παραγώγων (σε συμφωνία με αυτές που υπολογίζονται από την FI διατύπωση), αυξάνοντας όμως αρκετά το υπολογιστικό κόστος, το οποίο γίνεται ίσο με αυτό της FI διατύπωσης. Καθώς αυτό δεν είναι επιθυμητό, ιδιαίτερα σε μεγάλες βιομηχανικές εφαρμογές, προτείνεται μία νέα συζυγής διατύπωση, η Εμπλουτισμένη διατύπωση των Επιφανειακών Ολοκληρωμάτων (Enhanced Surface Integral -E-SI- adjoint formulation), η οποία έχει την ακρίβεια της FI διατύπωσης αλλά το υπολογιστικό κόστος της SI. Ακόμα, η συζυγής διατύπωση επεκτείνεται για ροές σε βαθμίδες στροβιλομηχανών, λαμβάνοντας υπόψη την αλληλεπίδραση κινητής και σταθερής πτερύγωσης. Χρησιμοποιείται το μοντέλο Πολλαπλών Συστημάτων Αναφοράς (Multiple Reference Frame) με την παραδοχή του ≪παγωμένου δρομέα≫. Η συζυγής διατύπωση παρουσιάζεται για στρωτές ροές καθώς η γενίκευσή της σε τυρβώδεις ροές δεν παρουσιάζει καμία δυσκολία, εάν υπάρχει διαθέσιμο το εν χρήση συζυγές μοντέλο τύρβης. Στοχεύοντας μέγιστο βαθμό απόδοσης, η προτεινόμενη συζυγής διατύπωση εφαρμόζεται στη βελτιστοποίηση φυγοκεντρικής αντλίας. Για την επιτάχυνση και σταθεροποίηση της αριθμητικής λύσης των εξισώσεων ροής καθώς και των συζυγών αυτών, αναπτύσσεται η Μέθοδος των Αναδρομικών Προβολών (Recursive Projection Method, RPM) σε περιβάλλον OpenFoam©. Ιδιαίτερη έμφαση δίνεται στη χρήση της μεθόδου για τη σταθεροποίηση της επίλυσης του συζυγούς προβλήματος σε ροές που παρουσιάζουν μικρής κλίμακας χρονικές μεταβολές (οι οποίες συνήθως οδηγούν σε απόκλιση του συζυγούς προβλήματος). Μέσω της αναγνώρισης και της κατάλληλης αντιμετώπισης της κυρίαρχης ιδιοτιμής του συστήματος προς επίλυση, εξασφαλίζεται η σύγκλιση. Χρησιμοποιώντας τις μεθόδους που παρουσιάστηκαν παραπάνω, χαράζονται χάρτες ευαισθησίας σε γεωμετρίες πραγματικών αυτοκινήτων. Συγκεκριμένα, εξετάζονται το πρωτότυπο αυτοκίνητο L1 της VW και το επιβατικό μοντέλο A7 της AUDI. Σε αυτά, οι χάρτες ευαισθησίας που χαράζονται δείχνουν στον σχεδιαστή ποιες περιοχές κάθε γεωμετρίας έχουν δυνατότητες βελτίωσης και πως πρέπει να διαμορφωθούν για τη μείωση της οπισθέλκουσας. Η πληροφορία αυτή βοηθά στην αεροδυναμική βελτιστοποίηση της μορφής του αυτοκινήτου, χωρίς όμως να υλοποιεί αναγκαστικά βρόχο βελτιστοποίησης. Επιπλέον, αναπτύσσεται η συζυγής μέθοδος για τις χρονικά μη-μόνιμες Navier–Stokes εξισώσεις, για προβλήματα βελτιστοποίησης μορφής καθώς και ελέγχου ροής με χρήση δεσμών ρευστού. Χρησιμοποιείται η μέθοδος των σταθμών ελέγχου (checkpoints) για την αντιμετώπιση της αντίθετης-στον-χρόνο ολοκλήρωσης των χρονικά μη-μόνιμων συζυγών εξισώσεων. Ως προς τον έλεγχο της ροής, εξετάζονται ροές γύρω από κυλινδρικές γεωμετρίες οι οποίες ελέγχονται μέσω παλλόμενων δεσμών ρευστού. Σχετικά με τη βελτιστοποίηση μορφής, η ϕυγοκεντρική αντλία, η οποία ήδη εξετάστηκε υπό την παραδοχή του ≪παγωμένου δρομέα≫, επανεξετάζεται χρησιμοποιώντας χρονικά μη-μόνιμους επιλύτες ροής και συζυγούς προβλήματος.


Sign in / Sign up

Export Citation Format

Share Document