Mesoscale description of an asymmetric lamellar phase

The relationship between the parameters in a description based on a mesoscale free energy functional for the concentration field and the macroscopic properties, such as the bending and compression moduli and the permeation constant, are examined for an asymmetric lamellar phase where the mass fractions of the hydrophobic and hydrophilic parts are not equal. The difference in the mass fractions is incorporated using a cubic term in the free energy functional, in addition to the usual quadratic and quartic terms in the Landau–Ginsburg formulation. The relationship between the coefficient of the cubic term and the difference in the mass fractions of the hydrophilic and hydrophobic parts is obtained. For a lamellar phase, it is important to ensure that the surface tension is zero due to symmetry considerations. The relationship between the parameters in the free energy functional for zero surface tension is derived. When the interface between the hydrophilic and hydrophobic parts is diffuse, it is found that the bending and compression moduli, scaled by the parameters in the free energy functional, do increase as the asymmetry in the bilayer increases. When the interface between the hydrophilic and hydrophobic parts is sharp, the scaled bending and compression moduli show no dependence on the asymmetry in the bilayer. The ratio of the permeation constant in between the water and bilayer in a molecular description and the Onsager coefficient in the mesoscale description is O(1) for both sharp and diffuse interfaces and it increases as the difference in the mass fractions is increased.

Structure of polyamidoamide dendrimers up to limiting generations: A mesoscale description

The polyamidoamide (PAMAM) class of dendrimers was one of the first dendrimers synthesized by Tomalia and co-workers at Dow. Since its discovery the PAMAMs have stimulated many discussions on the structure and dynamics of such hyperbranched polymers. Many questions remain open because the huge conformation disorder combined with very similar local symmetries have made it difficult to characterize experimentally at the atomistic level the structure and dynamics of PAMAM dendrimers. The higher generation dendrimers have also been difficult to characterize computationally because of the large size (294852 atoms for generation 11) and the huge number of conformations. To help provide a practical means of atomistic computational studies, we have developed an atomistically informed coarse-grained description for the PAMAM dendrimer. We find that a two-bead per monomer representation retains the accuracy of atomistic simulations for predicting size and conformational complexity, while reducing the degrees of freedom by tenfold. This mesoscale description has allowed us to study the structural properties of PAMAM dendrimer up to generation 11 for time scale of up to several nanoseconds. The gross properties such as the radius of gyration compare very well with those from full atomistic simulation and with available small angle x-ray experiment and small angle neutron scattering data. The radial monomer density shows very similar behavior with those obtained from the fully atomistic simulation. Our approach to deriving the coarse-grain model is general and straightforward to apply to other classes of dendrimers.

Polarimetric Signatures in the Stratiform Region of a Mesoscale Convective System

Four polarimetric measurands were collected in the stratiform region of a mesoscale convective system. The four are the reflectivity factor, the differential reflectivity, the correlation coefficient between orthogonal copolar echoes, and the differential propagation constant. Most striking is a signature of large aggregates (about 10 mm in size) seen in the differential phase through the melting layer. Another significant feature is an abrupt notch in the correlation coefficient that occurs towards the bottom of the bright band. Aircraft observations and a one-dimensional cloud model are used to explain some polarimetric measurements and to infer the presence of aggregates, graupel, and supercooled cloud water in the stratiform region. These unique observations and model data provide inferences concerning the presence of graupel and the growth of large aggregates in the melting layer.

Verifying scalings for bending rigidity of bilayer membranes using mesoscale models

The bending rigidity kappa of bilayer membranes was studied with coarse grained soft repulsive potentials using dissipative particle dynamics (DPD) simulations. Using a modified Andersen barostat to maintain the bilayers in a tensionless state, the bending rigidity was obtained from a Fourier analysis of the height fluctuations. From simulations carried out over a wide range of membrane thickness, the continuum scaling relation kappa proportional to d(2) was captured for both the L-alpha and L-beta phases. For membranes with 4 to 6 tail beads, the bending rigidity in the L-beta phase was found to be 10-15 times higher than that observed for the L-alpha phase. From the quadratic scalings obtained, a six fold increase in the area stretch modulus, k(A) was observed across the transition. The magnitude of increase in both kappa and k(A) from the L-alpha to the L-beta phase is consistent with current experimental observations in lipid bilayers and to our knowledge provides for the first time a direct evaluation of the mechanical properties in the L-beta phase.

Physically consistent simulation of mesoscale chemical kinetics: The non-negative FIS-alpha method

Biochemical pathways involving chemical kinetics in medium concentrations (i.e., at mesoscale) of the reacting molecules can be approximated as chemical Langevin equations (CLE) systems. We address the physically consistent non-negative simulation of the CLE sample paths as well as the issue of non-Lipschitz diffusion coefficients when a species approaches depletion and any stiffness due to faster reactions. The non-negative Fully Implicit Stochastic alpha (FIS alpha) method in which stopped reaction channels due to depleted reactants are deleted until a reactant concentration rises again, for non-negativity preservation and in which a positive definite Jacobian is maintained to deal with possible stiffness, is proposed and analysed. The method is illustrated with the computation of active Protein Kinase C response in the Protein Kinase C pathway. (C) 2011 Elsevier Inc. All rights reserved.

Design of a Portable Compliant Device for Estimating the Failure-Load of Mesoscale Cemented Sand Specimens

In this paper, we present the design and development of a portable, hand-operated composite compliant mechanism for estimating the failure-load of cm-sized stiff objects whose stiffness is of the order of 10 s of kN/m. The motivation for the design comes from the need to estimate the failure-load of mesoscale cemented sand specimens in situ, which is not possible with traditional devices used for large specimens or very small specimens. The composite compliant device, developed in this work, consists of two compliant mechanisms: a force-amplifying compliant mechanism (FaCM) to amplify sufficiently the force exerted by hand in order to break the specimen and a displacement-amplifying compliant mechanism (DaCM) to enable measurement of the force using a proximity sensor. The two mechanisms are designed using the selection-maps technique to amplify the force up to 100N by about a factor of 3 and measure the force with a resolution of 15 mN. The composite device, made using a FaCM, a DaCM, and a Hall effect-based proximity sensor, was tested on mesoscale cemented sand specimens that were 10mm in diameter and 20mm in length. The results are compared with those of a large commercial instrument. Through the experiments, it was observed that the failure-load of the cemented sand specimens varied from 0.95N to 24.33 N, depending on the percentage of cementation and curing period. The estimation of the failure-load using the compliant device was found to be within 1.7% of the measurements obtained using the commercial instrument and thus validating the design. The details of the design, prototyping, specimen preparation, testing, and the results comprise the paper.

Multiscale modeling of lamellar mesophases

The mesoscale simulation of a lamellar mesophase based on a free energy functional is examined with the objective of determining the relationship between the parameters in the model and molecular parameters. Attention is restricted to a symmetric lamellar phase with equal volumes of hydrophilic and hydrophobic components. Apart from the lamellar spacing, there are two parameters in the free energy functional. One of the parameters, r, determines the sharpness of the interface, and it is shown how this parameter can be obtained from the interface profile in a molecular simulation. The other parameter, A, provides an energy scale. Analytical expressions are derived to relate these parameters to r and A to the bending and compression moduli and the permeation constant in the macroscopic equation to the Onsager coefficient in the concentration diffusion equation. The linear hydrodynamic response predicted by the theory is verified by carrying out a mesoscale simulation using the lattice-Boltzmann technique and verifying that the analytical predictions are in agreement with simulation results. A macroscale model based on the layer thickness field and the layer normal field is proposed, and the relationship between the parameters in the macroscale model from the parameters in the mesoscale free energy functional is obtained.

Dynamics of edge dislocations in a sheared lamellar mesophase

The dynamics and interactions of edge dislocations in a nearly aligned sheared lamellar mesophase is analysed to provide insights into the relationship between disorder and rheology. First, the mesoscale permeation and momentum equations for the displacement field in the presence of external forces are derived from the model H equations for the concentration and momentum field. The secondary flow generated due to the mean shear around an isolated defect is calculated, and the excess viscosity due to the presence of the defect is determined from the excess energy dissipation due to the secondary flow. The excess viscosity for an isolated defect is found to increase with system size in the cross-stream direction as L-3/2 for an isolated defect, though this divergence is cut-off due to interactions in a defect suspension. As the defects are sheared past each other due to the mean flow, the Peach-Koehler force due to elastic interaction between pairs of defects is found to cause no net displacement relative to each other as they approach from large separation to the distance of closest approach. The equivalent force due to viscous interactions is found to increase the separation for defects of opposite sign, and decrease the separation for defects of same sign. During defect interactions, we find that there is no buckling instability due to dilation of layers for systems of realistic size. However, there is another mechanism, which is the velocity difference generated across a slightly deformed bilayer due to the mean shear, which could result in the creation of new defects. (C) 2013 AIP Publishing LLC.

Simulations and experiments in punching spring-steel devices with sub-millimeter features

This work demonstrates the feasibility of mesoscale (100 μm to mm) punching of multiple holes of intricate shapes in metals. Analytical modeling, finite element (FE)simulation, and experimentations are used in this work. Two dimensional FE simulations in ABAQUS were done with an assumed material modeling and plane-strain condition. A known analytical model was used and compared with the ABAQUS simulation results to understand the effects of clearance between the punch and the die. FE simulation in ABAQUS was done for different clearances and corner radii at punch, die, and holder. A set of punches and dies were used to punch out a miniature spring-steel gripper. Comparison of compliant grippers manufactured by wire-cut electro discharge machining(EDM) and punching shows that realizing sharp interior and re-entrant corners by punching is not easy to achieve. Punching of circular holes with 5 mm and 2.5 mm diameter is achieved. The possibility of realizing meso-scale parts with complicated shapes through punching is demonstrated in this work; and some strategies are suggested for improvement.

Secondary instability as a possible mechanism for clear-air turbulence: a case study

Many theories and mechanisms have been proposed to explain the phenomenon of clear-air turbulence (CAT), and some of them have been successful in predicting light, moderate and, in some cases, severe turbulence. It is only recently that skill in the forecasting of the severe form of CAT, which could lead to injuries to passengers and damage to aircraft, has improved. Recent observations and simulations suggest that some severe to extreme turbulence could be caused by horizontal vortex tubes resulting from secondary instabilities of regions of high shear in the atmosphere. We have conducted direct numerical simulations to understand the scale relationship between primary structures (larger-scale structures related to one of the causes mentioned above) and secondary structures (smaller-sized, shear structures of the size of aircraft). From shear layer simulations, we find that the ratio of sizes of primary and secondary vortices is of the right order to generate aircraft-scale vortex tubes from typical atmospheric shear layers. We have also conducted simulations with a mesoscale atmospheric model, to understand possible causes of turbulence experienced by a flight off the west coast of India. Our simulations show the occurrence of primary flow structures related to synoptic conditions around the time of the incident. The evidence presented for this mechanism also has implications for possible methods of detection and avoidance of severe CAT.

Investigations of aerosol black carbon from a semi-urban site in the Indo-Gangetic Plain region

Long-term (2009-2012) data from ground-based measurements of aerosol black carbon (BC) from a semi-urban site, Pantnagar (29.0 degrees N, 79.5 degrees E, 231 m amsl), in the Indo-Gangetic Plain (IGP) near the Himalayan foothills are analyzed to study the regional characterization. Large variations are seen in BC at both diurnal and seasonal scales, associated with the mesoscale and synoptic meteorological processes, and local/regional anthropogenic activities. BC diurnal variations show two peaks (morning and evening) arising from the combined effects of the atmospheric boundary layer (ABL) dynamics and local emissions. The diurnal amplitudes as well as the rates of diurnal evolution are the highest in winter season, followed by autumn, and the lowest in summer-monsoon. BC exhibits nearly an inverse relation with mixing layer depth in all seasons; being strongest in winter (R-2 = 0.89) and weakest (R-2 = 0.33) in monsoon (July-August). Unlike BC, co-located aerosol optical depths (AOD) and aerosol absorption are highest in spring over IGP, probably due to the presence of higher abundances of aerosols (including dust) above the ABL (in the free troposphere). AOD (500 nm) showed annual peak (>0.6) in May-June, dominated by coarse mode, while fine mode aerosols dominated in late autumn and early winter. Aerosols profiles from CALIPSO show highest values close to the surface in winter/autumn, similar to the feature seen in surface BC, whereas at altitudes > 2 km, the extinction is maximum in spring/summer. WRF-Chem model is used to simulate BC temporal variations and then compared with observed BC. The model captures most of the important features of the diurnal and seasonal variations but significantly underestimated the observed BC levels, suggesting improvements in diurnal and seasonal varying BC emissions apart from the boundary layer processes. (C) 2015 Elsevier Ltd. All rights reserved.

Structure-rheology relationship in a sheared lamellar fluid

The structure-rheology relationship in the shear alignment of a lamellar fluid is studied using a mesoscale model which provides access to the lamellar configurations and the rheology. Based on the equations and free energy functional, the complete set of dimensionless groups that characterize the system are the Reynolds number (rho gamma L-2/mu), the Schmidt number (mu/rho D), the Ericksen number (mu(gamma)/B), the interface sharpness parameter r, the ratio of the viscosities of the hydrophilic and hydrophobic parts mu(r), and the ratio of the system size and layer spacing (L/lambda). Here, rho and mu are the fluid density and average viscosity, (gamma) over dot is the applied strain rate, D is the coefficient of diffusion, B is the compression modulus, mu(r) is the maximum difference in the viscosity of the hydrophilic and hydrophobic parts divided by the average viscosity, and L is the system size in the cross-stream direction. The lattice Boltzmann method is used to solve the concentration and momentum equations for a two dimensional system of moderate size (L/lambda = 32) and for a low Reynolds number, and the other parameters are systematically varied to examine the qualitative features of the structure and viscosity evolution in different regimes. At low Schmidt numbers where mass diffusion is faster than momentum diffusion, there is fast local formation of randomly aligned domains with ``grain boundaries,'' which are rotated by the shear flow to align along the extensional axis as time increases. This configuration offers a high resistance to flow, and the layers do not align in the flow direction even after 1000 strain units, resulting in a viscosity higher than that for an aligned lamellar phase. At high Schmidt numbers where momentum diffusion is fast, the shear flow disrupts layers before they are fully formed by diffusion, and alignment takes place by the breakage and reformation of layers by shear, resulting in defects (edge dislocations) embedded in a background of nearly aligned layers. At high Ericksen number where the viscous forces are large compared to the restoring forces due to layer compression and bending, shear tends to homogenize the concentration field, and the viscosity decreases significantly. At very high Ericksen number, shear even disrupts the layering of the lamellar phase. At low Ericksen number, shear results in the formation of well aligned layers with edge dislocations. However, these edge dislocations take a long time to anneal; the relatively small misalignment due to the defects results in a large increase in viscosity due to high layer stiffness and due to shear localization, because the layers between defects get pinned and move as a plug with no shear. An increase in the viscosity contrast between the hydrophilic and hydrophobic parts does not alter the structural characteristics during alignment. However, there is a significant increase in the viscosity, due to pinning of the layers between defects, which results in a plug flow between defects and a localization of the shear to a part of the domain.