Numerical solution of steady and unsteady flow over a vibrating airfoil in a channel

The work deals with a numerical solution of 2D steady and unsteady inviscid incompressible flow over the profile NACA 0012 in a channel. The finite volume method (FVM) in a form of cell-centered explicit schemes at quadrilateral C-mesh is used. Governing system of equations is the system of incompressible Euler equations. The method of artificial compressibility and time dependent method is applied to steady computations. The small disturbance theory (SDT) applied to a numerical solution of flow over a rotated profile by a small angle only is mentioned. Brief introduction is given to the Arbitrary (Semi) Lagrangian-Eulerian (ALE) method used for unsteady computations. Some numerical results of unsteady flow over a vibrating profile achieved by both SDT and ALE method are presented. Unsteady flow is caused by prescribed oscillations of the profile (one degree of freedom) fixed in an elastic axis. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On a finite element method for a compressible viscous flow in a MOCVD reactor

We propose a stable finite element method for approximating the flow of a chemically reacting gas mixture in an MOCVD (metal‐organic chemical vapor deposition) reactor. The flow is governed by the full compressible Navier‐Stokes equations extended by transport equations for the chemical species, the energy equations and the equation of state, together with boundary conditions providing information on the reactor geometry and experimental conditions. The equations form a semilinear system with a constraint for which the corresponding pressure term is not the Lagrangian multiplier. An application of our method to a real world model of growth of GaAs shows the consistency with experimental data. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

Stagnation slip flow and heat transfer over a nonlinear stretching sheet

Stagnation slip flow and heat transfer characteristics of a viscous fluid over a nonlinear stretching surface has been investigated. The governing partial differential equations are transformed to nonlinear ordinary differential equations using similarity transformations. The analytical solution of the nonlinear system is obtained in series form using the very efficient homotopy analysis method (HAM). Convergence of the series solution is shown explicitly. Important features of flow and heat transfer characteristics are plotted and discussed. Comparison is made with existing numerical results when the stagnation‐point and slip effects are excluded. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

Lie group transformations for self‐similar shocks in a gas with dust particles

In this paper, Lie group of transformation method is used to investigate the self‐similar solutions for the system of partial differential equations describing one‐dimensional unsteady plane flow of an inviscid gas with large number of small dust particles. The forms of the drag force D and the heat transfer rate Q experienced by the particle not in equilibrium with the gas have been derived for which the system of equations admits self‐similar solutions. A particular solution to the problem in one case have been found out and is used to study the effect of the dust particles on the similarity exponent. Copyright © 2009 John Wiley & Sons, Ltd.     

Numerical solution of fractional differential equations with a Tau method based on Legendre and Bernstein polynomials

In this paper, we state and prove a new formula expressing explicitly the integratives of Bernstein polynomials (or B‐polynomials) of any degree and for any fractional‐order in terms of B‐polynomials themselves. We derive the transformation matrices that map the Bernstein and Legendre forms of a degree‐n polynomial on [0,1] into each other. By using their transformation matrices, we derive the operational matrices of integration and product of the Bernstein polynomials. These matrices together with the Tau method are then utilized to reduce the solution of this problem to the solution of a system of algebraic equations. The method is applied to solve linear and nonlinear fractional differential equations. Copyright © 2013 John Wiley & Sons, Ltd.     

Existence results for a degenerated nonlinear elliptic partial differential equation

The aim of this paper is to establish the existence of weak solutions to a steady state two-dimensional irrotational compressible flow around a thin profile. This flow is described by the small disturbance equations. If the speed of sound exceeds the fluid one, the governing equations remain elliptic. But when the fluid speed is beyond the sound one, the flow becomes locally hyperbolic and shock waves arise. For a modified elliptic model, using convexity arguments, we prove the existence of a solution which is the solution to the first model when the flow remains subsonic.     

The fully Sinc‐Galerkin method for time‐dependent boundary conditions

The fully Sinc‐Galerkin method is developed for a family of complex‐valued partial differential equations with time‐dependent boundary conditions. The Sinc‐Galerkin discrete system is formulated and represented by a Kronecker product form of those equations. The numerical solution is efficiently calculated and the method exhibits an exponential convergence rate. Several examples, some with a real‐valued solution and some with a complex‐valued solution, are used to demonstrate the performance of this method. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2004     

On global solutions of gravity driven nonhomogeneous viscous flows in a bounded domain

We consider an initial boundary value problem for nonhomogeneous Navier‐Stokes equations with a uniform gravitational field. For any given steady density profile whose derivatives are sufficiently close to a negative constant, we show that there exists a unique global solution if the initial perturbation with respect to the steady state is sufficiently small.     

Convergence of the solution for a domain decomposed,heterogeneously modeled flow past a concave profile problem

The free streamline problem investigated is that of fluid flow past a symmetric truncated concave‐shaped profile between walls. An open wake or cavity is formed behind the profile. Conformal mapping techniques are used to solve this problem. The problem formulated in the hodograph plane is decomposed into two nonoverlapping domains. Heterogeneous modeling is then used to describe the problems, i.e., a different governing differential equation in each domain. In one of these domains, a Baiocchi‐type transformation is used to obtain a fixed domain formulation for the part of the transformed problem containing an unknown boundary. In the other domain, the Baiocchi‐type transformation is extended across the boundary between the two domains, thus yielding a different problem formulation. This also assures that the dependent variables and their normal derivatives are continuous along this common boundary. The numerical solution scheme, a successive over‐relaxation approach, is applied over the whole problem domain with the use of a projection‐operation over only the fixed domain formulated part. Numerical results are obtained for the case of a truncated circular profile. These results are found to be in good agreement with another published result. The existence and uniqueness of the solution to the problem as a variational inequality is shown, and the convergence of the numerical solution using a domain decomposition method scheme is demonstrated by assuming some convergence property on the common interface of the two subdomains. © 2000 John Wiley & Sons, Inc. Numeer Methods Partial Differential Eq 16: 459–479, 2000     

Asymptotic stability of stationary solutions to the compressible bipolar Navier–Stokes–Poisson equations

In this paper, we consider the compressible bipolar Navier–Stokes–Poisson equations with a non‐flat doping profile in three‐dimensional space. The existence and uniqueness of the non‐constant stationary solutions are established when the doping profile is a small perturbation of a positive constant state. Then under the smallness assumption of the initial perturbation, we show the global existence of smooth solutions to the Cauchy problem near the stationary state. Finally, the convergence rates are obtained by combining the energy estimates for the nonlinear system and the L²‐decay estimates for the linearized equations. Copyright © 2017 John Wiley & Sons, Ltd.     

The tanh method: a tool for solving certain classes of non‐linear PDEs

The tanh (or hyperbolic tangent) method is a powerful technique to look for travelling waves when dealing with one‐dimensional non‐linear wave and evolution equations. In particular, this method is well suited for those problems where dispersion, convection and reaction–diffusion play an important role. To show the strength of this method we study a coupled set (the so‐called Boussinesq equations) which arises in the theory of non‐linear dispersive water waves. As a result, a solitary wave profile is found which generalizes an earlier result, the famous Korteweg‐de Vries solitary wave solution. Copyright © 2005 John Wiley & Sons, Ltd.     

The Cauchy problem of a fluid‐particle interaction model with external forces

In this paper, we consider the Cauchy problem of a fluid‐particle interaction model with external forces. We first construct the asymptotic profile of the system. The global existence and uniqueness theorem for the solution near the profile is given. Finally, optimal decay rate of the solution to the background profile is obtained by combining the decay rate analysis of a linearized equation with energy estimates for the nonlinear terms. The main method used in this paper is the energy method combining with the macro‐micro decomposition.     

Approximate solution of a complete singular integral equation of the first kind with a fixed hypersingularity by the overlapping method

We suggest a new method for the numerical solution of a singular integral equation of the first kind with a fixed hypersingularity, which arises in the problem on the flow past a profile with an ejector of the external flow. This method permits one to obtain a solution of the characteristic and complete integral equations with an interpolation degree of accuracy.     

Indirect boundary integral formulation for the solution of shear‐thinning flow inside single rotor mixers

An efficient indirect boundary integral formulation for the evaluation of inelastic non‐Newtonian shear‐thinning flows at low Reynolds number is presented in this article. The formulation is based on the solution of a homogeneous Stokes flow field and the use of a particular solution for the nonlinear non‐Newtonian terms that yields the complete solution to the problem. Matrix multiplications are reduced in comparison to other means of handling nonlinear terms in boundary integral formulations such as the dual reciprocity method. The iterative solution of the nonlinear system of equations has been performed with a modified Newton‐Raphson method obtaining accurate results for values of the power law index as low as 0.4 without domain partitioning. Geometries such as Couette flow and a typical industrial polymer mixer have been analyzed with the proposed method obtaining good results with a reduction in computational cost compared with other equivalent formulations. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27:1610–1627, 2011     

Analysis of a multistate control problem related to food technology

This paper is concerned with an optimal control problem related to the determination of an optimal profile for the steam temperature into the autoclave along the processing of canned foods. The problem studies a system coupling the evolution Navier-Stokes equations with the heat transfer equation by natural convection (the so-called Boussinesq equations), and with the microorganisms removal equation. The essential difficulties in the study of this multistate control problem arise from the lack of uniqueness for the solution of the state system. Here we obtain—after a careful analysis of the problem mathematical formulation—the uniqueness of part of the state, and the existence of optimal solutions.     

Towards a consistency proof for immiscible lattice BGK

Emanating from the incompressible Navier‐Stokes equations, a widely used mathematical model of immiscible twophase flow is rigorously derived. Once the model is fixed, it is analysed to what extend the immiscible lattice BGK method approximates the solution of the model equations. Using equivalent moment analysis and some elementary tools from differential geometry we are able to show weak consistency if certain well founded but unproven presumptions hold.     

Travelling waves of a predator–prey model with a spatiotemporal delay

This paper is concerned with the existence of traveling waves to a predator–prey model with a spatiotemporal delay. By analyzing the corresponding characteristic equations, the local stability of a positive steady state and each of boundary steady states are established, and the existence of Hopf bifurcation at the positive steady state is also discussed. By constructing a pair of upper–lower solutions and by using the cross‐iteration method as well as the Schauder's fixed‐point theorem, the existence of a traveling wave solution connecting the semi‐trivial steady state and the positive steady state is proved. Numerical simulations are carried out to illustrate the main theoretical results. Copyright © 2013 John Wiley & Sons, Ltd.     

Bivariate spline method for numerical solution of steady state Navier–Stokes equations over polygons in stream function formulation

We use the bivariate spline finite elements to numerically solve the steady state Navier–Stokes equations. The bivariate spline finite element space we use in this article is the space of splines of smoothness r and degree 3r over triangulated quadrangulations. The stream function formulation for the steady state Navier–Stokes equations is employed. Galerkin's method is applied to the resulting nonlinear fourth‐order equation, and Newton's iterative method is then used to solve the resulting nonlinear system. We show the existence and uniqueness of the weak solution in H²(Ω) of the nonlinear fourth‐order problem and give an estimate of how fast the numerical solution converges to the weak solution. The Galerkin method with C¹ cubic splines is implemented in MATLAB. Our numerical experiments show that the method is effective and efficient. © 2000 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 16: 147–183, 2000     

Using empirical Eigenfunctions and Galerkin method to two‐phase transport models

In this article, we consider a nonlinear partial differential system describing two‐phase transports and try to recover the source term and the nonlinear diffusion term when the state variable is known at different profile times. To this end, we use a POD‐Galerkin procedure in which the proper orthogonal decomposition technique is applied to the ensemble of solutions to derive empirical eigenfunctions. These empirical eigenfunctions are then used as basis functions within a Galerkin method to transform the partial differential equation into a set of ordinary differential equations. Finally, the validation of the used method has been evaluated by some numerical examples. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 23: 456–474, 2007     

Numerical solution of the free‐surface viscous flow on a horizontal rotating elliptical cylinder

The numerical solution of the free‐surface fluid flow on a rotating elliptical cylinder is presented. Up to the present, research has concentrated on the circular cylinder for which steady solutions are the main interest. However, for noncircular cylinders, such as the ellipse, steady solutions are no longer possible, but there will be periodic solutions in which the solution is repeated after one full revolution of the cylinder. It is this new aspect that makes the investigation of noncircular cylinders novel. Here we consider both the time‐dependent and periodic solutions for zero Reynolds number fluid flow. The numerical solution is expedited by first mapping the fluid film domain onto a rectangle such that the position of the free‐surface is determined as part of the solution. For the time‐dependent case a simple time‐marching method of lines approach is adopted. For the periodic solution the discretised nonlinear equations have to be solved simultaneously over a time period. The resulting large system of equations is solved using Newton's method in which the form of the Jacobian enables a straightforward decomposition to be implemented, which makes matrix inversion manageable. In the periodic case all derivatives have been approximated pseudospectrally with the time derivative approximated by a differentiation matrix which has been specially derived so that the weight of fluid is algebraically conserved. Of interest is the solution for which the weight of fluid is at its maximum possible value, and this has been obtained by increasing the weight until a consistency break‐down occurs. Time‐dependent solutions do not produce the periodic solution after a long time‐scale but have protuberances which are constantly appearing and disappearing. Periodic solutions exhibit spectral accuracy solutions and maximum supportable weight solutions have been obtained for ranges of eccentricity and angular velocity. The maximum weights are less than and approximately proportional to those obtained for the circular case. The shapes of maximum weight solutions is distinctly different from sub‐maximum weight solutions. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008