This paper presents an approach to the Kirchhoff-Love plate theory (KLPT) using the finite difference method (FDM). The KLPT 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 in the case of KLPT, the finite difference approximations are derived based on the fourth-order polynomial hypothesis (FOPH) and second-order polynomial hypothesis (SOPH) for the deflection surface. The FOPH is made for the fourth and third derivative of the deflection surface while the SOPH is made for its second and first derivative; this leads to a 13-point stencil for the governing equation. In addition, the boundary conditions and not the governing equations are applied at the plate edges. In this paper, the FOPH was made for all of the derivatives of the deflection surface; this led to a 25-point stencil for the governing equation. Furthermore, additional nodes were introduced at plate edges and at positions of discontinuity (continuous supports/hinges, incorporated beams, stiffeners, brutal change of stiffness, etc.), 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 the plate edges and to satisfy the boundary and continuity conditions. First-order analysis, second-order analysis, buckling analysis, and vibration analysis of plates were conducted with this model. Moreover, plates of varying thickness and plates with stiffeners were analyzed. 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. In first-order, second-order, buckling, and vibration analyses of rectangular plates, the results obtained in this paper were in good agreement with those of well-established methods, and the accuracy was increased through a grid refinement.

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.

This paper presents an approach to the Timoshenko beam theory (TBT) using the finite difference method (FDM). The Timoshenko beam theory covers cases associated with small deflections based on shear deformation and rotary inertia considerations. 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. The model developed in this paper consisted of formulating differential equations with finite differences and introducing additional points at the beam’s ends and at positions of discontinuity (concentrated loads or moments, supports, hinges, springs, brutal change of stiffness, spring-mass system, etc.). The introduction of additional points allowed us to apply the governing equations at the beam’s ends. Moreover, grid points with variable spacing were considered, the grid being uniform within beam segments. First-order, second-order, and vibration analyses of structures were conducted with this model. Furthermore, tapered beams were analyzed (element stiffness matrix, second-order analysis, vibration 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.

In this research paper, design and development of pyramidal horn antenna for J-band application is reported. It is particularly designed for 17 dB gain and half beam width about 25 degrees at 6.93 GHz. Horn aperture, horn axial length and distance from the throat of the antenna to aperture are the main design constraints which are calculated and used for the antenna design and simulation. Beam width in E-plane and H-plane horn is calculated and it is 19.18 dB and 22.86 dB respectively. The reported antenna design shows good performance for J-band in radiometry, satellite, and radar applications.

In this paper, we first discuss the research status and application progress of the finite element method and the smoothed particle method. By analyzing the advantages of the smoothed particle method and the finite element method, a new coupling algorithm, namely FEM-SPH algorithm, is proposed. By the method of comparison, it shows that finite element method and SPH method in the simulation of large deformation problems each have advantages and disadvantages, the finite element method smoothed particle coupling algorithm is effective to achieve the performance of high computational efficiency and can naturally simulate large deformation problems across. In the process of calculation, the large deformation unit can be freely into an algorithm to facilitate the calculation accuracy and efficiency of three methods of numerical simulation. Through the study found, FEM-SPH algorithm not only overcome the defect of smooth particle tensile instability, but also overcomes the problem of low efficiency of finite element computation. To further test the FEM-SPH algorithm has advantages in the practical engineering, we have carried out the actual test to the example of the super high speed collision, concluded that, since the target of most of the computational domain is always finite element, smoothed particle focused only in contact with the projectile and target of local area, particle number is not much, the whole calculation process just ten minutes, computational efficiency has been greatly improved, at the same time in the simulation of large deformation, the advantage is very obvious .This provides a criterion for the actual project, depending on the specific material deformation mode and choose a more appropriate conversion algorithm.

This paper presents an approach to the elastic analysis of beams subjected to Saint-Venant torsion using Green’s theorem and the finite difference method (FDM). The Saint-Venant torsion of beams, also called free torsion or unrestrained torsion, is characterized by the absence of axial stresses due to torsion; only shear stresses are developed. A solution to this torsion problem consists of finding a stress function that satisfies the governing equation and the boundary conditions. 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 nodes of the structure. In this paper, the beam’s cross-section was discretized using a two-dimensional grid and additional nodes were introduced on the boundaries. The introduction of additional nodes allowed us to apply the governing equations at boundary nodes and satisfy the boundary conditions. Beams with solid sections as well as multiply connected cross-sections were analyzed using this model; shear stresses and localized stresses at reentrant corners, torsion constant, and warping displacements were determined. Furthermore, beams with thin-walled closed sections, single-cell or multiple-cell, were analyzed using the Prandtl stress function whereby a linear distribution of the shear stresses over the thickness was considered; closed-form solutions for shear stresses and torsion constant were derived. The results obtained in this study showed good agreement with the exact results for rectangular cross-sections, and the accuracy was increased through a grid refinement. For thin-walled closed sections, the shear stresses obtained at the centerline using the closed-form solutions were in agreement with the values using Bredt’s analysis but the maximal values in the cross-section, which did not necessarily occur at the position with the smallest thickness, were higher; in addition, the results using the closed-form solutions were in good agreement with those using FDM.

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.

Although Eulerian approaches are standard in computational acoustics, they are less effective for certain classes of problems like bubble acoustics and combustion noise. A different approach for solving acoustic problems is to compute with individual particles following particle motion. In this paper, a Lagrangian approach to model sound propagation in moving fluid is presented and implemented numerically, using three meshfree methods to solve the Lagrangian acoustic perturbation equations (LAPE) in the time domain. The LAPE split the fluid dynamic equations into a set of hydrodynamic equations for the motion of fluid particles and perturbation equations for the acoustic quantities corresponding to each fluid particle. Then, three meshfree methods, the smoothed particle hydrodynamics (SPH) method, the corrective smoothed particle (CSP) method, and the generalized finite difference (GFD) method, are introduced to solve the LAPE and the linearized LAPE (LLAPE). The SPH and CSP methods are widely used meshfree methods, while the GFD method based on the Taylor series expansion can be easily extended to higher orders. Applications to modeling sound propagation in steady or unsteady fluids in motion are outlined, treating a number of different cases in one and two space dimensions. A comparison of the LAPE and the LLAPE using the three meshfree methods is also presented. The Lagrangian approach shows good agreement with exact solutions. The comparison indicates that the CSP and GFD method exhibit convergence in cases with different background flow. The GFD method is more accurate, while the CSP method can handle higher Courant numbers.

Polymer nanocomposites are composed of polymer materials reinforced with nano fillers. In the present study the effective thermal conductivity of the composites filled with nanofillers has been investigated using commercially available finite element software ABAQUS 6.11. The nanofillers used were alumina particles and multi-walled carbon nanotubes (MWNTs) and the matrix was considered to be made of epoxy. For the analysis 2D and 3D Representative Volume Elements (RVEs) were generated using Random Sequential Adsorption (RSA) algorithm using MATLAB and Python scripts. Thermal conductivity was found out for 2D and 3D RVEs for different area and weight fractions respectively. Two different shapes of the alumina nanoparticles were considered: spherical (circular) and ellipsoidal (elliptical) for 3D (2D) analysis. It was found that the thermal conductivity was increased with the addition of nanofillers. The increase in thermal conductivity was approximately same for both types of inclusions at corresponding area or weight fractions in 2D or 3D analysis. The results showed that addition of MWNTs to the composites lead to a significant increase in thermal conductivity than spherical or ellipsoidal inclusions.

This paper presents an approach to the vibration analysis of axially functionally graded non-prismatic Timoshenko beams (AFGNPTB) using the finite difference method (FDM). The characteristics (cross-sectional area, moment of inertia, elastic moduli, shear moduli, and mass density) of axially functionally graded beams vary along the longitudinal axis. The Timoshenko beam theory covers cases associated with small deflections based on shear deformation and rotary inertia considerations. 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 AFGNPTB 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 AFGNPTB, considering the damping. The results obtained in this study showed good agreement with those of other studies, and the accuracy was always increased through a grid refinement.

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.

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.

Combining modern thermodynamics theory branches, including finite time thermodynamics or entropy generation minimization, constructal theory and entransy theory, with metallurgical process engineering, this paper provides a new exploration on generalized thermodynamic optimization theory for iron and steel production processes. The theoretical core is to thermodynamically optimize performances of elemental packages, working procedure modules, functional subsystems, and whole process of iron and steel production processes with real finite-resource and/or finite-size constraints with various irreversibilities toward saving energy, decreasing consumption, reducing emission and increasing yield, and to achieve the comprehensive coordination among the material flow, energy flow and environment of the hierarchical process systems. A series of application cases of the theory are reviewed. It can provide a new angle of view for the iron and steel production processes from thermodynamics, and can also provide some guidelines for other process industries.

In this paper, a novel permanent magnet-assisted synchronous reluctance machine (PMASynRM) with rare-earth PMs and ferrite magnets is proposed. The performance of PMASynRM is discussed with respected to the different magnet ratio of rare-earth PMs and ferrite magnets. Some characteristics including the flux density, output torque, cogging torque, output power, power factor, torque ripple, loss, efficiency, and demagnetization are calculated by 2-D finite element analysis (FEA). The analysis results show that the excellent performance can be obtained by using hybrid magnet of rare-earth PMs and ferrite magnets with the suitable magnet ratio, and provide some desirable cost-performance trade-off.

The Wind Energy Conversion System (WECS) plays an inevitable role across the world. In particular, the attention for Permanent Magnet Synchronous Generators (PMSGs) connected with wind farm is popular. This paper deals with the literature review that describes the recent advances, progresses and innovatory trends on PMSGs for WECS. Comparison between geared and direct-driven conversion systems and the classification of electrical machines used in WECS are discussed. A detailed analysis on the design aspects considering various topologies of PMSGs are encompassed in the literature. The PMSG design and optimization problems are solved by field computation techniques and optimized by using Soft Computing (SC) techniques .The three-dimensional, finite element software platform for the analysis and design of PMSGs is discussed. This paper also deals with the interdisciplinary modeling, analysis, and optimization of PMSG using Finite Element Analysis (FEM) and SC techniques. Finally, PMSGs are reviewed and compared for further exploration.

The effects of CO₂-water-rock interactions on the mechanical properties of shale are essential for estimating the possibility of sequestrating CO₂ in shale reservoirs. In this study, uniaxial compressive strength (UCS) tests together with an acoustic emission (AE) system and SEM & EDS analysis were performed to investigate the mechanical properties and microstructural changes of black shales with different saturation times (10 days, 20 days and 30 days) in water dissoluted with sub-/super-critical CO₂. According to the experimental results, the values of UCS, Young’s modulus and brittleness index decrease gradually with increasing saturation time in water with sub-/super-critical CO₂. Compared to intact samples, 30-days’ saturation causes reductions of 56.43% in UCS and 54.21% in Young’s modulus for sub-critical saturated samples, and 66.05% in UCS and 56.32% in Young’s modulus for super-critical saturated samples, respectively. The brittleness index also decreases drastically from 84.3% for intact samples to 50.9% for samples saturated in water with sub-critical CO₂, to 47.9% for samples saturated in water with super-critical carbon dioxide (SC-CO₂). SC-CO₂ causes a greater reduction of shale’s mechanical properties. The crack propagation results obtained from the AE system show that longer saturation time produces higher peak cumulative AE energy. SEM images show that many pores occur when shale samples are saturated in water with sub-/super-critical CO₂. The EDS results show that CO₂-water-rock interactions increase the percentages of C and Fe and decrease the percentages of Al and K on the surface of saturated samples when compared to intact samples.

An effective method to inspect inaccessible nuclear power facility by interface wave which propagate along the shrink fit boundary of reactor head is proposed in this study. Reactor head is relatively thick to inspect from the outside of reactor by conventional ultrasonic testing. The proposed interface wave can propagate a long distance from the fixed transducer position. The inside of nuclear reactor is limited to access due to the high radiation, so transducers are located at outside of nuclear facility and interface wave propagates into the nuclear reactor for defect detection. The numerical simulation and experiments were carried out to validate the effectiveness of the proposed interface wave inspection method. Various defect cases simulating field failures are also presented with satisfactory detectability by the proposed technique with the features for defect classification.

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.

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.

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.

Due to the demands of personal travels and entertainments, light airplanes and small business aircrafts are developing rapidly. Light airplane structure is simple; however, it lacks crashworthiness design, especially the considerations on the impact of energy absorption. Therefore, in an event of accident, significant damage to passengers will be usually incurred. Airplanes made of composite materials structurally have high specific strength and good aerodynamic configuration. These materials have become the primary choice for new airplane development. This study mainly explores the topology optimization analysis of the light aircraft’s cockpit made of carbon fiber reinforced composites. This paper compares the compression amounts in the original models of composite material and aluminum alloy fuselages with the models after optimization during the crash-landing, in order to investigate the safety of fuselages made of different materials after structural optimization under the dynamic crashing. This study found that the energy absorbed by the aluminum alloy fuselage during crash-landing is still higher than that by the carbon fiber reinforced composites fuselage. On the other hand, the aluminum alloy fuselage after topology optimization could have an energy absorption capability enhanced by 40%, as compared to the that of the original model of aluminum alloy fuselage.

Structural deformations are one of the most significant factor that affects machine tool (MT) positioning accuracy. These induced errors are complex to be represented by a model, nevertheless they need to be evaluated and predicted in order to increase the machining performance. This paper presents a novel approach to calibrate a machine tool in real-time, analyzing the thermo-mechanical errors through Fibre Bragg Grating (FBG) sensors embedded in the MT frame. The proposed configuration consists of an adaptronic structure of passive materials, Carbon Fibre Reinforced Polymers (CFRP), equipped by FBG sensors that are able to measure in real-time the deformed conditions of the frame. By using a proper thermo-mechanical kinematic model, the displacement of the end effector may be predicted and corrected when it is subjected to external undesired factors. By starting from a set of FE simulations to develop a model able to describe the MT structure stresses, a prototype has been fabricated and tested. The scope was to compare the numerical model with the experimental tests using FBG sensors. The experimental campaign has been performed varying the structure temperature over time and measuring the tool tip point (TTP) positions. The obtained results showed a substantial matching between the real and the predicted position of TTP confirming the effectiveness of the proposed calibration system.

The aim of this work was to determine if the development of low-cost or no-cost finite element analysis (FEA) software has advanced to the point where it can be used in place of trusted commercial FEA packages for linear static structural analyses using isotropic material models. Nonlinear structural analysis will be covered in a separate paper. Several suitable packages were identified, these underwent a process of systematic elimination when they were unable to meet the minimum imposed qualitative criteria. Three packages were chosen to be subjected to performance benchmarking, namely: Code_Aster/Salome Meca; Mecway and Z88 Aurora. SimScale, a browser-based analysis package was included as well because it met all the baseline criteria and has the potential to offer a completely cloud-based approach to computer aided engineering, potentially reshaping the way an engineering business views its operational capabilities. This paper presents the test cases and simulation results for packages that fall under the linear static structural analysis type.

The modeling of flood-induced bend scour near embankment toes can provide important information for river engineering, embankment safety warnings, and emergency action management. During the rainy seasons, short-term general scour and bend scour are the most common causes for the failure of reinforced concrete embankments in Taiwan. To gain a deeper understanding of the scouring process near levee foundations, this study proposed a straightforward and practical method for bend scour simulation. The proposed simulation method is subdivided into three stages: two-dimensional flow simulation, general scour estimation, and bend scour estimation. A new bend-scour computation equation is proposed and incorporated into a two-dimensional hydraulic finite-volume model for simulating the evolution of bend scour depth around embankment toes. The proposed method is applied to simulate the temporal evolution of bend scouring near the Shuideliaw Embankment on the Cho-Shui River in Taiwan, where serious failure occurred during the June 2012 monsoon. Field data were gathered using the numbered-brick technique at the Shuideliaw Embankment to demonstrate the accuracy of the proposed method. The results of the bend scour simulations compared reasonably well with field measurements, indicating close agreement in terms of water levels and bend scour depths near the Shuideliaw Embankment. The proposed method was found to quickly estimate the maximum short-term general scour and bend scour depths for further enhancement of the safety of the embankment toe.

The production of lithium-ion battery cells is characterized by a high degree of complexit due to numerous cause-effect relationships between process characteristics. Knowledge about the multi-stage production is spread among several experts, rendering tasks such as failure analysis challenging. In this paper, a method is presented, which includes expert knowledge acquisition in production ramp-up by combining Failure Mode and Effects Analysis (FMEA) with a Bayesian Network. We show the effectiveness of this holistic method by building up a large scale, cross-process Bayesian Failure Network in lithium-ion battery production. Using this model, we are able to conduct root cause analyses as well as analyses of failure propagation. The former support operators in identifying root causes once a cell possesses a specific failure by calculating most-probable explanations matched to the individual battery cell data. The latter enable us to analyze propagation of failures and deviations in the production chain and thus provide support for placement of quality gates, leading to a significant reduction in scrap rate. Moreover, it gives an insight into which process steps are key drivers for which final product characteristics.

During HVDC earth return operation systems, a high magnitude current will be injected into soil through earth electrode, the potential on the surface would change widely and produce unfavorable effects on the AC systems around. This paper presents an effective finite element method (FEM) coupling electric field with thermal field to evaluate the electrical field induced by the injected DC current. Firstly, owe to the characteristic of FEM, this method can consider arbitrary soil and earth electrode structure. Secondly, by setting the electrical and thermal parameters of soil as a function of temperature at the same time, the dynamic coupling process of electric field and thermal field is simulated accurately. Thirdly, to deal with the singular point in FEM subdivision and the huge computation in traditional three-dimensional FEM, the FEM coupling 2-D earth electrode with 3-D soil based on "shell" theory is introduced. Finally, based on the suggested method, the effect of abnormal resistance region (ARR) near DC earth electrode on electric field distribution is analyzed.

Repair mortar is commonly used to rehabilitate reinforced concrete structures or components that exhibit a relatively high level of distresses. Yet, this repair mortar can be contaminated by salt from its service environment. This work employs a two-dimensional finite element model to investigate the transport behavior of ionic species in salt-contaminated and water-saturated repair mortar under an externally applied electric field. The model was experimentally validated and then utilized to evaluate the effectiveness of electrochemical chloride extraction (ECE) with or without electrical injection of corrosion inhibitor (EICI). In the case study, both the ECE alone and the ECE+EICI treatment was found effective in decontaminating the zone in front of the steel rebar. In both techniques, the magnitude of current density has a significant effect on removing chloride out of the mortar and increasing the pH of the pore solution near the rebar, whereas the treatment time any not have a significant effect under some scenarios. The injection of the organic corrosion inhibitor significantly slowed down the removal of chloride. Changes in the ionic distribution in the mortar were generally beneficial in reducing the corrosion risk of the steel rebar and thus extending the service life of the repair mortar.

We use Matlab 3D finite element fluid flow/transport modelling to simulate localized wellbore temperature events of order 0.05-0.1oC logged in Fennoscandia basement rock at ~ 1.5km depths. The temperature events are approximated as steady-state heat transport due to fluid draining from the crust into the wellbore via naturally occurring fracture-connectivity structures. Flow simulation is based on the empirics of spatially-correlated fracture-connectivity fluid flow widely attested by well-log, well-core, and well-production data. Matching model wellbore-centric radial temperature profiles to a 2D analytic expression for steady-state radial heat transport with Peclet number Pe ≡ r0φv0/D (r0 = wellbore radius, v0 = Darcy velocity at r0, φ = ambient porosity, D = rock-water thermal diffusivity), gives Pe ~ 10-15 for fracture-connectivity flow intersecting the well, and Pe ~ 0 for ambient crust. Darcy flow for model Pe ~ 10 at radius ~ 10 meters from the wellbore gives permeability estimate κ ~ 0.02Darcy for flow driven by differential fluid pressure between least principal crustal stress pore pressure and hydrostatic wellbore pressure. Model temperature event flow permeability κm ~ 0.02Darcy is related to well-core ambient permeability κ ~ 1µDarcy by empirical poroperm relation κm ~ κ exp(αmφ) for φ ~ 0.01 and αm ~ 1000. Our modelling of wellbore temperature events calibrates the concept of reactivating fossilized fracture-connectivity flow for EGS permeability stimulation of basement rock.

In the current study, the failure behavior of retrofitted steel structures was studied experimentally and theoretically with steel/CFRP double strap joints (DSJs) under quasi-static tensile loading. A series of DSJs with different bonding lengths are also considered and examined to experimentally assess the effective bond length. To predict the failure load values of the tested specimens, a new stress-based criterion, namely the point stress (PS) criterion is proposed. Although some theoretical predictive modelling for the strength between steel/CFRP joints under various loading conditions has been presented, in this work by using the new proposed approach, one can calculate rapidly and conveniently the failure loads of the steel/CFRP specimens. Furthermore, to assess the validity of the new proposed criterion, further experimental data on steel/CFRP DSJs available in the open literature are predicted using the PS criterion. Finally, it was found that a good agreement exists between the experimental results and the theoretical predictions based on the PS criterion.

A novel fiber optic sensing technology for high frequency dynamics detection is proposed in this paper, specifically tailored for structural health monitoring applications based on strain wave analysis, including both passive impact identification and active Lamb wave monitoring. The sensing solution relies on a fiber optic-based interferometric architecture associated to an innovative coherent detection scheme, which retrieves in a completely passive way the high-frequency phase information of the received optical signal. The optical fiber can be arranged into different configurations in order to meet the requirement of the specific application, the sensor sensitivity being maximized while still ensuring a limited gauge length if a local measure is required. For the active Lamb wave monitoring, this results in a sensing fiber arranged in multiple loops glued on an aluminum thin panel in order to increase the phase signal relative only to the sensing points of interest. Instead, for passive impact identification, the sensitivity is simply increased by exploiting a longer gauge length glued to the structure. The fiber optic coherent (FOC) sensor detects the strain waves emitted by a piezoelectric transducer placed on the aluminum panel or generated by an impulse hammer, respectively. The FOC sensor measurements have been compared with both a numerical model based on Finite Elements and traditional piezoelectric sensors, confirming a very good agreement between experimental and simulated results for both the active and passive impact monitoring scenarios.

Energy efficiency and leakage magnetic field (LMF) are two important issues in electric vehicle inductive chargers. In this work, the maximum achievable coil efficiency and the corresponding LMF strength are formulated as functions of hardware parameters, and figure of merits (FOM) are proposed for assessing the efficiency and LMF performance of the coil assembly pair. The impacts of the coil assemblies’ geometric parameters on both FOMs are examined with the aid of finite element analysis (FEA), and measures to improve the FOMs are extracted from FEA results. A coil assembly pair is manually optimized within given dimensional limits. Compared with the initial design, the optimized one achieves higher efficiency and lower LMF strength while consuming less copper. The performance improvement is verified by FEA results and experimental data measured on an 85 kHz electric vehicle inductive charger prototype. The key measures for coil assembly optimization are summarized.

We find that the degeneracies and many peculiarities of the DNA genetic code may be described thanks to two closely related (fivefold symmetric) finite groups. The first group has signature $G=\mathbb{Z}_5 \rtimes H$ where $H=\mathbb{Z}_2 . S_4\cong 2O$ is isomorphic to the binary octahedral group $2O$ and $S_4$ is the symmetric group on four letters/bases. The second group has signature $G=\mathbb{Z}_5 \rtimes GL(2,3)$ and points out a threefold symmetry of base pairings. For those groups, the representations for the $22$ conjugacy classes of $G$ are in one-to-one correspondence with the multiplets encoding the proteinogenic amino acids. Additionally, most of the $22$ characters of $G$ attached to those representations are informationally complete. The biological meaning of these coincidences is discussed.

To enable purposeful design and implementation of automated concrete technologies, precise assessment and prediction of the complex material flow at various stages of the process chain are necessary. This paper investigates the intermediate stage of the extrusion and deposition processes in extrusion-based 3D-concrete-printing, using a numerical model based on the Particle Finite Element Method (PFEM). In PFEM, due to the Lagrangian description of motion, remeshing algorithms and the alpha shape method are used to track the free surface during large deformation scenarios. The Bingham constitutive model was used for describing the rheological behaviour of fresh concrete. This model is validated by comparing the numerically predicted layer geometries with those obtained from laboratory 3D printing experiments. Extensive parametric studies were then conducted using the numerical simulation, delineating the influence of process and material parameters on the layer geometries, the dynamic surface forces generated under the extrusion nozzle and the inter-layer interactions.

The purpose of this work is to optimize angular extrusion for sever plastic deformation of tubular specimens. A finite element model (FEM) was built in Deform-3D^® for three extrusion dies and analyzed. The dies were ECAP (135°), TCAP with external groove of 90° (denoted as TCAP-e) and TCAP with internal groove of 90° (denoted as TCAP-i). The analysis for process parameters (such as coefficient of friction-μ, die angles-ѱ and φ, temperature - T, and radius ratio-R) showed that TCAP-i was the optimal die in processing Al6063 tubes based on strain and load distribution of the model. The TCAP-i die model was further optimized for different parameters namely die angles, coefficient of friction, deformation ratio and temperature. The results showed that at constant process temperature of 25 °C, the optimal TCAP-i has the following parameters: φ₂=80⁰, Ѱ₁=30⁰, Ѱ₂=80⁰, Ѱ₃=30⁰, μ =0.4 and R =1.5.

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

Long-range radio transmissions open new sensor application fields, in particular for environment monitoring. As an example, the {\sl LoRa} radio protocol enables to connect remote sensors at distance as long as ten kilometers in line-of-sight. However, the large area covered also bring several difficulties, such as the placement of sensing devices in regard to geography topology, or the variability of communication latency. Sensing the environment also carries constraints related to the interest of sensing points in relation with a physical phenomenon. Criteria for designs are thus evolving a lot from the existing methods, especially in complex terrains. This article describes simulation techniques based on geography analysis to compute long-range radio coverages and radio characteristics in these situations. As radio propagation is just a particular case of physical phenomena, it is shown how a unified approach also allows to characterize the behavior of potential physical risks. The case of heavy rainfall and flooding is investigated. Geography analysis is achieved using segmentation tools to produce cellular systems which are in turn translated into code for high-performance computations. The paper provides results from practical complex terrain experiments using LoRa, that confirm the accuracy of the simulation, scheduling characteristics for sample networks, and performance tables for simulations on middle range Graphics Processing Units (GPUs).

Parts produced by laser or electron-beam powder bed fusion (PBF) additive manufacturing are prone to residual stresses, deformations, and other defects linked to non-uniform temperature distribution during the manufacturing process. Several researchers have highlighted the important role scan sequence plays in achieving uniform temperature distribution in PBF. However, scan sequence continues to be determined offline based on trial-and-error or heuristics, which are neither optimal nor generalizable. To address these weaknesses, we have articulated a vision for an intelligent online scan sequence optimization approach to achieve uniform temperature distribution, hence reduced residual stresses and deformations, in PBF using physics-based and data-driven thermal models. This paper proposes SmartScan, our first attempt towards achieving our vision using a simplified physics-based thermal model. The conduction and convection dynamics of a single layer of the PBF process are modeled using the finite difference method and radial basis functions. Using the model, the next best feature (e.g., stripe or island) that minimizes a thermal uniformity metric is found using control theory. Simulations and experiments involving laser marking of a stainless steel plate are used to demonstrate the effectiveness of SmartScan in comparison to existing heuristic scan sequences for stripe and island scan patterns. In experiments, SmartScan yields up to 43% improvement in average thermal uniformity and 47% reduction in deformations (i.e., warpage) compared to existing heuristic approaches. It is also shown to be robust, and computationally efficient enough for online implementation.

Pure orbital blowout fractures occur within the confines of the internal orbital wall. Restoration of orbital form and volume is paramount to prevent functional and esthetic impairment. The anatomical peculiarity of the orbit has encouraged surgeons to develop implants with customized features to restore its architecture. This has resulted in worldwide clinical demand for patient-specific implants (PSIs) designed to fit precisely in the patient's unique anatomy. Fused filament fabrication (FFF) three-dimensional (3D) printing technology has enabled the fabrication of implant-grade polymers such as Polyetheretherketone (PEEK), paving the way for a more sophisticated generation of biomaterials. This study evaluates the FFF 3D printed PEEK orbital mesh customized implants with a metric considering the relevant design, biomechanical, and morphological parameters. The performance of the implants is studied as a function of varying thicknesses and porous design constructs through a finite element (FE) based computational model and a decision matrix based statistical approach. The maximum stress values achieved in our results predict the high durability of the implants, and the maximum deformation values were under one-tenth of a millimeter (mm) domain in all the implant profile configurations. The circular patterned implant (0.9 mm) had the best performance score. The study demonstrates that compounding multi-design computational analysis with 3D printing can be beneficial for the optimal restoration of the orbital floor.

Most of the pore-scale imaging and simulations of non-Newtonian fluid are based on the simplifying geometry of network modeling and overlook the fluid rheology and heat transfer. In the present paper, we developed a non-isothermal and non-Newtonian numerical model of the flow properties at pore-scale by direct simulation of the 3D micro-CT images using a Finite Volume Method (FVM). The numerical model is based on the resolution of the momentum and energy conservation equations. Owing to an adaptive meshing technique and appropriate boundary conditions, rock permeability and mobility are accurately computed. A temperature and concentration-dependent power-law viscosity model in line with the experimental measurement of the fluid rheology is adopted. The model is first applied at isothermal condition to 2 benchmark samples, namely Fontainebleau sandstone and Grosmont carbonate, and is found to be in good agreement with the Lattice Boltzmann method (LBM). Finally, at non-isothermal conditions, an effective mobility is introduced that enables to perform a numerical sensitivity study to fluid rheology, heat transfer, and operating conditions. While the mobility seems to evolve linearly with polymer concentration, the effect of the temperature seems negligible by comparison. However, a sharp contrast is found between carbonate and sandstone under the effect of a constant temperature gradient. Besides concerning the flow index and consistency factor, a good master curve is derived when normalizing the mobility for both the carbonate and the sandstone.

We develop a technique for realizing lithium niobate on insulator (LNOI) waveguides of a multi-centimeter-length with a propagation loss as low as 0.027 dB/cm. Our technique relies on patterning a chromium (Cr) thin film coated on the top surface of LNOI into a hard mask with a femtosecond laser followed by the chemo-mechanical polishing for structuring the LNOI into the waveguides. The surface roughness on the waveguides is determined to be 0.452 nm with an atomic force microscope (AFM). The approach is compatible with other surface patterning technologies such as optical and electron beam lithographies or laser direct writing, enabling high-throughput manufacturing of large-scale LNOI-based photonic integrated circuits.

Variations in the implant thread shape and occlusal load behavior may result in significant changes in the biological and mechanical properties of dental implants and surrounding bone tissue. Most previous studies consider a single implant thread design, an isotropic bone structure, and a static occlusal load. However, the effects of different thread designs, bone material properties, and loading conditions are important concerns in clinical practice. Accordingly, the present study performs Finite Element Analysis (FEA) simulations to investigate the static, quasi-static and dynamic response of the implant and implanted bone material under various thread designs and occlusal loading directions (buccal-lingual, mesiodistal and apical). The simulations focus specifically on the von Mises stress, displacement, shear stress, compressive stress and tensile stress within the implant and the surrounding bone. The results show that the thread design and occlusal loading rate have a significant effect on the stress distribution and deformation of the implant and bone structure during clinical applications. Overall, the results provide a useful insight into the design of enhanced dental implants for an improved load transfer efficiency and success rate.

Much remains to be understood for delayed ultraweak photon emission (UPE) from organism in association with oxidative burst following external perturbation, a phenomenon that has been experimented for over three decades. Delayed UPE often decays hyperbolically; yet, it is not uncommon to have delayed UPE that reveals first-order kinetic patterns characterized by single-exponential, double exponential, or multi-exponential changes. Some delayed UPEs also presented transient patterns that are characteristic of second-order responses. A soliton-based photon-storage model has addressed the hyperbolic decaying pattern of delayed UPE; however, there are questions outstanding regarding modeling other non-hyperbolic kinetics as well as the large range of temporal scales of delayed UPE that can vary from 8.5 microseconds to many hours. This work proposes an alternative, phenomenological model-framework for interpreting the various kinetic patterns of delayed UPE following stress. The delayed UPE is considered to be governed by two sequential phases: a stress-transfer phase that transforms the external stress to photo-genesis for emitting photons, and a photon-diffusion phase that transmits the photons emitted by the photo-genesis to the surface of organism for being detected as delayed UPE. Time-resolved photon diffusion analysis reveals that any delayed UPE in organism with a delay time >100ns cannot be addressed by the inherent temporal spread that realistic tissue scattering will cause. A slow stress-transfer phase is thus required to explain the delay time-scales of delayed UPE at a minimum of 8.5 µs as reported experimentally. The stress-transfer phase is hypothesized to carry the following types of responses: single or multiple 1st-order low-pass, single 2nd-order low-pass with various damping factors, and single 2nd-order band-pass with various damping factors. A single 1st–order low-pass response with a time-varying kinetic rate is also analyzed. The responses of these modeled pathways to bolus and step inputs demonstrate that a kinetic pattern other than the exact single-exponential one may have multiple causes.

Transient analysis of gas flow in pipeline was studied. Finite Element Method based on ANSYS FLOTRAN was used to account for changes in pressure, temperature and flow rate. Compressibility factor function of temperature and pressure was considered. For non-isothermal transient results, the pressure and the flowrate gave wave propagation as a result of slow transients created by demand condition at the outlet end of custody transfer. Results obtained were in agreement with the demand restrictions at the outlet end of custody transfer indicating that the predictions are accurate and reliable. The results demonstrated that the Finite Element Method gave accurate prediction of pressure, temperature and flowrate in transient gas studies. For steady state non-isothermal model, results showed that the magnitude of the average pressure drop was higher when pressure was predicted with constant compressibility factor, but the same average pressure drop was reduced when the pressure was predicted with variable compressibility factor, z. Since compressibility factor is a function of temperature and pressure, the above findings signifies that in the case when gas temperature does not stabilize, the prediction of pressure with isothermal model and constant compressibility factor will lead to significant errors.

A new approach is presented to defining the amount of information, in which information is understood as the data about a finite set as a whole, whereas the average length of an integrative code of elements serves as a measure of information. We obtain a formula for the syntropy of reflection, that is, the information that two intersecting finite sets reflect about each other. The specificities of reflection of discrete systems through the combination of their parts are considered, and it is shown that the additive syntropy and entropy of reflection are measures of structural order and chaos. Three information laws have been established: the law of conservation of the sum of chaos and order; information law of reflection; the law of conservation and transformation of information. An assessment of the structural organization and the level of development of discrete systems is presented. It is shown that various measures of information are structural characteristics of integrative codes of elements of discrete systems. The conclusion is made that, from the informational-genetic positions, the synergetic theory of information is primary in relation to the Hartley-Shannon information theory. In the appendix we consider the asymmetry of the mutual reflection of finite sets and the proportionality of arbitrary quantities.

In this study, we design a new refractive index sensor based on a metal–insulator–metal waveguide coupled with a notched ring and vertical rectangular resonators. We use the finite element method to study the propagation characteristics of the sensor. According to the calculation results, the transmission spectrum exhibits a typical Fano resonance shape. The phenomenon of Fano resonance is caused by the coupling between the wide spectrum and narrow spectrum. In the design, the wide spectrum signal is generated by the vertical rectangular resonator, while the narrow spectrum signal is generated by the notched ring resonator. In addition, we varied the structural parameters of the resonators and filled the structure with media of different refractive indices to study the sensing properties. The maximum achieved sensitivity of the sensor reached 1071.4 nm/RIU. The results reveal potential applications of the coupled system in the field of sensors.

A numerical application of linear-molecule symmetry properties, described by the D_∞h point group, is formulated in terms of lower-order symmetry groups D_nh with finite n. Character tables and irreducible representation transformation matrices are presented for D_nh groups with arbitrary n-values. These groups are subsequently used in the construction of symmetry-adapted ro-vibrational basis functions for solving the Schrödinger equations of linear molecules as part of the variational nuclear motion program TROVE. The TROVE symmetrisation procedure is based on a set of "reduced" vibrational eigenvalue problems with simplified Hamiltonians. The solutions of these eigenvalue problems have now been extended to include the classification of basis-set functions using ℓ, the eigenvalue (in units of ℏ) of the vibrational angular momentum operator ${\widehat{L}}_z$ . This facilitates the symmetry adaptation of the basis set functions in terms of the irreducible representations of D_nh. ¹²C₂H₂ is used as an example of a linear molecule of D_∞h point group symmetry to illustrate the symmetrisation procedure.

This paper presents a demonstration of the feasibility of fabricating micro-cantilever harvesters with extended stress distribution and enhanced bandwidth by exploiting an M-shaped two-degrees-of-freedom design. The measured mechanical response of the fabricated device displays the predicted dual resonance peak behavior with the fundamental peak at the intended frequency. This design has the features of high energy conversion efficiency in a miniaturized environment where the available vibrational energy varies in frequency. It makes such a design suitable for future large volume production of integrated self powered sensors nodes for the Internet-of-Things.

The meticulous study of finite automata has produced many important and useful results. Automata are simple yet efficient finite state machines that can be utilized in a plethora of situations. It comes, therefore, as no surprise that they have been used in classic game theory in order to model players and their actions. Game theory has recently been influenced by ideas from the field of quantum computation. As a result, quantum versions of classic games have already been introduced and studied. The PQ penny flip game is a famous quantum game introduced by Meyer in 1999. In this paper we investigate all possible finite games that can be played between the two players Q and Picard of the original PQ game. For this purpose we establish a rigorous connection between finite automata and the PQ game along with all its possible variations. Starting from the automaton that corresponds to the original game, we construct more elaborate automata for certain extensions of the game, before finally presenting a semiautomaton that captures the intrinsic behavior of all possible variants of the PQ game. What this means is that from the semiautomaton in question, by setting appropriate initial and accepting states, one can construct deterministic automata able to capture every possible finite game that can be played between the two players Q and Picard. Moreover, we introduce the new concepts of a winning automaton and complete automaton for either player.

Purpose: The aim of the present experimental 3D-finite element study was to evaluate the influence of an augmented sinus lift with an additional inserted bone graft block. The bone graft block stabilizes the implant in addition to conventional augmented bone. We placed the block in three different positions. The implants were loaded with axial force and forces secondary to laterotrusion and protrusion. Material and Methods: A simplified U-shaped 3D finite element model of the upper jaw and a more complex anatomical model of the left maxilla were created. The bone graft block was placed in three positions: in the lower third in contact with the sinus floor, the middle, and the upper third of the implant. Van Mises’ stress distribution was calculated and analyzed for the different models. We also compared the complex anatomical model with the simplified one. Results: The position of the bone graft block significantly influences the magnitude of stress distribution. A bone graft block positioned in the upper third or middle of the implant reduces the quantity of stress compared to the reference model without a bone graft block. The low bone graft block position is clearly associated with lower stress distribution in compact bone. We registered no significant differences in stress in compact bone with regard to laterotrusion or protrusion. Conclusions: Maximum values of von Mises stresses in compact bone can be reduced significantly by using a bone graft block. The reduction of stress is nearly the same for positions in the upper third and the middle of the implant. It is much more pronounced when the bone graft block is in the lower third of the implant near the sinus floor, which appeared to be the best position in the present study.

The fractional telegraph partial differential equation with fractional Atangana-Baleanu-Caputo (ABC) derivative is studied. Laplace method is used to find the exact solution of this equation. Stability inequalities are proved for the exact solution. Difference schemes for the implicit finite method are constructed. The implicit finite method is used to deal with modelling the fractional telegraph differential equation defined by Caputo fractional of Atangana-Baleanu (AB) derivative for different interval. Stability of difference schemes for this problem is proved by the matrix method. Numerical results with respect to the exact solution confirm the accuracy and effectiveness of the proposed method.

Lightning strike can cause a considerable damage in aircraft parts made from semiconducting materials such as Carbon Fiber Reinforced Plastics (CFRPs). Therefore, in recent years, the lightning strike phenomenon has attracted the interest of the academic community and the aircraft industry. Until now, the problem has been addressed mainly experimentally, while the reported numerical works are very limited. In the present work, a coupled electro-thermal FE model has been developed using the ANSYS commercial FE code to simulate the lightning strike damage in unidirectional CFRP laminates due to the Joule heat flux phenomenon. The model is based on the SOLID69 thermoelectric element and applies a non-linear, time-transient analysis. The main input to the model is the thermal-electrical properties of the composite material which vary with temperature. Using the model, a parametric study on the effect of mesh density and peak intensity on the thermal damage has been performed. Three electrical lightning strikes of low (10 kA), medium (30 kA) and high peak intensity (40 kA) have been applied according to the SAE ARP 5412 standard. The electro-thermal model has been validated against a numerical model from the literature. The numerical results reveal that the increase of peak intensity leads to the increase of the area and penetration depth of matrix thermal damage (pyrolysis) as well as to the increase of the area of fiber damage (deterioration and ablation). Through progressive damage modeling, the residual tensile strength of the CFRP plate after being subjected to lightning strike of different peak intensity has been predicted. Lightning strike initial damage has been simulated by translating the thermal field into degradation of elastic properties of the lamina. The results show an increase in the accumulated matrix damage and a decrease of tensile strength due to the initial lightning strike damage. For the maximum peak intensity of 40 kA, a decrease in tensile strength of 4.8% has been predicted

Quantitative dynamic strain measurements of the ground would be useful for engineering scale problems such as monitoring for natural hazards, soil-structure interaction studies and non-invasive site investigation using full waveform inversion (FWI). Distributed Acoustic Sensing (DAS), a promising technology for these purposes, needs to be better understood in terms of its directional sensitivity, spatial position, and amplitude for application to engineering-scale problems. This study investigates whether the physical measurements made using DAS are consistent with the theoretical reception patterns and experimental measurements of ground strain made by geophones. Results show that DAS and geophone measurements are consistent in both phase and amplitude for broadband (10s of Hz), high amplitude (10s of microstrain) and complex wavefields originating from different positions around the array when: (1) the DAS channels and geophone locations are properly aligned, (2) the DAS cable provides good deformation coupling to the internal optical fiber, (3) the cable is coupled to the ground through direct burial and compaction, and (4) laser phase noise is mitigated in the DAS measurements. The theoretical relationship between DAS-measured and point-wise strain for vertical and horizontal active sources is introduced using 3D elastic finite-difference simulations. The implications of using DAS strain measurements are discussed including directionality and magnitude differences between the actual and DAS-measured strain fields. A method for spatially aligning the DAS channels with the geophone locations at tolerances less than the spatial resolution of the DAS system is proposed.

This research aims to determine the similarities in groups of people to build a film recommender system for users. Users often have difficulty in finding suitable movies due to the increasing amount of movie information. The recommender system is very useful for helping customers choose a preferred movie with the existing features. In this study, the recommender system development is established by using several algorithms to obtain groupings, such as the K-Means algorithm, birch algorithm, mini-batch K-Means algorithm, mean-shift algorithm, affinity propagation algorithm, agglomerative clustering algorithm, and spectral clustering algorithm. We propose methods optimizing K so that each cluster may not significantly increase variance. We are limited to using groupings based on Genre and, Tags for movies. This research can discover better methods for evaluating clustering algorithms. To verify the quality of the recommender system, we adopted the mean square error (MSE), such as the Dunn Matrix and Cluster Validity Indices, and social network analysis (SNA), such as Degree Centrality, Closeness Centrality, and Betweenness Centrality. We also used Average Similarity, Computational Time, Association Rule with Apriori algorithm, and Clustering Performance Evaluation as evaluation measures to compare method performance of recommender systems using Silhouette Coefficient, Calinski-Harabaz Index, and Davies-Bouldin Index.

Evidence regarding the effect of the onlay preparation design for different CAD/CAM restorative materials considering the preservation of cusps is lacking. Molars were 3D modeled in four preparation designs for onlay restoration: traditional design with functional cusp coverage (TFC), non-retentive design with functional cusp coverage (NFC), traditional design with non-functional cusp coverage (TNFC), non-retentive design with non-functional cusp coverage (NNFC). Restorations were simulated with two CAD/CAM restorative materials: LD – lithium disilicate (IPS e.max CAD) and RC - resin composite (GrandioBloc). A 100 N axial load was applied to the occlusal surface simulating the centric contact point. Von Mises (VM) and maximum principal (Pmax) stresses were evaluated for restorations, cement layer and dental substrate. The non-retentive preparation design reduced the stress concentration in the tooth structure in comparison to the conventional retentive design. For LD onlays, the stress distribution on the restoration intaglio surface showed that the preparation design as well as the prepared cusp, influenced the stress magnitude. The non-retentive preparation design provided better load distribution in both restorative materials and more advantageous for tooth structure. The resin composite restoration on non-functional cusp is recommended when functional cusp is preserved, in order to associate conservative dentistry and low stress magnitude.

This paper deals with bistatic subsurface probing of a horizontally layered dielectric half-space by means of ultra-wideband electromagnetic pulses. A receiver collects reflections from the air-ground interface and from the gradients of dielectric permittivity in the half-space. This scenario is of interest for ground penetrating radar (GPR) applications. For the analytical description of the received signal, we developed and implemented a novel time-domain version of the coupled-wave Wentzel–Kramers–Brillouin approximation. Our solution is in very good agreement with finite-difference time-domain results, radically accelerates calculations, and effectively accounts for the protracted return signals observed in the lower part of the GPR spectrum. The paper includes results showing the application of the proposed technique to two case studies: in particular, the method was employed for the post-processing of experimental radargrams collected on Lake Chebarkul, in Russia, and to simulate GPR probing of the Moon surface, to detect smooth gradients of the dielectric permittivity in lunar regolith.

The Kummer surface was constructed in 1864. It corresponds to the desingularisation of the quotient of a 4-torus by 16 complex double points. Kummer surface is kwown to play a role in some models of quantum gravity. Following our recent model of the DNA genetic code based on the irreducible characters of the finite group G5:=(240,105)≅Z5⋊2O (with 2O the binary octahedral group), we now find that groups G6:=(288,69)≅Z6⋊2O and G7:=(336,118)≅Z7⋊2O can be used as models of the symmetries in hexamer and heptamer proteins playing a vital role for some biological functions. Groups G6 and G7 are found to involve the Kummer surface in the structure of their character table. An analogy between quantum gravity and DNA/RNA packings is suggested.

Linear elastic fracture modeling coupled with empirical material tension data result in good quantitative agreement with experimental measurements of fracture failure for both brittle and tough epoxy nanocomposites. The nanocomposites comprise diglycidyl ethers of bisphenol A cured with O,O’ bis (2-aminopropylpropylene glycol) (Jeffamine D230) and doped with rubber nanoparticles of varying concentrations. Toughness, critical load, and critical displacement in quasi-static single edge-notched three-point bending are predicted accurately using both surface-based cohesive zone (CZS) and extended finite element (XFEM) methods implemented in Abaqus software. Fracture initiation within a crack is taken at the yield stress from uniaxial tension data. Prediction of fracture processes using a generalized truncated linear traction-separation law was significantly improved by considering the case of a linear softening function. There are no adjustable parameters in the XFEM model. The CZS model requires only optimization of the element displacement at fracture parameter. Thus, these continuum methods describe these materials in mode I fracture with a minimum number of independent parameters.

Shot peening is an attractive technique for fatigue enhanced performance of metallic components, because it promotes crack initiation retardation and later crack growth. Engineering design based on fatigue crack propagation predictions applying the principles of fracture mechanics is commonly used in aluminum structures for aerospace engineering. The main purpose of present work was to analyze the effect of shot peening on the fatigue crack propagation of the 7475 aluminum alloy, under both constant amplitude loading and periodical overload blocks. The tests were performed on 4 and 8 mm thickness specimen's with stress ratios of 0.05 and 0.4. The analysis of the shot-peened surface showed a small increase of the micro-hardness values, due to the plastic deformations imposed by shot peening. The beneficial effect of surface peening on fatigue crack growth rates is quite limited to an increasing near the threshold. The specimens’ thickness has only marginal influence on the crack propagation, in opposite to the stress ratio. Periodic overload blocks of 300 cycles promotes a reduction of the fatigue crack growth rate for both intervals of 7,500 and 15,000 cycles.

Empirical equations representing, interpolating and smoothing groups of measurement results by regression methods are in widespread use in metrology and other fields of science and engineering. Standard equations for the propagation of measurement uncertainties to values computed from such equations are available but suffer from a lack of general acceptance and are only infrequently applied in practice. One reason for the slow uptake is the lack of clear methods to account for systematic error. In this paper, uncertainty propagation equations in terms of covariance matrices are generalized to allow for systematic errors in least-squares regression, and effects of the resulting uncertainties are investigated analytically. A stochastic ensemble model of systematic error is proposed for the computation of the non-diagonal elements of the weight matrix of the Generalized Least-Squares method (GLS) from the measurement uncertainty. A GLS projector formalism is described which separates the effects of measurement scatter on values calculated by the equation from those on the related “residuals”, i.e., the residual errors in the fitted data. The same projectors act similarly on the associated uncertainties and covariance matrices and permit the effects of systematic errors on the simulation covariance to be quantified. It is demonstrated that, in order to include systematic errors in the uncertainty estimates of GLS, covariance matrices may be substituted by novel, specifically defined dispersion matrices that are not specified yet by the GUM[1]. Systematic measurement errors may be estimated from various sources; a particularly easy way suggested here involves analyzing the structure in the regression residuals. It is demonstrated that GLS projectors exhibit inherent features of “error blindness” also with respect to systematic errors. [1] GUM: Guide to the Expression of Uncertainty in Measurement, http://www.bipm.org/en/publications/guides

This study presents a new form of velocity distribution in laminar liquid flow in rectangular microchannels using the eigenfunction expansion technique. Darcy friction factor and Poiseuille number are also obtained analytically. Due to the symmetry of the solutions, the effects of changing the aspect ratio from 0 to ∞ are also discussed. Using finite element method (FEM), the obtained analytical results are further compared with the 3D numerical simulations for the rectangular microchannels with different range of aspect ratio and pressure gradient, and excellent agreements were found. These findings provide additional insights in interpreting the results of the pressure-driven flows in finite aspect ratio microchannels, in which very precise comparison with the macroscale theory is crucial.

This paper presents a comprehensive and adaptive framework for orchid conservation programs illustrated with data from published and unpublished case studies. There is a specific focus on West Australian terrestrial orchids, but many of the approaches have universal relevance. Aspects of the framework include (1) setting appropriate objectives, (2) establishing effective collaborations between scientists, volunteers and regulators to fill knowledge and funding gaps, (3) use of survey and demographics data to determine extinction risks and management requirements for species, (4) effective habitat management to overcome threats such as grazing, (5) finding potential new habitats by modelling climate and site data, (6) investigating the effectiveness of pollinators and (7) using seed baiting to detect mycorrhizal fungi. The relative cost and effectiveness of different methods used to propagate orchids for translocation are compared. Methods known to be successful, in order of complexity, include placement of seed in situ, vegetative propagation, symbiotic germination in non-sterile organic matter, symbiotic germination in sterile culture, asymbiotic sterile germination and clonal division in tissue culture. These form a continuum of complexity, cost, time required, faculties needed, as well as the capacity to maintain genetic diversity and produce seedlings preadapted to survive in situ. They all start with seed collection and lead to seed storage, living collections used as tuber banks and seed orchards, as well as translocation for conservation. They could also lead to commercial availability and sustainable ecotourism, both of which are needed to reduce pressure on wild plants. Overall, there has been a strong preference to use relatively complex, expensive and time-consuming methods for orchid conservation, despite evidence that simpler approaches have also been successful. These simpler methods, which include in situ seed placement and non-sterile germination on inorganic substrates, should be trialled in combination with more complex orchid propagation methods as part of an adaptive management framework. It is essential that orchid conservation projects harness the unique biological features of orchids, such as abundant seed production and mycorrhizal fungi which are far more widespread than their hosts. This is necessary to increase the efficiency and coverage of recovery actions for the largest and most threatened plant family.

Slope failures pose a substantial threat to mining activity due to their destructive potential and high probability of occurrence on steep slopes close to limit equilibrium conditions, often found both in open pits and in waste and tailing disposal facilities. The development of slope monitoring and modeling programs usually entails the exploitation of in situ and remote sensing data together with the application of numerical modeling, and it plays an important role in the definition of prevention and mitigation measures aimed at minimizing the impact of slope failures in mining areas. Here we present a new methodology combining satellite radar interferometry and 2D finite element modeling for slope stability analysis at a regional scale, applied within slope unit polygons. We studied a former mining area in southeast Spain, and the method proved useful in detecting and characterizing a considerably large number of unstable slopes. Out of 1,959 slope units used for the spatial analysis of the radar interferometry data, 43 were unstable, with varying values of safety factor and landslide size. Out of the 43 active slope units, 21 exhibited line of sight velocities greater than the maximum error obtained through the validation analysis (2.5 cm/year). Eventually, this work discusses the possibility of using the results of the proposed approach to devise a proxy for landslide hazard. The proposed methodology can help to provide non-expert final users with intelligible, clear and easily comparable information to analyze slope instabilities in different settings, not limited to mining areas.

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.

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

Factorization reduces computational complexity and is therefore an important tool in statistical machine learning of high dimensional systems. Conventional molecular modeling, including molecular dynamics and Monte Carlo simulations of molecular systems, is a large research field based on approximate factorization of molecular interactions. Recently, the local distribution theory was proposed to factorize global joint distribution of a given molecular system into trainable local distributions. Belief propagation algorithms are a family of exact factorization algorithms for trees and are extended to approximate loopy belief propagation algorithms for graphs with loops. Despite the fact that factorization of probability distribution is their common foundation, computational research in molecular systems and machine learning studies utilizing belief propagation algorithms have been carried out independently with respective track of algorithm development. The connection and differences among these factorization algorithms are briefly presented in this perspective, with the hope to intrigue further development in factorization algorithms for physical modeling of complex molecular systems.

Every protein consists of a linear sequence over an alphabet of $20$ letters/amino acids. The sequence unfolds in the $3$-dimensional space through secondary (local foldings), tertiary (bonds) and quaternary (disjoint multiple) structures. The mere existence of the genetic code for the $20$ letters of the linear chain could be predicted with the (informationally complete) irreducible characters of the finite group $G_n:=\mathbb{Z}_n \rtimes 2O$ (with $n=5$ or $7$ and $2O$ the binary octahedral group) in our previous two papers. It turns out that some quaternary structures of protein complexes display $n$-fold symmetries. We propose an approach of secondary structures based on free group theory. Our results are compared to other approaches of predicting secondary structures of proteins in terms of $\alpha$ helices, $\beta$ sheets and coils, or more refined techniques. It is shown that the secondary structure of proteins shows similarities to the structure of some hyperbolic $3$-manifolds. The hyperbolic $3$-manifold of smallest volume --Gieseking manifold--, some other $3$ manifolds and the oriented hypercartographic group are singled out as tentative models of such secondary structures. For the quaternary structure, there are links to the Kummer surface.

Giant plasmonic surface enhancement has been observed in gold coated micron sized inverse pyramids entrapping a gold nanoparticle. The amplification of both surface enhanced Raman and photoluminescence signals was found to be dependent on the diameter of trapped gold nanoparticle and around 50-fold enhancement was detected for 250nm diameter sample relatively to the 50nm one. Finite differential time domain simulations, performed to determine the near-field distribution in the structure, showed that when the nanoparticle protrudes into the hotspot zone of the void, coupling of electromagnetic field occurs and the plasmon-related near-field enhancement is concentrated into the close vicinity of the nanoparticle, mainly into the close gaps around the tangential points of the curved sphere and the flat pyramid surface. This results in a more than 15 times increase of the near-field intensity, compared to the empty void.

A non-Newtonian fluid model is used to investigate the two-dimensional pulsatile blood flow through a tapered artery with a stenosis. The mixed convection effects of heat and mass transfer are also taken into account. By applying non-dimensionalization and radial coordinate transformation, we simplify the system in a tube. Under the finite difference scheme, numerical solutions are calculated for velocity, temperature concentration, resistance, impedance, wall shear stress and shearing stress at the stenosis throat. Finally, Quantitative analysis is carried out.

Fracture analysis of reinforced concrete deep beam strengthened with carbon fiber-reinforced polymer (CFRP) plates was carried out. The present research aims to find out whether the crack propagation in a strengthened deep beam follows linear elastic fracture mechanics (LEFM) theory or nonlinear fracture mechanics theory. To do so, a new energy release rate based on nonlinear fracture mechanics theory was formulated on the finite element method and the discrete cohesive zone model (DCZM) was developed in deep beams. To validate and compare with numerical models, three deep beams with rectangular cross-sections were tested. The code results based on nonlinear fracture mechanics models were compared with experimental results and ABAQUS results carried out based on LEFM. The predicted values of initial stiffness, yielding point and failure load, energy absorption, and compressive strain in the concrete obtained by the proposed model were very close to the experimental results. However, the ABAQUS software results have greater differences with the experimental results. For example, the predicted failure load for the shear-strengthened deep beam using the proposed model has only 6.3% differences compared to the experimental result. However, the predicted failure load using ABAQUS software based on LEFM has greater differences (25.1%) compared to the experimental result.

In this paper we survey some recent results due to the author concerning various inequalities and approximations for the finite Hilbert transform of a function belonging to several classes of functions, such as: Lipschitzian, monotonic, convex or with the derivative of bounded variation or absolutely continuous. More accurate estimates in the case that the higher order derivatives are absolutely continuous, are also provided. Some quadrature rules with error bounds are derived. They can be used in the numerical integration of the finite Hilbert transform and, due to the explicit form of the error bounds, enable the user to predict a priory the accuracy.

Background: Multifunctional nanoparticles are known for their wide range of biomedical applications. Controlling the magnetic properties of these nanoparticles is imperative for various applications, including therapeutic angiogenesis. The study was performed to evaluate the magnetic properties and their control mechanisms by the external magnetic field. Methods: A 100nm magnetic nanoparticle was placed in the magnetic field, and parametrically the magnet field strength and distance was evaluated. Various models of magnetic strength and disposition were evaluated. Magnetic flux density, force/weight, and magnetic gradient strength were the parameters evaluated in electromagnetic computational software. Results: The seven-coil method with three centrally placed coils as Halbach array, and each coil with a flux density of 7 Tesla, and with a coil dimension of 20cmx20cm (square model) of each coil showed a good magnetic strength and force/weight parameters in a distance of 15cm from the centrally placed coil. The particles were then evaluated for their motion characteristics in saline. It showed good displacement and acceleration properties. After that, the particles were theoretically assessed in a similar mathematical model after parametrically correcting the drag force. After application of high drag forces, the particles showed adequate motion characteristics. When the particle size was reduced further, the motion characteristics were preserved even with high drag forces. Conclusion: There is potential for a novel method of controlling multifunctional magnetic nanoparticles using high magnetic fields. Further studies are required to evaluate the motion characteristics of these particles in-vivo and invitro.

There is an increasing need for quantitative rockfall hazard and risk assessment that requires a precise definition of the terms and concepts used for this particular type of landslide. This paper suggests to use terms that appear to be the more logic and explicit as possible, and describes methods to derive some of the main hazard and risk descriptors. The terms and concepts presented concern the rockfall process (failure, propagation, fragmentation, modelling) and the hazard and risk descriptors, distinguishing the cases of localized hazards and diffused hazards. For a localized hazard, the failure probability of the considered rock compartment in a given period of time has to be assessed and the probability for a given element at risk to be impacted with a given energy must be derived combining the failure probability, the propagation probability and the exposure of the element. For a diffuse hazard that is characterized by a failure frequency, the number of rockfalls reaching the element at risk per unit of time and with a given energy (reach frequency) can be derived. However, when the element at risk is not replaced or repaired, the probability that it is impacted by at least one rockfall must be considered.

Expectation maximization (EM) algorithm is a powerful mathematical tool for estimating statistical parameter when data sample contains hidden part and observed part. EM is applied to learn finite mixture model in which the whole distribution of observed variable is average sum of partial distributions. Coverage ratio of every partial distribution is specified by the probability of hidden variable. An application of mixture model is soft clustering in which cluster is modeled by hidden variable whereas each data point can be assigned to more than one cluster and degree of such assignment is represented by the probability of hidden variable. However, such probability in traditional mixture model is simplified as a parameter, which can cause loss of valuable information. Therefore, in this research I propose a so-called conditional mixture model (CMM) in which the probability of hidden variable is modeled as a full probabilistic density function (PDF) that owns individual parameter. CMM aims to extend mixture model. I also propose an application of CMM which is called adaptive regressive model (ARM). Traditional regression model is effective when data sample is scattered equally. If data points are grouped into clusters, regression model tries to learn a unified regression function which goes through all data points. Obviously, such unified function is not effective to evaluate response variable based on grouped data points. The concept “adaptive” of ARM means that ARM solves the ineffectiveness problem by selecting the best cluster of data points firstly and then evaluating response variable within such best cluster. In order words, ARM reduces estimation space of regression model so as to gain high accuracy in calculation.

This study presents the development of a number of finite element (FE) models of the pelvis using different continuum and structural modelling approaches. Four FE models were developed using different modelling approaches: continuum isotropic, continuum orthotropic, hybrid isotropic and hybrid orthotropic. The models were subjected to an iterative adaptation process based on the Mechanostat principle. Each model was adapted to a number of common daily living activities (walking, stair ascent, stair descent, sit-to-stand and stand-to-sit) by applying onto it joint and muscle loads derived using a musculoskeletal modelling framework. The resulting models, along with a structural model previously developed by the authors, were compared visually in terms of bone architecture, and their response to a single load case was compared to a continuum FE model derived from CT imaging data. The main findings of this study were that the continuum orthotropic model was the closest to the CT derived model in terms of load response albeit having less total bone volume, suggesting that the role of material directionality in influencing the maximum orthotropic Young's modulus should be included in continuum bone adaptation models. In addition, the hybrid models, where trabecular and cortical bone were distinguished, had similar outcomes, suggesting that the approach to modelling trabecular bone is less influential when the cortex is modelled separately.

The possibility to test a radio frequency transceiver through the use of appropriate channel emulators allows to evaluate their performance under various operating conditions. Many systems are able to operate with a relatively limited instantaneous bandwidth by applying known statistical models. Sometimes it is necessary to evaluate the performance of a RF transceiver installed on high or very high speed platforms over predictable trajectories using optionally DEMs (Digital Elevation Models)data of the terrain to estimate the number of stimulated paths and their contributions during the flight. When the instantaneous bandwidth of the signal becomes high (over several hundred of MHz), one of the most important phenomenon to consider is the Doppler spread induced by the channel and traditional narrowband models become unuseful. In this paper there are presented some results when a time expansion is adopted in order to emulate the transceiver dynamic and the consequent Doppler spread with the aim to control the spectral purity of the emulated propagation channel.

Recently, the demand for electric vehicle is increasing worldwide due to eco-friendly policies and stricter emission regulations. As a traction motor for electric vehicle, interior permanent magnet synchronous motors are mainly used. For the design of the interior permanent magnet synchronous motor, the magnetic equivalent circuit method, which is a method of lumped constant circuit, and the finite element method, which is a method of distributed constant circuit, mainly are used. Magnetic equivalent circuit method is useful for simple design through fast and intuitive parameters, but it cannot derive the distribution of magnetic field. The finite element method can derive an accurate magnetic field distribution, but it takes a long time to analyze and it is difficult to analyze intuitive design parameters. In this paper, magnetic equivalent circuit method and Carter coefficient are mixed for rotor structure design. This design method will be called the hybrid magnetic equivalent circuit method. Intuitive design parameters are derived through this hybrid magnetic equivalent circuit method. We will derive the Air gap flux density distribution according to rotor shape, no-load induced voltage, and cogging torque, and compare and verify it with the finite element method.

Estimating the vehicle crashworthiness parameters experimentally is expensive and time consuming. For these reasons different modelling approaches are utilized to predict the vehicle behaviour and reduce the need for full-scale crash testing. The earlier numerical methods used for vehicle crashworthiness analysis were based on the use of lumped parameters models (LPM), a combination of masses and nonlinear springs interconnected in various configurations. Nowadays, the explicit nonlinear finite element analysis (FEA) is probably the most widely recognized modelling technique. Although informative, finite element models (FEM) of vehicle crash are expensive both in terms of man-hours put into assembling the model and related computational costs. A simpler analytical tool for early analysis of vehicle crashworthiness could greatly assist the modelling and save time. In this paper a simple piecewise LPM composed of a mass-spring-damper system, is used to estimate the vehicle crashworthiness parameters, focusing on the dynamic crush and the acceleration severity index (ASI). The model is first calibrated against a full-scale crash test and a FEM, post-processed with the LS-DYNA software, at an impact velocity of 56 km/h. The genetic algorithm is used to calibrate the model by estimating the piecewise lumped parameters (stiffness and damping of the front structure of the vehicle). After calibration, the LPM is applied to a range of velocities (40, 48, 64 and 72 km/h). The predictions for crashworthiness parameters from the LPM were compared with the predictions from the FEA and the results are much similar. It is shown that the LPM can assist in crash analysis, since LPM has some predictive capabilities and requires less computation time in comparison with the explicit nonlinear FEA.

Regeneration is usually regarded as a unique plant or some animal species process. In reality, regeneration is a ubiquitous process in all multicellular organisms. It ranges from response to wounding by healing the wounded tissue to whole body neoforming (remaking of the new body). In a larger context, regeneration is one facet of two reproduction schemes that dominate the evolution of life. Multicellular organisms can propagate their genes asexually or sexually. Here I present the view that the ability to regenerate tissue or whole-body regeneration is also determined by the sexual state of the multicellular organisms (from simple animals like hydra and planaria to plants and complex animals). The above idea is manifested here by showing evidence that many organisms, organs, or tissues show inhibited or diminished regeneration capacity when in reproductive status compared to the same organism organs or tissues in nonreproductive conditions or by exposure to sex hormones.

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.

One challenging problem is the representation of three-dimensional datasets that vary with time. These datasets can be though as a cloud of points that gradually deforms. But point-wise variations lack of information about the overall deformation pattern, and more importantly, about the extreme deformation locations inside the cloud. The present article applies a technique in computational mechanics to derive the strain-rate state of a time-dependent and three-dimensional data distribution, by which one can characterize its main trends of shift. Indeed, the tensorial analysis methodology is able to determine the global deformation rates in the entire dataset. With the use of this technique, one can characterize the significant fluctuations in a reduced multivariate description of an urban system and identify the possible causes of those changes: calculating the strain-rate state of a PCA-based multivariate description of an urban system, we are able to describe the clustering and divergence patterns between the districts of the city and to characterize the temporal rate in which those variations happen.

Network Simulators (NS) is typically used to study services and applications in complex scenarios due to the infeasibility of deploying real testbeds. This is the case of Mobile Ad hoc Networks (MANETs) which can be extensively used in emergency situations or to gather information without infraestructure. An important factor to get realiable results in the MANETs’ simulations is the wireless channel module of NSs. In this report, we explain how thw wireless channel is implemented in NS-3 and it can be configured.

Wave energy converter (WEC) arrays deployed in coastal regions may create physical disturbances potentially resulting in environmental stresses. Presently, limited information is available on the nature of these physical disturbance or the resultant effects. A quantitative Spatial Environmental Assessment Tool (SEAT) for evaluating potential effects of wave energy converter (WEC) arrays on nearshore hydrodynamics and sediment transport is presented for the central Oregon coast (USA) through coupled numerical model simulations of an array of WECs. Derived climatological wave conditions were used as inputs to the model to allow for the calculation of risk metrics associated with various hydrodynamic and sediment transport variables such as maximum shear stress, bottom velocity, and change in bed elevation. The risk maps provided simple, quantitative, and spatially-resolved means of evaluating physical changes in the vicinity of a hypothetical WEC array in response to varying wave conditions. Near-field risk of sediment mobility was determined to be moderate in the lee of the densely spaced array, where the potential for increased sediment deposition could result in benthic habitat alteration. Modifications to the nearshore sediment deposition and erosion patterns were observed near headlands and topographic features, which could have implications for littoral sediment transport. The results illustrate the benefits of a risk evaluation tool for facilitating coastal resource management at early market marine renewable energy sites.

In this paper, we investigated transducer placement strategies for detecting crack in primary aircraft structures using ultrasonic Structural Health Monitoring (SHM). The approach developed is for an expected damage location based on fracture mechanics, for example fatigue crack growth in a high stress location. To assess the performance of the developed approach, finite-element (FE) modelling of a damage tolerant aluminum fuselage has been performed by introducing an artificial crack at a rivet hole into the structural FE model and assessing its influence on the Lamb wave propagation, compared to a baseline measurement simulation. The efficient practical sensor position was determined from the largest change in area that is covered by reflected and missing wave scatter using an additive color model. Blob detection algorithms were employed in order to determine the boundaries of this area and to calculate the blob centroid. To demonstrate that the technique can be generalized, the results from different crack lengths and from tilted crack are also presented.

Detection of early stage corrosion on slender steel members is crucial for preventing buckling failures of steel structures. An active photoacoustic fiber optic sensors (FOS) system has been reported for early stage steel corrosion detection of steel plates and rebars using surface ultrasonic waves. The objective of this paper is to investigate the surface corrosion/rust detection problem on steel rods using numerically simulated surface ultrasonic waves. The finite element method (FEM) is applied in simulating the propagation of ultrasonic waves on steel rod models. Transmission mode of damage detection is adopted, in which one source (transmitter) and one sensor (receiver) are considered. In this research, radial displacements at the receiver were simulated and analyzed by short-time Fourier transform (STFT) for detecting, locating, and quantifying a surface rust located between the transmitter and the receiver. From our time domain and frequency domain analyses, it is found that the presence, location, and dimensions (length, width, and depth) of surface rust can be estimated by ultrasonic waves propagating through the surface rust.

Fatigue analysis plays a vital role in determining the structural integrity and life of a dental implant. With the use of such implants on the rise, there is a corresponding increase in the number of implant failures. As such, the aim of this research paper is to investigate the life of 3D-printed dental implants. The dental implants considered in this study were 3D printed according to the direct metal laser sintering (DMLS) method. Additionally, a finite element model was developed to study their performance, while fatigue life was predicted using Fe-Safe software®. The model was validated experimentally by performing fatigue tests. The life of the dental implants was analysed based on Normal strain and the Brown-Miller with Morrow mean correction factor algorithm. The model revealed that there was a strong correlation between the FEA and the experimental results. The clinical success of 3D-printed dental implant experimentally is 20.51 years and computationally under Normal strain is 19.89 years and Brown-Miller with Morrow mean correction factor is 26.82 years.

As long as the non-contact digital printing is not a common standard in the corrugated packaging industry, corrugated board crushing is a real issue that affects the load capacity of the boxes. Crushing mainly occurs during the converting of corrugated board (e.g. analog flexographic printing or laminating) and is a process that cannot be avoided. However, as show in this study, it can be controlled. In this work, extended laboratory tests were carried out on the crushing of double-walled corrugated board. The influence of fully controlled crushing (with a precision: ±10 μm) in the range from 10 to 70 % on different laboratory measurements was checked. Most of the typical mechanical tests were performed e.g. edge crush test, four-point bending test, shear stiffness test, torsional stiffness test, etc. on reference and crushed specimens. The residual thickness reduction of the crushed samples was also controlled. All empirical observations and performed measurements were the basis for building an analytical model of crushed corrugated board. The proven and verified model was then used to study the crushing effect of the selected corrugated board on the efficiency of simple packages with various dimensions.

This paper reveals the physical properties of Jeffery nanofluid flow past a moving plate embedded in porous medium under the existence of radiation and thermal diffusion. The analysis is carried out in three cases of moving plate, namely stationary plate λ = 0, forth-moving plate λ = 1, back-moving plate λ = −1. Finite difference method is applied to solve the governing equations of the flow and pointed out the variations in velocity, temperature and concentration with the use of graphical presentations. The impact of several parameters on local skin friction, Nusselt number and Sherwood number is also noticed and discussed. Enhancement of velocity is observed under the impact of Jeffery parameter for the cases of stationary plate and back-moving plate, whereas reverse nature is found in the case of forth-moving plate. The velocity enhances as the values of porosity parameter increases for the case of stationary plate and forth-moving plate but a reverse nature is noticed in the case of back-moving plate.

In this paper we discuss a problem of generalization of binomial distributed triangle, that is sequence A287326 in OEIS. The main property of A287326 that it returns a perfect cube n as sum of n-th row terms over k; 0 ≤ k ≤ n−1 or 1 ≤ k ≤ n, by means of its symmetry. In this paper we have derived a similar triangles in order to receive powers m = 5; 7 as row items sum and generalized obtained results in order to receive every odd-powered monomial n^2m+1; m ≥ 0 as sum of row terms of corresponding triangle. In other words, in this manuscript are found and discussed the polynomials D_m(n,k) and U_m(n,k), such that, when being summed up over k in some range with respect to m and n returns the monomial n^2m+1.

Applications in security and electronic surveillance require combination of excellent magnetic softness with good mechanical and anti-corrosive properties and low dimensionality. We overviewed the properties of soft magnetic glass-coated microwires and different post-processing making them quite attractive for electronic article surveillance and security applications. We studied the magnetostatic interaction between the microwires and routes to tune the non-linear hysteresis loops by using the arrays of different types of amorphous microwires. The presence of neighbouring microwire (either Fe or Co-based) significantly affects the hysteresis loop of the whole microwires array. In a microwires array containing magnetically bistable microwires we observed splitting of the initially rectangular hysteresis loop with a number of Barkhausen jumps correlated with the number of magnetically bistable microwires. Essentially, non-linear and irregular hysteresis loops have been observed in mixed arrays containing Fe and Co-rich microwires. The observed non-linearity in hysteresis loops allowed to increase the harmonics and tune their magnetic field dependencies. Non-linear hysteresis loops have been also observed upon devitrification of amorphous microwires. On the other hand, several routes allowing to tune the switching field by either post-processing or modifying the magnetoelastic anisotropy have been reviewed. Observed unique combination of magnetic properties together with thin dimensions and excellent mechanical and anti-corrosive properties provide excellent perspectives for the use of glass-coated microwires for security and electronic surveillance applications.

Propagation effects on the Narrow Bipolar Pulses (NBPs) or the radiation fields generated by compact cloud discharges as they propagate over finitely conducting ground are presented. The results are obtained using a sample of NBPs recorded with high time resolution from close thunderstorms in Sri Lanka. The results show that the peak amplitude and the temporal features such as the Full Width at Half Maximum (FWHM), zero crossing time and the time derivative of NBPs can be significantly distorted by propagation effects. For this reason the study of peak amplitudes and temporal features of NBPs and the remote sensing of current parameters of compact cloud discharges should be conducted using NBPs recorded under conditions where the propagation effects are minimal.

In this paper, we present a comprehensive 4D study of the early stage of plastic deformation in a polycrystalline binary AlLi alloy. The entire microstructure is mapped with X-ray diffraction contrast tomography and a set of bulk grains is further studied via X-ray topotomography during mechanical loading. The observed contrast is analyzed with respect to the slip system activation and the evolution of the orientation spread is measured as a function of applied strain. The experimental observations are augmented by the mechanical response predicted by crystal plasticity finite element simulations to analyze the onset of plasticity in details. Simulation results show a general agreement of the individual slip system activation during loading and that comparison with experiments at the length scale of the grains may be used to fine tune the constitutive model parameters.

This work deals with a numerical implementation of a fifth order CENTRAL WENO-TVD (\textit{Weighted Essentially Non-Oscillatory-Total Variation Dimimishing}) of Haschem (2006) scheme applied to the convective terms of some hyperbolic conservation laws problems, in a volume finite framework. The WENO-TVD scheme is used to solve the 1D advection and Burgers equations. For this case is implemented two different numerical fluxes: The Lax-Friedrichs and TVD fluxes. In the TVD fluxes the schemes applied are in flux-limiter form. The schemes implemented for this flux are: Van Albada-1 (van Albada et al.,1982), van Albada-2 (Kermani et al., 2003), van Leer (Hassanzadeh, 2009) and MINMOD (Hirsch, 2007). The WENO type schemes are characterized for their high order approximation, and do not produce spurious oscilations near discontinuities, shocks and higher gradients. A third order Runge-Kutta TVD for the temporal variable is used. Qualitative and quantitative comparison are presented. The numerical solutions are computed with an in-house computer code developed in MATLAB software. In future works, it will develope a paralelization of computer code for solving systems of conservation laws, e.g. Euler equations of gas dynamics.

Steel is one of the essential materials in the world's civilization. It is essential to produce many products such as pipelines, mechanical elements in machines, vehicles, profiles, and beam sections for buildings in many industries. Until the '50s of the 20th century, steel products required a complex process known as ingot casting; for years, steelmakers focused on developing and simplifying this process. The result was the con-tinuous casting process (CCP); it is the most productive method to produce steel. The CCP allows producing significant volumes of steel sections without interruption and is more productive than the formal ingot casting process. The CCP begins by transferring the liquid steel from the steel-ladle to a tundish. This tundish or vessel distributes the liquid steel, by flowing through its volume, to one or more strands having wa-ter-cooled copper molds. The mold is the primary cooling system, PCS, solidifying a steel shell to withstand a liquid core and its friction forces with the mold wall. Further down the mold, the rolls drive the steel section in the SCS. Here the steel section is cooled, solidifying the remaining liquid core, by sprays placed in every cooling segment all around the billet and along the curved section of the machine. Finally, the steel strand goes towards a horizontal-straight free-spray zone, losing heat by radiation mechanism, where the billet cools down further to total solidification. A moving torch cutting-scissor splits the billet to the desired length at the end of this heat-radiant zone.

Globally 30% of landslides occur in the northeastern part of India [1]. One of the major earthquake events in Sikkim, India occurred on 18th September 2011 (Mw 6.9) led to over 300 landslides and 122 human deaths [2]. These landslides not only controlled by natural disasters but initiated due to human activities. The present study considered Lungchok landslide occurred in south district of Sikkim due to 2011 seismic event. The study focused on the failure mechanism of the landslide based on finite element analysis by adopting eight different cases. The deformation characteristic was investigated for dry and saturated slope conditions under static and dynamic behavior considering vehicle loads using GeoStudio software. The FEM analysis has been carried out using load deformation and linear elastic. The analysis shows that the failure of the slope was not sudden due to the 2011 earthquake event, but progressive failure was observed with time and construction activity. The paper demonstrates that, an increase in infrastructure development including construction by hill cutting increased the initiation of landslide with soil erosion. The cracks developed after 2011 earthquake event led to further deformations during future disasters required effective stabilization measures.

Compressive sensing originates in the field of signal processing and has recently become a topic of energy-efficient data gathering in wireless sensor networks. In this paper, we propose an energy efficient distributed compressive sensing solution for sensor networks. Proposed solution utilizes sparsity distribution of signals to group sensor nodes into several coalitions and then implements localized compressive sensing inside coalitions. This solution improves data-gathering performance in terms of both data accuracy and energy consumption. The approach curbs both data-transmission costs and number of measurements. Coalition-based data gathering cuts transmission cost, and the number of measurements is reduced by scheduling sensor nodes and adjusting their sampling frequency. Our simulation showed that our approach enhances network performance by minimizing energy cost and improving data accuracy.

The effect of the different dental implants positioning region on the stress performance of the implant-supported prosthesis is not yet clear. This study evaluated the dental treatment with six dental implants in three different models and three different occlusal loading conditions, in terms of the biomechanical response of implants, prosthetic screw and maxilla, using three-dimensional finite element analysis. The finite element models were modelled containing external hexagon implants, as well as a Cobalt-Chromium superstructure. Three types of loads were applied: in the area of ​​the central incisors, first premolar and in the second molars. For the finite element simulations, the von-Mises stress peaks in the implant and in the surrounding cortical bone were analyzed. All recorded results reported higher values ​​for the implant-supported prosthesis in group C compared to the groups A and B. The highest stress values, ​​regardless the evaluated model, was in the prosthesis in group C and in screws, the smallest were in group A.

Long-term durability of the vacuum edge-seal plays a significant part in retrofitting triple vacuum glazing (TVG) to existing buildings in achieving towards zero-energy buildings (ZEB) target. Vacuum pressure decrement with respect to time between panes affect the thermal efficiency of TVG. This study reports a 3D finite element model, with validated mathematical methods and comparison, for the assessment of the influence of vacuum pressure diminution on the thermal transmittance (U value) of TVG. The centre-of-pane and total U values of TVG calculated to be 0.28 Wm⁻²K⁻¹ and 0.94 Wm⁻²K⁻¹ at the cavity vacuum pressure of 0.001 Pa. The results suggests that a rise in cavity pressure from 0.001 Pa to 100 kPa increases the centre-of-pane and total U values from 0.28 Wm⁻²K⁻¹ and 0.94 Wm⁻²K⁻¹ to 2.4 Wm⁻²K⁻¹ and 2.58 Wm⁻²K⁻¹, respectively. The temperature descent on the surfaces of TVG between hot and cold sides’ increases by decreasing the cavity vacuum pressure from 50 kPa to 0.001 Pa. To maintain the cavity vacuum pressure of 0.001 Pa for over 20 years of life span in the cavity of 10 mm wide edge sealed triple vacuum glazing, non-evaporable getters will maintain the cavity vacuum pressure that will enable the long-term durability to TVG.

Germplasm is a valuable natural resource in plant diversity that is crucial for its potential use. It provides knowledge about a species genetic composition. Germplasm protection strategies are not just planting hope threatened with extinction, they preserve medicinal and other essential plants on which survival rests. The successful use of genetic plant resources necessitates diligent collection, storage, analysis, documentation, and germplasm exchange. Slow growth cultures, cryopreservation, pollen and DNA banks, botanic gardens, genetic reserves and farmer’s fields are few conservation techniques. However, usage of an in vitro procedure with any chance of genetic instability leads to the destruction of the entire substance. Improved understanding of basic regeneration biology would, in turn, undoubtedly increase the capacity to regenerate plants from in vitro harvested explants, thus expanding selection possibilities. Germplasm conservation seeks to conserve endangered and vulnerable plant species worldwide for future proliferation and development; it is also the bedrock of agricultural production.

This study provides a contribution to the research field of 3D printed earthen buildings, focusing, for the first time, on the load-bearing capacity of these structures. The study involves the entire production and testing process of the earthen elements, from design, to the preparation of the mixture and the 3D printing, up to the uniaxial compression test on a wall segment. The results indicate that 3D printed earthen elements have a compressive strength of 2.32 MPa, comparable to that of rammed earth structures. The experimental data also made it possible to draw conclusions on the action of the infill, which seems to have the function of stopping the propagation of cracks. This has a positive effect on the overall behavior of 3D printed earthen elements, since it avoids the onset of dilative behavior in the final stages of the load test and maintains ultimate load values higher than 50% of the maximum load.

In this work, we have studied the propagation of ultrasonic waves of lysozyme solutions characterized by different degrees of aggregation and networking. The experimental investigation has been performed by means of the Transient Grating (TG) spectroscopy as a function of temperature; this technique enables to measure the ultrasonic acoustic proprieties over a wide time window, ranging from nanoseconds to milliseconds. The fitting of the measured TG signal allows the extraction of several dynamic properties, here we focused on the speed and the damping rate of sound. The temperature variation induces in the lysozyme solutions a series of processes: protein folding-unfolding, aggregation and sol-gel transition. Our TG investigation shows how these self-assembling phenomena modulate the sound propagation, affecting both the velocity and the damping rate of the ultrasonic waves. In particular, the damping of ultrasonic acoustic waves proves to be a dynamic property very sensitive to the protein conformational rearrangements and aggregation processes.