Detached-eddy Simulation for Time-dependent Turbulent Cavitating Flows

The Reynolds-averaged Navier-Stokes（RANS）,such as the original k-ω two-equation closures,have been very popular in providing good prediction for a wide variety of flows with presently available computational resource.But for cavitating flows,the above equations noticeably over-predict turbulent production and hence effective viscosity.In this paper,the detached eddy simulation（DES） method for time-dependent turbulent cavitating flows is investigated.To assess the state-of-the-art of computational capabilities,different turbulence models including the widely used RANS model and DES model are conducted.Firstly,in order to investigate the grid dependency in computations,different grid sizes are adopted in the computation.Furthermore,the credibility of DES model is supported by the unsteady cavitating flows over a 2D hydrofoil.The results show that the DES model can effectively reduce the eddy viscosities.From the experimental validations regarding the force analysis,frequency and the unsteady cavity visualizations,more favorable agreement with experimental visualizations and measurements are obtained by DES model.DES model is better able to capture unsteady phenomena including cavity length and the resulting hydrodynamic characteristics,reproduces the time-averaged velocity quantitatively around the hydrofoil,and yields more acceptable and unsteady dynamics features.The DES model has shown to be effective in improving the overall predictive capability of unsteady cavitating flows.

Detached-eddy simulation of flow around high-speed train on a bridge under cross winds

In order to describe an investigation of the flow around high-speed train on a bridge under cross winds using detached-eddy simulation(DES), a 1/8th scale model of a three-car high-speed train and a typical bridge model are employed, Numerical wind tunnel technology based on computational fluid dynamics(CFD) is used, and the CFD models are set as stationary models. The Reynolds number of the flow, based on the inflow velocity and the height of the vehicle, is 1.9×10~6. The computations are conducted under three cases, train on the windward track on the bridge(WWC), train on the leeward track on the bridge(LWC) and train on the flat ground(FGC). Commercial software FLUENT is used and the mesh sensitivity research is carried out by three different grids: coarse, medium and fine. Results show that compared with FGC case, the side force coefficients of the head cars for the WWC and LWC cases increases by 14% and 29%, respectively; the coefficients of middle cars for the WWC and LWC increase by 32% and 10%, respectively; and that of the tail car increases by 45% for the WWC whereas decreases by 2% for the LWC case. The most notable thing is that the side force and the rolling moment of the head car are greater for the LWC, while the side force and the rolling moment of the middle car and the tail car are greater for the WWC. Comparing the velocity profiles at different locations, the flow is significantly influenced by the bridge-train system when the air is close to it. For the three cases(WWC, LWC and FGC), the pressure on the windward side of train is mostly positive while that of the leeward side is negative. The discrepancy of train's aerodynamic force is due to the different surface area of positive pressure and negative pressure zone. Many vortices are born on the leeward edge of the roofs. Theses vortices develop downstream, detach and dissipate into the wake region. The eddies develop irregularly, leading to a noticeably turbulent flow at leeward side of train.

Detached-eddy simulation of supersonic turbulent flow over rearward facing step

Detailed formulation of the detached-eddy simulation （DES） based on shear stress transport model （SST） was deduced. Based on finite volume method, detached-eddy simulation of super- sonic turbulent flow over rearward facing step was conducted. Multi-block structured non-conform grid was designed. The obtained distribution of wall pressure coefficient along the downstream of the step was in good agreement with experiment. The structure of supersonic flow over rearward facing step is complicated and contains boundary layer-shock interaction, shock-shock interaction and sec- ondary flow separation. The predictions of the corner pressure and location of the reattachment shock are more exact than the Reynolds-averaged method. Results show that detached-eddy simulation can simulate the supersonic separated flows exactly and has advantages over Reynolds-averaged method. The partition of the detached-eddy simulation can be easily achieved by using multi-block non-conform grid, which reduces the requirement of the computational resources as well.

基于Plant Simulation仿真技术的装配生产线优化研究

【目的】优化装备生产线,缩短产品交付周期。【方法】基于Plant Simulation仿真技术,对装配生产线进行建模、编程、仿真、分析和优化,有效计算产品产量和成本,识别并优化装配生产线的瓶颈工位。【结果】优化后的三维产线仿真模型产量增幅接近10%,生产效率明显提升。【结论】通过将智能制造技能竞赛和科研教学活动相结合,能够凝练总结竞赛内容,促使教师紧盯前沿知识,创新改革教学内容,实现以赛促教、以赛促学、以赛促改、以赛促建的多重目标。

Analysis of Micromechanical Properties at the Interface of Pre-wet SBS Modified Asphalt Mixture Based on Molecular Simulation Technology

The pre-wetting of aggregate surface is a means to improve the interface performance of SBS modified asphalt and aggregate.The effect of pre-wetting technology on the interaction between SBS modified asphalt and aggregate was analyzed by molecular dynamics simulation.The diffusion coefficient and concentration distribution of SBS modified asphalt on aggregate surface are included.The simulation results show that the diffusion coefficient of the aggregate surface of SBS modified asphalt is increased by 47.6%and 70.5%respectively after 110#asphalt and 130#asphalt are pre-wetted.The concentration distribution of SBS modified asphalt on the aggregate surface after pre-wetting is more uniform.According to the results of interface energy calculation,the interface energy of SBS modified bitumen and aggregate can be increased by about 5%after pre-wetting.According to the results of molecular dynamics simulation,the pre-wetting technology can effectively improve the interface workability of SBS modified bitumen and aggregate,so as to improve the interface performance.

Molecular simulation study of the microstructures and properties of pyridinium ionic liquid[HPy][BF_(4)]mixed with acetonitrile

The microstructures and thermodynamic properties of mixed systems comprising pyridinium ionic liquid[HPy][BF_(4)]and acetonitrile at different mole fractions were studied using molecular dynamics simulation in this work.The following properties were determined:density,self-diffusion coefficient,excess molar volume,and radial distribution function.The results show that with an increase in the mole fraction of[HPy][BF_(4)],the self-diffusion coefficient decreases.Additionally,the excess molar volume initially decreases,reaches a minimum,and then increases.The rules of radial distribution functions(RDFs)of characteristic atoms are different.With increasing the mole fraction of[HPy][BF_(4)],the first peak of the RDFs of HA1-F decreases,while that of CT6-CT6 rises at first and then decreases.This indicates that the solvent molecules affect the polar and non-polar regions of[HPy][BF_(4)]differently.

An Algorithm for Cloud-based Web Service Combination Optimization Through Plant Growth Simulation

In order to improve the efficiency of cloud-based web services,an improved plant growth simulation algorithm scheduling model.This model first used mathematical methods to describe the relationships between cloud-based web services and the constraints of system resources.Then,a light-induced plant growth simulation algorithm was established.The performance of the algorithm was compared through several plant types,and the best plant model was selected as the setting for the system.Experimental results show that when the number of test cloud-based web services reaches 2048,the model being 2.14 times faster than PSO,2.8 times faster than the ant colony algorithm,2.9 times faster than the bee colony algorithm,and a remarkable 8.38 times faster than the genetic algorithm.

SolarDesign:An online photovoltaic device simulation and design platform

Solar Design(https://solardesign.cn/)is an online photovoltaic device simulation and design platform that provides engineering modeling analysis for crystalline silicon solar cells,as well as emerging high-efficiency solar cells such as organic,perovskite,and tandem cells.The platform offers user-updatable libraries of basic photovoltaic materials and devices,device-level multi-physics simulations involving optical–electrical–thermal interactions,and circuit-level compact model simulations based on detailed balance theory.Employing internationally advanced numerical methods,the platform accurately,rapidly,and efficiently solves optical absorption,electrical transport,and compact circuit models.It achieves multi-level photovoltaic simulation technology from“materials to devices to circuits”with fully independent intellectual property rights.Compared to commercial softwares,the platform achieves high accuracy and improves speed by more than an order of magnitude.Additionally,it can simulate unique electrical transport processes in emerging solar cells,such as quantum tunneling,exciton dissociation,and ion migration.

Numerical Simulation of Flow and Temperature Distribution in a Bottom-Blown Copper Bath

Smelting with oxygen bottom blowing is one of the main methods used in the frame of copper pyrometallurgy.With this approach,feed materials and oxygen-enriched air are introduced in reversed order to enhance multiphaseflow within the furnace.Understanding the flow structure and temperature distribution in this setup is crucial foroptimizing production.In this study,gas-liquid interactions,and temperature profiles under varying air-injectionconditions are examined by means of numerical simulation for a 3.2 m×20 m furnace.The results indicate that thehigh-velocity regions are essentially distributed near the lance within the reaction region and the flue gas outlet,while low-velocity regions are located close to the furnace walls on both side of the reaction region.Dead regionsappear in the sedimentation region,with gas velocities surpassing those of the molten phase.As the injection rateincreases from 0.50 to 0.80 Nm3/s,the stabilization time of the average liquid surface velocity decreases from 2.6 sto 1.9 s,exhibiting a similar trend to the gas holdup.During stabilization,the average liquid surface velocity risesfrom 0.505 to 0.702 m/s.The average turbulent kinetic energy(TKE)of the fluid in the molten bath increases from0.095 to 0.162 m^(2)/s^(2).The proportion of the area distribution with TKE greater than 0.10 m^(2)/s^(2) and the gas holdupat steady state both rise with an increase in the injection quantity.The maximum splashing height of the melt growsfrom approximately 0.756 to 1.154 m,with the affected area expanding from 14.239 to 20.498 m^(2).Under differentworking conditions with varying injection quantities,the average temperature changes in melt zone and flue gaszone of the furnace are small.The temperature in the melt and in the flue-gas zone spans the interval 1200℃–1257℃,and 1073℃–1121℃,respectively.The temperature distribution of the melt and flue gas reveals a patterncharacterized by elevated temperatures in the reaction zone,gradually transitioning to lower temperatures in thesedimentation region.

Numerical Simulation of Blood Flow Dynamics in a Stenosed Artery Enhanced by Copper and Alumina Nanoparticles

Nanotechnology holds immense importance in the biomedical field due to its ability to revolutionize healthcare on a molecular scale.Motivated by the imperative of enhancing patient outcomes,a comprehensive numerical simulation study on the dynamics of blood flow in a stenosed artery,focusing on the effects of copper and alumina nanoparticles,is conducted.The study employs a 2-dimensional Newtonian blood flow model infused with copper and alumina nanoparticles,considering the influence of a magnetic field,thermal radiation,and various flow parameters.The governing differential equations are first non-dimensionalized to facilitate analysis and subsequently solved using the 4th order collocation method,bvp4c module in MATLAB.This approach obtains velocity and temperature profiles,revealing the impact of relevant parameters crucial in the biomedical field.The findings of this study underscore the significance of understanding blood flow dynamics in stenosed arteries and the potential benefits of utilizing copper and alumina nanoparticles in treatment strategies.The incorporation of nanoparticles introduces novel avenues for enhancing therapeutic interventions,particularly in mitigating the effects of stenosis.The elucidation of velocity and temperature profiles provides valuable insights into the behavior of blood flow under different conditions,thereby informing the development of targeted biomedical applications.The arterial curvature flow parameter influences temperature profiles,with increased parameters promoting more efficient heat dissipation.The elevated values of Prandtl number and thermal radiation parameter showcase the diminished temperature profiles,indicating stronger dominance of momentum diffusion over thermal diffusion and radiative heat transfer mechanism.Sensitivity analysis of the pertinent physical parameters reveals that the Prandtl number has the most significant impact on blood flow dynamics.A statistical analysis of the present results and existing literature has also been included in the study.Overall,this research contributes to advancing our understanding of vascular health and lays the groundwork for innovative approaches in stenosis treatment and related biomedical fields.

Plastic deformation mechanism of γ-phase U–Mo alloy studied by molecular dynamics simulations

Uranium–molybdenum(U–Mo) alloys are critical for nuclear power generation and propulsion because of their superior thermal conductivity, irradiation stability, and anti-swelling properties. This study explores the plastic deformation mechanisms of γ-phase U–Mo alloys using molecular dynamics(MD) simulations. In the slip model, the generalized stacking fault energy(GSFE) and the modified Peierls–Nabarro(P–N) model are used to determine the competitive relationships among different slip systems. In the twinning model, the generalized plane fault energy(GPFE) is assessed to evaluate the competition between slip and twinning. The findings reveal that among the three slip systems, the {110}<111>slip system is preferentially activated, while in the {112}<111> system, twinning is favored over slip, as confirmed by MD tensile simulations conducted in various directions. Additionally, the impact of Mo content on deformation behavior is emphasized. Insights are provided for optimizing process conditions to avoid γ → α′′ transitions, thereby maintaining a higher proportion of γ-phase U–Mo alloys for practical applications.

Land use/cover change and ecological network in Gansu Province,China during 2000-2020 and their simulations in 2050

Land use/cover change(LUCC)constitutes the spatial and temporal patterns of ecological security,and the construction of ecological networks is an effective way to ensure ecological security.Exploring the spatial and temporal change characteristics of ecological network and analyzing the integrated relationship between LUCC and ecological security are crucial for ensuring regional ecological security.Gansu is one of the provinces with fragile ecological environment in China,and rapid changes in land use patterns in recent decades have threatened ecological security.Therefore,taking Gansu Province as the study area,this study simulated its land use pattern in 2050 using patch-generating land use simulation(PLUS)model based on the LUCC trend from 2000 to 2020 and integrated the LUCC into morphological spatial pattern analysis(MSPA)to identify ecological sources and extract the ecological corridors to construct ecological network using circuit theory.The results revealed that,according to the prediction results in 2050,the areas of cultivated land,forest land,grassland,water body,construction land,and unused land would be 63,447.52,39,510.80,148,115.18,4605.21,8368.89,and 161,752.40 km^(2),respectively.The number of ecological sources in Gansu Province would increase to 80,with a total area of 99,927.18 km^(2).The number of ecological corridors would increase to 191,with an estimated total length of 6120.66 km.Both ecological sources and ecological corridors showed a sparse distribution in the northwest and dense distribution in the southeast of the province at the spatial scale.The number of ecological pinch points would reach 312 and the total area would expect to increase to 842.84 km^(2),with the most pronounced increase in the Longdong region.Compared with 2020,the number and area of ecological barriers in 2050 would decrease significantly by 63 and 370.71 km^(2),respectively.In general,based on the prediction results,the connectivity of ecological network of Gansu Province would increase in 2050.To achieve the predicted ecological network in 2050,emphasis should be placed on the protection of cultivated land and ecological land,the establishment of ecological sources in desert areas,the reinforcement of the protection for existing ecological sources,and the construction of ecological corridors to enhance the stability of ecological network.This study provides valuable theoretical support and references for the future construction of ecological networks and regional land resource management decision-making.

LARGE-EDDY AND DETACHED-EDDY SIMULATIONS OF THE SEPARATED FLOW AROUND A CIRCULAR CYLINDER

The separated turbulent flow around a circular cylinder is investigated using Large-Eddy Simulation （LES）, Detached-Eddy Simulation （DES, or hybrid RANS/LES methods）, and Unsteady Reynolds-Averaged Navier-Stokes （URANS）. The purpose of this study is to examine some typical simulation approaches for the prediction of complex separated turbulent flow and to clarify the capability of applying these approaches to a typical case of the separated turbulent flow around a circular cylinder. Several turbulence models, i.e. dynamic Sub-grid Scale （SGS） model in LES, the DES-based Spalart-Allmaras （S-A） and κ-ω Shear-Stress- Transport （SST） models in DES, and the S-A and SST models in URANS, are used in the calculations. Some typical results, e.g., the mean pressure and drag coefficients, velocity profiles, Strouhal number, and Reynolds stresses, are obtained and compared with previous computational and experimental data. Based on our extensive calculations, we assess the capability and performance of these simulation approaches coupled with the relevant turbulence models to predict the separated turbulent flow.

Towards an entropy-based detached-eddy simulation

A concept of entropy increment ratio （ s- ） is introduced for compressible turbulence simulation through a series of direct nu- merical simulations （DNS）. s- represents the dissipation rate per unit mechanical energy with the benefit of independence of freestream Mach numbers. Based on this feature, we construct the shielding function f, to describe the boundary layer region and propose an entropy-based detached-eddy simulation method （SDES）. This approach follows the spirit of delayed de- tached-eddy simulation （DDES） proposed by Spalart et al. in 2005, but it exhibits much better behavior after their performanc- es are compared in the following flows, namely, pure attached flow with thick boundary layer （a supersonic fiat-plate flow with high Reynolds number）, fully separated flow （the supersonic base flow）, and separated-reattached flow （the supersonic cavity-ramp flow）. The Reynolds-averaged Navier-Stokes （RANS） resolved region is reliably preserved and the modeled stress depletion （MSD） phenomenon which is inherent in DES and DDES is partly alleviated. Moreover, this new hybrid strategy is simple and general, making it applicable to other models related to the boundary layer predictions.

Detached-eddy simulation of supersonic flow past a spike-tipped blunt nose

Space vehicle in atmosphere travels mostly at supersonic speed and generates a very strong bow shockwave around its blunt nose. Oblique shock and conical separated flow zone generated by a forward disk-tip spike significantly reduce the drag by reducing the high pressure area on the blunt nose. This study employs improved delayed detached eddy simulation to investigate the characteristic flow structures around a spike-tipped blunt nose at Mach number of 3 and Reynolds number（based on the blunt-body diameter） of 2.72x10;. The calculated time-averaged quantities agree well with experimental data. Characteristic frequencies in different flow regions are extracted using fast Fourier transform. It is found that two distinct instability modes exist: oscillation mode and pulsation mode. The former is related to the foreshock/turbulence interaction with nondimensional frequency at around 0.004. The latter corresponds to the interaction between turbulence and shock structures around the blunt nose, with a typical coherent structure shedding frequency at 0.092.

基于Plant Simulation的化纤自动落丝系统仿真实验分析

基于Plant Simulation软件,构建一种化纤自动落丝系统的3D数字化模型,对系统效率关键输入因子进行参数化设置,将系统各工位三维数字模型导入软件中建立层次化运动机构图形结构,并运用Simltalk语言实现三维动作仿真,通过基础物理参数设置,保证了仿真数字化模型与现实系统更具一致性。通过多级实验设计分析了系统影响因子对于系统效率的影响特性曲线,并进一步通过动态参数化实验方法,计算出双输入因子对系统效能影响的敏感度。该实验结果可为化纤自动落丝系统的建设成本控制与可行性分析提供指导,具有较好的工程应用价值。

基于Simulation X的多轴转向系统功能安全仿真分析

基于ISO 26262《道路车辆功能安全》的标准要求,对特种车的多轴电液转向系统进行分析,以提高系统的安全性和可靠性。运用Simulation X软件建立详细的多轴特种车仿真模型,通过模拟故障模式注入进行仿真试验。对仿真结果和数据进行分析,评估得到故障的严重性、暴露度和可控性,从而确定相应的汽车安全完整性等级。基于故障注入仿真的汽车功能安全分析方法,在系统早期设计阶段可以作为评估架构安全性的重要手段。

基于Plant Simulation的压气机叶片型线加工产线分析与优化

运用专业仿真软件Plant Simulation,根据压气机叶片型线机械加工工艺特点和物料运行流程建立生产线仿真模型,从产能、设备利用率及产线瓶颈等多方面进行仿真分析与优化。结果表明:Plant Simulation仿真平台的仿真可以找出规划设计方案中存在的问题并验证方案的合理性。该仿真结果为型线机械加工产线的优化设计提供了可靠依据,达到了节约投资成本和缩短设计周期的目的。

Exhaustive review of acceleration strategies for Monte Carlo simulations in photon transit

Monte Carlo simulation techniques have become the quintessence and a pivotal nexus of inquiry in the realm of simulating photon movement within biological fabrics.Through the stochastic sampling of tissue archetypes delineated by explicit optical characteristics,Monte Carlo simulations possess the theoretical capacity to render unparalleled accuracy in the depiction of exceedingly intricate phenomena.Nonetheless,the quintessential challenge associated with Monte Carlo simulation methodologies resides in their extended computational duration,which significantly impedes the refinement of their precision.Consequently,this discourse is specifically dedicated to exploring innovations in strategies and technologies aimed at expediting Monte Carlo simulations.It delves into the foundational concepts of various acceleration tactics,evaluates these strategies concerning their speed,accuracy,and practicality,and amalgamates a comprehensive overview and critique of acceleration methodologies for Monte Carlo simulations.Ultimately,the discourse envisages prospective trajectories for the employment of Monte Carlo techniques within the domain of tissue optics.

Global hybrid simulations of soft X-ray emissions in the Earth’s magnetosheath

Earth’s magnetopause is a thin boundary separating the shocked solar wind plasma from the magnetospheric plasmas,and it is also the boundary of the solar wind energy transport to the magnetosphere.Soft X-ray imaging allows investigation of the large-scale magnetopause by providing a two-dimensional(2-D)global view from a satellite.By performing 3-D global hybrid-particle-in-cell(hybrid-PIC)simulations,we obtain soft X-ray images of Earth’s magnetopause under different solar wind conditions,such as different plasma densities and directions of the southward interplanetary magnetic field.In all cases,magnetic reconnection occurs at low latitude magnetopause.The soft X-ray images observed by a hypothetical satellite are shown,with all of the following identified:the boundary of the magnetopause,the cusps,and the magnetosheath.Local X-ray emissivity in the magnetosheath is characterized by large amplitude fluctuations(up to 160%);however,the maximum line-of-sight-integrated X-ray intensity matches the tangent directions of the magnetopause well,indicating that these fluctuations have limited impact on identifying the magnetopause boundary in the X-ray images.Moreover,the magnetopause boundary can be identified using multiple viewing geometries.We also find that solar wind conditions have little effect on the magnetopause identification.The Solar wind Magnetosphere Ionosphere Link Explorer(SMILE)mission will provide X-ray images of the magnetopause for the first time,and our global hybrid-PIC simulation results can help better understand the 2-D X-ray images of the magnetopause from a 3-D perspective,with particle kinetic effects considered.