A NUMERICAL STUDY OF THE TRAIN-INDUCED UNSTEADY AIRFLOW IN A SUBWAY TUNNEL WITH NATURAL VENTILATION DUCTS USING THE DYNAMIC LAYERING METHOD

The objective of this study is to investigate numerically the characteristics of train-induced unsteady airflow in a subway tunnel with natural ventilation ducts.A three-dimensional numerical model using the dynamic layering method for the moving boundary of a train is first developed,and then it is validated against the model tunnel experimental data.With the tunnel and subway train geometries in the numerical model exactly the same as those in the model tunnel experimental test,but with the ventilation ducts being connected to the tunnel ceiling and a barrier placed at the tunnel outlet,the three-dimensional train-induced unsteady tunnel flows are numerically simulated.The computed distributions of the pressure and the air velocity in the tunnel as well as the time series of the mass flow rate at the ventilation ducts reveal the impact of the train motion on the exhaust and suction of the air through ventilation ducts and the effects of a barrier placed at the tunnel outlet on the duct ventilation performance.As the train approaches a ventilation duct,the air is pushed out of the tunnel through the duct.As the train passes the ventilation duct,the exhaust flow in the duct is changed rapidly to the suction flow.After the train passes the duct,the suction mass flow rate at the duct decreases with time since the air pressure at the opening of the duct is gradually recovered with time.A drastic change in the mass flow rate at a ventilation duct while a train passes the corresponding ventilation duct,causes a change in the exhaust mass flow rate at other ventilation ducts.Also,when a barrier is placed at the tunnel outlet,the air volume discharge rate at each ventilation duct is greatly increased,i.e.,the barrier placed at the tunnel outlet can improve remarkably the ventilation performance through each duct.

A Dynamic Adaptive Layered Multicast Congestion Control Mechanism

To solve the problem that most of existing layered multicast protocols cannot adapt todynamic network conditions because their layers are coarsely granulated and static,a new congestioncontrol mechanism for dynamic adaptive layered multicast(DALM) is presented.In this mechanism,anovel feedback aggregating algorithm is put forward,which can dynamically determine the number oflayers and the rate of each layer,and can efficiently improve network bandwidth utilization ratio.Additionally,because all layers is transmitted in only one group,the intricate and time-consuminginternet group management protocol(IGMP) operations,caused by receiver joining a new layer orleaving the topmost subscribed layer,are thoroughly eliminated.And this mechanism also avoids otherproblems resulted from multiple groups.Simulation results show that DALM is adaptive and TCPfriendly.

Design and Analysis of a New AUV's Sliding Control System Based on Dynamic Boundary Layer

The new AUV driven by multi-vectored thrusters not only has unique kinematic characteristics during the actual cruise but also exists uncertain factors such as hydrodynamic coefficients perturbation and unknown interference of tail fluid, which bring difficult to the stability of the AUV＇s control system. In order to solve the nonlinear term and unmodeled dynamics existing in the new AUV＇s attitude control and the disturbances caused by the external marine environment, a second-order sliding mode controller with double-loop structure that considering the dynamic characteristics of the rudder actuators is designed, which improves the robustness of the system and avoids the control failure caused by the problem that the design theory of the sliding mode controller does not match with the actual application conditions. In order to avoid the loss of the sliding mode caused by the amplitude and rate constraints of the rudder actuator in the new AUV＇s attitude control, the dynamic boundary layer method is used to adjust the sliding boundary layer thickness so as to obtain the best anti-chattering effects. Then the impacts of system parameters, rudder actuator＇s constraints and boundary layer on the sliding mode controller are computed and analyzed to verify the effectiveness and robustness of the sliding mode controller based on dynamic boundary layer. The computational results show that the original divergent second-order sliding mode controller can still effectively implement the AUV＇s attitude control through dynamically adjusting the sliding boundary layer thickness. The dynamic boundary layer method ensures the stability of the system and does not exceed the amplitude constraint of the rudder actuator, which provides a theoretical guidance and technical support for the control system design of the new AUV in real complex sea conditions.

Dynamics of Bottom Boundary Layers in the Yellow River Subaqueous Delta Based on Long-Term In-Situ Observations

Objective In geo-marine science,the generalized bottom boundary layer（BBL）represents a layer between sediments and seawater.The BBL plays an important role in geological,geobiochemical,geophysical and geotechnical research because it is the connection region of hydrosphere,

Influence of floc dynamic protection layer on alleviating ultrafiltration membrane fouling induced by humic substances

Cake layer formation is inevitable over time for ultrafiltration(UF)membrane-based drinking water treatment.Although the cake layer is always considered to cause membrane fouling,it can also act as a"dynamic protection layer",as it further adsorbs pollutants and dramatically reduces their chance of getting to the membrane surface.Here,the UF membrane fouling performance was investigated with pre-deposited loose flocs in the presence of humic acid(HA).The results showed that the floc dynamic protection layer played an important role in removing HA.The higher the solution pH,the more negative the floc charge,resulting in lower HA removal efficiency due to the electrostatic repulsion and large pore size of the floc layer.With decreasing solution pH,a positively charged floc dynamic protection layer was formed,and more HA molecules were adsorbed.The potential reasons were ascribed to the smaller floc size,greater positive charge,and higher roughness of the floc layer.However,similar membrane fouling performance was also observed for the negative and positive floc dynamic protection layers due to their strong looseness characteristics.In addition,the molecular weight(MW)distribution of HA also played an important role in UF membrane fouling behavior.For the small MW HA molecules,the chance of forming a loose cake layer was high with a negatively charged floc dynamic protection layer,while for the large MW HA molecules it was high with a positively charged floc dynamic protection layer.As a result,slight UF membrane fouling was induced.

THE WARMING MECHANISM IN THE SOUTHERN ARABIAN SEA DURING THE DEVELOPMENT OF INDIAN OCEAN DIPOLE EVENTS

This study aims to explore the relative role of oceanic dynamics and surface heat fluxes in the warming of southern Arabian Sea and southwest Indian Ocean during the development of Indian Ocean Dipole(IOD) events by using National Center for Environmental Prediction/National Center for Atmospheric Research(NCEP/NCAR) daily reanalysis data and Global Ocean Data Assimilation System(GODAS) monthly mean ocean reanalysis data from 1982 to2013,based on regression analysis,Empirical Orthogonal Function(EOF) analysis and combined with a 21/2layer dynamic upper-ocean model.The results show that during the initial stage of IOD events,warm downwelling Rossby waves excited by an anomalous anticyclone over the west Indian Peninsula,southwest Indian Ocean and southeast Indian Ocean lead to the warming of the mixed layer by reducing entrainment cooling.An anomalous anticyclone over the west Indian Peninsula weakens the wind over the Arabian Sea and Somali coast,which helps decrease the sea surface heat loss and shallow the surface mixed layer,and also contributes to the sea surface temperature(SST) warming in the southern Arabian Sea by inhibiting entrainment.The weakened winds increase the SST along the Somali coast by inhibiting upwelling and zonal advection.The wind and net sea surface heat flux anomalies are not significant over the southwest Indian Ocean.During the antecedent stage of IOD events,the warming of the southern Arabian Sea is closely connected with the reduction of entrainment cooling caused by the Rossby waves and the weakened wind.With the appearance of an equatorial easterly wind anomaly,the warming of the southwest Indian Ocean is not only driven by weaker entrainment cooling caused by the Rossby waves,but also by the meridional heat transport carried by Ekman flow.The anomalous sea surface heat flux plays a key role to damp the warming of the west pole of the IOD.

Chemical composition and formation mechanisms in the cathode-electrolyte interface layer of lithium manganese oxide batteries from reactive force field （ReaxFF） based molecular dynamics

Lithium manganese oxide （LiMn2O4） is a principal cathode material for high power and high energy density electrochemical storage on account of its low cost, non-toxicity, and ease of preparation relative to other cathode materials. However, there are well-documented problems with capacity fade of lithium ion batteries containing LiMn2O4. Experimental observations indicate that the manganese content of the electrolyte increases as an electrochemical cell containing LiMn2O4 ages, suggesting that active material loss by dissolution of divalent manganese from the LiMn2O4 surface is the primary reason for reduced cell life in LiMn2O4 batteries. To improve the retention of manganese in the active material, it is key to understand the reactions that occur at the cathode surface. Although a thin layer of electrolyte decomposition products is known to form at the cathode surface, the speciation and reaction mechanisms of Mn^2＋ in this interface layer are not yet well understood. To bridge this knowledge gap, reactive force field （ReaxFF） based molecular dynamics was applied to investigate the reactions occurring at the LiMn2O4 cathode surface and the mechanisms that lead to manganese dissolution. The ReaxFFMD simulations reveal that the cathode-electrolyte interface layer is composed of oxida- tion products of electrolyte solvent molecules including aldehydes, esters, alcohols, polycarbonates, and organic radicals. The oxidation reaction pathways for the electro- lyre solvent molecules involve the formation of surface hydroxyl species that react with exposed manganese atoms on the cathode surface. The presence of hydrogen fluoride （HF） induces formation of inorganic metal fluorides and surface hydroxyl species. Reaction products predicted by ReaxFF-based MD are in agreement with experimentally identified cathode-electrolyte interface compounds. An overall cathode-electrolyte interface reaction scheme is proposed based on the molecular simulation results.

Numerical Analyses of Caisson Breakwaters on Soft Foundations Under Wave Cyclic Loading

A caisson breakwater is built on soft foundations after replacing the upper soft layer with sand. This paper presents a dynamic finite element method to investigate the strength degradation and associated pore pressure development of the intercalated soft layer under wave cyclic loading. By combining the undrained shear strength with the empirical formula of overconsolidation clay produced by unloading and the development model of pore pressure, the dynamic degradation law that describes the undrained shear strength as a function of cycle number and stress level is derived. Based on the proposed dynamic degradation law and M-C yield criterion, a dynamic finite element method is numerically implemented to predict changes in undrained shear strength of the intercalated soft layer by using the general-purpose FEM software ABAQUS, and the accuracy of the method is verified. The effects of cycle number and amplitude of the wave force on the degradation of the undrained shear strength of the intercalated soft layer and the associated excess pore pressure response are investigated by analyzing an overall distribution and three typical sections underneath the breakwater. By comparing the undrained shear strength distributions obtained by the static method and the quasi-static method with the undrained shear strength distributions obtained by the dynamic finite element method in the three typical sections, the superiority of the dynamic finite element method in predicting changes in undrained shear strength is demonstrated.

Mechanisms associated with winter intraseasonal extreme sea ice extent in the Weddell Sea

Previous studies have shown evidence of atmospheric extratropical wave trains modulating sea ice area in the Weddell and Amundsen/Bellingshausen seas on intraseasonal time-scales(20–100 d). Here we investigate mechanisms relating intraseasonal extreme sea ice extent and Ekman layer dynamics with emphasis on the Weddell Sea. This study extends from 1989 to 2013 and focuses on the winter season. Wind stress τ is calculated with winds from the Climate Forecast System reanalysis(CFSR) to evaluate momentum transfer between the atmosphere and the Ekman layer. Lag-composites of the anomalies of Ekman transport and the Ekman pumping indicate that divergence of mass in the Ekman layer and upwelling lead the occurrence of extreme sea ice contraction on intraseasonal time-scales in the Weddell Sea. Opposite conditions(i.e., convergence of the mass and downwelling) lead extreme sea ice expansion on intraseasonal time-scales. This study suggests that the Ekman pumping resulting from the anomalous wind stress on intraseasonal time-scales can transport these warmer waters to the surface contributing to sea ice melting. Additionally, high resolution sea ice fraction and ocean currents obtained from satellite and in situ data are used to investigate in detail mechanisms associated with persistent extreme sea ice expansion and contraction on intraseasonal time-scales. These case studies reveal that atmospheric circumpolar waves on intraseasonal time-scales can induce contrasting anomalies of about ±20% in sea ice concentration at the Weddell and western Antarctica Peninsula margins within less than 30 d. This study shows that extreme anomalies in sea ice may lag between 5–25 d(1–5 pentads) the ocean-atmospheric forcing on intraseasonal time-scales.

Analysis on capabilities of density-based solvers within OpenFOAM to distinguish aerothermal variables in difusion boundary layer

Open source feld operation and manipulation（OpenFOAM）is one of the most prevalent open source computational fluid dynamics（CFD）software.It is very convenient for researchers to develop their own codes based on the class library toolbox within OpenFOAM.In recent years,several density-based solvers within OpenFOAM for supersonic/hypersonic compressible flow are coming up.Although the capabilities of these solvers to capture shock wave have already been verifed by some researchers,these solvers still need to be validated comprehensively as commercial CFD software.In boundary layer where diffusion is the dominant transportation manner,the convective discrete schemes＇capability to capture aerothermal variables,such as temperature and heat flux,is different from each other due to their own numerical dissipative characteristics and from viewpoint of this capability,these compressible solvers within OpenFOAM can be validated further.In this paper,frstly,the organizational architecture of density-based solvers within OpenFOAM is analyzed.Then,from the viewpoint of the capability to capture aerothermal variables,the numerical results of several typical geometrical felds predicted by these solvers are compared with both the outcome obtained from the commercial software Fastran and the experimental data.During the computing process,the Roe,AUSM＋（Advection Upstream Splitting Method）,and HLLC（Harten-Lax-van Leer-Contact）convective discrete schemes of which the spatial accuracy is 1st and 2nd order are utilized,respectively.The compared results show that the aerothermal variables are in agreement with results generated by Fastran and the experimental data even if the1st order spatial precision is implemented.Overall,the accuracy of these density-based solvers can meet the requirement of engineering and scientifc problems to capture aerothermal variables in diffusion boundary layer.

Numerical study of the effects of Planetary Boundary Layer structure on the pollutant dispersion within built-up areas

The effects of different Planetary Boundary Layer（PBL） structures on pollutant dispersion processes within two idealized street canyon configurations and a realistic urban area were numerically examined by a Computational Fluid Dynamics（CFD） model. The boundary conditions of different PBL structures/conditions were provided by simulations of the Weather Researching and Forecasting model. The simulated results of the idealized 2D and 3D street canyon experiments showed that the increment of PBL instability favored the downward transport of momentum from the upper flow above the roof to the pedestrian level within the street canyon. As a result, the flow and turbulent fields within the street canyon under the more unstable PBL condition are stronger. Therefore, more pollutants within the street canyon would be removed by the stronger advection and turbulent diffusion processes under the unstable PBL condition. On the contrary, more pollutants would be concentrated in the street canyon under the stable PBL condition. In addition, the simulations of the realistic building cluster experiments showed that the density of buildings was a crucial factor determining the dynamic effects of the PBL structure on the flow patterns. The momentum field within a denser building configuration was mostly transported from the upper flow, and was more sensitive to the PBL structures than that of the sparser building configuration. Finally, it was recommended to use the Mellor-Yamada-Nakanishi-Niino（MYNN） PBL scheme, which can explicitly output the needed turbulent variables, to provide the boundary conditions to the CFD simulation.