On dilatational plastic constitutive equation of ductile materials and plastic loading paths at bifurcation

The dilatational plastic constitutive equation presented in this paper is proved to be in a form of generality. Based on this equation, the constitutive behaviour of materials at the moment of bifurcation is demonstrated to follow a loading path with the response as "soft" as possible.

Evolution and propagation of density waves on galactic disk .1. General evolution equation of wavelet and discussion of its solution

This paper presents a general self-consistent theory of evolution and propagation of wavelets on the galactic disk. A simplified model for this theory, i. e. the thin transition-layer approximation is proposed.There are three types of solutions to the basic equation governing the evolution of wavelets on the disk: (ⅰ) normal propagating type; (ⅱ) swing type; (ⅲ) general evolving type. The results show that the first two types are applicable to a certain domain on the galactic disk and a certain region of the wave number of wavelets. The third is needed to join the other two types and to yield a coherent total picture of the wave motion. From the present theory, it can be seen that the well-known "swing theory" of the G-L sheet model holds only for a certain class of basic states of galaxies.

The wave model of mesothermal plasma near wakes and Korteweg-de Vries equation

The stationary two-dimensional (x, z) near wakes behind a flat-based projectile which moves at a constant mesothermal speed (V∞) along a z-axis in a rarefied, fully ionized, plasma is studied using the wave model previously proposed by one of the authors (VCL). One-fluid theory is used to depict the free expansion of ambient plasma into the vacuum produced behind a fast-moving projectile. This nonstationary, one-dimensional (x, t) flow which is approximated by the K-dV equation can be transformed, through substitution, t=z/V∞, into a stationary two-dimensional (x, z) near wake flow seen by an observer moving with the body velocity (V∞). The initial value problem of the K-dV equation in (x, t) variables is solved by a specially devised numerical method. Comparisons of the present numerical solution for the asymptotically small and large times with available analytical solutions are made and found in satisfactory agreements.

加卸载响应比在Poisson模型下的随机分布

本文在地震的发生时间服从Ｐｏｉｓｏｎ过程，而地震震级服从ＧｕｔｅｎｂｅｒｇＲｉｃｈｔｅｒ关系的前提下，对不同定义的加卸载响应比Ｙ值的随机分布进行了探讨。结果表明：当在计算窗口的地震发生的期望数目较大（＞４０）时，Ｙ１～Ｙ５值的分布基本稳定，出现高加卸载响应比的概率极低。然而当计算窗口的地震期望数目过小时，Ｙ２～Ｙ５值则变得不太稳定。也就是说，服从Ｐｏｉｓｏｎ过程的地震序列，在计算窗口的地震期望数目过小时，也可能产生Ｙ值较高的结果。为了使利用加卸载响应比预测地震更加可靠，文中给出了Ｙ１、Ｙ３在Ｐｏｉｓｏｎ模型下的９０％、９５％和９９％的置信区间，这对判别加卸载响应比异常是非常有用的。

Starting Nozzle Flow Simulation Using K-G Two-Equation Turbulence Model

The starting process of two-dimensional nozzle flows has been simulated with Euler, laminar and k - g two-equation turbulence Navier-Stokes equations. The flow solver is based on a combination of LUSGS subiteration implicit method and five spatial discretized schemes, which are Roe, HLLE, MHLLE upwind schemes and AUSM+, AUSMPW schemes. In the paper, special attention is for the flow differences of the nozzle starting process obtained from different governing equations and different schemes. Two nozzle flows, previously investigated experimentally and numerically by other researchers, are chosen as our examples. The calculated results indicate the carbuncle phenomenon and unphysical oscillations appear more or less near a wall or behind strong shock wave except using HLLE scheme, and these unphysical phenomena become more seriously with the increase of Mach number. Comparing the turbulence calculation, inviscid solution cannot simulate the wall flow separation and the laminar solution shows some different flow characteristics in the regions of flow separation and near wall.

A Semi-Empirical Equation of Penetration Depth on Concrete Target Impacted by Ogive-Nose Projectiles

In this paper, the penetration process of ogive-nose projectiles into the semi-infinite concrete target is investigated by the dimensional analysis method and FEM simulation. With the dimensional analysis, main non-dimensional parameters which control the penetration depth are obtained with some reasonable hypothesis. Then, a new semi-empirical equation is present based on the original work of Forrestal et al., has only two non-dimensional combined variables with definite physical meanings. To verify this equation, prediction results are compared with experiments in a wide variation region of velocity. Then, a commercial FEM code, LS-DYNA, is used to simulate the complex penetration process, that also show the novel semi-empirical equation is reasonable for determining the penetration depth in a concrete target.

Modified Simplified Rate Equation Model of Flowing Chemical Oxygen-Iodine Laser and its Application

A modified simplified rate equation (RE) model of flowing chemical oxygen-iodine laser (COIL), which is adapted to both the condition of homogeneous broadening and inhomogeneous broadening being of importance and the condition of inhomogeneous broadening being predominant, is presented for performance analyses of a COIL. By using the Voigt profile function and the gain-equal-loss approximation, a gain expression has been deduced from the rate equations of upper and lower level laser species. This gain expression is adapted to the conditions of very low gas pressure up to quite high pressure and can deal with the condition of lasing frequency being not equal to the central one of spectral profile. The expressions of output power and extraction efficiency in a flowing COIL can be obtained by solving the coupling equations of the deduced gain expression and the energy equation which expresses the complete transformation of the energy stored in singlet delta state oxygen into laser energy. By using these expressions, the RotoCOIL experiment is simulated, and obtained results agree well with experiment data. Effects of various adjustable parameters on the performances of COIL are also presented.

Numerical-Value Perturbation Method with Application of Solving Convective-Diffusion Integral Equation

The convective--diffusion equation is of primary importance in such fields as fluid dynamics and heat transfer hi the numerical methods solving the convective-diffusion equation, the finite volume method can use conveniently diversified grids (structured and unstructured grids) and is suitable for very complex geometry The disadvantage of FV methods compared to the finite difference method is that FV-methods of order higher than second are more difficult to develop in three-dimensional cases. The second-order central scheme (2cs) offers a good compromise among accuracy, simplicity and efficiency, however, it will produce oscillatory solutions when the grid Reynolds numbers are large and then very fine grids are required to obtain accurate solution. The simplest first-order upwind (IUW) scheme satisfies the convective boundedness criteria, however. Its numerical diffusion is large. The power-law scheme, QMCK and second-order upwind (2UW) schemes are also often used in some commercial codes. Their numerical accurate are roughly consistent with that of ZCS. Therefore, it is meaningful to offer higher-accurate three point FV scheme. In this paper, the numerical-value perturbational method suggested by Zhi Gao is used to develop an upwind and mixed FV scheme using any higher-order interpolation and second-order integration approximations, which is called perturbational finite volume (PFV) scheme. The PFV scheme uses the least nodes similar to the standard three-point schemes, namely, the number of the nodes needed equals to unity plus the face-number of the control volume. For instanc6, in the two-dimensional (2-D) case, only four nodes for the triangle grids and five nodes for the Cartesian grids are utilized, respectively. The PFV scheme is applied on a number of 1-D problems, 2~Dand 3-D flow model equations. Comparing with other standard three-point schemes, The PFV scheme has much smaller numerical diffusion than the first-order upwind (IUW) scheme, its numerical accuracy are also higher than the second-order central scheme (2CS), the power-law scheme (PLS), the QUICK scheme and the second-order upwind(ZUW) scheme.

非传播孤立波的数值模拟

非传播孤立波是近年来由中国学者发现的一种独特的孤立水波。本文通过数值求解非传播孤立波目前公认的控制方程-Miles导出的一个带共轭项的非线性立方Schrodinger方程，对非传播孤立波进行数值模拟。本文针对非传播孤立波的各种性质，作了大量的数值计算工作，并与实验观察的现象及人们对非传播孤立波的理论研究结果进行了比较和分析。为了得到稳定的非传播孤立波，本文讨论了Miles方程中的线性阻尼系数a的值，计算表明，线性阻尼a对能否形成稳定的非传播孤立波影响很大，在某些情况下，Laedke等人提出的Miles方程的非传播孤立波解的稳定性条件与我们对Miles方程的数值模拟的结果相当一致，a可在满足稳定性条件的区间内取值，但也发现在某些情况下Laedke等人的稳定性条件与我们的数值模拟不完全符合，证明Laedke等人关于非传播孤立波的稳定性条件只是一个必要条件，而不是充分条件。本文研究了两个非传播孤立波的相互作用，数值模拟表明，两个波的作用模式依赖于系统的参数，只有适当大小的外驱动频率和振幅及线性阻尼a可算出两个非传播孤立波周而复始的相互作用现象来，参数不合适时，两个波可能最终合并为一个非传播孤立波而不再分离，也可能彼此不发生作用，保持各自的独立。对不同的初始扰动及其演化的计算表明，要形成单个稳定的非传播孤立波，则初始扰动必须适当，否则扰动可能消失或发展成多个孤立波。关于形成非传播孤立波所需的外驱动条件，计算结果表明，只有适当大小的外驱动频率和振幅可形成稳定的非传播孤立波，数值结果可以描述驱动频率的上限和驱动振幅的上下限，但无法描述驱动频率的下限。我们的数值模拟工作说明Miles方程确实较好的描述了非传播孤立波的物理模型，该方程可以解释许多关于非传播孤立波的物理特性。但Miles方程无法对非传播孤立波的某些实验现象作出解释，因而有待于进一步研究改进。

Solute transport in nanochannels with roughness-like structures

Newfound attention has been given to solute transport in nanochannels. Because the electric double layer (EDL) thickness is comparable to characteristic channel dimensions, nanochannels have been used to separate ionic species with a constant charge-to-size ratio (i.e., electrophoretic mobility) that otherwise cannot be separated in electroosmotic or pressure- driven flow along microchannels. In nanochannels, the electrical fields within the EDL cause transverse ion distributions and thus yield charge-dependent mean ion speeds in the flow. Surface roughness is usually inevitable during microfabrication of microchannels or nanochannels. Surface roughness is usually inevitable during the fabrication of nanochannels. In the present study, we develop a numerical model to investigate the transport of charged solutes in nanochannels with hundreds of roughness-like structures. The model is based on continuum theory that couples Navier-Stokes equations for flows, Poisson-Boltzmann equation for electrical fields, and Nernst-Planck equation for solute transports. Different operating conditions are considered and the solute transport patterns in rough channels are compared with those in smooth channels. Results indicate that solutes move slower in rough nanochannels than in smooth ones for both pressure- driven and electroosmotic flows. Moreover, solute separation can be significantly improved by surface roughness under certain circumstances.

Verification of a scaling law in few-cycle laser pulses

The scaling law of photoionization in few-cycle laser pulses is verified in this paper. By means of numerical solution of time-dependent Schrodinger equation, the photoionization and the asymmetry degree of photoionization of atoms with different binding potential irradiated by various laser pulses are studied. We find that the effect of increasing pulse intensity is compensated by deepening the atomic binding potential. In order to keep the asymmetric photoionization unchanged, if the central frequency of the pulse is enlarged by k times, the atomic binding potential should also be enlarged by k times, and the laser intensity should be enlarged by k(3) times. (c) 2005 Optical Society of America.

A scaling law of photoionization in ultrashort pulses

By means of the numerical solution of time-dependant Schrodinger equation, we verify a scaling law of photoionization in ultrashort pulses. We find that for a given carrier-envelope phase and duration of the pulse, identical photoionizations are obtained provided that when the central frequency of the pulse is enlarged by k times, the atomic binding potential is enlarged by k times, and the laser intensity is enlarged by k(3) times. The scaling law allows us to reach a significant control over direction of photoemission and offers exciting prospects of reaching similar physical processes in different interacting systems which constitutes a novel kind of coherent control.

μ^-子催化核聚变中强脉冲激光对介原子μ^3He的电离

数值求解了一维含时的Schroedinger方程，研究了μ子催化核聚变反应中激光强度和波长对介原子μ^3He电离的影响．发现当激光强度为10^19-10^23W/cm^2量级时，介原子μ^3He有2．7％左右的电离率；当激光强度达到6．0×10^24W/cm^2时，对介原子μ^3He有显著的电离，并且电离率随着激光的强度、波长而递增，进而会有效提高μ子的催化效率．

Nonlinear X-wave formation and conical emission at different powers of a femtosecond laser pulse in water

Nonlinear X-wave formation at different pulse powers in water is simulated using the standard model of nonlinear Schrodinger equation (NLSE). It is shown that in near field X-shape originally emerges from the interplay between radial diffraction and optical Kerr effect. At relatively low power group-velocity dispersion (GVD) arrests the collapse and leads to pulse splitting on axis. With high enough power, multi-photon ionization (NIPI) and multi-photon absorption (MPA) play great importance in arresting the collapse. The tailing part of pulse is first defocused by MPI and then refocuses. Pulse splitting on axis is a manifestation of this process. Double X-wave forms when the split sub-pulses are self-focusing. In the far field, the character of the central X structure of conical emission (CE) is directly related to the single or double X-shape in the near field. (c) 2007 Elsevier B.V. All rights reserved.