We apply a previously developed asymptotic model (J. Fluid. Mech. 915, A133 (2021)) to study instabilities of free surface films of nanometric thickness on thermally conductive substrates in two and three spatial dimensions. While the specific focus
is on metal films exposed to laser heating, the model itself applies to any setup involving films on the nanoscale whose material parameters are temperature-dependent. For the particular case of metal films heated from above, an important aspect is that the considered heating is volumetric, since the absorption length of the applied laser pulse is comparable to the film thickness. In such a setup, absorption of thermal energy and film evolution are closely correlated and must be considered self-consistently. The asymptotic model allows for a significant simplification, which is crucial from both modeling and computational points of view, since it allows for asymptotically correct averaging of the temperature over the film thickness. We find that the properties of the thermally conductive substrate -- in particular its thickness and rate of heat loss -- play a critical role in controlling the film temperature and dynamics. The film evolution is simulated using efficient GPU-based simulations which, when combined with the developed asymptotic model, allow for fully nonlinear time-dependent simulations in large three-dimensional computational domains. In addition to uncovering the role of the substrate and its properties in determining the film evolution, one important finding is that, at least for the considered range of material parameters, strong in-plane thermal diffusion in the film results in negligible spatial variations of temperature, and the film evolution is predominantly influenced by temporal variation of film viscosity and surface tension (dictated by average film temperature), as well as thermal conductivity of the substrate.
Filtration of feed containing multiple species of particles is a common process in the industrial setting. In this work we propose a model for filtration of a suspension containing an arbitrary number of particle species, each with different affiniti
es for the filter membrane. We formulate a number of optimization problems pertaining to effective separation of desired and undesired particles in the special case of two particle species and we present results showing how properties such as feed composition affect the optimal filter design (internal pore structure). In addition, we propose a novel multi-stage filtration strategy, which provides a significant mass yield improvement for the desired particles, and, surprisingly, higher purity of the product as well.
We analyze the performance of membrane filters represented by pore networks using two criteria: 1) total volumetric throughput of filtrate over the filter lifetime and 2) accumulated foulant concentration in the filtrate. We first formulate the gover
ning equations of fluid flow on a general network, and we model transport and adsorption of particles (foulants) within the network by imposing an advection equation with a sink term on each pore (edge) as well as conservation of fluid and foulant volumetric flow rates at each pore junction (network vertex). Such a setup yields a system of partial differential equations on the network. We study the influence of three geometric network parameters on filter performance: 1) average number of neighbors of each vertex; 2) initial total void volume of the pore network; and 3) tortuosity of the network. We find that total volumetric throughput depends more strongly on the initial void volume than on average number of neighbors. Tortuosity, however, turns out to be a universal parameter, leading to almost perfect collapse of all results for a variety of different network architectures. In particular, the accumulated foulant concentration in the filtrate shows an exponential decay as tortuosity increases.
We consider a mathematical model that describes the flow of a Nematic Liquid Crystal (NLC) film placed on a flat substrate, across which a spatially-varying electric potential is applied. Due to their polar nature, NLC molecules interact with the (no
nuniform) electric field generated, leading to instability of a flat film. Implementation of the long wave scaling leads to a partial differential equation that predicts the subsequent time evolution of the thin film. This equation is coupled to a boundary value problem that describes the interaction between the local molecular orientation of the NLC (the director field) and the electric potential. We investigate numerically the behavior of an initially flat film for a range of film heights and surface anchoring conditions.
We consider a free surface thin film placed on a thermally conductive substrate and exposed to an external heat source in a setup where the heat absorption depends on the local film thickness. Our focus is on modeling film evolution while the film is
molten. The evolution of the film modifies local heat flow, which in turn may influence the film surface evolution through thermal variation of the films material properties. Thermal conductivity of the substrate plays an important role in determining the heat flow and the temperature field in the evolving film and in the substrate itself. In order to reach a tractable formulation, we use asymptotic analysis to develop a novel thermal model that is accurate, computationally efficient, and that accounts for the heat flow in both the in-plane and out-of plane directions. We apply this model to metal films of nanoscale thickness exposed to heating and melting by laser pulses, a setup commonly used for self and directed assembly of various metal geometries via dewetting while the films are in the liquid phase. We find that thermal effects play an important role, and in particular that the inclusion of temperature dependence in the metal viscosity modifies the time scale of the evolution significantly. On the other hand, in the considered setup the Marangoni (thermocapillary) effect turns out to be insignificant.
Pleated membrane filters, which offer larger surface area to volume ratios than unpleated membrane filters, are used in a wide variety of applications. However, the performance of the pleated filter, as characterized by a flux-throughput plot, indica
tes that the equivalent unpleated filter provides better performance under the same pressure drop. Earlier work (Sanaei & Cummings 2016) used a highly-simplified membrane model to investigate how the pleating effect and membrane geometry affect this performance differential. In this work, we extend this line of investigation and use asymptotic methods to couple an outer problem for the flow within the pleated structure to an inner problem that accounts for the pore structure within the membrane. We use our new model to formulate and address questions of optimal membrane design for a given filtration application.
Dense, stabilized, frictional particulate suspensions in a viscous liquid undergo increasingly strong continuous shear thickening (CST) as the solid packing fraction, $phi$, increases above a critical volume fraction, and discontinuous shear thickeni
ng (DST) is observed for even higher packing fractions. Recent studies have related shear thickening to a transition from mostly lubricated to predominantly frictional contacts with the increase in stress. The rheology and networks of frictional forces from two and three-dimensional simulations of shear-thickening suspensions are studied. These are analyzed using measures of the topology of the network, including tools of persistent homology. We observe that at low stress the frictional interaction networks are predominantly quasi-linear along the compression axis. With an increase in stress, the force networks become more isotropic, forming loops in addition to chain-like structures. The topological measures of Betti numbers and total persistence provide a compact means of describing the mean properties of the frictional force networks and provide a key link between macroscopic rheology and the microscopic interactions. A total persistence measure describing the significance of loops in the force network structure, as a function of stress and packing fraction, shows behavior similar to that of relative viscosity and displays a scaling law near the jamming fraction for both dimensionalities simulated.
We present the results of large scale simulations of 4th order nonlinear partial differential equations of dif- fusion type that are typically encountered when modeling dynamics of thin fluid films on substrates. The simulations are based on the alte
rnate direction implicit (ADI) method, with the main part of the compu- tational work carried out in the GPU computing environment. Efficient and accurate computations allow for simulations on large computational domains in three spatial dimensions (3D) and for long computational times. We apply the methods developed to the particular problem of instabilities of thin fluid films of nanoscale thickness. The large scale of the simulations minimizes the effects of boundaries, and also allows for simulating domains of the size encountered in published experiments. As an outcome, we can analyze the development of instabilities with an unprecedented level of detail. A particular focus is on analyzing the manner in which instability develops, in particular regarding differences between spinodal and nucleation types of dewetting for linearly unstable films, as well as instabilities of metastable films. Simulations in 3D allow for consideration of some recent results that were previously obtained in the 2D geometry (J. Fluid Mech. 841, 925 (2018)). Some of the new results include using Fourier transforms as well as topological invariants (Betti numbers) to distinguish the outcomes of spinodal and nucleation types of instabilities, describing in precise terms the complex processes that lead to the formation of satellite drops, as well as distinguishing the shape of the evolving film front in linearly unstable and metastable regimes. We also discuss direct comparison between simulations and available experimental results for nematic liquid crystal and polymer films.
We discuss instabilities of fluid films of nanoscale thickness, with a particular focus on films where the destabilising mechanism allows for linear instability, metastability, and absolute stability. Our study is motivated by nematic liquid crystal
films; however we note that similar instability mechanisms, and forms of the effective disjoining pressure, appear in other contexts, such as the well-studied problem of polymeric films on two-layered substrates. The analysis is carried out within the framework of the long-wave approximation, which leads to a fourth order nonlinear partial different equation for the film thickness. Within the considered formulation, the nematic character of the film leads to an additional contribution to the disjoining pressure, changing its functional form. This effective disjoining pressure is characterised by the presence of a local maximum for non-vanishing film thickness. Such a form leads to complicated instability evolution that we study by analytical means, including application of marginal stability criteria, and by extensive numerical simulations that help us develop a better understanding of instability evolution in the nonlinear regime. This combination of analytical and computational techniques allows us to reach novel understanding of relevant instability mechanisms, and of their influence on transient and fully developed fluid film morphologies.
When dense granular systems are exposed to external forcing, they evolve on the time scale that is typically related to the externally imposed one (shear or compression rate, for example). This evolution could be characterized by observing temporal e
volution of contact networks. However, it is not immediately clear whether the force networks, defined on contact networks by considering force interactions between the particles, evolve on a similar time scale. To analyze the evolution of these networks, we carry out discrete element simulations of a system of soft frictional disks exposed to compression that leads to jamming. By using the tools of computational topology, we show that close to jamming transition, the force networks evolve on the time scale which is much faster than the externally imposed one. The presentation will discuss the factors that determine this fast time scale.
