On the solution of Stokes equations for plane boundaries

In a recent paper(Li et al., Acta Mech. Sin. 31,32–44, 2015), the authors claimed that the general solution of steady Stokes flows can be compactly expressed using only two harmonic functions. They present two cases of a flat plate translating through a viscous fluid. The present paper shows that such a two-harmonic solution does not describe the rotation of a circular plate in an unbounded fluid and thus confirms that at least three independent harmonics are required to express the general solution of Stokes equations.     

A general solution for static piezothermoelastic problems of crystal class 6mm solids

In this paper, a general solution for three-dimensional static piezothermoelastic problems of crystal class 6mm solids is presented. The general solution involves four piezoelastic potential functions and a piezothermoelastic potential function, of which four piezoelastic potential functions are governed by weighted harmonic differential equations. Compared with the general solution given by Ashida et al., in which seven potential functions are introduced, the general solution proposed in the present paper is more rigorously derived. Moreover, it has a simple form and is convenient for application. As an illustrative example, the problem of a pyroelectric half-space subjected to axisymmetric heating is studied. Numerical results of displacements stresses, electric potential and electric displacements are obtained for a cadmium selenide half-space. Thermally induced displacements and stress distribution are compared with those obtained for the same material without piezoelectric and pyroelectric effects. Project supported by the National Natural Science Foundation of China (19872060).     

On Single Mode Forcing of the 2D-NSE

We examine how the global attractor of the 2D periodic Navier- Stokes equations projects in the energy–enstrophy-plane when the force is an eigenvector of the Stokes operator. We find sharper bounds than in Dascaliuc et al. (J Dyn Diff Equat 17:643–736, 2005) for a general force. We also find regions of the plane in which the energy and enstrophy must decrease, and calculate the curvature of the projected solution at certain initial data. Finally, we obtain restrictions on solutions that project onto a single point in the plane.     

Representing general solution of equations in theory of elasticity by harmonic functions

The general solution of the equations in the theory of elasticity is represented by seven harmonic functions, where there are only three harmonic functions independent of each other and every one of them has certain mechanics meaning. The examples applying the general solution to solve several simple inverse problems in elastostatics are presented.     

Solving the three-dimensional equations of the linear theory of elasticity

The three-dimensional Lamé equations are solved using Cartesian and curvilinear orthogonal coordinates. It is proved that the solution includes only three independent harmonic functions. The general solution of equations of elasticity for stresses is found. The stress tensor is expressed in both coordinate systems in terms of three harmonic functions. The general solution of the problem of elasticity in cylindrical coordinates is presented as an example. The three-dimensional stress–strain state of an elastic cylinder subjected, on the lateral surface, to arbitrary forces represented by a series of eigenfunctions is determined. An axisymmetric problem for a finite cylinder is solved numerically     

Boundary layer sensitivity and receptivitySensibilité et réceptivité d'une couche limite

The relation between the receptivity and the sensitivity of the incompressible flow in the boundary layer over a flat plate to harmonic perturbations is determined. Receptivity describes the birth of a disturbance, whereas sensitivity is a concept of larger breath, describing the modification incurred by the state of a system as a response to parametric variations. The governing equations ruling the system's state are the non-local stability equations. Receptivity and sensitivity functions can be obtained from the solution of the adjoint system of equations. An application to the case of Tollmien–Schlichting waves spatially developing in a flat plate boundary layer is studied. To cite this article: C. Airiau et al., C. R. Mecanique 330 (2002) 259–265.     

A general solution for Stokes flow and its application to the problem of a rigid plate translating in a fluid

A general solution for 3D Stokes flow is given which is different from, and more compact than the existing ones and more compact than them in that it involves only two scalar harmonic functions. The general solution deduced is combined with the potential theory method to study the Stokes flow induced by a rigid plate of arbitrary shape translating along the direction normal to it in an unbounded fluid.The boundary integral equation governing this problem is derived. When the plate is elliptic, exact analytical results are obtained not only for the drag force but also for the velocity distributions. These results include and complete the ones available for a circular plate. Numerical examples are provided to illustrate the main results for circular and elliptic plates. In particular, the elliptic eccentricity of a plate is shown to exhibit significant influences.     

Optimal Results for the Nonhomogeneous Initial-Boundary Value Problem for the Two-Dimensional Navier–Stokes Equations

In this work we study the fully nonhomogeneous initial boundary value problem for the two-dimensional time-dependent Navier–Stokes equations in a general open space domain in R² with low regularity assumptions on the initial and the boundary value data. We show that the perturbed Navier–Stokes operator is a diffeomorphism from a suitable function space onto its own dual and as a corollary we get that the Navier–Stokes equations are uniquely solvable in these spaces and that the solution depends smoothly on all involved data. Our source data space and solution space are in complete natural duality and in this sense, without any smallness assumptions on the data, we solve the equations for data with optimally low regularity in both space and time.     

Averaging method for the solution of non-linear differential equations with periodic non-harmonic solutions

While Krylov and Bogolyubov used harmonic functions in their averaging method for the approximate solution of weakly non-linear differential equations with oscillatory solution, we apply a similar averaging technique using Jacobi elliptic functions. These functions are also periodic and are exact solutions of strongly non-linear differential equations. The method is used to solve non-linear differential equations with linear and non-linear small dissipative terms and/or with time dependent parameters. It is also shown that quite general dissipative terms can be transformed into time-dependent parameters. As a special example, the Langevin (collisional) equation of motion of electrons in a neutralizing ion background under the influence of a time and space-dependent electric field is presented. The method may also be used for non-linear control theory, dynamic and parametric stabilization of non-linear oscillations in plasma physics, etc.     

3-D DYNAMIC RESPONSE OF TRANSVERSELY ISOTROPIC SATURATED SOILS

Introduction Thedynamicanalysisofsaturatedporoelasticmediaisofbroadmeaninginmanyaspectssuchasseismology,earthquakeengineering,soilmechanicsandgeophysics,etc.Thetransverselyisotropicsaturatedsoilsissolid fluidcoupledtwophasemediainwhichthe soilskeletondisp…     

Some methods of training radial basis neural networks in solving the Navier‐Stokes equations

The purpose of this research is to analyze the application of neural networks and specific features of training radial basis functions for solving 2‐dimensional Navier‐Stokes equations. The authors developed an algorithm for solving hydrodynamic equations with representation of their solution by the method of weighted residuals upon the general neural network approximation throughout the entire computational domain. The article deals with testing of the developed algorithm through solving the 2‐dimensional Navier‐Stokes equations. Artificial neural networks are widely used for solving problems of mathematical physics; however, their use for modeling of hydrodynamic problems is very limited. At the same time, the problem of hydrodynamic modeling can be solved through neural network modeling, and our study demonstrates an example of its solution. The choice of neural networks based on radial basis functions is due to the ease of implementation and organization of the training process, the accuracy of the approximations, and smoothness of solutions. Radial basis neural networks in the solution of differential equations in partial derivatives allow obtaining a sufficiently accurate solution with a relatively small size of the neural network model. The authors propose to consider the neural network as an approximation of the unknown solution of the equation. The Gaussian distribution is used as the activation function.     

Stokes flow with injection in a plane duct

Summary In this paper flow at low Reynolds number is studied in a plane duct with injection along its boundaries, by determining the Green function for the Stokes equation in the strip: in this way it is possible to calculate the fluid-dynamic field for any distribution of injection. To obtain the bi-harmonic function that governs this problem we write the stream function in terms of two harmonic functions f and g, represent these functions using the relevant Green function and write the integral equations that determine f and g. The solution is obtained by means of the Fourier transform technique. An asymptotic representation of this solution is also given.

---

Sommario In questo lavoro si studia il campo fluidodinamico determinato) dall'apporto di fluido in un condotto in cui scorre un fluido a piccoli numeri di Reynolds, determinando la funzione di Green per l'equazione di Stokes nella striscia: in tal modo si può calcolare il campo fluidodinamico per una qualsiasi distribuzione di iniezione. Per determinare la funzione biarmonica con derivata longitudinale infinita in alcuni punti della frontiera che governa questo problema, si scrive l'incognita in termini di funzioni armoniche, si utilizza la funzione di Green per rappresentare tali funzioni ausiliarie e si scrivono le equazioni integrali che le determinano. La soluzione di queste equazioni integrali viene ottenuta mediante la trasformata di Fourier. I risultati vengono presentati sia in forma numerica che con serie asintotiche.

     

Some useful hybrid approaches for predicting aerodynamic noise

In recent years, several numerical studies have shown the feasibility of Direct Noise Computation (DNC) where the turbulent flow and the radiated acoustic field are obtained simultaneously by solving the compressible Navier–Stokes equations. The acoustic radiation obtained by DNC can be used as reference solution to investigate hybrid methods in which the sound field is usually calculated as a by-product of the flow field obtained by a more conventional Navier–Stokes solver. A hybrid approach is indeed of practical interest when only the non-acoustic part of the aerodynamic field is available. In this review, some acoustic analogies or hybrid approaches are revisited in the light of CAA. To cite this article: C. Bailly et al., C. R. Mecanique 333 (2005).     

A Galerkin spectral method based on helical‐wave decomposition for the incompressible Navier–Stokes equations

This paper presents a global Galerkin spectral method for solving the incompressible Navier–Stokes equations in three‐dimensional bounded domains. The method is based on helical‐wave decomposition (HWD), which uses the vector eigenfunctions of the curl operator as orthogonal basis functions. We shall first review the general theory of HWD in an arbitrary simply connected domain, along with some new developments. We then employ the HWD to construct a Galerkin spectral method. The current method innovates the existing HWD‐based spectral method by (a) adding a series of auxiliary fields to the HWD of the velocity field to fulfill the no‐slip boundary condition and to settle the convergence problem of the HWD of the curl fields, and (b) providing a pseudo‐spectral method that utilizes a fast spherical harmonic transform algorithm and Gaussian quadrature to calculate the nonlinear term in the Navier–Stokes equations. The auxiliary fields are uniquely determined by solving the Stokes and Stokes‐like equations under adequate boundary conditions. The implementation of the method under the spherical geometry is presented in detail. Several numerical examples are provided to validate the proposed method. The method can be easily extended to other domains once the helical‐wave bases, which depend only on the geometry of the domains, are available. Copyright © 2015 John Wiley & Sons, Ltd.     

POINT FORCE SOLUTION FOR A TRANSVERSELY ISOTROPIC ELASTIC LAYER

ＰＯＩＮＴＦＯＲＣＥＳＯＬＵＴＩＯＮＦＯＲＡＴＲＡＮＳＶＥＲＳＥＬＹＩＳＯＴＲＯＰＩＣＥＬＡＳＴＩＣＬＡＹＥＲＰＯＩＮＴＦＯＲＣＥＳＯＬＵＴＩＯＮＦＯＲＡＴＲＡＮＳＶＥＲＳＥＬＹＩＳＯＴＲＯＰＩＣＥＬＡＳＴＩＣＬＡＹＥＲ￥ＤｉｎｇＨａｏｊｉａｎｇ（丁...     

DNS of premixed turbulent V-flame: coupling spectral and finite difference methods

To allow for a reliable examination of the interaction between velocity fluctuations, acoustics and combustion, a novel numerical procedure is discussed in which a spectral solution of the Navier–Stokes equations is directly associated to a high-order finite difference fully compressible DNS solver (sixth order PADE). Using this combination of high-order solvers with accurate boundary conditions, simulations have been performed where a turbulent premixed V-shape flame develops in grid turbulence. In the light of the DNS results, a sub-model for premixed turbulent combustion is analyzed. To cite this article: R. Hauguel et al., C. R. Mecanique 333 (2005).     

A reduced‐order approach for optimal control of fluids using proper orthogonal decomposition

In this article, a reduced‐order modeling approach, suitable for active control of fluid dynamical systems, based on proper orthogonal decomposition (POD) is presented. The rationale behind the reduced‐order modeling is that numerical simulation of Navier–Stokes equations is still too costly for the purpose of optimization and control of unsteady flows. The possibility of obtaining reduced‐order models that reduce the computational complexity associated with the Navier–Stokes equations is examined while capturing the essential dynamics by using the POD. The POD allows the extraction of a reduced set of basis functions, perhaps just a few, from a computational or experimental database through an eigenvalue analysis. The solution is then obtained as a linear combination of this reduced set of basis functions by means of Galerkin projection. This makes it attractive for optimal control and estimation of systems governed by partial differential equations (PDEs). It is used here in active control of fluid flows governed by the Navier–Stokes equations. In particular, flow over a backward‐facing step is considered. Reduced‐order models/low‐dimensional dynamical models for this system are obtained using POD basis functions (global) from the finite element discretizations of the Navier–Stokes equations. Their effectiveness in flow control applications is shown on a recirculation control problem using blowing on the channel boundary. Implementational issues are discussed and numerical experiments are presented. Copyright © 2000 John Wiley & Sons, Ltd.     

Stationary Stokes, Oseen and Navier–Stokes Equations with Singular Data

The concept of very weak solution introduced by Giga (Math Z 178:287–329, 1981) for the Stokes equations has hardly been studied in recent years for either the Navier–Stokes equations or the Navier–Stokes type equations. We treat the stationary Stokes, Oseen and Navier–Stokes systems in the case of a bounded open set, connected of class C^1,1{\mathcal{C}^{1,1}} of \mathbbR³{\mathbb{R}^3}. Taking up once again the duality method introduced by Lions and Magenes (Problèmes aus limites non-homogènes et applications, vols. 1 & 2, Dunod, Paris, 1968) and Giga (Math Z 178:287–329, 1981) for open sets of class C^￥{\mathcal{C}^{\infty}} [see also chapter 4 of Necas (Les méthodes directes en théorie des équations elliptiques. (French) Masson et Cie, éd., Paris; Academia, éditeurs, Prague, 1967), which considers the Hilbertian case p = 2 for general elliptic operators], we give a simpler proof of the existence of a very weak solution for stationary Oseen and Navier–Stokes equations when data are not regular enough, based on density arguments and a functional framework adequate for defining more rigourously the traces of non-regular vector fields. In the stationary Navier–Stokes case, the results will be valid for external forces not necessarily small, which lets us extend the uniqueness class of solutions for these equations. Considering more regular data, regularity results in fractional Sobolev spaces will also be discussed for the three systems. All these results can be extended to other dimensions.