Computational homogenization of materials with small deformation to determine configurational forces

In this work the mechanical boundary condition for the micro problem in a two-scaled homogenization using a FE² approach is discussed. The strain tensor is often used in the literature for small deformation problem to determine the boundary conditions for the boundary value problem on the micro level. This strain tensor based boundary condition gives consistent homogenized mechanical quantities, e.g. stress tensor and elasticity tensor, but the present work points out that it leads to unphysical homogenized configurational forces. Instead, we propose a displacement gradient based boundary condition for the micro problem. Results show that the displacement gradient based boundary condition can give not only the consistent homogenized mechanical quantities but also the appropriate homogenized configurational forces. The interpretation of the displacement gradient based boundary condition is discussed. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On Distributed $H^1$ Shape Gradient Flows in Optimal Shape Design of Stokes Flows: Convergence Analysis and Numerical Applications

We consider optimal shape design in Stokes flow using $H^1$ shape gradient flows based on the distributed Eulerian derivatives. MINI element is used for discretizations of Stokes equation and Galerkin finite element is used for discretizations of distributed and boundary $H^1$ shape gradient flows. Convergence analysis with a priori error estimates is provided under general and different regularity assumptions. We investigate the performances of shape gradient descent algorithms for energy dissipation minimization and obstacle flow. Numerical comparisons in 2D and 3D show that the distributed $H^1$ shape gradient flow is more accurate than the popular boundary type. The corresponding distributed shape gradient algorithm is more effective.     

A STUDY ON GRADIENT BLOW UP FOR VISCOSITY SOLUTIONS OF FULLY NONLINEAR, UNIFORMLY ELLIPTIC EQUATIONS

We investigate sharp conditions for boundary and interior gradient estimates of continuous viscosity solutions to fully nonlinear,uniformly elliptic equations under Dirichlet boundary conditions. When ...     

A gradient bound for free boundary graphs

We prove an analogue for a one‐phase free boundary problem of the classical gradient bound for solutions to the minimal surface equation. It follows, in particular, that every energy‐minimizing free boundary that is a graph is also smooth. The method we use also leads to a new proof of the classical minimal surface gradient bound. © 2010 Wiley Periodicals, Inc.     

基于简化的应变梯度理论下Kirchhoff板模型边值问题的提法及其应用北大核心CSCD

考虑应变梯度和速度梯度的影响,建立薄板控制微分方程及给出其边值问题的提法,修正了前人给出的薄板角点条件.采用Levy法,给出受分布力作用下简支板的挠度及自由振动频率的解析解.通过与文献中分子动力学数据对比,验证了该文模型的有效性并提出校核材料参数的一种方法.研究结果表明,增大弹性地基和应变梯度参数可以有效提高板的等效刚度,而速度梯度参数则相反.该文提出的板的边值问题为研究薄板在复杂支撑边界及外荷载等条件提供了理论依据.同时,有望为其有限元法、有限差分法和基于能量原理的Galerkin法等数值方法提供理论依据.     

Compact gradient tracking in shape optimization

In the present paper we consider the minimization of gradient tracking functionals defined on a compact and fixed subdomain of the domain of interest. The underlying state is assumed to satisfy a Poisson equation with Dirichlet boundary conditions. We proof that, in contrast to the situation of gradient tracking on the whole domain, the shape Hessian is not strictly H ^1/2-coercive at the optimal domain which implies ill-posedness of the shape problem under consideration. Shape functional and gradient require only knowledge of the Cauchy data of the state and its adjoint on the boundaries of the domain and the subdomain. These data can be computed by means of boundary integral equations when reformulating the underlying differential equations as transmission problems. Thanks to fast boundary element techniques, we derive an efficient algorithm to solve the problem under consideration.     

Solutions of thermal boundary layer equations when temperature gradient at the moving flat plate in parabolic flow is prescribed

Solutions of steady as well as unsteady three-dimensional incompressible thermal boundary layer equations are studied when the temperature gradient at the moving flat plate in parabolic flow is prescribed. A general analysis is made and different cases are studied by giving values to β and Cx which determine the gradient and curvature of the outer flow steam lines. The components of velocity in boundary layer are discussed by Sarma and Gupta and those results are used to analyse the thermal boundary layer equations. The response of temperature of the plate are studied for large and small times and curves are drawn representing the variation of temperature with time for various cases. The limiting time is also calculated.     

Control of phase boundary evolution in metal solidification for new thermodynamic parameters of the metal

The problem of controlling the phase boundary evolution in the course of solidification of metals with different thermodynamic properties is studied. The underlying mathematical model of the process is based on a three-dimensional nonstationary two-phase initial–boundary value problem of the Stefan type. The control functions are determined by optimal control problems, which are solved numerically with the help of gradient optimization methods. The gradient of the cost function is exactly computed by applying the fast automatic differentiation technique. The research results are described and analyzed. Some of them are illustrated.     

About the projection of kinematical quantities within the FE2-approach

Within the FE²-approach strain–like quantities are projected to the boundary of the attached microstructure yielding Dirichlet boundary conditions for the microscopic boundary value problem. By use of the divergence theorem one can show that the projection of gradient information to the boundary of the attached microstructure fulfills the condition that the volumetrical average of the strain–like quantities in the inner of the microstructure is exactly the same as the projected quantity, provided that the microstructure is continous. However, within the Cosserat continuum theory [1] the strain tensor depends directly on the extra rotation and not on the gradient of it. Thus, the projection can not be carried out in the same way as for gradient information. The present approach shows how to circumvent this problem by giving the extra–rotation within the Cosserat continuum theory a strong kinematical interpretation in the form of a displacement field, which shows the same effects as the extra–rotation itself. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Buckling analyses of three characteristic-lengths featured size-dependent gradient-beam with variational consistent higher order boundary conditions

In this work, a three characteristic-lengths featured size-dependent gradient-beam is constructed by adopting the modified nonlocal model, resulting in much more general constitutive equation with stress gradient up to four-order and strain gradient to two-order. The six-order differential governing equation for transverse displacement is formulated. All boundary conditions especially variational consistent higher order boundary conditions of the present model are derived with the aid of weighted residual approach. The closed-form solutions to critical buckling loads under different sets of boundary conditions are systematically formulated with higher order boundary conditions incorporated. The numerical results show that both nonlocal parameters have significant effect on the buckling behaviors. Meanwhile, if two nonlocal parameters are taken as same, the present results cannot always reduce to that from Eringen's nonlocal model. Due to its clear physical meaning, the present model is expected to be widely adopted in mechanical analyses of nano-structures.     

Asymptotic and numerical solutions of three‐dimensional boundary‐layer flow past a moving wedge

We consider a laminar boundary‐layer flow of a viscous and incompressible fluid past a moving wedge in which the wedge is moving either in the direction of the mainstream flow or opposite to it. The mainstream flows outside the boundary layer are approximated by a power of the distance from the leading boundary layer. The variable pressure gradient is imposed on the boundary layer so that the system admits similarity solutions. The model is described using 3‐dimensional boundary‐layer equations that contains 2 physical parameters: pressure gradient (β) and shear‐to‐strain‐rate ratio parameter (α). Two methods are used: a linear asymptotic analysis in the neighborhood of the edge of the boundary layer and the Keller‐box numerical method for the full nonlinear system. The results show that the flow field is divided into near‐field region (mainly dominated by viscous forces) and far‐field region (mainstream flows); the velocity profiles form through an interaction between 2 regions. Also, all simulations show that the subsequent dynamics involving overshoot and undershoot of the solutions for varying parameter characterizing 3‐dimensional flows. The pressure gradient (favorable) has a tendency of decreasing the boundary‐layer thickness in which the velocity profiles are benign. The wall shear stresses increase unboundedly for increasing α when the wedge is moving in the x‐direction, while the case is different when it is moving in the y‐direction. Further, both analysis show that 3‐dimensional boundary‐layer solutions exist in the range −1<α<∞. These are some interesting results linked to an important class of boundary‐layer flows.     

Micromechanical Simulation of the Hall-Petch Effect with a Crystal Gradient Theory including a Grain Boundary Yield Criterion

A viscoplastic strain gradient crystal plasticity theory based on the gradient of the equivalent plastic strain ∇γ_eq is proposed. A grain boundary yield condition is introduced. The microstructural explanation of the Hall-Petch effect, accounting for notch-like stress concentrations at the grain boundary as a result of discrete slip bands, is reviewed. Periodic tensile test FEM simulation results illustrate the prediction of the numerical model. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

一个新的边界层方程及其在形状控制中的应用

本文给出固壁边界上(即一个二维流形上) 的流体速度梯度和压力的二阶偏微分方程, 从而也给出边界上法向应力, 以及流体中运动物体所受的阻力和升力的计算公式. 本方法的创新在于边界上法向速度梯度不是通过在边界层内速度梯度的数值微分达到, 而是通过它与其他变量一起作为一组偏微分方程的解而得到, 证明边界层方程组的适定性问题, 并且给出解关于边界形状的Gâteaux 导数所满足的偏微分方程. 本文将本方法应用于飞机外形的形状最优控制, 给出阻力泛函关于形状第一变分的可计算形式. 数值例子表明, 用本方法得到的阻力精度比通用程序得到要高.     

紧致黎曼流形上一类非线性热方程的梯度估计

In this paper,we study gradient estimates for the nonlinear heat equation ut-△u =au log u,on compact Riemannian manifold with or without boundary.We get a Hamilton type gradient estimate for the positi...     

Parabolic variational inequality with parameter and gradient constraints

This paper concerns a singular control problem whose value function is governed by a time-dependent HJB equation with gradient constraints. The method is to transform a two-dimensional parabolic variational inequality with gradient constraints into a double obstacle problem with parameter involving two free boundaries that correspond to the investment and disinvestment policies. Moreover we analyze the behaviors of the free boundary surfaces. The main difficulties are to show the free boundary surfaces to be smooth with respect to time and to find the properties of free boundaries with respect to parameter.     

考虑二次梯度项影响的非线性不稳定渗流问题的精确解

考虑了二次梯度项影响的非线性径向流动问题的无限大地层和有界地层渗流模型．在井底定流量和定压生产时,对无限大地层及有界地层(包括封闭和定压地层)六种情况,利用广义Weber变换和广义Hankel变换求得了实空间的解析解,分析了非线性压力解与线性压力解的差异,发现在晚时段其差异可达8%以上．因此在试井长时要考虑二次梯度项的影响．     

A singular nonlinear elliptic equation with natural growth in the gradient

This paper is devoted to the existence and regularity of the homogenous Dirichlet boundary value problem for a singular nonlinear elliptic equation with natural growth in the gradient. By certain transformations, the problem can be transformed formally into either a Dirichlet problem or boundary blowup problems without gradient term, for which the corresponding existence results are also derived, which is a partial extension and supplement to the previous works. Copyright © 2008 John Wiley & Sons, Ltd.     

Well-posedness of the free boundary problem in compressible elastodynamics

We study the free boundary problem for the flow of a compressible isentropic inviscid elastic fluid. At the free boundary moving with the velocity of the fluid particles the columns of the deformation gradient are tangent to the boundary and the pressure vanishes outside the flow domain. We prove the local-in-time existence of a unique smooth solution of the free boundary problem provided that among three columns of the deformation gradient there are two which are non-collinear vectors at each point of the initial free boundary. If this non-collinearity condition fails, the local-in-time existence is proved under the classical Rayleigh–Taylor sign condition satisfied at the first moment. By constructing an Hadamard-type ill-posedness example for the frozen coefficients linearized problem we show that the simultaneous failure of the non-collinearity condition and the Rayleigh–Taylor sign condition leads to Rayleigh–Taylor instability.     

Control of hydrothermal waves in a thermocapillary flow using a gradient-based control strategy

The problem of controlling the hydrothermal waves in a thermocapillary flow is addressed using a gradient based control strategy. The control is a spatially and temporally varying heat flux boundary condition at the free surface. The control which minimizes the hydrothermal waves is found using a conjugate gradient method, where the gradient of the objective function with respect to the control variables is obtained from solving a set of adjoint equations. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

High-Resolution Color Image Reconstruction with Neumann Boundary Conditions

This paper studies the application of preconditioned conjugate gradient methods in high-resolution color image reconstruction problems. The high-resolution color images are reconstructed from multiple undersampled, shifted, degraded color frames with subpixel displacements. The resulting degradation matrices are spatially variant. To capture the changes of reflectivity across color channels, the weighted H ₁ regularization functional is used in the Tikhonov regularization. The Neumann boundary condition is also employed to reduce the boundary artifacts. The preconditioners are derived by taking the cosine transform approximation of the degradation matrices. Numerical examples are given to illustrate the fast convergence of the preconditioned conjugate gradient method.