In this paper, we propose a microscale liquid oscillator using electrowetting-on-dielectric (EWOD). Specifically, a mesoscale liquid bridge (LB) between two horizontal surfaces with EWOD is considered. When EWOD is applied, the solid surface becomes more hydrophilic, and hence the contact angle (CA) is reduced. Following the activation of EWOD, the LB can remain connected or it can break into either symmetric or asymmetric shapes depending on the initial liquid volume and wettability of the two surfaces. The LB dynamics activated by EWOD is studied using the multibody dissipative particle dynamics (MDPD) method. Our numerical results show that the behavior of an LB under EWOD can be interpreted via three modes. In the first mode, the LB does not break after applying EWOD. In the second mode, the LB breaks and does not reform. The third mode happens when, depending on the interplay of the volume of the liquid and CA manipulation, the LB continuously breaks, recoils, and reforms. For asymmetric cases, it was observed that the LB may completely detach from one surface and may not reform. It was also observed that decreasing the wettability of a surface, for cases with a continuous breaking and reformation behavior, increases the connecting time interval and decreases the breaking time interval in one break-reform cycle. The results provided in this investigation facilitate fundamental understanding of LB dynamics and their application for the design of microscale liquid oscillators using EWOD.

An explicit analytical solution to calculate the profiles after the shock collision with a planar contact surface is presented. The case when a shock is reflected after the incident shock refraction is considered. The goal of this work is to present explicit formulas to obtain the quantities behind the transmitted and reflected shocks valid for arbitrary initial preshock parameters.

Viscosity independence of lattice-Boltzmann methods is a crucial issue to ensure the physical relevancy of the predicted macroscopic flows over large ranges of physical parameters. The immersed-boundary (IB) method, a powerful tool that allows one to immerse arbitrary-shaped, moving, and deformable bodies in the flow, suffers from a boundary-slip error that increases as a function of the fluid viscosity, substantially limiting its range of application. In addition, low fluid viscosities may result in spurious oscillations of the macroscopic quantities in the vicinity of the immersed boundary. In this work, it is shown mathematically that the standard IB method is indeed not able to reproduce the scaling properties of the macroscopic solution, leading to a viscosity-related error on the computed IB force. The analysis allows us to propose a simple correction of the IB scheme that is local, straightforward and does not involve additional computational time. The derived method is implemented in a two-relaxation-time D2Q9 lattice-Boltzmann solver, applied to several physical configurations, namely, the Poiseuille flow, the flow around a cylinder towed in still fluid, and the flow around a cylinder oscillating in still fluid, and compared to a noncorrected immersed-boundary method. The proposed correction leads to a major improvement of the viscosity independence of the solver over a wide range of relaxation times (from 0.5001 to 50), including the correction of the boundary-slip error and the suppression of the spurious oscillations. This improvement may considerably extend the range of application of the IB lattice-Boltzmann method, in particular providing a robust tool for the numerical analysis of physical problems involving fluids of varying viscosity interacting with solid geometries.

We study the statistics of large deviations of the intensive work done in an interaction quench of a one-dimensional Bose gas with a large number $N$ of particles, system size $L$, and fixed density. We consider the case in which the system is initially prepared in the noninteracting ground state and a repulsive interaction is suddenly turned on. For large deviations of the work below its mean value, we show that the large-deviation principle holds by means of the quench action approach. Using the latter, we compute exactly the so-called rate function and study its properties analytically. In particular, we find that fluctuations close to the mean value of the work exhibit a marked non-Gaussian behavior, even though their probability is always exponentially suppressed below it as $L$ increases. Deviations larger than the mean value exhibit an algebraic decay whose exponent cannot be determined directly by large-deviation theory. Exploiting the exact Bethe ansatz representation of the eigenstates of the Hamiltonian, we calculate this exponent for vanishing particle density. Our approach can be straightforwardly generalized to quantum quenches in other interacting integrable systems.

We investigate the effect of dissipation in a Bose-Josephson junction (BJJ) coupled to baths of bosons at its two sites. Apart from the dynamical transition due to repulsive interactions, the BJJ undergoes a quantum phase transition by increasing the coupling strength with the bath modes. We analyze this system by mapping to an equivalent spin model coupled to the bosonic modes. The excitation energies and fluctuation of number imbalance are obtained within a Holstein-Primakoff approximation, which exhibits vanishing of the energy gap and enhanced quantum fluctuations at the critical point. We study the dynamics of BJJ using a time-dependent variational method and analyze stability of different types of steady states. As a special case we study in detail the phase space dynamics of BJJ coupled to a single mode, which reveals diffusive and incoherent behavior with increasing coupling to the bath mode. The dynamical steady states corresponding to the π oscillation and self-trapped state become unstable when their oscillation frequencies are in resonance with the bath modes. We study the Josephson dynamics in the presence of an Ohmic bath with Gaussian noise to incorporate the thermal fluctuations and obtain the Josephson oscillation frequency and damping analytically. We also observe the transition to the symmetry-broken state for strong coupling as well as decay of π oscillation and a self-trapped state to the ground state due to dissipation. Variation of the phase fluctuation with temperature of the bath shows similar behavior as observed in experiment. Finally we discuss the experimental setup to study the observable effects of dissipation in BJJ.

The autoresonant approach to excitation and control of large-amplitude uniformly precessing magnetization structures in finite-length easy axis ferromagnetic nanoparticles is suggested and analyzed within the Landau-Lifshitz-Gilbert model. These structures are excited by using a spatially uniform, oscillating, chirped frequency magnetic field, while the localization is imposed via boundary conditions. The excitation requires the amplitude of the driving oscillations to exceed a threshold. The dissipation effect on the threshold is also discussed. The autoresonant driving effectively compensates the effect of dissipation but lowers the maximum amplitude of the excited structures. Fully nonlinear localized autoresonant solutions are illustrated in simulations and described via an analog of a quasiparticle in an effective potential. The precession frequency of these solutions is continuously locked to that of the drive, while the spatial magnetization profile approaches the soliton limit when the length of the nanoparticle and the amplitude of the excited solution increase.

Absence epilepsy is characterized by a sudden paroxysmal loss of consciousness accompanied by oscillatory activity propagating over many brain areas. Although primary generalized absence seizures are supported by the global corticothalamic system, converging experimental evidence supports a focal theory of absence epilepsy. Here a physiology-based corticothalamic model is investigated with spatial heterogeneity due to focal epilepsy to unify global and focal aspects of absence epilepsy. Numeric and analytic calculations are employed to investigate the emergent spatiotemporal dynamics as well as their underlying dynamical mechanisms. They can be categorized into three scenarios: suppressed epilepsy, focal seizures, or generalized seizures, as summarized from a phase diagram vs focal width and characteristic axon range. The corresponding temporal frequencies and spatial extents of cortical waves in generalized seizures and focal seizures agree well with experimental observations of global and focal aspects of absence epilepsy, respectively. The emergence of the spatiotemporal dynamics corresponding to focal seizures provides a biophysical explanation of the temporally higher frequency but spatially more localized cortical waves observed in genetic rat models that display characteristics of human absence epilepsy. Predictions are also presented for further experimental test.

One prototypical instability in granular flows is the shear-banding instability, in which a uniform granular shear flow breaks into alternating bands of dense and dilute clusters of particles having low and high shear (shear stress or shear rate), respectively. In this work, the shear-banding instability in an arbitrarily inelastic granular shear flow is analyzed through the linear stability analysis of granular hydrodynamic equations closed with Navier–Stokes-level constitutive relations. It is shown that the choice of appropriate constitutive relations plays an important role in predicting the shear-banding instability. A parametric study is carried out to study the effect of the restitution coefficient, channel width, and mean density. Two global criteria relating the control parameters are found for the onset of the shear-banding instability.

A theory is developed that allows one to model the velocity field of acoustic microstreaming produced by nonspherical oscillations of an acoustically driven gas bubble. It is assumed that some of the bubble oscillation modes are excited parametrically and hence can oscillate at frequencies different from the driving frequency. Analytical solutions are derived in terms of complex amplitudes of oscillation modes, which means that the mode amplitudes are assumed to be known and serve as input data when the velocity field of acoustic microstreaming is calculated. No restrictions are imposed on the ratio of the bubble radius to the viscous penetration depth. The present paper is the first part of our study in which a general theory is developed and then applied to the case that acoustic microstreaming is generated by the interaction of the breathing mode (mode 0) with a mode of arbitrary order $m\ge 1$. Examples of numerical simulations and a comparison with experimental results are provided.

This paper continues a study that was started in our previous paper [A. A. Doinikov et al., Phys. Rev. E 100, 033104 (2019)]. The overall aim of the study is to develop a theory for modeling the velocity field of acoustic microstreaming produced by nonspherical oscillations of an acoustically driven gas bubble. In the previous paper, general equations were derived that describe the velocity field of acoustic microstreaming produced by modes $n$ and $m$ of bubble oscillations. In the present paper, the above equations are solved analytically in the case that acoustic microstreaming is the result of the interaction of the translational mode (mode 1) with a mode of arbitrary order $m\ge 1$. Solutions are expressed in terms of complex mode amplitudes, which means that the mode amplitudes are assumed to be known and serve as input data for the calculation of the velocity field of acoustic microstreaming. No restrictions are imposed on the ratio of the bubble radius to the viscous penetration depth. Analytical results are illustrated by numerical examples.

Magnetically driven artificial microswimmers have the potential to revolutionize many biomedical technologies, such as minimally invasive microsurgery, microparticle manipulation, and localized drug delivery. However, many of these applications will require the controlled dynamics of teams of these microrobots with minimal feedback. In this work, we study the motion and fluid dynamics produced by groups of artificial microswimmers driven by a torque induced through a uniform, rotating magnetic field. Through Stokesian dynamics simulations, we show that the swimmer motion produces a rotational velocity field in the plane orthogonal to the direction of the magnetic field's rotation, which causes two interacting swimmers to move in circular trajectories in this plane around a common center. The resulting overall motion is on a helical trajectory for the swimmers. We compare the highly rotational velocity field of the fluid to the velocity field generated by a rotlet, the point-torque singularity of Stokes flows, showing that this is a reasonable approximation on the time average. Finally, we study the motion of larger groups of swimmers, and we show that these groups tend to move coherently, especially when swimmer magnetizations are uniform. This coherence is achieved because the group center remains almost constant in the plane orthogonal to the net motion of the swimmers. The results in the paper will prove useful for controlling the ensemble dynamics of small collections of magnetic swimmers.

A simplified model is developed, which allows us to perform computer simulations of the particles transport in an evaporating droplet with a contact line pinned to a hydrophilic substrate. The model accounts for advection in the droplet, diffusion, and particle attraction by capillary forces. On the basis of the simulations, we analyze the physical mechanisms of forming of individual chains of particles inside the annular sediment. The parameters chosen correspond to the experiments of Park and Moon [Langmuir 22, 3506 (2006)], where an annular deposition and snakelike chains of colloid particles have been identified. The annular sediment is formed by advection and diffusion transport. We find that the close packing of the particles in the sediment is possible if the evaporation time exceeds the characteristic time of diffusion-based ordering. We show that the chains are formed by the end of the evaporation process due to capillary attraction of particles in the region bounded by a fixing radius, where the local droplet height is comparable to the particle size. At the beginning of the evaporation, the annular deposition is shown to expand faster than the fixing radius moves. However, by the end of the process, the fixing radius rapidly outreaches the expanding inner front of the ring. The snakelike chains are formed at this final stage when the fixing radius moves toward the symmetry axis.

Over the last decades, several types of collision models have been proposed to extend the validity domain of the lattice Boltzmann method (LBM), each of them being introduced in its own formalism. This article proposes a formalism that describes all these methods within a common mathematical framework, and in this way allows us to draw direct links between them. Here, the focus is put on single and multirelaxation time collision models in either their raw moment, central moment, cumulant, or regularized form. In parallel with that, several bases (nonorthogonal, orthogonal, Hermite) are considered for the polynomial expansion of populations. General relationships between moments are first derived to understand how moment spaces are related to each other. In addition, a review of collision models further sheds light on collision models that can be rewritten in a linear matrix form. More quantitative mathematical studies are then carried out by comparing explicit expressions for the post-collision populations. Thanks to this, it is possible to deduce the impact of both the polynomial basis (raw, Hermite, central, central Hermite, cumulant) and the inclusion of regularization steps on isothermal LBMs. Extensive results are provided for the D1Q3, D2Q9, and D3Q27 lattices, the latter being further extended to the D3Q19 velocity discretization. Links with the most common two and multirelaxation time collision models are also provided for the sake of completeness. This work ends by emphasizing the importance of an accurate representation of the equilibrium state, independently of the choice of moment space. As an addition to the theoretical purpose of this article, general instructions are provided to help the reader with the implementation of the most complicated collision models.

We study binary mixtures of hard particles, which exclude up to their $k\mathrm{th}$ nearest neighbors $\left(k\mathrm{NN}\right)$ on the simple cubic lattice and have activities $z_k$. In the first model analyzed, point particles (0NN) are mixed with 1NN ones. The grand-canonical solution of this model on a Husimi lattice built with cubes unveils a phase diagram with a fluid and a solid phase separated by a continuous and a discontinuous transition line which meet at a tricritical point. A density anomaly, characterized by minima in isobaric curves of the total density of particles against $z_0$ (or $z_1)$, is also observed in this system. Overall, this scenario is identical to the one previously found for this model when defined on the square lattice. The second model investigated consists of the mixture of 1NN particles with 2NN ones. In this case, a very rich phase behavior is found in its Husimi lattice solution, with two solid phases—one associated with the ordering of 1NN particles $\left(S1\right)$ and the other with the ordering of 2NN ones $\left(S2\right)$—beyond the fluid $\left(F\right)$ phase. While the transitions between $F\text{−}S2$ and $S1\text{−}S2$ phases are always discontinuous, the $F\text{−}S1$ transition is continuous (discontinuous) for small (large) $z_2$. The critical and coexistence $F\text{−}S1$ lines meet at a tricritical point. Moreover, the coexistence $F\text{−}S1,F\text{−}S2$, and $S1\text{−}S2$ lines meet at a triple point. Density anomalies are absent in this case.

In the present study of the nonequilibrium athermal random-field Ising model we focus on the behavior of the critical disorder $R_c\left(l\right)$ and the critical magnetic field $H_c\left(l\right)$ under different boundary conditions when the system thickness $l$ varies. We propose expressions for $R_c\left(l\right)$ and $H_c\left(l\right)$ as well as for the effective critical disorder $R_c^{\mathrm{eff}}(l,L)$ and effective critical magnetic field $H_c^{\mathrm{eff}}(l,L)$ playing the role of the effective critical parameters for the $L\times L\times l$ lattices of finite lateral size $L$. We support these expressions by the scaling collapses of the magnetization and susceptibility curves obtained in extensive simulations. The collapses are achieved with the two-dimensional (2D) exponents for $l$ below some characteristic value, providing thus a numerical evidence that the thin systems exhibit a 2D-like criticality which should be relevant for the experimental analyses of thin ferromagnetic samples.

We investigate the collective response of an ensemble of bistable oscillators to an external periodic signal, where the coupling strength between oscillators is diverse. We find that there exists an optimal level of coupling diversity, at which the collective response of the system can be largely improved, i.e., resonance induced by coupling diversity. We also observe that the system splits into three oscillation clusters when this resonance happens. We finally propose a reduced model based on the three oscillation clusters, which can well predict the collective response of the system.

A framework for data assimilation combining aspects of operator-theoretic ergodic theory and quantum mechanics is developed. This framework adapts the Dirac–von Neumann formalism of quantum dynamics and measurement to perform sequential data assimilation (filtering) of a partially observed, measure-preserving dynamical system, using the Koopman operator on the $L^2$ space associated with the invariant measure as an analog of the Heisenberg evolution operator in quantum mechanics. In addition, the state of the data assimilation system is represented by a trace-class operator analogous to the quantum mechanical density operator, and the assimilated observables by self-adjoint multiplication operators. An averaging approach is also introduced, rendering the spectrum of the assimilated observables discrete and thus amenable to numerical approximation. We present a data-driven formulation of the quantum mechanical data assimilation approach, utilizing kernel methods from machine learning and delay-coordinate maps of dynamical systems to represent the evolution and measurement operators via matrices in a data-driven basis. The data-driven formulation is structurally similar to its infinite-dimensional counterpart and shown to converge in a limit of large data under mild assumptions. Applications to periodic oscillators and the Lorenz 63 system demonstrate that the framework is able to naturally handle highly non-Gaussian statistics, complex state space geometries, and chaotic dynamics.

The dynamic susceptibility of concentrated ferrofluids of magnetite-kerosene type is studied experimentally to clarify the effect of interparticle interactions on the magnetization reversal dynamics and the ferrofluid relaxation time spectrum. We synthesize six ferrofluid samples, four of which have the same wide particle size distribution with a high (more than $2kT)$ average energy of magnetic dipole interactions. These samples differ in particle concentration and dynamic viscosity. The two remaining samples have a lower content of large particles and a moderate energy of magnetic dipole interactions. For all samples, we measure the dynamic susceptibility in the weak probing field at frequencies up to 160 kHz and the field amplitude dependence of the susceptibility at a frequency of 27 kHz. The results show that the susceptibility dispersion at frequencies up to 10 kHz is due to the rotational diffusion of colloidal particles and aggregates. Steric and hydrodynamic interparticle interactions are the main reason for the strong concentration dependence of the viscosity and so they also strongly influence the frequency dependence of the susceptibility. The influence of van der Waals and magnetic dipole interactions on the susceptibility is manifested indirectly, through the formation of multiparticle clusters. The contribution of clusters to the low-frequency susceptibility reaches 80%. Their large sizes (about 100 nm) shift the dispersion region to frequencies of 1–100 Hz, depending on the temperature and particle concentration. Experiments at 27 kHz demonstrate the increase in the dynamic susceptibility with increasing field amplitude. This growth is unexpected since all spectral amplitudes in the Debye function expansion of the dynamic susceptibility decrease monotonically with increasing field. To clarify the situation, the auxiliary problem of the magnetodynamics of a uniaxial particle in the alternating field is solved numerically. The Fokker-Planck-Brown rotational diffusion equation is used. It is shown that an increase in the field amplitude reduces the anisotropy barrier and the Néel relaxation time of particles and increases the dynamic susceptibility by one to two orders of magnitude compared to the weak-field limit. The calculation results are in qualitative agreement with the experimental data and allow us to propose a consistent interpretation of these data. We find that the increase in dynamic susceptibility with increasing amplitude is observed when two necessary conditions are met: (i) The suspension viscosity and the field frequency are high enough to cause the blocking of the rotational degrees of freedom of particles and aggregates and (ii) particles with a large magnetic anisotropy are present in the ferrofluid.

We investigate diffusion of colloidal particles carried by flow in tubes of variable diameter and under the influence of an external field. We generalize the method mapping the three-dimensional confined diffusion onto an effective one-dimensional problem to the case of nonconservative forces and use this mapping for the problem in question. We show that in the presence of hydrodynamic drag, the lowest approximation (the Fick-Jacobs approximation) may be insufficient, and inclusion of at least the first-order correction is desirable to obtain more reliable results. As a practical application, we use the method for investigation of separation of colloidal particles carried by a fluid flow according to their size, using flotation and centrifugation.

Inertial confinement fusion fuel suffers increased x-ray radiation losses when carbon from the capsule ablator mixes into the hot-spot. Here, we present one- and two-dimensional ion Vlasov-Fokker-Planck simulations that resolve hot-spot self-heating in the presence of a localized spike of carbon mix, totalling $1.9\%$ of the hot-spot mass. The mix region cools and contracts over tens of picoseconds, increasing its $\alpha$ particle stopping power and radiative losses. This makes a localized mix region more severe than an equal amount of uniformly distributed mix. There is also a purely kinetic effect that reduces fusion reactivity by several percent, since faster ions in the tail of the distribution are absorbed by the mix region. Radiative cooling and contraction of the spike induces fluid motion, causing neutron spectrum broadening. This artificially increases the inferred experimental ion temperatures and gives line of sight variations.