Thermal fluctuations and boundary conditions in the lattice Boltzmann method

The lattice Boltzmann method is a popular approach for simulating hydrodynamic interactions in soft matter and complex fluids. The solvent is represented on a discrete lattice whose nodes are populated by particle distributions that propagate on the discrete links between the nodes and undergo local collisions. On large length and time scales, the microdynamics leads to a hydrodynamic flow field that satisfies the Navier-Stokes equation. In this thesis, several extensions to the lattice Boltzmann method are developed. In complex fluids, for example suspensions, Brownian motion of the solutes is of paramount importance. However, it can not be simulated with the original lattice Boltzmann method because the dynamics is completely deterministic. It is possible, though, to introduce thermal fluctuations in order to reproduce the equations of fluctuating hydrodynamics. In this work, a generalized lattice gas model is used to systematically derive the fluctuating lattice Boltzmann equation from statistical mechanics principles. The stochastic part of the dynamics is interpreted as a Monte Carlo process, which is then required to satisfy the condition of detailed balance. This leads to an expression for the thermal fluctuations which implies that it is essential to thermalize all degrees of freedom of the system, including the kinetic modes. The new formalism guarantees that the fluctuating lattice Boltzmann equation is simultaneously consistent with both fluctuating hydrodynamics and statistical mechanics. This establishes a foundation for future extensions, such as the treatment of multi-phase and thermal flows. An important range of applications for the lattice Boltzmann method is formed by microfluidics. Fostered by the "lab-on-a-chip" paradigm, there is an increasing need for computer simulations which are able to complement the achievements of theory and experiment. Microfluidic systems are characterized by a large surface-to-volume ratio and, therefore, boundary conditions are of special relevance. On the microscale, the standard no-slip boundary condition used in hydrodynamics has to be replaced by a slip boundary condition. In this work, a boundary condition for lattice Boltzmann is constructed that allows the slip length to be tuned by a single model parameter. Furthermore, a conceptually new approach for constructing boundary conditions is explored, where the reduced symmetry at the boundary is explicitly incorporated into the lattice model. The lattice Boltzmann method is systematically extended to the reduced symmetry model. In the case of a Poiseuille flow in a plane channel, it is shown that a special choice of the collision operator is required to reproduce the correct flow profile. This systematic approach sheds light on the consequences of the reduced symmetry at the boundary and leads to a deeper understanding of boundary conditions in the lattice Boltzmann method. This can help to develop improved boundary conditions that lead to more accurate simulation results.

Hydrodynamic induced by an array of wave energy converters. Experimental and numerical analysis

The thesis analyses the hydrodynamic induced by an array of Wave energy Converters (WECs), under an experimental and numerical point of view. WECs can be considered an innovative solution able to contribute to the green energy supply and –at the same time– to protect the rear coastal area under marine spatial planning considerations. This research activity essentially rises due to this combined concept. The WEC under exam is a floating device belonging to the Wave Activated Bodies (WAB) class. Experimental data were performed at Aalborg University in different scales and layouts, and the performance of the models was analysed under a variety of irregular wave attacks. The numerical simulations performed with the codes MIKE 21 BW and ANSYS-AQWA. Experimental results were also used to calibrate the numerical parameters and/or to directly been compared to numerical results, in order to extend the experimental database. Results of the research activity are summarized in terms of device performance and guidelines for a future wave farm installation. The device length should be “tuned” based on the local climate conditions. The wave transmission behind the devices is pretty high, suggesting that the tested layout should be considered as a module of a wave farm installation. Indications on the minimum inter-distance among the devices are provided. Furthermore, a CALM mooring system leads to lower wave transmission and also larger power production than a spread mooring. The two numerical codes have different potentialities. The hydrodynamics around single and multiple devices is obtained with MIKE 21 BW, while wave loads and motions for a single moored device are derived from ANSYS-AQWA. Combining the experimental and numerical it is suggested –for both coastal protection and energy production– to adopt a staggered layout, which will maximise the devices density and minimize the marine space required for the installation.

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.

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⁻¹.

M-AARE - Coupling atmospheric, hydrological, hydrodynamic and damage models in the Aare river basin, Switzerland

The triggering mechanism and the temporal evolution of large flood events, especially of worst-case scenarios, are not yet fully understood. Consequently, the cumulative losses of extreme floods are unknown. To study the link between weather conditions, discharges and flood losses it is necessary to couple atmospheric, hydrological, hydrodynamic and damage models. The objective of the M-AARE project is to test the potentials and opportunities of a model chain that relates atmospheric conditions to flood losses or risks. The M-AARE model chain is a set of coupled models consisting of four main components: the precipitation module, the hydrology module, the hydrodynamic module, and the damage module. The models are coupled in a cascading framework with harmonized time-steps. First exploratory applications show that the one way coupling of the WRF-PREVAH-BASEMENT models has been achieved and provides promising new insights for a better understanding of key aspects in flood risk analysis.

Hydrodynamic modeling of mineral wool fiber suspensions in a two-dimensional flow

A consequence of a loss of coolant accident is that the local insulation material is damaged and maybe transported to the containment sump where it can penetrate and/or block the sump strainers. An experimental and theoretical study, which examines the transport of mineral wool fibers via single and multi-effect experiments is being performed. This paper focuses on the experiments and simulations performed for validation of numerical models of sedimentation and resuspension of mineral wool fiber agglomerates in a racetrack type channel. Three velocity conditions are used to test the response of two dispersed phase fiber agglomerates to two drag correlations and to two turbulent dispersion coefficients. The Eulerian multiphase flow model is applied with either one or two dispersed phases.

The hydrodynamic behaviour and mass transfer characteristics of single droplets in a pulsed sieve plate column

The literature relating to the performance of pulsed sieve plate liquid-liquid extraction columns and the relevant hydrodynamic phenomenon have been surveyed. Hydrodynamic behaviour and mass transfer characteristics of droplets in turbulent and non-turbulent conditions have also been reviewed. Hydrodynamic behaviour, i.e. terminal and characteristic velocity of droplets, droplet size and droplet breakup processes, and mass transfer characteristics of single droplets (d≤0.6 cm) were investigated under pulsed (mixer-settler & transitional regimes) and non-pulsed conditions in a 5.0 cm diameter, 100 cm high, pulsed sieve plate column with three different sieve plate types and variable plate spacing. The system used was toluene (displaced) - acetone - distilled water. Existing photographic techniques for following and recording the droplet behaviour, and for observing the parameters of the pulse and the pulse shape were further developed and improved. A unique illumination technique was developed by which a moving droplet could be photographed using cine or video photography with good contrast without using any dye. Droplet size from a given nozzle and droplet velocity for a given droplet diameter are reduced under pulsing condition, and it was noted that this effect is enhanced in the presence of sieve plate. The droplet breakup processes are well explained by reference to an impact-breakup mechanism. New correlations to predict droplet diameter based on this mechanism are given below.vskip 1.0cm or in dimensionless groups as follows:- (We)crit= 3.12 - 1.79 (Eo)crit A correlation based on the isotropic turbulence theory was developed to calculate droplet diameter in the emulsion regime.vskip 1.0cm Experimental results show that in the mixer-settler and transitional regimes, pulsing parameters had little effect on the overall dispersed phase mass transfer coefficient during the droplet formation and unhindered travel periods.

Non-Linear oil firm force coefficients for a journal bearing operating under aligned and misaligned conditions

It is well established that hydrodynamic journal bearings are responsible for self-excited vibrations and have the effect of lowering the critical speeds of rotor systems. The forces within the oil film wedge, generated by the vibrating journal, may be represented by displacement and velocity coefficient~ thus allowing the dynamical behaviour of the rotor to be analysed both for stability purposes and for anticipating the response to unbalance. However, information describing these coefficients is sparse, misleading, and very often not applicable to industrial type bearings. Results of a combined analytical and experimental investigation into the hydrodynamic oil film coefficients operating in the laminar region are therefore presented, the analysis being applied to a 120 degree partial journal bearing having a 5.0 in diameter journal and a LID ratio of 1.0. The theoretical analysis shows that for this type of popular bearing, the eight linearized coefficients do not accurately describe the behaviour of the vibrating journal based on the theory of small perturbations, due to them being masked by the presence of nonlinearity. A method is developed using the second order terms of Taylor expansion whereby design charts are provided which predict the twentyeight force coefficients for both aligned, and for varying amounts of journal misalignment. The resulting non-linear equations of motion are solved using a modified Newton-Raphson method whereby the whirl trajectories are obtained, thus providing a physical appreciation of the bearing characteristics under dynamically loaded conditions.

Concurrent multiscale modelling of atomistic and hydrodynamic processes in liquids

Fluctuations of liquids at the scales where the hydrodynamic and atomistic descriptions overlap are considered. The importance of these fluctuations for atomistic motions is discussed and examples of their accurate modelling with a multi-space-time-scale fluctuating hydrodynamics scheme are provided. To resolve microscopic details of liquid systems, including biomolecular solutions, together with macroscopic fluctuations in space-time, a novel hybrid atomistic-fluctuating hydrodynamics approach is introduced. For a smooth transition between the atomistic and continuum representations, an analogy with two-phase hydrodynamics is used that leads to a strict preservation of macroscopic mass and momentum conservation laws. Examples of numerical implementation of the new hybrid approach for the multiscale simulation of liquid argon in equilibrium conditions are provided. © 2014 The Author(s) Published by the Royal Society.

A CFD study of biomass pyrolysis in a downer reactor equipped with a novel gas-solid separator - I:hydrodynamic performance

This study presents the first part of a CFD study on the performance of a downer reactor for biomass pyrolysis. The reactor was equipped with a novel gas-solid separation method, developed by the co-authors from the ICFAR (Canada). The separator, which was designed to allow for fast separation of clean pyrolysis gas, consisted of a cone deflector and a gas exit pipe installed inside the downer reactor. A multi-fluid model (Eulerian-Eulerian) with constitutive relations adopted from the kinetic theory of granular flow was used to simulate the multiphase flow. The effects of the various parameters including operation conditions, separator geometry and particle properties on the overall hydrodynamics and separation efficiency were investigated. The model prediction of the separator efficiency was compared with experimental measurements. The results revealed distinct hydrodynamic features around the cone separator, allowing for up to 100% separation efficiency. The developed model provided a platform for the second part of the study, where the biomass pyrolysis is simulated and the product quality as a function of operating conditions is analyzed. Crown Copyright © 2014 Published by Elsevier B.V. All rights reserved.

Identification of a hydrodynamic threshold in karst rocks from the Biscayne Aquifer, south Florida, USA

A hydrodynamic threshold between Darcian and non-Darcian flow conditions was found to occur in cubes of Key Largo Limestone from Florida, USA (one cube measuring 0.2 m on each side, the other 0.3 m) at an effective porosity of 33% and a hydraulic conductivity of 10 m/day. Below these values, flow was laminar and could be described as Darcian. Above these values, hydraulic conductivity increased greatly and flow was non-laminar. Reynolds numbers (Re) for these experiments ranged from

The effect of a reduction of nitrogen and phosphorus to eutrophic and mesotrophic conditions towards silver nanoparticle (AgNP) on freshwater green algae species

We investigated the sensitivity of algae towards silver nanoparticles with OECD test medium and lower nutrient concentrations under standard test conditions to improve comparability and to exclude any other confounding factor aside nutrient levels. Two unicellular freshwater microalgae Desmodesmus subspicatus and Raphidocelis subcapitata were chosen due to their status as standard test organisms for the algae growth inhibition test and the response to changes in nutrient supply was compared. The original medium was used as the reference (standard). For the other four media, the amount of either nitrogen or phosphorus in the medium was lowered from half (50%) to one-fourth (25 %) of that of the OECD guideline, resulting in the following media: 50% N, 25% N, 50% P and 25% P medium. As test substance, the OECD reference material NM-300K was used. For this reason, the characterization of AgNP was done using DLS and Absorption spectra (UV/vis). Actual silver concentrations and ionic silver concentrations were measured at the highest test concentration used (100 µg Ag L-1) in R. subcapitata treatments only to reduce the number of samples. All tests were run according to the OECD guideline 201 with sterilized 50 mL cell culture flask. Each medium was tested using the test conditions for culturing with 3 replicates. Test concentrations for both algae species were 0, 25, 50 and 100 µg Ag L-1 for OECD, 50% P and 25% P while for both N reductions, the silver concentrations were 0, 10, 25 and 100 µg Ag L-1. Samples for determining the algal density were taken at every 24 h.

A Non-catalytic Deep Desulphurization Process using Hydrodynamic Cavitation

A novel approach is developed for desulphurization of fuels or organics without use of catalyst. In this process, organic and aqueous phases are mixed in a predefined manner under ambient conditions and passed through a cavitating device. Vapor cavities formed in the cavitating device are then collapsed which generate (in-situ) oxidizing species which react with the sulphur moiety resulting in the removal of sulphur from the organic phase. In this work, vortex diode was used as a cavitating device. Three organic solvents (n-octane, toluene and n-octanol) containing known amount of a model sulphur compound (thiophene) up to initial concentrations of 500 ppm were used to verify the proposed method. A very high removal of sulphur content to the extent of 100% was demonstrated. The nature of organic phase and the ratio of aqueous to organic phase were found to be the most important process parameters. The results were also verified and substantiated using commercial diesel as a solvent. The developed process has great potential for deep of various organics, in general, and for transportation fuels, in particular.