Search results

This study aims to discuss the numerical solutions of weakly singular Volterra and Fredholm integral equations, which are used to model the problems like heat conduction in engineering and the electrostatic potential theory, using the modified Lagrange polynomial interpolation technique combined with the biconjugate gradient stabilized method (BiCGSTAB). The framework for the existence of the unique solutions of the integral equations is provided in the context of the Banach contraction principle and Bielecki norm.

The authors have applied the modified Lagrange polynomial method to approximate the numerical solutions of the second kind of weakly singular Volterra and Fredholm integral equations.

Approaching the interpolation of the unknown function using the aforementioned method generates an algebraic system of equations that is solved by an appropriate classical technique. Furthermore, some theorems concerning the convergence of the method and error estimation are proved. Some numerical examples are provided which attest to the application, effectiveness and reliability of the method. Compared to the Fredholm integral equations of weakly singular type, the current technique works better for the Volterra integral equations of weakly singular type. Furthermore, illustrative examples and comparisons are provided to show the approach’s validity and practicality, which demonstrates that the present method works well in contrast to the referenced method. The computations were performed by MATLAB software.

The convergence of these methods is dependent on the smoothness of the solution, it is challenging to find the solution and approximate it computationally in various applications modelled by integral equations of non-smooth kernels. Traditional analytical techniques, such as projection methods, do not work well in these cases since the produced linear system is unconditioned and hard to address. Also, proving the convergence and estimating error might be difficult. They are frequently also expensive to implement.

There is a great need for fast, user-friendly numerical techniques for these types of equations. In addition, polynomials are the most frequently used mathematical tools because of their ease of expression, quick computation on modern computers and simple to define. As a result, they made substantial contributions for many years to the theories and analysis like approximation and numerical, respectively.

This work presents a useful method for handling weakly singular integral equations without involving any process of change of variables to eliminate the singularity of the solution.

To the best of the authors’ knowledge, the authors claim the originality and effectiveness of their work, highlighting its successful application in addressing weakly singular Volterra and Fredholm integral equations for the first time. Importantly, the approach acknowledges and preserves the possible singularity of the solution, a novel aspect yet to be explored by researchers in the field.

The purpose of this paper is to study large deviations for the solution processes of a stochastic equation incorporated with the effects of nonlocal condition.

A weak convergence approach is adopted to establish the Laplace principle, which is same as the large deviation principle in a Polish space. The sufficient condition for any family of solutions to satisfy the Laplace principle formulated by Budhiraja and Dupuis is used in this work.

Freidlin–Wentzell type large deviation principle holds good for the solution processes of the stochastic functional integral equation with nonlocal condition.

The asymptotic exponential decay rate of the solution processes of the considered equation towards its deterministic counterpart can be estimated using the established results.

This study aims to assess the robustness and accuracy of the face-centred finite volume (FCFV) method for the simulation of compressible laminar flows in different regimes, using numerical benchmarks.

The work presents a detailed comparison with reference solutions published in the literature –when available– and numerical results computed using a commercial cell-centred finite volume software.

The FCFV scheme provides first-order accurate approximations of the viscous stress tensor and the heat flux, insensitively to cell distortion or stretching. The strategy demonstrates its efficiency in inviscid and viscous flows, for a wide range of Mach numbers, also in the incompressible limit. In purely inviscid flows, non-oscillatory approximations are obtained in the presence of shock waves. In the incompressible limit, accurate solutions are computed without pressure correction algorithms. The method shows its superior performance for viscous high Mach number flows, achieving physically admissible solutions without carbuncle effect and predictions of quantities of interest with errors below 5%.

The FCFV method accurately evaluates, for a wide range of compressible laminar flows, quantities of engineering interest, such as drag, lift and heat transfer coefficients, on unstructured meshes featuring distorted and highly stretched cells, with an aspect ratio up to ten thousand. The method is suitable to simulate industrial flows on complex geometries, relaxing the requirements on mesh quality introduced by existing finite volume solvers and alleviating the need for time-consuming manual procedures for mesh generation to be performed by specialised technicians.

This paper aims to develop a novel grey Bernoulli model with memory characteristics, which is designed to dynamically choose the optimal memory kernel function and the length of memory dependence period, ultimately enhancing the model's predictive accuracy.

This paper enhances the traditional grey Bernoulli model by introducing memory-dependent derivatives, resulting in a novel memory-dependent derivative grey model. Additionally, fractional-order accumulation is employed for preprocessing the original data. The length of the memory dependence period for memory-dependent derivatives is determined through grey correlation analysis. Furthermore, the whale optimization algorithm is utilized to optimize the cumulative order, power index and memory kernel function index of the model, enabling adaptability to diverse scenarios.

The selection of appropriate memory kernel functions and memory dependency lengths will improve model prediction performance. The model can adaptively select the memory kernel function and memory dependence length, and the performance of the model is better than other comparison models.

The model presented in this article has some limitations. The grey model is itself suitable for small sample data, and memory-dependent derivatives mainly consider the memory effect on a fixed length. Therefore, this model is mainly applicable to data prediction with short-term memory effect and has certain limitations on time series of long-term memory.

In practical systems, memory effects typically exhibit a decaying pattern, which is effectively characterized by the memory kernel function. The model in this study skillfully determines the appropriate kernel functions and memory dependency lengths to capture these memory effects, enhancing its alignment with real-world scenarios.

Based on the memory-dependent derivative method, a memory-dependent derivative grey Bernoulli model that more accurately reflects the actual memory effect is constructed and applied to power generation forecasting in China, South Korea and India.

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.

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.

3D steady heat conduction analysis considering heat source is conducted on the fundamental of the fast multipole method (FMM)-accelerated line integration boundary element method (LIBEM).

Due to considering the heat source, domain integral is generated in the traditional heat conduction boundary integral equation (BIE), which will counteract the well-known merit of the BEM, namely, boundary-only discretization. To avoid volume discretization, the enhanced BEM, the LIBEM with dimension reduction property is introduced to transfer the domain integral into line integrals. Besides, owing to the unsatisfactory performance of the LIBEM when it comes to large-scale structures requiring massive computation, the FMM-accelerated LIBEM (FM-LIBEM) is proposed to improve the computation efficiency further.

Assuming N and M are the numbers of nodes and integral lines, respectively, the FM-LIBEM can reduce the time complexity from O(NM) to about O(N+ M), and a full discussion and verification of the advantage are done based on numerical examples under heat conduction.

(1) The LIBEM is applied to 3D heat conduction analysis with heat source. (2) The domain integrals can be transformed into boundary integrals with straight line integrals by the LIM. (3) A FM-LIBEM is proposed and can reduce the time complexity from O(NM) to O(N+ M). (4) The FM-LIBEM with high computational efficiency is exerted to solve 3D heat conduction analysis with heat source in massive computation successfully.

This paper aims to examine the uniform diffracted fields from a perfectly magnetic conductive (PMC) surface with the extended theory of boundary diffraction wave (BDW) approach.

Miyamoto and Wolf’s symbolic expression of the vector potential was used in the extended theory of BDW integral. This vector potential is applied to the problem, and the nonuniform field expression found was made uniform. Here, the expression is made uniform, using the detour parameter with the help of the asymptotic correlation of the Fresnel function. The BDW theory for the PMC surface extended the diffracted fields, and the uniform diffracted fields were calculated.

The field expressions obtained were interpreted with the graphs numerically for different aperture radii and observation distances. It has been shown that the BDW is continuous behind the diffracting aperture. There does not exist any discontinuity at the geometrically light-to-shadow transition boundary, as is required by the theory.

The results were graphically compared with diffracted fields for other surfaces. As far as we know, the uniform diffracted fields from the circular aperture on a PMC surface were calculated for the first time with the extended theory of the BDW approach.

The purpose of the present work is to introduce a wavelet method for the solution of linear and nonlinear psi-Caputo fractional initial and boundary value problem.

The authors have introduced the new generalized operational matrices for the psi-CAS (Cosine and Sine) wavelets, and these matrices are successfully utilized for the solution of linear and nonlinear psi-Caputo fractional initial and boundary value problem. For the nonlinear problems, the authors merge the present method with the quasilinearization technique.

The authors have drived the orthogonality condition for the psi-CAS wavelets. The authors have derived and constructed the psi-CAS wavelets matrix, psi-CAS wavelets operational matrix of psi-fractional order integral and psi-CAS wavelets operational matrix of psi-fractional order integration for psi-fractional boundary value problem. These matrices are successfully utilized for the solutions of psi-Caputo fractional differential equations. The purpose of these operational matrices is to make the calculations faster. Furthermore, the authors have derived the convergence analysis of the method. The procedure of implementation for the proposed method is also given. For the accuracy and applicability of the method, the authors implemented the method on some linear and nonlinear psi-Caputo fractional initial and boundary value problems and compare the obtained results with exact solutions.

Since psi-Caputo fractional differential equation is a new and emerging field, many engineers can utilize the present technique for the numerical simulations of their linear/non-linear psi-Caputo fractional differential models. To the best of the authors’ knowledge, the present work has never been introduced and implemented for psi-Caputo fractional differential equations.

The purpose of this study is to develop an efficient model to solve the electromagnetic forward problem using a novel semi-analytical approach to compute the electromagnetic fields because of the presence of a scatterer.

The proposed model involves a novel formulation of a complete orthonormal set of radiating/nonradiating polarization currents. Furthermore, an integral equation-based representation is derived, and the appropriate boundary conditions are imposed to get the scattered electromagnetic field. An error term is introduced to evaluate the obtained results.

The proposed model was tested using several examples at different frequencies. The results of this study show that the novel representation exhibits fast convergence behavior and achieves highly accurate results, when compared to the results provided by the transmission line method.

The derived formulations presented in this study are significant in the electromagnetic forward modelling field because of the meaningful physical representation they provide. This is an important aspect that leads to precise calculation of electromagnetic fields for various applications.