Recent developments in estimation of the survivability of a U.S. Navy transport barge in random seas are extended to improve accuracy. The single Degree-of-Freedom (DOF) model of a extreme roll response of a barge used in previous research is replaced by a 3-DOF roll-heave-sway model to include linear and nonlinear static and kinematic coupling between roll, sway and heave. The predominant nonlinearity in the model arises in an improved approximation of the roll righting moment and heave buoyant restoring force by coupling roll with heave. Kinematic coupling is introduced by allowing extreme displacements and rotations in the barge response. System coefficients in the 3-DOF roll-heave-sway model and a simpler 2-DOF roll-heave model are identified by comparing time domain simulations with measured physical model tests of barge motions. Predictions of the 3-DOF and 2-DOF models are compared to measured test data for the case of random waves. Monte Carlo simulations of the equations of motions are performed to predict the reliability of the barge in an operational sea state for a specified mission duration. Use of parallel computer processing is found to make this a viable option for stability estimations as we move into the next century. The stochastic nature of the ocean waves are modeled via filtered white noise. Estimations of the joint probability of the barge responses are presented after application of density estimation kernels. Both the 3-DOF roll-heave-sway model and 2-DOF roll-heave model are tested and compared. Last, examples are provided of some observed nonlinear behavior of the barge motions for variation in damping or ocean wave amplitude. Transient and intermittent chaotic responses are observed for deterministic input waves and quasiperiodic cases are illustrated. / Graduation date: 1996

Estrogen receptor α (ERα) is a member of a family of ligand-modulated transcription factors that have come to be known as nuclear receptors. ERα mediates the action of estrogens and plays an integral role in a wide range of physiological processes ranging from embryonic development and morphogenesis to reproduction to cardiovascular health. Not surprisingly, malfunction of the estrogen system is associated with a host of pathological conditions such as osteoporosis, heart disease and most notably breast cancer. Essential to its functioning as a transcription factor are specific protein-DNA interactions which are mediated by the binding of the DNA-binding (DB) domain of ERα to particular DNA sequences located within target gene promoters called estrogen response elements (EREs). Here, using a diverse array of biophysical techniques, including in particular isothermal titration calorimetry coupled with molecular modeling and semi-empirical analysis, I provide new insights into the ERα-DNA interaction in thermodynamic and structural terms. My data show that the binding of the DB domain of ERα to DNA is coupled to protonation at two specific amino acids, H196 and E203. Protonation of these residues is non-trivial and is required for high affinity binding. Amino acid sequence alignment of the DB domains of the NR family suggests that this may be a hallmark feature common to the functioning of all nuclear receptors. Furthermore, I demonstrate that the DB domain can tolerate all single nucleotide substitutions within the ERE and bind in the physiologically relevant nanomolar to micromolar range. Comparative thermodynamic analysis reveals that the DB domain binds to these ERE sequences utilizing a considerable range of energetic signatures such that any one thermodynamic component of binding is not predictive of associated affinity. In addition, it is shown that nucleotide substitution results in significant changes in secondary and three-dimensional features of the oligonucleotides and may impact binding affinity. Finally, I demonstrate that the zinc-finger of the DB domain of ERα is relatively promiscuous and can accommodate several heavy-metal divalent cations. Other than zinc, only DB domains reconstituted with cobalt, cadmium and mercury were capable of binding DNA. Incorporation of the metals resulted in a wide range of CD spectroscopic features which were found not to be predictive of DNA binding capacity. Thus, isostructure does not equate to isofunction in the case of metal reconstituted DB domain of ERα. This analysis suggests that metal coordination is not likely to be required for domain folding, but rather is required to bind DNA. Taken together, this thesis provides novel insights into the physicochemical basis of a key protein-DNA interaction essential to human health and disease. My studies bear the potential to impact the development of novel therapies harboring greater efficacy coupled with lower toxicity for the treatment of disease.

The experiments performed at the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Lab have discovered a state of matter called the strongly coupled quark-gluon plasma (sQGP). The strong coupling has limited the ability of the standard theory to describe such matter, namely Quantum Chromodynamics (QCD). However, string theory's anti-de Sitter/conformal field theory (AdS/CFT) correspondence has provided a new way to study the situation and in an analytical manner. So far, hydrodynamic properties of RHIC's plasma, such as elliptic flow and longitudinal expansion, have been seen to follow from classical supergravity calculations. In this dissertation I discuss some of the field's development as well as the research done by the author and collaborators.

When an electrically conducting fluid moves through a magnetic field, fluid mechanics and electromagnetism are coupled. This interaction is the object of magnetohydrodynamics, a discipline which covers a wide range of applications, from electromagnetic processing to plasma- and astro-physics. In this dissertation, the attention is restricted to turbulent liquid metal flows, typically encountered in steel and aluminium industries. Velocity measurements in such flows are extremely challenging because liquid metals are opaque, hot and often corrosive. Therefore, non-intrusive measurement devices are essential. One of them is the Lorentz force flowmeter. Its working principle is based on the generation of a force acting on a charge, which moves in a magnetic field. Recent studies have demonstrated that this technique can measure efficiently the mean velocity of a liquid metal. In the existing devices, however, the measurement depends on the electrical conductivity of the fluid. In this work, a novel version of this technique is developed in order to obtain measurements that are independent of the electrical conductivity. This is particularly appealing for metallurgical applications, where the conductivity often fluctuates in time and space. The study is entirely numerical and uses a flexible computational method, suitable for industrial flows. In this framework, the cost of numerical simulations increases drastically with the level of turbulence and the geometry complexity. Therefore, the simulations are commonly unresolved. Large eddy simulations are then very promising, since they introduce a subgrid model to mimic the dynamics of the unresolved turbulent eddies. The first part of this dissertation focuses on the quality and reliability of unresolved numerical simulations. The attention is drawn on the ambiguity that may arise when interpretating the results. Owing to coarse resolutions, numerical errors affect the performances of the discrete model, which in turn looses its physical meaning. In this work, a novel implementation of the turbulent strain rate appearing in the models is proposed. As opposed to its usual discretisation, the present strain rate is in accordance with the discrete equations of motion. Two types of flow are considered: decaying turbulence located far from boundaries, and turbulent flows between two parallel and infinite walls. Particular attention is given to the balance of resolved kinetic energy, in order to assess the role of the model. The second part of this dissertation deals with a novel version of Lorentz force flowmeters, consisting in one or two coils placed around a circular pipe. The forces acting on each coil are recorded in time as the liquid metal flows through the pipe. It is highlighted that the auto- or cross-correlation of these forces can be used to determine the flowrate. The reliability of the flowmeter is first investigated with a synthetic velocity profile associated to a single vortex ring, which is convected at a constant speed. This configuration is similar to the movement of a solid rod and enables a simple analysis of the flowmeter. Then, the flowmeter is applied to a realistic three-dimensional turbulent flow. In both cases, the influence of the geometrical parameters of the coils is systematically assessed.

In the context of cosmological structure formation sheets, filaments and eventually halos form due to gravitational instabilities. It is noteworthy, that at all times, the majority of the baryons in the universe does not reside in the dense halos but in the filaments and the sheets of the intergalactic medium. While at higher redshifts of z > 2, these baryons can be detected via the absorption of light (originating from more distant sources) by neutral hydrogen at temperatures of T ~ 10^4 K (the Lyman-alpha forest), at lower redshifts only about 20 % can be found in this state. The remain (about 50 to 70 % of the total baryons mass) is unaccounted for by observational means. Numerical simulations predict that these missing baryons could reside in the filaments and sheets of the cosmic web at high temperatures of T = 10^4.5 - 10^7 K, but only at low to intermediate densities, and constitutes the warm-hot intergalactic medium (WHIM). The high temperatures of the WHIM are caused by the formation of shocks and the subsequent shock-heating of the gas. This results in a high degree of ionization and renders the reliable detection of the WHIM a challenging task. Recent high-resolution hydrodynamical simulations indicate that, at redshifts of z ~ 2, filaments are able to provide very massive galaxies with a significant amount of cool gas at temperatures of T ~ 10^4 K. This could have an important impact on the star-formation in those galaxies. It is therefore of principle importance to investigate the particular hydro- and thermodynamical conditions of these large filament structures. Density and temperature profiles, and velocity fields, are expected to leave their special imprint on spectroscopic observations. A potential multiphase structure may act as tracer in observational studies of the WHIM. In the context of cold streams, it is important to explore the processes, which regulate the amount of gas transported by the streams. This includes the time evolution of filaments, as well as possible quenching mechanisms. In this context, the halo mass range in which cold stream accretion occurs is of particular interest. In order to address these questions, we perform particular hydrodynamical simulations of very high resolution, and investigate the formation and evolution of prototype structures representing the typical filaments and sheets of the WHIM. We start with a comprehensive study of the one-dimensional collapse of a sinusoidal density perturbation (pancake formation) and examine the influence of radiative cooling, heating due to an UV background, thermal conduction, and the effect of small-scale perturbations given by the cosmological power spectrum. We use a set of simulations, parametrized by the wave length of the initial perturbation L. For L ~ 2 Mpc/h the collapse leads to shock-confined structures. As a result of radiative cooling and of heating due to an UV background, a relatively cold and dense core forms. With increasing L the core becomes denser and more concentrated. Thermal conduction enhances this trend and may lead to an evaporation of the core at very large L ~ 30 Mpc/h. When extending our simulations into three dimensions, instead of a pancake structure, we obtain a configuration consisting of well-defined sheets, filaments, and a gaseous halo. For L > 4 Mpc/h filaments form, which are fully confined by an accretion shock. As with the one-dimensional pancakes, they exhibit an isothermal core. Thus, our results confirm a multiphase structure, which may generate particular spectral tracers. We find that, after its formation, the core becomes shielded against further infall of gas onto the filament, and its mass content decreases with time. In the vicinity of the halo, the filament's core can be attributed to the cold streams found in other studies. We show, that the basic structure of these cold streams exists from the very beginning of the collapse process. Further on, the cross section of the streams is constricted by the outwards moving accretion shock of the halo. Thermal conduction leads to a complete evaporation of the cold stream for L > 6 Mpc/h. This corresponds to halos with a total mass higher than M_halo = 10^13 M_sun, and predicts that in more massive halos star-formation can not be sustained by cold streams. Far away from the gaseous halo, the temperature gradients in the filament are not sufficiently strong for thermal conduction to be effective. / Im Rahmen der kosmologischen Strukturbildung entstehen durch Gravitationsinstabilitäten Flächen, Filamente und schließlich Halos. Interessanterweise befinden sich zu jedem Zeitpunkt der kosmologischen Entwicklung der Großteil der Baryonen nicht in den Halos, sondern in den Filamenten und Ebenen des intergalaktischen Mediums. Während diese Baryonen bei höheren Rotverschiebungen (z ~ 2) noch in Form durch die Absorbtion von Licht (von weit entfernteren Quellen) durch neutralen Wasserstoff bei einer Temperatur von T ~ 10^4 K beobachtbar sind (Lyman-Alpha Wald), gilt dies bei niedrigeren Rotverschiebungen für nur noch ca. 20 % der Baryonen. Der überwiegende Teil (ca. 50-70 % der gesamten baryonischen Masse) sind bisher noch nicht direkt beobachtbar. Numerische Simulationen sagen jedoch voraus, das sich diese Baryonen in den Filamenten und Flächen des kosmischen Netzes befinden. Die entsprechende Gasverteilung zeichnet sich durch hohe Temperaturen T = 10^5 - 10^7 K und geringe bis mittlere Dichten aus und wird als warm-heißes intergalaktisches Medium (WHIM) bezeichnet. Die hohen Temperaturen entstehen in Folge der Bildung von Stoßwellen und der darauf folgenden Erhitzung des Gases (shock-heating). Das WHIM ist daher hochgradig ionisiert und sein verlässlicher Nachweis stellt eine große Herausforderung für die beobachtende Kosmologie dar. Neuere hydrodynamische Simulationen zeigen, dass sich bei höheren Rotverschiebungen von z ~ 2 Gasströmungen entlang der Filamente bilden, die massive Galaxien mit erheblichen Mengen an relativ kaltem Gas (T ~ 10^4 K) versorgen können. Dies hätte einen erheblichen Einfluss auf die Sternentstehung in diesen Galaxien. Es ist daher von grundsätzlichem Interesse, die spezifischen hydro- und thermodynamischen Bedingungen in den Strukturen des WHIM zu untersuchen. Sowohl Dichte- und Temperaturprofile als auch Geschwindigkeitsfelder prägen spektroskopische Beobachtungen. Eine mögliche Mehrphasenstruktur des WHIM könnte daher als Indikator in beobachtenden Studien dienen. Im Zusammenhang mit den kalten Strömen ist es besonders interessant, Prozesse zu untersuchen die den Zufluss von kaltem Gas zu den Galaxien regulieren. Dies umfasst die Zeitentwicklung des Anteils an kaltem Gas in den Filamenten, sowie mögliche Mechanismen, die zum Versiegen des Zuflusses von kaltem Gas auf die Galaxienscheibe führen. Um diese Zusammenhänge zu erforschen, führen wir spezielle hydrodynamische Simulationen mit sehr hoher Auflösung durch, die zu ausgewählten, wohldefinierten Strukturen führen, die das WHIM charakterisieren. Wir beginnen mit einer ausführlichen Untersuchung des eindimensionalen Kollaps einer sinusförmigen Störung (pancake formation). Hierbei untersuchen wir den Einfluss von Strahlungkühlung, Heizung durch den intergalaktischen UV Hintergrund, Wärmeleitung, sowie von kleinskaligen Störungen, welche dem kosmologischen Störungsspektrum folgen. Wir benutzen hierbei eine Reihe von Simulationen, welche die Längenskala der anfänglichen Störung L als Parameter verwenden. Für L ~ 2 Mpc/h führt der Kollaps zur Ausbildung einer Stoßwelle. Zusätzlich entsteht als Folge der Strahlungskühlung und der Heizung durch den UV Hintergrund ein relativ dichter und kalter isothermer Kern. Mit ansteigendem L wird dieser Kern dichter und kompakter. Durch Wärmeleitung reduziert sich die räumliche Ausdehnung des Kerns. Für L ~ 30 Mpc/h führt dies zu einem Verschwinden des Kerns. Mit der Erweiterung unserer Methodik auf dreidimensionale Simulationen, entsteht nun eine Konfiguration, welche aus wohldefinierten Flächen, Filamenten und einem gasförmigen Halo besteht. Für L > 4 Mpc/h, erhalten wir Filamente, die vollständig durch Akkretionsschocks begrenzt sind. Wie in unseren eindimensionalen Simulationen weisen auch sie einen isothermen Kern auf. Dies legt nahe, dass das WHIM eine Mehrphasenstruktur besitzt und mögliche Spektralsignaturen erzeugen kann. Nach seiner Entstehung ist der Kern gegen weiteren Zufluss von Gas abgeschirmt und seine Masse reduziert sich mit der Zeit. In der direkten Umgebung des Halos entspricht der Kern des Filamentes den oben angesprochenen kalten Strömen. Unsere Untersuchung zeigt, dass diese während der gesamten Entwicklung des Halos existent sind. In der weiteren Entwicklung werden sie durch den expandierenden Akkretionsschock des Halos verengt. Ab einer Skala von L > 6 Mpc/h kann Wärmeleitung zu einem Verschwinden des Zustroms von kaltem Gas führen. Diese Skala entspricht Halos mit einer Gesamtmasse von M_halo = 10^13 M_sun. Galaxien, die sich in noch massiveren Halos bilden, können daher nicht durch kalte Ströme mit Gas für die Sternentstehung versorgt werden. Im Filament, weit außerhalb des gasförmigen Halos, sind die Temperaturgradienten zu klein, um effiziente Wärmeleitung zu ermöglichen.

This Thesis puts its focus on the physics of neutron stars and its description with methods of numerical relativity. In the first step, a new numerical framework the Whisky2D code will be developed, which solves the relativistic equations of hydrodynamics in axisymmetry. Therefore we consider an improved formulation of the conserved form of these equations. The second part will use the new code to investigate the critical behaviour of two colliding neutron stars. Considering the analogy to phase transitions in statistical physics, we will investigate the evolution of the entropy of the neutron stars during the whole process. A better understanding of the evolution of thermodynamical quantities, like the entropy in critical process, should provide deeper understanding of thermodynamics in relativity. More specifically, we have written the Whisky2D code, which solves the general-relativistic hydrodynamics equations in a flux-conservative form and in cylindrical coordinates. This of course brings in 1/r singular terms, where r is the radial cylindrical coordinate, which must be dealt with appropriately. In the above-referenced works, the flux operator is expanded and the 1/r terms, not containing derivatives, are moved to the right-hand-side of the equation (the source term), so that the left hand side assumes a form identical to the one of the three-dimensional (3D) Cartesian formulation. We call this the standard formulation. Another possibility is not to split the flux operator and to redefine the conserved variables, via a multiplication by r. We call this the new formulation. The new equations are solved with the same methods as in the Cartesian case. From a mathematical point of view, one would not expect differences between the two ways of writing the differential operator, but, of course, a difference is present at the numerical level. Our tests show that the new formulation yields results with a global truncation error which is one or more orders of magnitude smaller than those of alternative and commonly used formulations. The second part of the Thesis uses the new code for investigations of critical phenomena in general relativity. In particular, we consider the head-on-collision of two neutron stars in a region of the parameter space where two final states a new stable neutron star or a black hole, lay close to each other. In 1993, Choptuik considered one-parameter families of solutions, S[P], of the Einstein-Klein-Gordon equations for a massless scalar field in spherical symmetry, such that for every P > P⋆, S[P] contains a black hole and for every P < P⋆, S[P] is a solution not containing singularities. He studied numerically the behavior of S[P] as P → P⋆ and found that the critical solution, S[P⋆], is universal, in the sense that it is approached by all nearly-critical solutions regardless of the particular family of initial data considered. All these phenomena have the common property that, as P approaches P⋆, S[P] approaches a universal solution S[P⋆] and that all the physical quantities of S[P] depend only on |P − P⋆|. The first study of critical phenomena concerning the head-on collision of NSs was carried out by Jin and Suen in 2007. In particular, they considered a series of families of equal-mass NSs, modeled with an ideal-gas EOS, boosted towards each other and varied the mass of the stars, their separation, velocity and the polytropic index in the EOS. In this way they could observe a critical phenomenon of type I near the threshold of black-hole formation, with the putative solution being a nonlinearly oscillating star. In a successive work, they performed similar simulations but considering the head-on collision of Gaussian distributions of matter. Also in this case they found the appearance of type-I critical behaviour, but also performed a perturbative analysis of the initial distributions of matter and of the merged object. Because of the considerable difference found in the eigenfrequencies in the two cases, they concluded that the critical solution does not represent a system near equilibrium and in particular not a perturbed Tolmann-Oppenheimer-Volkoff (TOV) solution. In this Thesis we study the dynamics of the head-on collision of two equal-mass NSs using a setup which is as similar as possible to the one considered above. While we confirm that the merged object exhibits a type-I critical behaviour, we also argue against the conclusion that the critical solution cannot be described in terms of equilibrium solution. Indeed, we show that, in analogy with what is found in, the critical solution is effectively a perturbed unstable solution of the TOV equations. Our analysis also considers fine-structure of the scaling relation of type-I critical phenomena and we show that it exhibits oscillations in a similar way to the one studied in the context of scalar-field critical collapse. / Diese Arbeit legt seinen Schwerpunkt auf die Physik von Neutronensternen und deren Beschreibung mit Methoden der numerischen Relativitätstheorie. Im ersten Schritt wird eine neue numerische Umgebung, der Whisky2D Code entwickelt, dieser löst die relativistischen Gleichungen der Hydrodynamik in Axialymmetrie. Hierzu betrachten wir eine verbesserte Formulierung der sog. "flux conserved formulation" der Gleichungen. Im zweiten Teil wird der neue Code verwendet / um das kritische Verhalten zweier kollidierenden Neutronensternen zu untersuchen. In Anbetracht der Analogie, um Übergänge in der statistischen Physik Phase werden wir die Entwicklung der Entropie der Neutronensterne während des gesamten Prozesses betrachten. Ein besseres Verständnis der Evolution von thermodynamischen Größen, wie der Entropie in kritischer Prozess, sollte zu einem tieferen Verständnis der relativistischen Thermodynamik führen. Der Whisky2D Code, zur Lösung Gleichungen relativistischer Hydrodynamik wurde in einer „flux conserved form“ und in zylindrischen Koordinaten geschrieben. Hierdurch entstehen 1 / r singuläre Terme, wobei r der ist, die entsprechend behandelt werden müssen. In früheren Arbeiten, wird der Operator expandiert und die 1 / r spezifisch Therme auf die rechte Seite geschrieben, so dass die linke Seite eine Form annimmt, die identisch ist mit der kartesischen Formulierung. Wir nennen dies die Standard-Formulierung. Eine andere Möglichkeit ist, die Terme nicht zu expandieren, den und den 1/r Term in die Gleichung hinein zu ziehen. Wir nennen dies die Neue-Formulierung. Die neuen Gleichungen werden mit den gleichen Verfahren wie im kartesischen Fall gelöst. Aus mathematischer Sicht ist keine Unterschiede zwischen den beiden Formulierungen zu erwarten, erst die numerische Sicht zeigt die Unterschiede auf. Versuche zeigen, dass die Neue-Formulierung numerische Fehler um mehrere Größenordnungen reduziert. Der zweite Teil der Dissertation verwendet den neuen Code für die Untersuchung kritischer Phänomene in der allgemeinen Relativitätstheorie. Insbesondere betrachten wir die Kopf-auf-Kollision zweier Neutronensterne in einem Bereich des Parameter Raums, deren zwei mögliche Endzustände entweder einen neuen stabilen Neutronenstern oder ein Schwarzes Loch darstellen. Im Jahr 1993, betrachtete Choptuik Ein-Parameter-Familien von Lösungen, S [P], der Einstein-Klein-Gordon-Gleichung für ein masseloses Skalarfeld in sphärischer Symmetrie, so dass für jedes P> P ⋆, S[P] ein Schwarzes Loch enthalten ist und jedes P <P ⋆, S [P] eine Lösung ohne eine Singularität. Er studierte das numerisch Verhalten von S [P] als P → P ⋆ und stellte fest, dass die kritische Lösung, S[P ⋆], universell in dem Sinne ist, dass die Lösung nahe des kritischen Wertes unabhängig von den Anfangsdaten ist. Alle diese Phänomene haben die gemeinsame Eigenschaft, dass sobald sich P P ⋆ annähert, auch S[P] S [P ⋆] annähert und dass alle physikalischen Größen von S [P] nur noch von | P - P ⋆ | abhängen. Die erste Studie der kritischen Phänomene über den Frontalzusammenstoß von zwei Neutronensternen wurde von Jin und Suen im Jahr 2007 durchgeführt. Insbesondere untersuchten sie eine Reihe von Anfangsdaten gleicher Neutronensternmasse, mit einer idealen EOS, die aufeinander zu beschleunigt werden. Variiert wurden die Massen der Sterne, ihr Abstand, die Geschwindigkeit und die polytropen Index der EOS. Auf diese Weise konnten sie kritische Phänomen des Typ I beobachten. In weiteren Versuchen, führten sie ähnliche Simulationen frontal kollidierender Materie in Gauß Verteilungen durch. Auch in diesem Fall fanden sie Typ-I-kritisches Verhalten. Zudem führten sie eine störungstheoretische Analyse der Anfangsobjekte als auch der stabilen Endobjekte durch. Wegen der beträchtlichen Unterschiede in den Eigenfrequenzen in beiden Fallen, schlossen sie daraus, dass die kritische Lösung keine linear-gestörten Tolmann-Oppenheimer-Volkoff (TOV) Sterne im Gleichgewicht darstellen. In dieser Arbeit untersuchen wir die Dynamik der Frontalzusammenstoß zweier Neutronensterne gleicher Masse mit ähnlichem Setup wie oben besprochen. Während wir bestätigen, dass die erzeugten Objekte ebenfalls ein Typ-I-kritische Verhalten aufweisen, wiedersprechen wir der der Aussage, dass sich die kritische Lösung nicht als Gleichgewichtslösung dargestellt werden kann. In der Tat zeigen wir, dass die kritische Lösung als linear-gestörte instabile Lösung eines TOV-Sterns dargestellt werden kann. Unsere Analyse berücksichtigt auch die Feinstruktur der Skalenverhältnisse Typ-I-kritischer Phänomene und wir zeigen ebenfalle, dass hier Oszillationen auftreten, die bereits in ähnlichen Studien zum Kollaps kritischer Skalar-Feld gefunden wurden.

Liniengetriebene Winde werden durch Impulsübertrag von Photonen auf ein Plasma bei Absorption oder Streuung in zahlreichen Spektrallinien beschleunigt. Dieser Prozess ist besonders effizient für ultraviolette Strahlung und Plasmatemperaturen zwischen 10^4 K und 10^5 K. Zu den astronomischen Objekten mit liniengetriebenen Winden gehören Sterne der Spektraltypen O, B und A, Wolf-Rayet-Sterne sowie Akkretionsscheiben verschiedenster Größenordnung, von Scheiben um junge Sterne und in kataklysmischen Veränderlichen bis zu Quasarscheiben. Es ist bislang nicht möglich, das vollständige Windproblem numerisch zu lösen, also die Hydrodynamik, den Strahlungstransport und das statistische Gleichgewicht dieser Strömungen gleichzeitig zu behandeln. Die Betonung liegt in dieser Arbeit auf der Windhydrodynamik, mit starken Vereinfachungen in den beiden anderen Gebieten. <br /> Wegen persönlicher Beteiligung betrachte ich drei Themen im Detail. <br /> 1. Windinstabilität durch Dopplerde-shadowing des Gases. Die Instabilität bewirkt, dass Windgas in dichte Schalen komprimiert wird, die von starken Stoßfronten begrenzt sind. Schnelle Wolken entstehen im Raum zwischen den Schalen und stoßen mit diesen zusammen. Dies erzeugt Röntgenflashes, die die beobachtete Röntgenstrahlung heißer Sterne erklären können. <br /> 2. Wind runway durch radiative Wellen. Der runaway zeigt, warum beobachtete liniengetriebene Winde schnelle, kritische Lösungen anstelle von Brisenlösungen (oder shallow solutions) annehmen. Unter bestimmten Bedingungen stabilisiert der Wind sich auf masseüberladenen Lösungen, mit einem breiten, abbremsenden Bereich und Knicken im Geschwindigkeitsfeld. <br /> 3. Magnetische Winde von Akkretionsscheiben um Sterne oder in aktiven Galaxienzentren. Die Linienbeschleunigung wird hier durch die Zentrifugalkraft entlang korotierender poloidaler Magnetfelder und die Lorentzkraft aufgrund von Gradienten im toroidalen Feld unterstützt. Ein Wirbelblatt, das am inneren Scheibenrand beginnt, kann zu stark erhöhten Massenverlustraten führen. / Line driven winds are accelerated by the momentum transfer from photons to a plasma, by absorption and scattering in numerous spectral lines. Line driving is most efficient for ultraviolet radiation, and at plasma temperatures from 10^4 K to 10^5 K. Astronomical objects which show line driven winds include stars of spectral type O, B, and A, Wolf-Rayet stars, and accretion disks over a wide range of scales, from disks in young stellar objects and cataclysmic variables to quasar disks. It is not yet possible to solve the full wind problem numerically, and treat the combined hydrodynamics, radiative transfer, and statistical equilibrium of these flows. The emphasis in the present writing is on wind hydrodynamics, with severe simplifications in the other two areas. <br /> I consider three topics in some detail, for reasons of personal involvement. <br /> 1. Wind instability, as caused by Doppler de-shadowing of gas parcels. The instability causes the wind gas to be compressed into dense shells enclosed by strong shocks. Fast clouds occur in the space between shells, and collide with the latter. This leads to X-ray flashes which may explain the observed X-ray emission from hot stars. <br /> 2. Wind runaway, as caused by a new type of radiative waves. The runaway may explain why observed line driven winds adopt fast, critical solutions instead of shallow (or breeze) solutions. Under certain conditions the wind settles on overloaded solutions, which show a broad deceleration region and kinks in their velocity law. <br /> 3. Magnetized winds, as launched from accretion disks around stars or in active galactic nuclei. Line driving is assisted by centrifugal forces along co-rotating poloidal magnetic field lines, and by Lorentz forces due to toroidal field gradients. A vortex sheet starting at the inner disk rim can lead to highly enhanced mass loss rates.

Lipid molecules form thin biological membranes that envelop all living cells, and behave as two-dimensional liquid sheets immersed in bulk water. The interactions of such biomembranes with their environment lay the foundation of a plethora of biological processes rooted in the mesoscopic domain - length scales of 1-1000 nm and time scales of 1-1000 ns. Research in this intermediate regime has for a long time been out of reach for conventional experiments, but breakthroughs in computer simulation methods and scattering experimental techniques have made it possible to directly probe static and dynamic properties of biomembranes on these scales. Biomembranes are soft, with a relatively low energy cost of bending, and are thereby influenced by random, thermal fluctuations of individual molecules. Molecular dynamics simulations show how in-plane (density fluctuations) and out-of-plane (undulations) motions are intertwined in the bilayer in the mesoscopic domain. By novel methods, the fluctuation spectra of lipid bilayers can be calculated withdirect Fourier analysis. The interpretation of the fluctuation spectra reveals a picture where density fluctuations and undulations are most pronounced on different length scales, but coalesce in the mesoscopic regime. This analysis has significant consequences for comparison of simulation data to experiments. These new methods merge the molecular fluctuations on small wavelengths, with continuum fluctuations of the elastic membrane sheet on large wavelengths, allowing electron density profiles (EDP) and area per lipid to be extracted from simulations with high accuracy. Molecular dynamics simulations also provide insight on the small-wavelength dynamics of lipid membranes. Rapidly decaying density fluctuations can be described as propagating sound waves in the framework of linearized hydrodynamics, but there is a slow, dispersive, contribution that needs to be described by a stretched exponential over a broad range of length- and time scales - recent experiments suggest that this behavior can prevail even on micrometer length scales. The origin of this behavior is discussed in the context of fluctuations of the bilayer interface and the molecular structure of the bilayer itself. Connections to recent neutron scattering experiments are highlighted. / QC 20111014 / Modelling of biological membranes

The research focus is on the various modeling aspects of flow and sediment transport in water bodies and watersheds. The interaction of flow with a mobile bed involves a complex process in which various turbulent scales characterized by coherent structures cause a chaotic sediment motion. In many rivers and natural waterways secondary flows that are dominating flow struc-tures bring about more complications. In estuaries and open waterbodies thermal stratification and internal mixing control the flow structure besides the flow interaction with the mobile bed. To adequately model these processes 3D coupled flow and transport models are needed. The research is based on use and adaptation of open source codes for 3D hydrodynamic and sediment transport model known as Estuarine Coastal Ocean Model (ECOMSED) and the Soil and Water Assessment Tool (SWAT) model. A bed load transport model was developed and coupled to ECOMSED. The flow and sediment transport characteristics in a curved channel and a river reach were successfully captured by the model. Improvements in ECOMSED were made to study the effect of wind and basin bathymetry on mixing and flow exchange between two estuaries. Using spectral analysis the hydrological component of SWAT model was investigated for its applicability under limited data conditions in three Ethiopian catchments. / QC 20100827

Tsunamis are among the most terrifying and complex physical phenomena potentially affecting almost all coastal regions of the Earth. Tsunami waves propagate in the ocean over thousands of kilometres away from their generating source at considerable speeds. Among several other tsunamis that occurred during the past decade, the 2004 Indian Ocean Tsunami and the 2011 Tohoku Tsunami in Japan, considered to be the deadliest and costliest natural disasters in the history of mankind, respectively, have hit wide stretches of densely populated coastal areas. During these major events, severe destruction of inland structures resulted from the action of extreme hydrodynamic forces induced by tsunami flooding. Subsequent field surveys in which researchers from the University of Ottawa participated ultimately revealed that, in contrast to seismic forces, such hydrodynamic forces are not taken into proper consideration when designing buildings for tsunami prone areas. In view of these limitations, a novel interdisciplinary hydraulic-structural engineering research program was initiated at the University of Ottawa, in cooperation with the Canadian Hydraulic Centre of the National Research Council, to help develop guidelines for the sound design of nearshore structures located in such areas. The present study aims to simulate the physical laboratory experiments performed within the aforementioned research program using a single-phase three-dimensional weakly compressible Smoothed Particle Hydrodynamics (SPH) numerical model. These experiments consist in the violent impact of rapidly advancing tsunami-like hydraulic bores with individual slender structural elements. Such bores are emulated based on the classic dam-break problem. The quantitatively compared measurements include the time-history of the net base horizontal force and of the pressure distribution acting on columns of square and circular cross-sections, as well as flow characteristics such as bore-front velocity and water surface elevation. Good agreement was obtained. Results show that the magnitude and duration of the impulsive force at initial bore impact depend on the degree of entrapped air in the bore-front. The latter was found to increase considerably if the bed of the experimental flume is covered with a thin water layer of even just a few millimetres. In order to avoid large fluctuations in the pressure field and to obtain accurate simulations of the hydrodynamic forces, a Riemann solver-based formulation of the SPH method is utilized. However, this formulation induces excessive numerical diffusion, as sudden and large water surface deformations, such as splashing at initial bore impact, are less accurately reproduced. To investigate this particular issue, the small-scale physical experiment of Kleefsman et al. (2005) is also considered and modeled. Lastly, taking full advantage of the validated numerical model to better understand the underlying flow dynamics, the influence of the experimental test geometry and of the bed condition (i.e. dry vs. wet) is investigated. Numerical results show that when a bore propagates over a wet bed, its front is both deeper and steeper and it also has a lower velocity compared to when it propagates over a dry bed. These differences significantly affect the pressure distributions and resulting hydrodynamic forces acting on impacted structures.