SERPENT和SAFER密码算法的能量攻击

SERPENT和SAFER是AES的两个候选算法 ,本文使用能量攻击方法对它们进行了深入分析 ,结果表明 :对于 2 5 6、192和 12 8比特密钥的SERPENT算法 ,能量攻击平均需分别进行 2 159、2 119和 2 79次试验 .虽然所需的试验次数实际没法达到 ,但是此攻击方法大大地降低了SERPENT的密钥规模 ,并且发现对于能量攻击 ,SERPENT有许多弱密钥 .经过深入分析和穷尽搜索可知 :能量攻击可以获取SAFER的种子密钥 .文中还给出了两种抵抗能量攻击的SER PENT的改进密钥方案以及设计密钥方案时需注意的问题 .

对低轮SAFER++的差分-非线性密码分析

SAFER++是进入NESSIE第2轮评估的7个分组算法之一。采用差分密码分析和非线性密码分析相结合的方法对4轮、5轮和6轮SAFER++进行分析，结果表明：6轮SAFER++对这种攻击方法不免疫；攻击4轮和5轮SAFER++时，与已有结果相比，攻击复杂度大大减小。攻击对2\+{250}个256比特长度的密钥有效。

SAFER系列密码算法的设计与分析

SAFER系列密码算法的总体结构采用SP-网络,它的设计具有其独到的几个特色.分析SAFER系列密码算法的设计思想,沿着设计者对它们不断改进的思路,分别描述其混淆层、扩散层、密钥扩展算法的性质和对它们的攻击.最后提出几个尚需进一步考虑的问题.

Stress Wave Propagation in a Gradient Elastic Medium

The gradient elastic constitutive equation incorporating the second gradient of the strains is used to determine the monochromatic elastic plane wave propagation in a gradient infinite medium and thin rod. The equation of motion, together with the internal material length, has been derived. Various dispersion relations have been determined. We present explicit expressions for the relationship between various wave speeds, wavenumber and internal material length.

Two critical crack propagating velocities for PMMA fracture surface

Ultrasonic fractography and scanning electronic microscopy (SEM) are used to determine the direct relationship between the fracture surface morphology and the main crack velocity during the rapid rupture of polymethylmethacrylate (PMMA). Two critical crack velocities are found for the fracture. Quasi-parabolic markings will appear when the crack speed exceeds the first critical speed. Crack propagating at speed above the second critical speed leaves a thicket of small branches penetrating the surface behind them. Both critical speeds are functions of the thickness of the specimens.

Influence of Obstacle-Produced Turbulence on Development of Premixed Flames

An investigation into influence of obstructions on premixed flame propagation has been carried out in a semi-open tube. It is found that there exists flame acceleration and rising overpressure along the path of flame due to obstacles. According to the magnitude of flame speeds, the propagation of flame in the tube can be classified into three regimes: the quenching, the choking and the detonation regimes. In premixed flames near the flammability limits, the flame is observed first to accelerate and then to quench itself after propagating past a certain number of obstacles. In the choking regime, the maximum flame speeds are somewhat below the combustion product sound speeds, and insensitive to the blockage ratio. In the more sensitive mixtures, the transition to detonation (DDT) occurs when the equivalence ratio increases. The transition is not observed for the less sensitive mixtures. The dependence of overpressure on blockage ratio is not monotonous. Furthermore, a numerical study of flame acceleration and overpressure with the unsteady compressible flow model is performed, and the agreement between the simulation and measurements is good.

Application of Three-Dimensional Discrete Element Face-to-Face Contact Model with Fissure Water Pressure to Stability Analysis of Landslide in Panluo Iron Mine

Three-dimensional discrete element face-to-face contact model with fissure water pressure is established in this paper and the model is used to simulate three-stage process of landslide under fissure water pressure in the opencast mine, according to the actual state of landslide in Panluo iron mine where landslide happened in 1990 and was fathered in 1999. The calculation results show that fissure water pressure on the sliding surface is the main reason causing landslide and the local soft interlayer weakens the stability of slope. If the discrete element method adopts the same assumption as the limit equilibrium method, the results of two methods are in good agreement; while if the assumption is not adopted in the discrete element method, the critical phi numerically calculated is less than the one calculated by use of the limit equilibrium method for the same C. Thus, from an engineering point of view, the result from the discrete element model simulation is safer and has more widely application since the discrete element model takes into account the effect of rock mass structures.

Electromechanical influence of crack velocity at bifurcation for poled ferroelectric materials

The piezoelastodynamic field equations are solved to determine the crack velocity at bifurcation for poled ferroelectric materials where the applied electrical field and mechanical stress can be varied. The underlying physical mechanism, however, may not correspond to that assumed in the analytical model. Bifurcation has been related to the occurrence of a pair of maximum circumferential stress oriented symmetrically about the moving crack path. The velocity at which this behavior prevails has been referred to as the limiting crack speed. Unlike the classical approach, bifurcation will be identified with finite distances ahead of a moving crack. Nucleation of microcracks can thus be modelled in a single formulation. This can be accomplished by using the energy density function where fracture initiation is identified with dominance of dilatation in relation to distortion. Poled ferroelectric materials are selected for this study because the microstructure effects for this class of materials can be readily reflected by the elastic, piezoelectic and dielectric permittivity constants at the macroscopic scale. Existing test data could also shed light on the trend of the analytical predictions. Numerical results are thus computed for PZT-4 and compared with those for PZT-6B in an effort to show whether the branching behavior would be affected by the difference in the material microstructures. A range of crack bifurcation speed upsilon(b) is found for different r/a and E/sigma ratios. Here, r and a stand for the radial distance and half crack length, respectively, while E and a for the electric field and mechanical stress. For PZT-6B with upsilon(b) in the range 100-1700 m/s, the bifurcation angles varied from +/-6degrees to +/-39degrees. This corresponds to E/sigma of -0.072 to 0.024 V m/N. At the same distance r/a = 0.1, PZT-4 gives upsilon(b) values of 1100-2100 m/s; bifurcation angles of +/-15degrees to +/-49degrees; and E/sigma of -0.056 to 0.059 V m/N. In general, the bifurcation angles +/-theta(0) are found to decrease with decreasing crack velocity as the distance r/a is increased. Relatively speaking, the speed upsilon(b) and angles +/-theta(0) for PZT-4 are much greater than those for PZT-6B. This may be attributed to the high electromechanical coupling effect of PZT-4. Using upsilon(b)(0) as a base reference, an equality relation upsilon(b)(-) < upsilon(b)(0) < upsilon(b)(+) can be established. The superscripts -, 0 and + refer, respectively, to negative, zero and positive electric field. This is reminiscent of the enhancement and retardation of crack growth behavior due to change in poling direction. Bifurcation characteristics are found to be somewhat erratic when r/a approaches the range 10(-2)-10(-1) where the kinetic energy densities would fluctuate and then rise as the distance from the moving crack is increased. This is an artifact introduced by the far away condition of non-vanishing particle velocity. A finite kinetic energy density prevails at infinity unless it is made to vanish in the boundary value problem. Future works are recommended to further clarify the physical mechanism(s) associated with bifurcation by means of analysis and experiment. Damage at the microscopic level needs to be addressed since it has been known to affect the macrocrack speeds and bifurcation characteristics. (C) 2002 Published by Elsevier Science Ltd.

Analysis of asymmetrical cold rolling with varying coefficients of friction

In this research, asymmetrical cold rolling was produced by the difference in the coefficient of friction between rolls and sheets rather than the difference of roll radius or rotation speeds. The influence of friction coefficient ratio on the cross shear deformation, rolling pressure and torque was investigated using slab analysis. The results showed that the shear deformation zone length increased with the increase of the friction coefficient ratio. The rolling force decreased only under the condition that the friction coefficient ratio increased while the sum of the friction coefficients was held constant. As the reduction per pass was increased, the shear deformation zone length increased and the rolling force also increased. An increase of the front tension resulted in a decrease of the shear deformation zone length. An increase of back tension, however, led to an increase of the shear deformation zone length. The reduction of rolling torque for the work roll with higher surface roughness was greater than that for the work roll with lower surface roughness. (C) 2002 Elsevier Science B.V. All rights reserved.

Air flow over a mountain and the convective boundary-layer

In this paper, the analytical model coupling the convective boundary layer (CBL) with the free atmosphere developed by Qi and Fu (1992) is improved. And by this improved model, the interaction between airflow over a mountain and the CBL is further discussed. The conclusions demonstrate: (1) The perturbation potential temperatures in the free atmosphere can counteract the effect of orographic thermal forcing through entraining and mixing in the CBL. If u(M)BAR > u(F)BAR, the feedback of the perturbation potential temperatures in the free atmosphere is more important than orographic thermal forcing, which promotes the effect of interfacial waves. If u(M)BAR < u(F)BAR, orographic thermal forcing is more important, which makes the interfacial height and the topographic height identical in phase, and the horizontal speeds are a maximum at the top of the mountain. (2) The internal gravity waves propagating vertically in the free atmosphere cause a strong downslope wind to become established above the lee slope in the CBL and result in the hydraulic jump at the top of the CBL. (3) With the CBL deepening, the interfacial gravity waves induced by the potential temperature jump at the top of the CBL cause the airflow in the CBL to be subcritical.

Oscillatory Flow Induced Hydrodynamic Forces upon a Pipeline near Erosive Sandy Seabed

An experimental study of the properties of hydrodynamic forces upon a marine pipeline is presented in this paper, in the equilibrium scour conditions for various Keulegan-Carpenter numbers and various initial relative gaps between pipeline and the erosive sandy seabed. The tests are conducted in a U-shaped oscillatory water tunnel with a sand box located at the bottom of the test section. According to the experimental results, the maximum horizontal forces on the pipelines with an initial gap to seabed will decrease to some extent due to scouring process. For engineering appliances, it seems safer to estimate wave induced forces on pipelines under the assumption that seabed is plane. However, it should be noticed that great changes would be brought to the frequency properties of lift forces because of the sandy scour beneath the pipeline, which occurs for certain KC numbers.

Numerical Analyses of Transonic Fluid/Structure interaction in Aircraft Engineering

Through the coupling between aerodynamic and structural governing equations, a fully implicit multiblock aeroelastic solver was developed for transonic fluid/stricture interaction. The Navier-Stokes fluid equations are solved based on LU-SGS (lower-upper symmetric Gauss-Seidel) Time-marching subiteration scheme and HLLEW (Harten-Lax-van Leer-Einfeldt-Wada) spacing discretization scheme and the same subiteration formulation is applied directly to the structural equations of motion in generalized coordinates. Transfinite interpolation (TFI) is used for the grid deformation of blocks neighboring the flexible surfaces. The infinite plate spline (IPS) and the principal of virtual work are utilized for the data transformation between fluid and structure. The developed code was fort validated through the comparison of experimental and computational results for the AGARD 445.6 standard aeroelastic wing. In the subsonic and transonic range, the calculated flutter speeds and frequencies agree well with experimental data, however, in the supersonic range, the present calculation overpredicts the experimental flutter points similar to other computations. Then the flutter character of a complete aircraft configuration is analyzed through the calculation of the change of structural stiffness. Finally, the phenomenon of aileron buzz is simulated for the weakened model of a supersonic transport wing/body model at Mach numbers of 0.98 and l.05. The calculated unsteady flow shows, on the upper surface, the shock wave becomes stronger as the aileron deflects downward, and the flow behaves just contrary on the lower surface of the wing. Corresponding to general theoretical analysis, the flow instability referred to as aileron buzz is induced by a stronger shock alternately moving on the upper and lower surfaces of wing. For the rigid structural model, the flow is stable at all calculated Mach numbers as observed in experiment

Influence of Wind and Grain Size on Migration of Asymmetric Sand Waves

Large parts of shallow seas are covered by regular seabed patterns and sand wave is one kind of these patterns. The instability of the sedimentary structures may hazard pipelines and the foundations of offshore structures. In the last decade or so, it's a focus for engineers to investigate the movement mechanism of sand waves. Previous theoretical studies of the subject have developed a general model to predict the growth and migration of sand waves, which is based on the two-dimensional vertical shallow water equations and the bed-form deformation equations. Although the relation between wave-current flow and sand bed deformation has been established, the topographic influence has not been considered in the model. In this paper some special patterns, which are asymmetric and close to the reality, are represent as the perturbed seabed and the evolution of sand waves is calculated. The combination of a steady flow induced by wind and a sinusoidal tidal flow is considered as the basic flow. Finally the relations of some parameters (grain size, etc.) and sand waves' growth and migration are discussed, and the growth rate and migration speeds of asymmetric sand waves are carried out.