Industry 4.0 brings advanced technologies and automation, enabling greater productivity and efficiency. However, it also increases energy consumption and carbon emissions. Sustainable manufacturing in the era of Industry 4.0 is vital for addressing the environmental impact by minimizing resource consumption, optimizing energy efficiency, and adopting eco-friendly practices. CNC (Computer Numeric Control) 5-axis milling plays a significant role in the machining of precision molds and dies, aerospace, consumer electronics, etc. This research aims to explore the potential of the artificial intelligence (AI) technique in improving energy efficiency on CNC 5-axis milling process for sustainable manufacturing. The experiments were carried out using a CNC 5-axis milling machine with various machining parameters, toolpath planning, and dry cutting condition. Data on energy consumption of CNC 5-axis milling machine were collected, then the relationship between the machining parameters, toolpath planning, and energy consumption was analyzed and built. Furthermore, an AI-based optimization algorithm was developed that uses machine learning and optimization techniques to identify energy-efficient machining strategies. The developed algorithm was implemented within a simulation framework, the algorithm was rigorously assessed using machine learning analysis to effectively reduce energy consumption, all while ensuring the accuracy of machining result. Finally, verification experiments were carried out to evaluate the applicability and performance of the AI-based energy optimization approach on a CNC 5-axis milling machine. By integrating AI into the CNC machining process, manufacturers can achieve significant energy savings, leading to reduced carbon emissions and improved sustainability.

The layered oxyselenide BiCuSeO has garnered significant interest because of its ability to exhibit low thermal conductivity and a high Seebeck coefficient. This study involved the preparation of Bi1−xNaxCuSeO1−xFx(x=0, 0.05, 0.10, 0.15, and 0.20) ceramics through the processes of high energy ball milling and cold isostatic pressing. A systematic study was conducted to investigate the impact of co−doping Na/F on the thermoelectric performance of BiCuSeO ceramics. Replacing Bi3+ with Na+ results in the introduction of a considerable amount of holes, leading to a notable enhancement in the electrical conductivity and power factor. The conductivity was significantly increased from 9.10 S cm−1 in the pure BiCuSeO to 94.5 S cm−1 in Bi0.85Na0.15CuSeO0.85F0.15 at 323 K. At 823 K, the power factor of the Bi0.85Na0.15CuSeO0.85F0.15 sample achieved 44.8×10−5 W/ m K2. Moreover, the Bi1−xNaxCuSeO1−xFx ceramics exhibit a minimum thermal conductivity of 0.43 W m−1 K−1. As a result, the Bi0.85Na0.15CuSeO0.85F0.15 sample achieves a maximum ZT value of 0.78, which is 7.09 times greater than that of the pure BiCuSeO sample (0.11) .

Several studies have shown that combining TiO2 and Cu2O enhances the photocatalytic activity of the material by generating a heterojunction capable of extending the light absorption in the visible and reducing the electron-hole recombination rate. Ball milling has been chosen as an alternative methodology for photocatalyst preparation, among the several techniques documented in the literature review. The results of a previously reported investigation enabled the identification of the most effective photocatalyst that can be prepared for hydrogen generation by combining Cu2O and TiO2 (i.e., 1%wt Cu2O in TiO2 photocatalyst prepared by ball-milling method at 200 rpm and 1 min milling time). To optimize photocatalytic hydrogen generation in the presence of the greatest photocatalyst, the effects of (i) sacrificial species and their concentration, (ii) temperature, and (iii) pH of the system are taken into account, resulting in a light-to-chemical energy efficiency of 8% under the best-tested conditions. Last but not least, the possibility of using the present photocatalytic system under direct solar light irradiation is evaluated: the results indicate that nearly 60% of the hydrogen production recorded under sunlight can be attributed to the visible component of the solar spectrum, while the remaining 40% can be attributed to the UV component.

The extended finite element method (XFEM) is now widely used in crack simulations and the direction of propagation of the XFEM crack in SS-304 has been quite interesting. This paper discusses the direction,θ(+ve in anti-clockwise direction), of the XFEM crack (under mode-I) propagation based on the von Mises stress-field around the crack-tip for static loading under the light water reactors (LWRs) conditions. The Mises core region from the crack-tip is chosen for β=〖90〗^o, the angle between the load and the initial crack. Experimental data for LWRs’environment is obtained and the simulations have been carried out applying three different static loads and results for the direction of the crack are calculated by measuring the stress-field around the crack-tip. It has been found that the direction of the XFEM varies inversely as the stress-field around the crack tip is increased. The computational results have been compared with the numerical results based on the minimum strain energy density criterion[S-criterion].

Conducting computational stress-strain analysis using finite element methods (FEM) is a common approach when dealing with the complex geometries of atherosclerosis, which is a leading cause of global mortality and complex cardiovascular disease. The considerable expense linked to FEM analysis encourages the substitution of FEM with a considerably faster data-driven machine learning (ML) approach. This study investigated the potential of end-to-end deep learning tools as a more effective substitute for FEM in predicting stress-strain fields within 2D cross sections of arterial wall. We first proposed a U-Net based fully convolutional neural network (CNN) to predict the von Mises stress and strain distribution based on the spatial arrangement of calcification within arterial wall cross-sections. Further, we developed a conditional generative adversarial network (cGAN) to enhance, particularly from the perceptual perspective, the prediction accuracy of stress and strain field maps for arterial walls with various calcification quantities and spatial configurations. On top of U-Net and cGAN, we also proposed their ensemble approaches to further improve the prediction accuracy of field maps. Our dataset, consisting of input and output images, was generated by implementing boundary conditions and extracting stress-strain field maps. The trained U-Net models can accurately predict von Mises stress and strain fields, with structural similarity index scores (SSIM) of 0.854 and 0.830 and mean squared errors of 0.017 and 0.018 for stress and strain, respectively, on a reserved test set. Meanwhile, the cGAN models in a combination of ensemble and transfer learning techniques demonstrate high accuracy in predicting von Mises stress and strain fields, as evidenced by SSIM scores of 0.890 for stress and 0.803 for strain. Additionally, mean squared errors of 0.008 for stress and 0.017 for strain further support the model's performance on a designated test set. Overall, this study developed a surrogate model for finite element analysis, which can accurately and efficiently predict stress-strain fields of arterial walls regardless of complex geometries and boundary conditions.

The future of intelligent manufacturing machines involves autonomous selection of process parameters to maximize productivity while maintaining quality within specified constraints. To effectively optimize process parameters, these machines need to adapt to existing uncertainties in the physical system. This paper proposes a novel framework and methodology for feedrate optimization that is based on a physics-informed data-driven digital twin with quantified uncertainty. The servo dynamics are modeled using a digital twin, which incorporates the known uncertainty in the physics-based models and predicts the distribution of contour error using a data-driven model that learns the unknown uncertainty on-the-fly by sensor measurements. Using the quantified uncertainty, the proposed feedrate optimization maximizes productivity while maintaining quality under desired servo error constraints and stringency (i.e., the tolerance for constraint violation under uncertainty) using a model predictive control framework. Experimental results obtained using a 3-axis desktop CNC machine tool and a desktop 3D printer demonstrate significant cycle time reductions of up to 38% and 17% respectively, while staying close to the error tolerances compared to the existing methods.

The mechanical properties, structure, and morphology of milled and 3D printed Ti-6Al-4V alloys produced by selective laser melting were investigated in this study. The mechanical properties were investigated by means of nanoindentation, a tensile test, and a three-point bending test. The Atomic Force Microscope (AFM) was used to study the surface topography and roughness of both titanium alloys. The surface structure and phase analysis were studied by a scanning electron microscope (SEM) and through powder X-ray diffraction (XRD). The results from the nanoindentation experiments showed that the 3D-printed sample has higher indentation hardness and modulus than the milled one. The AFM observation of the surface topography of the samples showed that the milled sample has a higher roughness than the 3D-printed one. The tensile test results showed that the 3D printed sample by means of Selective Laser Melting (SLM) technology has about 26% higher tensile strength and smaller elongation than the milled one. The 3-point bending test revealed that the 3D printed Ti-6Al-4V sample has higher flexural strength than the milled one. It was found that the 3D-printed sample has a smaller crystal size than the milled one, which, according to the Hall-Petch relationship, leads to its higher indentation hardness.

Several studies have highlighted the importance and evolution of process safety leadership and culture in various industries. However, none has focused on the palm oil milling industry yet. This paper critically reviews the latest developments in the palm oil milling process and unit operations leading to process safety concerns. It also discusses the Principle of 3C that is applied to explain repeat accidents and the four-level safety culture in the palm oil milling industry. For this purpose, the author presents case studies of two key palm oil companies in Malaysia. Overall, this paper offers guidelines to leaders in the palm oil milling industry about the required process safety leadership and culture to be understood in order to improve their safety outcomes.

This study was aimed at investigating the effects of roasting and cryogenic milling on particle size, water solubility, and amount of bioactive components in ginseng root. Samples were pulverized by cryogenic milling, and one treatment condition was selected for each size range (10-50 μm, and >50 μm). The selected samples were roasted at different temperatures (160-200 °C) followed by cryogenic milling. Powdered samples were analyzed for their physicochemical characteristics. Results revealed that roasted samples exhibited significantly smaller particle size than controls (not roasted, p<0.05), and the particle size of roasted samples decreased with increase in roasting temperature until flocculation occurred around 180-190 °C. With decrease in particle size of ginseng; water solubility index, antioxidant activity, total polyphenol content, and total polysaccharide content were observed to increase. Ginseng samples showed an increasing trend in antioxidant activity with decreasing particle size (p>0.05) with significantly higher values for all roasted samples compared to the controls (p<0.05). Non-polar ginsenosides such as Rg2(S), Rg2(R), Rg3(S), Rg3(R), and Rh1(S) showed an increase in temperature-dependent manner. This study revealed that roasting at certain temperature range with cryogenic milling has positive effects on ginseng by reducing its particle size, and increasing water solubility and bioactive components.

This study aimed to decrease the particle size of ginseng by roasting and cryogenic milling to increase its water solubility and physiological activity. The samples were roasted for different times (9–21 min) and generated in different sizes (10–50, and > 50 μm). All roasted samples revealed significantly smaller particle sizes than did non-roasted samples based on Sauter mean diameter (D [3,2], p < 0.05). Further, the particle sizes of roasted samples decreased till roasting up to 15 min. In terms of the water solubility index (WSI), antioxidant activity, total polyphenol content (TPC), and total polysaccharides according to particle size, 10-20 μm-sized samples showed the highest values when compared with >50 μm-sized samples. Based on roasting time, WSI values of all samples roasted for up to 15 min were higher than those of the control (not roasted) (p < 0.05). Antioxidant activity and TPC also increased with increasing roasting time. Total polysaccharide content was the highest upon roasting for 15 min except for the 10-20 μm sample. Ginsenoside content of roasted samples >20 μm size was higher than that of the control (not roasted) except after 15 min of roasting. Therefore, roasting and cryogenic milling are effective in producing ginseng root powder.

The system design for Sustainable manufacturing is to consider both ecological and financial constraints. Manufacturing industry demands to advance in sustainable manufacturing by accounting in environment factors in it. All the factors that affect the environment need to be analyzed so that proper amendments or suggestions can be provided. To favour this, Computer aided life cycle inventory model has been presented for CNC milling and turning processes. Based on utilization of resources and stages, whole machining operation time can be divided into three phases named as process (milling or turning), idle and basic times. As parameters are different for evaluating the process times i.e. depth and width of cut in case of milling and initial and final diameters for turning, two different case studies has been presented one for each milling and turning process. Effect of material choice has been studied for different processes. Highly dense and hard materials takes more time in finishing the job due to low cutting speed and feed rates as compared to that of sot materials. In addition, face milling takes more time and consumes more power as compared to peripheral milling due to more retraction time caused by over travel distance and lower vertical transverse speeds than the horizontal transverse speed used in peripheral retraction process.

The microstructural evolution, thermodynamics and kinetics of Mo-21%Tm₂O₃ (mass fraction, %) powder mixtures during ball milling were investigated using X-ray diffraction and transmission electron microscopy. Ball milling induced Tm₂O₃ to be decomposed and then dissolved into Mo crystal. The supersaturated nanocrystalline solid solution of Mo (Tm, O) was obtained after 96 h of ball milling. The elements of Mo, Tm and O were distributed uniformly in the ball-milled particles. Based on the semi-experimental theory of Miedema, a thermodynamic model was developed to calculate the driving force of phase evolution. There was no chemical driving force to form a crystal solid solution of Tm atoms in Mo crystal or an amorphous phase because the Gibbs free energy for both processes was higher than zero. For Mo-21%Tm₂O₃, it was mechanical work, not negative heat of mixing, that provided the driving force to form supersaturated nanocrystalline Mo (Tm, O) solid solution.

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

Present research focus on the effect of different processing routes on physical and mechanical properties of nano Ti(CN) based cermet’s with metallic binders. Tungsten carbide (WC) is added as secondary carbide and Ni-Co are added as metallic binders in nano Ti(CN) based cermet processed via Conventional and Spark Plasma Sintering (SPS). Systematic comparison of composition and sintering conditions for different cermet’s’ systems is done to design the novel composition and sintering conditions. Nano TiCN powder was prepared by 30hr of ball milling. Highest density of >98.5% percent was achieved for SPS processed cermet’s sintered at 1200°C and 1250°C for 3 min. at 60 MPa pressure in comparison to conventionally sintered cermet’s at 1400°C for 1hr with two stage compaction process, uniaxially at 150 MPa and isostatically at 300 MPa pressure respectively. Comparison of XRD analysis of milled powders of different time intervals was done to understand the characteristics of the as received and milled powders. Peak broadening was observed after 5 hr. of ball milling, and it increased till 30hr. Also, peak broadening and refined carbide size was observed in XRD and SEM micrograph of SPS processed cermet. TEM analysis of milled powder showed internal structure having regular periodic arrangement of planes. SEM (BSE) images of all cermets primarily showed three major microstructural phases of core-rim-binder with black, grey and white contrast respectively. In the present sintering conditions high hardness of ~16 GPa and fracture toughness of ~9 MPa m1/2 was obtained for SPS processed cermets sintered at higher temperature.

The present work was aimed at the investigation of effect of High Energy Ball Milling (HEBM) time on the sintering kinetics, structure, and properties of heavy tungsten alloy (HTA) W-7%Ni-3%Fe. The HTA samples were obtained from the nanopowders (20-80 nm) by conventional liquid phase sintering (LPS) in hydrogen and by Spark Plasma Sintering (SPS) in vacuum. The HTA density was shown to depend on the HEBM time non-monotonously that originates from the formation of non-equilibrium solid solutions in the W-Ni-Fe systems during HEBM. The SPS kinetics of the HTA nanopowders was shown to have a two-stage character, the intensity of which depends on the Coble diffusion creep rate and on the intensity of diffusion of the tungsten atoms in the crystal lattice of the -phase. The kinetics of sintering of the initial submicron powders have a single-stage character originating from the intensity of the grain boundary diffusion in the -phase. The dependencies of the hardness and of the yield strength on the grain sizes were found to obey the Hall-Petch relation. The hardness, strength, and dynamic strength in the compression tests of the fine-grained tungsten alloys obtained by SPS and LPS were studied.

Heterojunctions based on metal halide perovskites (MHPs) are promising systems for photocatalytic hydrogen evolution reaction (HER). In this work, we coupled Cs3Bi2Br9 nanocrystals (NCs), obtained by wet ball milling synthesis, with g-C3N4 nanosheets (NSs) produced by thermal oxidation of bulk g-C3N4 in air. These methods are reproducible, inexpensive, and easy to scale up. Heterojunctions with different loadings of Cs3Bi2Br9 NCs were fully characterized and tested for HER. A relevant improvement of H2 production with respect to pristine carbon nitride was achieved at low NCs levels reaching values up to about 4600 µmol g-1 h-1. This work aims to provide insights on the synthesis of inexpensive and high performing heterojunctions using MHP for photocatalytic applications.

Palm oil production involves unit operations that lead to greenhouse gas emissions. A typical mill is estimated to produce greenhouse gas emissions of 637–1,131 kg CO2 eq/t crude palm oil. There has been a huge effort to reduce the carbon footprint of palm oil mills. However, the data from such research have not been consolidated. This paper reviews significant information, results, and conclusions derived from studies in the literature. The latest developments in palm oil milling operations and information on greenhouse gas emissions are presented. Current initiatives in reducing carbon footprint of palm oil mills are discussed along with an assessment of the technologies employed. These include the conventional method of capturing biogas from palm oil mill effluent and emerging approaches such as converting palm sludge oil into biodiesel, deploying enzyme-assisted oil extraction, and converting biomass into fuel for energy generation as an alternative to coal and other fossil fuels. The importance of self-sufficiency is deliberated because a self-sufficient palm oil mill is estimated to reduce emissions by 457 kg CO2 eq/t crude palm oil compared to a mill that requires an external power supply. Methods with the greatest positive effect on the carbon footprint are identified for further investigation.

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

The main purpose of the work was to determine the powder composition and properties of iron-based sinters used as a matrix in metallic diamond tools. The FeCuSnNi sinters obtained from a mixture of commercial powders were used for research. The material was subjected to milling in a ball mill and sinter have been performed with the hot-pressing technique in a graphite mold. The physical and mechanical properties of the sinters were tested. The investigations of obtained sinters concerned density, porosity, hardness, oxygen content, static tensile test and analysis of microstructure. The static tensile test was carried out with the use of a universal testing machine UTS-100. The study of the microstructure of sinters was performed using the JSM-7100F scanning electron microscope. The sinters are characterized by a high coefficient of thermal expansion and good mechanical properties. The FeCuSnNi sinter revealed a complex phase structure. Spot chemical analysis revealed the presence of Fe solution, Cu solution and the presence of iron oxides oxides. Nickel atoms were present throughout the sinter volume. The obtained results indicate that the presented material can be used as a matrix of a diamond-metal composite in diamond tools for the processing of natural stones and building materials.

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

In laser-assisted milling, higher temperature in shear zone softens the material potentially resulting in a shift of mean residual stress, which significantly affects the damage tolerance and fatigue performance of product. In order to guide the selection of laser and cutting parameters based on the preferred mean residual stress, inverse analysis is conducted by predicting residual stress based on guessed process parameters, which is defined as the forward problem, and applying iterative gradient search to find process parameters for next iteration, which is defined as the inverse problem. An analytical inverse analysis is therefore proposed for the mean residual stress in laser-assisted milling. The forward problem is solved by analytical prediction of mean residual stress after laser-assisted milling. The residual stress profile is predicted through the calculation of thermal stress, by treating laser beam as heat source, and plastic stress by first assuming pure elastic stress in loading process, then obtaining true stress with kinematic hardening followed by the stress relaxation. The variance-based recursive method is applied to solve inverse problem by updating process parameters to match the measured mean residual stress. Three cutting parameters including depth of cut, feed per tooth, and cutting speed, and two laser parameters including laser-tool distance and laser power, are updated with respected to the minimization of resulting residual stress and measurement in each iteration. Experimental measurements are referred on the laser-assisted milling of Ti-6Al-4V grade 5 and ELI. The percentage difference between experiments and predictions is less than 5% for both materials, and the selection is completed within 50 loops.

This paper presents a study of the 2D roughness profiles on a flat surface generated on a steel work piece, by ball nose end milling, with linear equidistant tool paths (pick-intervals). The milled surface exploration with a surface roughness tester (on pick and feed directions) produces usually 2D roughness profiles having periodical evolutions. These evolutions can be seen as time-dependent signals, describable as a sum of sinusoidal components (the wavelength of each component being regarded as a period). In order to detect a well approximated description of these components, two appropriate signal processing techniques are used in this work: first technique provides a direct mathematical (analytical) description and is based on computer aided curve (signal) fitting (more accurate), the second technique (providing an indirect and incomplete description) is based on the spectrum generated by fast Fourier transform (less accurate). This study can be seen as a way for a better understanding the interaction between the tool and the work piece or to achieve a mathematical description of the machined surface microgeometry in terms of roughness (e. g. its description being a collection of closely spaced 2D roughness profiles).

Printed circuit boards (PCBs) make up 3 to 5% of all electronic waste. The metal content of spent PCBs can reach 40%. They usually contain valuable metals - Ag, Au, Pd, as well as other metals such as Cu, Sn, Pb, Cd, Cr, Zn, Ni and Mn. However metallic part of whole PCB is 40-60% including Cu layers between the fiberglass-polymer layers. The paper describes economics of the the valuable metals (Ag, Au, Pd) containing concentrate preparation from raw PCB. Considered the influence of pre-treatment method of PCB before the extraction of valuable metals on the extraction self-cost change. The disintegration method is based on the high-energy impact of the particles of the material to be ground, thus causing the separation of the metallic components of the PCB. In the course of the work, single and double direct grinding by the method of disintegration was studied. For the calculation the test batch 10’000 kg of two types of PCB were taken for estimation of the self costs and potential profit in case of complete valuable metals (Ag, Au, Pd) plus Cu extraction. Is shown that from 10’000 kg of studied PCB is possible to obtain 1144 and 1644 kg of metal-rich concentrate, which further should be subjected to electro-hydrochlorination for metals leaching. The main findings shows, that despite of high content of Cu in studied PCB, the commercial value is insignificant in total income from the Ag, Au and Pd sale. Was established correlation of the self-cost decrease after separative disintegration of PCB by metal content increase (by specific metals such as Au, Ag, Pd and also Cu) with the metal potential yield after extraction.

There exists a great potential for application in soil mechanics with the use of mechanistic models such as the finite element model, and other non-local models, such as the peridynamic model, since the soil mass can be modelled as a group of particles that interact with each other. In this work, we determine the bearing capacity of deep foundations where stresses are transmitted to deep deposits, which generally present better characteristics in terms of compressibility and shear strength. One of the main elements are driven piles, which must be previously designed, made and tested before their final use. In this work, a comparison with different methods is presented. The finite element method and the peridynamic model are used. Because the effect of consolidation was not taken into account, in these examples, the assumption is made that the structures are cemented in inert soil. Likewise, it is assumed that the water table has no influence on the simulations. Failure envelopes were observed where bonds between particles present high gradients of deformation and fracture. Bearing capacities were estimated and compared with those obtained from the Terzaghi and Meyerhof methods.

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

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

Swarf of aluminum alloy with high corrosion resistance and ductility was successfully converted into fine hydroreactive powders via its ball milling with silver powder and either lithium chloride or gallium. The later substances significantly intensified particle size reduction, while silver formed ‘cathodic’ sites (Ag, Ag2Al) promoting Al corrosion in aqueous saline solutions with hydrogen generation. The diffraction patterns, microphotographs, and elemental analysis results demonstrated partial aluminum oxidation in the samples and their contamination with tungsten carbide from milling balls. Those effects were responsible for obtaining hydrogen yields lower than expected. For AlCl3 solution at 60 °C, Al–LiCl–Ag, Al–LiCl, Al–Ga–Ag, and Al–Ga composites delivered (84.6±0.2), (86.8±1.4), (80.2±0.5), and (76.7±0.7)% of the expected hydrogen respectively. Modification with Ag promoted Al oxidation thus providing higher hydrogen evolution rates. The samples with Ag were tested in CaCl2 solution as well, for which the reaction proceeded much slower. At a higher temperature (80 °C) after 3 h of experiment, the hydrogen yields for Al–LiCl–Ag and Al–Ga–Ag powders were (46.7±2.1) and (31.8±1.9)%. The tested Ag-modified composite powders were considered promising for hydrogen generation and had a potential for further improvement to provide higher performance.

Modeling and assessing the sustainability of machining systems has been considered to be a crucial approach to improving the environmental performance of the machining processes. As the most common machining system, the computer numerical control (CNC) milling system is a typical man-machine cooperative system where the activities of the machine tool and operator generate material and energy consumption. However, the energy consumption of the operator in the CNC milling system has often been ignored in most existing research. Therefore, existing methods fail to provide a comprehensive understanding of the sustainability of the CNC milling system. To fill this gap, an exergy loss assessment method is proposed to investigate the sustainability of the CNC milling system, where the energy consumption of the operator, the energy consumption of the machine tool, and material consumption are taken into consideration. The key performance indexes of the energy consumption of the operator, the energy consumption of the machine tool, the exergy loss, and the specific exergy loss (SEL) are analyzed and modeled. To demonstrate the feasibility of the proposed method, a case study was carried out on a three-axis machining center (XH714D), in which SEL was found to be 88.04J/mm3. The proposed method is effective to assess the sustainability of the CNC milling system, and the established exergy loss models build a good basis for exergy efficiency optimization.

Various metals and semiconductors containing Printed Circuit Boards (PCBs) are abundant in any electronic device equipped with controlling and computing features. These devices inevitably constitute E-waste after the end of service life. The typical construction of PCBs includes mechanically and chemically resistive materials, which significantly reduce the reaction rate or even avoid accessing chemical reagents (dissolvents) to target metals. Additionally, the presence of relatively reactive polymers and compounds from PCBs requires high energy consumption and reactive supply due to the formation of undesirable and sometimes environmentally hazardous reaction products. Preliminarily milling PCBs into powder is a promising method for increasing the reaction rate and avoiding liquid and gaseous emissions. Unfortunately, current state-of-the-art milling methods also lead to the presence of significantly more reactive polymers still adhered to milled target metal particles. This paper aims to find a novel single and two-stage disintegration-milling approach that can provide the formation of metal-rich particle size fractions. The morphology, particle fraction sizes, bulk density, and metal content in produced particles were measured and compared. Research results show the highest bulk density (up to 6.8 g·cm-3) and total metal content (up to 95.2 wt. %) in finest sieved fractions after the single-step milling of PCBs. Therefore, the concentrations of about half tested metallic elements are higher in the single milled specimen and with lower adhered plastics concentrations, as compared to double milled specimens.

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

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

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

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

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

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

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

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

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

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

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 paper, we aim to address the challenge of airflow interference during fault detection in high-speed train bogies by introducing a flow field and investigating the characteristics of the sound field distribution of critical components under its influence. This approach overcomes the limitation observed in previous studies that ignored the impact of the flow field. Furthermore, we evaluate and compare various layouts for inter-track acoustic sensor arrays. The proposed sensing methods hold significant practical engineering value, and the collected data serves as a valuable reference for acquiring acoustic signals from key bogie components in future research.

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

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

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.

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

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

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

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

Lithium titanates are used in various applications, such as anode materials for lithium intercalation (Li4Ti5O12) or breeding materials in fusion reactors (Li2TiO3). Here, we report the formation of nano-crystalline lithium titanates by a mechanochemical approach and present a deeper insight into their structural characteristics by PXRD and solid-state NMR spectroscopy. The compounds were synthesized in a high energy planetary ball mill with varying milling parameters and different grinding tools. NaCl type Li2TiO3 (α-Li2TiO3) was formed by dry milling of lithium hydroxide with titania (rutile or anatase) and by a milling induced structure transformation of monoclinic β Li2TiO3 or spinel type Li4Ti5O12. Heating of mechanochemical prepared α Li2TiO3 induces a phase transformation to the monoclinic phase similar to hydrothermal reaction products, but a higher thermal stability was observed for the mechanochemical formed product. Microstructure and crystallographic structure were characterized by PXRD via Rietveld analysis. Detailed phase analysis shows the formation of the cubic phase from the various educts. A set of two lattice parameters for α Li2TiO3 was refined, depending on the presence of OH- during the milling process. An average crystallite size of less than 15 nm was observed for the mechanochemical generated products. The local Li environment detected by 6Li SPE MAS NMR revealed Li defects in the form of tetrahedral instead of octahedral site occupation. Subsequent adjustment of the structural model for Rietveld refinement leads to better fits, supporting this interpretation.

Barriers to inventing electronic devices involve challenges of iterating electronic designs due to long lead times for professional circuit board milling or high-costs of commercial milling machines. To overcome these barriers this study provides open source (OS) designs for a low-cost circuit milling machine. First, design modifications for mechanical and electrical sub-systems of the OS D3D Robotics prototyping system are provided. Next, Copper Carve, an OS custom graphical user interface, is developed to enable circuit board milling by implementing backlash and substrate distortion compensation. The performance of the OS D3D circuit mill is then quantified and validated for: positional accuracy, cut quality, feature accuracy and distortion compensation. Finally, the return on investment is calculated for inventors using it. The results show by properly compensating for motion inaccuracies with Copper Carve, the machine achieves a motion resolution of 10 microns, which is more than adequate for most circuit designs. The mill is at least five times less expensive than all commercial alternatives and the material costs of the D3D mill are repaid from fabricating 20-43 boards. The results show that the OS circuit mill is of high-enough quality to enable rapid invention and distributed manufacturing of complex products containing custom electronics.

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

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

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

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

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

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

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.

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.

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

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

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

This article presents a case study on the structural assessment of a reinforced concrete (RC) foundation exposed to low temperatures. The foundation supports a 19,500 m³-capacity tank with low-temperature (-89°C) ethane. Icing and bubbling were observed on the tank’s surface soon after it started operations. Condensation was also observed at the bottom of the 0.8 m-depth RC slab, which raised concerns about the structural condition of the concrete. This study provides details of the field and analytical investigations conducted to assess the structural condition of the foundation. Heat transfer finite element (FE) analyses were performed to examine the concrete sections subjected to low temperatures. It was found that the ethane leakage produced a low temperature on the top side of the concrete foundation of +9.7°C. Overall, the temperatures calculated by the FE analyses were in good agreement with actual field measurements, within a ±5% accuracy. The simplified heat transfer equation for porous media used in this study was sufficiently accurate to model the effects of the ethane leakage in the concrete foundation, provided the ambient temperature at the site is taken into account in the analysis. The results also confirm that reinforcing bars can be neglected in the thermal analysis of massive concrete slabs. The results from the field measurements and FE analyses confirmed that the structural integrity of the RC foundation was never compromised. The approaches, methods and techniques discussed in this article are deemed suitable to solve the practical and scientific challenges involved in the structural assessment and repairs of large special structures. Accordingly, they can serve as useful reference and guidance for engineers and practitioners working in the field of forensic engineering.

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

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

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

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

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

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

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

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

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

Comparison between hydraulic and hydrologic computational methods is conducted in this study, regarding prismatic open channels under unsteady subcritical flow conditions. One-dimensional unsteady flow continuity and momentum equations are solved using explicit and implicit finite difference schemes for a symmetrical trapezoidal cross section, where the flow discharge and depth are the dependent variables. The results have been compared to those derived from Muskingum-Cunge hydraulic/hydrologic method as well as the commercial software HEC-RAS. The results from explicit and implicit code compare well to those from commercial software and hydraulic/hydrologic methods for long prismatic channels, thus directing the hydraulic engineer to quick preliminary design of prismatic open channels for unsteady flow with satisfactory accuracy.

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

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

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

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

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

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

In this paper, we analyze the convergence of a finite difference method with the implicit forward time central space (FTCS) scheme for solving the three-dimensional advection-diffusion-reaction equation (ADRE). It is proved that the method is unconditionally convergent. We apply the scheme to obtain numerical solutions for the transport of pollutants in street tunnel problems with various reaction coefficients and various rates of change of concentrations of sources or sinks of pollution.

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.

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

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

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

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

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

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

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

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

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

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

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

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

Leaf springs are suspension components constituted of several structural and assembly components. Once the leaves are assembled, a stress state due to the assembly process is developed. Thus, the quantification of the assembly stress level should be known in order to be accounted for in the structural design. The aim of this article is precisely to quantify the assembly stress level in parabolic leaf springs for railway freight wagons. An experimental method with the aid of extensometry is considered to evaluate elastic assembly stresses. Next, a numerical model of the leaf spring assembly process with the required boundary conditions is developed to support the experimental results and allow obtaining a distribution of surface stresses along the master leaf length. The development of this work resulted in a surface stress profile along the length of the master leaf generated by the assembly process, where the deviations of the numerical outcomes regarding the analysed spots with strain gauges are mostly of little relevance.

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

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

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

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

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