-
Fluctuation-Induced Transitions in Anisotropic Two-Dimensional Turbulence
Authors:
Lichuan Xu,
Adrian van Kan,
Chang Liu,
Edgar Knobloch
Abstract:
Two-dimensional (2D) turbulence features an inverse energy cascade that produces large-scale flow structures such as large-scale vortices (LSVs) and unidirectional jets. We investigate the dynamics of such structures using extensive direct numerical simulations (DNS) of randomly forced, viscously damped 2D turbulence within a periodic rectangular (Cartesian) domain $[0,L_x]\times[0,L_y]$. LSVs for…
▽ More
Two-dimensional (2D) turbulence features an inverse energy cascade that produces large-scale flow structures such as large-scale vortices (LSVs) and unidirectional jets. We investigate the dynamics of such structures using extensive direct numerical simulations (DNS) of randomly forced, viscously damped 2D turbulence within a periodic rectangular (Cartesian) domain $[0,L_x]\times[0,L_y]$. LSVs form and dominate the system when the domain aspect ratio $δ= L_x/L_y \approx 1$, while unidirectional jets predominate at $δ\gtrsim 1.1$. At intermediate $δ$, both structures are metastable, with noise-induced transitions between LSVs and jets. We derive and verify predictions for the dependence of kinetic energy E and flow polarity on the nondimensional control parameters. We further collect detailed statistics on the lifetimes of LSVs and jets from DNS runs that are up to 10738 viscous diffusive times long. The distribution of the lifetimes is consistent with that of a memoryless process. Our DNS show an exponential dependence of the mean lifetime on $δ$. Mean lifetimes depend sensitively on the Reynolds number Re: as Re increases, the energy gap between LSV (lower E) and jet states (higher E) arising from anisotropic dissipation increases, leading to an approximately exponential increase in lifetimes with Re for both LSVs and jets. Similarly, as the forcing scale decreases, transitions become less frequent. We study the transitions in detail, revealing that they occur in two stages: an initial, rapid redistribution of kinetic energy by nonlinear triadic interactions deforms LSVs into jets or vice versa. In the second stage, the energy of the newly formed structure slowly adjusts to its associated equilibrium value on a longer, viscous timescale, producing hysteresis. Our findings shed new light on the dynamics of coherent large-scale structures in anisotropic turbulence.
△ Less
Submitted 23 April, 2024; v1 submitted 13 November, 2023;
originally announced November 2023.
-
From a vortex gas to a vortex crystal in instability-driven two-dimensional turbulence
Authors:
A. van Kan,
B. Favier,
K. Julien,
E. Knobloch
Abstract:
We study structure formation in two-dimensional turbulence driven by an external force, interpolating between linear instability forcing and random stirring, subject to nonlinear damping. Using extensive direct numerical simulations, we uncover a rich parameter space featuring four distinct branches of stationary solutions: large-scale vortices, hybrid states with embedded shielded vortices (SVs)…
▽ More
We study structure formation in two-dimensional turbulence driven by an external force, interpolating between linear instability forcing and random stirring, subject to nonlinear damping. Using extensive direct numerical simulations, we uncover a rich parameter space featuring four distinct branches of stationary solutions: large-scale vortices, hybrid states with embedded shielded vortices (SVs) of either sign, and two states composed of many similar SVs. Of the latter, the first is a dense vortex gas where all SVs have the same sign and diffuse across the domain. The second is a hexagonal vortex crystal forming from this gas when the instability is sufficiently weak. These solutions coexist stably over a wide parameter range. The late-time evolution of the system from small-amplitude initial conditions is nearly self-similar, involving three phases: initial inverse cascade, random nucleation of SVs from turbulence and, once a critical number of vortices is reached, a phase of explosive nucleation of SVs, leading to a statistically stationary state. The vortex gas is continued in the forcing parameter, revealing a sharp transition towards the crystal state as the forcing strength decreases. This transition is analysed in terms of the diffusion of individual vortices and tools from statistical physics. The crystal can also decay via an inverse cascade resulting from the breakdown of shielding or insufficient nonlinear damping acting on SVs. Our study highlights the importance of the forcing details in two-dimensional turbulence and reveals the presence of nontrivial SV states in this system, specifically the emergence and melting of a vortex crystal.
△ Less
Submitted 6 January, 2024; v1 submitted 17 August, 2023;
originally announced August 2023.
-
Universal wrinkling of supported elastic rings
Authors:
Benjamin Foster,
Nicolás Verschueren,
Edgar Knobloch,
Leonardo Gordillo
Abstract:
An exactly solvable family of models describing the wrinkling of substrate-supported inextensible elastic rings under compression is identified. The resulting wrinkle profiles are shown to be related to the buckled states of an unsupported ring and are therefore universal. Closed analytical expressions for the resulting universal shapes are provided, including the one-to-one relations between the…
▽ More
An exactly solvable family of models describing the wrinkling of substrate-supported inextensible elastic rings under compression is identified. The resulting wrinkle profiles are shown to be related to the buckled states of an unsupported ring and are therefore universal. Closed analytical expressions for the resulting universal shapes are provided, including the one-to-one relations between the pressure and tension at which these emerge. The analytical predictions agree with numerical continuation results to within numerical accuracy, for a large range of parameter values, up to the point of self-contact.
△ Less
Submitted 29 August, 2022;
originally announced August 2022.
-
Stationary broken parity states in active matter models
Authors:
Tobias Frohoff-Hülsmann,
Max Philipp Holl,
Edgar Knobloch,
Svetlana V. Gurevich,
Uwe Thiele
Abstract:
We demonstrate that several nonvariational continuum models commonly used to describe active matter as well as other active systems exhibit nongeneric behavior: each model supports asymmetric but stationary localized states even in the absence of pinning at heterogeneities. Moreover, such states only begin to drift following a drift-transcritical bifurcation as the activity increases. Asymmetric s…
▽ More
We demonstrate that several nonvariational continuum models commonly used to describe active matter as well as other active systems exhibit nongeneric behavior: each model supports asymmetric but stationary localized states even in the absence of pinning at heterogeneities. Moreover, such states only begin to drift following a drift-transcritical bifurcation as the activity increases. Asymmetric stationary states should only exist in variational systems, i.e., in models with gradient structure. In other words, such states are expected in passive systems, but not in active systems where the gradient structure of the model is broken by activity. We identify a "spurious" gradient dynamics structure of these models that is responsible for this nongeneric behavior, and determine the types of additional terms that render the models generic, i.e., with asymmetric states that appear via drift-pitchfork bifurcations and are generically moving. We provide detailed illustrations of our results using numerical continuation of resting and steadily drifting states in both generic and nongeneric cases.
△ Less
Submitted 19 June, 2023; v1 submitted 28 May, 2022;
originally announced May 2022.
-
Pressure-driven wrinkling of soft inner-lined tubes
Authors:
Benjamin Foster,
Nicolás Verschueren,
Edgar Knobloch,
Leonardo Gordillo
Abstract:
Wrinkling of an inextensible elastic lining of an inner-lined tube under imposed pressure is considered. A simple equation modeling the elastic properties of the lining, the pressure, and the soft-substrate forces is derived. This equation aims to capture the wrinkling response of arterial endothelium to blood pressure changes. Numerical continuation is used to construct a bifurcation diagram as a…
▽ More
Wrinkling of an inextensible elastic lining of an inner-lined tube under imposed pressure is considered. A simple equation modeling the elastic properties of the lining, the pressure, and the soft-substrate forces is derived. This equation aims to capture the wrinkling response of arterial endothelium to blood pressure changes. Numerical continuation is used to construct a bifurcation diagram as a function of the imposed pressure for in-plane deformations, in excellent agreement with weakly nonlinear theory, which we also develop. Our approach explains how the wavelength and amplitude of the wrinkles are selected as a function of the parameters in compressed wrinkling systems and shows how localized folds and mixed-mode states form in secondary bifurcations from wrinkled states.
△ Less
Submitted 18 January, 2022; v1 submitted 9 April, 2021;
originally announced April 2021.
-
Two-dimensional localized states in an active phase-field-crystal model
Authors:
Lukas Ophaus,
Edgar Knobloch,
Svetlana V. Gurevich,
Uwe Thiele
Abstract:
The active phase-field-crystal (active PFC) model provides a simple microscopic mean field description of crystallization in active systems. It combines the PFC model (or conserved Swift-Hohenberg equation) of colloidal crystallization and aspects of the Toner-Tu theory for self-propelled particles. We employ the active PFC model to study the occurrence of localized and periodic active crystals in…
▽ More
The active phase-field-crystal (active PFC) model provides a simple microscopic mean field description of crystallization in active systems. It combines the PFC model (or conserved Swift-Hohenberg equation) of colloidal crystallization and aspects of the Toner-Tu theory for self-propelled particles. We employ the active PFC model to study the occurrence of localized and periodic active crystals in two spatial dimensions. Due to the activity, crystalline states can undergo a drift instability and start to travel while keeping their spatial structure. Based on linear stability analyses, time simulations and numerical continuation of the fully nonlinear states, we present a detailed analysis of the bifurcation structure of resting and traveling states. We explore, for instance, how the slanted homoclinic snaking of steady localized states found for the passive PFC model is modified by activity. The analysis is carried out for the model in two spatial dimensions. Morphological phase diagrams showing the regions of existence of various solution types are presented merging the results from all the analysis tools employed. We also study how activity influences the crystal structure with transitions from hexagons to rhombic and stripe patterns. This in-depth analysis of a simple PFC model for active crystals and swarm formation provides a clear general understanding of the observed multistability and associated hysteresis effects, and identifies thresholds for qualitative changes in behavior.
△ Less
Submitted 10 December, 2020;
originally announced December 2020.
-
Localized states in passive and active phase-field-crystal models
Authors:
Max Philipp Holl,
Andrew J. Archer,
Svetlana V. Gurevich,
Edgar Knobloch,
Lukas Ophaus,
Uwe Thiele
Abstract:
The passive conserved Swift-Hohenberg equation (or phase-field-crystal [PFC] model) corresponds to a gradient dynamics for a single order parameter field related to density. It provides a simple microscopic description of the thermodynamic transition between liquid and crystalline states. In addition to spatially extended periodic structures, the model describes a large variety of steady spatially…
▽ More
The passive conserved Swift-Hohenberg equation (or phase-field-crystal [PFC] model) corresponds to a gradient dynamics for a single order parameter field related to density. It provides a simple microscopic description of the thermodynamic transition between liquid and crystalline states. In addition to spatially extended periodic structures, the model describes a large variety of steady spatially localized structures. In appropriate bifurcation diagrams the corresponding solution branches exhibit characteristic slanted homoclinic snaking. In an active PFC model, encoding for instance the active motion of self-propelled colloidal particles, the gradient dynamics structure is broken by a coupling between density and an additional polarization field. Then, resting and traveling localized states are found with transitions characterized by parity-breaking drift bifurcations. Here, we first briefly review the snaking behavior of localized states in passive and active PFC models before discussing the bifurcation behavior of localized states in systems of (i) two coupled passive PFC equations described by common gradient dynamics, (ii) two coupled passive PFC where the coupling breaks the gradient dynamics structure, and (iii) a passive PFC coupled to an active PFC.
△ Less
Submitted 6 October, 2020;
originally announced October 2020.
-
First order phase transitions and the thermodynamic limit
Authors:
Uwe Thiele,
Tobias Frohoff-Hülsmann,
Sebastian Engelnkemper,
Edgar Knobloch,
Andrew J. Archer
Abstract:
We consider simple mean field continuum models for first order liquid-liquid demixing and solid-liquid phase transitions and show how the Maxwell construction at phase coexistence emerges on going from finite-size closed systems to the thermodynamic limit. The theories considered are the Cahn-Hilliard model of phase separation, which is also a model for the liquid-gas transition, and the phase fie…
▽ More
We consider simple mean field continuum models for first order liquid-liquid demixing and solid-liquid phase transitions and show how the Maxwell construction at phase coexistence emerges on going from finite-size closed systems to the thermodynamic limit. The theories considered are the Cahn-Hilliard model of phase separation, which is also a model for the liquid-gas transition, and the phase field crystal model of the solid-liquid transition. Our results show that states comprising the Maxwell line depend strongly on the mean density with spatially localized structures playing a key role in the approach to the thermodynamic limit.
△ Less
Submitted 29 January, 2020; v1 submitted 29 August, 2019;
originally announced August 2019.
-
Three-dimensional Phase Field Quasicrystals
Authors:
P. Subramanian,
A. J. Archer,
E. Knobloch,
A. M. Rucklidge
Abstract:
We investigate the formation and stability of icosahedral quasicrytalline structures using a dynamic phase field crystal model. Nonlinear interactions between density waves at two length scales stabilize three-dimensional quasicrystals. We determine the phase diagram and parameter values required for the quasicrystal to be the global minimum free energy state. We demonstrate that traits that promo…
▽ More
We investigate the formation and stability of icosahedral quasicrytalline structures using a dynamic phase field crystal model. Nonlinear interactions between density waves at two length scales stabilize three-dimensional quasicrystals. We determine the phase diagram and parameter values required for the quasicrystal to be the global minimum free energy state. We demonstrate that traits that promote the formation of two-dimensional quasicrystals are extant in three dimensions, and highlight the characteristics required for 3D soft matter quasicrystal formation.
△ Less
Submitted 17 February, 2016;
originally announced February 2016.
-
Soft-core particles freezing to form a quasicrystal and a crystal-liquid phase
Authors:
A. J. Archer,
A. M. Rucklidge,
E. Knobloch
Abstract:
Systems of soft-core particles interacting via a two-scale potential are studied. The potential is responsible for peaks in the structure factor of the liquid state at two different but comparable length scales, and a similar bimodal structure is evident in the dispersion relation. Dynamical density functional theory in two dimensions is used to identify two novel states of this system, the crysta…
▽ More
Systems of soft-core particles interacting via a two-scale potential are studied. The potential is responsible for peaks in the structure factor of the liquid state at two different but comparable length scales, and a similar bimodal structure is evident in the dispersion relation. Dynamical density functional theory in two dimensions is used to identify two novel states of this system, the crystal-liquid state, in which the majority of the particles are located on lattice sites but a minority remains free and so behaves like a liquid, and a 12-fold quasicrystalline state. Both are present even for deeply quenched liquids and are found in a regime in which the liquid is unstable with respect to modulations on the smaller scale only. As a result the system initially evolves towards a small scale crystal state; this state is not a minimum of the free energy, however, and so the system subsequently attempts to reorganize to generate the lower energy larger scale crystals. This dynamical process generates a disordered state with quasicrystalline domains, and takes place even when this large scale is linearly stable, i.e., it is a nonlinear process. With controlled initial conditions a perfect quasicrystal can form. The results are corroborated using Brownian dynamics simulations.
△ Less
Submitted 7 July, 2015;
originally announced July 2015.
-
Generation of defects and disorder from deeply quenching a liquid to form a solid
Authors:
A. J. Archer,
M. C. Walters,
U. Thiele,
E. Knobloch
Abstract:
We show how deeply quenching a liquid to temperatures where it is linearly unstable and the crystal is the equilibrium phase often produces crystalline structures with defects and disorder. As the solid phase advances into the liquid phase, the modulations in the density distribution created behind the advancing solidification front do not necessarily have a wavelength that is the same as the equi…
▽ More
We show how deeply quenching a liquid to temperatures where it is linearly unstable and the crystal is the equilibrium phase often produces crystalline structures with defects and disorder. As the solid phase advances into the liquid phase, the modulations in the density distribution created behind the advancing solidification front do not necessarily have a wavelength that is the same as the equilibrium crystal lattice spacing. This is because in a deep enough quench the front propagation is governed by linear processes, but the crystal lattice spacing is determined by nonlinear terms. The wavelength mismatch can result in significant disorder behind the front that may or may not persist in the latter stage dynamics. We support these observations by presenting results from dynamical density functional theory calculations for simple one- and two-component two-dimensional systems of soft core particles.
△ Less
Submitted 29 May, 2015;
originally announced May 2015.
-
Solidification in soft-core fluids: disordered solids from fast solidification fronts
Authors:
A. J. Archer,
M. C. Walters,
U. Thiele,
E. Knobloch
Abstract:
Using dynamical density functional theory we calculate the speed of solidification fronts advancing into a quenched two-dimensional model fluid of soft-core particles. We find that solidification fronts can advance via two different mechanisms, depending on the depth of the quench. For shallow quenches, the front propagation is via a nonlinear mechanism. For deep quenches, front propagation is gov…
▽ More
Using dynamical density functional theory we calculate the speed of solidification fronts advancing into a quenched two-dimensional model fluid of soft-core particles. We find that solidification fronts can advance via two different mechanisms, depending on the depth of the quench. For shallow quenches, the front propagation is via a nonlinear mechanism. For deep quenches, front propagation is governed by a linear mechanism and in this regime we are able to determine the front speed via a marginal stability analysis. We find that the density modulations generated behind the advancing front have a characteristic scale that differs from the wavelength of the density modulation in thermodynamic equilibrium, i.e., the spacing between the crystal planes in an equilibrium crystal. This leads to the subsequent development of disorder in the solids that are formed. For the one-component fluid, the particles are able to rearrange to form a well-ordered crystal, with few defects. However, solidification fronts in a binary mixture exhibiting crystalline phases with square and hexagonal ordering generate solids that are unable to rearrange after the passage of the solidification front and a significant amount of disorder remains in the system.
△ Less
Submitted 29 September, 2014; v1 submitted 26 July, 2014;
originally announced July 2014.
-
Quasicrystalline order and a crystal-liquid state in a soft-core fluid
Authors:
A. J. Archer,
A. M. Rucklidge,
E. Knobloch
Abstract:
A two-dimensional system of soft particles interacting via a two-length-scale potential is studied. Density functional theory and Brownian dynamics simulations reveal a fluid phase and two crystalline phases with different lattice spacing. Of these the larger lattice spacing phase can form an exotic periodic state with a fraction of highly mobile particles: a crystal liquid. Near the transition be…
▽ More
A two-dimensional system of soft particles interacting via a two-length-scale potential is studied. Density functional theory and Brownian dynamics simulations reveal a fluid phase and two crystalline phases with different lattice spacing. Of these the larger lattice spacing phase can form an exotic periodic state with a fraction of highly mobile particles: a crystal liquid. Near the transition between this phase and the smaller lattice spacing phase, quasicrystalline structures may be created by a competition between linear instability at one scale and nonlinear selection of the other.
△ Less
Submitted 18 September, 2013; v1 submitted 8 April, 2013;
originally announced April 2013.
-
Localized states in the conserved Swift-Hohenberg equation with cubic nonlinearity
Authors:
Uwe Thiele,
Andrew J. Archer,
Mark J. Robbins,
Hector Gomez,
Edgar Knobloch
Abstract:
The conserved Swift-Hohenberg equation with cubic nonlinearity provides the simplest microscopic description of the thermodynamic transition from a fluid state to a crystalline state. The resulting phase field crystal model describes a variety of spatially localized structures, in addition to different spatially extended periodic structures. The location of these structures in the temperature vers…
▽ More
The conserved Swift-Hohenberg equation with cubic nonlinearity provides the simplest microscopic description of the thermodynamic transition from a fluid state to a crystalline state. The resulting phase field crystal model describes a variety of spatially localized structures, in addition to different spatially extended periodic structures. The location of these structures in the temperature versus mean order parameter plane is determined using a combination of numerical continuation in one dimension and direct numerical simulation in two and three dimensions. Localized states are found in the region of thermodynamic coexistence between the homogeneous and structured phases, and may lie outside of the binodal for these states. The results are related to the phenomenon of slanted snaking but take the form of standard homoclinic snaking when the mean order parameter is plotted as a function of the chemical potential, and are expected to carry over to related models with a conserved order parameter.
△ Less
Submitted 18 January, 2013;
originally announced January 2013.
-
Solidification fronts in supercooled liquids: how rapid fronts can lead to disordered glassy solids
Authors:
A. J. Archer,
M. J. Robbins,
U. Thiele,
E. Knobloch
Abstract:
We determine the speed of a crystallisation (or more generally, a solidification) front as it advances into the uniform liquid phase after the system has been quenched into the crystalline region of the phase diagram. We calculate the front speed by assuming a dynamical density functional theory model for the system and applying a marginal stability criterion. Our results also apply to phase field…
▽ More
We determine the speed of a crystallisation (or more generally, a solidification) front as it advances into the uniform liquid phase after the system has been quenched into the crystalline region of the phase diagram. We calculate the front speed by assuming a dynamical density functional theory model for the system and applying a marginal stability criterion. Our results also apply to phase field crystal (PFC) models of solidification. As the solidification front advances into the unstable liquid phase, the density profile behind the advancing front develops density modulations and the wavelength of these modulations is a dynamically chosen quantity. For shallow quenches, the selected wavelength is precisely that of the crystalline phase and so well-ordered crystalline states are formed. However, when the system is deeply quenched, we find that this wavelength can be quite different from that of the crystal, so that the solidification front naturally generates disorder in the system. Significant rearrangement and ageing must subsequently occur for the system to form the regular well-ordered crystal that corresponds to the free energy minimum. Additional disorder is introduced whenever a front develops from random initial conditions. We illustrate these findings with results obtained from the PFC.
△ Less
Submitted 19 September, 2012; v1 submitted 5 June, 2012;
originally announced June 2012.
-
Modelling fluids and crystals using a two-component modified phase field crystal model
Authors:
Mark J. Robbins,
Andrew J. Archer,
Uwe Thiele,
Edgar Knobloch
Abstract:
A modified phase field crystal model in which the free energy may be minimised by an order parameter profile having isolated bumps is investigated. The phase diagram is calculated in one and two dimensions and we locate the regions where modulated and uniform phases are formed and also regions where localised states are formed. We investigate the effectiveness of the phase field crystal model for…
▽ More
A modified phase field crystal model in which the free energy may be minimised by an order parameter profile having isolated bumps is investigated. The phase diagram is calculated in one and two dimensions and we locate the regions where modulated and uniform phases are formed and also regions where localised states are formed. We investigate the effectiveness of the phase field crystal model for describing fluids and crystals with defects. We further consider a two component model and elucidate how the structure transforms from hexagonal crystalline ordering to square ordering as the concentration changes. Our conclusion contains a discussion of possible interpretations of the order parameter field.
△ Less
Submitted 9 December, 2011;
originally announced December 2011.