含节理岩块单轴受压试验三维离散元数值模拟

采用三维离散元方法，模拟含节理岩块的单轴受压试验。并针对节理空间分布，给出了解析解与数值计算比较的结果，验证了数值模拟的正确性。当节理正交时，离散元与有限元的计算结果一致。计算结果定量地说明了岩体的各向异性和尺寸效应。

超音速流动中侧向喷流干扰特性的数值模拟

通过与实验结果的比较验证Fluent软件可以用于超音速流动中侧向喷流干扰特性数值研究．使用Fluent软件对超音速流动中不同喷流压力和攻角下侧向喷流干扰特性进行的数值研究表明，超音速流动中，随喷流压力和负攻角增大，喷流前的高压区明显增大，喷流的控制效果更好．喷流包裹作用的影响范围在90°弹面内．喷流压力增大，喷流包裹作用加强，影响范围增大．

激光熔凝加工中瞬时温度场数值模拟

以激光熔凝表面强韧化处理为背景,应用空间弹塑性有限单元和高精度数值算法同时考虑材料组织性能的变化模拟工件的温度场,主要研究激光熔凝加工中瞬时温度场数值模拟,同时考虑相变潜热的影响,为第二步热应力场及残余应力的数值模拟做准备.用算例验证了模型的正确性,最后给出了激光熔凝加工不同时刻温度场分布.

气-液两相流界面迁移现象的数值模拟研究

借助于Level Set函数,建立了气-液两相的统一控制方程组,并在交错网格中进行离散.用两种格式,即Superbee-TVD格式和5阶WENO格式求解Level Set函数的输运方程,用SIMPLER算法的思想对主流场控制方程的求解方法进行改进.数值实验结果表明,在求解Level Set的控制方程时,5阶WENO方法比Superbee-TVD格式的结果更准确;用改进的数值算法可成功实现对密度比大于1 000/1的气-液两相流界面迁移问题的数值模拟.对几种典型大密度比气-液两相流问题的计算结果与实际问题的物理规律完全一致,验证了该方法的有效性和可靠性.

Oscillatory Instability of Rayleigh-Marangoni-Benard Convection in Two-Layer System

The convective instabilities in two or more superposed layers heated from below were studied extensively by many scientists due to several interfacial phenomena in nature and crystal growth application. Most works of them were performed mainly on the instability behaviors induced only by buoyancy force, especially on the oscillatory behavior at onset of convection (see Gershuni et. Al.(1982), Renardy et. Al. (1985,2000), Rasenat et. Al. (1989), and Colinet et. Al.(1994)) . But the unstable situations of multi-layer liquid convection will become more complicated and interesting while considering at the same time the buoyancy effect combined with thermocapillary effect. This is the case in the gravity reduced field or thin liquid layer where the thermocapillary effect is as important as buoyancy effect. The objective of this study was to investigate theoretically the interaction between Rayleigh-Bénard instability and pure Marangoni instability in a two-layer system, and more attention focus on the oscillatory instability both at the onset of convection and with increasing supercriticality. Oscillatory behavious of Rayleigh-Marangoni-Bénard convective instability (R-M-B instability) and flow patterns are presented in the two-layer system of Silicon Oil (10cSt) over Fluorinert (FC70) for a larger various range of two-layer depth ratios (Hr=Hupper/Hdown) from 0.2 to 5.0. Both linear instability analysis and 2D numerical simulation (A=L/H=10) show that the instability of the system depends strongly on the depth ratio of two-layer liquids. The oscillatory instability regime at the onset of R-M-B convection are found theoretically in different regions of layer thickness ratio for different two-layer depth H=12,6,4,3mm. The neutral stability curve of the system displaces to right while we consider the Marangoni effect at the interface in comparison with the Rayleigh-Bénard instability of the system without the Marangoni effect (Ma=0). The numerical results show different regimes of the developing of convection in the two-layer system for different thickness ratios and some differences at the onset of pure Marangoni convection and the onset of Rayleigh-Bénard convections in two-layer liquids. Both traveling wave and standing wave were detected in the oscillatory instability regime due to the competition between Rayleigh-Bénard instability and Marangoni effect. The mechanism of the standing wave formation in the system is presented numerically in this paper. The oscillating standing wave results in the competition of the intermediate Marangoni cell and the Rayleigh convective rolls. In the two-layer system of 47v2 silicone oil over water, a transition form the steady instability to the oscillatory instability of the Rayleigh-Marangoni-Bénard Convection was found numerically above the onset of convection for ε=0.9 and Hr=0.5. We propose that this oscillatory mechanism is possible to explain the experimental observation of Degen et. Al.(1998). Experimental work in comparison with our theoretical findings on the two-layer Rayleigh-Marangoni-Bénard convection with thinner depth for H<6mm will be carried out in the near future, and more attention will be paid to new oscillatory instability regimes possible in the influence of thermocapillary effects on the competition of two-layer liquids

Rayleigh-Marangoni-Bénard Convective Instability

The Rayleigh-Marangoni-Benard convective instability (R-M-B instability) and flow patterns in the two-layer system of silicon oil 10cSt and Fluorinert FC70 liquids are studied theoretically and experimentally. Both linear instability analysis and 2D numerical simulation (A=L/H=10) were performed to study the influence of thermocapillary force on the convective instability of the two-layer system. Time-dependent oscillations arising at the onset of convection were investigated in a larger various range of two-layer depth ratios (Hr=H1/H2) from 0.2 to 5.0 for different total depth less than 12mm. Our results are different from the previous study on the Rayleig-B閚ard instability and show the strong effects of thermocapillary force at the interface on the time-dependent oscillations at the onset of instability convection. Primary experimental results of the critical instability parameters and the convective structure in the R-M-B convection have been obtained by using the digital particle image velocimetry (DPIV) system, and a good agreement in comparison with the results of numerical simulation was obtained.

Transient Temperature Distribution in a Sheet Metal by Low-Energy Laser Scanning

Low-energy laser-heating techniques are widely used in engineering applications such as, thinfilm deposition, surface treatment, metal forming and micro-structural pattern formation. In this paper,under the conditions of ignoring the thermo-mechanical coupling, a numerical simulation on the spatialand temporal temperature distribution in a sheet metal produced by the laser beam scanning in virtue of thefinite element method is presented. Both the three-dimensional transient temperature field and thetemperature evolution as a function of heat penetrating depth in the metal sheet are calculated. Thetemperature dependence of material properties was taken into account. It was shown that, after taking thetemperature dependence of the material absorbance effect into consideration, the temperature change ratealong the scanning direction and the temperature maximum were both increased.

Optimization of Off-Null Ellipsometry for Optical Biosensor Applications

The optimization of off-null ellipsometry is described with emphasis on the improvement of resolution for visualizing biomolecule layers. For optical biosensor with layer thickness below 6.5 nm, a numerical simulation for the dependence of resolution on the azimuth settings of polarizer and analyzer is presented first. For comparison, three different resolutions are given at three azimuth settings which are near null and far away from null condition, respectively. Furthermore, the square or linear approximation relationship between the intensity and the layer thickness are also given at these settings. The difference among their accuracy is up to 100 times or so. Experimental results of the biosensor sample verify the optimization.

Propagation of Explosion Wave in Hard-Soft-Hard Layered Media

The paper presents a reasonable analysis for dynamic response and failure process of a plane multi-layered media, which are subjected to a blast loading. This blast loading is induced by a cylindric explosive put on the center of top surface of the layered media. With the help of numerical simulation technique provided by LS-DYNA software, the whole process of explosion wave propagation and attenuation can be revealed. The feature of local failure around the blasting site is also discussed in some detail. Our focus will be on the explosion wave attenuation for the hard-soft-hard sandwich layers. As seen in the paper, the computational results are delivered in a feasible way by comparing with experimental data.

Wave Propagation Analysis in a Pressure-Wave-Refrigerator

Pressure wave refrigerators (PWR) refrigerate the gas through periodical expansion waves. Due to its simple structure and robustness, PWR may have many potential applications if the efficiency becomes competitive with existing alternative devices. In order to improve the efficiency, the characteristics of wave propagation in a PWR are studied by experiment, numerical simulation and theoretical analysis. Based on the experimental results and numerical simulation, a simplified model is suggested, which includes the assumptions of flux-equilibrium and conservation of the free energy. This allows the independent analysis of the operation parameters and design specifics. Furthermore, the optimum operation condition can be deduced. Some considerations to improve the PWR efficiency are also given.

Shock-Induced Dust Cloud Over a Deposit Layer of Fine Particles

The present paper investigates particle density pro les of a dust cloud induced by a normal shock wave moving at a constant speed along a at surface deposited with ne particles. In shock-fixxed coordinates, numerical simulation of ow structures of the carrier- and dispersed- phases was performed for the M = 2 case. The neness and non-uniformity of the particle size are taken into account and their effcts on the dust cloud are discussed in detail.

Oscillatory Instabilities of Two-Layer Rayleigh-Marangoni-Benard Convection

The Rayleigh-Marangoni-Benard convective instability (R-M-B instability) in the two-layer systems such as Silicone oil (10cSt)/Fluorinert (FC70) and Silicone oil (2cSt)/water liquids are studied. Both linear instability analysis and nonlinear instability analysis (2D numerical simulation) were performed to study the influence of thermocapillary force on the convective instability of the two-layer system. The results show the strong effects of thermocapillary force at the interface on the time-dependent oscillations at the onset of instability convection. The secondary instability phenomenon found in the real two-layer system of Silicone oil over water could explain the difference in the comparison of the Degen's experimental observation with the previous linear stability analysis results of Renardy et al.

The induced flow field by internal solitary wave and its action on cylindrical piles in the stratified ocean

The induced flow fields by internal solitary waves and its actions on cylindrical piles in density stratified ocean with a basic density profile and a basic velocity profile are investigated. Some results, such as the time evolution of flow fields and hydrodynamic forces on the piles are yielded both by theoretical analysis and numerical calculation for general and specific cases. Several kinds of ambient sea conditions of the South China Sea are specified for numerical simulation. Moreover, the effects of relative density difference, depth ratio and wave steepness on maximal total force and total torque are analyzed.

Torsion fracture behavior of drawn pearlitic steel wires with different heat treatments

In this paper, torsion fracture behavior of drawn pearlitic steel wires with different heat treatments was investigated. Samples with different heat treatments was investigated. Samples with different heat treatment conditions were subjected to torsion and tensile tests. The shear strain along the torsion sample after fracture was measured. Fracture surface of wires was examined by Scanning Electron Microscopy. In addition, the method of Differential Scanning Calorimetry was used to characterize the thermodynamic process in the heat treatment. A numerical simulation via finite element method on temperature field evolution for the wire during heat treatment process was performed. The results show that both strain aging and recovery process occur in the material within the temperature range between room temperature and 435 degrees C. It was shown that the ductility measured by the number of twists drops at short heating times and recovers after further heating in the lead bath of 435 degrees C. On the other hand, the strenght of the wire increases at short heating times and decreases after further heating. The microstructure inhomogeneity due to short period of heat treatment, coupled with the gradient characteristics of shear deformation during torsion results in localized shear deformation of the wire. In this situation, shear cracks nucleate between lamella and the wire breaks with low number of twists.

Flow similariy applied in convection-based tilt sensor

The environment temperature has inevitable effects on property of the convect ion-based tilt sensors. It not only redefines the application, but also prevents the improvement of the sensor performance. Numerical simulation of the fluid flow in the chamber of a sensor was performed and the influence of the environment temperature was studied in this paper. At zero tilt angle, the temperature distribution along the perpendicular line cross the heat source at various environment temperatures was presented. It was found that the flow varied dramatically at different environment temperatures, which would cause the output signal vary accordingly, even when the tilt angle was kept at a constant, because this device works by sensing the change of flow. At the same condition, we present the numerical results when the temperature difference across the heat source and the environment was kept at the same, in those results, it was found that the temperature difference at every point along the perpendicular line cross the heat source keep the same, this result confirms the similarity principle of nature convection. Second, A method of eliminating environment temperature infect on property of convect ion-based tilt sensor, which is based on the theory of flow similarity, is proposed. It was found that a thermal transistance can be piped on the circuit of heat source to compensate the temperature of the heat source. A compensative circuit was specially designed which can keep flow similarity by changing heat source temperature in order to eliminate the influence of environment temperature. The experiment results show that above 70% temperature drift can be eliminated by this compensative circuit.