A numerical technique for solving IHCPs using Tikhonov regularization method

This study is intended to provide a numerical algorithm for solving a one-dimensional inverse heat conduction problem. The given heat conduction equation, the boundary conditions, and the initial condition are presented in a dimensionless form. The numerical approach is developed based on the use of the solution to the auxiliary problem as a basis function. To regularize the resultant ill-conditioned linear system of equations, we apply the Tikhonov regularization method to obtain the stable numerical approximation to the solution.     

Numerical solution of an equilibrium problem for an elastic body with a thin delaminated rigid inclusion

Under consideration is a 2D-problem of elasticity theory for a body with a thin rigid inclusion. It is assumed that there is a delamination crack between the rigid inclusion and the elastic matrix. At the crack faces, the boundary conditions are set in the form of inequalities providing mutual nonpenetration of the crack faces. Some numerical method is proposed for solving the problem, based on domain decomposition and the Uzawa algorithm for solving variational inequalities.We give an example of numerical calculation by the finite element method.     

Numerical simulation of equilibrium of an elastic two-layer structure with a through crack

In this paper, a problem of equilibrium of two elastic bodies pasted together along a curve is considered. It is assumed that there is a through crack on a part of the curve. Nonlinear boundary conditions providing mutual non-penetration between the crack faces are set. The main objective of the paper is to construct and test a numerical algorithm for solving the equilibrium problem. The algorithm is based on two approaches: a domain decomposition method and Uzawa method for solving variational inequalities. A numerical experiment illustrates the efficiency of the algorithm.     

Variational formulation of a modified couple stress theory and its application to a simple shear problem

A variational formulation is provided for the modified couple stress theory of Yang et al. by using the principle of minimum total potential energy. This leads to the simultaneous determination of the equilibrium equations and the boundary conditions, thereby complementing the original work of Yang et al. where the boundary conditions were not derived. Also, the displacement form of the modified couple stress theory, which is desired for solving many problems, is obtained to supplement the existing stress-based formulation. All equations/expressions are presented in tensorial forms that are coordinate-invariant. As a direct application of the newly obtained displacement form of the theory, a simple shear problem is analytically solved. The solution contains a material length scale parameter and can capture the boundary layer effect, which differs from that based on classical elasticity. The numerical results reveal that the length scale parameter (related to material microstructures) can have a significant effect on material responses.      

Variational formulation of a modified couple stress theory and its application to a simple shear problem

A variational formulation is provided for the modified couple stress theory of Yang et al. by using the principle of minimum total potential energy. This leads to the simultaneous determination of the equilibrium equations and the boundary conditions, thereby complementing the original work of Yang et al. where the boundary conditions were not derived. Also, the displacement form of the modified couple stress theory, which is desired for solving many problems, is obtained to supplement the existing stress-based formulation. All equations/expressions are presented in tensorial forms that are coordinate-invariant. As a direct application of the newly obtained displacement form of the theory, a simple shear problem is analytically solved. The solution contains a material length scale parameter and can capture the boundary layer effect, which differs from that based on classical elasticity. The numerical results reveal that the length scale parameter (related to material microstructures) can have a significant effect on material responses.     

Analytical solution to a generalized spectral problem in a method of recalculating boundary conditions for the biharmonic equation

An iterative algorithm with an efficient preconditioner for the numerical solution to an elasticity problem in the approximation of plate theory with mixed boundary conditions is proposed and substantiated. Exact constants of energy equivalence for optimization of the iterative method are obtained. Inversion of the preconditioner is equivalent to the double inversion of a discrete analog of the Laplace operator with Dirichlet boundary conditions.     

Asymptotic solution of a quasistatic thermoelasticity problem for a slender rod

An asymptotic expansion is constructed for solving a quasistatic thermo-elasticity problem for a slender cylindrical rod in the presence of mass forces and non-linear heat sources. The algorithm for constructing the asymptotic form, based on the method of boundary functions, is fairly simple and convenient for carrying out numerical calculations. A deduction is made on the basis of the asymptotic form constructed on how to select correctly a simplified one-dimensional model so as to obtain a better approximation for the solution of the initial two-dimensional problem. An existence theorem for the solution is proved under certain conditions.     

The cauchy problem for the radiative transfer equation with generalized conjugation conditions

The initial boundary problem for the nonstationary radiative transfer equation in a nonhomogeneous plane layer with generalized conjugation conditions on the material interface is studied. A generalized Monte Carlo algorithm is proposed for solving the problem, and numerical experiments are discussed.     

The use of the method of boundary states to analyse an elastic medium with cavities and inclusions

The analytical method of boundary states is developed and theoretically substantiated. A corollary of the Weierstrass theorem is proved according to which a function that is harmonic in a bounded, simply connected domain can be approximated by a series of homogeneous harmonic polynomials. A basis of the space of functions that are harmonic outside any neighbourhood of a point is constructed. An algorithm is developed for filling the basis of the space of the states of a multicavity elastic body. The method is used to solve a series of problems of determining of the stress-strain state of an unbounded elastic medium containing spherical cavities or inclusions with different boundary conditions: the boundary of the cavity is free (the Southwell problem), constrained or under conditions of contact with a rigid core. The effect of the width of the intercavity layer on the stress concentration is analysed in a non-axisymmetric problem with two cavities. The form of the relation between the mean-square discrepancy in the boundary conditions of the solution obtained and the number of elements in the basis is indicative of the numerical convergence of the solution of this problem.     

The Design of axisymmetric ducts for incompressible flow with vorticity.

In this paper a numerical algorithm is described for solving the boundary value problem associated with axisymmetric, inviscid, incompressible, rotational (and irrotational) flow in order to obtain duct wall shapes from prescribed wall velocity distributions. The governing equations are formulated in terms of the stream function and the function as independent variables where for irrotational flow can be recognized as the velocity potential function, for rotational flow ceases being the velocity potential function but does remain orthogonal to the stream lines. A numerical method based on finite differences on a uniform mesh is employed. The technique described is capable of tackling the so–called inverse problem where the velocity wall distributions are prescribed from which the duct wall shape is calculated, as well as the direct problem where the velocity distribution on the duct walls are calculated from prescribed duct wall shapes. The two different cases as outlined in this paper are in fact boundary value problems with Neumann and Dirichlet boundary conditions respectively. Even though both approaches are discussed, only numerical results for the case of the Dirichlet boundary conditions are given. A downstream condition is prescribed such that cylindrical flow, that is flow which is independent of the axial coordinate, exists. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

任意厚度具有自由边叠层板的精确解析解

自由边问题一直是三维弹性力学中的难题,通常很难满足自由边上一个正应力和两个剪应力都等于0．基于三维弹性力学基本方程和状态空间方法,引入自由边界位移函数并考虑全部弹性常数,建立了正交异性具有自由边单层和叠层板的状态方程．对状态方程中的变量以级数形式展开,通过边界条件的满足精确求解任意厚度具有自由边叠层板的位移和应力,此解满足层间应力和位移的连续条件．算例计算表明,采用引入的位移函数形式,简化了计算过程并且采用较少的级数项可以获得收敛解．与有限元方法计算结果进行了对比,可以得到较高精度的数值结果．其解可以作为其它数值方法和半解析方法的参考解．     

A level set approach for the solution of a state-constrained optimal control problem

Summary. State constrained optimal control problems for linear elliptic partial differential equations are considered. The corresponding first order optimality conditions in primal-dual form are analyzed and linked to a free boundary problem resulting in a novel algorithmic approach with the boundary (interface) between the active and inactive sets as optimization variable. The new algorithm is based on the level set methodology. The speed function involved in the level set equation for propagating the interface is computed by utilizing techniques from shape optimization. Encouraging numerical results attained by the new algorithm are reported on.Mathematics Subject Classification (1991): 35R35, 49K20, 49Q10, 65K10Revised version received March 19, 2003     

Numerical investigation of a model problem for the Poisson equation with inequality constraints in a domain with a cut

A model problem is considered for the Poisson equation in a two-dimensional domain with a cut. The Dirichlet and Neumann conditions are imposed on the exterior boundary of the domain together with the nonnegativity condition for the jump across the edges of the cut. In addition, the absolute value of the gradient inside the domain must be bounded by some constant. The boundary value problem turns into a variational problem, and the unknown function must yield the minimum of the energy functional on some convex set. After discretization of the problem by the finite element method, an Uzawa-type algorithm is used to find a solution. Some examples are included of solving the discrete problem.     

Mathematical and Numerical Simulation of Equilibrium of an Elastic Body Reinforced by a Thin Elastic Inclusion

A boundary value problem describing the equilibrium of a two-dimensional linear elastic body with a thin rectilinear elastic inclusion and possible delamination is considered. The stress and strain state of the inclusion is described using the equations of the Euler–Bernoulli beam theory. Delamination means the existence of a crack between the inclusion and the elastic matrix. Nonlinear boundary conditions preventing crack face interpenetration are imposed on the crack faces. As a result, problem with an unknown contact domain is obtained. The problem is solved numerically by applying an iterative algorithm based on the domain decomposition method and an Uzawa-type algorithm for solving variational inequalities. Numerical results illustrating the efficiency of the proposed algorithm are presented.     

An algorithm for maximizing a convex function over a simple set

The problem of maximizing a convex function on a so-called simple set is considered. Based on the optimality conditions [19], an algorithm for solving the problem is proposed. This numerical algorithm is shown to be convergent. The proposed algorithm has been implemented and tested on a variety of test problems.     

对称迭层矩形板的平面应力分析

常用的对称迭层板为各向异性板．根据平面应力问题的基本方程精确地用应力函数解法求得了各向异性板的一般解析解．推导出平面内应力和位移的一般公式,其中积分常数由边界条件来决定．一般解包括三角函数和双曲函数组成的解,它能满足4个边为任意边界条件的问题．还有代数多项式解,它能满足4个角的边界条件．因此一般解可用以求解任意边界条件下的平面应力问题．以4边承受均匀法向和切向载荷以及非均匀法向载荷的对称迭层方板为例,进行了计算和分析．     

New extragradient method for a class of equilibrium problems in Hilbert spaces

The paper proposes a new extragradient algorithm for solving strongly pseudomonotone equilibrium problems which satisfy a Lipschitz-type condition recently introduced by Mastroeni in auxiliary problem principle. The main novelty of the paper is that the algorithm generates the strongly convergent sequences in Hilbert spaces without the prior knowledge of Lipschitz-type constants and any hybrid method. Several numerical experiments on a test problem are also presented to illustrate the convergence of the algorithm.     

Saturation-free numerical scheme for computing the flow past a lattice of airfoils and the determination of separation points in a viscous fluid

A numerical method for computing the potential flow past a lattice of airfoils is described. The problem is reduced to a linear integrodifferential equation on the lattice contour, which is then approximated by a linear system of equations with the help of specially derived quadrature formulas. The quadrature formulas exhibit exponential convergence in the number of points on an airfoil and have a simple analytical form. Due to its fast convergence and high accuracy, the method can be used to directly optimize the airfoils as based on any given integral characteristics. The shear stress distribution and the separation points are determined from the velocity distribution at the airfoil boundary calculated by solving the boundary layer equations. The method proposed is free of laborious grid generation procedures and does not involve difficulties associated with numerical viscosity at high Reynolds numbers.     

A numerical solution to nonlinear multi‐point boundary value problems in the reproducing kernel space

In this paper, a new numerical algorithm is provided to solve nonlinear multi‐point boundary value problems in a very favorable reproducing kernel space, which satisfies all complex boundary conditions. Its reproducing kernel function is discussed in detail. The theorem proves that the approximate solution and its first‐ and second‐order derivatives all converge uniformly. The numerical experiments show that the algorithm is quite accurate and efficient for solving nonlinear multi‐point boundary value problems. Copyright © 2010 John Wiley & Sons, Ltd.