Sloshing and slamming oscillations in collapsible channel flow

We consider laminar high-Reynolds-number flow through a finite-length planar channel, where a portion of one wall is replaced by a thin massless elastic membrane that is held under longitudinal tension T and subject to an external pressure distribution. The flow is driven by a fixed pressure drop along the full length of the channel. We investigate the global stability of two-dimensional Poiseuille flow using a method of matched local eigenfunction expansions, which is compared to direct numerical simulations. We trace the neutral stability curve of the primary oscillatory instability of the system, illustrating a transition from high-frequency ‘sloshing’ oscillations at high T to vigorous ‘slamming’ motion at low T . Small-amplitude sloshing at high T can be captured using a low-order eigenmode truncation involving four surface-based modes in the compliant segment of the channel coupled to Womersley flow in the rigid segments. At lower tensions, we show that hydrodynamic modes contribute increasingly to the global instability and we demonstrate a change in the mechanism of energy transfer from the mean flow, with viscous effects being destabilising. Simulations of finite-amplitude oscillations at low T reveal a generic slamming motion, in which the the flexible membrane is drawn close to the opposite rigid wall before rapidly recovering. A simple model is used to demonstrate how fluid inertia in the downstream rigid channel segment, coupled to membrane curvature downstream of the moving constriction, together control slamming dynamics.

Fluidic flow delay by ionogel passive pumps in microfluidic paper-based analytical devices

[EN] A new concept for fluid flow manipulation in microfluidic paper-based analytical devices (m-PADs) is presented by introducing ionogel materials as passive pumps. m-PADs were fabricated using a new doubleside contact stamping process and ionogels were precisely photopolymerised at the inlet of the m-PADs.The ionogels remain mainly on the surface of the paper and get absorbed in the superficial paper-fibers allowing for the liquid to flow from the ionogel into the paper easily. As a proof of concept the fluid flowand mixing behaviour of two different ionogels mPADs were compared with the non-treated mPADs.It was demonstrated that both ionogels highly affect the fluid flow by delaying the flow due to their different physical and chemical properties and water holding capacities.

Low cost 3D global instability analysis and flow sensitivity based on dynamic mode decomposition and high-order numerical tools

We explore the recently developed snapshot-based dynamic mode decomposition (DMD) technique, a matrix-free Arnoldi type method, to predict 3D linear global flow instabilities. We apply the DMD technique to flows confined in an L-shaped cavity and compare the resulting modes to their counterparts issued from classic, matrix forming, linear instability analysis (i.e. BiGlobal approach) and direct numerical simulations. Results show that the DMD technique, which uses snapshots generated by a 3D non-linear incompressible discontinuous Galerkin Navier?Stokes solver, provides very similar results to classical linear instability analysis techniques. In addition, we compare DMD results issued from non-linear and linearised Navier?Stokes solvers, showing that linearisation is not necessary (i.e. base flow not required) to obtain linear modes, as long as the analysis is restricted to the exponential growth regime, that is, flow regime governed by the linearised Navier?Stokes equations, and showing the potential of this type of analysis based on snapshots to general purpose CFD codes, without need of modifications. Finally, this work shows that the DMD technique can provide three-dimensional direct and adjoint modes through snapshots provided by the linearised and adjoint linearised Navier?Stokes equations advanced in time. Subsequently, these modes are used to provide structural sensitivity maps and sensitivity to base flow modification information for 3D flows and complex geometries, at an affordable computational cost. The information provided by the sensitivity study is used to modify the L-shaped geometry and control the most unstable 3D mode.

Temporal Control of Ion Channel Activation

Ion channels are a large class of integral membrane proteins that allow for the diffusion of ions across a cellular membrane and are found in all forms of life. Pentameric ligand-gated ion channels (pLGICs) comprise a large family of proteins that include the nicotinic acetylcholine receptor (nAChR) and the γ-aminobutyric acid (GABA) receptor. These ion channels are responsible for the fast synaptic transmission that occurs in humans and as a result are of fundamental biological importance. pLGICs bind ligands (neurotransmitters), and upon ligand-binding undergo activation. The activation event causes an ion channel to enter a new physical state that is able to conduct ions. Ion channels allow for the flux of ions across the membrane through a pore that is formed upon ion channel activation. For pLGICs to function properly both ligand-binding and ion channel activation must occur. The ligand-binding event has been studied extensively over the past few decades, and a detailed mechanism of binding has emerged. During activation the ion channel must undergo structural rearrangements that allow the protein to enter a conformation in which ions can flow through. Despite this great and ubiquitous importance, a fundamental understanding of the ion channel activation mechanism and kinetics, as well as concomitant structural arrangements, remains elusive.

This dissertation describes efforts that have been made to temporally control the activation of ligand-gated ion channels. Temporal control of ion channel activation provides a means by which to activate ion channels when desired. The majority of this work examines the use of light to activate ion channels. Several photocages were examined in this thesis; photocages are molecules that release a ligand under irradiation, and, for the work described here, the released ligand then activates the ion channel. First, a new water-soluble photoacid was developed for the activation of proton-sensitive ion channels. Activation of acid-sensing ion channels, ASIC2a and GLIC, was observed only upon irradiation. Next, a variety of Ru²⁺ photocages were also developed for the release of amine ligands. The Ru²⁺ systems interacted in a deleterious manner with a representative subset of biologically essential ion channels. The rapid mixing of ion channels with agonist was also examined. A detection system was built to monitor ion channels activation in the rapid mixing experiments. I have shown that liposomes, and functionally-reconstituted ELIC, are not destroyed during the mixing process. The work presented here provides the means to deliver agonist to ligand-gated ion channels in a controlled fashion.

Flow around a hemisphere-cylinder at high angle of attack and low Reynolds number. Part II: POD and DMD applied to reduced domains

Three-dimensional direct numerical simulations (DNS) have been performed on a finite-size hemispherecylinder model at angle of attack AoA = 20◦ and Reynolds numbers Re = 350 and 1000. Under these conditions, massive separation exists on the nose and lee-side of the cylinder, and at both Reynolds numbers the flow is found to be unsteady. Proper orthogonal decomposition (POD) and dynamic mode decomposition (DMD) are employed in order to study the primary instability that triggers unsteadiness at Re = 350. The dominant coherent flow structures identified at the lower Reynolds number are also found to exist at Re = 1000; the question is then posed whether the flow oscillations and structures found at the two Reynolds numbers are related. POD and DMD computations are performed using different subdomains of the DNS computational domain. Besides reducing the computational cost of the analyses, this also permits to isolate spatially localized oscillatory structures from other, more energetic structures present in the flow. It is found that POD and DMD are in general sensitive to domain truncation and noneducated choices of the subdomain may lead to inconsistent results. Analyses at Re = 350 show that the primary instability is related to the counter rotating vortex pair conforming the three-dimensional afterbody wake, and characterized by the frequency St ≈ 0.11, in line with results in the literature. At Re = 1000, vortex-shedding is present in the wake with an associated broadband spectrum centered around the same frequency. The horn/leeward vortices at the cylinder lee-side, upstream of the cylinder base, also present finite amplitude oscillations at the higher Reynolds number. The spatial structure of these oscillations, described by the POD modes, is easily differentiated from that of the wake oscillations. Additionally, the frequency spectra associated with the lee-side vortices presents well defined peaks, corresponding to St ≈ 0.11 and its few harmonics, as opposed to the broadband spectrum found at the wake.

Real-time load-side control of electric power systems

Two trends are emerging from modern electric power systems: the growth of renewable (e.g., solar and wind) generation, and the integration of information technologies and advanced power electronics. The former introduces large, rapid, and random fluctuations in power supply, demand, frequency, and voltage, which become a major challenge for real-time operation of power systems. The latter creates a tremendous number of controllable intelligent endpoints such as smart buildings and appliances, electric vehicles, energy storage devices, and power electronic devices that can sense, compute, communicate, and actuate. Most of these endpoints are distributed on the load side of power systems, in contrast to traditional control resources such as centralized bulk generators. This thesis focuses on controlling power systems in real time, using these load side resources. Specifically, it studies two problems.

(1) Distributed load-side frequency control: We establish a mathematical framework to design distributed frequency control algorithms for flexible electric loads. In this framework, we formulate a category of optimization problems, called optimal load control (OLC), to incorporate the goals of frequency control, such as balancing power supply and demand, restoring frequency to its nominal value, restoring inter-area power flows, etc., in a way that minimizes total disutility for the loads to participate in frequency control by deviating from their nominal power usage. By exploiting distributed algorithms to solve OLC and analyzing convergence of these algorithms, we design distributed load-side controllers and prove stability of closed-loop power systems governed by these controllers. This general framework is adapted and applied to different types of power systems described by different models, or to achieve different levels of control goals under different operation scenarios. We first consider a dynamically coherent power system which can be equivalently modeled with a single synchronous machine. We then extend our framework to a multi-machine power network, where we consider primary and secondary frequency controls, linear and nonlinear power flow models, and the interactions between generator dynamics and load control.

(2) Two-timescale voltage control: The voltage of a power distribution system must be maintained closely around its nominal value in real time, even in the presence of highly volatile power supply or demand. For this purpose, we jointly control two types of reactive power sources: a capacitor operating at a slow timescale, and a power electronic device, such as a smart inverter or a D-STATCOM, operating at a fast timescale. Their control actions are solved from optimal power flow problems at two timescales. Specifically, the slow-timescale problem is a chance-constrained optimization, which minimizes power loss and regulates the voltage at the current time instant while limiting the probability of future voltage violations due to stochastic changes in power supply or demand. This control framework forms the basis of an optimal sizing problem, which determines the installation capacities of the control devices by minimizing the sum of power loss and capital cost. We develop computationally efficient heuristics to solve the optimal sizing problem and implement real-time control. Numerical experiments show that the proposed sizing and control schemes significantly improve the reliability of voltage control with a moderate increase in cost.

Vortex Shedding Dynamics in Long Aspect-Ratio Aerodynamics Bodies

The focus of the current dissertation is to study qualitatively the underlying physics of vortex-shedding and wake dynamics in long aspect-ratio aerodynamics in incompressible viscous flow through the use of the KLE method. We carried out a long series of numerical experiments in the cases of flow around the cylinder at low Reynolds numbers. The study of flow at low Reynolds numbers provides an insight in the fluid physics and also plays a critical role when applying to stalled turbine rotors. Many of the conclusions about the qualitative nature of the physical mechanisms characterizing vortex formation, shedding and further interaction analyzed here at low Re could be extended to other Re regimes and help to understand the separation of the boundary layers in airfoils and other aerodynamic surfaces. In the long run, it aims to provide a better understanding of the complex multi-physics problems involving fluid-structure-control interaction through improved mathematical computational models of the multi-physics process. Besides the scientific conclusions produced, the research work on streamlined and bluff-body condition will also serve as a valuable guide for the future design of blade aerodynamics and the placement of wind turbines and hydrakinetic turbines, increasing the efficiency in the use of expensive workforce, supplies, and infrastructure. After the introductory section describing the main fields of application of wind power and hydrokinetic turbines, we describe the main features and theoretical background of the numerical method used here. Then, we present the analysis of the numerical experimentation results for the oscillatory regime right before the onset of vortex shedding for circular cylinders. We verified the wake length of the closed near-wake behind the cylinder and analysed the decay of the wake at the wake formation region, and then studied the St-Re relationship at the Reynolds numbers before the wake sheds compared to the experimental data. We found a theoretical model that describes the time evolution of the amplitude of fluctuations in the vorticity field on the twin vortex wake, which accurately matches the numerical results in terms of the frequency of the oscillation and rate of decay. We also proposed a model based on an analog circuit that is able to interpret the concerning flow by reducing the number of degrees of freedom. It follows the idea of the non-linear oscillator and resembles the dynamics mechanism of the closed near-wake with a common configured sine wave oscillator. This low-dimensional circuital model may also help to understand the underlying physical mechanisms, related to vorticity transport, that give origin to those oscillations.

Development and Verification of Control and Protection Strategies in Hybrid AC/DC Power Systems for Smart Grid Applications

Modern power networks incorporate communications and information technology infrastructure into the electrical power system to create a smart grid in terms of control and operation. The smart grid enables real-time communication and control between consumers and utility companies allowing suppliers to optimize energy usage based on price preference and system technical issues. The smart grid design aims to provide overall power system monitoring, create protection and control strategies to maintain system performance, stability and security. This dissertation contributed to the development of a unique and novel smart grid test-bed laboratory with integrated monitoring, protection and control systems. This test-bed was used as a platform to test the smart grid operational ideas developed here. The implementation of this system in the real-time software creates an environment for studying, implementing and verifying novel control and protection schemes developed in this dissertation. Phasor measurement techniques were developed using the available Data Acquisition (DAQ) devices in order to monitor all points in the power system in real time. This provides a practical view of system parameter changes, system abnormal conditions and its stability and security information system. These developments provide valuable measurements for technical power system operators in the energy control centers. Phasor Measurement technology is an excellent solution for improving system planning, operation and energy trading in addition to enabling advanced applications in Wide Area Monitoring, Protection and Control (WAMPAC). Moreover, a virtual protection system was developed and implemented in the smart grid laboratory with integrated functionality for wide area applications. Experiments and procedures were developed in the system in order to detect the system abnormal conditions and apply proper remedies to heal the system. A design for DC microgrid was developed to integrate it to the AC system with appropriate control capability. This system represents realistic hybrid AC/DC microgrids connectivity to the AC side to study the use of such architecture in system operation to help remedy system abnormal conditions. In addition, this dissertation explored the challenges and feasibility of the implementation of real-time system analysis features in order to monitor the system security and stability measures. These indices are measured experimentally during the operation of the developed hybrid AC/DC microgrids. Furthermore, a real-time optimal power flow system was implemented to optimally manage the power sharing between AC generators and DC side resources. A study relating to real-time energy management algorithm in hybrid microgrids was performed to evaluate the effects of using energy storage resources and their use in mitigating heavy load impacts on system stability and operational security.

Movement effects on the flow physics and nutrient delivery in engineered valvular tissues

Mechanical conditioning has been shown to promote tissue formation in a wide variety of tissue engineering efforts. However the underlying mechanisms by which external mechanical stimuli regulate cells and tissues are not known. This is particularly relevant in the area of heart valve tissue engineering (HVTE) owing to the intense hemodynamic environments that surround native valves. Some studies suggest that oscillatory shear stress (OSS) caused by steady flow and scaffold flexure play a critical role in engineered tissue formation derived from bone marrow derived stem cells (BMSCs). In addition, scaffold flexure may enhance nutrient (e.g. oxygen, glucose) transport. In this study, we computationally quantified the i) magnitude of fluid-induced shear stresses; ii) the extent of temporal fluid oscillations in the flow field using the oscillatory shear index (OSI) parameter, and iii) glucose and oxygen mass transport profiles. Noting that sample cyclic flexure induces a high degree of oscillatory shear stress (OSS), we incorporated moving boundary computational fluid dynamic simulations of samples housed within a bioreactor to consider the effects of: 1) no flow, no flexure (control group), 2) steady flow-alone, 3) cyclic flexure-alone and 4) combined steady flow and cyclic flexure environments. We also coupled a diffusion and convention mass transport equation to the simulated system. We found that the coexistence of both OSS and appreciable shear stress magnitudes, described by the newly introduced parameter OSI-:τ: explained the high levels of engineered collagen previously observed from combining cyclic flexure and steady flow states. On the other hand, each of these metrics on its own showed no association. This finding suggests that cyclic flexure and steady flow synergistically promote engineered heart valve tissue production via OSS, so long as the oscillations are accompanied by a critical magnitude of shear stress. In addition, our simulations showed that mass transport of glucose and oxygen is enhanced by sample movement at low sample porosities, but did not play a role in highly porous scaffolds. Preliminary in-house in vitro experiments showed that cell proliferation and phenotype is enhanced in OSI-:τ: environments.^

Vortex evolution in the wake of accelerating low-aspect-ratio plates and three-dimensional bodies of revolution

To find examples of effecient locomotion and manoeuvrability, one need only turn to the elegant solutions natural flyers and swimmers have converged upon. This dissertation is specifically motivated by processes of evolutionary convergence, which have led to the propulsors and body shapes in nature that exhibit strong geometric collapse over diverse scales. These body features are abstracted in the studies presented herein using low-aspect-ratio at plates and a three-dimensional body of revolution (a sphere). The highly-separated vortical wakes that develop during accelerations are systematically characterized as a function of planform shape, aspect ratio, Reynolds number, and initial boundary conditions. To this end, force measurements and time-resolved (planar) particle image velocimetry have been used throughout to quantify the instantaneous forces and vortex evolution in the wake of the bluff bodies. During rectilinear motions, the wake development for the flat plates is primarily dependent on plate aspect ratio, with edge discontinuities and curvature playing only a secondary role. Furthermore, the axisymmetric case, i.e. the circular plate, shows strong sensitivity to Reynolds number, while this sensitivity quickly diminishes with increasing aspect ratio. For rotational motions, global insensitivity to plate aspect ratio has been observed. For the sphere, it has been shown that accelerations play an important role in the mitigation of flow separation. These results - expounded upon in this dissertation - have begun to shed light on the specific vortex dynamics that may be coopted by flying and swimming species of all shapes and sizes towards efficient locomotion.

Development of a continuous process for α-thio-β-chloroacrylamide synthesis with enhanced control of a cascade transformation

A continuous process strategy has been developed for the preparation of α-thio-β chloroacrylamides, a class of highly versatile synthetic intermediates. Flow platforms to generate the α-chloroamide and α-thioamide precursors were successfully adopted, progressing from the previously employed batch chemistry, and in both instances afford a readily scalable methodology. The implementation of the key α-thio-β-chloroacrylamide casade as a continuous flow reaction on a multi-gram scale is described, while the tuneable nature of the cascade, facilitated by continuous processing, is highlighted by selective generation of established intermediates and byproducts.

Water Delivery Quality and Automatic Control Modes in Irrigation Canals. A Case Study.

Regarding canal management modernization, water savings and water delivery quality, the study presents two automatic canal control approaches of the PI (Proportional and Integral) type: the distant and the local downstream control modes. The two PI controllers are defined, tuned and tested using an hydraulic unsteady flow simulation model, particularly suitable for canal control studies. The PI control parameters are tuned using optimization tools. The simulations are done for a Portuguese prototype canal and the PI controllers are analyzed and compared considering a demand-oriented-canal operation. The paper presents and analyzes the two control modes answers for five different offtake types – gate controlled weir, gate controlled orifice, weir with or without adjustable height and automatic flow adjustable offtake. The simulation results are compared using water volumes performance indicators (considering the demanded, supplied and the effectives water volumes) and a time indicator, defined taking into account the time during which the demand discharges are effective discharges. Regarding water savings, the simulation results for the five offtake types prove that the local downstream control gives the best results (no water operational losses) and that the distant downstream control presents worse results in connection with the automatic flow adjustable offtakes. Considering the water volumes and time performance indicators, the best results are obtained for the automatic flow adjustable offtakes and the worse for the gate controlled orifices, followed by the weir with adjustable height.

Special Issue “ Fluid Flow, Energy Transfer and Design”

Fluids are important because of their preponderance in our lives. Fluid mechanics touches almost every aspect of our daily lives, and it plays a central role in many branches of science and technology. Therefore, it is a challenging and exciting field of scientific activity due to the complexity of the subject studied and the breadth of the applications. The quest for advances in fluid mechanics, as in other scientific fields, emerge from analytical, computational (CFD) and experimental studies. The improvement in our ability to describe, predict and control the phenomena played (and plays) key roles in the technological breakthroughs. The present theme issue of “Fluid and Heat Flow: Simulation and Optimization” collects a selection of papers. selection of papers presented at Special Session “Fluid Flow, Energy Transfer and Design”