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.

Numerical Analyses of Transonic Fluid/Structure interaction in Aircraft Engineering

Through the coupling between aerodynamic and structural governing equations, a fully implicit multiblock aeroelastic solver was developed for transonic fluid/stricture interaction. The Navier-Stokes fluid equations are solved based on LU-SGS (lower-upper symmetric Gauss-Seidel) Time-marching subiteration scheme and HLLEW (Harten-Lax-van Leer-Einfeldt-Wada) spacing discretization scheme and the same subiteration formulation is applied directly to the structural equations of motion in generalized coordinates. Transfinite interpolation (TFI) is used for the grid deformation of blocks neighboring the flexible surfaces. The infinite plate spline (IPS) and the principal of virtual work are utilized for the data transformation between fluid and structure. The developed code was fort validated through the comparison of experimental and computational results for the AGARD 445.6 standard aeroelastic wing. In the subsonic and transonic range, the calculated flutter speeds and frequencies agree well with experimental data, however, in the supersonic range, the present calculation overpredicts the experimental flutter points similar to other computations. Then the flutter character of a complete aircraft configuration is analyzed through the calculation of the change of structural stiffness. Finally, the phenomenon of aileron buzz is simulated for the weakened model of a supersonic transport wing/body model at Mach numbers of 0.98 and l.05. The calculated unsteady flow shows, on the upper surface, the shock wave becomes stronger as the aileron deflects downward, and the flow behaves just contrary on the lower surface of the wing. Corresponding to general theoretical analysis, the flow instability referred to as aileron buzz is induced by a stronger shock alternately moving on the upper and lower surfaces of wing. For the rigid structural model, the flow is stable at all calculated Mach numbers as observed in experiment

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.

Donor-donor binding in In2O3: Engineering shallow donor levels

Using first-principles band structure methods, we investigate the interactions between different donors in In2O3. Through the formation energy and transition energy level calculations, we find that an oxygen-vacancy creates a deep donor level, while an indium-interstitial or a tin-dopant induces a shallow donor level. The coupling between these donor levels gives rise to even shallower donor levels and leads to a significant reduction in their formation energies. Based on the analysis of the PBE0-corrected band structure and the molecular-orbital bonding diagram, we demonstrate these effects of donor-donor binding. In addition, total energy calculations show that these defect pairs tend to be more stable with respect to the isolated defects due to their negative binding energies. Thus, we may design shallow donor levels to enhance the electrical conductivity via the donor donor binding.

Mechanism for Giycinebetaine to Improve Plant Salt Resistance and Its Research Advances in Genetic Engineering

The paper systematically discusses the mechanism for glycinebetaine to improve plant salt resistance and its research advances in genetic engineering at home and abroad as well as summarizing the research progresses about the key enzymes and their genetic engineering in glycinebetaine biosynthesis. It suggests that on the basis of further understanding the mechanism for glycinebetaine to improve plant salt resistance,the transformation of the genes relating to glycinebetaine biosynthesis should be carried out in major crops so that new plant varieties resistant to salt can be obtained.

Case Study: Influence Of Morphological Changes On Flooding In Jingjiang River

The middle reach of the Yangtze River, customarily called the Jingjiang River, together with its diversion channels and Dongting Lake, form a large complicated drainage system. In the last five decades, significant geomorphological changes have occurred in the drainage system, including the shrinkage of diversion channels, contraction of Dongting Lake, changes in the rating curve at the Luoshan station, and cutoffs of the lower Jingjiang River. These changes are believed to be the cause of the occurrence of abnormal floods in the Jingjiang River. Qualitative analyses suggest that the first three factors aggravate the flood situation in the lower Jingjiang River, while the last factor seems beneficial for flood prevention. To quantitatively evaluate these conclusions, a finite-volume numerical model was constructed. A series of numerical simulations were carried out to test the individual and combined effects of the aforementioned four factors, and these simulations showed that high flood stages in the Jingjiang River clearly are related to the geomorphological changes.

Impacts of socio-economic factors on sediment yield in the Upper Yangtze River

In recent years, the role of human activities in changing sediment yield has become more apparent for the construction of hydraulic engineering and water conservation projections in the Upper Yangtze River, but it has not been evaluated at the macro scale. Taking Sichuan Province and Chongqing City as an example, this paper studies the relationship between socio-economic factors and sediment yield in the Upper Yangtze River based on section data in 1989 and 2007. The results show that sediment yield is significantly correlated with population density and cultivated area, in which the former appears to be more closely related to sediment yield. Moreover, in the relation of sediment yield vs. population density, a critical value of population density exists, below which the sediment yield increases with the increase of population density and over which the sediment yield increases with the decrease of population density. The phenomenon essentially reflects the influence of natural factors, such as topography, precipitation and soil property, and some human activities on sediment yield. The region with a higher population density than critical value is located in the east of the study area and is characterized by plains, hills and low mountains, whereas the opposite is located in the west and characterized by middle and high mountains. In the eastern region, more people live on the lands with a low slope where regional soil erosion is slight; therefore, sediment yield is negatively related with population density. In contrast, in the western region, the population tends to aggregate in the areas with abundant soil and water resources which usually lead to a higher intensity of natural erosion, and in turn, high-intensity agricultural practices in these areas may further strengthen local soil erosion. It is also found that population tends to move from the areas with bad environment and high sediment yield to the areas with more comfortable environment and less sediment yield. The natural factors have greater influence on sediment yield of western region than that of eastern region. Generally, the natural factors play a dominant role on sediment yield in the Upper Yangtze River.

Simulation of Free-Surface Flow in a Tank Using the Navier-Stokes Model and Unstructured Finite Volume Method

Modelling free-surface flow has very important applications in many engineering areas such as oil transportation and offshore structures. Current research focuses on the modelling of free surface flow in a tank by solving the Navier-Stokes equation. An unstructured finite volume method is used to discretize the governing equations. The free surface is tracked by dynamically adapting the mesh and making it always surface conforming. A mesh-smoothing scheme based on the spring analogy is also implemented to ensure mesh quality throughout the computaiton. Studies are performed on the sloshing response of a liquid in an elastic container subjected to various excitation frequencies. Further investigations are also carried out on the critical frequency that leads to large deformation of the tank walls. Another numerical simulation involves the free-surface flow past as submerged obstacle placed in the tank to show the flow separation and vortices. All these cases demonstrate the capability of this numerical method in modelling complicated practical problems.

Validity ranges of interfacial wave theories in a two-layer fluid system

In the present paper, we endeavor to accomplish a diagram, which demarcates the validity ranges for interfacial wave theories in a two-layer system, to meet the needs of design in ocean engineering. On the basis of the available solutions of periodic and solitary waves, we propose a guideline as principle to identify the validity regions of the interfacial wave theories in terms of wave period T, wave height H, upper layer thickness d(1), and lower layer thickness d(2), instead of only one parameter-water depth d as in the water surface wave circumstance. The diagram proposed here happens to be Le Mehautes plot for free surface waves if water depth ratio r = d(1)/d(2) approaches to infinity and the upper layer water density rho(1) to zero. On the contrary, the diagram for water surface waves can be used for two-layer interfacial waves if gravity acceleration g in it is replaced by the reduced gravity defined in this study under the condition of sigma = (rho(2) - rho(1))/rho(2) -> 1.0 and r > 1.0. In the end, several figures of the validity ranges for various interfacial wave theories in the two-layer fluid are given and compared with the results for surface waves.

Numerical modeling of 3-D turbulent two-phase flow and coal combustion in a pulverized-coal combustor

In the present paper, a multifluid model of two-phase flows with pulverized-coal combustion, based on a continuum-trajectory model with reacting particle phase, is developed and employed to simulate the 3-D turbulent two-phase hows and combustion in a new type of pulverized-coal combustor with one primary-air jet placed along the wall of the combustor. The results show that: (1) this continuum-trajectory model with reacting particle phase can be used in practical engineering to qualitatively predict the flame stability, concentrations of gas species, possibilities of slag formation and soot deposition, etc.; (2) large recirculation zones can be created in the combustor, which is favorable to the ignition and flame stabilization.

A Mechanical Model For The Quantification Of The Effect Of Laser Quenching On Ctod In Steels

The technology of laser quenching is widely used to improve the surface properties of steels in surface engineering. Generally, laser quenching of steels can lead to two important results. One is the generation of residual stress in the surface layer. In general, the residual stress varies from the surface to the interior along the quenched track depth direction, and the residual stress variation is termed as residual stress gradient effect in this work. The other is the change of mechanical properties of the surface layer, such as the increases of the micro-hardness, resulting from the changes of the microstructure of the surface layer. In this work, a mechanical model of a laser-quenched specimen with a crack in the middle of the quenched layer is developed to quantify the effect of residual stress gradient and the average micro-hardness over the crack length on crack tip opening displacement (CTOD). It is assumed that the crack in the middle of the quenched layer is created after laser quenching, and the crack can be a pre-crack or a defect due to some reasons, such as a void, cavity or a micro-crack. Based on the elastic-plastic fracture mechanics theory and using the relationship between the micro-hardness and yield strength, a concise analytical solution, which can be used to quantify the effect of residual stress gradient and the average micro-hardness over the crack length resulting from laser quenching on CTOD, is obtained. The concise analytical solution obtained in this work, cannot only be used as a means to predict the crack driving force in terms of the CTOD, but also serve as a baseline for further experimental investigation of the effect after laser-quenching treatment on fracture toughness in terms of the critical CTOD of a specimen, accounting for the laser-quenching effect. A numerical example presented in this work shows that the CTOD of the quenched can be significantly decreased in comparison with that of the unquenched. (C) 2008 Elsevier B.V. All rights reserved.

Generalized Viscoplastic Modeling of Debris Flow

Various concepts have been proposed or used in the development of rheological models for debris flow. The earliest model developed by Bagnold was based on the concept of the “dispersive” pressure generated by grain collisions. Bagnold’s concept appears to be theoretically sound, but his empirical model has been found to be inconsistent with most theoretical models developed from non-Newtonian fluid mechanics. Although the generality of Bagnold’s model is still at issue, debris-flow modelers in Japan have generally accepted Takahashi’s formulas derived from Bagnold’s model. Some efforts have recently been made by theoreticians in non-Newtonian fluid mechanics to modify or improve Bagnold’s concept or model. A viable rheological model should consist both of a rate-independent part and a rate-dependent part. A generalized viscoplastic fluid (GVF) model that has both parts as well as two major rheological properties (i.e., the normal stress effect and soil yield criterion) is shown to be sufficiently accurate, yet practical, for general use in debris-flow modeling. In fact, Bagnold’s model is found to be only a particular case of the GVF model. Analytical solutions for (steady) uniform debris flows in wide channels are obtained from the GVF model based on Bagnold’s simplified assumption of constant grain concentration.

森林流域单元水文模型的构建与应用

目前全球缺水、水污染、洪涝灾害以及水土流失仍然非常严重，尤其在我国北方地区。流域水文模型可用来进行不同需水管理的情景分析，为解决我国水问题提供科学依据。分布式水文模型是流域水文模型的发展方向，具有显著特点：1）应用前景广泛，不仅可以模拟流域水文过程，还可以协助模拟泥沙或污染物的运移过程，为水利工程设计、水土保持、环境保护等领域提供技术支持；2）能够预测流域土地利用或气候变化下的流域水文响应过程变化，为管理部门提供决策支持；3）模型所需要的参数全部具有物理意义，可通过实际测量确定，适合模拟实测系列较短或是无观测流域的水文过程；4）对于目前国际水文界的前沿问题—水文尺度转换提供了一种有效的解决途径。 然而分布式水文模型还不完善，如1)真实性问题。对一些水文过程和边界条件还不确定。2)尺度转换问题。目前很少考虑尺度对参数有效性的影响。3)检验问题。还无法判断对有些难以测量的水文状态变量的模拟正确与否。4)计算时间和数据存储的问题。有些分布式水文模型虽然具有很强的水文物理基础和完善的模型结构，但是计算时间过长和(或)数据存储过大，难以应用。上述问题的核心就是对分布式水文模型的核心—单元水文模型的研究不够，需要为进一步完善单元水文模型进行研究。 本文采用饱和入渗理论、Saint-Venant方程、Richards方程、Penman-Monteith方程等等构建了以有限差分法求解的适用于森林流域的单元水文模型，并通过实验室模拟试验和坡地径流场资料进行了验证，主要结论为： 通过不同坡度和不同雨强下的室内坡面产汇流实验模拟，表明：该模型模拟的坡面流和壤中流过程与实测过程基本一致，峰现时间、径流历时、峰值流量、出流总量模拟值与实测值的相对误差均较小，基本小于10%。模型的模拟精度较高，实用性较强，为深入研究壤中流机制和改进流域降雨-径流模型提供了理论依据。 通过坡地径流观测场实测资料的验证，表明：该模型模拟的坡面流过程精度较高，累计流量的精度更高于小时过程的精度，离差系数、效率系数、确定系数均较理想，具有应用价值，有助于改善分布式水文模型在森林流域的模拟效果。

亚硝化—厌氧氨氧化处理猪场废水研究

猪场废水COD浓度高、氨氮浓度高、悬浮物浓度高，已成为农村面源污染的主要来源，并严重威胁到农村饮用水安全。猪场废水氨氮浓度高、处理难度大，如何采用经济高效的方法，去除氨氮使其达到排放标准，一直是猪场废水处理中面临的重要难题。 厌氧氨氧化是近年受到国内外水处理研究者广泛关注的新型生物脱氮技术，具有不需要外加有机碳源、节省供氧量、降低能耗等优点。虽然国内外研究者对厌氧氨氧化过程的脱氮机理、厌氧氨氧化菌的生理生化特性等进行了多方面的研究，但已有的报道大多以模拟废水为研究对象，以猪场废水为研究对象的报道，在国内外文献中极少有报导。 本论文以猪场废水为主要研究对象，考察了猪场废水的亚硝化过程、厌氧氨氧化的启动过程，并对亚硝化和厌氧氨氧化联合用于猪场废水脱氮进行了探索。 1．论文首先研究了猪场废水的亚硝化过程，考察了废水水质和主要运行条件对亚硝化过程的影响。实验表明：（1）亚硝化阶段反应时间为8到10h时，出水中氨氮和亚硝酸盐浓度比可达到1：1～1：1.23，满足厌氧氨氧化反应对二者比例的要求；达到前述要求时，氨氮去除率达到58.3～65.6 %，亚硝化率在整个过程均保持在97 %以上，COD去除率在59.2～68.6 %；（2）曝气量（溶解氧）对亚硝化过程影响显著，随着曝气量增大，达到厌氧氨氧化要求的氨氮与亚硝酸盐氮浓度比例所需水力停留时间τ越短，pH出现明显下降的时间越短；（3）τ对应的pH在7.8～8.1之间，无需进行pH调节即可满足厌氧氨氧化反应对pH的要求；（4）氨氮和COD降解过程遵循一级反应动力学，氨氮和COD降解的速率常数分别为0.0656～0.0724 1/h和0.0491～0.0664 1/h。 2．在进行亚硝化过程研究的同时，以模拟废水为试验对象，进行厌氧氨氧化启动研究。以反硝化污泥和养殖厂储水池厌氧底泥的混合污泥作为接种污泥，历时大约100天，培育出具有厌氧氨氧化活性的污泥，氨氮和亚硝酸盐氮最高进水浓度分别为223.8 mg/L和171.4 mg/L，去除率最高分别达48%和41.5%，此时二者消耗比例为1.33：1。 3．在猪场废水的亚硝化研究完成和厌氧氨氧化过程初步启动成功后，在模拟废水中逐步加入猪场废水的亚硝化处理出水，逐步实现亚硝化和厌氧氨氧化的组合。亚硝化出水添加到厌氧反应器后，厌氧氨氧化反应仍可继续进行，且去除效率逐步提高。研究发现添加的亚硝化出水中携带的亚硝化细菌在厌氧氨氧化菌膜外层生长并累积，增加了厌氧氨氧化反应基质的传质阻力，妨碍了厌氧氨氧化效率的提高。 4．亚硝化－厌氧氨氧化实际工程应用探索中，生物接触氧化池可在有效去除废水中的有机物的同时实现亚硝化，出水中氨氮和亚硝酸盐比例平均为1.10，可满足后续厌氧氨氧化的要求；在适宜的进水浓度和温度下，ABR池出现了厌氧氨氧化启动的迹象；研究同时发现，水质的波动和气温的变化是工程中影响厌氧氨氧化菌活性的重要因素。 论文的主要创新点在于：（1）以猪场废水为研究对象，以实现厌氧氨氧化为目标，对亚硝化过程进行了比较详细的考察，获得了亚硝化出水满足厌氧氨氧化要求的工艺条件，通过对其COD和氨氮降解过程的考察，得出亚硝化阶段COD降解和氨氮去除的动力学模型；（2）对亚硝化－厌氧氨氧化处理猪场废水进行了探索，确立了影响其污染物去除率稳定的重要因素。 论文的上述研究成果，为厌氧氨氧化技术的实用性研究提供理论依据。 Piggery wastewater, which is characterized by high concentration of COD、ammonium and suspend substance, has become a most important source of non-point source pollution and also severely threats drinking water security in rural area. How to discharge piggery wastewater with the ammonium concentration meeting standard by economical and effective method? This is the most urgent problem in piggery wastewater treatment. As a new biological nitrogen removal technology, Anammox process has been paid more and more attention by researchers all over the world. Anammox has advantages of no need of organic carbon addition, low oxygen consumption and energy consumption. Plenty of investigations have been carried out to the mechanism, physiological and biochemical characteristic of bacteria about Anammox. Most of researches focused on synthetic wastewater, there is rare report about its application in piggery wastewater. In this paper，experimental studies were performed to investigate Sharon process in treatment of piggery wastewater，the start up process of Annammox using synthetic wastewater were studied, the feasibility of applying Sharon-Anammox process in the nitrogen removal of piggery wastewater was evaluated. 1. Sharon process of piggery wastewater was firstly investigated to analyze the effects of water quality and main running parameters, which meet the NH4+-N to NO2--N ratio requirement of successive Anammox. Results showed: (1)During Sharon Process，after 8～10 hours’ reaction the NH4+-N to NO2--N ratio in effluent reached 1:1.0～1:1.23, when the removal percentage of NH4+-N was 58.3～65.6 %, a semi-nitration rate of above 97 % was achieved during the process; meanwhile 59.2～68.6 % of the COD was also removed. (2)The aeration rate（oxygen） had obvious effect on the hydraulic retention time(τ) which met the NH4+-N to NO2--N ratio requirement of Anammox. As aeration rate increased, the hydraulic retention time(τ) was shortened. (3) The pH corresponding to τ was between 7.8 and 8.1, thus it needed no artificial adjustment. (4) The reduction of ammonia and COD followed the first-order reaction kinetics. The velocity constants of ammonia and COD were 0.0656～0.0724 1/h and 0.0491～0.0664 1/h, respectively. 2. The startup of Anammox process using the artificial wastewater was performed simultaneously with Sharon. The aim was to investigate the running parameters of Anammox and make foundation for the combination stage. By using the mixture of denitrifying sludge and anaerobic sludge in tank of the breeding factory, sludge of Anammox activity was cultivated in UASB after 100 days. The removal percentage of NH4+-N and NO2-N were up to 48% and 41.5%, respectively, when the NH4+-N and NO2-N influent concentration were 223.8 mg/L and 171.4 mg/L, respectively, the NH4+-N and NO2-N removal rate was 1.33:1. 3. After investigation of Sharon and startup of Anammox, effluent of Sharon process was added into the synthetic wastewater to combine Sharon and Anammox step by step. It took some time after the addition of Sharon effluent that Anammox reaction continued and the removal rate kept increasing. It indicated that nitrifying bacteria were carried by the Sharon effluent cumulated in the outer layer of Anammox. This enhanced transfer resistance of Anammox reaction and the increasing removal rate was restrained. 4. In the bio-contact oxidation pond of practical project, Sharon process were carried out successfully and organic compounds were removed effectively. An average NO2-N/ NH4+-N rate of 1:1.0 was achieved in the effluent, which met the requirement of successive Anammox. Under condition of suitable influent concentration and temperature, there was evidence that Anammox could start up in ABR. The variety of wastewater and temperature had great affects on Anammox activity in practical engineering. Innovation of this paper: (1) The Sharon process for treating piggery wastewater was discussed in details. Technological parameters that met requirement of Anammox were obtained. The dynamic models of COD and ammonium removal in the process were educed. (2) Sharon-Ananmmox for treatment of piggery wastewater was investigated, and the primary influencing factors was studied. This paper could be a theoretical consult for research of Anammox utility.