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.

Classical molecular dynamics (MD) simulations were used to investigate how free surfaces, as well as supporting substrates, affect phase separation in a NiAg alloy. Bulk samples, droplets, and droplets deposited on a graphene substrate were investigated at temperatures that spanned regions of interest in the bulk NiAg phase diagram, i.e., miscible and immiscible liquid, liquid-crystal, and crystal-crystal regions. Using MD simulations to cool down a bulk sample from 3000 K to 800 K, it was found that phase separation below 2400 K takes place in agreement with the phase diagram. When free surface effects were introduced, phase separation was accompanied by a core-shell transformation: spherical droplets created from the bulk samples became core-shell nanoparticles with a shell made mostly of Ag atoms and a core made of Ni atoms. When such droplets were deposited on a graphene substrate, the phase separation was accompanied by Ni layering at the graphene interface and Ag at the vacuum interface. Thus, it should be possible to create NiAg core-shell and layer-like nanostructures by quenching liquid NiAg samples on tailored substrates. Furthermore, interesting bimetallic nanoparticle morphologies might be tuned via control of the surface and interface energies and chemical instabilities of the system.

We consider a free-surface thin film placed on a thermally conductive substrate and exposed to an external heat source in a set-up where the heat absorption depends on the local film thickness. Our focus is on modelling 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 film's 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 set-up 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 set-up the Marangoni (thermocapillary) effect turns out to be insignificant.

Metal films and other metal geometries of nanoscale thickness deposited on an insulating substrate, when exposed to laser irradiation, melt and evolve as fluids as long as their temperature is sufficiently high. This evolution often leads to pattern formation, which may be influenced strongly by material parameters that are temperature dependent. In addition, the laser heat absorption itself depends on the time-dependent metal thickness. Self-consistent modeling of evolving metal films shows that, by controlling the amount and geometry of the deposited metal, one can control the instability development. In particular, we demonstrate the "thermal crowding" effect: additional metal leads to elevated temperatures, which strongly influence the metal evolution, even if the metal geometries are disjoint. We demonstrate that the communication of disjoint metal domains occurs via heat diffusion through the underlying substrate. Fully self-consistent modeling focusing on the dominant effects, as well as accurate time-dependent simulations, allow us to describe the main features of thermal crowding and provide a route to control fluid instabilities and pattern formation on the nanoscale.

We extend the Lanchester equations of combat modeling by constructing a coupled system of partial differential equations with the goal of representing engagements between an attacking drone swarm and a layered defense. Our model includes an attacking force of drones traversing the battlefield (domain) in an attempt to reach its target. The defending force, which engages the incoming attackers, is made up of a static defense and a mobile drone force capable of intercepting and tracking the attackers. Drone swarm motion is modeled with advection–diffusion terms and engagement at range by nonlocal reaction terms. We solve the underlying model using a finite difference method and investigate how interaction range, swarm dispersion, and the allocation between static and mobile defenses affect attacker survivability. To account for expendable or single-use countermeasure platforms (e.g. kamikaze–style interceptors), we also examine how expendability affects the final engagement outcome. This study provides a foundation for modeling spatially distributed combat dynamics reflecting current and future combat scenarios.

In this paper we discuss the deformation of a flexible fiber clamped to a spherical body and immersed in a flow of fluid moving with a speed ranging between 0 and 50 cm/s by means of three dimensional numerical simulation developed in COMSOL . The effects of flow speed and initial configuration angle of the fiber relative to the flow are analyzed. A rigorous analysis of the numerical procedure is performed and our code is benchmarked against well established cases. The flow velocity and pressure are used to compute drag forces upon the fiber. Of particular interest is the behavior of these forces as a function of the flow speed and fiber orientation. The Vogel exponents, which characterize the rate of bending of a fiber in a flow, are found for the various configurations examined here and seem to display interesting variations. These exponents are then compared with our previously studied two-dimensional models.