周期性结构热动力时间-空间多尺度分析

研究一种时间-空间多尺度渐近均匀化分析方法，模拟不同的极端热和动力载荷下微尺度多 相周期性结构中热动力响应问题，并建立一个广义的波动函数场控制方程描述热动力响应. 通过引入一个放大空间尺度和两个缩小时间尺度，在不同时间尺度上获得由空间非均匀性引 起的波动效应和非局部效应. 根据高阶均匀化理论在空间和时间上进行均匀化，获得高阶非 局部函数场波动方程. 并进一步用C0连续修正了高阶非局部函数场波动方程的有限元近 似解，使问题的求解避免了对有限元离散的C1连续性要求. 并与经典的空间均匀化方法 相比较，指出了经典的空间均匀化方法的局限性，进一步以一维非傅立叶热传导和热动力问 题为例，讨论了各种情况下方法的正确性与有效性.     

Transport with a Very Low Density Contrast in Hele–Shaw Cell and Porous Medium: Evolution of the Mixing Zone

The optimal concentration of a blue dye solution with 'tracer' properties, enabling a pollutant to be marked was determined by the use of numerical, theoretical and experimental approaches. Experimental investigations were performed on a transparent Hele–Shaw cell and the concentration distribution was analyzed using an optical technique based on dye light absorption properties. The injected optimal concentration was established thanks to a theoretical and experimental study carried out on the output signal dynamics. Using the same experimental conditions, numerical simulations were performed. The very good agreement between the data (experimental and numerical) clarified that: (i) the choice of the blue dye optimal concentration was valid and (ii) the concentration-dependent density should not be neglected in flow and transport equations even if it concerns a so-called 'tracer'. Following this remark, a theoretical aspect was developed in order to determine the analogous conditions between a Hele–Shaw cell and a porous medium for the variable density transport phenomenon. The structure of the concentration-dependent dispersion tensor used in the numerical code was obtained by homogenizing the Stokes flow of a bi-component mixture. The numerical results show that, as long as the tracer density does not exceed a certain value, it is not necessary to take into account a density contrast in terms of the dispersion tensor. The classical form of the Taylor dispersion tensor can be used successfully.     

Non-gaussian diffusion modeling in composite porous media by homogenization: Tail effect

In the paper anomalous diffusion appearing in a porous medium composed of two porous components of considerably different diffusion characteristics is examined. The differences in diffusivities are supposed to result either from two medium types being present or from variations in pore size (double porosity media). The long-tail effect is predicted using the homogenization approach based on the application of multiple scale asymptotic developments. It is shown that, if the ratio of effective diffusion coefficients of two porous media is of the order of magnitude smaller or equal O( ²), where is a homogenization parameter, then the macroscopic behaviour of the composite may be affected by the presence of tail-effect. The results of the theoretical analysis were applied to a problem of diffusion in a bilaminate composite. Analytical calculations were performed to show the presence of the long-tail effect in two particular cases.Notations c _i the concentration of chemical species in water within the medium i - D _i the effective diffusion coefficient for the medium i - D _ij ^eff the macroscopic (or effective) diffusion tensor in the composite - ERV the elementary representative volume - h the thickness of the period - l a chracteristic length of the ERV or the periodic cell - L a characteristic macroscopic length - n the volumetric fraction of the material 2 - 1–n the volumetric fraction of the material 1 - N the unit vector normal to - t the time variable - x the macroscopic (or slow) space variable - y the microscopic (or fast) space variable - c _1c ,C _2c ,D _1c ,D _2c the characteristic quantities - T,T _1L ,T _2L ,T _1l ,T _2l the characteristic times - c ₁ ^* ,c ₂ ^* ,D ₁ ^* ,D ₂ ^* ,t ^* the non-dimensional variables - the homogenization parameter - ₁ the domain occupied by the material 1 - ₂ the domain occupied by the material 2 - the interface between the domains 1 and 2 - the total volume of the periodic cell - /x_i the gradient operator - the gradient operator     

Exact Results for the Homogenization of Elastic Fiber-Reinforced Solids at Finite Strain

This work is concerned with the homogenization of solids reinforced by aligned parallel continuous fibers or weakened by aligned parallel cylindrical pores and undergoing large deformations. By alternatively exploiting the nominal and material formulations of the corresponding homogenization problem and by applying the implicit function theorem, it is shown that locally homogeneous deformations can be produced in such inhomogeneous materials and form a differentiable manifold. For every macroscopic strain associated to a locally homogeneous deformation field, the effective nominal or material stress–strain relation is exactly determined and connections are also exactly established between the effective nominal and material elastic tangent moduli. These results are microstructure-independent in the sense that they hold irrespectively of the transverse geometry and distribution of the fibers or pores. A porous medium consisting of a compressible Mooney–Rivlin material with cylindrical pores is studied in detail to illustrate the general results. This work was the first time presented at the Euromech Colloqium 464 on “Fiber-reinforced Solids: Constitutive Laws and Instabilities”, September 28–October 1, 2004, Cantabria, Spain.     

HOMOGENIZATION—BASED TOPOLOGY DESIGN FOR PURE TORSION OF COMPOSITE SHAFTS

In conjunction with the homogenization theory and the finite element method, the mathematical models for designing the corss-section of composite shafts by maximizing the torsion rigidity are developed in this paper. To obtain the extremal torsion rigidity, both the cross-section of the macro scale shaft and the representative microstructure of the composite material are optimized using the new models. The micro scale computational model addresses the problem of finding the periodic microstructures with extreme shear moduli. The optimal microstructure obtained with the new model and the homogenization method can be used to improve and optimize natural or artificial materials. In order to be more practical for engineering applications, cellular materials rather than ranked materials are used in the optimal process in the existence of optimal bounds for the elastic properties. Moreover, the macro scale model is proposed to optimize the cross-section of the torsional shaft based on the tailared composites. The validating optimal results show that the models are very effective in obtaining composites with extreme elastic properties, and the cross-section of the composite shaft with the extremal torsion rigidity. The project supported by the National Natural Science Foundation of China (10172078 and 10102018)     

高速通道压裂裂缝的高导流能力分析及其影响因素研究

高速通道压裂是近年在非常规致密油气资源开采中出现的新工艺, 已在世界范围内推广实施, 并取得了良好的增产效果. 该技术可使支撑剂在人工压裂缝中形成簇团式分布, 从而形成油气高速流动通道, 提高裂缝的导流能力. 但目前对于高速通道压裂裂缝高导流能力的形成机理及其影响因素尚不清楚. 对此, 本文从流体力学理论出发, 首先将高速通道压裂裂缝内形成的支撑剂簇团视为渗流区域, 簇团间的大通道视为自由流动区域; 然后基于Darcy-Brinkman方程建立了裂缝内的流动数学模型, 采用均匀化理论对该流动数学模型进行了尺度升级, 推导得到了高速通道压裂裂缝的渗透率, 揭示了其高导流能力的形成机理; 并以此为基础, 分析了不同支撑剂簇团形状、大小以及分布方式等因素对其导流能力的影响, 可为高速通道压裂工艺参数设计与优化提供基础.     

Homogenization of the Signorini boundary‐value problem in a thick junction and boundary integral equations for the homogenized problem

We consider a mixed boundary‐value problem for the Poisson equation in a thick junction Ω_ε which is the union of a domain Ω₀ and a large number of ε—periodically situated thin cylinders. The non‐uniform Signorini conditions are given on the lateral surfaces of the cylinders. The asymptotic analysis of this problem is done as ε→0, i.e. when the number of the thin cylinders infinitely increases and their thickness tends to zero. We prove a convergence theorem and show that the non‐uniform Signorini boundary conditions are transformed in the limiting variational inequalities in the region that is filled up by the thin cylinders as ε→0. The convergence of the energy integrals is proved as well. The existence and uniqueness of the solution to this non‐standard limit problem is established. This solution can be constructed by using a penalty formulation and successive iteration. For some subclass, these problems can be reduced to an obstacle problem in Ω₀ and an appropriate postprocessing. The equations in Ω₀ finally are also treated with boundary integral equations. Copyright © 2010 John Wiley & Sons, Ltd.     

Dual-stage nested homogenization for rate-dependent anisotropic elasto-plasticity model of dendritic cast aluminum alloys

This paper proposes a nested dual-stage homogenization method for developing microstructure based continuum elasto-viscoplastic models for large secondary dendrite arm spacing or SDAS cast aluminum alloys. Microstructures of these alloys are characterized by extremely inhomogeneous distribution of inclusions along the dendrite cell boundaries. Traditional single-step homogenization methods are not suitable for this type of microstructure due to the size of the representative volume element (RVE) and the associated computations required for micromechanical analyses. To circumvent this limitation, two distinct RVE’s or statistically equivalent RVE’s are identified, corresponding to the inherent scales of inhomogeneity in the microstructure. The homogenization is performed in multiple stages for each of the RVE’s identified. The macroscopic behavior is described by a rate-dependent, anisotropic homogenization based continuum plasticity (HCP) model. Anisotropy and viscoplastic parameters in the HCP model are calibrated from homogenization of micro-variables for the different RVE’s. These parameters are dependent on microstructural features such as morphology and distribution of different phases. The uniqueness of the nested two-stage homogenization is that it enables evaluation of the overall homogenized model parameters of the cast alloy from limited experimental data, but also material parameters of constituents like inter-dendritic phase and pure aluminum matrix. The capabilities of the HCP model are demonstrated for a cast aluminum alloy AS7GU having a SDAS of 30 μm.     

复合材料周期结构数学均匀化方法的一种新型单胞边界条件

数学均匀化方法是计算周期复合材料结构的有效方法之一,单胞边界条件施加的合理性直接决定了影响函数控制方程的计算效率和精度,进而影响均匀化弹性参数和摄动位移的计算精度.本文首先将单胞影响函数作为虚拟位移处理,给出了单胞在结构中真实的边界条件,结果表明,四边固支适合作为二维结构单胞边界条件;其次,针对二维结构提出了超单胞周期边界条件,有效提高了影响函数的计算精度,并使用与虚拟位移相对应的虚拟势能泛函验证超单胞周期边界条件的有效性;最后,利用数值分析验证多尺度渐进展开方法的计算精度,强调了二阶摄动的必要性.