Transport equation with boundary conditions of hysteresis type

This paper is an advanced extension of the work reported in (Nonlinear Anal. 2005; 63 :1467–1473). A transport equation that describes the propagation of a substance in a moving fluid or gas is considered. The equation contains the transient, convection, and diffusion terms. The problem is formulated in a bounded domain provided with an inlet and an outlet for the fluid or gas flow. The crucial point of the problem setting is a hysteresis‐type condition posed on an active part of the boundary. This condition reflects the nondecreasing accumulation with saturation of the transported substance at each point of the active boundary part. We prove the existence and uniqueness of solutions to this problem, study the regularity properties of solutions, and perform numerical simulations that clarify the behavior of the model. Comparing with the results of (Nonlinear Anal. 2005; 63 :1467–1473), the advancement of this work consists in accounting for the motion of the fluid or gas and posing inlet and outlet boundary conditions. Copyright © 2009 John Wiley & Sons, Ltd.     

A tracking method for gas flow into vacuum based on the vacuum Riemann problem

In this paper, a tracking method is proposed for the expansion of gas flow into vacuum which may be combined with numerical methods for the equations of gas dynamics, the Euler equations. This tracking prevents the difficulties of the numerical approximation introduced by the vacuum as a region where the Euler equations are not valid due to the failure of the continuum assumption. The tracking algorithm is based on the exact or an approximate solution of the vacuum Riemann problem. This is the initial value problem with two constant states, one being the gas and the other the vacuum state, and a limit case of the usual Riemann problem. In this approach, the gas–vacuum boundary is sharply resolved within one mesh interval. For a test problem, the numerical results of gas flow into vacuum are presented which indicate that the gas vacuum boundary is captured very well.     

拟牛顿流的一种三变量域模型的有限元方法的数值分析

0．引言目前，涉及高温条件下材料蠕变性质的粘弹性流动问题已引起人们广泛的研究兴趣，不少文章讨论了如何对其进行数值求解（见［1］－－［41），首先，人们研究了较简单的仅以速度，压力两个变量来表述此现象的模型问题（如［1，2］）等．鉴于应力变量在材料性质方面的特殊重要性，最近J．Baxanzer等人在［3］中首次对应力满足幂函数规律的蠕变流研究了包含应力、速度和压力三种变量的模型问题的有限元逼近，当粘性的牛顿部分为零时（详见下述）在假定速度与应力、速度与压力有限元空间之间同时满足两种＊＊B条件以后，证明了有限元解…     

A study on preconditioning iterative methods and relativity of pressure in the numerical calculation of fluid flow

This article presents the effect of preconditioning iterative methods on boundary conditions of the pressure‐correction in the numerical computation of fluid flow with known velocity components on all boundaries using the SIMPLE algorithm. In such computation, a set of solutions of the pressure‐correction is indefinite, because only the Neumann condition is imposed on all boundaries. However, solutions become unique if the value of pressure‐correction is fixed at least on one boundary point, and the Dirichlet condition is additionally imposed. Though both conditions must give exactly the same velocity and temperature fields, this problem arises from the relativity of the pressure. The mathematical illustration for this problem is provided using the numerical computation of the natural convection in an enclosure. It is concluded that the preconditioner adopted and the condition that only the Neumann condition on all boundaries is given are effective to reduce the number of iterations in solving the linear system of equations of the pressure‐correction at the computation of the natural convection. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2004     

Boundary Value Problems for the Elliptic Sine-Gordon Equation in a Semi-strip

We study boundary value problems posed in a semistrip for the elliptic sine-Gordon equation, which is the paradigm of an elliptic integrable PDE in two variables. We use the method introduced by one of the authors, which provides a substantial generalization of the inverse scattering transform and can be used for the analysis of boundary as opposed to initial-value problems. We first express the solution in terms of a 2×2 matrix Riemann–Hilbert problem whose “jump matrix” depends on both the Dirichlet and the Neumann boundary values. For a well posed problem one of these boundary values is an unknown function. This unknown function is characterised in terms of the so-called global relation, but in general this characterisation is nonlinear. We then concentrate on the case that the prescribed boundary conditions are zero along the unbounded sides of a semistrip and constant along the bounded side. This corresponds to a case of the so-called linearisable boundary conditions, however, a major difficulty for this problem is the existence of non-integrable singularities of the function q _y at the two corners of the semistrip; these singularities are generated by the discontinuities of the boundary condition at these corners. Motivated by the recent solution of the analogous problem for the modified Helmholtz equation, we introduce an appropriate regularisation which overcomes this difficulty. Furthermore, by mapping the basic Riemann–Hilbert problem to an equivalent modified Riemann–Hilbert problem, we show that the solution can be expressed in terms of a 2×2 matrix Riemann–Hilbert problem whose “jump matrix” depends explicitly on the width of the semistrip L, on the constant value d of the solution along the bounded side, and on the residues at the given poles of a certain spectral function denoted by h(λ). The determination of the function h remains open.     

Far field boundary conditions for incompressible flows computation

Many far field boundary conditions are proposed in the literature to solve Navier-Stokes equations. It is necessary to distinguish the streamwise or outlet boundary conditions and the spanwise boundary conditions. In the first case the flow crosses the artificial frontier and it is required to avoid reflections that can change significantly the flow. In the second case the Navier-slip boundary condition is often used but if the frontier is not far enough the boundary is both inlet and outlet. Thus the Navier-slip boundary condition is not well suited as it imposes no flux through the frontier. The aim of this work is to compare some well-known boundary conditions, to quantify to which extend the artificial frontier can be close to the bodies in two- and three-dimensions and to take into account the flow rate through the spanwise directions.     

An Upwind Difference Scheme and Its Corresponding Boundary Scheme

An upwind difference scheme was given by the author in [5] for the numerical solution of steady-state problems. The present work studies this upwind scheme and its corresponding boundary scheme for the numerical solution of unsteady problems. For interior points the difference equations are approximations of the characteristic relations; for boundary points difference equatons are approximations of the characteristicrelations corresponding to the outgoing characteristics and the "non-reflecting" boundary conditions. Calculation of a Riemann problem in a finite computational region yields promising numerical results.     

On explicit,purely analytic solutions of off-centered stagnation flow towards a rotating disc by means of HAM

In this article, the off-centered stagnation flow towards a rotating disc is studied analytically. The governing non-linear equations and their associated boundary conditions are transformed into coupled non-linear ordinary differential equations. The series solution of the problem is obtained by utilizing the HAM. The convergence of the obtained series solutions is carefully checked. Graphical results are presented to investigate the influence of the rotation ratio on the flow field. An important point to note is that the non-alignment complicates the flow field and surface shear, but does not affect the torque. It is noted that the behavior of the HAM solution for velocity components is in good agreement with the numerical solution given in reference [C.Y. Wang, Off-centered stagnation flow towards a rotating disc, Int. J. Eng. Sci. 46 (2008) 391–396].     

层流管流入口压力突然升高引起的水力瞬变的理论分析

研究了层流状态下管道入口压力突然升高引起的水力瞬变过程,建立了瞬态压力分布的偏微分方程和初边值条件,用分离变量法求得了压力的理论解．根据压力和流量间的约束关系,得到了关于流量的偏微分方程和初边值条件,用分离变量求得了瞬变过程流量分布理论解．最后,用特征线法(MOC)对该问题进行了数值求解,理论解和数值解吻合很好．     

Solution for a problem of linear plane elasticity with mixed boundary conditions on an ellipse by the method of boundary integrals

A numerical boundary integral scheme is proposed for the solution of the system of field equations of plane, linear elasticity in stresses for homogeneous, isotropic media in the domain bounded by an ellipse under mixed boundary conditions. The stresses are prescribed on one half of the ellipse, while the displacements are given on the other half. The method relies on previous analytical work within the Boundary Integral Method [1], [2].The considered problem with mixed boundary conditions is replaced by two subproblems with homogeneous boundary conditions, one of each type, having a common solution. The equations are reduced to a system of boundary integral equations, which is then discretized in the usual way and the problem at this stage is reduced to the solution of a rectangular linear system of algebraic equations. The unknowns in this system of equations are the boundary values of four harmonic functions which define the full elastic solution inside the domain, and the unknown boundary values of stresses or displacements on proper parts of the boundary.On the basis of the obtained results, it is inferred that the tangential stress component on the fixed part of the boundary has a singularity at each of the two separation points, thought to be of logarithmic type. A tentative form for the singular solution is proposed to calculate the full solution in bulk directly from the given boundary conditions using the well-known Boundary Collocation Method. It is shown that this addition substantially decreases the error in satisfying the boundary conditions on some interval not containing the singular points.The obtained results are discussed and boundary curves for unknown functions are provided, as well as three-dimensional plots for quantities of practical interest. The efficiency of the used numerical schemes is discussed, in what concerns the number of boundary nodes needed to calculate the approximate solution.     

On the numerical solution of the flow between a rotating and a stationary disk

The flow between two co-axial, infinite disks, one rotating with constant angular velocity and one stationary is treated in this paper. The problem is reduced to that of finding the solution of a two-point boundary value for a sixth order nonlinear ordinary differential equation and three boundary conditions at each of a finite interval. The numerical solutions are obtained by using a fourth order Runge-Kutta integration scheme in modification due to Gill and in conjunction with a modified shooting method to correct the initial guesses at one boundary. The numerical calculations for different Reynolds numbers are carried out. The results obtained by this method are compared with available results. The comparison shows excellent agreement.     

Approaching Chaplygin pressure limit of solutions to the Aw–Rascle model

This paper studies the limit of solutions to the Aw–Rascle model as the pressure tends to the Chaplygin gas pressure. For concreteness, the pressure is taken as a modified Chaplygin gas pressure. Firstly, the Riemann problem for the Aw–Rascle model with the modified Chaplygin gas pressure is solved constructively. Secondly, it is shown that as the pressure tends to the Chaplygin gas pressure, some Riemann solutions containing a shock and a contact discontinuity tend to a delta-shock solution, whose propagation speed and strength are different from those of delta-shock solution to the Aw–Rascle model with a Chaplygin gas pressure. Besides, it is also proven that the rest Riemann solutions converge to a two-contact-discontinuity solution, which is exactly the solution to the Aw–Rascle model with a Chaplygin gas pressure. Thirdly, some numerical results are presented to exhibit the process of formation of delta-shocks.     

Inverse problem of hydrodynamics for doubly connected domain

We consider an inverse problem of hydrodynamics for flow past pair of aerofoils. We find a general form of its solution. The key part of problem’s solving is to determine numerical parameters defining flow domain and complex velocity in it up to conformal mapping (the parameters problem). The solvability of parameters problem is proved for various flow schemes. For that we essentially use the interpretation of the problem in terms of Riemann surface and the Riemann surfaces theory.     

Delta-shocks and vacuums in zero-pressure gas dynamics by the flux approximation

In this paper,firstly,by solving the Riemann problem of the zero-pressure flow in gas dynamics with a flux approximation,we construct parameterized delta-shock and constant density solutions,then we show that,as the flux perturbation vanishes,they converge to the delta-shock and vacuum state solutions of the zero-pressure flow,respectively.Secondly,we solve the Riemann problem of the Euler equations of isentropic gas dynamics with a double parameter flux approximation including pressure.Furthermore,we rigorously prove that,as the two-parameter flux perturbation vanishes,any Riemann solution containing two shock waves tends to a delta-shock solution to the zero-pressure flow;any Riemann solution containing two rarefaction waves tends to a two-contact-discontinuity solution to the zero-pressure flow and the nonvacuum intermediate state in between tends to a vacuum state.Finally,numerical results are given to present the formation processes of delta shock waves and vacuum states.     

On solution uniqueness of elliptic boundary value problems

In this paper, we consider the problem of solution uniqueness for the second order elliptic boundary value problem, by looking at its finite element or finite difference approximations. We derive several equivalent conditions, which are simpler and easier than the boundedness of the entries of the inverse matrix given in Yamamoto et al., [T. Yamamoto, S. Oishi, Q. Fang, Discretization principles for linear two-point boundary value problems, II, Numer. Funct. Anal. Optim. 29 (2008) 213–224]. The numerical experiments are provided to support the analysis made. Strictly speaking, the uniqueness of solution is equivalent to the existence of nonzero eigenvalues in the corresponding eigenvalue problem, and this condition should be checked by solving the corresponding eigenvalue problems. An application of the equivalent conditions is that we may discover the uniqueness simultaneously, while seeking the approximate solutions of elliptic boundary equations.     

A complete global solution to the pressure gradient equation

We study the domain of existence of a solution to a Riemann problem for the pressure gradient equation in two space dimensions. The Riemann problem is the expansion of a quadrant of gas of constant state into the other three vacuum quadrants. The global existence of a smooth solution was established in Dai and Zhang [Z. Dai, T. Zhang, Existence of a global smooth solution for a degenerate Goursat problem of gas dynamics, Arch. Ration. Mech. Anal. 155 (2000) 277-298] up to the free boundary of vacuum. We prove that the vacuum boundary is the coordinate axes.     

Two-parameter extremum problems of boundary control for stationary thermal convection equations

Two-parameter extremum problems of boundary control are formulated for the stationary thermal convection equations with Dirichlet boundary conditions for velocity and with mixed boundary conditions for temperature. The cost functional is defined as the root mean square integral deviation of the desired velocity (vorticity, or pressure) field from one given in some part of the flow region. Controls are the boundary functions involved in the Dirichlet condition for velocity on the boundary of the flow region and in the Neumann condition for temperature on part of the boundary. The uniqueness of the extremum problems is analyzed, and the stability of solutions with respect to certain perturbations in the cost functional and one of the functional parameters of the original model is estimated. Numerical results for a control problem associated with the minimization of the vorticity norm aimed at drag reduction are discussed.     

Riemann and Goursat step data problems for extensible strings

Riemann and Goursat step data problems for extensible nonlinear elastic strings are solved in the class of regulated functions. In the first paragraph, the solution to the simplest initial and boundary value problem, i.e., Goursat problem in strains, is constructed. This solution points out four vector-valued functions of a vector variable, which are the tools used in solving the Goursat problem in velocity and the Riemann problem.     

Meshless solution of two-dimensional incompressible flow problems using the radial basis integral equation method

The two-dimensional incompressible fluid flow problems governed by the velocity–vorticity formulation of the Navier–Stokes equations were solved using the radial basis integral (RBIE) equation method. The RBIE is a meshless method based on the multi-domain boundary element method with overlapping subdomains. It solves at each node for the potential and its spatial derivatives. This feature of the RBIE is advantageous in solving the velocity–vorticity formulation of the Navier–Stokes equations since the calculated velocity gradients can be used to compute the vorticity that is prescribed as a boundary condition to the vorticity transport equation. The accuracy of the numerical solution was examined by solving the test problem with known analytical solution. Two benchmark problems, i.e. the lid driven cavity flow and the thermally driven cavity flow were also solved. The numerical results obtained using the RBIE showed very good agreement with the benchmark solutions.