On the Semianalytical Two-Body Regularization in N-Body Simulations

A two-body regularization for N-body problem based on perturbation theory for Keplerian problem is discussed. We provide analytical estimations of accuracy and conduct N-body experiments in order to compare it with state-of-the-art Hermite integrator. It is shown that this regularization keeps some features that allow overcoming KS-regularization in some particular cases.

On a New Method of N-Body Simulations

The theoretical foundation of a new N-body simulation method for the dynamics of large numbers (N > 106) of gravitating bodies is described. The new approach is founded on the probability description of the physical parameters and a similarity method which permits a manifold reduction of the calculation time for the evolution of “large” systems. This is done by averaging the results of calculations over an ensemble of many “small” systems with total particle number in the ensemble equal to the number of stars in the large system. The method is valid for the approximate calculation of the evolution of large systems, including dissipative systems like AGN containing a supermassive black hole, accretion disc, and the surrounding stellar cluster.

Numerical Investigations of a New N-body Simulation Method

Numerical investigation of a new similarity method (the Aldar-Kose method) for N-body simulations is described. Using this method we have carried out numerical simulations for two tasks: 1) calculation of the temporal behavior of different physical parameters of active galactic nuclei (AGN) containing a super massive black hole (SMBH), an accretion disk, and a compact stellar cluster;2) calculation of the stellar capture rate to the central SMBH without accretion disk. The calculations show good perspectives for applications of the similarity method to optimize the evolution model calculations of large stellar systems and of AGN.

N-Body Simulation Inspired by Metaheuristics Optimization

The N-body problem in classical physics, is the calculation of force ofgravitational attraction of heavenly bodies towards each other. Solving this problem for many heavenly bodies has always posed a challenge to physicists andmathematicians. Large number of bodies, huge masses, long distances and exponentially increasing number of equations of motion of the bodies have been themajor hurdles in solving this problem for large and complex galaxies. Adventof high performance computational machines have mitigated the problem to muchextent, but still for large number of bodies it consumes huge amount of resourcesand days for computation. Conventional algorithms have been able to reduce thecomputational complexity from O n2 ð Þ to O nlogn ð Þ by splitting the space into atree or mesh network, researchers are still looking for improvements. In thisresearch work we propose a novel solution to N-body problem inspired by metaheuristics algorithms. The proposed algorithm is simulated for various time periods of selected heavenly bodies and analyzed for speed and accuracy. Theresults are compared with that of conventional algorithms. The outcomes showabout 50% time saving with almost no loss in accuracy. The proposed approachbeing a metaheuristics optimization technique, attempts to find optimal solution tothe problem, searching the entire space in a unique and efficient manner in a verylimited amount of time.

N-Body Simulations of Gas-Free Disc Galaxies with SMBH Seed in Binary Systems

We have shown the outcome of N-body simulations of the interactions of two disc galaxies without gas with the same mass. Both disc galaxies have halos of dark matter, central bulges and initial supermassive black hole (SMBH) seeds at their centers. The purpose of this work is to study the mass and dynamical evolution of the initial SMBH seed during a Hubble cosmological time. It is a complementation of our previous paper with different initial orbit conditions and by introducing the SMBH seed in the initial galaxy. The disc of the secondary galaxy has a coplanar or polar orientation in relation to the disc of the primary galaxy and their initial orbit are eccentric and prograde. The primary and secondary galaxies have mass and size of Milky Way with an initial SMBH seed. We have found that the merger of the primary and secondary discs can result in a final normal disc or a final warped disc. After the fusion of discs, the final one is thicker and larger than the initial disc. The tidal effects are very important, modifying the evolution of the SMBH in the primary and secondary galaxy differently. The mass of the SMBH of the primary galaxy has increased by a factor ranging from 52 to 64 times the initial seed mass, depending on the experiment. However, the mass of the SMBH of the secondary galaxy has increased by a factor ranging from 6 to 33 times the initial SMBH seed mass, depending also on the experiment. Most of the accreted particles have come from the bulge and from the halo, depleting their particles. This could explain why the observations show that the SMBH with masses of approximatelyis found in many bulgeless galaxies. Only a small number of the accreted particles has come from the disc. In some cases of final merging stage of the two galaxies, the final SMBH of the secondary galaxy was ejected out of the galaxy.

A Simulation of the Dependence of Tidal Interaction on Galaxy Type in Compact Groups

We have investigated the role that different galaxy types have in galaxy-galaxy interactions in compact groups. N-body simulations of 6 galaxies consisting of a differing mixture of galaxy types were run to compare the relative importance of galaxy population demographic on evolution. Three different groups with differing galaxy content were tested: all spiral, a single elliptical and 50% elliptical. Tidal interaction strength and duration were recorded to assess the importance of an interaction. A group with an equal number of spiral and elliptical galaxies has some of the longest and strongest interactions with elliptical-elliptical interactions being most significant. These elliptical-elliptical interactions are not dominated by a single large event but consist of multiple interactions. Elliptical galaxies tidally interacting with spiral galaxies, have the next strongest interaction events. For the case when a group only has a single elliptical, the largest magnitude tidal interaction is an elliptical on a spiral. Spirals interact with each other through many small interactions. For a spiral only group, the interactions are the weakest compared to the other group types. These spiral interactions are not dominated by any singular event that might be expected to lead to a merger but are more of an ongoing harassment. These results suggest that within a compact group, early type galaxies will not form via merger out of an assemblage of spiral galaxies but rather that compact groups, in effect form around an early type galaxy.

基于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.

The Probability Density Function Related to Shallow Cumulus Entrainment Rate and Its Influencing Factors in a Large-Eddy Simulation

The process of entrainment-mixing between cumulus clouds and the ambient air is important for the development of cumulus clouds.Accurately obtaining the entrainment rate(λ)is particularly important for its parameterization within the overall cumulus parameterization scheme.In this study,an improved bulk-plume method is proposed by solving the equations of two conserved variables simultaneously to calculateλof cumulus clouds in a large-eddy simulation.The results demonstrate that the improved bulk-plume method is more reliable than the traditional bulk-plume method,becauseλ,as calculated from the improved method,falls within the range ofλvalues obtained from the traditional method using different conserved variables.The probability density functions ofλfor all data,different times,and different heights can be well-fitted by a log-normal distribution,which supports the assumed stochastic entrainment process in previous studies.Further analysis demonstrate that the relationship betweenλand the vertical velocity is better than other thermodynamic/dynamical properties;thus,the vertical velocity is recommended as the primary influencing factor for the parameterization ofλin the future.The results of this study enhance the theoretical understanding ofλand its influencing factors and shed new light on the development ofλparameterization.

基于Simulation的薄壁桶搬运过程的形变研究

桶料的搬运方式主要有托底式、夹取式和挂取(吊取)式,文中针对薄壁桶的搬运方式进行模拟分析,通过模拟分析的具体数值来进行合理的选择。主要对挂取(吊取)和夹取2种薄壁桶的搬运方式进行模拟分析,基于Soildworks中的Simulation进行计算,并得出具体数值,选取合理的薄壁桶搬运方式。模拟分析表明:在同等桶壁厚的情况下,夹取式桶料搬运形式薄壁桶的变形程度较大,主要体现在薄壁桶的夹持部位,如特定场景需夹持薄壁桶的桶身进行搬运时,还需适当加厚薄壁桶的桶壁来抵抗夹持板的夹持力,使薄壁桶不会产生较大的形变量。

Microwave-induced thermoacoustic elastic imaging:A simulation study

Microwave-induced thermoacoustic imaging(MTI)has the advantages of high resolution,high contrast,non-ionization,and non-invasive.Recently,MTI was used in the¯eld of breast cancer screening.In this paper,based on the¯nite element method(FEM)and COMSOL Multiphysics software,a three-dimensional breast cancer model suitable for exploring the MTI process is proposed to investigate the in°uence of Young's modulus(YM)of breast cancer tissue on MTI.It is found that the process of electromagnetic heating and initial pressure generation of the entire breast tissue is earlier in time than the thermal expansion process.Besides,compared with normal breast tissue,tumor tissue has a greater temperature rise,displacement,and pressure rise.In particular,YM of the tumor is related to the speed of thermal expansion.In particular,the larger the YM of the tumor is,the higher the heating and contraction frequency is,and the greater the maximum pressure is.Di®erent Young's moduli correspond to di®erent thermoacoustic signal spectra.In MTI,this study can be used to judge di®erent degrees of breast cancer based on elastic imaging.In addition,this study is helpful in exploring the possibility of microwave-induced thermoacoustic elastic imaging(MTAE).

Simulation study of the performance of the Very Large Area gamma-ray Space Telescope

The Very Large Area gamma-ray Space Telescope(VLAST)is a mission concept proposed to detect gamma-ray photons through both Compton scattering and electron–positron pair production mechanisms,thus enabling the detection of photons with energies ranging from MeV to TeV.This project aims to conduct a comprehensive survey of the gamma-ray sky from a low-Earth orbit using an anti-coincidence detector,a tracker detector that also serves as a low-energy calorimeter,and a high-energy imaging calorimeter.We developed a Monte Carlo simulation application of the detector using the GEANT4 toolkit to evaluate the instrument performance,including the effective area,angular resolution,and energy resolution,and explored specific optimizations of the detector configuration.Our simulation-based analysis indicates that the current design of the VLAST is physically feasible,with an acceptance above 10 m^(2)sr which is four times larger than that of the Fermi-LAT,an energy resolution better than 2%at 10 GeV,and an angular resolution better than 0.2◦at 10 GeV.The VLAST project promises to make significant contributions to the field of gamma-ray astronomy and enhance our understanding of the cosmos.

Magnetopause properties at the dusk magnetospheric flank from global magnetohydrodynamic simulations,the kinetic Vlasov equilibrium,and in situ observations--Potential implications for SMILE

We derived the properties of the terrestrial magnetopause(MP)from two modeling approaches,one global–fluid,the other local–kinetic,and compared the results with data collected in situ by the Magnetospheric Multiscale 2(MMS2)spacecraft.We used global magnetohydrodynamic(MHD)simulations of the Earth’s magnetosphere(publicly available from the NASA-CCMC[National Aeronautics and Space Administration–Community Coordinated Modeling Center])and local Vlasov equilibrium models(based on kinetic models for tangential discontinuities)to extract spatial profiles of the plasma and field variables at the Earth’s MP.The global MHD simulations used initial solar wind conditions extracted from the OMNI database at the time epoch when the MMS2 observes the MP.The kinetic Vlasov model used asymptotic boundary conditions derived from the same in situ MMS measurements upstream or downstream of the MP.The global MHD simulations provide a three-dimensional image of the magnetosphere at the time when the MMS2 crosses the MP.The Vlasov model provides a one-dimensional local view of the MP derived from first principles of kinetic theory.The MMS2 experimental data also serve as a reference for comparing and validating the numerical simulations and modeling.We found that the MP transition layer formed in global MHD simulations was generally localized closer to the Earth(roughly by one Earth radius)from the position of the real MP observed by the MMS.We also found that the global MHD simulations overestimated the thickness of the MP transition by one order of magnitude for three analyzed variables:magnetic field,density,and tangential speed.The MP thickness derived from the local Vlasov equilibrium was consistent with observations for all three of these variables.The overestimation of density in the Vlasov equilibrium was reduced compared with the global MHD solutions.We discuss our results in the context of future SMILE(Solar wind Magnetosphere Ionosphere Link Explorer)campaigns for observing the Earth’s MP.

Numerical simulation on sand sedimentation and erosion characteristics around HDPE sheet sand barrier under different wind angles

For the safety of railroad operations,sand barriers are utilized to mitigate wind-sand disaster effects.These disasters,characterized by multi-directional wind patterns,result in diverse angles among the barriers.In this study,using numerical simulations,we examined the behavior of High Density Polyethylene(HDPE)sheet sand barriers under different wind angles,focusing on flow field distribution,windproof efficiency,and sedimentation erosion dynamics.This study discovered that at a steady wind speed,airflow velocity varies as the angle between the airflow and the HDPE barrier changes.Specifically,a 90°angle results in the widest low-speed airflow area on the barrier’s downwind side.If the airflow is not perpendicular to the barrier,it prompts a lateral airflow movement which decreases as the angle expands.The windproof efficiency correlates directly with this angle but inversely with the wind’s speed.Notably,with a wind angle of 90°,wind speed drops by 81%.The minimum wind speed is found at 5.1H(the sand barrier height)on the barrier’s downwind side.As the angle grows,the barrier’s windproof efficiency improves,extending its protective reach.Sedimentation is most prominent on the barrier’s downwind side,as the wind angle shifts from 30°to 90°,the sand sedimentation area on the barrier’s downwind side enlarges by 14.8H.As the angle grows,sedimentation intensifies,eventually overtakes the forward erosion and enlarges the sedimentation area.

Hybrid-Vlasov simulation of soft X-ray emissions at the Earth’s dayside magnetospheric boundaries

Solar wind charge exchange produces emissions in the soft X-ray energy range which can enable the study of near-Earth space regions such as the magnetopause,the magnetosheath and the polar cusps by remote sensing techniques.The Solar wind Magnetosphere Ionosphere Link Explorer(SMILE)and Lunar Environment heliospheric X-ray Imager(LEXI)missions aim to obtain soft Xray images of near-Earth space thanks to their Soft X-ray Imager(SXI)instruments.While earlier modeling works have already simulated soft X-ray images as might be obtained by SMILE SXI during its mission,the numerical models used so far are all based on the magnetohydrodynamics description of the space plasma.To investigate the possible signatures of ion-kinetic-scale processes in soft Xray images,we use for the first time a global hybrid-Vlasov simulation of the geospace from the Vlasiator model.The simulation is driven by fast and tenuous solar wind conditions and purely southward interplanetary magnetic field.We first produce global X-ray images of the dayside near-Earth space by placing a virtual imaging satellite at two different locations,providing meridional and equatorial views.We then analyze regional features present in the images and show that they correspond to signatures in soft X-ray emissions of mirrormode wave structures in the magnetosheath and flux transfer events(FTEs)at the magnetopause.Our results suggest that,although the time scales associated with the motion of those transient phenomena will likely be significantly smaller than the integration time of the SMILE and LEXI imagers,mirror-mode structures and FTEs can cumulatively produce detectable signatures in the soft X-ray images.For instance,a local increase by 30%in the proton density at the dayside magnetopause resulting from the transit of multiple FTEs leads to a 12%enhancement in the line-of-sight-and time-integrated soft X-ray emissivity originating from this region.Likewise,a proton density increase by 14%in the magnetosheath associated with mirror-mode structures can result in an enhancement in the soft X-ray signal by 4%.These are likely conservative estimates,given that the solar wind conditions used in the Vlasiator run can be expected to generate weaker soft X-ray emissions than the more common denser solar wind.These results will contribute to the preparatory work for the SMILE and LEXI missions by providing the community with quantitative estimates of the effects of small-scale,transient phenomena occurring on the dayside.

Flow field, sedimentation, and erosion characteristics around folded linear HDPE sheet sand fence: Numerical simulation study

Wind and sand hazards are serious in the Milan Gobi area of the Xinjiang section of the Korla Railway. In order to ensure the safe operation of railroads, there is a need for wind and sand protection in heavily sandy areas. The wind and sand flow in the region is notably bi-directional. To shield railroads from sand, a unique sand fence made of folded linear high-density polyethylene(HDPE) is used, aligning with the principle that the dominant wind direction is perpendicular to the fence. This study employed field observations and numerical simulations to investigate the effectiveness of these HDPE sand fences in altering flow field distribution and offering protection. It also explored how these fences affect the deposition and erosion of sand particles. Findings revealed a significant reduction in wind speed near the fence corner;the minimum horizontal wind speed on the leeward side of the first sand fence(LSF) decreased dramatically from 3 m/s to 0.64 m/s. The vortex area on the LSF markedly impacted horizontal wind speeds. Within the LSF, sand deposition was a primary occurrence. As wind speeds increased, the deposition zone shrank, whereas the positive erosion zone expanded. Close to the folded corners of the HDPE sand fence, there was a notable shift from the positive erosion zone to a deposition zone. Field tests and numerical simulations confirmed the high windproof efficiency(WE) and sand resistance efficiency(SE) in the HDPE sand fence. Folded linear HDPE sheet sand fence can effectively slow down the incoming flow and reduce the sand content, thus achieving good wind and sand protection. This study provides essential theoretical guidance for the design and improvement of wind and sand protection systems in railroad engineering.

基于Plant Simulation的生产车间仿真优化研究与应用

针对生产车间焊接阶段物料难以齐套引发的工序等待问题,通过使用Plant Simulation仿真软件构建模型,以车间实际生产数据为输入验证仿真模型的可信性,建立方便交互的可视化交互界面,方便模型复用和生产参数的调整.通过仿真模型结合启发式算法优化生产线,保证了生产车间的物料齐套性,缩短了焊接等待时间,提供了一种使用仿真软件进行生产线优化的可行方案.结果表明,所建立的仿真优化模型能够提供高效的排程方案,有效提高了企业生产效率.