Molecular/cluster statistical thermodynamics methods to simulate quasi-static deformations at finite temperature

The rapid evolution of nanotechnology appeals for the understanding of global response of nanoscale systems based on atomic interactions, hence necessitates novel, sophisticated, and physically based approaches to bridge the gaps between various length and time scales. In this paper, we propose a group of statistical thermodynamics methods for the simulations of nanoscale systems under quasi-static loading at finite temperature, that is, molecular statistical thermodynamics (MST) method, cluster statistical thermodynamics (CST) method, and the hybrid molecular/cluster statistical thermodynamics (HMCST) method. These methods, by treating atoms as oscillators and particles simultaneously, as well as clusters, comprise different spatial and temporal scales in a unified framework. One appealing feature of these methods is their “seamlessness” or consistency in the same underlying atomistic model in all regions consisting of atoms and clusters, and hence can avoid the ghost force in the simulation. On the other hand, compared with conventional MD simulations, their high computational efficiency appears very attractive, as manifested by the simulations of uniaxial compression and nanoindenation.     

A Continuum Theory of Deformable,Semiconducting Ferroelectrics

Ferroelectric solids, especially ferroelectric perovskites, are widely used as sensors, actuators, filters, memory devices, and optical components. While these have traditionally been treated as insulators, they are in reality wide-band-gap semiconductors. This semiconducting behavior affects the microstructures or domain patterns of the ferroelectric material and the interaction of ferroelectrics with electrodes, and is affected significantly by defects and dopants. In this paper, we develop a continuum theory of deformable, semiconducting ferroelectrics. A key idea is to introduce space charges and dopant density as field (state) variables in addition to polarization and deformation. We demonstrate the theory by studying oxygen vacancies in barium titanate. We find the formation of depletion layers, regions of depleted electrons, and a large electric field at the ferroelectric–electrode boundary. We also find the formation of a charge double layer and a large electric field across 90° domain walls but not across 180° domain walls. We show that these internal electric fields can give rise to a redistribution or forced diffusion of oxygen vacancies, which provides a mechanism for aging of ferroelectric materials.     

Nonlocal shell model for elastic wave propagation in single- and double-walled carbon nanotubes

This paper investigates the transverse and torsional wave in single- and double-walled carbon nanotubes (SWCNTs and DWCNTs), focusing on the effect of carbon nanotube microstructure on wave dispersion. The SWCNTs and DWCNTs are modeled as nonlocal single and double elastic cylindrical shells. Molecular dynamics (MD) simulations indicate that the wave dispersion predicted by the nonlocal elastic cylindrical shell theory shows good agreement with that of the MD simulations in a wide frequency range up to the terahertz region. The nonlocal elastic shell theory provides a better prediction of the dispersion relationships than the classical shell theory when the wavenumber is large enough for the carbon nanotube microstructure to have a significant influence on the wave dispersion. The nonlocal shell models are required when the wavelengths are approximately less than 2.36×10⁻⁹ and 0.95×10⁻⁹ m for transverse wave in armchair (15,15) SWCNT and torsional wave in armchair (10,10) SWCNT, respectively. Moreover, an MD-based estimation of the scale coefficient e₀ for the nonlocal elastic cylindrical shell model is suggested. Due to the small-scale effects of SWCNTs and the interlayer van der Waals interaction of DWCNTs, the phase difference of the transverse wave in the inner and outer tube can be observed in MD simulations in wave propagation at high frequency. However, the van der Waals interaction has little effect on the phase difference of transverse wave.     

A SIMPLE CONSTITUTIVE MODEL FOR FERROELECTRIC CERAMICS UNDER ELECTRICAL/MECHANICAL LOADING

A simple phenomenological model is developed for describing the macroscopic constitutive response of ferroelectric materials based on consideration of the fact that domain switching is a progressive evolution process with loading. The volume fraction of domain switching is taken as an internal variable, which is derived from the domain nucleation theory. The proposed theory can simulate the dielectric hysteresis, reversed butterfly hysteresis, nonlinear strain-stress hysteresis, as well as electric displacement-stress relation of ferroelectric materials. Its comparison with experimental results and two other theoretical models reveals that the model presented can well predict the nonlinear hysteresis of ferroelectrics under electrical or mechanical loading.     

Coupled Navier–Stokes—Molecular dynamics simulations using a multi‐physics flow simulation framework

Simulation of nano‐scale channel flows using a coupled Navier–Stokes/Molecular Dynamics (MD) method is presented. The flow cases serve as examples of the application of a multi‐physics computational framework put forward in this work. The framework employs a set of (partially) overlapping sub‐domains in which different levels of physical modelling are used to describe the flow. This way, numerical simulations based on the Navier–Stokes equations can be extended to flows in which the continuum and/or Newtonian flow assumptions break down in regions of the domain, by locally increasing the level of detail in the model. Then, the use of multiple levels of physical modelling can reduce the overall computational cost for a given level of fidelity. The present work describes the structure of a parallel computational framework for such simulations, including details of a Navier–Stokes/MD coupling, the convergence behaviour of coupled simulations as well as the parallel implementation. For the cases considered here, micro‐scale MD problems are constructed to provide viscous stresses for the Navier–Stokes equations. The first problem is the planar Poiseuille flow, for which the viscous fluxes on each cell face in the finite‐volume discretization are evaluated using MD. The second example deals with fully developed three‐dimensional channel flow, with molecular level modelling of the shear stresses in a group of cells in the domain corners. An important aspect in using shear stresses evaluated with MD in Navier–Stokes simulations is the scatter in the data due to the sampling of a finite ensemble over a limited interval. In the coupled simulations, this prevents the convergence of the system in terms of the reduction of the norm of the residual vector of the finite‐volume discretization of the macro‐domain. Solutions to this problem are discussed in the present work, along with an analysis of the effect of number of realizations and sample duration. The averaging of the apparent viscosity for each cell face, i.e. the ratio of the shear stress predicted from MD and the imposed velocity gradient, over a number of macro‐scale time steps is shown to be a simple but effective method to reach a good level of convergence of the coupled system. Finally, the parallel efficiency of the developed method is demonstrated. Copyright © 2009 John Wiley & Sons, Ltd.     

Phase-field simulation of the coupled evolutions of grain and twin boundaries in nanotwinned polycrystals

Nanotwinned polycrystals exhibit an excellent strength-ductility combination due to nanoscale twins and grains. However, nanotwin-assisted grain coarsening under mechanical loading reported in recent experiments may result in strength drop based on the Hall-Petch law. In this paper, a phase-field model is developed to investigate the effect of coupled evolutions of twin and grain boundaries on nanotwin-assisted grain growth. The simulation result demonstrates that there are three pathways for coupled motions of twin and grain boundaries in a bicrystal under the applied loading, dependent on the amplitude of applied loading and misorientation of the bicrystal. It reveals that a large misorientation angle and a large applied stress promote the twinning-driven grain boundary migration. The resultant twin-assisted grain coarsening is confirmed in the simulations for the microstructural evolutions in twinned and un-twinned polycrystals under a high applied stress.     

An energy-consistent fracture model for ferroelectrics

The fracture behavior of ferroelectrics has been intensively studied in recent decades, though currently a widely accepted fracture mechanism is still lacking. In this work, enlightened by previous experimental observations that crack propagation in ferroelectrics is always accompa-nied by domain switching, we propose a micromechanical model in which both crack propagation and domain switch-ing are controlled by energy-based criteria. Both electric energy and mechanical energy can induce domain switching, while only mechanical energy can drive crack propagation. Furthermore, constrained domain switching is considered in this model, leading to the gradient domain switching zone near the crack tip. Analysis results show that stress-induced ferroelastic switching always has a toughening effect as the mechanical energy release rate serves as the driving force for both fracture and domain switching. In compari-son, the electric-field-induced switching may have either a toughening or detoughening effect. The proposed model can qualitatively agree with the existing experimental results.     

A constrained particle dynamics for continuum-particle hybrid method in micro- and nano-fluidics

A hybrid method of continuum and particle dynamics is developed for micro- and nano-fluidics, where fluids are described by a molecular dynamics (MD) in one domain and by the Navier–Stokes (NS) equations in another domain. In order to ensure the continuity of momentum flux, the continuum and molecular dynamics in the overlap domain are coupled through a constrained particle dynamics. The constrained particle dynamics is constructed with a virtual damping force and a virtual added mass force. The sudden-start Couette flows with either non-slip or slip boundary condition are used to test the hybrid method. It is shown that the results obtained are quantitatively in agreement with the analytical solutions under the non-slip boundary conditions and the full MD simulations under the slip boundary conditions.The project supported by Chinese Academy of Sciences under the innovative project “Multi-scale modelling and simulation in complex system” (KJCX-SW-L08) and National Natural Science Foundation of China (10325211).     

Coupled Navier–Stokes/molecular dynamics simulations in nonperiodic domains based on particle forcing

This paper focuses on coupling methods for hybrid Navier–Stokes/molecular dynamics (MD) simulations. The computational domain is split in a continuum flow region, where a finite‐volume discretisation of the Navier–Stokes equations is used, and one or more particle domains, where molecular level modelling of the flow is employed. The domains are defined with a partial overlap, in which the flow states are coupled through an exchange of the velocity components. For the steady flows considered, an under‐relaxed Newton iteration method is used to drive the coupled system to convergence. The main focus of the present work is on methods to impose nonperiodic boundary conditions on the particle domain(s). A particle forcing is applied in the direction normal to the particle domain boundary to impose the boundary normal velocity component. A novel aspect of the present work is the extension of this method to more general nonplanar particle domain boundaries. The main contribution of the paper is the development of a particle forcing method in the direction tangential to the domain boundary, which is based on the equivalent continuum‐flow boundary shear stresses along with an iterative forcing strength adjustment based on the extrapolated particle boundary velocity. Furthermore, an adaptation scheme is presented, which uses the finite‐volume flux residuals of the particle bin averaged velocity field as a truncation criterion for the iterative force‐update scheme. It is demonstrated that by comparing the residual reduction for the momentum equation in the nonhomogeneous directions during the molecular dynamics simulations with that for a homogeneous direction, the forcing iteration at which the statistical noise in the velocity field dominates the uncertainty in the forcing strength can be determined. At this point the iteration can be truncated. It is shown that with adaptive schemes of this type, the total number of MD evaluations required in a coupled Navier–Stokes/MD simulation can be reduced relative to a hybrid scheme with a fixed number of forcing‐strength updates. Copyright © 2011 John Wiley & Sons, Ltd.     

On the physical interpretation of material parameters in phase field models for ferroelectrics

This publication intends to clarify the relations between the free energy expansion and the choice of the order parameter in phase field models for ferroelectrics. The free energy is expanded to include the piezoelectric coupling terms, reflecting the symmetry properties of the material below the Curie temperature. The order parameter can then be interpreted as the material polarization less than the dielectric and piezoelectric contributions. Furthermore, a 1d reduction of the problem is used to identify the physical meaning of the remaining model parameters. The theoretical results are verified by numerical simulations of 180^? and 90^? domain walls. Both the static case, where the emphasis lies on the specific energy and the interface width, and the dynamic case, where interfaces are subjected to an electric driving field, are considered in the simulations.     

An I-integral method for crack-tip intensity factor variation due to domain switching in ferroelectric single-crystals

In the present study, an I-integral method is established for solving the crack-tip intensity factors of ferroelectric single-crystals. The I-integral combined with the phase field model is successfully used to investigate crack-tip intensity factor variations due to domain switching in ferroelectricity subjected to electromechanical loadings, which exhibits several advantages over previous methods based on small-scale switching. First, the shape of the switching zone around a crack tip is predicted by the time-dependent Ginzburg–Landau equation, which does not require preset energy-based switching criterion. Second, the I-integral can directly solve the crack-tip intensity factors and decouple the crack-tip intensity factors of different modes based on superimposing an auxiliary state onto an actual state. Third, the I-integral is area-independent, namely, the I-integral is not affected by the integral area size, the polarization distributions, or domain walls. This makes the I-integral applicable to large-scale domain switching. To this end, the electro-elastic field intensity factors of an impermeable crack in PbTiO₃ ferroelectric single crystals are evaluated under electrical, mechanical, and combined loading. The intensity factors obtained by the I-integral agree well with those obtained by the extrapolation technique. From numerical results, the following conclusions can be drawn with respect to fracture behavior of ferroelectrics under large-scale switching. Under displacement controlled mechanical loading, the stress intensity factors (SIFs) decrease monotonically due to the domain switching process, which means a crack tip shielding or effective switching-induced toughening occurs. If an external electric field is applied, the electric displacement intensity factor (EDIF) increases in all cases, i.e., the formed domain patterns enhance the electric crack tip loading. The energy release rate, expressed by the crack-tip J-integral, is reduced by the domain switching in all examples, which underlines the switching-induced-toughening effect. In contrast, under stress controlled load, the SIF evolves due to large-scale switching to a stable value, which is higher than the non-switching initial value, i.e., fracture is promoted in this case.     

GIANT ANHYSTERETIC RESPONSE OF FERROELECTRIC SOLID SOLUTIONS WITH MORPHOTROPIC BOUNDARIES: THE ROLE OF POLAR ANISOTROPY

Computer modeling and simulation for the Pb(Zr 1 x Ti x )O 3 (PZT) system reveal the role of polar anisotropy on the giant anhysteretic response and structural properties of morphotropic phase boundary (MPB) ferroelectrics. It is shown that a drastic reduction of the composition- dependent polar anisotropy near the MPB flattens energy surfaces and thus facilitates reversible polarization rotation. It is further shown that the polar anisotropy favors formation of polar domains, promotes phase decomposition and results in a two-phase multidomain state, which response to applied electric field is anhysteretic when the polar domain reorientation is only caused by polarization rotation other than polar domain wall movement. This is the case for the decomposing ferroelectrics under a poling electric field with the formation of a two-phase multidomain microstructure, wherein most domain walls are pinned at the two-phase boundaries. Indication of the microstructure dependence of the anhysteretic strain response opens new avenues to improve the piezoelectric properties of these materials through the microstructure engineering.     

A continuum analysis of charge induced ferroelectric domain wall motions

A continuum model is presented for the motion of a domain wall in a plane 90°-domain configuration subjected to an isolated extrinsic charge near the surface of a ferroelectric single crystal. Local pinning is postulated for the kinetic law. Before the appearance of the extrinsic charge, all polarization surface charges are taken to be neutralized by environmental charges. The domain wall motion after the appearance of the extrinsic charge is assumed to proceed sufficiently fast without any significant conductive currents on the surface or in the interior of the crystal such that new surface and interface polarization charges remain unscreened and contribute to the ferroelectric anisotropy energy. A non-admissible divergence of the electric field and consequently of the local thermodynamic driving force and of the domain wall velocity appears in the model if the domain wall charged by interface polarization charges intersects the crystal surface charged by surface polarization charges under an arbitrary angle. The physically possible domain wall angle is identified using the condition of a non-divergent driving force. The ferroelectric anisotropy energy and an intrinsic surface energy of the domain wall, however, do not provide stability of the domain wall trajectory against an unlimited increase of its curvature at the surface. The problem has been solved conceptually by proper account of the domain wall bending energy. Numerical and dimensional analysis explain also why domain walls driven by extrinsic charges remain almost straight in soft ferroelectrics.     

A nonlocal continuum model for the buckling of carbon honeycombs

Using molecular dynamics (MD) simulations and Eringen’s nonlocal elasticity theory, in this paper we comprehensively study the small-scale effects on the buckling behaviours of carbon honeycombs (CHCs). The MD simulation results show that the small-scale effects stemming from the long-range van der Waals interaction between carbon atoms can significantly affect the buckling behaviours of CHCs. To incorporate the small-scale effects into the theoretical analysis of the buckling of CHCs, we develop a nonlocal continuum mechanics (CM) model by employing Eringen’s nonlocal elasticity theory. Our nonlocal CM model is found to fit MD simulations well by setting the nonlocal parameter in the nonlocal CM model as 0.67. It is shown in our MD-based nonlocal CM model that when the cell length of CHCs is smaller than 7.93 Å the influence of small-scale effects on the bucking of CHCs becomes unnegligible and the small-scale effects can greatly reduce the critical buckling stress of CHCs. This reduction in critical buckling stress induced by the small-scale effects becomes more significant as the length of the cell wall decreases. Moreover, CHCs are found to display two different buckling modes when they are under different states of loading. The critical condition for the transition between these two buckling modes of CHCs can be greatly affected by the small-scale effects when the vertical cell wall and the inclined cell wall of CHCs have different lengths.     

Feasibility study of symmetric solution of molecular argon flow inside microscale nozzles

The computational cost of numerical methods in microscopic-scales such as molecular dynamics (MD) is a deterrent factor that limits simulations with a large number of particles. Hence, it is desirable to decrease the computational cost and run time of simulations, especially for problems with a symmetrical domain. However, in microscopic-scales, implementation of symmetric boundary conditions is not straight-forward. Previously, the present authors have successfully used a symmetry boundary condition to solve molecular flows in constant-area channels. The results obtained with this approach agree well with the benchmark cases. Therefore, it has provided us with a sound ground to further explore feasibility of applying symmetric solutions of micro-fluid flows in other geometries such as variable-area ducts. Molecular flows are solved for the whole domain with and without the symmetric boundary condition. Good agreement has been reached between the results of the symmetric solution and the whole domain solution. To investigate robustness of the proposed method, simulations are conducted for different values of affecting parameters including an external force, a flow density, and a domain length. The results indicate that the symmetric solution is also applicable to variable-area ducts such as micro-nozzles.     

Modeling and simulation of large-scale ductile fracture in plates and shells

The work is concerned with the modeling and simulation of large scale ductile fracture in plate and shell structures. A meshfree method – the reproducing kernel particle method (RKPM) – is used in numerical computations in order to enact dynamic crack propagation without remeshing. There are several novelties in the present approach. First, we have developed a crack surface approximation and particle split algorithm for three-dimensional through-thickness cracks. Second, to represent evolving crack surface in 3D shell structures, a 3D parametric visibility condition algorithm is proposed, which re-constructs the local connectivity map for particles near the crack tip or crack surfaces, so that the meshfree interpolation field can represent physical material separation in the computational domain. Third, the constitutive update formulas in explicit time integration by different versions of Gurson models and the rate-dependent Johnson–Cook model are implemented for 3D computations. Finally, the performance of different Gurson-type models are investigated and compared with the experimental data of large scale in-plate tear process. Numerical simulations of crack propagation in stiffened plates and shells demonstrate that the proposed method provides an effective means to simulate ductile fracture in large scale plate/shell structures with engineering accuracy.     

A SIMPLIFIED MODEL FOR DOMAIN SWITCHING OF FERROELECTRIC CRYSTAL

I. INTRODUCTION Domain switching is the main source of nonlinear characteristics of ferroelectric materials. The trans-formation performance of domain is the basis for ferroelectric constitutive research. In many literatures[1??10], complete switching models were adopted, in which ferroelectric materials are considered tobe consist of several basic sorts of domains, which are independent of each other. Under the actionof su?ciently strong electric ?elds or mechanical stress, the orienta…     

An investigation of the combined size and rate effects on the mechanical responses of FCC metals

The molecular dynamics (MD) simulations are performed with single-crystal copper blocks under simple shear to investigate the size and strain rate effects on the mechanical responses of face-centered cubic (fcc) metals. It is shown that the yield stress decreases with the specimen size and increases with the strain rate. Based on the theory of dislocation nucleation, a modified power law is proposed to predict the scaled behavior of fcc metals, which agrees well with the numerical and experimental data ranging from nanoscale to macro-scale. In the MD simulations with different strain rates, a critical strain rate exists for each single-crystal copper block of given size, below which the yield stress is nearly insensitive to the strain rate. A hyper-surface is therefore formulated to describe the combined size and strain rate effects on the plastic yield strength of fcc metals. The preliminary results presented in this paper demonstrate the potential of the proposed simple procedure for engineering design at various spatial and temporal scales.     

MODELING AND DYNAMICS ANALYSIS OF SHELLS OF REVOLUTION BY PARTIALLY ACTIVE CONSTRAINED LAYER DAMPING TREATMENT

A new model for a smart shell of revolution treated with active constrained layer damping （ACLD） is developed, and the damping effects of the ACLD treatment are discussed. The motion and electric analytical formulation of the piezoelectric constrained layer are presented first. Based on the authors~ recent research on shells of revolution treated with passive constrained layer damping （PCLD）, the integrated first-order differential matrix equation of a shell of revolution partially treated with ring ACLD blocks is derived in the frequency domain. By virtue of the extended homogeneous capacity precision integration technology, a stable and simple numerical method is further proposed to solve the above equation. Then, the vibration responses of an ACLD shell of revolution are measured by using the present model and method. The results show that the control performance of the ACLD treatment is complicated and frequency-dependent. In a certain frequency range, the ACLD treatment can achieve better damping characteristics compared with the conventional PCLD treatment.     

Heat transfer and thermal management of electric vehicle batteries with phase change materials

This paper examines a passive thermal management system for electric vehicle batteries, consisting of encapsulated phase change material (PCM) which melts during a process to absorb the heat generated by a battery. A new configuration for the thermal management system, using double series PCM shells, is analyzed with finite volume simulations. A combination of computational fluid dynamics (CFD) and second law analysis is used to evaluate and compare the new system against the single PCM shells. Using a finite volume method, heat transfer in the battery pack is examined and the results are used to analyse the exergy losses. The simulations provide design guidelines for the thermal management system to minimize the size and cost of the system. The thermal conductivity and melting temperature are studied as two important parameters in the configuration of the shells. Heat transfer from the surroundings to the PCM shell in a non-insulated case is found to be infeasible. For a single PCM system, the exergy efficiency is below 50%. For the second case for other combinations, the exergy efficiencies ranged from 30–40%. The second shell content did not have significant influence on the exergy efficiencies. The double PCM shell system showed higher exergy efficiencies than the single PCM shell system (except a case for type PCM-1). With respect to the reference environment, it is found that in all cases the exergy efficiencies decreased, when the dead-state temperatures rises, and the destroyed exergy content increases gradually. For the double shell systems for all dead-state temperatures, the efficiencies were very similar. Except for a dead-state temperature of 302 K, with the other temperatures, the exergy efficiencies for different combinations are well over 50%. The range of exergy efficiencies vary widely between 15 and 85% for a single shell system, and between 30–80% for double shell systems.