High-order asymptotic field of tensile plane-strain nonlinear crack problems

In this paper, the governing equations and the analytical method of the secondorder asymptotic field for the plane-straln crack problems of mode I have been presented. The numerical calculation has been carried out. The amplitude coefficients k2 of the second term of the asymptotic field have been determined after comparing the present results with some fine results of the finite element calculation. The variation of coefficients k2 with changes of specimen geometry and developments of plastic zone have been analysed. It is shown that the second term of the asymptotic field has significant influence on the near-crack-tip field. Therefore, we may reasonably argue that both the J-integral and the coefficient k2 can beeome two characterizing parameters for crack initiation.

An asymptotic solution of hydrodynamic equations at the mid-plane of a plume in the earths mantle

From the partial differential equations of hydrodynamics governing the movements in the Earth's mantle of a Newtonian fluid with a pressure- and temperature-dependent viscosity, considering the bilateral symmetry of velocity and temperature distributions at the mid-plane of the plume, an analytical solution of the governing equations near the mid-plane of the plume was found by the method of asymptotic analysis. The vertical distribution of the upward velocity, viscosity and temperature at the mid-plane, and the temperature excess at the centre of the plume above the ambient mantle temperature were then calculated for two sets of Newtonian rheological parameters. The results obtained show that the temperature at the mid-plane and the temperature excess are nearly independent of the rheological parameters. The upward velocity at the mid-plane, however, is strongly dependent on the rheological parameters.

Variational approach to the closure problem of turbulence theory

A new method is proposed to solve the closure problem of turbulence theory and to drive the Kolmogorov law in an Eulerian framework. Instead of using complex Fourier components of velocity field as modal parameters, a complete set of independent real parameters and dynamic equations are worked out to describe the dynamic states of a turbulence. Classical statistical mechanics is used to study the statistical behavior of the turbulence. An approximate stationary solution of the Liouville equation is obtained by a perturbation method based on a Langevin-Fokker-Planck (LFP) model. The dynamic damping coefficient eta of the LFP model is treated as an optimum control parameter to minimize the error of the perturbation solution; this leads to a convergent integral equation for eta to replace the divergent response equation of Kraichnan's direct-interaction (DI) approximation, thereby solving the closure problem without appealing to a Lagrangian formulation. The Kolmogorov constant Ko is evaluated numerically, obtaining Ko = 1.2, which is compatible with the experimental data given by Gibson and Schwartz, (1963).

Perturbational finite volume scheme for the one-dimensional Navier-Stokes equations

Starting from the second-order finite volume scheme,though numerical value perturbation of the cell facial fluxes, the perturbational finite volume (PFV) scheme of the Navier-Stokes (NS) equations for compressible flow is developed in this paper. The central PFV scheme is used to compute the one-dimensional NS equations with shock wave.Numerical results show that the PFV scheme can obtain essentially non-oscillatory solution.

Basic equations for ferroelectrics

A set of new formula of energy functions for ferroelectrics was proposed, and then the new basic equations were derived in this paper. The finite element formulation based on the new basic equations was improved to avoid the equivalent nodal load produced by remnant polarization. With regard to the fundamentals of mathematics and physics, the new energy functions and basic equations are reasonable for the material element of ferroelectrics in finite element analysis.

Energy Functions And Basic Equations For Ferroelectrics

Energy functions (or characteristic functions) and basic equations for ferroelectrics in use today are given by those for ordinary dielectrics in the physical and mechanical communications. Based on these basic equations and energy functions, the finite element computation of the nonlinear behavior of the ferroelectrics has been carried out by several research groups. However, it is difficult to process the finite element computation further after domain switching, and the computation results are remarkably deviating from the experimental results. For the crack problem, the iterative solution of the finite element calculation could not converge and the solutions for fields near the crack tip oscillate. In order to finish the calculation smoothly, the finite element formulation should be modified to neglect the equivalent nodal load produced by spontaneous polarization gradient. Meanwhile, certain energy functions for ferroelectrics in use today are not compatible with the constitutive equations of ferroelectrics and need to be modified. This paper proposes a set of new formulae of the energy functions for ferroelectrics. With regard to the new formulae of the energy functions, the new basic equations for ferroelectrics are derived and can reasonably explain the question in the current finite element analysis for ferroelectrics.

A block SPP algorithm for multidimensional tridiagonal equations with optimal message vector length

A parallel strategy for solving multidimensional tridiagonal equations is investigated in this paper. We present in detail an improved version of single parallel partition (SPP) algorithm in conjunction with message vectorization, which aggregates several communication messages into one to reduce the communication cost. We show the resulting block SPP can achieve good speedup for a wide range of message vector length (MVL), especially when the number of grid points in the divided direction is large. Instead of only using the largest possible MVL, we adopt numerical tests and modeling analysis to determine an optimal MVL so that significant improvement in speedup can be obtained.

Work and energy equations and the principle of generalized effective stress for unsaturated soils

The effective stress principle has been efficiently applied to saturated soils in the soil mechanics and geotechnical engineering practice; however, its applicability to unsaturated soils is still under debate. The appropriate selection of stress state variables is essential for the construction of constitutive models for unsaturated soils. Owing to the complexity of unsaturated soils, it is difficult to determine the deformation and strength behaviors of unsaturated soils uniquely with the previous single-effective-stress variable theory and two-effective-stress-variable theory in all the situations. In this paper, based on the porous media theory, the specific expression of work is proposed, and the effective stress of unsaturated soils conjugated with the displacement of the soil skeleton is further derived. In the derived work and energy balance equations, the energy dissipation in unsaturated soils is taken into account. According to the derived work and energy balance equations, all of the three generalized stresses and the conjugated strains have effects on the deformation of unsaturated soils. For considering these effects, a principle of generalized effective stress to describe the behaviors of unsaturated soils is proposed. The proposed principle of generalized effective stress may reduce to the previous effective stress theory of single-stress variable or the two-stress variables under certain conditions. This principle provides a helpful reference for the development of constitutive models for unsaturated soils.

Diffraction integral formulas of the pulsed wave field in the temporal domain

To resolve the diffraction problems of the pulsed wave field directly in the temporal domain, we extend the Rayleigh diffraction integrals to the temporal domain and then discuss the approximation condition of this diffraction formula. (C) 1997 Optical Society of America.

平面交叉型微透镜阵列的制作及成像特性研究

介绍了采用光刻离子交换工艺制作平面交叉型微透镜阵列的方法。利用积分形式的光线方程式讨论了平面交叉型微透镜的近轴光学特性，研究了微透镜的光线轨迹方程式和一些重要的近轴成像特性，利用ABCD定理得到了平面交叉型微透镜像距、焦距、像高、横向放大率和主平面位置的数学表达式，焦距的理论计算结果和实验数据吻合得很好。

Vectorial solution of Kukhtarev equations for doubly doped crystals and optimal choice of recording directions in nonvolatile holographic storage

Vectorial Kukhtarev equations modified for the nonvolatile holographic recording in doubly doped crystals are analyzed, in which the bulk photovoltaic effect and the external electrical field are both considered. On the basis of small modulation approximation, both the analytic solution to the space-charge field with time in the recording phase and in the readout phase are deduced. The analytic solutions can be easily simplified to adapt the one-center model, and they have the same analytic expressions given those when the grating vector is along the optical axis. Based on the vectorial analyses of the band transport model an optimal recording direction is given to maximize the refractive index change in doubly doped LiNbO3:Fe: Mn crystals. (c) 2007 Optical Society of America.

Smoothing the side lobes of a focused pattern by spectral dispersion in the broadband laser

Starting from the Huygens-Fresnel diffraction integral, the propagation equations of a broadband laser passing through a dispersive lens and a dispersive wedge are derived. Smoothing effect on the side lobes of the focused pattern is achieved as the broadband laser passes through the lens because of the spectral dispersion of the lens. By inserting a dispersive wedge behind the lens, better smoothing effect is realized because a relative position shift between focused patterns of different frequency components is generated due to the spectral dispersion of the wedge. Smoothing effect on the side lobe is obtained even with small bandwidth of the broadband laser as the wedge is used. (c) 2006 Elsevier GmbH. All rights reserved.

Choice of calibration equations in calibrating microwave test fixtures

The error theory of linear equation system has been applied to the calibration procedure of microwave network analyser in this article. A new explanation for the choice of the linear calibration equations is proposed and a general principle for choosing calibration equations is presented. The method can also be used to predict the occurrence of the problem of frequency limitation at some periodic frequencies. This principle is employed to the thru-short-delay (TSD) method and the solution using the chosen equations gives the most accurate results. A good agreement between the theory and the experiment has been obtained.

A 16-term error model based on linear equations of voltage and current variables

Formulation of a 16-term error model, based on the four-port ABCD-matrix and voltage and current variables, is outlined. Matrices A, B, C, and D are each 2 x 2 submatrices of the complete 4 x 4 error matrix. The corresponding equations are linear in terms of the error parameters, which simplifies the calibration process. The parallelism with the network analyzer calibration procedures and the requirement of five two-port calibration measurements are stressed. Principles for robust choice of equations are presented. While the formulation is suitable for any network analyzer measurement, it is expected to be a useful alternative for the nonlinear y-parameter approach used in intrinsic semiconductor electrical and noise parameter measurements and parasitics' deembedding.