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.

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.

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.

To investigate the impact of a magnetic field on plaque development in a stenotic bifurcated artery, a finite element method is utilized. The blood flow is modelled as a stable, incompressible, Newtonian, biomagnetic, and laminar fluid. Furthermore, the arterial wall is assumed to be linear elastic. The Arbitrary Lagrangian Eulerian (ALE) method is employed to describe the hemodynamic flow in a bifurcated artery under the influence of an asymmetric magnetic field, taking into account two-way fluid-structure interaction coupling. A stable $P_2P_1$ finite element pair discretizes a nonlinear system of partial differential equations that requires a solution. The Newton-Raphson method is utilized to find a solution to the resulting nonlinear algebraic equation system. Numerical modelling is used to simulate the presence of magnetic fields, and the resulting displacement, velocity magnitude, pressure, and wall shear stresses are shown for a range of Reynolds numbers ($Re = 500$, $1000$, $1500$, and $2000$). The results of the numerical analysis demonstrate that the presence of a magnetic field has a significant effect not only on the magnitude of displacement but also on the velocity of the flow. The application of a magnetic field reduces flow separation, extends the recirculation area near the stenosis, and increases wall shear stress.

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.

The performance of a small reconfigurable unmanned aerial vehicle (UAV) is evaluated, 1 combining a multidisciplinary approach in computational analysis of additive manufactured struc- 2 tures, fluid dynamics, and experiments. Reconfigurable UAVs promise cost savings and efficiency 3 without sacrificing performance, while demonstrating versatility to fulfill different military mission 4 profiles. The use of computational fluid dynamics (CFD) in UAV design produces higher accuracy 5 aerodynamic data, which is particularly important for complex aircraft concepts such as blended 6 wing bodies. To address challenges relating to anisotropic materials, the Tsai-Wu failure criterion is 7 applied to structural analysis, using CFD solutions as load inputs. Aerodynamic performance results 8 show the low-speed variant attains an endurance of 1 hour, 48 minutes, whereas for its high-speed 9 counterpart, it is 29 minutes at 66.7% higher cruise speed. Each variant serves different aspects of 10 small UAS deployment in combat, with low-speed envisioned for close surveillance, and high-speed 11 for incursions. Experimental and simulation results suggest room for design iteration, in wing area 12 and geometry adjustments. Structural simulations demonstrated the need for airframe improvements 13 in the low-speed configuration. This paper highlights the potential of reconfigurable UAVs to disrupt 14 the industry, advocating for further research and design improvements.

Deformation behavior of MXene based polymer composites with bioinspired brick and mortar structures is analyzed. MXene/Polymer nanocomposites are modeled at microscale using bioinspired configurations of nacre-mimetic brick-and-mortar assembly structure. MXenes (brick) with polymer matrix (mortar) are modeled using classical analytical methods and numerical methods based on Finite Elements (FE). The analytical methods provide less accurate estimation of elastic properties compared to numerical one. MXene nanocomposite models analyzed with FE method provide estimates of elastic constants in the same order of magnitude as literature reported experimental results with good consistency. Bioinspired design of MXene nanocomposites results in the effective Young’s modulus of the nanocomposite increase by 25.1 % and the strength (maximum stress capacity within elastic limits) increase by 42.3 %. The brick and mortar structure of the nanocomposites leads to interlocking mechanism between MXene fillers in polymer matrix, resulting in effective load transfer, good strength, and damage resistance. This is demonstrated in this paper by numerical analysis of MXene nanocomposites subjected to quasi-static loads.

Variable Angle Tow (VAT) laminates offer a promising alternative to classical straight fiber composites in terms of design and performance. However, analyzing these structures can be more complex due to the introduction of new design variables. Carrera's Unified Formulation (CUF) has been successful in previous works for buckling, vibrational, and stress analysis of VAT plates. Typically, one-dimensional (1D) and two-dimensional (2D) CUF models are used, with a linear law describing the fiber orientation variation in the main plane of the structure. The objective of this article is to expand the CUF 2D plate finite elements family to perform free vibration analysis of composite laminated plate structures with curvilinear fibers. The primary contribution is the application of Reissner's Mixed Variational Theorem (RMVT) to a CUF finite element model. The Principle of Virtual Displacements (PVD) and RMVT are both used as variational statements for the study of monolayer and multilayer VAT plates dynamic behavior. The proposed approach is compared to Abaqus three-dimensional (3D) reference solutions, classical theories and literature results to investigate the effectiveness of the developed models. Results demonstrate that mixed theories provide the best approximation of the reference solution in all cases.

As the demand for liquefied petroleum gas (LPG) continues to increase worldwide, more LPG facilities are being constructed. However, the design and manufacture of pressure vessels used to store dangerous liquids or pressurised fluids, such as LPG, has been a major cause of hazards, including explosions and leakage. In this study, we aim to address this issue by designing and analysing a vertical pressure vessel with a capacity of 10 m3 of pressurised LPG in accordance with the American Society of Mechanical Engineers (ASME) code. Safety is the primary concern in designing pressure vessels due to the potential risks posed by accidents. Therefore, the main objective of this project is to design a pressure vessel that is safe from failure. The vessel is cylindrical with two elliptical heads, two nozzles, a manway and four leg supports, and we used Autodesk Inventor Professional 2023 for geometric modelling and Inventor Nastran for FEA analysis. To investigate the displacements, deflections and Von-Mises stresses in the pressure vessel, we used the finite element method, and Autodesk Inventor Nastran was employed for the analysis. Our results showed that changes to the structure of the pressure vessel were needed to reduce stress in the structure. Specifically, the displacement showed an inverse relationship with the tank section shell thickness. Moreover, we observed that the factor of safety increased linearly with shell thickness. We carefully considered permissible pressures and determined the required wall thickness to ensure acceptable maximum stresses. Our results indicate that the design is safe from failure, and the highest stresses are experienced by the manway followed by the shell, while the heads, nozzles and legs support experienced the lowest stresses. We also performed a theoretical calculation for the entire model and checked the results to ensure they were within the acceptable limits. In summary, our study highlights the importance of designing pressure vessels in accordance with ASME codes to ensure safety and prevent hazards caused by improper design and manufacturing. By using the finite element method, we identified potential stress points in the pressure vessel and were able to make modifications to improve its safety.

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.

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

This paper proposes a new structure of the dual-rotor hybrid-excited axial-flux permanent magnet vernier machine (DR-HEAFPMVM) with the modular stator and the consequent-pole PM (CPM) rotor for low-speed, high torque density applications such as in-wheel electric vehicles. The tooth-wound non-overlapping armature windings and direct current (DC) excitation windings are respectively arranged in stator main-teeth and split-teeth to obtain the modulated and adjustable air-gap flux densities, which result in the high torque density and outstanding flux-weakening capability. First, the design consideration, operation principle, and the air-gap flux density distributions of the proposed machine are elaborated based on the air-gap permeance function. Then, the influence of the pole ratios (PRs) and the DC excitation currents on the main electromagnetic performances of the DR-HEAFPMVM, such as the flux-weakening capability, back-electromotive force (back-EMF), on-load electromagnetic torque, loss distribution, and efficiencies, is investigated using the 3-D finite-element method (FEM). Results verify the feasibility of the flux adjustment of the DC excitation windings equipped in the split-tooth, and the design with a pole ratio of 8/1 tends to have higher torque density, higher machine efficiency, and considerable flux-weakening capability compared with the other two PRs.

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.

Wireless Sensor networks (WSNs) are gaining more and more research interest due to their ability to monitor large areas more or less independent. However, their reliability is an important issue, since their system reliability is influenced from a hardware, a data and an energetic point of view due to different factors like loading conditions, signal attenuation and battery lifetime. Regarding these issues, the positions of sensor nodes have to be chosen correctly to maximize the system reliability, which has to be considered during the development of new products equipped with WSNs. In this paper, an approach to estimate a WSNs system reliability during the development phase based on the analysis of the measurements on the example of strain measurements in finite elements (FE) models is presented. For that, the part under consideration is divided into regions with similar strains by the use of a Region Growing Algorithm (RGA). Afterwards, the WSN’s configuration is analyzed regarding its reliability in accordance to the data paths and measurement redundancy due to the sensor positions in the found measuring regions. This methodology is tested on an exemplary WSN configuration at an aircraft wing box under bending load and was found to work very well regarding the estimation of the hardware perspective on the system reliability. Therefore, the methodology and the developed algorithm show potentials for optimizations in sensor node positions for better reliability results.

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.

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.

Brushless Direct Current Motors are significantly improved for using many electrical applications. The design investigation is increased year by year and most of the motor design researchers find out the superior efficiency of the BLDC motor. The traditional motor consumed more electrical power than BLDC motor. The design of BLDC motors that are compact, small, easy to control, save energy and light weight. Due to this fact, the design of BLDC motor should be analyzed to improve the motor characteristics performance. The research focuses on the design of BLDC motor performance using Taguchi Method and Response Surface Methodology (RSM) for Finite Element Method (FEM). Firstly, the Taguchi Method is applied to identify the optimization point of design parameters and consider using the reference model of actual BLDC electric fan motor with a rated power of 26Watts design. Based on the Taguchi Method to create the RSM is used to two-step optimization design performance using in FEM analysis. The performances of motor design proposed approach to provide the effeteness for motor design investigator. In this investigation aims the better structural performances in output power, torque, efficiency and torque ripple. The final design of proposed model achieved an energy efficiency is 10% higher than the reference design of prototype modes, the output power result is 10W higher this reference design. We can reduce the loss of the proposed motor design. Furthermore, the maximum torque can also increase the 0.032Nm of that design. The design simulation results were produced and studied by the FEM analysis of JMAG software. From the evaluation results, the configurations of reference model design are compared with the proposed design by using the three methods of motor design optimization showed the best characteristics for performance of motor investigation design.

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

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.

The magnetic properties of 0.10 mm ultra-thin grain-oriented (UTGO) silicon steel and Fe-based amorphous (FBA) alloy under sinusoidal excitation were measured experimentally, and the magnetic field strength and iron loss of the two materials under different frequencies and magnetic densities were obtained. Based on the measured data, the magnetization and loss characteristics of the two materials were analyzed and compared. Furthermore, two Epstein square ring models of the same specification and different materials were designed, and the reliability of the models was verified. Then, the electromagnetic characteristics of the two Epstein square ring models at higher and lower frequencies were calculated by using the finite element method, and the iron loss was obtained and compared. The results show that the FBA alloy has good application characteristics at low frequency and low power, and the 0.10 mm UTGO silicon steel has good application characteristics at high frequency and high power. This research can provide important data basis to promote the application of these two materials in new energy equipment.

Eddy current based approaches have been investigated for a wide range of inspection applications. Dodd-Deeds model and the truncated region eigenfunction expansion (TREE) method are widely applied in various occasions, mostly for the cases that the sample is relatively larger than the radius of the sensor coil. The TREE method converts the integral expressions to the summation of many terms in the truncated region. In a recent work, the impedance of the co-axial air-cored sensor due to a plate of finite radius was calculated by the modified Dodd-Deeds analytical approach proposed by authors. In this paper, combining the modified analytical solution and the TREE method, a new finite region eigenfunction expansion (FREE) method is proposed. This method involves modifying its initial summation point from the first zero of the Bessel function to a value related to the radius of the plate, therefore makes it suitable for plate with finite dimensions. Experiments and simulations have been carried out and compared for the verification of the proposed method. Further, the planar size measurements of the metallic circular plate can be achieved by utilising the measured peak frequency feature.

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

Computing the rotation tensor is vital in the analysis of deformable bodies. This paper describes an explicit expression for the SO(3) rotation tensor R of the deformation gradient F, and successfully establishes an intrinsic relation between the exponential mapping Q = exp A and the deformation F. As an application, Truesdell's simple shear deformation is revisited.

Abstract—The static and fatigue analysis of plates with holes, made of 2D woven roving composites, is conducted. The parametrization of convex holes is proposed. Experimental results of specimens without holes and with different shapes of notches are discussed. The fatigue behaviour is considered with the use of the low cycle fatigue method. The analysis is supplemented by the numerical finite element modeling. The present work is an extension of the results discussed in the literature. The finite element modeling of plates with holes is presented in details in Ref [1-3] including also the problems of the accuracy. The experimental results and the form of specimens are shown in Ref [4]. Two parametrical description of static and fatigue behaviour of plates with holes characterizes the goal and the novelty of the paper. The fuzzy approach is employed to reduce the number of experimental data [5].

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.

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.

In the present article we describe a few simple and efficient finite volume type schemes on moving grids in one spatial dimension for the sake of simplicity. The underlying finite volume scheme is conservative and it is accurate up to the second order in space. The main novelty consists in the motion of the grid. This new dynamic aspect can be used to resolve better the areas with high solution gradients or any other special features. No interpolation procedure is employed, thus an unnecessary solution smearing is avoided. Thus, our method enjoys excellent conservation properties. The resulting grid is completely redistributed according the choice of the so-called monitor function. Several more or less universal choices of the monitor function are provided. Finally, the performance of the proposed algorithm is illustrated on several examples stemming from the simple linear advection to the simulation of complex shallow water waves. We believe that the techniques described in our paper can be beneficially used to model tsunami wave propagation and run-up.

This study investigates the magnetization of Cobalt (Co) nanoparticles with two different crystal structures: a centered cubic (fcc) and an hexagonal close packed (hcp), under the magnetic field of the Earth. For that reason, the hysteresis loops of the two material samples are exported and compared with the corresponding value by experimental procedure in the laboratory under the Earth’s magnetization. For both samples, the nanoparticle number variation and the nearest neighbor distance effect, in the Co hysteresis loop are investigated, structural information that is useful for the sample preparation. For the simulations the Finite Element Method (FEM) was used, while experiments are conducted using a measurement set-up based on fluxgate magnetometers. From the simulation, it is proven that, for both samples, the nanoparticle number effect does not change the hysteresis loop. In contrast, when the distance between the nanoparticles increases, the magnetization obtains smaller value. This study aims to build a novel, simple and economical measurement arrangement enabling in situ magnetic measurements for ferromagnetic nano-material characterization. As the need for continuous knowledge of the condition and the phase of ferromagnetic nanoparticles that is crucial in many applications, like catalytic processes, this work aims to introduce an alternative approach that will add another dimension to conventional chemical reaction engineering this of the magnetic state of the material resolved in time and space and to contribute in the field of the monitoring and the early fault diagnosis of nanoparticle materials.

Optimization of Magnetohydrodynamics parameters on the velocity profile and temperature distribution of incompressible fluid flow on a porous channel was evaluated. The fluid flow was considered to be unsteady, incompressible and flowing in 2-D in porous channel. The effects of magnetic parameter, Darcy’s number and fluid pressure on velocity profiles, and the effect of variable viscosity and Eckert number on temperature distribution in incompressible fluid flow were determined. The flow was considered to be in a channel running along the x-axis on which the magnetic field exist and finally along the y-axis on the porous part of the channel. The resulting model equations were solved by Finite Difference Method (FDM) in MATLAB software. Analysis of results indicated that increasing Darcy’s number and fluid pressure leads to increase in fluid velocity profile, while increasing magnetic parameter decreases fluid velocity profile. Also, it was observed that increase in both Eckert number and fluid viscosity lead to increase in temperature distribution. Optimization in temperature was achieved by increasing the magnetic field while viscosity was optimized by increasing the length of the porous part of the channel. This study will helpful to contribute alternative equations and methodology to engineering and in factories where getting the MHD parameters optimally is the main objective, particularly on temperature, velocity and pressure.

Cancer diseases lead to the second-highest death rate all over the world. The dynamics of invasion of cancer cells into the human body tissues and metastasis are the main causes of death in patients with cancer. This study deals with theoretical investigation of the dynamics of invasion of cancer cells for tumour growths in human body tissues using discretized Cahn-Hilliard, concentration and reaction-diffusion equations which were solved by Finite Difference Method with the aid of MATLAB computer software. A Crank-Nicolson numerical scheme was developed for the discretized model equations. The numerical result obtained was used to describe the dynamics of cancer invasion of tissues with respect to cancer cells density on tumour growth, turbulence and mobility and equilibrium between charge and discharge of cancer cells. The results of the study provide new insights into combating cancer disease by providing mitigating and intervention measures to this major health problem.

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.

Abstract-Digital Phase Plane is presented in this paper. For constant input second order digital systems with arbitrary initial conditions for arbitrary quantizaters, this digital phase plane shows, in plot form, the state trajectory and the asymptotic stability, including the existence and the position of arbitrary period limit cycles.

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

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

An estimation of the electric sources in the heart were conducted using a novel method, based on Huygens’ Principle, aiming at a direct estimation of equivalent bioelectric sources over the heart surface in real time. The main scope of this work was to establish a new, fast approach for the solution of the inverse electrocardiography problem. The study was based on recorded electrocardiograms (ECGs). Based on the Huygens’ principle, measurements obtained from a patient’s thorax surface were interpolated over the surface of the employed volume conductor model and considered as a secondary Huygens’ sources. These sources, being non-zero only over the surface under study, were employed to determine the weighting factors of the eigen-functions expansion describing the generated voltage distribution over the whole volume conductor volume. With the availability of the potential distribution stemming from measurements, the electromagnetics reciprocity theorem is applied once again to yield the equivalent sources over the pericardium. The methodology is self validated, since the surface potentials calculated from these equivalent sources are in very good agreement with ECG measurements. The ultimate aim of this effort is to create a tool providing the equivalent epicardial voltage or current sources in real time, i.e. during the ECG measurements with multiple electrodes.

Model-based approaches form the basis of the majority of damage detection and localization studies; however, built-in online detection technique can offer an alternative for future structural health monitoring technologies. Thus, in this paper, we present a dynamic localization algorithm for damage localization in functionally graded plates. The method depends on the creation of a grid matrix that contains the dynamic response of the structure over time. Subsequently, an optimization process is carried out via a linear equation representing the information contained within the grid to achieve accurate localization of damage. Due to the inherent sparsity of the localization nature, the FISTA algorithm has been employed as a problem solver. The approach is tested in the case of the functionally graded plate under clamped free boundary conditions. Several, damage scenarios were investigated for the case of damage signals on and off the grid. The results show that the proposed approach is capable of accurately predicting damage position and can be suitable for use in low size data systems.

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

Brain tumors are abnormal cells growth in the brain tissues that can be cancerous, or not. In any case, they could be a very aggressive disease that should be detected as early as possible. Usually, magnetic resonance imaging (MRI) is the main tool commonly adopted by neurologists and radiologists to identify and classify any possible anomalies present in the brain anatomy. In the present work, an automatic unsupervised method called Z2-$\gamma$, based on the use of adaptive finite-elements and suitable pre-processing and post-processing techniques, is introduced. In particular, the proposed methodology is able to automatically classify whether a given MR image represents a healthy or a diseased brain and in this latter case, is able to locate the tumor area.

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

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 study investigates the impact of specimen size on the mechanical properties of wood, specifically focusing on compression strength, elastic modulus, and Poisson's ratio. Compression tests were conducted using three different specimen sizes (20mm×20mm×30mm, 40mm×40mm×60mm, 60mm×90mm×90mm) in the longitudinal, radial, and tangential directions. Mechanical parameters, failure mechanisms, load-displacement curves, and stress-strain relationships were systematically analyzed for each size. The study also evaluated the influence of specimen size on the accuracy of finite element numerical simulations by utilizing the obtained mechanical parameters. The results reveal a significant correlation between compressive strength and specimen size, indicating a decrease in compressive strength with increasing specimen size. Conversely, elastic modulus and Poisson's ratio exhibit less sensitivity to variations in specimen size. Notably, parameters derived from small-sized specimens (20mm×20mm×30mm) exhibited substantial errors, while those obtained from medium-sized (40mm×40mm×60mm) and large-sized specimens (60mm×90mm×90mm) demonstrated greater reliability, providing precise results in finite element numerical simulations.

Lightweight component design is effectively achievable through the sandwich structures, many past research studies in the aerospace and racing sectors (since 1920s) has proven it. To extend their application into automotive and other transport industries, the manufacturing cycle times must be reduced. This can be achieved by thermoplastic sandwich materials of continuous fibre-reinforced thermoplastic (CFRTP) cover layers and thermoplastic honeycomb cores. To widespread the application of flat thermoplastic-based sandwich panels into complex parts, a novel forming technology was developed by the Fraunhofer Institute of Microstructure of Materials and Systems (IMWS). For an increase in reproducibility of the sandwich component production and a decrease in manufacturing defects, numerical modelling methods and experimental validation are needed. The finite element (FE) modelling of this thermoforming process is a thermo-mechanical assignment with all three: material, geometry, and boundary non-linearities. In the scope of this work, the material behavior of different sandwich components (cover layers, honeycomb core, interface) under thermo-mechanical loads is characterized in detail, and a complex numerical model is developed to describe the thermoforming process of thermoplastic-based sandwich laminate. The critical forming condition like wrinkling that occurs due to friction in a sandwich or poor intra-ply shear in cross-ply studied, validated, and optimized. A thermoplastic sandwich made of polypropylene (PP) honeycomb core and polypropylene glass fibre (PP/GF) cross-ply as cover layers was used for this study, and numerical modelling was executed in LS-DYNA software.

The insufficient toughness of concrete can lead to the occurrence of crack expansion, ultimately resulting in a decline in the physical and mechanical properties such as bearing capacity and deformation characteristics of structural components. In severe cases, this may lead to brittle failure, where the entire structure and system fail suddenly without significant deformation after being stressed. Additionally, during the process of shear and bending failure of reinforced concrete beams, the cracks generated on the side of the beam are mostly I-II composite cracks. Considering the material nonlinearity and geometric irregularity of reinforced concrete, it is of great significance to study the dynamic evolution law of crack expansion in reinforced concrete members. Based on experimental data, this paper aims to utilize the nonlinear finite element simulation software Abaqus to conduct finite element simulations of crack expansion in simple supported reinforced concrete beams. By accurately controlling the constitutive curves of concrete and reinforcement materials, the design of seismic performance of structural members can be enhanced, thereby making essential contributions to the preliminary analysis of actual crack control engineering. The nonlinear finite element results indicate that the joint cracks of the frame with the floor slab primarily manifest as flexural cracks, with the concrete damage degree of the floor slab being higher than that of the beam. This finding underscores the importance of considering the behavior of concrete and reinforcement materials in designing for seismic performance and crack control in actual engineering applications.

In order to examine the impact of geometric configuration of the human vocal tract on the redistribution of acoustic energy in the frequency spectrum, it is necessary to assemble computation models, which facilitate simple modification of the model’s geometric characteristics. This article describes how the human vocal tract model was constructed using so-called super-elements using the improved reduced system method. The mentioned procedure allows using such a model for demanding optimization calculations, for which the model is modified on the super-element level. Furthermore, the article presents an alternative approach, in which an isoparametric element of a higher degree polynomial was derived. Both approaches were compared and used to find the geometric configuration of the human vocal tract amplifying acoustic energy values on the frequency spectrum to around 3kHz.

In recent years, with the increase in maritime trade, the necessity of increasing the capacities of the ports has emerged. However, while it is planned to increase the capacities of the ports, it is important that the port continues to operate at the same time. In this respect, the old port structures should not be damaged during the capacity increase. In this study, the strengthening of a port in Guinea is discussed as a case study. In the study, the existing quay wall was evaluated, and geotechnical and structural alternatives of the new structure to be built for capacity increase were evaluated. A combined system was designed as a pile foundation and a reinforced foundation with plastic piles so as not to damage the existing quay wall. The pile capacities obtained as a result of the analyses were verified by loading tests.

This work introduces a novel numerical method designed to address three-dimensional unsteady free surface flows incorporating integral viscoelastic constitutive equations, specifically, the K-BKZ-PSM (Kaye–Bernstein, Kearsley, Zapas - Papanastasiou, Scriven, Macosko) model. The new proposed methodology employs a second-order finite difference approach along with the deformation fields method to solve the integral constitutive equation, and the marker-particle method (known as \emph{Marker-and-Cell}) to accurately capture the evolution of the fluid's free surface. The newly developed numerical method has proven its effectiveness in handling complex fluid flow scenarios, including confined flows and extrudate swell simulations of Boger fluids. Furthermore, a new semi-analytical solution for velocity and stress fields is derived, considering fully developed flows of a K-BKZ-PSM fluid in a pipe.

Crystallization often occurs in the laser welding of amorphous alloys, reducing the amorphous alloys’ properties. So, it is necessary to simulate the temperature field curve of the welding area by computer and then determine the power and laser moving speed of laser welding. In this paper, the temperature field curve of Zr41.2Ti13.8Cu12.5Ni10Be22.5(Vit1) amorphous alloy in laser welding is obtained by finite element analysis. The continuous heating curve (CHT) of Vit1 is fitted by the Vogel-Fulcher-Tammann (VFT) equation and the Kissinger equation. If the temperature field curve intersects with the CHT curve, crystallization occurs. The experiment results show that the VFT equation can be used to predict the crystallization of Vit1 better in laser welding. The tem-perature and welding time are increased by using low welding speed. Therefore, the temperature of the weld zone cannot fall in time, resulting intersection of the temperature field curve and the CHT curve. So, crystallization can be avoided if the welding speed is controlled within a rea-sonable range, and the highest temperature is kept under the CHT curve.

Wearable thermoelectric generators (w-TEGs) convert thermal energy into electrical energy to realize self-powering of intelligent electronic devices, thus reducing the burden of battery replacement and charging, and improving the usage time and efficiency of electronic devices. Through finite element simulation, this study successfully designed high-performance thermoelectric generator and made it into wearable thermoelectric module by adopting “rigid device - flexible connection” method. It was found that higher convective heat transfer coefficient on cold-end leads to larger effective temperature difference and better power generation performance of device in typical wearable scenario. Meanwhile, at same convective heat transfer coefficient on the cold-end, longer TE leg length leads to larger temperature difference established at both ends of device, larger device output power and open-circuit voltage. However, when the convective heat transfer coefficient increases to a certain level, optimization effect of increasing TE leg length on device power generation performance will gradually diminish. For devices with fixed temperature difference between two ends, longer TE leg length leads to higher resistance of TEG, resulting in lower device output power but slight increase in open-circuit voltage. Finally, sixteen 16×4×2 mm2 TEGs (L=1.38 mm, W=0.6 mm) and two modules were fabricated and tested. At hot end temperature Th=33 ℃ and cold end temperature Tc=30 ℃, the actual maximum output power Pout of TEG is about 0.2 mW, and the actual maximum output power Pout of TEG module is about 1.602 mW, which is highly consistent with the simulated value. This work brings great convenience to research and development of wearable thermoelectric modules and provides new, environmentally friendly and efficient power solution for wearable devices.

The link of fracture analysis of large ships is an essential part of the structural safety of ships. However, the finite element analysis of large structure according to the traditional method not only has high requirements on the computer, but also the influence of the mesh encryption at special locations makes further requirements on the accuracy of the ship. In this paper, based on the extended finite element method (XFEM) and the multi-point constraint method (MPC), the dynamic crack expansion process of the rectangular plate containing cracks under the multi-scale model is comparatively studied by using ABAQUS simulation software Finally, the applicability of the method is verified by applying the multi-scale model of ship nodes. The results show that the computational effort of the multi-scale model is significantly reduced with the same accuracy and allowable error, which proves the accuracy and efficiency of the method.

This paper presents the experimental and numerical studies of the vertical stress distribution beneath the centerline of circular footings resting on the sand. Vertical stress values due to circular load in medium-dense sand were measured using soil stress transducers. Reduced-scale laboratory model tests were carried out in a square-shaped test box and vertical stresses were measured at the centerline of the circular footing at different depths of the sand. In order to determine the effect of footing width on vertical stress distribution, laboratory experiments were performed on circular footings of different diameters. In addition, a series of finite element analyses was carried out on the test model to validate the findings of the experimental study. In the numerical analysis, the sand was modeled as a non-linear elastoplastic material. The analyses were performed using the finite element method with two-dimensional axisymmetric and three-dimensional conditions. The results obtained from laboratory model tests and numerical analysis were compared with the vertical stress values obtained from the theoretical solution in the literature. Vertical stress distribution in horizontal planes at any depth from the footing base shows the same characteristic form in all methods. However, some differences were observed between experimental, numerical, and theoretical results and discussed.

The Gol'denveizer problem of a torus was studied analytically by Audoly and Pomeau (2002), and the accuracy of the Audoly and Pomeau linear law was confirmed numerically by Sun (2021). However, the law does not include the major radius R of the torus. To find the influence of the major radius, we used finite element numerical simulation to simulate different cases, and we propose a modified Audoly and Pomeau linear law for vertical deformation, which includes R. A linear law of horizontal deformation is presented as well. Our studies show that the Audoly and Pomeau linear law has high accuracy. With modified vertical and horizontal deformation, a displacement-compatible relation between them is formulated.

Orofacial injuries are common occurrences during contact sports activities; however, there is an absence of data regarding the performance of hybrid occlusal splint mouthguards, especially during compressive loading. To evaluate the total deformation and stress concentration, a skull model was selected and duplicated to receive two different designs of mouthguard devices: one model received a conventional custom-made mouthguard (MG) with 4-mm thickness and the other received a novel hybrid occlusal splint-mouthguard (HMG) with the same thickness. Both models were subdivided into finite elements. The frictionless contacts were used, and a nonlinear analysis was performed simulating the compressive loading in occlusion. The results were presented in von-Mises stress maps (MPa) and Total Deformation (mm). A higher stress concentration in teeth was observed for the model with the conventional MG, while the HMG design displayed a promising mechanical response with lower stress magnitude. The HMG de-sign displayed a higher magnitude of stress on its occlusal portion than the MG design. The hybrid mouthguard (HMG) reduced (1) jaw displacement during chewing and (2) the generated stresses in maxil-lary and mandibular teeth.

Our aim is to implement and test some less known numerical methods that might be used to solve an advection equation. The description is restricted to information necessary to implement and test the presented methods and it is given only for one-dimensional case.

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.

Since large power transformers are custom-made, and their design process is a labor-intensive task, their design process is split into different parts. In tendering, the price calculation is based on the preliminary design of the transformer. Due to the complexity of this task, it belongs to the most general branch of discrete, non-linear mathematical optimization problems. Most of the published algorithms are using a copper filling factor based winding model to calculate the main dimensions of the transformer during this first, preliminary design step. Therefore, these cost optimization methods are not considering the detailed winding layout and the conductor dimensions. However, the knowledge of the exact conductor dimensions is essential to calculate the thermal behaviour of the windings and make a more accurate stray loss calculation. The paper presents a novel, evolutionary algorithm-based transformer optimization method which can determine the optimal conductor shape for the windings during this examined preliminary design stage. The accuracy of the presented FEM method was tested on an existing transformer design. Then the results of the proposed optimization method have been compared with a validated transformer design optimization algorithm.

Adaptable or adapted? Whether it is a question of physical, biological or even economic systems, this problem arises when all these systems are the location of matter and energy conversion. To this interdisciplinary question we propose a theoretical framework based on the two principles of thermodynamics. Considering a finite time linear thermodynamic approach, we show that non-equilibrium systems operating in quasi-static regime are quite deterministic as long as boundary conditions are correctly defined. The Novikov-Curzon-Ahlborn approach [1,2] applied to non-endoreversible systems then makes it possible to precisely determine the conditions for obtaining characteristic operating points. As a result, power maximization principle (MPP), entropy minimization principle(mEP), efficiency maximization, or waste minimization states are only specific modalities of system operation. We show that boundary conditions play a major role in defining operating points because they define the intensity of the feedback that ultimately characterizes the operation. Armed with these thermodynamic foundations, we show that the intrinsically most efficient systems are also the most constrained in terms of controlling the entropy and dissipation production. In particular, we show that the best figure of merit necessarily leads to a vanishing production of power. On the other hand, a class of systems emerges which, although they do not offer extreme efficiency or power, have a wide range of use and therefore marked robustness. It therefore appears that the number of degrees of freedom of the system leads to an optimization of the allocation of entropy production.

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.

In this paper we study the normal bases for Galois ring extension ${{R}} / {Z}_{p^r}$ where ${R}$ = ${GR}$(p^r, n). We present a criterion on normal basis for ${{R}} / {Z}_{p^r}$ and reduce this problem to one of finite field extension $\overline{R} / \overline{Z}_{p^r}=F_{q} / F_{p} (q=p^n)$ by Theorem 1. We determine all optimal normal bases for Galois ring extension.

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.

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 influence of applied loads on free vibrations analysis of plates and shells panels made of functionally graded materials is analyzed in the present work. Formulations for the static analysis considering geometrically nonlinear behavior, as well as linear buckling and free vibrations analyses are considered. The calculation of the through-thickness stress distribution is also performed. A finite element model based on a higher order shear deformation theory and using a non-conforming triangular flat plate/shell element with 3 nodes, and 8 degrees of freedom per node is used in the numerical implementation. The results obtained with this numerical model are presented, discussed, and compared with alternative solutions published by other authors in some benchmark applications.

Starting from the analysis of the lack of positivity of the Cattaneo heat equation, the article addresses the thermodynamic relevance of the positivity constraint in irreversible thermodynamics, than is at least as important as the entropic constraints. The fulfillment of this condition in hyperbolic models leads to the parametrization of the concentration fields with respect to internal variables associated with the microscopic dynamics. Using Brownian motion theory as a landmarking example for deriving macroscopic transport equations from the equations of motion at the particle/molecular level, we discuss two typical problems involving hydrodynamic interactions at microscale: surface chemical reactions at a solid interface of a diffusing reactant, and mass-balance equations in a complex viscoelastic fluid, in which the physics of the interaction leads either to the overcome the parabolic diffusion model or to consider the parametrization of the concentration with respect to the degrees of freedom associated with the relaxation dynamics of the solvent fluid.

As a new type of gear, bevoloid gears have the advantages of compensating axial error, smooth transmission, and eliminating turning error, and they are widely used in applications requiring high transmission accuracy and stability. However, research on calculating the time-varying mesh stiffness of bevoloid gears is still limited, and there is an urgent need to propose a method that can calculate the stiffness of bevoloid gears quickly and accurately. This paper first establishes the bevoloid gear tooth profile expressions, then assumes a pair of bevoloid gears meshing with the same rack and derives the contact line equations of parallel axis bevoloid gear pairs, and analyze the contact process of bevoloid gears. We propose an analytical algorithm that uses the slicing method to calculate the stiffness of helical gears, straight bevoloid gears, and helical bevoloid gears, change the parameters of helical bevoloid gears respectively, and analyze the influence of different parameters on stiffness. Finally, the finite element method is used to verify the analytical method, and the correctness of the analytical calculation results is verified, and the errors are analyzed.

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

One of the most important problems of geotechnical engineering is the design of pile foundations. Piles are generally designed for vertical and horizontal static loads. However, it is important to design piles against dynamic loads in regions with high seismicity. This study verified a centrifuge experiment obtained from the literature with a finite element program. With the kinematic interaction analyses made for this purpose, the effect of parameters such as layered soil, groundwater, pile stiffness, and earthquake acceleration were investigated numerically. According to the results obtained, it was understood that the earthquake magnitude is the most important parameter in kinematic interaction analyses. In addition, in layered soil conditions, the fact that 75% of the pile length is in the weak soil creates the most unfavorable situation. It has been observed that pile stiffness and groundwater level also have certain effects on the kinematic interaction.

The flavor, aroma and color of coffee may change due to mechanical damage, reducing its qual-ity. To measure the mechanical response of materials, compression tests can be performed, de-termining the effect of force application on the mechanical behavior of the fruit at different stages of ripeness. In this context, this study aimed to analyze the deformation, strain energy and equivalent von Mises stress of coffee fruits at mature, semi-mature and immature maturation stages and compressed by collapse forces. Compression in three directions (x, y and z) was sim-ulated on coffee fruit models using the finite element method, with a compression support in the direction opposite to the force application axis. The numerical simulation of the compression process allowed to verify that the more mature the fruit, the greater the associated mean defor-mation (2.20 mm mm-1, 0.78 mm mm-1 and 0.88 mm mm-1), the lower the mean strain energy (0.07 mJ, 0.21 mJ and 0. 34 mJ) and the lower the mean equivalent von Mises stress (0.25 MPa, 1.03 MPa and 1.25 MPa), corresponding to ripe, semi-ripe and immature fruits; and these proper-ties may be used to analyze how processes in the production chain affect the amount of mechan-ical damage.

River renaturation can be an effective management method for restoring the floodplain's natural capacity and minimizing the effects during high flow periods. A 1D-2D HEC-RAS model, in which the flood plain was considered as 2D and the main channel as 1D, was used to simulate flooding in the restored reach of the Spree River. When computing in this model, finite volume and finite difference approximations using the Preissmann approach are used for the 1D and 2D models, respectively. To comprehend the sensitivity of the parameters and model, several scenarios were simulated using different time steps and grid sizes. Additionally, dykes, dredging, and changes to the vegetation pattern have been used to simulate flood mitigation measures. The model predicted that flooding would occur mostly in the downstream portion of the channel in the majority of the scenarios without mitigation measures, whereas with mitigation measures, flooding in the floodplain would be greatly reduced. By preserving the natural balance on the channel's floodplain, the restored area needs to be kept in good condition. Therefore, mitigating measures that balance the area's economic and environmental aspects must be considered in light of the potential for floods.

Advanced manufacturing techniques, aimed at implants with high dependability, flexibility, and low manufacturing costs, are crucial in meeting the growing demand for high-quality products like biomedical implants. Incremental sheet forming is a promising flexible manufacturing approach for rapidly prototyping sheet metal components using low-cost tools. Titanium and its alloys are used to shape most biomedical implants because of its superior mechanical qualities, biocompatibility, low weight, and great structural strength. The poor formability of titanium sheets at room temperature, however, limits their widespread use. The goal of this research is to show that gradual sheet formation of a titanium biomedical implant is possible. The possibility of creative and cost-effective concepts for the manufacturing of such complicated shapes with significant wall angles is explored in this study. A numerical simulation based on finite element modeling as well as a design process tailored to metal forming is used to complete the development. The mean of uniaxial tensile tests with a constant strain rate was used to study the flow behavior of the studied material. To forecast the crack, the obtained flow behavior was modeled using the behavior model and failure model.

A connection between the phenomenon of wafer bifurcation and that of spontaneous symmetry breaking (SSB) has been established. Indeed, by developing an analytical approximation of the elastic energy of a wafer coated with a thin layer (e.g. a metal), it is shown as the elastic potential energy, interpreted as a quantity contributing to the thermodynamic free energy, can be investigated within the framework of the Landau theory of the second order phase transitions. The elastic energy of a bifurcated wafer is a complex function of the stress and the curvatures in the two perpendicular directions. In this work, it is shown as a translation of the coordinate system of the curvatures allows to gain a potential which has a “Mexican hat” shape. This is a distinctive trait of the phenomenon of spontaneous symmetry breaking (SSB). Moreover, it is shown as the values of the coordinates at the minimum of the SSB potential agrees with those provided by the theory. Bifurcation is hence a phenomenon that can also be interpreted as a spontaneous symmetry breaking where the rotational symmetry of a disc shaped wafer is broken. It occurs because of a lowering of the wafer energy in the SSB broken symmetry configuration.

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

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

With the increasingly stringent airworthiness standards, the noise generated during the rotorcraft flight is gradually attracting people’s attention. It widely operated helicopters at low altitudes because of their maneuverability. The way to reduce the noise caused by the complex airflow of the helicopter rotor system has progressively become a hot topic for researchers. Using a hybrid acoustic analysis method, this paper investigates the improvement of the noise and thrust of the helicopter’s tail rotor through the tail rotor structural parameters. For the basic model, the turbulence simulation is performed using an incompressible detached eddy simulation (DES) method, and the Lighthill acoustic analog equation is calculated using the finite element method (FEM). We verified the accuracy of the method through wind tunnel tests. We chose a series of structural parameters for sound simulation and fluid simulation calculations. The results indicate that the modified tail rotor noise reduced by 16.5 dBA and the total thrust increased by 19.9% from the prototype model. This work can enhance the duct tail rotor design to improve aerodynamic and aeroacoustic performance.

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.

Metasurfaces of gold (Au) nanoparticles on a SiO2-Si substrate were fabricated for the enhancement of second harmonic generation (SHG) using electron beam lithography and lift-off. Triangular Au nanoprisms which are non-centro-symmetric and support the second- order non-linearity were examined for SHG. The thickness of the SiO2 spacer is shown to be an efficient parameter to spectrally tune to maximise SHG. Electrical field enhancement at the fundamental wavelength was shown to define the intensity of the second harmonics. Numerical modeling of light enhancement was verified by experimental measurements of SHG and reflectivity spectra at the normal incidence. At the plasmonic resonance, SHG is enhanced up to ∼3.5×103 times for the optimised conditions.

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

In this paper, the results of finite element analyses of a single-layer cylindrical latticed shell under severe earthquake is presented. A 3D Finite Element model using fiber beam elements were used to investigate the collapse mechanism of this type of shell. The failure criteria of structural members are simulated based on the theory of damage accumulation. Severe earthquake of peak ground acceleration (PGA) of 500 gal was applied to the shell. The stress and defeomation of the shell were studied in detail. A three-stages collapse mechanism “double-diagonal -members-failure-belt” of this type of structure was discovered. Based on the analysis results, the measures to mitigate collapse of this type of structure is recomanded.

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.

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.

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

A novel compressible real gas 1D finite volume thermofluid network modelling methodology is proposed that can be applied in the analysis of supercritical CO₂ (sCO2) power cycles. It is applied to simulate centrifugal compressor performance via a mean line approach that directly solves the mass, energy, and momentum balance equations for the discretized flow channels within the machine along with the associated loss models, component characteristics, and applicable compressible real gas equation of state. This can be integrated seamlessly with the rest of the power cycle components like pipes, valves, heat exchangers, etc., thereby eliminating the need for static performance maps and the associated difficulties to ensure similarity over the whole range of operation. This direct simulation of the internal flows also makes it possible to do integrated cycle design optimization that includes the geometry of the compressors, and to simulate fast transients for dynamic analysis. The methodology was applied to simulate well-documented air and sCO2 centrifugal compressors and compared the generated results with that of conventional thermofluid network incompressible and compressible ideal gas approaches. The results indicate that the proposed compressible real gas approach can capture the stagnation pressure changes and isentropic efficiencies with good accuracy for both cases over the whole range of mass flow rates and shaft rotational speeds. In contrast, the conventional thermofluid network incompressible and compressible ideal gas approaches are only applicable in selected regions of the fluid property regime and suffer from discretization errors and inappropriate treatment of the relationship between stagnation and static pressures. Furthermore, comparative analyses between the newly proposed approach and the traditional mechanical energy balance method are performed, showing that the latter approach needs to be discretized more finely through the compression process to correctly capture the polytropic compression path shape, whereas the former more accurately captures the stagnation pressure changes through the analyzed compressors with only a single increment per compressor component.

Modal testing is a common step in aerostructure design, serving to validate the predicted natural frequencies and mode shapes obtained through computational methods. The strategic placement of sensors during testing is crucial to accurately measuring the intended natural frequencies. However, conventional methodologies for sensor placement are often time-consuming and involve iterative processes. This study explores the potential of machine learning techniques to enhance sensor selection methodologies. Three machine learning-based approaches are introduced and assessed, comparing their efficiency with established techniques. The evaluation of these methodologies is conducted using a numerical model of a beam to simulate real-world scenarios. The results offer insights into the efficacy of machine learning in optimizing sensor placement, presenting an innovative perspective on enhancing the efficiency and precision of modal testing procedures in aerostructure design.

It is commonly recognized that a ship’s ability to maintain stability is dependent upon the way its loads are distributed. This means making sure the center of gravity is sufficiently low and evenly distributed. Yet when a ship encounters turbulent water with extremely large wave amplitudes, load distribution is a critical factor in determining how stressed the ship becomes. In particular, this article examines several methods for increasing the ship's longitudinal strength to prevent ship structural damage during hogging and sagging. This research presents a number of proposed techniques to enhance the longitudinal strength of cargo ships in order to guarantee minimal stress and deformation during sagging and hogging. The first method is the stern-to-bow use of curved support structures to increase the longitudinal strength of a cargo ship. And another way to stiffen the lateral plate is the proposed use of different designs of the side plate, such as wavy curves, stars, etc. To further boost the strength, the lateral plate upper part is suggested to curve upward. A numerical simulation using ANSYS of the stiffened plate when subjected to uniform load distribution and fixed support on either end shows that there is a decrease in deformation of the plate when we use wavy curved plates compared with those of the other plate stiffening methods used, and furthermore, curving the plate on top increases its strength during sagging as well as hogging. Finally, the use of a curved support structure further enhanced the overall longitudinal strength of the vessel.

This article introduces a new physically based Lyapunov function difinition, which can be used as acoustic energy for finite amplitude pressure oscillations in reacting systems especially in combustion chambers. Reacting flow is seen as an open, non-equilibrium (parameters are distributed unevenly locally) thermodynamic system. This Lyapunov function is defined as the maximum mechanical work, which could be extracted by a heat engine from studied system, if this system would be disconnected from inlet and outlet and from any other surrounding environment, and the engine could transfer to the surrounding environment only mechanical work.

In this paper, we give all non splitting bi-unitary superperfect polynomials divisible by one or two irreducible polynomials over the prime field of two elements. We prove the nonexistence of odd bi-unitary superperfect polynomials over F2.

Precipitates are the primary source of strength for the Al-Mg-Si alloy. Aluminum alloy in the peak-aged state mainly contain β'' and β' precipitates. In the past, most of the literature only considered the strengthening effect of β''. Here we develop a single crystal intensity model including both precipitate enhancement effects for the first time. This model was subsequently implemented into a crystal plastic finite element method to model the uniaxial tensile process of a polycrystalline aggregate model of Al-Mg-Si alloy. The simulation results for uniaxial stretching are in good agreement with the experimental results, confirming that the constitutive parameters used for the single crystal strength model with two precipitates are based on realistic physical implications. Furthermore, by comparing the uniaxial tensile simulation results of a peak-aged alloy considering the actual precipitated phase composition of the alloy with those assuming that the precipitated phase is only the β'' phase, the predicted tensile strength of the former is around 5.65% lower than that of the latter, suggesting that the two kinds of precipitation should be considered separately when simulating the mechanical response of Al-Mg-Si alloy. It is highly expected that the present simulation strategy is not limited to the Al-Mg-Si alloys, and it can be equally applied to the other age-enhanced alloys.

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

According to the development needs of marine seismic station networking, it is necessary to develop a marine seismometer with a networking function to realize quasi-real-time and large-scale long-term observation of undersea seismic signals. This can lay a solid foundation for marine seismic research and research on earth tectonics, activities, and processes. Based on the advantages and disadvantages of the existing observation equipment, this study proposed to develop a new type of the Autonomous Glide Marine Seismometer(AGMS) that can be applied to the global sea area, focusing on the design and study of the main structure of the carrier of the AGMS. The ultimate strength of the ellipsoidal pressure-resistant spherical shell is investigated by using the finite element analysis method. A reasonable finite element analysis model is proposed through comparing the finite element parameter results with the theoretical formulas and relevant specifications. Finally, the design process and strength checking of the pressure-resistant shell of the AGMS are established. These relevant results can also provide a reference for the design of other deep-sea pressure-resistant chambers.

This study presents finite-time synchronization of chaotic fractional order systems with disturbance uncertainty and unknown time delay. First, a PID sliding surface is presented. Then, for the finite-time synchronization of the master and slave systems, a robust sliding-adaptive control method is provided. Update rules have been retrieved to estimate the system parameters by using the Lyapunov function while establishing the stability of the suggested mechanism and assuring the convergence of errors to zero. The proposed approach was used to solve the fractional order system of Gensou-Tsei with time-varying parameters, and the simulation results show that the presented approach performs well.

This research work aims to develop a mathematical modeling for production-inventory systems. Two mathematical models are present: An infinite dimensional, partial differential equation (PDE), production level modeling is present and a finite dimensional system for a coupled dynamic pricing, production rate and inventory level ordinary differential equations (ODE), integrating a proper Lyapunov stability analysis of the dynamical system and simulations.

Arthroplasty is commonly performed to treat advanced osteoarthritis or other degenerative joint conditions, but it can also be considered for young patients with severe joint damage that significantly limits their functionality and quality of life. Young patients are still at risk of aseptic mobilization and bone resorption due to the phenomenon of stress-shielding that causes an uneven distribution of tensions along the femoral contact surface prosthesis. This phenomenon can be limited by appropriately choosing the material of the prosthesis or by varying its stiffness, making sure that its mechanical behavior simulates that of the femur as much as possible. The aim of this study is to evaluate the mechanical behavior of a prosthesis optimized both in shape and material and in comparing the results with a standard titanium prosthesis. Methods: Through three-dimensional modeling and the use of software finite element method (FEM) such as ANSYS we simulate the mechanical behavior of traditional prosthesis and prosthesis optimized topologically respecting the ASTM F2996-13 standard Results: With topological optimization there is a stress reduction from 987 Mpa to 810 Mpa with a mass reduction of 30%. With the use of carbon fiber, it is possible to further reduce stress to 509 Mpa. Conclusions: The reduction of stress on the femoral stem allows an optimal distribution of the load on the cortical bone thus decreasing the problem of stress-shielding.

We consider an $N$ fermion system at low temperature $T$ in which we encounter special particle number values $N_m$ exhibiting special traits. These values arise in focusing attention upon the degree of mixture (DM) of the pertinent quantum states. Given the coupling constant of the Hamiltonian, the DMs stay constant for all $N$-values, but experience sudden jumps at the $N_m$. For a quantum state described by the matrix $\rho$, its purity is expressed by $Tr \rho^2$ and then the degree of mixture is given by $1 - Tr \rho^2$, a quantity that coincides with the entropy $S_q$ for $q=2$. Thus, Tsallis entropy of index two faithfully represents the degree of mixing of a state, that is, it measures the extent to which the state departs from maximal purity . Macroscopic manifestations of the degree of mixing can be observed through various physical quantities. Our present study is closely related to properties of many-fermion systemsn that are usually manipulated at zero temperature. Here we wish to study the subject at finite temperature. Gibbs' ensemble is appealed to. Some interesting insights are thereby gained.

Multi-chiplet technique is expected to be a promising solution to achieve high-density system integration with low power consumption and high usage ratio. This technique can be integrated with glass interposer to accomplish a competitive low fabrication cost compared with the silicon-based interposer architecture. In this study, process-oriented stress simulation is performed by the element activation and de activation technique in finite element analysis architecture. Submodeling technique is also utilized to mainly conquer the scale mismatch and difficulty in mesh gridding design. It is used as well to analyze the thermomechanical responses of glass interposer with chiplet arrangement and capped epoxy molding compound (EMC) during curing. A three-factor, three-level full factorial design is applied using the analysis of variance method to explore the significance of various structural design parameters on stress generation. Analytic results reveal that the maximum first principal stresses of 130.75 and 17.18 MPa are introduced on the sidewall of Cu-filled via and the bottom of glass interposer, respectively. Moreover, the EMC thickness and through glass via pitch are the dominant factors in the adopted vehicle. They significantly influence the stress magnitude during heating and cooling.

This paper proposed a time two-mesh (TT-M) finite difference numerical scheme to improve the efficiency of solving the symmetric regularized long wave (SRLW) equation. The TT-M Crank-Nicolson discretization and finite difference method are employed in time and space approximation respectively. The scheme involves three main steps: firstly, the time interval is divided into coarse and fine time meshes, then the nonlinear system is solved on the coarse time mesh; secondly, coarse numerical solutions on the fine time mesh are computed using an interpolation formula based on the solutions derived in the step one; lastly, the TT-M finite difference numerical solutions can be obtained through constructing the linearized fine time mesh system using Taylor’s formula. Compared to the currently existing TT-M numerical methods, the novelty of this study is that the nonlinear term including derivatives is linearized by Taylor’s formula for a function with three variables, whose error analysis is more complex. Finally, some numerical examples, including computational time and accuracy, preservation of conservation laws, are given to verify the efficiency of the scheme. By comparing it with the standard nonlinear finite difference scheme, this method can reduce CPU time without sacrificing accuracy.

Air conditioners are indispensable necessities in daily life and industrial production. However, the heat generation of their internal core power components can limit their performance release. In this work, we explore the application of thermoelectric coolers (TECs) in the field of power device heat dissipation through finite element simulation. First, we geometrically modeled the structure and typical operating conditions of core power devices in air conditioners. We compared the temperature fields in air cooling and TEC active cooling modes for high power consumption power devices in a 319 K operating environment. The simulation results show that in the single air cooling mode, the maximum temperature of the 173.8 W power device reaches 394.4 K and the average temperature reaches 373.9 K, which exceeds its rated operating temperature of 368.1 K. However, after adding TEC, the maximum temperature of the power device drops to 331.8 K at an operating current of 7.5 A and the average temperature of the device is 326.5 K. It indicates that TEC active cooling has a significant effect on the temperature control of the power device. At the same time, we also studied the effect of TEC working current on the temperature control effect of power devices. For TEC, there is a minimum working current of 5 A; when it is less than 5 A, TEC has no cooling effect; the cooling effect of TEC increases with an increase in working current. When the TEC working current is 10 A, the average temperature of the power device can be reduced to 292.9 K. This study has made a meaningful exploration of the application of TEC in chip temperature control and heat dissipation, providing a new solution for chip thermal management and accurate temperature control.

Rail impedance directly affects the transmission performance of track circuit . Considering the condition of earth stratification, for the difficult to calculate the rail impedance due to the semi-infinite integration interval and the oscillation of the integrand by using the Carson formula, The truncation method is proposed to divide the impedance formula is divided into definite integral and tail integral. The integral is approximated by the spline function, and the tail integral is calculated by using the exponential integral and Euler formula. Based on it, the rail impedance calculation formula of track circuit is obtained. The electromagnetic field model of track circuit with earth stratification is simulated by finite element method, and the correctness of the method is verified. Based on the formula, the influence of current frequency, soil depth and conductivity on rail impedance is studied. The relative error between the calculated results of rail impedance and the simulation results of finite element is within 5%. It can be seen that the formula has high accuracy and correctly reflects the law of rail impedance variation with current frequency, soil depth and resistivity. It provides a reliable reference for the theoretical calculation of rail impedance of track circuit.

The non-linear analysis of hollow-core concrete slabs requires the use of advanced numerical techniques, proper constitutive models both for concrete and steel as well as particular computational skills. If prestressing, cracking, crack opening, material softening, etc. are also to be taken into account, then the computational task can far exceed the capabilities of an ordinary engineer. In order for the calculations to be carried out in a traditional design office, simplified calculation methods are needed. Preferably based on the linear finite element (FE) method with a simple approach that takes into account material nonlinearities. In this paper the simplified analysis of hollow-core slabs based on the generalized nonlinear constitutive law is presented. In the proposed method a simple decomposition of the traditional iterative linear finite element analysis and the non-linear algebraic analysis of the plate cross-section is used. Through independent analysis of the plate cross-section in different deformation states, a degraded plate stiffness can be obtained, which allows iterative update of displacements and rotations in the nodes of the FE model. Which in turn allows to update the deformation state and then correct translations and rotations in the nodes again. The results obtained from the full detailed 3D nonlinear FEM model and from the proposed approach are compared for different slab cross-sections. The obtained results from both models are consistent.

This study investigates the failure of a roof with steel truss construction of a factory building in Tekirdag in North-western part of Turkey. The failure occurred under hefty weather conditions including thunderbolt, lightning strikes, heavy rain and fierce winds. In order to interpret the reason for the failure, the effects of different combinations of factors on the design and dimensioning of the roof are checked. Therefore, finite element analysis is performed several times under different assumptions and considering different factors aiming to determine the dominant ones that are responsible for the failure using the commercial software Abaqus (Dassault Systèmes, Vélizy-Villacoublay, France). Each loading condition gives out a characteristic form of failure. The scenario with the most similar form of failure to the real collapse is considered as the most likely scenario of failure. Also, the factors included in this scenario are expected to be the responsible factors for the partial collapse of the steel truss structure.

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.

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.