This paper presents an approach to the Timoshenko beam theory (TBT) using the finite difference method (FDM). The Timoshenko beam theory covers cases associated with small deflections based on shear deformation and rotary inertia considerations. The FDM is an approximate method for solving problems described with differential equations. It does not involve solving differential equations; equations are formulated with values at selected points of the structure. In addition, the boundary conditions and not the governing equations are applied at the beam’s ends. The model developed in this paper consisted of formulating differential equations with finite differences and introducing additional points at the beam’s ends and at positions of discontinuity (concentrated loads or moments, supports, hinges, springs, brutal change of stiffness, spring-mass system, etc.). The introduction of additional points allowed us to apply the governing equations at the beam’s ends. Moreover, grid points with variable spacing were considered, the grid being uniform within beam segments. First-order, second-order, and vibration analyses of structures were conducted with this model. Furthermore, tapered beams were analyzed (element stiffness matrix, second-order analysis, vibration analysis). Finally, a direct time integration method (DTIM) was presented; the FDM-based DTIM enabled the analysis of forced vibration of structures, with damping taken into account. The results obtained in this paper showed good agreement with those of other studies, and the accuracy was increased through a grid refinement. Especially in the first-order analysis of uniform beams, the results were exact for uniformly distributed and concentrated loads regardless of the grid.

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.

To obtain an accurate and reproducible experimental results in microbial fuel cell (MFC), it is important to know ‘anode maturation biofilm’ to produce a stable and maximum performance. For this purpose, four single chamber MFCs were tested in this study. The linear sweep voltammetry (LSV) polarization tests illustrated that maximum power densities of three MFCs became stable after 9 weeks. Although there were variations afterwards, such variations were negligible. Average maximum power densities from the 9th to the 17th week were 2,990 mW/m2 (MFC-4), 2,983 mW/m2 (MFC-2), 2,368 mW/m2 (MFC-3) and 837 mW/m2 (MFC-1). Polarization resistance shows that MFC-1 had much larger anode resistance (36.6-85.4 Ω) than the other MFCs (1.7-11.6 Ω). Anodic cyclic voltammetry (CV) shows that current production increased over time and MFC-1 had much smaller current production (24.4 mA) than the other MFCs (31.0-34.9 mA) at 17th week. The increased current production indicates anode biofilm became more mature over time, but overall cell performance did not increased accordingly. Possibly due to the bad inoculation, MFC-1 showed the lowest performance. However, its performance was restored to the initial performance and anode resistance was reduced by 47% at 17th week. This study shows that the optimum anode maturation time is 9 weeks and that bioanode performance is a key factor for MFC performance. This study also shows than LSV polarization and CV tests are accurate and non-destructive measurement methods for diagnosing anode performance.

Einstein’s publication of the Special Theory of Relativity in 1905 gave birth to the space-time model as a fundamental arena of the universe. Time has become 4th dimension of space and theoretically one could move backward and forward in time. Hypothetical time travel has been a hot topic of scientific discussion since Kurt Gödel published his article about closed time-like curves (CTC) in 1949. He was aware that his CTC model confirms that there is no physical time. Physicists understood his model as a theoretical frame for time travel. Recent research confirmed that universal space is time-invariant. Motion happens only in space and time is the duration of motion. Nobody can travel in time because time has no physical existence. Time as duration is an emergent physical quantity entering existence in the process of measurement of the duration of material change with clocks. The universal space is utterly timeless. There is no physical time in the universe.

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

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

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.

Pacing strategy refers to the distribution of effort and speed throughout the race to achieve optimal performance. The study aims to understand whether the choice of pacing strategy in swimming depends on the length of competitions and how sex, age, and performance level influence this strategy. Participants were the finalists of the 800-meter and 1500-meter freestyle events at the elite and junior world championships in 2022-2023. Race outcomes and pacing parameters were compared between the two distances and across different groups of swimmers. Swimmers in both distances did not break world records. Pacing strategy generally followed a U-shape with significant differences in the frequency and duration of speed changes between the two distances. The 800-meter exhibited more frequent changes in acceleration while the 1500-meter events generally followed a more consistent time series pattern. There were differences in pacing strategies between males and females and between junior and elite swimmers. Swimmers closer to world records showed more consistent pacing patterns compared to those farther from records. The study suggests that pacing strategies are influenced by race distance, sex, age, and performance level. The research highlights the complex interplay between physiological and psychological factors that shape a swimmer's decision-making during a race.

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.

Recently, there has been, on the one hand, an increase in the complexity of digital device designs and, on the other hand, an increase in the requirements for the development time and the reliability of the designs. One of the directions of solving this problem is developing new techniques for designing digital devices.This paper proposes a new technique for designing digital devices based on finite state machines with datapath (FSMD), when the functioning of the device is described in the form of an algorithm state machine with datapath (ASMD) charts. The new technique is called ASMD-FSMD. Different digital device design techniques are compared to each other using design examples of a synchronous multiplier on field programmable gate array (FPGA). The efficiency of the ASMD-FSMD technique compared to the traditional approach in terms of area and performance was investigated. The ASMD-FSMD technique, compared to the traditional one, reduces the area from 28.6% to 39.7% and increases the speed for some designs to 17.6%. In addition, using the ASMD-FSMD technique significantly reduces design time and increases design reliability. In conclusion, recommendations for using the ASMD-FSMD technique are made.

Research has shown that emotional states may influence the subjective experience of time. The present study involved 111 healthy adult participants and examined the impact of emotional stimuli on time perception, using a time estimation task with four conditions: in one of the conditions, participants estimated the exposure time of a black screen, wherein the duration could vary be-tween 2 to 7 s (condition without emotional interference); in the other conditions, participants es-timated the exposure time (also ranging from 2 to 7 s) of images that would induce pleasant, un-pleasant, or neutral emotions, to which participants were instructed to focus their attention (pleasant, unpleasant, or neutral emotions). The results showed that participants estimated “shorter” time intervals (underestimation) for pleasant and unpleasant images compared to the condition without interference, but no significant differences were found between the neutral condition and the condition without interference. It is suggested that increased attention to an ex-ternal stimuli, regardless of its emotional valence, produces an underestimated perception of time.

We consider the problem of scheduling $n$ jobs with identical processing times and given release as well as delivery times on $m$ uniform machines. The goal is to minimize the makespan, i.e., the maximum full completion time of any job. This problem is well-known to have an open complexity status even if the number of jobs is fixed. We present a polynomial-time algorithm for the problem which is based on the earlier introduced algorithmic framework blesscmore (``branch less and cut more''). We extend the analysis of the so-called behavior alternatives developed earlier for the version of the problem with identical parallel machines and show how the earlier used technique for identical machines can be extended to the uniform machine environment if a special condition on the job parameters is imposed. The time complexity of the proposed algorithm is $O(\gamma m^2 n\log n)$, where $\gamma$ can be either $n$ or the maximum job delivery time $q_{\max}$. This complexity can even be reduced further by using a smaller number $\kappa<n$ in the estimation describing the number of jobs of particular types. However, this number $\kappa$ becomes only known when the algorithm has terminated.

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

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

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 reliability of concrete cable-stayed bridge structure is time variant. This is due to strength deterioration and time-variant mechanical and environmental loadings. Demands for keeping concrete cable-stayed bridge structure in functional and reliable states require maintenance interventions. Limited financial resources require optimal maintenance strategies. These strategies have to consider the effects of uncertainties not only on the structural reliability deterioration but also on the maintenance interventions. By analyzing variation of concrete cable-stayed concrete cable-stayed bridge liabilities and benefit expectation through different maintenance programs, on the condition of satisfying the concrete cable-stayed bridge performance, with the maximum expected benefit as the target to discuss repair and reinforcement strategies for the concrete cable-stayed bridges in service. The method takes into account the present performance and future expected benefit, optimizing the performance of the concrete cable-stayed bridge and the expected benefit through determination policies of concrete cable-stayed concrete cable-stayed bridge maintenance program. The proposed method also includes an advanced cost model which takes into account the interaction between costs of maintenance interventions and their effect on reliability. Examples are presented demonstrating the application of the proposed method to a stock of existing deteriorating concrete cable-stayed bridges under various maintenance scenarios.

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

This paper presents an approach to the Euler-Bernoulli beam theory (EBBT) using the finite difference method (FDM). The EBBT covers the case of small deflections, and shear deformations are not considered. The FDM is an approximate method for solving problems described with differential equations. The FDM does not involve solving differential equations; equations are formulated with values at selected points of the structure. Generally, the finite difference approximations are derived based on fourth-order polynomial hypothesis (FOPH) and second-order polynomial hypothesis (SOPH) for the deflection curve; the FOPH is made for the fourth and third derivative of the deflection curve while the SOPH is made for its second and first derivative. In addition, the boundary conditions and not the governing equations are applied at the beam’s ends. In this paper, the FOPH was made for all of the derivatives of the deflection curve, and additional points were introduced at the beam’s ends and positions of discontinuity (concentrated loads or moments, supports, hinges, springs, etc.). The introduction of additional points allowed us to apply the governing equations at the beam’s ends and to satisfy the boundary and continuity conditions. Moreover, grid points with variable spacing were also considered, the grid being uniform within beam segments. First-order analysis, second-order analysis, and vibration analysis of structures were conducted with this model. Furthermore, tapered beams were analyzed (element stiffness matrix, second-order analysis). Finally, a direct time integration method (DTIM) was presented. The FDM-based DTIM enabled the analysis of forced vibration of structures, with damping taken into account. The results obtained in this paper showed good agreement with those of other studies, and the accuracy was increased through a grid refinement. Especially in the first-order analysis of uniform beams, the results were exact for uniformly distributed and concentrated loads regardless of the grid. Further research will be needed to investigate polynomial refinements (higher-order polynomials such as fifth-order, sixth-order…) of the deflection curve; the polynomial refinements aimed to increase the accuracy, whereby non-centered finite difference approximations at beam’s ends and positions of discontinuity would be used.

In Einstein’s Special Relativity(SR), time is treated as one of the four dimensions within spacetime and is relative in nature. Typically, a photon clock serves as the illustrative device in the time dilation thought experiment, featuring a single photon bouncing between two mirrors. However, the present paper identifies a fundamental flaw in this thought experiment, hence challenging its validity. Contrary to SR, this paper asserts that time is a consequence of motion. Space and time should be treated differently as proposed by classical mechanics. Importantly, this stance does not challenge the concept of curved space in General Relativity (GR), but rather questions the idea of relative time, as employed in both SR and GR. The absence of time dilation also dispels the twin paradox, prompting a necessary reassessment of the second postulate of relativity to ensure a more precise depiction of relativistic effects.

There are different ways to analyse energy absorbance (EA) in the human auditory system. In previous research, we developed a complete finite element model (FEM) of the human auditory system. In this mentioned work, the external auditory canal (EAC), middle ear, and inner ear (spiral cochlea, vestibule, and semi-circular canals) were modelled based on human temporal bone histological sections. Multiple acoustic, structure and fluid-coupled analyses were conducted using the FEM to perform harmonic analyses in the 0.1–10 kHz range. Once the FEM had been validated with published experimental data, the FEM numerical results were used to calculate the EA or energy reflected (ER) by the tympanic membrane. This EA was also measured in clinical audiology tests which were used as a diagnostic parameter. A mathematical approach was developed to calculate the EA and ER, with numerical and experimental results showing adequate correlation up to 1 kHz. Another published FEM had adapted its boundary conditions to replicate experimental results. Here, we recalculated those numerical results by applying the natural boundary conditions of human hearing and found that the results almost totally agreed with our FEM. This boundary problem is frequent and problematic in experimental hearing test protocols: the more invasive they are, the more the results are affected. One of the main objectives of using FEMs is to explore how the experimental test conditions influence the results. Further work will still be required to uncover the relationship between the middle ear structure and EA to clarify how to best use FEMs. Moreover, the FEM boundary conditions must be more representative in future work to ensure their adequate interpretation.

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

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

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

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

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

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

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

Time, traditionally viewed as a linear, non-dynamic parameter, is re-envisioned in this study as a Hyperdimensional concept. This paper conducts a cross-disciplinary examination, critically analysing the conceptualization of time in classical mechanics, quantum mechanics, and cosmology to propose a ground breaking reconceptualisation that extends beyond conventional frameworks. In classical mechanics, time is perceived as an absolute, continuously progressing backdrop, largely independent of events. Quantum mechanics treats time as a static parameter that does not influence quantum states but provides a framework for their evolution. In cosmology, time is considered a dimension that emerges from the Big Bang and serves as a measure for the universe's expansion, yet it does not interact with the structural dynamics of the cosmos. Challenging the transformative insights of Einstein’s relativity, which merges time with spatial dimensions under extreme conditions, we advocate for a perspective that views time as a Hyperdimensional and universal constant. This perspective posits that time, despite its unique and intrinsic properties, does not dynamically interact with or alter physical phenomena. Instead, it underpins our understanding of phenomena across different scales—from the minutiae of quantum states to the macroscopic dynamics of cosmology—without direct causation or change. A critical examination of time dilation and relativistic assumptions reveals significant discrepancies in traditional interpretations, particularly in how they are applied across different physical contexts. Our findings challenge the uniform applicability of time dilation, suggesting that observed phenomena often attributed to relativistic effects might better be explained through non-relativistic mechanisms such as phase shifts or changes in wavelength rather than temporal dilation. By synthesizing insights from various scientific domains, we advocate for a unified theory that recognizes time as a fundamental, universal dimension that is conceptual and non-interactive. Our goal is to bridge existing gaps between diverse scientific interpretations and promote a more integrated, profound understanding of time's autonomous and intrinsic nature.

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 &gt;12 h to a few seconds. All investigated methods showed excellent agreement with the golden standard by showing a &gt;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.

This paper presents an approach to the Kirchhoff-Love plate theory (KLPT) using the finite difference method (FDM). The KLPT covers the case of small deflections, and shear deformations are not considered. The FDM is an approximate method for solving problems described with differential equations. The FDM does not involve solving differential equations; equations are formulated with values at selected points of the structure. Generally in the case of KLPT, the finite difference approximations are derived based on the fourth-order polynomial hypothesis (FOPH) and second-order polynomial hypothesis (SOPH) for the deflection surface. The FOPH is made for the fourth and third derivative of the deflection surface while the SOPH is made for its second and first derivative; this leads to a 13-point stencil for the governing equation. In addition, the boundary conditions and not the governing equations are applied at the plate edges. In this paper, the FOPH was made for all of the derivatives of the deflection surface; this led to a 25-point stencil for the governing equation. Furthermore, additional nodes were introduced at plate edges and at positions of discontinuity (continuous supports/hinges, incorporated beams, stiffeners, brutal change of stiffness, etc.), the number of additional nodes corresponding to the number of boundary conditions at the node of interest. The introduction of additional nodes allowed us to apply the governing equations at the plate edges and to satisfy the boundary and continuity conditions. First-order analysis, second-order analysis, buckling analysis, and vibration analysis of plates were conducted with this model. Moreover, plates of varying thickness and plates with stiffeners were analyzed. Finally, a direct time integration method (DTIM) was presented. The FDM-based DTIM enabled the analysis of forced vibration of structures, with damping taken into account. In first-order, second-order, buckling, and vibration analyses of rectangular plates, the results obtained in this paper were in good agreement with those of well-established methods, and the accuracy was increased through a grid refinement.

The present work deals with the development of finite element methodology for obtaining the stress distributions in thick cylindrical HK40 stainless steel pipe that carry high temperature fluids. The material properties and loading are assumed to be random variables. Thermal stresses that are generated along radial, axial and tangential directions are computed generally using analytical expressions which are very complex. To circumvent such an issue, the probability theory and mathematical statistics have been applied to many engineering problems which allows to determine the safety both quantitatively and objectively based on the concepts of reliability. Monte Carlo simulation methodology is used to study the probabilistic characteristics of thermal stresses which is used for estimating the probabilistic distributions of stresses against the variations arising due to material properties and load. A 2-D Probabilistic finite element code is developed in MATLAB and the deterministic solution is compared with ABAQUS solutions. The values of stresses that are obtained from the variation of elastic modulus are found to be low as compared to the case where the load alone is varying. The probability of failure of the pipe structure is predicted against the variations in internal pressure and thermal gradient. These finite element framework developments are useful for the life estimation of piping structures in high temperature applications and subsequently quantifying the uncertainties in loading and material properties.

The present investigation focuses on the impact of jet nozzle orientation on mixing time in a cylindrical tank. The aim is to identify nozzle positions that improve mixing performance and to elucidate the governing parameters that influence it. A water tank was employed for the experiment. The vertical inclination angle (α) and the horizontal inclination angle (β) of the jet nozzle determined the nozzle positions. Mixing time was determined using an inert tracer and spectrophotometry measurements. The findings show that the mixing time is significantly influenced by the position of the jet nozzle position. The accuracy of existing jet turbulence and the circulation models for the prediction of mixing time was evaluated for the different nozzle positions. Our results indicate that both models provide accurate predictions for the conventional centrally aligned (β = 0°), upward-pointing jet nozzle positions only (α > 0). For the other nozzle positions where β > 0° and at varying α, the data follow the same trends as the jet turbulence and circulation models, however, the proportionality constants vary. Shorter mixing times can be attributed principally to longer jet path lengths and therefore higher fluid entrainment and circulation, as well as higher dissipation rates per jet length squared. However, it is suspected that the 3-dimensional nature of the flow pattern generated in the tank also plays a non-negligible role since mixing is hindered when the nozzle points more towards the tank wall.

This paper introduces a low-tech capacitive micromachined ultrasonic transducer (CMUT) designed with low environmental footprint materials. The fabrication process involves a copper plate as the fixed bottom electrode, a polymer-based adhesive as a dielectric material, and an aluminum foil as the top electrode.Finite element simulations include studies of displacement, mechanical stresses, and eigenfrequency. Experimental measurements validate the device's electromechanical behavior, showing an eigenfrequency of 88.6 kHz and a displacement of 22 pm. The low-tech CMUT demonstrates potential for applications such as ultrasonic actuation and energy harvesting, offering simplicity, biocompatibility, and low environmental impact. While not directly ready for applications, these transducers provide hands-on experience with technology similar to high-performance silicon-based implementations. These low-tech MUTs are perfect practical case studies for teaching purposes, combining simulation and experimental validation.

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.

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

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

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

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

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

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

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

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

Physicists have employed various approaches to solve the “paradox” and offered a consistent conclusion that the traveller twin is younger than the homebody twin. However, some authors attribute the cause of differential aging to acceleration, while others believe that acceleration should not be regarded as the source of differential aging, but frame switch is. So far there is no agreement on this issue. Actually, acceleration and frame switch should not be regarded as the source of differential aging. The author will introduce a thought experiment in which two very long spacecraft are moving in opposite directions at a speed close to the speed of light to discuss the real cause of differential aging. The author found that the real cause of differential aging in the twin “paradox” is because of the asymmetry of destination choice. People always choose their destinations in the universe rather than in spacecraft. The destinations are always stationary relative to the homebody twin and moving relative to the traveller twin. Once we choose a destination in a reference frame, the distance is the proper length in this reference frame, but the distance will contract in another moving reference frame, thus the recording time will be shorter.

This addendum contains clarifications and a sharpening of some of the results on the VIT transform framework developed in [1]. The focus is on the right-coefficient and left- coefficient forms of the transform, the extraction of a first-order term from a left polynomial fraction, and the application to linear time-varying systems.

We have introduced a sign operator of energy, analogous to the operator helicity, but in the direction of what we call energy vector. However, this energy vector needs a time vector. To give physical senses to the components of such a time vector, we try to explain the time dilation in special relativity and try to relate the components of the time vector to the tunneling times when an electron crosses a potential barrier.

Climate change such as sea level rise, change in temperature, precipitation, and storminess are expected to impact significantly coastal lagoons. The nature and magnitude of these impacts are uncertain. The objective of the research is to determine the climate change impacts on mixing and circulation at Songkhla lagoon, Thailand. Songkhla lagoon is the largest lagoonal water resource in Thailand and Southeast Asia. The lagoon is a combined freshwater and estuarine complex of high productivity which represents an extraordinary combination of environmental resources believed to be unique in the region. This work is part of a Climate Change impact assessment framework. It is the validation phase (step 5) of the framework applying a case study. Delft 3D was used to simulate CC scenarios in the climate downscaling models, part of the previous framework steps. These results were compared to the current conditions to determine the main changes in mixing and circulation in the coastal lagoon. Three indicators were applied to quantify the impacts: flushing time, salinity intrusion and stratification. The results suggest an increase in water velocities at the inlet in future scenarios and a decrease of flushing time. Salinity and stratification showed more complex changes in futures scenarios.

A number of experimental tests of time orientability are described as well as clear experimental signatures from non time orientability (time reversal). Some tests are well known, while others are based on more recent theoretical work. Surprisingly, the results all suggest that time is not orientable at a microscopic level; even definitive tests are positive. At a microscopic level the direction of time can reverse and a consistent forward time direction cannot be defined. That is the conclusion supported by a range of well-known experiments. The conflict between quantum theory and local realism; electrodynamics with electric charges; and spin half transformation properties of fermions; can all be interpreted as evidence of time reversal. While particle-antiparticle annihilation provides a definitive test. It offers both a new view of space-time and an novel interpretation of quantum theory with the potential to unify classical and quantum theories.

Can physical objects be in time-symmetry? Physical objects can only exist in a medium that has physical attributes, which means this medium is a type of energy. Is time energy? This article will show that time is not energy and there is no possibility that physical objects could be in time-symmetry. Physical objects only can be in symmetry in the time-invariant space in which they exist. In this perspective time measured with clocks is the result of the observer’s measurement in the time-invariant space. The time-symmetry model is flawed.

Capacity Expansion Models (CEMs) are optimization models used for long-term energy planning on national to continental scale. They are typically computationally demanding, thus in need of simplification, where one such simplification is to reduce the temporal representation. This paper investigates how using representative periods to reduce the temporal representation in CEMs distorts results compared to a benchmark model of a full chronological year. The test model is a generic CEM applied to Europe, equipped with a novel formulation for storage in model versions with reduced temporal representation. We test the performance of reduced models at penetration levels of wind and solar of 90%. Three measures for accuracy are used: (i) system cost, (ii) total capacity mix and (iii) regional capacity. We find that: (i) the system cost is well represented (~5% deviation from benchmark) with as few as ten representative days, (ii) the capacity mix is in general fairly well (~20% deviation) represented with 50 or more representative days, and (iii) the regional capacity mix displays large deviations (>50%) from benchmark for as many as 250 representative days. We conclude that modelers should be aware of the error margins when presenting results on these three aspects.

Today, due to recent developments in technology, children devote plenty of time for screen viewing. However, its harmful effects are not yet clear. The purpose of present study was to examine the associations among screen viewing and sleep duration, and body mass index (BMI) in under-five years old children. This cross-sectional study was conducted with 322 under-five healthy children that were selected using multistage stratified cluster sampling method in 2017. The data that were gathered by time-use diary method were analyzed using Kolmogorov-Smirnov test, Spearman correlation tests, multiple linear regression analysis, one-way ANCOVA, two-way ANCOVA. There was a negative correlation between screen time and sleep duration (rs = -0.42, p = 0.00), positive correlation between screen time and BMI (rs = 0.38, p = 0.00) and sleep duration negatively correlated with BMI (rs = -0.22, p = 0.00). screen viewing was a predictive factor for both sleep duration (β = -0.26, p = 0.00) and BMI (β = -0.26, p = 0.00). screen viewing had a significant impact on sleep duration (4, 314) = 5.02, P = 0.001) and BMI (F (4, 314) = 1.16, P=0.298). Results of this study indicated that screen viewing is related to sleep duration and BMI in under-five children. furthermore, screen time has an impact on sleep duration and BMI scores of children. findings of our study suggest that sleep duration negatively is associated with BMI in under-five-year-old children.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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.

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.

In order to reveal the mechanism and influencing factors of high pressure jet cavitation of nozzle in submerged environment, this study focused on the evolutionary process of cavitation bubbles and combined finite volume method and mixed multi-phase flow model to analyze the cavitation, velocity distribution and experimental cavitation intensity of fishing net cleaning equipment. Results show that the cavitation inception, growth and collapse mainly occurred in the peripheral region of the flow field. Ring-shaped cavitation erosion zone appeared on the test sample target. A lot of small dense erosion pits were densely distributed in the ring-shaped erosion zone, erosion marks were observed in the center. The cavitation erosion intensity was greatly affected by the nozzle structure. As the diameter of nozzle increased from 0.6 mm to 1 mm, the maximum gas volume fraction increased by 8.53%. The nozzle outlet enlargement angle greatly increased the cavitation intensity. The nozzle with an outlet angle of 30° exhibited the optimal cavitation erosion performance . The cavitation volume fraction of the nozzle with short necking structure was slightly larger than that of the nozzle with long necking structure at the same level in the necking length rang of 3mm to 7mm. In terms of the influence of nozzle structure on the cavitation erosion effect, the nozzle diameter D ranked the first, followed by the outlet angle α, and the necking length L was at the last.

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

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

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

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

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

Concrete-filled steel columns have been extensively used in the world due to having all suitable characteristics of concrete and steel, more ductility, increasing concrete confinement using steel wall, large energy-absorption capacity and appropriate fire behavior. In present paper, concrete-filled steel square columns have been simulated under the influence of blast load using ABAQUS software. These responses will be compared for scaled distances based on the distance to source and weight of explosive material. As result, it can be seen that although concrete deformation has been restricted using steel tube, but inner layer of concrete has been seriously damaged and column displacement will be decreased by increasing scaled distance. We also concluded that concrete-filled steel columns have high ductility and blast resistance.

The current work reveals the fine tuning between stability zones and randomness of Pointwise Stationary Fluid Flow Approximation (PSFFA) model of the non-stationary queueing D/M/1 system. More specifically, this clearly provides more insights into developing a contemporary PSFFA theory that unifies non-stationary queueing theory with chaos theory and fields in both theoretical physics and chaotic systems. This opens new grounds for stability analysis of non-stationary queueing systems.

In this paper, the concept of absolute airspace is further elaborated, and the time transform effect, the frequency transform relation of clock timing, and the principle of Doppler effect of special relativistic light are analyzed and demonstrated. Based on the concept of absolute airspace, a new and accurate Doppler effect relation of light (electromagnetic wave) is proposed. It is proved that there is a synchronous relationship between the time shift effect and the frequency change of the timing inside the clock. When the light (electromagnetic wave) source is in motion, the Doppler effect includes two physical mechanisms: the frequency change caused by time transformation effect and the frequency change caused by relative motion. When the light source is in motion, the Doppler effect of the light there are some differences in different directions due to the rotation of the earth. The relation of Doppler effect of special relativistic light is only a highly approximate relation, and its theoretical basis, special relativity, is flawed in principle.

On the basis of the findings of this paper, the theories of quantum mechanics and general relativity can be expressed in a comprehensive theory. The events were surveyed length of the time, in the Euclidean six-dimensional space-time (R6=x1+x2+x3+t1+t2+t3). Previous studies have not studied the relationship between the origin of matter's inherent properties, fundamental constants, and quantum mechanics phenomenon with space-time geometry. This theory prediction the with certainty obtained probabilistic results of quantum mechanics. This prediction is on the basis of metric resulting from the Eccentricity of the Ellipse due to the object density, in the length of time. The cause and nature of dark matter and dark energy are the findings of this article.

Noise model selection criteria has a significant impact on identifying the stochastic noise proper-ties of any GNSS daily coordinate time series. The low-frequency random walk noise existing in these time series could lead to overestimation of the tectonic rate, so it is of great significance to accurately detect the random walk component. This study focuses on noise model estimation cri-terion (BIC_tp) derived from the AIC and the BIC by introducing 2π factors. It is more sensitive to abnormal steps (random jumps). Using observation data from 72 GNSS stations from 1992 to 2022 and simulated data, four combined noise models are used to explore the impacts of time se-ries lengths (ranging from2 to 24 years) and data loss (between 2% and 30%) on noise models and velocity estimation. The results show that as the time length increased, the selected optimal noise model, and the estimated uncertainty of the tectonic trend with different data gap, gradually con-verge. When the time length is short (less than 8 years), it could lead to the FNRWWN, FNWN, and PLWN models being mistakenly estimated as GGMWN models, thereby affecting the accura-cy of determining the station velocity parameters. When the time length is 12 years, the RW noise component is more probably detected, As the time length increases, the impact of RW on velocity uncertainty is weakened. Finally, we conclude that for a time series with a minimum time length of 12 years, both the selection of the optimal stochastic noise model and the estimation of the ve-locity parameters are reliable.

The gelation time tg necessary for a solution of functional (associating) molecules to reach its gel point after a temperature jump, or a sudden concentration change, is theoretically calculated as a function of the concentration and temperature for the molecules carrying the number f of functional groups. Quite generally, tg is given by the product of the relaxation time tR and a thermodynamic factor Q. For pairwise cross-linking, and cross-linking with a fixed multiplicity k, precise values of tg are derived. For cross-linking with variable multiplicity, kinetic equation of stepwise association is derived for the reactivity pk of the functional groups, on the basis of which tg and tR are estimated. The thermodynamic factor Q is shown to generate a singularity of logarithmic divergence as the temperature (concentration) approaches the equilibrium gel point, while the relaxation time tR changes continuously across it. The retardation effect on the gelation time due to the reversibility of the cross-linking is explicitly calculated for some specific models of cross-linking. For a micellar cross-linking covering a wide range of the multiplicity, as seen in hydrophobically-modified water-soluble polymers, tR is shown to obey a formula similar to the Aniansson-Wall law.

We are using the book “Towards Quantum Gravity with an article by Claus Kiefer as to a quantum gravity interpretation of the density matrix in the early universe. The density matrix we are using is a one loop approximation, with inflaton value and potential terms, like V(phi) using the Padmanabhan values one can expect if the scale factor is a ~ a(Initial) times t ^ gamma , from early times . In doing so, we isolate out presuming a very small initial time step candidates initial time values which are from a polynomial for time values due to the Kiefer Density value.

In this paper, we study the Parisian time of a reflected Brownian motion with drift on a finite collection of rays. We derive the Laplace transform of the Parisian time using a recursive method, and provide an exact simulation algorithm to sample from the distribution of the Parisian time. The paper is motivated by the settlement delay in the real-time gross settlement (RTGS) system. Both the central bank and the participating banks in the system are concerned about the liquidity risk, and are interested in the first time that the duration of settlement delay exceeds a predefined limit, we reduce this problem to the calculation of the Parisian time. The Parisian time is also crucial in the pricing of Parisian type options; to this end, we will compare our results with the existing literature.

Future cyber-physical systems may extend over broad geographical areas, like cities or regions, thus requiring the deployment of large real-time networks. A strategy to guarantee predictable communication over such networks is to synthesize an offline time-triggered communication schedule. However, this synthesis problem is computationally hard (NP-complete), and existing approaches do not scale satisfactorily to the required network sizes. This article presents a segmented offline synthesis method which substantially reduces this limitation, being able to generate time-triggered schedules for large hybrid (wired and wireless) networks. We also present a series of algorithms and optimizations that increase the performance and compactness of the obtained schedules while solving some of the problems inherent to segmented approaches. We evaluate our approach on a set of realistic large-size multi-hop networks, significantly larger than those considered in the existing literature. The results show that our segmentation reduces the synthesis time up to two orders of magnitude.

In Newtonian mechanics, time as well as space are perceived as absolute entities. In 2 Einstein’s special relativity, time is frame dependent. Time is also affected by gravitational field 3 and as the field varies in space, time also varies throughout space. In the present article, a 4 thermodynamic-based time is investigated. The entity is called "irreversibility", which is generated 5 when availability (also known as exergy) is destroyed. Since each thermodynamic system may 6 generate different amount of irreversibility, this quantity is system dependent. The time’s arrow is 7 automatically satisfied, since irreversibility generation always proceeds in one direction (toward 8 future). We have demonstrated that, like common time, irreversibility is frame dependent, and 9 affected by gravity in the similar manner as the common time. For this reason, we propose to 10 assign the entity irreversibility of the system as thermodynamic time. A possible application of the 11 thermodynamic time is an interpretation and managing of the aging of biological systems. The 12 metabolic efficiency is related to the irreversibility of the chemical processes and affect the aging of 13 the system. Our sensation of time-flow may be attributed to the flow of availability and destruction 14 of it through the living system. It is shown by other authors that entropy generation (equivalent to 15 irreversibility) is a parameter for the human life span. Since the thermodynamic time is based on a 16 concept of thermodynamics that is universal, further applications to other subjects, such as biological 17 clock, telomere, and cosmology are possible.

With the great development of intelligent transportation systems (ITS), travel time prediction has attracted the attentions of many researchers and a large number of prediction methods have been developed. However, as an unavoidable topic, the predictability of travel time series is the basic premise for travel time prediction has received less attention than the methodology. Based on the analysis of the complexity of travel time series, this paper defines travel time predictability to express the probability of correct travel time prediction and proposes an entropy-based method to measure the upper bound of travel time predictability. Multiscale entropy is employed to quantify the complexity of travel time series, and the relationships between entropy and the upper bound of travel time predictability are presented. Empirical studies are made with vehicle trajectory data in an express road section. The effectiveness of time scales, tolerance, and series length to entropy and travel time predictability are analysis, and some valuable suggestions about the accuracy of travel time predictability are discussed. Finally, the comparisons between travel time predictability and actual prediction results from two prediction models, ARIMA and BPNN, are conducted. Experimental results demonstrate the validity and reliability of the proposed travel time predictability.

We explore the field equations in a 4-d complex space-time, in the same way, that general relativity does for our usual 4-d real space-time, forming this way, a new "general relativity" in C⁴ space-time, free of sources. Afterwards, by embedding our usual 4-d real space-time in C⁴ space-time, we describe geometrically the energy-momentum tensor T_μν as the lost geometric information of this embedding. We further give possible explanation of dark eld and dark energy.

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

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

Cycloid drives are generally used in precision machinery requiring high reduction ratios, such as robot joint (RV) reducers, in which the contact load and stress property of the internal meshing pairs greatly affect the dynamic transmission performance of the reducers. Therefore, in this study, Isogeometric Analysis (IGA) was employed to explore the multi-tooth contact problem of the cycloid pinwheel drive. Based on the NURBS curved surface generation method, the NURBS tooth profile of the cycloid gear was reconstructed. In addition, the NURBS surface of the cycloid gear - pin tooth - output pin was generated via the element splicing method. A geometrical analysis model of cycloid pinwheel drive was established to solve the contact force of the meshing pair under different input angles, and compared with the finite element method in terms of convergence, resultant accuracy, and solving timeliness. The results show that the Isogeometric analysis has higher accuracy and efficiency than the finite element method in calculating the contact stress and contact force. This method can improve the contact simulation accuracy of the cycloid drive and provides a reference for the design evaluation of RV reducer.

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

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

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

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

Emerging soft actuators with various soft materials and traveling-wave motion have been frequently discussed. Various configurations and the resulting performances have been proposed, and it is very challenging for soft actuators in realizing large strokes. This study presents an experimentally validated nonlinear finite element model to predict the deformation produced by a spiral pneumatic rubber actuator to generate traveling-wave motion. The actuator consists of a membrane mounted on a rubber substrate with three air chambers in a spiral configuration, simplifying the actuator's structure. The sequential deformations of the successive chambers interact with each other and produce radial traveling-waves on the membrane surface, driving the objects placed on the actuator. Finite element analysis using ANSYS computer software analyzed the elastic movement by considering the influence of different initial structural types. The structure with three different length ratios was introduced using the four parameters: membrane thickness, chamber width, chamber depth, and wall thickness. The simulation results indicated an optimal structure with specific ratios. A reasonable correlation was obtained during experimental validation; the predicted displacement values were approximately 17% smaller than the experimental values. Finally, the transportation performance of the prototype was tested and a velocity of 2.28 mm/s in the desired direction was achieved. We expect that our demonstration will expand the range of applications of the spiral pneumatic rubber actuator to include conveying or worm-like robots.

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

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

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