Damage analysis of tetragonal perovskite structure ceramics implicated by asymptotic field solutions and boundary conditions

Cracking of ceramics with tetragonal perovskite grain structure is known to appear at different sites and scale level. The multiscale character of damage depends on the combined effects of electromechanical coupling, prevailing physical parameters and boundary conditions. These detail features are exhibited by application of the energy density criterion with judicious use of the mode I asymptotic and full field solution in the range of r/a = 10(-4) to 10(-2) where r and a are, respectively, the distance to the crack tip and half crack length. Very close to the stationary crack tip, bifurcation is predicted resembling the dislocation emission behavior invoked in the molecular dynamics model. At the macroscopic scale, crack growth is predicted to occur straight ahead with two yield zones to the sides. A multiscale feature of crack tip damage is provided for the first time. Numerical values of the relative distances and bifurcation angles are reported for the PZT-4 ceramic subjected to different electric field to applied stress ratio and boundary conditions that consist of the specification of electric field/mechanical stress, electric displacement/mechanical strain, and mixed conditions. To be emphasized is that the multiscale character of damage in piezoceramics does not appear in general. It occurs only for specific combinations of the external and internal field parameters, elastic/piezoelectric/dielectric constants and specified boundary conditions. (C) 2002 Published by Elsevier Science Ltd.

Boundary conditions at the solid-liquid surface over the multiscale channel size from nanometer to micron

The boundary condition at the solid surface is one of the important problems for the microfluidics. In this paper we study the effects of the channel sizes on the boundary conditions (BC), using the hybrid computation scheme adjoining the molecular dynamics (MD) simulations and the continuum fluid mechanics. We could reproduce the three types of boundary conditions (slip, no-slip and locking) over the multiscale channel sizes. The slip lengths are found to be mainly dependent on the interfacial parameters with the fixed apparent shear rate. The channel size has little effects on the slip lengths if the size is above a critical value within a couple of tens of molecular diameters. We explore the liquid particle distributions nearest the solid walls and found that the slip boundary condition always corresponds to the uniform liquid particle distributions parallel to the solid walls, while the no-slip or locking boundary conditions correspond to the ordered liquid structures close to the solid walls. The slip, no-slip and locking interfacial parameters yield the positive, zero and negative slip lengths respectively. The three types of boundary conditions existing in "microscale" still occur in "macroscale". However, the slip lengths weakly dependent on the channel sizes yield the real shear rates and the slip velocity relative to the solid wall traveling speed approaching those with the no-slip boundary condition when the channel size is larger than thousands of liquid molecular diameters for all of the three types of interfacial parameters, leading to the quasi-no-slip boundary conditions.

A New Criterion Number for the Boundary Conditions at the Solid/Liquid Interface in Nanoscale

A simple, but important three-atom model was proposed at the solid/liquid interface, leading to a new criterion number, lambda, governing the boundary conditions (BCs) in nanoscale. The solid wall is considered as the face-centered-cubic (fcc) structure. The fluid is the liquid argon with the well-known LJ potential. Based on the concept, the two micro-systems have the same BCs if they have The same criterion number. The degree of the locking BCs is enhanced when lambda equals to 0.757. Such critical criterion number results in the substantial epitaxial ordering and one, two, or even three liquid layers are locked by the solid wall, depending on the coupling energy scale ratio of the solid and liquid atoms. With deviation from the critical criterion number, the flow approaches the slip BCs and there are little ordering structures within the liquid. Always at the same criterion number, the degree of the slip is decreased or the locking is enhanced with increasing the coupling energy scale ratio of the solid and liquid atoms. The above analysis is well confirmed by the molecular dynamics (MD) simulation. The slip length is well correlated in terms of the new criterion number. The future work is suggested to extend the present theory for other microstructures of the solid wall atoms and quasi-LJ potentials.

A Study of Composite Beam with Shape Memory Alloy Arbitrarily Embedded under Thermal and Mechanical Loadings

The constitutive relations and kinematic assumptions on the composite beam with shape memory alloy (SMA) arbitrarily embedded are discussed and the results related to the different kinematic assumptions are compared. As the approach of mechanics of materials is to study the composite beam with the SMA layer embedded, the kinematic assumption is vital. In this paper, we systematically study the kinematic assumptions influence on the composite beam deflection and vibration characteristics. Based on the different kinematic assumptions, the equations of equilibrium/motion are different. Here three widely used kinematic assumptions are presented and the equations of equilibrium/motion are derived accordingly. As the three kinematic assumptions change from the simple to the complex one, the governing equations evolve from the linear to the nonlinear ones. For the nonlinear equations of equilibrium, the numerical solution is obtained by using Galerkin discretization method and Newton-Rhapson iteration method. The analysis on the numerical difficulty of using Galerkin method on the post-buckling analysis is presented. For the post-buckling analysis, finite element method is applied to avoid the difficulty due to the singularity occurred in Galerkin method. The natural frequencies of the composite beam with the nonlinear governing equation, which are obtained by directly linearizing the equations and locally linearizing the equations around each equilibrium, are compared. The influences of the SMA layer thickness and the shift from neutral axis on the deflection, buckling and post-buckling are also investigated. This paper presents a very general way to treat thermo-mechanical properties of the composite beam with SMA arbitrarily embedded. The governing equations for each kinematic assumption consist of a third order and a fourth order differential equation with a total of seven boundary conditions. Some previous studies on the SMA layer either ignore the thermal constraint effect or implicitly assume that the SMA is symmetrically embedded. The composite beam with the SMA layer asymmetrically embedded is studied here, in which symmetric embedding is a special case. Based on the different kinematic assumptions, the results are different depending on the deflection magnitude because of the nonlinear hardening effect due to the (large) deflection. And this difference is systematically compared for both the deflection and the natural frequencies. For simple kinematic assumption, the governing equations are linear and analytical solution is available. But as the deflection increases to the large magnitude, the simple kinematic assumption does not really reflect the structural deflection and the complex one must be used. During the systematic comparison of computational results due to the different kinematic assumptions, the application range of the simple kinematic assumption is also evaluated. Besides the equilibrium study of the composite laminate with SMA embedded, the buckling, post-buckling, free and forced vibrations of the composite beam with the different configurations are also studied and compared.

A field model interpretation of crack initiation and growth behavior in ferroelectric ceramics: change of poling direction and boundary condition

Strain energy density expressions are obtained from a field model that can qualitatively exhibit how the electrical and mechanical disturbances would affect the crack growth behavior in ferroelectric ceramics. Simplification is achieved by considering only three material constants to account for elastic, piezoelectric and dielectric effects. Cross interaction of electric field (or displacement) with mechanical stress (or strain) is identified with the piezoelectric effect; it occurs only when the pole is aligned normal to the crack. Switching of the pole axis by 90degrees and 180degrees is examined for possible connection with domain switching. Opposing crack growth behavior can be obtained when the specification of mechanical stress sigma(infinity) and electric field E-infinity or (sigma(infinity), E-infinity) is replaced by strain e and electric displacement D-infinity or (epsilon(infinity), D-infinity). Mixed conditions (sigma(infinity),D-infinity) and (epsilon(infinity),E-infinity) are also considered. In general, crack growth is found to be larger when compared to that without the application of electric disturbances. This includes both the electric field and displacement. For the eight possible boundary conditions, crack growth retardation is identified only with (E-y(infinity),sigma(y)(infinity)) for negative E-y(infinity) and (D-y(infinity), epsilon(y)(infinity)) for positive D-y(infinity) while the mechanical conditions sigma(y)(infinity) or epsilon(y)infinity are not changed. Suitable combinations of the elastic, piezoelectric and dielectric material constants could also be made to suppress crack growth. (C) 2002 Published by Elsevier Science Ltd.

Compact difference approximation with consistent boundary condition

For simulating multi-scale complex flow fields it should be noted that all the physical quantities we are interested in must be simulated well. With limitation of the computer resources it is preferred to use high order accurate difference schemes. Because of their high accuracy and small stencil of grid points computational fluid dynamics (CFD) workers pay more attention to compact schemes recently. For simulating the complex flow fields the treatment of boundary conditions at the far field boundary points and near far field boundary points is very important. According to authors' experience and published results some aspects of boundary condition treatment for far field boundary are presented, and the emphasis is on treatment of boundary conditions for the upwind compact schemes. The consistent treatment of boundary conditions at the near boundary points is also discussed. At the end of the paper are given some numerical examples. The computed results with presented method are satisfactory.

Analysis of flexural vibration of viscoelastically damped sandwich plates

The simplified governing equations and corresponding boundary conditions of flexural vibration of viscoelastically damped unsymmetrical sandwich plates are given. The asymptotic solution of the equations is then discussed. If only the first terms of the asymptotic solution of all variables are taken as an approximate solution, the result is identical with that obtained from the Modal Strain Energy (MSE) Method. As more terms of the asymptotic solution are taken, the successive calculations show improved accuracy. With the natural frequencies and the modal loss factors of a damped sandwich plate known, one can calculate the response of the plate to various loads providing a reliable basis for engineering design.

A smoothing technique for discrete delta functions with application to immersed boundary method in moving boundary simulations.

The effects of complex boundary conditions on flows are represented by a volume force in the immersed boundary methods. The problem with this representation is that the volume force exhibits non-physical oscillations in moving boundary simulations. A smoothing technique for discrete delta functions has been developed in this paper to suppress the non-physical oscillations in the volume forces. We have found that the non-physical oscillations are mainly due to the fact that the derivatives of the regular discrete delta functions do not satisfy certain moment conditions. It has been shown that the smoothed discrete delta functions constructed in this paper have one-order higher derivative than the regular ones. Moreover, not only the smoothed discrete delta functions satisfy the first two discrete moment conditions, but also their derivatives satisfy one-order higher moment condition than the regular ones. The smoothed discrete delta functions are tested by three test cases: a one-dimensional heat equation with a moving singular force, a two-dimensional flow past an oscillating cylinder, and the vortex-induced vibration of a cylinder. The numerical examples in these cases demonstrate that the smoothed discrete delta functions can effectively suppress the non-physical oscillations in the volume forces and improve the accuracy of the immersed boundary method with direct forcing in moving boundary simulations.

One-dimensional quantum waveguide theory of Rashba electrons

The ballistic spin transport in one-dimensional waveguides with the Rashba effect is studied. Due to the Rashba effect, there are two electron states with different wave vectors for the same energy. The wave functions of two Rashba electron states are derived, and it is found that their phase depend on the direction of the circuit and the spin directions of two states are perpendicular to the circuit, with the +pi/2 and -pi/2 angles, respectively. The boundary conditions of the wave functions and their derivatives at the intersection of circuits are given, which can be used to investigate the waveguide transport properties of Rashba spin electron in circuits of any shape and structure. The eigenstates of the closed circular and square loops are studied by using the transfer matrix method. The transfer matrix M(E) of a circular arc is obtained by dividing the circular arc into N segments and multiplying the transfer matrix of each straight segment. The energies of eigenstates in the closed loop are obtained by solving the equation det[M(E)-I]=0. For the circular ring, the eigenenergies obtained with this method are in agreement with those obtained by solving the Schrodinger equation. For the square loop, the analytic formula of the eigenenergies is obtained first The transport properties of the AB ring and AB square loop and double square loop are studied using the boundary conditions and the transfer matrix method In the case of no magnetic field, the zero points of the reflection coefficients are just the energies of eigenstates in closed loops. In the case of magnetic field, the transmission and reflection coefficients all oscillate with the magnetic field; the oscillating period is Phi(m)=hc/e, independent of the shape of the loop, and Phi(m) is the magnetic flux through the loop. For the double loop the oscillating period is Phi(m)=hc/2e, in agreement with the experimental result. At last, we compared our method with Koga's experiment. (C) 2009 American Institute of Physics. [doi: 10.1063/1.3253752]

Micromixing enhanced by pulsating flows

We consider the micromixing enhancement by pulsating flows. Dimensionless governing equations and boundary conditions were developed for T-type micromixers with two inlet pulsating flows. The problem involves a set of parameters. Three key dimensionless parameters are identified: the Reynolds number, the Strouhal number, and the disturbance amplitude. Suitable Strouhal number or disturbance amplitude causes symmetrical meniscus-shape mixing interfaces, separating the whole mixing channel into a set of segments. Thus uniform exit species concentration can be reached. Too large or too small Strouhal number or disturbance amplitude yields the meniscus-shape mixing interfaces deviating from the centerline of the mixing channel, deteriorating the mixing performance. The optimized disturbance amplitude is increased with increases in Strouhal numbers. Low Reynolds number needs larger disturbance amplitude.

Ground state energy of He isoelectronic sequence treated variationally via Hylleraas-like wavefunction

In this study, the energy for the ground state of helium and a few helium-like ions (Z=1-6) is computed variationally by using a Hylleraas-like wavefunction. A four-parameters wavefunction, satisfying boundary conditions for coalescence points, is combined with a Hylleraas-like basis set which explicitly incorporates r12 interelectronic distance. The main contribution of this work is the introduction of modified correlation terms leading to the definition of integral transforms which provide the calculation of expectation value of energy to be done analytically over single-particle coordinates instead of Hylleraas coordinates.

黄土区考虑滞后作用的坡地水量转化模型

建立了考虑滞后作用的坡地水量转化数学模型，并对模型不同情况下的定解条件进行了分析。利用陕西长武王东沟小流域坡地天然降雨－入渗－产流的实测资料对模型进行了系统检验，结果表明该模型具有较高精度，考虑滞后作用时的模拟值更加接近实测值。在降雨入渗过程中，滞后作用对土壤含水率、入渗速率、径流过程、产流量都有相当大的影响。

A theoretical framework of tracer methods for marine sediment dynamics

A new theoretical framework of tracer methods is proposed in the present contribution, on the basis of mass conservation. This model is applicable for both artificial and natural tracers. It can be used to calculate the spatial distribution patterns of sediment transport rate, thus providing independent information and verification for the results derived from empirical formulae. For the procedures of the calculation, first, the tracer concentration and topographic maps of two times are obtained. Then, the spatial and temporal changes in the concentration and seabed elevation are calculated, and the boundary conditions required are determined by field observations (such as flow and bedform migration measurements). Finally, based upon eqs. (1) and (13), the transport rate is calculated and expressed as a function of the position over the study area. Further, appropriate modifications to the model may allow the tracer to have different densities and grain size distributions from the bulk sediment.

Observations on the nocturnal activity and feeding behavior of Anguilla japonica glass eels under laboratory conditions

Glass eels of the temperate anguillid species, Anguilla japonica, clearly showed a nocturnal activity rhythm under laboratory conditions. Light-dark cycle was a determinant factor affecting their photonegative behavior, nocturnal locomotor activity, and feeding behavior. Under natural light conditions, glass eels remained in shelters with little daytime feeding, but came out to forage during darkness. They moved and foraged actively in the following dark, and then their activity gradually declined possibly because of food satiation. They finally buried in the sand or stayed in tubes immediately after the lights came on. Under constant light, glass eels often came out of the shelters to forage in the lights but spent little time moving outside the shelters (e.g. swimming or crawling on the sand). Glass eels took shelter to avoid light and preferred tubes to sand for shelter possibly because tubes were much easier for them to take refuge in than sand. Feeding and locomotor activities of the glass eels were nocturnal and well synchronized. They appeared to depend on olfaction rather than vision to detect and capture prey in darkness. Feeding was the driving force for glass eels to come out of sand under constant light. However, in the dark, some glass eels swam or crept actively on sand even when they were fully fed. The lunar cycles of activity rhythms of glass eels that have been observed in some estuarine areas were not detected under these laboratory conditions.

Boundary Element Method (BEM) Analysis for Galvanic Corrosion of Hot Dip Galvanized Steel Immersed in Seawater

A numerical analysis of galvanic corrosion of hot-dip galvanized steel immersed in seawater was presented. The analysis was based on the boundary element methods (BEMs) coupled with Newton-Raphson iterative technique to treat the nonlinear boundary conditions, which were determined by the experimental polarization curves. Results showed that galvanic current density concentrates on the boundary of steel substrate and zinc coating, and the sacrificial protection of zinc coating to steel substrate results in overprotection of steel cathode. Not only oxygen reduction but also hydrogen reduction could occur as cathode reactions, which probably led up to the adsorption and absorption of hydrogen atoms. Flat galvanized steel tensile sample shows a brittle behavior similar to hydrogen embrittlement according to the SSRT (show strain rate test) in seawater.