Homogeneous isotropic superfluid turbulence in two dimensions: Inverse and forward cascades in the Hall-Vinen-Bekharevich-Khalatnikov model

We present the first direct-numerical-simulation study of the statistical properties of two-dimensional superfluid turbulence in the simplified, Hall-Vinen-Bekharevich-Khalatnikov two-fluid model. We show that both normalfluid and superfluid energy spectra can exhibit two power-law regimes, the first associated with an inverse cascade of energy and the second with the forward cascade of enstrophy. We quantify the mutual-friction-induced alignment of normal and superfluid velocities by obtaining probability distribution functions of the angle between them and the ratio of their moduli.

A two-dimensional model for crack closure effect in plates under bending

A two-dimensional model is proposed for taking into account the establishment of contact on the compression side of crack faces in plates under bending. An approximate but simple method is developed for evaluating reduction of stress intensity factor due to such ‘crack closure’. Analysis is first carried out permitting interference of the crack faces. Contact forces are then introduced on the crack faces and their magnitudes determined from the consideration that the interference is just eliminated. The method is based partly on finite element analysis and partly on a continuum analysis using Irwin's solution for point loads on the crack line.

Cooling of a composite slab in a two-fluid medium

A mixed boundary value problem associated with the diffusion equation that involves the physical problem of cooling of an infinite parallel-sided composite slab in a two-fluid medium, is solved completely by using the Wiener-Hopf technique. An analytical solution is derived for the temperature distribution at the quench fronts being created by two different layers of cold fluids having different cooling abilities moving on the upper surface of the slab at constant speedv. Simple expressions are derived for the values of the sputtering temperatures of the slab at the points of contact with the respective layers, assuming the front layer of the fluid to be of finite width and the back layer of infinite extent. The main problem is solved through a three-part Wiener-Hopf problem of a special type and the numerical results under certain special circumstances are obtained and presented in the form of a table.

Cooling of an infinite slab in a two-fluid medium

A mixed boundary-valued problem associated with the diffusion equation, that involves the physical problem of cooling of an infinite slab in a two-fluid medium, is solved completely by using the Wiener-Hopf technique. An analytical solution is derived for the temperature distribution at the quench fronts being created by two different layers of cold fluids having different cooling abilities moving on the upper surface of the slab at constant speed. Simple expressions are derived for the values of the sputtering temperatures of the slab at the points of contact with the respective layers, assuming one layer of the fluid to be of finite extent and the other of infinite extent. The main problem is solved through a three-part Wiener - Hopf problem of a special type, and the numerical results under certain special circumstances are obtained and presented in the form of a table.

Application of scaled particle theory to alkyl-aryl surfactants: a two zone model of micellar core

Scaled Particle Theory (SPT) has been applied to predict the total free energies of micellization of ionic as well as nonionic micellar systems containing an aryl ring. A modification of the previously developed model has been made, proposing a two-zone model of micellar core which corroborates with the structural information available for such systems. The results are in good agreement with experimental data and also confirm the dictating role of cavity forming free energies for such systems

The sputtering temperature of a cooling cylindrical rod without and with an insulated core in a two-fluid medium

Utilising Jones' method associated with the Wiener-Hopf technique, explicit solutions are obtained for the temperature distributions on the surface of a cylindrical rod without an insulated core as well as that inside a cylindrical rod with an insulated inner core when the rod, in either of the two cases, is allowed to enter, with a uniform speed, into two different layers of fluid with different cooling abilities. Simple expressions are derived for the values of the sputtering temperatures of the rod at the points of entry into the respective layers, assuming the upper layer of the fluid to be of finite depth and the lower of infinite extent. Both the problems are solved through a three-part Wiener-Hopf problem of special type and the numerical results under certain special circumstances are obtained and presented in tabular forms.

Reaction dynamics under confinement: an exact path integral treatment of a two-stage model of stochastic gene expression

Gene expression in living systems is inherently stochastic, and tends to produce varying numbers of proteins over repeated cycles of transcription and translation. In this paper, an expression is derived for the steady-state protein number distribution starting from a two-stage kinetic model of the gene expression process involving p proteins and r mRNAs. The derivation is based on an exact path integral evaluation of the joint distribution, P(p, r, t), of p and r at time t, which can be expressed in terms of the coupled Langevin equations for p and r that represent the two-stage model in continuum form. The steady-state distribution of p alone, P(p), is obtained from P(p, r, t) (a bivariate Gaussian) by integrating out the r degrees of freedom and taking the limit t -> infinity. P(p) is found to be proportional to the product of a Gaussian and a complementary error function. It provides a generally satisfactory fit to simulation data on the same two-stage process when the translational efficiency (a measure of intrinsic noise levels in the system) is relatively low; it is less successful as a model of the data when the translational efficiency (and noise levels) are high.

Solid-solid collapse transition in a two dimensional model molecular system

Solid-solid collapse transition in open framework structures is ubiquitous in nature. The real difficulty in understanding detailed microscopic aspects of such transitions in molecular systems arises from the interplay between different energy and length scales involved in molecular systems, often mediated through a solvent. In this work we employ Monte-Carlo simulation to study the collapse transition in a model molecular system interacting via both isotropic as well as anisotropic interactions having different length and energy scales. The model we use is known as Mercedes-Benz (MB), which, for a specific set of parameters, sustains two solid phases: honeycomb and oblique. In order to study the temperature induced collapse transition, we start with a metastable honeycomb solid and induce transition by increasing temperature. High density oblique solid so formed has two characteristic length scales corresponding to isotropic and anisotropic parts of interaction potential. Contrary to the common belief and classical nucleation theory, interestingly, we find linear strip-like nucleating clusters having significantly different order and average coordination number than the bulk stable phase. In the early stage of growth, the cluster grows as a linear strip, followed by branched and ring-like strips. The geometry of growing cluster is a consequence of the delicate balance between two types of interactions, which enables the dominance of stabilizing energy over destabilizing surface energy. The nucleus of stable oblique phase is wetted by intermediate order particles, which minimizes the surface free energy. In the case of pressure induced transition at low temperature the collapsed state is a disordered solid. The disordered solid phase has diverse local quasi-stable structures along with oblique-solid like domains. (C) 2013 AIP Publishing LLC.

Discrete Vortex Method-Based Model for Ground-Effect Studies

A discrete vortex method-based model has been proposed for two-dimensional/three-dimensional ground-effect prediction. The model merely requires two-dimensional sectional aerodynamics in free flight. This free-flight data can be obtained either from experiments or a high-fidelity computational fluid dynamics solver. The first step of this two-step model involves a constrained optimization procedure that modifies the vortex distribution on the camber line as obtained from a discrete vortex method to match the free-flight data from experiments/computational fluid dynamics. In the second step, the vortex distribution thus obtained is further modified to account for the presence of the ground plane within a discrete vortex method-based framework. Whereas the predictability of the lift appears as a natural extension, the drag predictability within a potential flow framework is achieved through the introduction of what are referred to as drag panels. The need for the use of the generalized Kutta-Joukowski theorem is emphasized. The extension of the model to three dimensions is by the way of using the numerical lifting-line theory that allows for wing sweep. The model is extensively validated for both two-dimensional and three-dimensional ground-effect studies. The work also demonstrates the ability of the model to predict lift and drag coefficients of a high-lift wing in ground effect to about 2 and 8% accuracy, respectively, as compared to the results obtained using a Reynolds-averaged Navier-Stokes solver involving grids with several million volumes. The model shows a lot of promise in design, particularly during the early phase.

Evolutionary games with two timescales

We consider a two timescale model of learning by economic agents wherein active or 'ontogenetic' learning by individuals takes place on a fast scale and passive or 'phylogenetic' learning by society as a whole on a slow scale, each affecting the evolution of the other. The former is modelled by the Monte Carlo dynamics of physics, while the latter is modelled by the replicator dynamics of evolutionary biology. Various qualitative aspects of the dynamics are studied in some simple cases, both analytically and numerically, and its role as a useful modelling device is emphasized.

Cell Dynamics of Model Proteins

We design rapidly folding sequences by assigning the strongest couplings to the contacts present in a target native state in a two dimensional model of heteropolymers. The pathways to folding and their dependence on the temperature are illustrated via a mapping of the dynamics into motion within the space of the maximally compact cells.

Multiphase models for simulating hot and slightly miscible DNAPL (dense non-aqueous phase fluids) in a saturated rock fracture under deformation

In the present study a two dimensional model is first developed to show the behaviour of dense non-aqueous phase liquids (DNAPL) within a rough fracture. To consider the rough fracture, the fracture is imposed with variable apertures along its plane. It is found that DNAPL follows preferential pathways. In next part of the study the above model is further extended for non-isothermal DNAPL flow and DNAPL-water interphase mass transfer phenomenon. These two models are then coupled with joint deformation due to normal stresses. The primary focus of these models is specifically to elucidate the influence of joint alteration due to external stress and fluid pressures on flow driven energy transport and interphase mass transfer. For this, it is assumed that the critical value for joint alteration is associated with external stress and average of water and DNAPL pressures in multiphase system and the temporal and spatial evolution of joint alteration are determined for its further influence on energy transport and miscible phase transfer. The developed model has been studied to show the influence of deformation on DNAPL flow. Further this preliminary study demonstrates the influence of joint deformation on heat transport and phase miscibility via multiphase flow velocities. It is seen that the temperature profile changes and shows higher diffusivity due to deformation and although the interphase miscibility value decreases but the lateral dispersion increases to a considerably higher extent.

Long-range Coulomb interactions and nanoscale electronic inhomogeneities in correlated oxides

Electronic, magnetic, or structural inhomogeneities ranging in size from nanoscopic to mesoscopic scales seem endemic and are possibly generic to colossal magnetoresistance manganites and other transition metal oxides. They are hence of great current interest and understanding them is of fundamental importance. We show here that an extension, to include long-range Coulomb interactions, of a quantum two-fluid l-b model proposed recently for manganites [Phys. Rev. Lett. 92, 157203 (2004)] leads to an excellent description of such inhomogeneities. In the l-b model two very different kinds of electronic states, one localized and polaronic (l) and the other extended or broad band (b) coexist. For model parameters appropriate to manganites and even within a simple dynamical mean-field theory (DMFT) framework, it describes many of the unusual phenomena seen in manganites, including colossal magnetoresistance (CMR), qualitatively and quantitatively. However, in the absence of long-ranged Coulomb interaction, a system described by such a model would actually phase separate, into macroscopic regions of l and b electrons, respectively. As we show in this paper, in the presence of Coulomb interactions, the macroscopic phase separation gets suppressed and instead nanometer scale regions of polarons interspersed with band electron puddles appear, constituting a kind of quantum Coulomb glass. We characterize the size scales and distribution of the inhomogeneity using computer simulations. For realistic values of the long-range Coulomb interaction parameter V-0, our results for the thresholds for occupancy of the b states are in agreement with, and hence support, the earlier approach mentioned above based on a configuration averaged DMFT treatment which neglects V-0; but the present work has features that cannot be addressed in the DMFT framework. Our work points to an interplay of strong correlations, long-range Coulomb interaction, and dopant ion disorder, all inevitably present in transition metal oxides as the origin of nanoscale inhomogeneities rather than disorder frustrated phase competition as is generally believed. As regards manganites, it argues against explanations for CMR based on disorder frustrated phase separation and for an intrinsic origin of CMR. Based on this, we argue that the observed micrometer (meso) scale inhomogeneities owe their existence to extrinsic causes, e.g., strain due to cracks and defects. We suggest possible experiments to validate our speculation.

Theory of the unusual doping and temperature dependence of photoemission spectra in manganites

A recent, major, puzzle in the core-level photoemission spectra of doped manganites is the observation of a 1–2 eV wide shoulder with intensity varying with temperature T as the square of the magnetization over a T scale of order 200 K, an order of magnitude less than electronic energies. This is addressed and resolved here, by extending a recently proposed two-fluid polaron–mobile electron model for these systems to include core-hole effects. The position of the shoulder is found to be determined by Coulomb and Jahn-Teller energies, while its spectral weight is determined by the mobile electron energetics which is strongly T and doping dependent, due to annealed disorder scattering from the polarons and the t2g core spins. Our theory accounts quantitatively for the observed T dependence of the difference spectra, and furthermore, explains the observed correspondence between spectral changes due to increasing doping and decreasing T.