Intermittent disturbance mechanical behavior and fractional deterioration mechanical model of rock under complex true triaxial stress paths

Mechanical excavation,blasting,adjacent rockburst and fracture slip that occur during mining excavation impose dynamic loads on the rock mass,leading to further fracture of damaged surrounding rock in three-dimensional high-stress and even causing disasters.Therefore,a novel complex true triaxial static-dynamic combined loading method reflecting underground excavation damage and then frequent intermittent disturbance failure is proposed.True triaxial static compression and intermittent disturbance tests are carried out on monzogabbro.The effects of intermediate principal stress and amplitude on the strength characteristics,deformation characteristics,failure characteristics,and precursors of monzogabbro are analyzed,intermediate principal stress and amplitude increase monzogabbro strength and tensile fracture mechanism.Rapid increases in microseismic parameters during rock loading can be precursors for intermittent rock disturbance.Based on the experimental result,the new damage fractional elements and method with considering crack initiation stress and crack unstable stress as initiation and acceleration condition of intermittent disturbance irreversible deformation are proposed.A novel three-dimensional disturbance fractional deterioration model considering the intermediate principal stress effect and intermittent disturbance damage effect is established,and the model predicted results align well with the experimental results.The sensitivity of stress states and model parameters is further explored,and the intermittent disturbance behaviors at different f are predicted.This study provides valuable theoretical bases for the stability analysis of deep mining engineering under dynamic loads.

Research progress on isotopic fractionation in the process of shale gas/coalbed methane migration

The research progress of isotopic fractionation in the process of shale gas/coalbed methane migration has been reviewed from three aspects: characteristics and influencing factors, mechanism and quantitative characterization model, and geological application. It is found that the isotopic fractionation during the complete production of shale gas/coalbed methane shows a four-stage characteristic of “stable-lighter-heavier-lighter again”, which is related to the complex gas migration modes in the pores of shale/coal. The gas migration mechanisms in shale/coal include seepage, diffusion, and adsorption/desorption. Among them, seepage driven by pressure difference does not induce isotopic fractionation, while diffusion and adsorption/desorption lead to significant isotope fractionation. The existing characterization models of isotopic fractionation include diffusion fractionation model, diffusion-adsorption/desorption coupled model, and multi-scale and multi-mechanism coupled model. Results of model calculations show that the isotopic fractionation during natural gas migration is mainly controlled by pore structure, adsorption capacity, and initial/boundary conditions of the reservoir rock. So far, the isotope fractionation model has been successfully used to evaluate critical parameters, such as gas-in-place content and ratio of adsorbed/free gas in shale/coal etc. Furthermore, it has shown promising application potential in production status identification and decline trend prediction of gas well. Future research should focus on:(1) the co-evolution of carbon and hydrogen isotopes of different components during natural gas migration,(2) the characterization of isotopic fractionation during the whole process of gas generation-expulsion-migration-accumulation-dispersion, and(3) quantitative characterization of isotopic fractionation during natural gas migration in complex pore-fracture systems and its application.

PLANE SURFACE SUDDENLY SET IN MOTION IN A VISCOELASTIC FLUID WITH FRACTIONAL MAXWELL MODEL

The fractional calculus approach in the constitutive relationship model of viscoelastic fluid is introduced.The flow near a wall suddenly set in mo- tion is studied for a non-Newtonian viscoelastic fluid with the fractional Maxwell model.Exact solutions of velocity and stress are obtained by using the discrete in- verse Laplace transform of the sequential fractional derivatives.It is found that the effect of the fractional orders in the constitutive relationship on the flow field is signif- icant.The results show that for small times there are appreciable viscoelastic effects on the shear stress at the plate,for large times the viscoelastic effects become weak.

Unsteady rotating flows of a viscoelastic fluid with the fractional Maxwell model between coaxial cylinders

The fractional calculus is used in the constitutive relationship model of viscoelastic fluid. A generalized Maxwell model with fractional calculus is considered. Based on the flow conditions described, two flow cases are solved and the exact solutions are obtained by using the Weber transform and the Laplace transform for fractional calculus.

A Comparative Study of Fractional Order Models on State of Charge Estimation for Lithium Ion Batteries

State of charge(SOC)estimation for lithium ion batteries plays a critical role in battery management systems for electric vehicles.Battery fractional order models(FOMs)which come from frequency-domain modelling have provided a distinct insight into SOC estimation.In this article,we compare five state-of-the-art FOMs in terms of SOC estimation.To this end,firstly,characterisation tests on lithium ion batteries are conducted,and the experimental results are used to identify FOM parameters.Parameter identification results show that increasing the complexity of FOMs cannot always improve accuracy.The model R(RQ)W shows superior identification accuracy than the other four FOMs.Secondly,the SOC estimation based on a fractional order unscented Kalman filter is conducted to compare model accuracy and computational burden under different profiles,memory lengths,ambient temperatures,cells and voltage/current drifts.The evaluation results reveal that the SOC estimation accuracy does not necessarily positively correlate to the complexity of FOMs.Although more complex models can have better robustness against temperature variation,R(RQ),the simplest FOM,can overall provide satisfactory accuracy.Validation results on different cells demonstrate the generalisation ability of FOMs,and R(RQ)outperforms other models.Moreover,R(RQ)shows better robustness against truncation error and can maintain high accuracy even under the occurrence of current or voltage sensor drift.

Modeling mechanism of a novel fractional grey mode based on matrix analysis

To fully display the modeling mechanism of the novelfractional order grey model （FGM （q,1））, this paper decomposesthe data matrix of the model into the mean generation matrix, theaccumulative generation matrix and the raw data matrix, whichare consistent with the fractional order accumulative grey model（FAGM （1,1））. Following this, this paper decomposes the accumulativedata difference matrix into the accumulative generationmatrix, the q-order reductive accumulative matrix and the rawdata matrix, and then combines the least square method, findingthat the differential order affects the model parameters only byaffecting the formation of differential sequences. This paper thensummarizes matrix decomposition of some special sequences,such as the sequence generated by the strengthening and weakeningoperators, the jumping sequence, and the non-equidistancesequence. Finally, this paper expresses the influences of the rawdata transformation, the accumulation sequence transformation,and the differential matrix transformation on the model parametersas matrices, and takes the non-equidistance sequence as an exampleto show the modeling mechanism.

Simulation of the Fate of Faecal Bacteria in Estuarine and Coastal Waters Based on A Fractionated Sediment Transport Model

A two-dimensional depth-integrated numerical model is refined in this paper to simulate the hydrodynamics, graded sediment transport process and the fate of faecal bacteria in estuarine and coastal waters. The sediment mixture is divided into several fractions according to the grain size. A bed evolution model is adopted to simulate the processes of the bed elevation change and sediment grain size sorting. The faecal bacteria transport equation includes enhanced source and sink terms to represent bacterial kinetic transformation and disappearance or reappearance due to sediment deposition or re-suspension. A novel partition ratio and dynamic decay rates of faecal bacteria are adopted in the numerical model. The model has been applied to the turbid water environment in the Bristol Channel and Severn estuary, UK. The predictions by the present model are compared with field data and those by non-fractionated model.

A Numerical Algorithm Based on Quadratic Finite Element for Two-Dimensional Nonlinear Time Fractional Thermal Diffusion Model

In this article,a high-order scheme,which is formulated by combining the quadratic finite element method in space with a second-order time discrete scheme,is developed for looking for the numerical solution of a two-dimensional nonlinear time fractional thermal diffusion model.The time Caputo fractional derivative is approximated by using the L2-1formula,the first-order derivative and nonlinear term are discretized by some second-order approximation formulas,and the quadratic finite element is used to approximate the spatial direction.The error accuracy O(h3+
t2)is obtained,which is verified by the numerical results.

A Theoretical Investigation of the SARS-CoV-2Model via Fractional Order Epidemiological Model

We propose a theoretical study investigating the spread of the novel coronavirus(COVID-19)reported inWuhan City of China in 2019.We develop a mathematical model based on the novel corona virus’s characteristics and then use fractional calculus to fractionalize it.Various fractional order epidemicmodels have been formulated and analyzed using a number of iterative and numerical approacheswhile the complications arise due to singular kernel.We use the well-known Caputo-Fabrizio operator for the purposes of fictionalization because this operator is based on the non-singular kernel.Moreover,to analyze the existence and uniqueness,we will use the well-known fixed point theory.We also prove that the considered model has positive and bounded solutions.We also draw some numerical simulations to verify the theoretical work via graphical representations.We believe that the proposed epidemic model will be helpful for health officials to take some positive steps to control contagious diseases.

A Fractional Drift Diffusion Model for Organic Semiconductor Devices

Because charge carriers of many organic semiconductors(OSCs)exhibit fractional drift diffusion(Fr-DD)transport properties,the need to develop a Fr-DD model solver becomes more apparent.However,the current research on solving the governing equations of the Fr-DD model is practically nonexistent.In this paper,an iterative solver with high precision is developed to solve both the transient and steady-state Fr-DD model for organic semiconductor devices.The Fr-DD model is composed of two fractionalorder carriers(i.e.,electrons and holes)continuity equations coupled with Poisson’s equation.By treating the current density as constants within each pair of consecutive grid nodes,a linear Caputo’s fractional-order ordinary differential equation(FrODE)can be produced,and its analytic solution gives an approximation to the carrier concentration.The convergence of the solver is guaranteed by implementing a successive over-relaxation(SOR)mechanism on each loop of Gummel’s iteration.Based on our derivations,it can be shown that the Scharfetter–Gummel discretization method is essentially a special case of our scheme.In addition,the consistency and convergence of the two core algorithms are proved,with three numerical examples designed to demonstrate the accuracy and computational performance of this solver.Finally,we validate the Fr-DD model for a steady-state organic field effect transistor(OFET)by fitting the simulated transconductance and output curves to the experimental data.

Equilibrium and kinetic Si isotope fractionation factors and their implications for Si isotope distributions in the Earth's surface environments

Several important equilibrium Si isotope fractionation factors among minerals,organic molecules and the H_4SiO_4 solution are complemented to facilitate the explanation of the distributions of Si isotopes in Earth's surface environments.The results reveal that,in comparison to aqueous H_4SiO_4,heavy Si isotopes will be significantly enriched in secondary silicate minerals.On the contrary,quadra-coordinated organosilicon complexes are enriched in light silicon isotope relative to the solution.The extent of ^(28)Si-enrichment in hyper-coordinated organosilicon complexes was found to be the largest.In addition,the large kinetic isotope effect associated with the polymerization of monosilicic acid and dimer was calculated,and the results support the previous statement that highly ^(28)Sienrichment in the formation of amorphous quartz precursor contributes to the discrepancy between theoretical calculations and field observations.With the equilibrium Si isotope fractionation factors provided here,Si isotope distributions in many of Earth's surface systems can be explained.For example,the change of bulk soil δ^(30)Si can be predicted as a concave pattern with respect to the weathering degree,with the minimum value where allophane completely dissolves and the total amount of sesquioxides and poorly crystalline minerals reaches their maximum.When,under equilibrium conditions,the well-crystallized clays start to precipitate from the pore solutions,the bulk soil δ^(30)Si will increase again and reach a constant value.Similarly,the precipitation of crystalline smectite and the dissolution of poorly crystalline kaolinite may explain the δ^(30)Si variations in the ground water profile.The equilibrium Si isotope fractionations among the quadracoordinated organosilicon complexes and the H_4SiO_4solution may also shed light on the Si isotope distributions in the Si-accumulating plants.

Numerical Analysis and Transformative Predictions of Fractional Order Epidemic Model during COVID-19 Pandemic: A Critical Study from Bangladesh

The COVID-19 pandemic is a curse and a threat to global health, development, the economy, and peaceful society because of its massive transmission and high rates of mutation. More than 220 countries have been affected by COVID-19. The world is now facing a drastic situation because of this ongoing virus. Bangladesh is also dealing with this issue, and due to its dense population, it is particularly vulnerable to the spread of COVID-19. Recently, many non-linear systems have been proposed to solve the SIR (Susceptible, Infected, and Recovered) model for predicting Coronavirus cases. In this paper, we have discussed the fractional order SIR epidemic model of a non-fatal disease in a population of a constant size. Using the Laplace Adomian Decomposition method, we get an approximate solution to the model. To predict the dynamic transmission of COVID-19 in Bangladesh, we provide a numerical argument based on real data. We also conducted a comparative analysis among susceptible, infected, and recovered people. Furthermore, the most sensitive parameters for the Basic Reproduction Number (R0) are graphically presented, and the impact of the compartments on the transmission dynamics of the COVID-19 pandemic is thoroughly investigated.

Analytical wave solutions of an electronically and biologically important model via two efficient schemes

We search for analytical wave solutions of an electronically and biologically important model named as the Fitzhugh–Nagumo model with truncated M-fractional derivative, in which the expafunction and extended sinh-Gordon equation expansion(ESh GEE) schemes are utilized. The solutions obtained include dark, bright, dark-bright, periodic and other kinds of solitons. These analytical wave solutions are gained and verified with the use of Mathematica software. These solutions do not exist in literature. Some of the solutions are demonstrated by 2D, 3D and contour graphs. This model is mostly used in circuit theory, transmission of nerve impulses, and population genetics. Finally, both the schemes are more applicable, reliable and significant to deal with the fractional nonlinear partial differential equations.

Model Identification and Control of Electromagnetic Actuation in Continuous Casting Process With Improved Quality

This paper presents an original theoretical framework to model steel material properties in continuous casting line process. Specific properties arising from non-Newtonian dynamics are herein used to indicate the natural convergence of distributed parameter systems to fractional order transfer function models. Data driven identification from a real continuous casting line is used to identify model of the electromagnetic actuator device to control flow velocity of liquid steel. To ensure product specifications, a fractional order control is designed and validated on the system. A projection of the closed loop performance onto the quality assessment at end production line is also given in this paper.

Letter to the Editor Re “Fractional Modeling and SOC Estimation of Lithium-ion Battery”

A Recent paper by Ma et al.,claims to estimate the state of charge of Lithium-ion batteries with a fractionalorder impedance model including a Warburg and a constant phase element(CPE)with a maximum error of 0.5%[1].The proposed equivalent circuit model from[1]is reproduced in Fig.1.

Regarding Deeper Properties of the Fractional Order Kundu-Eckhaus Equation and Massive Thirring Model

In this paper,the fractional natural decomposition method(FNDM)is employed to find the solution for the Kundu-Eckhaus equation and coupled fractional differential equations describing the massive Thirring model.Themassive Thirring model consists of a system of two nonlinear complex differential equations,and it plays a dynamic role in quantum field theory.The fractional derivative is considered in the Caputo sense,and the projected algorithm is a graceful mixture of Adomian decomposition scheme with natural transform technique.In order to illustrate and validate the efficiency of the future technique,we analyzed projected phenomena in terms of fractional order.Moreover,the behaviour of the obtained solution has been captured for diverse fractional order.The obtained results elucidate that the projected technique is easy to implement and very effective to analyze the behaviour of complex nonlinear differential equations of fractional order arising in the connected areas of science and engineering.

A New Fractional Model for the Falling Body Problem

Recently, a conformable fractional derivative has been proposed to calculate the derivative of non-integer order of time functions. It has been shown that this new fractional derivative definition obeys many advantages over the preceding definitions. For mathematical models in applied sciences and to preserve the dimensionality of the physical quantities, an auxiliary parameter （~r） which has the dimension of seconds should be implemented in the fractional derivative definition. We obtain analytic solutions for the resulting conformable fractional differential equations describing the vertical velocity and the height of the falling body. It is shown that the dimensions of velocity and height are always correct without any restrictions on the auxiliary parameter cr which contradicts with the results in the literature when applying the Caputo definition to the same problem. This may open the door for many future works either to describe the role of such an auxiliary parameter or to derive a more suitable definition for the fractional derivative.

Approximate Solutions for a Class of Fractional-Order Model of HIV Infection via Linear Programming Problem

In this paper, we provide a new approach to solve approximately a system of fractional differential equations (FDEs). We extend this approach for approximately solving a fractional-order differential equation model of HIV infection of CD4+T cells with therapy effect. The fractional derivative in our approach is in the sense of Riemann-Liouville. To solve the problem, we reduce the system of FDE to a discrete optimization problem. By obtaining the optimal solutions of new problem by minimization the total errors, we obtain the approximate solution of the original problem. The numerical solutions obtained from the proposed approach indicate that our approximation is easy to implement and accurate when it is applied to a systems of FDEs.

Implications for fault reactivation and seismicity induced by hydraulic fracturing

Evaluating the physical mechanisms that link hydraulic fracturing(HF) operations to induced earthquakes and the anticipated form of the resulting events is significant in informing subsurface fluid injection operations. Current understanding supports the overriding role of the effective stress magnitude in triggering earthquakes, while the impact of change rate of effective stress has not been systematically addressed. In this work, a modified critical stiffness was brought up to investigate the likelihood, impact,and mitigation of induced seismicity during and after hydraulic fracturing by developing a poroelastic model based on rate-and-state fraction law and linear stability analysis. In the new criterion, the change rate of effective stress was considered a key variable to explore the evolution of this criterion and hence the likelihood of instability slip of fault. A coupled fluid flow-deformation model was used to represent the entire hydraulic fracturing process in COMSOL Multiphysics. The possibility of triggering an earthquake throughout the entire hydraulic fracturing process, from fracturing to cessation, was investigated considering different fault locations, orientations, and positions along the fault. The competition between the effects of the magnitude and change rate of effective stress was notable at each fracturing stage. The effective stress magnitude is a significant controlling factor during fracturing events, with the change rate dominating when fracturing is suddenly started or stopped. Instability dominates when the magnitude of the effective stress increases(constant injection at each fracturing stage) and the change rate of effective stress decreases(the injection process is suddenly stopped). Fracturing with a high injection rate, a fault adjacent to the hydraulic fracturing location and the position of the junction between the reservoir and fault are important to reduce the Coulomb failure stress(CFS) and enhance the critical stiffness as the significant disturbance of stresses at these positions in the coupled process. Therefore,notable attention should be given to the injection rate during fracturing, fault position, and position along faults as important considerations to help reduce the potential for induced seismicity. Our model was verified and confirmed using the case of the Longmaxi Formation in the Sichuan Basin, China, in which the reported microseismic data were correlated with high critical stiffness values. This work supplies new thoughts of the seismic risk associated with HF engineering.

Pore-pressure and stress-coupled creep behavior in deep coal:Insights from real-time NMR analysis

Understanding the variations in microscopic pore-fracture structures(MPFS) during coal creep under pore pressure and stress coupling is crucial for coal mining and effective gas treatment. In this manuscript, a triaxial creep test on deep coal at various pore pressures using a test system that combines in-situ mechanical loading with real-time nuclear magnetic resonance(NMR) detection was conducted.Full-scale quantitative characterization, online real-time detection, and visualization of MPFS during coal creep influenced by pore pressure and stress coupling were performed using NMR and NMR imaging(NMRI) techniques. The results revealed that seepage pores and microfractures(SPM) undergo the most significant changes during coal creep, with creep failure gradually expanding from dense primary pore fractures. Pore pressure presence promotes MPFS development primarily by inhibiting SPM compression and encouraging adsorption pores(AP) to evolve into SPM. Coal enters the accelerated creep stage earlier at lower stress levels, resulting in more pronounced creep deformation. The connection between the micro and macro values was established, demonstrating that increased porosity at different pore pressures leads to a negative exponential decay of the viscosity coefficient. The Newton dashpot in the ideal viscoplastic body and the Burgers model was improved using NMR experimental results, and a creep model that considers pore pressure and stress coupling using variable-order fractional operators was developed. The model’s reasonableness was confirmed using creep experimental data. The damagestate adjustment factors ω and β were identified through a parameter sensitivity analysis to characterize the effect of pore pressure and stress coupling on the creep damage characteristics(size and degree of difficulty) of coal.