The effect of the earth pressure coefficients on the runout of granular material

In the framework of a better territory risk assessment and decision making, numerical simulation can provide a useful tool for investigating the propagation phase of phenomena involving granular material, like rock avalanches, when realistic geological contexts are considered.Among continuum mechanics models, the numerical model SHWCIN uses the depth averaged Saint Venant approach, in which the avalanche thickness (H) is very much smaller than its extent parallel to the bed (L). The material is assumed to be incompressible and the mass and the momentum equations are written in a depth averaged form.The SHWCIN code, based on the hypothesis of isotropy of normal stresses (σ_xx = σ_yy = σ_zz), has been modified (new code: RASH^3D) in order to allow for the assumption of anisotropy of normal stresses (σ_xx = Kσ_zz; σ_yy = Kσ_zz).A comparison among the results obtained by assuming isotropy or anisotropy is given through the back analysis of a set of laboratory experiments [Gray, J.M.N.T., Wieland, M., Hutter, K., 1999. Gravity-driven free surface flow of granular avalanches over complex basal topography. Proceedings of the Royal Society of London, Series A 455(1841)] and of a case history of rock avalanche (Frank slide, Canada).The carried out simulations have also underlined the importance of using a different earth pressure coefficient value (K) for directions of convergence and of divergence of the flux.     

Weyl Closure of a Linear Differential Operator

We study the Weyl closure Cl(L)  = K(x)〈∂〉L ∩ D for an operator L of the first Weyl algebra D = K〈x, ∂〉. We give an algorithm to compute Cl(L) and we describe its initial ideal under the order filtration. Our main application is an algorithm for constructing a Jordan–Hölder series for a holonomic D -module and a formula for its length. Using the closure, we also reproduce a result ofStrömbeck (1978), who described the initial ideals of left ideals of D under the order filtration, and a result ofCannings and Holland (1994), who described the isomorphism classes of right ideals of D.     

ANN based simulation and experimental verification of analytical four- and five-parameters models of PV modules

In this article, artificial neural network (ANN) is adopted to predict photovoltaic (PV) panel behaviors under realistic weather conditions. ANN results are compared with analytical four and five parameter models of PV module. The inputs of the models are the daily total irradiation, air temperature and module voltage, while the outputs are the current and power generated by the panel. Analytical models of PV modules, based on the manufacturer datasheet values, are simulated through Matlab/Simulink environment. Multilayer perceptron is used to predict the operating current and power of the PV module. The best network configuration to predict panel current had a 3–7–4–1 topology. So, this two hidden layer topology was selected as the best model for predicting panel current with similar conditions. Results obtained from the PV module simulation and the optimal ANN model has been validated experimentally. Results showed that ANN model provide a better prediction of the current and power of the PV module than the analytical models. The coefficient of determination (R²), mean square error (MSE) and the mean absolute percentage error (MAPE) values for the optimal ANN model were 0.971, 0.002 and 0.107, respectively. A comparative study among ANN and analytical models was also carried out. Among the analytical models, the five-parameter model, with MAPE = 0.112, MSE = 0.0026 and R² = 0.919, gave better prediction than the four-parameter model (with MAPE = 0.152, MSE = 0.0052 and R² = 0.905). Overall, the 3–7–4–1 ANN model outperformed four-parameter model, and was marginally better than the five-parameter model.     

Sylvester–Habicht Sequences and Fast Cauchy Index Computation

In this paper we show how Schönhage’s strategy for computing continued fractions (Schönhage 1971) can be combined with the theory of sub-resultants (Habicht, 1948; Collins, 1967; Brown, 1971; Brown and Traub, 1971; Loos, 1982; Gonzalez et al., 1990, 1994; Ducos, 1996; Ho and Yap, 1996; Lazard, 1998; Quitté, 1998) in order to compute the Cauchy index of a rational function or the signature of a non-singular Hankel matrix in a fast and also storage efficient way. Over the integers our algorithms have bit complexityO (M(d, σ) ·log(d)) withσ = O(dτ) where M(d,σ )  = O(dσ·log(dσ) ·loglog(dσ)) is Schönhage’s bound for multiplication of integer polynomials of degrees bounded by d and bit size bounded by σ in the multi-tape Turing machine model (Schönhage, 1982). Thus our bound isO(d²τ·log(dτ) ·loglog(dτ) ·log(d)).As a byproduct of the necessary analysis we obtain a refinement of the Sub-resultant Theorem. We present a new exact divisibility for sub-resultants in the defective case which extends the formulæ for the non-defective situation in a natural way. We also prove that the size of coefficients in the ordinary remainder sequence is quadratic in d.     

Computational assessment of several hydrogen-free high energy compounds

Tetrazino-tetrazine-tetraoxide (TTTO) is an attractive high energy compound, but unfortunately, it is not yet experimentally synthesized so far. Isomerization of TTTO leads to its five isomers, bond-separation energies were empolyed to compare the global stability of six compounds, it is found that isomer 1 has the highest bond-separation energy (1204.6 kJ/mol), compared with TTTO (1151.2 kJ/mol); thermodynamic properties of six compounds were theoretically calculated, including standard formation enthalpies (solid and gaseous), standard fusion enthalpies, standard vaporation enthalpies, standard sublimation enthalpies, lattice energies and normal melting points, normal boiling points; their detonation performances were also computed, including detonation heat (Q, cal/g), detonation velocity (D, km/s), detonation pressure (P, GPa) and impact sensitivity (h₅₀, cm), compared with TTTO (Q = 1311.01 J/g, D = 9.228 km/s, P = 40.556 GPa, h₅₀ = 12.7 cm), isomer 5 exhibites better detonation performances (Q = 1523.74 J/g, D = 9.389 km/s, P = 41.329 GPa, h₅₀ =  28.4 cm).     

Prediction of California bearing ratio (CBR) of fine grained soils by AI methods

Advances in field of artificial intelligence (AI) offers opportunities of utilizing new algorithms and models that enable researchers to solve the most complex systems. As in other engineering fields, AI methods have widely been used in geotechnical engineering. Unlikely, there seems quite insufficient number of research related to the use of AI methods for the estimation of California bearing ratio (CBR). There were actually some attempts to develop prediction models for CBR, but most of these models were essentially statistical correlations. Nevertheless, many of these statistical correlation equations generally produce unsatisfactory CBR values. However, this paper is likely one of the very first research which aims to investigate the applicability of AI methods for prediction of CBR. In this context, artificial neural network (ANN) and gene expression programming (GEP) were applied for the prediction of CBR of fine grained soils from Southeast Anatolia Region/Turkey. Using CBR test data of fine grained soils, some proper models are successfully developed. The results have shown that the both ANN and GEP are found to be able to learn the relation between CBR and basic soil properties. Additionally, sensitivity analysis is performed and it is found that maximum dry unit weight (γ_d) is the most effective parameter on CBR among the others such as plasticity index (PI), optimum moisture content (w_opt), sand content (S), clay + silt content (C + S), liquid limit (LL) and gravel content (G) respectively.     

Computer-based analysis of explicit approximations to the implicit Colebrook–White equation in turbulent flow friction factor calculation

The implicit Colebrook–White equation has been widely used to estimate the friction factor for turbulent fluid-flow in rough-pipes. In this paper, the state-of-the-art review for the most currently available explicit alternatives to the Colebrook–White equation, is presented. An extensive comparison test was established on the 20 × 500 grid, for a wide range of relative roughness (ε/D) and Reynolds number (R) values (1 × 10^?6 ? ε/D ? 5 × 10^?2; 4 × 10³ ? R ? 10⁸), covering a large portion of turbulent flow zone in Moody’s diagram. Based on the comprehensive error analysis, the magnitude points in which the maximum absolute and the maximum relative error are occurred at the pair of ε/D and R values, are observed. A limiting case of the most of these approximations provided friction factor estimates that are characterized by a mean absolute error of 5 × 10^?4, a maximum absolute error of 4 × 10^?3 whereas, a mean relative error of 1.3% and a maximum relative error of 5.8%, over the entire range of ε/D and R values, respectively. For practical purposes, the complete results for the maximum and the mean relative errors versus the 20 sets of ε/D value, are also indicated in two comparative figures. The examination results for error properties of these approximations gives one an opportunity to practically evaluate the most accurate formula among of all the previous explicit models; and showing in this way its great flexibility for estimating turbulent flow friction factor. Comparative analysis for the mean relative error profile revealed, the classification for the best-fitted six equations examined was in a good agreement with those of the best model selection criterion claimed in the recent literature, for all performed simulations.     

Computing Local Artin Maps,and Solvability of Norm Equations

Let L = K(α) be an Abelian extension of degree n of a number field K, given by the minimal polynomial of α over K. We describe an algorithm for computing the local Artin map associated with the extension L / K at a finite or infinite prime v of K. We apply this algorithm to decide if a nonzero a ∈ K is a norm from L, assuming that L / K is cyclic.     

Artificial neural network-based model for predicting VO2max from a submaximal exercise test

The goal of this study is to develop an accurate artificial neural network (ANN)-based model to predict maximal oxygen uptake (VO₂max) of fit adults from a single stage submaximal treadmill jogging test. Participants (81 males and 45 females), aged from 17 to 40 years, successfully completed a maximal graded exercise test (GXT) to determine VO₂max. The variables; gender, age, body mass, steady-state heart rate and jogging speed are used to build the ANN prediction model. Using 10-fold cross validation on the dataset, the average values of standard error of estimate (SEE), Pearson’s correlation coefficient (r) and multiple correlation coefficient (R) of the model are calculated as 1.80 ml kg^?1 min^?1, 0.95 and 0.93, respectively. Compared with the results of the other prediction models in literature that were developed using Multiple Linear Regression Analysis, the reported values of SEE, r and R in this study are considerably more accurate.     

Taylor approximation for hybrid systems

We propose a new approximation technique for Hybrid Automata. Given any Hybrid Automaton H, we call Approx(H, k) the Polynomial Hybrid Automaton obtained by approximating each formula ϕ in H with the formulae ϕ_k obtained by replacing the functions in ϕ with their Taylor polynomial of degree k. We prove that Approx(H, k) is an over-approximation of H. We study the conditions ensuring that, given any ϵ > 0, some k₀ exists such that, for all k > k₀, the “distance” between any vector satisfying ϕ_k and at least one vector satisfying ϕ is less than ϵ. We study also conditions ensuring that, given any ϵ > 0, some k₀ exists such that, for all k > k₀, the “distance” between any configuration reached by Approx(H, k) in n steps and at least one configuration reached by H in n steps is less than ϵ.     

Bounds for the Roots of Lacunary Polynomials

A polynomial P(X)  = X^d + a_d − 1X^d − 1 + ⋯ is called lacunary when a_d − 1 =  0. We give bounds for the roots of such polynomials with complex coefficients. These bounds are much smaller than for general polynomials.     

Simulation of solid thermal explosion and liquid thermal explosion of dicumyl peroxide using calorimetric technique

Dicumyl peroxide (DCPO), is produced by cumene hydroperoxide (CHP) process, is utilized as an initiator for polymerization, a prevailing source of free radicals, a hardener, and a linking agent. DCPO has caused several thermal explosion and runaway reaction accidents in reaction and storage zone in Taiwan because of its unstable reactive property. Differential scanning calorimetry (DSC) was used to determine thermokinetic parameters including 700 J g^–1 of heat of decomposition (ΔH_d), 110 °C of exothermic onset temperature (T₀), 130 kJ mol^–1 of activation energy (E_a), etc., and to analyze the runaway behavior of DCPO in a reaction and storage zone. To evaluate thermal explosion of DCPO with storage equipment, solid thermal explosion (STE) and liquid thermal explosion (LTE) of thermal safety software (TSS) were applied to simulate storage tank under various environmental temperatures (T_e). T_e exceeding the T₀ of DCPO can be discovered as a liquid thermal explosion situation. DCPO was stored under room temperature without sunshine and was prohibited exceeding 67 °C of self-accelerating decomposition temperature (SADT) for a tank (radius = 1 m and height = 2 m). SADT of DCPO in a box (width, length and height = 1 m, respectively) was determined to be 60 °C. The TSS was employed to simulate the fundamental thermal explosion behavior in a large tank or a drum. Results from curve fitting demonstrated that, even at the earlier stage of the reaction in the experiments, ambient temperature could elicit exothermic reactions of DCPO. To curtail the extent of the risk, relevant hazard information is quite significant and must be provided in the manufacturing process.     

Modal analysis and modified cascade neural networks in identification of geometrical parameters of circular arches

A hybrid computational system, composed of the finite element method (FEM) and cascade neural network system (CNNs), is applied to the identification of three geometrical parameters of elastic arches, i.e. span l, height f and cross-sectional thickness h. FEM is used in the direct (forward) analysis, which corresponds to the mapping α = {l, f, h} → {ω_j}, where: α – vector of control parameters, ω_j – arch eigenfrequencies. The reverse analysis is related to the identification procedure in which the reverse mapping is performed {ω_j} → {α_i}. For the identification purposes a recurrent, three level CNNs of structure (D^k-H^k-1)_s was formulated, where: k – recurrence step, s = I, II, III-levels of cascade system. The Semi-Bayesian approach is introduced for the design of CNNs applying the MML Maximum Marginal Likelihood) criterion. The computation of hyperparameters is performed by means of the Bayesian procedure evidence. The numerical analysis proves a great numerical efficiency of the proposed hybrid approach for both the perfect (noiseless) values of eigenfrequencies and noisy ones simulated by an added artificial noise.     

Dielectric and piezoelectric properties of [0 0 1] and [0 1 1]-poled relaxor ferroelectric PZN–PT and PMN–PT single crystals

The properties of PZN–PT and PMN–PT single crystals of varying compositions and orientations have been investigated. Among the various compositions studied, [0 0 1]-optimally poled PZN-(6–7)%PT and PMN-30%PT exhibit superior dielectric and piezoelectric properties, with K^T ≈ 6800–8000, d₃₃ ≈ 2800 pC/N, d₃₁ ≈ −(1200–1800) pC/N for PZN-(6–7)%PT; and K^T = 7500–9000, d₃₃ = 2200–2500 pC/N and d₃₁ = −(1100–1400) pC/N for PMN-30%PT. These two compositions are also fairly resistant to over-poling. The [0 0 1]-poled electromechanical coupling factors (k₃₃, k₃₁ and k_t) are relatively insensitive to crystal composition. [0 1 1]-optimally poled PZN-7%PT single crystal also exhibits extremely high d₃₁ values of up to −4000 pC/N with k₃₁ ≈ 0.90–0.96. While [0 1 1]-poled PZN-7%PT single crystal becomes over-poled with much degraded properties when poled at and above 0.6 kV/mm, PZN-6%PT crystal shows no signs of over-poling even when poled to 2.0 kV/mm. The presence of a certain amount (i.e., 10–15%) of orthorhombic phase in a rhombohedral matrix has been found to be responsible for the superior transverse piezoelectric properties of [0 1 1]-optimally poled PZN-(6–7)%PT. The present work shows that flux-grown PZN–PT crystals exhibit superior and consistent properties and improved over-poling resistance to flux-grown PMN–PT crystals and that, for or a given crystal composition, flux-grown PMN–PT crystals exhibit superior over-poling resistance to their melt-grown counterparts.     

Evaluation of response characteristics of resistive oxygen sensors based on porous cerium oxide thick film using pressure modulation method

In order to reduce the response time of resistive oxygen sensors using porous cerium oxide thick film, it is important to ascertain the factors controlling response. Pressure modulation method (PMM) was used to find the rate-limiting step of sensor response. This useful method measures the amplitude of sensor output (H(f)) for the sine wave modulation of oxygen partial pressure at constant frequency (f). In PMM, “break” response time, which is minimum period in which the sensor responds precisely, can be measured. Three points were examined: (1) simulated calculations of PMM were carried out using a model of porous thick film in which spherical particles are connected in a three-dimensional network; (2) sensor response speed was experimentally measured using PMM; and (3) the diffusion coefficient and surface reaction coefficient were estimated by comparison between experiment and calculation. The plot of log f versus log H(f) in the high f region was found to have a slope of approximately −0.5 for both porous thick film and non-porous thin film, when the rate-limiting step was diffusion. Calculations showed the response time of porous thick film was 1/20 that of non-porous thin film when the grain diameter of the porous thick film was the same as the thickness of non-porous thin film. At 973 K, “break” response time (t_b) of the resistive oxygen sensor was found by experiment to be 109 ms. It was concluded that the response of the resistive oxygen sensor prepared in this study was strongly controlled by diffusion at 923–1023 K, since the experiment revealed that the slope of plot of log f versus log H(f) in the high f region was approximately −0.5. At 923–1023 K, the diffusion coefficient of oxygen vacancy in porous ceria (D_V) was expressed as follows: D_V (m²s⁻¹) = 5.78 × 10⁻⁴ exp(−1.94 eV/kT). At 1023 K, the surface reaction coefficient (K) was found to exceed 10⁻⁴ m/s.     

Calibration and validation of hyperspectral indices for the estimation of broadleaved forest leaf chlorophyll content,leaf mass per area,leaf area index and leaf canopy biomass

This article aims at finding efficient hyperspectral indices for the estimation of forest sun leaf chlorophyll content (CHL, µg cm_leaf^? 2), sun leaf mass per area (LMA, g_{dry matter} m_leaf^? 2), canopy leaf area index (LAI, m²_leaf m_soil^? 2) and leaf canopy biomass (B_leaf, g_{dry matter} m_soil^? 2). These parameters are useful inputs for forest ecosystem simulations at landscape scale. The method is based on the determination of the best vegetation indices (index form and wavelengths) using the radiative transfer model PROSAIL (formed by the newly-calibrated leaf reflectance model PROSPECT coupled with the multi-layer version of the canopy radiative transfer model SAIL). The results are tested on experimental measurements at both leaf and canopy scales. At the leaf scale, it is possible to estimate CHL with high precision using a two wavelength vegetation index after a simulation based calibration. At the leaf scale, the LMA is more difficult to estimate with indices. At the canopy scale, efficient indices were determined on a generic simulated database to estimate CHL, LMA, LAI and Bleaf in a general way. These indices were then applied to two Hyperion images (50 plots) on the Fontainebleau and Fougères forests and portable spectroradiometer measurements. They showed good results with an RMSE of 8.2 µg cm^? 2 for CHL, 9.1 g m^? 2 for LMA, 1.7 m² m^? 2 for LAI and 50.6 g m^? 2 for Bleaf. However, at the canopy scale, even if the wavelengths of the calibrated indices were accurately determined with the simulated database, the regressions between the indices and the biophysical characteristics still had to be calibrated on measurements. At the canopy scale, the best indices were: for leaf chlorophyll content: ND_chl = (ρ₉₂₅ ? ρ₇₁₀)/(ρ₉₂₅ + ρ₇₁₀), for leaf mass per area: ND_LMA = (ρ₂₂₆₀ ? ρ₁₄₉₀)/(ρ₂₂₆₀ + ρ₁₄₉₀), for leaf area index: D_LAI = ρ₁₇₂₅ ? ρ₉₇₀, and for canopy leaf biomass: ND_Bleaf = (ρ₂₁₆₀ ? ρ₁₅₄₀)/(ρ₂₁₆₀ + ρ₁₅₄₀).     

Deadline-based scheduling in support of real-time data delivery

The use of deadline-based scheduling in support of real-time delivery of application data units (ADUs) in a packet-switched network is investigated. Of interest is priority scheduling where a packet with a smaller ratio of T/H (time until delivery deadline over number of hops remaining) is given a higher priority. We refer to this scheduling algorithm as the T/H algorithm. T/H has time complexity of O(log N) for a backlog of N packets and was shown to achieve good performance in terms of the percentage of ADUs that are delivered on-time. We develop a new and efficient algorithm, called T/H − p, that has O(1) time complexity. The performance difference of T/H, T/H − p and FCFS are evaluated by simulation. Implementations of T/H and T/H − p in high-speed routers are also discussed. We show through simulation that T/H − p is superior to FCFS but not as good as T/H. In view of the constant time complexity, T/H − p is a good candidate for high-speed routers when both performance and implementation cost are taken into consideration.