Search results

This paper aims to propose a new adaptive numerical method to find more accurate numerical solution for the heat source optimal control problem (OCP).

The main aim of this paper is to present an adaptive collocation approach based on the interpolating wavelets to solve an OCP for finding optimal heat source, in a two-dimensional domain. This problem arises when the domain is heated by microwaves or by electromagnetic induction.

This paper shows that combination of interpolating wavelet basis and finite difference method makes an accurate structure to design adaptive algorithm for such problems which usually have non-smooth solution.

The proposed numerical technique is flexible for different OCP governed by a partial differential equation with box constraint over the control or the state function.

This paper seeks to develop, propose and validate, through a series of presentable examples, a comprehensive high‐precision and ultra‐fast computing concept for solving stiff ordinary differential equations (ODEs) and partial differential equations (PDEs) with cellular neural networks (CNN).

The core of the concept developed in this paper is a straight‐forward scheme that we call “nonlinear adaptive optimization (NAOP)”, which is used for a precise template calculation for solving any (stiff) nonlinear ODEs through CNN processors.

One of the key contributions of this work (this is a real breakthrough) is to demonstrate the possibility of mapping/transforming different types of nonlinearities displayed by various classical and well‐known oscillators (e.g. van der Pol‐, Rayleigh‐, Duffing‐, Rössler‐, Lorenz‐, and Jerk‐ oscillators, just to name a few) unto first‐order CNN elementary cells, and thereby enabling the easy derivation of corresponding CNN‐templates. Furthermore, in case of PDEs solving, the same concept also allows a mapping unto first‐order CNN cells while considering one or even more nonlinear terms of the Taylor's series expansion generally used in the transformation of a PDEs in a set of coupled nonlinear ODEs. Therefore, the concept of this paper does significantly contribute to the consolidation of CNN as a universal and ultra‐fast solver of stiff differential equations (both ODEs and PDEs). This clearly enables a CNN‐based, real‐time, ultra‐precise, and low‐cost Computational Engineering. As proof of concept a well‐known prototype of stiff equations (van der Pol) has been considered; the corresponding precise CNN‐templates are derived to obtain precise solutions of this equation.

This paper contributes to the enrichment of the literature as the relevant state‐of‐the‐art does not provide a systematic and robust method to solve nonlinear ODEs and/or nonlinear PDEs using the CNN‐paradigm. Further, the “NAOP” concept developed in this paper has been proven to perform accurate and robust calculations. This concept is not based on trial‐and‐error processes as it is the case for various classes of optimization methods/tools (e.g. genetic algorithm, particle swarm, neural networks, etc.). The “NAOP” concept developed in this frame does significantly contribute to the consolidation of CNN as a universal and ultra‐fast solver of nonlinear differential equations (both ODEs and PDEs). An implantation of the concept developed is possible even on embedded digital platforms (e.g. field‐programmable gate array (FPGA), digital signal processing (DSP), graphics processing unit (GPU), etc.); this opens a broad range of applications. On‐going works (as outlook) are using NAOP for deriving precise templates for a selected set of practically interesting PDE models such as Navier Stokes, Schrödinger, Maxwell, etc.

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 purpose of this paper is to propose the use of the continuous adjoint method as a tool to identify the appropriate location and “type” (suction or blowing) of steady jets used in active flow control systems.

The method is based on continuous adjoint and covers both internal and external aerodynamics. The adjoint equations, including the adjoint to the SpalartAllmaras turbulence model and their boundary conditions are formulated. At the cost of solving the flow and adjoint equations just once, the sensitivity derivatives of the objective function with respect to hypothetical (normal) jet velocities at all wall nodes are computed. Comparisons of the computed sensitivities with finite differences and parametric studies to assess the present method are included.

Though the sensitivities are computed for zero jet velocities, they adequately support decision making on: the recommended location of jet(s), at boundary nodes with high absolute valued sensitivities; and the selection between suction or blowing jets, based on the sign of the computed sensitivities. Regarding adjoint methods, two important findings of this work are: the role of the adjoint pressure which proves to be an excellent sensor in flow control problems; and the prediction accuracy of the proposed adjoint method compared to the commonly made assumption of “frozen turbulence”.

First use of the continuous adjoint method using full differentiation of the turbulence model, in flow control optimization. A low‐cost design tool for recommending some of the most important jet characteristics.

The purpose of this paper is as follows: to significantly reduce the computation time (by a factor of 1,000 and more) compared to known numerical techniques for real-world problems with complex interfaces; and to simplify the solution by using trivial unfitted Cartesian meshes (no need in complicated mesh generators for complex geometry).

This study extends the recently developed optimal local truncation error method (OLTEM) for the Poisson equation with constant coefficients to a much more general case of discontinuous coefficients that can be applied to domains with different material properties (e.g. different inclusions, multi-material structural components, etc.). This study develops OLTEM using compact 9-point and 25-point stencils that are similar to those for linear and quadratic finite elements. In contrast to finite elements and other known numerical techniques for interface problems with conformed and unfitted meshes, OLTEM with 9-point and 25-point stencils and unfitted Cartesian meshes provides the 3-rd and 11-th order of accuracy for irregular interfaces, respectively; i.e. a huge increase in accuracy by eight orders for the new 'quadratic' elements compared to known techniques at similar computational costs. There are no unknowns on interfaces between different materials; the structure of the global discrete system is the same for homogeneous and heterogeneous materials (the difference in the values of the stencil coefficients). The calculation of the unknown stencil coefficients is based on the minimization of the local truncation error of the stencil equations and yields the optimal order of accuracy of OLTEM at a given stencil width. The numerical results with irregular interfaces show that at the same number of degrees of freedom, OLTEM with the 9-points stencils is even more accurate than the 4-th order finite elements; OLTEM with the 25-points stencils is much more accurate than the 7-th order finite elements with much wider stencils and conformed meshes.

The significant increase in accuracy for OLTEM by one order for 'linear' elements and by 8 orders for 'quadratic' elements compared to that for known techniques. This will lead to a huge reduction in the computation time for the problems with complex irregular interfaces. The use of trivial unfitted Cartesian meshes significantly simplifies the solution and reduces the time for the data preparation (no need in complicated mesh generators for complex geometry).

It has been never seen in the literature such a huge increase in accuracy for the proposed technique compared to existing methods. Due to a high accuracy, the proposed technique will allow the direct solution of multiscale problems without the scale separation.

This paper aims to discuss the approximate solution of the nonlinear thin film flow problems. A new analytic approximate technique for addressing nonlinear problems, namely, the optimal homotopy asymptotic method (OHAM), is proposed and used in an application to the nonlinear thin film flow problems.

This approach does not depend upon any small/large parameters. This method provides a convenient way to control the convergence of approximation series and to adjust convergence regions when necessary.

The obtained solutions show that the OHAM is more effective, simpler and easier than other methods. The results reveal that the method is explicit. By applying the method to nonlinear thin film flow problems, it was found to be simpler in applicability, and more convenient to control convergence. Therefore, the method shows its validity and great potential for the solution of nonlinear problems in science and engineering.

The proposed method is tested upon nonlinear thin film flow equation from the literature and the results are compared with the available approximate solutions including Adomian decomposition method (ADM), homotopy perturbation method, modified homotopy perturbation method and HAM. Moreover, the exact solution is compared with the available numerical solutions. The graphical representation of the solution is given by Maple and is physically interpreted.

The purpose of this paper is to present a numerical solution of a family of generalized fifth‐order Korteweg‐de Vries equations using a meshless method of lines. This method uses radial basis functions for spatial derivatives and Runge‐Kutta method as a time integrator and exhibits high accuracy as seen from the comparison with the exact solutions.

The study uses a meshless method of lines. This method uses radial basis functions for spatial derivatives and Runge‐Kutta method as a time integrator.

The paper reveals that this method exhibits high accuracy as seen from the comparison with the exact solutions.

This method is efficient method as it is easy to implement for the numerical solutions of PDEs.

This paper presents a novel framework for solving elliptic partial differential equations (PDEs) over irregularly spaced meshes on bounded domains.

Second‐generation wavelet construction gives rise to a powerful generalization of the traditional hierarchical basis (HB) finite element method (FEM). A framework based on piecewise polynomial Lagrangian multiwavelets is used to generate customized multiresolution bases that have not only HB properties but also additional qualities.

For the 1D Poisson problem, we propose – for any given order of approximation – a compact closed‐form wavelet basis that block‐diagonalizes the stiffness matrix. With this wavelet choice, all coupling between the coarse scale and detail scales in the matrix is eliminated. In contrast, traditional higher‐order (n>1) HB do not exhibit this property. We also achieve full scale‐decoupling for the 2D Poisson problem on an irregular mesh. No traditional HB has this quality in 2D.

Similar techniques may be applied to scale‐decouple the multiresolution finite element (FE) matrices associated with more general elliptic PDEs.

By decoupling scales in the FE matrix, the wavelet formulation lends itself particularly well to adaptive refinement schemes.

The paper explains second‐generation wavelet construction in a Lagrangian FE context. For 1D higher‐order and 2D first‐order bases, we propose a particular choice of wavelet, customized to the Poisson problem. The approach generalizes to other elliptic PDE problems.

Prediction of excess pore water pressure and estimation of soil parameters are the two key interests for consolidation problems, which can be mathematically quantified by a set of partial differential equations (PDEs). Generally, there are challenges in solving these two issues using traditional numerical algorithms, while the conventional data-driven methods require massive data sets for training and exhibit negative generalization potential. This paper aims to employ the physics-informed neural networks (PINNs) for solving both the forward and inverse problems.

A typical consolidation problem with continuous drainage boundary conditions is firstly considered. The PINNs, analytical, and finite difference method (FDM) solutions are compared for the forward problem, and the estimation of the interface parameters involved in the problem is discussed for the inverse problem. Furthermore, the authors also explore the effects of hyperparameters and noisy data on the performance of forward and inverse problems, respectively. Finally, the PINNs method is applied to the more complex consolidation problems.

The overall results indicate the excellent performance of the PINNs method in solving consolidation problems with various drainage conditions. The PINNs can provide new ideas with a broad application prospect to solve PDEs in the field of geotechnical engineering, and also exhibit a certain degree of noise resistance for estimating the soil parameters.

This study presents the potential application of PINNs for the consolidation of soils. Such a machine learning algorithm helps to obtain remarkably accurate solutions and reliable parameter estimations with fewer and average-quality data, which is beneficial in engineering practice.

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.