Asymptotic and numerical solutions of three‐dimensional boundary‐layer flow past a moving wedge

We consider a laminar boundary‐layer flow of a viscous and incompressible fluid past a moving wedge in which the wedge is moving either in the direction of the mainstream flow or opposite to it. The mainstream flows outside the boundary layer are approximated by a power of the distance from the leading boundary layer. The variable pressure gradient is imposed on the boundary layer so that the system admits similarity solutions. The model is described using 3‐dimensional boundary‐layer equations that contains 2 physical parameters: pressure gradient (β) and shear‐to‐strain‐rate ratio parameter (α). Two methods are used: a linear asymptotic analysis in the neighborhood of the edge of the boundary layer and the Keller‐box numerical method for the full nonlinear system. The results show that the flow field is divided into near‐field region (mainly dominated by viscous forces) and far‐field region (mainstream flows); the velocity profiles form through an interaction between 2 regions. Also, all simulations show that the subsequent dynamics involving overshoot and undershoot of the solutions for varying parameter characterizing 3‐dimensional flows. The pressure gradient (favorable) has a tendency of decreasing the boundary‐layer thickness in which the velocity profiles are benign. The wall shear stresses increase unboundedly for increasing α when the wedge is moving in the x‐direction, while the case is different when it is moving in the y‐direction. Further, both analysis show that 3‐dimensional boundary‐layer solutions exist in the range −1<α<∞. These are some interesting results linked to an important class of boundary‐layer flows.     

The Application of Neumann Type Systems for Solving Integrable Nonlinear Evolution Equations

An algorithm to obtain finite‐gap solutions of integrable nonlinear evolution equations (INLEEs) is provided by using the Neumann type systems in the framework of algebraic geometry. From the nonlinearization of Lax pairs, some INLEEs in 1+1 and 2+1 dimensions are reduced into a class of new Neumann type systems separating the spatial and temporal variables of INLEEs over a symplectic submanifold (M, ω²) . Based on the Lax representations of INLEEs, we deduce the Lax–Moser matrix for those Neumann type systems that yield the integrals of motion, elliptic variables, and a hyperelliptic curve of Riemann surface. Then, we attain the Liouville integrability for a hierarchy of Neumann type systems in view of a Lax equation on (M, ω²) and a set of quasi‐Abel–Jacobi variables. We also specify the relationship between Neumann type systems and INLEEs, where the involutive solutions of Neumann type systems give rise to the finite parametric solutions of INLEEs and the Neumann map cuts out a finite dimensional invariant subspace for INLEEs. Under the Abel–Jacobi variables, the Neumann type flows, the 1+1, and 2+1 dimensional flows are integrated with Abel–Jacobi solutions; as a result, the finite‐gap solutions expressed by Riemann theta functions for some 1+1 and 2+1 dimensional INLEEs are achieved through the Jacobi inversion with the aid of the Riemann theorem.     

Wrinkles as the Result of Compressive Stresses in an Annular Thin Film

It is well‐known that an elastic sheet loaded in tension will wrinkle and that the length scale of the wrinkles tends to 0 with vanishing thickness of the sheet [Cerda and Mahadevan, Phys. Rev. Lett. 90, 074302 (2003)]. We give the first mathematically rigorous analysis of such a problem. Since our methods require an explicit understanding of the underlying (convex) relaxed problem, we focus on the wrinkling of an annular sheet loaded in the radial direction [Davidovitch et al., Proc. Natl. Acad. Sci. 108 (2011), 18227]. Our main achievement is identification of the scaling law of the minimum energy as the thickness of the sheet tends to 0. This requires proving an upper bound and a lower bound that scale the same way. We prove both bounds first in a simplified Kirchhoff‐Love setting and then in the nonlinear three‐dimensional setting. To obtain the optimal upper bound, we need to adjust a naive construction (one family of wrinkles superimposed on a planar deformation) by introducing a cascade of wrinkles. The lower bound is more subtle, since it must be ansatz‐free. © 2014 Wiley Periodicals, Inc.     

Nonlinear Hamiltonian waves with constant frequency and surface waves on vorticity discontinuities

Waves with constant, nonzero linearized frequency form an interesting class of nondispersive waves whose properties differ from those of nondispersive hyperbolic waves. We propose an inviscid Burgers‐Hilbert equation as a model equation for such waves and give a dimensional argument to show that it models Hamiltonian surface waves with constant frequency. Using the method of multiple scales, we derive a cubically nonlinear, quasi‐linear, nonlocal asymptotic equation for weakly nonlinear solutions. We show that the same asymptotic equation describes surface waves on a planar discontinuity in vorticity in two‐dimensional inviscid, incompressible fluid flows. Thus, the Burgers‐Hilbert equation provides an effective equation for these waves. We describe the Hamiltonian structure of the Burgers‐Hilbert and asymptotic equations, and show that the asymptotic equation can also be derived by means of a near‐identity transformation. We derive a semiclassical approximation of the asymptotic equation and show that spatially periodic, harmonic traveling waves are linearly and modulationally stable. Numerical solutions of the Burgers‐Hilbert and asymptotic equations are in excellent agreement in the appropriate regime. In particular, the lifespan of small‐amplitude smooth solutions of the Burgers‐Hilbert equation is given by the cubically nonlinear timescale predicted by the asymptotic equation. © 2009 Wiley Periodicals, Inc.     

On the convergence of operator splitting for the Rosenau–Burgers equation

We present convergence analysis of operator splitting methods applied to the nonlinear Rosenau–Burgers equation. The equation is first splitted into an unbounded linear part and a bounded nonlinear part and then operator splitting methods of Lie‐Trotter and Strang type are applied to the equation. The local error bounds are obtained by using an approach based on the differential theory of operators in Banach space and error terms of one and two‐dimensional numerical quadratures via Lie commutator bounds. The global error estimates are obtained via a Lady Windermere's fan argument. Lastly, a numerical example is studied to confirm the expected convergence order.     

Numerical approximation of Quasi‐Newtonian flows by ALE‐FEM

The aim of this work is to investigate the numerical approximation of a nonlinear, time‐dependent quasi‐Newtonian flow problem formulated in the framework of Arbitrary Lagrangian Eulerian method. We present some stability results and convergence analysis of finite element solutions for semidiscrete and fully discretized problems, respectively. Numerical results supporting the derived error estimate are also presented. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2012     

A fully variational algorithmic formulation for wrinkling at finite strains

This contribution is concerned with an efficient novel algorithmic formulation for wrinkling at finite strains. In contrast to previously published numerical implementations, the advocated method is fully variational. More precisely, the parameters describing wrinkles or slacks, together with the unknown deformation mapping, are computed jointly by minimizing the potential energy of the considered mechanical system. Furthermore, the wrinkling criteria are naturally included within the presented variational framework. The presented approach allows to employ three–dimensional constitutive models directly, i.e., plane stress conditions characterizing membranes are variationally enforced by minimizing the potential energy with respect to the transversal strains. Since the proposed formulation for wrinkling in membranes is fully variational, it can be conveniently combined with other variational methods (based on energy minimization). As an example, a variationally consistent framework for finite strain plasticity theory is considered. More precisely, the minimization principle characterizing wrinkling in elastic membranes and that describing plasticity in inelastic solids are coupled leading to a novel variational approach for inelastic membranes. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

圆柱贮箱液体非线性晃动的多维模态分析方法

将Faltinsen等提出的多维模态理论应用到求解圆柱贮箱液体非线性晃动问题中．根据Narimanov-Moiseiev的三阶渐近假设关系,通过选取主导模态以及确定它们的阶次关系,将一般形式的无穷维模态系统降为五维渐近模态系统,即描述自由液面波高的广义坐标之间相互耦合的二阶非线性常微分方程组．通过对这个模态系统的数值积分,得到了与以前的理论分析和实验结果相吻合的非线性现象．研究结果表明,多维模态方法是用来求解液体非线性晃动动力学的一个很好的工具．在我们的下一步工作中,将继续发展这种方法,用来研究更为复杂的晃动问题．     

Web‐spline‐based finite element approximation of some quasi‐newtonian flows: Existence‐uniqueness and error bound

This article deals with the web‐spline‐based finite element approximation of quasi‐Newtonian flows. First, we consider the scalar elliptic p‐Laplace problem. Then, we consider quasi‐Newtonian flows where viscosity obeys power law or Carreau law. We prove well‐posedness at the continuous as well as the discrete level. We give some error bounds for the solution of quasi‐Newtonian flow problem based on the web‐spline method. Finally, we provide the numerical results for the p‐Laplace problem. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq31: 54–77, 2015     

The dual reciprocity boundary elements method for the linear and nonlinear two‐dimensional time‐fractional partial differential equations

In this paper, we apply the dual reciprocity boundary elements method for the numerical solution of two‐dimensional linear and nonlinear time‐fractional modified anomalous subdiffusion equations and time‐fractional convection–diffusion equation. The fractional derivative of problems is described in the Riemann–Liouville and Caputo senses. We employ the linear radial basis function for interpolation of the nonlinear, inhomogeneous and time derivative terms. This method is improved by using a predictor–corrector scheme to overcome the nonlinearity which appears in the nonlinear problems under consideration. The accuracy and efficiency of the proposed schemes are checked by five test problems. The proposed method is employed for solving some examples in two dimensions on unit square and also in complex regions to demonstrate the efficiency of the new technique. Copyright © 2016 John Wiley & Sons, Ltd.     

A preliminary investigation on an ELLAM scheme for linear transport equations

We present an Eulerian‐Lagrangian localized adjoint method (ELLAM) for linear advection‐reaction partial differential equations in multiple space dimensions. We carry out numerical experiments to investigate the performance of the ELLAM scheme with a range of well‐perceived and widely used methods in fluid dynamics including the monotonic upstream‐centered scheme for conservation laws (MUSCL), the minmod method, the flux‐corrected transport method (FCT), and the essentially non‐oscillatory (ENO) schemes and weighted essentially non‐oscillatory (WENO) schemes. These experiments show that the ELLAM scheme is very competitive with these methods in the context of linear transport PDEs, and suggest/justify the development of ELLAM‐based simulators for subsurface porous medium flows and other applications. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 22–43, 2003     

An efficient method for fractional nonlinear differential equations by quasi‐Newton's method and simplified reproducing kernel method

An efficient method for nonlinear fractional differential equations is proposed in this paper. This method consists of 2 steps. First, we linearize the nonlinear operator equation by quasi‐Newton's method, which is based on Fréchet derivative. Then we solve the linear fractional differential equations by the simplified reproducing kernel method. The convergence of the quasi‐Newton's method is discussed for the general nonlinear case as well. Finally, some numerical examples are presented to illustrate accuracy, efficiency, and simplicity of the method.     

A New Mechanism For Linear and Nonlinear Hydrodynamic Instability

A certain class of three dimensional disturbances is studied which exhibit a strong resonance and consequently lead to the rapid development of relatively large amplitudes. For parallel flows a nonlinear theory is developed and results in evolution equations quite different from those usually found. It is believed that the selection mechanism explored here is relevant to a variety of nonlinear wave and instability problems.     

MIXED MONOTONE ITERATIVE TECHNIQUES FOR SEMILINEAR EVOLUTION EQUATIONS IN BANACH SPACES

This paper is concerned with initial value problems for semilinear evolution equations in Banach spaces. The abstract iterative schemes are constructed by combining the theory of semigroups of linear operators and the method of mixed monotone iterations. Some existence results on minimal and maximal (quasi)solutions are established for abstract semilinear evolution equations with mixed monotone or mixed quasimonotone nonlinear terms. To illustrate the main results, applications to ordinary differential equations and partial differential equations are also given.     

Solving hyperbolic conservation laws using multiquadric quasi‐interpolation

In this article, we apply the univariate multiquadric (MQ) quasi‐interpolation to solve the hyperbolic conservation laws. At first we construct the MQ quasi‐interpolation corresponding to periodic and inflow‐outflow boundary conditions respectively. Next we obtain the numerical schemes to solve the partial differential equations, by using the derivative of the quasi‐interpolation to approximate the spatial derivative of the differential equation and a low‐order explicit difference to approximate the temporal derivative of the differential equation. Then we verify our scheme for the one‐dimensional Burgers' equation (without viscosity). We can see that the numerical results are very close to the exact solution and the computational accuracy of the scheme is ??(τ), where τ is the temporal step. We can improve the accuracy by using the high‐order quasi‐interpolation. Moreover the methods can be generalized to the other equations. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

A regular perturbation analytical‐numerical method for the evolution of precancerous lesions caused by the human papillomavirus

A numerical method is developed to analyze the behavior of the evolution of the lesions at the cervical cells caused by the human papillomavirus. The model to be solved consists in a one‐dimensional nonlinear advection–diffusion‐reaction equation. Such equation is approximated by a consistent explicit difference scheme which is based on regular perturbation theory. A constructive procedure for the numerical scheme is given and finally an illustrative example of the evolution of a mild dysplasia is included. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 847–855, 2015     

A fully discrete numerical scheme for weighted mean curvature flow

Summary. We analyze a fully discrete numerical scheme approximating the evolution of n–dimensional graphs under anisotropic mean curvature. The highly nonlinear problem is discretized by piecewise linear finite elements in space and semi–implicitly in time. The scheme is unconditionally stable und we obtain optimal error estimates in natural norms. We also present numerical examples which confirm our theoretical results. Received October 2, 2000 / Published online July 25, 2001     

Integrable Discrete Model for One‐Dimensional Soil Water Infiltration

We propose an integrable discrete model of one‐dimensional soil water infiltration. This model is based on the continuum model by Broadbridge and White, which takes the form of nonlinear convection–diffusion equation with a nonlinear flux boundary condition at the surface. It is transformed to the Burgers equation with a time‐dependent flux term by the hodograph transformation. We construct a discrete model preserving the underlying integrability, which is formulated as the self‐adaptive moving mesh scheme. The discretization is based on linearizability of the Burgers equation to the linear diffusion equation, but the naïve discretization based on the Euler scheme which is often used in the theory of discrete integrable systems does not necessarily give a good numerical scheme. Taking desirable properties of a numerical scheme into account, we propose an alternative discrete model that produces solutions with similar accuracy to direct computation on the original nonlinear equation, but with clear benefits regarding computational cost.     

The method of fundamental solutions for the Oseen steady‐state viscous flow past obstacles of known or unknown shapes

In this paper, the steady‐state Oseen viscous flow equations past a known or unknown obstacle are solved numerically using the method of fundamental solutions (MFS), which is free of meshes, singularities, and numerical integrations. The direct problem is linear and well‐posed, whereas the inverse problem is nonlinear and ill‐posed. For the direct problem, the MFS computations of the fluid flow characteristics (velocity, pressure, drag, and lift coefficients) are in very good agreement with the previously published results obtained using other methods for the Oseen flow past circular and elliptic cylinders, as well as past two circular cylinders. In the inverse obstacle problem the boundary data and the internal measurement of the fluid velocity are minimized using the MATLAB^© optimization toolbox lsqnonlin routine. Regularization was found necessary in the case the measured data are contaminated with noise. Numerical results show accurate and stable reconstructions of various star‐shaped obstacles of circular, bean, or peanut cross‐section.