Large-Eddy Simulation of the Three-Dimensional Experiment on Richtmyer-Meshkov Instability Induced Turbulence

A program MVFT3D of large-eddy simulation is developed and performed to solve the multi compressible Navier- Stokes equations. The SGS dissipation and molecular viscosity dissipation have been analyzed, and the former is much larger than the later. Our test shows that the SGS dissipation of Vreman model is smaller than the Smagorinsky model. We mainly simulate the experiment of fluid instability of shock-accelerated interface by Poggi in this paper. The decay of the turbulent kinetic energy before the first reflected shock wave–mixing zone interaction and its strong enhancement by re-shocks are presented in our numerical simulations. The computational mixing zone width under double re-shock agreement well with the experiment, and the decaying law of the turbulent kinetic energy is consistent with Mohamed and Larue’s investigation. Also, by using MVFT3D we give some simulation results of the inverse Chevron model from AWE. The numerical simulations presented in this paper allow us to characterize and better understand the Richtmyer-Meshkov instability induced turbulence, and the code MVFT3D is validated.

Quasi-3D Numerical Simulation of Tidal Hydrodynamic Field

Based on the 2D horizontal plane numerical model, a quasi-3D numerical model is established for coastal regions of shallow water. The characteristics of this model are that the velocity profiles;can be obtained at the same time when the equations of the value of difference between the horizontal current velocity and its depth-averaged velocity in the vertical direction are solved and the results obtained are consistent with the results of the 2D, model. The circulating flow in the rectangular area induced by wind is simulated and applied to the tidal flow field of the radial sandbanks in the South Yellow Sea. The computational results from this quasi-3D model are in good agreement with analytical results and observed data. The solution of the finite difference equations has been found to be stable, and the model is simple, effective and practical.

Efficient solution of large-scale matrix of acoustic wave equations in 3D frequency domain

In 3D frequency domain seismic forward and inversion calculation,the huge amount of calculation and storage is one of the main factors that restrict the processing speed and calculation efficiency.The frequency domain finite-difference forward simulation algorithm based on the acoustic wave equation establishes a large bandwidth complex matrix according to the discretized acoustic wave equation,and then the frequency domain wave field value is obtained by solving the matrix equation.In this study,the predecessor's optimized five-point method is extended to a 3D seven-point finite-difference scheme,and then a perfectly matched layer absorbing boundary condition(PML)is added to establish the corresponding matrix equation.In order to solve the complex matrix,we transform it to the equivalent real number domain to expand the solvable range of the matrix,and establish two objective functions to transform the matrix solving problem into an optimization problem that can be solved using gradient methods,and then use conjugate gradient algorithm to solve the problem.Previous studies have shown that in the conjugate gradient algorithm,the product of the matrix and the vector is the main factor that affects the calculation efficiency.Therefore,this study proposes a method that transform bandwidth matrix and vector product problem into some equivalent vector and vector product algorithm,thereby reducing the amount of calculation and storage.

Contrast between 2D inversion and 3D inversion based on 2D high-density resistivity data

The 2D data processing adopted by the high-density resistivity method regards the geological structures as two degrees, which makes the results of the 2D data inversion only an approximate interpretation;the accuracy and effect can not meet the precise requirement of the inversion. Two typical models of the geological bodies were designed, and forward calculation was carried out using finite element method. The forward-modeled profiles were obtained. 1% Gaussian random error was added in the forward models and then 2D and 3D inversions using a high-density resistivity method were undertaken to realistically simulate field data and analyze the sensitivity of the 2D and 3D inversion algorithms to noise. Contrast between the 2D and 3D inversion results of least squares inversion shows that two inversion results of high-density resistivity method all can basically reflect the spatial position of an anomalous body. However, the 3D inversion can more effectively eliminate the influence of interference from Gaussian random error and better reflect the distribution of resistivity in the anomalous bodies. Overall, the 3D inversion was better than 2D inversion in terms of embodying anomalous body positions, morphology and resistivity properties.

Clinical application of improved 2D computer-assisted fluoroscopic navigation through simulating a 3D vertebrae image to guide pedicle screw internal fixation

Objective To study the effect of using improved 2D computer-assisted fluoroscopic navigation through simulating 3D vertebrae image to guide pedicle screw internal fixation.Methods Posterior pedicle screw internal fixation,distraction

2D-3D registration for 3D analysis of lower limb alignment in a weight-bearing condition

X-ray imaging is the conventional method for diagnosing the orthopedic condition of a patient. Computerized Tomography(CT) scanning is another diagnostic method that provides patient’s 3D anatomical information. However, both methods have limitations when diagnosing the whole leg; X-ray imaging does not provide 3D information, and normal CT scanning cannot be performed with a standing posture. Obtaining 3D data regarding the whole leg in a standing posture is clinically important because it enables 3D analysis in the weight bearing condition.Based on these clinical needs, a hardware-based bi-plane X-ray imaging system has been developed; it uses two orthogonal X-ray images. However, such methods have not been made available in general clinics because of the hight cost. Therefore, we proposed a widely adaptive method for 2 D X-ray image and 3D CT scan data. By this method, it is possible to threedimensionally analyze the whole leg in standing posture. The optimal position that generates the most similar image is the captured X-ray image. The algorithm verifies the similarity using the performance of the proposed method by simulation-based experiments. Then, we analyzed the internal-external rotation angle of the femur using real patient data. Approximately 10.55 degrees of internal rotations were found relative to the defined anterior-posterior direction. In this paper, we present a useful registration method using the conventional X-ray image and 3D CT scan data to analyze the whole leg in the weight-bearing condition.

Multi-relaxation-time lattice Boltzmann simulations of lid driven flows using graphics processing unit

Large eddy simulation （LES） using the Smagorinsky eddy viscosity model is added to the two-dimensional nine velocity components （D2Q9） lattice Boltzmann equation （LBE） with multi-relaxation-time （MRT） to simulate incompressible turbulent cavity flows with the Reynolds numbers up to 1 × 10^7. To improve the computation efficiency of LBM on the numerical simulations of turbulent flows, the massively parallel computing power from a graphic processing unit （GPU） with a computing unified device architecture （CUDA） is introduced into the MRT-LBE-LES model. The model performs well, compared with the results from others, with an increase of 76 times in computation efficiency. It appears that the higher the Reynolds numbers is, the smaller the Smagorinsky constant should be, if the lattice number is fixed. Also, for a selected high Reynolds number and a selected proper Smagorinsky constant, there is a minimum requirement for the lattice number so that the Smagorinsky eddy viscosity will not be excessively large.

2D Coupled 3D:A New Numerical Model for Dual - Structured - Aquifer System

This paper points out that a successful numerical simulation is to construct a correct conceptional model which is very dose to the natural condition. A new model, two dimensional coupled three dimensional model (2D -3D ) is presented in the Present paper,which is the most suitable one for the dual - structured - aquifer system. An example of Wenyinghu area is shown.By using the 2D-3D model, a satisfied result of the simulated area is achieved.

Numerical modeling of 3D DC resistivity method in the mixed space-wavenumber domain

Forward modeling is the basis of inversion imaging and quantitative interpretation for DC resistivity exploration.Currently,a numerical model of the DC resistivity method must be finely divided to obtain a highly accurate solution under complex conditions,resulting in a long calculation time and large storage.Therefore,we propose a 3D numerical simulation method in a mixed space-wavenumber domain to overcome this challenge.The partial differential equation about abnormal potential is transformed into many independent ordinary differential equations with different wavenumbers using a 2D Fourier transform along the x axis and y axis direction.In this way,a large-scale 3D numerical simulation problem is decomposed into several 1D numerical simulation problems,which significantly reduces the computational and storage requirements.In addition,these ordinary 1D differential equations with different wavenumbers are independent of each other and high parallelelism of the algorithm.They are solved using a finite-element algorithm combined with a chasing method,and the obtained solution is modified using a contraction operator.In this method,the vertical direction is reserved as the spatial domain,then grid size can be determined flexibly based on the underground current density distribution,which considers the solution accuracy and calculation efficiency.In addition,for the first time,we use the contraction operator in the integral equation method to iterate the algorithm.The algorithm takes advantage of the high efficiency of the standard Fourier transform and chasing method,as well as the fast convergence of the contraction operator.We verified the accuracy of the algorithm and the convergence of the contraction operator.Compared with a volume integral method and goal-oriented adaptive finite-element method,the proposed algorithm has lower memory requirements and high computational efficiency,making it suitable for calculating a model with large-scale nodes.Moreover,different examples are used to verify the high adaptability and parallelism of the proposed algorithm.The findings show that the 3D numerical simulation method of DC resistivity method in a mixed space-wavenumber domain is highly efficient,precise,and parallel.

Assessment of Turbulence Models on Bridge-Pier Scour Using Flow-3D

The main purpose of this paper is to analyze the influence of different turbulence flow models on scouring pit of bridge-pier. Flow-3D software is applied in line with the purpose. The key motivation for this study is to contribute to the Flow-3D software by means of some modification and adjustment in the sediment scour model and shallow water model. An assessment of turbulence model adopted with the parameters of the Melville experiment to estimate the maximum scour-depth was performed. In the simulation results, the alternate eddy formation and shedding were repeated while the Karman vortex street formed behind the pier for the large eddy simulation LES turbulence model is more realistic in the flow phenomenon. The results of the scour development of large eddy simulation (LES) turbulence model were found to be more satisfied than the Renormalized group (RNG) turbulence model and close to the prior experiment results. The simulated scour results were significantly different with the observed data collected from previous literature in the reason of some unsuitability of meshing method in Flow-3D software.

LARGE EDDY SIMULATION OF VERTICAL JET IMPINGEMENT WITH A FREE SURFACE

The flow characteristics of a 2D slot jet vertically impinging on free surfaces are numerically investigated by the Large Eddy Simulation （LES） with a dynamic Sub-Grid Scale （SGS） model. The σ-coordinate transformation is introduced to map the depth-variable physical domain to a depth-uniform computational domain. The split-operator scheme, which splits the solution procedure into advection, diffusion and pressure propagation steps, is employed to solve the instantaneous velocity and pressure field. A fully nonlinear Lagrange-Euler method is used to compute the free surface elevation. The numerical results show that the jet retains good self-similarity in the Zone of Eestablished Flow （ZEF）, while the turbulence of jet shifts gradually from isotropic to anisotropic in the Zone of Surface Impingement （ZSI）. When the jet approaches the free surface, the centerline velocity decreases rapidly and the flow is deflected by the free surface. Two symmetrical surface jets are formed in the Zone of Horizontal Jets （ZHJ）. An approximate Guassian distribution of super-elevation is also formed on the free surface. The computed results are in good agreement with the available experimental data.

126kV真空灭弧室1/3线圈纵磁触头三维涡流场仿真

建立了126kV、1/3线圈纵磁触头三维电极结构模型,并利用有限元法对考虑了电弧影响的模型进行了三维涡流场仿真.在电流分别处于峰值和过零时,得到了纵向磁场分别在触头片表面和触头间隙中心平面上的分布,以及纵向磁场滞后时间沿径向路径的分布规律.

Prognostic analysis of prostaglandin D2 synthase in diffuse large B‑cell lymphoma

Purpose We aimed to analyse the correlation between prostaglandin D2 synthase(PTGDS)and a diffuse large B-cell lymphoma(DLBCL)prognosis.Methods We retrospectively collected two hundred paraffin-embedded tissue specimens that were pathologically diagnosed as DLBCL in the Fujian Tumour Hospital between January 2014 and December 2018.An abundance of paraffin-embedded tumour tissues were obtained.Twenty patients with lymphocyte-rich,benign,tissue-reactive,hypertrophic tonsillitis were selected as controls.Wax blocks were selected for primary cases and the controls were screened by professional pathologists.The levels of prostaglandin D2 synthase(PTGDS)and the EMT-related molecules,E-cadherin and vimentin,were detected by immunohistochemistry in clinical samples.A chi-square test revealed the correlations between PTGDS expression and clinicopathological characteristics,including age,sex,primary site,clinical stage,immunotyping,and International Prognostic Index(IPI)score.The Kaplan–Meier survival analysis was performed,and the diagnostic value was evaluated by receiver operating characteristic(ROC)curve analysis.Results A total of 138 cases(69%)were found to be PTGDS positive(>30%positive cells).PTGDS staining was negative(<30%positive cells)in 62 cases(31%).We collected the corresponding clinicopathological information and found that PTGDS expression was not significantly related to the patients’age,tumour stage,presence of extranodal invasion,or IPI score.According to the follow-up data,patients with low PTGDS expression had poor progression-free survival(PFS)and overall survival(OS),with 2-year PFS and OS rates of 41.7%and 50%,respectively.The 2-year PFS and OS rates of PTGDS-positive patients were 89.3%and 92.9%,respectively(P<0.0001),and the differences were significant.Conclusion We found that the expression level of PTGDS is significantly correlated with the prognosis of diffuse large B-cell lymphoma.

合成射流激励频率对D型体流动减阻控制的影响

采用合成射流作为主动控制手段,对自由空间中D型体尾流的控制机理和气动减阻进行了数值仿真.合成射流布置在D型体垂直背上下缘处并同时作用,射流方向与来流方向平行,频率从40 Hz变化到500 Hz.结果表明,在整个频率范围内都实现了减阻.在低频激励(65 Hz)作用下,尾流上下两侧剪切层趋于同步化,汇合位置延后,卡门涡街强度减弱;在高频激励(500 Hz)作用下,伴随着整流,剪切层向内偏转,尾流夹带现象减弱,回流区略延长.卡门涡街的抑制对于钝体尾流的形成以及减阻效果有着显著的影响.

A HYDRODYNAMIC LARGE-EDDY SIMULATION-MODEL USING DIRECTIONAL EDDY VISCOSITY METHOD IN VERTICAL DIRECTION

A non-linear turbulent hydrodynamic model was proposed by Zhan and Li (1993) using the large eddy simulation approach and directional eddy viscosity method. The model is simplified in this paper and is applied to the problem of a side discharge into open-channel flow,where a recirculation zone develops downstream of the discharge. The numerical results are compared with experimental results (Strazisar and Prahl (1973) and Mikhail, Chu, and Savage (1975)) and previous numerical results obtained using a k-model (McGuirk and Rodi, 1978). The agreement is good.

3D reconstruction of coal pore network and its application in CO_(2)-ECBM process simulation at laboratory scale

Three-dimensional(3D)reconstruction of the equivalent pore network model(PNM)using X-ray computed tomography(CT)data are of significance for studying the CO_(2)-enhanced coalbed methane recovery(CO_(2)-ECBM).The docking among X-ray CT technology,MATLAB,with COMSOL software not only can realize the 3D reconstruction of PNM,but also the CO_(2)-ECBM process simulation.The results show that the Median filtering algorithm enabled the de-noising of the original 2D CT slices,the image segmentation of all slices was realized based on the selected threshold,and the PNM can be constructed based on the Maximum Sphere algorithm.The mathematical model of CO_(2)-ECBM process fully coupled the expanded Langmuir equation.At the same time for CO_(2)injection,CH_(4)pressure tends to decrease with the increase of CO_(2)pressure,but its difference is not obvious.The CH_(4)pressure in the slice center changed a lot,while at the edge it changed a little under different CO_(2)pressures.The injected CO_(2)was transported to matrix along the macro and micro-fractures with continuous flow.The injected CO_(2)first replaced the adsorbed CH_(4)by covering the inner surface of macro-pores and meso-pores to form the single molecular layer adsorption of CO_(2).Then they migrated to micro-pores by Fick’s diffusion,sliding flow,and surface diffusion.Furthermore,the CO_(2)replaced CH_(4)adsorbed by volumetric filling in micro-pores,and formed the multi-molecular layer adsorption of CO_(2).The gas pressure and migration path between CO_(2)and CH_(4)are opposite.This study can provide a theoretical basis for studying digital rock physics technology and enrich the development of CO_(2)-ECBM technology.

竖向TMD用板式电涡流阻尼器磁路对比分析

为提升竖向TMD用板式电涡流阻尼器(PECD)的耗能效率,分别开展了典型竖向电涡流TMD样机的阻尼性能测试与有限元仿真分析,基于相邻永磁体有效耦合率性能指标,采用三维电磁场有限元分析方法分别优化了外置式与内置式PECD的磁路设计,并对比分析了各磁路的耗能效率、耐久性与占用空间等。研究结果表明:提出的相邻永磁体有效耦合率可以很好地表征永磁体相互作用对PECD等效阻尼系数的影响程度,有助于指导磁路优化设计;外置式与内置式PECD磁路优化设计方案在相邻永磁体磁极布置、水平与竖向适宜间距等方面均存在较大差别;各最优磁路永磁体用量一定时,内置式PECD的等效阻尼系数略低于外置式PECD(附加导磁钢板),但远高于外置式PECD(无导磁钢板);装配内置式PECD的竖向TMD结构相对紧凑,且耐久性也相对较高。

Numerical simulation of tomographic-SAR imaging and object reconstruction using compressive sensing with L_(1/2)-norm regularization

By making use of multiple acquisitions of synthetic aperture radar(SAR) observations over the same area, tomographic-SAR(tomo-SAR) technology can achieve three-dimensional(3-D) imaging of the objects of interest. The compressive sensing(CS) approach has been applied to deal with the sparseness of the elevation signals.Due to its sparsity and convexity, the L1-norm regularization, as an approximated L0-norm with an exact solution,has been employed in CS to reconstruct the reflectivity profile of the objects. In this paper, based on our studies on polarimetric scattering and SAR imaging simulations, we produce numerical multi-pass tomo-SAR observations of the terrain object. Then, we present the CS with novel L1/2-norm regularization to realize 3-D reconstruction. As a non-convex optimization problem, the L1/2-norm regularization is solved by an iterative algorithm. This numerical simulation of tomo-SAR imaging and 3-D reconstruction of the object modeling can be of great help for parameterized analysis of tomo-SAR imagery. As an example, a tomo-SAR image and 3-D reconstruction of the Beijing National Stadium model are presented.

126kV真空灭弧室1/2线圈纵磁触头三维磁场仿真

建立了126kV、1/2线圈纵磁触头三维电极结构模型,并利用有限元法对考虑了电弧影响的有限元模型进行了三维静磁场和涡流场仿真。结果表明静磁场下纵向磁场在触头片开槽处发生畸变,在触头间隙中心平面分布呈“平顶帐篷”形状。在涡流场下,当电流处于峰值时,纵向磁场在触头片上的分布和触头间隙中心平面上的分布与静磁场作用下相似:在电流过零时,纵向磁场在触头片中心形成一个最高峰,在最高峰周围形成6个次高峰;而在触头间隙中心平面分布呈“尖顶帐篷”形;沿路径(0,-60,110)～(0,60,110)纵向磁场的分布为开口向上的抛物线型,最大滞后时间0.36ms。

鄂尔多斯盆地长北气田山西组2段高分辨率层序构型与砂体预测

以钻井测井、岩芯描述和样品分析为基础,高分辨率层序分析为技术方法对鄂尔多斯盆地长北气田山西组2段进行了精细的沉积相和层序分析,将山2段确定为海相辫状河三角洲沉积体系,并将其划分为1个长期、2个中期、4个短期和9个超短期基准面旋回层序和2大类5种基准面旋回结构类型。分析了不同级别基准面旋回层序的叠加式样与砂体展布和演化规律,建立了以中期基准面升、降半旋回相域为等时地层单元的山2段地层格架和对单砂层进行了等时追踪对比,确定中期上升期是三角洲河道砂体连续叠置形成巨厚砂体的最有利时期,洪泛期为细粒沉积物沉积期,也是隔层和局域性盖层主要发育期,下降期受侵蚀影响河道砂体不发育。以高分辨率层序分析为基础,建立等时地层格架为约束条件,结合单砂体为等时地层单元编制砂岩等厚图、砂/泥比等值线图和沉积微相图,对长北气田山2段砂体分布开展序贯指示三维模拟和预测,为长北气田山2段气藏高效开发提供依据。