Reaching multi‐nanosecond timescales in combined QM/MM molecular dynamics simulations through parallel horsetail sampling

We report an enhanced sampling technique that allows to reach the multi‐nanosecond timescale in quantum mechanics/molecular mechanics molecular dynamics simulations. The proposed technique, called horsetail sampling, is a specific type of multiple molecular dynamics approach exhibiting high parallel efficiency. It couples a main simulation with a large number of shorter trajectories launched on independent processors at periodic time intervals. The technique is applied to study hydrogen peroxide at the water liquid–vapor interface, a system of considerable atmospheric relevance. A total simulation time of a little more than 6 ns has been attained for a total CPU time of 5.1 years representing only about 20 days of wall‐clock time. The discussion of the results highlights the strong influence of the solvation effects at the interface on the structure and the electronic properties of the solute. © 2017 Wiley Periodicals, Inc.     

Fictitious domain methods for two‐phase flow energy balance computations in nuclear components

This paper is dedicated to the numerical simulation of nuclear components (cores and steam generators) by fictitious domain methods. The fictitious domain approach consists in immersing the physical domain under study in a Cartesian domain, called the fictitious domain, and in performing the numerical resolution on this fictitious domain. The calculation times are then efficiently reduced by the use of fast solvers. In counterpart, one has to handle with an immersed boundary, generally non‐aligned with the Cartesian mesh, which can be non‐trivial. The two fictitious domain methods compared here on industrial simulations and developed by Ramière et al. deal with an approximate immersed interface directly derived from the uniform Cartesian mesh. All the usual immersed boundary conditions (Dirichlet, Robin, Neumann), possibly mixed, are handled through a unique formulation of the fictitious problem. This kind of approximation leads to first‐order methods in space that exhibit a good ratio of the precision of the approximate solution over the CPU time, which is very important for industrial simulations. After a brief recall of the fictitious domain method with spread interface (Ramière et al., CMAME 2007) and the fictitious domain method with immersed jumps (Ramière et al., JCP 2008), we will focus on the numerical results provided by these methods applied to the energy balance equation in a steam generator. The advantages and drawbacks of each method will be pointed out. Generally speaking, the two methods confirm their very good efficiency in terms of precision, convergence, and calculation time in an industrial context. Copyright © 2011 John Wiley & Sons, Ltd.     

Coupled Navier–Stokes/molecular dynamics simulations in nonperiodic domains based on particle forcing

This paper focuses on coupling methods for hybrid Navier–Stokes/molecular dynamics (MD) simulations. The computational domain is split in a continuum flow region, where a finite‐volume discretisation of the Navier–Stokes equations is used, and one or more particle domains, where molecular level modelling of the flow is employed. The domains are defined with a partial overlap, in which the flow states are coupled through an exchange of the velocity components. For the steady flows considered, an under‐relaxed Newton iteration method is used to drive the coupled system to convergence. The main focus of the present work is on methods to impose nonperiodic boundary conditions on the particle domain(s). A particle forcing is applied in the direction normal to the particle domain boundary to impose the boundary normal velocity component. A novel aspect of the present work is the extension of this method to more general nonplanar particle domain boundaries. The main contribution of the paper is the development of a particle forcing method in the direction tangential to the domain boundary, which is based on the equivalent continuum‐flow boundary shear stresses along with an iterative forcing strength adjustment based on the extrapolated particle boundary velocity. Furthermore, an adaptation scheme is presented, which uses the finite‐volume flux residuals of the particle bin averaged velocity field as a truncation criterion for the iterative force‐update scheme. It is demonstrated that by comparing the residual reduction for the momentum equation in the nonhomogeneous directions during the molecular dynamics simulations with that for a homogeneous direction, the forcing iteration at which the statistical noise in the velocity field dominates the uncertainty in the forcing strength can be determined. At this point the iteration can be truncated. It is shown that with adaptive schemes of this type, the total number of MD evaluations required in a coupled Navier–Stokes/MD simulation can be reduced relative to a hybrid scheme with a fixed number of forcing‐strength updates. Copyright © 2011 John Wiley & Sons, Ltd.     

A perturbed martingale approach to global optimization

A new global stochastic search, guided mainly through derivative-free directional information computable from the sample statistical moments of the design variables within a Monte Carlo setup, is proposed. The search is aided by imparting to the directional update term additional layers of random perturbations referred to as ‘coalescence’ and ‘scrambling’. A selection step, constituting yet another avenue for random perturbation, completes the global search. The direction-driven nature of the search is manifest in the local extremization and coalescence components, which are posed as martingale problems that yield gain-like update terms upon discretization. As anticipated and numerically demonstrated, to a limited extent, against the problem of parameter recovery given the chaotic response histories of a couple of nonlinear oscillators, the proposed method appears to offer a more rational, more accurate and faster alternative to most available evolutionary schemes, prominently the particle swarm optimization.     

线缆系统电磁脉冲效应测量系统研制和试验

分析了实验室瞬态X射线产生的系统电磁脉冲(SGEMP)效应测试所面临的技术问题,提出了解决方法、措施以及实验室模拟瞬态X射线的SGEMP模拟试验方法。通过电子屏蔽、电磁屏蔽、光电隔离、信号对称提取等特殊技术处理,解决了SGEMP效应模拟试验方法和测量系统抗X射线、抗电磁辐射等技术问题,并在大型瞬态X射线模拟源上,测出了瞬态X射线辐照时金属腔内线缆的SGEMP效应波形及幅值。     

NONLINEAR DYNAMICS AND SYNCHRONIZATION OF TWO COUPLED PIPES CONVEYING PULSATING FLUID

In this paper, the nonlinear dynamical behavior of two coupled pipes conveying pulsating fluid is studied. The connection between the two pipes is considered as a distributed linear spring. Based on this consideration, the equations of motion of the coupled two-pipe system are obtained. The two coupled nonlinear partial differential equations, discretized using the fourth- order Galerkin method, are solved by a fourth-order Runge-Kutta integration algorithm. Results show that the connection stiffness has a significant effect on the dynamical behavior of the coupled system. It is found that for some parameter values the motion types of the two pipes might be synchronous.     

Molecular dynamics simulations of displacement cascades in Fe-10%Cr systems

Molecular dynamics simulations of the displacement cascades in Fe-10%Cr systems are used to simulate the primary knocked-on atom events of the irradiation damage at temperatures 300,600,and 750 K with primary knockedon atom energies between 1 and 15 keV.The results indicate that the vacancies produced by the cascade are all in the central region of the displacement cascade.During the cascade,all recoil Fe and Cr atoms combine with each other to form Fe-Cr or Fe-Fe interstitial dumbbells as well as interstitial clusters.The number and the size of interstitial clusters increase with the energy of the primary knocked-on atom and the temperature.A few large clusters consist of a large number of Fe interstitials with a few Cr atoms,the rest are Fe-Cr clusters with small and medium sizes.The interstitial dumbbells of Fe-Fe and Fe-Cr are in the 111 and 110 series directions,respectively.     

Recursive dynamics simulator (ReDySim): A multibody dynamics solver

Recursive formulations have significantly helped in achieving real-time computations and model-based control laws. The recursive dynamics simulator (ReDySim) is a MATLAB-based recursive solver for dynamic analysis of multibody systems. ReDySim delves upon the decoupled natural orthogonal complement approach originally developed for serial-chain manipulators. In comparison to the commercially available software, dynamic analyses in ReDySim can be performed without creating solid model. The input parameters are specified in MATLAB environment. ReDySim has the capability to incorporate any control algorithm with utmost ease. In this work, the capabilities of ReDySim for solving open-loop and closed-loop systems are shown by examples of robotic gripper, KUKA KR5 industrial manipulator and four-bar mechanism. ReDySim can be downloaded for free from http://www.redysim.co.nr and can be used almost instantly.     

LED微阵列投影系统设计

为满足便携式投影仪的市场需求,设计了一种基于LED微型阵列的投影系统。该系统由显示单元和投影物镜构成。采用尺寸为12 mm×9 mm的自发光LED微型阵列作为系统的显示单元,利用光学设计软件设计了投影物镜。投影物镜采用反远距光学结构,全视场角为80°,焦距为8 mm,属于强光、广角镜头。在空间频率20 lp/mm处,该物镜的调制传递函数大于0. 85,畸变小于2%,符合投影系统的设计要求。该投影系统具有体积小,结构简单,投影效果好,易加工等诸多优势,可为第三代投影技术的发展提供参考。