Nonlinear preconditioning for efficient and accurate interface capturing in simulation of multicomponent compressible flows

Single fluid schemes that rely on an interface function for phase identification in multicomponent compressible flows are widely used to study hydrodynamic flow phenomena in several diverse applications. Simulations based on standard numerical implementation of these schemes suffer from an artificial increase in the width of the interface function owing to the numerical dissipation introduced by an upwind discretization of the governing equations. In addition, monotonicity requirements which ensure that the sharp interface function remains bounded at all times necessitate use of low-order accurate discretization strategies. This results in a significant reduction in accuracy along with a loss of intricate flow features. In this paper we develop a nonlinear transformation based interface capturing method which achieves superior accuracy without compromising the simplicity, computational efficiency and robustness of the original flow solver. A nonlinear map from the signed distance function to the sigmoid type interface function is used to effectively couple a standard single fluid shock and interface capturing scheme with a high-order accurate constrained level set reinitialization method in a way that allows for oscillation-free transport of the sharp material interface. Imposition of a maximum principle, which ensures that the multidimensional preconditioned interface capturing method does not produce new maxima or minima even in the extreme events of interface merger or breakup, allows for an explicit determination of the interface thickness in terms of the grid spacing. A narrow band method is formulated in order to localize computations pertinent to the preconditioned interface capturing method. Numerical tests in one dimension reveal a significant improvement in accuracy and convergence; in stark contrast to the conventional scheme, the proposed method retains its accuracy and convergence characteristics in a shifted reference frame. Results from the test cases in two dimensions show that the nonlinear transformation based interface capturing method outperforms both the conventional method and an interface capturing method without nonlinear transformation in resolving intricate flow features such as sheet jetting in the shock-induced cavity collapse. The ability of the proposed method in accounting for the gravitational and surface tension forces besides compressibility is demonstrated through a model fully three-dimensional problem concerning droplet splash and formation of a crownlike feature. (C) 2014 Elsevier Inc. All rights reserved.

A meshless adaptive multiscale method for fracture

The paper presents a multiscale method for crack propagation. The coarse region is modelled by the differential reproducing kernel particle method. Fracture in the coarse scale region is modelled with the Phantom node method. A molecular statics approach is employed in the fine scale where crack propagation is modelled naturally by breaking of bonds. The triangular lattice corresponds to the lattice structure of the (111) plane of an FCC crystal in the fine scale region. The Lennard-Jones potential is used to model the atom-atom interactions. The coupling between the coarse scale and fine scale is realized through ghost atoms. The ghost atom positions are interpolated from the coarse scale solution and enforced as boundary conditions on the fine scale. The fine scale region is adaptively refined and coarsened as the crack propagates. The centro symmetry parameter is used to detect the crack tip location. The method is implemented in two dimensions. The results are compared to pure atomistic simulations and show excellent agreement. (C) 2014 Elsevier B. V. All rights reserved.

Dielectric elastomers: Asymptotically-correct three-dimensional displacement field

Asymptotically-accurate dimensional reduction from three to two dimensions and recovery of 3-D displacement field of non-prestretched dielectric hyperelastic membranes are carried out using the Variational Asymptotic Method (VAM) with moderate strains and very small ratio of the membrane thickness to its shortest wavelength of the deformation along the plate reference surface chosen as the small parameters for asymptotic expansion. Present work incorporates large deformations (displacements and rotations), material nonlinearity (hyperelasticity), and electrical effects. It begins with 3-D nonlinear electroelastic energy and mathematically splits the analysis into a one-dimensional (1-D) through-the-thickness analysis and a 2-D nonlinear plate analysis. Major contribution of this paper is a comprehensive nonlinear through-the-thickness analysis which provides a 2-D energy asymptotically equivalent of the 3-D energy, a 2-D constitutive relation between the 2-D generalized strain and stress tensors for the plate analysis and a set of recovery relations to express the 3-D displacement field. Analytical expressions are derived for warping functions and stiffness coefficients. This is the first attempt to integrate an analytical work on asymptotically-accurate nonlinear electro-elastic constitutive relation for compressible dielectric hyperelastic model with a generalized finite element analysis of plates to provide 3-D displacement fields using VAM. A unified software package `VAMNLM' (Variational Asymptotic Method applied to Non-Linear Material models) was developed to carry out 1-D non-linear analysis (analytical), 2-D non-linear finite element analysis and 3-D recovery analysis. The applicability of the current theory is demonstrated through an actuation test case, for which distribution of 3-D displacements are provided. (C) 2014 Elsevier Ltd. All rights reserved.

Interphase anisotropy effects on lamellar eutectics: A numerical study

In directional solidification of binary eutectics, it is often observed that two-phase lamellar growth patterns grow tilted with respect to the direction z of the imposed temperature gradient. This crystallographic effect depends on the orientation of the two crystal phases alpha and beta with respect to z. Recently, an approximate theory was formulated that predicts the lamellar tilt angle as a function of the anisotropy of the free energy of the solid(alpha)-solid(beta) interphase boundary. We use two different numerical methods-phase field (PF) and dynamic boundary integral (BI)-to simulate the growth of steady periodic patterns in two dimensions as a function of the angle theta(R) between z and a reference crystallographic axis for a fixed relative orientation of alpha and beta crystals, that is, for a given anisotropy function (Wulff plot) of the interphase boundary. For Wulff plots without unstable interphase-boundary orientations, the two simulation methods are in excellent agreement with each other and confirm the general validity of the previously proposed theory. In addition, a crystallographic ``locking'' of the lamellae onto a facet plane is well reproduced in the simulations. When unstable orientations are present in the Wulff plot, it is expected that two distinct values of the tilt angle can appear for the same crystal orientation over a finite theta(R) range. This bistable behavior, which has been observed experimentally, is well reproduced by BI simulations but not by the PF model. Possible reasons for this discrepancy are discussed.

Effect of Surface Energy Anisotropy on the Stability of Growth Fronts in Multiphase Alloys

Eutectic growth offers a variety of examples for pattern formation which are interesting both for theoreticians as well as experimentalists. One such example of patterns is ternary eutectic colonies which arise as a result of instabilities during growth of two solid phases. Here, in addition to the two major components being exchanged between the solid phases during eutectic growth, there is an impurity component which is rejected by both solid phases. During progress of solidification, there develops a boundary layer of the third impurity component ahead of the solidification front of the two solid phases. Similar to Mullins-Sekerka type instabilities, such a boundary layer tends to make the global solidification envelope unstable to morphological perturbations giving rise to two-phase cells. This phenomenon has been studied numerically in two dimensions for the conditions of directional solidification, by Plapp and Karma (Phys Rev E 66:061608, 2002) using phase-field simulations. While, in the work by Plapp and Karma (Phys Rev E 66:061608, 2002) all interfaces are isotropic, in our presentation, we extend the phase-field model by considering interfacial anisotropy in the solid-solid and solid-liquid interfaces and characterize the role of interfacial anisotropy on the stability of the growth front through phase-field simulations in two dimensions.

Shape-Dependent Confinement in Ultrasmall Zero-, One-, and Two-Dimensional PbS Nanostructures

Spatial dimensionality affects the degree of confinement when an electron-hole pair is squeezed from one or more dimensions approaching the bulk exciton Bohr radius (alpha(B)) limit. The etectron-hole interaction in zero-dimensional (0D) dots, one-dimensional (1D) rods/wires, and two-dimensional (2D) wells/sheets should be enhanced by the increase in confinement dimensions in the order 0D > 1D > 2D. We report the controlled synthesis of PbS nanomateriats with 0D, 1D, and 2D forms retaining at least one dimension in the strongly confined regime far below alpha(B) (similar to 10 nm for PbS) and provide evidence through varying the exciton-phonon coupling strength that the degree of confinement is systematically weakened by the loss of confinement dimension. Geometry variations show distinguishable far-field optical polarizations, which could find useful applications in polarization-sensitive devices.

A solution framework for scheduling a BPM with non-identical job dimensions

A Batch Processing Machine (BPM) is one which processes a number of jobs simultaneously as a batch with common beginning and ending times. Also, a BPM, once started cannot be interrupted in between (Pre-emption not allowed). This research is motivated by a BPM in steel casting industry. There are three main stages in any steel casting industry viz., pre-casting stage, casting stage and post-casting stage. A quick overview of the entire process, is shown in Figure 1. There are two BPMs : (1) Melting furnace in the pre-casting stage and (2) Heat Treatment Furnace (HTF) in the post casting stage of steel casting manufacturing process. This study focuses on scheduling the latter, namely HTF. Heat-treatment operation is one of the most important stages of steel casting industries. It determines the final properties that enable components to perform under demanding service conditions such as large mechanical load, high temperature and anti-corrosive processing. In general, different types of castings have to undergo more than one type of heat-treatment operations, where the total heat-treatment processing times change. To have a better control, castings are primarily classified into a number of job-families based on the alloy type such as low-alloy castings and high alloy castings. For technical reasons such as type of alloy, temperature level and the expected combination of heat-treatment operations, the castings from different families can not be processed together in the same batch.

Experimental studies on a two stage pulse tube cryocooler reaching 2.5K

A two stage Pulse Tube Cryocooler (PTC) is designed and fabricated which reaches a no-load temperature of 2.5K in the second stage and similar to 60 K in the first stage respectively. The system provides a cooling power of similar to 250 mW at 5K in the second stage. Stainless steel meshes (size 200) and lead (Pb) granules are used as the first stage regenerator materials and combination of Pb with Er3Ni / HoCu2 are used as the second stage regenerator materials. The system operates at 1.6 Hz using a 6 kW water cooled helium compressor. Studies conducted by varying the dimensions of Pulse Tubes and regenerators show that the dimensions of the Pulse Tubes are more critical to the performance of the Cryocooler than those of the regenerators. Experimental studies show that the optimum volume ratios of Er3Ni to Pb and HoCu2 to Pb in the second stage regenerator should be 3:2 and 2:3 respectively for the best performance. Further, systems with HoCu2 performed better than those with Er3Ni. The theoretical analysis of the system has been carried out using a simple isothermal model. The experimentally measured cooling powers are in good agreement with the theoretical predictions.

Investigations of a two-stage pulse tube cryocooler operating down to 2.5K

A two-stage pulse tube cryocooler (PTC) which produces a no-load temperature of similar to 2.5 K in its second stage at an operating frequency of 1.6 Hz has been designed and fabricated. The second stage of the system provides a refrigeration power of similar to 250 mW at 5.0 K. The system uses stainless steel meshes (mesh size 200) along with lead (Pb) granules and combinations of Pb, Er3Ni, and HoCu2 as the first and second stage regenerator materials, respectively. Experimental studies have been carried out on different pulse tube configurations by varying the dimensions of the pulse tubes and regenerators to arrive at the best one, which leads to the lowest no-load second stage cold head temperature. Using this configuration, detailed experimental studies have been conducted by varying the volume percentage ratios of the second stage regenerator materials such as HoCu2, Er3Ni, and Pb (with an average grain size of similar to 250 mu m). This article presents the results of our experimental studies on cryocoolers with the regenerator material arranged in layered structures. Comparative studies have also been presented for specific cases where the regenerator materials are arranged as a homogeneous mixture in the second stage. The experimental results clearly indicate that the design of PTCs should use only layered structures of regenerator materials and not homogenous mixtures.

Non-periodic tilings in 2-dimensions with 4, 6, 8, 10 and 12-fold symmetries

The two dimensional plane can be filled with rhombuses, so as to generate non-periodic tilings with 4, 6, 8, 10 and 12-fold symmetries. Some representative tilings constructed using the rule of inflation are shown. The numerically computed diffraction patterns for the corresponding tilings are also shown to facilitate a comparison with possible X-ray or electron diffraction pictures.

Studies in biogas technology. Part II. Optimisation of plant dimensions

In this paper, the design basis of the conventional Khadi and Village Industries Commission biogas plants has been elucidated. It has been shown that minimisation of the cost of the gas holder alone leads to the narrow and deep digesters of conventional plants. If instead, the total capital cost of the gas holder plus digester is minimised, the optimisation leads to wide and shallow digesters, which are less expensive. To test this alternative, two prototype plants have been designed, constructed and operated. These plants are not only 25–40% cheaper, but their performance is actually slightly better than the conventional plants.

Colossal nonsaturating linear magnetoresistance in two-dimensional electron systems at a GaAs/(Al,Ga)As heterointerface

Engineering devices with a large electrical response to magnetic field is of fundamental importance for a range of applications such as magnetic field sensing and magnetic read heads. We show that a colossal nonsaturating linear magnetoresistance (NLMR) arises in two-dimensional electron systems hosted in a GaAs/AlGaAs heterostructure in the strongly insulating regime. When operated at high source-drain bias, the magnetoresistance of our devices increases almost linearly with magnetic field, reaching nearly 10 000% at 8 T, thus surpassing many known nonmagnetic materials that exhibit giant NLMR. The temperature dependence and mobility analysis indicate that the NLMR has a purely classical origin, driven by nanoscale inhomogeneities. A large NLMR combined with small device dimensions makes these systems an attractive candidate for on-chip magnetic field sensing.

Representing a cubic graph as the intersection graph of axis-parallel boxes in three dimensions

We show that every graph of maximum degree 3 can be represented as the intersection graph of axis parallel boxes in three dimensions, that is, every vertex can be mapped to an axis parallel box such that two boxes intersect if and only if their corresponding vertices are adjacent. In fact, we construct a representation in which any two intersecting boxes just touch at their boundaries. Further, this construction can be realized in linear time.

REPRESENTING A CUBIC GRAPH AS THE INTERSECTION GRAPH OF AXIS-PARALLEL BOXES IN THREE DIMENSIONS

We show that every graph of maximum degree 3 can be represented as the intersection graph of axis parallel boxes in three dimensions, that is, every vertex can be mapped to an axis parallel box such that two boxes intersect if and only if their corresponding vertices are adjacent. In fact, we construct a representation in which any two intersecting boxes touch just at their boundaries.

pentahexoctite: A new two-dimensional allotrope of carbon

The ability of carbon to exist in many forms across dimensions has spawned search in exploring newer allotropes consisting of either, different networks of polygons or rings. While research on various 3D phases of carbon has been extensive, 2D allotropes formed from stable rings are yet to be unearthed. Here, we report a new sp(2) hybridized two-dimensional allotrope consisting of continuous 5-6-8 rings of carbon atoms, named as ``pentahexoctite''. The absence of unstable modes in the phonon spectra ensures the stability of the planar sheet. Furthermore, this sheet has mechanical strength comparable to graphene. Electronically, the sheet is metallic with direction-dependent flat and dispersive bands at the Fermi level ensuring highly anisotropic transport properties. This sheet serves as a precursor for stable 1D nanotubes with chirality-dependent electronic and mechanical properties. With these unique properties, this sheet becomes another exciting addition to the family of robust novel 2D allotropes of carbon.