Particle dynamics in drying colloidal solution using discrete particle method

We study particle dynamics in drying colloidal solutions using the numerical simulation with discrete particle method (DPM). Simulations of two different systems were conducted; the drying dynamics of monodispersed and binary mixture of colloidal solution, and compared with those from the previous studies. In the monodispersed colloidal solution, the time evolution of particle concentration profile for varying Péclet number was simulated with the same initial particle concentration. In the binary colloidal solution, when the particle size ratio α is 3, three different stratification modes were observed varying Péclet number and initial particle concentration. By comparison, our method was in a good agreement with the existing methods. Additionally, because of the mesh-based Eulerian approach in our model, other various multi-physical phenomena, such as effect of thermal Marangoni or chemical reaction, can be included in an easy way. From the results, we expect that this work can provide a physical insight for predicting the quality of colloidal drying in a complicated situation.

Graphically interpreting how incision thresholds influence topographic and scaling properties of modeled landscapes

We examine the influence of incision thresholds on topographic and scaling properties of landscapes that follow a landscape evolution model (LEM) with terms for stream-power incision, linear diffusion, and uniform uplift. Our analysis uses three main tools. First, we examine the graphical behavior of theoretical relationships between curvature and the steepness index (which depends on drainage area and slope). These relationships plot as straight lines for the case of steady-state landscapes that follow the LEM. These lines have slopes and intercepts that provide estimates of landscape characteristic scales. Such lines can be viewed as counterparts of slope–area relationships, which follow power laws in detachment-limited landscapes but not in landscapes with diffusion. We illustrate the response of these curvature–steepness index lines to changes in the values of parameters. Second, we define a Péclet number that quantifies the competition between incision and diffusion, while taking the incision threshold into account. We examine how this Péclet number captures the influence of the incision threshold on the degree of landscape dissection. Third, we characterize the influence of the incision threshold using a ratio between it and the steepness index. This ratio is a dimensionless number in the case of the LEM that we use and reflects the fraction by which the incision rate is reduced due to the incision threshold; in this way, it quantifies the relative influence of the incision threshold across a landscape. These three tools can be used together to graphically illustrate how topography and process competition respond to incision thresholds.

Flow Rate Prediction for a Semi-permeable Membrane at Low Reynolds Number in a Circular Pipe

A concise and accurate prediction method is required for membrane permeability in chemical engineering and biological fields. As a preliminary study on this topic, we propose the concentration polarization model (CPM) of the permeate flux and flow rate under dominant effects of viscosity and solute diffusion. In this model, concentration polarization is incorporated for the solution flow through a semi-permeable membrane (i.e., permeable for solvent but not for solute) in a circular pipe. The effect of the concentration polarization on the flow field in a circular pipe under a viscous-dominant condition (i.e., at a low Reynolds number) is discussed by comparing the CPM with the numerical simulation results and infinitesimal Péclet number model (IPM) for the membrane permeability, strength of the osmotic pressure, and Péclet number. The CPM and IPM are confirmed to be a reasonable extension of the model for a pure fluid, which was proposed previously. The application range of the IPM is narrow because the advection of the solute concentration is not considered, whereas the CPM demonstrates superior applicability in a wide range of parameters, including the permeability coefficient, strength of the osmotic pressure, and Péclet number. This suggests the necessity for considering concentration polarization. Although the mathematical expression of the CPM is more complex than that of the IPM, the CPM exhibits a potential to accurately predict the permeability parameters for a condition in which a large permeate flux and osmotic pressure occur.

Exact solutions for the formation of stagnant caps of insoluble surfactant on a planar free surface

A class of exact solutions is presented describing the time evolution of insoluble surfactant to a stagnant cap equilibrium on the surface of deep water in the Stokes flow regime at zero capillary number and infinite surface Péclet number. This is done by demonstrating, in a two-dimensional model setting, the relevance of the forced complex Burgers equation to this problem when a linear equation of state relates the surface tension to the surfactant concentration. A complex-variable version of the method of characteristics can then be deployed to find an implicit representation of the general solution. A special class of initial conditions is considered for which the associated solutions can be given explicitly. The new exact solutions, which include both spreading and compactifying scenarios, provide analytical insight into the unsteady formation of stagnant caps of insoluble surfactant. It is also shown that first-order reaction kinetics modelling sublimation or evaporation of the insoluble surfactant to the upper gas phase can be incorporated into the framework; this leads to a forced complex Burgers equation with linear damping. Generalized exact solutions to the latter equation at infinite surface Péclet number are also found and used to study how reaction effects destroy the surfactant cap equilibrium.

Mass transfer from small spheroids suspended in a turbulent fluid

By coupling direct numerical simulation of homogeneous isotropic turbulence with a localised solution of the convection–diffusion equation, we model the rate of transfer of a solute (mass transfer) from the surface of small, neutrally buoyant, axisymmetric, ellipsoidal particles (spheroids) in dilute suspension within a turbulent fluid at large Péclet number, $\textit {Pe}$ . We observe that, at $\textit {Pe} = O(10)$ , the average transfer rate for prolate spheroids is larger than that of spheres with equivalent surface area, whereas oblate spheroids experience a lower average transfer rate. However, as the Péclet number is increased, oblate spheroids can experience an enhancement in mass transfer relative to spheres near an optimal aspect ratio $\lambda \approx 1/4$ . Furthermore, we observe that, for spherical particles, the Sherwood number $\textit {Sh}$ scales approximately as $\textit {Pe}^{0.26}$ over $\textit {Pe} = 1.4\times 10^{1}$ to $1.4\times 10^{4}$ , which is below the $\textit {Pe}^{1/3}$ scaling observed for inertial particles but consistent with available experimental data for tracer-like particles. The discrepancy is attributed to the diffusion-limited temporal response of the concentration boundary layer to turbulent strain fluctuations. A simple model, the quasi-steady flux model, captures both of these phenomena and shows good quantitative agreement with our numerical simulations.

Arterial pharmacokinetics in a patient-specific atherosclerotic artery-a simulation study

Of concern in the paper is a numerical study of endovascular drug delivery in a patient-specific atherosclerotic artery through a mathematical model in which the luminal flow is governed by an incompressible vis- cous Newtonian fluid, and the transport of luminal as well as tissue concentration is modeled as an unsteady convection-diffusion process. An image processing technique has been successfully adopted to detect the edges of the computational domain extracted from an asymmetric (about the centerline of the artery) patient-specific atherosclerotic artery. Considering each pixel as a control volume, the Marker and Cell (MAC) method has been leveraged to get a quantitative insight of the model considered by exploiting physiologically realistic initial, boundary as well as interface conditions. Simulated results reveal that the number as well as the length of separation zone does increase with increasing Re, and the near-wall velocity contour might be important for estimating the near-wall residence time for the pool of drug molecules available for tissue uptake. Results also show that the more the tissue porosity and interface roughness do not necessarily imply the more the effective- ness of delivery, even though they enhance the averaged concentration in the tissue domains, and also the area under concentration diminishes with increasing Peclet number. Thus, the tissue porosity, the Peclet number and various geometrical shapes (interface roughness) play a pivotal role in the dispersion and the effectiveness of drug delivery. GANITJ. Bangladesh Math. Soc.41.1 (2021) 62-77

Challenges in modelling diffusiophoretic transport

The methodology to simulate transport phenomena in bulk systems is well-established. In contrast, there is no clear consensus about the choice of techniques to model cross-transport phenomena and phoretic transport, mainly because some of the hydrodynamic descriptions are incomplete from a thermodynamic point of view. In the present paper, we use a unified framework to describe diffusio-osmosis(phoresis), and we report non-equilibrium molecular dynamics (NEMD) on such systems. We explore different simulation methods to highlight some of the technical problems that arise in the calculations. For diffusiophoresis, we use two NEMD methods: boundary-driven and field-driven. Although the two methods should be equivalent in the limit of very weak gradients, we find that finite Peclet-number effects are much stronger in boundary-driven flows than in the case where we apply fictitious color forces. Graphic abstract

Displacement washing of softwood pulp cooked to various levels of residual lignin content

This study investigates the influence of the degree of delignification of kraft spruce pulp cooked at seven different kappa numbers, ranging from 18.1 to 50.1, on the efficiency of displacement washing under laboratory conditions. Although the pulp bed is a polydispersive and heterogeneous system, the correlation dependence of the wash yield and bed efficiency on the Péclet number and the kappa number of the pulp showed that washing efficiency increased not only with an increasing Péclet number, but also with an increasing kappa number. The linear dependence between the mean residence time of the solute lignin in the bed and the space time, which reflects the residence time of the wash liquid in the pulp bed, was found for all levels of the kappa number. Washing also reduced the kappa number and the residual lignin content in the pulp fibers.