Showing papers in "arXiv: Fluid Dynamics in 2021"

Flow over an espresso cup: Inferring 3D velocity and pressure fields from tomographic background oriented schlieren videos via physics-informed neural networks

Abstract: Tomographic background oriented schlieren (Tomo-BOS) imaging measures density or temperature fields in 3D using multiple camera BOS projections, and is particularly useful for instantaneous flow visualizations of complex fluid dynamics problems. We propose a new method based on physics-informed neural networks (PINNs) to infer the full continuous 3D velocity and pressure fields from snapshots of 3D temperature fields obtained by Tomo-BOS imaging. PINNs seamlessly integrate the underlying physics of the observed fluid flow and the visualization data, hence enabling the inference of latent quantities using limited experimental data. In this hidden fluid mechanics paradigm, we train the neural network by minimizing a loss function composed of a data mismatch term and residual terms associated with the coupled Navier-Stokes and heat transfer equations. We first quantify the accuracy of the proposed method based on a 2D synthetic data set for buoyancy-driven flow, and subsequently apply it to the Tomo-BOS data set, where we are able to infer the instantaneous velocity and pressure fields of the flow over an espresso cup based only on the temperature field provided by the Tomo-BOS imaging. Moreover, we conduct an independent PIV experiment to validate the PINN inference for the unsteady velocity field at a center plane. To explain the observed flow physics, we also perform systematic PINN simulations at different Reynolds and Richardson numbers and quantify the variations in velocity and pressure fields. The results in this paper indicate that the proposed deep learning technique can become a promising direction in experimental fluid mechanics.

State-of-the-art SPH solver DualSPHysics: from fluid dynamics to multiphysics problems.

Abstract: DualSPHysics is a weakly compressible smoothed particle hydrodynamics (SPH) Navier-Stokes solver initially conceived to deal with coastal engineering problems, especially those related to wave impact with coastal structures. Since the first release back in 2011, DualSPHysics has shown to be robust and accurate for simulating extreme wave events along with a continuous improvement in efficiency thanks to the exploitation of hardware such as graphics processing units (GPUs) for scientific computing or the coupling with wave propagating models such as SWASH and OceanWave3D. Numerous additional functionalities have also been included in the DualSPHysics package over the last few years which allow the simulation of fluid-driven objects. The use of the discrete element method (DEM) has allowed the solver to simulate the interaction among different bodies (sliding rocks, for example), which provides a unique tool to analyse debris flows. In addition, the recent coupling with other solvers like Project Chrono or MoorDyn has been a milestone in the development of the solver. Project Chrono allows the simulation of articulated structures with joints, hinges, sliders and springs and MoorDyn allows simulating moored structures. Both functionalities make DualSPHysics one of the meshless model world leaders in the simulation of offshore energy harvesting devices. Lately, the present state of maturity of the solver goes beyond single phase simulations, allowing multi-phase simulations with gas-liquid and a combination of Newtonian and non-Newtonian models expanding further the capabilities and range of applications for the DualSPHysics solver. These advances and functionalities make DualSPHysics a state-of-the-art meshless solver with emphasis on free-surface flow modelling.

On closures for reduced order models $-$ A spectrum of first-principle to machine-learned avenues

Abstract: For over a century, reduced order models (ROMs) have been a fundamental discipline of theoretical fluid mechanics. Early examples include Galerkin models inspired by the Orr-Sommerfeld stability equation (1907) and numerous vortex models, of which the von Karman vortex street (1911) is one of the most prominent. Subsequent ROMs typically relied on first principles, like mathematical Galerkin models, weakly nonlinear stability theory, and two- and three-dimensional vortex models. Aubry et al. (1988) pioneered, with the proper orthogonal decomposition (POD) model of the turbulent boundary layer, a new data-driven paradigm, which has become the most popular avenue. In early POD modeling, available data was used to build an optimal basis, which was then utilized in a classical Galerkin procedure to construct the ROM. But data has made a profound impact on ROMs beyond the Galerkin expansion. In this paper, we take a modest step and illustrate the impact of data-driven modeling on one significant ROM area. Specifically, we focus on ROM closures, which are correction terms that are added to classical ROMs in order to model the effect of the discarded ROM modes in under-resolved simulations. Through simple examples, we illustrate the main modeling principles used to construct classical ROMs, motivate and introduce modern ROM closures, and show how data-driven modeling, artificial intelligence, and machine learning have changed the standard ROM methodology over the last two decades. Finally, we outline our vision on how state-of-the-art data-driven modeling can continue to reshape the field of reduced order modeling.

Physics-guided deep learning framework for predictive modeling of the Reynolds stress anisotropy

Abstract: Despite a cost-effective option in practical engineering, Reynolds-averaged Navier-Stokes simulations are facing the ever-growing demand for more accurate turbulence models. Recently, emerging machine learning techniques are making promising impact in turbulence modeling, but in their infancy for widespread industrial adoption. Towards this end, this work proposes a universal, inherently interpretable machine learning framework of turbulence modeling, which mainly consists of two parallel machine-learning-based modules to respectively infer the integrity basis and closure coefficients. At every phase of the model development, both data representing the evolution dynamics of turbulence and domain-knowledge representing prior physical considerations are properly fed and reasonably converted into modeling knowledge. Thus, the developed model is both data- and knowledge-driven. Specifically, a version with pre-constrained integrity basis is provided to demonstrate detailedly how to integrate domain-knowledge, how to design a fair and robust training strategy, and how to evaluate the data-driven model. Plain neural network and residual neural network as the building blocks in each module are compared. Emphases are made on three-fold: (i) a compact input feature parameterizing the newly-proposed turbulent timescale is introduced to release nonunique mappings between conventional input arguments and output Reynolds stress; (ii) the realizability limiter is developed to overcome under-constraint of modeled stress; and (iii) constraints of fairness and noisy-sensitivity are first included in the training procedure. In such endeavors, an invariant, realizable, unbiased and robust data-driven turbulence model is achieved, and does gain good generalization across channel flows at different Reynolds numbers and duct flows with various aspect ratios.

Uncovering near-wall blood flow from sparse data with physics-informed neural networks

Abstract: Near-wall blood flow and wall shear stress (WSS) regulate major forms of cardiovascular disease, yet they are challenging to quantify with high fidelity. Patient-specific computational and experimental measurement of WSS suffers from uncertainty, low resolution, and noise issues. Physics-informed neural networks (PINN) provide a flexible deep learning framework to integrate mathematical equations governing blood flow with measurement data. By leveraging knowledge about the governing equations (herein, Navier-Stokes), PINN overcomes the large data requirement in deep learning. In this study, it was shown how PINN could be used to improve WSS quantification in diseased arterial flows. Specifically, blood flow problems where the inlet and outlet boundary conditions were not known were solved by assimilating very few measurement points. Uncertainty in boundary conditions is a common feature in patient-specific computational fluid dynamics models. It was shown that PINN could use sparse velocity measurements away from the wall to quantify WSS with very high accuracy even without full knowledge of the boundary conditions. Examples in idealized stenosis and aneurysm models were considered demonstrating how partial knowledge about the flow physics could be combined with partial measurements to obtain accurate near-wall blood flow data. The proposed hybrid data-driven and physics-based deep learning framework has high potential in transforming high-fidelity near-wall hemodynamics modeling in cardiovascular disease.

From coarse wall measurements to turbulent velocity fields through deep learning

Abstract: This work evaluates the applicability of super-resolution generative adversarial networks (SRGANs) as a methodology for the reconstruction of turbulent-flow quantities from coarse wall measurements. The method is applied both for the resolution enhancement of wall fields and the estimation of wall-parallel velocity fields from coarse wall measurements of shear stress and pressure. The analysis has been carried out with a database of a turbulent open-channel flow with friction Reynolds number $Re_{\tau}=180$ generated through direct numerical simulation. Coarse wall measurements have been generated with three different downsampling factors $f_d=[4,8,16]$ from the high-resolution fields, and wall-parallel velocity fields have been reconstructed at four inner-scaled wall-normal distances $y^+=[15,30,50,100]$. We first show that SRGAN can be used to enhance the resolution of coarse wall measurements. If compared with direct reconstruction from the sole coarse wall measurements, SRGAN provides better instantaneous reconstructions, both in terms of mean-squared error and spectral-fractional error. Even though lower resolutions in the input wall data make it more challenging to achieve highly accurate predictions, the proposed SRGAN-based network yields very good reconstruction results. Furthermore, it is shown that even for the most challenging cases the SRGAN is capable of capturing the large-scale structures that populate the flow. The proposed novel methodology has great potential for closed-loop control applications relying on non-intrusive sensing.

Simulation of multi-species flow and heat transfer using physics-informed neural networks

Abstract: In the present work, single- and segregated-network PINN architectures are applied to predict momentum, species and temperature distributions of a dry air humidification problem in a simple 2D rectangular domain. The created PINN models account for variable fluid properties, species- and heat-diffusion and convection. Both the mentioned PINN architectures were trained using different hyperparameter settings, such as network width and depth to find the best-performing configuration. It is shown that the segregated-network PINN approach results in on-average 62% lower losses when compared to the single-network PINN architecture for the given problem. Furthermore, the single-network variant struggled to ensure species mass conservation in different areas of the computational domain, whereas, the segregated approach successfully maintained species conservation. The PINN predicted velocity, temperature and species profiles for a given set of boundary conditions were compared to results generated using OpenFOAM software. Both the single- and segregated-network PINN models produced accurate results for temperature and velocity profiles, with average percentage difference relative to the CFD results of approximately 7.5% for velocity and 8% for temperature. The mean error percentages for the species mass fractions are 9\% for the single-network model and 1.5% for the segregated-network approach. To showcase the applicability of PINNs for surrogate modelling of multi-species problems, a parameterised version of the segregated-network PINN is trained which could produce results for different water vapour inlet velocities. The normalised mean absolute percentage errors, relative to the OpenFOAM results, across three predicted cases for velocity and temperature are approximately 7.5% and 2.4% for water vapour mass fraction.

Convolutional neural networks for fluid flow analysis: toward effective metamodeling and low-dimensionalization

Abstract: We focus on a convolutional neural network (CNN), which has recently been utilized for fluid flow analyses, from the perspective on the influence of various operations inside it by considering some canonical regression problems with fluid flow data. We consider two types of CNN-based fluid flow analyses; 1. CNN metamodeling and 2. CNN autoencoder. For the first type of CNN with additional scalar inputs, which is one of the common forms of CNN for fluid flow analysis, we investigate the influence of input placements in the CNN training pipeline. As an example, estimation of force coefficients of an inclined flat plate and two side-by-side cylinders in laminar flows is considered. We find that care should be taken for the placement of additional scalar inputs depending on the problem setting and the complexity of flows that users handle. We then discuss the influence of various parameters and operations on the CNN performance, with the utilization of autoencoder (AE). A two-dimensional turbulence is considered for the demonstration of AE. The results of AE highly rely on the decaying nature. Investigation on the influence of padding operation at a convolutional layer is also performed. The zero padding shows reasonable ability compared to other methods which account for the boundary conditions assumed in the numerical data. Moreover, the effect of the dimensional reduction/extension methods inside CNN is also examined. The CNN model is robust against the difference in dimension reduction operations, while it is sensitive to the dimensional extension methods. The findings of this paper will help us better design a CNN architecture for practical fluid flow analysis

Towards extraction of orthogonal and parsimonious non-linear modes from turbulent flows.

Abstract: We propose a deep probabilistic-neural-network architecture for learning a minimal and near-orthogonal set of non-linear modes from high-fidelity turbulent-flow-field data useful for flow analysis, reduced-order modeling, and flow control. Our approach is based on $\beta$-variational autoencoders ($\beta$-VAEs) and convolutional neural networks (CNNs), which allow us to extract non-linear modes from multi-scale turbulent flows while encouraging the learning of independent latent variables and penalizing the size of the latent vector. Moreover, we introduce an algorithm for ordering VAE-based modes with respect to their contribution to the reconstruction. We apply this method for non-linear mode decomposition of the turbulent flow through a simplified urban environment, where the flow-field data is obtained based on well-resolved large-eddy simulations (LESs). We demonstrate that by constraining the shape of the latent space, it is possible to motivate the orthogonality and extract a set of parsimonious modes sufficient for high-quality reconstruction. Our results show the excellent performance of the method in the reconstruction against linear-theory-based decompositions. Moreover, we compare our method with available AE-based models. We show the ability of our approach in the extraction of near-orthogonal modes that may lead to interpretability.

Physics Informed Deep Learning for Flow and Transport in Porous Media

Abstract: We present our progress on the application of physics informed deep learning to reservoir simulation problems. The model is a neural network that is jointly trained to respect governing physical laws and match boundary conditions. The methodology is hereby used to simulate a 2-phase immiscible transport problem (Buckley-Leverett). The model is able to produce an accurate physical solution both in terms of shock and rarefaction and honors the governing partial differential equation along with initial and boundary conditions. We test various hypothesis (uniform and non-uniform initial conditions) and show that with the proper implementation of physical constraints, a robust solution can be trained within a reasonable amount of time and iterations. We revisit some of the limitations presented in previous work \cite{Fuks2020} and further the applicability of this method in a forward, pure hyperbolic setup. We also share some practical findings on the application of physics informed neural networks (PINN). We review various network architectures presented in the literature and show tips that helped improve their convergence and accuracy. The proposed methodology is a simple and elegant way to instill physical knowledge to machine-learning algorithms. This alleviates the two most significant shortcomings of machine-learning algorithms: the requirement for large datasets and the reliability of extrapolation. The principles presented can be generalized in innumerable ways in the future and should lead to a new class of algorithms to solve both forward and inverse physical problems.

Stable a posteriori LES of 2D turbulence using convolutional neural networks: Backscattering analysis and generalization to higher Re via transfer learning

Abstract: There is a growing interest in developing data-driven subgrid-scale (SGS) models for large-eddy simulation (LES) using machine learning (ML). In a priori (offline) tests, some recent studies have found ML-based data-driven SGS models that are trained on high-fidelity data (e.g., from direct numerical simulation, DNS) to outperform baseline physics-based models and accurately capture the inter-scale transfers, both forward (diffusion) and backscatter. While promising, instabilities in a posteriori (online) tests and inabilities to generalize to a different flow (e.g., with a higher Reynolds number, Re) remain as major obstacles in broadening the applications of such data-driven SGS models. For example, many of the same aforementioned studies have found instabilities that required often ad-hoc remedies to stabilize the LES at the expense of reducing accuracy. Here, using 2D decaying turbulence as the testbed, we show that deep fully convolutional neural networks (CNNs) can accurately predict the SGS forcing terms and the inter-scale transfers in a priori tests, and if trained with enough samples, lead to stable and accurate a posteriori LES-CNN. Further analysis attributes these instabilities to the disproportionally lower accuracy of the CNNs in capturing backscattering when the training set is small. We also show that transfer learning, which involves re-training the CNN with a small amount of data (e.g., 1%) from the new flow, enables accurate and stable a posteriori LES-CNN for flows with 16x higher Re (as well as higher grid resolution if needed). These results show the promise of CNNs with transfer learning to provide stable, accurate, and generalizable LES for practical use.

On the role of nonlinear correlations in reduced-order modeling

Abstract: A major goal for reduced-order models of unsteady fluid flows is to uncover and exploit latent low-dimensional structure. Proper orthogonal decomposition (POD) provides an energy-optimal linear basis to represent the flow kinematics, but converges slowly for advection-dominated flows and tends to overestimate the number of dynamically relevant variables. We show that nonlinear correlations in the temporal POD coefficients can be exploited to identify the underlying attractor, characterized by a minimal set of driving modes and a manifold equation for the remaining modes. By viewing these nonlinear correlations as an invariant manifold reduction, this least-order representation can be used to stabilize POD-Galerkin models or as a state space for data-driven model identification. In the latter case, we use sparse polynomial regression to learn a compact, interpretable dynamical system model from the time series of the active modal coefficients. We demonstrate this perspective on a quasiperiodic shear-driven cavity flow and show that the dynamics evolve on a torus generated by two independent Stuart-Landau oscillators. These results emphasize the importance of nonlinear dimensionality reduction to reveal underlying structure in complex flows.

Reinforcement-learning-based control of confined cylinder wakes with stability analyses

Abstract: This work studies the application of a reinforcement-learning-based (RL) flow control strategy to the flow past a cylinder confined between two walls in order to suppress vortex shedding. The control action is blowing and suction of two synthetic jets on the cylinder. The theme of this study is to investigate how to use and embed physical information of the flow in the RL-based control. First, global linear stability and sensitivity analyses based on the time-mean flow and the steady flow (which is a solution to the Navier-Stokes equations) are conducted in a range of blockage ratios and Reynolds numbers. It is found that the most sensitive region in the wake extends itself when either parameter increases in the parameter range we investigated here. Then, we utilise these physical results to help design RL-based control policies. We find that the controlled wake converges to the unstable steady base flow, where the vortex shedding can be successfully suppressed. A persistent oscillating control seems necessary to maintain this unstable state. The RL algorithm is able to outperform a gradient-based optimisation method (optimised in a certain period of time) in the long run. Furthermore, when the flow stability information is embedded in the reward function to penalise the instability, the controlled flow may become more stable. Finally, according to the sensitivity analyses, the control is most efficient when the probes are placed in the most sensitive region. The control can be successful even when few probes are properly placed in this manner.

Data-Driven Approach for Noise Reduction in Pressure-Sensitive Paint Data Based on Modal Expansion and Time-Series Data at Optimally Placed Points

Abstract: We propose a noise reduction method for unsteady pressure-sensitive paint (PSP) data based on modal expansion, the coefficients of which are determined from time-series data at optimally placed points. In this study, the proper orthogonal decomposition (POD) mode calculated from the time-series PSP data is used as a modal basis. Based on the POD modes, the points that effectively represent the features of the pressure distribution are optimally placed by the sensor optimization technique. Then, the time-dependent coefficient vector of the POD modes is determined by minimizing the difference between the time-series pressure data and the reconstructed pressure at the optimal points. Here, the coefficient vector is assumed to be a sparse vector. The advantage of the proposed method is a self-contained method, while existing methods use other data such as pressure tap data for the reduction of the noise. As a demonstration, we applied the proposed method to the PSP data measuring the Karman vortex street behind a square cylinder. The reconstructed pressure data is agreed very well with the pressures independently measured by pressure transducers.

Shock and shear layer interaction in a confined supersonic cavity flow

Abstract: The impinging shock of varying strengths on the free shear layer in a confined supersonic cavity flow is studied numerically using the detached-eddy simulation. The resulting spatiotemporal variations are analyzed between the different cases using unsteady statistics, $x-t$ diagrams, spectral analysis, and modal decomposition. A cavity of length to depth ratio $[L/D]=2$ at a freestream Mach number of $M_\infty = 1.71$ is considered to be in a confined passage. Impinging shock strength is controlled by changing the ramp angle ($\theta$) on the top-wall. The static pressure ratio across the impinging shock ($p_2/p_1$) is used to quantify the impinging shock strength. Five different impinging shock strengths are studied by changing the pressure ratio: $1.0,1.2,1.5,1.7$ and $2.0$. As the pressure ratio increases from 1.0 to 2.0, the cavity wall experiences a maximum pressure of 25% due to shock loading. At [$p_2/p_1]=1.5$, fundamental fluidic mode or Rossiter's frequency corresponding to $n=1$ mode vanishes whereas frequencies correspond to higher modes ($n=2$ and $4$) resonate. Wavefronts interaction from the longitudinal reflections inside the cavity with the transverse disturbances from the shock-shear layer interactions is identified to drive the strong resonant behavior. Due to Mach-reflections inside the confined passage at $[p_2/p_1]=2.0$, shock-cavity resonance is lost. Based on the present findings, an idea to use a shock-laden confined cavity flow in an enclosed supersonic wall-jet configuration as passive flow control or a fluidic device is also demonstrated.

Exact tensor closures for the three dimensional Jeffery's equation

Abstract: This paper presents an exact formula for calculating the fourth-moment tensor from the second-moment tensor for the three dimensional Jeffery's equation. Although this approach falls within the category of a moment tensor closure, it does not rely upon an approximation, either analytic or curve fit, of the fourth-moment tensor as do previous closures. This closure is orthotropic in the sense of \cite{cintra:95}, or equivalently, a natural closure in the sense of \cite{verleye:93}. The existence of these explicit formulae has been asserted previously, but as far as the authors know, the explicit forms have yet to be published. The formulae involve elliptic integrals, and are valid whenever fiber orientation was isotropic at some point in time. Finally, this paper presents the Fast Exact Closure (FEC), a fast and in principle exact method for solving Jeffery's equation, which does not require approximate closures, nor the elliptic integral computation.

Role of Jet Spacing and Strut Geometry on the formation of Large Scale Structures & Mixing Characteristics.

Abstract: The present study primarily focuses on the effect of jet spacing and strut geometry on the evolution and structure of the large-scale vortices which play a key role in mixing characteristics in turbulent supersonic flows. Numerically simulated results corresponding to varying parameters such as strut geometry and jet spacing (Xn=nDj such that n=2, 3 & 5) for square jet of height Dj=0.6 mm is presented in the current study, while the work also investigates the presence of the local quasi-two-dimensionality for the X2(2Dj) jet spacing; however, the same is not true for higher jet spacing. Further, the tapered strut (TS) section is modified to the straight strut (SS) for investigation, where the remarkable difference in flow physics is unfolded between the two configurations for similar jet spacing (X2: 2Dj). The instantaneous density and vorticity contours reveal the structures of varying scales undergoing different evolution for the different configurations. The effect of local spanwise rollers is clearly manifested in the mixing efficiency and the jet spreading rate. SS configuration exhibits excellent near field mixing behavior amongst all the arrangements. However, in case of TS cases, only X2(2Dj) configuration performs better due to the presence of local spanwise rollers. The qualitative and quantitative analysis reveals that near-field mixing is strongly affected by the two-dimensional rollers, while the early onset of wake mode is another crucial parameter to have improved mixing. Modal decomposition performed for SS arrangement sheds light into the spatial and temporal coherence of the structures, where the most dominant structures are found to be the Von-Karman street vortices in the wake region.

Comprehensive analysis of spherical bubble oscillations and shock wave emission in laser-induced cavitation.

Abstract: The dynamics of spherical laser-induced cavitation bubbles in water is investigated by plasma photography, time-resolved shadowgraphs, and single-shot probe beam scattering enabling to portray the transition from initial nonlinear to late linear oscillations. The frequency of late oscillations yields the bubble's gas content. Simulations with the Gilmore model using plasma size as input and oscillation times as fit parameter provide insights into experimentally not accessible bubble parameters and shock wave emission. The model is extended by a term covering the initial shock-driven acceleration of the bubble wall, an automated method determining shock front position and pressure decay, and an energy balance for the partitioning of absorbed laser energy into vaporization, bubble and shock wave energy, and dissipation through viscosity and condensation. These tools are used for analysing a scattering signal covering 102 oscillation cycles. The bubble was produced by a plasma with 1550 K average temperature and had 36 $\mu$m maximum radius. Predicted bubble wall velocities during expansion agree well with experimental data. Upon first collapse, most energy was stored in the compressed liquid around the bubble and then radiated away acoustically. The collapsed bubble contained more vapour than gas, and its pressure was 13.5 GPa. The pressure of the rebound shock wave initially decayed $\propto r^{-1.8}$, and energy dissipation at the shock front heated liquid near the bubble wall above the superheat limit. The shock-induced temperature rise reduces damping during late bubble oscillations. Bubble dynamics changes significantly for small bubbles with less than 10 $\mu$m radius.

High-fidelity reconstruction of turbulent flow from spatially limited data using enhanced super-resolution generative adversarial network

Abstract: In this study, a deep learning-based approach is applied with the aim of reconstructing high-resolution turbulent flow fields using minimal flow fields data. A multi-scale enhanced super-resolution generative adversarial network with a physics-based loss function is introduced as a model to reconstruct the high-resolution flow fields. The model capability to reconstruct high-resolution laminar flows is examined using data of laminar flow around a square cylinder. The results reveal that the model can accurately reproduce the high-resolution flow fields even when limited spatial information is provided. The case of turbulent channel flow is used to assess the ability of the model to reconstruct the high-resolution wall-bounded turbulent flow fields. The instantaneous and statistical results obtained from the model agree well with the ground truth data, indicating that the model can successfully learn to map the coarse flow fields to the high-resolution once. Furthermore, the computational cost of the proposed model, which is examined carefully, is found to be effectively low. This demonstrates that using high-fidelity training data with physics-guided generative adversarial network-based models can be practically efficient in reconstructing high-resolution turbulent flow fields from extremely coarse data.

A Quantitative Analysis of the Ignition Characteristics of Fine Iron Particles

Abstract: Ignition of iron particles in an oxidizing environment marks the onset of a self-sustained combustion. The objective of the current study is to quantitatively examine the ignition characteristics of fine iron particles governed by the kinetics of solid-phase iron oxidation. The oxidation rates are inversely proportional to the thickness of the oxide layer and calibrated using the experimentally measured growth of iron-oxide layers over time. Steady-state and unsteady analysis have been performed to probe the dependence of the critical gas temperature required to trigger a thermal runaway (namely, the ignition temperature $T_\mathrm{ign}$) on particle size, initial thickness of oxide layer, inert gas species, radiative heat loss, and the collective heating effect in a suspension of particles. Both analyses indicate that $T_\mathrm{ign}$ depends on $\delta_0$, i.e., the ratio between the initial oxide layer thickness and particle size, regardless of the absolute size of the particle. The unsteady analysis predicts that, for $\delta_0 \lesssim 0.003$, $T_\mathrm{ign}$ becomes independent of $\delta_0$. Under standard conditions in air, $T_\mathrm{ign}$ is approximately 1080 K for any particle size greater than 5 microns. Radiative heat loss has a minor effect on $T_\mathrm{ign}$. The collective effect of a suspension of iron particles in reducing $T_\mathrm{ign}$ is demonstrated. The transition behavior between kinetic-controlled and external-diffusion-controlled combustion regimes of an ignited iron particle is systematically examined. The influences of initial oxide-layer thickness and particle temperature on the ignition delay time, $\tau_\mathrm{ign}$, of iron particles are parametrically probed. A $d^2$-law scaling between $\tau_\mathrm{ign}$ and particle size is identified.

Passive control of the onset of vortex shedding in flow past a circular cylinder using Slit

Abstract: The present study investigates the passive flow control phenomena over a two-dimensional circular cylinder using numerical simulations in the laminar regime. The aim is to explore one of the passive control techniques, which involves the introduction of a slit to the geometry of the cylinder. The two parameters, slit width ratio S/D (slit width/diameter) and slit angle (measured with respect to the incoming flow direction), play an essential role in determining the trend of critical Reynolds number (Rec). Most of the analysis invokes flow visualization and saturation amplitude methods to obtain the critical Reynolds number (indicative of the onset of vortex shedding) for different cases. Further, Hopf bifurcation analysis using Stuart-Landau equation and global stability analysis confirm the accuracy and consistency in the predicted solutions. The additional amount of flow through slit increases the pressure downstream of the cylinder, which consequently leads to an increase in Rec of the modified cylinder. The critical Reynolds number increases with S/D of the modified cylinder at 0 deg slit angle (as an additional amount of flow grows with S/D). The critical Reynolds number shows an increasing trend with the slit angle in the range of S/D (0.05-0.15) as the fluctuation intensity reduces with slit angle in this range. For S/D= 0.15-0.25, the extra amount of flow through slit induces the instability in wake which causes a decrease in Rec with slit angle. A correlation is obtained that estimates the critical Reynolds number for a given S/D and slit angle.

Electrokinetic oscillatory flow and energy conversion of viscoelastic fluids in microchannels: a linear analysis

Abstract: We study the electrokinetic flow of viscoelastic fluids subjected to an oscillatory pressure gradient, and particularly focus on the resonance behaviors in the flow. The governing equations are restricted to linear regime so that the velocity and streaming potential fields can be solved analytically. Based on the interaction of viscoelastic shear waves, we explain the mechanism of resonance, and derive a critical Deborah number Dec = 1/4 which dictates the occurrence of resonance. Using the Maxwell fluid model, we show that the resonance enhances electrokinetic effects and results in a dramatic increase of electrokinetic energy conversion efficiency. However, by applying the Oldroyd-B fluid model it reveals that the amplification of efficiency is suppressed even for a very small Newtonian solvent contribution. This may be one of the reasons that experimental verification regarding the high efficiency predicted by Bandopadhyay & Chakraborty (Appl. Phys. Lett., vol. 101, 2012, 043905) is unavailable in the literature. Furthermore, the damping effect of solvent viscosity is more significant for higher-order resonances. Introducing the factor of multiple relaxation times, we show that the occurrence of resonances for the streaming potential field and the flow rate are still dominated by Dec. For the efficiency in the multi-mode case, the occurrence of resonance is dominated by the Deborah number De and the mode number N, and the resonance disappears for small De or large N. In addition, a new type of scaling relation between the streaming potential field and EDL thickness can be identified at large De.

Attention-Enhanced Neural Network Models for Turbulence Simulation

Abstract: Deep neural network models have shown a great potential in accelerating the simulation of fluid dynamic systems. Once trained, these models can make inference within seconds, thus can be extremely efficient. However, they suffer from a generalization problem when the flow becomes chaotic and turbulent. One of the most important reasons is that, existing models lack the mechanism to handle the unique characteristic of turbulent flow: multi-scale flow structures are non-uniformly distributed and strongly nonequilibrium. In this work, we address this issue with the concept of visual attention: intuitively, we expect the attention module to capture the nonequilibrium of turbulence by automatically adjusting weights on different regions. We benchmark the performance improvement with a state of the art neural network model, the Fourier Neural Operator (FNO), on two-dimensional (2D) turbulence prediction task. Numerical experiments show that the attention-enhanced neural network model can generalize well on higher Reynolds numbers flow, and can accurately reconstruct a variety of statistics and instantaneous spatial structures of turbulence. The attention mechanism provides 40% error reduction with 1% increase of parameters, at the same level of computational cost.

Role of corner flow separation in unsteady dynamics of hypersonic flow over a double wedge geometry

Abstract: This study investigates the origin and sustenance of self induced oscillations of shock structures in a hypersonic flow over a double wedge configuration. Previously, various researchers have considered the double wedge flow configuration for inviscid flow with variations of different inflows as well as geometric parameters such as inflow Mach number , wedge angles, and wedge lengths. Few recently published articles reveal an unsteady flow physics involved with the hypersonic viscous flow for double wedge configuration with large second wedge angles. However, the reason for such self sustained flow oscillations is not completely clear. The present work seeks out to investigate the origin of such oscillations in a low enthalpy hypersonic flow with different aft wedge angles and wedge length this http URL the current study, viscous flow over a double wedge at Mach 7 and fore wedge angle of 30° is considered. An improved version of rhoCentralFoam solver in OpenFOAM is used to investigate the double wedge flow over different grid resolutions in the separation region and shock shock interaction region. This study corroborates the observation from the previous literature with an improvement in the range of parameters which results in a self sustained periodic oscillation. The present study also suggests that the unsteadiness becomes possible when the incidence shock is in the proximity of the aft wedge expansion corner as a consequence of different wedge length ratios(L1/L2) or aft wedge angles ({\theta}2). Flow can still be steady at a large aft wedge angle if the incidence shock is far ahead of the aft wedge expansion corner.

Weakly nonlinear analysis of the viscoelastic instability in channel flow for finite and vanishing Reynolds numbers

Abstract: The recently-discovered centre-mode instability of rectilinear viscoelastic shear flow (Garg et al. Phy. Rev. Lett. 121, 024502, 2018) has offered an explanation for the origin of elasto-inertial turbulence (EIT) which occurs at lower Weissenberg ($Wi$) numbers. In support of this, we show using weakly nonlinear analysis that the subcriticality found in Page et al. (Phys. Rev. Lett. 125, 154501, 2020) is generic across the neutral curve with the instability only becoming supercritical at low Reynolds ($Re$) numbers and high $Wi$. We demonstrate that the instability can be viewed as purely elastic in origin even for $Re=O(10^3)$, rather than `elasto-inertial', as the underlying shear does not energise the instability. It is also found that the introduction of a realistic maximum polymer extension length, $L_{max}$, in the FENE-P model moves the neutral curve closer to the inertialess $Re=0$ limit at a fixed ratio of solvent-to-solution viscosities, $\beta$. In the dilute limit ($\beta \rightarrow 1$) with $L_{max} =O(100)$, the linear instability can brought down to more physically-relevant $Wi\gtrsim 110$ at $\beta=0.98$, compared with the threshold $Wi=O(10^3)$ at $\beta=0.994$ reported recently by Khalid et al. (arXiv: 2103.06794) for an Oldroyd-B fluid. Again the instability is subcritical implying that inertialess rectilinear viscoelastic shear flow is nonlinearly unstable - i.e. unstable to finite amplitude disturbances - for even lower $Wi$.

Characterization of turbulent supersonic flow over a backward facing step through POD

Abstract: The present work reports on the flow physics of turbulent supersonic flow over backward facing step (BFS) at Mach 2 using LES methodology where the dynamic Smagorinsky model is used for SGS modeling, while POD is invoked to identify the coherent structures present in the flow. The mean data obtained through the computations is in good agreement with the experimental measurements, while the iso-surfaces of Q-criterion at different time instants show the complex flow structures. The presence of counter rotating vortex pair in the shear layer along with the complex shock wave/boundary layer interaction leading to the separation of boundary layer is also evident from the contours of both Q and the modulus of vorticity. Further, the POD analysis reveals the presence of coherent structures, where the first and second modes confirm the vortical structures near the step as well as along the shear layer in the downstream region; while the second, third and fourth modes confirm the presence of vortices along the shear layer due to Kelvin-Helmholtz (K-H) instability. Moreover, POD as well as frequency analysis is extended at different planes to extract the detailed flow features.

Performance and wake characteristics of tidal turbines in an infinitely large array

Abstract: The efficiency of tidal-stream turbines in a large array depends on the balance between negative effects of turbine-wake interactions and positive effects of bypass-flow acceleration due to local blockage, both of which are functions of the layout of turbines. Here we investigate the hydrodynamics of turbines in an infinitely large array with aligned or staggered layouts for a range of streamwise and lateral turbine spacing. First, we present a theoretical analysis based on an extension of the Linear Momentum Actuator Disc Theory (LMADT) for perfectly aligned and staggered layouts, employing a hybrid inviscid-viscous approach to account for the local blockage effect within each turbine row and the viscous (turbulent) wake mixing behind each row in a coupled manner. We perform Large-Eddy Simulation (LES) of open-channel flow for 28 layouts of tidal turbines using an Actuator-Line Method with periodic boundary conditions. Both theoretical and LES results show that the efficiency of turbines (or the power of turbines for a given bulk velocity) in an aligned array decreases as we reduce the streamwise spacing, whereas that in a staggered array remains high and may even increase due to the positive local blockage effect (causing the local flow velocity upstream of each turbine to exceed the bulk velocity) if the lateral turbine spacing is sufficiently small. The LES results further reveal that the amplitude of wake meandering tends to decrease as we reduce the lateral spacing, which leads to a lower wake recovery rate in the near-wake region. These results will help to understand and improve the efficiency of tidal turbines in future large arrays, even though the performance of real tidal arrays may depend not only on turbine-to-turbine interactions but also on macro-scale interactions between the array and natural tidal currents, the latter of which are outside the scope of this study.

Modulation of homogeneous and isotropic turbulence in emulsions

Abstract: We present a numerical study of emulsions in homogeneous and isotropic turbulence at $Re_\lambda=137$. The problem is addressed via Direct Numerical Simulations (DNS), where the Volume of Fluid (VOF) is used to represent the complex features of the liquid-liquid interface. We consider a mixture of two iso-density fluids, where fluid properties are varied with the goal of understanding their role in turbulence modulation, in particular the volume fraction ($0.03<\alpha<0.5$), viscosity ratio ($0.01<\mu_d/\mu_c<100$) and large scale Weber number ($10.6

Experimental and numerical evidence of intensified non-linearity at the micro and nano scale: The lipid coated acoustic bubble.

Abstract: A lipid coated bubble (LCB) oscillator is a very interesting non-smooth oscillator with many important applications ranging from industry and chemistry to medicine. However, due to the complex behavior of the coating intermixed with the nonlinear behavior of the bubble itself, the dynamics of the LCB are not well understood. In this work, lipid coated Definity microbubbles (MBs) were sonicated with 25 MHz 30 cycle pulses with pressure amplitudes between 70kPa-300kPa. Here, we report higher order subharmonics in the scattered signals of single MBs at low amplitude high frequency ultrasound excitations. Experimental observations reveal the generation of period 2(P2), P3, and two different P4 oscillations at low excitation amplitude. Despite the reduced damping of the uncoated bubble system, such enhanced nonlinear oscillations has not been observed and can not be theoretically explained for the uncoated bubble. To investigate the mechanism of the enhanced nonlinearity, the bifurcation structure of the lipid coated MBs is studied for a wide range of MBs sizes and shell parameters. Consistent with the experimental results, we show that this unique oscillator can exhibit chaotic oscillations and higher order subharmonics at excitation amplitudes considerably below those predicted by the uncoated oscillator. Buckling or rupture of the shell and the dynamic variation of the shell elasticity causes the intensified non-linearity at low excitations. The simulated scattered pressure by single MBs are in good agreement with the experimental signals.

Effect of low-level jet height on wind farm performance

Abstract: Low-level jets (LLJs) are the wind maxima in the lowest 50 to 1000 m of atmospheric boundary layers. Due to their significant influence on the power production of wind farms it is crucial to understand the interaction between LLJs and wind farms. In the presence of an LLJ, there are positive and negative shear regions in the velocity profile. The positive shear regions of LLJs are continuously turbulent, while the negative shear regions have limited turbulence. We present large-eddy simulations of wind farms in which the LLJ is above, below, or in the middle of the turbine rotor swept area. We find that the wakes recover relatively fast when the LLJ is above the turbines. This is due to the high turbulence below the LLJ and the downward vertical entrainment created by the momentum deficit due to the wind farm power production. This harvests the jet's energy and aids wake recovery. However, when the LLJ is below the turbine rotor swept area, the wake recovery is very slow due to the low atmospheric turbulence above the LLJ. The energy budget analysis reveals that the entrainment fluxes are maximum and minimum when the LLJ is above and in the middle of the turbine rotor swept area, respectively. Surprisingly, we find that the negative shear creates a significant entrainment flux upward when the LLJ is below the turbine rotor swept area. This facilitates energy extraction from the jet, which is beneficial for the performance of downwind turbines.