On some applications of invariant manifolds

Taking transfer orbits of a collinear libration point probe, a lunar probe and an interplanetary probe as examples, some applications of stable and unstable invariant manifolds of the restricted three-body problem are discussed. Research shows that transfer energy is not necessarily conserved when invariant manifolds are used. For the cases in which the transfer energy is conserved, the cost is a much longer transfer time.     

Formation flying in elliptic orbits with the J2 perturbation

Relative dynamics between the chief satellite and the deputy ones in formation flying is crucial to maintaining the formation. A good choice of the formation usually requires a lower control frequency or less control energy. For formation flying missions in highly elliptic orbits, the well-known C-W equation is not accurate enough. Instead, Lawden’s equation is often used. First, the solution to Lawden’s equation with a very simple form is deduced. Then the J2 perturbation is added. It is found that Lawden’s solution is not necessarily valid when the J2 perturbation is considered. Completely discarding Lawden’s solution and borrowing the idea of mean orbit elements, two rules to initialize the formation are proposed. The deviation speed is greatly reduced. Different from previous studies on the J2 perturbation, except for the relatively simple expression for the semi-major axis, the tedious formulae of the long period terms and the short period terms of other orbital elements are not used. In addition, the deviation speed is further reduced by compensation of the nonlinear effects. Finally, a loose control strategy of the formation is proposed. To test the robustness of this strategy, a third body perturbation is added in numerical simulations.     

The role of mindanao dome in the variability of the Pacific North equatorial current bifurcation

The Mindanao Dome (MD) features prominent oceanic variability and is located geographically close to the bifurcation latitude Y _b of the Pacific North Equatorial Current. In this study, the role of the MD in the variability of Y _b is examined with 20 years of satellite altimetric sea surface height (SSH) data and a 1.5-layer linear Rossby wave model. It is shown that the seasonal variations of surface Y _b are related to not only the SSH fluctuations near the bifurcation point (bifurcation box; 125°–130°E, 12°–15°N) but also those outstanding in the MD region (MD box; 127°–132°E, 6°–9°N). The impact of the MD SSH changes is significant when the bifurcation point stays at southerly latitudes during February–September, which hinders (delays) the southward leap (northward retreat) of Y _b in April–May (July–August) and thus leads to the asymmetry of the mean Y _b seasonal cycle (with a positive skewness of γ = +0.64). Such asymmetry also shows year-to-year variations depending on yearly mean Y _b value. A southerly yearly mean Y _b involves larger contribution of the MD and thus causes larger asymmetry of Y _b seasonal cycle. At interannual and longer timescales, the MD acts to amplify the fluctuations of the bifurcation. It is responsible for about 20 % of the total low-frequency Y _b variances and plays an important role in the 0.12° year^?1 southward trend of Y _b in the past two decades. The impact of the MD on Y _b changes is becoming increasingly significant at various timescales such as the bifurcation point migrating southward in recent years.     

Updated Support Vector Machine for Seismic Liquefaction Evaluation Based on the Penetration Tests

Simplified techniques based on in situ testing methods are commonly used to predict liquefaction potential. Many of these simplified methods are based on finding the liquefaction boundary separating two categories (the occurrence or non-occurrence of liquefaction) through the analysis of liquefaction case histories. As the liquefaction classification problem is highly nonlinear in nature, it is difficult to develop a comprehensive model taking into account all the independent variables, such as the seismic and soil properties, using conventional modeling techniques. Hence, in many of the conventional methods that have been proposed, simplified assumptions have been made. In this study, an updated support vector machine (SVM) based on particle swarm optimization (PSO) is used to evaluate liquefaction potential in two separate case studies. One case is based on standard penetration test (SPT) data and the other is based on cone penetration test (CPT) data. The SVM model effectively explores the relationship between the independent and dependent variables without any assumptions about the relationship between the various variables. This study serves to demonstrate that the SVM can “discover” the intrinsic relationship between the seismic and soil parameters and the liquefaction potential. Comparisons indicate that the SVM models perform far better than the conventional methods in predicting the occurrence or non-occurrence of liquefaction.     

电活化硫酸盐处理四环素类废水的研究

本文构建以掺硼金刚石(Boron-doped diamond,BDD)为阳极、不锈钢为阴极、硫酸盐为电解质的电化学体系,考察了电流密度、pH值、硫酸盐浓度以及初始四环素浓度等四个因素对电化学氧化降解废水中四环素的影响,运用响应曲面法对运行参数进行优化;通过电子自旋共振检测技术分析电化学反应中产生的自由基,探究了间歇通电模式下电化学体系持续氧化机理。结果表明,四个因素对TOC去除率的影响大小次序为:电流密度>初始四环素浓度>初始pH值>硫酸盐浓度,其中初始pH值和硫酸盐浓度与电流密度和初始四环素浓度的交互作用对TOC去除率的影响较为显著;最佳运行参数为pH值为5,电流密度为100 mA·cm^-2,硫酸盐浓度为0.25 mol·L^-1,初始四环素浓度为1000 mg·L^-1;间歇通电模式下,BDD电极表面产生的SO₄^·-等高活性物质间相互转化提供了体系的可持续氧化能力。该研究结果为电化学氧化技术的实际应用提供了节省能耗的有效途径。     

半球形多向波聚合波道振荡水柱气室优化设计

为了解决振动水柱式波浪能转换装置收集多向波浪问题,本文设计了半球形多向聚合波道振荡水柱气室结构,以适合远海单点波浪能采集和发电。在规则波正向入射条件下,基于流体仿真分析软件(FLUENT)、流体动力学连续性假设和粘性不可压缩流体动量守恒的运动方程(Navier-Stokes方程)建立半球形振荡气室和三维数值波浪水槽模型。仿真结果表明:增设气室后壁,合理设计波道开口角度实现多向迎波捕获波浪能,优化前壁形状可降低波浪触底反射带来的能量耗散,同时提高了气室内空气压强和出气口速度,有效提升波浪能俘获效率,为后续发电的二次能量转换提供高效的空气动力。     

锦州湾表层海水微塑料分布特征

海洋微塑料是全球关注的新兴环境问题,海湾由于特殊的地理环境特征,成为微塑料分布研究的热点区域。本研究以锦州湾为研究海域,于2017年10月布设了11个点位开展表层海水微塑料样品采集,在实验室采用湿式氧化法开展样品前处理,应用傅立叶变换显微红外光谱仪分析鉴定微塑料成分。研究结果表明,锦州湾表层水体微塑料平均丰度为(0.93±0.59)个/m³,微塑料数量占全部塑料样品的96.2%。微塑料的主要成分为聚丙烯和聚乙烯,分别占55.0%和23.5%;线状和片状塑料的比例最高,分别占41.7%和26.2%;白色、蓝色和半透明微塑料分别占35.1%、26.0%和21.4%。受水动力条件和陆域河流输入等影响,锦州湾表层水体中微塑料的空间分布整体呈现北部偏高,向南部递减的趋势。     

Using stable isotopes to understand seasonal and interannual dynamics in moisture sources and atmospheric circulation in precipitation

Stable isotopes in precipitation are useful tracers to strengthen understanding of climate change and hydrological processes. In this study, the moisture sources of 190 precipitation events in Beijing were analysed using the Hybrid Single‐particle Lagrangian Integrated Trajectory model, based on which we studied the relation between variations in precipitation δ¹⁸O and dynamics in moisture sources and atmospheric circulation in seasonal and interannual timescales. Categorization of 7 groups of moisture sources was performed, among which oceanic moisture sources presented lower δ¹⁸O in precipitation than continental moisture sources. The results show that seasonal variations of precipitation δ¹⁸O were caused by changes of moisture sources. In summer, moisture from proximal oceans dominated vapour transport to Beijing due to increasing monsoon strength and resulted in a relatively small variation of precipitation δ¹⁸O. At the interannual timescale, the variations of δ¹⁸O in summer precipitation were related to dynamics in oceanic moistures, showing depleted values when the contribution of oceanic moistures, especially the proportion of long‐distance oceanic moisture, was high. Further analysis indicated that changes of oceanic moisture sources were controlled by the strength of summer monsoons. These findings address the complexity of moisture sources in midlatitude monsoon areas and suggest that isotopic signals in precipitation have the potential to deduce changes in moisture sources and atmospheric circulation and can therefore serve for palaeoclimate reconstruction.     

Why asymptotic limit of the acid dissolution capacity can lead to a sharp dissolution front in chemical dissolution of porous rocks?

The use of the asymptotic limit can greatly simplify the theoretical analysis of chemical dissolution front instabilities in fluid‐saturated rocks and therefore make it possible to obtain mathematical solutions, which often play a crucial role in understanding the propagation behavior of chemical dissolution fronts in chemical dissolution systems. However, there has been a debate in recent years that the asymptotic limit of the acid dissolution capacity (i.e., the acid dissolution capacity number approaching zero) alone cannot lead to a sharp dissolution front of the Stefan type in the acidization dissolution system, in which the dissolvable minerals of carbonate rocks are chemically dissolved by the injected acid flow. The acid dissolution capacity number is commonly defined as the ratio of the volume of the carbonate rock dissolved by an acid to that of the acid. In this paper, we use four different proof methods, including (i) direct use of the fundamental concepts; (ii) use of the mathematical governing equations of an acidization dissolution system; (iii) use of the different time scaling approach; and (iv) use of a moving coordinate system approach, to demonstrate that the asymptotic limit of the acid dissolution capacity can indeed lead to sharp dissolution fronts of the Stefan type in acidization dissolution systems on a much larger time scale (than the dissolution time scale). Our new finding is that on the reaction time scale, the condition of the conventional time derivative of porosity approaching zero alone can ensure that the acidization dissolution front has a sharp shape of the Stefan type. Copyright © 2017 John Wiley & Sons, Ltd.     

Generation of continuous surface soil moisture dataset using combined optical and thermal infrared images

Surface soil moisture (SSM) is a critical variable for understanding water and energy flux between the atmosphere and the Earth's surface. An easy to apply algorithm for deriving SSM time series that primarily uses temporal parameters derived from simulated and in situ datasets has recently been reported. This algorithm must be assessed for different biophysical and atmospheric conditions by using actual geostationary satellite images. In this study, two currently available coarse‐scale SSM datasets (microwave and reanalysis product) and aggregated in situ SSM measurements were implemented to calibrate the time‐invariable coefficients of the SSM retrieval algorithm for conditions in which conventional observations are rare. These coefficients were subsequently used to obtain SSM time series directly from Meteosat Second Generation (MSG) images over the study area of a well‐organized soil moisture network named REMEDHUS in Spain. The results show a high degree of consistency between the estimated and actual SSM time series values when using the three SSM dataset‐calibrated time‐invariable coefficients to retrieve SSM, with coefficients of determination (R²) varying from 0.304 to 0.534 and root mean square errors ranging from 0.020 m³/m³ to 0.029 m³/m³. Further evaluation with different land use types results in acceptable debiased root mean square errors between 0.021 m³/m³ and 0.048 m³/m³ when comparing the estimated MSG pixel‐scale SSM with in situ measurements. These results indicate that the investigated method is practical for deriving time‐invariable coefficients when using publicly accessed coarse‐scale SSM datasets, which is beneficial for generating continuous SSM dataset at the MSG pixel scale.