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.

MONTE CARLO SIMULATION OF RADIATION FIELD OPTIMIZATION FOR MEDICAL LINAC

A method for designing an X-ray flatness filter for medical electron linac is developed. It is used in the optimization process in the electron beam radiation system. Monte Carlo simulation method is used and two examples of real radiation system optimization processes for China-made medical electron linac are provided： 15 MV X- ray system of BJ-20 linac, and 12 MeV electron system of BJ-14. Results are verified by using the traditional method.

Arrhenius relationship and two-step scheme in AF hyperdynamics simulation of diffusion of Mg/Zn interface

The accelerating factor （AF） method is a simple and appropriate way to investigate the atomic long-time deep diffusion at solid-solid interface. In the framework of AF hyperdynamics （HD） simulation, the relationship between the diffusion coefficient along the direction of z-axis which is normal to the Mg/Zn interface and temperature was investigated, and the AF＇s impact on the diffusion constant （D0） and activation energy （Q^＊） was studied. Then, two steps were taken to simulate the atomic diffusion process and the formation of new phases： one for acceleration and the other for equilibration. The results show that： the Arrhenius equation works well for the description of Dz with different accelerating factors; the AF has no effect on the diffusion constant Do in the case of no phase transition; and the relationship between Q＊ and Q conforms to Q^＊=Q/A. Then, the new Arrhenius equation for AFHD is successfully constructed as Dz=Doexp[-Q/（ART）]. Meanwhile, the authentic equilibrium conformations at any dynamic moment can only be reproduced by the equilibration simulation of the HD-simulated configurations. Key words： accelerating factor method; Arrhenius equation; two-steps scheme; Mg/Zn interface; hyperdynamic simulation

Aptix硬件仿真器

随着设计复杂程度的增加,RTL级仿真验证时间大幅度增加,已经成为设计的瓶颈.使用硬件仿真器可以加快仿真速度,提高仿真效率,是设计百万门以上电路的不可或缺的工具.

PIC simulation of plasma properties in the discharge channel of a pulsed plasma thruster with flared electrodes

Plasma in the discharge channel of a pulsed plasma thruster(PPT) with flared electrodes is simulated by a self-developed two-dimensional code. The fully particle-in-cell method with Monte Carlo collision is employed to model the particle movement and collisions and investigate the plasma properties and acceleration process. Temporal and spatial variations of the electron density distribution and the ion velocity between electrodes are calculated and analyzed in detail.The computational results of the electron number density, which is in the order of 1023 m-3,show good agreements with experimental results of a PPT named ADD SIMP-LEX. The ion velocity distributions along the center line of the channel lead to a comprehensive understanding of ions accelerated by electromagnetic field. The electron distributions of PPT with discharge voltages varying from 1300 to 2000 V are compared. The diffusion of electrons presents strong dependency on discharge voltage and implies higher degree of ionization for higher voltage.

Design of a 200-MHz continuous-wave radio frequency quadrupole accelerator for boron neutron capture therapy

A high-intensity continuous-wave(CW) radio frequency quadrupole(RFQ) accelerator is designed for boron neutron capture therapy.The transmission efficiency of a 20-mA proton beam accelerated from 30 keV to 2.5 MeV can reach 98.7% at an operating frequency of 200 MHz.The beam dynamics have a good tolerance to errors.By comparing the high-frequency parameters of quadrilateral and octagonal RFQ cross sections,the quadrilateral structure of the four-vane cavity is selected owing to its multiple advantages,such as a smaller cross section at the same frequency and easy processing.In addition,tuners and undercuts are designed to tune the frequency of the cavity and achieve a flat electric field distribution along the cavity.In this paper,the beam dynamic simulation and electromagnetic design are presented in detail.

Two Improvements on Numerical Simulation of 2-DOF Vortex-Induced Vibration with Low Mass Ratio

Till now, there have been lots of researches on numerical simulation of vortex-induced vibration. Acceptable results have been obtained for fixed cylinders with low Reynolds number. However, for responses of 2-DOF vortex-induced vibration with low mass ratio, the accuracy is not satisfactory, especially for the maximum amplitudes. In Jauvtis and Williamson＇s work, the maximum amplitude of the cylinder with low mass ratio m＊=2.6 can reach as large as 1.5D to be called as the ＂super-upper branch＂, but from current literatures, few simulation results can achieve such value, even fail to capture the upper branch. Besides, it is found that the amplitude decays too fast in the lower branch with the RANS-based turbulence model. The reason is likely to be the defects of the turbulence model itself in the prediction of unsteady separated flows as well as the unreasonable setting of the numerical simulation parameters. Aiming at above issues, a modified turbulence model is proposed in this paper, and the effect of the acceleration of flow field on the response of vortex-induced vibration is studied based on OpenFOAM. By analyzing the responses of amplitude, phase and trajectory, frequency and vortex mode, it is proved that the vortex-induced vibration can be predicted accurately with the modified turbulence model under appropriate flow field acceleration.

QUICK ASSESSMENT METHODOLOGY FOR RELIABILITY OF SOLDER JOINTS IN BALL GRID ARRAY(BGA)ASSEMBLY——PART Ⅱ:RELIABILITY EXPERIMENT AND NUMERICAL SIMULATION

In the present study,a facility,i.e.,a mechanical deflection system (MDS),was established and applied to assess the long-term reliability of the solder joints in plastic ball grid array (BGA) assembly.It was found that the MDS not only quickly assesses the long-term reliability of solder joints within days,but can also mimic similar failure mechanisms in accelerated thermal cycling (ATC) tests. Based on the MDS and ATC reliability experiments,the acceleration factors (AF) were obtained for different reliability testing conditions.Furthermore,by using the creep constitutive relation and fatigue life model developed in part I,a numerical approach was established for the purpose of virtual life prediction of solder joints. The simulation results were found to be in good agreement with the test results from the MDS.As a result,a new reliability assessment methodology was established as an alternative to ATC for the evaluation of long-term reliability of plastic BGA assembly.

Numerical simulation of chorus-driving acceleration of relativistic electrons at extremely low L-shell during geomagnetic storms

During 2018 major geomagnetic storm,relativistic electron enhancements in extremely low L-shell regions(reaching L∼3)have been reported based on observations of ZH-1 and Van Allen probes satellites,and the storm is highly likely to be accelerated by strong whistler-mode waves occurring near very low L-shell regions where the plasmapause was suppressed.It is very interesting to observe the intense chorus-accelerated electrons locating in such low L-shells and filling into the slot region.In this paper,we further perform numerical simulation by solving the two-dimensional Fokker-Planck equation based on the bounce-averaged diffusion rates.Numerical results demonstrate the evolution processes of the chorus-driven electron flux and confirm the flux enhancement in low pitch angle ranges(20◦-50◦)after the wave-particle interaction for tens of hours.The simulation result is consistent with the observation of potential butterfly pitch angle distributions of relativistic electrons from both ZH-1 and Van Allen probes.

Numerical Simulations of the Equations of Particle Motion in the Gas Flow

Under the assumption of considering the gravity and without gravity, two different acceleration models to describe particle’ motion in the gas flow are formulated, respectively. The corresponding numerical simulations of these models do not only show the trend of the velocity of the particle in different density and particle diameter sizes, but also the relationship between the maximum particle velocity and its diameter size.

All-Atom Direct Folding Simulation for Proteins Using the Accelerated Molecular Dynamics in Implicit Solvent Model

We report the results of protein folding （219M, C34, N36, 2KES, 2KHK） by the method of accelerated molecular dynamics （aMD） at room temperature with the implicit solvent model. Starting from the linear structures, these proteins successfully fold to the native structure in a lO0-ns aMD simulation. In contrast, they are failed under the traditional MD simulation in the same simulation time. Then we find that the lowest root mean square deviations of helix structures from the native structures are 0.36 A, 0.63 A, 0.52 A, 1.1 A and 0.78 A. What is more, native contacts, cluster and free energy analyses show that the results of the aMD method are in accordance with the experiment very well. All analyses show that the aMD can accelerate the simulation process, thus we may apply it to the field of computer aided drug designs.

Energy and long-term hygrothermal performance of building enclosures based on dynamic environmental simulation test

Dynamic environmental testing is an effective means to study the energy and long-term hygrothermal performance of building enclosures. Southeast University is designing and building a large-scale dynamic environment simulation testing facility. It can simuhaneously and dynamically simulate temperature, relative humidity, infrared solar radiation, UV radiation, and precipitation. A transformation is needed to predict the energy and long-term hygrothermal performance of building enclosures under real service conditions using data obtained from accelerated tests.

Simulations of Stimulated Raman Scattering in Low-Density Plasmas

Stimulated Raman particle-in-cell （PIC） simulations scattering （SRS） in a low-density The backward stimulated Raman plasma slab is investigated by scattering （B-SRS） dominates initially and erodes the head of the pump wave, while the forward stimulated Raman scattering （F-SRS） subsequently develops and is located at the rear part of the slab. Two-stage electron acceleration may be more efficient due to the coexistence of these two instabilities. The B-SRS plasma wave with low phase velocities can accelerate the background electrons which may be further boosted to higher energies by the F-SRS plasma wave with high phase velocities. The simulations show that the peaks of the main components in both the frequency and wave number spectra occur at the positions estimated from the phase-matching conditions.

Implementation of a Particle Accelerator Beam Dynamics Code on Multi-Node GPUs

Particle accelerators play an important role in a wide range of scientific discoveries and industrial applications. The self-consistent multi-particle simulation based on the particle-in-cell (PIC) method has been used to study charged particle beam dynamics inside those accelerators. However, the PIC simulation is time-consuming and needs to use modern parallel computers for high-resolution applications. In this paper, we implemented a parallel beam dynamics PIC code on multi-node hybrid architecture computers with multiple Graphics Processing Units (GPUs). We used two methods to parallelize the PIC code on multiple GPUs and observed that the replication method is a better choice for moderate problem size and current computer hardware while the domain decomposition method might be a better choice for large problem size and more advanced computer hardware that allows direct communications among multiple GPUs. Using the multi-node hybrid architectures at Oak Ridge Leadership Computing Facility (OLCF), the optimized GPU PIC code achieves a reasonable parallel performance and scales up to 64 GPUs with 16 million particles.

Simulation-based Optimal Design of α-β-γ-δ Filter

The existing third-order tracker known as α-β-γ-δ filter has been used for target tracking and predicting for years. The filter can track the target＇s position and velocity, but not the acceleration. To extend its capability, a new fourth-order target tracker called α-β-γ-δ filter is proposed. The main objective of this study was to find the optimal set of filter parameters that leads to minimum position tracking errors. The tracking errors between using the α-β-γ-δ filter and the α-β-γ-δ filter are compared. As a result, the new filter exhibits significant improvement in position tracking accuracy over the existing third-order filter, but at the expense of computational time in search of the optimal filter. To reduce the computational time, a simulation-based optimization technique via Taguchi method is introduced.

MATHEMATICAL MODEL AND DIGITAL SIMULATION OF TURBOJET ACCELERATION

The acceleration performance of a turbojet is one of its important characteristics. In terms of control system, the engine is a controlled object. The acceleration performance not only depends on the engine, but also on the controller; therefore both the engine and the controller must be combined in a control system to make research for this performance. This paper presents a mathematical model of the turbojet acceleration control system and digital simulation method. Because the engine acceleration is a large variation transient process, the models of the engine and the controller are described by nonlinear equations. The methods of solving nonlinear equations and iterative techniques of calculating acceleration control system are discussed in this paper. The calculated results, as compared with test results, show that the simulation for this system is satisfactory.

Simulation of Plasma Density Gradient Generation by Wakefield and the Effect of Probe Pulse Time-Delay on the Photon Acceleration

Laser Wakefield is produced by ultra-high intensity laser pulse interacting with underdense plasma with special conditions for the laser wavelength and plasma density. In this mechanism, nonlinear forces appear due to the very high amplitudes of the electromagnetic wave and these forces evacuate plasma electrons from the path of the laser pulse leading to very high electron plasma density gradients. Due to the electrostatic forces which result from these density perturbations, the electrons move very fast in oscillatory manner to restore neutrality creating a wake of electron density perturbations behind the laser pulse. Detailed investigation has been dealt with the time-delay between the driver laser pulse and the probe pulse which can affect the production of high plasma gradients needed for photon acceleration process.

The Dynamic Simulation of Rotary Inertia on Light vehicle-Slope I

According to formula we can simulate their driven force and acceleration on the slope.The mechanical formula is used to obtain force and theoretical dynamics in the slope.The driven force decreases when rotation increases.When power increases the acceleration increases.it reduces when its weight raises.It is found that the a will decrease as slope becomes high from 5 to 11°to 22°,which fit the formula too.Meantime as the radius is high from 0.3m to 0.4m to 0.47m a will be low.The needed force will increase as the slope decline becomes big at the same power.

The Simulation on Dynamic of Rotary Inertia and Engine’s Inflamer in Light Vehicle

According to formula we can simulate their driven force and acceleration.The mechanical formula is used to obtain dynamics is used to simulate.The driven force increases when torque increases and tire diameter decreases.We need torque to increase so this is our plan.Acceleration raises when torque raises and it reduces when its weight raises.With the decreasing of radius of road the centripetal acceleration is increasing in the condition of light vehicle.It is that it decreases sluggishly before 0.35m/s2 then it maintains a steep decline to 0.62m/s2 and at last becomes sluggish again.It is valued that the economical efficiency about consumed fuel under different power.In the time of 0.2hr the fuel inflamer inclines sharply first then turns stable.It is the smallest value.Beyond it the fuel maintains a high value all the time.The discharged pollution gas decreases with the decreasing initial temperature.The low initial temperature is good to fuel gas.Meantime the smallest incline range is 300~350K which explains that it is the most save one.

Simulation of high-intense beam transport in electrostatic accelerating column

An electrostatic accelerating column was designed and fabricated by Lanzhou University for an intense DT/D-D neutron generator. In order to achieve a neutron yield of 5.0 × 10^(12)n/s, a deuteron beam of 30 mA,accelerated to 400 kV, and transported in the electrostatic accelerating column smoothly are required. One particle-in-cell code BEAMPATH was used to simulate the beam transport, and the IONB 1.0 code was used to simulate the intense beam envelopes. Emittance growths due to space charge effect and spherical aberration were analyzed. The simulation results show that the accelerating column can transport deuteron beam of 30 mA smoothly and the requirement for the neutron generator is satisfied.