No abstract

In this paper, we review the existing interpolation functions and introduce a finite element interpolation function to be used in the immersed boundary and finite element methods. This straightforward finite element interpolation function for unstructured grids enables us to obtain a sharper interface that yields more accurate interfacial solutions. The solution accuracy is compared with the existing interpolation functions such as the discretized Dirac delta function and the reproducing kernel interpolation function. The finite element shape function is easy to implement and it naturally satisfies the reproducing condition. They are interpolated through only one element layer instead of smearing to several elements. A pressure jump is clearly captured at the fluid-solid interface. Two example problems are studied and results are compared with other numerical methods. A convergence test is thoroughly conducted for the independent fluid and solid meshes in a fluid-structure interaction system. The required mesh size ratio between the fluid and solid domains is obtained.

The thermostats in molecular dynamics (MD) simulations of highly confined channel flow may have significant influences on the fidelity of transport phenomena. In this study, we exploit non-equilibrium MD simulations to generate Couette flows with different combinations of thermostat algorithms and strategies. We provide a comprehensive analysis on the effectiveness of three thermostat algorithms Nosé-Hoover chain (NHC), Langevin (LGV) and dissipative particle dynamics (DPD) when applied in three thermostat strategies, thermostating either walls (TW) or fluid (TF), and thermostating both the wall and fluid (TWTF). Our results of thermal and mechanical properties show that the TW strategy more closely resembles experimental conditions. The TF and TWTF systems also produce considerably similar behaviors in weakly sheared systems, but deviate the dynamics in strongly sheared systems due to the isothermal condition. The LGV and DPD thermostats used in the TF and TWTF systems provide vital ways to yield correct dynamics in coarse-grained systems by tuning the fluid transport coefficients. Using conventional NHC thermostat to thermostat fluid only produces correct thermal behaviors in weakly sheared systems, and breaks down due to significant thermal inhomogeneity in strongly sheared systems.

In this paper, nanoscale wetting on groove-patterned surfaces is thoroughly studied using molecular dynamics simulations. The results are compared with Wenzel's and Cassie's predictions to determine whether these continuum theories are still valid at the nanoscale for both hydrophobic and hydrophilic types of surfaces when the droplet size is comparable to the groove size. A system with a liquid mercury droplet and grooved copper substrate is simulated. The wetting properties are determined by measuring contact angles of the liquid droplet at equilibrium states. Correlations are established between the contact angle, roughness factor r, and surface fraction f. The results show that, for hydrophobic surfaces, the contact angle as a function of roughness factor and surface fraction on nanogrooved surfaces obeys the predictions from Wenzel's theory for wetted contacts and Cassie's theory for composite contacts. However, slight deviations occur in composite contacts when a small amount of liquid penetration is observed. The contact angle of this partial wetting cannot be accurately predicted using either Cassie's or Wenzel's theories. For hydrophilic surfaces, only wetted contacts are observed. In most cases, the resulting contact angles are found to be higher than Wenzel's predictions. At the nanoscale, high surface edge density plays an important role, which results in contact line pinning near plateau edges. For both hydrophobic and hydrophilic surfaces, substantial amount of anistropic spreading is found in the direction that is parallel to the grooves, especially at wetted or partially wetted contacts.

Through nanomedicine, game-changing methods are emerging to deliver drug molecules directly to diseased areas. One of the most promising of these is the targeted delivery of drugs and imaging agents via drug carrier-based platforms. Such drug delivery systems can now be synthesized from a wide range of different materials, made in a number of different shapes, and coated with an array of different organic molecules, including ligands. If optimized, these systems can enhance the efficacy and specificity of delivery compared with those of non-targeted systems. Emerging integrated multiscale experiments, models and simulations have opened the door for endless medical applications. Current bottlenecks in design of the drug-carrying particles are the lack of knowledge about the dispersion of these particles in the microvasculature and of their subsequent internalization by diseased cells (Bao et al . 2014 J. R. Soc. Interface 11 , 20140301 ( doi:10.1098/rsif.2014.0301 )). We describe multiscale modelling techniques that study how drug carriers disperse within the microvasculature. The immersed molecular finite-element method is adopted to simulate whole blood including blood plasma, red blood cells and nanoparticles. With a novel dissipative particle dynamics method, the beginning stages of receptor-driven endocytosis of nanoparticles can be understood in detail. Using this multiscale modelling method, we elucidate how the size, shape and surface functionality of nanoparticles will affect their dispersion in the microvasculature and subsequent internalization by targeted cells.