High-Resolution Finite Compact Difference Schemes for Hyperbolic Conservation Laws

A finite compact (FC) difference scheme requiring only bi-diagonal matrix inversion is proposed by using the known high-resolution flux. Introducing TVD or ENO limiters in the numerical flux, several high-resolution FC-schemes of hyperbolic conservation law are developed, including the FC-TVD, third-order FC-ENO and fifth-order FC-ENO schemes. Boundary conditions formulated need only one unknown variable for third-order FC-ENO scheme and two unknown variables for fifth-order FC-ENO scheme. Numerical test results of the proposed FC-scheme were compared with traditional TVD, ENO and WENO schemes to demonstrate its high-order accuracy and high-resolution.

Small-dynamic-angle measurement using tangent relationship and phase-modulating technique

A method using two prisms for measurement of small dynamic angles is proposed in which the measurement is based on a simple tangent equation and a phase-modulating interferometer with a laser diode to measure dynamic optical path differences with higher accuracy. Owing to the simple tangent equation, the symmetry requirement on the two prisms in the optical configuration is eliminated, and easy measurement of the separations between two parallel beams with a position-sensitive detector is achieved. Small-dynamic-angle measurements are experimentally demonstrated with high accuracy. (C) 2007 Society of Photo-Optical Instrumentation Engineers.

Small-dynamic-angle measurement using tangent relationship and phase-modulating technique

A method using two prisms for measurement of small dynamic angles is proposed in which the measurement is based on a simple tangent equation and a phase-modulating interferometer with a laser diode to measure dynamic optical path differences with higher accuracy. Owing to the simple tangent equation, the symmetry requirement on the two prisms in the optical configuration is eliminated, and easy measurement of the separations between two parallel beams with a position-sensitive detector is achieved. Small-dynamic-angle measurements are experimentally demonstrated with high accuracy. (C) 2007 Society of Photo-Optical Instrumentation Engineers.

A new hardening law for strain gradient plasticity

A new hardening law of the strain gradient theory is proposed in this paper, which retains the essential structure of the incremental version of conventional J(2) deformation theory and obeys thermodynamic restrictions. The key feature of the new proposal is that the term of strain gradient plasticity is represented as an internal variable to increase the tangent modulus. This feature which is in contrast to several proposed theories, allows the problem of incremental equilibrium equations to be stated without higher-order stress, higher-order strain rates or extra boundary conditions. The general idea is presented and compared with the theory given by Fleck and Hutchinson (Adv. in Appl. Mech. (1997) 295). The new hardening law is demonstrated by two experimental tests i.e. thin wire torsion and ultra-thin beam bending tests. The present theoretical results agree well with the experiment results.

双曲守恒型方程的二阶摄动有限差分格式

对双曲守恒型方程，将其一阶迎风格式空间差商的常系数摄动展开为时间步长和空间步长的幂级数，通过确定幂级数系数而获得二阶精度的摄动有限差分（PFD）格式。进而从双曲守恒型方程的通量分裂型一阶迎风格式出发，通过娄似的摄动展开方法，获得空间精度为二阶的通量分裂形式的摄动有限差分（FPFD）格式。这两类格式保留了一阶守恒迎风格式的简洁结构形式，使用三节点即可达到二阶精度，又避免了三点二阶格式的非物理数值振荡。并将这两类格式推广应用到双曲守恒型方程组，最后通过模型方程和一维激波管流动的数值算例验证了格式的高精度、高分辨率性质。

A Study of Size-Dependent Microindentation

We recently proposed a strain gradient theory to account for the size dependence of plastic deformation at micron and submicron length scales. The strain gradient theory includes the effects of both rotation gradient and stretch gradient such that the rotation gradient influences the material character through the interaction between the Cauchy stresses and the couple stresses; the stretch gradient measures explicitly enter the constitutive relations through the instantaneous tangent modulus. Indentation tests at scales on the order of one micron have shown that measured hardness increases significantly with decreasing indent size. In the present paper, the strain gradient theory is used to model materials undergoing small-scale indentations. A strong effect of including strain gradients in the constitutive description is found with hardness increasing by a factor of two or more over the relevant range behavior. Comparisons with the experimental data for polycrystalline copper and single crystal copper indeed show an approximately linear dependence of the square of the hardness, H 2, on the inverse of the indentation depth, 1/h, I.e., H-2 proportional to 1/h, which provides an important self-consistent check of the strain gradient theory proposed by the authors earlier.

Finite Element Solutions for Plane Strain Mode I Crack with Strain Grading Effect

In this paper, a new phenomenological theory with strain gradient effects is proposed to account for the size dependence of plastic deformation at micro- and submicro-length scales. The theory fits within the framework of general couple stress theory and three rotational degrees of freedom omega(i) are introduced in addition to the conventional three translational degrees of freedom mu(i). omega(i) is called micro-rotation and is the sum of material rotation plus the particles' relative rotation. While the new theory is used to analyze the crack tip field or the indentation problems, the stretch gradient is considered through a new hardening law. The key features of the theory are that the rotation gradient influences the material character through the interaction between the Cauchy stresses and the couple stresses; the term of stretch gradient is represented as an internal variable to increase the tangent modulus. In fact the present new strain gradient theory is the combination of the strain gradient theory proposed by Chen and Wang (Int. J. Plast., in press) and the hardening law given by Chen and Wang (Acta Mater. 48 (2000a) 3997). In this paper we focus on the finite element method to investigate material fracture for an elastic-power law hardening solid. With remotely imposed classical K fields, the full field solutions are obtained numerically. It is found that the size of the strain gradient dominance zone is characterized by the intrinsic material length l(1). Outside the strain gradient dominance zone, the computed stress field tends to be a classical plasticity field and then K field. The singularity of stresses ahead of the crack tip is higher than that of the classical field and tends to the square root singularity, which has important consequences for crack growth in materials by decohesion at the atomic scale. (C) 2002 Elsevier Science Ltd. All rights reserved.

Crack Tip Field and J-Integral with Strain Gradient Effect

The mode I plane strain crack tip field with strain gradient effects is presented in this paper based on a simplified strain gradient theory within the framework proposed by Acharya and Bassani. The theory retains the essential structure of the incremental version of the conventional J_2 deformation theory No higher-order stress is introduced and no extra boundary value conditions beyond the conventional ones are required. The strain gradient effects are considered in the constitutive relation only through the instantaneous tangent modulus. The strain gradient measures are included into the tangent modulus as internal parameters. Therefore the boundary value problem is the same as that in the conventional theory Two typical crack Problems are studied: (a) the crack tip field under the small scale yielding condition induced by a linear elastic mode-I K-field and (b) the complete field for a compact tension specimen. The calculated results clearly show that the stress level near the crack tip with strain gradient effects is considerable higher than that in the classical theory The singularity of the strain field near the crack tip is nearly equal to the square-root singularity and the singularity of the stress field is slightly greater than it. Consequently, the J-integral is no longer path independent and increases monotonically as the radius of the calculated circular contour decreases.

An improvement and proof of OGY method

OGY method is the most important method of controlling chaos. It stabilizes a hyperbolic periodic orbit by making small perturbations for a system parameter. This paper improves the method of choosing parameter, and gives a mathematics proof of it.

Direct Numerical Simulation Of Hypersonic Boundary-Layer Transition Over A Blunt Cone

Direct numerical simulation of transition How over a blunt cone with a freestream Mach number of 6, Reynolds number of 10,000 based on the nose radius, and a 1-deg angle of attack is performed by using a seventh-order weighted essentially nonoscillatory scheme for the convection terms of the Navier-Stokes equations, together with an eighth-order central finite difference scheme for the viscous terms. The wall blow-and-suction perturbations, including random perturbation and multifrequency perturbation, are used to trigger the transition. The maximum amplitude of the wall-normal velocity disturbance is set to 1% of the freestream velocity. The obtained transition locations on the cone surface agree well with each other far both cases. Transition onset is located at about 500 times the nose radius in the leeward section and 750 times the nose radius in the windward section. The frequency spectrum of velocity and pressure fluctuations at different streamwise locations are analyzed and compared with the linear stability theory. The second-mode disturbance wave is deemed to be the dominating disturbance because the growth rate of the second mode is much higher than the first mode. The reason why transition in the leeward section occurs earlier than that in the windward section is analyzed. It is not because of higher local growth rate of disturbance waves in the leeward section, but because the growth start location of the dominating second-mode wave in the leeward section is much earlier than that in the windward section.

A Multimoment Finite-Volume Shallow-Water Model On The Yin-Yang Overset Spherical Grid

A numerical model for shallow-water equations has been built and tested on the Yin-Yang overset spherical grid. A high-order multimoment finite-volume method is used for the spatial discretization in which two kinds of so-called moments of the physical field [i.e., the volume integrated average ( VIA) and the point value (PV)] are treated as the model variables and updated separately in time. In the present model, the PV is computed by the semi-implicit semi-Lagrangian formulation, whereas the VIA is predicted in time via a flux-based finite-volume method and is numerically conserved on each component grid. The concept of including an extra moment (i.e., the volume-integrated value) to enforce the numerical conservativeness provides a general methodology and applies to the existing semi-implicit semi-Lagrangian formulations. Based on both VIA and PV, the high-order interpolation reconstruction can only be done over a single grid cell, which then minimizes the overlapping zone between the Yin and Yang components and effectively reduces the numerical errors introduced in the interpolation required to communicate the data between the two components. The present model completely gets around the singularity and grid convergence in the polar regions of the conventional longitude-latitude grid. Being an issue demanding further investigation, the high-order interpolation across the overlapping region of the Yin-Yang grid in the current model does not rigorously guarantee the numerical conservativeness. Nevertheless, these numerical tests show that the global conservation error in the present model is negligibly small. The model has competitive accuracy and efficiency.

Analysis and application of ellipticity of stability equations on fluid mechanics

By using characteristic analysis of the linear and nonlinear parabolic stability equations (PSE), PSE of primitive disturbance variables are proved to be parabolic intotal. By using sub-characteristic analysis of PSE, the linear PSE are proved to be elliptical and hyperbolic-parabolic for velocity U, in subsonic and supersonic, respectively; the nonlinear PSE are proved to be elliptical and hyperbolic-parabolic for relocity U + u in subsonic and supersonic, respectively. The methods are gained that the remained ellipticity is removed from the PSE by characteristic and sub-characteristic theories, the results for the linear PSE are consistent with the known results, and the influence of the Mach number is also given out. At the same time, the methods of removing the remained ellipticity are further obtained from the nonlinear PSE.

Singularity of dynamic stress fields around an interface crack between viscoelastic bodies

The singular nature of the dynamic stress fields around an interface crack located between two dissimilar isotropic linearly viscoelastic bodies is studied. A harmonic load is imposed on the surfaces of the interface crack. The dynamic stress fields around the crack are obtained by solving a set of simultaneous singular integral equations in terms of the normal and tangent crack dislocation densities. The singularity of the dynamic stress fields near the crack tips is embodied in the fundamental solutions of the singular integral equations. The investigation of the fundamental solutions indicates that the singularity and oscillation indices of the stress fields are both dependent upon the material constants and the frequency of the harmonic load. This observation is different from the well-known -1/2 oscillating singularity for elastic bi-materials. The explanation for the differences between viscoelastic and elastic bi-materials can be given by the additional viscosity mismatch in the case of viscoelastic bi-materials. As an example, the standard linear solid model of a viscoelastic material is used. The effects of the frequency and the material constants (short-term modulus, long-term modulus and relaxation time) on the singularity and the oscillation indices are studied numerically.

Fatigue-crack closure measurements on 2024-t3 sheet specimens

A computer-controlled procedure has been developed for automatic measurement of the crack opening stress S-op during fatigue tests. A crack opening displacement gauge (GOD meter) is used to obtain digital data on the load versus COD curves. Three methods for deriving S-op from the data sets are compared: (1) a slope method, (2) a tangent lines intersecting method, and (3) a tangent point method. The effect of the position of the COD meter with respect to the crack tip on S-op is studied in tests of 2024-T3 specimens. Results of crack growth and S-op are presented for CA loading with an overload, and with an overload followed by an underload.

Discussion on the geometric structure of strange attractor

We discuss the transversal heteroclinic cycle formed by hyperbolic periodic pointes of diffeomorphism on the differential manifold. We point out that every possible kind of transversal heteroclinic cycle has the Smalehorse property and the unstable manifolds of hyperbolic periodic points have the closure relation mutually. Therefore the strange attractor may be the closure of unstable manifolds of a countable number of hyperbolic periodic points. The Henon mapping is used as an example to show that the conclusion is reasonable.