Transonic Aeroelastic Numerical Simulation in Aeronautical Engineering

A lower-upper symmetric Gauss-Seidel (LU-SGS) subiteration scheme is constructed for time-marching of the fluid equations. The Harten-Lax-van Leer-Einfeldt-Wada (HLLEW) scheme is used for the spatial discretization. The same subiteration formulation is applied directly to the structural equations of motion in generalized coordinates. Through subiteration between the fluid and structural equations, a fully implicit aeroelastic solver is obtained for the numerical simulation of fluid/structure interaction. To improve the ability for application to complex configurations, a multiblock grid is used for the flow field calculation and transfinite interpolation (TFI) is employed for the adaptive moving grid deformation. The infinite plate spline (IPS) and the principal of virtual work are utilized for the data transformation between the fluid and structure. The developed code was first validated through the comparison of experimental and computational results for the AGARD 445.6 standard aeroelastic wing. Then, the flutter character of a tail wing with control surface was analyzed. Finally, flutter boundaries of a complex aircraft configuration were predicted.

Integration property of a wedge-elliptic cone waverider with scramjet engine

The hypersonic waverider forebody is designed in this paper. For the present waverider, the undersurface is carved out as a stream surface of a hypersonic inviscid flow field around wedge-elliptic cone, and the upper surface is assumed to be a freestream surface. A finite-volume code is used to generate the three-dimensional flow field. The leading edge is determined by satisfying the condition that the lip is situated at the intersection line of shocks.

Development of a Cell Size Relaxed Scheme for the Direct Simulation Monte Carlo Method

The conventional direct simulation Monte Carlo (DSMC) method has a strong restriction on the cell size because simulated particles are selected randomly within the cell for collisions. Cells with size larger than the molecular mean free path are generally not allowed in correct DSMC simulations. However, the cell-size induced numerical error can be controlled if the gradients of flow properties are properly involved during collisions. In this study, a large cell DSMC scheme is proposed to relax the cell size restriction. The scheme is applied to simulate several test problems and promising results are obtained even when the cell size is greater than 10 mean free paths of gas molecules. However, it is still necessary, of course, that the cell size be small with respect to the flow field structures that must be resolved.

Information Preservation Modeling of Rayleigh-Benard Transition from Thermal Conduction to Convection

Onset and evolution of the Rayleigh-Benard (R-B) convection are investigated using the Information Preservation (IP) method. The information velocity and temperature are updated using the Octant Flux Splitting (OFS) model developed by Masters & Ye based on the Maxwell transport equation suggested by Sun & Boyd. Statistical noise inherent in particle approaches such as the direct simulation Monte Carlo (DSMC) method is effectively reduced by the IP method, and therefore the evolutions from an initial quiescent fluid to a final steady state are shown clearly. An interesting phenomenon is observed: when the Rayleigh number (Ra) exceeds its critical value, there exists an obvious incubation stage. During the incubation stage, the vortex structure clearly appears and evolves, whereas the Nusselt number (Nu) of the lower plate is close to unity. After the incubation stage, the vortex velocity and Nu rapidly increase, and the flow field quickly reaches a steady, convective state. A relation of Nu to Ra given by IP agrees with those given by DSMC, the classical theory and experimental data.

Subgrid scale fluid velocity timescales seen by inertial particles in large-eddy simulation of particle-laden turbulence

Large-eddy simulation (LES) has emerged as a promising tool for simulating turbulent flows in general and, in recent years,has also been applied to the particle-laden turbulence with some success (Kassinos et al., 2007). The motion of inertial particles is much more complicated than fluid elements, and therefore, LES of turbulent flow laden with inertial particles encounters new challenges. In the conventional LES, only large-scale eddies are explicitly resolved and the effects of unresolved, small or subgrid scale (SGS) eddies on the large-scale eddies are modeled. The SGS turbulent flow field is not available. The effects of SGS turbulent velocity field on particle motion have been studied by Wang and Squires (1996), Armenio et al. (1999), Yamamoto et al. (2001), Shotorban and Mashayek (2006a,b), Fede and Simonin (2006), Berrouk et al. (2007), Bini and Jones (2008), and Pozorski and Apte (2009), amongst others. One contemporary method to include the effects of SGS eddies on inertial particle motions is to introduce a stochastic differential equation (SDE), that is, a Langevin stochastic equation to model the SGS fluid velocity seen by inertial particles (Fede et al., 2006; Shotorban and Mashayek, 2006a; Shotorban and Mashayek, 2006b; Berrouk et al., 2007; Bini and Jones, 2008; Pozorski and Apte, 2009).However, the accuracy of such a Langevin equation model depends primarily on the prescription of the SGS fluid velocity autocorrelation time seen by an inertial particle or the inertial particle–SGS eddy interaction timescale (denoted by $\delt T_{Lp}$ and a second model constant in the diffusion term which controls the intensity of the random force received by an inertial particle (denoted by C_0, see Eq. (7)). From the theoretical point of view, dTLp differs significantly from the Lagrangian fluid velocity correlation time (Reeks, 1977; Wang and Stock, 1993), and this carries the essential nonlinearity in the statistical modeling of particle motion. dTLp and C0 may depend on the filter width and particle Stokes number even for a given turbulent flow. In previous studies, dTLp is modeled either by the fluid SGS Lagrangian timescale (Fede et al., 2006; Shotorban and Mashayek, 2006b; Pozorski and Apte, 2009; Bini and Jones, 2008) or by a simple extension of the timescale obtained from the full flow field (Berrouk et al., 2007). In this work, we shall study the subtle and on-monotonic dependence of $\delt T_{Lp}$ on the filter width and particle Stokes number using a flow field obtained from Direct Numerical Simulation (DNS). We then propose an empirical closure model for $\delta T_{Lp}$. Finally, the model is validated against LES of particle-laden turbulence in predicting single-particle statistics such as particle kinetic energy. As a first step, we consider the particle motion under the one-way coupling assumption in isotropic turbulent flow and neglect the gravitational settling effect. The one-way coupling assumption is only valid for low particle mass loading.

Aerodynamic characteristics research on wide-speed range waverider configuration

Waverider generated from a given flow field has a high lift-to-drag ratio because of attached bow shock on leading edge. However, leading edge blunt and off-design condition can make bow shock off leading edge and have unfavorable influence on aerodynamic characteristics. So these two problems have always been concerned as important engineering science issues by aeronautical engineering scientists. In this paper, through respectively using low speed and high speed waverider design principles, a wide-speed rang vehicle is designed, which can level takeoff and accelerate to hypersonic speed for cruise. In addition, sharp leading edge is blunted to alleviated aeroheating. Theoretical study and wind tunnel test show that this vehicle has good aerodynamic performance in wide-speed range of subsonic, transonic, supersonic and hypersonic speeds.

硅生长过程中尖形磁场对流场及提拉速率的影响

通过使用基于非正交网格有限体积法的FLUTRAPP(fluid flow and transport phenomena pro-gram)程序模拟了工业8 in单晶硅提拉法生长过程.数值模拟结果表明,流场及提拉速率在单晶硅生长过程中具有振荡特性,提拉速率的振荡周期大约为2 min.尖形磁场的引入能够抑制坩埚中熔体流动的振荡,减小提拉速率的振幅,从而有利于提高所生长单晶的质量.

Spatial resolved temperature measurement based on absorption spectroscopy using a single tunable diode laser

A novel method based on wavelength-multiplexed line-of-sight absorption and profile fitting for non-uniform flow field measurement is reported. A wavelength scanning combing laser temperature and current modulation WMS scheme is used to implement the wavelength-multiplexed-profile fitting method. Second harmonic (2f) signal of eight H2O transitions features near 7,170 cm(-1) are measured in one period using a single tunable diode laser. Spatial resolved temperature distribution upon a CH4/air premixed flat flame burner is obtained. The result validates the feasibility of strategy for non-uniform flow field diagnostics by means of WMS-2f TDLAS.

Experimental Investigation of Performance of a Miniature Direct Methanol Fuel Cell in Short-Term Microgravity

Experimental study of a liquid fed direct methanol fuel cell has been conducted in different gravity environments. A small single cell with 5 cm x 5 cm active area has single serpentine channel on the graphite cathode polar plate and 11 parallel straight channels on the graphite anode flow bed. Cell voltage and current have been measured and two-phase flow in anode channels has been in situ visually observed. The experimental results indicate that the effect of gravity on power performance of the direct methanol fuel cell is large when the concentration polarization governs fuel cells operation. Gravitational effect becomes larger at higher current density. Increasing methanol feeding molarity is conducive to weaken the influence of gravity on performance of liquid fed direct methanol fuel cells. Increasing feeding flow rate of methanol solution from 6 to 15 ml/min could reduce the size of carbon dioxide bubbles, while the influence of gravity still exist. Transport phenomena inside direct methanol fuel cells in microgravity is also analyzed and discussed.

A study on water film in saturated sand

Water film can serve as a sliding surface and cause landslides on gentle slopes. The development of "water film" in saturated sand is analyzed numerically and theoretically based on a quasi-three-phase model. It is shown that stable water films initiate and grow if the choking state (where the fluid velocity decreases to near zero) remains steady in a liquefied sand column. Discontinuity can occur in pore water velocity, grain velocity and pore pressure after the initiation of a water film. However, the discontinuity and water film can disappear once the choking state is changed. The key to the formation of water film is the choking in the sand column caused by eroded fine grains.

Numerical Study on the Distribution Characteristics of TP in Taihu Lake, China

A depth-averaged two-dimensional hydrodynamics model and a two-dimensional water quality model, bases on the alternating direction implicit (ADI) method, is developed to study the distribution characteristics of total phosphorus (TP) in Taihu Lake. Wind stress, the pollution source from the inflow rivers, releasing rate of TP from the bottom sediment and water diversion from Yangtze River are the effecting factors of TP distribution. By using the model proposed in this paper, the concentration field of the total phosphorus was simulated, which leads to the conclusion that the flow field has a great influence on the spatial and temporal distribution of TP in the Taihu Lake.

Analysis of current induced by long internal solitary waves in stratified ocean

An approximate theoretical expression for the current induced by long internal solitary waves is presented when the ocean is continuously or two-layer stratified. Particular attention is paid to characterizing velocity fields in terms of magnitude, flow components, and their temporal evolution/spatial distribution. For the two-layer case, the effects of the upper/lower layer depths and the relative layer density difference upon the induced current are further studied. The results show that the horizontal components are basically uniform in each layer with a shear at the interface. In contrast, the vertical counterparts vary monotonically in the direction of the water depth in each layer while they change sign across the interface or when the wave peak passes through. In addition, though the vertical components are generally one order of magnitude smaller than the horizontal ones, they can never be neglected in predicting the heave response of floating platforms in gravitationally neutral balance. Comparisons are made between the partial theoretical results and the observational field data. Future research directions regarding the internal wave induced flow field are also indicated.

渗流场流线构造方法研究

流线是分析渗流场流体流动特性的一个非常重要的工具,其构造方法的研究已受到油藏工程师们的普遍关注. 给出了基于约束Delaunay三角网格的等值线和流线生成算法. 提出了一种新的计算渗流场压力梯度的方法,借助于"翼边"数据结构,设计了"指南针"算法加速流线的跟踪过程. 实例表明,算法具有良好的时空效率,结果能够较为直观地反映出油藏流体在注采井间的运动轨迹

Numerical Study on the Spatial and Temporal Characteristics of Water Quality in the Taihu Lake of China

Taihu Lake is the third largest fresh water lake in China. With the fast economic development, abundant industrial and agricultural waste water has been discharged into Taihu Lake, causing the eutrophication of the water quality, which greatly affected the water utility. In the past decades, the treatment of Taihu Lake has witnessed limited success. Therefore, it is practically and theoretically significant to study the eutrophication of Taihu Lake. This research has focused on the issue of water quality including the characteristics of spatial and temporal distributions, and the rules of nutrient diffusion in the Taihu lake area. Based on the monitoring data, the basis distribution characteristics of water quality in Taihu Lake are analyzed. Comparing Taihu Lake with other Lakes shows that one important reason for Taihu eutrophication is the long period of water retention. A transporting and diffusing model of Taihu nutrient is developed by combining with the hydrodynamics model. Using the model, the concentration field of the total phosphorus (TP) and the influence of wind-driven current are numerically investigated, which leads to the conclusion that the flow field has a great influence on the spatial and temporal distributions of TP in Taihu Lake. Furthermore, the effect for improving the water quality by the project of water diversion from the Yangtze River to Taihu Lake was analyzed by simulation. The results demonstrate that short-term water diversion cannot improve the water quality of the heavily-polluted Meiliang Bay and the western bank areas of Taihu Lake.

迷宫流道转角对灌水器抗堵塞性能的影响

以转角分别为45.0°、60.0°、67.5°和75.0°的齿形流道灌水器为研究对象,应用CFD流场速度数值分析、PIV颗粒运动轨迹线和速度观测对比以及浑水抗堵塞测试相结合的方法,研究了转角对灌水器水力性能和抗堵塞能力的影响。结果表明转角与流量系数及流态指数均呈负相关关系,而灌水器的抗堵塞能力随着转角的增加呈下降趋势。综合分析转角对水力性能和抗堵塞性能的影响,提出迷宫流道结构灌水器的合理转角为60.0°。