Consistency between independence theorems and generalized self-consistent method

Recently Zheng & Hwang established a series of independence theorems concerning with planar effective elastic properties. It is manifested that the estimation of the effective elastic properties of microcracked solids through the generalized self-consistent method (GSCM) contradicts with these independence theorems. In this paper it is shown that such contradiction is actually caused by the approximate algorithm adopted, while the exact solution of GSCM is consistent with these rigorously established independence theorems. Since only an approximate algorithm in GCSM is available in dealing with problems involving non-circular inclusions or holes, an intrinsic GSCM is proposed, which can be performed based on an approximate algorithm and the corresponding estimations are consistent with the independence theorems.     

Unique Orbital Normal Form for Vector Fields of Hopf-Zero Singularity

Normal forms for vector fields of Hopf-zero singularity in R³ are studied. Multiple Lie bracket method is used to give unique normal forms under both conjugacy and orbital equivalence for such vector fields with a generic quadratic part.Dedicated to Professor Shui-Nee Chow on the occasion of his 60th birthday.     

A Grobman–Hartman Theorem for Control Systems

We consider the problem of locally linearizing a control system via topological transformations. According to [2,3], there is no naive generalization of the classical Grobman–Hartman theorem for ODEs to control systems: a generic control system, when viewed as a set of under-determined differential equations parametrized by the control, cannot be linearized using pointwise transformations on the state and the control values. However, if we allow the transformations to depend on the control at a functional level (open loop transformations), we are able to prove a version of the Grobman–Hartman theorem for control systems.     

强迫谐振动下连续体结构拓扑优化

应用结构拓扑优化ICM(独立连续映射)方法,对强迫谐振动下结构拓扑优化问题建立了以重量极小为目标,位移幅值为约束的优化模型.位移幅值采用一阶泰勒展式近似,由于拓扑优化中设计变量数目通常很多,对强迫谐振动位移幅值的敏度分析推导了伴随法公式,使得一次敏度分析可以计算出对所有设计变量的偏导数,克服了采用直接法敏度分析中一次只能计算出对一个设计变量的偏导数的不足.算例表明用伴随法分析敏度在结构拓扑优化中可以大幅提高计算效率,ICM方法采用独立于截面及形状参数的拓扑优化设计变量更清晰地反映了拓扑优化的本质.     

Analytic solutions of a class of nonlinear partial differential equations

An approach is presented for computing the adjoint operator vector of a class of nonlinear （that is, partial-nonlinear） operator matrices by using the properties of conjugate operators to generalize a previous method proposed by the author. A unified theory is then given to solve a class of nonlinear （partial-nonlinear and including all linear） and non-homogeneous differential equations with a mathematical mechanization method. In other words, a transformation is constructed by homogenization and triangulation, which reduces the original system to a simpler diagonal system. Applications are given to solve some elasticity equations.     

A continuous adjoint formulation for the computation of topological and surface sensitivities of ducted flows

Topology optimization of fluid dynamic systems is a comparatively young optimal design technique. Its central ingredient is the computation of topological sensitivity maps. Whereas, for finite element solvers, implementations of such sensitivity maps have been accomplished in the past, this study focuses on providing this functionality within a professional finite volume computational fluid dynamics solver. On the basis of a continuous adjoint formulation, we derive the adjoint equations and the boundary conditions for typical cost functions of ducted flows and present first results for two‐ and three‐dimensional geometries. Emphasis is placed on the versatility of our approach with respect to changes in the objective function. We further demonstrate that surface sensitivity maps can also be computed with the implemented functionality and establish their connection with topological sensitivities. Copyright © 2008 John Wiley & Sons, Ltd.     

胶凝原油蠕变的时间-温度-应力等效性

黏弹性蠕变是胶凝原油重要的流变特性,要描述胶凝原油的黏弹性蠕变行为,必须同时考虑 应力、时间和温度等参数. 于是对胶凝原油进行了不同应力和不同温度下的蠕变实验,利用 时间-温度-应力等效原理分析了胶凝原油的非线性蠕变行为,得到了参考应力和参考温度下 的蠕变主曲线,并用黏弹性流变模型拟合得到了蠕变主曲线的表达式. 从而可以分别通过较 高温度和应力水平下的短期蠕变行为来预测较低温度和应力水平下的长期蠕变行为.     

Efficient computation of operator‐type response sensitivities for uncertainty quantification and predictive modeling: illustrative application to a spent nuclear fuel dissolver model

This work honors the 75th birthday of Professor Ionel Michael Navon by presenting original results highlighting the computational efficiency of the adjoint sensitivity analysis methodology for function‐valued operator responses by means of an illustrative paradigm dissolver model. The dissolver model analyzed in this work has been selected because of its applicability to material separations and its potential role in diversion activities associated with proliferation and international safeguards. This dissolver model comprises eight active compartments in which the 16 time‐dependent nonlinear differential equations modeling the physical and chemical processes comprise 619 scalar and time‐dependent model parameters, related to the model's equation of state and inflow conditions. The most important response for the dissolver model is the time‐dependent nitric acid in the compartment furthest away from the inlet, where measurements are available at 307 time instances over the transient's duration of 10.5 h. The sensitivities to all model parameters of the acid concentrations at each of these instances in time are computed efficiently by applying the adjoint sensitivity analysis methodology for operator‐valued responses. The uncertainties in the model parameters are propagated using the above‐mentioned sensitivities to compute the uncertainties in the computed responses. A predictive modeling formalism is subsequently used to combine the computational results with the experimental information measured in the compartment furthest from the inlet and then predict optimal values and uncertainties throughout the dissolver. This predictive modeling methodology uses the maximum entropy principle to construct an optimal approximation of the unknown a priori distribution for the a priori known mean values and uncertainties characterizing the model parameters and the computed and experimentally measured model responses. This approximate a priori distribution is subsequently combined using Bayes' theorem with the “likelihood” provided by the multi‐physics computational models. Finally, the posterior distribution is evaluated using the saddle‐point method to obtain analytical expressions for the optimally predicted values for the parameters and responses of both multi‐physics models, along with corresponding reduced uncertainties. This work shows that even though the experimental data pertains solely to the compartment furthest from the inlet (where the data were measured), the predictive modeling procedure used herein actually improves the predictions and reduces the predicted uncertainties for the entire dissolver, including the compartment furthest from the measurements, because this predictive modeling methodology combines and transmits information simultaneously over the entire phase‐space, comprising all time steps and spatial locations. Copyright © 2016 John Wiley & Sons, Ltd.     

基于拓扑优化的声学结构材料分布设计

本文针对结构的声学设计问题进行研究,通过优化两种不同的材料在结构设计域内的拓扑分 布来最小化谐振结构所产生的声场中指定参考面/参考域内的声压。在研究中假定结构为线弹性小变 形结构,材料阻尼为Rayleigh阻尼,声学介质为无粘、可压缩、小扰动流体。对结构响应采用有限 元格式进行计算,对声场采用基于Helmholtz积分的边界元格式进行计算,由于声场在无穷远自由边 界的无反射条件在边界积分中能自动得到满足,该格式特别适合于具有开放边界的声场计算。建立 了结构有限元-声场边界元格式的耦合系统拓扑优化模型,导出了耦合系统敏感度分析的一般格式及 伴随格式。数值算例验证了所提出的结构-声学耦合系统优化方法的有效性和可靠性,并揭示了基于 声学准则的拓扑优化结果的有关特性 关键词边界积分,结构声学耦合系统,拓扑优化,敏感度分析,伴随方法     

On initial,boundary conditions and viscosity coefficient control for Burgers' equation

In order to use the optimal control techniques in models of geophysical flow circulation, an application to a 1D advection–diffusion equation, the so-called Burgers' equation, is described. The aim of optimal control is to find the best parameters of the model which ensure the closest simulation to the observed values. In a more general case, the continuous problem and the corresponding discrete form are formulated. Three kinds of simulation are realized to validate the method. Optimal control processes by initial and boundary conditions require an implicit discretization scheme on the first time step and a decentered one for the non-linear advection term on boundaries. The robustness of the method is tested with a noised dataset and random values of the initial controls. The optimization process of the viscosity coefficient as a time- and space-dependent variable is more difficult. A numerical study of the model sensitivity is carried out. Finally, the numerical application of the simultaneous control by the initial conditions, the boundary conditions and the viscosity coefficient allows a possible influence between controls to be taken into account. These numerical experiments give methodological rules for applications to more complex situations. © 1998 John Wiley & Sons, Ltd.