Molecular dynamics simulations of oscillatory Couette flows with slip boundary conditions

The effect of interfacial slip on steady-state and time-periodic flows of monatomic liquids is investigated using non-equilibrium molecular dynamics simulations. The fluid phase is confined between atomically smooth rigid walls, and the fluid flows are induced by moving one of the walls. In steady shear flows, the slip length increases almost linearly with shear rate. We found that the velocity profiles in oscillatory flows are well described by the Stokes flow solution with the slip length that depends on the local shear rate. Interestingly, the rate dependence of the slip length obtained in steady shear flows is recovered when the slip length in oscillatory flows is plotted as a function of the local shear rate magnitude. For both types of flows, the friction coefficient at the liquid-solid interface correlates well with the structure of the first fluid layer near the solid wall.Comment: 31 pages, 11 figure

The impact of hydraulic retention time on the performance of two configurations of anaerobic pond for municipal sewage treatment

Anaerobic ponds have the potential to contribute to low carbon wastewater treatment, however are currently restricted by long hydraulic residence time (HRT) which leads to large land requirements. A two-stage anaerobic pond (SAP) design was trialled against a single-stage control (CAP) over four HRTs down to 0.5 days, to determine the lowest HRT at which the ponds could operate effectively. No statistical differences were observed in particulate removal between the ponds over all four HRTs, suggesting solids loading is not a critical factor in AP design. Significantly higher biogas production rates were observed in the SAP than the CAP at 1.5 d and 1.0 d HRT, and microbial community profiling suggests the two-stage design may be facilitating spatial separation of the anaerobic digestion process along reactor length. Hydrogenotrophic methanogensis dominated over aceticlastic, with acetate oxidisation a likely degradation pathway. Experimental tracer studies were compared to CFD simulations, with the SAP showing greater hydraulic efficiency, and differences more pronounced at shorter HRTs. Greater flow recirculation between baffles was observed in CFD velocity profiles, demonstrating baffles can dissipate preferential flow patterns and increase effective pond volume, especially at high flow rates. The study demonstrates the potential of APs to be operated at shorter HRTs in psychrophilic conditions, presenting an opportunity for use as pre-treatments (in place of septic tanks) and primary treatment for full wastewater flows. Two-stage designs should be investigated to separate the stages of the anaerobic digestion process by creating preferential conditions along the pond length

Multiscale molecular dynamics/hydrodynamics implementation of two dimensional “Mercedes Benz” water model

A multiscale Molecular Dynamics/Hydrodynamics implementation of the 2D Mercedes Benz (MB or BN2D) [1] water model is developed and investigated. The concept and the governing equations of multiscale coupling together with the results of the two-way coupling implementation are reported. The sensitivity of the multiscale model for obtaining macroscopic and microscopic parameters of the system, such as macroscopic density and velocity fluctuations, radial distribution and velocity autocorrelation functions of MB particles, is evaluated. Critical issues for extending the current model to large systems are discussed

Multiscale simulation of nanofluidic networks of arbitrary complexity

We present a hybrid molecular-continuum method for the simulation of general nanofluidic networks of long and narrow channels. This builds on the multiscale method of Borg et al. (Microfluid Nanofluid 15(4):541–557, 2013; J Comput Phys 233:400–413, 2013) for systems with a high aspect ratio in three main ways: (a) the method has been generalised to accurately model any nanofluidic network of connected channels, regardless of size or complexity; (b) a versatile density correction procedure enables the modelling of compressible fluids; (c) the method can be utilised as a design tool by applying mass-flow-rate boundary conditions (and then inlet/outlet pressures are the output of the simulation). The method decomposes the network into smaller components that are simulated individually using, in the cases in this paper, molecular dynamics micro-elements that are linked together by simple mass conservation and pressure continuity relations. Computational savings are primarily achieved by exploiting length-scale separation, i.e. modelling long channels as hydrodynamically equivalent shorter channel sections. In addition, these small micro-elements reach steady state much quicker than a full simulation of the network does. We test our multiscale method on several steady, isothermal network flow cases and show that it converges quickly (within three iterations) to good agreement with full molecular simulations of the same cases

Hybrid molecular-continuum simulations of water flow through carbon nanotube membranes of realistic thickness

We present new hybrid molecular-continuum simulations of water flow through filtration membranes. The membranes consist of aligned carbon nanotubes (CNTs) of high aspect ratio, where the tube diameters are ~1–2 nm and the tube lengths (i.e. the membrane thicknesses) are 2–6 orders of magnitude larger than this. The flow in the CNTs is subcontinuum, meaning standard continuum fluid equations cannot adequately model the flow; also, full molecular dynamics (MD) simulations are too computationally expensive for modelling these membrane thicknesses. However, various degrees of scale separation in both time and space in this problem can be exploited by a multiscale method: we use the serial-network internal-flow multiscale method (SeN-IMM). Our results from this hybrid method compare very well with full MD simulations of flow cases up to a membrane thickness of 150 nm, beyond which any full MD simulation is computationally intractable. We proceed to use the SeN-IMM to predict the flow in membranes of thicknesses 150 nm–2 μm, and compare these results with both a modified Hagen–Poiseuille flow equation and experimental results for the same membrane configuration. We also find good agreement between experimental and our numerical results for a 1-mm-thick membrane made of CNTs with diameters around 1.1 nm. In this case, the hybrid simulation is orders of magnitude quicker than a full MD simulation would be

Advances in modelling of biomimetic fluid flow at different scales

The biomimetic flow at different scales has been discussed at length. The need of looking into the biological surfaces and morphologies and both geometrical and physical similarities to imitate the technological products and processes has been emphasized. The complex fluid flow and heat transfer problems, the fluid-interface and the physics involved at multiscale and macro-, meso-, micro- and nano-scales have been discussed. The flow and heat transfer simulation is done by various CFD solvers including Navier-Stokes and energy equations, lattice Boltzmann method and molecular dynamics method. Combined continuum-molecular dynamics method is also reviewed

Molecular dynamics pre-simulations for nanoscale computational fluid dynamics

We present a procedure for using molecular dynamics (MD) simulations to provide essential fluid and interface properties for subsequent use in computational fluid dynamics (CFD) calculations of nanoscale fluid flows. The MD pre-simulations enable us to obtain an equation of state, constitutive relations, and boundary conditions for any given fluid/solid combination, in a form that can be conveniently implemented within an otherwise conventional Navier–Stokes solver. Our results demonstrate that these enhanced CFD simulations are then capable of providing good flow field results in a range of complex geometries at the nanoscale. Comparison for validation is with full-scale MD simulations here, but the computational cost of the enhanced CFD is negligible in comparison with the MD. Importantly, accurate predictions can be obtained in geometries that are more complex than the planar MD pre-simulation geometry that provides the nanoscale fluid properties. The robustness of the enhanced CFD is tested by application to water flow along a (15,15) carbon nanotube, and it is found that useful flow information can be obtained

Fluid simulations with atomistic resolution : a hybrid multiscale method with field-wise coupling

We present a new hybrid method for simulating dense fluid systems that exhibit multiscale behaviour, in particular, systems in which a Navier–Stokes model may not be valid in parts of the computational domain. We apply molecular dynamics as a local microscopic refinement for correcting the Navier–Stokes constitutive approximation in the bulk of the domain, as well as providing a direct measurement of velocity slip at bounding surfaces. Our hybrid approach differs from existing techniques, such as the heterogeneous multiscale method (HMM), in some fundamental respects. In our method, the individual molecular solvers, which provide information to the macro model, are not coupled with the continuum grid at nodes (i.e. point-wise coupling), instead coupling occurs over distributed heterogeneous fields (here referred to as field-wise coupling). This affords two major advantages. Whereas point-wise coupled HMM is limited to regions of flow that are highly scale-separated in all spatial directions (i.e. where the state of non-equilibrium in the fluid can be adequately described by a single strain tensor and temperature gradient vector), our field-wise coupled HMM has no such limitations and so can be applied to flows with arbitrarily-varying degrees of scale separation (e.g. flow from a large reservoir into a nano-channel). The second major advantage is that the position of molecular elements does not need to be collocated with nodes of the continuum grid, which means that the resolution of the microscopic correction can be adjusted independently of the resolution of the continuum model. This in turn means the computational cost and accuracy of the molecular correction can be independently controlled and optimised. The macroscopic constraints on the individual molecular solvers are artificial body-force distributions, used in conjunction with standard periodicity. We test our hybrid method on the Poiseuille flow problem for both Newtonian (Lennard-Jones) and non-Newtonian (FENE) fluids. The multiscale results are validated with expensive full-scale molecular dynamics simulations of the same case. Very close agreement is obtained for all cases, with as few as two micro elements required to accurately capture both the Newtonian and non-Newtonian flowfields. Our multiscale method converges very quickly (within 3–4 iterations) and is an order of magnitude more computationally efficient than the full-scale simulation

Multiscale simulation of heat transfer in a rarefied gas

We present a new hybrid method for dilute gas flows that couples a continuum- fluid description to the direct simulation Monte Carlo (DSMC) technique. Instead of using a domain-decomposition framework, we adopt a heterogeneous approach with micro resolution that can capture non-equilibrium or non-continuum fluid behaviour both close to bounding walls and in the bulk. A continuum-fluid model is applied across the entire domain, while DSMC is applied in spatially-distributed micro regions. Using a field-wise coupling approach, each micro element provides a local correction to a continuum sub- region, the dimensions of which are identical to the micro element itself. Interpolating this local correction between the micro elements then produces a correction that can be applied over the entire continuum domain. Key advantages of this method include its suitability for flow problems with varying degrees of scale separation, and that the location of the micro elements is not restricted to the nodes of the computational mesh. Also, the size of the micro elements adapts dynamically with the local molecular mean free path. We demonstrate the method on heat transfer problems in dilute gas flows, where the coupling is performed through the computed heat fluxes. Our test case is micro Fourier flow over a range of rarefaction and temperature conditions: this case is simple enough to enable validation against a pure DSMC simulation, and our results show that the hybrid method can deal with both missing boundary and constitutive information