Finite-size scaling analysis of the avalanches in the three-dimensional Gaussian random-field Ising model with metastable dynamics

A numerical study is presented of the third-dimensional Gaussian random-field Ising model at T=0 driven by an external field. Standard synchronous relaxation dynamics is employed to obtain the magnetization versus field hysteresis loops. The focus is on the analysis of the number and size distribution of the magnetization avalanches. They are classified as being nonspanning, one-dimensional-spanning, two-dimensional-spanning, or three-dimensional-spanning depending on whether or not they span the whole lattice in different space directions. Moreover, finite-size scaling analysis enables identification of two different types of nonspanning avalanches (critical and noncritical) and two different types of three-dimensional-spanning avalanches (critical and subcritical), whose numbers increase with L as a power law with different exponents. We conclude by giving a scenario for avalanche behavior in the thermodynamic limit.

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.

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.

Lämmönsiirron numeerinen laskenta profiloidulle kanavalle.

Työn tavoitteena oli tutkia aaltomaisen profiloinnin vaikutusta suorien jäähdytyskanavien lämmönsiirtoon ja painehäviöön. Erilaisia profiileja oli kymmenen kappaletta ja ne olivat 5mm leveitä ja 30cm pitkiä kukin. Ne laskettiin kolmeulotteisina tapauksina FINFLO-virtausratkaisijalla kolmella eri Reynoldsin luvulla, jotka vastasivat laminaarista, osittain turbulenttista ja lähes kokonaan turbulenttista virtausta. Lämmönsiirtoaine oli kuiva +30°C ilma ja profiloinnin toteutustapa oli toisiaan sivuavat ympyräkaaret kolmella erilaisella säteen arvolla ja kolmella erilaisella aallonpituuden arvolla. Lisäksi laskettiin saman levyisen tasokanavan arvot jokaisella Reynoldsin luvulla kaksiulotteisina tapauksina. Näitä profiloimattomia kanavia pidettiin referenssitapauksina. Tuloksena havaittiin että profiloimalla saadaan yksiselitteisesti suurempi lämpöteho ulos samasta tilavuudesta. Lämmönsiirtokerroin kasvaa profiloinnin avulla parhaimmillaan n. 20% käytetystä turbulenssimallista tai lämmönsiirtokertoimen määritelmästä riippumatta. Painehäviö kasvaa myös aina, mutta kitkakerroin voi hieman pienentyä. Profiilin varsinaisena hyvyyskriteerinä pidettiin lämmönsiirtokertoimen ja kitkakertoimen suhdetta h/f. Se osoittautui riippuvaksi Reynoldsin luvusta ja turbulenssimallista; ASM ja Chien k-έ -mallit ennustavat transitioetäisyyden eri tavalla. Laminaarisilla virtauksilla h/f :n vaihtelu oli vähäistä; suhde vaihteli vain ±5% eri profiilien kesken. ASM-mallilla havaittiin sekundääripyörteilyä, ehkä siksi että se mallintaa anisotrooppisen turbulenssin. Chien k-έ malli ennusti suuremman ja aikaisemmin alkavan turbulenttisuuden kuin ASM. Lisäksi havaittiin mm. että tietyillä profiileilla muodostuu kanavan kapeimpaan kohtaan selvä virtausnopeuden paikallinen minimi seinämän läheisyyden takia.

Radiaalikompressorin spiraalivirtauksen numeerinen mallinnus

Diplomityössä mallinnetaan numeerisesti radiaalikompressorin spiraalin virtaus. Spiraalin tarkoituksena radiaalikompressorissa on kerätä tasaisesti virtaus diffuusorin kehältä. Spiraaliin on viime aikoina kiinnitetty enemmän huomiota, koska on havaittu, että kompressorin hyötysuhdetta voidaan parantaa spiraalia optimoimalla. Spiraalin toimintaa tarkastellaan kolmella eri massavirralla. Työn alussa käsitellään spiraalin toimintaperiaatteita. Numeerisena ratkaisijana käytetään Teknillisessä korkeakoulussa kehitettyä FIN-FLO -koodia. FINFLO -laskentaohjelmassa ratkaistaan Navier-Stokes yhtälöt kolmeulotteiselle laskenta-alueelle. Diskretointi perustuu kontrollitilavuus menetelmään. Työssä käsitellään laskentakoodin toimintaperiaatteitta. Turbulenssia mallinnetaan algebrallisella Baldwin-Lomaxin ja kahden yhtälön Chienin k-e turbulenssimalleilla. Laskentatuloksia verrataan Lappeenrannan teknillisessä korkeakoulussa tehtyihin mittauksiin kyseessä olevasta kompressorin spiraalista. Myös eri turbulenssimalleilla ja hilatasoilla saatuja tuloksia verrataan keskenään. Laskentatuloksien jälkikäsittelyä varten ohjelmoitiin neljä eri tietokoneohjelmaa. Laskennalla pyritään saamaan lisäselvyyttä virtauksen käyttäytymiseen spiraalissa ja erityisesti ns. kielen alueella. Myös kahden eri turbulenssimallin toimivuutta kompressorin numeerisessa mallinnuksessa tutkitaan.

A Farewell to Flat Biology. Three-dimensional Cell Culture Models in Cancer Drug Target Identification and Validation

Cells of epithelial origin, e.g. from breast and prostate cancers, effectively differentiate into complex multicellular structures when cultured in three-dimensions (3D) instead of conventional two-dimensional (2D) adherent surfaces. The spectrum of different organotypic morphologies is highly dependent on the culture environment that can be either non-adherent or scaffold-based. When embedded in physiological extracellular matrices (ECMs), such as laminin-rich basement membrane extracts, normal epithelial cells differentiate into acinar spheroids reminiscent of glandular ductal structures. Transformed cancer cells, in contrast, typically fail to undergo acinar morphogenic patterns, forming poorly differentiated or invasive multicellular structures. The 3D cancer spheroids are widely accepted to better recapitulate various tumorigenic processes and drug responses. So far, however, 3D models have been employed predominantly in the Academia, whereas the pharmaceutical industry has yet to adopt a more widely and routine use. This is mainly due to poor characterisation of cell models, lack of standardised workflows and high throughput cell culture platforms, and the availability of proper readout and quantification tools. In this thesis, a complete workflow has been established entailing well-characterised 3D cell culture models for prostate cancer, a standardised 3D cell culture routine based on high-throughput-ready platform, automated image acquisition with concomitant morphometric image analysis, and data visualisation, in order to enable large-scale high-content screens. Our integrated suite of software and statistical analysis tools were optimised and validated using a comprehensive panel of prostate cancer cell lines and 3D models. The tools quantify multiple key cancer-relevant morphological features, ranging from cancer cell invasion through multicellular differentiation to growth, and detect dynamic changes both in morphology and function, such as cell death and apoptosis, in response to experimental perturbations including RNA interference and small molecule inhibitors. Our panel of cell lines included many non-transformed and most currently available classic prostate cancer cell lines, which were characterised for their morphogenetic properties in 3D laminin-rich ECM. The phenotypes and gene expression profiles were evaluated concerning their relevance for pre-clinical drug discovery, disease modelling and basic research. In addition, a spontaneous model for invasive transformation was discovered, displaying a highdegree of epithelial plasticity. This plasticity is mediated by an abundant bioactive serum lipid, lysophosphatidic acid (LPA), and its receptor LPAR1. The invasive transformation was caused by abrupt cytoskeletal rearrangement through impaired G protein alpha 12/13 and RhoA/ROCK, and mediated by upregulated adenylyl cyclase/cyclic AMP (cAMP)/protein kinase A, and Rac/ PAK pathways. The spontaneous invasion model tangibly exemplifies the biological relevance of organotypic cell culture models. Overall, this thesis work underlines the power of novel morphometric screening tools in drug discovery.

Ring Dark Solitons in Toroidal Bose-Einstein Condensates

In this Thesis, we study various aspects of ring dark solitons (RDSs) in quasi-two-dimensional toroidally trapped Bose-Einstein condensates, focussing on atomic realisations thereof. Unlike the well-known planar dark solitons, exact analytic expressions for RDSs are not known. We address this problem by presenting exact localized soliton-like solutions to the radial Gross-Pitaevskii equation. To date, RDSs have not been experimentally observed in cold atomic gases, either. To this end, we propose two protocols for their creation in experiments. It is also currently well known that in dimensions higher than one, (ring) dark solitons are susceptible, in general, to an irreversible decay into vortex-antivortex pairs through the snake instability. We show that the snake instability is caused by an unbalanced quantum pressure across the soliton's notch, linking the instability to the Bogoliubov-de Gennes spectrum. In particular, if the angular symmetry is maintained (or the toroidal trapping is restrictive enough), we show that the RDS is stable (long-lived with a lifetime of order seconds) in two dimensions. Furthermore, when the decay does take place, we show that the snake instability can in fact be reversible, and predict a previously unknown revival phenomenon for the original (many-)RDS system: the soliton structure is recovered and all the point-phase singularities (i.e. vortices) disappear. Eventually, however, the decay leads to an example of quantum turbulence; a quantum example of the laminar-to-turbulent type of transition.

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.

Response curves of deterministic and probabilistic cellular automata in one and two dimensions

One of the most important problems in the theory of cellular automata (CA) is determining the proportion of cells in a specific state after a given number of time iterations. We approach this problem using patterns in preimage sets - that is, the set of blocks which iterate to the desired output. This allows us to construct a response curve - a relationship between the proportion of cells in state 1 after niterations as a function of the initial proportion. We derive response curve formulae for many two-dimensional deterministic CA rules with L-neighbourhood. For all remaining rules, we find experimental response curves. We also use preimage sets to classify surjective rules. In the last part of the thesis, we consider a special class of one-dimensional probabilistic CA rules. We find response surface formula for these rules and experimental response surfaces for all remaining rules.

Finite-size scaling analysis of the avalanches in the three-dimensional Gaussian random-field Ising model with metastable dynamics

A numerical study is presented of the third-dimensional Gaussian random-field Ising model at T=0 driven by an external field. Standard synchronous relaxation dynamics is employed to obtain the magnetization versus field hysteresis loops. The focus is on the analysis of the number and size distribution of the magnetization avalanches. They are classified as being nonspanning, one-dimensional-spanning, two-dimensional-spanning, or three-dimensional-spanning depending on whether or not they span the whole lattice in different space directions. Moreover, finite-size scaling analysis enables identification of two different types of nonspanning avalanches (critical and noncritical) and two different types of three-dimensional-spanning avalanches (critical and subcritical), whose numbers increase with L as a power law with different exponents. We conclude by giving a scenario for avalanche behavior in the thermodynamic limit.