The successful handling of large semiconductor wafers is crucial for scaling up their production. Early-stage warpage control allows to prevent undesirable asymmetric warpage, known as wafer bifurcation or buckling. Indeed, even in a gravity-free environment, thinning an 8" or 12" semiconductor wafer can result in warpage and bifurcation. To mitigate this issue, the taiko method, which involves creating a thicker ring region around the rim of the wafer, has been widely used. Previous research has focused on the theoretical factors affecting the warpage of a backside metalized taiko wafer. This work extends the case to a front-side metalized taiko wafer and introduces the concept of the equivalent thickness of a taiko wafer. The equivalent thickness of a taiko wafer, influenced by the ring region, lies in between the thickness of the central region and that of the annular region. Because of the limited number of taiko wafers that can be produced on a production line, modelling can be beneficial. In this work we compared the results of a developed analytical model with those obtained from a finite element analysis (FEA) approach with ANSYSY® software to model the equivalent thickness of a taiko wafer. We investigated the curvature as a function of the stress of the metal layer, considering key design factors such as the substrate region thickness, the thickness of the thin metal film, the step height, and the width of the ring region. By systematically varying the thickness of the central region of the taiko wafer, we explored the curvature as a function of stress induced by thermal load in the linear regime and determined the slopes in the linear region of the curvature vs. stress curves. The aim of this study is to identify regularities and similarities with the Stoney equation and investigate the validity of the analytical approach for the case of a taiko substrate. It results that there is a good agreement between the analytical model of a taiko wafer and the numerical analysis gained by the FEA methods.

This paper presented a two-dimensional, well-balanced hydrodynamic and sediment transport model based on the solutions of variable density shallow water equations (VDSWEs) for sediment-laden flows, and the Exner equation for bed changes. Those equations are solved in a coupled way by the first-order Godunov-type finite volume method. The Harten-Lax-van Leer-Contact (HLLC) Riemann solver is extended to find the local Riemann fluxes in order to maintain the exact balance between the momentum term and the bed slope term. The advantage of a well-balanced scheme over an unbalanced scheme is demonstrated by the synthetic standing contact-discontinuity test case. Then, the model is employed to simulate two laboratory experiments. At last, a field case, the 1996 Lake Ha! Ha! flood event (Canada), is simulated. Results of cross sectional geometries and profiles of longitudinal thalweg agree well with measurements. The accuracy and simplicity of the numerical model, together with the robust implementation, make the model a good candidate for practical engineering applications.

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.

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.

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.

The inverse Finite Element Method (iFEM) is a model-based technique to compute the displacement (and then the strain) field of a structure from strain measurements and a geometrical discretization of the same. Different literature works exploit the error between the numerically reconstructed strains and the experimental measurements to perform damage identification in a Structural Health Monitoring framework. However, only damage detection and localization are performed, without attempting a proper damage size estimation. The latter could be based on machine learning techniques, however, an a priori definition of the damage conditions would be required. To overcome these limitations, the present work proposes a new approach in which the damage is systematically introduced in the iFEM model to minimize its discrepancy with respect to the physical structure. This is performed with a maximum likelihood estimation framework, where the most accurate damage scenario is selected among a series of different models. The proposed approach is experimentally verified on an aluminum plate subjected to fatigue crack propagation, which enables the creation of a Digital-Twin of the structure itself. The strain field fed to the iFEM routine is experimentally measured with an Optical Backscatter Reflectometry fiber and the methodology is validated with independent observations of lasers and the Digital Image Correlation.

In this paper, we present the investigation results of radiation-induced effects in metal-oxide-semiconductor (MOS) doped with moderate amounts of Zinc Oxide (ZnO) as a potential candidate for space-borne application. The samples were fabricated via the sputtering method at a working pressure of 3mTorr and a deposition temperature of 300oC. The ZnO samples were exposed to 1.25-MeV gamma-ray utilizing Co60 source, and their electronic response was measured at ionizing doses ranging from 10 kGy to 300 kGy. A comparative work was performed through finite element method to simulate the electronic response of the PN junction diode due to ionizing radiation. The results indicate that the ideality factor of the MOS diode increases as the ionizing dose increases, rendering it unsuitable for use as a diode. The degradation of the electrical parameters was also simulated, showing the increase in hole concentration. These findings suggest that the ejection of electrons occurred, which agrees with the gamma radiation effects trend. Furthermore, as the intensity of radiation increases, the spatial charge that arises from the separation of hole-electron pairs results in a substantial reduction of the electric field in the central portion of the n-type region. These findings provide insights into the degradation of electrical parameters in MOS devices under gamma radiation and have implications for their use in space-borne applications.

H-adaptivity is an effective tool to introduce local mesh refinement in FEM-based numerical simulation of crack propagation. The implementation of h-adaptivity could benefit the numerical simulation of fatigue or accidental load scenarios involving large structures such as ship hulls. In engineering applications, the element deletion method is frequently used to represent cracks. However, the element deletion method has some drawbacks such as strong mesh dependency and loss of mass or energy. In order to mitigate this problem, the element splitting method could be applied. In this study, a numerical method called ‘h-adaptive element splitting’ (h-AES) is introduced. The h-AES method is applied in FEM programs by combining h-adaptivity with the element splitting method. Two examples using the h-AES method to simulate cracks in large structures under linear-elastic fracture mechanics scenario are presented. The numerical results are verified against analytical solutions. Based on the examples, the h-AES method is proven to be able to introduce mesh refinement in large-scale numerical models that consist of structured coarse meshes. By employing the mesh refinement introduced in this paper, very small cracks are well represented in large structures.

In this study, a numerical simulation of fluid flow through railway ballast in the time domain is presented, providing a model for unsteady-state flow. It is demonstrated that the position of the free surface with respect to time can also be used to solve the steady flow case. The effect of ballast fouling is included in the model to capture the realistic behavior of railway ballast, which is critical to understanding the impact of flooding. A thorough comparison with a range of previous studies, including theoretical and experimental approaches, is made and very close agreement is obtained. The significant impact of ballast fouling on fluid flow and its potential consequences for railway infrastructure are highlighted by the simulation. Valuable insights into the behavior of water flow through porous media and its relevance to railway ballast management are offered by this study.

To understand the enhanced protection mechanism of CoCrNiAlY-YSZ-LaMgAl11O19 dou-ble-layer ceramic coating with aluminum plating, a finite element simulation method was used to simulate the distribution of thermal stress in the coating in all directions. The results show that in the air exposure of the un-aluminized coating, high temperature causes a large radial thermal stress on the surface of the LaMgAl11O19 (LMA) layer, and it increases with the increase in temperature, which is the main reason for the initiation of axial cracks. After arc aluminum plating, the aluminum plating layer effectively inhibited the volume shrinkage of the coating through good adhesion to the coating and internal diffusion, the thermal stress of the coating was considerably reduced, and the CoCrNiAlY-YSZ-LMA coating had an effective enhancement and protection effect; however, there was still a certain amount of shear thermal stress inside the LMA layer, the top of the crack, and the bottom of the crack. This thermal stress caused the initi-ation of radial microcracks in the LMA layer, which also becomes a risk point for the failure of the aluminum coating.

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.

In this work we explore the utilization of plasmonic resonance (PR) in silver nanowires to enhance the performance of organic solar cells. Plasmonic resonance is a phenomenon in which nanoscale conductive materials exhibit oscillation of conduction electrons, resulting in the creation of an electric field. Enhancing light absorption is crucial for improving organic solar cell efficiency and incorporating metallic nanostructures to induce surface plasmon resonance (SPR) shows promise in achieving this goal. We discuss the two key mechanisms of plasmonic effects: far-field scattering and near-field resonance modes. Far-field scattering extends the optical path of incident light, while near-field plasmonic effects involve localized surface plasmon resonance (LSPR) and plasmonic cavity modes, enhancing absorption by strengthening electric fields near the nanostructures. Silver nanowires are the focus of this study, and finite-difference time-domain (FDTD) simulation software is used to investigate their plasmonic resonance behavior in a ZnO/Silver nanowires/ZnO (ZAZ) electrode structure. The simulations reveal the dominance of LSPR in this configuration, with intense electric fields inside the nanowire and propagation into the surrounding medium, offering opportunities for enhanced light absorption in the organic solar cell's active layer.

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.

The main objective of the paper is to investigate the wave propagation behaviour in composites using numerical and experimental methods. Specifically, the study focuses on the actuation of waves by two PZT transducers placed on the surface and embedded at different positions within the composites, namely the quarter and middle positions. The numerical analysis provides simulated results, while the experimental analysis involves practical testing. The paper demonstrates that the numerical results are consistent with the experimental findings. It is observed that when the PZT transducers are placed at deeper positions within the composites, the peak amplitude of the first wave packet for the A0 mode decreases. Overall, the paper contributes to the understanding of wave propagation in composites by comparing the effects of different positions of PZT transducers, both on the surface and embedded within the material. The findings highlight the influence of transducer placement on wave characteristics, particularly the reduction in peak amplitude with deeper transducer positioning.

Accurate cloud quantification is essential in climate change research. In this work, we construct an automated computer vision framework by synergistically incorporating deep neural networks and finite sector clustering to achieve robust whole sky image-based cloud classification, adaptive segmentation, and recognition under intricate illumination dynamics. A bespoke YOLOv8 architecture attains over 95% categorical precision across four archetypal cloud varieties curated from extensive annual observations (2020) at a Tibetan highland station. Tailor-made segmentation strategies adapted to distinct cloud configurations, allied with illumination-invariant image enhancement algorithms, effectively eliminate solar interference and substantially boost quantitative performance even in illumination-adverse analysis scenarios. Compared with the traditional threshold analysis method, the cloud quantification accuracy calculated within the framework of this paper is significantly improved. Collectively, the methodological innovations provide an advanced solution to markedly escalate cloud quantification precision levels imperative for climate change research, while offering a paradigm for cloud analytics transferable to various meteorological stations.

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.

A novel phase-field method based numerical approach to modeling the compaction process of sedi- ments is presented. Water-logged rock or soil sediments are deposited on the water basins over time and the increasing volume of the sediment compacts under its own weight and external pressure. Coupled evolution of mass conservation, Darcy flow, and the viscoelastic constitutive response, in conjunction with the evolution of porosity and permeability, make this problem highly non-linear and involve moving boundaries. We adopt a phase-field approach to represent the moving sediment interface, and the underlying multiphasic model permits for studying compaction in dynamically evolving sediments. This approach does not necessitate a change of domain sizes as is the case of existing traditional models of sedimentation but rather treats sedimentation growth as a problem of interface motion. We first model a classical compaction problem in 1D to compare with existing results in the literature, and then extend the framework to 2D to model the compaction process taking place under the influence of gravity. The model is then extensively applied to understand the effect of the sediment initial state and the sediment material properties on the compaction process and its spatiotemporal evolution. Lastly, a strong validation of our phase-field treatment results against predictions from traditional methods of solving open boundary sediment compaction problems, and a good agreement between the conventional understanding and the numerical predictions of porosity dependence on the fluid pore pressure for various cases of geological interest, are used to demonstrate the applicability and scalability of this novel numerical framework to model more advanced sediment compaction problems and geometries.

Long-span steel structure trusses are widely used in factory buildings, but with the growth of service time and the increase of dynamic load fatigue, a considerable part of the long-span truss with dynamic load appears serious transverse cracks at the bottom of the middle span and oblique deformation of the abdomen in the operation process. The U-shaped cracks at the bottom and belly, as well as the mid-span down deflection of the main truss, reduced the functional function of the factory building truss structure, and had to limit the original crane load, which affected the normal safety and durability of the structure. Therefore, the application principle of the variable axial force cable system in the long-span factory building truss structure and 3D3S software modeling [1] were used. Analyzing and studying the reinforcement method of large-span powerhouse truss can provide practical experience for subsequent similar projects. In view of the above phenomenon, the large-span powerhouse trusses of Hongcheng Powerhouse 1 and No.2 located in Tonglu, Zhejiang Province are used as the research object, and the variable axial force cable method is proposed to strengthen and lift the load. Considering the span of large-span powerhouse truss, the cable system with 22m controlling force of 400kN is proposed to be selected for powerhouse 1, and the cable system with variable axial force of 24m is proposed to be selected for Powerhouse 2. The force model of large-span truss is established by using the finite element method commonly used to analyze the force of truss. Under two working conditions, The influence of reinforcement effect is analyzed and compared from three aspects: stiffness, bearing capacity and stability. And the phenomenon of uneven stress distribution is analyzed. The stress distribution characteristics of each node are understood by simulating the most disadvantageous node plates with the greatest internal force before and after reinforcement.

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.

The combination of interstitial hyperthermia treatment (IHT) with high dose rate brachytherapy (HDR-BT) for prostate cancer treatment and has the potential to improve clinical outcome, since it highly enhances the efficiency of cell kill, especially when applied simultaneously. Therefore, we have developed the ThermoBrachy applicators. To effectively apply optimal targeted IHT, treatment planning is considered essential. However, treatment planning in IHT is rarely applied since it is regarded difficult to accurately calculate the deposited energy in the tissue in a short enough time for clinical practice. In this study, we investigated various time-efficient methods for fast computation of the electromagnetic (EM) energy deposition resulting from the ThermoBrachy applicators. Initially, we investigated the use of an electro-quasistatic solver. Next, we extended our investigation to the application of geometric simplifications. Furthermore, we investigated the validity of the superpositioning principle, which can enable adaptive treatment plan optimization without the need for continuous recomputation of the EM field. Finally, we evaluated the accuracy of the methods by comparing them to the golden standard Finite-Difference Time-Domain calculation method using gamma-index analysis. The simplifications considerably reduced the computation time needed, improving from >12 h to a few seconds. All investigated methods showed excellent agreement with the golden standard by showing a >99% passing rate with 1%/0.5 mm Dose Difference and Distance-to-Agreement criteria. These results allow the proposed electromagnetic simulation method to be used for fast and accurate adaptive treatment planning.

The splitting method developed in a previous study for linear, multi – component, convection-diffusion mass transfer equations is extended to non-linear, multi-component, convection-diffusion mass transfer equations. The diffusion term is non-linear (the diffusion coefficients depend on chemical species concentrations). Theoretical results about the splitting error are presented, as well as numerical simulations for mixtures with n ≤ 5. The numerical experiments show second-order accuracy and stability for the present splitting algorithm.

The calculation of the collapse load of spherical domes is addressed using a semi-analytical approach under the hypotheses of small displacements and perfect plasticity. The procedure is based on the numerical approximation of the self-stress that represents the projection of the balance equilibrium null space on a finite dimensional manifold. The so obtained self-equilibrated stress span is superimposed to a finite element linear elastic solution to the prescribed loads yielding to the statically admissible set accordingly to Melan’s theorem. The compatibility of the stress with the constitutive law of the material has been enforced using linearized limit domain in terms of generalized stress, namely axial force and bending moment along the local spherical curvilinear coordinates. The procedure has been tested with reference to numerical and experimental data from the literature confirming the accuracy of the proposed method. The comparison with the literature confirms that the buckling load is much greater than the plastic collapse loads both calculated through the proposed procedure and reported in the quoted literature.

Biosensors detect specific bio-analytes by generating a measurable signal from the interaction between the sensing element and the target molecule. Surface acoustic wave (SAW) biosensors offer unique advantages due to their high sensitivity, real-time response capability, and label-free detection. The typical SAW modes are the Rayleigh mode and the shear-horizontal mode. Both present pros and cons for biosensing applications and generally need different substrates and device geometries to be efficiently generated. This study investigates and characterizes ultra-high-frequency (UHF-) SAW resonator biosensors. It reveals the simultaneous presence of the two typical SAW modes, clearly separated in frequency, called slow and fast. The two modes are studied by numerical simulations and biosensing experiments with the glial-fibrillary-acidic-protein (GFAP) biomarker. The slow mode is generally more sensitive to changes in surface properties, such as temperature and mass changes, by a factor of about 1.4 with respect to the fast mode.

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.

Several sun models suggest the radioactive balance where the concentration of greenhouse gases and the albedo effect are related to the Earth's surface temperature. There is a considerable increment of greenhouse gases due to anthropogenic activities. Climate change correlates with this alteration in the atmosphere and an increase in surface temperature. Efficient forecasting of climate change and its impacts of 1.5°C global warming above pre-industrial levels could be helpful to respond to the threat of c.c. and develop sustainably. Many studies have predicted the temperature change in the coming years. The global community has to create a model that can realize good predictions to ensure the best way to deal with the warming. Thus, we propose a finite-time thermodynamic (FTT) approach in the present work. The FTT can solve problems such as the faint young sun paradox. In addition, we use different machine learning models to evaluate our method and compare the experimental prediction and results.

This study investigates the transition from surface to underground mining at the Sepon Gold mine. The stability of surface slopes is assessed prior to commencing underground operations. Stope mining is adopted based on ore body characteristics and geological features. Finite element numerical analysis, employing the Generalized Hoek-Brown criterion, evaluates slope stability for surface slopes and opening stopes. The pit design comprises a 70° slope angle, 20 m height, and 10-15 m safety berm, meeting stability requirements with a factor of safety of 2.46. Underground mining design includes a 62° ore body dip, 50 m thick crown pillar to prevent surface subsidence, and 5 m wide, 5 m high stopes. Sill pillars are left for support after each level excavation. Rock bolts are employed in specific stope areas where is necessary, with continuous monitoring for surface deformation. The study analyzes the influence of stope sizes on pit wall and pit bottom stability, identifying slope failures near the hanging wall close to the pit bottom during underground mining. A slight increase in the displacement zone is observed on the upper crest of the top footwall. Overall, the findings demonstrate successful stability in the transition from surface to underground mining at Sepon Gold mine.

Nickel superalloys occupy a reasonably broad niche in the industry. One of this group's most common and well-known alloys is Inconel 718. Parts from Inconel 718 are used in aerospace, energy, automotive, and some other vital industries, which creates a demand for research and improvement of the machining conditions of this alloy. This article is devoted to the systematization and generalization of the accumulated experience of machining of Inconel 718 at turning operation. The research methodology is based on the Pareto distribution. In the study, more than 50 articles devoted to turning Inconel 718 were analyzed. The tool materials, methods, and types of coating used in turning Inconel 718 were analyzed. Finite element method modeling of processing processes has been considered. The trends in the selection of cutting parameters and the geometry of the cutting tool and finite element method modeling of processing processes have been considered. Lacuna for promising future research was formulated.

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

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.

The discovered light modulation capabilities of diatom silicious valves make them an excellent toolkit for photonic devices and applications. In this work, a reproducible surface-enhanced Raman-scattering (SERS) enhancement was achieved with hybrid substrates employing diatom silica valves coated with an ultrathin uniform gold film. Three structurally different hybrid substrates, based on the valves of three dissimilar diatom species, have been compared to elucidate the structural contribution to SERS enhancement. The comparative analysis of obtained results showed that substrates containing cylindrical Aulacoseira sp. valves achieved the highest enhancement up to 14-fold. Numerical analysis based on the frequency domain finite element method was carried out to supplement the experimental results. Our results demonstrate that diatom valves of different shapes can enhance the SERS signal, offering a toolbox for SERS-based sensors, where the magnitude of the enhancement depends on valve geometry and ultrastructure.

We construct a microscopic model to study discrete randomness in bistable systems coupled to an environment comprising many degrees of freedom. A quartic double well is bilinearly coupled to a finite number N of harmonic oscillators. Solving the time-reversal invariant Hamiltonian equations of motion numerically, we show that for N=1, the system exhibits a transition with increasing coupling strength from integrable to chaotic motion, following the KAM scenario. Raising N to values of the order of 10 and higher, the dynamics crosses over to a quasi-relaxation, approaching either one of the stable equilibria at the two minima of the potential. We corroborate the irreversibility of this relaxation on other characteristic timescales of the system by recording the time dependences of autocorrelation, partial entropy, and the frequency of jumps between the wells as functions of N and other parameters. Preparing the central system in the unstable equilibrium at the top of the barrier and the bath in a random initial state drawn from a Gaussian distribution, symmetric under spatial reflection, we demonstrate that the decision whether to relax into the left or the right well is determined reproducibly by residual asymmetries in the initial positions and momenta of the bath oscillators. This result reconciles the randomness and spontaneous symmetry breaking of the asymptotic state with the conservation of entropy under canonical transformations and the manifest symmetry of potential and initial condition of the bistable system.

Large values and gradients of stresses and deformations, triggering concentrations of stresses and deformations, arise in the corner areas of a structure. The action of forced deformations, leading to the finite rupture of the contact between the elements of a structure, also triggers the concentration of stresses, while the rupture reaches an irregular point, a line on the area boundary. The theoretical analysis of the stress-strain state (SSS) of areas with angular cutouts in the boundary under the action of discontinuous forced deformations is reduced to the study of singular solutions to the homogeneous problem of the elasticity theory that has power-related features. The calculation of stress concentration coefficients in the domain of a singular solution to the elastic problem makes no sense. It is experimentally proven that the zone, that is close to the vertex of the angular cutout in the area boundary, has substantial deformations, rotations, and it corresponds to rising values of the first and second derivatives of displacements along the radius in cases of sufficiently small radii in the neighbourhood of the irregular point of the boundary. For such areas, it is necessary to consider the plane problem of the elasticity theory, taking into account the geometric nonlinearity under the action of forced deformations. This will allow analyzing the effect of relations between orders of values of deformations, rotations, and forced deformations on the form of the equation of equilibrium. The purpose of this work is to analyze the effect of relations of deformation orders, rotations, forced deformations on the form of the equilibrium equation in the polar coordinate system for a V-shaped area under the action of forced temperature-induced deformations with regard for the geometrical non-linearity and physical linearity.

In this paper, we design two inertial-type subgradient extragradient algorithms with line-search process for solving the pseudomonotone variational inequality problems (VIPs) and common fixed-point problem (CFPP) of finite Bregman relatively nonexpansive mapping and a Bregman relatively demicontractive mapping in p-uniformly convex and uniformly smooth Banach spaces, which are more general than Hilbert spaces. Under mild conditions, we derive weak and strong convergence of the suggested algorithms to a common solution of the VIPs and CFPP, respectively. Additionally, an illustrated example is furnished to back up the feasibility and implementability of our proposed methods.

Digital twin, an innovative technology propelled by data and models, plays a seminal role in the digital transformation and intelligent upgrade of ships. This study introduces a digital twin methodology for real-time monitoring of ship structure deformation fields, based on finite discrete strain data, and a visualization tool framework is developed using virtual reality technology. First, the inverse Finite Element Method (iFEM) is employed to derive the deformation field of the ship structure in real time via sensor strain data. Secondly, the deformation field data obtained based on the iFEM algorithm is converted into general visualization data conducive to interpretation within VR applications. Lastly, a digital twin software tool is built to enable synchronous response and interaction between the virtual scene and the physical scene, directly superposing particular virtual objects (data acquired by sensors, CAD virtual models, and deformation field cloud images) onto the physical scene in real-time. The digital twin tool embodies a virtual reality visualization system framework integrating physical data measurement, reconstruction, analysis, expression, storage, rendering and interaction of deformation field data. Through practical application, the flexibility, effectiveness and compatibility of the developed prototype tool are verified. According to the results, the system can enhance the efficiency of scientific communication, model validation, and interdisciplinary sharing during the analysis and evaluation of mechanical properties of ship structures.

We define the notions of strong and strict recurrency for actions of countable ordered groups on $\sigma$-finite non atomic measure spaces with quasi-invariant measures. We show that strong recurrency is equivalent to non existence of weakly wandering sets of positive measure. We also show that for certain p.m.p ergodic actions the system is not strictly recurrent, which shows that strong and strict recurrency are not equivalent.

(1) The article is devoted to the development of a theoretical and algorithmic basis for numerical modeling of the spontaneous potential method (SP) as applied to the study of sandy-argillaceous reservoirs. (2) In terms of coupled flows, we consider a physical-mathematical model of SP signals from an electrochemical source, with regard to the case of fluid-saturated shaly sandstone. (3) An algorithm for 2D finite-element modeling of SP signals was developed and implemented in software, along with its internal and external testing with analytical solutions. The numerical SP modeling was carried out, with determining the dependences on the reservoir thickness and porosity, the amount of argillaceous material and the type of minerals. We performed a comparative analysis of the simulated and field SP data, using the results of laboratory core examinations taken from wells in a number of fields in the Latitudinal Ob Region of Western Siberia. (4) The results of the study may be used either for the development of the existing SP techniques, by providing them with a consistent computational model, or for the design of new experimental approaches.

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.

Alternating current field measurement (ACFM) testing is one of promising techniques in the field of non-destructive testing with advantages of the non-contact capability and the reduction of lift-off effects. In this paper, a novel crack detection approach is proposed to reduce the effect of the angled crack (cack orientation) by using rotated ACFM techniques. The sensor probe is composed of an excitation coil and two receiving coils. Two receiving coils are orthogonally placed in the centre of the excitation coil where the magnetic field is measured. It is found that the change of the x component and the peak value of the z component of the magnetic field when the sensor probe rotates around a crack follows a sine wave shape. A customised accelerated finite element method solver programmed in MATLAB is adopted to simulate the performance of the designed sensor probe which can significantly improve the computation efficiency due to the small crack perturbation. The experiments have also been carried out to validate the simulations. It is found that the ratio between the z and x components of the magnetic field remains stable under various rotation angles. It shows the potential to estimate the depth of the crack from the ratio detected by combining the magnetic fields from both receiving coils (i.e., the x and z components of the magnetic field) using the rotated ACFM technique.

Peri-implant bone resorption has been reported around some implants after loading and this could create problems for the long-term stability of peri-implant soft and hard tissues. The causes are still not completely known, but a relevant importance could be assumed by the presence of a bacterial contamination at the micro-gap level of the implant-abutment junction. In this regard, external and internal implant-abutment assemblies have been shown to be much more permeable to bacterial colonization than Cone-Morse or conical connections. A subcrestal implant placement could have aesthetic advantages allowing a better prosthetic emergence profile. In literature controversial experimental and clinical results have been reported about bone resorption around implants placed equicrestally and subcrestally. Interestingly, Finite Element Analysis (FEA) studies revealed to be extremely useful for assessing the peri-implant bone strain and stress. Thus, the aim of this study was a FEA evaluation of implants with a Cone-Morse implant-abutment assembly (Implacil De Bortoli, São Paulo, Brazil) inserted into a bone block model (width: 10 mm, vertical height: 17.5 mm, thickness of the cortical area: 1.5 mm) mimicking equicrestal and subcrestal placements (1 and 2 mm). The results demonstrated that maximum stresses were observed within the cortical bone around the equicrestally placed implants, the lowest in the implant placed 2-mm subcrestally and intermediate stresses within the implant placed 1-mm subcrestally. Cortical bone was more stressed under lateral loads than axial loads. In conclusion, this FEA study suggested that implant placement at -1 mm could be recommended.

In this work, we incorporate a thixotropic-viscoelastic model into the widely used Computational Fluid Dynamics (CFD) software OpenFOAM, along with the rheoTool library. The model we implement is known as the Modified-Bautista-Manero (MBM), and effectively describes the rheological behavior of worm-like micellar solutions in extensional flows. We provide a detailed explanation of the numerical implementation of the model, specifically using the log-conformation tensor approach. Unlike previous works focused on this kind of fluids, we simulate inertial flows while considering convective terms in the governing equations, thus obtaining a more realistic behavior on the calculated results. The MBM model implementation is validated through numerical simulations on two different industrial-relevant geometries: the planar 4:1 contraction and the 4:1:4 contraction-expansion configurations. Furthermore, we investigate the influence of inertial, viscoelastic, and thixotropic effects on various flow field variables. These variables include velocity, viscosity, normal stresses, and corner vortex size. Our analysis encompasses both transient and steady solutions of corner vortexes across a range of Deborah and Reynolds numbers. Our results are also directly compared with simulations obtained using the non-thixotropic rubber network-based exponential Phan-Thien-Tanner (EPTT) model. From our planar 4:1 contraction results, we found that vortex-enhancement is seen when high elasticity is coupled with quick structural reformation and very low inertial effects. From our planar 4:1:4 contraction-expansion simulations, we show that an increase in inertia leads both to vortex-inhibition in the upstream channel and slight vortex-enhancement in the downstream channel. Lastly, we show the strong effect of the convection of fluidity into the fluidity profiles and into the upstream/downstream corner vortex sizes.

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.

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.

Lamb waves occur in thin-walled structures in two wave modes, the symmetric and antisymmetric mode. Their oscillation on the structures‘ surfaces is either in phase (symmetric) or shifted by a phase angle of π (antisymmetric). In this work, a method is developed to compare the surfaces‘ oscillation phase relation. It is based on the evaluation of time signals regarding the instantaneous phase angle using the continuous wavelet transformation and as a comparative method the short-time Fourier transformation. For this purpose, numerical simulations utilizing the finite element method provide time signals from the top and bottom surface of different thin-walled structures. They differ with respect to their material settings and laminate configurations. The numerically obtained time signals are evaluated by the developed methods. The occurring oscillation phase differences on the top and bottom surface are studied and both methods are compared. Subsequently, the oscillation phase is evaluated experimentally for the wave propagation in a fiber metal laminate. It is shown that the method based on the continuous wavelet transformation is suitable for the evaluation of oscillation phase relations in time signals. Additionally, it is proven that fiber metal laminates show only two phase relations which indicates the occurrence of Lamb waves.

The article presents a simplified method for determining the strength of a corrugated board packaging subjected to dynamic transport loads. The proposed algorithm consists of several calculation steps: (1) a static analysis of the compressive strength of the package, (2) an analysis of random vibrations in the frequency domain, used to determine the resonance frequencies, and (3) a dynamic analysis of the package loaded with computed resonant frequencies. For this purpose, numerical models of the static compression test of the packaging before and after the dynamic analysis of the package subjected to general transport loads were developed. In order to validate the model, the laboratory packaging compression tests were also performed for samples of boxes with 3-layer cardboard. Thanks to this, it was possible to verify the results of numerical simulations of compression tests for several box geometries. Which, in turn, allowed for the development of a method based on dynamic and post-dynamic (static) numerical analysis, allowing to determine with high accuracy the resistance of selected packaging to vibrations and dynamic loads. The results of the (validated experimentally) numerical analysis prove the usefulness of the simplified method presented here for precise estimation of the load capacity of various packages dynamically loaded during transport.

This paper applied the space-time generalized finite difference scheme for solving the nonlinear dispersive shallow water waves as the modified Camassa−Holm equation, modified Degasperis-Procesi equation, Fornberg-Whitham equation, and its modified form. The proposed meshless numerical scheme was composed of the space-time generalized finite difference method, the two-step Newton-Raphson method, and the time-marching method. The space-time approach can treat the temporal derivative as one of the spatial derivatives. This numerical technique enables all the partial derivatives in the governing equation can be discretized by a spatial discretization method, and the mixed derivative can efficiently deal with using the proposed meshless numerical scheme. The space-time generalized finite difference method is advanced from the Taylor series expansion and the moving-least square method. The numerical discretization process is only related to functional data and weighting coefficients on the central node and its nearby nodes. Thus, the matrix system composed of nonlinear algebraic equations will be a sparse matrix and can be efficiently solved by the two-step Newton-Raphson method. Furthermore, the time-marching method was utilized to proceed with the space-time domain along the time axis. In this paper, several numerical examples were tested to verify the capability of the proposed space-time generalized finite difference scheme.

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.

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.

Abstract: This paper outlines the methodology and results for a two species finite volume scalar computational drug transport model developed for simulating the mass transport of Poly(lactic-co-glycolic acid (PLGA) ) from a half-embedded single strut implanted in a coronary arterial vessel wall. The mathematical drug transport model incorporates the convection-diffusion equation in scalar form (dimensionless) with a two species (free-drug and bound-drug) mass transport setup including reversible equilibrium reaction source terms for the free and bound-drug states to account for the pharmaco-kinetic reactions in the arterial wall. The relative reaction rates of the added source terms control the interconversion of drug between the free and bound-drug states. The model is solved by a 2-D finite-volume method for discretizing and solving the free and bound drug transport equations with anisotropic vascular drug diffusivities. This model is an improvement over our previously developed model using the finite-difference method. A dimensionless characteristic scaling pre-analysis was conducted a priori to evaluate the significance of implementing the reaction source terms in the transport equations. This paper reports the findings of an investigation of the interstitial flow profile into the arterial wall and the free and bound drug diffusion profiles with a parametric study of varying polymer drug concentration (low and high), tortuosity, porosity, and Peclet and DamKöhler numbers over the course of 400 hours (16.67 days). The results also reveal how a single species drug delivery model that neglects both a reversible binding reaction source term and the porosity and tortuosity of the arterial wall cannot accurately predict the distribution of both the free and bound drug.

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

Microstructure evolution in metal additive manufacturing (AM) is a complex multi-physics and multi-scale problem. Understanding the impact of AM process conditions on the microstructure evolution and the resulting mechanical properties of the printed part is an active area of research. At the meltpool scale, the thermo-fluidic governing equations have been extensively modeled in the literature to understand the meltpool conditions and the thermal gradients in its vicinity. In many phenomena governed by partial differential equations, dimensional analysis and identification of important dimensionless numbers can provide significant insights into the process dynamics. In this context, a novel strategy using dimensional analysis and the method of linear least squares regression to numerically investigate the thermo-fluidic governing equations of the Laser Powder Bed Fusion AM process is presented in this work. First, the governing equations are solved using the Finite Element Method, and the model predictions are validated by comparing with experimentally estimated cooling rates, and with numerical results from the literature. Then, through dimensional analysis, an important dimensionless quantity - interpreted as a measure of heat ab- sorbed by the powdered material and the meltpool, is identified. This dimensionless measure of heat absorbed, along with classical dimensionless quantities such as Peclet, Marangoni, and Stefan numbers, is used to investigate advective transport in the meltpool for different alloys. Further, the framework is used to study the variations of thermal gradients and the solidification cooling rate. Important correlations linking meltpool morphology and microstructure evolution related variables with classical dimensionless numbers are the key contribution of this work.

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.

Nitrile butadiene rubber (NBR) and its various composites are extensively studied as matrix materials for oil-free lubricated mechanical friction surfaces. While numerous macroscopic friction experiments can compare the frictional performance of different NBR composites, the complex material preparation and low experimental efficiency pose challenges. Nanoscale molecular dynamics (MD) simulations can help develop materials with improved friction properties, but they cannot directly provide the friction coefficients of actual macroscopic friction surfaces. Therefore, it is important to integrate the advantages of macroscopic and nanoscale friction studies and perform a synergistic analysis to modify the friction properties of composite materials. In this study, a multi-scale approach is proposed to simulate the frictional characteristics of carbon nanotube (CNT)-reinforced NBR by combining the MD friction simulation， micromechanical bridging model, and macroscopic finite element analysis (FEA) methods. The results of the multi-scale friction simulation of copper (Cu)-CNT/NBR composites show that the addition of CNTs significantly improves the frictional properties of NBR, and the extent of the improvement is related to the mass fraction of CNTs. As the mass fraction of CNTs increases (0%, 1.25%, 2.5%, 5%), the coefficient of friction decreases from 0.50 to 0.38. The results are validated by ring-block friction experiments at the same scale as the simulations.

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.

Cross-field innovation in faba bean breeding has added significant value to the traditional faba bean planting industry by allowing for the selection of faba bean varieties with a dual purpose as ornaments and food. In order to improve the selection efficiency, we determined 17 phenotypes related to ‘dual purpose’ faba beans and evaluated 37 new breeding lines, which were bred for 10 years by our team, and two traditional faba bean varieties by means of an Analytic Hierarchy Process. The lines and varieties were evaluated and sorted into five grades, from which two breeding lines with excellent comprehensive phenotypes and dual ornamental and food functions were selected. The new line GS5 passed the identification of Crop Identification Committee of Chongqing Seed Industry Association on November 25, 2021 and was renamed "Doumei 1 hao" with the identification number of " 202111". The results confirmed that the comprehensive evaluation system constructed in this study can be used for the effective selection of dual-purpose faba bean lines, providing an empirical approach for strategic decision-making by breeders.

Traumatic Brain Injury represents a significant public health concern due to its role in traumatic death and disability, often caused by head impacts or rapid accelerations. Despite advancements in Finite Element Head Models, existing models primarily represent adult males, neglecting anatomical and physiological differences between genders and age groups. Additionally, the mechanistic roles of specific intracranial tissues, such as the dura mater, have conventionally been overlooked. To better understand the mechanical behavior of the brain and its injured regions in the context of TBI, including the dura mater, this study proposes the development of a Female Finite Element Head Model. This model, representative of a middle-aged female subject, was developed using medical image-derived geometry and finite element modeling techniques. The validation results demonstrate a similarity between the numerical displacement curves and those obtained experimentally. In the model with dura mater, the numerical results show analogous behavior to the experimental results, despite minimal variations in the amplitude of the curves. These advancements highlight the significance of including the dura mater in biomechanical brain modeling. Incorporating the dura mater in Finite Element Head Models produces more realistic results than models without it. This emphasizes the importance of enhancing biomechanical models to better represent anatomical complexity.

In this study, a new structure is proposed based on the body-centered cubic (BCC) lattice structure: adding a cubic truss in the center of the BCC structure, which is denoted by TLC (Truss-lattice-cube). The different dimensions of the central cube can affect the mechanical properties of the lattice structure notably. With a fixed length (15mm) of a unit cell, the optimal size for the central cube is determined to be 5 mm. Quasi-static compressive tests are performed on specimens made of Polylactic Acid (PLA) using additive manufacturing technology. The deformation characteristics of the new structure is analyzed in detail by experiments and numerical simulations. Compared to the BCC structure, the mechanical properties of the TLC structure are significantly improved. The initial flow stress of the TLC has been increased by 122% at a strain of 0.1; the specific strength has been enhanced by 293% at a strain of 0.5; and the specific energy absorption has been improved by 312% at the strain of 0.6. Printing defects in the lattice structure could deteriorate the mechanical property remarkably. In this work, incorporation of microcracks into the finite element model makes the simulation capture the influence of printing defects and improve the predictive accuracy of simulation significantly.

This paper introduces an innovative approach to numerically model Structure-Soil-Structure Interaction (SSSI) by integrating the Boundary Element Method (BEM) and the Finite Element Method (FEM) in a coupled manner. To assess the accuracy of the proposed method, a comparative study is undertaken, comparing its outcomes with those generated by the conventional FEM technique. Alongside accuracy, the computational efficiency aspect is crucial for the analysis of large-scale SSSI problems. Hence, the computational performance of the coupled BEM-FEM method undergoes a thorough examination and is compared with that of the standalone FEM method. The results from these comparisons illustrate the superior capabilities of the proposed method in comparison to the FEM method. The novel approach provides more reliable results compared to traditional FEM methods, serving as a valuable tool for engineers and researchers involved in structural analysis and design.

Recent studies have shown that the mechanical properties of bioinspired periodic composite materials can be strongly influenced by finite deformation effects, leading to highly nonlinear static and dynamic behaviors at multiple length scales. For instance, in porous periodic nacre-like microstructures, microscopic and macroscopic instabilities may occur for a given uniaxial loading process and, as a consequence, wave attenuation properties may evolve as a function of the microstructural evolution, giving it the designation of metamaterials. The numerical outcomes provide new opportunities to design bioinspired soft composite metamaterials characterized by high deformability and enhanced elastic wave attenuation capabilities given by the insertion of voids and lead cores.

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.

The inverse finite element method (iFEM) is a novel method for reconstructing the full-field displacement of structures by discrete measurement strain. In practical engineering applications, the accuracy of iFEM is reduced due to the positional offset of strain sensors during installation, and errors in structural installation. Therefore, coarse and fine two-stage calibration (CFTSC) method is proposed to enhance the accuracy of the reconstruction of structures. Firstly, the coarse calibration is based on a single-objective particle swarm optimization algorithm to optimize the displacement-strain transformation matrix related to the sensor position. Secondly, as selecting different training data can affect the training effect of self-constructed fuzzy networks, this paper proposes to screen the appropriate training data based on residual analysis. Finally, the experiments of the wing integrated antenna structure verify the efficient of the method on the reconstruction accuracy of the structural body displacement field.

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.

This work aims to describe the dielectric behaviours of four main tissues in the human forearm using mathematical modelling, including fat, muscle, blood, and bone. Determining the electrical response of multiple tissues can help develop physiological monitoring of an organ or a body’s section by filtering out the impedance contributions from the surrounding tissues to blood flow-induced impedance variations. Multi-frequency bioimpedance analysis (MF-BIA) was initially performed using the finite element method (FEM) with a 3D forearm model followed by a pilot study to characterise the response of actual forearm tissues from 1 kHz to 349 kHz. Both the simulation and experimental results were fitted to a single-dispersion Cole model (SDCM) and a multi-dispersion Cole model (MDCM) to determine the Cole parameters for each tissue. The correlation analysis and Bland-Altman plot indicated a good fit between raw and fitted impedance values using both SDCM and MDCM. Overall, MDCM exhibited better performance in fitting and estimation of the Cole parameters with correlation coefficient (R2) of 0.99 and 0.97, root mean squared error (RMSE) of 0.09 Ω and 0.14 Ω, and mean difference (mean±standard deviation) of 0.00±0.09 Ω and -0.03±0.14 Ω for the real part and imaginary part of impedance, respectively.

Theoretical approaches were applied to study the effect of magnetic field and heat transfer on the flow of blood plasma through an asymmetric arterial segment. The plasma was considered to be unsteady, laminar and an incompressible fluid through non-uniform arterial segment in a two-dimensional flow. Axial velocity, pressure gradient, stream function and pressure shift per wavelength were evaluated in blood plasma flow in arteries by use of coupled linear partial differential equations, solved with the help of Finite difference method. The corresponding initial and boundary conditions were obtained by discretizing the flow channel. The effect of magnetic field on blood plasma in arteries was determined by creating a magnetic field gradient through the application of a varying strength of magnetic field on the flow. Heat transfer characteristics due to applied magnetic field and viscosity of blood was obtained by use of linear partial differential equations to determine varying temperature conditions in heat transfer characteristics. Numerical results for velocity profiles, magnetic profiles and temperature profiles were obtained to characterize blood plasma viscosity by use of various non-dimensional parameters. The results were graphically presented using MATLAB software. The results obtained helped in analyzing theoretically the effects of magnetic field and heat transfer in arterial plasma flow.

The aim of this paper is to optimize thickness variation function of simply supported and cantilever beams, in terms of maximizing gaps between chosen neighboring frequencies and analyze obtained results. The optimization results are examined in terms of achieving the objective function (related to eigenproblem), but also in terms of their dynamics stiffness (forced vibrations excited by point harmonic load). In the optimization process, a Genetic Algorithm was used. Problems related to structural dynamics were solved by FEM implementation into the algorithm. Sample results were presented, and the developed algorithm was analyzed in terms of result convergence by examining several variable parameters. The authors demonstrated the validity of applying the described optimization tool to the presented problems. Conclusions were drawn regarding the correlation between stiffness and mass distribution in the optimized beams and the natural frequency modes in terms of which they were optimized.

This study explores the significance of active suspension systems for vehicles with lower chassis compared to conventional ones, aiming at the development of future automobiles. Conventional linear MR (Magneto-Rheological) dampers were found inadequate in ensuring sufficient vibration control because the vehicle's chassis becomes lowered in the unmanned vehicles or purposed based vehicles. As an alternative, a rotary type of MR damper is proposed in this work. The proposed damper is designed based on pre-specified design parameters through mathematical modeling and magnetic field analyses. Subsequently, a prototype of the rotary MR damper identical to the design is fabricated, and effectiveness is shown through experimental investigations. In configuring the experiments, a proportional-integral (PI) controller is employed for current control to reduce the response time of the damper. The results presented in this work provides useful guidelines to develop a new type of MR damper applicable to various types of future vehicles’ suspension systems with low distance from the tire to the body floor.

Bio-inspired gyroid triply periodic minimum surface (TPMS) lattice structures have been the focus of research in automotive engineering because they can absorb a lot of energy and have wider plateau ranges. The main challenge is determining the optimal energy absorption capacity and accurately capturing plastic plateau areas using finite element analysis (FEA). Using nTop's Boolean subtraction method, this study combined walled TPMS gyroid structures with normal TPMS gyroid lattice. This made a composite TPMS-gyroid lattice with relative densities ranging from 14% to 54%. Using the ideaMaker software and the fused deposition modeling (FDM) Raise3D Pro 2 3D printer to print polylactic acid (PLA) bioplastics in 1.75 mm filament made it possible to slice computer aided design (CAD) models and fabricate 36 lattices samples precisely layer by layer technique. Shimadzu 100kN testing equipment was utilized for the mechanical compression experiments. And the validation using finite element analysis (FEA). Further, CTG was examined using a field emission scanning electron microscope (FE-SEM) before and after compression testing. The composite TPMS gyroid lattice showed potential as shock absorbers for vehicles with relative densities of 33%, 38%, and 54%. The Gibson-Ashby model showed that the composite TPMS gyroid lattice deformed mainly by bending, and the size effect was seen when the relative densities were less than 15%. The lattice's relative density had a significant impact on its ability to absorb energy. The research also explored the use of these innovative foam-like composite TPMS-gyroid lattices in high-speed crash box scenarios, potentially enhancing vehicle safety and performance.

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 presented research article shows a functionally graded lattice as an infill structure in designing a 3D printed mechanical element—a bicycle crank arm. The authors want to answer the fundamental question: whether the functionally graded lattice structures can efficiently improve structural performance and how they can be implemented in real complex components. Intuitively, everyone senses that using such elements must bring profit. A good example might be nature, which commonly reaches for such a solution in the skeletons of animals and humans. Two aspects determine their realisations: the lack of adequate construction and analysis methods and the limitations of existing manufacturing methods. Thus, the authors used relatively simple crank arm design and design exploration methods for structural analysis. The design exploration method enabled us to find the optimum solution efficiently. A developed prototype was created utilising additive manufacturing. One of the 3D printing methods - fused filament fabrication technology for metals, allowed the prototype of a crank arm with the optimised infill to be made. As a result, the authors developed a lightweight and manufacturable crank arm showing the new construction and analysis method implementable in similar 3D printed elements.

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

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.

The design of water Vessel fender system can minimize impact energy, which could result in severe vessel or offshore structure damage during collision. In order to protect Offshore structures and vessels from damage after impact, a Fender design approach is proposed which uses Soft materials like rubber, sand, inflated rubber tube, etc. to reduce damage on impacted vessels as well as offshore structures like wind turbines, bridges, oil drilling platforms in water etc. The fender is designed in such a way that it can handle vessels made of different material, sizes, and weights.The proposed solution is a fender consisting of a number of layers of soft sandwich materials, each layer getting softer as it moves away from the structure i.e. Hard-Soft fender design which can absorb impact energy upon collision. To handle lighter vessels made of soft or brittle materials, the outer layer is made of inflatable or Very soft foam rubber; however, the inner layer covering the offshore structure is made of hard material to handle heavy vessels made of hard material such as steel. Furthermore, the use of sand as a natural fender along waterways such as narrow channels is also proposed to reduce the cost of fender installation and catastrophic vessel impact damage on rocky ground since it produces less deformation just like rubber fender as demonstrated in the simulation. Finite Element Analysis using ANSYS software, Explicit Dynamic Autodyn, was performed. The results showed a considerable decrease in the total deformation and internal energy of the vessel and fended structure when collision occurs. Furthermore, the result showed that Hard-Soft fender layer performs much better than Soft-Hard-Soft followed by Soft-Hard Fender layer design.

In this study, the autofrettage method, which is used to increase the durability of cylindrical parts, was applied to a barrel made of steel material according to the AISI 4340 standard. Assuming that the barrel works continuously under the pressure of 400 MPa, optimum autofrettage was performed. The stresses occurring on the barrel before and after autofretage were calculated. Additionally, analytical equivalent stress values (Von-Mises) formed on the part's wall were compared with numerical analysis results obtained with a commercial finite element program ANSYS®, and their compatibility was examined. As a result of the calculations and analysis performed, it was observed that autofretage has positive effect on increasing the pressure carrying capacity of the barrel. After autofrettage the cylindrical part, it was determined that the Von-Mises stress values formed in the part wall changed and these stresses took the lowest value at a point (elastic-plastic interface). It has been evaluated that this point gives the optimal autofretage pressure value for a certain operating pressure. After the autofrettage application, it was determined that the Von-Mises Stress value in the barrel wall decreased by 25.27 % and the pressure carrying capacity of the barrel increased without any design or dimensional change. It has been seen that the analytical and numerical analysis results are consistent with each other, and it is expected that this study will guide and contribute to the studies to be conducted on increasing the strength or lightening the weight of the parts.

For enhancing the efficiency of cantilever casting the construction in the bridge, a novel rhombus traveling track of hanging basket was developed and optimized. The civil software titled Midas was applied to analyse the mechanical properties of the hanging basket with full load and control load. The strength and distortion of the walking mechanism were within the specified range when the maximum beam unit internal force of 149.69Mpa, which was less than the allowable stress of steel 175Mpa, and the anti-overturning safety factor of the hanging basket was 2.5 to meet the requirements. Through the comparative analysis of the key components and the finite element calculation, it was found that there was about 30% redundancy in the structure performance. Therefore, further optimization of each structure was carried out, and the front elevation of the main truss of 20° was achieved. It obtained the best performance in all aspects.

Introducing the Moving Least-Squares Aided Finite Element Method (MLS-FEM), we present a robust computational technique that effectively evaluates the development of mixing in polymeric mixing devices. Our approach uses the Moving Least-Squares (MLS) interpolant to enhance the classic finite element shape functions. As a result, it becomes possible to treat fixed or moving components, such as rotors of a mixer, as discontinuities within the flow domain. This allows analysis of the flow field with a background mesh independent of the rotors' shape, orientation, and position. We compared the MLS-FEM results with the classical boundary-fitted mesh FEM, which comprises the mesh regeneration steps by the movement of the rotors. We then showed that the new method provides accurate results for the flow field parameters (i.e., the velocity components and the pressure values) and the mixing evolution without requiring mesh regeneration and without compromising accuracy.

Cretan lyra is a stringed instrument very popular on the island of Crete, Greece and an important part of its musical tradition. For stringed musical instruments, the air mode resonance plays a vital part in their sound, especially in the low frequency range. For this study, the air mode resonance of a Cretan lyra is investigated with the use of finite element method (FEM). Two different FEM acoustic models were utilized: first, a pressure acoustics model with the Cretan lyra body treated as rigid was used to provide an approximate result; secondly, an acoustic-structure interaction model was applied for a more accurate representation. In addition, acoustic measurements were performed to identify the air mode resonance frequency. Results of this study reveal that the acoustic-structure interaction model has a 3.7% difference regarding the actual measurements of the resonance frequency. In contrast, the pressure acoustics solution is approximately 13.8% too high compared with the actual measurements. Taken together, findings of this study support the idea that utilizing FEM acoustic-structure interaction models could possibly predict more accurately the vibroacoustic behavior of musical instruments, which in turn can enable the determination of key aspects that can be used to control the instrument’s tone and sound quality.

This paper proposes an analytical model of the double-stator spoke-type permanent-magnet vernier machine (DSSTVM) using the subdomain model (SDM), which can be used to calculate the mag-netic field distribution and corresponding electromagnetic parameters of the DSSTVM. The whole field domain is divided into several subdomains according to the magnetic characteristics of each region, within which Laplace's and Poisson's equations are solved accordingly in terms of magnetic vector potential (MVP). Then, the magnetic flux density distribution, back electromotive force (EMF), and electromagnetic torque of the DSSTVM can be obtained correspondingly. Ultimately, finite element analysis (FEA) is adopted to validate the proposed analytical model's effectiveness for quickly predicting the no-load and on-load performances of the DSSTVM.

Nickel-based Superalloy Inconel 718 is widely used in the aerospace industry for its excellent high-temperature strength and thermal stability. However, milling Inconel 718 presents challenges due to significantly increased cutting force and vibration, which is a typical difficult-to-machine material. This paper focuses on the milling process of Inconel 718, establishing a milling force model to analyze the force trends under various processing parameters. Finite element analysis is employed to study the stress and temperature fields during milling. Dynamic equations for milling Inconel 718 are developed, and stability lobe diagrams are generated based on modal experiments. Milling experiments on Inconel 718 validate the milling force model and finite element analysis results. The fmincon optimization algorithm is utilized to identify the optimal machining parameters for Inconel 718. Through this research, valuable insights into enhancing the efficiency and quality of Inconel 718 machining are provided. This study contributes to a deeper understanding of Inconel 718 milling behavior, offering crucial guidance for more efficient machining processes.

One of the uses of Helmholtz resonators is as sound absorbers for room acoustic applications, especially for the low frequency range. Their efficiency is centered around their resonance frequency which is mainly depended on elements of their geometry such as the resonator volume and neck dimensions. Incorporating additional necks on the body of a Helmholtz resonator (depending on whether they are open or closed) has been found to alter the resulted resonance frequency. For this study, tunable Helmholtz resonators to multiple resonance frequencies, are proposed and investigated utilizing additional necks. The resonance frequencies of various multi-neck Helmholtz resonators are first modeled with the use of finite element method (FEM), then calculated with the use of an analytical approach and the results of the two approaches are finally compared. Results of this study show that Helmholtz resonators with multiple resonances at desired frequencies are achievable with the use of additional necks, while FEM and analytical methods can be used for the estimation of the resonance frequencies. Analytical and FEM approach results show a good agreement in cases of small number of additional necks, while the increasing differences in cases of higher neck additions, were attributed to the change in effective length of the necks as demonstrated by FEM. The proposed approach can be useful for tunable sound absorbers for room acoustics applications according to the needs of a space. Also, this approach can be applied in cases of additional tunable air resonances of acoustic instruments (e.g. string instruments).

With simple but rigorous arguments we prove that the ordinary second quantization of bosonic and fermionc fields is formally equivalent to constraining the elementary degrees of freedom of the classical fields to have intrinsically periodic dynamics in time. This result confirms and extends the formal equivalence obtained in previous papers for general Hamiltonian systems in terms of canonical quantization, for elementary relativistic particles also in terms of the Feynman Path Integrals, as well as other remarkable correspondences of both phenomenological and theoretical fundamental interest.

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.

A model of a single-server queuing-inventory system (QIS) with catastrophes and a finite waiting room for consumer customers (c-customers) has been developed. When a catastrophe occurs, all inventory in the system's warehouse is destroyed, but c-customers in the system are still waiting for restocking. In addition to c-customers, negative customers (n-customers) are also taken into account, which displaces one c-customers (if any). The policy (?, ?) is used to replenish stocks. If upon arrival of the c-customer the inventory level is zero, then according to Bernoulli’s scheme this c-customer is either lost or joining the queue. The mathematical model of the studied QIS is constructed in the form of a two-dimensional continuous-time Markov chain (2D CTMC). Exact and approximate methods for calculating the steady-state probabilities of constructed 2D CTMCs are developed and closed-form expressions are derived for calculating the performance measures. The results of numerical experiments are presented, demonstrating the high accuracy of the developed approximate formulas, as well as the behavior of performance indicators depending on the reorder point. In addition, an optimization problem is solved to obtain the optimal reorder point to minimize the expected total cost.

The ultrasonic vibration assisted pressing process can improve the fluidity, the uneven distribution of density and particle size of WC-Co powder. However, the microscopic mechanism of ultrasonic vibration on the powder remains unclear. In this paper, WC and Co particles are simulated by three-dimensional spherical models of two particle sizes with the aide of the secondary development of Python. At the same time, the forming process of powder at mesoscale is simulated in the virtue of the finite element analysis software ABAQUS. The influence of vibration amplitude on the fluidity, the filling density and the stress distribution of WC-Co powder is investigated when the ultrasonic vibration is applied to the conventional pressing process. From simulation results, the ultrasonic vibration amplitude has great influences on the density of the compact. With the increase of the ultrasonic amplitude, the compact density also increases gradually. At the initial stage of the compaction, the particles move violently under the action of ultrasonic vibration, the fluidity is obviously enhanced, the arch bridge effect between particles is destroyed, the fine particles quickly fill the pores between the large particles, which increase the density, and the residual stress in the billet decreases after the compaction. From the experimental results, the size distribution of the billet is more uniform, the elastic aftereffect is reduced, the dimensional instability is improved, and the density curves obtained by experiment and simulation are within a reasonable error range.

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.

The tension cable-supported power transmission structure (TC-PTS) is a new type of power transmission structure suitable for mountainous terrain, which is sensitive to wind load. In this regard, a nonlinear finite element analysis model of wind-induced vibration is proposed for the TC-PTS, and the wind-induced vibration response of the structure is analyzed. Firstly, the tangent stiffness matrix of the three-dimensional truss element for the supporting suspension cable and transmission line, considering the geometric nonlinearity of structures, is derived through the relationship between the element elastic energy and its displacement. Subsequently, the element mass matrix and damping matrix of the supporting suspension cable and transmission line, as well as the element nodal load vector obtained from wind load equivalence are given. Then, based on the nonlinear finite element theory, the nonlinear dynamic equation of wind-induced vibration is established for the TC-PTS and solved by Newmark-β method combined with Newton-Raphson iterative method. Furthermore, the rain-flow counting method and Miner's linear fatigue cumulative damage theory were used for wind induced fatigue damage assessment. Finally, a two-span TC-PTS is selected as an example, and the wind-induced nonlinear vibration and fatigue damage assessment are analyzed through the proposed model. The results show that the proposed model has high computational accuracy and efficiency. With the increase of wind speed and wind direction angle, the maximum lateral displacement and tension of the supporting suspension cable and transmission line increase, and their degree of increase shows a nonlinear trend. In terms of the wind-induced fatigue analysis results of TC-PTS, the fatigue damage at the end of the supporting-conductor suspension cable is greater than the fatigue damage at its midpoint. Compared to the fatigue damage at the midpoint of the conductor, the fatigue damage at the end of the conductor is less affected by wind direction angle, and both are more significantly affected by the wind speed.

Transcutaneous nerve stimulation has several neuromodulatory applications, including pain management, sensory restoration, and motor recovery. Transcutaneous stimulation delivered via electrode arrays offers ease of use and reconfigurability through non-invasive means. Also, employing electrode arrays allows regulation of the permeated electric fields for greater control over the stimulation performance. Field shaping involves the regulation of permeated electric fields and is a well-established practice for deep brain stimulation; however, its applicability to peripheral nerve stimulation still needs exploration. Adopting such techniques can help to overcome current limitations with transcutaneous stimulation, like poor targetability, early onset of fatigue, and perceived discomfort. This article provides a perspective on the applicability of field-shaping techniques in improving the overall stimulation outcome.

The paper presents the numerical and experimental results of research aimed at determining the influence of hardness in the range of 50-90 Shore A of layered tools made of elastomeric materials on the possibility of forming Inconel 625 nickel-based alloy sheets. A stamping die made of 90MnCrV8 steel (hardness 60HRC) was designed for forming embosses in drawpieces, ensuring various stress states on the cross-section of the formed element. The principle of operating the stamping die was based on the Guerin method. Finite element-based numerical modelling of the forming process for various configurations of polyurethane inserts was also carried out. The drawpieces obtained by sheet forming were subjected to geometry tests using optical 3D scanning. The results confirm that in the case of forming difficult-to-deform Inconel 625 Ni-based alloy sheets, the hardness of the polyurethane inserts significantly affects the geometric quality of the obtained drawpieces. The assumptions made in numerical simulations were verified in experimental studies. Based on the test results, it was concluded that the selection of polyurethane hardness should be determined by the shape of the formed element. Significant non-uniform sheet metal deformation was also found, which may pose a problem in the process of designing forming tools and technology of plastic forming of Inconel 625 Ni-based alloy sheets.

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

Fatigue failure remains a critical concern in structural engineering and material science, prompting extensive research to understand and predict the behaviour of materials under cyclic loading conditions. The present study aims to investigate the fatigue life of carbon steel specimens containing opposite semicircular edge notches through a comprehensive experimental and numerical analysis. In this study, stress concentration factors (SCF, K_t) of rectangular plate with opposite semicircular notches are considered under uniform tensile stress to analyse the notch deformation because of stretching of plate. Furthermore, the research focuses on quantifying stress concentration factors (SCFs) for these notches based on S-N curves of carbon steel. The study employs a combination of experimental and numerical techniques to understand the influence of these notches on the fatigue performance of carbon steel structures. A plate with opposite semicircular edge single notches under the axial load creates stress concentration near the notch and it is much larger than the average stress on the plate. Both analytical and finite element methods are used to calculate the maximum stress around the notch. SOLIDWORKS Premium Student Edition 2023 has been employed for modelling and SOLIDWORKS Simulation Premium Student Edition 2023 has been used for stress analysis and fatigue notch factor of rectangular plate of size 31 mm x 25.4 mm x 6.35 mm. The uniform tensile load with a magnitude of 20195 N is applied on one sides of rectangular plate normal to the sides of notches with ratio h/r = 1, for the semicircular notch. The result obtained on both analytical and finite element methods are compared and the percentage of error has been evaluated. Subsequently, these specimens undergo fatigue testing under varying loading conditions to capture their fatigue behaviour. The acquired fatigue data is then plotted against stress amplitude to construct S-N curves, forming the foundation for assessing the fatigue life of the notched specimens. To complement the experimental findings and to gain a deeper understanding of the stress concentration phenomenon, numerical simulations are conducted using advanced finite element analysis (FEA) techniques. The finite element models are carefully calibrated against the experimental results to ensure their accuracy and reliability. The FEA simulations enable the determination of stress concentration factors at critical locations within the notched specimens, further validating the experimental observations. The investigation reveals crucial insights into the effect of opposite semicircular edge notches on the fatigue life of carbon steel structures. The obtained S-N curves allow engineers and designers to predict the fatigue life of components with similar notches, aiding in the development of reliable and durable structures in practical applications. Moreover, the stress concentration factors determined from the numerical simulations provide valuable data to assess the potential failure modes and to optimise designs, effectively mitigating fatigue-related failures. The combination of experimental and numerical approaches ensures a comprehensive and rigorous analysis of the fatigue behaviour in notched specimens, offering a reliable basis for making informed engineering decisions. The comparison between the analytical method and the Finite Element Method (FEM) demonstrated good agreement, with an error percentage of 4.272%. The analysis revealed that the specimen would experience failure after approximately 2882 cycles, with a maximum stress of 395.914 MPa. This research study enhances the understanding of fatigue life in carbon steel structures containing opposite semicircular edge notches and contributes valuable data to the field of fatigue mechanics. The outcomes serve as a valuable resource for professionals engaged in structural engineering, material science, and design optimisation, ultimately leading to safer and more durable industrial components in critical applications. The findings of this research contribute to the understanding of fatigue behaviour in carbon steel components with stress concentration effects caused by semicircular notches. Moreover, the validated numerical simulations and data curves facilitate the prediction of fatigue life and aid in determining the critical conditions leading to fatigue failure. In conclusion, this research highlights the significance of combining experimental testing with numerical simulations to comprehensively analyse the fatigue life of carbon steel specimens with opposite semicircular edge notches. The obtained stress concentration factors provide crucial information for structural integrity assessment and offer potential for further optimising design criteria to mitigate fatigue-related failures. The findings of this study could play a vital role in enhancing the reliability and safety of carbon steel structures subjected to cyclic loading conditions. The comprehensive experimental and numerical analyses establish a foundation for future studies exploring other materials and geometries with notches, fostering advancements in fatigue life prediction and structural integrity assessment.

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.

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.

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.

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.

In the field of computational fluid dynamics (CFD) for rotorcraft, two significant challenges are resolving the complex vortex structures in rotor wakes, and representing the moving rotor blades in the ambient airflow. This study addresses the first challenge by using a third-order unstructured finite volume solver, which, compared to its second-order counterpart, has substantially lower dissipation. Consequently, even relatively coarse meshes are capable of resolving small vortices. With this background flow field solver, the second challenge is addressed by modeling each rotor as an Actuator Disk Model (ADM), or describing individual rotor blades as actuator lines, designated as the Actuator Line Model (ALM). Both of these models are equipped with an improved correction for aerodynamic losses at blade tips, which is thoroughly presented in the methodology section. The numerical experiments section centers on analyzing errors linked to various sampling approaches. Additionally, the article discusses comparisons between vortex theory and ALM, specifically regarding calculations for fixed-wing aircraft. Furthermore, high-order precision and parallel efficiency are exemplified in scenarios encompassing rotors engaged in both hovering and forward flight rotors. The results in this paper demonstrate that the combination of the ALM/ADM with the improved tip loss correction and the third-order finite volume solver presents a new way of developing efficient tools for the aerodynamic analysis of helicopter rotors.

Compared with ordinary asphalt pavement, the application of modified asphalt will improve the performance of asphalt pavement. However, using the same modification method for the whole road section can cause modification waste or result in some parts not meeting performance standards. This paper proposes a whole road section modification method. First of all, a finite element model was established using Abqus according to the actual road pavement material to analyze the mechanical response of each layer under different working conditions and establish the grading standard under comprehensive working conditions. Secondly, different modification methods were analyzed according to investigation and their modification effects were graded through clustering analysis method. Finally, different working conditions and different modification methods were reasonably matched and a section of highway was selected for verification. The results show that before modification, the calculated rutting is 6.23mm and after modification, it is reduced to 2.45mm which indicates that reasonable modification has a better effect and can achieve the purpose of rutting equalization control on the whole road section. It is of great significance to carry out reasonable modification on different conditions to control rutting equalization of asphalt pavements.

We present a modified characteristic finite element method that exhibits second-order spatial accuracy for solving convection-reaction-diffusion equations on surfaces. The temporal direction adopted backward-Euler method, while the spatial direction employed surface finite element method. In contrast to regular domains, it is observed that the point in the characteristic direction traverses the surface only once within a brief time. Thus, the good approximation of the solution in the characteristic direction holds significant importance for the numerical scheme. In this regard, Taylor expansion is employed to reconstruct the solution beyond the surface in the characteristic direction. The stability of our scheme is then proved. A comparison is made with an existing characteristic finite element method based on face mesh. Numerical examples are provided to validate the effectiveness of our proposed method.

Scanning underwater areas, using magnetometers, in search of unexploded ordnance is a difficult challenge, where machine learning methods can find a significant application. However, this requires the creation of a set enabling the training of prediction models. Such a task is difficult and costly due to the limited availability of relevant data. To meet this challenge in the article, we propose the use of numerical modeling to solve this task. The conducted experiments allow us to conclude that it is possible to obtain high compliance of the numerical model with the results of physical tests. In addition, the paper discusses the methodology of simplifying the computational model, allowing for almost three times reduction of the calculation time. In addition, in the work we present the methodology of creating an appropriate data set, enabling the generation of any number of training samples.

Due to the ongoing global warming on the Earth, permafrost degradation has been extensively taking place, which poses a substantial threat to civil and industrial facilities and infrastructure elements, as well as to the utilization of natural resources in the Arctic and high-latitude regions. In order to prevent the negative consequences of permafrost thawing under the foundations of constructions, various geophysical techniques for monitoring permafrost have been proposed and applied so far: temperature, electrical, seismic and many others. At the same time, non-stationary electromagnetic methods seem to have significant potential and a number of important practical advantages. We propose a cross-borehole exploration system for the highest localization of target objects in the cryolithozone. A novel mathematical apparatus for three-dimensional modeling of transient electromagnetic (TEM) signals by the vector finite element method has been developed. The original combination of the latter, the Sumudu integral transform and artificial neural networks makes it possible to examine spatially heterogeneous objects of the cryolithozone with a high contrast of geoelectric parameters, significantly reducing computational costs. We have studied the sensitivity of the TEM signals to the boundary between thawed and frozen rocks depending on the inter-borehole distance and spatial orientations of the transmitters and receivers. Numerical simulation results of the TEM monitoring of industrial facilities located on permafrost are considered. The formation of a talik has been shown to significantly manifest itself in the measured electromagnetic responses, which enables timely prevention of industrial disasters and environmental catastrophes.

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