Concentrators of deformations and fracture of plastic bodies

An approach to the investigation of shape discontinuity regions as strain concentrators is proposed. The near-concentrator strain fields are determined on the basis of the theory of ideal rigid-plastic body; under the condition of plane deformation, their determination is reduced to integration of ordinary differential equations. The deformation as a function of the location of the plastic region and its shape evolution in the process of plastic flow is studied. The plastic flow is demonstrated to be not unique (within the framework of solution completeness). A deformation criterion for the choice of the preferred plastic flow is suggested. The fracture of a V-notched strip is considered. On the basis of the solutions obtained, an approach to the investigation of the fracture processes for more complicated models is formulated.     

Flow sharing and bankruptcy games

We discuss the relationship between bankruptcy problems and flow sharing problems, respectively, and show that the latter type of problem can be interpreted as a generalization of the former. The corresponding flow sharing game is convex, hence we can use the converse reduced game property to characterize the nucleolus and the constrained egalitarian solution. Our main contribution is thus to introduce the nucleolus as an alternative solution concept for flow sharing problems, and to offer a game-theoretic interpretation of the traditional egalitarian solution.     

POD and CVT Galerkin reduced-order modeling of the flow around a cylinder

Reduced-order models are applied to the laminar vortex-shedding flow around a circular cylinder. The models rely on approximations of the solutions in the space spanned by a set of POD or CVT reduced basis vectors, which are computed from snapshots of a numerical solution. To enable the computation of drag and lift forces, the modeling of the velocity and pressure is realized via a one-way coupling. A comparison of the results of the CVT and POD reduced-order models is presented. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

带有线性外力场的双曲平均曲率流Cauchy问题经典解的生命跨度

本文提出并研究带有线性外力场的双曲平均曲率流,通过凸曲线的支撑函数,导出一个双曲型Monge-Ampère 方程并将其转化成Riemann 不变量满足的拟线性双曲方程组。利用拟线性双曲方程组Cauchy 问题的局部解理论,讨论带有线性外力场的双曲平均曲率流Cauchy 问题经典解的生命跨度（即局部解存在的最大时间区间）。     

A class of methods for linear programming

A class of methods is presented for solving standard linear programming problems. Like the simplex method, these methods move from one feasible solution to another at each iteration, improving the objective function as they go. Each such feasible solution is also associated with a basis. However, this feasible solution need not be an extreme point and the basic solution corresponding to the associated basis need not be feasible. Nevertheless, an optimal solution, if one exists, is found in a finite number of iterations (under nondegeneracy). An important example of a method in the class is the reduced gradient method with a slight modification regarding selection of the entering variable.     

Elastoplastic torsion of a cylindrical rod for finite deformations

A solution of the problem of the torsion of a cylindrical rod was obtained in /1/ for a general, isotropic, incompressible elastic material. The present paper gives an analytical solution of the elastoplastic torsion problem for finite deformations, written in terms of quadratures of elliptic functions. The non-linear kinematics of elastoplastic deformation is introduced into the defining equations with the help of a multiplicative decomposition of the deformation gradient into elastic and plastic components /2, 3/. The elastic deformation and rate of plastic deformation are related to the state of stress of the body, in accordance with the defining Mooney-Rivlin equations /4/ and the law of flow for finite deformations associated with the Tresca yield condition /5/. A non-linear first-order partial differential equation and the initial data at the elastoplastic boundary are obtained in order to determine the angle of rotation within the plastic zone of the basis formed from the eigenvectors of the stress tensor, relative to the radial direction. The integration of the resulting equation is reduced to determining the general integral of the Ricatti equation with right-hand side determined from the angular velocity of flow of the material within the plastic zone. It is shown that neglecting the finiteness of the deformation leads to too high an estimate of the rigidity of the rod.     

Three-dimensional flow over a stretching surface in a viscoelastic fluid

This article looks at the hydrodynamic elastico-viscous fluid over a stretching surface. The equations governing the flow are reduced to ordinary differential equations, which are analytically solved by applying an efficient technique namely the homotopy analysis method (HAM). The solutions for the velocity components are computed. The numerical values of wall skin friction coefficients are also tabulated. The present HAM solution is compared with the known exact solution for the two-dimensional flow and an excellent agreement is found.     

Creating Advanced Bases For Large Scale Linear Programs Exploiting Embedded Network Structure

In this paper, we investigate how an embedded pure network structure arising in many linear programming (LP) problems can be exploited to create improved sparse simplex solution algorithms. The original coefficient matrix is partitioned into network and non-network parts. For this partitioning, a decomposition technique can be applied. The embedded network flow problem can be solved to optimality using a fast network flow algorithm. We investigate two alternative decompositions namely, Lagrangean and Benders. In the Lagrangean approach, the optimal solution of a network flow problem and in Benders the combined solution of the master and the subproblem are used to compute good (near optimal and near feasible) solutions for a given LP problem. In both cases, we terminate the decomposition algorithms after a preset number of passes and active variables identified by this procedure are then used to create an advanced basis for the original LP problem. We present comparisons with unit basis and a well established crash procedure. We find that the computational results of applying these techniques to a selection of Netlib models are promising enough to encourage further research in this area.     

Application of the multigrid approach for solving 3D Navier-Stokes equations on hexahedral grids using the discontinuous Galerkin method

The Galerkin method with discontinuous basis functions is adapted for solving the Euler and Navier-Stokes equations on unstructured hexahedral grids. A hybrid multigrid algorithm involving the finite element and grid stages is used as an iterative solution method. Numerical results of calculating the sphere inviscid flow, viscous flow in a bent pipe, and turbulent flow past a wing are presented. The numerical results and the computational cost are compared with those obtained using the finite volume method.     

Solution of problems of oblique penetration of axisymmetric projectiles into soft soil based on local interaction models

A comparative analysis of the solutions of the three-dimensional problem of the oblique penetration of a rigid body into soft soil is carried out arsing interaction models based on one-dimensional solutions of the problem of the spherical cavity expansion. Both the well-known self-similar analytical solutions for an incompressible medium as well as the generalized solution for a compressible elastoplastic medium with separation of the shock wave which arises are considered. Use of the incompressible medium hypothesis, disregarding flow separation, in estimating the maximum values of the resistive forces leads to large errors. Taking account of compressibility enables the resistive forces to be refined appreciably and enables a satisfactory estimate of the deviation of the trajectories of bodies from the initial direction of motion to the obtained. In the proposed method of solving oblique penetration problems, a three-dimensional problem is reduced, on the basis of the plane sections hypothesis and disregarding peripheral mass and momentum flows, to the combined solution of a number of axisymmetric problems for each meridional section. It is shown that, with well-known local interaction models, this approach enables the reliability of the calculation of both the force and the kinematic characteristics of the penetration process to be increased considerably due to the fact that the dynamics of the free surface and cavitation effects of the covitating flow are taken into account.     

Model order reduction by projection applied to the universal Reynolds equation

The approach presented in this paper yields a reduced order solution to the universal Reynolds equation for incompressible fluids, which is valid in lubrication as well as in cavitation regions, applied to oil-film lubricated journal bearings in internal combustion engines. The extent of cavitation region poses a free boundary condition to the problem and is determined by an iterative spatial evaluation of a superposed modal solution. Using a Condensed Galerkin and Petrov–Galerkin method, the number of degrees of freedom of the original grid is reduced to obtain a fast but still accurate short-term prediction of the solution. Based on the assumption that a detailed solution of a previous combustion cycle is available, a basis and an optimal test space for the Galerkin method is generated. The resulting reduced order model is efficiently exploited in a time-saving evaluation of the Jacobian matrix describing the elastohydrodynamic coupling in a multi-body dynamics simulation using flexible components. Finally, numerical results are presented for a single crankshaft main bearing of typical dimensions.     

A mathematical model of Helmholtz type for a parachute profile in the presence of gravity

A model of Helmholtz type for a plane inviscid incompressible and potential fluid flow past a curvilinear obstacle of parachute in the presence of gravity is considered. Assuming that the “attack” (wind) flow is unsteady, it is shown that a bounded cavity zone should occur behind the obstacle. The determination of the fluid flow is reduced to a boundary value problem of Volterra type, for a half plane whose solution is explicitly set up, once the unknown separation (jet) lines are found under some approximation hypotheses.     

A problem in the theory of mixing layers

The classical problem of the free steady mixing layer which is formed as the result of the interaction between two parallel homogeneous flows which move with different velocities and come into contact in a certain section is considered. Subject to the additional condition that the first derivative of the solution in a class of self-similar functions is positive, a boundary-value problem is studied, for values of the self-similarity index m > 0, which describes the mixing of two viscous streams of the same fluid for m = 1 [1] and for m = 2 [2]. The method of investigation used [3–5] enables the third-order non-linear equation to be reduced to a first-order equation and enables the corresponding solutions (Gz) to be constructed in a parametric form as a function of the values of m. A knowledge of the behaviour of the velocity profile of the main stream can be used to investigate the flow stability. The results obtained form the basis of the subsequent construction of the solution of Lock's problem [6] and the investigation of the uniqueness of the solutions obtained.     

幂律流体在分形介质中不定常椭圆渗流

以椭圆渗流模型为基础,得到了分形油藏中,垂直裂缝井的幂律型非牛顿流体的定常渗流的压力分布公式和产量公式;并用数值差分的方法研究了分形油藏中,垂直裂缝井的不定常渗流,分析了分形参数对不定常压力的影响;同时从平均质量守恒方程出发,得到了相应的不定常渗流的近似解析解。分析表明:用椭圆渗流模型研究垂直裂缝井的渗流,大大简化了渗流问题的复杂性。     

A capable neural network model for solving the maximum flow problem

This paper presents an optimization technique for solving a maximum flow problem arising in widespread applications in a variety of settings. On the basis of the Karush–Kuhn–Tucker (KKT) optimality conditions, a neural network model is constructed. The equilibrium point of the proposed neural network is then proved to be equivalent to the optimal solution of the original problem. It is also shown that the proposed neural network model is stable in the sense of Lyapunov and it is globally convergent to an exact optimal solution of the maximum flow problem. Several illustrative examples are provided to show the feasibility and the efficiency of the proposed method in this paper.     

Blasius flow of thixotropic fluids: A numerical study

In the present work, laminar, two-dimensional flow of an incompressible thixotropic fluid obeying Harris rheological model is investigated above a fixed semi-infinite plate- the so-called Blasius flow. Assuming that the flow is occurring at high Reynolds number, use will be made of the boundary layer theory to simplify the equations of motion. The equations so obtained are then reduced to a single fourth-order ODE using a suitable similarity variable. It is shown that Harris fluids do not render themselves to a self-similar solution in Blasius flow. A local similarity solution is found which enabled investigating the effects of the model parameters on the velocity profile and wall shear stress at a given location above the plate. Numerical results show that for the Harris model to represent thixotropic fluids, the sign and magnitude of the material parameters appearing in this fluid model cannot be arbitrary.     

A free boundary problem of elliptic equation arising in supersonic flow past a conical body

We study free boundary value problems of elliptic equation caused by a supersonic flow past a non-symmetric conical body. The flow is described by the potential flow equation. In the self-similar coordinate system the problem can be reduced to a boundary value problem of second order nonlinear elliptic equation with a free boundary. Applying the partial hodograph transformation and the method of nonlinear alternative iteration we proved the existence of solution to this boundary value problem. Consequently, we also proved the conclusion that for the problem of supersonic flow past a conical body, if the conical body is slightly different from a circular cone with its vertex angle less than a given value determined by the parameters of the coming flow, then there exists a weak entropy solution with an attached conical shock.     

Nonisothermal gas flow in a long channel analyzed on the basis of the kinetic <Emphasis Type="Italic">S</Emphasis>-Model

The nonisothermal steady rarefied gas flow driven by a given pressure gradient (Poiseuille flow) or a temperature gradient (thermal creep) in a long channel (pipe) of an arbitrary cross section is studied on the basis of the linearized kinetic S-model. The solution is constructed using a high-order accurate conservative method. The numerical computations are performed for a circular pipe and for a cross section in the form of a regular polygon inscribed in a circle. The basic characteristic of interest is the gas flow rate through the channel. The solutions are compared with previously known results. The flow rates computed for various cross sections are also compared with the corresponding results for a circular pipe.     

Indirect boundary integral formulation for the solution of shear‐thinning flow inside single rotor mixers

An efficient indirect boundary integral formulation for the evaluation of inelastic non‐Newtonian shear‐thinning flows at low Reynolds number is presented in this article. The formulation is based on the solution of a homogeneous Stokes flow field and the use of a particular solution for the nonlinear non‐Newtonian terms that yields the complete solution to the problem. Matrix multiplications are reduced in comparison to other means of handling nonlinear terms in boundary integral formulations such as the dual reciprocity method. The iterative solution of the nonlinear system of equations has been performed with a modified Newton‐Raphson method obtaining accurate results for values of the power law index as low as 0.4 without domain partitioning. Geometries such as Couette flow and a typical industrial polymer mixer have been analyzed with the proposed method obtaining good results with a reduction in computational cost compared with other equivalent formulations. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27:1610–1627, 2011     

Plane potential flow past a cylinder with porous surface

In the paper two models for the potential flow past a cylinder with porous or perforated surface are reduced to nonlinear Riemann-Hilbert problems for the complex velocity function and general existence and uniqueness theorems for the solution of the problems are derived. Further the problems are reduced to nonlinear singular integral equations of Hilbert type whose solvability is investigated by means of the Schauder and Banach fixed point theorem.