Numerical simulation of a stratigraphic model

A multi-lithology diffusive stratigraphic model is considered, which simulates at large scales in space and time the infill of sedimentary basins governed by the interaction between tectonics displacements, eustatic variations, sediment supply, and sediment transport laws. The model accounts for the mass conservation of each sediment lithology resulting in a mixed parabolic, hyperbolic system of partial differential equations (PDEs) for the lithology concentrations and the sediment thickness. It also takes into account a limit on the rock alteration velocity modeled as a unilaterality constraint. To obtain a robust, fast, and accurate simulation, fully and semi-implicit finite volume discre tization schemes are derived for which the existence of stable solutions is proved. Then, the set of nonlinear equations is solved using a Newton algorithm adapted to the unilaterality constraint, and preconditioning strategies are defined for the solution of the linear system at each Newton iteration. They are based on an algebraic approximate decoupling of the sediment thickness and the concentration variables as well as on a proper preconditioning of each variable. These algorithms are studied and compared in terms of robustness, scalability, and efficiency on two real basin test cases.     

风机基础TLP平台在波浪中运动性能的CFD数值分析

随着海上风能的开发向深水发展,支撑风机的载体平台越来越受到关注。在经济性与安全性、稳定性的多重要求下,张力腿平台（TLP）在海洋风能资源的开发中体现出了重要地位。采用基于开源平台OpenFOAM开发的计算流体动力学（CFD）水动力学求解器naoe-FOAM-SJTU对一座处于中等水深下的风机基础水下TLP（STLP）的运动响应进行了数值模拟与研究。文中使用弹簧锚链模型模拟STLP的垂向系泊锁链系统,模拟该平台在不同波浪环境下的运动响应情况。首先将STLP单自由度自由衰减CFD模拟结果与已有全耦合时域分析结果进行对比,验证了naoe-FOAM-SJTU求解器及使用弹簧模型模拟STLP系泊系统的准确性与可靠性。随后在考虑非线性波浪载荷的情况下研究极端海况下与一般作业海况下STLP的运动响应情况,计算工况中的风机基础所受弯矩及锚链受力情况,并详细展示流场、速度场信息,分析高阶波浪成分、不同海况等条件对于STLP运动性能的影响。研究结果表明,TLP在中等水深中具有良好的运动性能,naoe-FOAM-SJTU求解器可以有效模拟水中生产平台在波浪环境下的水动力问题,并可以对整个流场进行可视化展示与分析。     

New finite-gap solutions for the coupled Burgers equations engendered by the Neumann systems

On the tangent bundle TS^N-1 of the unit sphere S^N-1, this paper reduces the coupled Burgers equations to two Neumann systems by using the nonlinearization of the Lax pair, whose Liouville integrability is displayed in the scheme of the r-matrix technique. Based on the Lax matrix of the Neumann systems, the Abel--Jacobi coordinates are appropriately chosen to straighten out the restricted Neumann flows on the complex torus, from which the new finite-gap solutions expressed by Riemann theta functions for the coupled Burgers equations are given in view of the Jacobi inversion.     

Investigation of highly efficient algorithms for solving linear equations in the discontinuous deformation analysis method

Large‐scale engineering computing using the discontinuous deformation analysis (DDA) method is time‐consuming, which hinders the application of the DDA method. The simulation result of a typical numerical example indicates that the linear equation solver is a key factor that affects the efficiency of the DDA method. In this paper, highly efficient algorithms for solving linear equations are investigated, and two modifications of the DDA programme are presented. The first modification is a linear equation solver with high efficiency. The block Jacobi (BJ) iterative method and the block conjugate gradient with Jacobi pre‐processing (Jacobi‐PCG) iterative method are introduced, and the key operations are detailed, including the matrix‐vector product and the diagonal matrix inversion. Another modification consists of a parallel linear equation solver, which is separately constructed based on the multi‐thread and CPU‐GPU heterogeneous platforms with OpenMP and CUDA, respectively. The simulation results from several numerical examples using the modified DDA programme demonstrate that the Jacobi‐PCG is a better iterative method for large‐scale engineering computing and that adoptive parallel strategies can greatly enhance computational efficiency. Copyright © 2015 John Wiley & Sons, Ltd.     

非结构网格上的三维浅水流动数值模型

针对当前复杂环境水流模拟的需求,建立了新型的基于特征型高分辨率数值算法的三维非结构网格浅水动力模型。模型采用有限体积法离散sigma坐标下的三维浅水方程,运用Roe黎曼近似解评估水平界面通量。模型网格拟合边界能力强,可根据需要局部加密;格式数值性能优良,具有守恒性、单调迎风性、高数值分辨率等特性。同时,应用干湿判别法处理动边界,以适应浅滩地形漫/露过程模拟的需要。封闭水池内部风生环流、干河床上溃坝过程和长江口实际潮流场的模拟从不同侧面展示了模型的特点,结果表明它能够准确地预测水流的三维流动结构,而且计算简单高效,具有良好的数值稳定性。     

Numerical and Experimental Analysis of A Ducted Propeller Designed by A Fully Automated Optimization Process Under Open Water Condition

A fully automated optimization process is provided for the design of ducted propellers under open water conditions, including 3D geometry modeling, meshing, optimization algorithm and CFD analysis techniques. The developed process allows the direct integration of a RANSE solver in the design stage. A practical ducted propeller design case study is carried out for validation. Numerical simulations and open water tests are fulfilled and proved that the optimum ducted propeller improves hydrodynamic performance as predicted.     

Optimizing dewatering design for a metro station on the Chengdu Metro Line 7

At present, most calculation results regarding foundation pit dewatering are ideal values, making construction resources prone to being wasted. In order to optimize the traditional pipe well design of large wells, the linear programming solution module in Excel is used, with the total water inflow taken as the objective function, the water level drawdown used as the constraint and test condition, and a station project on the Chengdu Metro Line 7 serving as the subject of this study. The total water inflow of the traditional pipe well design is optimized by the simplex method, producing a total water inflow of 4 040.65 m~3/d, which, compared with 4 829.79 m~3/d, the total water inflow calculated by means of the traditional design optimization method, engenders a reduction of roughly 16% per day. The feasibility of the optimization methodology is verified by the drawdown constraint, which reveals the decrease of construction costs and the diminution of the influence that the lowered groundwater level has on the surroundings of the metro station. Finally, references are provided as to optimizing the dewatering designs for other metro stations in similar engineering and hydrogeological conditions.     

基于二次场的可控源电磁法三维有限元-无限元数值模拟

本文提出了一种新的混合有限元-无限元三维可控源电磁法(CSEM)问题快速高精度正演模拟算法.首先从电场双旋度方程出发,推导了水平电偶极子源的二次场边值问题,采用无限元代替截断边界条件和有限元离散内部计算区域的新策略,达到减小计算区域的目的,基于并行直接求解技术,实现多源CSEM问题的快速精确求解.其次,通过层状解析模型测试,一方面验证了新算法的正确性,另一方面通过与其他三种已知CSEM问题求解策略进行对比,表明了本文提出的基于二次场有限元-无限元算法具有离散区域小、求解速度快和计算精度高等优点.最后,通过3D模型计算,清晰直观地模拟了场源阴影效应,为野外数据的处理与解释提供指导.     

扎鲁特地区无烟煤分子结构特征和模型构建

为了构建内蒙古扎鲁特地区无烟煤的分子结构模型,采用工业分析、元素分析、核磁共振碳谱、高分辨率透射电镜和X射线光电子能谱等测试方法对该地区无烟煤的分子结构进行研究。结果表明无烟煤分子结构中以芳香碳为主,脂肪碳以短烷基侧链和环烷烃为主。煤分子结构中氧原子主要以酚羟基和醚氧基的形式存在,氮原子主要以吡咯型氮和吡啶型氮的形式存在,硫原子主要的赋存形态是硫醇硫酚。结合由核磁共振碳谱测得的煤结构参数和由高分辨率透射电镜测得的芳香环尺寸和含量,构建了无烟煤大分子平均结构模型。对构建的模型进行了结构优化与能量分析,发现煤分子结构的芳香片层趋于平行排列,非键结势能中的范德华力是保持煤结构稳定的主要因素。本研究为采用分子动力学模拟从分子尺度研究扎鲁特地区无烟煤在石墨化过程中,其芳香碳层的拼叠过程及其反应路径提供了模型基础。      

The chemical changes of DOM from black waters to coastal marine waters by HPLC combined with ultrahigh resolution mass spectrometry

How dissolved organic matter (DOM) undergoes chemical changes during its transit from river to ocean remains a challenge due to its complex structure. In this study, DOM along a river transect from black waters to marine waters is characterized using an offline combination of reversed-phase high performance liquid chromatography (RP-HPLC) coupled to electrospray ionization Fourier transform ion cyclotron resonance mass spectrometry (ESI-FTICR-MS), as well as tandem ESI-FTICR-MS. In addition, a water extract from degraded wood that mainly consists of lignins is used for comparison to the DOM from this transect. The HPLC chromatograms of all DOM samples and the wood extract show two major well-separated components; one is hydrophilic and the other is hydrophobic, based on their elution order from the C₁₈ column. From the FTICR-MS analysis of the HPLC fractions, the hydrophilic components mainly contain low molecular weight compounds (less than 400 Da), while the hydrophobic fractions contain the vast majority of compounds of the bulk C₁₈ extracted DOM. The wood extract and the DOM samples from the transect of black waters to coastal marine waters show strikingly similar HPLC chromatograms, and the FTICR-MS analysis further indicates that a large fraction of molecular formulas from these samples are the same, existing as lignin-like compounds. Tandem mass spectrometry experiments show that several representative molecules from the lignin-like compounds have similar functional group losses and fragmentation patterns, consistent with modified lignin structural entities in the wood extract and these DOM samples. Taken together, these data suggest that lignin-derived compounds may survive the transit from the river to the coastal ocean and can accumulate there because of their refractory nature.