Search results

The main aim of the current paper is to find a numerical plan for hydraulic fracturing problem with application in extracting natural gases and oil.

First, time discretization is accomplished via Crank-Nicolson and semi-implicit techniques. At the second step, a high-order finite element method using quadratic triangular elements is proposed to derive the spatial discretization. The efficiency and time consuming of both obtained schemes will be investigated. In addition to the popular uniform mesh refinement strategy, an adaptive mesh refinement strategy will be employed to reduce computational costs.

Numerical results show a good agreement between the two schemes as well as the efficiency of the employed techniques to capture acceptable patterns of the model. In central single-crack mode, the experimental results demonstrate that maximal values of displacements in x- and y- directions are 0.1 and 0.08, respectively. They occur around both ends of the line and sides directly next to the line where pressure takes impact. Moreover, the pressure of injected fluid almost gained its initial value, i.e. 3,000 inside and close to the notch. Further, the results for non-central single-crack mode and bifurcated crack mode are depicted. In central single-crack mode and square computational area with a uniform mesh, computational times corresponding to the numerical schemes based on the high order finite element method for spatial discretization and Crank-Nicolson as well as semi-implicit techniques for temporal discretizations are 207.19s and 97.47s, respectively, with 2,048 elements, final time T = 0.2 and time step size τ = 0.01. Also, the simulations effectively illustrate a further decrease in computational time when the method is equipped with an adaptive mesh refinement strategy. The computational cost is reduced to 4.23s when the governed model is solved with the numerical scheme based on the adaptive high order finite element method and semi-implicit technique for spatial and temporal discretizations, respectively. Similarly, in other samples, the reduction of computational cost has been shown.

This is the first time that the high-order finite element method is employed to solve the model investigated in the current paper.

This study aims to derive a novel spatial numerical method based on multidimensional local Taylor series representations for solving high-order advection-diffusion (AD) equations.

The parabolic AD equations are reduced to the nonhomogeneous elliptic system of partial differential equations by utilizing the Chebyshev spectral collocation method (ChSCM) in the temporal variable. The implicit-explicit local differential transform method (IELDTM) is constructed over two- and three-dimensional meshes using continuity equations of the neighbor representations with either explicit or implicit forms in related directions. The IELDTM yields an overdetermined or underdetermined system of algebraic equations solved in the least square sense.

The IELDTM has proven to have excellent convergence properties by experimentally illustrating both h-refinement and p-refinement outcomes. A distinctive feature of the IELDTM over the existing numerical techniques is optimizing the local spatial degrees of freedom. It has been proven that the IELDTM provides more accurate results with far fewer degrees of freedom than the finite difference, finite element and spectral methods.

This study shows the derivation, applicability and performance of the IELDTM for solving 2D and 3D advection-diffusion equations. It has been demonstrated that the IELDTM can be a competitive numerical method for addressing high-space dimensional-parabolic partial differential equations (PDEs) arising in various fields of science and engineering. The novel ChSCM-IELDTM hybridization has been proven to have distinct advantages, such as continuous utilization of time integration and optimized formulation of spatial approximations. Furthermore, the novel ChSCM-IELDTM hybridization can be adapted to address various other types of PDEs by modifying the theoretical derivation accordingly.

This paper aims to comprehensively explore techniques for reducing solution time in finite element analysis (FEA), addressing the critical need for expediting computations to facilitate agile design exploration within project timelines.

Drawing from a wide array of literature sources, this paper synthesizes and analyzes various methodologies used to enhance the efficiency of FEA. Techniques are scrutinized in terms of their applicability, effectiveness and potential limitations.

The review signifies application of linear assumptions across multiple facets of analysis and delves into matrix order reduction strategies, geometry simplification, symmetry exploitation, submodeling and mesh attribute control. It reveals how these techniques can effectively reduce computational burdens while maintaining acceptable levels of accuracy.

While this review provides a comprehensive overview of existing efficiency enhancement techniques in FEA, it acknowledges inherent limitations of any synthesis-based study. Future research should focus on refining these methodologies.

The insights provided in this paper offer practical guidance for structural engineers and researchers seeking to optimize FEA workflows. By implementing these techniques, practitioners can expedite solution times and enhance their ability to explore design alternatives efficiently ultimately leading to cost savings and more robust structures.

This review contributes to the existing literature by offering a comprehensive synthesis of efficiency enhancement techniques in FEA. By highlighting the originality and value of each discussed methodology, this paper provides a roadmap for future research and practical implementation in the field of structural engineering.

A finite element computational methodology on a curved boundary using an efficient subparametric point transformation is presented. The proposed collocation method uses one-side curved and two-side straight triangular elements to derive exact subparametric shape functions.

Our proposed method builds upon the domain discretization into linear, quadratic and cubic-order elements using subparametric spaces and such a discretization greatly reduces the computational complexity. A unique subparametric transformation for each triangle is derived from the unique parabolic arcs via a one-of-a-kind relationship between the nodal points.

The novel transformation derived in this paper is shown to increase the accuracy of the finite element approximation of the boundary value problem (BVP). Our overall strategy is shown to perform well for the BVP considered in this work. The accuracy of the finite element approximate solution increases with higher-order parabolic arcs.

The proposed collocation method uses one-side curved and two-side straight triangular elements to derive exact subparametric shape functions.

The paper aims to the size-dependent analysis of functionally graded materials in thermal environment based on the modified couple stress theory using finite element method.

The element formulation is developed within the framework of the penalty unsymmetric finite element method (FEM) in that the C¹ continuity requirement is satisfied in weak sense and thus, C⁰ continuous interpolation enhanced by independent nodal rotation is employed as the test function. Meanwhile, the trial function is designed based on the stress functions and the weighted residual method. Besides, the special Gauss quadrature scheme is employed for integrals of matrices in accordance with the graded variation of the material properties.

The numerical results reveal that in thermal environment, functionally graded materials exhibit better bending performance compared to homogeneous materials, Moreover, the findings also indicate that with an increase in MLSP, the natural frequencies of out-of-plane modes gradually increase, while the natural frequencies of in-plane modes show much less variation, leading to a mode switch phenomenon.

The work provides an efficient numerical tool for analyzing and designing the functionally graded structures in thermal environment in practical engineering applications.

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.

The influence of the temperature discrepancy parameter and higher order of the time-derivative is discussed. Classical coupled and generalized hygrothermoelasticity models are recovered by considering the various special cases and illustrated graphically.

The theory of integral transformations has been used to study a new hygrothermal model that includes higher-order time derivatives with three-phase-lags and memory-dependent derivatives (MDD). This model considers the microscopic structure’s influence on a non-simple hygrothermoelastic infinitely long cylinder. The generalized Fourier and Fick’s law was adopted to derive the linearly coupled partial differential equations with higher-order time-differential with the two-phase lag model, including memory-dependent derivatives for the hygrothermal field. The investigation of microstructural interactions and the subsequent hygrothermal change has been undertaken as a result of the delay time and relaxation time translations.

These two-phase-lag models are also practically applicable in modeling nanoscale heat and moisture transport problems applied to almost all important devices. This work will enable future investigators to gain insight into non-simple hygrothermoelasticity with different phase delays of higher order in detail.

To the best of my knowledge, and after completing an intensive search of the relevant literature, the author could not learn any published research that presents a general solution for a higher-order time-fractional three-phase-lag hygrothermoelastic infinite circular cylinder with memory memory-dependent derivative.

The ability to characterise highly conducting objects, that may also be highly magnetic, by the complex symmetric rank–2 magnetic polarizability tensor (MPT) is important for metal detection applications including discriminating between threat and non-threat objects in security screening, identifying unexploded anti-personnel landmines and ordnance and identifying metals of high commercial value in scrap sorting. Many everyday non-threat items have both a large electrical conductivity and a magnetic behaviour, which, for sufficiently weak fields and the frequencies of interest, can be modelled by a high relative magnetic permeability. This paper aims to discuss the aforementioned idea.

The numerical simulation of the MPT for everyday non-threat highly conducting magnetic objects over a broad range of frequencies is challenging due to the resulting thin skin depths. The authors address this by employing higher order edge finite element discretisations based on unstructured meshes of tetrahedral elements with the addition of thin layers of prismatic elements. Furthermore, computer aided design (CAD) geometrical models of the non-threat and threat object are often not available and, instead, the authors extract the geometrical features of an object from an imaging procedure.

The authors obtain accurate numerical MPT characterisations that are in close agreement with experimental measurements for realistic physical objects. The assessment of uncertainty shows the impact of geometrical and material parameter uncertainties on the computational results.

The authors present novel computations and measurements of MPT characterisations of realistic objects made of magnetic materials. A novel assessment of uncertainty in the numerical predictions of MPT characterisations for uncertain geometry and material parameters is included.

In structural health monitoring, localization of multiple slight damage without baseline data is significant and difficult. The purpose of this paper is to discuss these issues.

Damage in the structure causes singularities of displacement modes, which in turn reveals damage. Methods based on the displacement modes may fail to accurately locate the slight damage because the slight damage in engineering structure results in a relatively small variation of the displacement modes. In comparison with the displacement modes, the strain modes are more sensitive to the slight damage because the strain is the derivative of the displacement. As a result, the slight variation in displacement data will be magnified by the derivative, leading to a significant variation of the strain modes. A novel method based on strain modes is proposed for the purpose of accurately locating the multiple slight damage.

In the two bay beam and steel fixed-fixed beams, the numerical simulations and the experimental cases, respectively, illustrate that the proposed method can achieve more accurate localization in comparison with the one based on the displacement modes.

The paper offers a practical approach for more accurate localization of multiple slight damage without baseline data. And the robustness to measurement noise of the proposed method is evaluated for increasing levels of artificially added white Gaussian noise until its limit is reached, defining its range of practical applicability.

The purpose of this study is to develop a new numerical algorithm to simulate the phase-field model.

First, the derivative of the temporal direction is discretized by a second-order linearized finite difference scheme where it conserves the energy stability of the mathematical model. Then, the isogeometric collocation (IGC) method is used to approximate the derivative of spacial direction. The IGC procedure can be applied on irregular physical domains. The IGC method is constructed based upon the nonuniform rational B-splines (NURBS). Each curve and surface can be approximated by the NURBS. Also, a map will be defined to project the physical domain to a simple computational domain. In this procedure, the partial derivatives will be transformed to the new domain by the Jacobian and Hessian matrices. According to the mentioned procedure, the first- and second-order differential matrices are built. Furthermore, the pseudo-spectral algorithm is used to derive the first- and second-order nodal differential matrices. In the end, the Greville Abscissae points are used to the collocation method.

In the numerical experiments, the efficiency and accuracy of the proposed method are assessed through two examples, demonstrating its performance on both rectangular and nonrectangular domains.

This research work introduces the IGC method as a simulation technique for the phase-field crystal model.