A High-resolution Simulation of Supertyphoon Rammasun(2014)--Part Ⅰ:Model Verification and Surface Energetics Analysis

A 72-h high-resolution simulation of Supertyphoon Rammasun （2014） is performed using the Advanced Research Weather Research and Forecasting model. The model covers an initial 18-h spin-up, the 36-h rapid intensification （RI） period in the northern South China Sea, and the 18-h period of weakening after landfall. The results show that the model reproduces the track, intensity, structure of the storm, and environmental circulations reasonably well. Analysis of the surface energetics under the storm indicates that the storm＇s intensification is closely related to the net energy gain rate （eg）, defined as the difference between the energy production （PD） due to surface entropy flux and the energy dissipation （Ds） due to surface friction near the radius of maximum wind （RMW）. Before and during the RI stage, the ~：g is high, indicating sufficient energy supply for the storm to intensify. However, the Sg decreases rapidly as the storm quickly intensifies, because the Ds increases more rapidly than the PD near the RMW. By the time the storm reaches its peak intensity, the Ds is about 20% larger than the PD near the RMW, leading to a local energetics deficit under the eyewall. During the mature stage, the PD and Ds can reach a balance within a radius of 86 km from the storm center （about 2.3 times the RMW）. This implies that the local PD under the eyewall is not large enough to balance the Ds, and the radially inward energy transport from outside the eyewall must play an important role in maintaining the storm＇s intensity, as well as its intensification.

High-Resolution Simulation and Analysis of the Mature Structure of a Polar Low over the Sea of Japan on 21 January 1997

This paper presents a high-resolution simulation of a remarkable polar low observed over the Sea of Japan on 21 January 1997 by using a 5-km mesh non-hydrostatic model MRI-NHM (Meteorological Research Institute Non-Hydrostatic Model). A 24-hour simulation starting from 0000 UTC 21 January 1997 successfully reproduced the observed features of the polar low such as the wrapping of western part of an initial E-W orientation vortex, the spiral-shaped bands, the cloud-free 'eye', and the warm core structure at its mature stage. The 'eye' of the simulated polar low was relatively dry, and was associated with a strong downdraft. A thermodynamic budget analysis indicates that the 'warm core' in the 'eye' region was mainly caused by the adiabatic warming associated with the downdraft. The relationship among the condensational diabatic heating, the vertical velocity, the convergence of the moisture flux, and the circulation averaged within a 50 km×50 km square area around the polar low center shows that they form a positive feedback loop, and this loop is not inconsistent with the CISK (Conditional Instability of the Second Kind) mechanism during the developing stage of the polar low.

High-resolution simulation of upper-ocean submesoscale variability in the South China Sea:Spatial and seasonal dynamical regimes

Submesoscale processes in marginal seas usually have complex generating mechanisms,highly dependent on the local background flow and forcing.This numerical study investigates the spatial and seasonal differences of submesoscale activities in the upper ocean of the South China Sea(SCS)and the different dynamical regimes for sub-regions.The spatial and seasonal variations of vertical vorticity,horizontal convergence,lateral buoyancy gradient,and strain rate are analyzed to compare the submesoscale phenomenon within four sub-regions,the northern region near the Luzon Strait(R1),the middle ocean basin(R2),the western SCS(R3),and the southern SCS(R4).The results suggest that the SCS submesoscale processes are highly heterogeneous in space,with different seasonalities in each sub-region.The submesoscale activities in the northern sub-regions(R1,R2)are active in winter but weak in summer,while there appears an almost seasonal anti-phase in the western region(R3)compared to R1 and R2.Interestingly,no clear seasonality of submesoscale features is shown in the southern region(R4).Further analysis of Ertel potential vorticity reveals different generating mechanisms of submesoscale processes in different sub-regions.Correlation analyses also show the vertical extent of vertical velocity and the role of monsoon in generating submesoscale activities in the upper ocean of sub-regions.All these results suggest that the sub-regions have different regimes for submesoscale processes,e.g.,Kuroshio intrusion(R1),monsoon modulation(R2),frontal effects(R3),topography wakes(R4).

A Comparison Study of Tropical Pacific Ocean State Estimation:Low-Resolution Assimilation vs.High-Resolution Simulation

A comparison study is performed to contrast the improvements in the tropical Pacific oceanic state of a low-resolution model respectively via data assimilation and by an increase in horizontal resolution. A low resolution model (LR) (1°lat by 2°lon) and a high-resolution model (HR) (0.5°lat by 0.5°lon) are employed for the comparison. The authors perform 20-yr numerical experiments and analyze the annual mean fields of temperature and salinity. The results indicate that the low-resolution model with data assimilation behaves better than the high-resolution model in the estimation of ocean large-scale features. From 1990 to 2000, the average of HR's RMSE (root-mean-square error) relative to independent Tropical Atmosphere Ocean project (TAO) mooring data at randomly selected points is 0.97℃ compared to a RMSE of 0.56℃ for LR with temperature assimilation. Moreover, the LR with data assimilation is more frugal in computation. Although there is room to improve the high-resolution model, the low-resolution model with data assimilation may be an advisable choice in achieving a more realistic large-scale state of the ocean at the limited level of information provided by the current observational system.

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

Relative Impacts of Sea Ice Loss and Atmospheric Internal Variability on the Winter Arctic to East Asian Surface Air Temperature Based on Large-Ensemble Simulations with NorESM2

To quantify the relative contributions of Arctic sea ice and unforced atmospheric internal variability to the “warm Arctic, cold East Asia”(WACE) teleconnection, this study analyses three sets of large-ensemble simulations carried out by the Norwegian Earth System Model with a coupled atmosphere–land surface model, forced by seasonal sea ice conditions from preindustrial, present-day, and future periods. Each ensemble member within the same set uses the same forcing but with small perturbations to the atmospheric initial state. Hence, the difference between the present-day(or future) ensemble mean and the preindustrial ensemble mean provides the ice-loss-induced response, while the difference of the individual members within the present-day(or future) set is the effect of atmospheric internal variability. Results indicate that both present-day and future sea ice loss can force a negative phase of the Arctic Oscillation with a WACE pattern in winter. The magnitude of ice-induced Arctic warming is over four(ten) times larger than the ice-induced East Asian cooling in the present-day(future) experiment;the latter having a magnitude that is about 30% of the observed cooling. Sea ice loss contributes about 60%(80%) to the Arctic winter warming in the present-day(future) experiment. Atmospheric internal variability can also induce a WACE pattern with comparable magnitudes between the Arctic and East Asia. Ice-lossinduced East Asian cooling can easily be masked by atmospheric internal variability effects because random atmospheric internal variability may induce a larger magnitude warming. The observed WACE pattern occurs as a result of both Arctic sea ice loss and atmospheric internal variability, with the former dominating Arctic warming and the latter dominating East Asian cooling.

Verification and Validation of High-Resolution Inviscid and Viscous Conical Nozzle Flows

Capturing elaborated flow structures and phenomena is required for well-solved numerical flows.The finite difference methods allow simple discretization of mesh and model equations.However,they need simpler meshes,e.g.,rectangular.The inverse Lax-Wendroff(ILW)procedure can handle complex geometries for rectangular meshes.High-resolution and high-order methods can capture elaborated flow structures and phenomena.They also have strong mathematical and physical backgrounds,such as positivity-preserving,jump conditions,and wave propagation concepts.We perceive an effort toward direct numerical simulation,for instance,regarding weighted essentially non-oscillatory(WENO)schemes.Thus,we propose to solve a challenging engineering application without turbulence models.We aim to verify and validate recent high-resolution and high-order methods.To check the solver accuracy,we solved vortex and Couette flows.Then,we solved inviscid and viscous nozzle flows for a conical profile.We employed the finite difference method,positivity-preserving Lax-Friedrichs splitting,high-resolution viscous terms discretization,fifth-order multi-resolution WENO,ILW,and third-order strong stability preserving Runge-Kutta.We showed the solver is high-order and captured elaborated flow structures and phenomena.One can see oblique shocks in both nozzle flows.In the viscous flow,we also captured a free-shock separation,recirculation,entrainment region,Mach disk,and the diamond-shaped pattern of nozzle flows.