Rate effect and coupled evolution of atomic motions and potential landscapes

Since rate effect of materials plays a key role in impact engineering, the microscopic mechanism of rate effect is investigated at molecular level in this paper. The results show that rate effect on the strength of atomic system is closely related to the coupled evolution of atomic motions and potential landscapes. Accordingly, it becomes possible to develop a new algorithm of molecular simulation, which could properly and efficiently demonstrate strain rate effect under a wide range of loading rates and unveil the mecha- nisms underlying the strain rate effects.     

Calculating particle pair potentials from fluid‐state pair correlations: Iterative ornstein–zernike inversion

An iterative Monte Carlo inversion method for the calculation of particle pair potentials from given particle pair correlations is proposed in this article. The new method, which is best referred to as Iterative Ornstein–Zernike Inversion, represents a generalization and an improvement of the established Iterative Boltzmann Inversion technique (Reith, Pütz and Müller‐Plathe, J. Comput. Chem. 2003, 24, 1624). Our modification of Iterative Boltzmann Inversion consists of replacing the potential of mean force as an approximant for the pair potential with another, generally more accurate approximant that is based on a trial bridge function in the Ornstein–Zernike integral equation formalism. As an input, the new method requires the particle pair correlations both in real space and in the Fourier conjugate wavenumber space. An accelerated iteration method is included in the discussion, by which the required number of iterations can be greatly reduced below that of the simple Picard iteration that underlies most common implementations of Iterative Boltzmann Inversion. Comprehensive tests with various pair potentials show that the new method generally surpasses the Iterative Boltzmann Inversion method in terms of reliability of the numerical solution for the particle pair potential. © 2018 Wiley Periodicals, Inc.     

(Adiabatic) phase boundaries in a bistable chain with twist and stretch

Mass–spring chains with only extensional degrees of freedom have provided insights into the behavior of crystalline solids, including those capable of phase transitions. Here we add rotational degrees of freedom to the masses in a chain and study the dynamics of phase boundaries across which both the twist and stretch can jump. We solve impact and Riemann problems in the chain by numerical integration of the equations of motion and show that the solutions are analogous to those in a phase transforming rod whose stored energy function depends on both twist and stretch. From the dynamics of phase boundaries in the chain we extract a kinetic relation whose form is familiar from earlier studies involving chains with only extensional degrees of freedom. However, for some combinations of parameters characterizing the energy landscape of our springs we find propagating phase boundaries for which the rate of dissipation, as calculated using isothermal expressions for the driving force, is negative. This suggests that we cannot neglect the energy stored in the oscillations of the masses in the interpretation of the dynamics of mass–spring chains. Keeping this in mind we define a local temperature of our chain and show that it jumps across phase boundaries, but not across sonic waves. Hence, impact problems in our mass–spring chains are analogous to those on continuum thermoelastic bars with Mie–Gruneisen type constitutive laws. At the end of the paper we use our chain to shed some light on experiments involving yarns that couple twist and stretch to perform useful work in response to various stimuli.     

Grid‐independent depth‐averaged simulations with a collocated unstructured finite volume scheme

In this article, the depth‐averaged transport equations are written in a new way so that it is possible to solve the transport equations for very small water depths. Variables are interpolated into the cell face with two different schemes and, the schemes are compared in terms of computational cost and accuracy. The bed source terms are computed using two different assumptions. The effect of these assumptions on numerical simulations is then investigated. Solutions of transport equations on different types of unstructured triangular grids are compared and, an appropriate choice of grid is suggested. Copyright © 2011 John Wiley & Sons, Ltd.     

Correlation of thermodynamic modeling and molecular simulations for liquid-liquid equilibrium of ternary polymer mixtures based on a phenomenological scaling method

In our previous study [S.Y. Oh, Y.C. Bae, J. Phys. Chem. B 114 (2010) 8948-8953], we presented a new method to predict liquid-liquid equilibria in ternary simple liquid mixtures by using a combination of a thermodynamic model and molecular simulations. As a continuation of that effort, we extend our previously developed method to ternary polymer systems. In the simulations, we used the dummy atoms to calculate the pair interaction energy values between the polymer segments and the solvent molecules. Furthermore, a thermodynamic model scaling concept is introduced to consider the chain length dependence of the energy parameters. This method was applied to ternary mixtures incorporating low to high molecular weight polymers. The method presented here well described the experimental observations using one or no adjustable parameters.     

A boundary condition for arbitrary shaped inlets in lattice–Boltzmann simulations

We introduce a mass‐flux‐based inlet boundary condition for the lattice‐Boltzmann method. The proposed boundary condition requires minimal amount of boundary data, it produces a steady‐state velocity field which is accurate close to the inlet even for arbitrary inlet geometries, and yet it is simple to implement. We demonstrate its capability for both simple and complex inlet geometries by numerical experiments. For simple inlet geometries, we show that the boundary condition provides very accurate inlet velocities when Re?1. Even with moderate Reynolds number, the inlet velocities are accurate for practical purposes. Furthermore, the potential of our boundary condition to produce inlet velocities which convincingly adapt to complex inlet geometries is highlighted with two specific examples. Copyright © 2009 John Wiley & Sons, Ltd.     

Internal flow numerical simulation of a cross‐flow fan in refrigerant operating condition using user‐defined functions

In the paper, a cross‐flow fan in refrigerant operating condition is systematically simulated using user‐defined functions. Three‐dimensional simulations are acquired with Navier–Stokes equations coupled with k–ε turbulence model, and internal flow characteristics of an indoor split‐type air conditioner are obtained, which is mainly composed of cross‐flow fan and heat exchanger. It has systematically been simulated in the isothermal flow condition that the performance of cross‐flow fan may be reduced easily with dry or humid air, and in the refrigerant operating condition in which user‐defined functions are applied to the humid air, considered as a mixture of dry air and vapor. A density‐modulated function is adopted to deal with the condensation of the vapor at the heat‐transfer region approximately. The results show flow mechanism of the two gas‐phase flow, including phase‐vary process. The distribution of the parameters is not uniform at the inlet of the machine, the intensity and position of pressure and velocity vary along the axial direction of the fan, the distribution of vapor volume fraction and turbulent intensity in heat‐transfer region is obtained, and the external characteristic data of the indoor machine are obtained and analyzed. Compared with the experimental data, the calculated characteristic curves and designed parameters are on target. © British Crown Copyright 2010/MOD. Reproduced with permission. Published by John Wiley & Sons, Ltd.     

Comparison of non‐reflective boundary conditions for a free‐rising turbulent axisymmetric plume

The stability and accuracy of radiation type non‐reflective outflow boundary conditions, as well as the standard Neumann boundary condition with zero normal derivative, have been compared for the numerical simulation of a turbulent axisymmetric plume with Reynolds number of 7700 and Prandtl number of 0.71. Comparison of the performance of the boundary conditions with respect to each other, and to the results obtained for an extended domain, shows that a one‐dimensional scheme in which advection and diffusion terms are included in the radiation equation is the optimum approach for the plume simulation. Copyright © 2013 John Wiley & Sons, Ltd.     

Electric field assisted annealing effects on microstructure and ionic conductivity in ceria/YSZ oxide heterostructures

The effect of electric field assisted annealing on the microstructure, composition and ionic conductivity properties in CeO₂/YSZ oxide heterostructures have been investigated using molecular dynamics simulations. Amorphization–recrystallization steps were performed with and without external electric field of strength 10?MV/cm along three different orientations: in-plane (YZ), normal (X) and 45° resultant (XY) with respect to the oxide heterointerfaces. The microstructural and compositional differences at the interfaces and in the interior of the oxide heterolayers were evaluated and were found to show a clear correlation with the orientations of the applied field. In particular, the XY configuration displayed a compressive lattice strain of ～9% along with a reduced oxygen vacancy concentration when compared to the others. Ionic density profiles suggest pronounced segregation (～60% higher compared to the average value in the interior) of yttrium ions closer to the YSZ/CeO₂ interface for the XY configuration. Other configurations exhibit minimal to no such variations. These microstructural differences are found to affect the number of mobile charge carriers and the activation barriers associated with ionic migration through the oxide lattice and consequently, influence the ionic conductivity.     

Simulation of metastable zone width and induction time for a seeded aqueous solution of potassium sulfate

The metastable zone width (MSZW, ΔT_m) and induction time (t_ind) were determined with computer simulation for seeded batch crystallization of potassium sulfate from aqueous solution. The MSZW and induction time determined with simulation showed the same behavior as experimental values reported in the literature; log (ΔT_m) increased linearly with an increase in log R (R: cooling rate) and t_ind decreases in proportion to (ΔT)⁻ⁿ (ΔT: supercooling, n: nucleation order in the secondary rate expression of B=k_n(ΔT)ⁿ). The secondary nucleation parameters (k_n and n) were deduced both from the simulated MSZW and induction times by using the previously proposed model [J. Cryst. Growth, 2010, 312, 548–554]. The secondary nucleation rate calculated with the deduced parameters was in agreement with that calculated with the parameters input for simulation.