Instability Analysis of Incompressible Open Cavity Flows

El problema del flujo sobre una cavidad abierta ha sido estudiado en profundidad en la literatura, tanto por el interés académico del problema como por sus aplicaciones prácticas en gran variedad de problemas ingenieriles, como puede ser el alojamiento del tren de aterrizaje de aeronaves, o el depósito de agua de aviones contraincendios. Desde hace muchos a˜nos se estudian los distintos tipos de inestabilidades asociadas a este problema: los modos bidimensionales en la capa de cortadura, y los modos tridimensionales en el torbellino de recirculación principal dentro de la cavidad. En esta tesis se presenta un estudio paramétrico completo del límite incompresible del problema, empleando la herramienta de estabilidad lineal conocida como BiGlobal. Esta aproximación permite contemplar la estabilidad global del flujo, y obtener tanto la forma como las características de los modos propios del problema físico, sean estables o inestables. El estudio realizado permite caracterizar con gran detalle todos los modos relevantes, así como la envolvente de estabilidad en el espacio paramétrico del problema incompresible (Mach nulo, variación de Reynolds, espesor de capa límite incidente, relación altura/profundidad de la cavidad, y longitud característica de la perturbación en la dirección transversal). A la luz de los resultados obtenidos se proponen una serie de relaciones entre los parámetros y características de los modos principales, como por ejemplo entre el Reynolds crítico de un modo, y la longitud característica del mismo. Los resultados numéricos se contrastan con una campaña experimental, siendo la principal conclusión de dicha comparación que los modos lineales están presentes en el flujo real saturado, pero que existen diferencias notables en frecuencia entre las predicciones teóricas y los experimentos. Para intentar determinar la naturaleza de dichas diferencias se realiza una simulación numérica directa tridimensional, y se utiliza un algoritmo de DMD (descomposición dinámica de modos) para describir el proceso de saturación. ABSTRACT The problem of the flow over an open cavity has been studied in depth in the literature, both for being an interesting academical problem and due to the multitude of industrial applications, like the landing gear of aircraft, or the water deposit of firefighter airplanes. The different types of instabilities appearing in this flow studied in the literature are two: the two-dimensional shear layer modes, and the three-dimensional modes that appear in the main recirculating vortex inside the cavity. In this thesis a parametric study in the incompressible limit of the problem is presented, using the linear stability analysis known as BiGlobal. This approximation allows to obtain the global stability behaviour of the flow, and to capture both the morphological features and the characteristics of the eigenmodes of the physical problem, whether they are stable or unstable. The study presented here characterizes with great detail all the relevant eigenmodes, as well as the hypersurface of instability on the parameter space of the incompressible problem (Mach equal to zero, and variation of the Reynolds number, the incoming boundary layer thickness, the length to depth aspect ratio of the cavity and the spanwise length of the perturbation). The results allow to construct parametric relations between the characteristics of the leading eigenmodes and the parameters of the problem, like for example the one existing between the critical Reynolds number and its characteristic length. The numerical results presented here are compared with those of an experimental campaign, with the main conclusion of said comparison being that the linear eigenmode are present in the real saturated flow, albeit with some significant differences in the frequencies of the experiments and those predicted by the theory. To try to determine the nature of those differences a three-dimensional direct numerical simulation, analyzed with Dynamic Mode Decomposition algorithm, was used to describe the process of saturation.

Linear instability analysis of incompressible flow over a cuboid cavity

Direct numerical simulations are performed to analyze the three-dimensional instability of flows over three-dimensional cavities. The flow structures at different Reynolds numbers are investigated by using the spectral-element solver nek5000. As the Reynolds number increasing, the lateral wall effects become more important, the recirculation zone shrinks, the front vortex increases and the flow structure inside of the cavity becomes more complex. Results show that the flow bifurcates from a steady state to an oscillatory regime beyond a value of Reynolds number Re = 1100.

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.

GroEL-GroES-mediated protein folding requires an intact central cavity

The chaperonin GroEL is an oligomeric double ring structure that, together with the cochaperonin GroES, assists protein folding. Biochemical analyses indicate that folding occurs in a cis ternary complex in which substrate is sequestered within the GroEL central cavity underneath GroES. Recently, however, studies of GroEL “minichaperones” containing only the apical substrate binding subdomain have questioned the functional importance of substrate encapsulation within GroEL-GroES complexes. Minichaperones were reported to assist folding despite the fact that they are monomeric and therefore cannot form a central cavity. Here we compare directly the folding activity of minichaperones with that of the full GroEL-GroES system. In agreement with earlier studies, minichaperones assist folding of some proteins. However, this effect is observed only under conditions where substantial spontaneous folding is also observed and is indistinguishable from that resulting from addition of the nonchaperone protein α-casein. By contrast, the full GroE system efficiently promotes folding of several substrates under conditions where essentially no spontaneous folding is observed. These data argue that the full GroEL folding activity requires the intact GroEL-GroES complex, and in light of previous studies, underscore the importance of substrate encapsulation for providing a folding environment distinct from the bulk solution.

Evidence that a prominent cavity in the coiled coil of HIV type 1 gp41 is an attractive drug target

Synthetic C peptides, corresponding to the C helix of the HIV type 1 (HIV-1) gp41 envelope protein, are potent inhibitors of HIV-1 membrane fusion. One such peptide is in clinical trials. The crystal structure of the gp41 core, in its proposed fusion-active conformation, is a trimer of helical hairpins in which three C helices pack against a central coiled coil. Each C helix shows especially prominent contacts with one of three symmetry-related, hydrophobic cavities on the surface of the coiled coil. We show that the inhibitory activity of the C peptide C34 depends on its ability to bind to this coiled-coil cavity. Moreover, examining a series of C34 peptide variants with modified cavity-binding residues, we find a linear relationship between the logarithm of the inhibitory potency and the stability of the corresponding helical-hairpin complexes. Our results provide strong evidence that this coiled-coil cavity is a good drug target and clarify the mechanism of C peptide inhibition. They also suggest simple, quantitative assays for the identification and evaluation of analogous inhibitors of HIV-1 entry.

A GIS-based methodology to quantitatively define an Adjacent Protected Area in a shallow karst cavity: The case of Altamira cave

Different types of land use are usually present in the areas adjacent to many shallow karst cavities. Over time, the increasing amount of potentially harmful matter and energy, of mainly anthropic origin or influence, that reaches the interior of a shallow karst cavity can modify the hypogeal ecosystem and increase the risk of damage to the Palaeolithic rock art often preserved within the cavity. This study proposes a new Protected Area status based on the geological processes that control these matter and energy fluxes into the Altamira cave karst system. Analysis of the geological characteristics of the shallow karst system shows that direct and lateral infiltration, internal water circulation, ventilation, gas exchange and transmission of vibrations are the processes that control these matter and energy fluxes into the cave. This study applies a comprehensive methodological approach based on Geographic Information Systems (GIS) to establish the area of influence of each transfer process. The stratigraphic and structural characteristics of the interior of the cave were determined using 3D Laser Scanning topography combined with classical field work, data gathering, cartography and a porosity–permeability analysis of host rock samples. As a result, it was possible to determine the hydrogeological behavior of the cave. In addition, by mapping and modeling the surface parameters it was possible to identify the main features restricting hydrological behavior and hence direct and lateral infiltration into the cave. These surface parameters included the shape of the drainage network and a geomorphological and structural characterization via digital terrain models. Geological and geomorphological maps and models integrated into the GIS environment defined the areas involved in gas exchange and ventilation processes. Likewise, areas that could potentially transmit vibrations directly into the cave were identified. This study shows that it is possible to define a Protected Area by quantifying the area of influence related to each transfer process. The combined maximum area of influence of all the processes will result in the new Protected Area. This area will thus encompass all the processes that account for most of the matter and energy carried into the cave and will fulfill the criteria used to define the Protected Area. This methodology is based on the spatial quantification of processes and entities of geological origin and can therefore be applied to any shallow karst system that requires protection.

Tight-Binding Approximations in 1D and 2D Coupled-Cavity Photonic Crystal Structures

Light confinement and controlling an optical field has numerous applications in the field of telecommunications for optical signals processing. When the wavelength of the electromagnetic field is on the order of the period of a photonic microstructure, the field undergoes reflection, refraction, and coherent scattering. This produces photonic bandgaps, forbidden frequency regions or spectral stop bands where light cannot exist. Dielectric perturbations that break the perfect periodicity of these structures produce what is analogous to an impurity state in the bandgap of a semiconductor. The defect modes that exist at discrete frequencies within the photonic bandgap are spatially localized about the cavity-defects in the photonic crystal. In this thesis the properties of two tight-binding approximations (TBAs) are investigated in one-dimensional and two-dimensional coupled-cavity photonic crystal structures We require an efficient and simple approach that ensures the continuity of the electromagnetic field across dielectric interfaces in complex structures. In this thesis we develop \textrm{E} -- and \textrm{D} --TBAs to calculate the modes in finite 1D and 2D two-defect coupled-cavity photonic crystal structures. In the \textrm{E} -- and \textrm{D} --TBAs we expand the coupled-cavity \overrightarrow{E} --modes in terms of the individual \overrightarrow{E} -- and \overrightarrow{D} --modes, respectively. We investigate the dependence of the defect modes, their frequencies and quality factors on the relative placement of the defects in the photonic crystal structures. We then elucidate the differences between the two TBA formulations, and describe the conditions under which these formulations may be more robust when encountering a dielectric perturbation. Our 1D analysis showed that the 1D modes were sensitive to the structure geometry. The antisymmetric \textrm{D} mode amplitudes show that the \textrm{D} --TBA did not capture the correct (tangential \overrightarrow{E} --field) boundary conditions. However, the \textrm{D} --TBA did not yield significantly poorer results compared to the \textrm{E} --TBA. Our 2D analysis reveals that the \textrm{E} -- and \textrm{D} --TBAs produced nearly identical mode profiles for every structure. Plots of the relative difference between the \textrm{E} and \textrm{D} mode amplitudes show that the \textrm{D} --TBA did capture the correct (normal \overrightarrow{E} --field) boundary conditions. We found that the 2D TBA CC mode calculations were 125-150 times faster than an FDTD calculation for the same two-defect PCS. Notwithstanding this efficiency, the appropriateness of either TBA was found to depend on the geometry of the structure and the mode(s), i.e. whether or not the mode has a large normal or tangential component.

DESIGN OF A MINIATURIZED SOP INTEGRATED GNSS RF FRONT-END MODULE

Navigation devices used to be bulky and expensive and were not widely commercialized for personal use. Nowadays, all useful electronic devices are turning into being handheld so that they can be conveniently used anytime and anywhere. One can claim that almost any mobile phone, used today, has quite strong navigational capabilities that can efficiently work anywhere in the globe. No matter where you are, you can easily know your exact location and make your way smoothly to wherever you would like to go. This couldn’t have been made possible without the existence of efficient and small microwave circuits responsible for the transmission and reception of high quality navigation signals. This thesis is mainly concerned with the design of novel highly miniaturized and efficient filtering components working in the Global Navigational Satellite Systems (GNSS) frequency band to be integrated within an efficient Radio Frequency (RF) front-end module (FEM). A System-on-Package (SoP) integration technique is adopted for the design of all the components in this thesis. Two novel miniaturized filters are designed, where one of them is a wideband filter targeting the complete GNSS band with a fractional bandwidth of almost 50% at a center frequency of 1.385 GHz. This filter utilizes a direct inductive coupling topology to achieve the required wide band performance. It also has very good out-of-band rejection and low IL. Whereas the other dual band filter will only cover the lower and upper GNSS bands with a rejection notch in between the two bands. It has very good inter band rejection. The well-known “divide and conquer” design methodology was applied for the design of this filter to help save valuable design and optimization time. Moreover, the performance of two commercially available ultra-Low Noise Amplifiers (LNAs) is studied. The complete RF FEM showed promising preliminary performance in terms of noise figure, gain and bandwidth, where it out performed other commercial front-ends in these three aspects. All the designed circuits are fabricated and tested. The measured results are found to be in good agreements with the simulations.

Genetically engineered silk and genipin-enhanced fibrin hydrogel for annulus fibrosus repair

INTRODUCTION: Around 80% of people are affected by low back pain at least once in their life, often caused by trauma provoking intervertebral disc (IVD) herniation and/or IVD degeneration. Apart from some promising approaches for nucleus pulposus repair, so far no treatment or repair is available for the outer fibrous tissue, annulus fibrosus (AF). We aimed for sealing and repairing an AF injury in a bovine IVD organ culture model in vitro over 14 days under different loading conditions. For this purpose, a silk fleece composite from Bombyx mori silk was combined with genipin-enhanced fibrin hydrogel [1]. METHODS: Bovine IVDs of 12-17 months old animals were isolated by first removing all surrounding tissue, followed by cutting out the IVDs [2]. Culturing of discs occurred in high glucose Dulbecco's Modified Eagle Medium (HG-DMEM) supplemented with 5% serum as previously described. On the next day, injury was induced using a 2mm biopsy punch (Polymed, Switzerland). The formed cavity was filled with (0.4%) genipin-enhanced human based fibrin hydrogel (35- 55mg/mL human fibrinogen, Baxter, Austria) and sealed with a silk fleece-membrane composite (Spintec Engineering, Germany). Different culture conditions were applied: free swelling, static diurnal load of 0.2MPa for 8h/d and complex loading at 0.2MPa compression combined with ± 2° torsion at 0.2Hz for 8h/d. Complex loading was applied by a custom built 2 degree of freedom bioreactor [3]. After 14 days of culture cell activity was determined with resazurin assay. Additionally, glycosaminoglycan (dimethyl-methylene blue), DNA (Hoechst) and collagen content (hydroxy-proline) were determined. Finally, real-time qPCR of major IVD marker genes was performed. RESULTS: The silk seal closing the injury site could successfully withstand the forces of all three loading conditions with no misplacement over the two weeks’ culture. Nevertheless, disc height of the repaired discs did not significantly differ from the injured group. The disc phenotype could be maintained as demonstrated by biochemical analysis of gene expression, cell activity, DNA-, collagen- and GAG content. The silk itself was evaluated to be highly biocompatible for hMSC, as revealed by cytotoxicity assays. DISCUSSION & CONCLUSIONS: The silk can be considered a highly-elastic and biocompatible material for AF closure and the genipin-enhanced fibrin hydrogel has also good biomechanical properties. However, the cyto-compatibility of genipin seems rather poor and other hydrogels and/or cross-linkers should be looked into. REFERENCES: 1 C.C. Guterl et al. (2014) Characterization of Mechanics and Cytocompatibility of Fibrin Genipin Annulus Fibrosus Sealant with the Addition of Cell Adhesion Molecules, Tissue Eng Part A 2 S.C. Chan, B. Gantenbein-Ritter (2012) Preparation of intact bovine tail intervertebral discs for organ culture, J Vis Exp 3 B Gantenbein et al. (2015) Organ Culture Bioreactors - Platforms to Study Human Intervertebral Disc Degeneration and Regenerative Therapy, Curr Stem Cell Res Ther [epub ahead of print] ACKNOWLEDGEMENTS: This project is supported by the Gebert Rüf Stiftung project # GRS-028/13.

Maintenance of coelomic fluid pH in sea urchins exposed to elevated CO2: the role of body cavity epithelia and stereom dissolution

Experimental ocean acidification leads to a shift in resource allocation and to an increased [HCO3-] within the perivisceral coelomic fluid (PCF) in the Baltic green sea urchin Strongylocentrotus droebachiensis. We investigated putative mechanisms of this pH compensation reaction by evaluating epithelial barrier function and the magnitude of skeleton (stereom) dissolution. In addition, we measured ossicle growth and skeletal stability. Ussing chamber measurements revealed that the intestine formed a barrier for HCO3- and was selective for cation diffusion. In contrast, the peritoneal epithelium was leaky and only formed a barrier for macromolecules. The ossicles of 6 week high CO2-acclimatised sea urchins revealed minor carbonate dissolution, reduced growth but unchanged stability. On the other hand, spines dissolved more severely and were more fragile following acclimatisation to high CO2. Our results indicate that epithelia lining the PCF space contribute to its acid-base regulation. The intestine prevents HCO3- diffusion and thus buffer leakage. In contrast, the leaky peritoneal epithelium allows buffer generation via carbonate dissolution from the surrounding skeletal ossicles. Long-term extracellular acid-base balance must be mediated by active processes, as sea urchins can maintain relatively high extracellular [HCO3-]. The intestinal epithelia are good candidate tissues for this active net import of HCO3- into the PCF. Spines appear to be more vulnerable to ocean acidification which might significantly impact resistance to predation pressure and thus influence fitness of this keystone species.

Shell size, body mass and average pO2 in mantle cavity water of different marine mollusc species

Marine invertebrates with open circulatory system establish low and constant oxygen partial pressure (Po2) around their tissues. We hypothesized that as a first step towards maintenance of low haemolymph and tissue oxygenation, the Po2 in molluscan mantle cavity water should be lowered against normoxic (21 kPa) seawater Po2, but balanced high enough to meet the energetic requirements in a given species. We recorded Po2 in mantle cavity water of five molluscan species with different lifestyles, two pectinids (Aequipecten opercularis, Pecten maximus), two mud clams (Arctica islandica, Mya arenaria), and a limpet (Patella vulgata). All species maintain mantle cavity water oxygenation below normoxic Po2. Average mantle cavity water Po2 correlates positively with standard metabolic rate (SMR): highest in scallops and lowest in mud clams. Scallops show typical Po2 frequency distribution, with peaks between 3 and 10 kPa, whereas mud clams and limpets maintain mantle water Po2 mostly <5 kPa. Only A. islandica and P. vulgata display distinguishable temporal patterns in Po2 time series. Adjustment of mantle cavity Po2 to lower than ambient levels through controlled pumping prevents high oxygen gradients between bivalve tissues and surrounding fluid, limiting oxygen flux across the body surface. The patterns of Po2 in mantle cavity water correspond to molluscan ecotypes.

Greenland and Antarctic ice sheet topography, cavity geometry, and global bathymetry (RTopo-2), links to NetCDF files

The ocean plays an important role in modulating the mass balance of the polar ice sheets by interacting with the ice shelves in Antarctica and with the marine-terminating outlet glaciers in Greenland. Given that the flux of warm water onto the continental shelf and into the sub-ice cavities is steered by complex bathymetry, a detailed topography data set is an essential ingredient for models that address ice-ocean interaction. We followed the spirit of the global RTopo-1 data set and compiled consistent maps of global ocean bathymetry, upper and lower ice surface topographies and global surface height on a spherical grid with now 30-arc seconds resolution. We used the General Bathymetric Chart of the Oceans (GEBCO, 2014) as the backbone and added the International Bathymetric Chart of the Arctic Ocean version 3 (IBCAOv3) and the Interna- tional Bathymetric Chart of the Southern Ocean (IBCSO) version 1. While RTopo-1 primarily aimed at a good and consistent representation of the Antarctic ice sheet, ice shelves and sub-ice cavities, RTopo-2 now also contains ice topographies of the Greenland ice sheet and outlet glaciers. In particular, we aimed at a good representation of the fjord and shelf bathymetry sur- rounding the Greenland continent. We corrected data from earlier gridded products in the areas of Petermann Glacier, Hagen Bræ and Sermilik Fjord assuming that sub-ice and fjord bathymetries roughly follow plausible Last Glacial Maximum ice flow patterns. For the continental shelf off northeast Greenland and the floating ice tongue of Nioghalvfjerdsfjorden Glacier at about 79°N, we incorporated a high-resolution digital bathymetry model considering original multibeam survey data for the region. Radar data for surface topographies of the floating ice tongues of Nioghalvfjerdsfjorden Glacier and Zachariæ Isstrøm have been obtained from the data centers of Technical University of Denmark (DTU), Operation Icebridge (NASA/NSF) and Alfred Wegener Institute (AWI). For the Antarctic ice sheet/ice shelves, RTopo-2 largely relies on the Bedmap-2 product but applies corrections for the geometry of Getz, Abbot and Fimbul ice shelf cavities.