Anisotropic Flow in Event-by-Event Ideal Hydrodynamic Simulations of root s(NN) 200 GeV Au+Au Collisions

We simulate top-energy Au + Au collisions using ideal hydrodynamics in order to make the first comparison to the complete set of midrapidity flow measurements made by the PHENIX Collaboration. A simultaneous calculation of nu(2), nu(3), nu(4), and the first event-by-event calculation of quadrangular flow defined with respect to the nu(2) event plane (nu(4){Psi(2)}) gives good agreement with measured values, including the dependence on both transverse momentum and centrality. This provides confirmation that the collision system is indeed well described as a quark-gluon plasma with an extremely small viscosity and that correlations are dominantly generated from collective effects. In addition, we present a prediction for nu(5).

Experimental model identification of open-frame underwater vehicles

Most of the works published on hydrodynamic parameter identification of open-frame underwater vehicles focus their attention almost exclusively on good coherence between simulated and measured responses, giving less importance to the determination of “actual values” for hydrodynamic parameters. To gain insight into hydrodynamic parameter experimental identification of open-frame underwater vehicles, an experimental identification procedure is proposed here to determine parameters of uncoupled and coupled models. The identification procedure includes: (i) a prior estimation of actual values of the forces/torques applied to the vehicle, (ii) identification of drag parameters from constant velocity tests and (iii) identification of inertia and coupling parameters from oscillatory tests; at this stage, the estimated values of drag parameter obtained in item (ii) are used. The procedure proposed here was used to identify the hydrodynamic parameters of LAURS—an unmanned underwater vehicle developed at the University of São Paulo. The thruster–thruster and thruster–hull interactions and the advance velocity of the vehicle are shown to have a strong impact on the efficiency of thrusters appended to open-frame underwater vehicles, especially for high advance velocities. Results of tests with excitation in 1-DOF and 3-DOF are reported and discussed, showing the feasibility of the developed procedure.

Performance Optimization of a Floating Breakwater Model Using SPH Method, with a Practical Application

Máster en Oceanografía

Novel simulation methods for Coulomb and hydrodynamic interactions

This thesis presents new methods to simulate systems with hydrodynamic and electrostatic interactions. Part 1 is devoted to computer simulations of Brownian particles with hydrodynamic interactions. The main influence of the solvent on the dynamics of Brownian particles is that it mediates hydrodynamic interactions. In the method, this is simulated by numerical solution of the Navier--Stokes equation on a lattice. To this end, the Lattice--Boltzmann method is used, namely its D3Q19 version. This model is capable to simulate compressible flow. It gives us the advantage to treat dense systems, in particular away from thermal equilibrium. The Lattice--Boltzmann equation is coupled to the particles via a friction force. In addition to this force, acting on {it point} particles, we construct another coupling force, which comes from the pressure tensor. The coupling is purely local, i.~e. the algorithm scales linearly with the total number of particles. In order to be able to map the physical properties of the Lattice--Boltzmann fluid onto a Molecular Dynamics (MD) fluid, the case of an almost incompressible flow is considered. The Fluctuation--Dissipation theorem for the hybrid coupling is analyzed, and a geometric interpretation of the friction coefficient in terms of a Stokes radius is given. Part 2 is devoted to the simulation of charged particles. We present a novel method for obtaining Coulomb interactions as the potential of mean force between charges which are dynamically coupled to a local electromagnetic field. This algorithm scales linearly, too. We focus on the Molecular Dynamics version of the method and show that it is intimately related to the Car--Parrinello approach, while being equivalent to solving Maxwell's equations with freely adjustable speed of light. The Lagrangian formulation of the coupled particles--fields system is derived. The quasi--Hamiltonian dynamics of the system is studied in great detail. For implementation on the computer, the equations of motion are discretized with respect to both space and time. The discretization of the electromagnetic fields on a lattice, as well as the interpolation of the particle charges on the lattice is given. The algorithm is as local as possible: Only nearest neighbors sites of the lattice are interacting with a charged particle. Unphysical self--energies arise as a result of the lattice interpolation of charges, and are corrected by a subtraction scheme based on the exact lattice Green's function. The method allows easy parallelization using standard domain decomposition. Some benchmarking results of the algorithm are presented and discussed.

Quantum hydrodynamic models for semiconductors with and without collisions

In dieser Arbeit werden Quantum-Hydrodynamische (QHD) Modelle betrachtet, die ihren Einsatz besonders in der Modellierung von Halbleiterbauteilen finden. Das QHD Modell besteht aus den Erhaltungsgleichungen für die Teilchendichte, das Momentum und die Energiedichte, inklusive der Quanten-Korrekturen durch das Bohmsche Potential. Zu Beginn wird eine Übersicht über die bekannten Ergebnisse der QHD Modelle unter Vernachlässigung von Kollisionseffekten gegeben, die aus ein­em Schrödinger-System für den gemischten-Zustand oder aus der Wigner-Glei­chung hergeleitet werden können. Nach der Reformulierung der eindimensionalen QHD Gleichungen mit linearem Potential als stationäre Schrö­din­ger-Gleichung werden die semianalytischen Fassungen der QHD Gleichungen für die Gleichspannungs-Kurve betrachtet. Weiterhin werden die viskosen Stabilisierungen des QHD Modells be­rück­sich­tigt, sowie die von Gardner vorgeschlagene numerische Viskosität für das {sf upwind} Finite-Differenzen Schema berechnet. Im Weiteren wird das viskose QHD Modell aus der Wigner-Glei­chung mit Fokker-Planck Kollisions-Ope­ra­tor hergeleitet. Dieses Modell enthält die physikalische Viskosität, die durch den Kollision-Operator eingeführt wird. Die Existenz der Lösungen (mit strikt positiver Teilchendichte) für das isotherme, stationäre, eindimensionale, viskose Modell für allgemeine Daten und nichthomogene Randbedingungen wird gezeigt. Die dafür notwendigen Abschätzungen hängen von der Viskosität ab und erlauben daher den Grenzübergang zum nicht-viskosen Fall nicht. Numerische Simulationen der Resonanz-Tunneldiode modelliert mit dem nichtisothermen, stationären, eindimensionalen, viskosen QHD Modell zeigen den Einfluss der Viskosität auf die Lösung. Unter Verwendung des von Degond und Ringhofer entwickelten Quanten-Entropie-Minimierungs-Verfahren werden die allgemeinen QHD-Gleichungen aus der Wigner-Boltzmann-Gleichung mit dem BGK-Kollisions-Operator hergeleitet. Die Herleitung basiert auf der vorsichtige Entwicklung des Quanten-Max­well­ians in Potenzen der skalierten Plankschen Konstante. Das so erhaltene Modell enthält auch vertex-Terme und dispersive Terme für die Ge­schwin­dig­keit. Dadurch bleibt die Gleichspannungs-Kurve für die Re­so­nanz-Tunnel­diode unter Verwendung des allgemeinen QHD Modells in einer Dimension numerisch erhalten. Die Ergebnisse zeigen, dass der dispersive Ge­schwin­dig­keits-Term die Lösung des Systems stabilisiert.

Model reduction techniques in flexible multibody dynamics with application to engine cranktrain simulation

The development of a multibody model of a motorbike engine cranktrain is presented in this work, with an emphasis on flexible component model reduction. A modelling methodology based upon the adoption of non-ideal joints at interface locations, and the inclusion of component flexibility, is developed: both are necessary tasks if one wants to capture dynamic effects which arise in lightweight, high-speed applications. With regard to the first topic, both a ball bearing model and a journal bearing model are implemented, in order to properly capture the dynamic effects of the main connections in the system: angular contact ball bearings are modelled according to a five-DOF nonlinear scheme in order to grasp the crankshaft main bearings behaviour, while an impedance-based hydrodynamic bearing model is implemented providing an enhanced operation prediction at the conrod big end locations. Concerning the second matter, flexible models of the crankshaft and the connecting rod are produced. The well-established Craig-Bampton reduction technique is adopted as a general framework to obtain reduced model representations which are suitable for the subsequent multibody analyses. A particular component mode selection procedure is implemented, based on the concept of Effective Interface Mass, allowing an assessment of the accuracy of the reduced models prior to the nonlinear simulation phase. In addition, a procedure to alleviate the effects of modal truncation, based on the Modal Truncation Augmentation approach, is developed. In order to assess the performances of the proposed modal reduction schemes, numerical tests are performed onto the crankshaft and the conrod models in both frequency and modal domains. A multibody model of the cranktrain is eventually assembled and simulated using a commercial software. Numerical results are presented, demonstrating the effectiveness of the implemented flexible model reduction techniques. The advantages over the conventional frequency-based truncation approach are discussed.

Finite Element Modeling of Shell Shape in the Freshwater Turtle Pseudemys concinna Reveals a Trade-Off between Mechanical Strength and Hydrodynamic Efficiency

Aquatic species can experience different selective pressures on morphology in different flow regimes. Species inhabiting lotic regimes often adapt to these conditions by evolving low-drag (i.e., streamlined) morphologies that reduce the likelihood of dislodgment or displacement. However, hydrodynamic factors are not the only selective pressures influencing organismal morphology and shapes well suited to flow conditions may compromise performance in other roles. We investigated the possibility of morphological trade-offs in the turtle Pseudemys concinna. Individuals living in lotic environments have flatter, more streamlined shells than those living in lentic environments; however, this flatter shape may also make the shells less capable of resisting predator-induced loads. We tested the idea that ‘‘lotic’’ shell shapes are weaker than ‘‘lentic’’ shell shapes, concomitantly examining effects of sex. Geometric morphometric data were used to transform an existing finite element shell model into a series of models corresponding to the shapes of individual turtles. Models were assigned identical material properties and loaded under identical conditions, and the stresses produced by a series of eight loads were extracted to describe the strength of the shells. ‘‘Lotic’’ shell shapes produced significantly higher stresses than ‘‘lentic’’ shell shapes, indicating that the former is weaker than the latter. Females had significantly stronger shell shapes than males, although these differences were less consistent than differences between flow regimes. We conclude that, despite the potential for many-to-one mapping of shell shape onto strength, P. concinna experiences a trade-off in shell shape between hydrodynamic and mechanical performance. This trade-off may be evident in many other turtle species or any other aquatic species that also depend on a shell for defense. However, evolution of body size may provide an avenue of escape from this trade-off in some cases, as changes in size can drastically affect mechanical performance while having little effect on hydrodynamic performance.

Hydrodynamic Characterization of Bile Salt Host-Guest Complexes by Pulsed Field Gradient NMR Spectroscopy

The work described herein is aimed at understanding primary and secondary aggregation of bile salt micelles and how micelles can perform chiral recognition of binapthyl analytes. Previous work with cholate and deoxycholate using micellar electrokinetic chromatography (MEKC) and nuclear magnetic resonance (NMR) has provided insightinto cholate and deoxycholate micelle formation, especially with respect to the critical micelle concentration (CMC). Chiral separations of the model analyte, 1,1â??-binaphthyl-2,2â??-diyl hydrogen phosphate (BNDHP), via cholate (C) and deoxycholate (DC) mediated MEKC separataions previously have shown the DC CMC to be 7-10 mM andthe cholate CMC at 14 mM at ph 12. A second model analyte,1,1â??-binaphthol (BN), was also previously investigated to probe micellar structure, but the MEKC data for this analyte implied a higher CMC, which may be interpreted as secondary aggregation. Thiswork extends the investigation of bile salts to include pulsed field gradient spin echo (PFGSE) NMR experiments being used to gain information about the size and degree of polydispersity of cholate and deoxycholate micelles. Concentrations of cholate below 10mM show a large variation in effective radius likely due to the existence of transient preliminary aggregates. The onset of the primary micelle shows a dramatic increase in effective radius of the micelle in cholate and deoxycholate. In the region of expectedsecondary aggregation a gradual increase of effective radius was observed with cholate; deoxycholate showed a persistent aggregate size in the secondary micelle region that is modulated by the presence of an analyte molecule. Effective radii of cholate anddeoxycholate (individually) were compared with and without R- and S-BNDHP in order to observe the effective radius difference of micelles with and without analyte present. The presence of S-BNDHP consistently resulted in a larger effective aggregate radius incholate and deoxycholate, confirming previous data of the S-BNDHP interacting more with the micelle than R-BNDHP. In total, various NMR techniques, like diffusion NMR can be used to gain a greater understanding of the bile salt micellization process and chiral resolution.

Energy harvesting through arterial wall deformation: design considerations for a magneto-hydrodynamic generator

As the complexity of active medical implants increases, the task of embedding a life-long power supply at the time of implantation becomes more challenging. A periodic renewal of the energy source is often required. Human energy harvesting is, therefore, seen as a possible remedy. In this paper, we present a novel idea to harvest energy from the pressure-driven deformation of an artery by the principle of magneto-hydrodynamics. The generator relies on a highly electrically conductive fluid accelerated perpendicularly to a magnetic field by means of an efficient lever arm mechanism. An artery with 10 mm inner diameter is chosen as a potential implantation site and its ability to drive the generator is established. Three analytical models are proposed to investigate the relevant design parameters and to determine the existence of an optimal configuration. The predicted output power reaches 65 μW according to the first two models and 135 μW according to the third model. It is found that the generator, designed as a circular structure encompassing the artery, should not exceed a total volume of 3 cm3.

A Model Study of the Circulation in the Pearl River Estuary (PRE) and Its Adjacent Coastal Waters: 1. Simulations and Comparison with Observations

The Princeton Ocean Model is used to study the circulation features in the Pearl River Estuary and their responses to tide, river discharge, wind, and heat flux in the winter dry and summer wet seasons. The model has an orthogonal curvilinear grid in the horizontal plane with variable spacing from 0.5 km in the estuary to 1 km on the shelf and 15 sigma levels in the vertical direction. The initial conditions and the subtidal open boundary forcing are obtained from an associated larger-scale model of the northern South China Sea. Buoyancy forcing uses the climatological monthly heat fluxes and river discharges, and both the climatological monthly wind and the realistic wind are used in the sensitivity experiments. The tidal forcing is represented by sinusoidal functions with the observed amplitudes and phases. In this paper, the simulated tide is first examined. The simulated seasonal distributions of the salinity, as well as the temporal variations of the salinity and velocity over a tidal cycle are described and then compared with the in situ survey data from July 1999 and January 2000. The model successfully reproduces the main hydrodynamic processes, such as the stratification, mixing, frontal dynamics, summer upwelling, two-layer gravitational circulation, etc., and the distributions of hydrodynamic parameters in the Pearl River Estuary and coastal waters for both the winter and the summer season.

Application of a computational model for complex fluvial ecosystems: The population dynamics of zebra mussel Dreissena polymorpha as a case study

The potential and adaptive flexibility of population dynamic P-systems (PDP) to study population dynamics suggests that they may be suitable for modelling complex fluvial ecosystems, characterized by a composition of dynamic habitats with many variables that interact simultaneously. Using as a model a reservoir occupied by the zebra mussel Dreissena polymorpha, we designed a computational model based on P systems to study the population dynamics of larvae, in order to evaluate management actions to control or eradicate this invasive species. The population dynamics of this species was simulated under different scenarios ranging from the absence of water flow change to a weekly variation with different flow rates, to the actual hydrodynamic situation of an intermediate flow rate. Our results show that PDP models can be very useful tools to model complex, partially desynchronized, processes that work in parallel. This allows the study of complex hydroecological processes such as the one presented, where reproductive cycles, temperature and water dynamics are involved in the desynchronization of the population dynamics both, within areas and among them. The results obtained may be useful in the management of other reservoirs with similar hydrodynamic situations in which the presence of this invasive species has been documented.

Gravitational wave background from Standard Model physics: qualitative features

Because of physical processes ranging from microscopic particle collisions to macroscopic hydrodynamic fluctuations, any plasma in thermal equilibrium emits gravitational waves. For the largest wavelengths the emission rate is proportional to the shear viscosity of the plasma. In the Standard Model at 0T > 16 GeV, the shear viscosity is dominated by the most weakly interacting particles, right-handed leptons, and is relatively large. We estimate the order of magnitude of the corresponding spectrum of gravitational waves. Even though at small frequencies (corresponding to the sub-Hz range relevant for planned observatories such as eLISA) this background is tiny compared with that from non-equilibrium sources, the total energy carried by the high-frequency part of the spectrum is non-negligible if the production continues for a long time. We suggest that this may constrain (weakly) the highest temperature of the radiation epoch. Observing the high-frequency part directly sets a very ambitious goal for future generations of GHz-range detectors.

Comparison of 3D kinetic and hydrodynamic models to ROSINA-COPS measurements of the neutral coma of 67P/Churyumov-Gerasimenko

67P/Churyumov-Gerasimenko (67P) is a Jupiter-family comet and the object of investigation of the European Space Agency mission Rosetta. This report presents the first full 3D simulation results of 67P’s neutral gas coma. In this study we include results from a direct simulation Monte Carlo method, a hydrodynamic code, and a purely geometric calculation which computes the total illuminated surface area on the nucleus. All models include the triangulated 3D shape model of 67P as well as realistic illumination and shadowing conditions. The basic concept is the assumption that these illumination conditions on the nucleus are the main driver for the gas activity of the comet. As a consequence, the total production rate of 67P varies as a function of solar insolation. The best agreement between the model and the data is achieved when gas fluxes on the night side are in the range of 7% to 10% of the maximum flux, accounting for contributions from the most volatile components. To validate the output of our numerical simulations we compare the results of all three models to in situ gas number density measurements from the ROSINA COPS instrument. We are able to reproduce the overall features of these local neutral number density measurements of ROSINA COPS for the time period between early August 2014 and January 1 2015 with all three models. Some details in the measurements are not reproduced and warrant further investigation and refinement of the models. However, the overall assumption that illumination conditions on the nucleus are at least an important driver of the gas activity is validated by the models. According to our simulation results we find the total production rate of 67P to be constant between August and November 2014 with a value of about 1 × 10²⁶ molecules s⁻¹.

Validation of a Spatially Continuous EDEN Water-Surface Model for the Everglades, Florida

The Everglades Depth Estimation Network (EDEN) is an integrated network of realtime water-level monitoring, ground-elevation modeling, and water-surface modeling that provides scientists and managers with current (2000-present), online water-stage and water-depth information for the entire freshwater portion of the Greater Everglades. Continuous daily spatial interpolations of the EDEN network stage data are presented on grid with 400-square-meter spacing. EDEN offers a consistent and documented dataset that can be used by scientists and managers to: (1) guide large-scale field operations, (2) integrate hydrologic and ecological responses, and (3) support biological and ecological assessments that measure ecosystem responses to the implementation of the Comprehensive Everglades Restoration Plan (CERP) (U.S. Army Corps of Engineers, 1999). The target users are biologists and ecologists examining trophic level responses to hydrodynamic changes in the Everglades. The first objective of this report is to validate the spatially continuous EDEN water-surface model for the Everglades, Florida developed by Pearlstine et al. (2007) by using an independent field-measured data-set. The second objective is to demonstrate two applications of the EDEN water-surface model: to estimate site-specific ground elevation by using the validated EDEN water-surface model and observed water depth data; and to create water-depth hydrographs for tree islands. We found that there are no statistically significant differences between model-predicted and field-observed water-stage data in both southern Water Conservation Area (WCA) 3A and WCA 3B. Tree island elevations were derived by subtracting field water-depth measurements from the predicted EDEN water-surface. Water-depth hydrographs were then computed by subtracting tree island elevations from the EDEN water stage. Overall, the model is reliable by a root mean square error (RMSE) of 3.31 cm. By region, the RMSE is 2.49 cm and 7.77 cm in WCA 3A and 3B, respectively. This new landscape-scale hydrological model has wide applications for ongoing research and management efforts that are vital to restoration of the Florida Everglades. The accurate, high-resolution hydrological data, generated over broad spatial and temporal scales by the EDEN model, provides a previously missing key to understanding the habitat requirements and linkages among native and invasive populations, including fish, wildlife, wading birds, and plants. The EDEN model is a powerful tool that could be adapted for other ecosystem-scale restoration and management programs worldwide.

(Table 5) Experimental hydrodynamic parameters for different soil types obtained in the Lüchow-Dannenberg district, Lower Saxony

In three typical sandy soils of Northern Germany the mobility of radioactive fission products of technetium, iodine, ruthenium and zirconium have been investigated in dependence of the hydrodynamic and physico-chemical soil properties. The laboratory experiments, which simulated fall-out events, used soil columns (1 m length, 30 cm diameter) taken as undisturbed as possible. By measurements of the breakthrough curves in the percolate and of the depth distribution of radionuclides in the soil columns after 6 months the average transport velocity could be determined. These values could be compared with the average water velocity measured by 3H tagging. Three qualitative mobility relations were observed: Ranker: Tc > Ru > I > Zr; Podsol: Tc > Ru > I > Zr; Brown forest soil: Tc = Ru > I > Zr. Relations between some physico-chemical soil properties and the retardation of radionuclides due to adsorption could be observed (eg. retardation of iodine and technetium by organic substances). The average retardation factors of the radionuclides and the hydrodynamic soil parameters are used in a model which gives a quantitative assessment of the hazard of groundwater contamination by a fall-out event in areas covered with comparable soils.