Search results

This paper aims to investigate the performance of two novel numerical methods, the face-based smoothed finite element method (FS-FEM) and the edge-based smoothed finite element method (ES-FEM), which employ linear tetrahedral elements, for the purpose of strength assessment of a high-speed train hollow axle.

The calculation of stress for the wheelset, comprising an axle and two wheels, is facilitated through the application of the European axle strength design standard. This standard assists in the implementation of loading and boundary conditions and is exemplified by the typical CRH2 high-speed train wheelset. To evaluate the performance of these two methods, a hollow cylinder cantilever beam is first used as a benchmark to compare the present methods with other existing methods. Then, the strength analysis of a real wheelset model with a hollow axle is performed using different numerical methods.

The results of deflection and stress show that FS-FEM and ES-FEM offer higher accuracy and better convergence than FEM using linear tetrahedral elements. ES-FEM exhibits a superior performance to that of FS-FEM using linear tetrahedral elements, showing accuracy and convergence close to FEM using hexahedral elements.

This study channels the novel methods (FS-FEM and ES-FEM) in the static stress analysis of a railway wheelset. Based on the careful testing of FS-FEM and ES-FEM, both methods hold promise as more efficient tools for the strength analysis of complex railway structures.

The purpose of this study is to investigate the interaction between magnetotactic bacteria and Fe₃O₄–water nanofluid (NF) in a wavy enclosure in the presence of 2D natural convection flow.

Uniform magnetic field (MF), Brownian and thermophoresis effects are also contemplated. The dimensionless, time-dependent equations are governed by stream function, vorticity, energy, nanoparticle concentration and number of bacteria. Radial basis function-based finite difference method for the space derivatives and the second-order backward differentiation formula for the time derivatives are performed. Numerical outputs in view of isolines as well as average Nusselt number, average Sherwood number and flux density of microorganisms are presented.

Convective mass transfer rises if any of Lewis number, Peclet number, Rayleigh number, bioconvection Rayleigh number and Brownian motion parameter increases, and the flux density of microorganisms is an increasing function of Rayleigh number, bioconvection Rayleigh number, Peclet number, Brownian and thermophoresis parameters. The rise in buoyancy ratio parameter between 0.1 and 1 and the rise in Hartmann number between 0 and 50 reduce all outputs average Nusselt, average Sherwood numbers and flux density of microorganisms.

This study implies the importance of the presence of magnetotactic bacteria and magnetite nanoparticles inside a host fluid in view of heat transfer and fluid flow. The limitation is to check the efficiency on numerical aspect. Experimental observations would be more effective.

In practical point of view, in a heat transfer and fluid flow system involving magnetite nanoparticles, the inclusion of magnetotactic bacteria and MF effect provide control over fluid flow and heat transfer.

This is a scientific study. However, this idea may be extended to sustainable energy or biofuel studies, too. This means that a better world may create better social environment between people.

The presence of magnetotactic bacteria inside a Fe₃O₄–water NF under the effect of a MF is a good controller on fluid flow and heat transfer. Since the magnetotactic bacteria is fed by nanoparticles Fe₃O₄ which has strong magnetic property, varying nanoparticle concentration and Brownian and thermophoresis effects are first considered.

This paper focuses on the unconditionally optimal error estimates of a linearized second-order scheme for a nonlocal nonlinear parabolic problem. The first step of the scheme is based on Crank–Nicholson method while the second step is the second-order BDF method.

A rigorous error analysis is done, and optimal L² error estimates are derived using the error splitting technique. Some numerical simulations are presented to confirm the study’s theoretical analysis.

Optimal L² error estimates and energy norm.

The goal of this research article is to present and establish the unconditionally optimal error estimates of a linearized second-order BDF finite element scheme for the reaction-diffusion problem. An optimal error estimate for the proposed methods is derived by using the temporal-spatial error splitting techniques, which split the error between the exact solution and the numerical solution into two parts, that is, the temporal error and the spatial error. Since the spatial error is not dependent on the time step, the boundedness of the numerical solution in L∞-norm follows an inverse inequality immediately without any restriction on the grid mesh.

The space fractional PDEs (SFPDEs) play an important role in the fractional calculus field. Proposing a high-order, stable and flexible numerical procedure for solving SFPDEs is the main aim of most researchers. This paper devotes to developing a novel spectral algorithm to solve the FitzHugh–Nagumo models with space fractional derivatives.

The fractional derivative is defined based upon the Riesz derivative. First, a second-order finite difference formulation is used to approximate the time derivative. Then, the Jacobi spectral collocation method is employed to discrete the spatial variables. On the other hand, authors assume that the approximate solution is a linear combination of special polynomials which are obtained from the Jacobi polynomials, and also there exists Riesz fractional derivative based on the Jacobi polynomials. Also, a reduced order plan, such as proper orthogonal decomposition (POD) method, has been utilized.

A fast high-order numerical method to decrease the elapsed CPU time has been constructed for solving systems of space fractional PDEs.

The spectral collocation method is combined with the POD idea to solve the system of space-fractional PDEs. The numerical results are acceptable and efficient for the main mathematical model.

This study aims to examine the effects of cross-flow and multiple jet impingement on conductive panel cooling performance when subjected to uniform magnetic field effects. The cooling system has double rotating cylinders.

Cross-flow ratios (CFR) ranging from 0.1 to 1, magnetic field strength (Ha) ranging from 0 to 50 and cylinder rotation speed (Rew) ranging from −5,000 to 5,000 are the relevant parameters that are included in the numerical analysis. Finite element method is used as solution technique. Radial basis networks are used for the prediction of average Nusselt number (Nu), average surface temperature of the panel and temperature uniformity effects when varying the impacts of cross-flow, magnetic field and rotations of the double cylinder in the cooling channel.

The effect of CFR on cooling efficiency and temperature uniformity is favorable. By raising the CFR to the highest value under the magnetic field, the average Nu can rise by up to 18.6%, while the temperature drop and temperature difference are obtained as 1.87°C and 3.72°C. Without cylinders, magnetic field improves the cooling performance, while average Nu increases to 4.5% and 8.8% at CR = 0.1 and CR = 1, respectively. When the magnetic field is the strongest with cylinders in channel at CFR = 1, temperature difference (ΔT) is obtained as 2.5 °C. The rotational impacts on thermal performance are more significant when the cross-flow effects are weak (CFR = 0.1) compared to when they are substantial (CFR = 1). Cases without a cylinder have the worst performance for both weak and severe cross-flow effects, whereas using two rotating cylinders increases cooling performance and temperature uniformity for the conductive panel. The average surface temperature lowers by 1.2°C at CFR = 0.1 and 0.5°C at CFR = 1 when the worst and best situations are compared.

The outcomes are relevant in the design and optimization-based studies for electric cooling, photo-voltaic cooling and battery thermal management.

This study aims to understand how the texture shape, number of textures and addition of nanoparticle additives in lubricants impact the dynamic characteristics of journal bearing by comparing six different texture shapes like triangle, chevron, arc, circle, rectangle and elliptical applied in pressure-increasing region under various geometrical and operating conditions.

The finite element method approach has been employed to solve governing Reynold’s equation, assuming iso-viscous Newtonian fluid, for computation of performance parameters like stiffness and damping coefficient, threshold speed, etc. By using a regression model, the impact of adding nanoparticles Al₂O₃ and CuO to the base lubricant on viscosity variation is calculated for selected temperature ranges and weight fractions of nanoparticles.

The arc-shaped texture with an area density of 28.27%, eccentricity ratio of 0.2 and texture depth of 0.6 exhibited 35.22% higher direct stiffness and 41.4% higher damping coefficient compared to the lowest value in the circle-shaped texture. Increasing the number of arc-shaped textures on the bearing surface with low area density led to declining stiffness and damping parameters. However, with nanoparticle additives, the arc-shaped texture further showed 10.75% and 8.11% improvement in stiffness and 9.99% and 4.87% enhancement in damping coefficient for Al₂O₃ and CuO, respectively, at 90 °C temperature and 0.5% weight fraction.

By understanding the influence of texture shapes on the dynamic characteristics, engineers can design bearings that exhibit improved stability and enhance overall performance.

This paper aims to focus on the finite element method in the frequency domain (FD-FEM) for the transient electric field in the non-sinusoidal steady state under the non-sinusoidal periodic voltage excitation.

Firstly, the boundary value problem of the transient electric field in the frequency domain is described, and the finite element equation of the FD-FEM is derived by Galerkin’s method. Secondly, the constrained electric field equation on the boundary in the frequency domain (FD-CEFEB) is also derived, which can solve the electric field intensity on the boundary and the dielectric interface with high accuracy. Thirdly, the calculation procedures of the FD-FEM with FD-CEFEB are introduced in detail. Finally, a numerical example of the press-packed insulated gate bipolar transistor under the working condition of the repetitive turn-on and turn-off is given.

The FD-CEFEB improves numerical accuracy of electric field intensity on the boundary and interfacial charge density, which can be achieved by modifying the existing FD-FEMs’ code in appropriate steps. Moreover, the proposed FD-FEM and the FD-CEFEB will only increase calculation costs by a little compared with the traditional FD-FEMs.

The FD-CEFEB can directly solve the electric field intensity on the boundary and the dielectric interface with high accuracy. This paper provides a new FD-FEM for the transient electric field in the non-sinusoidal steady state with high accuracy, which is suitable for combined insulation structure with a long time constant.

This study establishes a method for predicting ground vibrations caused by railway tunnels in unsaturated soils with spatial variability.

First, an improved 2.5D finite-element-method-perfect-matching-layer (FEM-PML) model is proposed. The Galerkin method is used to derive the finite element expression in the ub-pl-pg format for unsaturated soil. Unlike the ub-v-w format, which has nine degrees of freedom per node, the ub-pl-pg format has only five degrees of freedom per node; this significantly enhances the calculation efficiency. The stretching function of the PML is adopted to handle the unlimited boundary domain. Additionally, the 2.5D FEM-PML model couples the tunnel, vehicle and track structures. Next, the spatial variability of the soil parameters is simulated by random fields using the Monte Carlo method. By incorporating random fields of soil parameters into the 2.5D FEM-PML model, the effect of soil spatial variability on ground vibrations is demonstrated using a case study.

The spatial variability of the soil parameters primarily affected the vibration acceleration amplitude but had a minor effect on its spatial distribution and attenuation over time. In addition, ground vibration acceleration was more affected by the spatial variability of the soil bulk modulus of compressibility than by that of saturation.

Using the 2.5D FEM-PML model in the ub-pl-pg format of unsaturated soil enhances the computational efficiency. On this basis, with the random fields established by Monte Carlo simulation, the model can calculate the reliability of soil dynamics, which was rarely considered by previous models.

Fused Filament Fabrication (FFF) is an extrusion-based manufacturing process using fused thermoplastics. Despite its low cost, the FFF is not extensively used in high-value industrial sectors mainly due to parts' anisotropy (related to the deposition strategy) and residual stresses (caused by successive heating cycles). Thus, this study aims to investigate the process improvement and the optimization of the printed parts.

In this work, a meshless technique – the Radial Point Interpolation Method (RPIM) – is used to numerically simulate the viscoplastic extrusion process – the initial phase of the FFF. Unlike the FEM, in meshless methods, there is no pre-established relationship between the nodes so the nodal mesh will not face mesh distortions and the discretization can easily be modified by adding or removing nodes from the initial nodal mesh. The accuracy of the obtained results highlights the importance of using meshless techniques in this field.

Meshless methods show particular relevance in this topic since the nodes can be distributed to match the layer-by-layer growing condition of the printing process.

Using the flow formulation combined with the heat transfer formulation presented here for the first time within an in-house RPIM code, an algorithm is proposed, implemented and validated for benchmark examples.

1) The OE-MLPG penalty meshfree method is developed to solve cracked structure.(2) Smartening the numerical meshfree method by combining the particle swarm optimization (PSO) optimization algorithms and Voronoi computational geometric algorithm. (3). Selection of base functions, finding optimal penalty factor and distribution of appropriate nodal points to the accuracy of calculation in the meshless local Petrov–Galekrin (MLPG) meshless method.

Using appropriate shape functions and distribution of nodal points in local domains and sub-domains and choosing an approximation or interpolation method has an effective role in the application of meshless methods for the analysis of computational fracture mechanics problems, especially problems with geometric discontinuity and cracks. In this research, computational geometry technique, based on the Voronoi diagram (VD) and Delaunay triangulation and PSO algorithm, are used to distribute nodal points in the sub-domain of analysis (crack line and around it on the crack plane).

By doing this process, the problems caused by too closeness of nodal points in computationally sensitive areas that exist in general methods of nodal point distribution are also solved. Comparing the effect of the number of sentences of basic functions and their order in the definition of shape functions, performing the mono-objective PSO algorithm to find the penalty factor, the coefficient, convergence, arrangement of nodal points during the three stages of VD implementation and the accuracy of the answers found indicates, the efficiency of V-E-MLPG method with Ns = 7 and ß = 0.0037–0.0075 to estimation of 3D-stress intensity factors (3D-SIFs) in computational fracture mechanics.

The present manuscript is a continuation of the studies (Ref. [33]) carried out by the authors, about; feasibility assessment, improvement and solution of challenges, introduction of more capacities and capabilities of the numerical MLPG method have been used. In order to validate the modeling and accuracy of calculations, the results have been compared with the findings of reference article [34] and [35].