Source author record

Burkhard Duenweg

Burkhard Duenweg appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

cond-mat.soft cond-mat.stat-mech physics.comp-ph physics.chem-ph physics.flu-dyn

Catalog footprint

What is connected

7works

5topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2020arXiv

High-precision estimate of the hydrodynamic radius for self-avoiding walks

The universal asymptotic amplitude ratio between the gyration radius and the hydrodynamic radius of self-avoiding walks is estimated by high-resolution Monte Carlo simulations. By studying chains of length of up to $N = 2^{25} \approx 34 \times 10^6$ monomers, we find that the ratio takes the value $R_{\mathrm{G}}/R_{\mathrm{H}} = 1.5803940(45)$, which is several orders of magnitude more accurate than the previous state of the art. This is facilitated by a sampling scheme which is quite general, and which allows for the efficient estimation of averages of a large class of observables. The competing corrections to scaling for the hydrodynamic radius are clearly discernible. We also find improved estimates for other universal properties that measure the chain dimension. In particular, a method of analysis which eliminates the leading correction to scaling results in a highly accurate estimate for the Flory exponent of $ν= 0.58759700(40)$.

preprint2020arXiv

Semi-automatic construction of Lattice Boltzmann models

A crucial step in constructing a Lattice Boltzmann model is the definition of a suitable set of lattice velocities, and the correct assignment of the associated weights. For high-order models, the solution of this problem requires a non-trivial effort. The paper outlines the functioning of a publicly available Python script which has been written to assist researchers in that task. The speed of sound $c_s$ is considered as a parameter, which can, within limits, be chosen at will. Under this premise, the Maxwell-Boltzmann constraint equations are a system of linear equations to determine the weights, and hence amenable to numerical solution by standard linear algebra library routines. By suitable contractions, the tensor equations are mapped to a set of equivalent scalar equations, which simplifies the treatment significantly. For a user-supplied set of velocity shells, the software first checks if a solution for the weights exists, and returns it if it also happens to be unique. In such a case, the software also calculates the range of $c_s$ values that yield positive weights. Standard models like D3Q19 with a well-defined special $c_s$ value then result as limiting cases where one of the weights vanishes. In case of an infinite set of solutions, the user may find one particular solution by supplying a $c_s$ value, and then minimizing one or several weights within the framework of standard linear programming. Some examples illustrate the feasibility and usefulness of the approach. A number of models that have been discussed in the literature are nicely reproduced, while the software has also been able to find some new models of even higher order.

preprint2018arXiv

Polymer Solutions

The article provides a brief general introduction into the concepts of scaling, universality, and crossover scaling, plus the blob concept that provides an intuitive picture of crossover phenomena. We present the most important static and dynamic scaling laws for unentangled uncharged polymer solutions, together with their test and refinement by careful computer simulations. A hoard of simulation methods has been developed for these systems, and these will be briefly discussed as well.

preprint2014arXiv

Consistent two-phase Lattice Boltzmann model for gas-liquid systems

A new lattice Boltzmann method for simulating multiphase flows is developed theoretically. The method is adjusted such that its continuum limit is the Navier-Stokes equation, with a driving force derived from the Cahn-Hilliard free energy. In contrast to previous work, however, the bulk and interface terms are decoupled, the former being incorporated into the model through the local equilibrium populations, and the latter through a forcing term. We focus on gas-liquid phase equilibria with the possibility to implement an arbitrary equation of state. The most novel aspect of our approach is a systematic Chapman-Enskog expansion up to the third order. Due to the third-order gradient in the interface forcing term, this is needed for full consistency with both hydrodynamics and thermodynamics. Our construction of a model that satisfies all conditions is based upon previous work by Chen, Goldhirsch, and Orszag (J. Sci. Comp. 34, 87 (2008)), and implies 59 and 21 velocities in three and two dimensions, respectively. Applying the conditions of positivity of weights, existence of a two-phase region in the phase diagram, and positivity of the bulk viscosity, we find substantial restrictions on the permitted equation of state, which can only be lifted by an even more refined model. Moreover, it turns out that it is necessary to solve a self-consistent equation for the hydrodynamic flow velocity, in order to enforce the identity of momentum density and mass current on the lattice. The analysis completely identifies all spurious terms in the Navier-Stokes equation, and thus shows how to systematically eliminate each of them, by constructing a suitable collision operator. (continued - see article PDF)

preprint2010arXiv

Collapse dynamics of copolymers in a poor solvent: Influence of hydrodynamic interactions and chain sequence

We investigate the dynamics of the collapse of a single copolymer chain, when the solvent quality is suddenly quenched from good to poor. We employ Brownian dynamics simulations of a bead-spring chain model and incorporate fluctuating hydrodynamic interactions via the Rotne-Prager-Yamakawa tensor. Various copolymer architectures are studied within the framework of a two-letter HP model, where monomers of type H (hydrophobic) attract each other, while all interactions involving P (polar or hydrophilic) monomers are purely repulsive. The hydrodynamic interactions are found to assist the collapse. Furthermore, the chain sequence has a strong influence on the kinetics and on the compactness and energy of the final state. The dynamics is typically characterised by initial rapid cluster formation, followed by coalescence and final rearrangement to form the compact globule. The coalescence stage takes most of the collapse time, and its duration is particularly sensitive to the details of the architecture. Long blocks of type P are identified as the main bottlenecks to find the globular state rapidly.

preprint2009arXiv

Implicit and explicit solvent models for the simulation of a single polymer chain in solution: Lattice Boltzmann vs Brownian dynamics

We present a comparative study of two computer simulation methods to obtain static and dynamic properties of dilute polymer solutions. The first approach is a recently established hybrid algorithm based upon dissipative coupling between Molecular Dynamics and lattice Boltzmann (LB), while the second is standard Brownian Dynamics (BD) with fluctuating hydrodynamic interactions. Applying these methods to the same physical system (a single polymer chain in a good solvent in thermal equilibrium) allows us to draw a detailed and quantitative comparison in terms of both accuracy and efficiency. It is found that the static conformations of the LB model are distorted when the box length L is too small compared to the chain size. Furthermore, some dynamic properties of the LB model are subject to an $L^{-1}$ finite size effect, while the BD model directly reproduces the asymptotic $L \to \infty$ behavior. Apart from these finite size effects, it is also found that in order to obtain the correct dynamic properties for the LB simulations, it is crucial to properly thermalize all the kinetic modes. Only in this case, the results are in excellent agreement with each other, as expected. Moreover, Brownian Dynamics is found to be much more efficient than lattice Boltzmann as long as the degree of polymerization is not excessively large.

preprint2002arXiv

Corrections to Scaling in the Hydrodynamic Properties of Dilute Polymer Solutions

We discuss the hydrodynamic radius $R_H$ of polymer chains in good solvent, and show that the leading order correction to the asymptotic law $R_H \propto N^ν$ ($N$ degree of polymerization, $ν\approx 0.59$) is an ``analytic'' term of order $N^{-(1 - ν)}$, which is directly related to the discretization of the chain into a finite number of beads. This result is further corroborated by exact calculations for Gaussian chains, and extensive numerical simulations of different models of good--solvent chains, where we find a value of $1.591 \pm 0.007$ for the asymptotic universal ratio $R_G / R_H$, $R_G$ being the chain's gyration radius. For $Θ$ chains the data apparently extrapolate to $R_G / R_H \approx 1.44$, which is different from the Gaussian value 1.5045, but in accordance with previous simulations. We also show that the experimentally observed deviations of the initial decay rate in dynamic light scattering from the asymptotic Benmouna--Akcasu value can partly be understood by similar arguments.