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.

This paper presents a new numerical method for multiscale modeling of composite materials. The new numerical model, called DECM, consists in a DEM (Discrete Element Method) approach of the Cell Method (CM) and combines the main features of both the DEM and the CM. In particular, it offers the same degree of detail as the CM, on the microscale, and manages the discrete elements individually such as the DEM—allowing finite displacements and rotations—on the macroscale. Moreover, the DECM is able to activate crack propagation until complete detachment and automatically recognizes new contacts. Unlike other DEM approaches for modeling failure mechanisms in continuous media, the DECM does not require prior knowledge of the failure position. Furthermore, the DECM solves the problems in the space domain directly. Therefore, it does not require any dynamic relaxation techniques to obtain the static solution. For the sake of example, the paper shows the results offered by the DECM for axial and shear loading of a composite two-dimensional domain with periodic round inclusions. The paper also offers some insights into how the inclusions modify the stress field in composite continua.

An in-depth study of the failure of granular materials, which is known as a mechanism to generate defects, can reveal the facts about the origin of the imperfections such as cracks in the carbon anodes. The initiation and propagation of the cracks in the carbon anode, especially the horizontal cracks below the stub-holes, reduce the anode efficiency during the electrolysis process. In order to avoid the formation of cracks in the carbon anodes, the failure analysis of coke aggregates can be employed to determine the appropriate recipe and operating conditions. In this paper, it will be shown that a particular failure mode can be responsible for the crack generation in the carbon anodes. The second-order work criterion is employed to analyze the failure of the coke aggregate specimens and the relationships between the second-order work, the kinetic energy, and the instability of the granular material are investigated. In addition, the coke aggregates are modeled by exploiting the discrete element method (DEM) to reveal the micro-mechanical behavior of the dry coke aggregates during the compaction process. The optimal number of particles required for the failure analysis in the DEM simulations is determined. The effects of the confining pressure and the strain rate as two important compaction process parameters on the failure are studied. The results reveal that increasing the confining pressure enhances the probability of the diffusing mode of the failure in the specimen. On the other hand, the increase of strain rate augments the chance of the strain localization mode of the failure in the specimen.

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.

Aiming at the problems of high cost of manual pulling, low reliability of existing pulling devices, high breaking rate and high leakage rate in the process of cotton stalk reuse after leaving the field in Xinjiang cotton area, a loosening soil and root-cutting cotton stalk pulling the setting machine has been researched and designed, established the root-cutting force model, calculated and optimized the key parameters of the V-shaped root-cutting knife, and determined the range of slide cutting angle, cutting edge angle and soil entry angle. The shoveling process simulation of V-shaped root-cutting knife and root-soil complex was constructed, and the working mechanism of V-shaped root-cutting knife was clarified. In order to verify the reliability and operation performance of the V-shaped root-cutting knife, the slide cutting angle, cutting edge angle and soil entry angle were used as the test factors, the response surface test of three-factors and three-levels was carried out with root-breaking force and mean value of the cutting resistance as the test index. The test results were analyzed by variance analysis, the significant of influencing factors on the root-breaking force are in a descending order as: slide cutting angle, cutting edge angle, soil entry angle; the degree of influence on the mean value of the cutting resistance are as follows: slide cutting angle, soil entry angle, cutting edge angle. In order to make the V-shaped root-cutting knife achieve the optimal working state, the parameters of the test index were optimized, and the optimal design parameters of the V-shaped root-cutting knife were set as follows: the slide cutting angle was 48.3°, cutting edge angle was 43.4°, and soil entry angle was 26.2°. The field uprooting test showed that the average pass rate of root breakage was 94.8% and the average pull-out rate of the cotton stalk was 93.2%. This study provides theoretical guidance for the development of root breaking mechanism of cotton straw harvester.

The paper presents the problem of damage detection in steel girders. Static displacements at the selected point of the structure play the role of measured variables. Structural response signal decomposition is performed according to the Mallat pyramid algorithm, which is used to perform the discrete wavelet transform (DWT). This procedure allows us to quite well determine the location of structural damage. The geometry and the placement of any defective part of the structure may have a random character. It can be assumed that the random processes occurring in the broadly understood structure mechanics are Gaussian in nature. The first four probabilistic moments are estimated using three approaches independent: semi-analytical (SAM), perturbative (SPT), and Monte-Carlo simulations (MCS). The semi-analytical random approach seems to be the most optimal due to the necessary computation time. The incorporation of the mathematical stochastic apparatus into the classical (deterministic) analysis of the statics of the structure makes it possible to estimate the reliability measures of the analyzed girder.

Computation of Discrete Fourier Transform (DFT) is a challenging task. Especially, on computational machines/embedded systems where resources are limited. The importance of Fourier Transform (FT) cannot be denied in the field of signal processing. This paper proposes a technique that can compute Discrete Fourier Transforms for a matrice or vector with the help of matrix multiplication. Moreover, this paper discusses the trivial methods used for computation of DFT along with methods based on matrix multiplication used to compute discrete Fourier Transform in addition to the shortcomings. The proposed method can help in the calculation of a Discrete Fourier Transformation matrix by truncation of values from the proposed generic method which can help in computing DFT of varying lengths of vectors. On legacy computing machines and programming environments, having support for matrix multiplication, the proposed methodology can be implemented.

With the rapid development of earth observation, satellite navigation, mobile communication and other technologies, the order of magnitude of the spatial data we acquire and accumulate is increasing, and higher requirements are put forward for the application and storage of spatial data. Under this circumstance, a new form of spatial data organization emerged-the global discrete grid. This form of data management can be used for the efficient storage and application of large-scale global spatial data, which is a digital multi-resolution the geo-reference model that helps to establish a new model of data association and fusion. It is expected to make up for the shortcomings in the organization, processing and application of current spatial data. There are different types of grid system according to the grid division form, including global discrete grids with equal latitude and longitude, global discrete grids with variable latitude and longitude, and global discrete grids based on regular polyhedrons. However, there is no accuracy evaluation index system for remote sensing images expressed on the global discrete grid to solve this problem. This paper is dedicated to finding a suitable way to express remote sensing data on discrete grids, and establishing a suitable accuracy evaluation system for modeling remote sensing data based on hexagonal grids to evaluate modeling accuracy. The results show that this accuracy evaluation method can evaluate and analyze remote sensing data based on hexagonal grids from multiple levels, and the comprehensive similarity coefficient of the images before and after conversion is greater than 98%, which further proves that the availability hexagonal grid-based remote sensing data of remote sensing images. And among the three sampling methods, the image obtained by the nearest interpolation sampling method has the highest correlation with the original image.

An approach to estimating the objective function value of minimization maximum lateness problem is proposed. It is shown how to use transformed instances to define a new continuous objective function. After that, using this new objective function, the approach itself is formulated. We calculate the objective function value for some polynomially solvable transformed instances and use them as interpolation nodes to estimate the objective function of the initial instance. What is more, two new polynomial cases, that are easy to use in the approach, are proposed. In the end of the paper numeric experiments are described and their results are provided.

In this work, numerical simulations for the absorption and scattering efficiencies of spheroid core–shell nanoparticles (CSNs) were conducted and studied using the discrete-dipole approximation method. The characteristics of surface plasmon resonances (SPR) depend upon shell thickness, the compositions of the core and shell materials, and the aspect ratio of the constructed CSNs. We used different core@shell compositions, specifically Au@SiO₂, Ag@SiO₂, Au@TiO₂, Ag@TiO₂, Au@Ag, and Ag@Au, for extinction spectra analysis. We also investigated coupled resonance mode wavelengths by adjusting the composition’s layer thickness and aspect ratio. In this study, we show that the extinction efficiency of the Ag@TiO₂ core–shell nanoparticles (CSNPs) was higher than that of the others, and we examined the impact of TiO₂ shell thickness and Ag core radius on SPR peak positions. From the extinction spectra we found that the Ag@TiO₂ nanoparticle had better refractive index sensitivity and figure of merit when the aspect ratio was set to 0.3. All of the experimental results proved that the tunability of these plasmonic resonances was highly dependent on the material used, the layer thickness, and the aspect ratio of the core@shell CSNPs.

In this paper we introduce a computational model for the simulation of hypervelocity impact (HVI) phenomena which is based on the Discrete Element Method (DEM). Our paper constitutes the first application of DEM to the modeling and simulating of impact events for velocities beyond 5 kms⁻¹. We present here the results of a systematic numerical study on HVI of solids. For modeling the solids, we use discrete spherical particles that interact with each other via potentials. In our numerical investigations we are particularly interested in the dynamics of material fragmentation upon impact. We model a typical HVI experiment configuration where a sphere strikes a thin plate and investigate the properties of the resulting debris cloud. We provide a quantitative computational analysis of the resulting debris cloud caused by impact and a comprehensive parameter study by varying key parameters of our model. We compare our findings from the simulations with recent HVI experiments performed at our institute. Our findings are that the DEM method leads to very stable, energy–conserving simulations of HVI scenarios that map the experimental setup where a sphere strikes a thin plate at hypervelocity speed. Our chosen interaction model works particularly well in the velocity range where the local stresses caused by impact shock waves markedly exceed the ultimate material strength.

In today’s world face detection is the most important task. Due to the chromosomes disorder sometimes a human face suffers from different abnormalities. For example, one eye is bigger than the other, cliff face, different chin-length, variation of nose length, length or width of lips are different, etc. For computer vision currently this is a challenging task to detect normal and abnormal face and facial parts from an input image. In this research paper a method is proposed that can detect normal or abnormal faces from a frontal input image. This method used Fast Fourier Transformation (FFT) and Discrete Cosine Transformation of frequency domain and spatial domain analysis to detect those faces.

Discrete Double Hilbert Exponential sums along polynomials. This is preprint version (not peer reviewed) And we now check the proof of (ii) of Lemma 3.2 whether it is really true.

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.

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.

High Efficiency Video Coding (HEVC) uses an 8-point filter and a 7-point filter, which are based on the discrete cosine transform (DCT), for the 1/2-pixel and 1/4-pixel interpolations, respectively. In this paper, discrete sine transform (DST)-based interpolation filters (IF) are proposed. The first proposed DST-based IFs (DST-IFs) use 8-point and 7-point filters for the 1/2-pixel and 1/4-pixel interpolations, respectively. The final proposed DST-IFs use 12-point and 11-point filters for the 1/2-pixel and 1/4-pixel interpolations, respectively. These DST-IF methods are proposed to improve the motion-compensated prediction in HEVC. The 8-point and 7-point DST-IF methods showed average BD-rate reductions of 0.7% and 0.3% in the random access (RA) and low delay B (LDB) configurations, respectively. The 12-point and 11-point DST-IF methods showed average BD-rate reductions of 1.4% and 1.2% in the RA and LDB configurations for the Luma component, respectively.

In this paper we present two algorithms for computing discrete logarithms. We present some properties and show how we break the Diffie-Hellman protocol, which is considered to be very secure.

Several factors influence traffic congestion and overall traffic dynamics. Simulation modeling has been utilized to understand the traffic performance parameters during traffic congestions. This paper focuses on driver behavior of route selection by differentiating three distinguishable decisions, which are shortest distance routing, shortest time routing and less crowded road routing. This research generated 864 different scenarios to capture various traffic dynamics under collective driving behavior of route selection. Factors such as vehicle arrival rate, behaviors at system boundary and traffic light phasing were considered. The simulation results revealed that shortest time routing scenario offered the best solution considering all forms of interactions among the factors. Overall, this routing behavior reduces traffic wait time and total time (by 69.5% and 65.72%) compared to shortest distance routing.

Four Fourier transforms are usually defined, the Integral Fourier transform, the Discrete-Time Fourier transform (DTFT), the Discrete Fourier transform (DFT) and the Integral Fourier transform for periodic functions. However, starting from their definitions, we show that all four Fourier transforms can be reduced to actually only one Fourier transform, the Fourier transform in the distributional sense.

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.

The current theories of consciousness do not address the two major problems:What is the solution to the hard problem of consciousness or What is the experience of self-existence? and How is conscious perception computed? We address these questions using our hypothesized theory:The Self Simulation theory and propose a model of computation of consciousness.The components of the model are an interconnected network of auxiliary processors or Auxiliary Neural Networks(AN) which receives input from the external sensory modalities,computes a representation of the information and relays it to the NCC(Neural Correlates of Consciousness), the NCC – a processing unit which can be neuro anatomically distributed in location but computationally connected to each other and to the Auxiliary Networks,and a working memory which the NCC references for further computation of information.Computation of conscious perception is achieved via the decoding of the content of the memory by the NCC.The self simulation theory accounts not only for the various phenomenal experiences but also experimentation paradigms conducted in the past.

Drilling and blasting is the conventional method used for rock fragmentation in open pit mining. Blast-induced damage can reduce the level of stability of benches and pit slopes. To develop an optimal blast design, an adequate knowledge of the rock properties and in-situ fractures is needed. Fractures are generally the paths of least resistance for explosive energy and can affect the intensity of blast-induced damage. Discrete Fracture Networks (DFN) are 3D representations of joint systems used for estimating the distribution of in-situ fractures in a rock mass. The combined finite/discrete element method (FDEM) can be used to simulate the complex rock breakage process during a blast. The objective of this paper is to develop a methodology for assessing the influence of in-situ joints on the post-blast fracturing and the associated wall damage in 2D bench blasts scenarios. First, a simple one-blasthole scenario is analyzed with the FDEM software Irazu 2D and calibrated based on a laboratory-scale blasting experiment available from previous literature. Secondly, more complex scenarios consisting of one-blasthole models at the bench-scale were simulated. A bench blast without DFN (base case) and one with DFN, were numerically simulated. The model with DFN demonstrated that the growth path and intensity of blast-induced fractures were governed by pre-existing fractures which led to smaller wall damage area. The damage intensity for the base case scenario is about 42% higher than for the blast model with DFN included, which highlights the significance of in-situ fractures in the resulting blast damage intensity. The methodology for developing the DFN-included blasting simulation provides a more realistic modeling process for blast-induced wall damage assessment. This results in a better characterization of the blast damage zone and can lead to improved slope stability analyses.

The theory of the continuous two-dimensional (2D) Fourier Transform in polar coordinates has been recently developed but no discrete counterpart exists to date. In the first part of this two-paper series, we proposed and evaluated the theory of the 2D discrete Fourier Transform (DFT) in polar coordinates. The theory of the actual manipulated quantities was shown, including the standard set of shift, modulation, multiplication, and convolution rules. In this second part of the series, we address the computational aspects of the 2D DFT in polar coordinates. Specifically, we demonstrate how the decomposition of the 2D DFT as a DFT, Discrete Hankel Transform (DHT) and inverse DFT sequence can be exploited for efficient code. We also demonstrate how the proposed 2D DFT can be used to approximate the continuous forward and inverse Fourier transform in polar coordinates in the same manner that the 1D DFT can be used to approximate its continuous counterpart.

The theory of the continuous two-dimensional (2D) Fourier Transform in polar coordinates has been recently developed but no discrete counterpart exists to date. In this paper, we propose and evaluate the theory of the 2D discrete Fourier Transform (DFT) in polar coordinates. This discrete theory is shown to arise from discretization schemes that have been previously employed with the 1D DFT and the discrete Hankel Transform (DHT). The proposed transform possesses orthogonality properties, which leads to invertibility of the transform. In the first part of this two-part paper, the theory of the actual manipulated quantities is shown, including the standard set of shift, modulation, multiplication, and convolution rules. Parseval and modified Parseval relationships are shown, depending on which choice of kernel is used. Similar to its continuous counterpart, the 2D DFT in polar coordinates is shown to consist of a 1D DFT, DHT and 1D inverse DFT.

In this research article we solve the problem of synchronization and anti-synchronization of chaotic systems described by discrete and time-delayed variable fractional order differential equations. To guarantee the synchronization and anti-synchronization of these systems, we use the well-known PID control theory and the Lyapunov-Krasovskii stability theory for discrete systems of variable fractional order.We illustrate the results obtained through simulation with examples, in which it can be seen that our results are satisfactory, thus achieving synchronization and anti-synchronization of chaotic systems of variable fractional order with discrete time delay.

In my paper, a divisibility rule named the Queen function is de- scribed, which has been generalized for every integer and encompasses all pre- viously established divisibility rules. Subsequently, the paper discusses other areas where this function is instrumental. These include simplifying ratios, generating prime numbers, and more. Towards the end of the paper, a set of hypotheses is also presented. Among these, the most significant is the novel approach involving the application of the Queen function for prime number generation.

In this work, the hypothesis that the universe is made up of 4-dimensional spheres of space, whose diameter is the Planck length, allows us to calculate most of the constants used in current physical theories. Calculated constants include: elementary charge, fine structure constant, electron mass, Planck constant, gravitational constant, electrical constant, Boltzmann constant, mass and charge of up and down quarks, muon mass, and Higgs boson. All of these depend on the speed of light and the Planck length, which shows that all the constants are related and not due to chance.

By applying a running average (with a window-size= d), we could transform Discrete data to broad-range, Continuous values. When we have more than 2 columns and one of them is containing data about the tags of classification (Class Column), we could compare and sort the features (Non-class Columns) based on the R2 coefficient of the regression for running averages. The parameters tuning could help us to select the best features (the non-class columns which have the best correlation with the Class Column). “Window size” and “Ordering” could be tuned to achieve the goal. this optimization problem is hard and we need an Algorithm (or Heuristics) for simplifying this tuning. We demonstrate a novel heuristics, Called Simulated Distillation (SimulaD), which could help us to gain a somehow good results with this optimization problem.

Selenium and iodine are essential micronutrients for humans. They are often deficient in food supply due to low phytoavailable concentrations in soil. Agronomic biofortification of food crops is one approach to overcome micronutrient malnutrition. This study focused on German consumers’ willingness to purchase selenium- and/or iodine-biofortified apples. For this purpose, an online survey was carried out. In this context, consumers were asked to choose their most preferred apple product from a set card of product alternatives (Discrete Choice Experiment). The multinomial logit model results demonstrated that German consumers’ have a preference in particular for iodine-biofortified apples. Furthermore, apple choice was mainly influenced by price, health claims, and plastic-free packaging material. Viewed individually, selenium did not exert an effect on product choice whereas positive interactions between both micronutrients exist.

The basis of the discrete-impulse energy supply (DIES) concept is the efficient use of supplied energy. The references describe in detail the general principles of DIES, examine the energy and thermodynamic aspects and the main mechanisms of intensification that can be initiated on the basis of this principle. DIES mechanisms conveniently can be divided into hard and soft ones. The former should be used to stimulate hydromechanical processes, and the latter to accelerate the processes of phase heat and mass transfer, or for the purpose of intensive mixing of multicomponent media. The authors have studied the possibility of using DIES to intensify the hydromechanical processes, in particular emulsification of milk fat (homogenization of milk, preparation of spreads), processing of cream cheese masses. Objects of research were whole non-homogenized milk, fat emulsions, cream cheese mass. In order to evaluate the efficiency of milk homogenization the homogenization coefficient change was studied, which was determined by centrifugation method as the most affordable and accurate one. Emulsions were evaluated according to the degree of destabilization, resistance and dispersion of the fat phase. The rheological characteristics of cheese masses were evaluated by the effective viscosity change.

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 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%.

We formulate the statistics of the discrete multicomponent fragmentation event using a methodology borrowed from statistical mechanics. We generate the ensemble of all feasible distributions that can be formed when a single integer multicomponent mass is broken into fixed number of fragments and calculate the combinatorial multiplicity of all distributions in the set. We define random fragmentation by the condition that the probability of distribution be proportional to its multiplicity and obtain the partition function and the mean distribution in closed form. We then introduce a functional that biases the probability of distribution to produce in a systematic manner fragment distributions that deviate to any arbitrary degree from the random case. We corroborate the results of the theory by Monte Carlo simulation and demonstrate examples in which components in sieve cuts of the fragment distribution undergo preferential mixing or segregation relative to the parent particle.

This paper deals with the service performance analysis and improvement using discrete event simulation has been used. The simulation of the heath care has been done by arena master development 14-version software. The performance measurement for this study are patients output, service rate, service efficiency and it is directly related to waiting time of patients in each service station, work in progress, resource utilization.Simulation model was building for Bahir Dar clinic and then, prepared the proposed model for the system. Based on the simulation model run result, the output of the existing healthcare service system is low due to presence of bottlenecks on the service system. Moreover, the station with the largest queue and high resource utilization are identified as a bottleneck. The bottlenecks, which have identified are reduced by using reassigning the existing resources and add new resources and merging the similar services, which has under low resource utilization (nurses). Finally, the researchers have proposed a developed model from different scenarios. Moreover, the best scenario is developed by combining scenario 2 and 3. And then, service efficiency of the healthcare has increased by 9.86 percent, the work in progress (WIP) are reduced by 3 patients from the system and the service capacity of the system is increased 34 to 40 patients per day due to the reduction of bottleneck stations.

Cheating prevention in online exams is often hard and costly to tackle with proctoring, and it even sometimes involves privacy issues, especially in social distancing due to the pandemic of COVID-19. Here we propose a low-cost and privacy-preserving anti-cheating scheme by programmatically minimizing the cheating gain. A novel anti-cheating scheme we developed theoretically ensures that the cheating gain of all students can be controlled below a desired level aided by the prior knowledge of students’ abilities and a proper assignment of question sequences. Furthermore, a heuristic greedy algorithm we developed can refine an assignment of questions from a cyclic pool of question sequences to efficiently reduce the cheating gain. Compared to the integer linear programming and min-max matching methods in a small-scale simulation, our heuristic algorithm provides results close to the optimal solutions offered by the two standard discrete optimization methods. Hence, our anti-cheating approach could potentially be a cost-effective solution to the well-known cheating problem even without proctoring.

A common reactor type in the chemical and process industry is the fixed-bed reactor. Accurate modeling can be achieved with particle-resolved Computational Fluid Dynamic (CFD) simulations. However, the underlying bed morphology plays a paramount role. Synthetic bed-generation methods are much more flexible and faster than image-based approaches. In this study, we look critically at the two different bed generation methods: Discrete Element Method (DEM) (in the commercial software STAR-CCM+) and the rigid-body model (in the open-source software Blender). The two approaches are compared in terms of synthetically generated beds with experimental data of overall and radial porosity, particle orientation, as well as radial velocities. Both models show accurate agreement for the porosity. However, only Blender shows similar particle orientation than the experimental results. The main drawback of the DEM is the long calculation time and the shape approximation with composite particles.

Background: Gait analysis systems are crucial for rehabilitation but are often time-consuming and less accessible in low and middle-income countries like Brazil. GAITWell^® was developed to address these challenges with its modular and portable design for automated gait data extraction and analysis. This study evaluates its methodological properties. Methods: GAITWell^® uses discrete binary sensors on interconnected plates to collect gait data, analyzed via the DBSCAN algorithm. Each plate measures 44 cm x 37 cm and features an 11 x 7 discrete sensor array with 4 cm spacing. Test-retest reliability was assessed using the intraclass correlation coefficient (ICC), standard error of the mean (SEM), and Bland-Altman plots. Concurrent validity was evaluated by comparing the GAITWell^® measurements with those obtained from the Qualisys Pro-Reflex through correlations. Results: Participated 38 healthy young adults, including 18 females (47.4%). The sample had an average age of 33.2 years (SD 13.0) and an average BMI of 26.6 kg/m² (SD 5.4). Correlations between GAITWell^® and Qualisys ranged from moderate to very high for most gait variables, except stride length (r = .360, p < .05). ICC_2,1 showed moderate to good agreement for most variables, except for double support times, cadence, and base of support. Preliminary analysis suggests that increasing the sensor resolution to 1.2 cm x 1.2 cm could potentially decrease the measurement error from 2 cm to 0.6 cm, representing a theoretical reduction of 70%. Conclusions: While low resolution affected GAITWell^® reliability, it remains a promising tool. Future research should focus on enhancing sensor accuracy to improve reliability.

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.

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.

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.

This addendum contains clarifications and a sharpening of some of the results on the VIT transform framework developed in [1]. The focus is on the right-coefficient and left- coefficient forms of the transform, the extraction of a first-order term from a left polynomial fraction, and the application to linear time-varying systems.

Both drag and lift forces impact an inclined plane when it is dragged through a granular bed. In this paper, the following results have been obtained: the drag and lift forces grow with the velocity of motion; when the immersion depth is constant, the inclination angle has no effect on drag force, however, the lift force increases linearly with this inclination angle; the ratio of drag and lift forces is exactly equal to the tangent value of the inclined angle. In order to describe this physical process macroscopically, a continuum wedge model based on the Coulomb model is established to predict drag and lift forces. Particularly，the dynamic friction angle in the assumed shear band is predicted as a function of both inclined angle and moving velocity.

Given the Income enabling nature of education as stipulated by human capital theory, it can be postulated that “ceteris paribus”, households with formally educated heads experience less food insecurity and monetary poverty than those with uneducated heads. We test this claim in the case of Burkina Faso, using the 2014 National Survey on Households Living Conditions, along with semi-parametric modeling techniques. In its design the study uses households “willingness and ability” to spend annually on food consumption a per-capita amount above the food poverty line of 102,040 CFA Franc to characterize “household food security”, and households “unwillingness or inability” to spend above the overall poverty line of 153,530 CFA Franc to characterize “monetary poverty”. In addition, the study relies not only on single equation multivariate probit and logit specifications, but also on both fully parametric and semi-parametric bivariate probit representations of food insecurity and monetary poverty. The results show that relaxing the linearity and independence assumptions through joint semi-parametric bivariate modeling captures better the true effects of heads of households’ educational attainment on households' food insecurity and monetary poverty. In fact, compared to households headed by someone with no education, those headed by someone with a primary, secondary or higher education are respectively 19.8% , 49.7% and 118.9% less likely to experience food insecurity, and respectively 40.1%, 77.0% and 172.3% less likely to experience monetary poverty in Burkina Faso. In addition, the experience of food insecurity and monetary poverty are highly correlated at 92.7%, suggesting that educational policies that alleviate poverty in Burkina Faso should also impact positively food security in the country.

The inability to delineate a unified physical ontology that accounts simultaneously for the laws of special relativity and the results of quantum experiments has been a defining problem in physics for more than 100 years. This analysis addresses the problem by positing a novel mixed ontology composed of a “discrete” 4D spacetime and a physical, ultra-high dimensional (3 x N) “Planck Space.” Together, Planck Space and the three spatial dimensions of 4D spacetime form a tightly integrated ((3 x N) + 3) hyperspace (the “Dual Ontology”).Critically, the Dual Ontology’s structure replaces 1) the continuous, differentiable manifold of 4D spacetime with a discrete, background-independent 4D spacetime that is intrinsically curved and 2) mathematical spaces with a physical (3 x N) Planck Space. Moreover, the Dual Ontology is structurally and dynamically predicated on the one-to-one mapping and identity between the discrete spatial units that simultaneously form 4D spacetime and Planck Space. The one-to-one mapping and identity of the discrete spatial units support an integrated quantum dynamics based upon the dynamic evolution of single and N-body quantum states in 4D spacetime in full compliance with the laws of special relativity and the instantaneous collapse of all quantum states in an ontic Planck Space, where special and general relativity and more generally, 4D spacetime’s laws of physics, do not apply.

The paper presents two procedures for creating 3D discrete models of apple using the finite element method (FEM). The study focused on a new Polish apple variety, Chopin at development, ripening, and senescence maturity stages. The analysis of the results allowed for the determination that, depending on the desired outcomes, the parenchymal tissue can be considered as a material with either elastic or elastoplastic properties. Model validation was performed using displacement-force relationships and changes in contact area as a function of loading force obtained from compression tests. Elastoplastic models reflect empirical studies much better. The fit of the elastic models, based on force vs. displacement curves, ranged from 79% to 83%, while the fit of the elastoplastic models to empirical data ranged from 91% to 96%. Elastic models, on the other hand, are well-suited to the modified force vs. displacement curves, which reflect the elastic behavior of apples under compression test. Contour plots of pressures for elastoplastic models align much better with the test results than those for elastic models. 4. Similarly, in the case of contact area values, the elastoplastic models outperformed the elastic models, achieving a fit between 92% and 96%, while the elastic models reached only 82% to 83% compliance.

We consider a system of two coupled parallel queues, queue 1 and queue 2, with infinite waiting rooms. The time setting is discrete. In either queue, the service of a customer requires exactly one discrete time slot. Arrivals of new customers occur independently from slot to slot, but the numbers of arrivals of both types within a slot may be mutually interdependent. Their joint probability generating function (pgf) is indicated as A(z1, z2) and characterizes the whole model. It is well-known that, in general, determining the steady-state joint probability mass function (pmf) u(m, n), m, n ≥ 0 or the corresponding joint pgf U(z1, z2) of the system contents (i.e., numbers of customers present) in both queues is a formidable task. Only in a number of isolated cases, for very specific choices of the arrival pgf A(z1, z2), explicit results are known in the literature. In this paper, we identify a multiparameter, generic class of arrival pgfs A(z1, z2), for which we can explicitly determine the system-content pgf U(z1, z2). We find that, for arrival pgfs of this class, U(z1, z2) has a denominator which is a product, say r1(z1)r2(z2) of two univariate functions. This property allows a straightforward inversion of U(z1, z2), resulting in a pmf u(m, n) which can be expressed as a finite linear combination of bivariate geometric terms. We also observe that our generic model encompasses most of the previously known results as special cases.

The dynamic evolution and variation of electrical loads is now a priority for optimum management and, above all, forecasting. Indeed, these dynamic load variations require computer tools able to implement optimal load forecasting models. Scientific research into automated models for forecasting electrical loads therefore represents a challenge for scientific researchers. Several research studies have therefore been carried out. These include machine learning approaches such as: LSTM (Long Short-Term Memory), Support Vector Machine (SVM) and others. This article proposes a comparison study between a hybrid model based on the wavelet transform coupled with the adaboost method, and LSTM, for forecasting electricity consumption loads. Python 3.10 was used for the studies. The first results obtained consisted in the study of electricity load forecasts over four (04) data sets of the year (January, April, August and December). These results are very satisfactory, showing a good correlation between the real results of the electrical load consumption data, and the results of the implemented model. The electrical load consumption data are collected from a renewable energy production source: photovoltaic solar energy. The model test results (R2) obtained range from 0.919 to 0.958. These results are representative of real data and reflect the model's performance. The second set of results concerns a comparative study between the methods used by Y. Xie et al (CNN, MLP, LSTM-AM and LSTM-AM-MLP), and our approach. The prediction results obtained between the approaches show a prediction gap ranging from 0.15 to 0.78 for RMSE and 0.83 for MAPE. These results show significant differences between the methods. The minimization of errors by our model, therefore, reflects the model's performance in terms of accuracy, which is necessary for the optimal management of consumer electricity load forecasts in order to ensure balance between supply and demand.

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.

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.

Trypanosoma cruzi, the protozoan causative of Chagas disease (ChD), exhibits striking genetic and phenotypic intraspecific diversity, along with ecoepidemiological complexity. Human-pathogen interactions lead to distinct clinical presentations of ChD. In 2009, an international consensus classified T. cruzi strains into six discrete typing units (DTUs), TcI to TcVI, later including TcBat, and proposed reproducible genotyping schemes for DTU identification. This article aims to review the impact of classifying T. cruzi strains into DTUs on our understanding of biological, ecoepidemio-logical, and pathogenic aspects of T. cruzi. We will explore the likely origin of DTUs and the intrinsic characteristics of each group of strains concerning genome organization, genomics, and suscepti-bility to drugs used in ChD treatment. We will also provide an overview of the association of DTUs with mammalian reservoirs, and summarize the geographic distribution, and the clinical implica-tions, of prevalent specific DTUs in ChD patients. Throughout this review, we will emphasize the crucial roles of both parasite and human genetics in defining ChD pathogenesis and chemotherapy outcome.

The rapid growth of Earth's global geospatial data necessitates an efficient system for organizing the data, facilitating data fusion from diverse sources, and promoting interoperability. Mapping the spheroidal surface of the planet presents significant challenges, as it involves balancing distortion and splitting the surface into multiple partitions. The distortion decreases as the number of partitions increases, but at the same time the complexity of data processing increases, since each partition represents a separate data set and is defined in its own local coordinate system. In this paper, we propose the dual orthogonal equidistant cylindrical projection method to mitigate distortion and reduce the number of partitions. Additionally, we use the rotation of projection cylinders to effectively minimize average angular and areal distortions of the Earth’s landmass and reduce the interruption of continental plates caused by partition edges. By incorporating auxiliary latitudes and proposing an approximate authalic latitude, we further enhance the mapping of the ellipsoid onto the sphere, simplifying calculations. Experimental results demonstrate a substantial reduction in distortion and interruption of continental plates. With only two partitions, an average landmass angular distortion of less than 3.56 degrees and an average normalized surface distortion of less than 1.07 were achieved.

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.

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.

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.

In this paper, we examine two strategies for boosting the performance of ensembles of Siamese networks (SNNs) for image classification using two loss functions (Triplet and Binary Cross Entropy) and two methods for building the dissimilarity spaces (FULLY and DEEPER). With FULLY, the distance between a pattern and a prototype is calculated by comparing two images using the fully connected layer of the Siamese network. With DEEPER, each pattern is described using a deeper layer combined with dimensionality reduction. The basic design of the SNNs takes advantage of supervised k-means clustering for building the dissimilarity spaces that train a set of support vector machines, which are then combined by sum rule for a final decision. The robustness and versatility of this approach are demonstrated on several cross-domain image data sets, including a portrait data set, two bioimage and two animal vocalization data sets. Results show that the strategies employed in this work to increase the performance of dissimilarity image classification using SNN is closing the gap with standalone CNNs. Moreover, when our best system is combined with an ensemble of CNNs, the resulting performance is superior to an ensemble of CNNs, demonstrating that our new strategy is extracting additional information.

Today, maritime transport is responsible for moving approximately 80 percent of the volume of world trade, which has favored Colombia, which, due to its geographical position, has managed to become one of the most competitive and dynamic economies in the world. South America. For this reason, this investigation is carried out in a port of their country, with the aim of identifying and analyzing restrictions or critical points that may cause delays in the unloading of goods, these would also cause delays in the loading of vehicles. To do this, it began with a diagnosis of the current situation through interviews with the operational personnel of the port, then a simulated model was designed in the ARENA software, in which it was observed that the weighing activities are where the queues were formed, this will cause a delay in the download. As possible solutions, it is recommended to senior management that a proper verification of the operation of the weighing equipment be carried out in order to unload goods more quickly. In addition, a replacement weighing equipment must be purchased. With this change, it would be possible to reduce inefficiencies or cost overruns caused by download delays, which are reducing the fluidity of the goods and the competitiveness of the organization

In this paper, Hybrid Automata, which is a formal model for hybrid systems, has been introduced. A summary of its theory is presented and some of its special and important classes are listed and some properties that can be studied and checked for it are mentioned. Finally, the purposes of use, the most widely used areas, and the tools that provide H.A. Support are addressed.

We present a generic and open-source framework for the numerical modeling of the expected transport and storage mechanisms in unconventional gas reservoirs. These unconventional reservoirs typically contain natural fractures at multiple scales. Considering the importance of these fractures in shale gas production, we perform a rigorous study on the accuracy of different fracture models. The framework is validated against an industrial simulator and is used to perform a history-matching study on the Barnett shale. This work presents an open-source code that leverages cutting-edge numerical modeling capabilities like automatic differentiation, stochastic fracture modeling, multi-continuum modeling and other explicit and discrete fracture models. We modified the conventional mass balance equation to account for the physical mechanisms that are unique to organic-rich source rocks. Some of these include the use of an adsorption isotherm, a dynamic permeability-correction function, and an embedded discrete fracture model (EDFM) with fracture-well connectivity. We explore the accuracy of the EDFM for modeling hydraulically-fractured shale-gas wells, which could be connected to natural fractures of finite or infinite conductivity, and could deform during production. Simulation results indicates that although the EDFM provides a computationally efficient model for describing flow in natural and hydraulic fractures, it could be inaccurate under these three conditions: 1. when the fracture conductivity is very low. 2. when the fractures are not orthogonal to the underlying Cartesian grid blocks, and 3. when sharp pressure drops occur in large grid blocks with insufficient mesh refinement. Each of these results are very significant considering that most of the fluids in these ultra-low matrix permeability reservoirs get produced through the interconnected natural fractures, which are expected to have very low fracture conductivities. We also expect sharp pressure drops near the fractures in these shale gas reservoirs, and it is very unrealistic to expect the hydraulic fractures or complex fracture networks to be orthogonal to any structured grid. In conclusion, this paper presents an open-source numerical framework to facilitate the modeling of the expected physical mechanisms in shale-gas reservoirs. The code was validated against published results and a commercial simulator. We also performed a history-matching study on a naturally-fractured Barnett shale-gas well considering adsorption, gas slippage & diffusion and fracture closure as well as proppant embedment, using the framework presented. This work provides the first open-source code that can be used to facilitate the modeling and optimization of fractured shale-gas reservoirs. To provide the numerical flexibility to accurately model stochastic natural fractures that are connected to hydraulically-fractured wells, it is built atop other related open-source codes. We also present the first rigorous study on the accuracy of using EDFM to model both hydraulic fractures and natural fractures that may or may not be interconnected.

Up-conversion nanoparticles (UCNP) under near-infrared (NIR) light irradiation have been well investigated in the field of bio-imaging. However, the low up-conversion luminescence (UCL) intensity limits applications. Plasmatic modulation has been proposed as an effective tool to adjust the luminescence intensity and lifetime. In this study discrete dipole approximation (DDA) was explored concerning guiding the design of UCNP@mSiO2-Au structures with enhanced UCL intensity. The extinction effects of gold shells could be changed by adjusting the distance between the UCNPs and the Au NPs by synthesized tunable mesoporous silica (mSiO2) spacers. Enhanced UCL was obtained under 808 nm irradiation. Theoretical predictions could not be demonstrated to full extend by experimental data, indicating that better models for simulation need to better take into account inhomogeneities in particle morphologies.

This article takes an approach to explaining the behavioral manifestations of the decision making in US companies’ offer of health insurance that is grounded not only on their cost minimizing behavior, but also in a humanness dimension based on the African concept of Ubuntu. In this way, we define an Ubuntu based Random Utility modeling framework, describing the choice process as a tripartite decision making, and implemented using a nationally representative random sample of 1,061 American companies from the Dunn and Bradstreet Business data, supplied by Survey Sampling International to the Associated Press-NORC Center for Public Affairs Research. The results from the three sequentially implemented specifications showed that the relationship between management culture and health plan offering strategy is dependent on other relevant co-variates, which when left out, leads to the problem of omitted variables bias. However, when all variables are included but assumed to enter the relationship exogenously, this results in management culture not having any statistically significant effect on companies' decisions about scope of health plan offering. When the exogeneity assumption is relaxed through a recursively Bivariate Probit model, the system of two equations produces a highly significant management culture effect. In fact, in this later case we see that companies with groups and formal committee management culture are 1.58 times less likely to choose a multiple plan strategy over a single plan strategy, hence failing to show the more wholesome plan offering that would theoretically prevail under Ubuntu style management.

In this paper, we proposed the fusion of two projection based face recognition algorithms: local binary Patterns in DCT domain and singular value decomposition (SVD) characterized by its simplicity and efficiently. Standard databases ORL are used to test the experimental results which prove that proposed system achieves more accurate face recognition as compared to individual method.

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.

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.

Well placement optimization in a given geological setting for a fractured geothermal reservoir is a prerequisite for enhanced geothermal operations. High computational cost associated in the framework of fully coupled thermo-hydraulic-mechanical (THM) processes in a fractured reservoir simulation, makes the well positioning as a missing point in developing a field scale investigation. Here, in this study, we shed light on this topic through examining different injection-production well (doublet) position in a given real fracture network. Water and CO2 are used as working fluids for geothermal operations and importance of well positions are examined using coupled THM numerical simulations for both the fluids. Results of this study are examined through the thermal breakthrough time, mass flux and the energy extraction potential to assess the impact of well position in a two-dimensional reservoir framework. Almost ten times of the difference between the final amount of heat extraction is observed for different well position but with the same well spacing and geological characteristics. Furthermore, stress field is be a strong function of well position that is important with respect to the possibility of unwanted stress development. As part of the MEET project, this study recommends to perform similar well placement optimization study for each fracture set in a fully coupled THM manner before a field well drilling.

The normal and Poisson distribution assumptions in the normal-Poisson mixed effects regression model are often too restrictive for many real count data. Several works have independently relaxed the Poisson conditional distribution assumption for counts or the normal distribution assumption for random effects. This work couples some recent advances in these two regards to develop a skew t–discrete gamma regression model in which the count outcomes have full dispersion flexibility and random effets can be skewed and heavy tailed. Inference in the model is achieved by maximum likelihood using pseudo-adaptive Gaussian quadature. The use of the proposal is demonstrated on a popular owl sibling negotiation data. It appears that, for this example, the proposed approach outperforms models based on normal random effects and the Poisson or negative binomial count distribution.

Establishing global high-rate secure communications is a potential application of continuous-variable quantum key distribution (CVQKD) but also challenging for long-distance transmissions in metropolitan areas. The discrete modulation(DM) can make up for the shortage of transmission distance that has a unique advantage against all side-channel attacks, however its further performance improvement requires source preparation in the presence of noise and loss. Here, we consider the effects of photon catalysis (PC) on the DM-involved source preparation for lengthening the maximal transmission distance of the CVQKD system. We address a zero-photon catalysis (ZPC)-based source preparation for enhancing the DM-CVQKD system. The statistical fluctuation due to the finite length of data is taken into account for the practical security analysis. Numerical simulations show that the ZPC-based DM-CVQKD system can not only achieve the extended maximal transmission distance, but also contributes to the reasonable increase of the secret key rate. This approach enables the DM-CVQKD to tolerate lower reconciliation efficiency, which may promote the practical implementation solutions compatible with classical optical communications using state-of-the-art technology.

In this paper, using the discrete time model, we consider the average age of all files for a cached-files-updating system where a server generates N files and transmits them to a local cache. In order that the cached files are fresh, in each time slot the server updates files with certain probabilities. The age of one file or its age of information (AoI) is defined as the time the file stays in cache since it was last time sent to cache. Assume that each file in cache has corresponding request popularity. In this paper, we obtain the distribution function of the popularity-weighted average age over all files, which gives a complete description of this average age. For the random age of single file, both the mean and its distribution have been derived before by establishing a simple Markov chain. Using the same idea, we show that an N dimensional stochastic process can be constituted to characterize the changes of N file ages simultaneously. By solving the steady-state of the resulting process, we obtain the explicit expression of stationary probability for an arbitrary state-vector. Then, the distribution function of the popularity-weighted average age can be derived by mergering a proper set of stationary probabilities. For the possible applications, the distribution function can be utilized to calculate the probability that the average age violates certain statistical guarantee.

This paper presents a discrete particle swarm optimization (DPSO) algorithm with heterogeneous (non-uniform) parameter values for solving the dynamic travelling salesman problem (DTSP). The DTSP can be modelled as a sequence of static sub-problems, each of which is an instance of the TSP. We present a method for automatically setting the values of the DPSO parameters without three parameters, which can be defined based on the size of the problem, the size of the particle swarm, the number of iterations, and the particle neighbourhood size. We show that the diversity of parameter values has a positive effect on the quality of the generated results. We compare the performance of the proposed heterogeneous DPSO with two ant colony optimization (ACO) algorithms. The proposed algorithm outperforms the base DPSO and is competitive with the ACO.

Complex human behavior, including interlimb and interpersonal coordination, has been studied from a dynamical system perspective. We review the applications of a dynamical system approach to a sporting activity, which includes continuous, discrete, and switching dynamics. Continuous dynamics identified switching between in- and anti-phase synchronization, controlled by an interpersonal distance of 0.1 m during expert kendo matches, using a relative phase analysis. As discrete dynamics, return map analysis was applied to the time series of movements during kendo matches, and six coordination patterns were classified. Furthermore, state transition probabilities were calculated based on the two states, which clarified the coordination patterns and switching behavior. We introduced switching dynamics with temporal inputs to clarify the simple rules underlying the complex behavior corresponding to switching inputs in a striking action as a non-autonomous system. As a result, we determined that the time evolution of the striking action was characterized as fractal-like movement patterns generated by a simple Cantor set rule with rotation. Finally, we propose a switching hybrid dynamics to understand both court-net sports, as strongly coupled interpersonal competition, and weakly coupled sports, such as martial arts.

The purpose of this work was analysing the compaction of a cohesive material using different DEM simulators to determine the equivalent contact models and identify how some parameters of the simulations affect the compaction results (maximum force and compacts appearance) and computational costs. For that purpose, three cohesion contact models were tested (‘linear cohesion’ in EDEM; ‘SJKR’ and ‘SJKR2’ in LIGGGHTS). The influence of the particle size distribution (PSD) on the results was also investigated. Further assessments were performed on the effect of selecting different timesteps, using distinct conversion tolerances for exporting the 3D models to STL files and moving the punch with different speeds. Consequently, it was possible to determine that a timestep equal to a 10% Rayleigh timestep, a conversion tolerance of 0.01 mm and a punch speed of 0.2 m/s are adequate for simulating the compaction process using the contact models in this work. In addition, the results determined that the maximum force was influenced by the PSD because of the rearrangement of the particles. The PSD was also related to the computational cost because of the number of simulated particles and their sizes. Finally, an equivalence was found between the linear cohesion and SJKR2 contact models

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.

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.

Hydraulic fracturing for oil-gas and geothermal reservoir stimulation is closely related to the generation and propagation of Mode I crack. Nonlinear deformation due to rock heterogeneity occurs at such crack tip, which causes the fracture process zone (FPZ) to form before the crack propagates unsteadily. However, the relationship between the FPZ characteristics and rock heterogeneity still remain elusive. We used three rock types common in reservoir for experimental investigation, and each of them includes two subtypes with different heterogeneity due to grain size or microstructural characteristics. Drawing on the experiment results, we calculated the FPZ size in each cracked chevron notched Brazilian disk, and we reproduced the formation process of the FPZ in marble by discrete element method. We showed that strong heterogeneity is favorable to large FPZ size, can enhance the ability of crack generation and complicate crack morphology. Coupling the Weibull distribution with fracture mechanics, the dependence of the FPZ size on heterogeneity degree can be theoretically explained, which suggests that inherent heterogeneity of rocks set physical foundation for formation of FPZ. These findings can improve our recognition to formation mechanism of the Mode I crack and provide useful guidelines for evaluating reservoir fracability.

It has been known that the performance of the High-Pressure Grinding Rolls (HPGR) varies as a function of method used to confine laterally the rolls, their diameter/length (aspect) ratio as well as their condition, if new or worn. However, quantifying these effects through direct experimentation in machines with reasonably large dimensions is not straightforward given the challenge, among others, of guaranteeing that the feed material remains unchanged. The present work couples the discrete element method (DEM) to multibody dynamics (MBD) and a novel particle replacement model (PRM) to simulate the performance of pilot-scale HPGRs grinding pellet feed. It shows that rotating side plates, in particular when fitted with studs, allow reaching more uniform forces along the bed, which also translates in a more constant product size along the rolls as well as higher throughput. It also shows that the edge effect is relatively constant with roll length, leading to substantially larger proportional edge regions for high-aspect ratio rolls. On the other hand, the product from the center region of such rolls was found to be finer when pressed at identical specific forces. Finally, rolls were found to have higher throughput, but generate a coarser product when worn following the commonly observed trapezoidal profile. The approach used in industry to compensate for roller wear by increasing the specific force and roll speed has then been demonstrated to be effective to maintain and potentially even increase product fineness and throughput, as long as the minimum safety gap is not reached.

This paper addresses the study of the stress field in composites continua with the multiscale approach of the DECM (Discrete Element modeling with the Cell Method). The analysis focuses on composites consisting of a matrix with inclusions of various shapes, to investigate whether and how the shape of the inclusions changes the stress field. The purpose is to provide a numerical explanation for some of the main failure mechanisms of concrete, which is precisely a composite consisting of a cement-based matrix and aggregates of various shapes. Actually, while extensive experimental campaigns detailed the shape effect of concrete aggregates in the past, so far it has not been possible to model the stress field within the inclusions and on the interfaces accurately. The reason for this lies in the limits of the differential formulation, which is the basis of the most commonly used numerical methods. The Cell Method (CM), on the contrary, is an algebraic method that provides descriptions up to the micro-scale, independently of the presence of rheological discontinuities or concentrated sources. This makes the CM useful for describing the shape effect of the inclusions, on the micro-scale. When used together with a multiscale approach, it also models the macro-scale behavior of periodic composite continua, without losing accuracy on the micro-scale. The DECM uses discrete elements precisely to provide the CM with a multiscale approach.

Protein-DNA interactions are critical for the successful functioning of all natural systems. The key role in these interactions is played by processes of protein search for specific sites on DNA. Although it has been studied for many years, only recently microscopic aspects of these processes became more clear. In this work, we present a review on current theoretical understanding of the molecular mechanisms of the protein target search. A comprehensive discrete-state stochastic method to explain the dynamics of the protein search phenomena is introduced and explained. Our theoretical approach utilizes a first-passage analysis and it takes into account the most relevant physical-chemical processes. It is able to describe many fascinating features of the protein search, including unusually high effective association rates, high selectivity and specificity, and the robustness in the presence of crowders and sequence heterogeneity.

This paper presents a multi-swarm PSO algorithm for the Quadratic Assignment Problem (QAP) implemented on the OpenCL platform. Our work was motivated by results of time efficiency tests performed for single-swarm algorithm implementation that showed clearly that the benefits of a parallel execution platform can be fully exploited provided the processed population is large. The described algorithm can be executed in two modes: with independent swarms or with migration. We discuss the algorithm construction as well as we report results of tests performed on several problem instances from the QAPLIB library. During the experiments the algorithm was configured to process large populations. This allowed us to collect statistical data related to values of goal function reached by individual particles. We use them to demonstrate on two test cases that although single particles seem to behave chaotically during the optimization process, when the whole population is analyzed, the probability that a particle will select a near-optimal solution grows.

Alkali Silica Reaction (ASR) is known to be a serious problem for concrete worldwide, especially in high humidity and high temperature regions. ASR is a slow process that develops over years to decades and it is influenced by changes in environmental and loading conditions of the structure. The problem becomes even more complicated if one recognizes that other phenomena like creep and shrinkage are coupled with ASR. This results in synergistic mechanisms that can not be easily understood without a comprehensive computational model. In this paper, coupling between creep, shrinkage and ASR is modeled within the Lattice Discrete Particle Model (LDPM) framework. In order to achieve this, a multi-physics formulation is used to compute the evolution of temperature, humidity, cement hydration, and ASR in both space and time, which is then used within physics-based formulations of cracking, creep and shrinkage. The overall model is calibrated and validated on the basis of experimental data available in the literature. Results show that even during free expansions (zero macroscopic stress), a significant degree of coupling exists because ASR induced expansions are relaxed by meso-scale creep driven by self-equilibriated stresses at the meso-scale. This explains and highlights the importance of considering ASR and other time dependent aging and deterioration phenomena at an appropriate length scale in coupled modeling approaches.

This paper studies the duration patterns of unemployment spells for immigrants and the determinants of unemployment’s completion into one of a number of studied labour market states in Finland. We estimate a duration model for unemployment with competing risks of its terminating into employment, labour market training or economic inactivity. Taking into account the wide period of observation and opportunities to analyse processes of labour market integration during various periods of economic development in Finland, in combination with the individualistic character of the labour careers of immigrants, this research is beneficial owing to the many various findings concerning labour market integration of immigrants. The approach undertaken in this research has a dualistic “descriptive-dynamic” character under which integration is understood as a never-ending process, which is conditioned by a time period of long-term existence and a context of solitary action. We find that transitions out of unemployment spells have a cyclical character; after every new “cycle” in unemployment, the probability of terminating unemployment decreases further. We also find that ascriptive factors make sense in the process of job-placement of immigrants from unemployment. Therefore, the gender, education and age of immigrants, as well as the effect of the period in which first unemployment occurred, potentially predict transitions out of unemployment and further labour market integration of immigrants.

This work addresses the analysis and characterization of deadlocks in discrete-event systems modeled by labeled Petri nets (LPNs) with undistinguishable and unobservable transitions. To provide a solution for the notorious problem, it is essential to present an effective characterization such that deadlock control and synthesis are technically and methodologically possible. To this end, we introduce the notion of dangerous implicit vectors (DIVs) which implicitly threaten the system deadlock-freedom. The set of dead markings is divided into two subsets: dead basis markings (DBMs) and dangerous implicit markings (DIMs). An algorithm is designed to compute the sets of DIVs and DIMs at a given basis state of a system. Moreover, by virtue of linear algebraic equations, we formulate sufficient conditions to identify the existence of blocking markings in an LPN. Finally, an algorithm is developed to construct an observed graph that is a compact representation of the reachability graph of a Petri net with respect to the existence of dead reaches. At the end of this paper, experiment results that illustrate the correctness and effectiveness of the reported solution are presented.

Analysis of non-stationary and nonlinear sound signals obtained from dynamical processes is one of the greatest challenges in signal processing. Turning machine operation is a highly dynamic process persuaded by many events such as dynamical responses, chip formations and operational conditions of machining. Traditional and widely used fast Fourier transformation and spectrogram is not suitable for processing sound signals acquired from dynamical systems as its result have significant deficiencies because of stationary assumptions and having priori basis. The relatively new technique, Discrete Wavelet Transform, which uses Wavelet decomposition and the recently developed technique, Hilbert Huang Transform (HHT) which uses Empirical Mode Decomposition have notably better properties in the analysis of nonlinear and non-stationary sound signals. This paper presents a comparative study on the decomposition of multicomponent sound signals using Empirical Mode Decomposition and Wavelet Decomposition to highlight the suitability found in HHT to analyze tool-emitted sound signals received from turning processes. Apart from the short mathematical and theoretical foundations of the transformations, this paper demonstrates their decomposition strength using an experimental case study of tool flank wear monitoring in turning. It also concludes HHT is more suitable than DWT to analyze tool-emitted sound signals received from turning processes.

This study investigates the fixed-time consensus-based attitude control problem for large telescope on-orbit construction subject to discrete-time communication, unknown tiny disturbance torque, model uncertainty, and actuator saturation under an undirected communication graph. Motivated by the need to save wireless communication resources, the communication scheme was designed to exploit interference, which can save communication resources in proportion to the number of modular unit mirrors. Then, a fixed-time observer with discrete-time communication was proposed, which also handles the unknown channel attenuation problem when exploiting information interference. Moreover, based on the backstepping procedure and the suitably defined adaptive compensation law, the fixed-time attitude consensus control law was proposed to handle unknown tiny disturbance torque, model uncertainty, and actuator saturation. Finally, comparative simulation results were provided to verify the effectiveness of the fixed-time observer and proposed fixed-time control law.

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.

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.

Magnetorheological fluids (MRFs) are smart materials consisting of micro-scale magnetizable particles suspended in a carrier fluid. The rheological properties of a MRF can be changed from a fluid-state to a solid-state upon the application of an external magnetic field. This study reports the development of a particle-level simulation code for magnetic solid spheres moving through an incompressible Newtonian carrier fluid. The numerical algorithm is implemented within an open-source finite-volume solver coupled with an immersed boundary method (FVM-IBM) to perform fully-resolved simulations. The particulate phase of the MRF is modeled using the discrete element method (DEM). The resultant force acting on the particles due to the external magnetic field (i.e., magnetostatic polarization force) is computed based on the Clausius-Mossotti relationship. The fixed and mutual dipole magnetic models are then used to account for the magnetic (MAG) interactions between particles. Several benchmark flows were simulated using the newly-developed FVM-IBM-DEM-MAG algorithm to assess the accuracy and robustness of the calculations. First, the sedimentation of two spheres in a rectangular duct containing a Newtonian fluid is computed without the presence of an external magnetic field, mimicking the so-called drafting-kissing-tumbling (DKT) phenomenon. The numerical results obtained for the DKT case study are verified against published data from the scientific literature. Second, we activate both the magnetostatic polarization and the dipole-dipole forces and resultant torques between the spheres for the DKT case study. Next, we study the robustness of the FVM-IBM-DEM-MAG solver by computing multi-particle chaining (i.e., particle assembly) in a two-dimensional (2D) domain for area volume fractions of 20% (260 particles) and 30% (390 particles) under vertical and horizontal magnetic fields. Finally, the fourth computational experiment describes the multi-particle chaining in a three-dimensional (3D) domain allowing to study fully-resolved MRF simulations of 580 magnetic particles under vertical and horizontal magnetic fields.

In this paper the D 3 dihedral logistic map of fractional order is introduced. The map 1 presents a dihedral symmetry D 3 . It is numerically shown that the construction and interpretation 2 of the bifurcation diagram versus the fractional order require special attention. The system stability 3 is determined and the problem of hidden attractors is analyzed. Also, analytical and numerical 4 results show that the chaotic attractor of integer order, with D 3 symmetries, looses its symmetry 5 in the fractional-order variant.

Sub-Saharan social-ecological systems are undergoing changes in environmental conditions, including modifications in rainfall pattern and biodiversity loss. Consequences of such changes depend on complex causal chains which call for integrated management strategies whose efficiency could benefit from ecosystem dynamic modelling. However, ecosystem models often require lots of quantitative information for estimating parameters, which is often unavailable. Alternatively, qualitative modelling frameworks have proved useful for explaining ecosystem response to perturbations, while requiring fewer information and providing more general predictions. However, current qualitative methods have some shortcomings which may limit their utility for specific issues. In this paper, we propose the Ecological Discrete-Event Network (EDEN), an innovative qualitative dynamic modelling framework based on "if-then" rules which generates many alternative event sequences (trajectories). Based on expert knowledge, observations and literature, we use this framework to assess the effect of permanent changes in surface water and herbivores diversity on vegetation and socio-economic transitions in an East African savanna. Results show that water availability drives changes in vegetation and socio-economic transitions, while herbivore functional groups had highly contrasted effects depending on the group. This first use of EDEN in a savanna context is promising for bridging expert knowledge and ecosystem modelling.

Most existing flexible count regression models allow only approximate inference. Balanced discretization is a simple method to produce a mean-parametrizable flexible count distribution starting from a continuous probability distribution. This makes easy the definition of flexible count regression models allowing exact inference under various types of dispersion (equi-, under- and overdispersion). This study describes maximum likelihood (ML) estimation and inference in count regression based on balanced discrete gamma (BDG) distribution and introduces a likelihood ratio based latent equidispersion (LE) test to identify the parsimonious dispersion model for a particular dataset. A series of Monte Carlo experiments were carried out to assess the performance of ML estimates and the LE test in the BDG regression model, as compared to the popular Conway-Maxwell-Poisson model (CMP). The results show that the two evaluated models recover population effects even under misspecification of dispersion related covariates, with coverage rates of asymptotic 95% confidence interval approaching the nominal level as the sample size increases. The BDG regression approach, nevertheless, outperforms CMP regression in very small samples (n = 15 − 30), mostly in overdispersed data. The LE test proves appropriate to detect latent equidispersion, with rejection rates converging to the nominal level as the sample size increases. Two applications on real data are given to illustrate the use of the proposed approach to count regression analysis.

State legitimacy and effectiveness could be seen by the way to deliver welfare to citizens to mitigate social grievances, that could eventually lead to conflicts (Kivimäki, 2021). Social security systems in Maghreb countries are quite similar in their architecture and aims to provide social insurance to all the workers in the labor market. However, they suffer from the same main problem: the low rate of enrollment of workers. Many workers (employees and self-employed) work informally without any social security coverage. The issue of whether informal jobs are chosen voluntarily by workers or as a strategy of last resort is controversial. Many authors recognize that the informal sector is heterogeneous and it is made up of workers who voluntary choose it and others who are pushed inside because of entry barriers to the formal sector (Günther & Launov, 2012). Using the SAHWA survey and discrete choice models, this article confirms the heterogeneity of the informal labor market in three Maghreb countries: Algeria, Morocco, and Tunisia. Furthermore, this article highlights the profiles of workers who voluntarily choose informality, which is missing from previous studies. Finally, this article proposes policy recommendations in order to extend social security to informal workers and to include them in the formal labour market.

This study introduces a new discrete distribution which is a weighted version of Poisson-Lindley distribution. The weighted distribution is obtained using the negative binomial weight function and can be fitted to count data with over-dispersion. The p.m.f., p.g.f. and simulation procedure of the new weighted distribution, namely weighted negative binomial Poisson-Lindley (WNBPL), are provided. The maximum likelihood method for parameter estimation is also presented. The WNBPL distribution is fitted to several insurance datasets, and is compared to the Poisson and negative binomial distributions in terms of several statistical tests.

A basic limitation of the periodogram as frequency estimator is that any of its significant peaks may be produced by a diffuse (or spread) frequency component instead of a pure one. Diffuse components are common in applications such as channel estimation in which a given periodogram peak reveals the presence of a complex multipath distribution (unresolvable propagation paths or diffuse scattering, for example). We present a method to detect the presence of a diffuse component in a given peak based on analyzing the projection of the data vector onto the span of the signature’s derivatives up to a given order. Fundamentally, a diffuse component is detected if the energy in the derivatives’ subspace is too high at the peak’s frequency, and its spread is estimated as the ratio between this last energy and the peak’s energy. The method is based on exploiting the signature’s Vandermonde structure through the properties of discrete Chebyshev polynomials. We also present an efficient numerical procedure for computing the data component in the derivatives’ span based on barycentric interpolation. The paper contains a numerical assessment of the proposed estimator and detector.

As an important civil and military infrastructure, airport runway pavement is faced up with threats from cluster munitions, since it is vulnerable to projectile impact with internal explosion. Aiming at the damage assessment of island airport runway pavement under impact, this work dealt with discrete modeling of rigid projectile penetration into concrete pavement and calcareous sand subgrade multi-layer structure. First, Discrete Element Method (DEM) is introduced to model concrete and calcareous sand granular material features like cohesive fracture and strain hardening due to compression, with mesoscale constitutive laws governing the normal and shear interactions between adjacent particles. Second, sequent DEM simulations of uniaxial and triaxial compression were performed to calibrate the DEM parameters for pavement concrete as well as subgrade calcareous sand. Piror to multi-layer structure investigations, penetration into sole concrete or calcareous sand is validated in terms of projectile deceleration and depth of penetration (DOP) with relative error

Abstract: As more and more numerical and analytical solutions to the linear neutron transport equation become available, verification of numerical results is increasingly important. This presentation concerns the development of another benchmark for the linear neutron transport equation. There are numerous ways of solving the transport equation, such as the Wiener-Hopf method based on analyticity, method of singular eigenfunctions, Laplace and Fourier transforms and analytical discrete ordinates, which is arguably one of the most straightforward, to name a few. Another potential method is the PN method, where the solution is expanded in terms of full range orthogonal Legendre polynomials and with orthogonality and truncation, the moments form a set of second order ODEs. Because of the half-range boundary conditions for incoming particles however, full range Legendre expansions are inaccurate near material discontinuities. For this reason, a double PN (DPN) expansion is more appropriate, where the incoming and exit-ing flux distributions are expanded separately to preserve the discontinuity at material interfaces. Here, a new method of solution for the DPN equations is proposed and demonstrated for an iso-tropically scattering medium.

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.

Gauss quadrature integral approximation is extended to include integrals with a measure consisting of a continuous as well as a discrete component. That is, we give an approximation for the integral of a function plus its sum over a discrete weighted set.

With growing city tourism, there is an increasing need for urban travel demand models to consider traffic generated by visitors. Existing research has concentrated on socio-demographic and journey-related factors to determine what influences the mode choice of tourists. In contrast, revealed preference data, like travel time, is almost never considered. In this article, we present the results of discrete choice modeling of city tourists’ mode choice based on revealed preference data from a survey we conducted in Kassel, Germany. We used multinomial logit models and determined the model parameters using maximum likelihood estimations. Surprisingly, travel time played a smaller role in mode choice than understood from previously established knowledge about everyday mobility. In the final model, travel time was only significant for the alternative walking. Also, most other sociodemographic and journey-related variables showed no significant influence. The final model reproduced the mode choice, but the goodness of fit was lower than expected from other research. We conclude that modeling the travel behavior of tourists is more complex than everyday mobility. An alternative approach that we suggest would be to model trip chains rather than single trips.

The high-power Asymmetric half-bridge Converter (AHBC) LED constant current driver controlled by digital current mode is a fourth-order system. Static operating point, parasitic resistance, load characteristics, sampling effect, modulation mode and loop delay will have great influence on its dynamic performance. In this paper, the small-signal pulse transfer function of the driver is established by the discrete-time modeling method for the two operating points corresponding to the three modulation modes of the trailing edge, leading edge and double edge. And, the effects of parasitic parameters, delay effect, sampling effect and load effect are fully considered in modeling. For a large number of complex exponential matrix operations, the first order Taylor formula is used for approximate calculation after the coefficient matrix is obtained by substituting the data. Then, Matlab software is used to compare and analyze the discrete-time model and the discrete-average model. The results show that the proposed discrete-time model can more accurately characterize the resonant peak and high-frequency dynamic characteristics, and is very suitable for the design of high frequency digital controller.

Investigation into geological storage of CO2 is underway at the Technology Development Plant (TDP) at Hontomín (Burgos, Spain), the only current onshore injection site in the European Union. The storage reservoir is a deep saline aquifer located within Low Jurassic Formations (Lias and Dogger), formed by fractured carbonates with low matrix permeability. Understanding the processes involved in CO2 migration within this kind of low-primary permeability carbonates influenced by fractures and faults is key to ensure safe operation and reliable plume prediction. During the hydraulic characterization tests, 2300 tons of liquid CO2 and 14000 m3 of synthetic brine were co-injected on site in various sequences to characterize the pressure response of the seal-storage pair [de Dios et al, 2017] The injection tests were analyzed with a compositional dual media model which accounts for both temperature effects (as the CO2 is liquid at the bottom of the wellbore) and multiphase flow hysteresis (to effectively simulate the alternating brine and CO2 injection tests that were performed). The pressure and temperature responses of the storage formation were history-matched mainly through the petrophysical characteristics of the fracture network [Le Gallo et al, 2017]. The dynamic characterization of the fracture network dominates the CO2 migration while the matrix does not appear to significantly contribute to the storage capacity. This initial modeling approach was improved using the characterization workflow developed within the European FP7 CO2ReMove project for sandstone fractured reservoirs [Ringrose et al., 2011; Deflandre et al., 2011]. Fractured reservoirs are challenging to handle because of their high level of heterogeneity that conditions the reservoir behaviour during the injection. In particular, natural fractures have a significant impact on well performance [Ray et al, 2012]. Furthermore, the understanding of the processes involved in CO2 migration within relatively low-permeability storage influenced by fractures and faults is essential for enabling safe storage operation [Iding and Ringrose, 2010]. As part of the European H2020 ENOS project, the site geological model is updated by integration of the recently acquired data such as the image log interpretations from injection and observation wells. The geological model is generated through the analysis and integration of data including borehole images and well test data. Following a methodology developed for naturally fractured hydrocarbon reservoirs [Ray et al., 2012], the image log analysis identified two sets of diffuse fractures. A Discrete Fracture Network [Bourbiaux et al., 2005] was built around both wells which encompass the caprock, storage and underburden formations. The fracture characteristics of the two sets of diffuse fractures, such as orientations, densities and conductivities, are calibrated upon the interpretation of the injection tests [Le Gallo et al, 2017]. For each facies, the DFN characteristics were then upscaled and propagated to the full-field reservoir simulation model as 3D fracture properties (fracture porosity, fracture permeability and equivalent block size).

Experimental recordings of the collective activity of interacting spiking neurons exhibit random behavior and memory effects, thus the stochastic process modeling the spiking activity is expected to show some degree of time irreversibility. We use the thermodynamic formalism to build a framework, in the context of spike train statistics, to quantify the degree of irreversibility of any parametric maximum entropy measure under arbitrary constraints, and provide an explicit formula for the information entropy production of the inferred Markov maximum entropy process. We provide examples to illustrate our results and discuss the importance of time irreversibility for modeling the spike train statistics.

Human number understanding is thought to rely on the analogue number system (ANS), working according to Weber’s law. We propose an alternative account, suggesting that symbolic mathematical knowledge is based on a discrete semantic system (DSS), a representation that stores values in a semantic network, similar to the mental lexicon or to a conceptual network. Here, focusing on the phenomena of numerical distance and size effects in comparison tasks, first we discuss how a DSS model could explain these numerical effects. Second, we demonstrate that DSS model can give quantitatively as appropriate a description of the effects as the ANS model. Finally, we show that symbolic numerical size effect is mainly influenced by the frequency of the symbols, and not by the ratios of their values. This last result suggests that numerical distance and size effects cannot be caused by the ANS, while the DSS model might be the alternative approach that can explain the frequency-based size effect.