Elastodynamic stress intensity factor history for a semi-infinite crack under three-dimensional transient loading

The dynamic stress intensity factor history for a semi-infinite crack in an otherwise unbounded elastic body is analyzed. The crack is subjected to a pair of suddenly-applied point loadings on its faces at a distance L away from the crack tip. The exact expression for the mode I stress intensity factor as a function of time is obtained. The method of solution is based on the direct application of integral transforms, the Wiener-Hopf technique and the Cagniard-de Hoop method. Due to the existence of the characteristic length in loading this problem was long believed a knotty problem. Some features of the solutions are discussed and graphical result for numerical computation is presented.

爆轰驱动激波管缝合激波马赫数计算

采用数值求解的方法得出了氢氧爆轰驱动激波管的缝合状态参数。以空气为试验气体时，缝合激波马赫数随着H2摩尔浓度的增加而增加，H_2摩尔浓度达到90％左右时达到最大。当缝合马赫数较高时，需要考虑高温真实气体效应的影响，缝合激波马赫数较理想气体的高。以氢空气混合气为试验气体时，缝合激波马赫数较以空气为试验气体的小。通过调整驱动气体与被驱动气体的初始参数，可以得到即能恰好消除Taylor波又能缝合的运行状态。

用钛宝石再生放大器产生高重复率啁啾脉冲列

针对啁啾脉冲放大技术建成的钛宝石激光装置，提出一种获得高重复率激光脉冲列的方法．通过改变钛宝石再生放大器中泡克耳斯盒电光开关的传统工作模式，使得腔内放大的脉冲从某特定时刻起，每当在腔内往返一次就以一定的倒出比例（倒出率）倒出腔内脉冲能量的一部分，从而可以在有限的时间段内产生高重复率的啁啾激光脉冲列．基于Franz-Nodvik放大理论，建立了该高重复率再生放大器的理论模型，通过数值计算，系统地分析了初始增益、倒出时刻、倒出率对输出的脉冲序列的影响．在抽运功率为35mJ、倒出率为1／2的实验条件下，通过腔外

Charge displacement in ultraintense laser-plasma interaction

The distribution of optical held and charge density in the interaction between ultraintense ultrashort pulse laser and plasma is studied by numerical computation. The plasma considered has an exponential density profile. which corresponds to isothermal expanding. Our calculation shows that electrons are pushed forward by the incident laser, but ions, due to their much greater inertia, remain stationary. The resulting charge displacement forms a strong electrostatic field in the plasma. After the interaction of laser pulse and plasma. electrostatic energy still exists even after the laser pulse and will be absorbed by the plasma finally. This serves as an explanation to the mechanism of laser energy deposited into plasma.

MODIFIED DIRECT TWOS-COMPLEMENT PARALLEL ARRAY MULTIPLICATION ALGORITHM FOR COMPLEX MATRIX OPERATION

A direct twos-complement parallel array multiplication algorithm is introduced and modified for digital optical numerical computation. The modified version overcomes the problems encountered in the conventional optical twos-complement algorithm. In the array, all the summands are generated in parallel, and the relevant summands having the same weights are added simultaneously without carries, resulting in the product expressed in a mixed twos-complement system. In a two-stage array, complex multiplication is possible with using four real subarrays. Furthermore, with a three-stage array architecture, complex matrix operation is straightforwardly accomplished. In the experiment, parallel two-stage array complex multiplication with liquid-crystal panels is demonstrated.

A mobilidade eletroforética e o perfil de potencial da membrana celular

O objetivo do presente trabalho foi estudar o comportamento dos potenciais superficiais e do perfil de potencial atraves da membrana de eritr ocito em func ao da forca i onica e das cargas superficiais, usando um modelo que leva em conta as cargas el etricas do glicoc alix e das proteınas citoplasm aticas, al em das cargas superficiais da bicamada lipıdica e os efeitos dos eletr olitos divalentes. Programas especıficos em linguagem C foram elaborados para o c alculo desses potenciais, tomando como dados num ericos resultados experimentais de medidas de mobilidade eletrofor etica de eritr ocitos para diferentes valores de forca i onica. Neste c alculo, o metodo para tratamento dos dados eletrofor eticos indicado por Hsu et al.[57] foi incluıdo em nosso modelo. A equac ao de Poisson-Boltzmann nao linear foi resolvida por computac ao num erica, usando o metodo de Runge-Kutta de quarta ordem, obtendo-se os perfis de potencial. Os resultados mostraram que a estimativa da densidade de carga el etrica na superfıcie de c elulas usando a equac ao cl assica de Helmholtz-Smoluchowski conduz a valores que nao conseguem refletir as forcas que regem o comportamento eletrofor etico das mesmas. O presente modelo gerou valores de potenciais superficiais e perfis de potencial para a membrana do eritr ocito bem distintos daqueles obtidos anteriormente para um modelo descrito por uma equac ao de Poisson-Boltzmann linear. Nossos resultados confirmam que a avaliac ao de parametros el etricos superficiais da membrana de eritr ocito, envolvendo dados oriundos de eletroforese, deve incluir c alculos hidrodin amicos al em de eletroest aticos, como sugerido por Hsu et al. [57].

Estudo do comportamento caótico e determinação de dimensão fractal em modelos pré-inflacionários não compactos

O caos determinístico é um dos aspectos mais interessantes no que diz respeito à teoria moderna dos sistemas dinâmicos, e está intrinsecamente associado a pequenas variações nas condições iniciais de um dado modelo. Neste trabalho, é feito um estudo acerca do comportamento caótico em dois casos específicos. Primeiramente, estudam-se modelos préinflacionários não-compactos de Friedmann-Robertson-Walker com campo escalar minimamente acoplado e, em seguida, modelos anisotrópicos de Bianchi IX. Em ambos os casos, o componente material é um fluido perfeito. Tais modelos possuem constante cosmológica e podem ser estudados através de uma descrição unificada, a partir de transformações de variáveis convenientes. Estes sistemas possuem estruturas similares no espaço de fases, denominadas centros-sela, que fazem com que as soluções estejam contidas em hipersuperfícies cuja topologia é cilíndrica. Estas estruturas dominam a relação entre colapso e escape para a inflação, que podem ser tratadas como bacias cuja fronteira pode ser fractal, e que podem ser associadas a uma estrutura denominada repulsor estranho. Utilizando o método de contagem de caixas, são calculadas as dimensões características das fronteiras nos modelos, o que envolve técnicas e algoritmos de computação numérica, e tal método permite estudar o escape caótico para a inflação.

环形非稳腔倒向波抑制镜抑制效果的数值模拟

对环形非稳腔激光器做数值计算模拟，给带有倒向波抑制(RWS)镜的环形非稳腔激光器的设计提供了数值参考。通过自洽方程，用Matlab计算带有倒向波抑制镜的环形非稳腔正、倒向波的稳定过程。改变腔长、放大率等参数来计算，发现使用倒向波抑制镜能明显提高正向波的输出，并能抑制倒向波。气体分子的多普勒频移导致正、倒向波的频率变化。频率分开的越小，竞争越大，抑制效果越好，正向波输出功率则有所下降。对于特定的结构，存在最佳的放大率，使正向波的输出最大。计算结果表明，倒向波抑制镜的引入能增大环形非稳腔输出功率。计算结果对实

飞秒激光用多层介质膜脉宽压缩光栅的设计

以薄膜光学的干涉理论和衍射光学的傅里叶模式理论为基础，给出了0．8μm飞秒激光器用多层介质膜脉宽压缩光栅的理论设计；设计采用H3L（HL）^9H0．5L2．4H的多层介质膜为基底，当刻蚀后表面浮雕结构的占宽比为0．35，线密度为1480线／mm，槽深为0．2μm，顶层HfO2的剩余厚度为0．15μm时，对于Littrow角度（36．7°）和TE波模式入射的衍射光栅其-1级衍射效率达到95％以上．

Approximate analytical optical transfer function for a wavefront-coded imaging system with two-dimensional rectangularly separable phase masks

Wavefront coding can be used to extend the depth of field of incoherent imaging systems and is a powerful system-level technique. In order to assess the performance of a wavefront-coded imaging system, defocused optical transfer function (OTF) is the metric frequently used. Unfortunately, to the best of our knowledge, among all types of phase masks, it is usually difficult to obtain the analytical OTF except the cubic one. Although numerical computation seems good enough for performance evaluation, the approximate analytical OTF is still indispensable because it can reflect the relationship between mask parameters and system frequency response in a clearer way. Thus, a method is proposed to derive the approximate analytical OTF for two-dimensional rectangularly separable phase masks. The analytical results are well consistent with the direct numerical computations, but the proposed method can be accepted only from engineering point of view and needs rigorous proof in future. (c) 2010 Society of Photo-Optical Instrumentation Engineers. [DOI: 10.1117/1.3485759]

Effects of Fracture Seepage on the Stability of Landslide during Reservoir Water Level Fluctuation

The hydraulic conductivity function of fractures is a key scientific question to describe and reveal the process and the role of water seepage reasonably. In this paper, the generation technology of random fracture network and the latest numerical computation method for equivalent permeability tensor of fracture network are applied to analyze the landslide located at Wangjiayuanzi in Wanzhou District of Chongqing by simulating the changes of the seepage field caused by the running of the Three Gorges Reservoir. The influences of the fracture seepage on the seepage field and stability of the landslide were discussed with emphasis. The results show that the fractures existing in the soil increase the permeability coefficient of the landslide body and reduce the delay time of the underground water level in the landslide which fluctuates relative to the water level of reservoir，that causes the safe coefficient of the slope changes more gently than that of the same slope without fractures. It means, if only water level fluctuating condition is concerned, the fractures existing in the soil plays a positive role to the stability of slopes.

Analysis of chemical calorific effect during reactive extrusion processes for free radical polymerization

In the reactive extrusion process for polymerization, the chemical calorific effect has a great influence on the temperature. In order to quantitatively analyze the polymerization trend and optimize the processing conditions, the phenomena of the chemical calorific effect during reactive extrusion processes for free radical polymerization were analyzed. Numerical computation expressions of the heat of chemical reaction and the reactive calorific intensity were deduced, and then a numerical simulation of the reactive extrusion process for the polymerization of n-butyl methacrylate was carried out. The evolutions of the heat of chemical reaction and the reactive calorific intensity along the! axial direction of the extruder are presented, on the basis of which reactive processing conditions can be optimized.

Computer simulation of reactive extrusion processes for free radical polymerization

The reactive extrusion for polymerization is an integrated polymer processing technology. A new semi-implicit iterative algorithm was proposed to deal with the complicated relationships among the chemical reaction, the macromolecular structure and the chemorheological property. Then the numerical computation expressions of the average molecular weight, the monomer conversion, and the initiator concentration were deduced, and the computer simulation of the reactive extrusion process for free radical polymerization was carried out, on basis of which reactive processing conditions can be optimized.

Free radical grafting of polyethylene with vinyl monomers by reactive extrusion

The free radical grafting of polyethylene with vinyl monomers by reactive extrusion was studied numerically. Numerical computation expressions of key variables, such as the concentrations of the initiator and polymer, grafting degree, average molecular weight and apparent viscosity, were deduced. The evolutions of the above variables were predicted by means of an uncoupled semi-implicit iterative algorithm. The monomer conversion monotonically increases with decreasing throughput or increasing initial initiator concentration; with increasing barrel temperature, the monomer conversion first increases then decreases. The simulated results are nearly in good agreement with the experimental results.

An extended variable-grid global ocean circulation model and its preliminary results of the equatorial Pacific circulation

To investigate the interaction between the tropical Pacific and China seas a variable-grid global ocean circulation model with fine grid[(1/6)degrees] covering the area from 20degreesS to 50degreesN and from 99degrees to 150degreesE is developed. Numerical computation of the annually cyclic circulation fields is performed. The results of the annual mean zonal currents and deep to abyssal western boundary currents in the equatorial Pacific Ocean are reported. The North Equatorial Current,the North Equatorial Countercurrent, the South Equatorial Current and the Equatorial Undercurrent are fairly well simulated. The model well reproduces the northward flowing abyssal western boundary current. From the model results a lower deep western boundary current east of the Bismarck-Solomon-New Hebrides Island chain at depths around 2 000 in has been found. The model results also show that the currents in the equatorial Pacific Ocean have multi-layer structures both in zonal currents and western boundary currents, indicating that the global ocean overturning thermohaline circulation appears of multi-layer pattern.