Search results

This study aims to treat a novel system of Volterra integro-differential equations with multiple delays and variable bounds, constituting a generic numerical method based on the matrix equation and a combinatoric-parametric Charlier polynomials. The proposed method utilizes these polynomials for the matrix relations at the collocation points.

Thanks to the combinatorial eligibility of the method, the functional terms can be transformed into the generic matrix relations with low dimensions, and their resulting matrix equation. The obtained solutions are tested with regard to the parametric behaviour of the polynomials with $\alpha$, taking into account the condition number of an outcome matrix of the method. Residual error estimation improves those solutions without using any external method. A calculation of the residual error bound is also fulfilled.

All computations are carried out by a special programming module. The accuracy and productivity of the method are scrutinized via numerical and graphical results. Based on the discussions, one can point out that the method is very proper to solve a system in question.

This paper introduces a generic computational numerical method containing the matrix expansions of the combinatoric Charlier polynomials, in order to treat the system of Volterra integro-differential equations with multiple delays and variable bounds. Thus, the method enables to evaluate stiff differential and integral parts of the system in question. That is, these parts generates two novel components in terms of unknown terms with both differentiated and delay arguments. A rigorous error analysis is deployed via the residual function. Four benchmark problems are solved and interpreted. Their graphical and numerical results validate accuracy and efficiency of the proposed method. In fact, a generic method is, thereby, provided into the literature.

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.

The purpose of this study is to present the effectiveness and robustness of a numerical algorithm based on the block Newton method for the nonlinear kinematic hardening rules adopted in modeling ductile materials.

Elastoplastic problems can be defined as a coupled problem of the equilibrium equation for the overall structure and the yield equations for the stress state at every material point. When applying the Newton method to the coupled residual equations, the displacement field and the internal variables, which represent the plastic deformation, are updated simultaneously.

The presented numerical scheme leads to an explicit form of the hardening behavior, which includes the evolution of the equivalent plastic strain and the back stress, with the internal variables. The features of the present approach allow the displacement field and the hardening behavior to be updated straightforwardly. Thus, the scheme does not have any local iterative calculations and enables us to simultaneously decrease the residuals in the coupled boundary value problems.

A pseudo-stress for the local residual and an algebraically derived consistent tangent are applied to elastic-plastic boundary value problems with nonlinear kinematic hardening. The numerical procedure incorporating the block Newton method ensures a quadratic rate of asymptotic convergence of a computationally efficient solution scheme. The proposed algorithm provides an efficient and robust computation in the elastoplastic analysis of ductile materials. Numerical examples under elaborate loading conditions demonstrate the effectiveness and robustness of the numerical scheme implemented in the finite element analysis.

The purpose of this study is to propose a precise and standardized strategy for numerically simulating vehicle aerodynamics.

Error sources in computational fluid dynamics were analyzed. Additionally, controllable experiential and discretization errors, which significantly influence the calculated results, are expounded upon. Considering the airflow mechanism around a vehicle, the computational efficiency and accuracy of each solution strategy were compared and analyzed through numerous computational cases. Finally, the most suitable numerical strategy, including the turbulence model, simplified vehicle model, calculation domain, boundary conditions, grids and discretization scheme, was identified. Two simplified vehicle models were introduced, and relevant wind tunnel tests were performed to validate the selected strategy.

Errors in vehicle computational aerodynamics mainly stem from the unreasonable simplification of the vehicle model, calculation domain, definite solution conditions, grid strategy and discretization schemes. Using the proposed standardized numerical strategy, the simulated steady and transient aerodynamic characteristics agreed well with the experimental results.

Building upon the modified Low-Reynolds Number k-e model and Scale Adaptive Simulation model, to the best of the authors’ knowledge, a precise and standardized numerical simulation strategy for vehicle aerodynamics is proposed for the first time, which can be integrated into vehicle research and design.

The fractional order HIV model has an important role in biological science. To study the HIV model in a better way, the model is presented with the help of Atangana- Baleanu operator which is in Caputo sense. Also, the characteristics of the solutions are described briefly with the help of the advance numerical techniques for the different values of fractional order derivatives. This paper aims to discuss the aforementioned objectives.

In this work, Adams-Bashforth method and Euler method are used to get the solution of the HIV model. These are the important numerical methods. The comparison results also are described with the physical meaning of the solutions of the model.

HIV model is analyzed under the view of fractional and AB derivative in Atangana-Baleanu-Caputo sense. The uniqueness of the solution is proved by using Banach Fixed point. The solution is derived with the help of Sumudu transform. Further, the authors employed fractional Adam-Bashforth method and Euler method to enumerate numerical results. The authors have used several values of fractional orders to present the outcomes graphically. The above calculations have been done with the help of MATLAB (R2016a). The numerical scheme used in the proposed study is valid and fruitful, and the same can be used to explore other real issues.

This investigation can be done for the real data sets.

This paper aims to express the solution of the HIV model in a better way with the effect of non-locality, this work is very useful.

In this work, HIV model is developed with the help of Atangana- Baleanu operator in Caputo sense. By using Banach Fixed point, the authors proved that the solution is unique. Also, the solution is presented with the help of Sumudu transform. The behaviors of the solutions are checked for different values of fractional order derivatives with the physical meaning with help of the Adam-Bashforth method and the Euler method.

The purpose of this study is to introduce a new numerical scheme with no stability condition and high-order accuracy for the solution of two-dimensional coupled groundwater flow and transport simulation problems with regular and irregular geometries and compare the results with widely acceptable programs such as Modular Three-Dimensional Finite-Difference Ground-Water Flow Model (MODFLOW) and Modular Three-Dimensional Multispecies Transport Model (MT3DMS).

The newly proposed numerical scheme is based on the method of lines (MOL) approach and uses high-order approximations both in space and time. Quintic B-spline (QBS) functions are used in space to transform partial differential equations, representing the relevant physical phenomena in the system of ordinary differential equations. Then this system is solved with the DOPRI5 algorithm that requires no stability condition. The obtained results are compared with the results of the MODFLOW and MT3DMS programs to verify the accuracy of the proposed scheme.

The results indicate that the proposed numerical scheme can successfully simulate the two-dimensional coupled groundwater flow and transport problems with complex geometry and parameter structures. All the results are in good agreement with the reference solutions.

To the best of the authors' knowledge, the QBS-DOPRI5 method is used for the first time for solving two-dimensional coupled groundwater flow and transport problems with complex geometries and can be extended to high-dimensional problems. In the future, considering the success of the proposed numerical scheme, it can be used successfully for the identification of groundwater contaminant source characteristics.

Reynolds-averaged Navier–Stokes (RANS) models often perform poorly in shock/turbulence interaction regions, resulting in excessive wall heat load and incorrect representation of the separation length in shockwave/turbulent boundary layer interactions. The authors suggest that this can be traced back to inadequate numerical treatment of the inviscid fluxes. The purpose of this study is an extension to the well-known Harten, Lax, van Leer, Einfeldt (HLLE) Riemann solver to overcome this issue.

It explicitly takes into account the broadening of waves due to the averaging procedure, which adds numerical dissipation and reduces excessive turbulence production across shocks. The scheme is derived based on the HLLE equations, and it is tested against three numerical experiments.

Sod’s shock tube case shows that the scheme succeeds in reducing turbulence amplification across shocks. A shock-free turbulent flat plate boundary layer indicates that smooth flow at moderate turbulence intensity is largely unaffected by the scheme. A shock/turbulent boundary layer interaction case with higher turbulence intensity shows that the added numerical dissipation can, however, impair the wall heat flux distribution.

The proposed scheme is motivated by implicit large eddy simulations that use numerical dissipation as subgrid-scale model. Introducing physical aspects of turbulence into the numerical treatment for RANS simulations is a novel approach.

The numerical methods are of great importance for approximating the solutions of a system of nonlinear singular ordinary differential equations. In this paper, the authors present the biorthogonal flatlet multiwavelet collocation method (BFMCM) as a numerical scheme for a class of system of Lane–Emden equations with initial or boundary or four-point boundary conditions.

The approach is involved in combining the biorthogonal flatlet multiwavelet (BFM) with the collocation method. The authors investigate the properties and procedure of the BFMCM for first time on this class of equations. By using the BFM and the collocation points, the method is constructed and it transforms the nonlinear differential equations problem into a system of nonlinear algebraic equations. The unknown coefficients of the assuming solution are determined by solving the obtained system. Additionally, convergence analysis and numerical stability of the suggested method are provided.

According to the attained results, the proposed BFMCM has more accurate results in comparison with results of other methods. The maximum absolute errors are calculated by using the BFMCM for comparison purposes provided.

The key desirable properties of BFMCM are its efficiency, simple applicability and minimizes errors. Therefore, the proposed method can be used to solve nonlinear problems or problems with singular points.

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.

The paper aims to expand the scope of application of the lattice Boltzmann method (LBM), especially in the field of aircraft engineering. The traditional LBM is usually applied to incompressible flows at a low Reynolds number, which is not sufficient to satisfy the needs of aircraft engineering. Devoted to tackling the defect, the paper proposes a developed LBM combining the subgrid model and the multiple relaxation time (MRT) approach. A multilayer adaptive Cartesian grid method to improve the computing efficiency of the traditional LBM is also employed.

The subgrid model and the multilayer adaptive Cartesian grid are introduced into MRT-LBM for simulations of incompressible flows at a high Reynolds number. Validated by several typical flow simulations, the numerical methods in this paper can efficiently study the flows under high Reynolds numbers.

Some numerical simulations for the lid-driven flow of cavity, flow around iced GLC305, LB606b and ONERA-M6 are completed. The paper presents the investigation results, indicating that the methods are accurate and effective for the separated flow after icing.

LBM is developed with the addition of the subgrid model and the MRT method. A numerical strategy is proposed using a multilayer adaptive Cartesian grid method and its treatment of boundary conditions. The paper refers to innovative algorithm developments and applications to the aircraft engineering, especially for iced wing simulations with flow separations.