Experimental and numerical studies of micro PEM fuel cell

A single micro proton exchange membrane fuel cell (PEMFC) has been produced using Micro-electromechanical systems (MEMS) technology with the active area of 2.5 cm 2 and channel depth of about 500 μ m.A theoretical analysis is performed in this study for a novel MEMS-based design of a micro PEMFC.The model consists of the conservation equations of mass,momentum,species and electric current in a fully integrated finite-volume solver using the CFD-ACE+ commercial code.The polarization curves of simulation are well correlated with experimental data.Three-dimensional simulations are carried out to treat prediction and analysis of micro PEMFC temperature,current density and water distributions in two different fuel flow rates (15 cm 3 /min and 40 cm 3 /min).Simulation results show that temperature distribution within the micro PEMFC is affected by water distribution in the membrane and indicate that low and uniform temperature distribution in the membrane at low fuel flow rates leads to increased membrane water distribution and obtains superior micro PEMFC current density distribution under 0.4 V operating voltage.Model predictions are well within those known for experimental mechanism phenomena.     

Experimental and numerical studies of weak blast waves in air

Effects of viscosity and vibrational nonequilibrium on the profile of a weak, spherical N-wave in air are experimentally and numerically studied. Weak blast waves were generated, in a quiescent air dome, by spark discharges and exploding wires and observed by high frequency response microphones over 40 meters. Some similarity relationships were obtained from the blast wave experiments. For observed N-waves having less than 100 Pa peak overpressure, the peak overpressure p _f and the duration of the positive phaset _d+ are found to vary with the radial distance from the sourcer as p _f r ^–1.38 andt _d + r ^0.19, whilst the rise time of the blast wave t _f linearly increases with distance. Similar trends were also found for the negative phase of the blast wave. Numerical simulations were carried out to compare with the blast wave data. The Navier-Stokes equations for spherical symmetric flows were solved by coupling with a relaxation equation for vibrational excitation of oxygen using the random choice method (RCM) adapted to supercomputing with an operator splitting technique. The resultant N-wave profiles are in good agreement with the experimental results. The numerical results clearly indicate that the wave-easing process due to the dispersive effect of vibrational relaxation plays a dominant role in determining the rise time of the N-wave.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1990.     

Experimental and numerical studies of rotational relaxation behind a strong shock wave in air

This paper describes the experimental and numerical investigations of unknown characteristics of the rotational nonequilibrium phenomena behind a strong shock wave in air. Experiments were carried out using a piston-driven shock tube with helium as driving gas and air as driven (test) gas, operated as a two-stage shock tube. In the experiments, emission spectra of NO were measured to evaluate the rotational temperature behind a strong shock wave. The numerical calculations use the computational code for the thermal and chemical nonequilibrium flow behind a strong shock wave developed by the present author's group, where 11 chemical species (N, O, NO, N, O, N, O, NO, N, O, e) and 48 chemical reactions of high-temperature air are considered. The thermal nonequilibrium is expressed by introducing an 8 temperature model composed of translational temperature, rotational and vibrational temperatures for N, O, NO, and electron temperature. The coupling of a rotation, vibration and dissociation (CRVD) model was incorporated to take sufficiently into account the rotational nonequilibrium. The calculations were conducted for the same conditions as the experimental ones. From the calculated flow properties, emission spectra were re-constructed using the code for computing spectra of high temperature air “SPRADIAN”. Furthermore, rotational and vibrational temperatures of NO (0,1) were determined from a curve fitting method and compared with the computed results. Received 12 September 2001 / Accepted 18 February 2002     

无氧铜层裂的实验与数值研究

利用一级气体炮对国产的无氧铜（OFHC）进行了平板撞击致层裂实验,3 mm厚的OFHC飞片撞击6 mm厚的OFHC靶,采用锰铜应力计记录试样/有机玻璃（PMMA）界面附近的应力历史,获得试样发生层裂的信息。软回收试样,观测回收试样的层裂片。采用基于空穴聚集的层裂模型,数值模拟这些平板撞击致层裂实验。数值模拟的应力剖面以及试样层裂片厚度与实验结果基本一致。此外,对于国外相关的OFHC层裂实验,采用基于空穴聚集的层裂模型也作了相应的数值模拟,并进行了比较。     

粉尘空气混合物层流燃烧速度的实验测定

给出了计算粉尘层流火焰速度的直接方法,此法简便易行,在粉尘浓度较低时,计算精度较高。实验结果表明:粉尘浓度对火焰传播速度和燃烧速度有很大的影响,粉尘浓度过大时,粒子运动轨迹就难观测,用直接法计算层流火焰速度的误差增大;管径大小也对燃烧速度有很大影响;小管径中的所得值比大管径中的所得值约低８％。     

Microgravity experiments and numerical studies on ethanol/air spray flames

Spray flames are known to exhibit amazing features in comparison with single-phase flames. The weightless situation offers the conditions in which the spray characteristics can be well controlled before and during combustion. The article reports on a joint experimental/numerical work that concerns ethanol/air spray flames observed in a spherical chamber using the condensation technique of expansion cooling (based on the Wilson cloud chamber principle), under microgravity.We describe the experimental set-up and give details on the creation of a homogeneous and nearly monosized aerosol. Different optical diagnostics are employed successfully to measure the relevant parameters of two-phase combustion. A classical shadowgraphy system is used to track the flame speed propagation and allow us to observe the flame front instability. The complete characterization of the aerosol is performed with a laser diffraction particle size analyser by measuring the droplet diameter and the droplet density number, just before ignition. A laser tomography device allows us to measure the temporal evolution of the droplet displacement during flame propagation, as well as to identify the presence of droplets in the burnt gases. The numerical modelling is briefly recalled. In particular, spray-flame propagation is schematized by the combustion spread in a 2-D lattice of fuel droplets surrounded by an initial gaseous mixture of fuel vapour and air.In its spherical expansion, the spray flame presents a corrugated front pattern, while the equivalent single-phase flame does not. From a numerical point of view, the same phenomena of wrinkles are also observed in the simulations. The front pattern pointed out by the numerical approach is identified as of Darrieus–Landau (DL) type. The droplets are found to trigger the instability. Then, we quantitatively compare experimental data with numerical predictions on spray-flame speed. The experimental results show that the spray-flame speed is of the same order of magnitude as that of the single-phase premixed flame. On the other hand, the numerical results exhibit the role played by the droplet radius in spray-flame propagation, and retrieve the experiments only when the droplets are small enough and when the Darrieus–Landau instability is triggered. A final discussion is developed to interpret the various patterns experimentally observed for the spray-flame front.     

Experimental studies on the yield behavior of ductile and brittle aluminum foams

Specimen size effects are a major cause of the unreliability of foam models in finite element codes. Here, the modified Arcan apparatus is used to investigate the biaxial yielding of ductile and brittle Al foams. This apparatus subjects a central section of a “butterfly-shaped” specimen to a uniform state of plane stress. The stresses have local maxima at the central section, thus ensuring that yielding occurs there. A yield envelope, which directly relates to the crushing process, can then be determined. Size effects are introduced when using conventional methods such as tri-axial or plate-shear tests. In such tests, averages of stress and strain are measured. These measures do not represent the actual yield event, because foam's internal structure is inhomogeneous and so is the deformation field. Strain localization and failure can occur at any weak layer of cells in the bulk. In this study, we have performed a series of biaxial tests on isotropic Alporas and anisotropic Hydro closed-cell Al foams of approximately equal densities. Alporas failed locally by a ductile phenomenon of progressive crushing of cells. It also possessed uniaxial strength asymmetry. Hydro specimens parallel and perpendicular to ‘foam rise’ were investigated. The Hydro foam developed a local, characteristic brittle crack at loads in the vicinity of the yield point. Phenomenological yield surfaces, which incorporate these features are obtained for the foams, and show dependence on both the deviatoric and hydrostatic stresses. We also provide expressions for the shear and hydrostatic strengths in terms of the uniaxial strengths. Finally, the size-independence of the yield surface is verified using the uniaxial compression of tapered specimens.     

Experimental and numerical studies of creasing of paperboard

A laboratory creasing device to capture the most important properties of a commercial rotary creasing tool was designed. Finite element analysis of the creasing of a multiply paperboard in the laboratory crease device was presented. The multiply paperboard was modeled as a multilayered structure with cohesive softening interface model connecting the paperboard plies. The paperboard plies were modeled by an anisotropic elastic–plastic material model. The purpose of the analysis of the laboratory creasing device was to present material models that represent paperboard, and to investigate how well the analysis captured the multiply paperboard behavior during laboratory creasing. And to increase the understanding of what multiply paperboard properties that influence the laboratory crease operation. The result of the simulations showed very good correlations with the experimental obtained results. The results indicated that the paperboard properties that have the most influence is the out-of-plane shear, out-of-plane compression and the friction between the laboratory creasing device and the paperboard.     

Experimental and numerical studies of impact on filament-wound composite cylinder

This study focused on the impact behavior of carbon-fiber-wrapped composite cylinders subjected to impact from flat-ended,hemispherical-nosed and conical-nosed impactors.Damage morphologies of the cylinders and mechanisms of the damage were analyzed.Change laws of the maximum impact forces,durations of impact processes and energies absorbed by the cylinders after impact with different impactors and impact energies were obtained.A finite element model was developed and the simulation results were in reasonable agreement with the tests.Finally,taking the flat-ended impactor as an example,stress distributions of the cylinders under pressurization and impact were discussed.     

钙质砂场地爆炸成坑实验与数值模拟研究

为了建立钙质砂场地爆炸成坑效应的计算方法,首先在开挖出的钙质砂模型场地开展了不同当量、不同埋深的野外爆炸实验,然后基于有限元与光滑粒子流耦合算法建立了适用于钙质砂爆炸成坑计算的数值模型,并分析了炸药形状和土体参数对爆坑形态的影响,最后建立了适用于钙质砂场地中的爆坑计算公式。结果表明：埋置爆下,钙质砂场地爆坑尺寸大于硅质砂土中的爆坑尺寸;光滑粒子流算法能较好地揭示钙质砂场地中爆坑轮廓的形成机理;炸药形状和土体密实度等参数对于钙质砂爆坑形态具有不同程度的影响,拟合得到的钙质砂场地接触爆和埋置爆抛掷型爆坑尺寸计算公式,可较好地预测不同爆炸当量作用下的爆炸成坑尺寸。

     

Experimental and numerical studies of shock wave focusing in water

The focusing of plane shock waves and spherical blast waves in water by parabolic and ellipsoidal reflectors is studied experimentally and numerically. The numerical computations are performed by a second order method based on an extension of Godunov's method, which is called piecewise-linear method. Tests of this method and a comparison with pressure measurements show a strong dependence of the maximum pressure attained in the focal region on the mesh size. Outside this focal region, however, good agreement between the pressure measurements and the result of the numerical computation is found. Furthermore it is shown by numerical computation, that non-linear effects in the focusing process may be partly avoided by an optimized reflector contour which depends on the strength of the incident shock wave. This modification provides that the maximum pressure is attained at the geometric focus and not at a location between reflector surface and geometrical focus.     

泡沫铝合金动态力学性能实验研究

利用分离式霍布金森压杆(SHPB)实验技术和MTS材料实验机对两组不同孔径、不同密度的开孔泡沫铝合金进行了准静态和动态压缩实验研究。实验结果表明：泡沫铝合金的静态和动态变形过程均具有泡沫材料变形的三个阶段特征。开孔泡沫铝合金的变形是均匀变化过程,并不出现局部的变形带。与相对密度对力学性能的影响相比,孔径大小的影响可以忽略不计。在考察的应变率范围内,屈服应力对应变率并不很敏感。     

Experimental and numerical studies of underwater shock wave attenuation

The attenuation of an underwater shock wave by a thin porous layer is studied both experimentally and numerically. The shock waves are generated by exploding 10 mg silver azide pellets and the pressures at different distances from the explosion center are measured. Measurements are also carried out with a gauze layer placed between the explosion source and the pressure gauge. The results with and without the gauze layer are compared evaluating the shock wave attenuation. Numerical simulations of the phenomenon are also carried out for a simple wave attenuation model. The results are compared with the experimental data. Despite the simple mathematical model of wave attenuation, the agreement between the experimental and numerical results is reasonable.Received: 22 October 2002, Accepted: 17 June 2003, Published online: 5 August 2003PACS: 47.11.+j, 47.40.Nm, 47.55.Mh     

Experimental and numerical studies of the pressure field generated by DDT events

Results of experimental and numerical study of maximum pressures and impulses developed in closed vessel as a function of DDT location are presented. Experiments were carried out in the detonation tube with stoichiometric hydrogen-oxygen mixture diluted with nitrogen. The effect of two variables was investigated: distance between tube closed end and DDT origin, and the ratio of burned and unburned mixture volume prior to the DDT moment. It was found that impulse value is proportional to the total length of tube and does not depend significantly on DDT location. Results of numerical simulation are in good agreement with experimental data. Received March 10, 1995 / Accepted March 12, 1995     

小尺寸管道内二氧化碳抑制甲烷爆炸效果的实验及数值模拟

为了有效防治矿井瓦斯爆炸事故, 以瓦斯的主要成分甲烷作为模拟气体, 运用自主设计改装的XKWB-S型小尺寸石英玻璃管道实验系统, 结合高速摄影仪, 并采用FLACS数值模拟软件, 研究惰性气体抑爆条件下甲烷燃烧爆炸特性, 进行体积分数为6%~27%的CO₂抑制体积分数为9%CH₄爆炸的实验及数值模拟, 结果表明:各组分混合气体在爆炸传播过程中, 爆炸压力、火焰锋面速度和气体运动速度均呈现一定程度的波动, 且压力和速度没有同时达到最大值; CO₂的加入有效抑制了甲烷/空气反应, 且添加CO₂体积分数越大, 抑爆效果越明显, 模拟结果与实验结果基本吻合。     

Experimental and numerical studies of thermo-hydrous transfers in concrete exposed to high temperature

The compressive strength of concrete can be as high as 80 MPa at 28 days. High strength concrete (HSC) can be obtained by decreasing porosity and lowering permeability. Concrete, especially HSC, performs poorly when subjected to fire. This is attributed to high thermal stresses and water vapour pressure. High thermal gradient induces high thermo-mechanical stresses in the concrete system. Low permeability prevents water from escaping and induces high water vapour pressure causing cracking and spalling. The aim of this study is both experimentally and numerically study the coupled heat and mass transfers in concrete exposed to elevated temperature. Five concrete mixtures with various cement contents and water cement ratios of a constant aggregate content were studied before and after heating–cooling cycles. The concrete cylindrical specimens were subjected to several tests: compression and splitting tensile tests, measurement of modulus of elasticity, heating–cooling cycles, thermal field and mass loss during the heating–cooling cycles, and permeability tests. Comparisons between the numerical and experimental results on the thermo-hydrous behaviour were reported. Parametric analyses were carried out in order to underline main parameters involved in concrete behaviour at high temperature. The numerical and experimental results included thermal gradient, water vapour pressure, relative humidity, concrete mass losses due to dehydration, and water content for concrete elements heated from 20 to 600°C. The results show the degrees of damage due to the concrete chemical transformations at high temperature.     

Experimental and numerical investigations on slot air jet impingement cooling of a heated cylindrical concave surface

Experimental and numerical studies have been carried out for slot air jet impingement on a heated concave surface of a partially opened-top horizontal cylinder of length L = 20 cm. The slot jet is situated at the symmetry line of the partially opened-top cylinder along the gravity vector and impinges to the bottom of the cylinder which is designated as θ = 0°. The width of the opening at the top of the horizontal cylinder is W = 3 cm which corresponds to a circumferential angle Δθ = 50.8°. The experiments are performed by a Mach–Zehnder interferometer which enables to measure the local convection heat transfer coefficient. Also, a finite volume method based on the SIMPLE algorithm and non-orthogonal grid discretization scheme is used to solve the continuity, momentum, and energy equations. The Poisson equations are solved for (x, y) to find the grid points which are distributed in a non-uniform manner with higher concentration close to the solid regions. The effects of jet Reynolds number (Re _j) in the range from 190 to 1,600 and the ratio of spacing between nozzle and cylinder surface to the jet width from H = 1.5 to H = 10.7 on the local and average Nusselt numbers are examined. It is observed that maximum Nusselt number occurs at the stagnation point at (θ = 0°) and the local heat transfer coefficient decreases on the circumferential surface of the cylinder with increase of θ as a result of thermal boundary layer thickness growth. Also results show that the local and average heat transfer coefficients are raised by increasing the jet Reynolds number and by decreasing the nozzle-to-surface spacing.     

装甲钢/UHPC复合靶体抗侵彻性能试验与数值模拟研究

装甲钢/超高性能混凝土（UHPC）复合防护结构在重点工程中抵抗弹体的高速侵彻作用具有广泛的应用前景。为评估该复合结构的抗侵彻性能,对两种复合靶体开展侵彻试验与数值模拟研究。首先,开展了12发30 mm口径30CrMnSiNi2A弹体372~646 m/s速度侵彻复合靶试验。随后通过一系列静动态力学性能试验标定装甲钢材料的本构模型参数,并建立三维有限元模型对上述试验开展数值模拟分析。通过对比试验和数值模拟得到的弹体侵彻深度、残余弹体长度和装甲钢板的失效模式,验证了装甲钢本构模型参数的可靠性。进一步基于弹道效益系数对复合靶抗侵彻性能进行了定量评估。最后,确定了不同装甲钢板厚度复合靶体的临界贯穿速度,并对弹体侵彻复合靶的弹、靶失效模式进行了讨论。