Search results

This work compares the performance of the three boundary element techniques for solving Helmholtz problems: dual reciprocity, multiple reciprocity and direct interpolation. All techniques transform domain integrals into boundary integrals, despite using different principles to reach this purpose.

Comparisons here performed include the solution of eigenvalue and response by frequency scanning, analyzing many features that are not comprehensively discussed in the literature, as follows: the type of boundary conditions, suitable number of degrees of freedom, modal content, number of primitives in the multiple reciprocity method (MRM) and the requirement of internal interpolation points in techniques that use radial basis functions as dual reciprocity and direct interpolation.

Among the other aspects, this work can conclude that the solution of the eigenvalue and response problems confirmed the reasonable accuracy of the dual reciprocity boundary element method (DRBEM) only for the calculation of the first natural frequencies. Concerning the direct interpolation boundary element method (DIBEM), its interpolation characteristic allows more accessibility for solving more elaborate problems. Despite requiring a greater number of interpolating internal points, the DIBEM has presented higher-quality results for the eigenvalue and response problems. The MRM results were satisfactory in terms of accuracy just for the low range of frequencies; however, the neglected higher-order primitives impact the accuracy of the dynamic response as a whole.

There are safe alternatives for solving engineering stationary dynamic problems using the boundary element method (BEM), but there are no suitable comparisons between these different techniques. This paper presents the particularities and detailed comparisons approaching the accuracy of the three important BEM techniques, aiming at response and frequency evaluation, which are not found in the specialized literature.

Hexahedral meshing is one of the most important steps in performing an accurate simulation using the finite element analysis (FEA). However, the current hexahedral meshing method in the industry is a nonautomatic and inefficient method, i.e. manually decomposing the model into suitable blocks and obtaining the hexahedral mesh from these blocks by mapping or sweeping algorithms. The purpose of this paper is to propose an almost automatic decomposition algorithm based on the 3D frame field and model features to replace the traditional time-consuming and laborious manual decomposition method.

The proposed algorithm is based on the 3D frame field and features, where features are used to construct feature-cutting surfaces and the 3D frame field is used to construct singular-cutting surfaces. The feature-cutting surfaces constructed from concave features first reduce the complexity of the model and decompose it into some coarse blocks. Then, an improved 3D frame field algorithm is performed on these coarse blocks to extract the singular structure and construct singular-cutting surfaces to further decompose the coarse blocks. In most modeling examples, the proposed algorithm uses both types of cutting surfaces to decompose models fully automatically. In a few examples with special requirements for hexahedral meshes, the algorithm requires manual input of some user-defined cutting surfaces and constructs different singular-cutting surfaces to ensure the effectiveness of the decomposition.

Benefiting from the feature decomposition and the 3D frame field algorithm, the output blocks of the proposed algorithm have no inner singular structure and are suitable for the mapping or sweeping algorithm. The introduction of internal constraints makes 3D frame field generation more robust in this paper, and it can automatically correct some invalid 3–5 singular structures. In a few examples with special requirements, the proposed algorithm successfully generates valid blocks even though the singular structure of the model is modified by user-defined cutting surfaces.

The proposed algorithm takes the advantage of feature decomposition and the 3D frame field to generate suitable blocks for a mapping or sweeping algorithm, which saves a lot of simulation time and requires less experience. The user-defined cutting surfaces enable the creation of special hexahedral meshes, which was difficult with previous algorithms. An improved 3D frame field generation method is proposed to correct some invalid singular structures and improve the robustness of the previous methods.

This paper aims to provide SIF calculation method for engineering application.

In this paper, the stress intensity factors (SIFs) calculation method is applied to the anisotropic Ni-based single crystal film cooling holes (FCHs) structure.

Based on contour integral, the anisotropic SIFs analysis finite element method (FEM) in Ni-based single crystal is proposed. The applicability and mesh independence of the method is assessed by comparing the calculated SIFs using mode of plate with an edge crack. Anisotropic SIFs can be calculated with excellent accuracy using the finite element contour integral approach. Then, the effect of crystal orientation and FCHs interference on the anisotropic SIFs is clarified. The SIFs of FCH edge crack in the [011] orientated Ni-based single crystal increases faster than the other two orientations. And the SIF of horizontal interference FCHs edge crack is also larger than that of the inclined interference one.

The SIFs of the FCH edge crack in the turbine air-cooled blade are innovatively computed using the sub-model method. Both the Mode I and II SIFs of FCHs edge crack in blade increase with crack growing.

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.

The ionization of the air surrounding the phase conductor in high-voltage transmission lines results in a phenomenon known as the Corona effect. To avoid this, Corona rings are used to dampen the electric field imposed on the insulator. The purpose of this study is to present a fast and intelligent surrogate model for determination of the electric field imposed on the surface of a 120 kV composite insulator, in presence of the Corona ring.

Usually, the structural design parameters of the Corona ring are selected through an optimization procedure combined with some numerical simulations such as finite element method (FEM). These methods are slow and computationally expensive and thus, extremely reducing the speed of optimization problems. In this paper, a novel surrogate model was proposed that could calculate the maximum electric field imposed on a ceramic insulator in a 120 kV line. The surrogate model was created based on the different scenarios of height, radius and inner radius of the Corona ring, as the inputs of the model, while the maximum electric field on the body of the insulator was considered as the output.

The proposed model was based on artificial intelligence techniques that have high accuracy and low computational time. Three methods were used here to develop the AI-based surrogate model, namely, Cascade forward neural network (CFNN), support vector regression and K-nearest neighbors regression. The results indicated that the CFNN has the highest accuracy among these methods with 99.81% R-squared and only 0.045468 root mean squared error while the testing time is less than 10 ms.

To the best of the authors’ knowledge, for the first time, a surrogate method is proposed for the prediction of the maximum electric field imposed on the high voltage insulators in the presence Corona ring which is faster than any conventional finite element method.

Two-dimensional (2D) problems are governed by unsteady anisotropic modified-Helmholtz equation of time–space dependent coefficients are considered. The problems are transformed into a boundary-only integral equation which can be solved numerically using a standard boundary element method (BEM). Some examples are solved to show the validity of the analysis and examine the accuracy of the numerical method.

The 2D problems which are governed by unsteady anisotropic modified-Helmholtz equation of time–space dependent coefficients are solved using a combined BEM and Laplace transform. The time–space dependent coefficient equation is reduced to a time-dependent coefficient equation using an analytical transformation. Then, the time-dependent coefficient equation is Laplace transformed to get a constant coefficient equation, which can be written as a boundary-only integral equation. By utilizing a BEM, this integral equation is solved to find numerical solutions to the problems in the frame of the Laplace transform. These solutions are then inversely transformed numerically to obtain solutions in the original time–space frame.

The main finding of this research is the derivation of a boundary-only integral equation for the solutions of initial-boundary value problems governed by a modified-Helmholtz equation of time–space dependent coefficients for anisotropic functionally graded materials with time-dependent properties.

The originality of the research lies on the time dependency of properties of the functionally graded material under consideration.

The assembly simulation in tolerance analysis is one of the most important steps for the tolerance design of mechanical products. However, most assembly simulation methods are based on the rigid body assumption, and those assembly simulation methods considering deformation have a poor efficiency. This paper aims to propose a novel efficient and precise tolerance analysis method based on stable contact to improve the efficiency and reliability of assembly deformation simulation.

The proposed method comprehensively considers the initial rigid assembly state, the assembly deformation and the stability examination of assembly simulation to improve the reliability of tolerance analysis results. The assembly deformation of mating surfaces was first calculated based on the boundary element method with optimal initial assembly state, then the stability of assembly simulation results was assessed by the density-based spatial clustering of applications with noise algorithm to improve the reliability of tolerance analysis. Finally, combining the small displacement torsor theory, the tolerance scheme was statistically analyzed based on sufficient samples.

A case study of a guide rail model demonstrated the efficiency and effectiveness of the proposed method.

The present study only considered the form error when generating the skin model shape, and the waviness and the roughness of the matching surface were not considered.

To the best of the authors’ knowledge, the proposed method is original in the assembly simulation considering stable contact, which can effectively ensure the reliability of the assembly simulation while taking into account the computational efficiency.

The need for precise synthesis of customized designs has resulted in the development of advanced coating processes for modern nanomaterials. Achieving accuracy in these processes requires a deep understanding of thermophysical behavior, rheology and complex chemical reactions. The manufacturing flow processes for these coatings are intricate and involve heat and mass transfer phenomena. Magnetic nanoparticles are being used to create intelligent coatings that can be externally manipulated, making them highly desirable. In this study, a Keller box calculation is used to investigate the flow of a coating nanofluid containing a viscoelastic polymer over a circular cylinder.

The rheology of the coating polymer nanofluid is described using the viscoelastic model, while the effects of nanoscale are accounted for by using Buongiorno’s two-component model. The nonlinear PDEs are transformed into dimensionless PDEs via a nonsimilar transformation. The dimensionless PDEs are then solved using the Keller box method.

The transport phenomena are analyzed through a comprehensive parametric study that investigates the effects of various emerging parameters, including thermal radiation, Biot number, Eckert number, Brownian motion, magnetic field and thermophoresis. The results of the numerical analysis, such as the physical variables and flow field, are presented graphically. The momentum boundary layer thickness of the viscoelastic polymer nanofluid decreases as fluid parameter increases. An increase in mixed convection parameter leads to a rise in the Nusselt number. The enhancement of the Brinkman number and Biot number results in an increase in the total entropy generation of the viscoelastic polymer nanofluid.

Intelligent materials rely heavily on the critical characteristic of viscoelasticity, which displays both viscous and elastic effects. Viscoelastic models provide a comprehensive framework for capturing a range of polymeric characteristics, such as stress relaxation, retardation, stretching and molecular reorientation. Consequently, they are a valuable tool in smart coating technologies, as well as in various applications like supercapacitor electrodes, solar collector receivers and power generation. This study has practical applications in the field of coating engineering components that use smart magnetic nanofluids. The results of this research can be used to analyze the dimensions of velocity profiles, heat and mass transfer, which are important factors in coating engineering. The study is a valuable contribution to the literature because it takes into account Joule heating, nonlinear convection and viscous dissipation effects, which have a significant impact on the thermofluid transport characteristics of the coating.

The momentum boundary layer thickness of the viscoelastic polymer nanofluid decreases as the fluid parameter increases. An increase in the mixed convection parameter leads to a rise in the Nusselt number. The enhancement of the Brinkman number and Biot number results in an increase in the total entropy generation of the viscoelastic polymer nanofluid. Increasing the strength of the magnetic field promotes an increase in the density of the streamlines. An increase in the mixed convection parameter results in a decrease in the isotherms and isoconcentration.

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.

In the present work, we focus on developing an in-house parallel meshless local Petrov-Galerkin (MLPG) code for the analysis of heat conduction in two-dimensional and three-dimensional regular as well as complex geometries.

The parallel MLPG code has been implemented using open multi-processing (OpenMP) application programming interface (API) on the shared memory multicore CPU architecture. Numerical simulations have been performed to find the critical regions of the serial code, and an OpenMP-based parallel MLPG code is developed, considering the critical regions of the sequential code.

Based on performance parameters such as speed-up and parallel efficiency, the credibility of the parallelization procedure has been established. Maximum speed-up and parallel efficiency are 10.94 and 0.92 for regular three-dimensional geometry (343,000 nodes). Results demonstrate the suitability of parallelization for larger nodes as parallel efficiency and speed-up are more for the larger nodes.

Few attempts have been made in parallel implementation of the MLPG method for solving large-scale industrial problems. Although the literature suggests that message-passing interface (MPI) based parallel MLPG codes have been developed, the OpenMP model has rarely been touched. This work is an attempt at the development of OpenMP-based parallel MLPG code for the very first time.