Some Closed-Form Solutions of Burgers’ Equation

This paper shows how to use the method of quasisolutions to construct exact solutions to Burgers’ equation. A function υ=υ(x, y) is called a quasisolution of a PDE in case there exists a function φ (not a constant function) of one variable so that u(x, y)=φ(υ(x, y)) is a solution of the equation. We prove a theorem giving necessary and sufficient conditions for υ to be a quasisolution to Burgers’ equation. A function φ can then be found explicitly so that u=φ(υ) is an actual solution. Combining this technique with similarity methods, we find a continuum of solutions to Burgers' equation.     

Explicit hybrid six–step,sixth order,fully symmetric methods for solving y ″ = f (x,y)

A family of explicit, fully symmetric, sixth order, six‐step methods for the numerical solution of y′′ = f(x,y) is studied. This family wastes two function evaluations per step and can be derived through interpolation techniques. An interval of periodicity is possessed and the phase lag is of high order. Numerical instabilities usually present in such type of multistep methods were circumvented. We conclude with extended numerical tests over a set of problems justifying our effort of dealing with the new methods.     

On the Behavior of a Nonstationary Poiseuille Solution as <Emphasis Type="Italic">t</Emphasis> → ∞

A nonstationary Poiseuille solution describing the flow of a viscous incompressible fluid in an infinite cylinder is defined as a solution to an inverse problem for the heat equation. The behavior as t → ∞ of the nonstationary Poiseuille solution corresponding to the prescribed flux F(t) ofthe velocity field is studied. In particular, it is proved that if the flux F(t) tends exponentially to a constant flux F _* then the nonstationary Poiseuille solution tends exponentially as t → ∞ to the stationary Poiseuille solution having the flux F _*.Original Russian Text Copyright © 2005 Pileckas K.__________Translated from Sibirskii Matematicheskii Zhurnal, Vol. 46, No. 4, pp. 890–900, July–August, 2005.     

Generalized Dirichlet‐to‐Neumann Map in Time‐Dependent Domains

We study the heat, linear Schrödinger (LS), and linear KdV equations in the domain l(t) < x < ∞ , 0 < t < T , with prescribed initial and boundary conditions and with l(t) a given differentiable function. For the first two equations, we show that the unknown Neumann or Dirichlet boundary value can be computed as the solution of a linear Volterra integral equation with an explicit weakly singular kernel. This integral equation can be derived from the formal Fourier integral representation of the solution. For the linear KdV equation we show that the two unknown boundary values can be computed as the solution of a system of linear Volterra integral equations with explicit weakly singular kernels. The derivation in this case makes crucial use of analyticity and certain invariance properties in the complex spectral plane. The above Volterra equations are shown to admit a unique solution.     

Data‐Driven Optimal Transport

The problem of optimal transport between two distributions ρ(x) and μ(y) is extended to situations where the distributions are only known through a finite number of samples {x_i} and {y_j}. A weak formulation is proposed, based on the dual of the Kantorovich formulation, with two main modifications: replacing the expected values in the objective function by their empirical means over the {x_i} and {y_j}, and restricting the dual variables u(x) and v(y) to a suitable set of test functions adapted to the local availability of sample points. A procedure is proposed and tested for the numerical solution of this problem, based on a fluidlike flow in phase space, where the sample points play the role of active Lagrangian markers. © 2016 Wiley Periodicals, Inc.     

Phase‐fitted,six‐step methods for solving x′′ = f(t,x)

An explicit, six‐step method of sixth order is presented and tuned for the numerical solution of x′′ = f(t,x). This method is explicit, hybrid, and uses two function evaluations (stages) per step. Its coefficients are varied and depend on the step size. This variance comes from the demand of the method to nullify the phase errors produced when solving the standard simple oscillator. The first and second derivative of this error vanish also. Numerical tests in a set of relevant problems illustrate the efficiency of the newly derived method.     

Fréchet differentiability for a damped sine-Gordon equation

We consider a damped sine-Gordon equation with a variable diffusion coefficient. The goal is to derive necessary conditions for the optimal set of parameters minimizing the objective function J. First, we show that the solution map is continuous under a weak assumption on the topology of the admissible set P. Then the solution map is shown to be weakly Gâteux differentiable on P, implying the Gâteux differentiability of the objective function. Finally we show the Fréchet differentiability of J. The optimal set of parameters is shown to satisfy a bang–bang control law.     

Well‐posedness for the incompressible magneto‐hydrodynamic system

This paper is concerned with well‐posedness of the incompressible magneto‐hydrodynamics (MHD) system. In particular, we prove the existence of a global mild solution in BMO^?1 for small data which is also unique in the space C([0, ∞); BMO^?1). We also establish the existence of a local mild solution in bmo^?1 for small data and its uniqueness in C([0, T); bmo^?1). In establishing our results an important role is played by the continuity of the bilinear form which was proved previously by Kock and Tataru. In this paper, we give a new proof of this result by using the weighted L^p‐boundedness of the maximal function. Copyright © 2006 John Wiley & Sons, Ltd.     

Global existence to a reaction–diffusion equation through a self‐similar‐like upper solution

A demonstration method is presented, which will ensure the existence of positive global solutions in time to the reaction–diffusion equation ?u_t+Δu+u^p=0 in ?ⁿ×[0, ∞), for exponents p?3 and space dimensions n?3. This method does not require the initial value to have a specific uniform smallness condition, but rather to satisfy a bell‐like form. The method is based on a specific upper solution, which models the diffusion process of the heat equation. The upper solution is not self‐similar, but does have a self‐similar‐like form. After transforming the reaction–diffusion problem into an equivalent one, whose initial value is uniformly very small, a local solution is obtained in the time interval [0, 1] by the use of this upper solution. This local solution is then extended to [0, ∞) through an infinite sequence of extensions. At each step, an appropriate change of variables will transform the extension into a problem nearly identical to the local problem in [0, 1]. These transformations exploit the diffusive and self‐similar‐like nature of the upper solution. Copyright © 2009 John Wiley & Sons, Ltd.     

A <Emphasis Type="Italic">d</Emphasis>-person Differential Game with State Space Constraints

We consider a network of d companies (insurance companies, for example) operating under a treaty to diversify risk. Internal and external borrowing are allowed to avert ruin of any member of the network. The amount borrowed to prevent ruin is viewed upon as control. Repayment of these loans entails a control cost in addition to the usual costs. Each company tries to minimize its repayment liability. This leads to a d -person differential game with state space constraints. If the companies are also in possible competition a Nash equilibrium is sought. Otherwise a utopian equilibrium is more appropriate. The corresponding systems of HJB equations and boundary conditions are derived. In the case of Nash equilibrium, the Hamiltonian can be discontinuous; there are d interlinked control problems with state constraints; each value function is a constrained viscosity solution to the appropriate discontinuous HJB equation. Uniqueness does not hold in general in this case. In the case of utopian equilibrium, each value function turns out to be the unique constrained viscosity solution to the appropriate HJB equation. Connection with Skorokhod problem is briefly discussed.     

A best approximation for the solution of one‐dimensional variable‐coefficient Burgers' equation

In this article, an iterative method for the approximate solution to one‐dimensional variable‐coefficient Burgers' equation is proposed in the reproducing kernel space W_(3,2). It is proved that the approximation w_n(x,t) converges to the exact solution u(x,t) for any initial function w₀(x,t) ε W_(3,2), and the approximate solution is the best approximation under a complete normal orthogonal system . Moreover the derivatives of w_n(x,t) are also uniformly convergent to the derivatives of u(x,t).© 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

On the ‘one and one-half dimensional’ relativistic Vlasov–Maxwell system

The time evolution of a collisionless plasma is studied in the case when the Viasov density ? is a function of the time, one space variable and two velocity variables. The electromagnetic fields E, B also have a special structure, and the magnetic field B is non-trivial. It is shown that smooth, consistent initial values generate a uniquc smooth global solution.     

An application of <Emphasis Type="Italic">H</Emphasis>-differentiability to nonnegative and unrestricted generalized complementarity problems

This paper deals with nonnegative nonsmooth generalized complementarity problem, denoted by GCP(f,g). Starting with H-differentiable functions f and g, we describe H-differentials of some GCP functions and their merit functions. We show how, under appropriate conditions on H-differentials of f and g, minimizing a merit function corresponding to f and g leads to a solution of the generalized complementarity problem. Moreover, we generalize the concepts of monotonicity, P ₀-property and their variants for functions and use them to establish some conditions to get a solution for generalized complementarity problem. Our results are generalizations of such results for nonlinear complementarity problem when the underlying functions are C ¹, semismooth, and locally Lipschitzian.     

A new smoothing Broyden-like method for solving nonlinear complementarity problem with a <Emphasis Type="Italic">P</Emphasis><Subscript>0</Subscript>-function

In this paper, we propose a new smoothing Broyden-like method for solving nonlinear complementarity problem with P ₀ function. The presented algorithm is based on the smoothing symmetrically perturbed minimum function φ(a, b) = min{a, b} and makes use of the derivative-free line search rule of Li et al. (J Optim Theory Appl 109(1):123–167, 2001). Without requiring any strict complementarity assumption at the P ₀-NCP solution, we show that the iteration sequence generated by the suggested algorithm converges globally and superlinearly under suitable conditions. Furthermore, the algorithm has local quadratic convergence under mild assumptions. Some numerical results are also reported in this paper.     

Pairs of positive solutions for nonlinear elliptic equations with the <Emphasis Type="Italic">p</Emphasis>-Laplacian and a nonsmooth potential

In this paper we examine a nonlinear elliptic problem driven by the p-Laplacian differential operator and with a potential function which is only locally Lipschitz, not necessarily C¹ (hemivariational inequality). Using the nonsmooth critical point theory of Chang, we obtain two strictly positive solutions. One solution is obtained by minimization of a suitable modification of the energy functional. The second solution is obtained by generalizing a result of Brezis-Nirenberg about the local C¹₀-minimizers versus the local H¹₀-minimizers of a C¹-functional. Mathematics Subject Classification (2000) 35J50, 35J85, 35R70     

A new non-interior continuation method for <Emphasis Type="Italic">P</Emphasis><Subscript>0</Subscript>-NCP based on a SSPM-function

In this paper, we consider a new non-interior continuation method for the solution of nonlinear complementarity problem with P ₀-function (P ₀-NCP). The proposed algorithm is based on a smoothing symmetric perturbed minimum function (SSPM-function), and one only needs to solve one system of linear equations and to perform only one Armijo-type line search at each iteration. The method is proved to possess global and local convergence under weaker conditions. Preliminary numerical results indicate that the algorithm is effective.     

Une résolution élémentaire de l'équation diophantienne x = y − 1

A complete solution in integers of the well-known diophantine equation y² = x³ + 1 is given by using elementary arguments based on the Pell-Fermat equation X² − 3Y² = 1.     

Duality for <Emphasis Type="Italic">ε</Emphasis>-Variational Inequality

In this paper, following the method in the proof of the composition duality principle due to Robinson and using some basic properties of the ε-subdifferential and the conjugate function of a convex function, we establish duality results for an ε-variational inequality problem. Then, we give Fenchel duality results for the ε-optimal solution of an unconstrained convex optimization problem. Moreover, we present an example to illustrate our Fenchel duality results for the ε-optimal solutions. The authors thank the referees for valuable suggestions and comments. This work was supported by Grant No. R01-2003-000-10825-0 from the Basic Research Program of KOSEF.     

Analysis of positive solution and Hyers‐Ulam stability for a class of singular fractional differential equations with p‐Laplacian in Banach space

This paper deals with 2 core aspects of fractional calculus including existence of positive solution and Hyers‐Ulam stability for a class of singular fractional differential equations with nonlinear p‐Laplacian operator in Caputo sense. For these aims, the suggested problem is converted into an integral equation via Green function , for ε∈(n−1,n], where n≥4. Then, the Green function is examined whether it is increasing or decreasing and positive or negative function. After these properties, some classical fixed‐point theorems are used for the existence of positive solution. Hyers‐Ulam stability of the proposed problem is also considered. For the application of the results, an expressive example is included.     

Well-posedness in <Emphasis Type="Italic">BV</Emphasis><Subscript><Emphasis Type="Italic">t</Emphasis></Subscript> and convergence of a difference scheme for continuous sedimentation in ideal clarifier-thickener units

Summary. We consider a scalar conservation law modeling the settling of particles in an ideal clarifier-thickener unit. The conservation law has a nonconvex flux which is spatially dependent on two discontinuous parameters. We suggest to use a Krukov-type notion of entropy solution for this conservation law and prove uniqueness (L¹ stability) of the entropy solution in the BV_t class (functions W(x,t) with _tW being a finite measure). The existence of a BV_t entropy solution is established by proving convergence of a simple upwind finite difference scheme (of the Engquist-Osher type). A few numerical examples are also presented.Mathematics Subject Classification (2000):35L65, 35R05, 65M06, 76T20