Big Data Application Simulation Platform Design for Onboard Distributed Processing of LEO Mega-Constellation Networks

Due to the restricted satellite payloads in LEO mega-constellation networks(LMCNs),remote sensing image analysis,online learning and other big data services desirably need onboard distributed processing(OBDP).In existing technologies,the efficiency of big data applications(BDAs)in distributed systems hinges on the stable-state and low-latency links between worker nodes.However,LMCNs with high-dynamic nodes and long-distance links can not provide the above conditions,which makes the performance of OBDP hard to be intuitively measured.To bridge this gap,a multidimensional simulation platform is indispensable that can simulate the network environment of LMCNs and put BDAs in it for performance testing.Using STK's APIs and parallel computing framework,we achieve real-time simulation for thousands of satellite nodes,which are mapped as application nodes through software defined network(SDN)and container technologies.We elaborate the architecture and mechanism of the simulation platform,and take the Starlink and Hadoop as realistic examples for simulations.The results indicate that LMCNs have dynamic end-to-end latency which fluctuates periodically with the constellation movement.Compared to ground data center networks(GDCNs),LMCNs deteriorate the computing and storage job throughput,which can be alleviated by the utilization of erasure codes and data flow scheduling of worker nodes.

Simulation of Fracture Process of Lightweight Aggregate Concrete Based on Digital Image Processing Technology

The mechanical properties and failure mechanism of lightweight aggregate concrete(LWAC)is a hot topic in the engineering field,and the relationship between its microstructure and macroscopic mechanical properties is also a frontier research topic in the academic field.In this study,the image processing technology is used to establish a micro-structure model of lightweight aggregate concrete.Through the information extraction and processing of the section image of actual light aggregate concrete specimens,the mesostructural model of light aggregate concrete with real aggregate characteristics is established.The numerical simulation of uniaxial tensile test,uniaxial compression test and three-point bending test of lightweight aggregate concrete are carried out using a new finite element method-the base force element method respectively.Firstly,the image processing technology is used to produce beam specimens,uniaxial compression specimens and uniaxial tensile specimens of light aggregate concrete,which can better simulate the aggregate shape and random distribution of real light aggregate concrete.Secondly,the three-point bending test is numerically simulated.Thirdly,the uniaxial compression specimen generated by image processing technology is numerically simulated.Fourth,the uniaxial tensile specimen generated by image processing technology is numerically simulated.The mechanical behavior and damage mode of the specimen during loading were analyzed.The results of numerical simulation are compared and analyzed with those of relevant experiments.The feasibility and correctness of the micromodel established in this study for analyzing the micromechanics of lightweight aggregate concrete materials are verified.Image processing technology has a broad application prospect in the field of concrete mesoscopic damage analysis.

Numerical Simulation of Dam-Break Flows Using Radial Basis Functions: Application to Urban Flood Inundation

Dam-break flows pose significant threats to urban areas due to their potential for causing rapid and extensive flooding. Traditional numerical methods for simulating these events struggle with complex urban landscapes. This paper presents an alternative approach using Radial Basis Functions to simulate dam-break flows and their impact on urban flood inundation. The proposed method adapts a new strategy based on Particle Swarm Optimization for variable shape parameter selection on meshfree formulation to enhance the numerical stability and convergence of the simulation. The method’s accuracy and efficiency are demonstrated through numerical experiments, including well-known partial and circular dam-break problems and an idealized city with a single building, highlighting its potential as a valuable tool for urban flood risk management.

MicroMagnetic.jl:A Julia package for micromagnetic and atomistic simulations with GPU support

MicroMagnetic.jl is an open-source Julia package for micromagnetic and atomistic simulations.Using the features of the Julia programming language,MicroMagnetic.jl supports CPU and various GPU platforms,including NVIDIA,AMD,Intel,and Apple GPUs.Moreover,MicroMagnetic.jl supports Monte Carlo simulations for atomistic models and implements the nudged-elastic-band method for energy barrier computations.With built-in support for double and single precision modes and a design allowing easy extensibility to add new features,MicroMagnetic.jl provides a versatile toolset for researchers in micromagnetics and atomistic simulations.

Numerical Simulation and Parallel Computing of Acoustic Wave Equation in Isotropic-Heterogeneous Media

In this paper,we consider the numerical implementation of the 2D wave equation in isotropic-heterogeneous media.The stability analysis of the scheme using the von Neumann stability method has been studied.We conducted a study on modeling the propagation of acoustic waves in a heterogeneous medium and performed numerical simulations in various heterogeneous media at different time steps.Developed parallel code using Compute Unified Device Architecture(CUDA)technology and tested on domains of various sizes.Performance analysis showed that our parallel approach showed significant speedup compared to sequential code on the Central Processing Unit(CPU).The proposed parallel visualization simulator can be an important tool for numerous wave control systems in engineering practice.

Molecular simulation study of the stabilization process of NEPE propellant

In this reported study, the density functional theory(DFT) was used at the(U)B3LYP/6-311G(d,p) level to investigate the stabilization process of the nitrate ester plasticized polyether propellant(NEPE). Molecular simulations were conducted of the reaction that generates NO_(2), the autocatalytic and aging reaction triggered by the NO_(2), and the nitrogen dioxide absorption reaction of the stabilizers during the propellent stabilization process. These simulations were derived using the transition-state theory(TST)and variational transition-state theory(VTST). The simulation results suggested that the stabilization of the NEPE propellant consisted of three stages. First, heat and NO_(2) were generated during the denitrification reaction of nitroglycerine(NG) and 1,2,4-butanetriol trinitrate(BTTN) in the NEPE propellant.Second, nitroso products were generated by the reactions of N-Methyl-4-nitroaniline(MNA) and 2-nitrodiphenylamine(2NDPA) with NO_(2). Third, the stabilizers were exhausted and the autocatalytic reaction of NG and BTTN and the aging reaction of polyethylene glycol(PEG) were triggered by the heat and NO_(2)generated in the first stage. By comparing the energy barriers of the various reactions, it was found that the NO_(2)generated from the denitrification reaction significantly reduced the reaction energy barrier to 105.56-126.32 kJ/mol, also increased the reaction rate constant, and decreased the thermal stability and energetic properties of the NEPE propellant. In addition, the NO_(2)also weakened the mechanical properties of the NEPE propellant by attacking the-CH2groups and the O atoms in the PEG molecular chain. The energy barriers of the reactions of MNA and 2NDPA with NO_(2)(94.61-133.61 k J/mol) were lower than those of the autocatalytic and decomposition reactions of NG, BTTN, and the aging reactions of PEG(160.30-279.46 kJ/mol). This indicated that, by eliminating NO_(2), the stabilizer in the NEPE propellant can effectively prevent NO_(2)from reacting with the NG, BTTN, and PEG in the NEPE propellant. Consequently, this would help maintain the energy and mechanical properties of the NEPE propellant, thereby improving its thermal stability.

Three-dimensional finite element simulation and reconstruction of jointed rock models using CT scanning and photogrammetry

The geometry of joints has a significant influence on the mechanical properties of rocks.To simplify the curved joint shapes in rocks,the joint shape is usually treated as straight lines or planes in most laboratory experiments and numerical simulations.In this study,the computerized tomography (CT) scanning and photogrammetry were employed to obtain the internal and surface joint structures of a limestone sample,respectively.To describe the joint geometry,the edge detection algorithms and a three-dimensional (3D) matrix mapping method were applied to reconstruct CT-based and photogrammetry-based jointed rock models.For comparison tests,the numerical uniaxial compression tests were conducted on an intact rock sample and a sample with a joint simplified to a plane using the parallel computing method.The results indicate that the mechanical characteristics and failure process of jointed rocks are significantly affected by the geometry of joints.The presence of joints reduces the uniaxial compressive strength (UCS),elastic modulus,and released acoustic emission (AE) energy of rocks by 37%–67%,21%–24%,and 52%–90%,respectively.Compared to the simplified joint sample,the proposed photogrammetry-based numerical model makes the most of the limited geometry information of joints.The UCS,accumulative released AE energy,and elastic modulus of the photogrammetry-based sample were found to be very close to those of the CT-based sample.The UCS value of the simplified joint sample (i.e.38.5 MPa) is much lower than that of the CT-based sample (i.e.72.3 MPa).Additionally,the accumulative released AE energy observed in the simplified joint sample is 3.899 times lower than that observed in the CT-based sample.CT scanning provides a reliable means to visualize the joints in rocks,which can be used to verify the reliability of photogrammetry techniques.The application of the photogrammetry-based sample enables detailed analysis for estimating the mechanical properties of jointed rocks.

Chemical simulation teaching system based on virtual reality and gesture interaction

Background Most existing chemical experiment teaching systems lack solid immersive experiences,making it difficult to engage students.To address these challenges,we propose a chemical simulation teaching system based on virtual reality and gesture interaction.Methods The parameters of the models were obtained through actual investigation,whereby Blender and 3DS MAX were used to model and import these parameters into a physics engine.By establishing an interface for the physics engine,gesture interaction hardware,and virtual reality(VR)helmet,a highly realistic chemical experiment environment was created.Using code script logic,particle systems,as well as other systems,chemical phenomena were simulated.Furthermore,we created an online teaching platform using streaming media and databases to address the problems of distance teaching.Results The proposed system was evaluated against two mainstream products in the market.In the experiments,the proposed system outperformed the other products in terms of fidelity and practicality.Conclusions The proposed system which offers realistic simulations and practicability,can help improve the high school chemistry experimental education.

Image processing based three-dimensional model reconstruction for cross-platform numerical simulation

Numerical simulation is the most powerful computational and analysis tool for a large variety of engineering and physical problems.For a complex problem relating to multi-field,multi-process and multi-scale,different computing tools have to be developed so as to solve particular fields at different scales and for different processes.Therefore,the integration of different types of software is inevitable.However,it is difficult to perform the transfer of the meshes and simulated results among software packages because of the lack of shared data formats or encrypted data formats.An image processing based method for three-dimensional model reconstruction for numerical simulation was proposed,which presents a solution to the integration problem by a series of slice or projection images obtained by the post-processing modules of the numerical simulation software.By means of mapping image pixels to meshes of either finite difference or finite element models,the geometry contour can be extracted to export the stereolithography model.The values of results,represented by color,can be deduced and assigned to the meshes.All the models with data can be directly or indirectly integrated into other software as a continued or new numerical simulation.The three-dimensional reconstruction method has been validated in numerical simulation of castings and case studies were provided in this study.

Simulation Research on the Effect of Spreading Process Parameters on the Quality of Lunar Regolith Powder Bed in Additive Manufacturing

Lunar surface additive manufacturing with lunar regolith is a key step in in-situ resource utilization.The powder spreading process is the key process,which has a major impact on the quality of the powder bed and the precision of molded parts.In this study,the discrete element method(DEM)was adopted to simulate the powder spreading process with a roller.The three powder bed quality indicators,including the molding layer offset,voidage fraction,and surface roughness,were established.Besides,the influence of the three process parameters,which are roller’s translational speed,rotational speed,and powder spreading layer thickness on the powder bed quality indicators was also analyzed.The results show that with the reduction of the powder spreading layer thickness and the increase of the rotational speed,the offset increased significantly;when the translational speed increased,the offset first increased and then decreased,which resulted in an extreme value;with the increase of the layer thickness and the decrease of the translational speed,the values for voidage fraction and surface roughness significantly reduced.The powder bed quality indicators were adopted as the optimization objective,and the multi-objective parameter optimization was carried out.The predicted optimal powder spreading parameters and powder bed quality indicators were then obtained.Moreover,the optimal values were then verified.This study can provide informative guidance for in-situ manufacturing at the moon in future deep space exploration missions.

Analysis of seismic disaster failure mechanism and dam-break simulation of high arch dam

Based on a Chinese national high arch dam located in a meizoseismal region, a nonlinear numerical analysis model of the damage and failure process of a dam-foundation system is established by employing a 3-D deformable distinct element code(3DEC) and its re-development functions. The proposed analysis model considers the dam-foundation-reservoir coupling effect, infl uence of nonlinear contact in the opening and closing of the dam seam surface and abutment rock joints during strong earthquakes, and radiation damping of far fi eld energy dissipation according to the actual workability state of an arch dam. A safety assessment method and safety evaluation criteria is developed to better understand the arch dam system disaster process from local damage to ultimate failure. The dynamic characteristics, disaster mechanism, limit bearing capacity and the entire failure process of a high arch dam under a strong earthquake are then analyzed. Further, the seismic safety of the arch dam is evaluated according to the proposed evaluation criteria and safety assessment method. As a result, some useful conclusions are obtained for some aspects of the disaster mechanism and failure process of an arch dam. The analysis method and conclusions may be useful in engineering practice.

Photolithography Process Simulation for Integrated Circuits and Microelectromechanical System Fabrication

Simulations of photoresist etching,aerial image,exposure,and post-bake processes are integrated to obtain a photolithography process simulation for microelectromechanical system（MEMS） and integrated circuit（IC） fabrication based on three-dimensional （3D） cellular automata（CA）. The simulation results agree well with available experimental results. This indicates that the 3D dynamic CA model for the photoresist etching simulation and the 3D CA model for the post-bake simulation could be useful for the monolithic simulation of various lithography processes. This is determined to be useful for the device-sized fabrication process simulation of IC and MEMS.

Material driven workability simulation by FEM including 3D processing maps for magnesium alloy

The three-dimensional （3D） processing maps considering strain based on the two-dimensional （2D） processing maps proposed by PRASAD can describe the distribution of the efficiency of power dissipation and flow instability regions at various temperatures, strain rates and strains, which exhibit intrinsic workability related to material itself. Finite element （FE） simulation can obtain the distribution of strain, strain rate, temperature and die filling status, which indicates state-of-stress （SOS） workability decided by die shape and different processing conditions. On the basis of this, a new material driven analysis method for hot deformation was put forward by the combination of FE simulation with 3D processing maps, which can demonstrate material workability of the entire hot deformation process including SOS workability and intrinsic workability. The hot forging process for hard-to-work metal magnesium alloy was studied, and the 3D thermomechanical FE simulation including 3D processing maps of complex hot forging spur bevel gear was first conducted. The hot forging experiments were carried out. The results show that the new method is reasonable and suitable to determine the aoorooriate nrocess narameters.

Molecular Dynamic Simulation on the Absorbing Process of Isolating and Coating of α-olefin Drag Reducing Polymer

The absorbing process in isolating and coating process of α-olefin drag reducing polymer was studied by molecular dynamic simulation method, on basis of coating theory of α-olefin drag reducing polymer particles with polyurethane as coating material. The distributions of sodium laurate, sodium dodeeyl sulfate, and sodium dodeeyl benzene sulfonate on the surface of α-olefin drag reducing polymer particles were almost the same, but the bending degrees of them were obviously different. The bending degree of SLA molecules was greater than those of the other two surfactant molecules. Simulation results of absorbing and accumulating structure showed that, though hydrophobie properties of surfactant molecules were almost the same, water density around long chain sulfonate sodium was bigger than that around alkyl sulfate sodium. This property goes against useful absorbing and accumulating on the surface of α-olefin drag reducing polymer particles; simulation results of interactions of different surfactant and multiple hydroxyl compounds on surface of particles showed that, interactions of different surfaetant and one kind of multiple hydroxyl compound were similar to those of one kind of surfaetant and different multiple hydroxyl compounds. These two contrast types of interactions also exhibited the differences of absorbing distribution and closing degrees to surface of particles. The sequence of closing degrees was derived from simulation; control step of addition polymerization interaction in coating process was absorbing mass transfer process, so the more closed to surface of particle the multiple hydroxyl compounds were, the easier interactions With isoeyanate were. Simulation results represented the compatibility relationship between surfactant and multiple hydroxyl compounds. The isolating and coating processes of α-olefin drag reducing polymer were further understood on molecule and atom level through above simulation research, and based on the simulation, a referenced theoretical basis was provided for practical optimal selection and experimental preparation of α-olefin drag reducing polymer particles suspension isolation agent.

基于Plant Simulation的航空发动机维修生产线工序优化

针对航空发动机维修生产线生产效率较低的问题,考虑航空发动机大修工序间关系复杂度高、动态性及工时不确定等特点,采用Plant Simulation生产系统仿真软件对大修生产线建立仿真模型。在对工序流程、维修区域和工时等分析的基础上,通过仿真模型对26个维修工序工位从堵塞时间、工作时间和等待时间3个方面进行分析,识别得出叶轮性能恢复、涡轮转子性能恢复和试车3个瓶颈工序。对不同影响因素开展仿真优化与敏感性分析,得出瓶颈工序的优化措施及效果评价。

NUMERICAL SIMULATION OF THE THERMO-MECHANICAL PROCESS FOR BEAM BLANK CONTINUOUS CASTING

The aim of this study was to simulate the solidification process of beam blank continuous casting, and then find the reasons for the typical defects of the beam blank. A two-dimensional transient coupled finite element model has been developed to compute the temperature and stress profile in beam blank continuous casting. The enthalpy method was used in the heat conduction equation. The thermo-mechanical property in the mushy zone was taken into consideration in this calculation. It is shown that at the mold exit the thickness of the shell had its maximum value at the flange tip and its minimum value at the fillet. The temperature had a great fluctuation on the surface of the beam blank in the secondary cooling zone. At the unbending point, the surface temperature of the web was in the brittleness temperature range under the present condition. To ensure the quality, it is necessary to weaken the intensity of secondary cooling. At the mold exit the equivalent stress and strain have higher values at the flange tip and at the web. From the spray 1 to the unbending point, the maximum values of stress and strain gradually moved to the internal section of the flange tip and the web. However, whenever, there were bigger stress and strain values near the flange tip and the web than in the other parts, it must be very easy to generate cracks at those positions. Now, online verification of this simulation has been developed, which has proved to be very useful and efficient to instruct the practical production of beam blank continuous casting.

A simulation study of microstructure evolution during solidification process of liquid metal Ni

A molecular dynamics simulation study has been performed for the microstructure evolution in a liquid metal Ni system during crystallization process at two cooling rates by adopting the embedded atom method （EAM） model potential. The bond-type index method of Honeycutt-Andersen （HA） and a new cluster-type index method （CTIM-2） have been used to detect and analyse the microstructures in this system. It is demonstrated that the cooling rate plays a critical role in the microstructure evolution： below the crystallization temperature Tc, the effects of cooling rate are very remarkable and can be fully displayed. At different cooling rates of 2.0 × 10^13 K·s^-1 and 1.0 × 10^12 K·s^-1, two different kinds of crystal structures are obtained in the system. The first one is the coexistence of the hcp （expressed by （12 0 0 0 6 6） in CTIM-2） and the fcc （12 0 0 0 12 0） basic clusters consisting of 1421 and 1422 bond-types, and the hcp basic cluster becomes the dominant one with decreasing temperature, the second one is mainly the fcc （12 0 0 0 12 0） basic clusters consisting of 1421 bond-type, and their crystallization temperatures Tc would be 1073 and 1173 K, respectively.

A wind tunnel simulation of the dynamic processes involved in sand dune formation on the western coast of Hainan Island

The western coast of Hainan Island exhibits a savanna landscape. Many types of sand dunes, including transverse dune ridges, longitudinal dune ridges, elliptical dunes, coppice dunes, and climbing dunes, are widely distributed in the coastal zone. In winter, high-frequency and high-energy NE winds （dominant winds） are prevalent, with a resultant drift direction （RDD） of S35.6°W. In spring, low-frequency and low-energy SW secondary winds prevail, with a RDD of N25.1°E. Wind tunnel simulations revealed that the airflow over the dune surface is the main factor controlling the erosion and deposition patterns of dune surfaces and the morphological development of dunes. In the region＇s bidirectional wind environment, with two seasonally distinct energy levels, the airflow over the surface of elliptical dunes, barchan dunes, and transverse dune ridges will exhibit a transverse pattern, whereas the airflow over longitudinal dunes ridges exhibits a lateral pattern and that over climbing dunes exhibits a climbing-circumfluent pattern. These patterns represent different dynamic processes. The coastal dunes on the western coast of Hainan Island are influenced by factors such as onshore winds, sand sources, coastal slopes, rivers, and forest shelter belts. The source of the sand that supplements these dunes particularly influences the development pattern： when there is more sand, the pattern shows positive equilibrium deposition between dune ridges and dunes; otherwise, it shows negative equilibrium deposition. The presence or absence of forest shelter belts also influences deposition and dune development patterns and transformation of dune forms. Coastal dunes and inland desert dunes experience similar dynamic processes, but the former have more diversified shapes and more complex formation mechanisms.

Modeling, Simulation and Optimization of a Whole Industrial Catalytic Naphtha Reforming Process on Aspen Plus Platform

A new 18-lump kinetic model for naphtha catalytic reforming reactions is discussed. By developing this model as a user module, a whole industrial continuous catalytic reforming process is simulated on Aspen plus plat-form. The technique utilizes the strong databases, complete sets of modules, and flexible simulation tools of the Aspen plus system and retains the characteristics of the proposed kinetic model. The calculated results are in fair agreement with the actual operating data. Based on the model of the whole reforming process, the process is opti-mized and the optimization results are tested in the actual industrial unit for about two months. The test shows that the process profit increases about 1000yuan·h-1 averagely, which is close to the calculated result.

CFD simulation on the generation of turbidites in deepwater areas: a case study of turbidity current processes in Qiongdongnan Basin, northern South China Sea

Turbidity currents represent a major agent for sediment transport in lakes, seas and oceans. In particu-lar, they formulate the most significant clastic accumulations in the deep sea, which become many of the world＇s most important hydrocarbon reservoirs. Several boreholes in the Qiongdongnan Basin, the north-western South China Sea, have recently revealed turbidity current deposits as significant hydrocarbon res-ervoirs. However, there are some arguments for the potential provenances. To solve this problem, it is es-sential to delineate their sedimentary processes as well as to evaluate their qualities as reservoir. Numerical simulations have been developed rapidly over the last several years, offering insights into turbidity current behaviors, as geologically significant turbidity currents are difficult to directly investigate due to their large scale and often destructive nature. Combined with the interpretation of the turbidity system based on high-resolution 3D seismic data, the paleotophography is acquired via a back-stripping seismic profile integrated with a borehole, i.e., Well A, in the western Qiongdongnan Basin; then a numerical model is built on the basis of this back-stripped profile. After defining the various turbidity current initial boundary conditions, includ-ing grain size, velocity and sediment concentration, the structures and behaviors of turbidity currents are investigated via numerical simulation software ANSYS FLUENT. Finally, the simulated turbidity deposits are compared with the interpreted sedimentary bodies based on 3D seismic data and the potential provenances of the revealed turbidites by Well A are discussed in details. The simulation results indicate that a sedimen-tary body develops far away from its source with an average grain size of 0.1 mm, i.e., sand-size sediment. Taking into account the location and orientation of the simulated seismic line, the consistence between normal forward simulation results and the revealed cores in Well A indicates that the turbidites should have been transported from Vietnam instead of Hainan Island. This interpretation has also been verified by the planar maps of sedimentary systems based on integration of boreholes and seismic data. The identification of the turbidity provenance will benefit the evaluation of extensively distributed submarine fans for hydro-carbon exploration in the deepwater areas.