Numerical simulation of microwave-induced cracking and melting of granite based on mineral microscopic models

This study introduces a coupled electromagnetic–thermal–mechanical model to reveal the mechanisms of microcracking and mineral melting of polymineralic rocks under microwave radiation.Experimental tests validate the rationality of the proposed model.Embedding microscopic mineral sections into the granite model for simulation shows that uneven temperature gradients create distinct molten,porous,and nonmolten zones on the fracture surface.Moreover,the varying thermal expansion coefficients and Young's moduli among the minerals induce significant thermal stress at the mineral boundaries.Quartz and biotite with higher thermal expansion coefficients are subjected to compression,whereas plagioclase with smaller coefficients experiences tensile stress.In the molten zone,quartz undergoes transgranular cracking due to theα–βphase transition.The local high temperatures also induce melting phase transitions in biotite and feldspar.This numerical study provides new insights into the distribution of thermal stress and mineral phase changes in rocks under microwave irradiation.

Numerical and theoretical study of large-scale failure of strata overlying sublevel caving mines with steeply dipping discontinuities

The deformation and fracture evolution mechanisms of the strata overlying mines mined using sublevel caving were studied via numerical simulations.Moreover,an expression for the normal force acting on the side face of a steeply dipping superimposed cantilever beam in the surrounding rock was deduced based on limit equilibrium theory.The results show the following:(1)surface displacement above metal mines with steeply dipping discontinuities shows significant step characteristics,and(2)the behavior of the strata as they fail exhibits superimposition characteristics.Generally,failure first occurs in certain superimposed strata slightly far from the goaf.Subsequently,with the constant downward excavation of the orebody,the superimposed strata become damaged both upwards away from and downwards toward the goaf.This process continues until the deep part of the steeply dipping superimposed strata forms a large-scale deep fracture plane that connects with the goaf.The deep fracture plane generally makes an angle of 12°-20°with the normal to the steeply dipping discontinuities.The effect of the constant outward transfer of strata movement due to the constant outward failure of the superimposed strata in the metal mines with steeply dipping discontinuities causes the scope of the strata movement in these mines to be larger than expected.The strata in the metal mines with steeply dipping discontinuities mainly show flexural toppling failure.However,the steeply dipping structural strata near the goaf mainly exhibit shear slipping failure,in which case the mechanical model used to describe them can be simplified by treating them as steeply dipping superimposed cantilever beams.By taking the steeply dipping superimposed cantilever beam that first experiences failure as the key stratum,the failure scope of the strata(and criteria for the stability of metal mines with steeply dipping discontinuities mined using sublevel caving)can be obtained via iterative computations from the key stratum,moving downward toward and upwards away from the goaf.

Geopolymer-based modification of blasting sealing materials and optimization of blasting block size in coal seams of open pit mines

This research proposes the utilization of a geopolymer-based blasting sealing material to improve the profitability of coal sales and reduce the rate of coal fragmentation during blasting in open pit mines.The study first focused on optimizing the strength of the sealant material and reducing curing time.This was achieved by regulating the slag doping and sodium silicate solution modulus.The findings demonstrated that increasing slag content and improving the material resulted in an early rise in strength while increasing the modulus of the sodium silicate solution extended the curing time.The slag doping level was fixed at 80 g,and the sodium silicate solution modulus was set at 1.5.To achieve a strength of 3.12 MPa,the water/gel ratio was set at 0.5.The initial setting time was determined to be 33 min,meeting the required field test duration.Secondly,the strength requirements for field implementation were assessed by simulating the action time and force destruction process of the sealing material during blasting using ANSYS/LS-DYNA software.The results indicated that the modified material meets these requirements.Finally,the Shengli Open Pit Coal Mine served as the site for the field test.It was observed that the hole-sealing material’s hydration reaction created a laminated and flocculated gel inside it.This enhanced the density of the modified material.Additionally,the pregelatinized starch,functioning as an organic binder,filled the gaps between the gels,enhancing the cohesion and bonding coefficient of the material.Upon analyzing the post-blasting shooting effect diagram using the Split-Desktop software,it was determined that the utilization of the modified blast hole plugging material resulted in a decrease in the rate of coal fragmentation from 33.2%to 21.1%.This reduction exhibited a minimal error of 1.63%when compared to the field measurement,thereby providing further confirmation of the exceptional plugging capabilities of the modified material.This study significantly contributes to establishing a solid theoretical basis for enhancing the blasting efficiency of open pit mines and,in turn,enhancing their economic advantages.

Safety analysis of Sormeh underground mine to improve sublevel stoping stability

In underground mines,sublevel stoping is used among a variety of different methods for mining an orebody,which creates large underground openings.In this case,the stability of these openings is affected by a number of factors,including the geometrical characteristics of the rock and mining-induced stresses.In this study,a sensitivity analysis was conducted with the numerical,squat pillar,and Mathews stability methods using the Taguchi technique to properly understand the influence of geometric parameters and stress on stope stability according to Sormeh underground mine data.The results show a full factorial analysis is more reliable since stope stability is a complex process.Furthermore,the numerical results indicate that overburden stress has the most impact on stope stability,followed by stope height.However,the results obtained with Mathews and squat pillar methods show that stope height has the greatest impact,followed by overburden stress and span.It appears that these methods overestimate the impact of stope height.Therefore,it is highly recommended that Mathews and squat pillar methods should not be used in high stope that is divided with several sill pillars.Nonetheless,Mathews method cannot accurately predict how the sill pillar impacts the stope stability.In addition,numerical analysis shows that all geometric parameters affect the roof safety factor,whereas the sill pillar has no significant influence on the safety factor of the hanging wall,which is primarily determined by the stope height–span ratio.

Cancelling ore pillars in large-scale coastal gold deposit: A case study in Sanshandao gold mine, China

According to the actual conditions in Sanshandao Gold Mine, the cancelling ore pillars mining method was researched. Firstly, a series of tests for the physical and mechanical characteristics of rock mass were carried out and a quality classification system of rock mass applied in coastal metal deposit was established. Secondly, the reasonable demarcation depth of cancelling ore pillars was simulated using the finite element method, and the simulation results show that the ore pillars can be cancelled below the level of-555 m. Thirdly, a novel layer-backfill mining method of room-pillar alternation was designed to reduce the disturbance and settlement of terrain in mining area. Engineering practice shows that the new mining method enhanced the mining output and relieved rock disturbance. Furthermore, the settlement of the roof strata was small and no disaster occurred. The new mining technology effectively controlled the deformation of the terrain, indicating that the mining of the large-scale gold coastal deposit in Sanshandao Gold Mine was achieved safely, efficiently, and with a low loss rate.

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.

Strategic equilibrium price analysis and numerical simulation of preponderant high-tech metal mineral resources

Based on exploitation compensation value system of preponderant high-tech metal mineral resources and distortion of pricing mechanism, social utility function was constructed to modify decision utility function of developers, and was extended to Stackelberg production decision model of preponderant high-tech metal mineral resources development. Analyzing the influences on market monopoly, output and price decision-making exerted by altruistic preferences, inequity aversion and sequential reciprocity fairness belief equilibrium, game fairness equilibrium which is significant in experimental economics can be obtained and verified by numerical simulation. In process of strategic pricing, method that uses the variation of producer surplus to measure strategic value from psychological preferences was proposed for the first time and technical support to improve exploitation compensation value system of preponderant metal mineral resources was available.

Effects of seepage pressure on the mechanical behaviors and microstructure of sandstone

Surrounding rocks of underground engineering are subjected to long-term seepage pressure,which can deteriorate the mechanical properties and cause serious disasters.In order to understand the impact of seepage pressure on the mechanical property of sandstone,uniaxial compression tests,P-wave velocity measurements,and nuclear magnetic resonance(NMR)tests were conducted on saturated sandstone samples with varied seepage pressures(i.e.0 MPa,3 MPa,4 MPa,5 MPa,6 MPa,7 MPa).The results demonstrate that the mechanical parameters(uniaxial compressive strength,peak strain,elastic modulus,and brittleness index),total energy,elastic strain energy,as well as elastic strain energy ratio,decrease with increasing seepage pressure,while the dissipation energy and dissipation energy ratio increase.Moreover,as seepage pressure increases,the micro-pores gradually transform into meso-pores and macro-pores.This increases the cumulative porosity of sandstone and decreases P-wave velocity.The numerical results indicate that as seepage pressure rises,the number of tensile cracks increases progressively,the angle range of microcracks is basically from 50-120to 80-100,and as a result,the failure mode transforms to the tensile-shear mixed failure mode.Finally,the effects of seepage pressure on mechanical properties were discussed.The results show that decrease in the effective stress and cohesion under the action of seepage pressure could lead to deterioration of strength behaviors of sandstone.

Research on the influences of motion characteristics of jetting projectile charge under water

Shaped charge warhead is important for enhancing the damage performance of underwater weapons.This paper used finite element analysis software and based on JPC water penetration experiments to examine the influence of liner parameters(wall thickness,material),charge aspect ratio,and stand-off distance on the movement characteristics of JPC in water.The findings reveal that the head diameter of the JPC increases and experiences significant erosion after entering the water,the effective length of the JPC in water undergoes two distinct phases:a growth phase and a decrease phase,with the velocity of the JPC decaying exponentially.Increasing the liner thickness,stand-off distance and the charge aspect ratio can improve the erosion resistance and the velocity retention capacity of the JPC.The optimal ranges for liner thickness and stand-off distance are 0.0363D_(k) to 0.0545D_(k)(D_(k) is the charge diameter),the stand-off distance should be within 1.0D_(k).After the charge aspect ratio higher than 1.25,the charge structure exerts minimal influence the movement characteristics of the JPC in water.Material density plays a crucial role in the velocity decay pattern of the JPC during penetration.JPC with higher densities exhibit superior velocity retention capabilities in water,with the velocity decay pattern converging if the densities are similar.Consequently,copper,tantalum and tungsten liners are deemed appropriate for underwater shaped charge warhead.Finally,the results will provide an important reference for the design of underwater shaped charge warhead.

Recent progresses in the development of tubular segmented-in-series solid oxide fuel cells:Experimental and numerical study

Solid oxide fuel cells(SOFCs)have attracted a great deal of interest because they have the highest efficiency without using any noble metal as catalysts among all the fuel cell technologies.However,traditional SOFCs suffer from having a higher volume,current leakage,complex connections,and difficulty in gas sealing.To solve these problems,Rolls-Royce has fabricated a simple design by stacking cells in series on an insulating porous support,resulting in the tubular segmented-in-series solid oxide fuel cells(SIS-SOFCs),which achieved higher output voltage.This work systematically reviews recent advances in the structures,preparation methods,perform-ances,and stability of tubular SIS-SOFCs in experimental and numerical studies.Finally,the challenges and future development of tubular SIS-SOFCs are also discussed.The findings of this work can help guide the direction and inspire innovation of future development in this field.

Evaluating the stability and volumetric flowback rate of proppant packs in hydraulic fractures using the lattice Boltzmann-discrete element coupling method

The stability and mobility of proppant packs in hydraulic fractures during hydrocarbon production are numerically investigated by the lattice Boltzmann-discrete element coupling method(LB-DEM).This study starts with a preliminary proppant settling test,from which a solid volume fraction of 0.575 is calibrated for the proppant pack in the fracture.In the established workflow to investigate proppant flowback,a displacement is applied to the fracture surfaces to compact the generated proppant pack as well as further mimicking proppant embedment under closure stress.When a pressure gradient is applied to drive the fluid-particle flow,a critical aperture-to-diameter ratio of 4 is observed,above which the proppant pack would collapse.The results also show that the volumetric proppant flowback rate increases quadratically with the fracture aperture,while a linear variation between the particle flux and the pressure gradient is exhibited for a fixed fracture aperture.The research outcome contributes towards an improved understanding of proppant flowback in hydraulic fractures,which also supports an optimised proppant size selection for hydraulic fracturing operations.

Thermal safety boundary of lithium-ion battery at different state of charge

Thermal runaway(TR)is a critical issue hindering the large-scale application of lithium-ion batteries(LIBs).Understanding the thermal safety behavior of LIBs at the cell and module level under different state of charges(SOCs)has significant implications for reinforcing the thermal safety design of the lithium-ion battery module.This study first investigates the thermal safety boundary(TSB)correspondence at the cells and modules level under the guidance of a newly proposed concept,safe electric quantity boundary(SEQB).A reasonable thermal runaway propagation(TRP)judgment indicator,peak heat transfer power(PHTP),is proposed to predict whether TRP occurs.Moreover,a validated 3D model is used to quantitatively clarify the TSB at different SOCs from the perspective of PHTP,TR trigger temperature,SOC,and the full cycle life.Besides,three different TRP transfer modes are discovered.The interconversion relationship of three different TRP modes is investigated from the perspective of PHTP.This paper explores the TSB of LIBs under different SOCs at both cell and module levels for the first time,which has great significance in guiding the thermal safety design of battery systems.

Influence of water coupling coefficient on the blasting effect of red sandstone specimens

This study investigates the impact of different water coupling coefficients on the blasting effect of red sandstone.The analysis is based on the theories of detonation wave and elastic wave,focusing on the variation in wall pressure of the blasting holes.Using DDNP explosive as the explosive load,blasting tests were conducted on red sandstone specimens with four different water coupling coefficients:1.20,1.33,1.50,and 2.00.The study examines the morphologies of the rock specimens after blasting under these different water coupling coefficients.Additionally,the fractal dimensions of the surface cracks resulting from the blasting were calculated to provide a quantitative evaluation of the extent of rock damage.CT scanning and 3D reconstruction were performed on the post-blasting specimens to visually depict the extent of damage and fractures within the rock.Additionally,the volume fractal dimension and damage degree of the post-blasting specimens are calculated.The findings are then combined with numerical simulation to facilitate auxiliary analysis.The results demonstrate that an increase in the water coupling coefficient leads to a reduction in the peak pressure on the hole wall and the crushing zone,enabling more of the explosion energy to be utilized for crack propagation following the explosion.The specimens exhibited distinct failure patterns,resulting in corresponding changes in fractal dimensions.The simulated pore wall pressure–time curve validated the derived theoretical results,whereas the stress cloud map and explosion energy-time curve demonstrated the buffering effect of the water medium.As the water coupling coefficient increases,the buffering effect of the water medium becomes increasingly prominent.

Toward understanding the microstructure characteristics,phase selection and magnetic properties of laser additive manufactured Nd-Fe-B permanent magnets

Nd-Fe-B permanent magnets play a crucial role in energy conversion and electronic devices.The essential magnetic properties of Nd-Fe-B magnets,particularly coercivity and remanent magnetization,are significantly infuenced by the phase characteristics and microstructure.In this work,Nd-Fe-B magnets were manufactured using vacuum induction melting(VIM),laser directed energy deposition(LDED)and laser powder bed fusion(LPBF)technologies.Themicrostructure evolution and phase selection of Nd-Fe-B magnets were then clarified in detail.The results indicated that the solidification velocity(V)and cooling rate(R)are key factors in the phase selection.In terms of the VIM-casting Nd-Fe-B magnet,a large volume fraction of theα-Fe soft magnetic phase(39.7 vol.%)and Nd2Fe17Bxmetastable phase(34.7 vol.%)areformed due to the low R(2.3×10-1?C s-1),whereas only a minor fraction of the Nd2Fe14B hard magnetic phase(5.15 vol.%)is presented.For the LDED-processed Nd-Fe-B deposit,although the Nd2Fe14B hard magnetic phase also had a low value(3.4 vol.%)as the values of V(<10-2m s-1)and R(5.06×103?C s-1)increased,part of theα-Fe soft magnetic phase(31.7vol.%)is suppressed,and a higher volume of Nd2Fe17Bxmetastable phases(47.5 vol.%)areformed.As a result,both the VIM-casting and LDED-processed Nd-Fe-B deposits exhibited poor magnetic properties.In contrast,employing the high values of V(>10-2m s-1)and R(1.45×106?C s-1)in the LPBF process resulted in the substantial formation of the Nd2Fe14B hard magnetic phase(55.8 vol.%)directly from the liquid,while theα-Fe soft magnetic phase and Nd2Fe17Bxmetastable phase precipitation are suppressed in the LPBF-processed Nd-Fe-B magnet.Additionally,crystallographic texture analysis reveals that the LPBF-processedNd-Fe-B magnets exhibit isotropic magnetic characteristics.Consequently,the LPBF-processed Nd-Fe-B deposit,exhibiting a coercivity of 656 k A m-1,remanence of 0.79 T and maximum energy product of 71.5 k J m-3,achieved an acceptable magnetic performance,comparable to other additive manufacturing processed Nd-Fe-B magnets from MQP(Nd-lean)Nd-Fe-Bpowder.

A methodology for damage evaluation of underground tunnels subjected to static loading using numerical modeling

We have proposed a methodology to assess the robustness of underground tunnels against potential failure.This involves developing vulnerability functions for various qualities of rock mass and static loading intensities.To account for these variations,we utilized a Monte Carlo Simulation(MCS)technique coupled with the finite difference code FLAC^(3D),to conduct two thousand seven hundred numerical simulations of a horseshoe tunnel located within a rock mass with different geological strength index system(GSIs)and subjected to different states of static loading.To quantify the severity of damage within the rock mass,we selected one stress-based(brittle shear ratio(BSR))and one strain-based failure criterion(plastic damage index(PDI)).Based on these criteria,we then developed fragility curves.Additionally,we used mathematical approximation techniques to produce vulnerability functions that relate the probabilities of various damage states to loading intensities for different quality classes of blocky rock mass.The results indicated that the fragility curves we obtained could accurately depict the evolution of the inner and outer shell damage around the tunnel.Therefore,we have provided engineers with a tool that can predict levels of damages associated with different failure mechanisms based on variations in rock mass quality and in situ stress state.Our method is a numerically developed,multi-variate approach that can aid engineers in making informed decisions about the robustness of underground tunnels.

Surrogate modeling for unsaturated infiltration via the physics and equality-constrained artificial neural networks

Machine learning(ML)provides a new surrogate method for investigating groundwater flow dynamics in unsaturated soils.Traditional pure data-driven methods(e.g.deep neural network,DNN)can provide rapid predictions,but they do require sufficient on-site data for accurate training,and lack interpretability to the physical processes within the data.In this paper,we provide a physics and equalityconstrained artificial neural network(PECANN),to derive unsaturated infiltration solutions with a small amount of initial and boundary data.PECANN takes the physics-informed neural network(PINN)as a foundation,encodes the unsaturated infiltration physical laws(i.e.Richards equation,RE)into the loss function,and uses the augmented Lagrangian method to constrain the learning process of the solutions of RE by adding stronger penalty for the initial and boundary conditions.Four unsaturated infiltration cases are designed to test the training performance of PECANN,i.e.one-dimensional(1D)steady-state unsaturated infiltration,1D transient-state infiltration,two-dimensional(2D)transient-state infiltration,and 1D coupled unsaturated infiltration and deformation.The predicted results of PECANN are compared with the finite difference solutions or analytical solutions.The results indicate that PECANN can accurately capture the variations of pressure head during the unsaturated infiltration,and present higher precision and robustness than DNN and PINN.It is also revealed that PECANN can achieve the same accuracy as the finite difference method with fewer initial and boundary training data.Additionally,we investigate the effect of the hyperparameters of PECANN on solving RE problem.PECANN provides an effective tool for simulating unsaturated infiltration.

Fractal Fractional Order Operators in Computational Techniques for Mathematical Models in Epidemiology

New fractional operators, the COVID-19 model has been studied in this paper. By using different numericaltechniques and the time fractional parameters, the mechanical characteristics of the fractional order model areidentified. The uniqueness and existence have been established. Themodel’sUlam-Hyers stability analysis has beenfound. In order to justify the theoretical results, numerical simulations are carried out for the presented methodin the range of fractional order to show the implications of fractional and fractal orders.We applied very effectivenumerical techniques to obtain the solutions of themodel and simulations. Also, we present conditions of existencefor a solution to the proposed epidemicmodel and to calculate the reproduction number in certain state conditionsof the analyzed dynamic system. COVID-19 fractional order model for the case of Wuhan, China, is offered foranalysis with simulations in order to determine the possible efficacy of Coronavirus disease transmission in theCommunity. For this reason, we employed the COVID-19 fractal fractional derivative model in the example ofWuhan, China, with the given beginning conditions. In conclusion, again the mathematical models with fractionaloperators can facilitate the improvement of decision-making for measures to be taken in the management of anepidemic situation.

Numerical manifold method for thermo-mechanical coupling simulation of fractured rock mass

As a calculation method based on the Galerkin variation,the numerical manifold method(NMM)adopts a double covering system,which can easily deal with discontinuous deformation problems and has a high calculation accuracy.Aiming at the thermo-mechanical(TM)coupling problem of fractured rock masses,this study uses the NMM to simulate the processes of crack initiation and propagation in a rock mass under the influence of temperature field,deduces related system equations,and proposes a penalty function method to deal with boundary conditions.Numerical examples are employed to confirm the effectiveness and high accuracy of this method.By the thermal stress analysis of a thick-walled cylinder(TWC),the simulation of cracking in the TWC under heating and cooling conditions,and the simulation of thermal cracking of the SwedishÄspöPillar Stability Experiment(APSE)rock column,the thermal stress,and TM coupling are obtained.The numerical simulation results are in good agreement with the test data and other numerical results,thus verifying the effectiveness of the NMM in dealing with thermal stress and crack propagation problems of fractured rock masses.

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.

Construction of a Computational Scheme for the Fuzzy HIV/AIDS Epidemic Model with a Nonlinear Saturated Incidence Rate

This work aimed to construct an epidemic model with fuzzy parameters.Since the classical epidemic model doesnot elaborate on the successful interaction of susceptible and infective people,the constructed fuzzy epidemicmodel discusses the more detailed versions of the interactions between infective and susceptible people.Thenext-generation matrix approach is employed to find the reproduction number of a deterministic model.Thesensitivity analysis and local stability analysis of the systemare also provided.For solving the fuzzy epidemic model,a numerical scheme is constructed which consists of three time levels.The numerical scheme has an advantage overthe existing forward Euler scheme for determining the conditions of getting the positive solution.The establishedscheme also has an advantage over existing non-standard finite difference methods in terms of order of accuracy.The stability of the scheme for the considered fuzzy model is also provided.From the plotted results,it can beobserved that susceptible people decay by rising interaction parameters.