Advances in Hybrid Molecular/Continuum Methods for Micro and Nano Flows

This paper was presented at the 4th Micro and Nano Flows Conference (MNF2014), which was held at University College, London, UK. The conference was organised by Brunel University and supported by the Italian Union of Thermofluiddynamics, IPEM, the Process Intensification Network, the Institution of Mechanical Engineers, the Heat Transfer Society, HEXAG - the Heat Exchange Action Group, and the Energy Institute, ASME Press, LCN London Centre for Nanotechnology, UCL University College London, UCL Engineering, the International NanoScience Community, www.nanopaprika.eu.Next generation uid ow systems are likely to depend on micro or nano scale dynamics that make the system behaviour multiscale in both space and time. There may be strong or weak separation between the length scales and between the time scales in di erent parts of the ow, and these scale-separations may also vary in space and time. In this paper we discuss a practical approach to improving the e ciency of hybrid particle/continuum models of such multiscale ows. Our focus is on adapting the solution method to the local scale-separation conditions, in order to balance compu- tational e ciency with accuracy. We compare results from our new hybridisation in space and time with full molecular simulations of benchmark nanoscale ows

A lumped-parameter model of the cerebrospinal system for simulating syringomyelia

This paper was presented at the 2nd Micro and Nano Flows Conference (MNF2009), which was held at Brunel University, West London, UK. The conference was organised by Brunel University and supported by the Institution of Mechanical Engineers, IPEM, the Italian Union of Thermofluid dynamics, the Process Intensification Network, HEXAG - the Heat Exchange Action Group and the Institute of Mathematics and its Applications.Syringomyelia is a disease in which high-pressure fluid-filled cavities, called syrinxes, form in the spinal cord (SC) which can cause progressive loss of sensory and motor functions. Poor treatment outcomes have led to myriad hypotheses for its pathogenesis, which unfortunately are often based on small numbers of patients due to the relative rarity of the disease. However, accumulating evidence in the last decade from animal studies implicates arterial pulsations in syrinx formation. In particular, it has been suggested that a phase difference between the pressure pulse in the spinal subarachnoid space and the perivascular spaces, due to a pathologically disturbed blood supply, could result in a net influx of cerebrospinal fluid (CSF) into the SC. A lumped-parameter model is developed of the cerebrospinal system to investigate this conjecture. It is found that although this phase-lag mechanism may operate, it requires the SC to have an intrinsic storage capacity due to the collapsibility of the contained venous reservoir. If this storage requirement is met then the results presented here suggest that, on mechanical grounds, a syringo-subarachnoid shunt may be a better surgical treatment option than a subarachnoid bypass for posttraumatic syringomyelia

Gas dynamics at the micro-scale: A review of progress in hydrodynamic modelling

This paper was presented at the 2nd Micro and Nano Flows Conference (MNF2009), which was held at Brunel University, West London, UK. The conference was organised by Brunel University and supported by the Institution of Mechanical Engineers, IPEM, the Italian Union of Thermofluid dynamics, the Process Intensification Network, HEXAG - the Heat Exchange Action Group and the Institute of Mathematics and its Applications.We review some recent developments in the modelling of non-equilibrium (rarefied) gas flows at the micro- and nano-scale using extended hydrodynamic models. Following a brief exposition of the challenges that non-equilibrium poses in micro- and nano-scale gas flows, we outline the field of extended hydrodynamics, describing the effective abandonment of Burnett-type models in favour of high-order regularised moment equations. We then review the boundary conditions required if the conventional Navier-Stokes-Fourier (NSF) fluid dynamic model is applied at the micro scale, describing how 2nd-order Maxwelltype conditions can be used to compensate for some of the non-equilibrium flow behaviour near solid surfaces. While extended hydrodynamics is not yet widely-used for real flow problems because of its inherent complexity, we finish with an outline of recent ‘phenomenological extended hydrodynamics’ (PEH) techniques — essentially the NSF equations scaled to incorporate non-equilibrium behaviour close to solid surfaces — which offer promise as engineering models.This work is funded in the UK by the Engineering and Physical Sciences Research Council through grants EP/F002467/1, EP/D07455X/1, EP/D007488/1 and EP/F028865/1

Coupled continuum hydrodynamics and molecular dynamics method for multiscale simulation

This paper was presented at the 3rd Micro and Nano Flows Conference (MNF2011), which was held at the Makedonia Palace Hotel, Thessaloniki in Greece. The conference was organised by Brunel University and supported by the Italian Union of Thermofluiddynamics, Aristotle University of Thessaloniki, University of Thessaly, IPEM, the Process Intensification Network, the Institution of Mechanical Engineers, the Heat Transfer Society, HEXAG - the Heat Exchange Action Group, and the Energy Institute.We present a new hybrid methodology for carrying out multiscale simulations of flow problems lying between continuum hydrodynamics and molecular dynamics, where macro/micro lengthscale separation exists only in one direction. Our multiscale method consists of an iterative technique that couples mass and momentum flux between macro and micro domains, and is tested on a converging/diverging nanochannel case containing flow of a simple Lennard-Jones liquid. Comparisons agree well with a full MD simulation of the same test case.EPSRC Programme Grant EP/I011927/

A New Heterogeneous Multiscale Technique for Microscale Gas Flows

This paper was presented at the 4th Micro and Nano Flows Conference (MNF2014), which was held at University College, London, UK. The conference was organised by Brunel University and supported by the Italian Union of Thermofluiddynamics, IPEM, the Process Intensification Network, the Institution of Mechanical Engineers, the Heat Transfer Society, HEXAG - the Heat Exchange Action Group, and the Energy Institute, ASME Press, LCN London Centre for Nanotechnology, UCL University College London, UCL Engineering, the International NanoScience Community, www.nanopaprika.eu.We present a new hybrid method for dilute gas flows that heterogeneously couples a continuumfluid description to the direct simulation Monte Carlo (DSMC) method. 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, DSMC sub-domains of any size can be placed at any location. The sub-domain arrangement can therefore be adjusted for each problem to capture non-equilibrium behaviour both close to bounding walls and in the bulk. We demonstrate our method on a test case of high-speed micro Couette flow. With large differences in wall velocity, significant viscous heating is present, and so our coupling considers the transfer of both momentum and heat. Our hybrid results are validated against a pure DSMC simulation, and the results show that the method can deal with missing boundary and constitutive information

Exploiting timescale separation in micro and nano flows

This paper was presented at the 3rd Micro and Nano Flows Conference (MNF2011), which was held at the Makedonia Palace Hotel, Thessaloniki in Greece. The conference was organised by Brunel University and supported by the Italian Union of Thermofluiddynamics, Aristotle University of Thessaloniki, University of Thessaly, IPEM, the Process Intensification Network, the Institution of Mechanical Engineers, the Heat Transfer Society, HEXAG - the Heat Exchange Action Group, and the Energy Institute.In this paper we describe how timescale separation in micro/nano flows can be exploited for computational acceleration. A modified version of the seamless heterogenous multiscale method (SHMM) is proposed: a multi-step SHMM. This maintains the main advantages of SHMM (e.g., re-initialisation of micro data is not required; temporal gearing (computational speed-up) is easily controlled; and it is applicable to full and intermediate degrees of timescale separation) while improving on accuracy and greatly reducing the number of macroscopic computations and micro/macro coupling instances required. The improved accuracy of the multi-step SHMM is demonstrated for two canonical one-dimensional transient flows (oscillatory Poiseuille and oscillatory Couette flow) and for rarefied-gas oscillatory Poiseuille flow.This research is financially supported by the EPSRC Programme Grant EP/I011927/1

Simulation of the Head-disk Interface Gap Using a Hybrid Multi-scale Method

We present a hybrid multi-scale method that provides a capability to capture the disparate scales associated with modelling flow in micro- and nano-devices. Our model extends the applicability of an internal-flow multi-scale method by providing a framework to couple the internal (small scale) flow regions to the external (large scale) flow regions. We demonstrate the application of both the original methodology and the new hybrid approach to model the flow field in the vicinity of the head-disk interface gap of a hard disk drive enclosure. The internal flow regions within the gap are modelled by an extended internal-flow multi-scale method that utilises a finite-difference scheme for non-uniform grids. Our proposed hybrid multi-scale method is then employed to couple the internal micro-flow region to the flow external to the gap, to capture entrance/exit effects. We also demonstrate the successful application of the method in capturing other localised phenomena (e.g. those due to localised wall heating)

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

Radiation Campaign of HPK Prototype LGAD sensors for the High-Granularity Timing Detector (HGTD)

We report on the results of a radiation campaign with neutrons and protons of Low Gain Avalanche Detectors (LGAD) produced by Hamamatsu (HPK) as prototypes for the High-Granularity Timing Detector (HGTD) in ATLAS. Sensors with an active thickness of 50~$\mu$m were irradiated in steps of roughly 2$\times$ up to a fluence of $3\times10^{15}~\mathrm{n_{eq}cm^{-2}}$. As a function of the fluence, the collected charge and time resolution of the irradiated sensors will be reported for operation at $-30^{\circ}$