Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userThe influence of thermal effects on the breakup of thin films of nanometric thickness
97
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
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 alternate 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.
A charged droplet can be electrodynamically levitated in the air using a quadrupole trap by typically applying a sinusoidal electric field. When a charged drop is levitated it exhibits surface oscillations simultaneously building charge density due to continuous evaporation and subsequently undergoes breakup due to Rayleigh instability. In this work, we examined large-amplitude surface oscillations of a sub-Rayleigh charged drop and its subsequent breakup, levitated by various applied signals such as sine, square and ramp waveform at various imposed frequencies, using high-speed imaging (recorded at 100-130 thousand Frames Per Second (fps)). It is observed that the drop surface oscillates in sphere-prolate-sphere-oblate (SPSO) mode and seldom in the sphere-prolate-sphere (SPS) mode depending on the intricate interplay of various forces due to charge(q), the intensity of applied field ($Lambda$) and shift of the droplet from the geometric center of the trap ($z_{shift}$). The Fast Fourier Transformation (FFT) analysis shows that the droplet oscillates with the forced frequency irrespective of the type of the applied waveform. While in the sinusoidal case, the nonlinearities are significant, in the square and ramp potentials, there is an admittance of all the harmonic frequencies of the applied potential. Interestingly, the breakup characteristics of a critically charged droplet is found to be unaffected by the type of the applied waveform. The experimental observations are validated with an analytical theory as well as with the Boundary Integral (BI) simulations in the potential flow limit and the results are found to be in a reasonable agreement.
Recent experiments of thin films flowing down a vertical fiber with varying nozzle diameters present a wealth of new dynamics that illustrate the need for more advanced theory. We present a detailed analysis using a full lubrication model that includes slip boundary conditions, nonlinear curvature terms, and a film stabilization term. This study brings to focus the presence of a stable liquid layer playing an important role in the full dynamics. We propose a combination of these physical effects to explain the observed velocity and stability of traveling droplets in the experiments and their transition to isolated droplets. This is also supported by stability analysis of the traveling wave solution of the model.
Turbulence structure resulting from multi-fluid or multi-species, variable-density isotropic turbulence interaction with a Mach 2 shock is studied using turbulence-resolving shock-capturing simulations and Eulerian (grid) and Lagrangian (particle) methods. The complex roles density play in the modification of turbulence by the shock wave are identified. Statistical analyses of the velocity gradient tensor (VGT) show that the density variations significantly change the turbulence structure and flow topology. Specifically, a stronger symmetrization of the joint probability density function (PDF) of second and third invariants of the anisotropic velocity gradient tensor, PDF$(Q^ast, R^ast)$, as well as the PDF of the vortex stretching contribution to the enstrophy equation, are observed in the multi-species case. Furthermore, subsequent to the interaction with the shock, turbulent statistics also acquire a differential distribution in regions having different densities. This results in a nearly symmetrical PDF$(Q^ast, R^ast)$ in heavy fluid regions, while the light fluid regions retain the characteristic tear-drop shape. To understand this behavior and the return to standard turbulence structure as the flow evolves away from the shock, Lagrangian dynamics of the VGT and its invariants are studied by considering particle residence times and conditional particle variables in different flow regions. The pressure Hessian contributions to the VGT invariants transport equations are shown to be not only affected by the shock wave, but also by the density in the multi-fluid case, making them critically important to the flow dynamics and turbulence structure.
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.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
