Computersimulationen zur Dynamik von verdünnten und halbverdünnten Polymerlösungen

Diese Arbeit legt eine neue Methode zur Simulation derDynamik vonPolymeren in verdünnter und halbverdünnterLösung vor. Die Effizienz der Methode und derAnstieg der Computerleistung in den letzten Jahren erlaubenes, weitaus komplexere Systeme als bisher zu betrachten.Die neue Methode modelliert die Polymere als Kugel-Feder-Ketten, die mittels Molekulardynamik simuliertwerden. Die Flüssigkeit wird durch die numerischeLösung der Kontinuitätsgleichungund der Navier-Stokes-Gleichung mit derLattice-Boltzmann-Methodemodelliert. Die Flüssigkeit wird über eineReibungskraft an die Monomere des Kugel-Feder-Modellsgekoppelt. Die Methode wird auf das Problem einer flexiblen EinzelketteimLösungsmittel angewendet. Der Vergleich derErgebnisse mit einer existierenden reinenMolekulardynamik-Simulationergibt Übereinstimmung innerhalb weniger Prozent,während die neueMethode um etwa einen Faktor 20 weniger CPU-Zeitbenötigt. Eine semiflexible Kette zeigt völliganderes Verhalten: Die Hydrodynamik spielt im Gegensatz zur flexiblen Ketteeineuntergeordnete Rolle. Simulationen von halbverdünntenLösungen flexibler Kettenbestehend aus insgesamt 50000 Monomeren zeigen zum erstenMal direkt dieAbschirmung sowohl der Volumenausschluss-Wechselwirkung alsauch derHydrodynamik.

Modellierung von Porenraumgeometrien und Transport in korngestützten porösen Medien

Poröse Medien spielen in der Hydrosphäre eine wesentliche Rolle bei der Strömung und beim Transport von Stoffen. In diesem Raum finden komplexe Prozesse statt: Advektion, Kon-vektion, Diffusion, hydromechanische Dispersion, Sorption, Komplexierung, Ionenaustausch und Abbau. Die strömungsmechanischen- und die Transportverhältnisse in porösen Medien werden direkt durch die Geometrie des Porenraumes selbst und durch die Eigenschaften der transportierten (oder strömenden) Medien bestimmt. In der Praxis wird eine Vielzahl von empirischen Modellen verwendet, die die Eigenschaften des porösen Mediums in repräsentativen Elementarvolumen wiedergeben. Die Ermittlung der in empirischen Modellen verwendeten Materialparameter erfolgt über Labor- oder Feldbestimmungsmethoden. Im Rahmen dieser Arbeit wurde das Computer-modell PoreFlow entwickelt, welches die hydraulischen Eigenschaften eines korngestützten porösen Mediums aus der mikroskopischen Modellierung des Fluidflusses und Transportes ableitet. Das poröse Modellmedium wird durch ein dreidimensionales Kugelpackungsmodell, zusam-mengesetzt aus einer beliebigen Kornverteilung, dargestellt. Im Modellporenraum wird die Strömung eines Fluids basierend auf einer stationären Lösung der Navier-Stokes-Gleichung simuliert. Die Ergebnisse der Modellsimulationen an verschiedenen Modellmedien werden mit den Ergebnissen von Säulenversuchen verglichen. Es zeigt sich eine deutliche Abhängigkeit der Strömungs- und Transportparameter von der Porenraumgeometrie sowohl in den Modell-simulationen als auch in den Säulenexperimenten.

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.

Scherverhalten kolloidaler Lagenphasen in Rohrströmungen

Kolloidale Suspensionen aus identischen kugelförmigen, geladenen Partikeln in wässrigen Medien stellen ein ideales Modellsystem zur Untersuchung des Gleichgewichtsverhaltens, aber auch des Nicht-Gleichgewichtsverhaltens Weicher Materie dar. So bilden derartige Systeme bei hinreichend starker und langreichweitiger elektrostatischer Repulsion fluid und kristallin geordnete Strukturen aus, die wegen der weitreichenden Analogie zu atomar kondensierter Materie als kolloidale Fluide und Kristalle bezeichnet werden. Von großem Vorteil ist dabei die Möglichkeit zur kontrollierten Einstellung der Wechselwirkung und die gute optische Zugänglichkeit für Mikroskopie und Lichtstreuung sowie die Weichheit der Materialien, aufgrund derer sich auch Zustände fernab des mechanischen Gleichgewichts gezielt präparieren lassen. Themenstellung der vorliegenden Arbeit ist die Untersuchung des Phasenverhaltens und der Fließmechanismen kolloidaler Kristalle in einer Rohrströmung. Im ersten Teil der Arbeit wird gezeigt, dass beim Fluss durch eine zylindrische Röhre Mehrphasenkoexistenz auftritt, wobei ein polykristalliner Kern von einer isotropen Scherschmelze umgeben ist. Zusätzlich treten an der Grenze zwischen diesen Phasen und an der Rohrwand Phasen hexagonal geordneter übereinander hinweggleitender Lagen auf. Der Vergleich zwischen auf der Basis der Navier-Stokes-Gleichung theoretisch berechneten und gemessenen Geschwindigkeitsprofilen zeigt, dass jede dieser Phasen für sich Newtonsches Fließverhalten aufweist. Die Gesamtviskosität ist hingegen durch die mit dem Durchsatz veränderliche Phasenzusammensetzung Nicht-Newtonsch. Damit gelang es, die erstmalig von Würth beschriebene Scherverdünnung auf eine Veränderung der Phasenzusammensetzung zurückzuführen. Im zweiten Teil der Arbeit wurde erstmals das Fließverhalten der Lagenphasen mittels Lichtstreuung und Korrelationsanalyse untersucht. Dafür wurde ein im Prinzip einfacher, aber leistungsstarker Aufbau realisiert, der es erlaubt, die zeitliche Veränderung der Bragg-Reflexe der Lagenphase in radialer und azimutaler Richtung zu verfolgen und mittels Fourieranalyse zu analysieren. In Abhängigkeit vom Durchsatz geht die zunächst rastend gleitende Lagenphase in eine frei gleitende Lagenphase über, wobei charakteristische Veränderungen der Spektren sowie der Korrelationsfunktionen auftreten, die detailliert diskutiert werden. Der Übergang im Gleitmechanismus ist mit einem Verlust der Autokorrelation der Rotationskomponente der periodischen Intra-Lagenverzerrung verbunden, während die Kompressionskomponente erhalten bleibt. Bei hohen Durchflüssen lassen die Reflexbewegungen auf das Auftreten einer Eigenschwingung der frei gleitenden Lagen schließen. Diese Schwingung lässt sich als Rotationsbewegung, gekoppelt mit einer transversalen Auslenkung in Vortexrichtung, beschreiben. Die Ergebnisse erlauben eine detaillierte Diskussion von verschiedenen Modellvorstellungen anderer Autoren.

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.

Numerical modelling of fluid-structure interaction with application in hemodynamics

The thesis deals with numerical algorithms for fluid-structure interaction problems with application in blood flow modelling. It starts with a short introduction on the mathematical description of incompressible viscous flow with non-Newtonian viscosity and a moving linear viscoelastic structure. The mathematical model consists of the generalized Navier-Stokes equation used for the description of fluid flow and the generalized string model for structure movement. The arbitrary Lagrangian-Eulerian approach is used in order to take into account moving computational domain. A part of the thesis is devoted to the discussion on the non-Newtonian behaviour of shear-thinning fluids, which is in our case blood, and derivation of two non-Newtonian models frequently used in the blood flow modelling. Further we give a brief overview on recent fluid-structure interaction schemes with discussion about the difficulties arising in numerical modelling of blood flow. Our main contribution lies in numerical and experimental study of a new loosely-coupled partitioned scheme called the kinematic splitting fluid-structure interaction algorithm. We present stability analysis for a coupled problem of non-Newtonian shear-dependent fluids in moving domains with viscoelastic boundaries. Here, we assume both, the nonlinearity in convective as well is diffusive term. We analyse the convergence of proposed numerical scheme for a simplified fluid model of the Oseen type. Moreover, we present series of experiments including numerical error analysis, comparison of hemodynamic parameters for the Newtonian and non-Newtonian fluids and comparison of several physiologically relevant computational geometries in terms of wall displacement and wall shear stress. Numerical analysis and extensive experimental study for several standard geometries confirm reliability and accuracy of the proposed kinematic splitting scheme in order to approximate fluid-structure interaction problems.

Implementation of monotonic higher order upwind scheme in KIVA 4

KIVA is a FORTRAN code developed by Los Alamos national lab to simulate complete engine cycle. KIVA is a flow solver code which is used to perform calculation of properties in a fluid flow field. It involves using various numerical schemes and methods to solve the Navier-Stokes equation. This project involves improving the accuracy of one such scheme by upgrading it to a higher order scheme. The numerical scheme to be modified is used in the critical final stage calculation called as rezoning phase. The primitive objective of this project is to implement a higher order numerical scheme, to validate and verify that the new scheme is better than the existing scheme. The latest version of the KIVA family (KIVA 4) is used for implementing the higher order scheme to support handling the unstructured mesh. The code is validated using the traditional shock tube problem and the results are verified to be more accurate than the existing schemes in reference with the analytical result. The convection test is performed to compare the computational accuracy on convective transfer; it is found that the new scheme has less numerical diffusion compared to the existing schemes. A four valve pentroof engine, an example case of KIVA package is used as application to ensure the stability of the scheme in practical application. The results are compared for the temperature profile. In spite of all the positive results, the numerical scheme implemented has a downside of consuming more CPU time for the computational analysis. The detailed comparison is provided. However, in an overview, the implementation of the higher order scheme in the latest code KIVA 4 is verified to be successful and it gives better results than the existing scheme which satisfies the objective of this project.

Noninvasive pressure difference mapping derived from 4D flow MRI in patients with unrepaired and repaired aortic coarctation.

PURPOSE To develop a method for computing and visualizing pressure differences derived from time-resolved velocity-encoded three-dimensional phase-contrast magnetic resonance imaging (4D flow MRI) and to compare pressure difference maps of patients with unrepaired and repaired aortic coarctation to young healthy volunteers. METHODS 4D flow MRI data of four patients with aortic coarctation either before or after repair (mean age 17 years, age range 3-28, one female, three males) and four young healthy volunteers without history of cardiovascular disease (mean age 24 years, age range 20-27, one female, three males) was acquired using a 1.5-T clinical MR scanner. Image analysis was performed with in-house developed image processing software. Relative pressures were computed based on the Navier-Stokes equation. RESULTS A standardized method for intuitive visualization of pressure difference maps was developed and successfully applied to all included patients and volunteers. Young healthy volunteers exhibited smooth and regular distribution of relative pressures in the thoracic aorta at mid systole with very similar distribution in all analyzed volunteers. Patients demonstrated disturbed pressures compared to volunteers. Changes included a pressure drop at the aortic isthmus in all patients, increased relative pressures in the aortic arch in patients with residual narrowing after repair, and increased relative pressures in the descending aorta in a patient after patch aortoplasty. CONCLUSIONS Pressure difference maps derived from 4D flow MRI can depict alterations of spatial pressure distribution in patients with repaired and unrepaired aortic coarctation. The technique might allow identifying pathophysiological conditions underlying complications after aortic coarctation repair.

Noninvasive 4D pressure difference mapping derived from 4D flow MRI in patients with repaired aortic coarctation: comparison with young healthy volunteers

To assess spatial and temporal pressure characteristics in patients with repaired aortic coarctation compared to young healthy volunteers using time-resolved velocity-encoded three-dimensional phase-contrast magnetic resonance imaging (4D flow MRI) and derived 4D pressure difference maps. After in vitro validation against invasive catheterization as gold standard, 4D flow MRI of the thoracic aorta was performed at 1.5T in 13 consecutive patients after aortic coarctation repair without recoarctation and 13 healthy volunteers. Using in-house developed processing software, 4D pressure difference maps were computed based on the Navier-Stokes equation. Pressure difference amplitudes, maximum slope of pressure amplitudes and spatial pressure range at mid systole were retrospectively measured by three readers, and twice by one reader to assess inter- and intraobserver agreement. In vitro, pressure differences derived from 4D flow MRI showed excellent agreement to invasive catheter measurements. In vivo, pressure difference amplitudes, maximum slope of pressure difference amplitudes and spatial pressure range at mid systole were significantly increased in patients compared to volunteers in the aortic arch, the proximal descending and the distal descending thoracic aorta (p < 0.05). Greatest differences occurred in the proximal descending aorta with values of the three parameters for patients versus volunteers being 19.7 ± 7.5 versus 10.0 ± 2.0 (p < 0.001), 10.9 ± 10.4 versus 1.9 ± 0.4 (p = 0.002), and 8.7 ± 6.3 versus 1.6 ± 0.9 (p < 0.001). Inter- and intraobserver agreements were excellent (p < 0.001). Noninvasive 4D pressure difference mapping derived from 4D flow MRI enables detection of altered intraluminal aortic pressures and showed significant spatial and temporal changes in patients with repaired aortic coarctation.

Numerical simulation of magnetite segregation in a dense medium cyclone

Numerical simulations of turbulent driven flow in a dense medium cyclone with magnetite medium have been conducted using Fluent. The predicted air core shape and diameter were found to be close to the experimental results measured by gamma ray tomography. It is possible that the Large eddy simulation (LES) turbulence model with Mixture multi-phase model can be used to predict the air/slurry interface accurately although the LES may need a finer grid. Multi-phase simulations (air/water/medium) are showing appropriate medium segregation effects but are over-predicting the level of segregation compared to that measured by gamma-ray tomography in particular with over prediction of medium concentrations near the wall. Further, investigated the accurate prediction of axial segregation of magnetite using the LES turbulence model together with the multi-phase mixture model and viscosity corrections according to the feed particle loading factor. Addition of lift forces and viscosity correction improved the predictions especially near the wall. Predicted density profiles are very close to gamma ray tomography data showing a clear density drop near the wall. The effect of size distribution of the magnetite has been fully studied. It is interesting to note that the ultra-fine magnetite sizes (i.e. 2 and 7 mu m) are distributed uniformly throughout the cyclone. As the size of magnetite increases, more segregation of magnetite occurs close to the wall. The cut-density (d(50)) of the magnetite segregation is 32 gm, which is expected with superfine magnetite feed size distribution. At higher feed densities the agreement between the [Dungilson, 1999; Wood, J.C., 1990. A performance model for coal-washing dense medium cyclones, Ph.D. Thesis, JKMRC, University of Queensland] correlations and the CFD are reasonably good, but the overflow density is lower than the model predictions. It is believed that the excessive underflow volumetric flow rates are responsible for under prediction of the overflow density. (c) 2006 Elsevier Ltd. All rights reserved.

Study of flow behaviour in a three-product cyclone using computational fluid dynamics

Simplicity in design and minimal floor space requirements render the hydrocyclone the preferred classifier in mineral processing plants. Empirical models have been developed for design and process optimisation but due to the complexity of the flow behaviour in the hydrocyclone these do not provide information on the internal separation mechanisms. To study the interaction of design variables, the flow behaviour needs to be considered, especially when modelling the new three-product cyclone. Computational fluid dynamics (CFD) was used to model the three-product cyclone, in particular the influence of the dual vortex finder arrangement on flow behaviour. From experimental work performed on the UG2 platinum ore, significant differences in the classification performance of the three-product cyclone were noticed with variations in the inner vortex finder length. Because of this simulations were performed for a range of inner vortex finder lengths. Simulations were also conducted on a conventional hydrocyclone of the same size to enable a direct comparison of the flow behaviour between the two cyclone designs. Significantly, high velocities were observed for the three-product cyclone with an inner vortex finder extended deep into the conical section of the cyclone. CFD studies revealed that in the three-product cyclone, a cylindrical shaped air-core is observed similar to conventional hydrocyclones. A constant diameter air-core was observed throughout the inner vortex finder length, while no air-core was present in the annulus. (c) 2006 Elsevier Ltd. All rights reserved.

Reconstruction of an unsteady flow from incomplete boundary data:Part I. Theory

Purpose – To propose and investigate a stable numerical procedure for the reconstruction of the velocity of a viscous incompressible fluid flow in linear hydrodynamics from knowledge of the velocity and fluid stress force given on a part of the boundary of a bounded domain. Design/methodology/approach – Earlier works have involved the similar problem but for stationary case (time-independent fluid flow). Extending these ideas a procedure is proposed and investigated also for the time-dependent case. Findings – The paper finds a novel variation method for the Cauchy problem. It proves convergence and also proposes a new boundary element method. Research limitations/implications – The fluid flow domain is limited to annular domains; this restriction can be removed undertaking analyses in appropriate weighted spaces to incorporate singularities that can occur on general bounded domains. Future work involves numerical investigations and also to consider Oseen type flow. A challenging problem is to consider non-linear Navier-Stokes equation. Practical implications – Fluid flow problems where data are known only on a part of the boundary occur in a range of engineering situations such as colloidal suspension and swimming of microorganisms. For example, the solution domain can be the region between to spheres where only the outer sphere is accessible for measurements. Originality/value – A novel variational method for the Cauchy problem is proposed which preserves the unsteady Stokes operator, convergence is proved and using recent for the fundamental solution for unsteady Stokes system, a new boundary element method for this system is also proposed.

Extended self-similarity works for the Burgers equation and why

Extended self-similarity (ESS), a procedure that remarkably extends the range of scaling for structure functions in Navier-Stokes turbulence and thus allows improved determination of intermittency exponents, has never been fully explained. We show that ESS applies to Burgers turbulence at high Reynolds numbers and we give the theoretical explanation of the numerically observed improved scaling at both the IR and UV end, in total a gain of about three quarters of a decade: there is a reduction of subdominant contributions to scaling when going from the standard structure function representation to the ESS representation. We conjecture that a similar situation holds for three-dimensional incompressible turbulence and suggest ways of capturing subdominant contributions to scaling.

Chapman-Enskog analysis for finite-volume formulation of lattice Boltzmann equation

The classical Chapman-Enskog expansion is performed for the recently proposed finite-volume formulation of lattice Boltzmann equation (LBE) method D.V. Patil, K.N. Lakshmisha, Finite volume TVD formulation of lattice Boltzmann simulation on unstructured mesh, J. Comput. Phys. 228 (2009) 5262-5279]. First, a modified partial differential equation is derived from a numerical approximation of the discrete Boltzmann equation. Then, the multi-scale, small parameter expansion is followed to recover the continuity and the Navier-Stokes (NS) equations with additional error terms. The expression for apparent value of the kinematic viscosity is derived for finite-volume formulation under certain assumptions. The attenuation of a shear wave, Taylor-Green vortex flow and driven channel flow are studied to analyze the apparent viscosity relation.

A generalized wave action conservative equation for the dissipative dynamical system in the nearshore region

To describe the various complex mechanisms of the dissipative dynamical system between waves, currents, and bottoms in the nearshore region that induce typically the wave motion on large-scale variation of ambient currents, a generalized wave action equation for the dissipative dynamical system in the nearshore region is developed by using the mean-flow equations based on the Navier-Stokes equations of viscous fluid, thus raising two new concepts: the vertical velocity wave action and the dissipative wave action, extending the classical concept, wave action, from the ideal averaged flow conservative system into the real averaged flow dissipative system (that is, the generalized conservative system). It will have more applications.