Der vorliegende Bericht gibt einen zusammenfassenden Überblick der im Vorhaben erreichten Ergebnisse. Ziel war die Qualifikation von CFD-Methoden für Zweiphasenströmungen mit Phasenüber¬gang. Dafür werden neuartige experimentelle Daten benötigt. Diese können an der TOPFLOW-Anlage des HZDR generiert werden, da die Anlage Experimente in für die Reaktorsicher-heits¬forschung relevanten Skalen und Parametern mit innovativen Messtechniken verbindet. Die experimentellen Arbeiten umfassen Untersuchungen zu Strömungen in vertikalen Rohren mit Hilfe der ultraschnellen Röntgentomographie, zu Strömungen mit und ohne Phasenübergang in einem Testbassin sowie zur Gegenstrombegrenzung in einem Heißstrangmodell. Diese werden im vorliegenden Bericht nur kurz dargestellt, da es zu allen 3 Versuchsserien ausführliche Dokumentationen in separaten Berichten gibt. Ein wichtiges Ergebnis der Arbeiten zur CFD-Qualifizierung ist der Erstellung des Baseline-Modellkonzepts sowie die Erstellung des Baseline-Modells für polydisperse Blasenströmungen. Damit wird ein wesentlicher Beitrag zur Erhöhung der Vorhersagefähigkeit von CFD-Codes auf Basis des Zwei- oder Mehr-Fluid-Modells erreicht. Das innovative Generalized Two-Phase Flow Konzept (GENTOP) zielt hingegen auf eine Erweiterung der Einsatzmöglichkeiten der Zweiphasen-CFD. In vielen Strömungen treten unterschiedlicher Morphologien der Phasen bzw. Strömungsformen parallel in einer Strömungsdomäne auf. Außerdem gibt es Übergänge zwischen diesen Morphologien. Mit dem GENTOP-Konzept wurde erstmals ein Rahmen geschaffen der die Simulation solcher Strömungen auf konsistente Art und Weise ermöglicht. Spezielle Modellentwicklungen erfolgten mit dem Ziel einer besseren Modellierung des Phasenübergangs.

No description available.

京都大学 / 0048 / 新制・課程博士 / 博士(工学) / 甲第20413号 / 工博第4350号 / 新制||工||1674(附属図書館) / 京都大学大学院工学研究科化学工学専攻 / (主査)教授 宮原 稔, 教授 田門 肇, 教授 山本 量一 / 学位規則第4条第1項該当 / Doctor of Philosophy (Engineering) / Kyoto University / DGAM

No description available.

In this thesis, I study fluctuations and rare events of time-additive observables of discrete-time Markov chains on finite state spaces. The observable of interest is the mean node connectivity visited by a random walk running on instances of an Erdős-Rényi (ER) random graph. I implement and analyze the Adaptive Power Method (APM) which converges to the driven process, a biased random walk defined through a control parameter that simulates trajectories corresponding to rare events of the observable in the original dynamics. The APM demonstrates good convergence and accurately produces the desired quantities from a single trajectory. Due to the bulk-dangling-chain structure in the ER graph, the driven process seems to undergo a dynamical phase transition (DPT) for infinitely large graphs, meaning the behavior of the trajectories changes abruptly as the control parameter is varied. Observations show that the random walk visits two distinct phases, being de-localized in the bulk or localized in the chain. Through two simpler models capturing the bulk-dangling-chain property of the ER graph I study how the DPT occurs as the graph size increases. I observe that the trajectories of the driven process near the transition show intermittent behavior between the two phases. The diverging time scale of the DPT is found to be the average time that the random walk spends in a phase before it transitions to the other one. On the ER graph the trajectories are also intermittent but the form of the time scaling remains open due to computational limits on the graph size.

Patman is a fluorescent membrane probe related to Laurdan. The structural distinctions between the two probes are the lengths of the aliphatic tails (eleven carbons in Laurdan and fifteen in Patman) and the presence of a trimethylammonium group on Patman that produces a positively-charged head. Preliminary studies exploring Patman as a probe to detect membrane properties during apoptosis revealed that the fluorescence intensity of two edges of the emission spectrum (435 and 500 nm) stabilizes at different rates as the probe binds to the cell membrane. To test whether these differences represent dissimilarities in probe binding to ordered and disordered domains, experiments were conducted to monitor Patman equilibration with bilayers composed of various mixtures of saturated and unsaturated phosphatidylcholines at temperatures above, at, and below the main thermotropic phase transition. In general, Patman equilibrated more rapidly with bilayers in the liquid-disordered phase than in the solid-ordered phase. With solid phase membranes, the fluorescence stabilized faster at 500 nm than at 435 nm. Similar, yet more subtle, results occurred in the lipid disordered phase. In contrast, the situation was reversed at the phase transition temperature; equilibration was faster at 435 nm than at 500 nm. To determine whether these results reflected specific properties of Patman, the experiments were repeated with Laurdan, and several distinctions were observed. First, equilibration with solid phase lipids was faster than for Patman and not different from equilibration with the fluid phase. Second, differences in rates between the two wavelengths were less than with Patman for solid phase membranes but greater than with Patman for melted bilayers. Third, at the phase transition temperature, the difference in equilibration rates was the opposite of the result obtained with Patman. Computer simulations were used to assist with interpretation of these results. The data suggest that both probes bind superficially to the membrane before incorporating among the lipid molecules. Once within the membrane, Patman localizes to at least two distinct depths within the bilayer. Probe molecules in the shallow, more hydrated position favor 500 nm emission and those occupying a deeper, dehydrated position emit primarily at 435 nm. Laurdan's equilibration additionally represents movement of the probe between leaflets and multiple bilayers.

We investigated the antiferromagnetic phase transition in the frustrated and multiferroic hexagonal manganites h-YMnO3 (YMO) and h-(Y0.98Eu0.02)MnO3 (YEMO). Elastic neutron scattering was used to study, in detail, the phase transition in YMO and YEMO under zero pressure and in YMO under a hydrostatic pressure of 1.5 GPa. Under conditions of zero pressure, we found critical temperatures of TN = 71.3(1) K and 72.11(5) K and the critical exponent 0.22(2) and b = 0.206(3), for YMO and YEMO, respectively. This is in agreement with earlier work by Roessli et al. Under an applied hydrostatic pressure of 1.5 GPa, the ordering temperature increased to TN = 75.2(5) K, in agreement with earlier reports, while b was unchanged. Inelastic neutron scattering was used to determine the size of the anisotropy spin wave gap close to the phase transition. From spin wave theory, the gap is expected to close with a critical exponent, b0, identical to the order parameter b. Our results indicate that the gap in YEMO indeed closes at TN = 72.4(3) K with b0 = 0.24(2), while the in-pressure gap in YMO closes at 75.2(5) K with an exponent of b0 = 0.19(3). In addition, the low temperature anisotropy gap was found to have a slightly higher absolute value under pressure. The consistent values obtained for b in the two systems support the likelihood of a new universality class for triangular, frustrated antiferromagnets.

This dissertation develops a computational method for simulating laser-induced cavitation and investigates the mechanism behind the formation of non-spherical bubbles induced by long-pulsed lasers. The proposed computational method accounts for the laser emission and absorption, phase transition, and the dynamics and thermodynamics of a two-phase fluid flow. In this new method, the model combines the Navier-Stokes (NS) equations for a compressible inviscid two-phase fluid flow, a new laser radiation equation, and a novel local thermodynamic model of phase transition. The Navier-Stokes equations are solved using the FInite Volume method with Exact two-phase Riemann solvers (FIVER). Following this method, numerical fluxes across phase boundaries are computed by constructing and solving one-dimensional bi-material Riemann problems. The new laser radiation equation is derived by customizing the radiative transfer equation (RTE) using the special properties of laser, including monochromaticity, directionality, high intensity, and a measurable focusing or diverging angle. An embedded boundary finite volume method is developed to solve the laser radiation equation on the same mesh created for the NS equations. The fluid mesh usually does not resolve the boundary and propagation directions of the laser beam, leading to the challenges of imposing the boundary conditions on the laser domain. To overcome this challenge, ghost nodes outside the laser domain are populated by mirroring and interpolation techniques. The existence and uniqueness of the solution are proved for the two-dimensional case, leveraging the special geometry of the laser domain. The method is up to second-order accuracy, which is also proved, and verified using numerical tests. A method of latent heat reservoir is developed to predict the onset of vaporization, which accounts for the accumulation and release of latent heat. In this work, the localized level set method is employed to track the bubble surface. Furthermore, the continuation of phase transition is possible in laser-induced cavitation problems, especially for long-pulsed lasers. A method of local correction and reinitialization is developed to account for continuous phase transitions. Several numerical tests are presented to verify the convergence of these methods. This multiphase laser-fluid coupled computational model is employed to simulate the formation and expansion of bubbles with different shapes induced by different long-pulsed lasers. The simulation results show that the computational method can capture the key phenomena in the laser-induced cavitation problems, including non-spherical bubble expansion, shock waves, and the ``Moses effect''. Additionally, the observed complex non-spherical shapes of vapor bubbles generated by long-pulsed laser reflect some characteristics (e.g., direction, width) of the laser beam. The dissertation also investigates the relation between bubble shapes and laser parameters and explores the transition between two commonly observed shapes -- namely, a rounded pear-like shape and an elongated conical shape -- using the proposed computational model. Two laboratory experiments are simulated, in which Holmium:YAG and Thulium fiber lasers are used respectively to generate bubbles of different shapes. In both cases, the predicted bubble nucleation and morphology agree reasonably well with the experimental observation. The full-field results of laser radiance, temperature, velocity, and pressure are analyzed to explain bubble dynamics and energy transmission. It is found that due to the lasting energy input, the vapor bubble's dynamics is driven not only by advection, but also by the continued vaporization at its surface. Vaporization lasts less than 1 microsecond in the case of the pear-shaped bubble, compared to over 50 microseconds for the elongated bubble. It is thus hypothesized that the bubble's morphology is determined by a competition between the speed of bubble growth due to advection and continuous vaporization. When the speed of advection is higher than that of vaporization, the bubble tends to grow spherically. Otherwise, it elongates along the laser beam direction. To test this hypothesis, the two speeds are defined analytically using a model problem and then estimated for the experiments using simulation results. The results support the hypothesis and also suggest that when the laser's power is fixed, a higher laser absorption coefficient and a narrower beam facilitate bubble elongation. / Doctor of Philosophy / Laser-induced cavitation is a process where laser beams create bubbles in a liquid. This phenomenon is widely applied in research and microfluidic applications for precise control of bubble dynamics. It also naturally occurs in various laser-based processes involving liquid environments. Understanding laser-induced cavitation is important for enhancing the effectiveness and safety of related technologies. However, experimental studies encounter limitations, highlighting the development of numerical methods to advance the understanding of laser-induced cavitation. The laser-induced cavitation can be roughly described as localized boiling through thermal radiation. The detailed physics involves the absorption of laser light by a liquid, the formation of vapor bubbles due to localized heating, and the dynamics of both the bubbles and the surrounding liquid. The first part of the dissertation introduces a new computational method for modeling these phenomena. The dynamics of the two-phase flow are modeled by the Navier-Stokes equations, which are solved using the FInite Volume method with Exact two-phase Riemann solvers (FIVER). The absorption of the laser light is modeled by a new laser radiation equation, which is derived from laser energy conservation and special properties of the laser. An embedded boundary finite volume method is developed to solve this equation on the same mesh created for the NS equations. Additionally, a method of latent heat reservoir is developed to predict the onset of vaporization. In this work, the level set method is employed to track the bubble surface, and a method of local correction and reinitialization is developed to account for possible continuous phase transitions. After developing this new method, several test cases are simulated. The simulation results show that the method can capture the key phenomena in the laser-induced cavitation problems, including the absorption of laser light, non-spherical bubble expansion, and shock waves. When the laser pulse is comparable to or longer than the acoustic time scale (long-pulsed laser), vapor bubbles generated often have complex non-spherical shapes. The bubble shapes reflect some characteristics (e.g., direction, width) of the laser beam. The second part of the dissertation investigates the relation between bubble shapes and laser parameters. Two laboratory experiments are simulated, in which two different lasers are used to generate bubbles of different shapes, namely, a rounded pear-like shape and an elongated conical shape. In both cases, the simulated bubbles exhibit shapes and sizes that reasonably match the experimental results. The simulation results of temperature, pressure, and velocity fields are analyzed to explain bubble dynamics and energy transmission. The analysis shows that the expansion of bubbles induced by long-pulsed lasers is determined not only by advection but also by the continued vaporization at its surface. Vaporization lasts less than $1$ microsecond in the case of the pear-shaped bubble, compared to over $50$ microseconds for the elongated bubble. It is thus hypothesized that the bubble expansion is determined by a competition between the speed of bubble growth due to advection and continuous vaporization. When the speed of advection is higher than that of vaporization, the bubble tends to grow spherically. Otherwise, it elongates along the laser beam direction. To test this hypothesis, the two speeds are defined analytically using a model problem and then estimated for the experiments using simulation results. The results support the hypothesis and also suggest that when the laser's power is fixed, a higher laser absorption coefficient and a narrower beam facilitate bubble elongation.

This dissertation is the collection of a progressive research on the topic of topological quantum computation and information with the focus on the error threshold of the well-known models such as the unpaired Majorana, the toric code, and the planar code. We study the basics of quantum computation and quantum information, and in particular quantum error correction. Quantum error correction provides a tool for enhancing the quantum computation fidelity in the noisy environment of a real world. We begin with a brief introduction to stabilizer codes. The stabilizer formalism of the theory of quantum error correction gives a well-defined description of quantum codes that is used throughout this dissertation. Then, we turn our attention to a quite new subject, namely, topological quantum codes. Topological quantum codes take advantage of the topological characteristics of a physical many-body system. The physical many-body systems studied in the context of topological quantum codes are of two essential natures: they either have intrinsic interaction that self-corrects errors, or are actively corrected to be maintained in a desired quantum state. Examples of the former are the toric code and the unpaired Majorana, while an example for the latter is the surface code. A brief introduction and history of topological phenomena in condensed matter is provided. The unpaired Majorana and the Kitaev toy model are briefly explained. Later we introduce a spin model that maps onto the Kitaev toy model through a sequence of transformations. We show how this model is robust and tolerates local perturbations. The research on this topic, at the time of writing this dissertation, is still incomplete and only preliminary results are represented. As another example of passive error correcting codes with intrinsic Hamiltonian, the toric code is introduced. We also analyze the dynamics of the errors in the toric code known as anyons. We show numerically how the addition of disorder to the physical system underlying the toric code slows down the dynamics of the anyons. We go further and numerically analyze the presence of time-dependent noise and the consequent delocalization of localized errors. The main portion of this dissertation is dedicated to the surface code. We study the surface code coupled to a non-interacting bosonic bath. We show how the interaction between the code and the bosonic bath can effectively induce correlated errors. These correlated errors may be corrected up to some extend. The extension beyond which quantum error correction seems impossible is the error threshold of the code. This threshold is analyzed by mapping the effective correlated error model onto a statistical model. We then study the phase transition in the statistical model. The analysis is in two parts. First, we carry out derivation of the effective correlated model, its mapping onto a statistical model, and perform an exact numerical analysis. Second, we employ a Monte Carlo method to extend the numerical analysis to large system size. We also tackle the problem of surface code with correlated and single-qubit errors by an exact mapping onto a two-dimensional Ising model with boundary fields. We show how the phase transition point in one model, the Ising model, coincides with the intrinsic error threshold of the other model, the surface code.

No description available.