Volume effects and associations in liquid alloys

The regular associated solution model for binary systems has been modified by incorporating the size of the complex as an explicit variable. The thermodynamic properties of the liquid alloy and the interactions between theA ?B type of complex and the unassociated atoms in anA-B binary have been evaluated as a function of relative size of the complex using the activity coefficients at infinite dilution and activity data at one other composition in the binary. The computational procedure adopted for determining the concentration of clusters and interaction energies in the associated liquid is similar to that proposed by Lele and Rao. The analysis has been applied to the thermodynamic mixing functions of liquid Al-Ca alloys believed to contain Al2Ca associates. It is found that the size of the cluster significantly affects the interaction energies between the complex and the unassociated atoms, while the equilibrium constant and enthalpy change for the association reaction exhibit only minor variation, when the equations are fitted to experimental data. The interaction energy between unassociated free atoms remains virtually unaltered as the size of the complex is varied between extreme values. Accurate data on free energy, enthalpy, and volume of mixing at the same temperature on alloy systems with compound forming tendency would permit a rigorous test of the proposed model.

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.

Model-free simulations for compressible mixing layer

Confined supersonic mixing layer is explored through model-free simulations. Both two- and three-dimensional spatio-temporal simulations were carried out employing higher order finite difference scheme as well as finite volume scheme based on open source software (OpenFOAM) to understand the effect of three-dimensionality on the development of mixing layer. It is observed that although the instantaneous structures exhibit three-dimensional features, the average pressure and velocities are predominantly two-dimensional. The computed wall pressures match well with experimental results fairly well, although three-dimensional simulation underpredicts the wall pressure in the downstream direction. The self-similarity of the velocity profiles is obtained within the duct length for all the simulations. Although the mixing layer thicknesses differ among different simulations, their growth rate is nearly the same. Significant differences are observed for species and temperature distribution between two- and three-dimensional calculations, and two-dimensional calculations do not match the experimental observation of smooth variations in species mass fraction profiles as reported in literature. Reynolds stress distribution for three-dimensional calculations show profiles with less peak values compared to two-dimensional calculations; while normal stress anisotropy is higher for three-dimensional case.

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.

Isotropic finite volume discretization

Finite volume methods traditionally employ dimension by dimension extension of the one-dimensional reconstruction and averaging procedures to achieve spatial discretization of the governing partial differential equations on a structured Cartesian mesh in multiple dimensions. This simple approach based on tensor product stencils introduces an undesirable grid orientation dependence in the computed solution. The resulting anisotropic errors lead to a disparity in the calculations that is most prominent between directions parallel and diagonal to the grid lines. In this work we develop isotropic finite volume discretization schemes which minimize such grid orientation effects in multidimensional calculations by eliminating the directional bias in the lowest order term in the truncation error. Explicit isotropic expressions that relate the cell face averaged line and surface integrals of a function and its derivatives to the given cell area and volume averages are derived in two and three dimensions, respectively. It is found that a family of isotropic approximations with a free parameter can be derived by combining isotropic schemes based on next-nearest and next-next-nearest neighbors in three dimensions. Use of these isotropic expressions alone in a standard finite volume framework, however, is found to be insufficient in enforcing rotational invariance when the flux vector is nonlinear and/or spatially non-uniform. The rotationally invariant terms which lead to a loss of isotropy in such cases are explicitly identified and recast in a differential form. Various forms of flux correction terms which allow for a full recovery of rotational invariance in the lowest order truncation error terms, while preserving the formal order of accuracy and discrete conservation of the original finite volume method, are developed. Numerical tests in two and three dimensions attest the superior directional attributes of the proposed isotropic finite volume method. Prominent anisotropic errors, such as spurious asymmetric distortions on a circular reaction-diffusion wave that feature in the conventional finite volume implementation are effectively suppressed through isotropic finite volume discretization. Furthermore, for a given spatial resolution, a striking improvement in the prediction of kinetic energy decay rate corresponding to a general two-dimensional incompressible flow field is observed with the use of an isotropic finite volume method instead of the conventional discretization. (C) 2014 Elsevier Inc. All rights reserved.

An improved free surface modeling for incompressible SPH

In Incompressible Smooth Particle Hydrodynamics (ISPH), a pressure Poisson equation (PPE) is solved to obtain a divergence free velocity field. When free surfaces are simulated using this method a Dirichlet boundary condition for pressure at the free surface has to be applied. In existing ISPH methods this is achieved by identifying free surface particles using heuristically chosen threshold of a parameter such as kernel sum, density or divergence of the position, and explicitly setting their pressure values. This often leads to clumping of particles near the free surface and spraying off of surface particles during splashes. Moreover, surface pressure gradients in flows where surface tension is important are not captured well using this approach. We propose a more accurate semi-analytical approach to impose Dirichlet boundary conditions on the free surface. We show the efficacy of the proposed algorithm by using test cases of elongation of a droplet and dam break. We perform two dimensional simulations of water entry and validate the proposed algorithm with experimental results. Further, a three dimensional simulation of droplet splash is shown to compare well with the Volume-of-Fluid simulations. (C) 2014 Elsevier Ltd. All rights reserved.

Sphere to ring morphological transformation in drying nanofluid droplets in a contact-free environment

Understanding the transients of buckling in drying colloidal suspensions is pivotal for producing new functional microstructures with tunable morphologies. Here, we report first observations and elucidate the buckling instability induced morphological transition (sphere to ring structure) in an acoustically levitated, heated nanosuspension droplet using dynamic energy balance. Droplet deformation featuring the formation of symmetric cavities is initiated by capillary pressure that is two to three orders of magnitude greater than the acoustic radiation pressure, thus indicating that the standing pressure field has no influence on the buckling front kinetics. With an increase in heat flux, the growth rate of surface cavities and their post-buckled volume increase while the buckling time period reduces, thereby altering the buckling pathway and resulting in distinct precipitate structures. However, irrespective of the heating rate, the volumetric droplet deformation exhibits a linear time dependence and the droplet vaporization is observed to deviate from the classical D-2-law.

Volume conservation issues in incompressible smoothed particle hydrodynamics

A divergence-free velocity field is usually sought in numerical simulations of incompressible fluids. We show that the particle methods that compute a divergence-free velocity field to achieve incompressibility suffer from a volume conservation issue when a finite time-step position update scheme is used. Further, we propose a deformation gradient based approach to arrive at a velocity field that reduces the volume conservation issues in free surface flows and maintains density uniformity in internal flows while retaining the simplicity of first order time updates. (C) 2015 Elsevier Inc. All rights reserved.

Interrelationship between Interphase Boundaries and Phase Contents near the Critical Compositions of Lead-Free Ferroelectric (Na0.5Bi0.5)TiO3-BaTiO3

Heterophase structures in lead-free perovskite-type ferroelectric solid solutions of (1 - z)(Na0.5Bi0.5)TiO3 - zBaTiO(3) are analysed for a few critical compositions near the morphotropic phase boundary (z = 0.05-0.07). Examples of the phase coexistence and elastic matching of the phases from different symmetry groups are considered to find optimum volume fractions of specific domain types and coexisting phases at the complete stress relief in two-phase samples. Some interrelations between these volume fractions are described using variants of the domain arrangement at changes in the composition and unit-cell parameters. The evaluated room-temperature volume fractions of the ferroelectric monoclinic (Cm symmetry) and tetragonal (P4mm symmetry) phases near the morphotropic phase boundary are in agreement with experimental data.

Limited-view light sheet fluorescence microscopy for three dimensional volume imaging

We propose and demonstrate a limited-view light sheet microscopy (LV-LSM) for three dimensional (3D) volume imaging. Realizing that longer and frequent image acquisition results in significant photo-bleaching, we have taken limited angular views (18 views) of the macroscopic specimen and integrated with maximum likelihood (ML) technique for reconstructing high quality 3D volume images. Existing variants of light-sheet microscopy require both rotation and translation with a total of approximately 10-fold more views to render a 3D volume image. Comparatively, LV-LSM technique reduces data acquisition time and consequently minimizes light-exposure by many-folds. Since ML is a post-processing technique and highly parallelizable, this does not cost precious imaging time. Results show noise-free and high contrast volume images when compared to the state-of-the-art selective plane illumination microscopy. (C) 2015 AIP Publishing LLC.

Microstructural characterization of ultrafine-grain interstitial-free steel by X-ray diffraction line profile analysis

This paper highlights the microstructural features of commercially available interstitial free (IF) steel specimens deformed by equal channel angular pressing (ECAP) up to four passes following the route A. The microstructure of the samples was studied by different techniques of X-ray diffraction peak profile analysis as a function of strain (epsilon). It was found that the crystallite size is reduced substantially already at epsilon=2.3 and it does not change significantly during further deformation. At the same time, the dislocation density increases gradually up to epsilon=4.6. The dislocation densities estimated from X-ray diffraction study are found to correlate very well with the experimentally obtained yield strength of the samples.

Singularity-free path planning for the Stewart platform manipulator

This paper addresses the problem of singularity-free path planning for the six-degree-of-freedom parallel manipulator known as the Stewart platform manipulator. Unlike serial manipulators, the Stewart platform possesses singular configurations within the workspace where the manipulator is uncontrollable. An algorithm has been developed to construct continuous paths within the workspace of the manipulator by avoiding singularities and ill-conditioning. Given two end-poses of the manipulator, the algorithm finds out safe (well-conditioned) via points and plans a continuous path from the initial pose to the final one. When the two end-poses belong to different branches and no singularity-free path is possible, the algorithm indicates the impossibility of a valid path. A numerical example has also been presented as illustration of the path planning strategy.

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.

Unsteady free convection flow in the stagnation-point region of a rotating sphere

The unsteady free convection boundary-layer flow in the forward stagnation-point region of a sphere, which is rotating with time-dependent angular velocity in an ambient fluid, has been studied. Both constant wall temperature and constant hear flux conditions have been considered. The non-linear coupled parabolic partial differential equations governing the flow have been solved numerically using an implicit finite-difference scheme. The skin friction and the heat transfer are enhanced by the buoyancy force. The effect of the buoyancy force is found to be more pronounced for smaller Prandtl numbers than for larger Prandtl numbers. For a given buoyancy force, the heat transfer increases with an increase in Prandtl number, but the skin friction decreases.

Unsteady free convection flow in the stagnation-point region of a three-dimensional body

The unsteady free convection flow in the stagnation-point region of a heated three-dimensional body placed in an ambient fluid is studied under boundary layer approximations. We have considered the case where there is an initial steady state that is perturbed by a step-change in the wall temperature. The non-linear coupled partial differential equations governing the free convection flow are solved numerically using a finite difference scheme. The presented results show the temporal development of the momentum and thermal boundary layer characteristics.