基于mixed Lagrangian/Eulerian方法的轮轨滚动接触特性分析

为了缩短有限元方法显式计算轮轨滚动接触的时间,采用mixed Lagrangian/Eulerian方法建立了轮轨滚动接触有限元模型.应用该模型对轮轨接触区的单元进行细化,计算了列车在启动、运行和制动工况下轮轨的接触特性.计算结果表明:不同工况下,轮轨滚动接触区最大Mises应力、最大接触应力和接触斑面积等法向特性变化幅度均在2%以内,但接触区纵向截面中Mises应力分布及纵向剪切应力分布有较大变化;启动和制动工况下,最大Mises应力和最大纵向剪切应力位置均比自由滚动时更接近于轮轨表面;不同工况下,摩擦力大小和方向发生变化,在列车牵引和制动工况中,摩擦力达到极限时轮轨间出现完全滑动,摩擦力方向与滑动方向相反,且不同速度等级下的纵向摩擦力变化幅度也在2%以内.

高温气流内顺/逆流喷雾的Eulerian-Lagrangian数值仿真(英文)

用Eulerian-Lagrangian方法对高温气流内顺/逆流喷雾过程进行了数值仿真。气相方程用Eulerian方法建立,用k-ε模型封闭;液滴相方程用Lagrangian方法建立,采用了随机轨道模型;考虑了气体和液滴间的双向耦合。结果表明,高温气流内喷雾后,气相流场内产生了速度亏损和温度亏损,气体出口压力变化不大,液滴运动决定于气流湍流强度和自身惯性力,液滴在气道内的停留时间决定于液滴运动特性和气流湍流强度;逆流喷雾对气流的影响比顺流喷雾时大。仿真结果与文献结果基本一致。

基于level set的Eulerian-Lagrangian耦合方法及其应用

提出了一种基于水平集的Eulerian-Lagrangian耦合方法,其中Lagrangian方法采用相容显式有限元拉氏方法,Eulerian方法采用基于近似Riemann解的有限体积Eulerian方法,多介质界面处理采用新的水平集和Ghost方法计算.给出了若干数值算例,包括激波管问题以及金属和气体的运动界面及其大变形问题,并分别与精确解和相容显式有限元拉氏方法的计算结果进行了对比.数值结果表明,该方法计算结果正确,精度较高,能够准确捕捉物质界面,适用于处理大变形问题.

基于Lagrangian-Eulerian耦合的放射性物质大气扩散模型研究

针对中尺度扩散下拉格朗日模型和欧拉模型在计算效率和准确度方面的矛盾,建立了适用于模拟核事故下放射性物质大气扩散的基于拉格朗日模型与欧拉模型耦合的预测模型,以国内某核电站为例进行了放射性物质扩散的虚拟算例模拟,对模型进行验证和分析。结果表明:拉格朗日-欧拉耦合模型可以较好地模拟放射性物质在大气中的扩散,不仅可以提升中尺度和大尺度下的计算效率,还可以有效减小欧拉模型初始条件所带来的系统误差。

Efficient modeling of cable-pulley system with friction based on arbitrary-Lagrangian-Eulerian approach

In conventional modeling of a cable-pulley system, the cable must be finely meshed with Lagrangian elements for valid contact detections with pulleys, leading to extremely low efficiency. The sliding joint method based on the arbitrary-Lagrangian- Eulerian （ALE） formulation still lacks an efficient cable element, and in particular, modeling of friction between a sliding joint and the cable has not been studied. This paper presents efficient multi-body modeling of a cable-pulley system with friction. A variable- length cable element with a node movable along the cable, which is described with ALE, is developed to mesh the cable. A transitional cable element is then proposed to model the contact part of the cable by fixing its two nodes to the two corresponding locations of the pulley. Friction of the cable-pulley is derived as a simple law of tension decay and embedded in the multi-body system modeling. It is simplified as a generalized friction force acting only on the arc-length coordinate. This approach can use a rough mesh on the cable, and is free of contact detections, thus significantly saving computation time. Several examples are presented to validate the proposed method, and show its effectiveness in real engineering applications.

A hybrid Lagrangian-Eulerian numerical model for sea-ice dynamics

A hybrid Lagrangian - Eulerian （HLE） method is developed for sea ice dynamics, which combines the high computational efficiency of finite difference method （FDM） with the high numerical accuracy of smoothed particle hydrodynamics （SPH）. In this HLE model, the sea ice cover is represented by a group of Lagrangian ice particles with their own thicknesses and concentrations. These ice variables are interpolated to the Eularian gird nodes using the Gaussian interpolation function. The FDM is used to determine the ice velocities at Eulerian grid nodes, and the velocities of Lagrangian ice particles are interpolated from these grid velocities with the Gaussian function also. The thicknesses and concentrations of ice particles are determined based on their new locations. With the HLE numerical model, the ice ridging process in a rectangular basin is simulated, and the simulated results are validated with the analytical solution. This method is also applied to the simulation of sea ice dynamics in a vortex wind field. At last, this HLE model is applied to the Bohai Sea, and the simulated concentration, thickness and velocity match the satellite images and the field observed data well.

New Eulerian-Lagrangian Method for Salinity Calculation

A difference scheme in curvilinear coordinates is put forward for calculation of salinity in estuaries and coastal waters, which is based on Eulerian-Lagrangian method. It combines first-order and second-order Lagrangian interpolation to reduce numerical dispersion and oscillation. And the length of the curvilinear grid is also considered in the interpolation. Then the scheme is used in estuary, coast and ocean model, and several numerical experiments for the Yangtze Estuary and the Hangzhou Bay are conducted to test it. These experiments show that it is suitable for simulations of salinity in estuaries and coastal waters with the models using curvilinear coordinates.

Simulation of sheet metal extrusion processes with Arbitrary Lagrangian-Eulerian method

An Arbitrary Lagrangian-Eulerian(ALE) method was employed to simulate the sheet metal extrusion process,aiming at avoiding mesh distortion and improving the computational accuracy.The method was implemented based on MSC/MARC by using a fractional step method,i.e.a Lagrangian step followed by an Euler step.The Lagrangian step was a pure updated Lagrangian calculation and the Euler step was performed using mesh smoothing and remapping scheme.Due to the extreme distortion of deformation domain,it was almost impossible to complete the whole simulation with only one mesh topology.Therefore,global remeshing combined with the ALE method was used in the simulation work.Based on the numerical model of the process,some deformation features of the sheet metal extrusion process,such as distribution of localized equivalent plastic strain,and shrinkage cavity,were revealed.Furthermore,the differences between conventional extrusion and sheet metal extrusion process were also analyzed.

Comparisons of some difference forms for compressible flow in cylindrical geometry on arbitrary Lagrangian and Eulerian framework

The study of cylindrically symmetric compressible fluid is interesting from both theoretical and numerical points of view. In this paper, the typical spherical sym- metry properties of the numerical schemes are discussed, and an area weighted scheme is extended from a Lagrangian method to an arbitrary Lagrangian and Eulerian （ALE） method. Numerical results are presented to compare three discrete configurations, i.e., the control volume scheme, the area weighted scheme, and the plane scheme with the addition of a geometrical source. The fact that the singularity arises from the geometri- cal source term in the plane scheme is illustrated. A suggestion for choosing the discrete formulation is given when the strong shock wave problems are simulated.

Hybrid N-order Lagrangian Interpolation Eulerian-Lagrangian Method for Salinity Calculation

The Eulerian？Lagrangian method（ELM） has been used by many ocean models as the solution of the advection equation,but the numerical error caused by interpolation imposes restriction on its accuracy.In the present study,hybrid N-order Lagrangian interpolation ELM（Li ELM） is put forward in which the N-order Lagrangian interpolation is used at first,then the lower order Lagrangian interpolation is applied in the points where the interpolation results are abnormally higher or lower.The calculation results of a step-shaped salinity advection model are analyzed,which show that higher order（N=3？8） Li ELM can reduce the mean numerical error of salinity calculation,but the numerical oscillation error is still significant.Even number order Li ELM makes larger numerical oscillation error than its adjacent odd number order Li ELM.Hybrid N-order Li ELM can remove numerical oscillation,and it significantly reduces the mean numerical error when N is even and the current is in fixed direction,while it makes less effect on mean numerical error when N is odd or the current direction changes periodically.Hybrid odd number order Li ELM makes less mean numerical error than its adjacent even number order Li ELM when the current is in the fixed direction,while the mean numerical error decreases as N increases when the current direction changes periodically,so odd number of N may be better for application.Among various types of Hybrid N-order Li ELM,the scheme reducing N-order directly to 1st-order may be the optimal for synthetic selection of accuracy and computational efficiency.

FINITE ELEMENT ANALYSIS FOR CHIP FORMATION IN HIGH SPEED TURNING OPERATIONS BY ARBITRARY LAGRANGIAN EULERIAN METHOD

A two-dimensional finite element （FE） model for the high speed turning operations when orthogonally machining AISI H13 tool steel at 49HRC using poly crystalline cubic boron nitride （PCBN） is described. An arbitrary Lagrangian Eulerian （ALE） method has been adopted which does not need any chip separation criteria as opposed to the traditional Lagrangian approach. Through FE simulations temperature and stresses distributions are presented that could be helpful in predicting tool life and improving process parameters. The results show that high temperatures are generated along the tool rake face as compared to the shear zone temperatures due to high thermal conductivity of PCBN tools.

A multifractal model for linking Lagrangian and Eulerian velocity structure functions

A multifractal model is developed to connect the Lagrangian multifractal dimensions with their Eulerian counterparts. We propose that the characteristic time scale of a Lagrangian quantity should be the Lagrangian time scale, and it should not be the Eulerian time scale which was widely used in previous studies on Lagrangian statistics. Using the present model, we can obtain the scaling exponents of Lagrangian velocity structure functions from the existing data or models of scaling exponents of Eulerian velocity structure functions. This model is validated by comparing its prediction with the results of experiments, direct numerical simulations, and the previous theoretical models. The comparison shows that the proposed model can better predict the scaling exponents of Lagrangian velocity structure functions, especially for orders larger than 6.

COMPARISON OF THE EULERIAN AND LAGRANGIAN TIDAL RESIDUALS IN THE BOHAI SEA

Tidal residual is very important to the transport of water particles, nutrients, plankton, etc. in the coastal sea. Eulerian scheme and Lagrangian scheme are two different ways to get the time averaged residual. Solution of the Bohai Sea’s hydrodynamic system using a semi implicit layer averaged numerical model yielded different direction Eulerian and Lagrangian tidal residuals. The latter were stronger than the former in most sea areas. Their different directions produced different circulation pattern in some areas. Compared with the Eulerian residual, the Lagrangian residual seemed to be more in accord with the observation.

A hybrid subcell-remapping algorithm for staggered multi-material arbitrary Lagrangian-Eulerian methods

A new flux-based hybrid subcell-remapping algorithm for staggered multimaterial arbitrary Lagrangian-Eulerian (MMALE) methods is presented. This new method is an effective generalization of the original subcell-remapping method to the multi-material regime (LOUBERE, R. and SHASHKOV,M. A subcell remapping method on staggered polygonal grids for arbitrary-Lagrangian-Eulerian methods. Journal of Computational Physics, 209, 105–138 (2005)). A complete remapping procedure of all fluid quantities is described detailedly in this paper. In the pure material regions, remapping of mass and internal energy is performed by using the original subcell-remapping method. In the regions near the material interfaces, remapping of mass and internal energy is performed with the intersection-based fluxes where intersections are performed between the swept regions and pure material polygons in the Lagrangian mesh, and an approximate approach is then introduced for constructing the subcell mass fluxes. In remapping of the subcell momentum, the mass fluxes are used to construct the momentum fluxes by multiplying a reconstructed velocity in the swept region. The nodal velocity is then conservatively recovered. Some numerical examples simulated in the full MMALE regime and several purely cyclic remapping examples are presented to prove the properties of the remapping method.

基于耦合Lagrangian-Eulerian算法的流体机械数值模拟

采用耦合的Lagrangian-Eulerian算法,对流体机械旋转工况进行数值模拟,计算结果表明,液位在叶片扰动之后的瞬时分布与实际情况基本相符,可以准确的计算流动对叶片的作用压力,可同时用于校核旋转工况下的叶片强度。

Coupling of 3D Eulerian and Lagrangian Spray Approaches in Industrial Combustion Engine Simulations

The Lagrangian DDM （discrete droplet model） is state-of-the-art for CFD （computational fluid dynamics） simulations of mixture formation and combustion in industrial engines. A commonly known drawback of the DDM approach is the attenuated validity in the dense spray, where the bulk liquid disintegrates into droplets. There the assumption of single droplets surrounded by a homogenous gas field is not reasonable. In this region, the Eulerian-Eulerian multi-phase approach performs better because instead of parcels the spray is represented by the volume fractions of one bulk liquid and several droplet size class phases. A further drawback of the DDM approach is that increasing the spatial resolution of the computational grid leads to a reduced statistical convergence, since the number of spray parcels per computational cell becomes smaller. It is desirable to combine the benefits of both spray approaches in coupled CFD simulations. Therefore, the dense spray region is simulated separately with the Eulerian spray approach on a highly resolved mesh covering only the region close to the nozzle orifice. The entire engine domain with combustion and emission models is simulated with the Eulerian-Lagrangian spray approach for the dilute spray region. The two simulations are coupled through exchange of boundary conditions and model source terms. An on-line coupling interface manages the data transfer between the two simulation clients, i.e., Eulerian spray and engine client. The aim of this work is to extend the coupled spray approach in terms of exchanging combustion related heat and species sources, and consequently creating the link between Eulerian spray and combustion models. The results show mixture formation and combustion in real-case engine simulations, and demonstrate the feasibility of spray model combination in engineering applications.

Lagrangian time scales and its relationship to Eulerian equivalents in turbulent channel flow

Lagrangian and Eulerian time scales were obtained from the direct numerical simulation of turbulent channel flow at two Reynolds numbers based on the friction velocity and channel half-height, Rer= 80, 100. The Lagrangian integral time scales and time microscales were compared to their Eulerian equivalents. It is found that the ratio of Lagrangian to TL Eulerian integral time scales is given by TE/TiE= 1 ＋ 0.1y＋ for y＋ ≤ 10, and that the ratios between the Lagrangian to theEulerian time microscales are almost the same irrespective of the components. Those increase with y＋ are approximated by ≈ 2.75 - 1.75 exp （-v＋/a） . These results also show that these expressions are independent of the Reynolds number.

Numerical simulation of magnetic drug targeting with Eulerian-Lagrangian model and effect of viscosity modification due to diabetics

The targeted drug delivery and targeted drug therapy transport a drug directly to the center of the disease under various conditions and thereby treat it deliber- ately without effects on the body. This paper studies the magnetic drug targeting （MDT） technique by particle tracking in the presence of magnetic field in bifurcation vessels of a healthy person and a diabetes patient. The Lagrangian particle tracking is performed to estimate particle behavior under effects of imposed magnetic field gradients along the bifurcation. The results show that the magnetic field increases the volume fraction of particle in the target region, and the efficiency of MDT on a patient with the diabetes disease is better than a healthy person. Also, for the higher magnetic numbers, the flow in the upper branch is strongly affected by the magnetic field.

Lagrangian coherent structure analysis on transport of Acetes chinensis along coast of Lianyungang,China

Spatial heterogeneity or“patchiness”of plankton distributions in the ocean has always been an attractive and challenging scientific issue to oceanographers.We focused on the accumulation and dynamic mechanism of the Acetes chinensis in the Lianyungang nearshore licensed fishing area.The Lagrangian frame approaches including the Lagrangian coherent structures theory,Lagrangian residual current,and Lagrangian particle-tracking model were applied to find the transport pathways and aggregation characteristics of Acetes chinensis.There exist some material transport pathways for Acetes chinensis passing through the licensed fishing area,and Acetes chinensis is easy to accumulate in the licensed fishing area.The main mechanism forming this distribution pattern is the local circulation induced by the nonlinear interaction of topography and tidal flow.Both the Lagrangian coherent structure analysis and the particle trajectory tracking indicate that Acetes chinensis in the licensed fishing area come from the nearshore estuary.This work contributed to the adjustment of licensed fishing area and the efficient utilization of fishery resources.

对抗条件下基于SAC-Lagrangian的UAV智能规划

无人机因其低成本、可消耗、分布部署、敏捷灵活的优势,在多个民用领域大放异彩。但受其智能化程度限制,如何在复杂对抗条件下自主安全完成任务仍存在巨大挑战。针对目前无人机任务规划存在的智能性和安全性问题,提出一种基于安全强化学习算法的无人机智能规划方法(SAC-Lagrangian)。考虑了雷达威胁、禁飞区安全约束和地导对抗条件,将任务规划问题建模为约束马尔可夫决策过程(CMDP),通过拉格朗日乘子法变为对偶问题,采用最大熵柔性行动者-评论家(SAC)算法近似求解最优策略,保证了智能体在遵守安全约束条件下最大化期望回报。仿真结果表明,与其他基线算法相比,所提方法能在保证任务性能的同时确保安全性,适应动态变化的场景,任务完成率达到96%,因此,具有高效、鲁棒和安全的优势。