A semi‐implicit scheme for large Eddy simulation of piston engine flow and combustion

A semi‐implicit scheme is presented for large eddy simulation of turbulent reactive flow and combustion in reciprocating piston engines. First, the governing equations in a deforming coordinate system are formulated to accommodate the moving piston. The numerical scheme is made up of a fourth‐order central difference for the diffusion terms in the transport equations and a fifth‐order weighted essentially nonoscillatory (WENO) scheme for the convective terms. A second‐ order Adams–Bashforth scheme is used for time integration. For higher density ratios, it is combined with a predictor–corrector scheme. The numerical scheme is explicit for time integration of the transport equations, except for the continuity equation which is used together with the momentum equation to determine the pressure field and velocity field by using a Poisson equation for the pressure correction field. The scheme is aimed at the simulation of low Mach number flows typically found in piston engines. An efficient multigrid method that can handle high grid aspect ratio is presented for solving the pressure correction equation. The numerical scheme is evaluated on two test engines, a laboratory four‐stroke engine with rectangular‐shaped engine geometry where detailed velocity measurements are available, and a modified truck engine with practical cylinder geometry where lean ethanol/air mixture is combusted under a homogeneous charge compression ignition (HCCI) condition. Copyright © 2012 John Wiley & Sons, Ltd.     

In-cylinder flows of a motored four-stroke engine with flat-crown and slightly concave-crown pistons

The temporal and spatial evolution processes of the in-cylinder flow structures and turbulence intensities in the symmetry and offset planes of a motored four-valve, four-stroke engine during the intake and compression strokes are diagnosed by using a particle image velocimeter. Two pistons of different crown shapes (flat-crown and slightly concave-crown pistons) are studied. The inception, establishment, and evolution of the tumbling vortical flow structures during the intake and compression strokes are clearly depicted. Quantitative strengths of the rotating vortical flow motions are presented by a dimensionless parameter, the tumble ratio, which can represent the mean angular velocity of the vortices in the target plane. The turbulence intensity of the in-cylinder flow is also calculated by using the measured time-varying velocity data. The results show that the flat-crown piston induces higher bulk-averaged tumble ratio and turbulence intensity than the slightly concave-crown piston does because the tumble ratio and turbulence generated by the flat-crown piston in the offset planes during the compression stroke are particularly large. The engine with the flat-crown piston also presents larger torque and power outputs and lower hydrocarbon emission than that with the slightly concave-crown piston. This might be caused by the enhanced combustion in the engine cylinder due to the stronger tumble ratio and turbulence intensity.     

基于分子动力学模拟的单晶硅冲击压缩相变研究

晶体硅具有复杂的相变机制,在相图研究中受到广泛关注,其在动载荷下的变形机制是当前研究热点。为揭示晶体硅在强动加载下的变形和相变行为特征,基于分子动力学方法,采用平板冲击加载方式,模拟研究了单晶硅在初始环境温度为300 K时分别沿[001]、[110]和[111]晶向的不同强度下的冲击压缩行为,冲击粒子速度为0.3～3.2 km/s。研究发现,随着冲击粒子速度的增加,单晶硅剪切应力在逐渐增加后由于结构相变发生急剧下降,相变阈值和相变机制均呈现各向异性。其中,沿[001]晶向冲击压缩下观察到多种固-固相变以及固-液相变,并观察到与最新文献的实验高度一致的固-液共存现象。研究结果可为动加载下晶体硅的相变研究提供纳米尺度的结果支撑。     

Holographic analysis of the deformation of a diesel engine piston

Holography has been used to measure the deformation of a diesel-engine piston heated through a low temperature rise in order to confirm the finite-element modeling of an isoformic piston being developed by Wellworthy Ltd. Previous work relating to the study indicated that a body movement had occurred on the piston and one of the main features of the current study was to measure the absolute movement at diametrically opposite positions on the piston. A single holographic plate was used for this purpose and produced the required results which, when averaged, compare favorably with finite-element results.Paper was presented at the 1988 SEM Fall Conference on Experimental Mechanics held on November 6–8 in Indianapolis, IN.     

Damper modules with adapted stiffness ratio

A mechanism for the excitation of piston rod vibrations in automotive damper modules is discussed by a simple model. An improved nonlinear model based on elasticity effects leads to good simulation results. It is shown theoretically and experimentally that the adaptation of the stiffness of the piston rod bushing to the “stiffness” of the damper force characteristic can eliminate the piston rod oscillations completely.     

高径比对砂岩单轴压缩声发射特性的影响研究

岩石尺寸效应对声发射特性的影响对于利用声发射技术监测岩石失稳破坏具有重要的实践意义.本文选择不同高径比的长石细砂岩试样开展单轴压缩声发射试验,探讨尺寸效应对岩石声发射事件数、振铃计数、峰值频率等声发射特性参数的影响.研究结果表明:根据岩石试件撞击数的增长速率,将试件破坏阶段分为突变期、平静期、爆发期三个阶段,随着岩石试...     

Comparative study of in-cylinder tumble flows in an internal combustion engine using different piston shapes—an insight using particle image velocimetry

This paper deals with the experimental investigations of the in-cylinder tumble flows in a single-cylinder engine with five different piston crown shapes at an engine speed of 1,000 rev/min., during suction and compression strokes under motoring conditions using particle image velocimetry. Two-dimensional in-cylinder tumble flow measurements and analysis are carried out in combustion space on a vertical plane passing through cylinder axis. Ensemble average velocity vectors are used to analyze the tumble flow structure. Tumble ratio and average turbulent kinetic energy are evaluated and used to characterize the tumble flows. From results, it is found that at end of compression, pentroof-offset-bowl piston shows about 41 and 103% improvement in tumble ratio and average turbulent kinetic energy respectively, compared to that of flat piston. The present study will be useful in understanding effect of piston crown shapes on nature of the in-cylinder fluid tumble flows under real engine conditions.     

A contribution to film coefficient estimation in piston cooling galleries

The need to reduce fuel consumption and exhaust emissions in internal combustion engines has been drastically increased during last years. One of the most important processes affecting these parameters is heat transfer from the in-cylinder gas to the surrounding walls, as this mechanism has a direct influence on the combustion process. Regarding the different walls (liner, cylinder head and piston surfaces), heat flow to the piston is especially important, as it is essential to avoid excessively high temperatures that could result in material damage and/or oil cracking. With this purpose different cooling strategies are used, among which the improvement of the piston cooling system by using oil galleries is preferred. In this work, the heat flow through the oil gallery in a Diesel piston was investigated on a dedicated test bench. This bench consists of a controlled heat source and a piston oil cooling system in which different test conditions were evaluated in order to obtain a correlation for the film coefficient associated with piston oil cooling. These experimental results were then incorporated into a lumped model for engine heat transfer. Finally, in order to evaluate the accuracy of this model and the effects of the correlation for oil gallery coefficient on engine heat flows, results obtained on a conventional engine test bench equipped with a Diesel engine, in which two piston temperatures had been measured, were used. The results show an improvement in piston temperature predictions when compared with those obtained using a previously reported expression for the calculation of the oil film coefficient.     

Three-dimensional photoelastic study of the load-carrying capacity/face width ratio of Wildhaber-Novikov gears for automotive applications

Wildhaber-Novikov gears are becoming more and more popular for heavy load applications. They have been tried as speed-reducing gears, in spite of heavy noise generation, in aircraft jet engines, marine engines and agricultural machinery. These circular-arc gears, though stronger than involute geras, have the disadvantage of needing larger face width. Axial face width required can be minimized only at the expense of load-carrying capacity. This is not a serious disadvantage with stationary engines. But in the case of automotive applications, this fact limits the load-carrying capacity of Wildhaber-Novikov gears, since space problem is very critical in these applications and large face widths, therefore, cannot be used for such gear-reduction units. The face width is determined by the helix angle and, therefore, a thorough investigation of the dependence of the load-carrying capacity on helix angle is necessary, if these special types of gears have to make headway into the automotive field. In this experimental investigation, three-dimensional photoelastic technique has been employed to study the load-carrying capacity/face width ratio. Three-dimensional gear models made of epoxy castings (Araldite B) were loaded in a specially built gear-loading fixture and were stress frozen. Gear cutters of Wildhaber-Novikov gears developed at the laboratory were used for the preparation of these three-dimensional gear models. The stress-frozen models were analyzed using the conventional slicing technique to study the contact and bending-stress distribution along the face width. The load-carrying capacity in terms of contact stress and bending stress has been studied for different face widths and helix angles. Twenty-deg pressure angle, 14-mm module gears of alladdendum type of Wildhaber-Novikov gears with different helix angles up to 40 deg were tested. The helical-overlap ratio used was 1.0. The conventional profile parameters were employed for the manufacture of gear cutters of end-mill type. The decrease in contact and bending strength with increase in the helix angle or decrease in the face width, as obtained from this photoelastic method of stress analysis has been compared with existing theoretical results.     

影响内燃机活塞环—缸套擦伤的因素及防擦伤的摩擦学设计

阐述了内燃机活塞环－缸套擦伤的机理。分析了磨合，缸套表面形貌，活塞环外形，活塞组件的热变形与机械变形，润滑油，材料和表面覆层等因素对活塞环擦伤的影响。指出可以采用内燃机强化系数ｐｅｖｍ作为活塞环防擦伤摩擦学设计的判别特征值，从而为从设计上防止活塞环擦伤提供了思路。     

Entrance flow in eccentric annular ducts

A control-volume-based solution of the complete set of Navier-Stokes equations for the laminar, three-dimensional developing flow in straight, eccentric, cylindrical annular ducts is described. Numerical results for velocity and pressure development, pressure defect and entrance lengths are presented for a wide range of duct parameters, i.e. relative eccentricity ? and radius ratio γ. The present results match very well with earlier numerical solutions for the limiting cases of developing flow in concentric ducts and fully developed flow in eccentric ducts. Comparison with earlier approximate results for developing flow in eccentric ducts indicates that the approximate model predicts the velocity and pressure development with an error of about 10%. However, the development length predicted by the approximate model is grossly in error. The pressure defect and development length in eccentric ducts are very high compared with their counterparts in concentric ducts. The pressure defect, development length and maximum velocity increase with the radius ratio for eccentric ducts, while the reverse is true for concentric ducts. Also, the apparent friction factor decreases as the eccentricity increases.     

A novel peristaltic micropump with low compression ratios

It is common for peristaltic micropumps to have large compression ratios. In the limit, the chamber of the pump is completely blocked by the membrane to prevent back flow. Different from this kind of pump, a micropump with small compression ratios is proposed in this study. With small oscillation amplitudes the membrane of the pump can reciprocate at high frequencies to improve its pumping flow. Both the multidimensional method and the lumped‐element method are employed for analysis. For this kind of peristaltic micropump the working fluid is allowed to flow freely in the forward and backward directions. Therefore, the operating sequences for the high‐compression ratio type of pumps are not appropriate. It is shown that the theoretical net flow rate is zero for the four‐phase and six‐phase modes of sequence and becomes negative for the three‐phase mode unless regulators, such as the nozzle/diffusers, are incorporated to rectify the flow. However, this pump becomes very attractive by reversing the operating sequence of the three‐phase mode because positive net flow is yielded. It is seen that with the reversed three‐phase mode and the nozzle/diffuser as connecting channels, the pumping effectiveness is greatly enhanced. The pumps with both two chambers and three chambers are under consideration in the study. Copyright ©2011 John Wiley & Sons, Ltd.     

Dynamic effects of the power-conversion module in a reciprocating engine

The forced response characteristics of piston, connecting rod and their assembly, henceforth called power-conversion module, is studied subjecting a forced response model of such a module to combustion characteristics in order to investigate clattering noise characteristics brought with compression ignition excitation. Existing research either focused on the piston or the connecting rod solely. As demonstrated by the modal analysis of the whole power-conversion module, it is revealed that the natural frequencies of the entire module dominate the noise-characteristics of clattering noise even when using a linear model. A subsequent parametric study applying different combustion characteristics with different pressure rise rates, but similar peak pressures on the modal-model of the power-conversion module delivered novel insights into the root cause of clattering noise characteristics. Moreover, the approach delivers an amended understanding of disturbing noises occurring in knock control systems of internal combustion engines. The reason for empirically elaborated limits of the maximum cylinder pressure rise rate to achieve smooth engine acoustics, published first in the late 1920s, was revealed.     

高焓膨胀管中压缩管最佳长度设计分析

中国空气动力研究与发展中心超高速所即将建成的高焓喷胀管采用自由活塞驱动,为了能够使该设备有效运行试验时间达到最长,需要对设备中的压缩管长度进行优化。针对压缩管中的自由活塞运动,本文通过数值求解单一组分和混合驱动气体条件下的活塞运动准一维常微分控制方程组,分析了不同试验设计参数和不同驱动气体介质对压缩管定压驱动时间的影响。采用数值解结合脊线(局部极值曲线)理论获得了定压驱动时间局部极值曲面,基于压缩管四种摩尔比例驱动气体(100%H2,85%H2+15%N2,100%He,85%He+15%Ar)和三组压缩比(λ=60,100,140)下的极值曲面,研究发现,定压驱动时间脊线高度随着压缩管的有效长度(初始活塞头部到压缩管末端膜片的距离)和内径之比L/D增加而增加。本文研究最终获得了给定参数条件下满足最长定压驱动时间要求的最佳压缩管长度。     

Effect of aspect ratio on triaxial compression of multi-sphere ellipsoid assemblies simulated using a discrete element method

Here, we present a numerical investigation of the mechanical behavior of ellipsoids under triaxial compression for a range of aspect ratios. Our simulations use a multi-sphere approach in a three-dimensional discrete element method. All assemblies were prepared at their densest condition, and triaxial compression tests were performed up to extremely large strains, until a critical state was reached. The stress–strain relationship and the void ratio–strain behavior were evaluated. We found that the stress–dilatancy relationship of ellipsoids with different aspect ratios could be expressed as a linear equation. In particular, the aspect ratio influenced the position of the critical state lines for these assemblies. Particle-scale characteristics at the critical state indicate that particles tend to be flat lying, and the obstruction of particle rotation that occurs with longer particles affects their contact mechanics. Lastly, anisotropic coefficients related to aspect ratio were investigated to probe the microscopic origins of the macroscopic behavior. A detailed analysis of geometrical and mechanical anisotropies revealed the microscopic mechanisms underlying the dependency of peak and residual strengths on aspect ratio.     

Spatio-temporal structure and cycle to cycle variations of an in-cylinder tumbling flow

The aim of this paper is to make use of PIV and high-speed PIV in a research engine of moderate tumbling ratio in order to analyze both the spatial structure of the flow and its temporal evolution during series of consecutive cycles. Appropriate analyzing tools are introduced, and four different points are addressed: (1) the chain of events driving the generation of the three-dimensional mean tumbling motion is investigated; (2) a Lagrangian analysis of the roll-up of the tumbling jet in individual cycles demonstrates a strong cycle to cycle variation during the compression phase (the rms of the position of the jet front being approximately 10% of the piston stroke); (3) focussing on the “breakdown” phase, phase invariant proper orthogonal decomposition enables us to distinguish cycles according to their structure near top dead center (TDC). We show that when the coherent energy of the flow is conserved, there is no increase in the fluctuating kinetic energy; (4) finally, the phase-averaged Reynolds stresses is decomposed into a contribution of the in-cycle coherence and the turbulence carried by the flow states. Approximately 30% of the fluctuating kinetic energy is due to cycle to cycle fluctuations in this chamber near TDC.     

FRP约束钢管混凝土柱轴力-弯矩相关关系研究

为了建立FRP(Fiber Reinforced Polymer)约束钢管混凝土压弯构件轴力-弯矩相关关系曲线,分析组合结构中各组份对柱体压弯性能的影响,建立了FRP约束钢管混凝土截面分析计算模型,考虑不同FPR-钢复合约束程度和荷载偏心率对混凝土本构的影响。与现有CFRP(Carbon Fiber Reinforced Polymer)-钢管混凝土压弯试验结果对比,验证了模型的合理性。对模型开展参数分析,结果表明,FRP约束混凝土压弯性能主要受环向FRP影响,环向FRP层数的增加在提高柱体轴压承载力的同时可有效降低其界限偏心率,提高大偏心状态下柱体的轴压比;纵向FRP层数仅会影响柱体承载力,对其界限偏心率及大偏心受压状态下的轴压比影响不显著;随着环向FRP约束程度的提高,柱体对P-Δ效应的敏感性增强。     

驱动气体的密度梯度对弹丸发射速度的影响

为进一步提升轻气炮的发射能力,提出采用梯度气体替代单一氢气或氦气作为驱动气体的方法,通过对等直径发射器进行分析,建立了弹丸在梯度气体驱动下的加速运动模型,对比了氖-氦梯度气体驱动与单一氦气驱动的发射能力差异,分析了梯度气体参数对发射性能的影响。结果表明,与单一氦气驱动相比,氖-氦梯度气体驱动能够提升0.4～1.4 km/s的发射速度或降低0.2～0.9 GPa的发射过载;气体的密度和活塞的运动速度对发射速度和过载的影响最大,气体压力和多方气体指数的影响次之;梯度气体中,高密度气体应选择多方气体指数和密度较高的气体（如氖气、氩气等）;梯度气体界面位置（高密度气体占比）对发射速度的影响不大,但高密度气体占比少有利于降低弹底压力。     

Heat transfer during transient compression: Measurements and simulations

Experiments have been performed to assess the utility of unsteady one-dimensional heat conduction modelling for the calculation of heat losses during a free piston compression process. Heat transfer measurements have been obtained within a gun tunnel barrel using surface junction thermocouple instrumentation. The gun tunnel was operated with a relatively heavy piston such that the shock waves induced by the piston motion were weak. Peak heat transfer values are estimated reasonably well by the unsteady one-dimensional model. However, overall quantitative agreement between the measurements and calculations has not been achieved at this stage, principally because the development of turbulent heat transport was not properly modelled. Received 21 September 2001 / Accepted 11 March 2002 – Published online 11 June 2002     

气液活塞式脉冲液体射流泵装置稳定性的理论研究

气液活赛式脉冲液体射流泵装置用于输送有毒、高温或放射性液体,装置在运行时,气液活塞筒内的气液交界面必须稳定在活塞筒内,形成稳定的脉冲运动,才能保证装置正常的稳定运行,本文在分析影响气液活塞式脉冲液体射流泵装置稳定性的主要因素基础上,运用流体力学的基本理论,导出了装置稳定性的基本方程组以及简化方程组,通过试验对上述理论进行了验证。