Dynamics of pulled desorption with effects of excluded-volume interaction: The p-Laplacian diffusion equation and its exact solution

We analyze the dynamics of desorption of a polymer molecule which is pulled at one of its ends with force f, trying to desorb it. We assume a monomer to desorb when the pulling force on it exceeds a critical value f(c). We formulate an equation for the average position of the n-th monomer, which takes into account excluded-volume interaction through the blob-picture of a polymer under external constraints. The approach leads to a diffusion equation with a p-Laplacian for the propagation of the stretching along the chain. This has to be solved subject to a moving boundary condition. Interestingly, within this approach, the problem can be solved exactly in the trumpet, stem-flower and stem regimes. In the trumpet regime, we get tau = tau(0)n(d)(2), where n(d) is the number of monomers that have desorbed at the time tau. tau(0) is known only numerically, but for f close to f(c), it is found to be tau(0) similar to f(c)/(f(2/3) - f(c)(2/3)) If one used simple Rouse dynamics, this result would change to tau similar to f(c)n(d)(2)/(f - f(c)). In the other regimes too, one can find exact solution, and interestingly, in all regimes tau similar to n(d)(2). Copyright (C) EPLA, 2011

Critical phenomena in polymer solutions: Scaling of the free energy

The thermodynamics of monodisperse solutions of polymers in the neighborhood of the phase separation temperature is studied by means of Wilson’s recursion relation approach, starting from an effective ϕ4 Hamiltonian derived from a continuum model of a many‐chain system in poor solvents. Details of the chain statistics are contained in the coefficients of the field variables ϕ, so that the parameter space of the Hamiltonian includes the temperature, coupling constant, molecular weight, and excluded volume interaction. The recursion relations are solved under a series of simplifying assumptions, providing the scaling forms of the relevant parameters, which are then used to determine the scaling form of the free energy. The free energy, in turn, is used to calculate the other singular thermodynamic properties of the solution. These are characteristically power laws in the reduced temperature and molecular weight, with the temperature exponents being the same as those of the 3d Ising model. The molecular weight exponents are unique to polymer solutions, and the calculated values compare well with the available experimental data.

Computer simulation of liquid crystals

In this article we review the current status in the modelling of both thermotropic and lyotropic Liquid crystal. We discuss various coarse-graining schemes as well as simulation techniques such as Monte Carlo (MC) and Molecular dynamics (MD) simulations.In the area of MC simulations we discuss in detail the algorithm for simulating hard objects such as spherocylinders of various aspect ratios where excluded volume interaction enters in the simulation through overlap test. We use this technique to study the phase diagram, of a special class of thermotropic liquid crystals namely banana liquid crystals. Next we discuss a coarse-grain model of surfactant molecules and study the self-assembly of the surfactant oligomers using MD simulations. Finally we discuss an atomistically informed coarse-grained description of the lipid molecules used to study the gel to liquid crystalline phase transition in the lipid bilayer system.

Confinement and viscoelastic effects on chain closure dynamics

Chemical reactions inside cells are typically subject to the effects both of the cell's confining surfaces and of the viscoelastic behavior of its contents. In this paper, we show how the outcome of one particular reaction of relevance to cellular biochemistry - the diffusion-limited cyclization of long chain polymers - is influenced by such confinement and crowding effects. More specifically, starting from the Rouse model of polymer dynamics, and invoking the Wilemski-Fixman approximation, we determine the scaling relationship between the mean closure time t(c) of a flexible chain (no excluded volume or hydrodynamic interactions) and the length N of its contour under the following separate conditions: (a) confinement of the chain to a sphere of radius d and (b) modulation of its dynamics by colored Gaussian noise. Among other results, we find that in case (a) when d is much smaller than the size of the chain, t(c) similar to Nd-2, and that in case (b), t(c) similar to N-2/(2 (2H)), H being a number between 1/2 and 1 that characterizes the decay of the noise correlations. H is not known a priori, but values of about 0.7 have been used in the successful characterization of protein conformational dynamics. At this value of H (selected for purposes of illustration), t(c) similar to N-3.4, the high scaling exponent reflecting the slow relaxation of the chain in a viscoelastic medium. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4729041]

Molecular dynamic simulations of elastomer structure and its influence on anisotropic stress under time-varying strain

We investigate the evolution of polymer structure and its influence on uniaxial anisotropic stress under time-varying uniaxial strain, and the role of external control variables such as temperature, strain rate, chain length, and density, using molecular dynamics simulation. At temperatures higher than glass transition, stress anisotropy in the system is reduced even though the bond stretch is greater at higher temperatures. There is a significant increase in the stress level with increasing density. At higher densities, the uncoiling of the chains is suppressed and the major contribution to the deformation is by internal deformation of the chains. At faster rates of loading stress anisotropy increases. The deformation mechanism is mostly due to bond stretch and bond bending rather than overall shape and size. Stress levels increase with longer chain length. There is a critical value of the functionality of the cross-linkers beyond which the uniaxial stress developed increases caused primarily by bond stretching due to increased constraint on the motion of the monomers. Stacking of the chains in the system also plays a dominant role in the behaviour in terms of excluded volume interactions. Low density, high temperature, low values of functionality of cross-linkers, and short chain length facilitate chain uncoiling and chain slipping in cross-linked polymers.

Dynamics of the reaction between the free end of a tethered self-avoiding polymer and a flat penetrable surface: A renormalization group study

The average time tau(r) for one end of a long, self-avoiding polymer to interact for the first time with a flat penetrable surface to which it is attached at the other end is shown here to scale essentially as the square of the chain's contour length N. This result is obtained within the framework of the Wilemski-Fixman approximation to diffusion-limited reactions, in which the reaction time is expressed as a time correlation function of a ``sink'' term. In the present work, this sink-sink correlation function is calculated using perturbation expansions in the excluded volume and the polymer-surface interactions, with renormalization group methods being used to resum the expansion into a power law form. The quadratic dependence of tau(r) on N mirrors the behavior of the average time tau(c) of a free random walk to cyclize, but contrasts with the cyclization time of a free self-avoiding walk (SAW), for which tau(r) similar to N-2.2. A simulation study by Cheng and Makarov J. Phys. Chem. B 114, 3321 (2010)] of the chain-end reaction time of an SAW on a flat impenetrable surface leads to the same N-2.2 behavior, which is surprising given the reduced conformational space a tethered polymer has to explore in order to react. (C) 2014 AIP Publishing LLC.

Finite volume evolution Galerkin method for hyperbolic conservation laws with spatially varying flux functions

We present a generalization of the finite volume evolution Galerkin scheme [M. Lukacova-Medvid'ova,J. Saibertov'a, G. Warnecke, Finite volume evolution Galerkin methods for nonlinear hyperbolic systems, J. Comp. Phys. (2002) 183 533-562; M. Luacova-Medvid'ova, K.W. Morton, G. Warnecke, Finite volume evolution Galerkin (FVEG) methods for hyperbolic problems, SIAM J. Sci. Comput. (2004) 26 1-30] for hyperbolic systems with spatially varying flux functions. Our goal is to develop a genuinely multi-dimensional numerical scheme for wave propagation problems in a heterogeneous media. We illustrate our methodology for acoustic waves in a heterogeneous medium but the results can be generalized to more complex systems. The finite volume evolution Galerkin (FVEG) method is a predictor-corrector method combining the finite volume corrector step with the evolutionary predictor step. In order to evolve fluxes along the cell interfaces we use multi-dimensional approximate evolution operator. The latter is constructed using the theory of bicharacteristics under the assumption of spatially dependent wave speeds. To approximate heterogeneous medium a staggered grid approach is used. Several numerical experiments for wave propagation with continuous as well as discontinuous wave speeds confirm the robustness and reliability of the new FVEG scheme.

Prediction of compressive strength of SCC and HPC with high volume fly ash using ANN

An artificial neural network (ANN) is presented to predict a 28-day compressive strength of a normal and high strength self compacting concrete (SCC) and high performance concrete (HPC) with high volume fly ash. The ANN is trained by the data available in literature on normal volume fly ash because data on SCC with high volume fly ash is not available in sufficient quantity. Further, while predicting the strength of HPC the same data meant for SCC has been used to train in order to economise on computational effort. The compressive strengths of SCC and HPC as well as slump flow of SCC estimated by the proposed neural network are validated by experimental results.

Role of free volume in strain softening of as-cast and annealed bulk metallic glass

Plasticity in amorphous alloys is associated with strain softening, induced by the creation of additional free volume during deformation. In this paper, the role of free volume, which was a priori in the material, on work softening was investigated. For this, an as-cast Zr-based bulk metallic glass (BMG) was systematically annealed below its glass transition temperature, so as to reduce the free volume content. The bonded-interface indentation technique is used to generate extensively deformed and well defined plastic zones. Nanoindentation was utilized to estimate the hardness of the deformed as well as undeformed regions. The results show that the structural relaxation annealing enhances the hardness and that both the subsurface shear band number density and the plastic zone size decrease with annealing time. The serrations in the nanoindentation load-displacement curves become smoother with structural relaxation. Regardless of the annealing condition, the nanohardness of the deformed regions is similar to 12-15% lower, implying that the prior free volume only changes the yield stress (or hardness) but not the relative flow stress (or the extent of strain softening). Statistical distributions of the nanohardness obtained from deformed and undeformed regions have no overlap, suggesting that shear band number density has no influence on the plastic characteristics of the deformed region.

Finite volume TVD formulation of lattice Boltzmann simulation on unstructured mesh

A numerical scheme is presented for accurate simulation of fluid flow using the lattice Boltzmann equation (LBE) on unstructured mesh. A finite volume approach is adopted to discretize the LBE on a cell-centered, arbitrary shaped, triangular tessellation. The formulation includes a formal, second order discretization using a Total Variation Diminishing (TVD) scheme for the terms representing advection of the distribution function in physical space, due to microscopic particle motion. The advantage of the LBE approach is exploited by implementing the scheme in a new computer code to run on a parallel computing system. Performance of the new formulation is systematically investigated by simulating four benchmark flows of increasing complexity, namely (1) flow in a plane channel, (2) unsteady Couette flow, (3) flow caused by a moving lid over a 2D square cavity and (4) flow over a circular cylinder. For each of these flows, the present scheme is validated with the results from Navier-Stokes computations as well as lattice Boltzmann simulations on regular mesh. It is shown that the scheme is robust and accurate for the different test problems studied.

Critical Evaluation of Determining Swelling Pressure by Swell-Load Method and Constant Volume Method

For any construction activity in expansive soils, determination of swelling pressure/heave is an essential step. Though many attempts have been made to develop laboratory procedures by using the laboratory one-dimensional oedometer to determine swelling pressure of expansive soils, they are reported to yield varying results. The main reason for these variations could be heterogeneous moisture distribution of the sample over its thickness. To overcome this variation the experimental procedure should be such that the soil gets fully saturated. Attempts were made to introduce vertical sand drains in addition to the top and bottom drains. In this study five and nine vertical sand drains were introduced to experimentally find out the variations in the swell and swelling pressure. The variations in the moisture content at middle, top, and bottom of the sample in the oedometer test are also reported. It is found that swell-load method is better as compared to zero-swell method. Further, five number of vertical sand drains are found to be sufficient to obtain uniform moisture content distribution.

Effectiveness factor for zeroth-order catalytic reactions with volume change and nonuniform catalyst activity profiles

Analytical solutions are presented for the effectiveness factor of a zeroth-order reaction with volume change and nonuniform catalyst activity profile in slab, cylinder and spherical pellets. The possibility of shape normalization is considered for a variety of activity profiles and pellet shapes. When the catalyst activity at the external surface of the pellet is non-zero, shape normalization is obtained, which makes the asymptotic behavior of the effectiveness factor identical for small and large values of Thiele modulus, however, the normalization can lead to significant errors, particularly for the case of activity profiles decreasing towards the outer surface of the catalyst.

R-parameter: A local truncation error based adaptive framework for finite volume compressible flow solvers

A residual-based strategy to estimate the local truncation error in a finite volume framework for steady compressible flows is proposed. This estimator, referred to as the -parameter, is derived from the imbalance arising from the use of an exact operator on the numerical solution for conservation laws. The behaviour of the residual estimator for linear and non-linear hyperbolic problems is systematically analysed. The relationship of the residual to the global error is also studied. The -parameter is used to derive a target length scale and consequently devise a suitable criterion for refinement/derefinement. This strategy, devoid of any user-defined parameters, is validated using two standard test cases involving smooth flows. A hybrid adaptive strategy based on both the error indicators and the -parameter, for flows involving shocks is also developed. Numerical studies on several compressible flow cases show that the adaptive algorithm performs excellently well in both two and three dimensions.

Rocket Propellant Combustion Studies in a Constant Volume Bomb

A constant volume window bomb has been used to measure the characteristic velocity (c*) of rocket propellants. Analysis of the combustion process inside the bomb including heat losses has been made. The experiments on double base and composite propellants have revealed some (i) basic heat transfer aspects inside the bomb and (ii) combustion characteristics of Ammonium Perchlorate-Polyester propellants. It has been found that combustion continues even beyond the peak pressure and temperature points. Lithium Fluoride mixed propellants do not seem to indicate significant differences in c*) though the low pressure deflagration limit is increased with percentage of Lithium Fluoride.

Vitrification, relaxation and free volume in glycerol-water binary liquid mixture:Spin probe ESR studies

Glass transition and relaxation of the glycerol-water (G-W) binary mixture system have been studied over the glycerol concentration range of 5-85 mol% by using the highly sensitive technique of electron spin resonance (ESR). For the water rich mixture the glass transition,sensed by the dissolved spin probe, arises from the vitrified mesoscopic portion of the binary system. The concentration dependence of the glass transition temperature manifests a closely related molecular level cooperativity in the system. A drastic change in the mesoscopic structure of the system at the critical concentration of 40 mol is confirmed by an estimation of the spin probe effective volume in a temperature range where the tracer reorientation is strongly coupled to the system dynamics.