Three-dimensional aspects of fluid flows in channels. II. Effects of meniscus and thin film regimes on viscous fingers

We perform a three-dimensional study of steady state viscous fingers that develop in linear channels. By means of a three-dimensional lattice-Boltzmann scheme that mimics the full macroscopic equations of motion of the fluid momentum and order parameter, we study the effect of the thickness of the channel in two cases. First, for total displacement of the fluids in the channel thickness direction, we find that the steady state finger is effectively two-dimensional and that previous two-dimensional results can be recovered by taking into account the effect of a curved meniscus across the channel thickness as a contribution to surface stresses. Second, when a thin film develops in the channel thickness direction, the finger narrows with increasing channel aspect ratio in agreement with experimental results. The effect of the thin film renders the problem three-dimensional and results deviate from the two-dimensional prediction.

Three-dimensional aspects of fluid flows in channels. II. Effects of meniscus and thin film regimes on viscous fingers

We perform a three-dimensional study of steady state viscous fingers that develop in linear channels. By means of a three-dimensional lattice-Boltzmann scheme that mimics the full macroscopic equations of motion of the fluid momentum and order parameter, we study the effect of the thickness of the channel in two cases. First, for total displacement of the fluids in the channel thickness direction, we find that the steady state finger is effectively two-dimensional and that previous two-dimensional results can be recovered by taking into account the effect of a curved meniscus across the channel thickness as a contribution to surface stresses. Second, when a thin film develops in the channel thickness direction, the finger narrows with increasing channel aspect ratio in agreement with experimental results. The effect of the thin film renders the problem three-dimensional and results deviate from the two-dimensional prediction.

Approximation and support theorem for a wave equation in two space dimensions

We prove a characterization of the support of the law of the solution for a stochastic wave equation with two-dimensional space variable, driven by a noise white in time and correlated in space. The result is a consequence of an approximation theorem, in the convergence of probability, for equations obtained by smoothing the random noise. For some particular classes of coefficients, approximation in the Lp-norm for p¿1 is also proved.

Electrophoretic mobility of supercoiled, catenated and knotted DNA molecules.

We systematically varied conditions of two-dimensional (2D) agarose gel electrophoresis to optimize separation of DNA topoisomers that differ either by the extent of knotting, the extent of catenation or the extent of supercoiling. To this aim we compared electrophoretic behavior of three different families of DNA topoisomers: (i) supercoiled DNA molecules, where supercoiling covered the range extending from covalently closed relaxed up to naturally supercoiled DNA molecules; (ii) postreplicative catenanes with catenation number increasing from 1 to ∼15, where both catenated rings were nicked; (iii) knotted but nicked DNA molecules with a naturally arising spectrum of knots. For better comparison, we studied topoisomer families where each member had the same total molecular mass. For knotted and supercoiled molecules, we analyzed dimeric plasmids whereas catenanes were composed of monomeric forms of the same plasmid. We observed that catenated, knotted and supercoiled families of topoisomers showed different reactions to changes of agarose concentration and voltage during electrophoresis. These differences permitted us to optimize conditions for their separation and shed light on physical characteristics of these different types of DNA topoisomers during electrophoresis.

Stochastic inversion of tracer test and electrical geophysical data to estimate hydraulic conductivities.

Quantifying the spatial configuration of hydraulic conductivity (K) in heterogeneous geological environments is essential for accurate predictions of contaminant transport, but is difficult because of the inherent limitations in resolution and coverage associated with traditional hydrological measurements. To address this issue, we consider crosshole and surface-based electrical resistivity geophysical measurements, collected in time during a saline tracer experiment. We use a Bayesian Markov-chain-Monte-Carlo (McMC) methodology to jointly invert the dynamic resistivity data, together with borehole tracer concentration data, to generate multiple posterior realizations of K that are consistent with all available information. We do this within a coupled inversion framework, whereby the geophysical and hydrological forward models are linked through an uncertain relationship between electrical resistivity and concentration. To minimize computational expense, a facies-based subsurface parameterization is developed. The Bayesian-McMC methodology allows us to explore the potential benefits of including the geophysical data into the inverse problem by examining their effect on our ability to identify fast flowpaths in the subsurface, and their impact on hydrological prediction uncertainty. Using a complex, geostatistically generated, two-dimensional numerical example representative of a fluvial environment, we demonstrate that flow model calibration is improved and prediction error is decreased when the electrical resistivity data are included. The worth of the geophysical data is found to be greatest for long spatial correlation lengths of subsurface heterogeneity with respect to wellbore separation, where flow and transport are largely controlled by highly connected flowpaths.

Dilaton quantum cosmology in two dimensions

We consider a renormalizable two-dimensional model of dilaton gravity coupled to a set of conformal fields as a toy model for quantum cosmology. We discuss the cosmological solutions of the model and study the effect of including the back reaction due to quantum corrections. As a result, when the matter density is below some threshold new singularities form in a weak-coupling region, which suggests that they will not be removed in the full quantum theory. We also solve the Wheeler-DeWitt equation. Depending on the quantum state of the Universe, the singularities may appear in a quantum region where the wave function is not oscillatory, i.e., when there is not a well-defined notion of classical spacetime.

Supergravity models of (3+1)-dimensional QCD

The most general black M5-brane solution of eleven-dimensional supergravity (with a flat R4 spacetime in the brane and a regular horizon) is characterized by charge, mass and two angular momenta. We use this metric to construct general dual models of large-N QCD (at strong coupling) that depend on two free parameters. The mass spectrum of scalar particles is determined analytically (in the WKB approximation) and numerically in the whole two-dimensional parameter space. We compare the mass spectrum with analogous results from lattice calculations, and find that the supergravity predictions are close to the lattice results everywhere on the two dimensional parameter space except along a special line. We also examine the mass spectrum of the supergravity Kaluza-Klein (KK) modes and find that the KK modes along the compact D-brane coordinate decouple from the spectrum for large angular momenta. There are however KK modes charged under a U(1)×U(1) global symmetry which do not decouple anywhere on the parameter space. General formulas for the string tension and action are also given.

Energy-momentum tensor for scalar fields coupled to the dilaton in two dimensions

We clarify some issues related to the evaluation of the mean value of the energy-momentum tensor for quantum scalar fields coupled to the dilaton field in two-dimensional gravity. Because of this coupling, the energy-momentum tensor for matter is not conserved and therefore it is not determined by the trace anomaly. We discuss different approximations for the calculation of the energy-momentum tensor and show how to obtain the correct amount of Hawking radiation. We also compute cosmological particle creation and quantum corrections to the Newtonian potential.

Single acquisition electrical property mapping based on relative coil sensitivities: A proof-of-concept demonstration.

PURPOSE: All methods presented to date to map both conductivity and permittivity rely on multiple acquisitions to compute quantitatively the magnitude of radiofrequency transmit fields, B1+. In this work, we propose a method to compute both conductivity and permittivity based solely on relative receive coil sensitivities ( B1-) that can be obtained in one single measurement without the need to neither explicitly perform transmit/receive phase separation nor make assumptions regarding those phases. THEORY AND METHODS: To demonstrate the validity and the noise sensitivity of our method we used electromagnetic finite differences simulations of a 16-channel transceiver array. To experimentally validate our methodology at 7 Tesla, multi compartment phantom data was acquired using a standard 32-channel receive coil system and two-dimensional (2D) and 3D gradient echo acquisition. The reconstructed electric properties were correlated to those measured using dielectric probes. RESULTS: The method was demonstrated both in simulations and in phantom data with correlations to both the modeled and bench measurements being close to identity. The noise properties were modeled and understood. CONCLUSION: The proposed methodology allows to quantitatively determine the electrical properties of a sample using any MR contrast, with the only constraint being the need to have 4 or more receive coils and high SNR. Magn Reson Med, 2014. © 2014 Wiley Periodicals, Inc.

Separation of ammonia and carbon dioxide by adsorption

Tässä työssä on tutkittu ammoniakin ja hiilidioksidin erottamista adsorptio prosessilla ja suunniteltiin paineen muunteluun perustuvan adsorptioprosessin (PSA) käyttöä. Työn tarkoituksena oli laskea adsorptioon perustuvan prosessin kannattavuus melamiinitehtaan poistokaasujen erotuksessa. Tätä varten työssä suunniteltiin tehdasmitta-kaavainen prosessi ja arvioitiin sen kannattavuus. Työssä mitattiin adsorptiotasapainot, joiden perusteella sovitettiin sopiva kokeellinen adsorptioisotermi. Adsorptioisotermi lisättiin simulointiohjelmaan, jonka avulla suunniteltiin kaksi vaihtoehtoista pilot laitteistoa kaasujen erottamiseksi. Toisella pilot laitteistolla saadaan mitattua vain läpäisykäyrät, mutta paremmalla versiolla saadaan myös tietoa erotettujen komponenttien puhtaudesta. Suunnittelun tärkeimpiä lähtökohtia on molempien komponenttien mahdollisimman korkea puhtaus ja talteenottoaste. Täysimittakaavainen tehdas suunniteltiin simulointiohjelmiston avulla kahdelle eri kapasiteetille ja arvioitiin niiden kustannukset ja kannattavuus. Adsorptioprosessit osoittautuivat kannattaviksi kaasuseoksen erottamisessa kummassakin tapauksessa

Kalman Filtering in Online Paper Quality Assessment

Controlling the quality variables (such as basis weight, moisture etc.) is a vital part of making top quality paper or board. In this thesis, an advanced data assimilation tool is applied to the quality control system (QCS) of a paper or board machine. The functionality of the QCS is based on quality observations that are measured with a traversing scanner making a zigzag path. The basic idea is the following: The measured quality variable has to be separated into its machine direction (MD) and cross direction (CD) variations due to the fact that the QCS works separately in MD and CD. Traditionally this is done simply by assuming one scan of the zigzag path to be the CD profile and its mean value to be one point of the MD trend. In this thesis, a more advanced method is introduced. The fundamental idea is to use the signals’ frequency components to represent the variation in both CD and MD. To be able to get to the frequency domain, the Fourier transform is utilized. The frequency domain, that is, the Fourier components are then used as a state vector in a Kalman filter. The Kalman filter is a widely used data assimilation tool to combine noisy observations with a model. The observations here refer to the quality measurements and the model to the Fourier frequency components. By implementing the two dimensional Fourier transform into the Kalman filter, we get an advanced tool for the separation of CD and MD components in total variation or, to be more general, for data assimilation. A piece of a paper roll is analyzed and this tool is applied to model the dataset. As a result, it is clear that the Kalman filter algorithm is able to reconstruct the main features of the dataset from a zigzag path. Although the results are made with a very short sample of paper roll, it seems that this method has great potential to be used later on as a part of the quality control system.

Critical behavior of a water monolayer under hydrophobic confinement

The properties of water can have a strong dependence on the confinement. Here, we consider a water monolayer nanoconfined between hydrophobic parallel walls under conditions that prevent its crystallization. We investigate, by simulations of a many-body coarse-grained water model, how the properties of the liquid are affected by the confinement. We show, by studying the response functions and the correlation length and by performing finite-size scaling of the appropriate order parameter, that at low temperature the monolayer undergoes a liquid-liquid phase transition ending in a critical point in the universality class of the two-dimensional (2D) Ising model. Surprisingly, by reducing the linear size L of the walls, keeping the walls separation h constant, we find a 2D-3D crossover for the universality class of the liquid-liquid critical point for L/h=~50, i.e. for a monolayer thickness that is small compared to its extension. This result is drastically different from what is reported for simple liquids, where the crossover occurs for , and is consistent with experimental results and atomistic simulations. We shed light on these findings showing that they are a consequence of the strong cooperativity and the low coordination number of the hydrogen bond network that characterizes water.

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.

Large-eddy simulation of wind flows over complex terrains for wind energy applications

Wind energy has obtained outstanding expectations due to risks of global warming and nuclear energy production plant accidents. Nowadays, wind farms are often constructed in areas of complex terrain. A potential wind farm location must have the site thoroughly surveyed and the wind climatology analyzed before installing any hardware. Therefore, modeling of Atmospheric Boundary Layer (ABL) flows over complex terrains containing, e.g. hills, forest, and lakes is of great interest in wind energy applications, as it can help in locating and optimizing the wind farms. Numerical modeling of wind flows using Computational Fluid Dynamics (CFD) has become a popular technique during the last few decades. Due to the inherent flow variability and large-scale unsteadiness typical in ABL flows in general and especially over complex terrains, the flow can be difficult to be predicted accurately enough by using the Reynolds-Averaged Navier-Stokes equations (RANS). Large- Eddy Simulation (LES) resolves the largest and thus most important turbulent eddies and models only the small-scale motions which are more universal than the large eddies and thus easier to model. Therefore, LES is expected to be more suitable for this kind of simulations although it is computationally more expensive than the RANS approach. With the fast development of computers and open-source CFD software during the recent years, the application of LES toward atmospheric flow is becoming increasingly common nowadays. The aim of the work is to simulate atmospheric flows over realistic and complex terrains by means of LES. Evaluation of potential in-land wind park locations will be the main application for these simulations. Development of the LES methodology to simulate the atmospheric flows over realistic terrains is reported in the thesis. The work also aims at validating the LES methodology at a real scale. In the thesis, LES are carried out for flow problems ranging from basic channel flows to real atmospheric flows over one of the most recent real-life complex terrain problems, the Bolund hill. All the simulations reported in the thesis are carried out using a new OpenFOAM® -based LES solver. The solver uses the 4th order time-accurate Runge-Kutta scheme and a fractional step method. Moreover, development of the LES methodology includes special attention to two boundary conditions: the upstream (inflow) and wall boundary conditions. The upstream boundary condition is generated by using the so-called recycling technique, in which the instantaneous flow properties are sampled on aplane downstream of the inlet and mapped back to the inlet at each time step. This technique develops the upstream boundary-layer flow together with the inflow turbulence without using any precursor simulation and thus within a single computational domain. The roughness of the terrain surface is modeled by implementing a new wall function into OpenFOAM® during the thesis work. Both, the recycling method and the newly implemented wall function, are validated for the channel flows at relatively high Reynolds number before applying them to the atmospheric flow applications. After validating the LES model over simple flows, the simulations are carried out for atmospheric boundary-layer flows over two types of hills: first, two-dimensional wind-tunnel hill profiles and second, the Bolund hill located in Roskilde Fjord, Denmark. For the twodimensional wind-tunnel hills, the study focuses on the overall flow behavior as a function of the hill slope. Moreover, the simulations are repeated using another wall function suitable for smooth surfaces, which already existed in OpenFOAM® , in order to study the sensitivity of the flow to the surface roughness in ABL flows. The simulated results obtained using the two wall functions are compared against the wind-tunnel measurements. It is shown that LES using the implemented wall function produces overall satisfactory results on the turbulent flow over the two-dimensional hills. The prediction of the flow separation and reattachment-length for the steeper hill is closer to the measurements than the other numerical studies reported in the past for the same hill geometry. The field measurement campaign performed over the Bolund hill provides the most recent field-experiment dataset for the mean flow and the turbulence properties. A number of research groups have simulated the wind flows over the Bolund hill. Due to the challenging features of the hill such as the almost vertical hill slope, it is considered as an ideal experimental test case for validating micro-scale CFD models for wind energy applications. In this work, the simulated results obtained for two wind directions are compared against the field measurements. It is shown that the present LES can reproduce the complex turbulent wind flow structures over a complicated terrain such as the Bolund hill. Especially, the present LES results show the best prediction of the turbulent kinetic energy with an average error of 24.1%, which is a 43% smaller than any other model results reported in the past for the Bolund case. Finally, the validated LES methodology is demonstrated to simulate the wind flow over the existing Muukko wind farm located in South-Eastern Finland. The simulation is carried out only for one wind direction and the results on the instantaneous and time-averaged wind speeds are briefly reported. The demonstration case is followed by discussions on the practical aspects of LES for the wind resource assessment over a realistic inland wind farm.

Multicellular spheroids of bone marrow stromal cells: a three-dimensional in vitro culture system for the study of hematopoietic cell migration

Cell fate decisions are governed by a complex interplay between cell-autonomous signals and stimuli from the surrounding tissue. In vivo cells are connected to their neighbors and to the extracellular matrix forming a complex three-dimensional (3-D) microenvironment that is not reproduced in conventional in vitro systems. A large body of evidence indicates that mechanical tension applied to the cytoskeleton controls cell proliferation, differentiation and migration, suggesting that 3-D in vitro culture systems that mimic the in vivo situation would reveal biological subtleties. In hematopoietic tissues, the microenvironment plays a crucial role in stem and progenitor cell survival, differentiation, proliferation, and migration. In adults, hematopoiesis takes place inside the bone marrow cavity where hematopoietic cells are intimately associated with a specialized three 3-D scaffold of stromal cell surfaces and extracellular matrix that comprise specific niches. The relationship between hematopoietic cells and their niches is highly dynamic. Under steady-state conditions, hematopoietic cells migrate within the marrow cavity and circulate in the bloodstream. The mechanisms underlying hematopoietic stem/progenitor cell homing and mobilization have been studied in animal models, since conventional two-dimensional (2-D) bone marrow cell cultures do not reproduce the complex 3-D environment. In this review, we will highlight some of the mechanisms controlling hematopoietic cell migration and 3-D culture systems.