Nonlinear Boundary Value Problems on Semi-Infinite Intervals using Weighted Spaces: An Upper and Lower Solution Approach

This paper presents a lower and upper solution approach for singular second order boundary value problems on the half line and establishes the existence of positive, unbounded and monotone solutions. The project is supported by the fund of National Natural Science(10571111) and the fund of Natural Science of Shandong Province(Y2005A07).     

Algorithms for vector field generation in mass consistent models

Diagnostic models in meteorology are based on the fulfillment of some time independent physical constraints as, for instance, mass conservation. A successful method to generate an adjusted wind field, based on mass conservation equation, was proposed by Sasaki and leads to the solution of an elliptic problem for the multiplier. Here we study the problem of generating an adjusted wind field from given horizontal initial velocity data, by two ways. The first one is based on orthogonal projection in Hilbert spaces and leads to the same elliptic problem but with natural boundary conditions for the multiplier. We derive from this approach the so called E–algorithm. An innovative alternative proposal is obtained from a second approach where we consider the saddle–point formulation of the problem—avoiding boundary conditions for the multiplier— and solving this problem by iterative conjugate gradient methods. This leads to an algorithm that we call the CG–algorithm, which is inspired from Glowinsk's approach to solve Stokes–like problems in computational fluid dynamics. Finally, the introduction of new boundary conditions for the multiplier in the elliptic problem generates better adjusted fields than those obtained with the original boundary conditions. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

A finite difference method for singularly perturbed differential-difference equations with layer and oscillatory behavior

In this paper, we present a finite difference method for singularly perturbed linear second order differential-difference equations of convection–diffusion type with a small shift, i.e., where the second order derivative is multiplied by a small parameter and the shift depends on the small parameter. Similar boundary value problems are associated with expected first-exit times of the membrane potential in models of neurons. Here, the study focuses on the effect of shift on the boundary layer behavior or oscillatory behavior of the solution via finite difference approach. An extensive amount of computational work has been carried out to demonstrate the proposed method and to show the effect of shift parameter on the boundary layer behavior and oscillatory behavior of the solution of the problem.     

Unsteady MHD flow of a non-Newtonian fluid on a porous plate

This paper deals with the study of the MHD flow of non-Newtonian fluid on a porous plate. Two exact solutions for non-torsionally generated unsteady hydromagnetic flow of an electrically conducting second order incompressible fluid bounded by an infinite non-conducting porous plate subjected to a uniform suction or blowing have been analyzed. The governing partial differential equation for the flow has been established. The mathematical analysis is presented for the hydromagnetic boundary layer flow neglecting the induced magnetic field. The effect of presence of the material constants of the second order fluid on the velocity field is discussed.     

广义二阶流体管内轴向流动

在流体的本构关系中引入分数阶导数运算，对于介于粘性与弹性之间的流体的描述更具有合理性。本文将这种关系引入二阶流体，研究其管内轴向流动。我们先求出了１／２阶导数的解析解，用以验证Ｌａｐｌａｃｅ数值反演的ＣＲＵＭＰ方法的有效性。然后用ＣＲＵＭＰ法分析二阶流体管内轴向流动的特征。分析表明粘弹性特征越明显的流体，其速度与应力对分数导数的阶数越具有敏感性。     

Rare events in the Boussinesq system with fluctuating dynamical boundary conditions

The Boussinesq system models various phenomena in geophysical and climate dynamics. It is a coupled system of the Navier-Stokes equations and the salinity transport equation. Due to uncertainty in salinity flux on fluid boundary, this system is subject to random fluctuations on the boundary. This stochastic Boussinesq system can be transformed into a random dynamical system. Rare events, or small probability events, are investigated in the context of large deviations. A large deviations principle is established via a weak convergence approach based on a recently developed variational representation of functionals of infinite dimensional Brownian motion.     

Exact solutions for some oscillating flows of a second grade fluid with a fractional derivative model

In the present work, the exact analytic solutions for some oscillating flows of a generalized second grade fluid are investigated using Fourier sine and Laplace transforms. A more appropriate model is presented for fluid material between viscous and elastic to introduce the fractional calculus approach into the constitutive relationship. This paper employs the fractional calculus approach to study second grade fluid flows. In order to avoid lengthy calculations of residues and contour integrals, the discrete inverse Laplace transform method has been used. Similar solutions for second grade fluid appear as the limiting cases of our solutions. The influence of pertinent parameters on the flows is delineated and appropriate conclusions are drawn.     

Joint Distributions for Interacting Fluid Queues

Motivated by recent traffic control models in ATM systems, we analyse three closely related systems of fluid queues, each consisting of two consecutive reservoirs, in which the first reservoir is fed by a two-state (on and off) Markov source. The first system is an ordinary two-node fluid tandem queue. Hence the output of the first reservoir forms the input to the second one. The second system is dual to the first one, in the sense that the second reservoir accumulates fluid when the first reservoir is empty, and releases fluid otherwise. In these models both reservoirs have infinite capacities. The third model is similar to the second one, however the second reservoir is now finite. Furthermore, a feedback mechanism is active, such that the rates at which the first reservoir fills or depletes depend on the state (empty or nonempty) of the second reservoir.The models are analysed by means of Markov processes and regenerative processes in combination with truncation, level crossing and other techniques. The extensive calculations were facilitated by the use of computer algebra. This approach leads to closed-form solutions to the steady-state joint distribution of the content of the two reservoirs in each of the models.     

Interacting Wave Fronts and Rarefaction Waves in a Second Order Model of Nonlinear Thermoviscous Fluids

A wave equation including nonlinear terms up to the second order for a thermoviscous Newtonian fluid is proposed. In the lossless case this equation results from an expansion to third order of the Lagrangian for the fundamental non-dissipative fluid dynamical equations. Thus it preserves the Hamiltonian structure, in contrast to the Kuznetsov equation, a model often used in nonlinear acoustics. An exact traveling wave front solution is derived from a generalized traveling wave assumption for the velocity potential. Numerical studies of the evolution of a number of arbitrary initial conditions as well as head-on colliding and confluent wave fronts exhibit several nonlinear interaction phenomena. These include wave fronts of changed velocity and amplitude along with the emergence of rarefaction waves. An analysis using the continuity of the solutions as well as the boundary conditions is proposed. The dynamics of the rarefaction wave is approximated by a collective coordinate approach in the energy balance equation.     

A cellular-automata method for studying porous media properties

A cellular-automata (CA) approach for investigating properties of porous media with tortuous channels and different smoothness of pore walls is proposed. This approach is aimed at combining two different CA models: the first one is intended for constructing the morphology of a porous material; the second, for simulating a fluid flow through it. The porous media morphology is obtained as a result of evolution of a cellular automaton, forming a “steady pattern.” The result is then used for simulating a fluid flow through a porous medium by applying the Lattice Gas CA model. The method has been tested on a small fragment of a porous material and implemented for investigating a carbon electrode of a hydrogen fuel cell on a multiprocessor cluster.     

On the solvability of a nonlinear pseudoparabolic problem

This paper is concerned with the study of an initial boundary value problem for a nonlinear second order pseudoparabolic equation arising in the unidirectional flow of a thermodynamic compatible third grade fluid. We establish some a priori bounds for the solution and prove its existence.     

Fluid-structure interaction of an elastic pipe immersed in a coaxial rigid pipe: a model for the Tullio Phenomenon

An axisymmetric, elastic pipe is filled with an incompressible fluid and is immersed in a second, coaxial rigid pipe which contains the same fluid. A pressure pulse in the outer fluid annulus deforms the elastic pipe which invokes a fluid motion in the fluid core. It is the aim of this study to investigate streaming phenomena in the core which may originate from such a fluid-structure interaction. This work presents a numerical solver for such a configuration. It was developed in the OpenFOAM software environment and is based on the Arbitrary Lagrangian Eulerian (ALE) approach for moving meshes. The solver features a monolithic integration of the one-dimensional, coupled system between the elastic structure and the outer fluid annulus into a dynamic boundary condition for the moving surface of the fluid core. Results indicate that our configuration may serve as a mechanical model of the Tullio Phenomenon (sound-induced vertigo). (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Modified equation based mesh adaptation algorithm for evolutionary scalar partial differential equations

It is well known that on uniform mesh classical higher order schemes for evolutionary problems yield an oscillatory approximation of the solution containing discontinuity or boundary layers. In this article, an entirely new approach for constructing locally adaptive mesh is given to compute nonoscillatory solution by representative “second” order schemes. This is done using modified equation analysis and a notion of data dependent stability of schemes to identify the solution regions for local mesh adaptation. The proposed algorithm is applied on scalar problems to compute the solution with discontinuity or boundary layer. Presented numerical results show underlying second order schemes approximate discontinuities and boundary layers without spurious oscillations.     

On non-Newtonian p-fluids. The pseudo-plastic case

In the following we study a class of stationary Navier-Stokes equations with shear dependent viscosity, under the non-slip (Dirichlet) boundary condition. We consider pseudo-plastic fluids. A fluid is said pseudo-plastic, or shear thinning, if in Eq. (1.1) below one has p<2. We are interested in global (i.e., up to the boundary) regularity results, in dimension n=3, for the second order derivatives of the velocity and the first order derivatives of the pressure. We consider a cubic domain Ω and impose the non-slip boundary condition only on two opposite faces. On the other faces we assume periodicity, as a device to avoid effective boundary conditions. This choice is made so that we work in a bounded domain Ω and simultaneously with a flat boundary.     

A second‐order linearized difference scheme for a strongly coupled reaction‐diffusion system

This article deals with the numerical solution to some models described by the system of strongly coupled reaction–diffusion equations with the Neumann boundary value conditions. A linearized three‐level scheme is derived by the method of reduction of order. The uniquely solvability and second‐order convergence in L₂‐norm are proved by the energy method. A numerical example is presented to demonstrate the accuracy and efficiency of the proposed method. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

A new stabilizing technique for boundary integral methods for water waves

Boundary integral methods to compute interfacial flows are very sensitive to numerical instabilities. A previous stability analysis by Beale, Hou and Lowengrub reveals that a very delicate balance among terms with singular integrals and derivatives must be preserved at the discrete level in order to maintain numerical stability. Such balance can be preserved by applying suitable numerical filtering at certain places of the discretization. While this filtering technique is effective for two-dimensional (2-D) periodic fluid interfaces, it does not apply to nonperiodic fluid interfaces. Moreover, using the filtering technique alone does not seem to be sufficient to stabilize 3-D fluid interfaces.

Here we introduce a new stabilizing technique for boundary integral methods for water waves which applies to nonperiodic and 3-D interfaces. A stabilizing term is added to the boundary integral method which exactly cancels the destabilizing term produced by the point vortex method approximation to the leading order. This modified boundary integral method still has the same order of accuracy as the point vortex method. A detailed stability analysis is presented for the point vortex method for 2-D water waves. The effect of various stabilizing terms is illustrated through careful numerical experiments.

     

On some unsteady flows of a non-Newtonian fluid

Some properties of unsteady unidirectional flows of a fluid of second grade are considered for flows impulsively started from rest by the motion of a boundary or two boundaries or by sudden application of a pressure gradient. Flows considered are: unsteady flow over a plane wall, unsteady Couette flow, flow between two parallel plates suddenly set in motion with the same speed, flow due to one rigid boundary moved suddenly and one being free, unsteady Poiseuille flow and unsteady generalized Couette flow. The results obtained are compared with those of the exact solutions of the Navier–Stokes equations. It is found that the stress at time zero on the stationary boundary for the flows generated by impulsive motion of a boundary or two boundaries is finite for a fluid of second grade and infinite for a Newtonian fluid. Furthermore, it is shown that for unsteady Poiseuille flow the stress at time zero on the boundary is zero for a Newtonian fluid, but it is not zero for a fluid of second grade.     

Existence of solutions to the nonlinear,singular second order Bohr boundary value problems

This paper investigates the existence of solutions for nonlinear systems of second order, singular boundary value problems (BVPs) with Bohr boundary conditions. A key application that arises from this theory is the famous Thomas–Fermi equations for the model of the atom when it is in a neutral state. The methodology in this paper uses an alternative and equivalent BVP, which is in the class of resonant singular BVPs, and thus this paper obtains novel results by implementing an innovative differential inequality, Lyapunov functions and topological techniques. This approach furnishes new results in the area of singular BVPs for a priori bounds and existence of solutions, where the BVP has unrestricted growth conditions and subject to the Bohr boundary conditions. In addition, the results can be relaxed and hold for the non-singular case too.     

On the Parametric Cubic Spline Approach for Solving Second Order Boundary Value Problems

In this article some comments on the paper “parametric cubic spline approach to the solution of a system of second order boundary value problems” in (Khan and Aziz, J. Optim. Theory Appl. 118:45–54, 2003) are given. This paper concerns with a numerical method for solving a second order boundary value problem associated with obstacle, unilateral and contact problems. Corrections are given for the convergence analysis of the numerical method and the computational experiments.     

Exterior electromagnetic shaping using wavelet BEM

The present paper is devoted to exterior electromagnetic shaping in two dimensions. We model the conductors by regular densities which leads to a finite objective and allows a line‐search. In order to compute the surface pressure we optimize an Augmented Lagrangian by a Newton method using a second‐order approach for the Lagrange multiplier. Since the underlying state function satisfies an exterior boundary value problem, we compute first and second order derivatives of its boundary data by boundary integral equations which are numerically solved by a fast wavelet Galerkin scheme. Numerical results prove that we succeeded in finding a fast and robust algorithm for solving the considered class of problems. Copyright © 2004 John Wiley & Sons, Ltd.