Computational fluid dynamics-discrete element method simulation of stirred tank reactor for graphene production

Liquid phase exfoliation(LPE)process for graphene production is usually carried out in stirred tank reactor and the interactions between the solvent and the graphite particles are important as to improve the production efficiency.In this paper,these interactions were revealed by computational fluid dynamics–discrete element method(CFD-DEM)method.Based on simulation results,both liquid phase flow hydrodynamics and particle motion behavior have been analyzed,which gave the general information of the multiphase flow behavior inside the stirred tank reactor as to graphene production.By calculating the threshold at the beginning of graphite exfoliation process,the shear force from the slip velocity was determined as the active force.These results can support the optimization of the graphene production process.

Review:Recent Development of High⁃Order⁃Spectral Method Combined with Computational Fluid Dynamics Method for Wave⁃Structure Interactions

The present paper reviews the recent developments of a high⁃order⁃spectral method(HOS)and the combination with computational fluid dynamics(CFD)method for wave⁃structure interactions.As the numerical simulations of wave⁃structure interaction require efficiency and accuracy,as well as the ability in calculating in open sea states,the HOS method has its strength in both generating extreme waves in open seas and fast convergence in simulations,while computational fluid dynamics(CFD)method has its advantages in simulating violent wave⁃structure interactions.This paper provides the new thoughts for fast and accurate simulations,as well as the future work on innovations in fine fluid field of numerical simulations.

Computational Fluid Dynamics(CFD) Analysis and Optimization of Reconstructed Intake System of Cylinder Head Based on Slicing Reverse Method

To find out and improve the flow characteristics inside the intake system of cylinder head,the application of computational fluid dynamics(CFD)in the evaluation and optimization of the reconstructed intake system based on slicing reverse method was proposed.The flow characteristics were found out through CFD,and the velocity vector field,pressure field and turbulent kinetic energy field for different valve lifts were discussed,which were in good agreement with experimental data,and the quality of reconstruction was evaluated.In order to improve its flow characteristic,an optimization plan was proposed.The results show that the flow characteristics after optimization are obviously improved.The results can provide a reference for the design and optimization of the intake system of cylinder head.

Analysis of Temperature Rise in High-Speed Permanent Magnet Synchronous Traction Motors by Coupling the Equivalent Thermal Circuit Method and Computational Fluid Dynamics

To solve the problem of temperature rise caused by the high power density of high-speed permanent magnet synchronous traction motors,the temperature rise of various components in the motor is analyzed by coupling the equivalent thermal circuit method and computational fluid dynamics.Also,a cooling strategy is proposed to solve the problem of temperature rise,which is expected to prolong the service life of these devices.First,the theoretical bases of the approaches used to study heat transfer and fluid mechanics are discussed,then the fluid flow for the considered motor is analyzed,and the equivalent thermal circuit method is introduced for the calculation of the temperature rise.Finally,the stator,rotor loss,motor temperature rise,and the proposed cooling method are also explored through experiments.According to the results,the stator temperature at 50,000 r/min and 60,000 r/min at no-load operation is 68℃ and 76℃,respectively.By monitoring the temperature of the air outlets inside and outside the motor at different speeds,it is also found that the motor reaches a stable temperature rise after 65 min of operation.Coupling of the thermal circuit method and computational fluid dynamics is a strategy that can provide the average temperature rise of each component and can also comprehensively calculate the temperature of each local point.We conclude that a hybrid cooling strategy based on axial air cooling of the inner air duct of the motor and water cooling of the stator can meet the design requirements for the ventilation and cooling of this type of motors.

Computational fluid dynamics simulations of respiratory airflow in human nasal cavity and its characteristic dimension study

To study the airflow distribution in human nasal cavity during respiration and the characteristic parameters of nasal structure, three-dimensional, anatomically accurate representations of 30 adult nasal cavity models were recons- tructed based on processed tomography images collected from normal people. The airflow fields in nasal cavities were simulated by fluid dynamics with finite element software ANSYS. The results showed that the difference of human nasal cavity structure led to different airflow distribution in the nasal cavities and variation of the main airstream passing through the common nasal meatus. The nasal resistance in the regions of nasal valve and nasal vestibule accounted for more than half of the overall resistance. The characteristic model of nasal cavity was extracted on the basis of characteristic points and dimensions deduced from the original models. It showed that either the geometric structure or the airflow field of the two kinds of models was similar. The characteristic dimensions were the characteristic parameters of nasal cavity that could properly represent the original model in model studies on nasal cavity.

Effect of a two-phase wedge-sliding model on the ingredient drift of a stable mixed fluid and its computing method

A two-phase wedge-sliding model is developed based on the micro-cellular structure and minimum entropy theory of a stable system, and it is used to describe the ingredient distribution of a mixed fluid in a non-uniform stress field and to analyse its phase drift phenomenon. In the model, the drift-inhibition angle and the expansion-inhibition angle are also deduced and used as evaluating indexes to describe the drifting trend of different ingredients among the mixed fluids. For solving above two indexes of the model, a new calculation method is developed and used to compute the phase distributions of multiphase fluid at peak stress and gradient area stress, respectively. As an example, the flow process of grease in a pipe is analysed by simulation method and used to verify the validity of the model.

A COMPUTATIONAL METHOD FOR THE NAVIER-STOKES EQUATIONS AT ALL SPEEDS

A PLU-SGS method based on a time-derivative preconditioning algorithm and LU-SGS method is developed in order to calculate the Navier-Stokes equations at all speeds. The equations were discretized using A USMPW scheme in conjunction with the third-order MUSCL scheme with Van Leer limiter. The present method was applied to solve the multidimensional compressible Navier-Stokes equations in curvilinear coordinates. Characteristic boundary conditions based on the eigensystem of the preconditioned equations were employed. In order to examine the performance of present method, driven-cavity flow at various Reynolds numbers and viscous flow through a convergent-divergent nozzle at supersonic were selected to rest this method. The computed results were compared with the experimental data or the other numerical results available in literature and good agreements between them are obtained. The results show that the present method is accurate, self-adaptive and stable for a wide range of flow conditions from low speed to supersonic flows.

Deep Learning Applied to Computational Mechanics:A Comprehensive Review,State of the Art,and the Classics

Three recent breakthroughs due to AI in arts and science serve as motivation:An award winning digital image,protein folding,fast matrix multiplication.Many recent developments in artificial neural networks,particularly deep learning(DL),applied and relevant to computational mechanics(solid,fluids,finite-element technology)are reviewed in detail.Both hybrid and pure machine learning(ML)methods are discussed.Hybrid methods combine traditional PDE discretizations with ML methods either(1)to help model complex nonlinear constitutive relations,(2)to nonlinearly reduce the model order for efficient simulation(turbulence),or(3)to accelerate the simulation by predicting certain components in the traditional integration methods.Here,methods(1)and(2)relied on Long-Short-Term Memory(LSTM)architecture,with method(3)relying on convolutional neural networks.Pure ML methods to solve(nonlinear)PDEs are represented by Physics-Informed Neural network(PINN)methods,which could be combined with attention mechanism to address discontinuous solutions.Both LSTM and attention architectures,together with modern and generalized classic optimizers to include stochasticity for DL networks,are extensively reviewed.Kernel machines,including Gaussian processes,are provided to sufficient depth for more advanced works such as shallow networks with infinite width.Not only addressing experts,readers are assumed familiar with computational mechanics,but not with DL,whose concepts and applications are built up from the basics,aiming at bringing first-time learners quickly to the forefront of research.History and limitations of AI are recounted and discussed,with particular attention at pointing out misstatements or misconceptions of the classics,even in well-known references.Positioning and pointing control of a large-deformable beam is given as an example.

CVBEM and FVM Computational Model Comparison for Solving Ideal Fluid Flow in a 90-Degree Bend

While finite volume methodologies (FVM) have predominated in fluid flow computations, many flow problems, including groundwater models, would benefit from the use of boundary methods, such as the Complex Variable Boundary Element Method (CVBEM). However, to date, there has been no reporting of a comparison of computational results between the FVM and the CVBEM in the assessment of flow field characteristics. In this work, the CVBEM is used to develop a flow field vector outcome of ideal fluid flow in a 90-degree bend which is then compared to the computational results from a finite volume model of the same situation. The focus of the modelling comparison in the current work is flow field trajectory vectors of the fluid flow, with respect to vector magnitude and direction. Such a comparison is necessary to validate the development of flow field vectors from the CVBEM and is of interest to many engineering flow problems, specifically groundwater modelling. Comparison of the CVBEM and FVM flow field trajectory vectors for the target problem of ideal flow in a 90-degree bend shows good agreement between the considered methodologies.

Stabilization for Equal-Order Polygonal Finite Element Method for High Fluid Velocity and Pressure Gradient

This paper presents an adapted stabilisation method for the equal-order mixed scheme of finite elements on convex polygonal meshes to analyse the high velocity and pressure gradient of incompressible fluid flows that are governed by Stokes equations system.This technique is constructed by a local pressure projection which is extremely simple,yet effective,to eliminate the poor or even non-convergence as well as the instability of equal-order mixed polygonal technique.In this research,some numerical examples of incompressible Stokes fluid flow that is coded and programmed by MATLAB will be presented to examine the effectiveness of the proposed stabilised method.

MIXED COMPATIBLE ELEMENT AND MIXED HYBRID INCOMPATIBLE ELEMENT VARIATIONAL METHODS IN DYNAMICS OF VISCOUS BAROTROPIC FLUIDS

This paper presents and proves the mixed compatible finite element variationalprinciples in dynamics of viscous barotropic fluids. When the principles are proved, itis found that the compatibility conditions of stress can be naturally satisfied. The gene-rallzed variational principles with mixed hybrid incompatible finite elements are alsopresented and proved, and they can reduce the computation of incompatible elements indynamics of viscous barotropic flows.

Application of the Hierarchical Functions Expansion Method for the Solution of the Two Dimensional Navier-Stokes Equations for Compressible Fluids in High Velocity

This work presents a new application for the Hierarchical Function Expansion Method for the solution of the Navier-Stokes equations for compressible fluids in two dimensions and in high velocity. This method is based on the finite elements method using the Petrov-Galerkin formulation, know as SUPG (Streamline Upwind Petrov-Galerkin), applied with the expansion of the variables into hierarchical functions. To test and validate the numerical method proposed as well as the computational program developed simulations are performed for some cases whose theoretical solutions are known. These cases are the following: continuity test, stability and convergence test, temperature step problem, and several oblique shocks. The objective of the last cases is basically to verify the capture of the shock wave by the method developed. The results obtained in the simulations with the proposed method were good both qualitatively and quantitatively when compared with the theoretical solutions. This allows concluding that the objectives of this work are reached.

My Way: A Computational Autobiography

In this paper,the author recounts his forty-year plus struggle to find a sound basis for understanding the computational fluid dynamics of compressible flow.

A Generic Computational Solution of a Natural Convection Flow past an Infinite Vertical Porous Plate

A simulation was carried out on an unsteady flow of a viscous, incompressible and electrically conducting fluid past an infinite vertical porous plate. A generic computer program using the Galerkin finite element method is employed to solve the coupled non-linear differential equations for velocity and temperature fields. The diffusion equation, the energy equation, the momentum equations and other relevant parameters are transformed into interpretable postfix codes. Numerical calculations are carried out on the flow fields both in the presence of cooling and heating of the plate by free convection currents. The effects of the dimensionless parameters, namely, the Prandtl number, the Eckert number, the modified Grashof number, the Schmidt number and the time on the temperature and velocity distributions are discussed.

恶劣海况下船舶砰击颤振响应特性数值计算与试验研究

[目的]针对恶劣海况下船舶所受砰击颤振响应现象,探究船舶非线性波浪载荷与瞬态高幅值砰击载荷的耦合作用。[方法]采用计算流体动力学与有限元方法(CFD-FEM)相结合的双向流固耦合方法对S175集装箱船进行数值仿真计算,并与试验结果及切片理论计算结果进行对比验证;采用分段变截面弹性龙骨梁模型开展船舶的砰击颤振特性模型试验,基于CFD-FEM双向流固耦合方法开展船艏砰击载荷及高频非线性砰击颤振响应特性分析,并与模型试验结果进行对比验证。[结果]结果显示,波浪砰击载荷对船艏颤振响应的影响不可忽视,6级海况下由砰击颤振诱发的二阶高频成分分量占低频波浪弯矩的59.86%。[结论]采用基于CFD-FEM的双向流固耦合方法可准确计算船首砰击颤振响应;在高海况下船舶所受非线性波浪载荷及结构动态响应易受船首瞬态砰击载荷的影响,在船舶结构设计与安全评估中需考虑高频砰击颤振的情况。

嫦娥五号发动机降落羽流扬尘特性研究

在月球探测器着陆过程中,发动机羽流与月面相互作用后溅起的月尘是月面环境危害的主要来源。本文以嫦娥五号任务测试数据作为仿真入口条件,采用计算流体动力学(CFD)两阶段法建立了嫦娥五号任务中使用的喷管1∶1模型和真空羽流扩散侵蚀模型,研究了喷管在不同降落高度下的侵蚀速率,并计算了发动机距离月面高度为0.5~2.0 m范围时月尘颗粒的运动轨迹、扬尘角和速度特性。结果表明,基于剪切应力得到的最大侵蚀速率为8.83 kg/m^(2)s,随着高度增加,侵蚀速率降低,与嫦娥五号降落相机相同高度下的分析结果一致。粒径为1、70ìm的月尘颗粒最大扬尘高度分别为0.72、0.36 m,最大速度分别为2520、1010 m/s。不同粒径月尘的扬尘角范围为1.44°~2.27°,计算的扬尘角与Apollo探月任务中的结果相近。

孔的排布方式对孔型阻尼密封泄漏特性和动力学性能影响研究

为揭示孔的排布方式对孔型阻尼密封性能的影响规律以及探究最佳排布方式,基于多频椭圆轨迹涡动模型非定常计算流体力学摄动数值方法,研究了不同孔径、孔排布方式下孔型阻尼密封的泄漏特性和动力学性能。采用孔型阻尼密封试件试验数据,验证了数值方法的可靠性和准确性,通过计算分析了不同孔径(3.175、12.268 mm)、孔排布方式(轴向顺排、周向顺排)、进口预旋(0、0.5)对孔型阻尼密封的泄漏量、动力特性系数和腔室动态压力的影响规律。数值结果表明:大孔径孔型阻尼密封动力学性能较优,具有较高的有效阻尼和较低的穿越频率,但封严性能较差,对应的泄漏量增加了6%~10%;小孔径孔型阻尼密封的泄漏特性和动力学性能对孔的排布方式不敏感,泄漏量偏差小于0.5%,动力特性系数偏差小于7%;对于大孔径孔型阻尼密封,采用周向顺排的方式可获得更好的封严性能以及更优的动力学性能,泄漏量可降低4%,有效阻尼的穿越频率也得到显著减小。该研究可为孔型阻尼密封的设计提供参考。

基于CFD-DEM的流态化速冻影响因素研究

利用Fluent多孔介质理论建立了速冻流化床等效流域模型;以表观气速、入口温度、脉动频率为试验因素,根据CFD-DEM数值模拟方法对典型颗粒形状物料的流态化和冻结进行单因素试验分析,探究了各试验因素对物料冻结的影响规律;基于单因素试验结果,以平均降温速率为评价指标,进行三因素、四水平的正交试验,试验结果表明,在进风量一定的情况下,影响程度由强到弱依次为入口温度、脉动频率、表观气速。以期分析结果对脉冲气流式速冻工程实践具有一定的参考作用。

基于计算流体力学-离散单元法的U形管冲蚀磨损数值模拟研究

采用计算流体力学-离散单元法(CFD-DEM)和切向撞击能模型(SIEM)来研究气固两相流中高速煤粉颗粒冲击导致的U形弯管冲蚀磨损问题。分析了颗粒粒径、气体流速、固体颗粒形状对其冲蚀磨损特性的影响,并验证粗粒化模型应用于此磨损过程的准确性。结果表明,U形弯头的冲蚀磨损位置主要有3处;磨损率随粒径的增大呈现整体减小的趋势;U形管拱背磨损率随气速的增加而增大,且前2处磨损率峰值位置随气速的增加向入口方向偏移;长椭球对U形管造成的磨损最严重,扁椭球次之,球形最轻;粗粒化CFD-DEM能够基本准确预测U形管冲蚀磨损,且显著降低计算时间。

考虑跨洋特征及碎冰对快速性影响的极地探险邮轮型线优化

[目的]全球减碳背景下,为应对极地船舶设计建造的环保要求,需开展跨洋航行时碎冰水域对极地探险邮轮船体型线优化设计的影响规律研究,以获得最佳节能型线。[方法]针对极地探险邮轮的跨洋航行特征,采用航区权重方法进行量化评估,分析碎冰对阻力和推进效率的影响。应用计算流体力学耦合离散元法(CFD-DEM)来分析螺旋桨的碎冰水域性能,建立以联合自航功率为目标的优化模型,进而对全船参数化模型开展设计航速下的优化计算。[结果]计算结果表明,优化后的船型可以满足排水量要求,有效降低了2种水域下的航行功率,其联合自航功率降低了9.71%。[结论]研究成果验证了基于权重优化方法的可行性和合理性,可为极地探险邮轮的型线和推进器优化设计提供参考。