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.

We present the results on photothermal (PT) and heat conductive properties of nanogranular silicon (Si) films synthesized by evaporation of colloidal droplets (drop-casting) of 100 ± 50 nm sized crystalline Si nanoparticles (NP) deposited on glass substrates. Finite difference time domain (FDTD) and finite element mesh (FEM) modeling of absorbed light intensity and photo-induced spatial temperature distribution across the Si NP films were well correlated with the local temperatures measured by micro-Raman spectroscopy and used for determination of heat conductivities in the films of various thicknesses. Cubic-to-hexagonal phase transition in these films caused by laser heating was found to be heavily influenced by the film thickness and heat conductive properties of glass substrate, on which the films were deposited. Heat conductivities across the drop-casted Si nanogranular films were found to be in the range of lowest heat conductivities of other types of nanostructurely voided Si films due to enhanced phonon scattering across inherently voided topology, weak NP-NP and NP-substrate interface bonding within nanogranular Si films.

This paper presents the current lexical, morphological, syntactic and rule-based machine translation work for Erzya and Moksha that can and should be used in the development of a roadmap for Mordvin linguistic research. We seek to illustrate and outline initial problem types to be encountered in the construction of an Apertium-based shallow-transfer machine translation system for the Mordvin language forms. We indicate reference points within Mordvin Studies and other parts of Uralic studies, as a point of departure for outlining a linguistic studies with a means for measuring its own progress and developing a roadmap for further studies.

Cogging torque is a pulsating, parasitic and undesired torque ripple intrinsic of the design of a permanent magnet synchronous generator (PMSG), which should be minimized due to its adverse effects: vibration and noise. In addition, as aerodynamic power is low during start-up at low wind speeds in small wind energy systems, the cogging torque must be as low as possible to achieve a low cut-in speed. A novel mitigation technique using compound pre-slotting, based on a combination of magnetic and non-magnetic materials, is investigated. The finite element technique is used to calculate the cogging torque of a real PMSG design for a small wind turbine, with and without using compound pre-slotting. The results show that cogging torque can be reduced by a factor of 48% with this technique, while avoiding the main drawback of the conventional pre-slotting technique: the reduction of induced voltage due to leakage flux between stator teeth. Furthermore, through a combination of pre-slotting and other cogging torque optimization techniques, 84%, cogging torque can be eliminated for a given design.

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

The effect of thermal radiation on the two – dimensional, steady-state, conjugate heat transfer from a circular cylinder with an internal heat source in steady laminar crossflow is investigated in this work. P₀ (Rosseland) and P₁ approximations were used to model the radiative transfer. The mathematical model equations were solved numerically. Qualitatively, P₀ and P₁ approximations show the same effect of thermal radiation on conjugate heat transfer; the increase in the radiation – conduction parameter decreases the cylinder surface temperature and increases the heat transfer rate. Quantitatively, there are significant differences between the results provided by the two approximations.

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.

Here we investigate how the local properties of particles in a thermal bath may influence the thermodynamics of the bath, and consequently alter the statistical mechanics of subsystems that comprise the bath. We are guided by the theory of small-system thermodynamics, which is based on two primary postulates: that small systems can be treated self-consistently by coupling them to an ensemble of similarly small systems, and that a large ensemble of small systems forms its own thermodynamic bath. We adapt this “nanothermodynamics” to investigate how a large system may subdivide into an ensemble of smaller subsystems, causing internal heterogeneity across multiple size scales. For the semi-classical ideal gas, maximum entropy favors subdividing a large system of “atoms” into an ensemble of “regions” of variable size. The mechanism of region formation could come from quantum exchange symmetry that makes atoms in each region indistinguishable, while decoherence between regions allows atoms in separate regions to be distinguishable by their distinct locations. Combining regions reduces the total entropy, as expected when distinguishable particles become indistinguishable, and as required by a theorem in quantum mechanics for sub-additive entropy. Combining large volumes of small regions gives the usual entropy of mixing for a semi-classical ideal gas, resolving Gibbs paradox without invoking quantum symmetry for atoms that may be meters apart. Other models presented here are based on Ising-like spins, which are solved analytically in one dimension. Focusing on the bonds between the Ising-like spins we find similarity in the equilibrium properties of a two-state model in the nanocanonical ensemble and a three-state model in the canonical ensemble. Thus, emergent phenomena may alter the thermal behavior of microscopic models, and the correct ensemble is necessary for fully-accurate predictions. Another result using Ising-like spins involves simulations that include a nonlinear correction to Boltzmann’s factor, which mimics the statistics of indistinguishable states by imitating the dynamics of spin exchange on intermediate lengths. These simulations exhibit 1/f-like noise at low frequencies (f), and white noise at higher f, similar to equilibrium thermal fluctuations found in many materials.

In this paper, we present the extension of the pressure-based solver designed for the simulation of compressible and/or incompressible two-phase flows of viscous fluids. The core of the numerical scheme is based on the hybrid Kurganov— Noele — Petrova/PIMPLE algorithm. The governing equations are discretized in the conservative form and solved for velocity and pressure, with the density evaluated by an equation of state. The acoustic-conservative interface discretization technique helps to prevent the unphysical instabilities on the interface. The solver was validated on various cases in wide range of Mach number, both for single-phase and two-phase flows. The numerical algorithm was implemented on the basis of the well-known open-source Computational Fluid Dynamics library OpenFOAM in the solver called interTwoPhaseCentralFoam. The source code and the pack of test cases are available on GitHub: https://github.com/unicfdlab/hybridCentralSolvers

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

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

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

The variable-length arbitrary Lagrange-Euler (ALE)-ANCF finite element, which employ nonrational interpolating polynomials, cannot exactly describe the rational cubic Bezier curves such as conic and circular curves. The rational absolute nodal coordinate formulation (RANCF) finite element, whose reference length (undeformed length) is constant, can exactly represent the rational cubic Bezier curves. A new variable-length finite element called the ALE-RANCF finite element, which is capable of accurately describe the rational cubic Bezier curves, is proposed by combining the desirable features of the ALE-ANCF and RANCF finite element. In order to control the reference length of ALE-RANCF element within a suitable range, element segmentation and merging schemes are proposed. It is demonstrated that exact geometry and mechanic is maintained after the ALE-RANCF element is divided into two shorter ones, and compared with the ALE-ANCF elements, there are smaller deviations and oscillations after two ALE-RANCF elements are merged into a longer one. Numerical examples are presented and the feasibility and advantages of the ALE-RANCF finite element are demonstrated.

The fatigue life of welded joints under cyclic loading is a complex process that can be roughly divided into fatigue crack initiation, crack propagation and final fracture. Usually, these phases are evaluated separately. Several methods are available to estimate the extent of the different phases in welded components, such as stress-based fatigue concepts for crack initiation or linear elastic fracture mechanics for crack propagation; however, distinguishing these phases in welded components is a complex procedure and often relies on assumptions such as transition crack sizes. The objective of this study is to investigate the relationship between fatigue crack initiation and propagation in welded joints using artificially notched specimens with welded joints characteristics of different notch acuity (different radii and opening angle). The experiments show that the investigated relationship basically depends on the notch acuity, the load level and the stress ratio.

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.

Corrugated cardboard is an ecological material, mainly because, in addition to virgin cellulose fibers also the fibers recovered during recycling process are used in its production. However, the use of recycled fibers causes slight deterioration of the mechanical properties of the corrugated board. In addition, converting processes such as printing, die-cutting, lamination, etc. cause micro-damage in the corrugated cardboard layers. In this work, the focus is precisely on the crushing of corrugated cardboard. A series of laboratory experiments were conducted, in which the different types of single-walled corrugated cardboards were pressed in a fully controlled manner to check the impact of the crush on the basic material parameters. The amount of crushing (with a precision of 10 micrometers) was controlled by a precise FEMat device, for crushing the corrugated board in the range from 10 to 70 % of its original thickness. In this study, the influence of crushing on bending, twisting and shear stiffness as well as a residual thickness and edge crush resistance of corrugated board was investigated. Then, a procedure based on a numerical homogenization, taking into account a partial delamination in the corrugated layers to determine the degraded material stiffness was proposed. Finally, using the empirical-numerical method, a simplified calculation model of corrugated cardboard was derived, which satisfactorily reflects the experimental results.

Steels studs are an inevitable part of drywall construction as they are lightweight and offer the required structural stability. However, the studs act as sound bridges between the plasterboards reducing the overall sound insulation of the wall. Overcoming this often calls for wider cavity walls and complex stud decoupling fixtures that increase the installation cost while reducing the floor area. As an alternative approach, this research reveals the potential of perforated studs to improve the acoustic insulation of drywall partitions. The acoustic and structural performance is characterized using a validated finite element model that acted as a prediction tool in reducing the number of physical tests required. The results established that an acoustic numerical model featuring fluid-structure-interaction can predict the weighted sound reduction index of a stud wall assembly at an accuracy of ±1 dB. The model was used to analyze six perforated stud designs and found to outperform the sound insulation of non-perforated drywall partitions by reducing the sound bridging. Overall, the best performing perforated stud design was found to offer improvements in acoustic insulation of up to 4 dB, while being structurally compliant.

The steady-state, coupled heat and mass transfer from a fluid flow to a sphere accompanied by an exothermal catalytic chemical reaction on the surface of the sphere is analysed taking into consideration the effect of thermal radiation. The flow past the sphere is considered steady, laminar and incompressible. The radiative transfer is modeled by P₀ and P₁ approximations. The mathematical model equations were discretized by the finite difference method. The discrete equations were solved by the defect correction – multigrid method. The influence of thermal radiation on the sphere surface temperature, concentration and reaction rate was analysed for three parameter sets of the dimensionless reaction parameters. The numerical results show that only for very small values of the Prater number the effect of thermal radiation on the surface reaction is not significant.

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

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

Background: Cervical laminoplasty is a useful for treatment of cervical myelopathy. However, this procedure has limitations for kyphotic cervical alignments. We used the finite element (FE) analysis and investigated the biomechanical changes in an intact and laminoplasty models with lordosis, straight, and kyphosis cervical alignments. Methods: A three-dimensional FE model of the cervical spine (C2-C7) with ligaments was created from computed tomography. The model was modified with the following cobb angles and the C3-C6 laminoplasty was conducted; a) laminoplasty-lordotic model (LM-L; C2-C7 angle: -10°), b) laminoplasty-straight model (LM-S; C2-C7 angle: 0°), and c) laminoplasty-kyphotic model (LM-K; C2-C7 angle: 10°). A pure moment with a compressive follower load of 100N to represent the weight of the head/cranium and cervical muscle stabilization was applied to these models. The range of motion (ROM), annular stress, nucleus stress and facet forces were analyzed. Results: ROM of LM-K increased when compared to the other models except for flexion. The LM-K had the highest mobility with 49% increase in ROM observed under extension, compared to the intact model. In all motion except for flexion, LM-L models’ ROM decreased by more than 10%, and LM-S and LM-K models’ ROM increased by more than 10% at C2-C7 compared to the intact model. The annular stresses and nucleus stresses in LM-K were higher compared to the other models. The maximum increase in annular stresses was about 194% in LM-K compared to the intact model, observed at the C3-C4 segment. The facet contact forces were lowest in the LM-K, compared to the other models.Conclusions: Patients with a cervical kyphosis alignment are at a disadvantage of increased kyphosis compared to cases with lordosis or straight alignment and should be treated with caution.

A new two-dimensional nanomaterial – Titanium Carbide MXene (Ti₃C₂-MXene) – was reported in 2011. In this work, the microscale models of Ti₃C₂-MXene nanomaterial are considered with polymer matrix. The nanocomposites are modeled using nacre-mimetic brick-and-mortar assembly configurations due to enhanced mechanical properties and interlocking mechanism between the Ti₃C₂-MXene (brick) and polymer matrices (mortar). The polymer matrix material (Epoxy-resin) is modeled with elastic and viscoelastic behavior (Kelvin-Voigt Model). The Finite Element Method is used for numerical analysis of the microscale models with the multi-point constraint method to include Ti₃C₂-MXene fillers in the polymer matrix. Ti₃C₂-MXenes are considered as thick plate elements with transverse shear effects. The response of elastic and viscoelastic models of polymer matrix are studied. Finally, a tensile and compressive load is applied at the microscale and the effective load transfer due to nacre-mimetic configuration is discussed. This paper provides nacre-mimetic models to pre-design the nanocomposite for optimal performance with damage resistance and enhanced strength.

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.

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

The formability and failure behavior of TRIP steel blanks were investigated through various stress states. The forming limit diagram (FLD) at fracture is constructed both experimentally and numerically. Numerical studies are performed to evaluate the applicability of different damage criteria in predicting the FLD as well as complex cross-die deep drawing process. The fracture surface and numerical results revealed that the material failed in a different mode for different strain path. Therefore, Tresca model which is based on shear stress accurately predicted the conditions where shear had the profound effect on the damage initiation, whereas Situ localized necking criterion was able to calculate the conditions which localization was dominant.

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.

Gol'denveizer's problem of a torus has been analyzed by Audoly and Pomeau (2002) and Sun (2021). However, all of the investigations of Gol'denveizer's problem of an elastic torus have been linear. In this paper, the finite element method is used to more accurately address this problem. Furthermore, Sun (2021) cannot be solved by nonlinear analysis. We research the nonlinear mechanical properties of Gol'denveizer's problem of circular and elliptic tori, and relevant nephograms are given. We study the buckling of Gol'denveizer's problem of an elastic torus, and propose failure patterns and force-displacement curves of tori in the nonlinear range. Investigations reveal that circular tori have more rich buckling phenomena as the parameter a increases. Gol'denveizer's problem of the buckling of an elliptic torus is analyzed, and we find a new buckling phenomenon called a "skirt." As a/b increases, the collapse load of an elliptic torus of the Gol'denveizer problem is enhanced gradually.

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.

Vibration suppression in a field of frequencies and the creation of band gaps which do not allow the propagation of the waves is studied in this paper by means of microstructures designed through topology optimization. Topology optimization is formulated in frequency domain. Band gap design is based on Floquet-Bloch theory and genetic optimization. The result is appealing in view of modern 3D printing techniques.

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.

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

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.

Magnetometric Resistivity(MMR) Method, which can use the power supply method of the traditional apparent resistivity method to measure the magnetic field. At present, the application research abroad is relatively extensive, but the domestic(China) research on the application of the MMR method is very few, and it is not even well-known. This paper is based on the MMR theoretical method under the point source DC condition, combined with the abnormal potential method and the modified Biot-Savart law, and using the three-dimensional numerical calculation method of the finite difference method to calculate the abnormal potential field, electric field, and magnetic field on the matlab2018a platform. Calculate, realize the multi-physics simulation in the electromagnetic field through Matlab platform programming, and verify the correctness of the algorithm by the spherical anomaly model with an analytical solution. Through the visual simulation of the three-dimensional data volume of the multi-physics field in the electromagnetic field, we can better understand the response mechanism of the electromagnetic field under DC conditions and grasp their three-dimensional spatial distribution rules. It is hoped that the research in this article can help the research of MMR Personnel better use this method for exploration.

Zero-rate output (ZRO) drift induces deteriorated micro-electromechanical system (MEMS) gy-roscope performances, severely limiting its practical applications. Hence, it is vital to explore an effective method toward ZRO drift reduction. In this work, we conduct an elaborate investigation on the impacts of the internal and packaging stresses on the ZRO drift at the thermal start-up stage, and propose a temperature-induced stress release method to reduce the duration and magnitude of ZRO drift. Self-developed high-Q dual mass tuning fork gyroscopes (TFGs) are adopted to study the correlations between temperature, frequency and ZRO drift. Furthermore, a rigorous finite element simulation model is built based on the actual device and packaging structure, revealing the temperature and stresses distribution inside TFGs. Meanwhile, the relationship between temperature and stresses are deeply explored. Moreover, we introduce a temperature-induced stress release process to generate thermal stresses and reduce the temperature-related device sensitivity. By this way, the ZRO drift duration is drastically reduced from ~2000 s to ~890 s, and the drift magnitude decreases from ~0.4 °/s to ~0.23 °/s. This stress release method achieves a small bias instability (BI) of 7.903 °/h and a low angle random walk (ARW) of 0.792 °/√h, and the long-term bias performance is significantly improved.

It is well known that the Continuous Galerkin Finite Element (CGFE) method is globally consistent with respect to the first law of thermodynamics. This means that, for any mesh, all obtained discrete solutions will conserve total energy. One might expect, that the method is, also, globally consistent with respect to the second law of thermodynamics. In this paper, we formally study if such conjecture is true. The heat conduction equation is used as the physical model for this analysis. In the present study it is proved that the conjecture is false: at least, for standard piecewise linear (1D and 2D) elements, the CGFE method is not always globally consistent with respect to the second law of thermodynamics. In other words, some obtained discrete solutions can violate the global postulate of the second law, which asserts that total entropy can never decrease.

Despite significant advances in metallurgy and the potential to create high and very high-strength steel, all international specifications for steel design provided little information about the limits of slenderness for high-strength steel sections (HSSs) and didn’t give anything about the design of very high-strength steel sections (VHSSs). The American structural steel design specification covers circular hollow sections (CHS) of steel grades up to S690, while it isn't in ANSI/AISC 360-16 and ASTM A514/A514M. Euro code EN 1993-1-1 can be applicable only to steel grades up to S460. EN 1993-1-3 and EN 1993-1-12 contain supplementary rules for cold-formed members and HSS (S460- S700), respectively. Though ductility requirements for (HSS) are specifically established in EN 1993-1-12, it is largely a continuation of the normal-strength-steel design procedures standardized in EN 1993-1-1, where (HSSs) are largely treated as normal-strength steels in these rules, and the same design approaches are used. The behavior of high and very high-strength steel circular hollow sections under axial compression load is studied in this research. A total of 16 nonlinear finite-element (FE) models were generated to replicate stub column tests that were experimentally tested by others in previous research. Hence, a parametric study is conducted using forty FE models developed to investigate the local buckling behavior under various slenderness ratios comprehensively. The developed models cover slenderness ratios ranging from 20 to 1226 and steel grades S460 and S1100 with yield stress equal 460 MPa and 1152 MPa, respectively. The FE results were combined with 105 previously collected experimental results to assess the applicability of existing codified design methodologies in the Euro code and the North American codes of cold-formed CHS. Based on the results of this study, new cross-sectional slenderness limits and new design equations for more efficient simple designs were presented for high and very-high strength steel circular hollow sections and Compared with the results of experimental tests and FE models.

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.

To investigate the effect of Polyvinylpyrrolidone (PVP) addition and consequently porosity, two different sets of membranes are manufactured, since PVP is a widely used poring agent which has an impact on the mechanical properties of the membrane material. The first set (PAN 1) includes PAN and the necessary solvent while the second set (PAN 2) is made of PAN and PVP. These membranes are put through several characterisation processes including tensile testing. The obtained data are used to model the static behaviour of the membranes with different geometries, but similar loading and boundary conditions that represent their operating conditions. This modelling process is undertaken by using finite element method. The main idea is to investigate how geometry affects the load-carrying capacity of the membranes. Alongside membrane modelling, their materials are modelled with representative elements with hexagonal and rectangular pore arrays (RE) to understand the impact of porosity on the mechanical properties. Exploring the results, the best geometry is found as the elliptic membrane with the aspect ratio 4 and the better RE as the hexagonal array which can predict the elastic properties with an approximate error of 12%.

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

The cuttlebone-like structure is a complex porous bionic structure with asymmetric sinusoidal S-shaped wall structure connecting laminar septa, and studies have shown that the cuttlebone-like structure has light weight, high strength and excellent energy absorption capability. It has become an important research to investigate the mechanical mechanism of cuttlebone-like structure and to design a bionic structure beyond nature structure based on it. In order to investigate the mechanical properties of the static compression and shearing processes of the cuttlebone-like structure, this paper establishes theoretical formulas and analyzes the influence of dimensionless parameters such as height-to-thickness ratio $\lambda$, period-to-thickness ratio $\xi$ and amplitude to period ratio $\eta$ on the compression and shearing processes with the help of ABAQUS finite element software.The parameter sensitivity analysis method was used to compare the important influence degree of each dimensionless parameter on the mechanical properties of the structure and to determine the relative optimal parameters of the structural mechanical properties under compression and shearing. Based on the relatively optimal parameter of the cuttlebone-like structure, the cross-sectional geometry of the structure is improved to build a new cuttlebone-like structure.The experimental results showed that the new cuttlebone-like structure with revised cross-section improved the mechanical properties of the cuttlebone-like structure.The study results are of theoretical guidance to enhance the mechanical properties and geometric design of cuttlebone-like structures.

Fire is a significant threat to human life and civil infrastructures. Builders and architects are hankering for safer and sustainable alternatives of concrete that do not compromise with their design intent or fire safety requirements. The aim of the present work is to improve the residual compressive performance of concrete in the post-fire exposure, by incorporating by-products from urban residues. Based on sustainability and circular economy motivations, the attention is focused on rubber tire fly ash, aged brick powder, and plastic (PET) bottle residuals used as partial sand replacement. The selected by-products from urban residues are used for the preparation of Cement-Based Composites (CBCs) in two different proportions (10 % and 15 %). The thermal CBC behaviour is thus investigated under realistic fire scenarios (i.e., Direct Flame (DF) for 1 hour (1h)), by following the ISO 834 standard provisions, but necessarily resulting in non-uniform thermal exposure for the cubic specimens. The actual thermal exposure is further explored with a Finite Element (FE) model, giving evidence of thermal boundaries effects. The post-fire residual compressive strength of heated concrete and CBC samples is hence experimentally derived, and compared to unheated specimens in ambient conditions. From the experimental study, the enhanced post-fire performance of CBCs with PET bottle residual is generally found superior to other CBCs or concrete. The structure-property relation is also established, with the support of Scanning Electron Microscopy (SEM) micrographs. Based on existing empirical models of literature for the prediction of the compressive or residual compressive strength of standard concrete, newly developed empirical relations for both concrete and CBCs are assessed.

Rapid post-earthquake damage diagnosis of bridges can guide decision-making for emergency response management and recovery. This can be facilitated using digital technologies to remove the barriers of manual post-event inspections. Prior mechanics-based Finite Element (FE) models can be used for post-event response simulation using the measured ground motions at nearby stations; however, the damage assessment outcomes would suffer from uncertainties in structural and soil material properties, input excitations, etc. For instrumented bridges, these uncertainties can be reduced by integrating sensory data with prior models through a model updating approach. This study presents a sequential Bayesian model updating technique, through which a linear/nonlinear FE model, including soil-structure interaction effects, and the foundation input motions are jointly identified from measured acceleration responses. The efficacy of the presented model updating technique is first examined through a numerical verification study. Then, seismic data recorded from the San Rogue Canyon Bridge in California are used for a real-world case study. Comparison between the free-field and the foundation input motions reveals valuable information regarding the soil-structure interaction effects at the bridge site. Moreover, the reasonable agreement between the recorded and estimated bridge responses shows the potentials of the presented model updating technique for real-world applications. The updated FE model is considered as the digital twin of the bridge and can be used to analyze the bridge and monitor the structural response at element, section, and fiber levels to diagnose the location and severity of any potential damage mechanism.

In the process of metal cutting, the anti-wear performance of the tool determines the life of the tool and affects the surface quality of the workpiece. The finite element simulation method can directly show the tool wear state and morphology, but due to the limitations of the simulation time and complex boundary conditions, it has not been commonly used in tool life prediction. Based on this, a tool wear model was established on the platform of a finite element simulation software for the cutting process of titanium alloy TC4 by end milling. The key technique is to embed different types of tool wear models into the finite element model in combination with the consequent development technology. The effectiveness of the tool wear model is validated by comparing the experimental results with the simulation results. At the same time, in order to quickly predict the tool life, an empirical prediction formula of tool wear was established, and the change course of tool wear under time change was obtained.

Corrosion is considered as one of the main factors in the structural performance deterioration of steel members. In this study, experimental and numerical methods were used to assess the reduction in compressive strength of short tubular steel columns with local corrosion damage. The corrosion damage was varied with different depths (0, 1.5, 2, 3, 4, 4.5, and 6 mm), height (0, 20, 40, 60, 80, 100, 120, 140, 160, and 180 mm), circumference (0, 90, 180, 270, and 360°), and location along the column. A parametric numerical study was performed to establish a correlation between the residual compressive strength and the severity of corrosion damage. The results showed that as the corrosion depth, height and circumference increased, the compressive strength decreased linearly. As for the corrosion height, the residual compressive strength became constant after decreasing linearly when the corrosion height was greater than the half-wavelength of buckling of the short columns. An equation is presented to evaluate the residual compressive strength of short columns with local corrosion wherein the volume of the corrosion damage was used as a reduction factor in calculating the compressive strength. The percentage error using the presented equation was found to be within 11.4%.

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.

In this paper a new monolithic Eulerian formulation in the framework of non-classical continuum is presented for the analysis of fluid-strucutre interaction problems. In this respect, Cosserat continuum theory taking into account the micro-rotational degrees of freedom of the particles is considered. Continuum description of the model and variational formulation of the governing flow dynamics for non-classical -fluid-structure interaction nCFSI is presented. The model is analyzed by computing a well known benchmark problem by Hecht and Pironneau [16]. The algorithmic description is presented and implemented using FreeFEM++. Code is validated with the benchmark solution of Turek and Hron [38] in case of flow around a flag attached with cylinder. New microstructral behavior of the solution is studied and numerical simulations and results are shown in the form of figures. Some interesting feature of the flow is observed and microstructural characteristics are discussed.

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.

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.

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.

This paper studies an effective finite difference scheme for solving two-dimensional Heston stochastic volatility option pricing model problems. A dynamically balanced up-downwind strategy for approximating the cross-derivative is implemented and analyzed. Semi-discretized and spatially nonuniform platforms are utilized. The numerical method comprised is simple, straightforward with reliable first order overall approximations. The spectral norm is used throughout the investigation and numerical stability is proven. Simulation experiments are given to illustrate our results.

This work presents a linear smoothing scheme over high-order triangular elements in the framework of a cell-based strain smoothed finite element method for two-dimensional nonlinear problems. The main idea behind the proposed linear smoothing scheme for strain-smoothed finite element method (S-FEM) is no subdivision of finite element cells to sub-cells while the classical S-FEM needs sub-cells. Since the linear smoothing function is employed, S-FEM is able to use quadratic triangular or quadrilateral elements. The modified smoothed matrix obtained node-wise is evaluated. In the same manner with the computation of the strain-displacement matrix, the smoothed stiffness matrix and deformation graident are obtained over smoothing domains. A series of benchmark tests are investigated to demonstrate validity and stability of the proposed scheme. The validity and accuracy are confirmed by comparing the obtained numerical results with the standard FEM using 2nd-order triangular element and the exact solutions.

Polydimethylsiloxane (PDMS) is one of the most popular elastomers and has been used in different fields, especially in biomechanics research. Among the many interesting features of this material, its hyperelastic behaviour stands out, which allows the use of PDMS in various applications, like the ones that mimic soft tissues. However, the hyperelastic behaviour is not linear and needs detailed analysis, especially the characterization of shear strain. In this work, two approaches, numerical and experimental, were proposed to characterize the effect of shear strain on PDMS. The experimental method was implemented as a simple shear testing associated with 3D digital image correlation and was made using two specimens with two thicknesses of PDMS (2 and 4 mm). A finite element software was used to implement the numerical simulations, in which four different simulations using the Mooney-Rivlin, Yeoh, Gent, and Polynomial hyperelastic constitutive models were performed. These approaches showed that the maximum value of shear strain occurred in the central region of the PDMS, and higher values emerged for the 2 mm PDMS thickness. Qualitatively, in the central area of the specimen, the numerical and experimental results have similar behaviours and the values of shear strain are close. For higher values of displacement and thicknesses, the numerical simulation results move further away from experimental values.

EIT-MESHER (https://github.com/EIT-team/Mesher) is C++ software, based on the CGAL library, which generates high quality Finite Element Model tetrahedral meshes from binary masks of 3D volume segmentations. Originally developed for biomedical applications in Electrical Impedance Tomography (EIT) to address the need for custom, non-linear refinement in certain areas (e.g. around electrodes), EIT-MESHER can also be used in other fields where custom FEM refinement is required, such as Diffuse Optical Tomography (DOT).

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.

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.

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.

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.

Computation of Common Middle Point seismic sections and their subsequent time migration and diffraction imaging provides very important knowledge about the internal structure of 3D heterogeneous geological media and are key elements for successive geological interpretation. Full-scale numerical simulation, that computes all single shot seismograms, provides a full understanding of how the features of the image reflect the properties of the subsurface prototype. Unfortunately, this kind of simulations of 3D seismic surveys for realistic geological media needs huge computer resources, especially for simulation of seismic waves’ propagation through multiscale media like cavernous fractured reservoirs. In order to significantly reduce the query of computer resources we propose to model these 3D seismic cubes directly rather than the shot-by-shot simulation with subsequent CMP stacking. To do that we modify the well known "exploding reflectors principle" for 3D heterogeneous multiscale media by use of the finite-difference technique on the base of grids locally refined in time and space. To be able to simulate realistic models and acquisition we develop scalable parallel software, which needs reasonable computational costs. Numerical results for simulation of Common Middle Points sections and their time migration are presented and discussed.

Rubble stone masonry walls are widely diffused in most of the cultural and architectural heritage of historical cities. The mechanical response of such material is rather complicated to predict due to their composite nature. Vertical compression tests, diagonal compression tests, and shear-compression tests are usually adopted to experimentally investigate the mechanical properties of stone masonries. However, further tests are needed for the safety assessment of these ancient structures. Since the relation between normal and shear stresses plays a major role in the shear behavior of masonry joints, governing the failure mode, triplet test configuration was here investigated. First, the experimental tests carried out at the laboratory (LPMS) of the University of L'Aquila on stone masonry specimens were presented. Then, the triplet test was simulated by using the Total Strain Crack Model, which reflects all the ultimate states of quasi-brittle material such as cracking, crushing and shear failure. The goal of the numerical investigation was to evaluate the shear mechanical parameters of the masonry sample, including strength, dilatancy, normal and shear deformations. Furthermore, the effect of (i) confinement pressure and (ii) bond behavior at the sample-plates interfaces were investigated, showing that they can strongly influence the mechanical response of the walls.

A simplified nonlinear dispersive system of BBM-type, initially derived by D. Mitsotakis, is employed here in order to model the generation and propagation of surface water waves over variable bottom. The simplification consists in applying the so-called Boussinesq approximation. Using the finite element method and the FreeFem++ software, we solve numerically this system for three different complexities for the bathymetry function: a flat bottom case, a variable bottom in space, and a variable bottom both in space and in time. The last case is illustrated with the Java 2006 tsunami event. This article is designed rather as a tutorial paper even if it contains the description of completely new adaptation techniques.

The surface enhanced fluorescence（SEF）detection bases by plasmonic nanopillars array with nanoparticles has opened up a new gate in the application of biological imaging and sensing. The fluorescence enhancement of the probe molecule depends on its position in equilibrium, which is close to the hot spot leading to the electromagnetic field enhancement, but not too close to the metal surface resulting in quenching. Here, a large scale SiO₂-Ag-cicada wing SEF substrate was fabricated by magnetron sputtering with correction enhancement factor of 797.6. Thereinto the cicada wing provides the skeleton of the nanopillars array structure, the deposited Ag constructs two kinds of hot spots, and SiO₂ forms a separation layer to prevent quenching. Moreover, the substrate exhibited good reproducibility, high sensitivity with low limits of detection (LOD) and high stability for oxidation resistance. We propose that SEF substrate with modification of SiO₂ can not only improve the enhancement performance, but also expanding its application in the biological investigations.

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 paper is dedicated to the topical problem of examining permafrost state and the processes of its geocryological changes by means of geophysical methods. To monitor the cryolithozone, we propose and scientifically substantiate a new technique of pulsed electromagnetic cross-well sounding. Based on the vector finite-element method, we created a mathematical model of the cross-well sounding process with a pulsed source in a three-dimensional spatially heterogeneous medium. A high-performance parallel computing algorithm was developed and verified. Through realistic geoelectric models of permafrost with a talik under a highway, constructed following the results of electrotomography field data interpretation, we numerically simulated the pulsed sounding on the computing resources of the Siberian Supercomputer Center of SB RAS. The simulation results suggest the proposed system of pulsed electromagnetic cross-well monitoring to be characterized by a high sensitivity to the presence and dimensions of the talik. The devised approach can be oriented to addressing a wide range of issues related to monitoring permafrost rocks under civil and industrial facilities, buildings and constructions.

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 this article we explore the limiting behavior of the universal prior distribution obtained when applied over multiple meta-level hierarchy of programs and output data of a computational automata model. We were motivated to alleviate the effect of Solomonoff's assumption that all computable functions or hypotheses of the same length are equally likely, by weighing each program in turn by the algorithmic probability of their description number encoding. In the limiting case, the set of all possible program strings of a fixed-length converges to a distribution of self-replicating quines and quine-relays - having the structure of a constructor. We discuss how experimental algorithmic information theory provides insights towards understanding the fundamental metrics proposed in this work and reflect on the significance of these result in digital physics and the constructor theory of life.

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.

Cervical fusion has been a standard procedure for treating the abnormalities associated with the cervical spine. However, the reliability of anterior cervical discectomy and fusion (ACDF) has become arguable due to its adverse effects on the biomechanics of adjacent segments. One of the drawbacks associated with ACDF is adjacent segment degeneration (ASD) which has served as the base for the development of dynamic stabilization systems (DSS) and total disc replacement (TDR) devices for cervical spine. However, the hybrid surgical technique has also gained popularity recently but their effect on the biomechanics of cervical spine is not well researched. Thus, the objective of this FE study was to draw the comparison among single, bi-level and hybrid surgery with DCI implant with traditional fusion. Reduction in range of motion (ROM) for all the implanted models was observed for all the motions except extension, compared to intact model. The maximum increase in ROM of 42% was observed at C5-C6 level in Hybrid-DCI model. The maximum increase in adjacent segment’s ROM of 8.7% was observed in multilevel fusion model. The maximum von Mises stress in the implant was highest for the multilevel DCI model. Our study also showed that the shape of DCI implant permits flexion/extension relatively more compared to lateral bending and axial rotation.

The human mandible's cancellous bone, which is characterized by its unique porosity and directional sensitivity to external forces, is crucial for sustaining biting stress. Traditional CAD models fail to fully represent bone's anisotropic structure and thus depend on simple isotropic assumptions. For our research, we use the latest versions of nTopology and Creo Parametric software to make biomimetic Voronoi lattice models that accurately reflect the complex geometry and mechanical properties of trabecular bone. The porosity of human cancellous bone is accurately modeled in this work using biomimetic Voronoi lattice models. The porosities range from 70% to 95%, which can be achieved by changing the pore sizes to 1.0 mm, 1.5 mm, 2.0 mm, and 2.5 mm. Finite element analysis (FEA) was used to examine the displacements, stresses, and strains acting on dental implants with a buttress thread, abutment, retaining screw, and biting load surface. The results show that the Voronoi model accurately depicts the complex anatomy of the trabecular bone in the human jaw, compared to standard solid block models. The ideal pore size for biomimetic Voronoi lattice trabecular bone models is 2 mm, taking into account both the von Mises stress distribution over the dental implant, screw retention, cortical bone, cancellous bone, and micromotions. This pore size displayed balanced performance by successfully matching natural bone's mechanical characteristics. Advanced finite element analysis (FEA) improves biomechanical understanding of how bones and implants interact by creating more accurate models of biological problems and dynamic loading situations. This makes biomechanical engineering better.

Transition between free-surface and pressurized flows is an crucial phenomenon in many hydraulic systems, including water distribution systems, urban drainage systems, etc. During the transition, the force exerted on the structures changes drastically, thus it is meaningful to simulate this process. However, severe numerical oscillations are widely observed behind filling-bores, causing unphysical pressure variations and even computation failure. In this paper, some oscillation-suppressing approaches are reviewed and evaluated on a benchmark model. Then a new oscillation-suppressing approach is proposed to admit numerical viscosity when the water surface is at proximity of conduct roof which has first order accuracy. This approach adds numerical viscosity when water surface is at the proximity of conduct roof. It can sufficiently suppress numerical oscillations under an acoustic wave speed of 1000m/s and is simple to apply. In comparison with two experiments, the simulation results of this method show good agreement and little numerical oscillations. The results in this paper can help readers to choose an appropriate oscillation-suppressing method to improve the robustness and accuracy of flow regime transition simulations.

The aim of this work is to introduce a mathematical model representing the evolution of the temperature in a vegetation cover and the ground underneath it. Vegetation, and its interaction with soil, plays a very important role in the protection of soil surface from the action of sun and precipitations. A reduction in the vegetated mass increase the risk of desertification, soil erosion or surface runoffs which which can give rise to soil loss and sediment retention. These processes can favour climate change and global warming, which are major concerns nowadays. The mathematical model presented takes into account the main processes involved in vegetation cover and the interaction with the soil, among which, we can mention the Leaf Area Index, which is a dimensionless quantity defined as the one-sided green leaf area per unit ground surface area, or albedo and co-albedo which are clearly influenced by the vegetation. It is also considered a nonlinear heat capacity in the soil which incorporates the latent heat of fusion, when the phase change takes place. The numerical technique used to solve the mathematical model is based on a finite volume scheme with Weighted Essentially Non Oscillatory technique for spatial reconstruction and the third order Runge-Kutta Total Variation Diminishing numerical scheme is used for time integration. Some numerical examples are solved to obtain the distribution of temperature both in the vegetation cover and the soil.

We present a finite element model updating technique for soil-structure system identification of the Millikan Library building using the seismic data recorded during the 2002 Yorba Linda earthquake. A detailed finite element (FE) model of the Millikan Library building is developed in OpenSees and updated using a sequential Bayesian estimation approach for joint parameter and input identification. A two-step system identification approach is devised. First, the fixed-base structural model is updated to estimate the structural model parameters (including effective elastic modulus of structural components, distributed floor mass, and Rayleigh damping parameters) and some uncertain components of the foundation-level motion. Then, the identified structural model is used for soil-structure model updating wherein the Rayleigh damping parameters, the stiffness and viscosity of the soil subsystem (modeled using a substructure approach), and the foundation input motions (FIMs) are estimated. The identified model parameters are compared with state-of-practice recommendations. While a specific application is made for the Millikan Library, the present work offers a framework for integrating large-scale FE models with measurement data for model inversion. By utilizing this framework for different civil structures and earthquake records, key structural model parameters can be estimated from the real-world recorded data, which can subsequently be used for assessing and improving, as necessary, state-of-the-art seismic analysis and structural modeling techniques. This paper presents an effort towards using real-world measurements for large-scale FE model updating in the soil and structure, uniform soil time domain for joint parameter and input estimation, and thus paves the way for future applications in system identification, health monitoring and diagnosis of civil structures.

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 the present study, a mathematical model of non-steady partial differential equation from the process of oxygen mass transport in the human pulmonary circulation is proposed. Mathematical modelling of this kind of problems lead to a non-steady partial differential equation and for its numerical simulation, we have used finite differences. The aim of the process is the exact numerical analysis of the study, wherein consistency, stability and convergence is proposed. The necessity of doing the process is that, we would like to increase the order of numerical solution to a higher order scheme. An increment in the order of numerical solution makes the numerical simulation more accurate, also makes the numerical simulation being more complicated. In addition, the process of numerical analysis of the study in this order of solution needs more research work.

The fifth Industrial revolution (I5.0) prioritizes resilience and sustainability, integrating cognitive cyber-physical systems and advanced technologies to enhance machining processes. Numerous research studies have been conducted to optimize machining operations by identifying and reducing the sources of uncertainty and estimating the optimal cutting parameters. Virtual modeling and Tool Condition Monitoring (TCM) methodologies have been developed to assess the cutting states during machining processes. With a precise estimation of cutting states, the safety margin necessary to deal with the uncertainties can be reduced, resulting in improved process productivity. This paper reviews the recent advances in high-performance machining systems with a focus on the cyber-physical models developed for the cutting operation of difficult-to-cut materials using cemented carbide tools. An overview of the literature and background on the advances in the offline and online process optimization approaches have been presented. Process optimization objectives such as tool life utilization, dynamic stability, enhanced productivity, improved machined part quality, and reduced energy consumption and carbon emissions are independently investigated for the offline and online optimization methods. Addressing the critical objectives and constraints prevalent in industrial applications, the paper explores the challenges and opportunities inherent in developing a robust cyber-physical optimization system

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.

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

Abstract Background Context Pelvic ring fractures are becoming more common in the aging population and can prove to be fatal, having mortality rates between 10% and 16%. Stabilization of these fractures is challenging and often require immediate internal fixation. Therefore, it is necessary to have a biomechanical understanding of the different fixation techniques for pelvic ring fractures. Methods A previously validated three-dimensional finite element model of the lumbar spine, pelvis, and femur was used for this study. A unilateral pelvic ring fracture was simulated by resecting the left side of the sacrum and pelvis. Five different fixation techniques were used to stabilize the fracture. A compressive follower load and pure moment was applied to compare different biomechanical parameters including range of motion (contralateral sacroiliac joint, L1-S1 segment, L5-S1 segment), and stresses (L5-S1 nucleus stresses, instrument stresses) between different fixation techniques. Results Transiliac -trans sacral screw fixation at S1 and S2 showed the highest stabilization for horizontal and vertical displacement at the sacral fracture site and reduction of contralateral sacroiliac joint for bending and flexion range of motion by 165% and 121%, respectively. DTSF (Double transiliac rod and screw fixation) model showed highest stabilization in horizontal displacement at the pubic rami fracture site, while the L5_PF_W_CC (L5-Ilium posterior screw fixation with cross connectors) and L5_PF_WO_CC ( L5-Ilium posterior screw fixation without cross connectors) showed higher rod stresses, reduced L1-S1 (approximately 28%), and L5-S1 (approximately 90%) range of motion. Conclusions Longer sacral screw fixations were superior in stabilizing sacral and contralateral sacroiliac joint range of motion. Lumbopelvic fixations displayed a higher degree of stabilization in the horizontal displacement compared to vertical displacement of pubic rami fracture, while also indicating the highest rod stresses. When determining the surgical approach for pelvic ring fractures, patient-specific factors should be accounted for to weigh the advantages and disadvantages for each technique.

The advancement of technology and the availability of specialized digital signal processing chips have made digital filter design and implementation more feasible in a variety of fields, including biomedical engineering. This paper makes two key contributions. First, it uses a genetic algorithm to optimize the coefficients of Finite Impulse Response (FIR) filters. Second, it conducts a case study on using genetic algorithms to optimize FIR filters for electrocardiogram (ECG) biomedical signal noise removal. The goal of the proposed filter design approach is to achieve the desired signal bandwidth while minimizing the side lobe level and eliminating unwanted signals using a genetic algorithm. The results of the comprehensive analysis of the impact of different genetic operators show that the genetic algorithm-based filter outperforms conventional filter designs.

The authors previously found a model of universal quantum computation by making use of the coset structure of subgroups of a free group $G$ with relations. A valid subgroup $H$ of index $d$ in $G$ leads to a \lq magic' state $\left|\psi\right\rangle$ in $d$-dimensional Hilbert space that encodes a minimal informationally complete quantum measurement (or MIC), possibly carrying a finite \lq contextual' geometry. In the present work, we choose $G$ as the fundamental group $\pi_1(V)$ of an exotic $4$-manifold $V$, more precisely a \lq small exotic' (space-time) $R^4$ (that is homeomorphic and isometric, but not diffeomorphic to the Euclidean $\mathbb{R}^4$). Our selected example, due to to S. Akbulut and R.~E. Gompf, has two remarkable properties: (a) it shows the occurence of standard contextual geometries such as the Fano plane (at index $7$), Mermin's pentagram (at index $10$), the two-qubit commutation picture $GQ(2,2)$ (at index $15$) as well as the combinatorial Grassmannian Gr$(2,8)$ (at index $28$) , (b) it allows the interpretation of MICs measurements as arising from such exotic (space-time) $R^4$'s. Our new picture relating a topological quantum computing and exotic space-time is also intended to become an approach of \lq quantum gravity'.

Functionally graded materials (FGM) have attracted considerable attention in the field of composite materials and rekindled interest in research on composite materials, due to their unique mechanical response achieved through material design and optimization. Compared to conventional composites, FGMs offer several advantages and exceptional properties, including improved deformation resistance, improved toughness, lightness properties, and excellent recoverability. This study focused on the production of functionally graded (FG) polymer materials by additive manufacturing (AM) method. FG structures were produced by the Fused Deposition Modeling (FDM) method using Acrylonitrile benzidine styrene (ABS) and Polylactic acid (PLA) materials and tensile tests were performed according to ASTM D638. The effects of different layer thicknesses, volume ratios, and total thicknesses on mechanical behavior were investigated. In addition, tensile tests were carried out by producing single-layer samples at certain volume ratios to create a numerical model with the Finite element method to verify the experimental data. As a result of this study, it is presented that the FG structure produced with FDM improves mechanical behavior.

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

This study employs a high-fidelity numerical framework to determine the plastic material flow patterns and temperature distributions that lead to void formation during friction stir welding (FSW), and to relate the void morphologies to the underlying alloy material properties and process conditions. Three aluminum alloys, viz., 6061-T6, 7075-T6, and 5053-H18 were investigated under varying traverse speeds. The choice of aluminum alloys enables investigation of a wide range of thermal and mechanical properties. The numerical simulations were validated using experimental observations of void morphologies in these three alloys. Temperatures, plastic strain rates, and material flow patterns are considered. The key results from this study are: (1) The predicted stir zone and void morphology are in good agreement with the experimental observations, (2) The temperature and plastic strain-rate maps in the steady-state process conditions show a strong dependency on the alloy type and traverse speeds, (3) The material velocity contours provide a good insight into the material flow in the stir zone for the FSW process conditions that result in voids as well as those that do not result in voids. The numerical model and the ensuing parametric studies presented in this work provide a framework for understanding material flow under different process conditions in aluminum alloys, and potentially in other alloys. Furthermore, the utility of the numerical model for making quantitative predictions and investigating different process parameters to reduce void formation is demonstrated.

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.

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

This paper presents a new approach of eddy current methods for determining the size of the co-axial hole in the metallic circular disk. In recent decades, for the air-cored sensor probe, the impedance change due to the presence of an infinite metal plate can be calculated by the Dodd-Deeds model. However, in practical measurements, the sample cannot match with the condition required - ‘infinite’, thus the Dodd-Deeds model could not be applied to the disk with finite size and certainly not a co-axial hole in the center. In this paper, a dual-constraint analytical method is proposed. That is, the upper and lower limits of the integration are substituted with specific values instead of the original 0 and . Besides, it is found that, once the outer radius of the disk is fixed (i.e. the lower limit of integration is fixed), the upper limit reduces linearly as the size of the coaxial hole increases. Both the FEM simulation and experiments have been carried out to validate this method. The radius of the hole can be estimated based on the dual-constraint integration feature.

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

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 paper presents an algorithm for the fully automatic mesh generation for the finite element analysis of ships and offshore structures. The quality requirements on the mesh generator are imposed by the acceptance criteria of the classification societies as well as the need to avoid shear locking when using low degree shell elements. The meshing algorithm will be generating quadrilateral dominated meshes (consisting of quads and triangles) and the mesh quality requirements mandate that quadrilaterals with internal angles close to $90\degree$ are to be preferred. The geometry is described by a dictionary containing points, rods, surfaces and openings. The first part of the proposed method consist of an algorithm to automatically clean the geometry. The corrected geometry is then meshed by the frontal Delaunay mesh generator as implemented in the gmsh package. We present a heuristic method to precondition the cross field of the frontal quadrilateral mesher. Also the influence of the order in which the plates are meshed will be explored as a preconditioning step.

We overview recent findings achieved in the field of model-driven development of additively manufactured porous materials for development of a new generation of bioactive implants for orthopedic applications. Porous structures produced of biocompatible titanium alloys by selective laser melting can present a promising material to design scaffolds with regulated mechanical properties and with capacity to be loaded with pharmaceutical products. Adjusting pore geometry, one could control elastic modulus and strength/fatigue properties of the engineered structures to be compatible with bone tissues, thus preventing the stress shield effect when replacing a diseased bone fragment. Adsorption of medicals by internal spaces would make it possible to emit the antibiotic and anti-tumor agents into surrounding tissues. We critically analyze the recent advances in the field featuring model design approaches, virtual testing of the designed structures, capabilities of additive printing of porous structures, biomedical issues of the engineered scaffolds and so on. A special attention is paid to highlight the current troubles in the field and the ways of their solutions.

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

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

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

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

Recently, more wind turbine systems are being installed in deep waters far from the coast. Several concepts of floating wind turbine systems (FWTSs) are developed, among which the semisubmersible platform, due to its applicability in different water depths, good hydrodynamic performance, and facility in the installation process, constitutes the most explored technology compared to the others. However, a significant obstacle to the industrialization of this technology is a design of a cost-effective FWTS, which can be achieved by optimizing the geometry, size and weight of the floating platform, along with the mooring system. This is only possible by selecting a method capable of accurately analyzing the FWTS coupled hydro-aero-structural dynamics at each design stage. Accordingly, this paper aims at providing a detailed overview of the most common coupled numerical and physical methods, including their basic assumptions, formulations, limitations, and costs, used for analyzing the dynamics of FWTSs, mainly those supported by a semisubmersible, to assist the choice of the most suitable method at each design phase of FWTSs. Finally, the article discusses possible future research to address challenges in modeling FWTSs dynamics that persist to date.

This investigation introduces an efficient method for assessing the reliability of semi-rigid steel frame connections in existing buildings, explicitly targeting the Top Angle and Seat Pad (TA-SP) semi-rigid connection commonly used in Australian cities' high-rise buildings. The moment capacity of these connections was examined through parameter power model testing, meticulously considering the bending, shear, and prying forces actions to determine the weakest elements in the connection system. A plastic hinge model was applied to ensure computational efficiency and accuracy in assessing the connections' ultimate strength. Comparatively, the finite element (F.E) method M−θ curves were introduced to scrutinise the nonlinear behaviour of the connection, providing validation for the predicted analytical approach. Further emphasis was placed on identifying the reliability of the weakest elements in the connection system, employing the First-Order and Second-Order Reliability Methods models to estimate their reliability index. The correlation between the model's failure probability and a percentile error was also explored, utilising Monte Carlo simulation reference for refined validation. Novel formulas and methodologies were introduced to compute the locked-chain-series system reliability, which accounts for the inherent material and geometric variability. The findings highlight the critical significance of the weakest element in determining the overall reliability of the existing TA-SP connection.

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

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.

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

This study presents a comprehensive investigation of entropy generation during natural convection in a porous medium with Casson fluid. The governing equations, including the momentum, energy, and entropy balance equations, are solved numerically using finite element method. Through the analysis, various parameters affecting entropy generation are investigated, such as the Casson fluid parameter, Radiation, Prandtl number, and Rayleigh number. The results indicate that the Casson fluid parameter significantly influences the flow and heat transfer characteristics, while the Darcy number and Rayleigh number control the intensity of natural convection. Moreover, the Prandtl number determines the relative significance of heat transfer compared to viscous effects.

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.