skip to main content

Title: Phase-field model for a weakly compressible soft layered material: morphological transitions on smectic–isotropic interfaces
A coupled phase-field and hydrodynamic model is introduced to describe a two-phase, weakly compressible smectic (layered phase) in contact with an isotropic fluid of different density. A non-conserved smectic order parameter is coupled to a conserved mass density in order to accommodate non-solenoidal flows near the smectic–isotropic boundary arising from density contrast between the two phases. The model aims to describe morphological transitions in smectic thin films under heat treatment, in which arrays of focal conic defects evolve into conical pyramids and concentric rings through curvature dependent evaporation of smectic layers. The model leads to an extended thermodynamic relation at a curved surface that includes its Gaussian curvature, non-classical stresses at the boundary and flows arising from density gradients. The temporal evolution given by the model conserves the overall mass of the liquid crystal while still allowing for the modulated smectic structure to grow or shrink. A numerical solution of the governing equations reveals that pyramidal domains are sculpted at the center of focal conics upon a temperature increase, which display tangential flows at their surface. Other cases investigated include the possible coalescence of two cylindrical stacks of smectic layers, formation of droplets, and the interactions between focal conic domains more » through flow. « less
; ;
Award ID(s):
Publication Date:
Journal Name:
Soft Matter
Page Range or eLocation-ID:
6140 to 6159
Sponsoring Org:
National Science Foundation
More Like this
  1. The present work uses a reduced-order model to study the motion of a buoyant vortex ring with non-negligible core size. Buoyancy is considered in both non-Boussinesq and Boussinesq situations using an axisymmetric contour dynamics formulation. The density of the vortex ring differs from that of the ambient fluid, and both densities are constant and conserved. The motion of the ring is calculated by following the boundary of the vortex core, which is also the interface between the two densities. The velocity of the contour comes from a combination of a specific continuous vorticity distribution within its core and a vortex sheet on the core boundary. An evolution equation for the vortex sheet is derived from the Euler equation, which simplifies considerably in the Boussinesq limit. Numerical solutions for the coupled integro-differential equations are obtained. The dynamics of the vortex sheet and the formation of two possible singularities, including singularities in the curvature and the shock-like profile of the vortex sheet strength, are discussed. Three dimensionless groups, the Atwood, Froude and Weber numbers, are introduced to measure the importance of physical effects acting on the motion of a buoyant vortex ring.
  2. Inspired by the recent realization of a two-dimensional (2-D) chiral fluid as an active monolayer droplet moving atop a 3-D Stokesian fluid, we formulate mathematically its free-boundary dynamics. The surface droplet is described as a general 2-D linear, incompressible and isotropic fluid, having a viscous shear stress, an active chiral driving stress and a Hall stress allowed by the lack of time-reversal symmetry. The droplet interacts with itself through its driven internal mechanics and by driving flows in the underlying 3-D Stokes phase. We pose the dynamics as the solution to a singular integral–differential equation, over the droplet surface, using the mapping from surface stress to surface velocity for the 3-D Stokes equations. Specializing to the case of axisymmetric droplets, exact representations for the chiral surface flow are given in terms of solutions to a singular integral equation, solved using both analytical and numerical techniques. For a disc-shaped monolayer, we additionally employ a semi-analytical solution that hinges on an orthogonal basis of Bessel functions and allows for efficient computation of the monolayer velocity field, which ranges from a nearly solid-body rotation to a unidirectional edge current, depending on the subphase depth and the Saffman–Delbrück length. Except in the near-wall limit,more »these solutions have divergent surface shear stresses at droplet boundaries, a signature of systems with codimension-one domains embedded in a 3-D medium. We further investigate the effect of a Hall viscosity, which couples radial and transverse surface velocity components, on the dynamics of a closing cavity. Hall stresses are seen to drive inward radial motion, even in the absence of edge tension.« less
  3. Abstract
    This dataset contains monthly average output files from the iCAM6 simulations used in the manuscript "Enhancing understanding of the hydrological cycle via pairing of process-oriented and isotope ratio tracers," in review at the Journal of Advances in Modeling Earth Systems. A file corresponding to each of the tagged and isotopic variables used in this manuscript is included. Files are at 0.9° latitude x 1.25° longitude, and are in NetCDF format. Data from two simulations are included: 1) a simulation where the atmospheric model was "nudged" to ERA5 wind and surface pressure fields, by adding an additional tendency (see section 3.1 of associated manuscript), and 2) a simulation where the atmospheric state was allowed to freely evolve, using only boundary conditions imposed at the surface and top of atmosphere. Specific information about each of the variables provided is located in the "usage notes" section below. Associated article abstract: The hydrologic cycle couples the Earth's energy and carbon budgets through evaporation, moisture transport, and precipitation. Despite a wealth of observations and models, fundamental limitations remain in our capacity to deduce even the most basic properties of the hydrological cycle, including the spatial pattern of the residence time (RT) of water inMore>>
  4. The discovery of oxide electronics is of increasing importance today as one of the most promising new technologies and manufacturing processes for a variety of electronic and optoelectronic applications such as next-generation displays, batteries, solar cells, and photodetectors. The high potential use seen in oxide electronics is due primarily to their high carrier mobilities and their ability to be fabricated at low temperatures. However, since the majority of oxide semiconductors are n-type oxides, current applications are limited to unipolar devices, eventually developing oxide-based bipolar devices such as p-n diodes and complementary metal-oxide semiconductors. We have contributed to wide range of oxide semiconductors and their electronics and optoelectronic device applications. Particularly, we have demonstrated n-type oxide-based thin film transistors (TFT), integrating In2O3-based n-type oxide semiconductors from binary cation materials to ternary cation species including InZnO, InGaZnO (IGZO), and InAlZnO. We have suggested channel/metallization contact strategies to achieve stable TFT performance, identified vacancy-based native defect doping mechanisms, suggested interfacial buffer layers to promote charge injection capability, and established the role of third cation species on the carrier generation and carrier transport. More recently, we have reported facile manufacturing of p-type SnOx through reactive magnetron sputtering from a Sn metal target. The fabricatedmore »p-SnOx was found to be devoid of metallic phase of Sn from x-ray photoelectron spectroscopy and demonstrated stable performance in a fully oxide based p-n heterojunction together with n-InGaZnO. The oxide-based p-n junctions exhibited a high rectification ratio greater than 103 at ±3 V, a low saturation current of ~2x10-10, and a small turn-on voltage of -0.5 V. With all the previous achievements and investigations about p-type oxide semiconductors, challenges remain for implementing p-type oxide realization. For the implementation of oxide-based p-n heterojunctions, the performance needs to be further enhanced. The current on/off ration may be limited, in our device structure, due to either high reverse saturation current (or current density) or non-ideal performance. In this study, two rational strategies are suggested to introduce an “intrinsic” layer, which is expected to reduce the reverse saturation current between p-SnOx and n-IGZO and hence increase the on/off ratio. The carrier density of n-IGZO is engineered in-situ during the sputtering process, by which compositionally homogeneous IGZO with significantly reduced carrier density is formed at the interface. Then, higher carrier density IGZO is formed continuously on the lower carrier density IGZO during the sputtering process without any exposure of the sample to the air. Alternatively, heterogeneous oxides of MgO and SiO2 are integrated into between p-SnOx and n-IGZO, by which the defects on the surface can be passivated. The interfacial properties are thoroughly investigated using transmission electron microscopy and atomic force microscopy. The I-V characteristics are compared between the set of devices integrated with two types of “intrinsic” layers. The current research results are expected to contribute to the development of p-type oxides and their industrial application manufacturing process that meets current processing requirements, such as mass production in p-type oxide semiconductors.« less
  5. Abstract. Polycyclic aromatic hydrocarbons (PAHs) are carcinogenic air pollutants. The dispersion of PAHs in the atmosphere is influenced by gas–particle partitioning and chemical loss. These processes are closely interlinked and may occur at vastly differing timescales, which complicates their mathematical description in chemical transport models. Here, we use a kinetic model that explicitly resolves mass transport and chemical reactions in the gas and particle phases to describe and explore the dynamic and non-equilibrium interplay of gas–particle partitioning and chemical losses of PAHs on soot particles. We define the equilibration timescale τeq of gas–particle partitioning as the e-folding time for relaxation of the system to the partitioning equilibrium. We find this metric to span from seconds to hours depending on temperature, particle surface area, and the type of PAH. The equilibration time can be approximated using a time-independent equation, τeq≈1kdes+kads, which depends on the desorption rate coefficient kdes and adsorption rate coefficient kads, both of which can be calculated from experimentally accessible parameters. The model reveals two regimes in which different physical processes control the equilibration timescale: a desorption-controlled and an adsorption-controlled regime. In a case study with the PAH pyrene, we illustrate how chemical loss can perturb the equilibrium particulatemore »fraction at typical atmospheric concentrations of O3 and OH. For the surface reaction with O3, the perturbation is significant and increases with the gas-phase concentration of O3. Conversely, perturbations are smaller for reaction with the OH radical, which reacts with pyrene on both the surface of particles and in the gas phase. Global and regional chemical transport models typically approximate gas–particle partitioning with instantaneous-equilibration approaches. We highlight scenarios in which these approximations deviate from the explicitly coupled treatment of gas–particle partitioning and chemistry presented in this study. We find that the discrepancy between solutions depends on the operator-splitting time step and the choice of time step can help to minimize the discrepancy. The findings and techniques presented in this work not only are relevant for PAHs but can also be applied to other semi-volatile substances that undergo chemical reactions and mass transport between the gas and particle phase.« less