Decoding complex flow phenomena via Dynamic Mode Decomposition

In this thesis, the viability of the Dynamic Mode Decomposition (DMD) as a technique to analyze and model complex dynamic real-world systems is presented. This method derives, directly from data, computationally efficient reduced-order models (ROMs) which can replace too onerous or unavailable high-fidelity physics-based models. Optimizations and extensions to the standard implementation of the methodology are proposed, investigating diverse case studies related to the decoding of complex flow phenomena. The flexibility of this data-driven technique allows its application to high-fidelity fluid dynamics simulations, as well as time series of real systems observations. The resulting ROMs are tested against two tasks: (i) reduction of the storage requirements of high-fidelity simulations or observations; (ii) interpolation and extrapolation of missing data. The capabilities of DMD can also be exploited to alleviate the cost of onerous studies that require many simulations, such as uncertainty quantification analysis, especially when dealing with complex high-dimensional systems. In this context, a novel approach to address parameter variability issues when modeling systems with space and time-variant response is proposed. Specifically, DMD is merged with another model-reduction technique, namely the Polynomial Chaos Expansion, for uncertainty quantification purposes. Useful guidelines for DMD deployment result from the study, together with the demonstration of its potential to ease diagnosis and scenario analysis when complex flow processes are involved.

Polyethylene flow prediction with a differential multi-mode Pom-Pom model

We report the first steps of a collaborative project between the University of Queensland, Polyflow, Michelin, SK Chemicals, and RMIT University; on simulation, validation and application of a recently introduced constitutive model designed to describe branched polymers. Whereas much progress has been made on predicting the complex flow behaviour of many - in particular linear - polymers, it sometimes appears difficult to predict simultaneously shear thinning and extensional strain hardening behaviour using traditional constitutive models. Recently a new viscoelastic model based on molecular topology, was proposed by McLeish and Larson (1998). We explore the predictive power of a differential multi-mode version of the pom-pom model for the flow behaviour of two commercial polymer melts: a (long-chain branched) low-density polyethylene (LDPE) and a (linear) high-density polyethylene (HDPE). The model responses are compared to elongational recovery experiments published by Langouche and Debbaut (1999), and start-up of simple shear flow, stress relaxation after simple and reverse step strain experiments carried out in our laboratory.

Simulation of the velocity field in compound channel flow using different closure models

1st European IAHR Congress,6-4 May, Edinburg, Scotland

Quantification of aortic flow by phase-contrast magnetic resonance in patients with bicuspid aortic valve.

AIMS: Bicuspid aortic valve (BAV) causes complex flow patterns in the ascending aorta (AAo), which may compromise the accuracy of flow measurement by phase-contrast magnetic resonance (PC-MR). Therefore, we aimed to assess and compare the accuracy of forward flow measurement in the AAo, where complex flow is more dominant in BAV patients, with flow quantification in the left ventricular outflow tract (LVOT) and the aortic valve orifice (AV), where complex flow is less important, in BAV patients and controls. METHODS AND RESULTS: Flow was measured by PC-MR in 22 BAV patients and 20 controls at the following positions: (i) LVOT, (ii) AV, and (iii) AAo, and compared with the left ventricular stroke volume (LVSV). The correlation between the LVSV and the forward flow in the LVOT, the AV, and the AAo was good in BAV patients (r = 0.97/0.96/0.93; P < 0.01) and controls (r = 0.96/0.93/0.93; P < 0.01). However, in relation with the LVSV, the forward flow in the AAo was mildly underestimated in controls and much more in BAV patients [median (inter-quartile range): 9% (4%/15%) vs. 22% (8%/30%); P < 0.01]. This was not the case in the LVOT and the AV. The severity of flow underestimation in the AAo was associated with flow eccentricity. CONCLUSION: Flow measurement in the AAo leads to an underestimation of the forward flow in BAV patients. Measurement in the LVOT or the AV, where complex flow is less prominent, is an alternative means for quantifying the systolic forward flow in BAV patients.

Estimating traveltimes and groundwater flow patterns using 3D time-lapse crosshole ERT imaging of electrical resistivity fluctuations induced by infiltrating river water

The infiltration of river water into aquifers is of high relevance to drinking-water production and is a key driver of biogeochemical processes in the hyporheic and riparian zone, but the distribution and quantification of the infiltrating water are difficult to determine using conventional hydrological methods (e.g., borehole logging and tracer tests). By time-lapse inverting crosshole ERT (electrical resistivity tomography) monitoring data, we imaged groundwater flow patterns driven by river water infiltrating a perialpine gravel aquifer in northeastern Switzerland. This was possible because the electrical resistivity of the infiltrating water changed during rainfall-runoff events. Our time-lapse resistivity models indicated rather complex flow patterns as a result of spatially heterogeneous bank filtration and aquifer heterogeneity. The upper part of the aquifer was most affected by the river infiltrate, and the highest groundwater velocities and possible preferential flow occurred at shallow to intermediate depths. Time series of the reconstructed resistivity models matched groundwater electrical resistivity data recorded on borehole loggers in the upper and middle parts of the aquifer, whereas the resistivity models displayed smaller variations and delayed responses with respect to the logging data. in the lower part. This study demonstrated that crosshole ERT monitoring of natural electrical resistivity variations of river infiltrate could be used to image and quantify 3D bank filtration and aquifer dynamics at a high spatial resolution.

Numerical Investigation of Flow Past a Circular Cylinder and in a Staggered Tube Bundle Using Various Turbulence Models

Transitional flow past a three-dimensional circular cylinder is a widely studied phenomenon since this problem is of interest with respect to many technical applications. In the present work, the numerical simulation of flow past a circular cylinder, performed by using a commercial CFD code (ANSYS Fluent 12.1) with large eddy simulation (LES) and RANS (κ - ε and Shear-Stress Transport (SST) κ - ω! model) approaches. The turbulent flow for ReD = 1000 & 3900 is simulated to investigate the force coefficient, Strouhal number, flow separation angle, pressure distribution on cylinder and the complex three dimensional vortex shedding of the cylinder wake region. The numerical results extracted from these simulations have good agreement with the experimental data (Zdravkovich, 1997). Moreover, grid refinement and time-step influence have been examined. Numerical calculations of turbulent cross-flow in a staggered tube bundle continues to attract interest due to its importance in the engineering application as well as the fact that this complex flow represents a challenging problem for CFD. In the present work a time dependent simulation using κ – ε, κ - ω! and SST models are performed in two dimensional for a subcritical flow through a staggered tube bundle. The predicted turbulence statistics (mean and r.m.s velocities) have good agreement with the experimental data (S. Balabani, 1996). Turbulent quantities such as turbulent kinetic energy and dissipation rate are predicted using RANS models and compared with each other. The sensitivity of grid and time-step size have been analyzed. Model constants sensitivity study have been carried out by adopting κ – ε model. It has been observed that model constants are very sensitive to turbulence statistics and turbulent quantities.

Application of Computational Fluid Dynamics in Modeling Blood Flow in Human Thoracic Aorta

The aim of this study was to simulate blood flow in thoracic human aorta and understand the role of flow dynamics in the initialization and localization of atherosclerotic plaque in human thoracic aorta. The blood flow dynamics in idealized and realistic models of human thoracic aorta were numerically simulated in three idealized and two realistic thoracic aorta models. The idealized models of thoracic aorta were reconstructed with measurements available from literature, and the realistic models of thoracic aorta were constructed by image processing Computed Tomographic (CT) images. The CT images were made available by South Karelia Central Hospital in Lappeenranta. The reconstruction of thoracic aorta consisted of operations, such as contrast adjustment, image segmentations, and 3D surface rendering. Additional design operations were performed to make the aorta model compatible for the numerical method based computer code. The image processing and design operations were performed with specialized medical image processing software. Pulsatile pressure and velocity boundary conditions were deployed as inlet boundary conditions. The blood flow was assumed homogeneous and incompressible. The blood was assumed to be a Newtonian fluid. The simulations with idealized models of thoracic aorta were carried out with Finite Element Method based computer code, while the simulations with realistic models of thoracic aorta were carried out with Finite Volume Method based computer code. Simulations were carried out for four cardiac cycles. The distribution of flow, pressure and Wall Shear Stress (WSS) observed during the fourth cardiac cycle were extensively analyzed. The aim of carrying out the simulations with idealized model was to get an estimate of flow dynamics in a realistic aorta model. The motive behind the choice of three aorta models with distinct features was to understand the dependence of flow dynamics on aorta anatomy. Highly disturbed and nonuniform distribution of velocity and WSS was observed in aortic arch, near brachiocephalic, left common artery, and left subclavian artery. On the other hand, the WSS profiles at the roots of branches show significant differences with geometry variation of aorta and branches. The comparison of instantaneous WSS profiles revealed that the model with straight branching arteries had relatively lower WSS compared to that in the aorta model with curved branches. In addition to this, significant differences were observed in the spatial and temporal profiles of WSS, flow, and pressure. The study with idealized model was extended to study blood flow in thoracic aorta under the effects of hypertension and hypotension. One of the idealized aorta models was modified along with the boundary conditions to mimic the thoracic aorta under the effects of hypertension and hypotension. The results of simulations with realistic models extracted from CT scans demonstrated more realistic flow dynamics than that in the idealized models. During systole, the velocity in ascending aorta was skewed towards the outer wall of aortic arch. The flow develops secondary flow patterns as it moves downstream towards aortic arch. Unlike idealized models, the distribution of flow was nonplanar and heavily guided by the artery anatomy. Flow cavitation was observed in the aorta model which was imaged giving longer branches. This could not be properly observed in the model with imaging containing a shorter length for aortic branches. The flow circulation was also observed in the inner wall of the aortic arch. However, during the diastole, the flow profiles were almost flat and regular due the acceleration of flow at the inlet. The flow profiles were weakly turbulent during the flow reversal. The complex flow patterns caused a non-uniform distribution of WSS. High WSS was distributed at the junction of branches and aortic arch. Low WSS was distributed at the proximal part of the junction, while intermedium WSS was distributed in the distal part of the junction. The pulsatile nature of the inflow caused oscillating WSS at the branch entry region and inner curvature of aortic arch. Based on the WSS distribution in the realistic model, one of the aorta models was altered to induce artificial atherosclerotic plaque at the branch entry region and inner curvature of aortic arch. Atherosclerotic plaque causing 50% blockage of lumen was introduced in brachiocephalic artery, common carotid artery, left subclavian artery, and aortic arch. The aim of this part of the study was first to study the effect of stenosis on flow and WSS distribution, understand the effect of shape of atherosclerotic plaque on flow and WSS distribution, and finally to investigate the effect of lumen blockage severity on flow and WSS distributions. The results revealed that the distribution of WSS is significantly affected by plaque with mere 50% stenosis. The asymmetric shape of stenosis causes higher WSS in branching arteries than in the cases with symmetric plaque. The flow dynamics within thoracic aorta models has been extensively studied and reported here. The effects of pressure and arterial anatomy on the flow dynamic were investigated. The distribution of complex flow and WSS is correlated with the localization of atherosclerosis. With the available results we can conclude that the thoracic aorta, with complex anatomy is the most vulnerable artery for the localization and development of atherosclerosis. The flow dynamics and arterial anatomy play a role in the localization of atherosclerosis. The patient specific image based models can be used to diagnose the locations in the aorta vulnerable to the development of arterial diseases such as atherosclerosis.

Experimental analysis of the flow between stay and guide vanes of a pump-turbine in pumping mode

An experimental study was conducted in a pump-turbine model in pumping mode, in order to characterize the flow field structure in the region between stay and guide vanes, using mainly the laser-Doppler anemometry in a two-color and back-scattered light-based system. The structure of the steady and unsteady flow was analyzed. The measurements were performed at three operation points. The obtained data provide appropriate boundary conditions and a good base of validation for numerical codes, and for the understanding of main loss mechanisms of this complex flow.

The Inertio-Elastic Planar Entry Flow of Low-Viscosity Elastic Fluids in Micro-fabricated Geometries

The non-Newtonian flow of dilute aqueous polyethylene oxide (PEO) solutions through microfabricated planar abrupt contraction-expansions is investigated. The contraction geometries are fabricated from a high-resolution chrome mask and cross-linked PDMS gels using the tools of soft-lithography. The small length scales and high deformation rates in the contraction throat lead to significant extensional flow effects even with dilute polymer solutions having time constants on the order of milliseconds. The dimensionless extra pressure drop across the contraction increases by more than 200% and is accompanied by significant upstream vortex growth. Streak photography and videomicroscopy using epifluorescent particles shows that the flow ultimately becomes unstable and three-dimensional. The moderate Reynolds numbers (0.03 ≤ Re ≤ 44) associated with these high Deborah number (0 ≤ De ≤ 600) microfluidic flows results in the exploration of new regions of the Re-De parameter space in which the effects of both elasticity and inertia can be observed. Understanding such interactions will be increasingly important in microfluidic applications involving complex fluids and can best be interpreted in terms of the elasticity number, El = De/Re, which is independent of the flow kinematics and depends only on the fluid rheology and the characteristic size of the device.

Thermal performance of multipass parallel and counter-cross-flow heat exchangers

A thorough study of the thermal performance of multipass parallel cross-flow and counter-cross-flow heat exchangers has been carried out by applying a new numerical procedure. According to this procedure, the heat exchanger is discretized into small elements following the tube-side fluid circuits. Each element is itself a one-pass mixed-unmixed cross-flow heat exchanger. Simulated results have been validated through comparisons to results from analytical solutions for one- to four-pass, parallel cross-flow and counter-cross-flow arrangements. Very accurate results have been obtained over wide ranges of NTU (number of transfer units) and C* (heat capacity rate ratio) values. New effectiveness data for the aforementioned configurations and a higher number of tube passes is presented along with data for a complex flow configuration proposed elsewhere. The proposed procedure constitutes a useful research tool both for theoretical and experimental studies of cross-flow heat exchangers thermal performance.

A new approach for thermal performance calculation of cross-flow heat exchangers

A new numerical methodology for thermal performance calculation in cross-flow heat exchangers is developed. Effectiveness-number of transfer units (epsilon-NTU) data for several standard and complex flow arrangements are obtained using this methodology. The results are validated through comparison with analytical solutions for one-pass cross-flow heat exchangers with one to four rows and with approximate series solution for an unmixed-unmixed heat exchanger, obtaining in all cases very small errors. New effectiveness data for some complex configurations are provided. (c) 2005 Elsevier Ltd. All rights reserved.

Effectiveness-ntu computation with a mathematical model for cross-flow heat exchangers

Due to the wide range of design possibilities, simple manufactured, low maintenance and low cost, cross-flow heat exchangers are extensively used in the petroleum, petrochemical, air conditioning, food storage, and others industries. In this paper a mathematical model for cross-flow heat exchangers with complex flow arrangements for determining epsilon -NTU relations is presented. The model is based on the tube element approach, according to which the heat exchanger outlet temperatures are obtained by discretizing the coil along the tube fluid path. In each cross section of the element, tube-side fluid temperature is assumed to be constant because the heat capacity rate ratio C*=Cmin/Cmax tends toward zero in the element. Thus temperature is controlled by effectiveness of a local element corresponding to an evaporator or a condenser-type element. The model is validated through comparison with theoretical algebraic relations for single-pass cross-flow arrangements with one or more rows. Very small relative errors are obtained showing the accuracy of the present model. epsilon -NTU curves for several complex circuit arrangements are presented. The model developed represents a useful research tool for theoretical and experimental studies on heat exchangers performance.

Global instability analysis and control of three-dimensional long open cavity flow

The three-dimensional wall-bounded open cavity may be considered as a simplified geometry found in industrial applications such as leading gear or slotted flats on the airplane. Understanding the three-dimensional complex flow structure that surrounds this particular geometry is therefore of major industrial interest. At the light of the remarkable former investigations in this kind of flows, enough evidences suggest that the lateral walls have a great influence on the flow features and hence on their instability modes. Nevertheless, even though there is a large body of literature on cavity flows, most of them are based on the assumption that the flow is two-dimensional and spanwise-periodic. The flow over realistic open cavity should be considered. This thesis presents an investigation of three-dimensional wall-bounded open cavity with geometric ratio 6:2:1. To this aim, three-dimensional Direct Numerical Simulation (DNS) and global linear instability have been performed. Linear instability analysis reveals that the onset of the first instability in this open cavity is around Recr 1080. The three-dimensional shear layer mode with a complex structure is shown to be the most unstable mode. I t is noteworthy that the flow pattern of this high-frequency shear layer mode is similar to the observed unstable oscillations in supercritical unstable case. DNS of the cavity flow carried out at different Reynolds number from steady state until a nonlinear saturated state is obtained. The comparison of time histories of kinetic energy presents a clearly dominant energetic mode which shifts between low-frequency and highfrequency oscillation. A complete flow patterns from subcritical cases to supercritical case has been put in evidence. The flow structure at the supercritical case Re=1100 resembles typical wake-shedding instability oscillations with a lateral motion existed in the subcritical cases. Also, This flow pattern is similar to the observations in experiments. In order to validate the linear instability analysis results, the topology of the composite flow fields reconstructed by linear superposition of a three-dimensional base flow and its leading three-dimensional global eigenmodes has been studied. The instantaneous wall streamlines of those composited flows display distinguish influence region of each eigenmode. Attention has been focused on the leading high-frequency shear layer mode; the composite flow fields have been fully recognized with respect to the downstream wave shedding. The three-dimensional shear layer mode is shown to give rise to a typical wake-shedding instability with a lateral motions occurring downstream which is in good agreement with the experiment results. Moreover, the spanwise-periodic, open cavity with the same length to depth ratio has been also studied. The most unstable linear mode is different from the real three-dimensional cavity flow, because of the existence of the side walls. Structure sensitivity of the unstable global mode is analyzed in the flow control context. The adjoint-based sensitivity analysis has been employed to localized the receptivity region, where the flow is more sensible to momentum forcing and mass injection. Because of the non-normality of the linearized Navier-Stokes equations, the direct and adjoint field has a large spatial separation. The strongest sensitivity region is locate in the upstream lip of the three-dimensional cavity. This numerical finding is in agreement with experimental observations. Finally, a prototype of passive flow control strategy is applied.

Constriction flows of monodisperse linear entangled polymers: Multiscale modeling and flow visualization

We explore both the rheology and complex flow behavior of monodisperse polymer melts. Adequate quantities of monodisperse polymer were synthesized in order that both the materials rheology and microprocessing behavior could be established. In parallel, we employ a molecular theory for the polymer rheology that is suitable for comparison with experimental rheometric data and numerical simulation for microprocessing flows. The model is capable of matching both shear and extensional data with minimal parameter fitting. Experimental data for the processing behavior of monodisperse polymers are presented for the first time as flow birefringence and pressure difference data obtained using a Multipass Rheometer with an 11:1 constriction entry and exit flow. Matching of experimental processing data was obtained using the constitutive equation with the Lagrangian numerical solver, FLOWSOLVE. The results show the direct coupling between molecular constitutive response and macroscopic processing behavior, and differentiate flow effects that arise separately from orientation and stretch. (c) 2005 The Society of Rheology.

Variations of local heat transfer coefficient in piped flow of viscous liquids

Measurements were carried out to determine local coefficients of heat transfer in short lengths of horizontal pipe, and in the region of an discontinuity in pipe diameter. Laminar, transitional and turbulent flow regimes were investigated, and mixtures of propylene glycol and water were used in the experiments to give a range of viscous fluids. Theoretical and empirical analyses were implemented to find how the fundamental mechanism of forced convection was modified by the secondary effects of free convection, temperature dependent viscosity, and viscous dissipation. From experiments with the short tube it was possible to determine simple empirical relationships describing the axial distribution of the local 1usselt number and its dependence on the Reynolds and Prandtl numbers. Small corrections were made to account for the secondary effects mentioned above. Two different entrance configurations were investigated to demonstrate how conditions upstream could influence the heat transfer coefficients measured downstream In experiments with a sudden contraction in pipe diameter the distribution of local 1u3se1t number depended on the Prandtl number of the fluid in a complicated way. Graphical data is presented describing this dependence for a range of fluids indicating how the local Nusselt number varied with the diameter-ratio. Ratios up to 3.34:1 were considered. With a sudden divergence in pipe diameter, it was possible to derive the axial distribution of the local Nusse1t number for a range of Reynolds and Prandtl numbers in a similar way to the convergence experiments. Difficulty was encountered in explaining some of the measurements obtained at low Reynolds numbers, and flow visualization techniques wore used to determine the complex flow patterns which could lead to the anomalous results mentioned. Tests were carried out with divergences up to 1:3.34 to find the way in which the local Nusselt number varied with the diameter ratio, and a few experiments were carried out with very large ratios up .to 14.4. A limited amount of theoretical analysis of the 'divergence' system was carried out to substantiate certain explanations of the heat transfer mechanisms postulated.