Search results

The purpose of this paper is to develop an efficient method for solving the magneto-hydrodynamic (MHD) boundary layer flow of an upper-convected Maxwell (UCM) fluid over a porous isothermal stretching sheet.

The paper applied a collocation approach based on rational Legendre functions for solving the third-order non-linear boundary value problem, describing the MHD boundary layer flow of an UCM fluid over a porous isothermal stretching sheet. This method solves the problem on the semi-infinite domain without transforming domain of the problem to a finite domain.

This approach reduces the solution of a problem to the solution of a system of algebraic equations. The numerical values of the skin friction coefficient are presented and analyzed for various parameters of interest in the problem. The authors also compare the results of this work with some recent results and show that the new method is efficient and applicable.

The method solves this problem without use of discrete variables and linearization or small perturbation. Also it was confirmed by the theorem and figure of absolute coefficients that this approach has exponentially convergence rate.

The finite element method (FEM) is used to calculate the two-dimensional anti-plane dynamic response of structure embedded in D’Alembert viscoelastic multilayered soil on the rigid bedrock. This paper aims to research a time-domain absorbing boundary condition (ABC), which should be imposed on the truncation boundary of the finite domain to represent the dynamic interaction between the truncated infinite domain and the finite domain.

A high-order ABC for scalar wave propagation in the D’Alembert viscoelastic multilayered media is proposed. A new operator separation method and the mode reduction are adopted to construct the time-domain ABC.

The derivation of the ABC is accurate for the single layer but less accurate for the multilayer. To achieve high accuracy, therefore, the distance from the truncation boundary to the region of interest can be zero for the single layer but need to be about 0.5 times of the total layer height of the infinite domain for the multilayer. Both single-layered and multilayered numerical examples verify that the accuracy of the ABC is almost the same for both cases of only using the modal number excited by dynamic load and using the full modal number of infinite domain. Using the ABC with reduced modes can not only reduce the computation cost but also be more friendly to the stability. Numerical examples demonstrate the superior properties of the proposed ABC with stability, high accuracy and remarkable coupling with the FEM.

A high-order time-domain ABC for scalar wave propagation in the D’Alembert viscoelastic multilayered media is proposed. The proposed ABC is suitable for both linear elastic and D’Alembert viscoelastic media, and it can be coupled seamlessly with the FEM. A new operator separation method combining mode reduction is presented with better stability than the existing methods.

The purpose of this paper is the derivation and efficient implementation of surface impedance boundary conditions (SIBCs) for nonlinear magnetic conductors.

An approach based on perturbation theory is proposed, which expands to nonlinear problems the methods already developed by the authors for linear problems. Differently from the linear case, for which the analytical solution of the diffusion equation in the semi-infinite space for the magnetic field is available, in the nonlinear case the corresponding nonlinear diffusion equation must be solved numerically. To this aim, a suitable smooth map is defined to reduce the semi-infinite computational domain to a finite one; then the diffusion equation is solved by a Galerkin method relying on basis functions constructed via the push-forward of a Lagrangian polynomial basis whose degrees of freedom are collocated at Gauss–Lobatto nodes. The use of such basis in connection with a suitable under-integration naturally leads to mass-lumping without impacting the order of the method. The solution of the diffusion equation is coupled with a boundary element method formulation for the case of parallel magnetic conductors in terms of E and B fields.

The results are validated by comparison with full nonlinear finite element method simulations showing very good accordance at a much lower computational cost.

Limitations of the method are those arising from perturbation theory: the introduced small parameter must be much less than one. This implies that the penetration depth of the magnetic field into the magnetic and conductive media must be much smaller than the characteristic size of the conductor.

The efficient implementation of a nonlinear SIBC based on a perturbation approach is proposed for an electric and magnetic field formulation of the two-dimensional problem of current driven parallel solid conductors.

This paper aims to solve the Falkner–Skan problem over an isothermal moving wedge using the combination of the quasilinearization method and the fractional order of rational Chebyshev function (FRC) collocation method on a semi-infinite domain.

The quasilinearization method converts the equation into a sequence of linear equations, and then by using the FRC collocation method, these linear equations are solved. The governing nonlinear partial differential equations are reduced to the nonlinear ordinary differential equation by similarity transformations.

The entropy generation and the effects of the various parameters of the problem are investigated, and various graphs for them are plotted.

Very good approximation solutions to the system of equations in the problem are obtained, and the convergence of numerical results is shown by using plots and tables.

This paper aims to propose a new numerical method for the solution of the Blasius and magnetohydrodynamic (MHD) Falkner-Skan boundary-layer equations. The Blasius and MHD Falkner-Skan equations are third-order nonlinear boundary value problems on the semi-infinite domain.

The approach is based upon modified rational Bernoulli functions. The operational matrices of derivative and product of modified rational Bernoulli functions are presented. These matrices together with the collocation method are then utilized to reduce the solution of the Blasius and MHD Falkner-Skan boundary-layer equations to the solution of a system of algebraic equations.

The method is computationally very attractive and gives very accurate results.

Many problems in science and engineering are set in unbounded domains. One approach to solve these problems is based on rational functions. In this work, a new rational function is used to find solutions of the Blasius and MHD Falkner-Skan boundary-layer equations.

The purpose of this paper is to present a new approximate analytical procedure to obtain dual solutions of nonlinear differential equations arising in mixed convection flow in a semi-infinite domain. This method, which is based on Padé-approximation and homotopy–Padé technique, is applied to a model of magnetohydrodynamic Falkner–Skan flow as well. These examples indicate that the method can be successfully applied to solve nonlinear differential equations arising in science and engineering.

Homotopy–Padé method.

The main focus of the paper is on the prediction of the multiplicity of the solutions, however we have calculated multiple (dual) solutions of the model problem namely, mixed convection heat transfer in a porous medium.

The authors conjecture here that the combination of traditional–Pade and Hankel–Pade generates a useful procedure to predict multiple solutions and to calculate prescribed parameter with acceptable accuracy as well. Validation of this conjecture for other further examples is a challenging research opportunity.

Dual solutions of nonlinear differential equations arising in mixed convection flow in a semi-infinite domain.

In this study, the authors are using two modified methods.

The purpose of this paper is to propose a Tau method for solving nonlinear Blasius equation which is a partial differential equation on a flat plate.

The operational matrices of derivative and product of modified generalized Laguerre functions are presented. These matrices together with the Tau method are then utilized to reduce the solution of the Blasius equation to the solution of a system of nonlinear equations.

The paper presents the comparison of this work with some well‐known results and shows that the present solution is highly accurate.

This paper demonstrates solving of the nonlinear Blasius equation with an efficient method.

The purpose of this paper is to propose a stable high-order absorbing boundary condition (ABC) based on new continued fraction for scalar wave propagation in 2D and 3D unbounded layers.

The ABC is obtained based on continued fraction (CF) expansion of the frequency-domain dynamic stiffness coefficient (DtN kernel) on the artificial boundary of a truncated infinite domain. The CF which has been used to the thin layer method in [69] will be applied to the DtN method to develop a time-domain high-order ABC for the transient scalar wave propagation in 2D. Furthermore, a new stable composite-CF is proposed in this study for 3D unbounded layers by nesting the above CF for 2D layer and another CF.

The ABS has been transformed from frequency to time domain by using the auxiliary variable technique. The high-order time-domain ABC can couple seamlessly with the finite element method. The instability of the ABC-FEM coupled system is discussed and cured.

This manuscript establishes a stable high-order time-domain ABC for the scalar wave equation in 2D and 3D unbounded layers, which is based on the new continued fraction. The high-order time-domain ABC can couple seamlessly with the finite element method. The instability of the coupled system is discussed and cured.

The response of the Pine Flat dam–water–foundation rock system is studied by a new described approach (i.e. FE-(FE-TE)-FE). The initial part of study is focused on the time harmonic analysis. In this part, it is possible to compare the transfer functions against corresponding responses obtained by the FE-(FE-HE)-FE approach (referred to as exact method which employs a rigorous fluid hyper-element). Subsequently, the transient analysis is carried out. In that part, it is only possible to compare the results for low and high normalized reservoir length cases. Therefore, the sensitivity of results is controlled due to normalized reservoir length values.

In the present study, dynamic analysis of a typical concrete gravity dam–water–foundation rock system is formulated by the FE-(FE-TE)-FE approach. In this technique, dam and foundation rock are discretized by plane solid finite elements while, water domain near-field region is discretized by plane fluid finite elements. Moreover, the H-W (i.e. Hagstrom–Warburton) high-order condition is imposed at the reservoir truncation boundary. This task is formulated by employing a truncation element at that boundary. It is emphasized that reservoir far-field is excluded from the discretized model.

High orders of H-W condition, such as O5-5 considered herein, generate highly accurate responses for both possible excitations under both types of full reflective and absorptive reservoir bottom conditions. It is such that transfer functions are hardly distinguishable from corresponding exact responses obtained through the FE-(FE-HE)-FE approach in time harmonic analyses. This is controlled for both low and high normalized reservoir length cases (L/H = 1 and 3). Moreover, it can be claimed that transient analysis leads practically to exact results (in numerical sense) when one is employing high order H-W truncation element. In other words, the results are not sensitive to reservoir normalized length under these circumstances.

Dynamic analysis of concrete gravity dam–water–foundation rock systems is formulated by a new method. The salient aspect of the technique is that it utilizes H-W high-order condition at the truncation boundary. The method is discussed for all types of excitation and reservoir bottom conditions.

The combination method, combined finite element‐boundary element approach, is suitable for unbounded field problems. Although this technique attains a high degree of accuracy, the matrix of the discretized system equation is not banded but sometimes densely or sparsely populated. We reported the development of an infinite boundary element for 2‐D Laplace problems, with which the bandwidth of the discretized system matrix does not increase beyond that of the finite element region. In this paper, we extend this approach and propose another infinite boundary element for 2‐D Helmholtz problems. To illustrate the validity of the proposed technique, some numerical examples are given and their results are compared with those of other methods.