Numerical solution for curved crack problem in elastic half-plane using hypersingular integral equation

A hypersingular integral equation for the curved crack problems of an elastic half-plane is introduced. Formulation of the equation is based on the usage of a modified complex potential. The potential is generally expressed in the form of a Cauchy-type integral. The modified complex potential is composed of the principal part and the complementary part. The principal part of the complex potential is actually equivalent to the original complex potential for the curved crack in an infinite plate. The role of the complementary part is to eliminate the boundary traction along the boundary of the half-plane caused by the principal part. From the assumed boundary traction condition, a hypersingular integral equation is obtained for the curved crack problems of an elastic half-plane. The curve length coordinate method is used to obtain a final solution. Several numerical examples are presented that prove the efficiency of the suggested method.     

Microstructure evolution during severe plastic deformation

Radiotracer diffusion studies of severely deformed, ultra-fine grained materials have revealed the presence of ultra-fast transport paths, which include “non-equilibrium” grain boundaries and free volume. Under some experimental conditions, percolating porosity is produced even in pure copper. Micro-cracks may form in metals, if the local maximum shear stress exceeds the shear yield stress. However, their growth and propagation is postponed till late in the deformation process owing to the ductility of metals, the hydrostatic component of the stress system and/or dynamic recovery/recrystallization. In other words, crack growth and propagation is present only when the scope for further deformation is highly restricted. Using this approach, the load required for equal channel angular pressing, the change in the slope of the Hall–Petch plot with decreasing grain size and the theoretical limit for the smallest grain size attainable in a metal in a severe plastic deformation process are predicted and validated by experimental results. Experimentally successful prevention of percolated crack formation by the superposition of a hydrostatic pressure is also accounted for using this model.     

CONSERVATION LAWS FOR HEATED LAMINAR RADIAL LIQUID AND FREE JETS

The conserved quantities for the heated radial liquid jet and the heated radial free jet are established by using conservation laws. The flow in a heated radial jet is described by Prandtl's momentum boundary layer equation, the continuity equation and the energy equation. Viscous dissipation is neglected. The multiplier approach is used to derive the conservation laws for the system of three equations for the velocity components and the temperature and three conserved vectors are obtained. The conservation laws for the system of two partial differential equations for the stream function formulation are also computed by the multiplier approach and three conserved vectors are obtained. One of these is a non-local conserved vector for the system. The conserved quantities for the heated radial liquid jet and the heated radial free jet, emitted into a stationary fluid of uniform temperature θ_∞, are derived by integrating the conservation laws across the jet.     

搅拌器内部流动的无网格法三维数值模拟

本文将移动粒子半隐式法(MPS)的基本算法由二维扩展至三维。将圆柱坐标系引入到初场粒子的布置中,避免了在笛卡儿坐标系下处理不规则形状(如斜边或曲边)问题时粒子初场布置困难和精确度较低的问题,改善了对计算边界条件表达的精确性。引入移动边界模型,对直叶片搅拌器的内部流动进行了三维数值模拟。还提出了一种新的初始粒子布置简易方法,明显简化粒子初始布置时的复杂程度,提高了对三维复杂几何形状问题的可操作性。     

A Solvable N-body Problem of Goldfish Type Featuring N2 Arbitrary Coupling Constants

A N-body problem “of goldfish type” is introduced, the Newtonian (“acceleration equal force”) equations of motion of which describe the motion of N pointlike unit-mass particles moving in the complex z-plane. The model—for arbitrary N—is solvable, namely its configuration (positions and velocities of the N “particles”) at any later time t can be obtained from its configuration at the initial time by algebraic operations. It features specific nonlinear velocity-dependent many-body forces depending on N² arbitrary (complex) coupling constants. Sufficient conditions on these constants are identified which cause the model to be isochronous—so that all its solutions are then periodic with a fixed period independent of the initial data. A variant with twice as many arbitrary coupling constants, or even more, is also identified.     

A Finite Volume Method Preserving Maximum Principle for the Conjugate Heat Transfer Problems with General Interface Conditions

In this paper, we present a unified finite volume method preserving discrete maximum principle (DMP) for the conjugate heat transfer problems with general interface conditions. We prove the existence of the numerical solution and the DMP-preserving property. Numerical experiments show that the nonlinear iteration numbers of the scheme in [24] increase rapidly when the interfacial coefficients decrease to zero. In contrast, the nonlinear iteration numbers of the unified scheme do not increase when the interfacial coefficients decrease to zero, which reveals that the unified scheme is more robust than the scheme in [24]. The accuracy and DMP-preserving property of the scheme are also veried in the numerical experiments.     

Inverse conductivity Problem with Internal Data

This paper concerns the reconstruction of a scalar coefficient of a second-order elliptic equation in divergence form posed on a bounded domain from internal data. This problem finds applications in multi-wave imaging, greedy methods to approximate parameter-dependent elliptic problems, and image treatment with partial differential equations. We first show that the inverse problem for smooth coefficients can be rewritten as a linear transport equation. Assuming that the coefficient is known near the boundary, we study the well-posedness of associated transport equation as well as its numerical resolution using discontinuous Galerkin method. We propose a regularized transport equation that allow us to derive rigorous convergence rates of the numerical method in terms of the order of the polynomial approximation as well as the regularization parameter. We finally provide numerical examples for the inversion assuming a lower regularity of the coefficient, and using synthetic data.     

Uniform Convergence of Multigrid V-Cycle on Adaptively Refined Finite Element Meshes for Elliptic Problems with Discontinuous Coefficients

The multigrid V-cycle methods for adaptive finite element discretizations of two-dimensional elliptic problems with discontinuous coefficients are considered. Under the conditions that the coefficient is quasi-monotone up to a constant and the meshes are locally refined by using the newest vertex bisection algorithm, some uniform convergence results are proved for the standard multigrid V-cycle algorithm with Gauss-Seidel relaxations performed only on new nodes and their immediate neighbours. The multigrid V-cycle algorithm uses $\mathcal{O}(N)$ operations per iteration and is optimal.     

半线性抛物最优控制问题全离散插值系数有限元方法的收敛性分析

本文考虑全离散插值系数有限元方法求解半线性抛物最优控制问题,其中控制变量用分片常数函数逼近,状态变量和对偶状态变量用分片线性函数逼近.对于方程中的半线性项,先用插值系数技巧处理,再用牛顿迭代法求解.通过引入一些辅助变量和投影算子,并利用有限元空间的逼近性质,得到半线性抛物最优控制问题插值系数有限元方法的收敛性结果;数值算例结果验证了理论结果的正确性.