经过修正的层流流动的流动稳定性问题(Ⅰ)——基本概念和理论

本文提出了经过修正的层流流动的流动稳定性理论,并在文中给出平行剪切流中平均速度的一类修正剖面,使这种理论可用于研究平行剪切流的流动稳定性,指出了流动失稳的一条新的可能途径.     

Model reduction on inertial manifolds for N–S equations approached by multilevel finite element method

Approximate Inertial Manifolds (AIMs) is approached by multilevel finite element method, which can be referred to as a Post-processed nonlinear Galerkin finite element method, and is applied to the model reduction for fluid dynamics, a typical kind of nonlinear continuous dynamic system from viewpoint of nonlinear dynamics. By this method, each unknown variable, namely, velocity and pressure, is divided into two components, that is the large eddy and small eddy components. The interaction between large eddy and small eddy components, which is negligible if standard Galerkin algorithm is used to approach the original governing equations, is considered essentially by AIMs, and consequently a coarse grid finite element space and a fine grid incremental finite element space are introduced to approach the two components. As an example, the flow field of incompressible flows around airfoil is simulated numerically and discussed, and velocity and pressure distributions of the flow field are obtained accurately. The results show that there exists less essential degrees-of-freedom which can dominate the dynamic behaviors of the discretized system in comparison with the traditional methods, and large computing time can be saved by this efficient method. In a sense, the small eddy component can be captured by AIMs with fewer grids, and an accurate result can also be obtained.     

Two-dimensional pursuit-evasion game with penalty on turning rates

An approach, both analytical and numerical, is used to solve a two-dimensional pursuit-evasion game characterized by a difficulty level intermediate between that of thesimple motion game (with freely and instantaneously oriented velocities) and that of thegame of two cars (with lower bounds on curvature radii). Each player's velocity has a constant modulus. The maneuvers are penalized by introducing, in the performance index, an integral term for the squared velocity turning rate.The local problem solution is relatively easy to find: the equations of motion and the adjoint equations can be integrated by means of elliptic functions and integrals. The global problem is more delicate to solve, because of the existence of a dispersal singular surface requiring an important numerical search to be determined. Thesynthesis problem (how to express the optimal strategies as functions of state) is not explicitly solvable, but a numerical approach using successive approximations can be developed. Illustrative interception trajectories are given.The authors are grateful to Mr. J. P. Peltier, Head, Guidance Group, Aerospace Mechanics Division, Systems Department, ONERA, Châtillon, France, for his suggestions and his efficient assistance in the numerical aspect of this study.     

Viscous flow over a shrinking sheet with an arbitrary surface velocity

In this paper, an analytical solution in a closed form for the boundary layer flow over a shrinking sheet is presented when arbitrary velocity distributions are applied on the shrinking sheet. The solutions with seven typical velocity profiles are derived based on a general closed form expression. Such flow is usually not self-similar and the solution can only be implemented when the mass transfer at the wall is prescribed and determined by the moving velocity of the wall. The characteristics of the flows with the typical velocity distributions are discussed and compared with previous similarity solutions. The flow is observed to have quite different behavior from that of the self-similar flow reported in the literature and the results demonstrate distinctive momentum and energy transport characteristics. Some plots of the stream functions are also illustrated to show the difference in flow field between the shrinking sheet and the stretching sheet. An integral approach to solve boundary layer flow over a shrinking or stretching sheet with uncoupled arbitrary surface velocity and wall mass transfer velocity is outlined and the effectiveness of this approach is discussed.     

Numerical solution of the fully non-linear weakly dispersive serre equations for steep gradient flows

We demonstrate a numerical approach for solving the one-dimensional non-linear weakly dispersive Serre equations. By introducing a new conserved quantity the Serre equations can be written in conservation law form, where the velocity is recovered from the conserved quantities at each time step by solving an auxiliary elliptic equation. Numerical techniques for solving equations in conservative law form can then be applied to solve the Serre equations. We demonstrate how this is achieved. The system of conservation equations are solved using the finite volume method and the associated elliptic equation for the velocity is solved using a finite difference method. This robust approach allows us to accurately solve problems with steep gradients in the flow, such as those generated by discontinuities in the initial conditions.The method is shown to be accurate, simple to implement and stable for a range of problems including flows with steep gradients and variable bathymetry.     

Semi-implicit finite volume level set method for advective motion of interfaces in normal direction

In this paper a semi-implicit finite volume method is proposed to solve the applications with moving interfaces using the approach of level set methods. The level set advection equation with a given speed in normal direction is solved by this method. Moreover, the scheme is used for the numerical solution of eikonal equation to compute the signed distance function and for the linear advection equation to compute the so-called extension speed [1]. In both equations an extrapolation near the interface is used in our method to treat Dirichlet boundary conditions on implicitly given interfaces. No restrictive CFL stability condition is required by the semi-implicit method that is very convenient especially when using the extrapolation approach. In summary, we can apply the method for the numerical solution of level set advection equation with the initial condition given by the signed distance function and with the advection velocity in normal direction given by the extension speed. Several advantages of the proposed approach can be shown for chosen examples and application. The advected numerical level set function approximates well the property of remaining the signed distance function during whole simulation time. Sufficiently accurate numerical results can be obtained even with the time steps violating the CFL stability condition.     

An immersed interface method for the Navier-Stokes equations on irregular domains

An augmented method based on a Cartesian grid is proposed for the incompressible Navier Stokes equations on an irregular domain. The irregular domain is embedded into a rectangular one so that a fast Poisson solver can be used in the projection method. Unlike several methods suggested in the literature that set the force strengths as unknowns, which often results an ill-conditioned linear system, we set the jump in the normal derivative of the velocity as the augmented variable. The new approach improve the condition number of the system for the augmented variable significantly. Using the immersed interface method, we achieve second order accuracy for the velocity. Numerical results and comparisons are given to validate the new method. Some interesting new numerical experiments results are also presented. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Wind identification along a flight trajectory,part 3: 2D-dynamic approach

This paper deals with the identification of the wind profile along a flight trajectory by means of a two-dimensional dynamic approach. In this approach, the wind velocity components are computed as the difference between the inertial velocity components and the airspeed components. The airspeed profile as well as the nominal thrust, drag, and lift profiles are obtained from the available DFDR measurements. The actual values of the thrust, drag, and lift are assumed to be proportional to the respective nominal values via multiplicative parameters, called the thrust, drag, and lift factors. The thrust, drag, and lift factors plus the inertial velocity components at impact are determined by matching the flight trajectory computed from DFDR data with the flight trajectory available from ATCR data. This leads to a least-square problem which is solved analytically under the additional requirement of closeness of the multiplicative factors to unity. Application of the 2D-dynamic approach to the case of Flight Delta 191 shows that, with reference to the last 180 sec before impact, the values of the multiplicative factors were 1.09, 0.84, and 0.89; this implies that the actual values of the thrust, drag, and lift were 9% above, 16% below, and 11% below their respective nominal values. For the last 60 sec before impact, the aircraft was subject to severe windshear, characterized by a horizontal wind velocity difference of 123 fps and a vertical wind velocity difference of 80 fps. The 2D-dynamic approach is applicable to the analysis of windshear accidents in take-off or landing, especially for the case of older-generation, shorter-range aircraft which do not carry the extensive instrumentation of newer-generation, longer-range aircraft. The same methodology can be extended to the investigation of aircraft accidents originating from causes other than windshear (e.g., icing, incorrect flap position, engine malfunction), above all if its precision is further increased by combining the 2D-dynamic approach and the 2D-kinematic approach.     

An approach based on level set method for void identification of continuum structure with time-domain dynamic response

An approach based on the level set method has been developed to identify the position and geometry of voids in continuum structure using time-domain dynamic response. The level set method is employed in the proposed approach to represent the boundary of the voids implicitly. The voids are identified by solving an optimization problem which minimizes an objective function about the displacement error. The boundary of the voids is evolved by updating the level set function. The shape derivative of the objective function for the time-domain dynamic response is derived and used to construct the velocity field. Then, the level set function is updated through the velocity field. The proposed approach has been applied to several numerical examples of void identification in continuum structure. The results indicate that the proposed approach based on the level set method can identify voids effectively and accurately with time-domain dynamic response. Moreover, the effects of measure points, excitation force, noise, void distribution, numerical error, element size and boundary conditions on the approach are studied. Meanwhile, the computational costs of some examples are provided.     

A DAE formulation for geared rotor dynamics including frictional contact between the teeth

A DAE approach is presented for geared rotor dynamics simulations with rigid helical evolvent gears. It includes the normal contact force between the teeth as well as tangential components. Given the evolvent tooth flank geometry of gear 1 and gear 2 [1], the contact line and the velocity difference in the contact are found. The requirement of no penetration of the teeth yields a non-holonomic constraint and the contact normal force. The friction caused force and moment are obtained by applying Coulomb's friction model. This approach is used to investigate the dynamics of two ideal rotors with translational DoFs, which are connected by gears to one another. The driving rotor has a given angular speed, while the driven rotates unrestrainedly and is connected to a rotational damper. Because of the periodic friction terms, the solution is periodic. A direct time integration or a harmonic approach can be used for the numerical computation. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On numerical implementations of a new iterative method with boundary condition splitting for solving the nonstationary stokes problem in a strip with periodicity condition

Based on finite-difference approximations in time and a bilinear finite-element approximation in spatial variables, numerical implementations of a new iterative method with boundary condition splitting are constructed for solving the Dirichlet initial-boundary value problem for the nonstationary Stokes system. The problem is considered in a strip with a periodicity condition along it. At each iteration step of the method, the original problem splits into two much simpler boundary value problems that can be stably numerically approximated. As a result, this approach can be used to construct new effective and stable numerical methods for solving the nonstationary Stokes problem. The velocity and pressure are approximated by identical bilinear finite elements, and there is no need to satisfy the well-known difficult-to-verify Ladyzhenskaya-Brezzi-Babuska condition, as is usually required when the problem is discretized as a whole. Numerical iterative methods are constructed that are first- and second-order accurate in time and second-order accurate in space in the max norm for both velocity and pressure. The numerical methods have fairly high convergence rates corresponding to those of the original iterative method at the differential level (the error decreases approximately 7 times per iteration step). Numerical results are presented that illustrate the capabilities of the methods developed.     

The lift on a small sphere in a linear shear flow near the interface of two immiscible fluids

Analytical expressions have been derived which predict, to lowest order, the inertial lift and the lateral migration velocity of a rigid sphere translating and rotating in a linear shear flow field near the flat interface of two immiscible fluids. This asymptotic analysis is primarily based on the assumption that the two Reynolds numbers defined by the gap width between the interface and the sphere, the shear rate and the translational slip velocity with which the spherical particle moves parallel to the interface are small. Furthermore, the radius of the sphere is assumed to be small compared to the gap width. To leading order in this creeping flow regime, the linear Stokes equations are obtained and a symmetry argument can be used to show that the Stokes solution does not predict any lift force. The transverse force experienced by the sphere and its migration velocity are due to the small but finite inertial terms in the Navier-Stokes equations, which can be studied by perturbation techniques. By applying a Green's function approach and matched asymptotic methods, which also incorporate the effects of the outer Oseen-like flow regime, the three components comprising the lift velocity have been calculated in closed form: the one induced by the shear rate only, the purely slip induced one and the one due to the interaction of the slip velocity with the shear flow field. The thus obtained expressions for the case of two immiscible fluids with arbitrary density and viscosity ratios extend the results that already exist in the literature for other flow configurations, such as an unbounded shear flow field [1] or a wall-bounded one, where the wall lies either within the leading order Stokes region [2] or in the outer Oseen region [3]. (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Robust aircraft conflict resolution under trajectory prediction uncertainty

In air traffic control, aircraft velocity may be perturbed due to weather effects or measurements errors and affect trajectory prediction. We address the aircraft conflict resolution problem in the presence of such perturbations. We consider polyhedral uncertainty sets on aircraft velocity, develop a budgeted robust optimization approach and show that it is amenable to mixed-integer optimization. Numerical experiments reveal that perturbations of the order of 5% on aircraft velocities can be accounted for without significantly impacting performance.     

圆管Poiseuille流动中平均速度的一种修正剖面及其稳定性研究

本文给出了圆管Poiseuille流动中Hagen-Poiseuille速度剖面的一种修正剖面.这种剖面可看作是轴对称扰动各谐波分量非线性相互作用对平均流影响的一般体现.通过对这种速度剖面的稳定性研究,本文首次得到轴对称扰动造成失稳的结果,提出了Hagen-Poiseuille流动一种新的产生失稳的可能途径.     

Numerical simulation of discharged waste heat and contaminants into the South Estuary of the Yangtze River

In this paper, a two-dimensional depth-averaged mathematical model based on the finite volume approach is formed, which can be used to compute the flow fields and contaminant movements driven by tidal flows in estuaries with variations in the river bed. The paper also presents a set of appropriate outflow boundary conditions for immersed outlets in water, when these outlets are important to the entire flow region. As an example, the distributions and variations of velocity, temperature, and concentration caused by discharging contaminants from seven sources, two of which are submerged, near the south estuary of the Yangtze river have been simulated over a full tidal cycle. The results of the simulation and prediction are presented, and the effects of the outflow boundary conditions are also discussed.     

Mathematical tools of the kinetic theory of active particles with some reasoning on the modelling progression and heterogeneity

The mathematical approach proposed in this paper refers to the modelling, and related mathematical problems, of large systems of interacting entities whose microscopic state includes not only geometrical and mechanical variables (typically position and velocity), but also peculiar functions or specific activities. The number of the above entities is sufficiently large for describing the overall state of the system using a suitable probability distribution over the microscopic state. The first part of the paper is devoted to the derivation of suitable mathematical structures which can be properly used to model a variety of models in different fields of applied sciences. Then some research perspectives are analyzed, focussed on applications to biological systems.     

二维气体动力学中压差方程的特征分解和简单波

利用直接的方法讨论了在自相似平面上气体动力学中二维压差方程的特征分解理论,得到了压强P和特征值A±的特征分解.进一步地,若流动来自常状态,还可得到速度（u,v）的特征分解.由此,可以得到与常状态流动相邻的流动是简单波,并说明了简单波的流动区域是被一族直线所覆盖,且沿着每条直线, （P,u,v）为常数.     

A multi-mesh finite element method for Lagrange elements of arbitrary degree

We consider within a finite element approach the usage of different adaptively refined meshes for different variables in systems of nonlinear, time-depended PDEs. To resolve different solution behaviors of these variables, the meshes can be independently adapted. The resulting linear systems are usually much smaller, when compared to the usage of a single mesh, and the overall computational runtime can be more than halved in such cases. Our multi-mesh method works for Lagrange finite elements of arbitrary degree and is independent of the spatial dimension. The approach is well defined, and can be implemented in existing adaptive finite element codes with minimal effort. We show computational examples in 2D and 3D ranging from dendritic growth to solid–solid phase-transitions. A further application comes from fluid dynamics where we demonstrate the applicability of the approach for solving the incompressible Navier–Stokes equations with Lagrange finite elements of the same order for velocity and pressure. The approach thus provides an easy way to implement alternative to stabilized finite element schemes, if Lagrange finite elements of the same order are required.     

Vortex Perturbation Dynamics

An initial value approach is used to examine the dynamics of perturbations introduced into a vortex under strain. Both the basic vortex considered and the perturbations are taken as three-dimensional. An explicit solution for the time evolution of the vorticity perturbations is given for arbitrary initial vorticity. Analytical solutions for the resulting velocity components are found when the initial vorticity is assumed to be localized. For more general initial vorticity distributions, the velocity components are determined numerically. It is found that the variation in the radial direction of the initial vorticity disturbance is the most important factor influencing the qualitative behavior of the solutions. Transient growth in the magnitude of the velocity components is found to be directly attributable to the compactness of the initial vorticity.