Integrability and Self-Similarity in Transient Stimulated Raman Scattering

Summary. The phenomenon of stimulated Raman scattering (SRS) can be described by three coupled PDEs which define the pump electric field, the Stokes electric field, and the material excitation as functions of distance and time. In the transient limit these equations are integrable, i.e., they admit a Lax pair formulation. Here we study this transient limit. The relevant physical problem can be formulated as an initial-boundary value (IBV) problem where both independent variables are on a finite domain. A general method for solving IBV problems for integrable equations has been introduced recently. Using this method we show that the solution of the equations describing the transient SRS can be obtained by solving a certain linear integral equation. It is interesting that this equation is identical to the linear integral equation characterizing the solution of an IBV problem of the sine-Gordon equation in light-cone coordinates. This integral equation can be solved uniquely in terms of the values of the pump and Stokes fields at the entry of the Raman cell. The asymptotic analysis of this solution reveals that the long-distance behavior of the system is dominated by the underlying self-similar solution which satisfies a particular case of the third Painlevé transcendent. This result is consistent with both numerical simulations and experimental observations. We also discuss briefly the effect of frequency mismatch between the pump and the Stokes electric fields. Received December 10, 1996; second revision received October 10, 1997; final revision received January 20, 1998     

Semi‐invariant solutions of the Navier‐Stokes equations

We present new exact solutions and reduced differential systems of the Navier‐Stokes equations of incompressible viscous fluid flow. We apply the method of semi‐invariant manifolds, introduced earlier as a modification of the Lie invariance method. We show that many known solutions of the Navier‐Stokes equations are, in fact, semi‐invariant and that the reduced differential systems we derive using semi‐invariant manifolds generalize previously obtained results that used ad hoc methods. Many of our semi‐invariant solutions solve decoupled systems in triangular form that are effectively linear. We also obtain several new reductions of Navier‐Stokes to a single nonlinear partial differential equation. In some cases, we can solve reduced systems and generate new analytic solutions of the Navier‐Stokes equations or find their approximations, and physical interpretation.     

A variational approach to estimate incompressible fluid flows

A variational approach is used to recover fluid motion governed by Stokes and Navier–Stokes equations. Unlike previous approaches where optical flow method is used to track rigid body motion, this new framework aims at investigating incompressible flows using optical flow techniques. We formulate a minimization problem and determine conditions under which unique solution exists. Numerical results using finite element method not only support theoretical results but also show that Stokes flow forced by a potential are recovered almost exactly.     

Direct solution of Navier–Stokes equations by radial basis functions

The pressure–velocity formulation of the Navier–Stokes (N–S) equation is solved using the radial basis functions (RBF) collocation method. The non-linear collocated equations are solved using the Levenberg–Marquardt method. The primary novelty of this approach is that the N–S equation is solved directly, instead of using an iterative algorithm for the primitive variables. Two flow situations are considered: Couette flow with and without pressure gradient, and 2D laminar flow in a duct with and without flow obstruction. The approach is validated by comparing the Couette flow results with the analytical solution and the 2D results with those obtained using the well-validated CFD-ACE™ commercial package.     

Influence of eddies on heat transfer in Couette flow over undulated substrates

In Couette flows over undulated substrates eddies can be generated under creeping flow conditions. In contrast to free surface flows on undulated substrates even smooth bottom undulations allow for eddy generation due to the kinematical constraints. The subject of our paper is how these flow patterns interact with the temperature field in non–isothermal flows. Our analysis of the thermo–mechanical coupling is focused on the two dominant effects, namely convection and thermoviscosity, whereas dissipation heat, buoyancy and temperature–dependence of the remaining material parameters are neglected. We solve the problem in two steps: First, the influence of the eddies on the convective heat transfer is considered by solving the heat conduction equation with convection. For the velocity field we take the solution resulting analytically from Reynolds' lubrication approximation for the isothermal flow. The thermoviscous feedback of the resulting temperature field to the flow is considered in forthcoming papers. For the construction of the solution an analytical approach based on a nonorthogonal series representation of the fundamental fields and a variational formulation of the field equations is used. The results are visualised and the physical effects they reveal are discussed. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Homogenization of the Stokes Problem With Non-homogeneous Slip Boundary Conditions

We study the Stokes system with non-homogeneous Fourier boundary conditions depending on a parameter, in a domain with periodic inclusions of the size of the period. Following the values of this parameter, we obtain at the limit a Darcy's law, a Brinkmann type equation or a Stokes type equation. We also present a physical model to which the results apply. This model describes the flow of an incompressible viscous fluid through a porous medium under the action of an exterior electric field.     

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 second‐order backward difference time‐stepping scheme for penalized Navier–Stokes equations modeling filtration through porous media

We investigate the stability and convergence of a fully implicit, linearly extrapolated second‐order backward difference time‐stepping scheme for the penalized Navier–Stokes equations modeling filtration through porous media. In the penalization approach, an extended Navier–Stokes equation is used in the entire computational domain with suitable resistance terms to mimic the presence of porous medium. It is widely used as an alternative to the heterogeneous approach in which different types of partial differential equations (PDEs) are used in fluid and porous subregions along with suitable continuity conditions at the interface. However, the introduction of extra resistance terms makes the penalized Navier–Stokes equations more nonlinear. We prove that the linearly extrapolated scheme is unconditionally stable and derive optimal order error estimates without any stability condition. To show feasibility and applicability of the approach, it is used to numerically solve a passive control problem in which flow around a solid body is controlled by adding porous layers on the surface. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 681–705, 2016     

Coupling nonlinear Stokes and Darcy flow using mortar finite elements

We study a system composed of a nonlinear Stokes flow in one subdomain coupled with a nonlinear porous medium flow in another subdomain. Special attention is paid to the mathematical consequence of the shear-dependent fluid viscosity for the Stokes flow and the velocity-dependent effective viscosity for the Darcy flow. Motivated by the physical setting, we consider the case where only flow rates are specified on the inflow and outflow boundaries in both subdomains. We recast the coupled Stokes–Darcy system as a reduced matching problem on the interface using a mortar space approach. We prove a number of properties of the nonlinear interface operator associated with the reduced problem, which directly yield the existence, uniqueness and regularity of a variational solution to the system. We further propose and analyze a numerical algorithm based on mortar finite elements for the interface problem and conforming finite elements for the subdomain problems. Optimal a priori error estimates are established for the interface and subdomain problems, and a number of compatibility conditions for the finite element spaces used are discussed. Numerical simulations are presented to illustrate the algorithm and to compare two treatments of the defective boundary conditions.     

A new approach to the classical Stokes flow problem: Part I Methodology and first-order analytical results

The problem of determining the axisymmetric Stokes flow past an arbitrary body, the boundary shape of which can be represented by an analytic function, is examined by developing an exact method. An appropriate nonorthogonal coordinate system is introduced, and it is shown that the Hilbert space to which the stream function belongs is spanned by the set of Gegenbauer polynomials based on the physical argument that the drag on a body should be finite. The partial differential equation of the original problem is then reduced to two simultaneous vector differential equations. By the truncation of this infinite-dimensional system to the one-dimensional subspace, an explicit analytic solution to the Stokes equation valid for all bodies in question is obtained as a first approximation.     

New Exact Solutions of the Navier–Stokes Equations for Axisymmetric Automodel Flows of Fluids

We investigate self-similar solutions of the Navier–Stokes equations for the axisymmetric flow of a viscous incompressible fluid. The original equations are transformed by the Slezkin method. On the basis of analysis of physical properties of the flow and the Slezkin general equation, we show that, in parallel with the known solutions of this equation, there exist several other solutions with physical meaning. We consider the simplest case of irrotational flows for which current lines may be circles, ellipses, parabolas, and hyperbolas. Unlike the Landau and Squire solutions, these flows are interpreted as nonjet flows of fluid flowing into and out of a homogeneous porous axially symmetric body.     

On the uniqueness of the solution of the two‐dimensional Navier–Stokes equation with a Dirac mass as initial vorticity

We propose two different proofs of the fact that Oseen's vortex is the unique solution of the two‐dimensional Navier–Stokes equation with a Dirac mass as initial vorticity. The first argument, due to C. E. Wayne and the second named author, is based on an entropy estimate for the vorticity equation in self‐similar variables. The second proof is new and relies on symmetrization techniques for parabolic equations. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Convergence and conditioning of a Nyström method for Stokes flow in exterior three-dimensional domains

Convergence and conditioning results are presented for the lowest-order member of a family of Nyström methods for arbitrary, exterior, three-dimensional Stokes flow. The flow problem is formulated in terms of a recently introduced two-parameter, weakly singular boundary integral equation of the second kind. In contrast to methods based on product integration, coordinate transformation and singularity subtraction, the family of Nyström methods considered here is based on a local polynomial correction determined by an auxiliary system of moment equations. The polynomial correction is designed to remove the weak singularity in the integral equation and provide control over the approximation error. Here we focus attention on the lowest-order method of the family, whose implementation is especially simple. We outline a convergence theorem for this method and illustrate it with various numerical examples. Our examples show that well-conditioned, accurate approximations can be obtained with reasonable meshes for a range of different geometries.     

Supercomputer simulation of viscous noncompressible flow past bodies of complex shape

The article reports the results of numerical analysis of viscous thermally conducting noncompressible flow past bodies of complex shape. The effect of physico-chemical processes is ignored. Flow past bodies with various nose cone designs is examined. The process is described using a system of Navier–Stokes equations augmented with the energy equation. The continuity equation is used to control the numerical accuracy. The simulation results are reported in the form of vector fields and surface components of the velocity vector in various channel sections. The results are analyzed for various body configurations.     

Triple-deck regularization of weak laminar hydraulic jumps in a two layer shallow water flow

The inviscid, weakly nonlinear, shallow-water limit of the Navier–Stokes equations leads to a hyperbolic conservation law. In certain cases of a two-layer flow the flux-function is non-convex, thus leading to the possibility of shocks (in physical terms hydraulic jumps) violating the Oleinik-entropy criteria, so called non-classical shocks. To rule out the inadmissible shocks their internal structure is studied, based on an asymptotic approach consistent with the Navier–Stokes equations. This leads to a triple-deck problem with a novel, non-linear interaction equation in the form of a forced, extended KdV-equation. The limit of vanishing and weak influence of the displacement effect is studied analytically, and in addition representative numerical solutions of the full problem are presented. Of particular interest is a solution, which has a pronounced, almost vanishing minimum in he wall shear. Its local structure is studied in some detail. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Singular Functions and Characterizations of Field Concentrations: a Survey

In the presence of closely located inclusions of the extreme material property, the physical fields, such as the electric field and the stress tensor, may be concentrated and arbitrarily large in the narrow region between two inclusions. Recently there has been significant progress on the quantitative characterization of the field concentration in the contexts of electrostatics(Laplace equation), linear elasticity(Lamé system), and viscous flow(Stokes system). This paper is to review such progress in a coherent way.     

Dual hybrid methods for the elasticity and the Stokes problems: a unified approach

Summary. A unified approach to construct finite elements based on a dual-hybrid formulation of the linear elasticity problem is given. In this formulation the stress tensor is considered but its symmetry is relaxed by a Lagrange multiplier which is nothing else than the rotation. This construction is linked to the approximations of the Stokes problem in the primitive variables and it leads to a new interpretation of known elements and to new finite elements. Moreover all estimates are valid uniformly with respect to compressibility and apply in the incompressible case which is close to the Stokes problem. Received June 20, 1994 / Revised version received February 16, 1996     

Numerical solution of the time-dependent incompressible Navier–Stokes equations by piecewise linear finite elements

In this paper we present a new method to solve the 2D generalized Stokes problem in terms of the stream function and the vorticity. Such problem results, for instance, from the discretization of the evolutionary Stokes system. The difficulty arising from the lack of the boundary conditions for the vorticity is overcome by means of a suitable technique for uncoupling both variables. In order to apply the above technique to the Navier–Stokes equations we linearize the advective term in the vorticity transport equation as described in the development of the paper. We illustrate the good performance of our approach by means of numerical results, obtained for benchmark driven cavity problem solved with classical piecewise linear finite element.     

Rapid Solution of Stokes Flow using Multiscale Galerkin BEM

Two issues that arise when solving Stokes flow problems with the hydrodynamical single‐layer potential are addressed. First, the resulting boundary integral equation is singular, and second, discretizations lead to dense matrices. We discuss a well‐posed modified equation which is equivalent for zero net‐flux. Furthermore, we describe a multiscale basis that lead to sparse stiffness matrices. This approach is suitable for complicated geometries and is an extension of our previous work for the Laplace equation.