Numerical simulations of earthquake rupture induced by pressure perturbation

The subsurface fluid injection can cause pressure increase within faults,leading to earthquake occurrences.However,the factors controlling earthquake rupture due to pressure perturbation remain poorly understood.To resolve this problem,we simulate the physical processes of earthquake nucleation and rupture on strike-slip faults perturbated by pressure migration based on the slip-weakening law.Multiple kinds of factors,including background stress,fluid injection rates,the area of the pressurized region,fault geometry,and fault friction coefficients,are considered in our simulations.Our simulation results reveal that the ratio of shear stress to normal stress rather than their absolute values controls the rupture behavior.With the large stress ratios,high injection rates,and large pressurized areas,earthquakes are prone to propagate as runaway ruptures.Additionally,faults with large aspect ratios of length to width are also favorable for causing runaway ruptures.In contrast,the factors of fault strike,dip angles and friction coefficients have minor influence on rupture behavior.

Gyro-Landau-fluid simulations of impurity effects on ion temperature gradient driven turbulence transport

The effects of impurities on ion temperature gradient(ITG)driven turbulence transport in tokamak core plasmas are investigated numerically via global simulations of microturbulence with carbon impurities and adiabatic electrons.The simulations use an extended fluid code(ExFC)based on a four-field gyro-Landau-fluid(GLF)model.The multispecies form of the normalized GLF equations is presented,which guarantees the self-consistent evolution of both bulk ions and impurities.With parametric profiles of the cyclone base case,well-benchmarked ExFC is employed to perform simulations focusing on different impurity density profiles.For a fixed temperature profile,it is found that the turbulent heat diffusivity of bulk ions in a quasi-steady state is usually lower than that without impurities,which is contrary to the linear and quasilinear predictions.The evolutions of the temperature gradient and heat diffusivity exhibit a fast relaxation process,indicating that the destabilization of the outwardly peaked impurity profile is a transient state response.Furthermore,the impurity effects from different profiles can obviously influence the nonlinear critical temperature gradient,which is likely to be dominated by linear effects.These results suggest that the improvement in plasma confinement could be attributed to the impurities,most likely through adjusting both heat diffusivity and the critical temperature gradient.

Damage constitutive model of lunar soil simulant geopolymer under impact loading

Lunar base construction is a crucial component of the lunar exploration program,and considering the dynamic characteristics of lunar soil is important for moon construction.Therefore,investigating the dynamic properties of lunar soil by establishing a constitutive relationship is critical for providing a theoretical basis for its damage evolution.In this paper,a split Hopkinson pressure bar(SHPB)device was used to perform three sets of impact tests under different pressures on a lunar soil simulant geopolymer(LSSG)with sodium silicate(Na_(2)SiO_(3))contents of 1%,3%,5%and 7%.The dynamic stressestrain curves,failure modes,and energy variation rules of LSSG under different pressures were obtained.The equation was modified based on the ZWT viscoelastic constitutive model and was combined with the damage variable.The damage element obeys the Weibull distribution and the constitutive equation that can describe the mechanical properties of LSSG under dynamic loading was obtained.The results demonstrate that the dynamic compressive strength of LSSG has a marked strain-rate strengthening effect.Na_(2)SiO_(3) has both strengthening and deterioration effects on the dynamic compressive strength of LSSG.As Na_(2)SiO_(3) grows,the dynamic compressive strength of LSSG first increases and then decreases.At a fixed air pressure,5%Na_(2)SiO_(3) had the largest dynamic compressive strength,the largest incident energy,the smallest absorbed energy,and the lightest damage.The ZWT equation was modified according to the stress response properties of LSSG and the range of the SHPB strain rate to obtain the constitutive equation of the LSSG,and the model’s correctness was confirmed.

Gyrokinetic simulation of low-n Alfvénic modes in tokamak HL-2A plasmas

The turbulence characteristics of plasmas with internal transport barriers in the HL-2A tokamak are analyzed by means of linear gyrokinetic simulations. It is found that turbulence is dominated by the ion temperature gradient(ITG)mode together with large-scale modes characterized by high-frequency electromagnetic fluctuation, which are destabilized by the steep ion temperature gradient in the weak magnetic shear regime. Comparison with solutions of analytical dispersion relations shows that their linear features match well with the beta-induced Alfvén eigenmode branch of the shear Alfvénic spectrum. It is further clarified that the large population of fast ions in these plasmas plays a stabilization role through the dilution mechanism in high-n ITG mode regimes.

Validation of the current and pressure coupling schemes with nonlinear simulations of TAE and analysis on the linear stability of tearing mode in the presence of energetic particles

Both current and pressure coupling schemes have been adopted in the hybrid kinetic–magnetohydrodynamic code CLT-K recently.Numerical equivalences between these two coupling schemes are strictly verified under different approximations.First,when considering only the perturbed distribution function of energetic particles(EPs),the equivalence can be proved analytically.Second,when both the variations of the magnetic field and the EP distribution function are included,the current and pressure coupling schemes numerically produce the same result in the nonlinear simulations.On this basis,the influences of co-/counter-passing and trapped EPs on the linear stabilities of the m/n=2/1 tearing mode(TM)have been investigated(where m and n represent the poloidal and toroidal mode numbers,respectively).The results of scanningβh of EPs show that the co-passing and trapped EPs are found to stabilize the TM,while the counter-passing EPs tend to destabilize the TM.The behind(de)stabilization mechanisms of the TM by EPs are carefully analyzed.Furthermore,after exceeding critical EP betas,the same branch of the high-frequency mode is excited by co-/counterpassing and trapped EPs,which is identified as the m/n=2/1 energetic particle mode.

Simulation Mode Application in Winter Wheat North-moving

[Objective] The aim was to study the simulation mode application in winter wheat north moving.[Method] Based on meteorological data of many years collected in north of Hebei province,and arrangement of multi-location field experiment,mode against cold was constructed considering cold-resistance characteristics of wheat.[Result] Wheat overwintering was forecasted with the help of working mode,results were verified,and the safe northern boundary for wheat overwintering in Hebei was determined in regions below dam in the south of 41°N.[Conclusion] This study provided reference for exploration of dynamic rules of wheat growth,and growth influence by environment and plant technique.

Model test and numerical simulation on the dynamic stability of the bedding rock slope under frequent microseisms

Shake table testing was performed to investigate the dynamic stability of a mid-dip bedding rock slope under frequent earthquakes. Then, numerical modelling was established to further study the slope dynamic stability under purely microseisms and the influence of five factors, including seismic amplitude, slope height, slope angle, strata inclination and strata thickness, were considered. The experimental results show that the natural frequency of the slope decreases and damping ratio increases as the earthquake loading times increase. The dynamic strength reduction method is adopted for the stability evaluation of the bedding rock slope in numerical simulation, and the slope stability decreases with the increase of seismic amplitude, increase of slope height, reduction of strata thickness and increase of slope angle. The failure mode of a mid-dip bedding rock slope in the shaking table test is integral slipping along the bedding surface with dipping tensile cracks at the slope rear edge going through the bedding surfaces. In the numerical simulation, the long-term stability of a mid-dip bedding slope is worst under frequent microseisms and the slope is at risk of integral sliding instability, whereas the slope rock mass is more broken than shown in the shaking table test. The research results are of practical significance to better understand the formation mechanism of reservoir landslides and prevent future landslide disasters.

An Integrated Dynamic Model of Ocean Mining System and Fast Simulation of Its Longitudinal Reciprocating Motion

An integrated dynamic model of China＇s deep ocean mining system is developed and the fast simulation analysis of its longitudinal reciprocating motion operation processes is achieved. The seafloor tracked miner is built as a three-dimensional single-body model with six-degree-of-freedom. The track-terrain interaction is modeled by partitioning the track-terrain interface into a certain number of mesh elements with three mutually perpendicular forces, including the normal force, the longitudinal shear force and the lateral shear force, acting on the center point of each mesh element. The hydrodynamic force of the miner is considered and applied. By considering the operational safety and collection efficiency, two new mining paths for the miner on the seafloor are proposed, which can be simulated with the established single-body dynamic model of the miner. The pipeline subsystem is built as a three-dimensional multi-body discrete element model, which is divided into rigid elements linked by flexible connectors. The flexible connector without mass is represented by six spring-damper elements. The external hydrodynamic forces of the ocean current from the longitudinal and lateral directions are both considered and modeled based on the Morison formula and applied to the mass center of each corresponding discrete rigid element. The mining ship is simplified and represented by a general kinematic point, whose heave motion induced by the ocean waves and the longitudinal and lateral towing motions are considered and applied. By integrating the single-body dynamic model of the miner and the multi-body discrete element dynamic model of the pipeline, and defining the kinematic equations of the mining ship, the integrated dynamic model of the total deep ocean mining system is formed. The longitudinal reciprocating motion operation modes of the total mining system, which combine the active straight-line and turning motions of the miner and the ship, and the passive towed motions of the pipeline, are proposed and simulated with the developed 3D dynamic model. Some critical simulation results are obtained and analyzed, such as the motion trajectories of key subsystems, the velocities of the buoyancy modules and the interaction forces between subsystems, which in a way can provide important theoretical basis and useful technical reference for the practical deep ocean mining system analysis, operation and control.

An agent-based simulation with NetLogo platform to evaluate forward osmosis process(PRO Mode)

Forward osmosis(FO), as an emerging technology, is influenced by different factors such as operating conditions,module characteristics, and membrane properties. The general aim of this study was to develop a suitable(flexible,comprehensive, and convenient to use) computational tool which is able to simulate osmosis through an asymmetric membrane oriented in pressure retarded osmosis(PRO) mode in a wide variety of scenarios. For this purpose, an agent-based model was created in NetLogo platform, which is an easy-to-use application environment with graphical visualization abilities and well suited for modeling a complex system evolving over time. The simulation results were validated with empirical data obtained from literature and a great agreement was observed. The effect of various parameters on process performance was investigated in terms of temperature,cross-flow velocity, length of the module, pure water permeability coefficient, and structural parameter of the membrane. Results demonstrated that the increase in all parameters, except structural parameter of the membrane and the length of module led to the increase of average water flux. Moreover, nine different draw solutes were selected in order to assess the influence of net bulk osmotic pressure difference between the draw solution(DS) and feed solution(FS)(known as the driving force of FO process) on water flux. Based on the findings of this paper, the performance of FO process(PRO mode) can be efficiently evaluated using the NetL ogo platform.

Simulation research of energy management strategy for dual mode plug-in hybrid electrical vehicles

In this paper, a plug-in hybrid electrical vehicle(PHEV) is taken as the research object, and its dynamic performance and economic performance are taken as the research goals. Battery charge-sustaining(CS) period is divided into power mode and economy mode. Energy management strategy designing methods of power mode and economy mode are proposed. Maximum velocity, acceleration performance and fuel consumption are simulated during the CS period in the AVL CRUISE simulation environment. The simulation results indicate that the maximum velocity and acceleration time of the power mode are better than those in the economy mode. Fuel consumption of the economy mode is better than that in the power mode. Fuel consumption of PHEV during the CS period is further improved by using the methods proposed in this paper, and this is meaningful for research and development of PHEV.

An Improved Multi-Objective Hybrid Genetic-Simulated Annealing Algorithm for AGV Scheduling under Composite Operation Mode

This paper presents an improved hybrid algorithm and a multi-objective model to tackle the scheduling problem of multiple Automated Guided Vehicles(AGVs)under the composite operation mode.The multi-objective model aims to minimize the maximum completion time,the total distance covered by AGVs,and the distance traveled while empty-loaded.The improved hybrid algorithm combines the improved genetic algorithm(GA)and the simulated annealing algorithm(SA)to strengthen the local search ability of the algorithm and improve the stability of the calculation results.Based on the characteristics of the composite operation mode,the authors introduce the combined coding and parallel decoding mode and calculate the fitness function with the grey entropy parallel analysis method to solve the multi-objective problem.The grey entropy parallel analysis method is a combination of the grey correlation analysis method and the entropy weighting method to solve multi-objective solving problems.A task advance evaluation strategy is proposed in the process of crossover and mutation operator to guide the direction of crossover and mutation.The computational experiments results show that the improved hybrid algorithm is better than the GA and the genetic algorithm with task advance evaluation strategy(AEGA)in terms of convergence speed and solution results,and the effectiveness of the multi-objective solution is proved.All three objectives are optimized and the proposed algorithm has an optimization of 7.6%respectively compared with the GA and 3.4%compared with the AEGA in terms of the objective of maximum completion time.

Experimental and simulation studies on similitude design method for shock responses of beam-plate coupled structure

The similitude theory helps to understand the physical behaviors of large structures through scaled models. Several papers have studied the similitude of shock issues. However, the dynamic similitude for shock responses of coupled structures is rarely incorporated in open studies. In this paper, scaling laws are derived for the shock responses and spectra of coupled structures. In the presented scaling laws, the geometric distortion and energy loss are considered. The ability of the proposed scaling laws is demonstrated in the simulation and experimental cases. In both cases, the similitude prediction for the prototype's time-domain waveform and spectrum is conducted with the scaled model and scaling laws. The simulation and experimental cases indicate that the predicted shock responses and spectra agree well with those of the prototype, which verifies the proposed scaling laws for predicting shock responses.

Nonlinear evolution and secondary island formation of the double tearing mode in a hybrid simulation

Double tearing modes(DTMs),induced by double current sheet configurations or two neighboring rational surfaces with the same safety factor in tokamaks,are widely observed in solar,space,and fusion plasmas.In this paper,the evolution of DTMs without a guide field is investigated numerically using a hybrid model(electron fluid+ion PIC).The overall evolution processes of DTMs are qualitatively consistent with previous works using other models.The particle dynamics during the evolution of DTMs is analyzed in detail.Behaviors of ions and electrons present different characteristics around the reconnection region which gives rise to Hall effects producing the out-of-plane quadrupole magnetic field.In the explosive reconnection process with interactions between two DTMs islands,the asymmetric drive and the thin current layer feature lead to the emergence of secondary magnetic islands which develop with the late evolution of the DTMs.

A SIMULATION STUDY ON DYNAMIC-SAMPLING-MODE FOR FLUORESCENCE MOLECULAR TOMOGRAPHY

Fluorescence tomography can obtain a sufficient dataset and optimal three-dimensional imageswhen projections are captured over 360◦ by CCD camera. Herein a non-stop dynamic samplingmode for fluorescence tomography is proposed in an attempt to improve the optical measurementspeed of the traditional imaging system and stability of the object to be imaged. A series ofsimulations are carried out to evaluate the accuracy of dataset acquired from the dynamicsampling mode. Reconstruction with the corresponding data obtained in the dynamic-modeprocess is also performed with the phantom. The results demonstrate the feasibility of suchan imaging mode when the angular velocity is set to the appropriate value, thus laying thefoundation for real experiments to verify the superiority in performance of this new imagingmode over the traditional one.

Calibration of a rule-based intelligent network simulation model

This paper is focused on calibration of an intelligent network simulation model （INS1M） with reallife transportation network to analyse the INSIM＇s feasibility in simulating commuters＇ travel choice behaviour under the influence of real-time integrated multimodal traveller information （IMTI）. A transportation network model for the central and western areas of Singapore was simulated in PARAMICS and integrated with INSIM expert system by means of an application programming interface to form the INSIM. Upon calibration, INSIM was able to realistically present complicated scenarios in which real-time IMTI was provided to commuters and the network performance measures being recorded.

Research and Simulation of Parking Resource Sharing Mode in Old Community—Taking Yingmenkou Community as an Example

In recent years, the number of permanent residents in Chengdu, Sichuan Province, China has continued to grow, and the number of private cars keeps rising. A series of problems caused by difficult parking and random parking in old communities are common. This paper will take the parking problem of the old community in Yingmenkou community, Chengdu, China as an example to summarize the current situation and parking problems of the private car parking facilities in these old communities. Based on Chengdu’s current old community parking management policy, parking facility renovation and expansion policy, this paper researches and establishes a parking resource sharing model suitable for old communities. The simulation results show that this model has a small footprint and is easy to install and disassemble. It can be freely spliced and combined units according to the size of the old community and the needs of residents. It is suitable for old communities with no centralized parking lot and high parking demand by residents.

Density of excess modes below the first phonon mode in four-dimensional glasses

Glasses are known to possess low-frequency excess modes beyond the Debye prediction.For decades,it has been assumed that evolution of low-frequency density of excess modes D(ω) with frequency ω follows a power-law scaling:D(ω)~ω~γ.However,it remains debated on the value of γ at low frequencies below the first phonon-like mode in finitesize glasses.Early simulation studies reported γ=4 at low frequencies in two-(2D),three-(3D),and four-dimensional(4D)glasses,whereas recent observations in 2D and 3D glasses suggested γ=3.5 in a lower-frequency regime.It is uncertain whether the low-frequency scaling of D(ω)~ω^(3.5) could be generalized to 4D glasses.Here,we conduct numerical simulation studies of excess modes at frequencies below the first phonon-like mode in 4D model glasses.It is found that the system size dependence of D(ω) below the first phonon-like mode varies with spatial dimensions:D(ω) increases in2D glasses but decreases in 3D and 4D glasses as the system size increases.Furthermore,we demonstrate that the ω^(3.5)scaling,rather than the ω~4 scaling,works in the lowest-frequency regime accessed in 4D glasses,regardless of interaction potentials and system sizes examined.Therefore,our findings in 4D glasses,combined with previous results in 2D and 3D glasses,suggest a common low-frequency scaling of D(ω)~ ω^3.5) below the first phonon-like mode across different spatial dimensions,which would inspire further theoretical studies.

STATIC CORRECTION MODES IN DYNAMICAL SIMULATION OF FLEXIBLE MULTIBODY SYSTEM WITH CLOSED LOOPS

By use of the recursive formulation of flexible multibody system with closed loops,this paper studies the dynamic simulation of complicated multibody system.The conception of static correction modes in structural dynamics is introduced to present the local deformation effects due to reaction forces at kinematic joints.A space four bar mechanism is utilized here as an example to describe the method of selecting static correction modes.Compared with vibration normal modes, static correction modes are demonstrated to be effective in numerical simulation.

Studies on fluid model for numerical simulation of gas discharges in color plasma displays

The fluid models of gas discharge in alternating current plasma display panel (AC PDP) cell are discussed. From the Boltzmann equation, the hydrodynamic equations are derived, but this model consumes much computa- tional time for simulation. The drift-diffusion approximation model and the local field approximation model are ob- tained to simplify the numerical computation, and the approximation conditions of these two models are discussed in detail. The drift-diffusion approximation model gives more satisfactory result for PDP simulation, and the expression of energy balance equation is given completely in this model.

GENERAL FAILURE PROBABILITY SIMULATION AND APPLICATION FOR MULTI-MODE

A general failure probability simulation and deviation evaluation methods were presented for fuzzy safety state and fuzzy failure state. And the corresponding number integral method was simultaneously established. As the distribution of state variable and the membership of the state variable to the fuzzy safety set were normal, the general failure probability of the single failure mode had precise analytic solution, which was used to verify the precision of the presented methods. The results show that the evaluation of the simulation method convergences to the analytic solution with the number increase of the sampling. The above methods for the single failure mode was extended to the multi-mode by the expansion and probability principles. The presented methods were applied to the engineering problem. For the number of significant mode is not too many, the high precision solution can be given by the presented number simulation and number integral methods, which is illustrated by the engineering examples. In addition, the application scope of the methods was discussed.