Stress intensity factors calculation in anti-plane fracture problem by orthogonal integral extraction method based on femol

For an anti-plane problem, the differential operator is self-adjoint and the corresponding eigenfunctions belong to the Hilbert space. The orthogonal property between eigenfunctions (or between the derivatives of eigenfunctions) of anti-plane problem is exploited. We developed for the first time two sets of radius-independent orthogonal integrals for extraction of stress intensity factors (SIFs), so any order SIF can be extracted based on a certain known solution of displacement (an analytic result or a numerical result). Many numerical examples based on the finite element method of lines (FEMOL) show that the present method is very powerful and efficient.

Extremum principle of crystal plasticity

An analytical method for determining slip shear rate under prescribed stress rate or prescribed strain rate has been presented on the basis of the incremental theory of crystal plasticity. The problem has been reduced to a quadric convex programming.In order to analyse the plastic response of crystals subjected to external load, two new extremum principles are proposed. They are equivalent to the boundary-value problem of crystal plasticity. By the new extremum principles, the slip shear rates are independent function which can be obtained from the variational equation.

Magnetostatic relations including fluctuation fields - the local expansion

Using the approach of local expansion, we analyze the magnetostatic relations in the case of conventional turbulence. The turbulent relations are obtained consisten tly for themomentum equation and induction equation of both the average and fluctuation relations.In comparison with the magnetostatic relations as discussed usually, turbulent fluctuationfields produce forces, one of which 1/(4π)(α1×B0)×B0 may have parallel and perpendicular components in the direction of magnetic field, the other of which 1/(4π)K×B0 is introduced by the boundary value of turbulence and is perpendicular to the magnetic field. In the case of 2-dimensional configuration of magnetic field, the basic equation will be reduced into a second-order elliptic equation, which includes some linear and nonlinear terms introduced by turbulent fluctuation fields. Turbulent fields may change the configuration of magnetic field and even shear it non-uniformly. The study on the influence of turbulent fields is significant since they are observed in many astrophysical environments.

Low temperature annealing effects on the structure and optical properties of ZnO films grown by pulsed laser deposition

ZnO thin films were deposited on glass substrates at room temperature (RT) similar to 500 degrees C by pulsed laser deposition (PLD) technique and then were annealed at 150-450 degrees C in air. The effects of annealing temperature on the microstructure and optical properties of the thin films deposited at each substrate temperature were investigated by XRD, SEM, transmittance spectra, and photoluminescence (PL). The results showed that the c-axis orientation of ZnO thin films was not destroyed by annealing treatments: the grain size increased and stress relaxed for the films deposited at 200-500 degrees C, and thin films densified for the films deposited at RT with increasing annealing temperature. The transmittance spectra indicated that E-g of thin films showed a decreased trend with annealing temperature. From the PL measurements, there was a general trend, that is UV emission enhanced with lower annealing temperature and disappeared at higher annealing temperature for the films deposited at 200-500 degrees C; no UV emission was observed for the films deposited at RT regardless of annealing treatment. Improvement of grain size and stoichiometric ratio with annealing temperature can be attributed to the enhancement of UV emission, but the adsorbed oxygen species on the surface and grain boundary of films are thought to contribute the annihilation of UV emission. It seems that annealing at lower temperature in air is an effective method to improve the UV emission for thin films deposited on glass substrate at substrate temperature above RT.

Excitonic optical absorption in semiconductors under intense terahertz radiation

The excitonic optical absorption of GaAs bulk semiconductors under intense terahertz (THz) radiation is investigated numerically. The method of solving initial-value problems, combined with the perfect matched layer technique, is used to calculate the optical susceptibility. In the presence of a driving THz field, in addition to the usual exciton peaks, 2p replica of the dark 2p exciton and even-THz-photon-sidebands of the main exciton resonance emerge in the continuum above the band edge and below the main exciton resonance. Moreover, to understand the shift of the position of the main exciton peak under intense THz radiation, it is necessary to take into consideration both the dynamical Franz-Keldysh effect and ac Stark effect simultaneously. For moderate frequency fields, the main exciton peak decreases and broadens due to the field-induced ionization of the excitons with THz field increasing. However, for high frequency THz fields, the characteristics of the exciton recur even under very strong THz fields, which accords with the recent experimental results qualitatively.

Analysis for temperature and pressure fields in process of hydrate dissociation by depressurization

An improved axisymmetric mathematic modeling is proposed for the process of hydrate dissociation by depressurization around vertical well. To reckon in the effect of latent heat of gas hydrate at the decomposition front, the energy balance equation is employed. The semi-analytic solutions for temperature and pressure fields are obtained by using Boltzmann-transformation. The location of decomposition front is determined by solving initial value problem for system of ordinary differential equations. The distributions of pressure and temperature along horizontal radiate in the reservoir are calculated. The numeric results indicate that the moving speed of decomposition front is sensitively dependent on the well pressure and the sediment permeability. Copyright (C) 2010 John Wiley & Sons, Ltd.

Dielectric responses of anisotropic graded granular composites having arbitrary inclusion shapes

The transformation field method (TFM) originated from Eshelby's transformation field theory is developed to estimate the effective permittivity of an anisotropic graded granular composite having inclusions of arbitrary shape and arbitrary anisotropic grading profile. The complicated boundary-value problem of the anisotropic graded composite is solved by introducing an appropriate transformation field within the whole composite region. As an example, the effective dielectric response for an anisotropic graded composite with inclusions having arbitrary geometrical shape and arbitrary grading profile is formulated. The validity of TFM is tested by comparing our results with the exact solution of an isotropic graded composite having inclusions with a power-law dielectric grading profile and good agreement is achieved in the dilute limit. Furthermore, it is found that the inclusion shape and the parameters of the grading profile can have profound effect on the effective permittivity at high concentrations of the inclusions. It is pointed out that TFM used in this paper can be further extended to investigate the effective elastic, thermal, and electroelastic properties of anisotropic graded granular composite materials.

Transformation field method for calculating the effective properties of isotropic graded composites having arbitrary shapes

A method of transformation field is developed to estimate the effective properties of graded composites whose inclusions have arbitrary shapes and gradient profiles by means of a periodic cell model. The boundary-value problem of graded composites having arbitrary inclusion shapes is solved by introducing the transformation field into the inclusion region. As an example, the effective dielectric response of isotropic graded composites having arbitrary shapes and gradient profiles is handled by the transformation field method (TFM). Moreover, TFM results are validated by the exact solutions of isotropic graded spherical inclusions having a power-law profile and good agreement is obtained in the dilute limit. Furthermore, it is found that the inclusion shapes and the parameters of the gradient profiles can have profound effect on the effective properties of composite systems at high concentration of inclusions.

化学冷铁——焓变涂料的研究

焓变涂料是以在一定的温度区间内可发生吸热化学反应的物质为基料配制的,将其涂在砂型(芯)表面,可代替常规冷铁.它对解决难于放置冷铁的串皮芯部位的激冷问题,效果明显.焓变涂料激冷作用的大小,取决于基料物质比焓变的大小以及涂层厚度等因素。

扩展CSS仿真程序及其应用

本文介绍了一个连续-离散复合仿真程序。该程序的连续部分基本取自美国的 CSS 程序。经作者扩充后成为复合仿真程序。扩充的模块包括采样系统仿真模块、两点边值问题求解以及参量最优化等。作者提出的单纯形搜索法改进策略可以显著加快参量最优化过程的收敛速度。

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.

Application of enthalpy method for moving boundary problem in laser surface melting

A numerical analysis was carried out to study the moving boundary problem in the physical process of pulsed Nd-YAG laser surface melting prior to vaporization. The enthalpy method was applied to solve this two-phase axisymmetrical melting problem Computational results of temperature fields were obtained, which provide useful information to practical laser treatment processing. The validity of enthalpy method in solving such problems is presented.

Boundary knot method for 2D and 3D Helmholtz and convection-diffusion problems under complicated geometry

The boundary knot method (BKM) of very recent origin is an inherently meshless, integration-free, boundary-type, radial basis function collocation technique for the numerical discretization of general partial differential equation systems. Unlike the method of fundamental solutions, the use of non-singular general solution in the BKM avoids the unnecessary requirement of constructing a controversial artificial boundary outside the physical domain. The purpose of this paper is to extend the BKM to solve 2D Helmholtz and convection-diffusion problems under rather complicated irregular geometry. The method is also first applied to 3D problems. Numerical experiments validate that the BKM can produce highly accurate solutions using a relatively small number of knots. For inhomogeneous cases, some inner knots are found necessary to guarantee accuracy and stability. The stability and convergence of the BKM are numerically illustrated and the completeness issue is also discussed.

Improving Nucleation in the Fabrication of High-Quality 3D Macro-Porous Copper Film through The Surface-Modification of A Polystyrene Colloid-Assembled Template

A new approach is developed to the fabrication of high-quality three-dimensional macro-porous copper films. A highly-ordered macroporous copper film is successfully produced on a polystyrene sphere (PS) template that has been modified by sodium dodecyl sulfate (SDS). It is shown that this procedure can change a hydrophobic surface of PS template into a hydrophilic surface. The present study is devoted to the influence of the electrolyte solution transport on the nucleation process. It is demonstrated that the permeability of the electrolyte solution in the nanochannels of the PS template plays an important role in the chemical electrodeposition of high-quality macroporous copper film. The permeability is drastically enhanced in our experiment through the surface modi. cation of the PS templates. The method could be used to homogeneously produce a large number of nucleations on a substrate, which is a key factor for the fabrication of the high-quality macroporous copper film.

Molecular dynamics studies on the dislocation gliding near a tilt boundary

The gliding behavior of edge dislocation near a grain boundary(QB) in copper under pure shear stresses is simulated by using molecular dynamics(MD) method. Many-body potential incorporating the embedded atom method (EAM) is used. The critical shear stresses for a single disocation to pass across GB surface are obtained at values of sigma(c)=23MPa similar to 68 MPa and 137 MPa similar to 274 MPa for Sigma=165 small angle tilt GB at 300 K and 20 K, respectively. The first result agrees with the experimental yield stress sigma(y)(=42 MPa) quite well. It suggests that there might be one of the reasons of initial plastic yielding caused by single dislocation gliding across GB. In addition, there might be possibility to obtain yield strength from microscopic analysis. Moreover, the experimental value of sigma(y) at low temperature is generally higher than that at room temperature. So, these results are in conformity qualitatively with experimental fact. On the other hand, the Sigma=25 GB is too strong an obstacle to the dislocation. In this case, a dislocation is able to pass across GB under relatively low stress only when it is driven by other dislocations. This is taken to mean that dislocation pile-up must be built up in front of this kind of GB, if this GB may take effect on the process of plastic deformation.