Remoulded shear strength at plastic and semi-solid states

The undrained shear strength of remoulded soils is of great concern in geotechnical engineering applications. This study aims to develop a reliable approach for determining the undrained shear strength of remoulded fine-grained soils, through the use of index test results, at both the plastic and semi-solid states of consistency. Experimental investigation and subsequent analysis involving a number of fine-grained soils of widely varying plasticity and geological origin have led to a two-parameter linear model of the relationship between logarithm of remoulded undrained shear strength and liquidity index. The numerical values of the parameters are found to be dependent to a lesser extent on the soil group and to a greater extent on the soil state. Based on the values of regression coefficient, ranking index and ranking distance, it seems that the relationship represents the experimental results well. It may be pointed out that the possibility of such a relationship in the semi-solid state of a soil has not been explored in the past. It is also shown that the shear strength at the plastic limit is about 32-34 times that at the liquid limit.

Discussion: Remoulded shear strength at plastic and semi-solid states

The undrained shear strength of remoulded soils is of great concern in geotechnical engineering applications. This study aims to develop a reliable approach for determining the undrained shear strength of remoulded fine-grained soils, through the use of index test results, at both the plastic and semi-solid states of consistency. Experimental investigation and subsequent analysis involving a number of fine-grained soils of widely varying plasticity and geological origin have led to a two-parameter linear model of the relationship between logarithm of remoulded undrained shear strength and liquidity index. The numerical values of the parameters are found to be dependent to a lesser extent on the soil group and to a greater extent on the soil state. Based on the values of regression coefficient, ranking index and ranking distance, it seems that the relationship represents the experimental results well. It may be pointed out that the possibility of such a relationship in the semi-solid state of a soil has not been explored in the past. It is also shown that the shear strength at the plastic limit is about 32–34 times that at the liquid limit.

An enthalpy-based model of dendritic growth in a convecting binary alloy melt

Purpose-In the present work, a numerical method, based on the well established enthalpy technique, is developed to simulate the growth of binary alloy equiaxed dendrites in presence of melt convection. The paper aims to discuss these issues. Design/methodology/approach-The principle of volume-averaging is used to formulate the governing equations (mass, momentum, energy and species conservation) which are solved using a coupled explicit-implicit method. The velocity and pressure fields are obtained using a fully implicit finite volume approach whereas the energy and species conservation equations are solved explicitly to obtain the enthalpy and solute concentration fields. As a model problem, simulation of the growth of a single crystal in a two-dimensional cavity filled with an undercooled melt is performed. Findings-Comparison of the simulation results with available solutions obtained using level set method and the phase field method shows good agreement. The effects of melt flow on dendrite growth rate and solute distribution along the solid-liquid interface are studied. A faster growth rate of the upstream dendrite arm in case of binary alloys is observed, which can be attributed to the enhanced heat transfer due to convection as well as lower solute pile-up at the solid-liquid interface. Subsequently, the influence of thermal and solutal Peclet number and undercooling on the dendrite tip velocity is investigated. Originality/value-As the present enthalpy based microscopic solidification model with melt convection is based on a framework similar to popularly used enthalpy models at the macroscopic scale, it lays the foundation to develop effective multiscale solidification.

Influence of Spatial Variation of Hydraulic Conductivity of Municipal Solid Waste on Performance of Bioreactor Landfill

The current study analyzes the leachate distribution in the Orchard Hills Landfill, Davis Junction, Illinois, using a two-phase flow model to assess the influence of variability in hydraulic conductivity on the effectiveness of the existing leachate recirculation system and its operations through reliability analysis. Numerical modeling, using finite-difference code, is performed with due consideration to the spatial variation of hydraulic conductivity of the municipal solid waste (MSW). The inhomogeneous and anisotropic waste condition is assumed because it is a more realistic representation of the MSW. For the reliability analysis, the landfill is divided into 10 MSW layers with different mean values of vertical and horizontal hydraulic conductivities (decreasing from top to bottom), and the parametric study is performed by taking the coefficients of variation (COVs) as 50, 100, 150, and 200%. Monte Carlo simulations are performed to obtain statistical information (mean and COV) of output parameters of the (1) wetted area of the MSW, (2) maximum induced pore pressure, and (3) leachate outflow. The results of the reliability analysis are used to determine the influence of hydraulic conductivity on the effectiveness of the leachate recirculation and are discussed in the light of a deterministic approach. The study is useful in understanding the efficiency of the leachate recirculation system. (C) 2013 American Society of Civil Engineers.

Supersolid in a one-dimensional model of hard-core bosons

We study a system of hard-core boson on a one-dimensional lattice with frustrated next-nearest-neighbor hopping and nearest-neighbor interaction. At half filling, for equal magnitude of nearest- and next-nearest-neighbor hopping, the ground state of this system exhibits a first-order phase transition from a bond-ordered solid to a charge-density-wave solid as a function of the nearest- neighbor interaction. Moving away from half filling we investigate the system at incommensurate densities, where we find a supersolid phase which has concurrent off-diagonal long-range order and density-wave order which is unusual in a system of hard-core bosons in one dimension. Using the finite-size density-matrix renormalization group method, we obtain the complete phase diagram for this model.

Transition from Cassie to impaled state during drop impact on groove-textured solid surfaces

Liquid drops impacted on textured surfaces undergo a transition from the Cassie state characterized by the presence of air pockets inside the roughness valleys below the drop to an impaled state with at least one of the roughness valleys filled with drop liquid. This occurs when the drop impact velocity exceeds a particular value referred to as the critical impact velocity. The present study investigates such a transition process during water drop impact on surfaces textured with unidirectional parallel grooves referred to as groove-textured surfaces. The process of liquid impalement into a groove in the vicinity of drop impact through de-pinning of the three-phase contact line (TPCL) beneath the drop as well as the critical impact velocity were identified experimentally from high speed video recordings of water drop impact on six different groove-textured surfaces made from intrinsically hydrophilic (stainless steel) as well as intrinsically hydrophobic (PDMS and rough aluminum) materials. The surface energy of various 2-D configurations of liquid-vapor interface beneath the drop near the drop impact point was theoretically investigated to identify the locally stable configurations and establish a pathway for the liquid impalement process. A force balance analysis performed on the liquid-vapor interface configuration just prior to TPCL de-pinning provided an expression for the critical drop impact velocity, U-o,U-cr, beyond which the drop state transitions from the Cassie to an impaled state. The theoretical model predicts that Uo, cr increases with the increase in pillar side angle, a, and intrinsic hydrophobicity whereas it decreases with the increase in groove top width, w, of the groove-textured surface. The quantitative predictions of the theoretical model were found to show good agreement with the experimental measurements of U-o,U-cr plotted against the surface texture geometry factor in our model, {tan(alpha/2)/w}(0.5).

Study of ground state phases for spin-1/2 Falicov-Kimball model on a triangular lattice

The spin dependent Falicov-Kimball model (FKM) is studied on a triangular lattice using numerical diagonalization technique and Monte-Carlo simulation algorithm. Magnetic properties have been explored for different values of parameters: on-site Coulomb correlation U, exchange interaction J and filling of electrons. We have found that the ground state configurations exhibit long range Neel order, ferromagnetism or a mixture of both as J is varied. The magnetic moments of itinerant (d) and localized U) electrons are also studied. For the one-fourth filling case we found no magnetic moment from d- and f-electrons for U less than a critical value. `.2014 Elsevier Ltd. All rights reserved.

Model for triplet state engineering in organic light emitting diodes

Engineering the position of the lowest triplet state (T-1) relative to the first excited singlet state (S-1) is of great importance in improving the efficiencies of organic light emitting diodes and organic photovoltaic cells. We have carried out model exact calculations of substituted polyene chains to understand the factors that affect the energy gap between S-1 and T-1. The factors studied are backbone dimerisation, different donor-acceptor substitutions, and twisted geometry. The largest system studied is an 18 carbon polyene which spans a Hilbert space of about 991 x 10(6). We show that for reverse intersystem crossing process, the best system involves substituting all carbon sites on one half of the polyene with donors and the other half with acceptors. (C) 2014 AIP Publishing LLC.

Exact Solution of the PPP Model for Correlated Electronic States of Tetracene and Substituted Tetracene

Tetracene is an important conjugated molecule for device applications. We have used the diagrammatic valence bond method to obtain the desired states, in a Hilbert space of about 450 million singlets and 902 million triplets. We have also studied the donor/acceptor (D/A)-substituted tetracenes with D and A groups placed symmetrically about the long axis of the molecule. In these cases, by exploiting a new symmetry, which is a combination of C-2 symmetry and electron-hole symmetry, we are able to obtain their low-lying states. In the case of substituted tetracene, we find that optically allowed one-photon excitation gaps reduce with increasing D/A strength, while the lowest singlet triplet gap is only wealdy affected. In all the systems we have studied, the excited singlet state, S-i, is at more than twice the energy of the lowest triplet state and the second triplet is very close to the S-1 state. Thus, donor-acceptor-substituted tetracene could be a good candidate in photovoltaic device application as it satisfies energy criteria for singlet fission. We have also obtained the model exact second harmonic generation (SHG) coefficients using the correction vector method, and we find that the SHG responses increase with the increase in D/A strength.

An extended finite element model coupled with level set method for stable pitting corrosion

Mass balance between metal and electrolytic solution, separated by a moving interface, in stable pit growth results in a set of governing equations which are solved for concentration field and interface position (pit boundary evolution), which requires only three inputs, namely the solid metal concentration, saturation concentration of the dissolved metal ions and diffusion coefficient. A combined eXtended Finite Element Model (XFEM) and level set method is developed in this paper. The extended finite element model handles the jump discontinuity in the metal concentrations at the interface, by using discontinuous-derivative enrichment formulation for concentration discontinuity at the interface. This eliminates the requirement of using front conforming mesh and re-meshing after each time step as in conventional finite element method. A numerical technique known as level set method tracks the position of the moving interface and updates it over time. Numerical analysis for pitting corrosion of stainless steel 304 is presented. The above proposed method is validated by comparing the numerical results with experimental results, exact solutions and some other approximate solutions.

The generalized Onsager model for a binary gas mixture

The Onsager model for the secondary flow field in a high-speed rotating cylinder is extended to incorporate the difference in mass of the two species in a binary gas mixture. The base flow is an isothermal solid-body rotation in which there is a balance between the radial pressure gradient and the centrifugal force density for each species. Explicit expressions for the radial variation of the pressure, mass/mole fractions, and from these the radial variation of the viscosity, thermal conductivity and diffusion coefficient, are derived, and these are used in the computation of the secondary flow. For the secondary flow, the mass, momentum and energy equations in axisymmetric coordinates are expanded in an asymptotic series in a parameter epsilon = (Delta m/m(av)), where Delta m is the difference in the molecular masses of the two species, and the average molecular mass m(av) is defined as m(av) = (rho(w1)m(1) + rho(w2)m(2))/rho(w), where rho(w1) and rho(w2) are the mass densities of the two species at the wall, and rho(w) = rho(w1) + rho(w2). The equation for the master potential and the boundary conditions are derived correct to O(epsilon(2)). The leading-order equation for the master potential contains a self-adjoint sixth-order operator in the radial direction, which is different from the generalized Onsager model (Pradhan & Kumaran, J. Fluid Mech., vol. 686, 2011, pp. 109-159), since the species mass difference is included in the computation of the density, viscosity and thermal conductivity in the base state. This is solved, subject to boundary conditions, to obtain the leading approximation for the secondary flow, followed by a solution of the diffusion equation for the leading correction to the species mole fractions. The O(epsilon) and O(epsilon(2)) equations contain inhomogeneous terms that depend on the lower-order solutions, and these are solved in a hierarchical manner to obtain the O(epsilon) and O(epsilon(2)) corrections to the master potential. A similar hierarchical procedure is used for the Carrier-Maslen model for the end-cap secondary flow. The results of the Onsager hierarchy, up to O(epsilon(2)), are compared with the results of direct simulation Monte Carlo simulations for a binary hard-sphere gas mixture for secondary flow due to a wall temperature gradient, inflow/outflow of gas along the axis, as well as mass and momentum sources in the flow. There is excellent agreement between the solutions for the secondary flow correct to O(epsilon(2)) and the simulations, to within 15 %, even at a Reynolds number as low as 100, and length/diameter ratio as low as 2, for a low stratification parameter A of 0.707, and when the secondary flow velocity is as high as 0.2 times the maximum base flow velocity, and the ratio 2 Delta m/(m(1) + m(2)) is as high as 0.5. Here, the Reynolds number Re = rho(w)Omega R-2/mu, the stratification parameter A = root m Omega R-2(2)/(2k(B)T), R and Omega are the cylinder radius and angular velocity, m is the molecular mass, rho(w) is the wall density, mu is the viscosity and T is the temperature. The leading-order solutions do capture the qualitative trends, but are not in quantitative agreement.

Estimation of the components of municipal solid waste settlement

Estimation of the municipal solid waste settlements and the contribution of each of the components are essential in the estimation of the volume of the waste that can be accommodated in a landfill and increase the post-usage of the landfill. This article describes an experimental methodology for estimating and separating primary settlement, settlement owing to creep and biodegradation-induced settlement. The primary settlement and secondary settlement have been estimated and separated based on 100% pore pressure dissipation time and the coefficient of consolidation. Mechanical creep and biodegradation settlements were estimated and separated based on the observed time required for landfill gas production. The results of a series of laboratory triaxial tests, creep tests and anaerobic reactor cell setups were conducted to describe the components of settlement. All the tests were conducted on municipal solid waste (compost reject) samples. It was observed that biodegradation accounted to more than 40% of the total settlement, whereas mechanical creep contributed more than 20% towards the total settlement. The essential model parameters, such as the compression ratio (C-c'), rate of mechanical creep (c), coefficient of mechanical creep (b), rate of biodegradation (d) and the total strain owing to biodegradation (E-DG), are useful parameters in the estimation of total settlements as well as components of settlement in landfill.

Optimal Linear Glauber Model

Contrary to the actual nonlinear Glauber model, the linear Glauber model (LGM) is exactly solvable, although the detailed balance condition is not generally satisfied. This motivates us to address the issue of writing the transition rate () in a best possible linear form such that the mean squared error in satisfying the detailed balance condition is least. The advantage of this work is that, by studying the LGM analytically, we will be able to anticipate how the kinetic properties of an arbitrary Ising system depend on the temperature and the coupling constants. The analytical expressions for the optimal values of the parameters involved in the linear are obtained using a simple Moore-Penrose pseudoinverse matrix. This approach is quite general, in principle applicable to any system and can reproduce the exact results for one dimensional Ising system. In the continuum limit, we get a linear time-dependent Ginzburg-Landau equation from the Glauber's microscopic model of non-conservative dynamics. We analyze the critical and dynamic properties of the model, and show that most of the important results obtained in different studies can be reproduced by our new mathematical approach. We will also show in this paper that the effect of magnetic field can easily be studied within our approach; in particular, we show that the inverse of relaxation time changes quadratically with (weak) magnetic field and that the fluctuation-dissipation theorem is valid for our model.

Deformation field evolution in indentation of a porous brittle solid

An experimental study of plane strain wedge indentation of a model porous brittle solid has been made to understand the effect of indentation parameters on the evolution of the deformation field and the accompanying volume change. A series of high-speed, high-resolution images of the indentation region and simultaneous measurements of load response were captured through the progression of the indentation process. Particle image velocimetry analysis of the images facilitated in situ measurement of the evolution of the resulting plastic zone in terms of incremental material displacement (velocity), strain rate, strain and volume change (e.g., local pore collapse). These measurements revealed initiation and propagation of flow localizations and fractures, as well as enabled estimate of volume changes occurring in the deformation zone. The results were directly compared with theoretical estimates of indentation pressure and deformation zone geometry and were used to validate a modified cavity expansion solution that incorporates effects of volume changes in the plastic zone. (C) 2015 Elsevier Ltd. All rights reserved.

Effects of heat exchanger design on the performance of a solid state hydrogen storage device

Heat exchanger design plays a significant role in the performance of solid state hydrogen storage device. In the present study, a cylindrical hydrogen storage device with an embedded annular heat exchanger tube with radial circular copper fins, is considered. A 3-D mathematical model of the storage device is developed to investigate the sorption performance of metal hydride (MH). A prototype of the device is fabricated for 1 kg of MH alloy, LaNi5, and tested at constant supply pressure of hydrogen, validating the simulation results. Absorption characteristics of storage device have been examined by varying different operating parameters such as hydrogen supply pressure and cooling fluid temperature and velocity. Absorption process is completed in 18 min when these parameters are 15 bar, 298 K and 1 m/s respectively. A study of geometric parameters of copper fins (such as perforation, number and thickness of fin) has been carried out to investigate their effects on absorption process. Copyright (C) 2015, Hydrogen Energy Publications, LLC. Published by Elsevier Ltd. All rights reserved.