Law of the wall in an unstably stratified turbulent channel flow

We perform direct numerical simulations of an unstably stratified turbulent channel flow to address the effects of buoyancy on the boundary layer dynamics and mean field quantities. We systematically span a range of parameters in the space of friction Reynolds number ($Re_{\tau}$) and Rayleigh number ($Ra$). Our focus is on deviations from the logarithmic law of the wall due to buoyant motion. The effects of convection in the relevant ranges are discussed providing measurements of mean profiles of velocity, temperature and Reynolds stresses as well as of the friction coefficient. A phenomenological model is proposed and shown to capture the observed deviations of the velocity profile in the log-law region from the non-convective case

Airframe Noise Simulations of a Full-Scale Aircraft

Computational results for a full-scale simulation of a Gulfstream G-III aircraft are presented. In support of a NASA airframe noise flight test campaign, Exa Corporations lattice Boltzmann PowerFLOW solver was used to perform time-accurate simulations of the flow around a highly detailed, full-scale aircraft model. Free-air boundary conditions were used at a Mach number of 0.23 and a Reynolds number of 10.5 10(exp 6) based on mean aerodynamic chord. This paper documents the simulation campaign for the baseline aircraft configuration at several flight conditions, including multiple flap deflections and main landing gear deployed or retracted. The high-fidelity, synthetic data were post-processed using a Ffowcs-Williams and Hawkings integral approach to estimate farfield acoustic behavior, with pressures on the model solid surface or a permeable surface enveloping the acoustic near field used as input. The numerical approach, simulation attributes, and the effects of grid resolution, gear deployment, and multiple flap deflections, are discussed as well

Finite Volume vs.vs. Streaming-based Lattice Boltzmann algorithm for fluid-dynamics simulations: a one-to-one accuracy and performance study

A new finite volume (FV) discretisation method for the Lattice Boltzmann (LB) equation which combines high accuracy with limited computational cost is presented. In order to assess the performance of the FV method we carry out a systematic comparison, focused on accuracy and computational performances, with the standard $streaming$ (ST) Lattice Boltzmann equation algorithm. To our knowledge such a systematic comparison has never been previously reported. In particular we aim at clarifying whether and in which conditions the proposed algorithm, and more generally any FV algorithm, can be taken as the method of choice in fluid-dynamics LB simulations. For this reason the comparative analysis is further extended to the case of realistic flows, in particular thermally driven flows in turbulent conditions. We report the first successful simulation of high-Rayleigh number convective flow performed by a Lattice Boltzmann FV based algorithm with wall grid refinement.Comment: 15 pages, 14 figures (discussion changes, improved figure readability

Synthetic turbulence generator for the wall-modeled LES lattice Boltzmann method

The synthetic turbulence generator (STG) lies at the interface of the Reynolds averaged Navier-Stokes (RANS) simulation and large eddy simulation (LES). This paper presents a STG for the multiple-relaxation-time(MRT) lattice Boltzmann method(LBM) framework at high friction Reynolds numbers, with consideration of near wall modeling. The Reichardt wall law, in combination with a force-based method, is used to model the near wall field. The STG wall-modeled(STG-WM) LES results are compared with turbulent channel flow simulations at $Re_{\tau}=1000,2000,5200$ at different resolutions. The results demonstrate good agreement with DNS, with the adaptation length of 6 to 8 boundary layer thickness. This method has a wide range of potentials for hybrid RANS/LES-LBM related applications at high friction Reynolds numbers

Lattice Boltzmann Methods for Turbulent Flows – Application to Coriolis Mass Flowmeter

Komplexe Strömungsphänomene machen es schwierig Ingenieursanwendungen so detailliert und genau zu simulieren, dass eine Charakterisierung und Verbesserung ihres Funktionsprinzips möglich ist. Diese Arbeit zeigt, dass die Lattice-Boltzmann-Methode (LBM) sehr gut für diesen Zweck geeignet ist. Im Vordergrund stehen hierbei die Simulation und Modellierung von turbulenten Strömungen. Diese lassen sich auf Grund der hervorragenden Parallelisierbarkeit der LBM mit Large-eddy Simulationen an Stelle von Reynolds-gemittelten Navier--Stokes Modellen, die im industriellen Umfeld üblich sind, berechnen. Somit können komplexe transiente turbulente Strömungen simulativ untersucht werden. Die daraus gewonnenen Erkenntnisse dienen insbesondere der Auslegung und Optimierung von Bauteilen und Prozessen. Alle beschriebenen LBM Simulationen werden mit der Open Source Software OpenLB durchgeführt. Dazu wird OpenLB erweitert, um eine Validierung von implementierten Turbulenzmodellen mittels kanonischer Strömungsformen zu ermöglichen. Des Weiteren wird ein Framework für die Simulation von Fluid-Struktur Interaktion (FSI) geschaffen. Anfangs werden die Kollisionsoperatoren Bhatnagar--Gross--Krook (BGK), Entropic Lattice Boltzmann (ELB), Two-Relaxation-Time (TRT), Regularized Lattice Boltzmann (RLB) und Multiple-Relaxation-Time (MRT) in der Taylor-Green Vortex Strömung, einem klassischen Beispiel für abklingende homogene isotrope Turbulenz (DHIT), untersucht. Hierbei liegt der Fokus auf Stabilität, Konsistenz und Genauigkeit der verwendeten Schemata. Die Studie beinhaltet den Vergleich der turbulenten kinetischen Energie, der Dissipationsrate der Energie und dem Energiespektrum zu einer Referenzlösung. Drei unterschiedliche Reynoldszahlen, $\mathrm{Re}=800$, $\mathrm{Re}=1600$ und $\mathrm{Re}=3000$, werden sowohl unter Verwendung einer akustischen als auch einer diffusiven Skalierung betrachtet, um den Einfluss der Lattice Machzahl zu charakterisieren. In stark unteraufgelösten Gitterkonfigurationen zeigt das BGK Schema ein instabiles Verhalten. Divergierende Simulationen unter der Verwendung des MRT Schemas sind auf eine starke Abhängigkeit von der Lattice Machzahl zurückzuführen. Obwohl ELB die Viskosität verändert, kann kein Verhalten, das einem Wirbelviskositätsmodell entspricht, gefunden werden. Bei geringen Lattice Machzahlen zeigt das RLB Schema sehr geringe Energielevel bei hohen Wellenzahlen. Der ,,magic parameter" des TRT Schemas wird bestimmt im Hinblick auf den Energieeintrag. Trotzdem wird keine erhöhte Stabilität im Vergleich zum BGK Schema festgestellt. Insgesamt sollte die Lattice Machzahl bezüglich des verwendeten Kollisonsschemas gewählt werden, um die Stabilität zu gewährleisten und die Genauigkeit zu verbessern. Für die Realisierung eines wandmodellierten Large-Eddy Simulation (NWM-LES) Ansatzes wird der BGK Kollisionsoperator ausgewählt. Das Smagorinsky Wirbelviskositätsmodell kommt hierbei zum Einsatz und wird in der turbulenten Grenzschicht mit der van Driest\u27schen Dämpfungsfunktion verwendet. Der Einfluss verschiedener Implementierungen von Geschwindigkeitsrandbedingungen und Wandfunktionen wird in einer biperiodischen, voll ausgebildeten turbulenten Kanalströmung für Schubspannungs-Reynoldszahlen von $\mathrm{Re}_\tau=1000$, $\mathrm{Re}_\tau=2000$ und $\mathrm{Re}_\tau=5200$ untersucht. Die Validierung erfolgt mittels Daten einer direkten numerischen Simulation (DNS) für Turbulenzstatistiken erster und zweiter Ordnung. Die Anwendung dieses Ansatzes auf einen Coriolis Massendurchflussmesser (CMF) zeigt, dass der Druckverlust bis zu einer Reynoldszahl $\mathrm{Re}=127800$ beschrieben werden kann. Des Weiteren wird der entwickelte NWM-LES LBM Ansatz mit OpenFOAM, einer Open Source Implementierung der finititen Volumen Methode (FVM) für komplexe turbulente Strömungen, die relevant für Verbrennungsmotoren sind, verglichen. Der zuvor entwickelte und validierte LBM Ansatz wird mit einer Geschwindigkeitsrandbedingung für gekrümmte Ränder erweitert. Die Ergebnisse beider Strömungslöser werden mit Daten eines Particle Image Velocimetry (PIV) Experiments verglichen. Die Validierung umfasst sowohl die zeitgemittelten als auch die quadratisch gemittelten (RMS) Geschwindigkeitsfelder. Zusätzlich wird sowohl die Laufzeit der Simulation als auch die Dauer der unterschiedlichen Gittergenerierungsprozesse bestimmt. Die Performanceanalyse der getesteten Konfiguration zeigt, dass OpenLB 32-mal schneller ist als OpenFOAM. Folglich ist der entwickelte NWM-LES LBM Ansatz dazu in der Lage, komplexe turbulente Strömungen in einer Ingenieursanwendung akkurat und mit einem verringerten Rechenaufwand zu beschreiben. Wirbel induzierte Vibrationen (VIV) sind ein weiterer wichtiger Anwendungsfall für Ingenieursapplikationen. Für die Untersuchung dieser werden verschiedene Fluid-Struktur Ansätze für LBM implementiert, verglichen und evaluiert. Die zwei untersuchten Klassen sind die Moving Boundary Methods (MBM) und die Partially Saturated Methods (PSM). Als erstes wird die Galiläische Invarianz von aerodynamischen Koeffizienten für die einzelnen Schemata untersucht. Dazu wird das BGK Schema verwendet, um einen exzentrisch positionierten Zylinder in einer Couette Strömung zu simulieren. Überdies werden verschiedene Volumenapproximationsmethoden für PSM und Auffüllmechanismen für MBM verglichen. Sowohl die Gitterkonvergenz als auch die Konvergenz der Galiläischen Invarianz werden betrachtet. Die Studie der VIV-Phänomene umfasst einen transvers oszillierenden Zylinder in einem Freistrom bei einer Reynoldszahl von $\mathrm{Re}=100$. Dabei werden freie und erzwungene Oszillation betrachtet, um bekannte Phänomene, wie Lock-in und Lock-out Zonen, zu untersuchen. Die Ergebnisse zeigen, dass sowohl MBM als auch PSM eine gute Übereinstimmung zu Literaturdaten aufweisen, womit die Eignung für VIV-Simulationen bestätigt werden kann. Schließlich wird ein Fluid-Struktur Interaktionsansatz unter der Verwendung eines MBM Ansatzes für die Simulation eines CMFs realisiert. Hierbei wird OpenLB mit Elmer, einer Open Source Implementierung der Finite-Elemente-Methode, gekoppelt, um auch die Strukturdynamik zu beschreiben. Ein gestaffelter Kopplungsansatz zwischen den beiden Softwarepaketen wird präsentiert. Das Finite-Elemente-Gitter wird durch das Gittergenerierungstool Gmsh erstellt, um einen kompletten Open Source Workflow zu garantieren. Zunächst werden die Eigenmoden des CMFs berechnet und mit Messdaten verglichen. Die daraus bestimmte Anregungsfrequenz wird zur Bestimmung des Phasenshifts in einer partitionierten voll gekoppelten FSI Simulation verwendet. Der berechnete Phasenshift zeigt eine gute Übereinstimmung mit den Messdaten und bestätigt, dass dieses Modell in der Lage ist, das Funktionsprinzip eines CMFs zu beschreiben. Die durchgeführten Studien zeigen das große Potential der LBM für die Simulation von Ingenieursapplikationen, insbesondere wenn turbulente Strömungen betrachtet werden