Level‐set topology optimization with many linear buckling constraints using an efficient and robust eigensolver

Linear buckling constraints are important in structural topology optimization for obtaining designs that can support the required loads without failure. During the optimization process, the critical buckling eigenmode can change; this poses a challenge to gradient‐based optimization and can require the computation of a large number of linear buckling eigenmodes. This is potentially both computationally difficult to achieve and prohibitively expensive. In this paper, we motivate the need for a large number of linear buckling modes and show how several features of the block Jacobi conjugate gradient (BJCG) eigenvalue method, including optimal shift estimates, the reuse of eigenvectors, adaptive eigenvector tolerances and multiple shifts, can be used to efficiently and robustly compute a large number of buckling eigenmodes. This paper also introduces linear buckling constraints for level‐set topology optimization. In our approach, the velocity function is defined as a weighted sum of the shape sensitivities for the objective and constraint functions. The weights are found by solving an optimization sub‐problem to reduce the mass while maintaining feasibility of the buckling constraints. The effectiveness of this approach in combination with the BJCG method is demonstrated using a 3D optimization problem. Copyright © 2016 John Wiley & Sons, Ltd.     

MEMS器件虚拟运行中运动规律求解器的设计与实现

针对虚拟运行（在器件动态模型的基础上对MEMS器件运动性能进行评测的设计工具,其作用在于建立MEMS设计参数与器件运行之间的直接联系）中运动规律的联立、求解和实时仿真等问题,分析了几种典型MEMS器件的动态模型,提出了运动规律求解器的设计原理和体系结构,其中包括三项关键技术：物理运动规律在部件库系统中的形式化表示,物理问题自动求解的推理算法,以及离线编译的优化技术。最后,设计和实现了运动规律求解器原型系统,并在该系统中测试了微泵的虚拟运行,达到了通过修改阀片设计参数和膜片驱动参数实时观察虚拟运行中状态相应变化的目标,即在微泵的设计参数与其运行之间建立了直接的关联关系。     

Hybrid time domain solvers for the Maxwell equations in 2D

We present a new approach to time domain hybrid schemes for the Maxwell equations. By combining the classical FD‐TD scheme with two unstructured solvers, one explicit finite volume solver and one implicit finite element solver, we achieve a very efficient and flexible second‐order scheme. The second‐order accuracy of the hybrid scheme is verified through convergence studies on perfectly conducting as well as dielectric and diamagnetic circular cylinders. The numerical results also show its superiority to the FD‐TD scheme. Copyright © 2002 John Wiley & Sons, Ltd.     

半功率带宽法与INV阻尼计法求阻尼比的研究

半功率带宽法是国内外常用的测试阻尼的方法，笔者发现半功率带宽法由于受到信号处理中的①采样频率SF②频率分辨率△f③采样频率与信号频率之比K④采样分析的长度T和⑤阻尼比D的大小等因素的影响使得测得的阻尼比值波动很大，除少数状态是可信的外，多数状态是不稳定的，有时误差惊人，得到的测试数据不能使用，易导致错误的结论。文内提出了一种创新的INV阻尼计法，它具有高精度、高准确率、操作方便等特点，使阻尼比测试有了一个突破。文内把两种方法与时域波形衰减法求阻尼比进行了对比分析，以供同行专家参考。     

Domain decomposition and multigrid solvers for flow simulation in porous media on distributed memory parallel processors

Realistic simulations of fluid flow in oil reservoirs have been proven to be computationally intensive. In this work, techniques for solving large sparse systems of linear equations that arise in simulating these processes are developed for parallel computers such as INTEL hypercubes iPSC/2 and iPSC/860. This solver is based on a combined multigrid and domain decomposition approach. The Algorithm uses line corrections solved using a multigrid method, line Jacobi and block incomplete domain decomposition as an overall preconditioner for a conjugate gradient-like acceleration method, ORTHOMIN (k). This is shown to be a factor of ten times faster on a 32-processor hypercube compared to widely used sequential solvers. Three test problems are used to validate the results which include implicit wells and faults: The first is based on highly heterogeneous two-phase flow, the second on the SPE Third Comparative Solution and the third on real production compositional data.     

一类非主型方程古沙及柯西问题的解

应用文献提出的双曲型方程的Riemann函数法,并利用超几何函数的特殊性质对一类非主型方程的古沙问题和柯西问题提供求解的方法,给出了解的表达式。     

Implementation of a parallel unstructured Euler solver on shared- and distributed-memory architectures

An efficient three-dimensional unstructured Euler solver is parallelized on a CRAY Y-MP C90 shared-memory computer and on an Intel Touchstone Delta distributed-memory computer. This paper relates the experiences gained and describes the software tools and hardware used in this study. Performance comparisons between the two differing architectures are made.This work was sponsored in part by ARPA (NAG-1-1485) and by NASA Contract No. NAS1-19480 while authors Mavriplis, Saltz and Das were in residence at ICASE, NASA Langley Research Center, Hampton, Virginia. This research was performed in part using the Intel Touchstone Delta System operated by Caltech on behalf of the Concurrent Supercomputing Consortium. Access to this fecility was provided by NASA Langley Research Center and the Center for Research in Parallel Processing. The content of the information does not necessarily reflect the position or the policy of the Government and no official endorsement should be inferred.     

A comparison of the different extensions of a weak formulation of an approximate riemann solver for supercritical flows and their relationship to existing schemes

A weak formulation of the Roe linearised Riemann solver proposed recently is examined in the context of the equations governing steady, supercritical flows. Two possible alternatives present themselves, and it is shown that these are equivalent and, moreover, are equivalent to existing schemes.     

A numerical study of combined forced and natural convection heat transfer in a vertical cylinder bundle

Combined forced and natural convection in a bundle of vertical cylinders have been studied numerically with the aid of Riemann geometry. The cylinders heated uniformly are arranged in an equilateral triangular pattern of three bundle spacings S/D = 1.1, 1.5, and 1.9. Numerical calculations are made of the mixed convection of air (Pr = 0.7) in a subchannel of the bundle at ratios of modified Grashof number to Reynolds number Gr*/Re from 1 to 10⁵. The results show that the reverse flows appear along the center line in the subchannel with increasing Gr*/Re. Local Nusselt numbers for the thermally fully‐developed region are well correlated with Gr*/Re. © 2003 Wiley Periodicals, Inc. Heat Trans Asian Res, 32(7): 639–649, 2003; Published online in Wiley InterScience ( www.interscience.wiley.com ). DOI 10.1002/htj.10118     

针对PRESENT分组密码算法的代数分析*

本文研究针对PRESENT分组密码的代数分析。通过使用S盒的表达式形式,构建出多轮PRESENT加密中的代数方程组。这种构建方程的方法被推广到具有小型S盒的典型SPN型分组密码算法的方程构建问题中。文中还对简化的PRESENT算法进行了攻击实验。采用MiniSAT作为攻击过程中的求解工具,对4轮、6轮PRESENT加密进行实际攻击。可以在一分钟之内恢复4轮加密的所有密钥,数小时内恢复6轮加密的密钥。并且通过引入了差分思想,首次将有效攻击轮数提高到8轮。