Simulation studies on the performance of thermoacoustic prime movers and refrigerator

The prime movers and refrigerators based on thermoacoustics have gained considerable importance toward practical applications in view of the absence of moving components, reasonable efficiency, use of environmental friendly working fluids, etc. Devices such as twin Standing Wave ThermoAcoustic Prime Mover (SWTAPM), Traveling Wave ThermoAcoustic Prime Mover (TWTAPM) and thermoacoustically driven Standing Wave ThermoAcoustic Refrigerator (SWTAR) have been studied by researchers. The numerical modeling and simulation play a vital role in their development. In our efforts to build the above thermoacoustic systems, we have carried out numerical analysis using the procedures of CFD on the above systems. The results of the analysis are compared with those of DeltaEC (freeware from LANL, USA) simulations and the experimental results wherever possible. For the CFD analysis commercial code Fluent 6.3.26 has been used along with the necessary boundary conditions for different working fluids at various average pressures. The results of simulation indicate that choice of the working fluid and the average pressure are critical to the performance of the above thermoacoustic devices. Also it is observed that the predictions through the CFD analysis are closer to the experimental results in most cases, compared to those of DeltaEC simulations. (C) 2015 Elsevier Ltd. All rights reserved.

Implicit scheme for meshless compressible Euler solver

In this paper, an implicit scheme is presented for a meshless compressible Euler solver based on the Least Square Kinetic Upwind Method (LSKUM). The Jameson and Yoon's split flux Jacobians formulation is very popular in finite volume methodology, which leads to a scalar diagonal dominant matrix for an efficient implicit procedure (Jameson & Yoon, 1987). However, this approach leads to a block diagonal matrix when applied to the LSKUM meshless method. The above split flux Jacobian formulation, along with a matrix-free approach, has been adopted to obtain a diagonally dominant, robust and cheap implicit time integration scheme. The efficacy of the scheme is demonstrated by computing 2D flow past a NACA 0012 airfoil under subsonic, transonic and supersonic flow conditions. The results obtained are compared with available experiments and other reliable computational fluid dynamics (CFD) results. The present implicit formulation shows good convergence acceleration over the RK4 explicit procedure. Further, the accuracy and robustness of the scheme in 3D is demonstrated by computing the flow past an ONERA M6 wing and a clipped delta wing with aileron deflection. The computed results show good agreement with wind tunnel experiments and other CFD computations.

Thermal Modeling of Organic Light-Emitting Diode Display Panels

Power densities required to operate active-matrix organic-light-emitting diode (AMOLED) based displays for high luminance applications, lead to temperature rise due to self heating. Temperature rise leads to significant degradation and consequent reduction in life time. In this work numerical techniques based computational fluid dynamics (CFD) is used to determine the temperature rise and its distribution for an AMOLED based display for a given power density and size. Passive cooling option in form of protruded rectangular fins is implemented to reduce the display temperature.

Numerical Investigations for Leakage and Windage Heating in Straight-Through Labyrinth Seals

The ability to quantify leakage flow and windage heating for labyrinth seals with honeycomb lands is critical in understanding gas turbine engine system performance and predicting its component life. Variety of labyrinth seal configurations (number of teeth, stepped or straight, honeycomb cell size) are in use in gas turbines, and for each configuration, there are many geometric factors that can impact a seal's leakage and windage characteristics. This paper describes the development of a numerical methodology aimed at studying the effect of honeycomb lands on leakage and windage heating. Specifically, a three-dimensional computational fluid dynamics (CFD) model is developed utilizing commercial finite volume-based software incorporating the renormalization group (RNG) k-epsilon turbulence model with modified Schmidt number. The modified turbulence model is benchmarked and fine-tuned based on several experiments. Using this model, a broad parametric study is conducted by varying honeycomb cell size, pressure ratio (PR), and radial clearance for a four-tooth straight-through labyrinth seal. The results show good agreement with available experimental data. They further indicate that larger honeycomb cells predict higher seal leakage and windage heating at tighter clearances compared to smaller honeycomb cells and smooth lands. However, at open seal clearances larger honeycomb cells have lower leakage compared to smaller honeycomb cells.

COUPLED OPTICAL-THERMAL-FLUID MODELING OF A DIRECTLY HEATED TUBULAR SOLAR RECEIVER FOR SUPERCRITICAL CO2 BRAYTON CYCLE

Recent studies have evaluated closed-loop supercritical carbon dioxide (s-CO2) Brayton cycles to be a higher energy density system in comparison to conventional superheated steam Rankine systems. At turbine inlet conditions of 923K and 25 MPa, high thermal efficiency (similar to 50%) can be achieved. Achieving these high efficiencies will make concentrating solar power (CSP) technologies a competitive alternative to current power generation methods. To incorporate a s-CO2 Brayton power cycle in a solar power tower system, the development of a solar receiver capable of providing an outlet temperature of 923 K (at 25 MPa) is necessary. The s-CO2 will need to increase in temperature by similar to 200 K as it passes through the solar receiver to satisfy the temperature requirements of a s-CO2 Brayton cycle with recuperation and recompression. In this study, an optical-thermal-fluid model was developed to design and evaluate a tubular receiver that will receive a heat input similar to 2 MWth from a heliostat field. The ray-tracing tool SolTrace was used to obtain the heat-flux distribution on the surfaces of the receiver. Computational fluid dynamics (CFD) modeling using the Discrete Ordinates (DO) radiation model was used to predict the temperature distribution and the resulting receiver efficiency. The effect of flow parameters, receiver geometry and radiation absorption by s-CO2 were studied. The receiver surface temperatures were found to be within the safe operational limit while exhibiting a receiver efficiency of similar to 85%.

轴对称近似等熵压缩流场的乘波前体优化设计

以升阻比为优化目标，在来流马赫数Ma＝2～4及飞行高度H＝20km～24km条件下，进行了轴对称近似等熵压缩流场的乘波前体优化设计，通过CFD验证M_∞＝4优化乘波体的气动特性，并研究了Ma＝3优化乘波前体在非设计条件下的气动特性。结果表明：近似等熵压缩下表面的乘波前体在设计条件下具有良好的气流压缩效果，可满足机体／发动机一体化设计的需要；乘波前体升阻比在1.5～1.9之间，纵向压心位置靠后；非设计条件下，压缩波不聚焦，小于设计马赫数升阻比时降低，大于设计马赫数时升阻比略大。

轴对称射流气动声场的数值模拟

构造了Kirchhoff积分和CFD计算相结合的算法，采用高精度差分格式对不同喷口马赫数的轴对称射流进行数值模拟，并以此作为近场声源，运用Kirchhoff积分方法研究远场气动噪声。数值结果表明，远场噪声具有方向性，噪声声压在离开对称轴20℃处达到最大值。随着传播距离增大，噪声方向性逐渐减弱。

COIL亚声速段横向喷流混合流场数值分析

COIL是一个气体动力学、化学反应动力学以及光学相互耦合的复杂过程.作为高总压COIL系统研究的第一步,利用三维CFD技术对传统的COIL亚声速段横向射流混合过程进行了数值分析,讨论了包括压力梯度驱动项的分子扩散机制,得到了横向射流的较精细的结构,如马蹄形射流界面、逆旋涡对以及射流剪切层.结果表明,压力梯度驱动项对重、轻组分的扩散作用相反,重组分沿压力梯度正向扩散,轻组分逆压力梯度方向扩散.在射流穿透不足的情况下,仍然在喷管出口得到了约为0.01 cm-1的小信号增益系数.

The hemodynamic effects of in-tandem carotid artery stenosis: implications for carotid endarterectomy.

OBJECTIVES: It remains controversial whether patients with severe disease of the internal carotid artery and a coexisting stenotic lesion downstream would benefit from a carotid endarterectomy (CEA) of the proximal lesion. The aim of this study was to simulate the hemodynamic and wall shear effects of in-tandem internal carotid artery stenosis using a computational fluid dynamic (CFD) idealized model to give insight into the possible consequences of CEA on these lesions. METHODS: A CFD model of steady viscous flow in a rigid tube with two asymmetric stenoses was introduced to simulate blood flow in arteries with multiple constrictions. The effect of varying the distance between the two stenoses, and the severity of the upstream stenosis on the pressure and wall shear stress (WSS) distributions on the second plaque, was investigated. The influence of the relative positions of the two stenoses was also assessed. RESULTS: The distance between the plaques was found to have minimal influence on the overall hemodynamic effect except for the presence of a zone of low WSS (range -20 to 30 dyne/cm2) adjacent to both lesions when the two stenoses were sufficiently close (<4 times the arterial diameter). The upstream stenosis was protective if it was larger than the downstream stenosis. The relative positions of the stenoses were found to influence the WSS but not the pressure distribution. CONCLUSIONS: The geometry and positions of the lesions need to be considered when considering the hemodynamic effects of an in-tandem stenosis. Low WSS is thought to cause endothelial dysfunction and initiate atheroma formation. The fact that there was a flow recirculation zone with low WSS in between the two stenoses may demonstrate how two closely positioned plaques may merge into one larger lesion. Decision making for CEA may need to take into account the hemodynamic situation when an in-tandem stenosis is found. CFD may aid in the risk stratification of patients with this problem.

Three-Dimensional Propagation of the Transmitted Shock Wave in a Square Cross-Sectional Chamber

An investigation into the three-dimensional propagation of the transmitted shock wave in a square cross-section chamber was described in this paper, and the work was carried out numerically by solving the Euler equations with a dispersion-controlled scheme. Computational images were constructed from the density distribution of the transmitted shock wave discharging from the open end of the square shock tube and compared directly with holographic interferograms available for CFD validation. Two cases of the transmitted shock wave propagating at different Mach numbers in the same geometry were simulated. A special shock reflection system near the corner of the square cross-section chamber was observed, consisting of four shock waves: the transmitted shock wave, two reflection shock waves and a Mach stem. A contact surface may appear in the four-shock system when the transmitted shock wave becomes stronger. Both the secondary shock wave and the primary vortex loop are three-dimensional in the present case due to the non-uniform flow expansion behind the transmitted shock.

Pulverized-coal combustion in Wall-Protecting-Jets combustor (WPJC)

A new type of pulverized-coal combustor, called "Wall-Protecting-Jets Combustor" (hereafter, WPJC has been proposed, designed and studied with both CFD (Computational Fluid Dynamics) and experimental methods. The WPJC is based on a novel concept in which all inlet jets are along the combustor wall. Pilot combustion experiments were conducted to investigate the combustion performance of WPJC. Two-phase flows and pulverized-coal combustion were simulated to study the mechanism of),WPJC using the commercial software FLUENT. The results show that the WPJC has many remarkable advantages: wall-protection by the cold jets without the use of refractory materials; low-temperature and three-stage combustion with low NOx emission; negligible ash/slag-deposition; multiple functions with convenient switching between them; effective adjustment of the combustion intensity and the ignition position.

一种易于并行求解Euler方程组的分区技术

在具有复杂边界的计算区域内,求解偏微分方程组时,经常需要分区和并行计算,分区方法直接关系到数值计算的并行化程度,本文在应用时间算子分裂方法求解Euler方程组的过程中,提出了一种非常容易实现并行化计算的分区技术.

Characterisation of cell adhesion to substrate materials and the resistance to enzymatic and mechanical cell-removal

Cell-implant adhesive strength is important for prostheses. In this paper, an investigation is described into the adhesion of bovine chondrocytes to Ti6Al4V-based substrates with different surface roughnesses and compositions. Cells were cultured for 2 or 5 days, to promote adhesion. The ease of cell removal was characterised, using both biochemical (trypsin) and mechanical (accelerated buoyancy and liquid flow) methods. Computational fluid dynamics (CFD) modelling has been used to estimate the shear forces applied to the cells by the liquid flow. A comparison is presented between the ease of cell detachment indicated using these methods, for the three surfaces investigated. © 2008 Materials Research Society.

Large eddy simulations for fan-OGV broadband noise prediction

The work in this paper forms part of a project on the use of large eddy simulation (LES) for broadband rotor-stator interaction noise prediction. Here we focus on LES of the flow field near a fan blade trailing edge. The first part of the paper aims to evaluate LES suitability for predicting the near-field velocity field for a blunt NACA-0012 airfoil at moderate Reynolds numbers (2× 10 5 and 4× 10 5). Preliminary computations of turbulent mean and root-mean-square velocities, as well as energy spectra at the trailing edge, are compared with those from a recent experiment.1 The second part of the paper describes preliminary progress on an LES calculation of the fan wakes on a fan rig. 2 The CFD code uses a mixed element unstructured mesh with a median dual control volume. A wall-adapting local eddy-viscosity sub-grid scale model is employed. A very small amount of numerical dissipation is added in the numerical scheme to keep the compressible solver stable. Further results for the fan turbulentmean and RMS velocity, and especially the aeroacoustics field will be presented at a later stage. Copyright © 2008 by Qinling LI, Nigel Peake & Mark Savill.

用FLUENT软件计算化学氧碘激光流场

用计算流体动力学(CFD)软件FLUENT,数值研究化学氧碘激光喷管形线,碘喷孔位置,副流入口压力,O2(1△)初始产额,水蒸汽含量以及稀释气体对平均小信号增益系数沿流动方向分布的影响;并对大连化学物理研究所超音速化学氧碘激光的两次实验进行数值模拟,计算结果显示平均小信号增益系数与实验测试的摩尔功率趋势一致.