Simulating the evolution of focused waves by a two-layer Boussinesq-type model

Accurate simulation of the evolution of freak waves by the wave phase focusing method requires accurate linear and nonlinear properties,especially in deep-water conditions.In this paper,we analyze the ability to simulate deep-water focused waves of a two-layer Boussinesq-type(BT)model,which has been shown to have excellent linear and nonlinear performance.To further improve the numerical accuracy and stability,the internal wavegenerated method is introduced into the two-layer Boussinesq-type model.Firstly,the sensitivity of the numerical results to the grid resolution is analyzed to verify the convergence of the model;secondly,the focused wave propagating in two opposite directions is simulated to prove the symmetry of the numerical results and the feasibility of the internal wave-generated method;thirdly,the limiting focused wave condition is simulated to compare and analyze the wave surface and the horizontal velocity of the profile at the focusing position,which is in good agreement with the measured values.Meanwhile the simulation of focused waves in very deep waters agrees well with the measured values,which further demonstrates the capability of the two-layer BT model in simulating focused waves in deep waters.

A new model simulating the development of gas condensate reservoirs

A new simulation model for the development of gas condensate reservoirs is introduced based on the influence that phase change,non-Darcy flow,and capillary pressure have on the production of gas condensates.The model predicts well performance,including bottom-hole pressure,oil/gas production rate,oil/gas recovery,gaseoil ratio,and the change in produced fluid composition.It also calculates dynamic characters,such as the change of pressure field and oil/gas saturation field during the development of gas condensate reservoirs.The model is applicable to different boundary conditions(both constant-pressure and sealed boundary)and different production modes(both constant-pressure and constant-volume production modes).Model validation attempted using numerical simulation results for sealed boundary conditions with constant-pressure production mode has shown a relatively good match,proving its validity.For constant-pressure boundary conditions with constant-volume production mode,four stages are defined according to the dynamic behavior of production performance in the development of gas condensate reservoirs.

Absorbing Boundary Conditions for Simulating Water Waves Near Solid Bodies

The objective of this paper is to present a new method for designing absorbing or non-reflective boundary conditions (ABC) or (NRBC), illustrated by the case study of the modelling of a solid body in water, specifically the capillary gravity waves generated by its motion at the surface. The study analyses the flow of an inviscid, barotropic, and compressible fluid around the stationary solid body. The dynamic behaviour of the fluid is analysed using a two-dimensional coupled Neumann-Kelvin model extended with capillarity and inertia terms. For computational purposes, it is necessary to truncate the unbounded spatial domain with artificial boundaries and then introduce appropriate absorbing boundary conditions. The propagation of short wavelength waves in a convective fluid medium with significant differences in properties between the interior and the surface of the fluid presents a number of difficulties in the design of these conditions. The results are illustrated numerically and commented upon.

An approach for simulating the air brake system of long freight trains based on fluid dynamics

Air brake systems are critical equipment for railway trains, which affects the running safety of the trains significantly. To study air braking characteristics of long freight trains, an approach for simulating air brake systems based on fuid dynamics theory was proposed. The structures and working mechanisms of locomotive and wagon air brakes are introduced, and mathematical models of the pipes, brake valves, reservoirs or chambers, cylinders, etc., are presented.Besides, the dynamic motions of parts in the main valve are considered. The simulation model of the whole air brake system is then formulated, and the solving method based on the finite-difference method is used. New efficient pipe boundary conditions without iterations are developed for brake pipes and branch pipes, which can achieve higher computational efficiency. The proposed approach for simulating the air brake system is validated by comparing with published measured data. Simulation results of different train formations indicate that models that consider the dynamic behavior of brake pipes are recommended for predicting the characteristics of long trains under service braking conditions.

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

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

Simulating the strong ground motion of the 2022 MS6.8 Luding,Sichuan,China Earthquake

Stochastic finite-fault simulations are effective for simulating ground motions and are widely used in engineering to determine the impacts of ground motion and develop relevant predictive equations.In this study,the source,path,and site amplification coefficient of western Sichuan Province,China,and stochastic finite-fault simulations were used to simulate the acceleration time series,Fourier amplitude spectra,and 5%damped response spectra of 28 strong-motion stations with rupture distances within 300 km of the 2022 MS6.8 Luding earthquake.The simulation results of 14 stations at rupture distances of 45-185 km match the observation.However,the simulation results of 3 near-and 6 far-field stations at rupture distances of 12-36 km and 222-286 km,respectively,were obviously deviated from the observations.Simulation results of the near-field stations are larger than the observations at high frequencies(>6 Hz).The discrepancy likely comes from the nonlinear site effect of near-field stations,which reduced the site amplification at high frequencies.Simulation result of the far-field stations is smaller than the observation at frequencies above 1 Hz.As these stations are located close to the Longmenshan Fault Zone(LFZ),thus,we obtained a new quality factor(Q)from data of historical events and stations located around LFZ.Using the new Q value,the discrepancies of the high-frequency simulation results of the far-field stations were corrected.This result indicated that the laterally varying Q values can be used to address the impact of strong crustal lateral heterogeneity on simulation.

Evaluation of Performance of Polar WRF Model in Simulating Precipitation over Qinghai-Tibet Plateau

Considering the complex topographic forcing and large cryosphere concentration,the present study utilized the polar-optimized WRF model(Polar WRF)to conduct downscaling simulations over the Qinghai-Tibet Plateau(TP)and its surrounding regions.Multi-group experiments with the 10 km horizontal resolution are used to evaluate the modeling of precipitation.Firstly,on the basis of the model ground surface properties upgrade and the optimized Noah-MP,the“better-performing”configuration suite for modeling precipitation is comprehensively examined.Various model parameters such as nudging options,five cumulus parameterization schemes,two planetary boundary layer schemes,and six microphysics schemes are investigated to further refine the Polar WRF configuration.Moreover,the precipitation simulation for a full calendar year is compared with multiple reanalyses and observations.The simulations demonstrate that the Polar WRF model successfully captures the general features of precipitation over this region and is sensitive to model parameters.Based on the results,it is recommended to use grid nudging with q intensity coefficient of 0.0002,the multi-scale kain-fritsch cumulus parameterization,the Yonsei University boundary layer scheme,and the Morrison 2-mom microphysics with reduced default droplet concentration value of 100 cm-3.Overall,the model performance is better than the ERA-interim and TRMM 3b42.It is comparable to,and in some cases slightly better than,the CRA-Land,especially in the prediction for the western part of the plateau where in situ observations are limited,and the cryosphere-atmosphere interaction is more pronounced.

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

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

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

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

基于CarSim/Simulink的计算机控制系统综合实验设计

针对计算机控制系统课程中控制算法设计与仿真实现,设计了基于CarSim/Simulink的控制综合实验方案。以4轮独立驱动小车为对象,要求学生熟悉其基本工作原理,并在非线性动力学模型基础上,进行局部线性近似,设计反馈控制器,模拟小车控制。使用CarSim仿真软件可直观地模拟调整4轮独立驱动小车的质量、质心位置和转动惯量等,实现个性化选择,搭建车辆模型;与Simulink接口互通,在Simulink中完成控制算法。实验设计旨在使学生夯实计算机控制系统的理论基础,并逐步过渡到控制工程,理论联系实际,培养其解决复杂工程问题的能力。

基于Simulink的2PSK调制与解调系统的仿真分析

分析2PSK调制和解调的原理,采用Matlab环境下的Simulink仿真工具箱实现2PSK调制及解调系统的仿真建模与分析。仿真实验表明,仿真结果和理论分析内容一致,验证了仿真模型的正确性。该方法有助于相关研究者直观地理解相移调制与解调系统的组成及工作原理,对实际工程应用也有一定的借鉴作用。

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).

电动汽车电控液压制动系统CarSim/Simulink联合仿真研究

针对某款国产电动汽车的电控液压制动系统,文章基于滑移率的比例-积分-微分(PID)控制提出了防抱死制动系统(ABS)模型,以改善车辆制动性能,提高在对开路面上的制动有效性和安全性。首先,使用Simulink软件建立液压电控制动系统的动力学模型;然后,基于滑移率设计PID控制器,并通过1/4车辆模型验证其有效性;最后,利用S函数将ABS控制器的Simulink模型输入CarSim平台中,开展CarSim/Simulink双平台联合控制仿真。结果表明,基于滑移率的PID控制的ABS对比无ABS的车辆,制动性能更加优越,在对开路面上制动稳定性更加明显。

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.