During few last years, climate change including global warming which is attributed to human activities and also its long-term adverse effects on the planet’s functions have been identified as the most challenging discussion topics which have arisen many concerns and efforts to find the possible solutions. Since the warmth arising from Earth’s landscapes affects the world’s weather and climate patterns, we decided to study the changes in the Land Surface Temperature (LST) patterns in different seasons through non-linear methods. Here, we particularly want to estimate the non-integer dimension and fractal structure of the land surface temperature. For this study, the (LST) data has been obtained during the daytime by the Moderate Resolution Imaging Spectroradiometer (MODIS) on NASA’s Terra satellite. Depending on what time of the year data has been collected, temperatures change in different ranges. Since equatorial regions remain warm, and Antarctica and Greenland remain cold, and also because altitude affects temperature, we selected Riley County in the U.S. state of Kansas, which does not belong to any of this type locations and we are interested to observe the seasonal changes in temperature in this county. The results of the present study show that the Land Surface Temperature (LST) belongs to the class of fractal process with non-integer dimension.

The study considers building a dynamic model of profit maximization for the car-sharing system and its verification on the example of implementing such a system in Astana (Republic of Kazakhstan). The region, bounded by the administrative boundaries of Astana, was divided into subregions by covering the region with regular hexagons placed side by side. A dataset was built with information on 1168 trips to Astana from January to March 2023. The Kepler visualization service constructed maps of the beginning and end of trips in the region and a map of trips binding to the hexagonal grid cells. Each cell of the grid corresponds to a specific sub-region, for which the quantitative parameters necessary for solving the profit maximization problem in the car-sharing system are calculated. Stations with cars will be placed in the cells of the grid, which will be available to car-sharing service customers. Based on the collected data, a dynamic (four periods per day) and static profit maximization model in the car-sharing system was built. Modeling was carried out based on the built models in the example of Astana. It was established that using a dynamic model of profit maximization in the car-sharing system increases the profit by 3.7%.

Home care centers face both an increase in demand and many variations during the execution of routes, compromising the routes initially planned: robust solutions are not effective enough, it is necessary to move on to resilient approaches. We create a close to reality use case supported by interviews of staff at home health care centers, where caregivers are faced with unexpected events that compromise their initial route. We model, analyze and compare three resilient approaches to deal with these disruptions: a baseline approach without any collaboration, a distributed collaborative approach, and a centralized collaborative approach, where we propose a centralization and sharing of information to improve local decision-making. The latter reduces the number of late arrivals by 11% and the total time of late arrival by 21%, and also halves the number of routes exceeding the end of work time (contrary to the distributed collaborative approach, due to the time wasted reaching colleagues). The use of a device, such as a smartphone application, to centralize and share information thus allows better mutual assistance between caregivers. Moreover, we highlight several possible openings like the coupling of simulation and optimization to propose a more resilient approach.

Our study essentially concerns the dynamic behavior of an SIRS epidemic model in discrete time. Two equilibrium points are obtained; one is disease-free while the other is endemic. We are interested in the endemic fixed point as well as its asymptotic stability. Depending on the parameters which are estimated using the data from US Department of health and SIRS modelling with optimization, two Flip and Transcritical bifurcations appear. We illustrate their diagrams, as well as their bifurcation curves using the method of Carcasses \cite{carcasses1993determination,carcasses1995singularities} for the Flip bifurcation and by an implicit function deduced from such an equation for the Transcritical bifurcation. We use the scanning of the parametric plane to have a global view of the behavior of the model and to highlight the zones of stability of the existing singularities. A superposition of the bifurcation curves with the parametric plane can show the overlap of the curves with the boundaries of the stability domains, which confirms the smooth running of the simulation and its correspondence with the theory, we finish this article with constrained optimal control applied to infection rate and recruitment rate for an SIRS discrete epidemic model. Pontryagin's maximum principle is used to determine these optimal controls. Finally using COVID-19 data in the USA, we obtain results that demonstrate the effectiveness of the proposed control strategy to mitigate the spread of the pandemic.

A stochastic dominance (SD) relation can be defined by two different perspectives: One from the view of distributions, and the other one from the view of expected utilities. In early days, Fishburn investigated SD from the view of distributions and we refer this perspective as Fishburn’s SD. One of his many results was the development of fractional-order SD for continuous distributions. However, discrete fractional-order SD may not be generalized directly since some properties of fractional calculus do not have a discrete counterpart. In this paper, we develop a discrete analogue of fractional-order SD from the view of distributions. We generalize the order of SD by Lizama’s fractional delta operator, show the preservation of SD hierarchy, and formulate the utility classes that are congruent with our SD relations. This work brings a message that some results of discrete SD cannot be generalized directly from continuous SD. We characterize the difference between discrete and continuous fractional-order SD, as well as the way to handle them.

Discrete distribution estimation is a fundamental statistical tool, which is widely used to perform data analysis tasks in various applications involving sensitive personal information. Due to privacy concerns, individuals may not always provide their raw information, which leads to unpredictable biases in the final results of estimated distribution. Local Differential Privacy (LDP) is an advanced technique for privacy protection of discrete distribution estimation. Currently, typical LDP mechanisms provide same protection for all items in the domain, which imposes unnecessary perturbation on less sensitive items and thus degrades the utility of final results. Although, several recent works try to alleviate this problem, the utility can be further improved. In this paper, we propose a novel notion called Item-Oriented Personalized LDP (IPLDP), which independently perturbs different items with different privacy budgets to achieve personalized privacy protection. Furthermore, to satisfy IPLDP, we propose the Item-Oriented Personalized Randomized Response (IPRR) based on the observation that the sensitivity of data shows an inverse relationship with the population size of respective individuals. Theoretical analysis and experimental results demonstrate that our method can provide fine-grained privacy protection and improve data utility simultaneously.

This paper presents three types of tensor Conjugate-Gradient methods for solving large-scale linear discrete ill-posed problems based on the t-product between third-order tensors. An automatical determination strategy of a suitable regularization parameter is proposed for the tensor conjugate gradient (tCG) method. A truncated version and a preprocessed verion of the tCG method are further presented. The discrepancy principle is employed to determine a suitable regularization parameter. Several numerical examples are given to show the effectiveness of the proposed tCG methods in image and video restoration.

Positive invariant set is an important concept of dynamic systems. The purpose of this paper is to study the sufficient and necessary conditions that the set of ellipsoids and the Lorenz cone are positive invariant sets of discrete time systems. By means of nonlinear programming and induced norm, the problem of positive invariance is formulated as an optimization problem, and the equivalent dual form optimization is also presented using the dual optimization method. Our results provide more alternative methods for determining the positive invariance of quadratic form convex sets from the point of view of optimization and dual optimization. The effectiveness of this method is demonstrated by numerical examples.

This paper is concerned with the optimal linear quadratic Gaussian (LQG) control problem for discrete time-varying system with multiple input delays and multiplicative noises. The main contributions are two-fold. Firstly, when the state variables can be observed exactly, we obtain a necessary and suffcient condition for the multiple-delays system in terms of the non-homogeneous relationship between the state and costate, which is the solution to the coupled forward and backward stochastic difference equations. Secondly, when the state variables are partially observed, we derive a suboptimal linear output feedback controller for the discrete-time system based on the obtained results of the optimal LQG control. Numerical examples are shown to illustrate the proposed algorithm.

We offered the circuit of a discrete analog filter (DAF) low frequencies of second-order on switchable capacitors, which differs from the known ones in that three capacitors with one terminal are connected to a common bus, and the frequency of the DAF pole depends not only on the ratio of the capacitances of the frequency-setting capacitors, but also on the ratio of resistances two additional resistors in the feedback. This allows at constant values ​​of the capacitances and a fixed switching frequency of electronic switches, control the pole frequency by changing the resistance ratio of additional resistors. Computer simulation of the developed circuit was performed in the Micro-Cap environment. It has been established that the DAF has a transmission coefficient close to minus one at very low frequencies and a transmission coefficient close to zero at higher frequencies.

We propose a promising model of discrete spacetime based on nonassociative geometry and complex networks. Our approach treats space as a simplicial 3-complex (or complex network), built from "atoms" of spacetime and entangled states forming n-dimensional simplices (n=1,2,3). At large scales, a highly connected network is a coarse, discrete representation of a smooth spacetime. We show that for high temperatures, the network describes disconnected discrete space. At the Planck temperature, the system experiences the phase transition, and for low temperatures, the space becomes a triangulated discrete space. We show that the cosmological constant depends on the universe’s topology. The "foamy" structure, analogous to Wheeler’s "spacetime foam," significantly contributes to the effective cosmological constant, which is determined by the Euler characteristic of the universe.

The paper delves into the critical issue of damage detection within the guy cable of a truss steel mast, a pivotal component in structural integrity. It introduces a model wherein damage manifests as a localized cross-section reduction in a single element of the cable. Employing Discrete Wavelet Transform (DWT), a pioneering methodology, the study scrutinizes the behavior of the affected elements through static structural analyses. Signal decomposition via the Mallat pyramid algorithm facilitates comprehensive examination. Static displacements at the connection point between the cable and the truss mast serve as the measured variables. Through systematic investigation, the paper evaluates the impact of damage size and location, external loading force, and cable tension force on the efficacy of the proposed approach. Utilizing data derived from Finite Element Method (FEM) computations for wavelet analysis the authors substantiate the findings with numerical examples, thus offering valuable insights into damage detection strategies for structural health monitoring and engineering applications.

Polymeric composites such as Poly-ether-ether-ketone (PEEK)/carbon fiber (CF) have been widely utilized due to outstanding performances such as high specific strength and specific modulus. The PEEK/CF components via powder-bed fusion additive manufacturing usually show brittle fracture behaviors induced by their poor interfacial affinity and inner voids. These defects are strongly associated with powder packing quality upon deposition. The particle dynamic model has been widely employed to study the interactions of particle motions. Powder property, bulk material property, and interfacial features of composite powders are key factors in the particle dynamic model. In this work, an efficient and systematic material evaluation is developed for composite powders to investigate their deposition mechanism. The discrete element method is utilized to simulate the dynamic behaviors of PEEK/CF composite powders. The powder properties, bulk material properties, and interfacial features of powders are calibrated and justified by experimental measurement, numerical simulation, and design of experiments. The particle dynamic model can well explain the powder flow behaviors and interactions. It reveals that the addition of short CF particles can assist the flow of PEEK powders and improve the packing quality of the composite powders.

Numerous reservoirs that play a significant role in worldwide petroleum production and reserves contain fractures. Typically, the fractures must form a connected network for a reservoir to be classified as naturally fractured. Characterizing the reservoir with a focus on its fracture network is crucial for modeling and predicting production performance. To simplify the solution, dual continuum modeling techniques are commonly employed. However, to use continuum-scale approaches, properties such as average aperture, permeability, and matrix fracture interaction parameters must be assigned, making it necessary to improve fracture depiction and modeling methods. This study investigates a fractured reservoir with a low matrix permeability and a well-connected fracture network. The focus is on the impact of the hyrachical fracture network on the production performance of gas-based enhanced oil recovery methods. The Discrete Fracture Model (DFN) was utilized to create comprehensive two-dimensional models for three processes: Gas Injection (GI), Water Alternating Gas (WAG), and Foam Assisted Water Alternating Gas (FAWAG). Moreover, dimensionless numbers were employed to establish connections between properties across the entire fracture hierarchy, spanning from minor to major fractures and encompassing the fracture intensity. The results indicate that the FAWAG is more sensitive to fracture types and networks than the WAG and GI processes. Hence, the sensitivity of the individual EOR method to the fracture network requires a respective depth of description of the fracture network. However, other factors, such as reservoir fluids properties and fracture properties, might influence the recovery when the minor fracture networks are excluded. This study has determined that, among the enhanced oil recovery (EOR) techniques examined, the significance of the hyrachical depth of fracture networks diminishes as the ratio of major (primary fracture) aperture to the aperture of medium and minor fractures increases. Additionally, the impact of the assisted-gravity drainage method was greater with increased reservoir height; however, as the intensity ratio increased, the relative importance of medium and minor fracture networks decreased.

The celebrated geodesic congruence equation of Raychaudhuri, together with the resulting singularity theorems of Penrose and Hawking that it enabled, yield a highly general set of conditions under which a spacetime (or, more generically, a pseudo-Riemannian manifold) is expected to become geodesically incomplete. However, the proofs of these theorems traditionally depend upon a collection of assumptions about the continuum spacetime (and, in the physical case, the stress-energy distribution defined over it), including its global structure, its energy conditions, the existence of trapped null surfaces and the various volume/intersection properties of geodesic congruences, that are inherently difficult to translate to the case of discrete spacetime formalisms, such as causal set theory or the Wolfram model. Some, such as the discrete analog of Raychaudhuri's equation for the volumes of geodesic congruences, are subtle to formulate due to intrinsic differences in the behavior of discrete vs. continuous geodesics; for others, such as the definition of a trapped null surface, no appropriate translation is known for the general discrete case due to a lack of a priori coordinate information. It is therefore a non-trivial question to ask whether (and to what extent) there exist equivalently general conditions under which one expects discrete spacetimes to become geodesically incomplete, and how these conditions might differ from those in the continuum. This article builds upon previous work, in which the conformal and covariant Z4 (CCZ4) formulation of the Cauchy problem for the Einstein field equations, with constraint-violation damping, was defined in terms of Wolfram model evolution over discrete (spatial) hypergraphs for the case of vacuum spacetimes, and proceeds to consider a minimal extension to the non-vacuum case by introducing a massive scalar field distribution, defined in either spherical or axial symmetry. Under appropriate assumptions, this scalar field distribution admits a physical interpretation as a collapsing (and, in the axially-symmetric case, uniformly rotating) dust, and we are able to show, through a combination of rigorous mathematical analysis and explicit numerical simulation, that the resulting discrete spacetimes converge asymptotically to either non-rotating Schwarzschild black hole solutions or maximally-rotating (extremal) Kerr black hole solutions, respectively. Although the assumptions used in obtaining these preliminary results are very strong, they nevertheless offer hope that a more general, perhaps ultimately ``Penrose-like'', singularity theorem may be provable in the discrete spacetime case too.

In this study, we focus on an abstract form of summation equations, which can be regarded as a particular class of nonlinear difference equations. We investigate sufficient conditions for the existence of discrete almost automorphic solutions by appealing to the fixed point theory and establish a linkage between the existence of discrete almost automorphic solutions and the existence of bounded solutions for the nonlinear difference equations in the spotlight. Consequentially, we propose a discrete counterpart of the Bohr-Neugebauer theorem for the qualitative theory of difference equations.

Marine vessels operating on the Arctic Sea route are constantly prone to collisions and friction with ice. This study discusses the wear of the hull plate caused by the collision of ice against vessels operating in Arctic Sea routes. The abrasive wear of the hull due to ice impact was numerically assessed based on both the incident behavior of ice particles interacting with the flow around the hull and the wear loss of the hull surface caused by the contact force of ice particles. A multi-phase approach was adopted to account for the behavior of ice particles continuously affected by the fluid force around the hull. The fluid force acting on the ice floe was evaluated using computational fluid dynamics (CFD) and the dynamic motion of the drift ice was evaluated using the discrete element method (DEM). The motion of the floating ice particles was updated in real time by iteratively coupling the fluid force and the motion of the ice floe at each time step of the numerical simulation. The results of the wear simulation models were presented in terms of the shape change of the hull surface due to wear. At first, the wear was evaluated for cases in which only the surface paint of the hull was damaged. Thereafter, a computation model considering the shape change of the hull surface experiencing long-term friction of ice particles was introduced. Finally, the numerical procedures to predict the abrasive wear of the hull surface by ice impact were discussed.

The vibration response of soil is a key property in the field of agricultural soil tillage. Vibration components of tillage machinery are generally used to reduce tillage resistance and improve work efficiency, the pressure variation under low-frequency vibration will affect the fragmentation and dispersion of farmland soil. However, the gradient of pressure variation, frequency domain response, and effective transmission range are unclear. A new method based on DEM (Discrete Element Method) is presented to study the vibration response and pressure transmission under low-frequency vibration. Bench test results showed that peak pressure positively correlates with the vibration frequency and attenuates rapidly at the vibration distance of 100 to 250 mm. The results data were also selected to determine the simulation model parameters. Amplitude, vibration frequency, and soil depth were used as test factors in single-factor simulation tests, and their effects on the peak pressure, frequency domain response, and effective transmission distance were analyzed. The results showed a positive relationship between the peak pressure and the test factors. The peak pressure increases with a maximum gradient of 19.02 kPa/mm at a vibration distance of 50 mm. The amplitude, vibration frequency, and soil depth positively correlate with the dominant frequency amplitude. The main frequency is independent of amplitude and soil depth. At the vibration distance of 250 mm, the dominant frequency is approximately twice the vibration frequency at 7–11 Hz and approximately equal to the vibration frequency at 13–15 Hz. Multiple exponential functions were used to fit the peak pressure attenuation function, obtaining an effective transmission distance range of 347.15 to 550.37 mm for the 5kPa cut-off pressure. For a soil depth of 300 mm, the vertical shear wave diffusion angle is greater than the horizontal shear wave diffusion angle. The study clarifies the vibration response of soil under low-frequency vibration, which helps to design vibration type soil-engaging components of tillage machinery and match vibration parameters for energy saving and resistance reduction in soil tillage.

The paper presents an error model of a measurement chain containing link performing discrete wavelet transform algorithm which is most often the last stage of measurement signal processing. The goal is to determine the uncertainty budget of the input quantities of the wavelet transform algorithm. The error model takes into account parts of analog, analog-to-digital and digital processing, describing the properties of subsequent fragments of the chain using their transmittance and processing functions. The proposed model enables the description of both deterministic and non-deterministic properties of signal errors. The proposed model was validated using an example measurement chain created for this purpose.

In industrial settings, gears play a crucial role by assisting various machinery functions such as speed control, torque manipulation, and altering motion direction.The malfunction or failure of these gear components can have serious repercussions, resulting in production halts and financial losses. As a result, there is an increasing requirement to monitor the state of these components in order to avoid such issues from occurring. To address this need, research efforts have focused on early defect detection in gears in order to reduce the impact of possible failures. This study focused on analyzing vibration and thermal datasets from two extruder machine gearboxes using an autoencoder Long Short-Term Memory (LSTM) model. The major goal is to implement an outlier detection approach to detect and classify defects. The results of this study highlighted the extraordinary performance of the Autoencoder LSTM model, which achieved an impressive accuracy rate of 94.42% in recognizing malfunctioning gearboxes within the extruder machine system. Furthermore, the study used a thorough global metrics evaluation methodology to further test the model’s dependability and efficacy, consequently substantiating the proposed approach’s validity.

The deployment of a manufacturing execution system (MES) holds promising potential in facilitating the accumulation of a substantial amount of inspection data. Low quality levels in discrete manufacturing environments are the result of multi-factor coupling and failure to find quality issues in a timely manner within manufacturing settings may trigger the propagation of defects downstream. Currently, most of the inspection quality methods are direct measurement followed by manual judgment. The integration of deep learning methods provides a feasible way to identify defects in a timely manner, thus improving the acceptance rate of factories. This paper focuses on the design of a data-driven quality prediction and control model around discrete manufacturing characteristics, and use fuzzy theory to evaluate the quality level of production stages. Building Multivariate Long and short-term memory learning hidden quality representations to extract predictions from multi-level inspection data in manufacturing systems. Finally, by validating the data of actual produced water dispensers according to three evaluation indexes, RMSE, MAE, MAPE, the results show that Multivariate Long and short-term memory has better prediction performance.

Numerical simulations of standard situations in the transition region from gas kinetics to fluid dynamics at small Mach numbers indicate a clear dependence of the simulation results on the underlying kinetic model (here: nonlinear and linearized Boltzmann collision operator vs. BGK relaxation model). We develop an improved mathematical framework (trace theory) to explain these differences. In particular we reveal certain deficiencies for the classical BKG system as well as for the standard Navier Stokes approach.

We study the problem of a spherically-symmetric distribution of a perfect relativistic fluid accreting onto a (potentially spinning) black hole within a fully discrete spacetime setting. This problem has previously been studied extensively in the context of continuum spacetimes, beginning with the purely analytic work of Bondi in the spherically-symmetric Newtonian case, Michel in the spherically-symmetric general relativistic case, and Petrich, Shapiro and Teukolsky in the axially-symmetric general relativistic case relevant for spinning black holes. However, the purpose of the present work is to determine the effect of discretization of the underlying spacetime upon the mass/energy and momentum accretion rates, the overall morphology and characteristics of the accretion flow, and the drag force exerted on the black hole in the case of non-zero spin. In order to achieve this, we first develop a novel formulation of the equations of general relativistic hydrodynamics that is more directly amenable to rigorous analysis within a discrete spacetime setting, and we then proceed to implement this formulation into the \textsc{Gravitas} computational general relativity framework. Through a combination of mathematical analysis and explicit numerical simulation in \textsc{Gravitas}, we discover that the mass/energy and momentum accretion rates both decrease monotonically as functions of the underlying spacetime discretization scale, with this effect becoming more pronounced for higher values of the black hole spin parameter, higher fluid temperatures, and stiffer equation of state parameters. We also find that the exerted drag force is highly sensitive to the value of the underlying discretization scale in the case of spinning black hole spacetimes, with certain instabilities becoming significantly more pronounced at certain critical values of the discretization parameter. We discuss some potentially observable consequences of these results, as well as some directions for future theoretical investigation.

The terrain of Changbai Mountain has great influence on the distribution of atmospheric flows and the occurrence and development of precipitation. This study quantitatively analyzes the regional characteristic of topographic perturbations and the relationship between terrain, wind and precipitation in Changbai Mountain region, by using a spectral analysis of the two-dimensional discrete cosine transform. Three domains with relatively heavy summer precipitation are selected as the study region. The results indicate that the overall terrain of Changbai Mountain region exhibits anisotropic characteristics. The terrain spectra of domain B are less than those of domain A and C across the whole wavelength (λ) bands, indicating that the large-scale topographic perturbations of domain B are relatively weak. The largest topographic spectral peak of domain C shows the most pronounced undulation of terrain among the three domains. The dominant wavelengths of terrain height variance for domain A and C, both close to the respective maximum wavelengths, indicate more prominent large-scale topographic perturbations. For domain A, the variation of the precipitation spectra is consistent with that of the wind spectra at the wavelength bands of λ < 390 km, showing a high correlation between wind field and the occurrence of rainfall. The inverse relationship at larger wavelengths indicates that multiple factors contribute to the occurrence of rainfall. For domain B, there is consistency in the fluctuations of terrain spectra, precipitation spectra, and wind spectra at the wavelength bands of λ < 278.3 km, implying that the smaller-scale terrain has an important effect on the occurrence of summer precipitation. For domain C, the variations of terrain spectra, precipitation spectra, and wind spectra are almost consistent across the whole wavelength bands, indicating that the large-scale terrain and minor terrain both play a crucial role in atmospheric uplift and the occurrence and development of summer rainfall.

To assure the opacity, an optimal problem is to retain the occurring events sequences as many as possible. Contrary to such problem, the other optimal problem is to preserve the minimal occurring events sequences. In the paper, based on the choosing cost, an optimal opacity-enforcing model is presented to minimize the sum of discount choosing cost, which is not only to preserve the opacity, but also to retain the secret to the maximum. To solve the model, two scenarios on opacity are considered. For the two scenarios, some algorithms are presented to achieve the optimal solution for the model by using the method of dynamic programming. And then, the solutions produced by the algorithms are proved correctly by theoretical proof. Finally, some illustrations for the algorithms are given.

The establishment of an isolation layer in submarine mining has long presented an enduring challenge. In the context of this research, we conducted a similarity simulation test to preliminarily ascertain the interplay between the thickness and extent of the isolation layer. Subsequently, we introduce an innovative approach that integrates fractal theory and the bonded block model (BBM) to simulate undersea isolation layer mining. The validation of this method hinges on on-site borehole scanning and displacement monitoring, which serve to depict the intricate fractal evolution of fractures and predict the optimal thickness of the isolation layer. Our findings affirm the robustness and validity of this method. Evaluation of the fractal dimensions of fractures reveals that a critical threshold of 1.7 exists to prevent structural failure of the isolation layer, while a limit of 1.5 is necessary to avert significant water ingress. Remarkably, the correlation dimension of the settlement time series closely aligns with the fractal dimension of the fractures, underscoring the feasibility of ensuring the safety of isolation layer mining through real-time settlement monitoring.

The accompanying mechanical behavior during the flow of ore-rock dispersion is an important factor leading to the instability and failure of the orepass in the ore storage section. It is of great significance to accurately understand the instability characteristics of the well wall in the ore section, during the ore drawing process by understanding the flow characteristics and internal mechanical transfer mechanism of ore-rock dispersion. Using the discrete element method to analyze the flow characteristics, contact compactness, stress distribution characteristics and contact force probability distribution of the ore-rock dispersion, it realizes the quantitative characterization of the damage degree of ore-rock flow, and reveals the damage mechanism of mine storage wall. The results show that: (1) The friction between ore-rock particles, between ore-rock and well wall is an important reason for the change of the macroscopic flow form of ore-rock from “one” to “V” type. (2) The compactness and stress concentration of ores are increasing with the increase of storage depth, which is shown that with the increase of storage depth, the greater the probability of well wall instability. (3) The probability distribution of contact force strength between particles decreases exponentially in the whole process of ore drawing, in which the strong chain plays a major role in the stability of the dispersion system. (4) The number of overpressure and the overpressure coefficient can be used to quantitatively characterize the wall damage degree under the action of ore-rock flow. The dynamic side pressure decreases periodically in the form of exponential fluctuations, and the dynamic load formed by the ore-rock flow mainly played a role in the lower area of the ore storage section.

In India, where the transport sector contributes significantly to CO2 emissions, there has been a growing shift towards electric vehicles (EVs), particularly in the three-wheeler segment. Traditional auto-rickshaws, a staple in urban transportation, are often powered by internal combustion engines, leading to pollution and high fuel consumption. Electric three-wheelers, or e-rickshaws, have emerged as a cleaner, more efficient alternative, with zero tailpipe emissions. However, the successful adoption of electric three-wheelers depends on understanding the preferences and challenges faced by auto-rickshaw drivers, who are key to the Intermediate Public Transport (IPT) system. Based on this motivation, the study contributed to the literature by analyzing stated preference survey data from 2,403 respondents who purchased electric tricycles in 2023. We provide insights into the relationship between consumer demographics (such as age, income, gender, education), Vehicle operations, Vehicle financing and Electrification perspective and their preferences for specific electric and fuel tricycle segments. Finally, we estimate consumers' preference for electric tricycles and fuel tricycles. To do this, we use a discrete choice model and estimate the impact factor of Indian consumers in purchasing electric tricycles to obtain market penetration strategies for future electric tricycles.