Euler sums are alternating (or level two) extension of multiple zeta values (MZVs). Kaneko and Tsumura initiated the study of multiple T-values (MTVs), another level two generalization, by restricting the summation indices in the definition of MZVs to a fixed parity pattern. In this paper, we shall study finite MTVs and their alternating versions which are level two and level four variations of finite MZVs, respectively. We conjecture that all finite MZVs are in the Q-span of finite MTVs which in turn apparently lie in the span of finite Euler sums, and the inclusions are both proper. We shall first provide some structural results for Euler sums of small weights, guided by the author's previous conjecture that the finite Euler sum space of weight w is isomorphic to a quotient Euler sum space of weight w. Then, by utilizing some well-known properties of the classical alternating MTVs, we shall derive a few important Q-linear relations among the finite alternating MTVs, including the reversal, linear shuffle and sum relations. We then compute the upper bound for the dimension of the Q-span of weight w finite (alternating) MTVs for w<9, both rigorously using the newly discovered relations and numerically aided by computers.

This paper presents an approach to the vibration analysis of axially functionally graded (AFG) non-prismatic Euler-Bernoulli beams using the finite difference method (FDM). The characteristics (cross-sectional area, moment of inertia, elastic moduli, and mass density) of AFG beams vary along the longitudinal axis. The FDM is an approximate method for solving problems described with differential equations. It does not involve solving differential equations; equations are formulated with values at selected points of the structure. In addition, the boundary conditions and not the governing equations are applied at the beam’s ends. In this paper, differential equations were formulated with finite differences, and additional points were introduced at the beam’s ends and at positions of discontinuity (supports, hinges, springs, concentrated mass, spring-mass system, etc.). The introduction of additional points allowed us to apply the governing equations at the beam’s ends and to satisfy the boundary and continuity conditions. Moreover, grid points with variable spacing were also considered, the grid being uniform within beam segments. Vibration analysis of AFG non-prismatic Euler-Bernoulli beams was conducted with this model, and natural frequencies were determined. Finally, a direct time integration method (DTIM) was presented. The FDM-based DTIM enabled the analysis of forced vibration of AFG non-prismatic Euler-Bernoulli beams, considering the damping. The results obtained in this paper showed good agreement with those of other studies, and the accuracy was always increased through a grid refinement.

This paper treats the connection problem of expressing sums of finite products of Chebyshev polynomials of the third and fourth kinds in terms of five classical orthogonal polynomials. In fact, by carrying out explicit computations each of them are expressed as linear combinations of Hermite, generalized Laguerre, Legendre, Gegenbauer, and Jacobi polynomials which involve some terminating hypergeometric functions ${}_2 F_0, {}_2 F_1$, and ${}_3 F_2$.

Here we consider sums of finite products of Chebyshev polynomials of the third and fourth kinds. Then we represent each of those sums of finite products as linear combinations of the four kinds of Chebyshev polynomials which involve the hypergeometric function ₃F₂.

Two-photon imaging of light stimulus-evoked neuronal activity has been used to study all neuron classes in the vertebrate retina, from the photoreceptors to the retinal ganglion cells. Clearly, the ability to study retinal circuits down to the level of single synapses or zoomed out at the level of complete populations of neurons, has been a major asset in our understanding of this beautiful circuit. In this chapter, we discuss the possibilities and pitfalls of using an all-optical approach in this highly light-sensitive part of the brain.

This paper presents an approach to the Euler-Bernoulli beam theory (EBBT) using the finite difference method (FDM). The EBBT covers the case of small deflections, and shear deformations are not considered. The FDM is an approximate method for solving problems described with differential equations. The FDM does not involve solving differential equations; equations are formulated with values at selected points of the structure. Generally, the finite difference approximations are derived based on fourth-order polynomial hypothesis (FOPH) and second-order polynomial hypothesis (SOPH) for the deflection curve; the FOPH is made for the fourth and third derivative of the deflection curve while the SOPH is made for its second and first derivative. In addition, the boundary conditions and not the governing equations are applied at the beam’s ends. In this paper, the FOPH was made for all of the derivatives of the deflection curve, and additional points were introduced at the beam’s ends and positions of discontinuity (concentrated loads or moments, supports, hinges, springs, etc.). The introduction of additional points allowed us to apply the governing equations at the beam’s ends and to satisfy the boundary and continuity conditions. Moreover, grid points with variable spacing were also considered, the grid being uniform within beam segments. First-order analysis, second-order analysis, and vibration analysis of structures were conducted with this model. Furthermore, tapered beams were analyzed (element stiffness matrix, second-order analysis). Finally, a direct time integration method (DTIM) was presented. The FDM-based DTIM enabled the analysis of forced vibration of structures, with damping taken into account. The results obtained in this paper showed good agreement with those of other studies, and the accuracy was increased through a grid refinement. Especially in the first-order analysis of uniform beams, the results were exact for uniformly distributed and concentrated loads regardless of the grid. Further research will be needed to investigate polynomial refinements (higher-order polynomials such as fifth-order, sixth-order…) of the deflection curve; the polynomial refinements aimed to increase the accuracy, whereby non-centered finite difference approximations at beam’s ends and positions of discontinuity would be used.

In this paper, we study sums of finite products of Chebyshev polynomials of the first kind and Lucas polynomials and represent each of them in terms of Chebyshev polynomials of all kinds. Here the coefficients involve terminating hypergeometric functions ₂F₁ and these representations are obtained by explicit computations.

This paper is concerned with representing sums of finite products of Chebyshev polynomials of the second kind and of Fibonacci polynomials in terms of several classical orthogonal polynomials. Indeed, by explicit computations each of them is expressed as linear combinations of Hermite, generalized Laguerre, Legendre, Gegenbauer and Jacobi polynomials which involve the hypergeometric functions ${}_1 F_1$ and ${}_2 F_1$.

We evaluate various sums involving the powers of Fibonacci and Lucas numbers.

In this paper, we give explicit asymptotic formulas for some sums over primes involving generalized alternating hyperharmonic numbers of types I, II and III. Analogous results for numbers with $k$-prime factors will also be considered.

We evaluate some binomial cubic Fibonacci sums.

We evaluate various binomial sums involving the powers of Fibonacci and Lucas numbers.

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.

We study on the sums of powers of consequtive integers and alternating sums of power of consequtive integers. We derive many identities and correlations including Bernoulli, Euler and Genocchi polynomials and numbers.

Let {ξ1, ξ2. . . .} be a sequence of independent possibly differently distributed random variables defined on a probability space (Ω, F, P) with distribution functions {F_ξ1., F_ξ2 , . . .}. Let η be a counting random variable independent of sequence {ξ1, ξ2. . . .}. In this paper, we find conditions under which distribution function of randomly stopped sum Sη = ξ1 + ξ2 + . . . + ξη belongs to the class of generalized subexponential distributions.

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.

Frames are more stable as compared to bases under the action of a bounded linear operator. Sums of different frames under the action of a bounded linear operator are studied with the help of analysis, synthesis and frame operators. A simple construction of frames from the existing ones under the action of such an operator is presented here. It is shown that a frame can be added to its alternate dual frames, yielding a new frame. It is also shown that the sum of a pair of orthogonal frames is a frame. This provides an easy construction of a frame where the frame bounds can be computed easily. Moreover, for a pair of orthogonal frames, the necessary and sufficient condition is presented for their alternate dual frames to be orthogonal. This allows for an easy construction of a large number of new frames.

Partial sums of analytic univalent functions and partial sums of starlike have been investigated extensively by several researchers. In this paper, we investigate a partial sums of convex harmonic functions that are univalent and sense preserving in the open unit disk.

In this paper, we investigate the recently introduced degenerate Bernstein polynomials and operators and derive some of their properties. Also, we give some properties of the degenerate Euler numbers and polynomials and their connection with the degenerate Euler polynomials.

In the present paper, we introduce (p, q)-extension of Apostol type Frobenius-Euler polynomials and numbers and investigate some basic identities and properties for these polynomials and numbers, including addition theorems, difference equations, derivative properties, recurrence relations and so on. Then, we provide integral representations, explicit formulas and relations for these polynomials and numbers. Moreover, we discover (p, q)-extensions of Carlitz's result [L. Carlitz, Mat. Mag., 32 (1959), 247–260] and Srivastava and Pintér addition theorems in [H. M. Srivastava, A. Pinter, Appl. Math. Lett., 17 (2004), 375–380].

In this paper, we mainly show that generalized hyperharmonic number sums with reciprocal binomial coefficients can be expressed in terms of classical (alternating) Euler sums, zeta values and generalized (alternating) harmonic numbers.

New q-series in the spirit of Jacobi have been found in a publication rst published in 1884 written in Russian and translated into English in 1928. This work was found by chance and appears to be almost totally unknown. From these entirely new q-series, fresh lattice sums have been discovered and are presented here.

Motivated by Kim-Kim [19] introduced the new type of degenerate poly- Bernoulli polynomials by means of the degenerate polylogarithm function. In this paper, we define the degenerate poly-Frobenius-Euler polynomials, called the new type of degenerate poly-Frobenius-Euler polynomials, by means of the degenerate polylogarithm function. Then, we derive explicit expressions and some identities of those numbers and polynomials.

In this paper, we construct the degenerate poly-Frobenius-Genocchi polynomials, called the type 2 degenerate poly-Frobenius-Euler polynomials, by means of polyexponential function. We derive explicit expressions and some identities of those polynomials. In the last section, we introduce type 2 degenerate unipoly-Frobenius-Genocchi polynomials by means of unipoly function and derive explicit multifarious properties.

The simplest elasticity model of foundation underlying a slender beam under flexure was conceived by Winkler, requiring local proportionality between soil reactions and beam deflection. Such an approach leads to well-posed elastostatic and elastodynamic problems, but, as highlighted by Wieghardt, it provides elastic responses which are not technically significant for a wide variety of engineering applications. Thus, Winkler's model was replaced by Wieghardt himself by assuming that the beam deflection is the convolution integral between soil reaction field and an averaging kernel. Due to conflict between constitutive and kinematic compatibility requirements, the corresponding elastic problem of an inflected beam resting on Wieghardt foundation results to be ill-posed. Modifications of the original Wieghardt model were proposed by introducing fictitious boundary concentrated forces of constitutive type, which are physically questionable, being significantly influenced on prescribed kinematic boundary conditions. Inherent difficulties and issues are overcome in the present research using a displacement-driven nonlocal integral strategy got by swapping input and output fields involved in Wieghardt's original formulation. That is, nonlocal soil reaction fields are output of integral convolutions of beam deflection fields with an averaging kernel. Equipping the displacement-driven nonlocal integral law with the bi-exponential averaging kernel, an equivalent nonlocal differential problem, supplemented with non-standard constitutive boundary conditions involving nonlocal soil reactions, is established. As a key implication, the integro-differential equations governing the elastostatic problem of an inflected elastic slender beam resting on displacement-driven nonlocal integral foundation are replaced with much simpler differential equations supplemented with kinematic, static and new constitutive boundary conditions. The proposed nonlocal approach is illustrated by examining and analytically solving exemplar problems of structural engineering. Benchmark solutions for numerical analyses are also detected.

In the paper, by a very simple approach, the author establishes an expression in terms of a lower Hessenberg determinant for the Euler polynomials. By the determinantal expression, the author finds a recurrence relation for the Euler polynomials. By the way, the author derives the corresponding expression and recurrence relation for the Euler numbers.

Jacobian theta function identities are related to the Addition formula for the Weierstrass 1 Sigma function. In this paper, from the Addition formula for Weierstrass Sigma function, new and 2 simple proofs of four square sums theorem are given. Using the addition formula of the theta func- 3 tions, we also give simple proofs of two square sums and eight square sums theorem.

In this paper, we consider a new class of polynomials which is called the multi-poly-Euler polynomials. Then, we investigate their some properties and relations. We provide that the type 2 degenerate multi-poly-Euler polynomials equals a linear combination of the degenerate Euler polynomials of higher order and the degenerate Stirling numbers of the first kind. Moreover, we provide an addition formula and a derivative formula. Furthermore, in a special case, we acquire a correlation between the type 2 degenerate multi-poly-Euler polynomials and degenerate Whitney numbers.

In the present paper, the (p,q)-Hermite based Apostol type Frobenius-Euler polynomials and numbers are firstly considered and then diverse basic identities and properties for the mentioned polynomials and numbers, including addition theorems, difference equations, integral representations, derivative properties, recurrence relations. Moreover, we provide summation formulas and relations associated with the Stirling numbers of the second kind.

In this paper, closed forms of the summation formulas for generalized Fibonacci numbers are presented. As special cases, we give summation formulas of Fibonacci, Lucas, Pell, Pell-Lucas, Jacobsthal, Jacobsthal-Lucas numbers.

The Finite size Lyapunov exponent (FSLE) has been used extensively since the late 1990’s to diagnose turbulent regimes from Lagrangian experiments and to detect Lagrangian coherent structures in geophysical flows and two-dimensional turbulence. Historically, the FSLE was defined in terms of its computational method rather than via a mathematical formulation, and the behavior of the FSLE in the turbulent inertial ranges is based primarily on scaling arguments. Here we propose an exact definition of the FSLE based on conditional averaging of the finite amplitude growth rate (FAGR) of the particle pair separation. With this new definition, we show that the FSLE is a close proxy for the inverse structural time, a concept introduced a decade before the FSLE. The (in)dependence of the FSLE on initial conditions is also discussed, as well as the links between the FAGR and other relevant Lagrangian metrics, such as the finite time Lyapunov exponent and the second order velocity structure function.

Kim-Kim [12] introduced the new type of degenerate Bernoulli numbers and polynomials arising from the degenerate logarithm function. In this paper, we introduce a new type of degenerate poly-Euler polynomials and numbers, are called degenerate poly-Euler polynomials and numbers, by using the degenerate polylogarithm function and derive several properties on the degenerate poly-Euler polynomials and numbers. In the last section, we also consider the degenerate unipoly-Euler polynomials attached to an arithmetic function, by using the degenerate polylogarithm function and investigate some identities of those polynomials. In particular, we give some new explicit expressions and identities of degenerate unipoly polynomials related to special numbers and polynomials.

In the paper, by virtue of the Fa`a di Bruno formula, some properties of the Bell polynomials of the second kind, and the inversion formulas of binomial numbers and the Stirling numbers of the first and second kinds, the authors simplify meaningfully and significantly coefficients in two families of ordinary differential equations associated with higher order Frobenius–Euler numbers.

The objective of this study was to analyze the in vitro stability and toxicity of liposomes containing guarana in skin cell lines. The liposomes were produced by the reverse phase evaporation method containing 1 mg/mL guarana. The stability of the liposomes was evaluated by physical-chemical parameters for up to 90 d using three different storage conditions. The cytotoxicity of guarana (GL), liposomes (B-Lip), and guarana-loaded liposomes (G-Lip) was evaluated on spontaneously immortalized human keratinocyte cell lines (HaCaT), murine swiss albino fibroblasts (3T3), and human fibroblast (1BR.3.G). The evaluation was performed using cellular viability analysis. The techniques used were 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) and neutral red capturing (NRU) and the analyses were conducted after 24, 48, and 72 h exposure of these cells to the different treatments. The G-Lip exhibited physical-chemical stability for 60 d when the samples were stored in a refrigerator. The GL, B-Lip, and G-Lip demonstrated low cytotoxicity in the three cell cultures tested since a low reduction in cell viability was only observed at the highest concentrations. In addition, greater cell damage was observed for B-Lip; however, guarana protects cells from this damage. Thus, G-Lip structures can be considered as promising systems for topical applications.

Traditional discrimination techniques for Euler deconvolution use only the color spectrums of structural indexes, without considering the spatial distribution characteristics and inherent relationships among the Euler solutions to separate adjacent causative sources. In the present study, a new approach was developed for discriminating uncorrelated Euler solutions from coherent solutions based on the focusing levels indicated by the probability density distributions generated using multivariate kernel density estimations (KDE). A novel multiple coverage technique was proposed by using a series of different sized moving windows over gridded gravity data, which formed tight clusters of Euler solutions for different sized causative sources. The results of the probability density distributions were obtained using a 3-D KDE method for the Euler solution subsets {x, y, z} of synthetical models, and real data from a survey conducted in British Columbia (Canada) which had successfully established more credible and meaningful geological models when compared with three other subsets.

We evaluate the nested sum $\sum_{a_{n - 1} = c}^{a_n } {\sum_{a_{n - 2} = c}^{a_{n - 1} } { \cdots \sum_{a_0 = c}^{a_1 } {x^{a_0 } } } }$ where $a_n$ and $c$ are any integers and $x$ is a real or complex variable. Consequently, we evaluate multiple sums involving the terms of a general second order sequence, the Horadam sequence (W_j (a, b; p, q)), defined for all non-negative integers j by the recurrence relation W₀ = a, W₁ = b; W_j = pW_j−1 − qW_j−2 (j≥2); where a, b, p and q are arbitrary complex numbers, with p ≠ 0, q ≠ 0.

In statistics and probability theory, one of the most important statistics is the sums of random variables. After introducing a probability distribution, determining the sum of n independent and identically distributed random variables is one of the interesting topics for authors. This paper presents the probability density functions for the sum of n independent and identically distributed random variables such as Shanker, Akash, Ishita, Pranav, Rani, and Ram Awadh. In order to determine all aforementioned distributions, the problem-solving methods are applied which is based on the change-of-variables technique.

Industrial manipulators and parallel robots are often used for tasks like drilling or milling, that require three translational, but only two rotational degrees of freedom (“3T2R”). While kinematic models for specific mechanisms for these tasks exist, a general kinematic model for parallel robots is still missing. This paper presents the definition of the rotational component of kinematic constraints equations for parallel robots based on two reciprocal sets of Euler angles for the end-effector orientation and the orientation residual. The method allows to completely remove the redundant coordinate in 3T2R tasks and to solve the inverse kinematics for general serial and parallel robots with the gradient-descent algorithm. The functional redundancy of robots with full mobility is exploited using nullspace projection.

In the paper, we develop some numerical schemes including finite difference scheme and finite volume scheme for the fractional Laplacian operator, and apply the resulted numerical schemes to solve some fractional diffusion equations. First, the fractional Laplacian operator can be characterized as the weak singular integral by a integral operator with zero boundary condition. Second, because of the solution of fractional diffusion equation is usually singular near the boundary, we apply finite difference scheme to discrete fractional Laplacian operator and fractional diffusion equation by fractional interpolation function. Moreover, it is find that the differential matrix of the above scheme is symmetric matrix, and strictly row-wise diagonally dominant in special fractional interpolation function. Third, we show finite volume scheme to discrete fractional diffusion equation by fractional interpolation function, and analyze the properties of differential matrix. Finally, the numerical experiments are given, and verify the correctness of theoretical results and the efficiency of the schemes.

In this paper, we present sums and partial sums of Horadam sequences via generating functions which extends a recent result of Prodinger.

The two-parameter Wright special function is an interesting mathematical object that arises in the theory of the space and time-fractional diffusion equations. Moreover, many other special functions are particular instantiations of the Wright function. The article demonstrates finite representations of the Wright function in terms of sums of generalized hypergeometric functions, which in turn provide connections with the theory of the Gaussian, Airy, Bessel, and Error functions, etc. The main application of the presented results is envisioned in computer algebra for testing numerical algorithms for the evaluation of the Wright function.

We derive a concentration inequality for the uncertainty in the mean computed by stratified random sampling, and provide an online sampling method based on this inequality. Our concentration inequality is versatile and considers a range of factors including: the data ranges, weights, sizes of the strata, the number of samples taken, the estimated sample variances, and whether strata are sampled with or without replacement. Sequentially choosing samples to minimize this inequality leads to a online method for choosing samples from a stratified population. We evaluate and compare the effectiveness of our method against others for synthetic data sets, and also in approximating the Shapley value of cooperative games. Results show that our method is competitive with the performance of Neyman sampling with perfect variance information, even without having prior information on strata variances. We also provide a multidimensional extension of our inequality and discuss future applications.

This paper presents the simplified and quick way to obtain an experimental measure of the modulus of elasticity(E) of a piece of metal rod using the Smartphone microphone and through the experimental configuration in mode of unsupport beam Euler-Bernoulli(E-B). To understand it, it is required to have knowledge about flexural vibration of E-B beams and to know the analysis of the spectogram of the audio signal. The methodology consists in measuring the flexural resonance frequency through an app installed on a Smartphone and using the E formula deduced from the E-B theoretical framework. The results of the experimental measurements have a discrepancy of less than 5% with respect to the factory value, for a carbon steel rod with circular cross section using the impact excitation technique. It describes the precision and accuracy of the experimental measurement. It is useful, because it helps to classify the materials in the Faculty of Engineering and it helps to make a more efficient quality control.

In this work we study the intersection properties of a finite disk system in the euclidean space. We accomplish this by utilizing subsets of spheres with varying dimensions and analyze specific points within them, referred to as poles. Additionally, we introduce two applications: estimating the common scale factor for the radii that makes the re-scaled disks intersects in a single point, this is the \v{C}ech scale, and constructing the minimal Axis-Aligned Bounding Box (AABB) that encloses the intersection of all disks in the system.

Drawing inspiration from a recent construction of a polyhedral structure associated with an icosahedrally-symmetric map on the Riemann sphere, the article shows how to build such "dynamical polyhedra" for other icosahedral maps. First, icosahedral algebra is used to determine a special family of maps with 60 periodic critical points. The topological behavior of each map is worked out and results in a geometric algorithm that constructs a system of edges---the dynamical polyhedron---in natural correspondence to a map's topology. It turns out that the maps' descriptions fall into classes the presentation of which concludes the paper.

Besides the Lagrangian and the Eulerian descriptions, the motion of a body can also be expressed relative to the present configuration of the body, known as the relative motion description. It is interesting to consider such a relative motion description in general to formulate the basic system of field equations for solid bodies. In doing so, when the time increment from the present state is small enough, the nonlinear constitutive equations can be linearized relative to the present state so that the resulting system becomes linear. This will be done for thermoelastic materials with a brief comment on the exploitation of entropy principle in general. Relative Lagrangian formulation is based on the well-known ``small-on-large'' idea, and can be implemented for solving problems with large deformation in successive incremental manner. Some applications of such a formulation in numerical simulations are briefly reviewed.

Enhancement in chemisorption is one of the active research areas in carbon materials. To remedy the thermally degraded chemisorption at high temperatures, we report a comprehensive study of kink structures in free-standing monoatomic carbon nanowires upon heating. Our Monte Carlo simulation considers multi-monoatomic carbon chains laterally interacted by the Van der Waals force. Our study reveals that the carbon nanowires maintain their linearity regardless of chain length at low temperatures, but this is not the case at high temperatures. Disordered kink structure is observed in the short carbon chains, especially above the Peierls transition temperature. The severe kink structure might increase the possibility of attaching negatively charged atoms, thereby contributing to the development of next-generation materials for chemisorption at high temperatures. We have also provided an important inspection that any physical property of the finite-length carbon chain predicted by ab-initio calculation should reconsider the atomic rearrangement due to the thermal instability at high temperatures.

The cosmological constant problem is examined by taking an Einstein--scalar with a Higgs-type potential and scrutinizing the infrared structure induced by finite temperature effects. A variant optimal perturbation theory is implemented in the recently proposed quantum-gravitational framework. The optimized renormalized mass, i.e., the renormalized mass determined by the variant optimal perturbation theory, of the scalar field turns out to be on the order of the temperature. This shifts the cosmological constant problem to compatibility of the consequent perturbative analysis. The compatibility is guaranteed essentially by renormalization group invariance of physical quantities. We point out the resummation behind the invariance.

We employ another approach to quantize electromagnetic fields in the coordinate space, instead of the mode (or Fourier) space, such that local features of photons can be efficiently, physically, and more intuitively described. To do this, coordinate-ladder operators are defined from mode-ladder operators via the unitary transformation of systems involved in arbitrary inhomogeneous dielectric media. Then, one can expand electromagnetic field operators through the coordinate-ladder operators weighted by non-orthogonal and spatially-localized bases, which are propagators of initial quantum electromagnetic (complex-valued) field operators. Here, we call them QEM-CV-propagators. However, there are no general closed form solutions available for them. This inspires us to develop a quantum finite-difference time-domain (Q-FDTD) scheme to numerically time evolve QEM-CV-propagators. In order to check the validity of the proposed Q-FDTD scheme, we perform computer simulations to observe the Hong-Ou-Mandel effect resulting from the destructive interference of two photons in a 50/50 quantum beam splitter.

In our publication from 8 years ago1 we calculated RKKY interaction between two magnetic impurities in graphene. The consideration was based on the perturbation theory for the thermodynamic potential in the imaginary time representation and direct evaluation of real space spin susceptibility. Only the case of zero temperature was considered. We show in this short notice that the approach can be easily generalized to the case of finite temperature.

In this paper, by extending some results of informational genomics, we present a new randomness test based on the empirical entropy of strings and some properties of repeatability and unrepeatability of substrings of certain lengths. We give the theoretical motivations of our method and some experimental results of its application to a wide class of strings: decimal representations of real numbers, roulette outcomes, logistic maps, linear congruential generators, quantum measurements, natural language texts, and genomes. It will be evident that the evaluation of randomness resulting from our tests does not distinguish among the different sources of randomness (natural, or pseudo-casual).

The process of establishing secure communication between devices in the Internet of Things (IoT) can be approached by using a One-time-Pad (OTP) protocol. We propose using a known secure multiparty sum protocol for generating a One-time-Pad key using true (physical) random numbers in each device (party). We implemented the proposal using ZeroC-Ice, a middleware for distributed computing. The protocol security properties were analyzed under the assumptions of the Dolev-Yao threat model.

This paper focuses on the problem of finite-time stabilization of homogeneous, non-Lipschitz systems with dilations. A key contribution of this paper is the design of a virtual recursive Hölder, non-Lipschitz state feedback which renders the non-Lipschitz systems in the special case dominated by a lower-triangular nonlinear system, finite-time stable. The proof is based on a recursive design algorithm developed recently, to construct the virtual Hölder continuous, finite-time stabilizer as well as a C¹ positive definite and proper Lyapunov function that guarantees finite-time stability of the non-Lipschitz non-linear systems.

Sums of M consecutive squared integers (a+i)^2 equaling squared integers (for a<=1, 0<= i<= M-1) yield remarkable regular linear features when plotting values of M in function of a. These features correspond to groupings of pairs of a values for successive same values of M around straight lines of equation mu*M = 2a and are characterized in this paper for rational values of mu.

The current work aims to develop an approximation of the slice of Minkowski sum of finite number of ellipsoids, sliced up by an arbitrarily oriented plane in Euclidean space $\mathbb{R}^3$ that, to the best of the author's knowledge, has not been addressed yet. This approximation of the actual slice is in a closed form of an explicit parametric equation in the case that the slice is not passing through those zones of the Minkowski surface with high curvatures, namely the "corners". At corners with high curvatures an alternative computational algorithm is introduced in which a family of ellipsoidal inner and outer bounds of the Minkowski's sum are used to construct a "narrow strip" for the actual slice of Minkowski sum. This strip can narrow persistently for a few more number of constructing bounds to precisely coincide on the actual slice of Minkowski sum. This algorithm is also applicable to the cases of slice of Minkowski sum of ellipsoids with high aspect ratio. In line with the main goal, some ellipsoidal inner and outer bounds of the Minkowski sum are reviewed, including the so called "Kurzhanski's" bounds. Also some ellipsoidal approximations are suggested for the Minkowski sum, which they can be used, as well as the inner and outer bounds, in calculation of the suggested parametric approximation of the slice of Minkowski sum.

The object of study is a transversely bent triangular plate made of orthotropic material fixed along the edges of the plate under the action of a uniformly distributed load. Fourth-order differential equilibrium equations with variable orthotropy parameters were used. The equations were approximated by finite differences for a grid of scalene triangles. This type of mesh accurately describes the boundary contour of the triangular-shaped plates. The boundary conditions for the mesh consider the orthotropy of the plate material. Seven standard finite-difference equations have been developed considering the boundary conditions along the three edges of the plate and the presence of three angles of an irregularly shaped triangle. A finite difference matrix is obtained. The matrix structure makes calculating a triangular plate at different base angles possible. The boundary conditions were varied in the form of rigid or hinged support of triangular plates. The calculation method considers the orthotropic parameters of the material in two mutually perpendicular planes. The adaptation of the numerical method to the calculation of orthotropic plates of arbitrary shape is described. Relationships are given for determining the rigidity characteristics of orthotropic materials. An algorithm for simple engineering calculation of triangular orthotropic plates is proposed, which makes it possible to perform accurate calculations in variant design.

A simple optimization strategy for the computation of 3D finite-differencing kernels on many-cores architectures is proposed. The 3D finite-differencing computation is split direction-by-direction and exploits two level of parallelism: in-core vectorization and multi-threads shared-memory parallelization. The main application of this method is to accelerate the high-order stencil computations in numerical relativity codes.

The statistical bundle is the set of couples $(Q,W)$ of a probability density Q and a random variable W such that $\mathbb{E}$_Q [W] = 0. On a finite state space, we assume Q to be a probability density with respect to the uniform probability and give an affine atlas of charts such that the resulting manifold is a model for Information Geometry. Velocity and accelleration of a one-dimensional statistical model are computed in this set up. The Euler-Lagrange equations are derived from the Lagrange action integral. An example of Lagrangian using minus the entropy as potential energy is briefly discussed.

Recently, we have applied the generalized Littlewood theorem concerning contour integrals of the logarithm of analytical function to find the sums over inverse powers of zeroes for the incomplete Gamma- and Riemann zeta- functions, polygamma functions, and elliptical functions. Here, the same theorem is applied to study such sums for the zeroes of the Hurwitz zeta-function , including the sum over the inverse first power of its appropriately defined non-trivial zeroes. We also study some related properties of the Hurwitz zeta-function zeroes. In particular, we show that for any natural N and small real epsilon, when z tends to n=0, -1, -2… we can find at least N zeroes of zeta in the - vicinity of 0 for sufficiently small epsilon, as well as one simple zero tending to 1, etc.

Cycloid drives are generally used in precision machinery requiring high reduction ratios, such as robot joint (RV) reducers, in which the contact load and stress property of the internal meshing pairs greatly affect the dynamic transmission performance of the reducers. Therefore, in this study, Isogeometric Analysis (IGA) was employed to explore the multi-tooth contact problem of the cycloid pinwheel drive. Based on the NURBS curved surface generation method, the NURBS tooth profile of the cycloid gear was reconstructed. In addition, the NURBS surface of the cycloid gear - pin tooth - output pin was generated via the element splicing method. A geometrical analysis model of cycloid pinwheel drive was established to solve the contact force of the meshing pair under different input angles, and compared with the finite element method in terms of convergence, resultant accuracy, and solving timeliness. The results show that the Isogeometric analysis has higher accuracy and efficiency than the finite element method in calculating the contact stress and contact force. This method can improve the contact simulation accuracy of the cycloid drive and provides a reference for the design evaluation of RV reducer.

We study the atom-field interaction at finite temperature and in the dispersive regime. We show that the master equation for this sytem may be solved with the use of superoperator techniques. We calculate the linear entropy in case the field is initially in a coherent state and the atom in a superposition of its ground and excited states.

The aim of this work is to broaden the theory of computation by making use of fuzzy sets and multisets, which are generalization of ordinary sets. We study fuzzy multiset finite automata (FMFA), which are characterized by the incorporation of vagueness when it comes to the choice of the next state and multiple occurrences of the same input symbol. The resulting automaton is a genaralization of both classical automata and fuzzy automata. We discuss the algebraic characteristics of FMFA sub-automata—separatedness, strongly connectedness, cyclic, directable, and triangular FMFA in terms of equivalence classes induced by an equivalence relation defined on state set of FMFA. Intrestingly, for a given FMFA, a new FMFA can be constructed that is a homomorphic image of the original FMFA. We define different posets structures that can be associated with a given FMFA and show that some of them are upper semilattices and discuss the inter-relationship of a given poset, a finite upper semilattice (FUSL), a FMFA and a posets/FUSL associated with a given FMFA. Finally, we introduce the concept of decomposition of a FMFA in two different ways and characterize the strongly connected, triangular and directable FMFA.

In this paper, we consider the fourth-order parabolic equations with integer and fractional order time derivatives with Neumann boundary conditions. The integer order time derivatives are approximated by backward Euler difference quotients, and the fractional order time derivatives are approximated by L1 interpolation. We propose the block-centered finite difference scheme for fourth-order parabolic equations of integer and fractional order time derivatives. We prove the stability of the block-centered finite difference scheme and the second-order convergence of the discrete L2 norms of the approximate solution and its derivatives of every order. Numerical examples are given to verify the effectiveness of the block-centered finite difference scheme.

In the era of big data with increasingly complex data structures and ever-larger data scales, matrix-type data are becoming highly valued and their applications in the fields of medicine, industry, education, geography, and astronomy are growing in extent. In recent years, significant progress has been made in the practical use of matrix variable t-distribution finite mixture models for handling data in order to address the issues of multi-subgroup structures and long data tails. In this paper, the expectation-maximization (EM) algorithm with penalized maximum likelihood is proposed to resolve the problem of the unbounded nature of the likelihood function applied to the model by considering the degeneracy of the variance-covariance matrix of this model. Our data were analyzed through simulations and real data, and the results demonstrate that our model is effective in both preventing the likelihood function from being unbounded and in ensuring the accuracy of the estimated parameters of the EM algorithm.

We give three new authentication schemes without secrecy. The first two on finite fields and Galois rings, using Gray map for this link. The third construction is given on Galois rings. The main achievement in this work is to obtain optimal impersonation and substitution probabilities in the schemes. Additionally, in the first and second scheme, we simplify the source space and bring a better relationship between the size of the message space and the key space than the given in [8]. Finally, we provide a third scheme on Galois rings, which generalizes the scheme over finite fields constructed in [9].

This paper describes a new approach that can be used to determine the mechanical properties of unknown materials and complex material systems. The approach uses inverse finite element modelling (FEM) accompanied with a designed algorithm to obtain the modulus of elasticity, yield stress and strain hardening material constants of an isotropic hardening material model, as well as the material constants of the Drucker-Prager material model (modulus of elasticity, cap yield stress and angle of friction). The algorithm automatically feeds the input material properties data to finite element software and automatically runs simulations to establish a convergence between the numerical loading-unloading curve and the target data obtained from continuous indentation tests using common indenter geometries. A further module was developed to optimise convergence using an inverse FEM analysis interfaced with a non-linear MATLAB algorithm. A sensitivity analysis determined that the dual Spherical and Berkovich (S&B) approach delivered better results than other dual indentation methods such as Berkovich and Vickers (B&V) and Vickers and Spherical (V&S). It was found that better convergence values can be achieved despite a large variation in the starting parameter values and / or material constitutive model and such behaviour reflects the uniqueness of the dual S&B indentation in predicting complex material systems. The study has shown that a robust optimization method based on a non-linear least-squares curve fitting function (LSQNONLIN) within MATLAB and ABAQUS can be used to accurately predict a unique set of elastic plastic material properties and Drucker-Prager material properties. This is of benefit to the scientific investigation of properties of new materials or obtaining the material properties at different location of a part which might be not be similar due to manufacturing processes (e.g. different heating and cooling rates at different locations).

The paper aims in assessing risks of mandible fractures consequent to impacts or sport accidents. The role of the structural stiffness of mandible, related to disocclusion state, is evaluated through numerical simulations using the finite element method (FEM). It has been assumed that the quasi-static stress field, due to distributed forces developed during accidents, could explain the common types of mandibular fractures. Geometric model of adolescent mandible was built, upon the basis of medical imaging, in CAD software with distinction between cortical layer and inner spongy bone. The finite element model of disoccluded mandible was next created. Mandibular condyles were supposed jammed in the maxillary fossae. The total force of 700 N, simulating an impact on mandible, has been sequentially applied in three distinct areas: centrally, at canine zone and at the mandibular angle. Clinically most frequent fractures of mandible were recognized through the analysis of maximal principal stress and maximal principal strain fields. Mandibular fracture during accidents can be analyzed at satisfactory level using linear quasi-static FE models for designing protections in sport and transport. The proposed approach can be improved by introducing more realistic interactions between condylar processes and fossae.

The dislocation density tensors of thin elastic shells have been formulated explicitly in terms of the Riemann curvature tensor. The formulation reveals that the dislocation density of the shells is proportional to KA3=2, where K is the Gauss curvature and A is the determinant of metric tensor of the middle surface.

In this paper, we gave an attack on RSA when Euler function has small multiplicative inverse modulo "e" and the prime sum p+q is of the form p+q=2^nk_0+k_1 where n is a given positive integer and k_0 and k_1 are two suitably small unknown integers using sublattice reduction techniques and Coppersmith's methods for finding small roots of modular polynomial equations. When we compare this method with an approach using lattice based techniques, this procedure slightly improves the bound and reduces the lattice dimension.

To investigate the variation patterns of spectral characteristics of main branches in fruit trees due to uneven distribution of growth, pruning, and fruit quality, as well as to determine the optimal fruit vi-bration harvesting, a vibration differential model of tree branches was constructed using the concen-trated mass method. A three-dimensional model of the fruit tree was built using finite element software. The spectral characteristics of each branch were obtained using the linear sweep frequency method. The method of changing the mass of each branch was employed to simulate the variation in branch mass. The frequency of occurrence, mean acceleration, and coefficient of variation of the inherent frequencies were used as evaluation indicators. The results revealed that when the vibration frequency of the fruit tree was at the peak frequency in the tree trunk spectrum, the frequency of occurrence and mean accel-eration were optimal. And the coefficient of variation of branch acceleration was less than 1, indicating a moderate dispersion. Additionally, reducing the branch mass resulted in an increase in the peak fre-quency of the tree trunk spectrum. Therefore, the growth, pruning, and uneven distribution of the branches in fruit trees can alter the spectral characteristics. When harvesting fruit by vibration, selecting the peak frequency in the trunk vibration spectrum as the excitation force frequency yields the best fruit detachment.

Most of the damage to Vleg cars was caused by material fatigue. The objectives of the fatigue analysis on Xenia 2011 car wheels with the Finite Element Method using Solidworks are: 1) Knowing the fatigue analysis process on car wheels, 2) Knowing the results of the fatigue analysis on car wheels. The research method used is the Finite Element Method (FEM) method; the re-search steps include model data retrieval from the Wheel code, 3D model creation process, static simulation process according to SAE J 328 standard, mesh variations, fatigue simulation, and finally, standardization of results fatigue simulation regarding the SAE J 328 standard. The re-sults of the fatigue analysis on the Xenia 2011 car wheels obtained an age of 1000,000 cycles when compared with the minimum average of SAE J 328 of 600,000 cycles, this product design is clas-sified as a high cycle, which means it will experience a prolonged failure, and this design has a maximum alternating stress value of 627.4702 MPa, a strain of 0.007105, a displacement of 7.148 mm and a SOF of 0.087890 with a load factor used of 2.5, so overall this Xenia car can be loaded with a load of 1775 kg, if it is more than that, then Wheels will experience cracking.

A complete model was developed to simulate the behavior of a circular clamped axisymmetric fluid-coupled Piezoelectric Micromachined Ultrasonic Transducer (PMUT). Combining Finite Difference and Boundary Element Matrix (FD-BEM), this model is based on the discretization of the partial differential equation used to translate the mechanical behavior of a PMUT. In the model, both the axial and the transverse displacements are preserved in the equation of motion and used to properly define the neutral line position. To introduce fluid coupling, a Green’s function dedicated to axisymmetric circular radiating sources is employed. The resolution of the behavioral equations is used to establish the equivalent electroacoustic circuit of a PMUT that preserves the average particular velocity, the mechanical power, and the acoustic power. Particular consideration is given to verifying the validity of certain assumptions that are usually made across various steps of previously reported analytical models. In this framework, the advantages of the membrane discretization performed in the FD-BEM model are highlighted through accurate simulations of the first vibration mode and especially the cutoff frequency that many other models do not predict. This high cutoff frequency corresponds to cases where the spatial average velocity of the plate is null and is of great importance for PMUT design because it defines the upper limit above which the device is considered to be mechanically blocked. These modeling results are compared with electrical and dynamic membrane displacement measurements of AlN-based PMUTs in air and fluid. This complete PMUT model using the FD-BEM approach is shown to be very efficient in terms of computation time and accuracy.

Developing detailed Finite Element Models related to aerospace, automobile, civil, material science etc could be a very challenging task, especially generating a mesh that could produce accurate results. This review article provides an overview of different types of structural mechanics problems which are of interest to researchers and recent developments in meshing methodology related to each type of problems. The article also describes an automatic 2D mesh generation framework leverag- ing the capabilities of commercial FEA software MSC.PATRAN

The Schwarzschild metric describes a non-rotating and charge-free celestial body, and it results in things stopping at the event horizon of a black hole and spending infinite time across the event horizon by the observer far away from the black hole. The analysis of the particle’s behavior at the event horizon tells us that this solution predicts an un-expanded black hole which violates the astronomical observations and our knowledge about the black hole. Although some alternative metrics have been proposed, the singularity problem is still unsolved. In this research, the degenerate Fermi electron gas is used to reveal that the Fermi electron gas cannot shrink to a point no matter how large energy it obtains, so the singularity exists at the center very unreasonably. In order to avoid these problems, a finite-size nucleus of the black hole is proposed and reasonably explained by the behaviors of the Fermi electron gas and the Fermi neutron gas there. On the other hand, the Kerr-Newman metric is the one describing the rotating and charged black hole and the equation of the light velocity at each space point can be obtained. It tells us that there are two real and non-imaginary solutions for the radial speed of light at the position larger than the Schwarzschild radius.

We present a unified approach to the thermodynamics of hadron-quark-gluon matter at finite temperatures on the basis of a quark cluster expansion in the form of a generalized Beth-Uhlenbeck approach with a generic ansatz for the hadronic phase shifts that fulfills the Levinson theorem. The change in the composition of the system from a hadron resonance gas to a quark-gluon plasma takes place in the narrow temperature interval of 150−185 MeV where the Mott dissociation of hadrons is triggered by the dropping quark mass as a result of the restoration of chiral symmetry. The deconfinement of quark and gluon degrees of freedom is regulated by the Polyakov loop variable that signals the breaking of the Z(3) center symmetry of the color SU(3) group of QCD. We suggest a Polyakov-loop quark-gluon plasma model with O(αs) virial correction and solve the stationarity condition of the thermodynamic potential (gap equation) for the Polyakov loop. The resulting pressure is in excellent agreement with lattice QCD simulations up to high temperatures.

Infiltration based stormwater best management practices bring considerable economic, social and ecological benefits. Controlling stormwater quantity and quality are primarily important to prevent urban flooding and minimizing loads of pollutants to the receiving waters. However, there have been growing concerns about how the traditional design approach contributes to the failure of infiltration based BMP’s that have caused flooding, ponding, prolonged movement of surface water, and frequent clogging, etc. Many of these problems were due to the fact that the current design approaches of stormwater BMP’s only focus on surface hydrology and give little or no attention to the underline subsoil permeability rate and other constraints during the design and sizing process. As a result, we are exhibiting many newly constructed infiltration based BMP’s are failing to function well. This paper presents and demonstrates a new paradigm shift in designing infiltration-based stormwater BMP’s by combining subsurface hydrology and undelaying native soil constraints to establish acceptable criteria for sizing infiltration based BMPs.

Every finite simple group P can be generated by two of its elements. Pairs of generators forP are available in the Atlas of finite group representations as (not necessarily minimal) permutation representations P . It is unusual but significant to recognize that a P is a Grothendieck’s dessin d’enfant D and that a wealth of standard graphs and finite geometries G - such as near polygons and their generalizations - are stabilized by a D. In our paper, tripods P − D − G of rank larger than two,corresponding to simple groups, are organized into classes, e.g. symplectic, unitary, sporadic, etc (as in the Atlas). An exhaustive search and characterization of non-trivial point-line configurationsdefined from small index representations of simple groups is performed, with the goal to recognize their quantum physical significance. All the defined geometries G's have a contextuality parameterclose to its maximal value 1.

Knowing the mechanical properties of fibrous composite materials used in aeronautical industry from the design phase is a primary objective for designers. There are many methods for determining them, but they are laborious and require a long calculation time. And some of the methods are very approximate, giving only upper and lower bounds for these values. Experimental measurements also consume time and others resources. A method for sufficiently accurate and quickly obtained estimation of the engineering constants of the homogenized material is presented in the work. FEM is used to determine the natural frequencies of a standard bar, for which there are sufficiently precise classical methods for expressing these constants according to the properties of the homogenized material. In this way, these sought constants can be determined. In the work, Young's modulus will be determined for such a material, using the relationships that provide the proper pulsations for the longitudinal vibrations. Within the adopted model, transverse and torsional vibrations will be eliminated by blocking the nodes on the surface of the bars. In this way, more longitudinal self-pulsations can be obtained, so the precision in calculating the Young's modulus will increase.

We extend Ziv and Lempel's model of finite-state encoders to the realm of lossy compression of individual sequences. In particular, the model of the encoder includes a finite-state reconstruction codebook followed by an information lossless finite-state encoder that compresses the reconstruction codeword with no additional distortion. We first derive two different lower bounds to the compression ratio that depend on the number of states of the lossless encoder. Both bounds are asymptotically achievable by conceptually simple coding schemes. We then show that when the number of states of the lossless encoder is large enough in terms of the reconstruction block-length, the performance can be improved, sometimes significantly so. In particular, the improved performance is achievable using a random-coding ensemble that is universal, not only in terms of the source sequence, but also in terms of the distortion measure.

This paper presents a finite element analysis of the bi-directional orthogonal model, which incorporates individual crack strain separation and tracking. The objective of the research is to expand the current shear friction model to manage bi-directional cracking at any angle, allowing for a more universal model that can be applied to intricate structures and non-proportional loading cases. The proposed model was initially developed as a total strain-based model, with the assumption that crack strains are equivalent to total strains, but was subsequently recalculated to improve accuracy by separating crack strains from total strains. Furthermore, a separate crack strain formulation was created to account for strains in the concrete's uncracked portions and locked-in crack strains. The article then discusses the testing of various convergence methods and loading programs to achieve high convergence. Comparative analyses of the generalized shear friction model with other models for crack orientation and loading cases similar to those of a reinforced concrete membrane are also presented. The MATLAB program successfully applied the bi-directional cracking model for one finite element under a uniform cyclical strain state, using a secant stiffness formulation.

In this paper, the tensile behavior of unidirectional carbon/epoxy prepreg is experimentally analyzed at varying loads and temperatures, focusing on the considerable nonlinearity at the beginning of the stress-strain curve. The high viscosity of this material presented difficulties in securely holding specimens during tensile testing. As a result, modifications were made to the conventional gripping method to enable the acquisition of reliable test data. Additionally, a longer gauge length was chosen to minimize the impact of slippage on the elastic modulus measurement. To support the experiment, a micromechanical model of a prepreg tow with fiber waviness is proposed. An RVE model of periodically distributed unidirectional waved cylindrical fibers embedded within the matrix is developed to predict the effective material stiffness parameters. Numerical results are presented for different amplitude-to-wavelength ratios which indicate the fiber waviness reduces the tensile modulus of the composite. The simulation results are in good agreement with the uniaxial tensile test of the prepreg tow.

Music-domain visual programming languages (VPLs) have shown to be Turing complete. However, the common lack of built-in flow control structures can obstruct using VPLs implementing general-purpose algorithms, which harms the direct use of algorithms and algorithm theory in art creation processes using VPLs. In this article, we show how to systematically implement general-purpose algorithms in music-domain visual languages by using the Finite State Machines with Datapath computation model. The results expose a finite state machine and a set of internal state variables that walk paths whose speed can be controlled using a metronome ticks and whose path depends on the initial conditions of the algorithm. These elements can be further mapped to music elements according to the musician's intentions. We demonstrate this technique by implementing Euclid's Greatest Common Divider algorithm and using it to control high-level music elements in an implementation of Terry Riley's In C, and to control audio synthesis parameters in a FM synthesizer.

In this study, we have considered for numerical solution of a Poisson equation in a domain. we have presented a first and second order finite difference method for solving the system of the intial boundary value. We have presented the numerical solution of the Poisson equations in a two-dimensional finite region. Numerical Solution and exact solution for the presented method has been compared and the figures presented. The present error analysis tables demonstrate the efficiency of the method.

This works is related to the modeling of cellulose pyrolysis with a pseudo equilibrium approach. The objective is to model the kinetics of the cellulose pyrolysis with a semi-global mechanism obtained from the literature, in order to obtain the yield and the rate of formation of char. The pseudo equilibrium sense consists in the supposition that the solid phase devolatilization can be described kinetically - at finite rate - , preserving the competitive characteristic between the production of char and tar, while the gas phase can be described by means of chemical equilibrium. A set of ordinary, linear and non linear, differential equations was obtained and solved numerically with a simple but consistent scheme using the Totally Implicit Euler method. Chemical equilibrium was solved using CANTERA coupled with a code written in Matlab. Results showed that the scheme preserve the tar-gas competitive characteristic for cellulose pyrolysis. The gas phase, is defined as a mixture of CO₂, CO, H₂O, CH₄, H₂ and N₂ showed similarly composition compared with models from literature. Finally, the extension of the model to biomass in general, is straightforward, in order to include the hemicellulose and lignin.

The estimation of the Hubble constant in the past few decades has increasingly become more accurate with the advance of new techniques. But its value seems to depend on the epoch at which the measurements are made. The Planck estimate of the Hubble constant from the observations of the CMBR in the early Universe is about 67.

This paper presents a steady-state and transient heat conduction analysis framework using the polygonal scaled boundary finite element method (PSBFEM) with polygon/quadtree meshes. The PSBFEM is implemented with commercial finite element code Abaqus by the User Element Subroutine (UEL) feature. The detailed implementation of the framework, defining the UEL element, and solving the stiffness/mass matrix by the eigenvalue decomposition are presented. Several benchmark problems from heat conduction are solved to validate the proposed implementation. Results show that the PSBFEM is reliable and accurate for solving heat conduction problems. Not only can the proposed implementation help engineering practitioners analyze the heat conduction problem using polygonal mesh in Abaqus, but it also provides a reference for developing the UEL to solve other problems using the scaled boundary finite element method.

We provide an axiomatic base for the set of natural numbers, that has been proposed as a canonical construction, and use this definition of $\mathbb N$ to find several results on finite group theory. Every finite group $G$, is well represented with a natural number $N_G$; if $N_G=N_H$ then $H,G$ are in the same isomorphism class. We have a linear order, on the quotient space of isomorphism classes of finite groups, that is well behaved with respect to cardinality. If $H,G$ are two finite groups such that $|H|=m<n=|G|$, then $H<\mathbb Z_n\leq G\leq\mathbb Z_{p_1}^{n_1}\oplus\mathbb Z_{p_2}^{n_2}\oplus\cdots\oplus\mathbb Z_{p_k}^{n_k}$ where $n=p_1^{n_1}p_2^{n_2}\cdots p_{k}^{n_k}$ is the prime factorization of $n$. We find a canonical order for the objects of $G$ and define equivalent objects of $G$, thus finding the automorphisms of $G$. The Cayley table of $G$ takes canonical block form, and we are provided with a minimal set of independent equations that define the group. We show how to find all groups of order $n$, and order them. We give examples using all groups with order smaller than $10$, and we find the canonical block form of the symmetry group $\Delta_4$. In the next section, we extend our results to the infinite case, which defines a real number as an infinite set of natural numbers. A real function is a set of real numbers, and a sequence of real functions $f_1,f_2,\ldots$ is well represented by a set of real numbers, as well. In general, we represent mathematical objects using the smallest possible data-type. In the last section, mathematical objects are well assigned to tree structures. We conclude with brief comments on type theory and future work on computational aspects of these representations.

4-component relativistic atomic and molecular calculations are typically performed within the no-pair approximation where negative-energy solutions are discarded, hence the symmetry between electronic and positronic solutions is not considered. These states are however needed in QED calculations, where furthermore charge conjugation symmetry becomes an issue. In this work we shall discuss the realization of charge conjugation symmetry of the Dirac equation in a central field within the finite basis approximation. Three schemes for basis set construction are considered: restricted, inverse and dual kinetic balance. We find that charge conjugation symmetry can be realized within the restricted and inverse kinetic balance prescriptions, but only with a special form of basis functions that does not obey the right boundary conditions of the radial wavefunctions. The dual kinetic balance prescription is on the other hand compatible with charge conjugation symmetry without restricting the form of the radial basis functions. However, since charge conjugation relates solutions of opposite value of the quantum number κ, this requires the use of basis sets chosen according to total angular momentum j rather than orbital angular momentum ` . As a special case, we consider the free-particle Dirac equation, where the solutions of opposite sign of energy are related by charge conjugation symmetry. We note that there is additional symmetry in those solutions of the same value of κ come in pairs of opposite energy.

A finite-difference approach with non-uniform meshes was presented for simulating magnetotelluric responses in 2D structures. We presented the formulation of this scheme and gave some sights into its successful implementation, and compared finite-difference solution with known numerical results and simple analytical solutions. First, a homogeneous half-space model was tested and the finite-difference approach can provide very good accuracy for 2D magnetotelluric modeling. Then we compared to the analytical solutions for the two-layered model, the relative errors of the apparent resistivity and the impedance phase were both increased when the frequency was increased. In the end, we compare our finite-difference simulation results with COMMEMI 2D-0 model with the finite-element solutions. Both results are in close agreement to each other. These comparisons confirm the validity and reliability of our finite-difference algorithm.

Equal-Channel Angular Pressing (ECAP) is a method used to introduce severe plastic deformation into a metallic billet without changing its geometry. In special cases strain localization occurs and a pattern consisting of regions with high and low deformation (so-called shear and matrix bands) can emerge. This paper studies this phenomenon numerically adopting two-dimensional finite element simulations of one ECAP pass. The mechanical behavior of aluminum is modeled using phenomenological plasticity theory with isotropic or kinematic hardening. The effects of the two different strain hardening types are investigated numerically by systematic parameter studies: While isotropic hardening only causes minor fluctuations in the plastic strain fields, a material with high initial hardening rate and sufficient strain hardening capacity can exhibit pronounced localized deformation after ECAP. The corresponding finite element simulation results show a regular pattern of shear and matrix bands. This result is confirmed experimentally by ECAP-processing of AA6060 material in a severely cold worked condition, where microstructural analysis also reveals the formation of shear and matrix bands. Excellent agreement is found between the experimental and numerical results in terms of shear and matrix band width and length scale. The simulations provide additional insights regarding the evolution of the strain and stress states in shear and matrix bands.

Concrete-filled steel columns have been extensively used in the world due to having all suitable characteristics of concrete and steel, more ductility, increasing concrete confinement using steel wall, large energy-absorption capacity and appropriate fire behavior. In present paper, concrete-filled steel square columns have been simulated under the influence of blast load using ABAQUS software. These responses will be compared for scaled distances based on the distance to source and weight of explosive material. As result, it can be seen that although concrete deformation has been restricted using steel tube, but inner layer of concrete has been seriously damaged and column displacement will be decreased by increasing scaled distance. We also concluded that concrete-filled steel columns have high ductility and blast resistance.

In the current paper, a residual generation filter with finite memory structure is proposed for sensor fault detection. The proposed finite memory residual generation filter provides the residual by real-time filtering of fault vector using only the most recent finite observations and inputs on the window. It is shown that the residual given by the proposed residual generation filter provides the exact fault for noise-free systems. The proposed residual generation filter is specified to the digital filter structure for the amenability to hardware implementation. Finally, to illustrate the capability of the proposed residual generation filter, numerical examples are performed for the discretized DC motor system having the multiple sensor faults.

The ANSYS 2024R1 student edition was used to create a finite element model of a basketball rim and backboard. This finite element model included the use of steel for the rim and its mount, a tempered glass backboard, and an aluminum frame behind the backboard. After a mesh was created, fixed support boundary conditions were applied to the four corners of the aluminum frame. The subsequently calculated mode shapes and frequencies were compared to empirical modal analysis previously done at the United States Military Academy at West Point, New York. Four mode shapes and frequencies agreed rather well between the theoretical finite element analysis and previously published empirical modal analysis, specifically where the rim was vibrating in the vertical direction, which was the direction that the accelerometer was aligned for modal analysis. However, three mode shapes missed by the empirical modal analysis were found where the vibration of the rim was confined to the horizontal plane, which was orthogonal to the orientation of our accelerometer.

Evaluation of structural behaviour of a masonry Venetian Church with pointed vault, is presented in this paper through an analysis following the neccessary steps of a monument study. Having a detailed geometric model and the material estimation, the finite element method was used in order the influence of specific structural parts of the structure like masonry butresses and walls connections to the structural behaviour are examined. The comparison of the eigenfrequencies which are estimated by in situ measurements and finite element modal analysis is used to model identification. The response spectrum analysis, the static analysis for subsidence of some parts and the transient analysis of a specific seismic excitation are used to evaluation of the church strength. Comparison of modal analysis results with the structural pathology is helpful to this research about the structural behaviour of the church.

Financial time series and other human-driven, non-natural processes are known to exhibit fat-tailed outcome distributions. That is, such processes demonstrate a greater tendency for extreme outcomes than the normal distribution or other natural distributional processes would predict. We examine the mathematical expectation, or simply "expectation," traditionally the probability-weighted outcome, regarded since the seventeenth century as the mathematical definition of "expectation." However, when considering the "expectation" of an individual confronted with a finite sequence of outcomes, particularly existential outcomes (e.g., a trader with a limited time to perform or lose his position in a trading operation), we find this individual "expects" the median terminal outcome over those finite trials, with the classical seventeenth-century definition being the asymptotic limit as trials increase to this. Since such finite-sequence "expectations" often differ in values from the classic one, so do the optimal allocations (e.g., growth-optimal). We examine these for fat-tailed distributions. The focus is on implementation, and the techniques described can be applied to all distributional forms.

New families of direct serendipity and direct mixed finite elements on general planar, strictly convex polygons were recently defined by the authors. The finite elements of index r are H1 and H(div) conforming, respectively, and approximate optimally to order r+1 while using the minimal number of degrees of freedom. The shape function space consists of the full set of polynomials defined directly on the element and augmented with a space of supplemental functions. The supplemental functions were constructed as rational functions, which can be difficult to integrate accurately using numerical quadrature rules when the index is high. This can result in a loss of accuracy in certain cases. In this work, we propose alternative ways to construct the supplemental functions on the element as continuous piecewise polynomials. One approach results in supplemental functions that are in Hp for any p≥1. We prove the optimal approximation property for these new finite elements. We also perform numerical tests on them, comparing results for the original supplemental functions and the various alternatives. The new piecewise polynomial supplements can be integrated accurately, and therefore show better robustness with respect to the underlying meshes used.

In this paper, to make the missile intercept the maneuvering target, the parallel approaching guidance law is developed. In order to estimate the target maneuver more accurately and reduce that’s influence on guidance accuracy, the distance-scalar disturbance observer is employed. Specifically, the estimation accuracy of the observer designed is regardless of the relative distance. The finite-time prescribed performance is employed to ensure that the line-of-sight angular rate is capable of converging to a predesigned small region in the specified finite time. All signals of interception system can guarantee ultimately uniformly boundedness proved by the Lyapunov stability theory. Finally, the availability of the parallel approaching guidance law is demonstrated by the numerical simulation.

The purpose of this study was to investigate the functional load capacity of the periodontal ligament (PDL) in a full arch maxilla and mandible model using numerical simulation. The goal of was to determine the functional load pattern in multi- and single rooted teeth with full and reduced periodontal support. CBCT data were used to create 3D-models of a maxilla and mandible. The DICOM dataset was used to create a CAD-model. For a precise description of the surfaces of each structure (enamel, dentin, cementum, pulp, PDL, gingiva, bone), each tooth was segmented separately, and the biomechanical characteristics were considered. A Finite-Element-Analysis (FEA) software computed the biomechanical behavior of stepwise increased force of 700N in cranial and 350N in ventral direction of the muscle approach of the Masseter muscle. The periodontal attachment (cementum-PDL-bone contact) was subsequently reduced in 1mm increments and the simulation repeated. Quantitative (pressure, tension, and deformation) and qualitative (color-coded images) data were recorded and descriptively analyzed. The teeth with the highest load capacities were the upper and lower molars (0.4-0.6MPa), followed by the premolars (0.4-0.5MPa) and canines (0.3-0.4MPa) when vertically loaded. Qualitative data showed that the area with the highest stress in the PDL were for single rooted teeth the cervical and apical area and for molars additionally the furcation roof. In both, single- and multi rooted teeth the gradual reduction of the bone levels caused an increase of the load on the remaining PDL. Cervical and apical areas as well as the furcation roof are the zones with the highest functional stress. The more bone loss, the higher the mechanical load on the residual periodontal supporting structures.

Emerging soft actuators with various soft materials and traveling-wave motion have been frequently discussed. Various configurations and the resulting performances have been proposed, and it is very challenging for soft actuators in realizing large strokes. This study presents an experimentally validated nonlinear finite element model to predict the deformation produced by a spiral pneumatic rubber actuator to generate traveling-wave motion. The actuator consists of a membrane mounted on a rubber substrate with three air chambers in a spiral configuration, simplifying the actuator's structure. The sequential deformations of the successive chambers interact with each other and produce radial traveling-waves on the membrane surface, driving the objects placed on the actuator. Finite element analysis using ANSYS computer software analyzed the elastic movement by considering the influence of different initial structural types. The structure with three different length ratios was introduced using the four parameters: membrane thickness, chamber width, chamber depth, and wall thickness. The simulation results indicated an optimal structure with specific ratios. A reasonable correlation was obtained during experimental validation; the predicted displacement values were approximately 17% smaller than the experimental values. Finally, the transportation performance of the prototype was tested and a velocity of 2.28 mm/s in the desired direction was achieved. We expect that our demonstration will expand the range of applications of the spiral pneumatic rubber actuator to include conveying or worm-like robots.

This paper presents an approach to the two-dimensional analysis of elastic isotropic deep beams using the finite difference method (FDM). Deep beams are subjected to in-plane loading and present a shear span to height ratio of less than 2.50; consequently, Euler-Bernoulli beam theory and Timoshenko beam theory do not apply. Deep beams analysis is generally conducted using numerical methods such as the finite element method and to a lesser extent the FDM; the strut-and-tie model and the stress field method are also widely utilized. Analytical approaches usually make use of the Airy stress function, where stresses are formulated in terms of the stress function; however, the exact solution of this function satisfying all of the boundary conditions can hardly be found, even for simple cases. In this paper, deep beams were analyzed using the FDM. The FDM is an approximate method for solving problems described with differential equations. The FDM does not involve solving differential equations; equations are formulated with values at selected nodes of the structure. Therefore, the deep beam was discretized with a two-dimensional grid, and additional nodes were introduced at the boundaries and at positions of discontinuity (openings, brutal change of material properties, non-uniform grid spacing), the number of additional nodes corresponding to the number of boundary conditions at the node of interest. The introduction of additional nodes allowed us to apply the governing equations at boundary nodes and satisfy the boundary and continuity conditions. An Airy stress function approach and a displacement potential function approach were considered in this study whereby strong formulations of equations (equilibrium, kinematic, and constitutive) were set. Stress and stability analyses were carried out with this model; furthermore, deep beams of varying stiffness, layered beams, and beams having openings were analyzed. For slender beams, the results obtained with the Airy stress function approach showed good agreement with those of the Euler-Bernoulli beam theory, and for deep beams, the shapes of stress distributions were in good agreement with a proper understanding of the behavior of structures. On the other hand, the displacement potential function approach delivered unsatisfactory results, probably due to the use of an inefficient equation solver; a more powerful tool will be needed in future research for this purpose.