Large-scale computer simulation of fully developed turbulent channel flow with heat transfer

Recently, with the advent of supercomputers, there has been considerable interest in the use of direct numerical simulation to obtain information about turbulent shear flow at low Reynolds number. This paper presents a pseudospectral technique to solve the full three-dimensional time-dependent Navier-Stokes and advection-diffusion equations without the use of subgrid-scale modelling. The technique has not been previously used for fully developed turbulent channel flow simulation and is based on methods applied in other contexts. The emphasis of this paper is to provide a reasonably detailed account of how the simulation is done rather than to present new calculations of turbulence. The details of an algorithm for turbulent channel flow simulation and the grid and time step sizes needed to integrate through transient behaviour to steady state turbulence have not been published before and are presented here. Results from a Cray-2 simulation of fully developed turbulent flow in a channel with heat transfer are presented along with a critical comparison between experiment and computation. The first- and second-order moments agree well with experimental measurements; the agreement is poor for higher-order moments such as the skewness and flatness near the walls of the channel. Detailed information given about the effects of spatial grid resolution on a computed results is important for estimating the size of the computation required to study various aspects of a turbulent flow.     

Uniform random number generation

In typical stochastic simulations, randomness is produced by generating a sequence of independent uniform variates (usually real-valued between 0 and 1, or integer-valued in some interval) and transforming them in an appropriate way. In this paper, we examine practical ways of generating (deterministic approximations to) such uniform variates on a computer. We compare them in terms of ease of implementation, efficiency, theoretical support, and statistical robustness. We look in particular at several classes of generators, such as linear congruential, multiple recursive, digital multistep, Tausworthe, lagged-Fibonacci, generalized feedback shift register, matrix, linear congruential over fields of formal series, and combined generators, and show how all of them can be analyzed in terms of their lattice structure. We also mention other classes of generators, like non-linear generators, discuss other kinds of theoretical and empirical statistical tests, and give a bibliographic survey of recent papers on the subject.     

Standing waves, clustering, and phase waves in 1D simulations of kinetic relaxation oscillations in NO+NH3 on Pt(1 0 0) coupled by diffusion

The Lombardo–Imbihl–Fink (LFI) ODE model of the NO+NH₃ reaction on a Pt(1 0 0) surface shows stable relaxation oscillations with very sharp transitions for temperatures T between 404 and 433 K. Here we study numerically the effect of linear diffusive coupling of these oscillators in one spatial dimension. Depending on the parameters and initial conditions we find a rich variety of spatio-temporal patterns which we group into four main regimes: bulk oscillations (BOs), standing waves (SW), phase clusters (PC), and phase waves (PW). Two key ingredients for SW and PC are identified, namely the relaxation type of the ODE oscillations and a nonlocal (and nonglobal) coupling due to relatively fast diffusion of the kinetically slaved variables NH₃ and H. In particular, the latter replaces the global coupling through the gas phase used to obtain SW and PC in models of related surface reactions. The PW exist only under the assumption of (relatively) slow diffusion of NH₃ and H.     

Validation and verification of social processes within agent-based computational organization models

The use of simulation modeling in computational analysis of organizations is becoming a prominent approach in social science research. However, relying on simulations to gain intuition about social phenomena has significant implications. While simulations may give rise to interesting macro-level phenomena, and sometimes even mimic empirical data, the underlying micro and macro level processes may be far from realistic. Yet, this realism may be important to infer results that are relevant to existing theories of social systems and to policy making. Therefore, it is important to assess not only predictive capability but also explanation accuracy of formal models in terms of the degree of realism reflected by the embedded processes. This paper presents a process-centric perspective for the validation and verification (V&V) of agent-based computational organization models. Following an overview of the role of V&V within the life cycle of a simulation study, emergent issues in agent-based organization model V&V are outlined. The notion of social contract that facilitates capturing micro level processes among agents is introduced to enable reasoning about the integrity and consistency of agent-based organization designs. Social contracts are shown to enable modular compositional verification of interaction dynamics among peer agents. Two types of consistency are introduced: horizontal and vertical consistency. It is argued that such local consistency analysis is necessary, but insufficient to validate emergent macro processes within multi-agent organizations. As such, new formal validation metrics are introduced to substantiate the operational validity of emergent macro-level behavior. Levent Yilmaz is Assistant Professor of Computer Science and Engineering in the College of Engineering at Auburn University and co-founder of the Auburn Modeling and Simulation Laboratory of the M&SNet. Dr. Yilmaz received his Ph.D. and M.S. degrees from Virginia Polytechnic Institute and State University (Virginia Tech). His research interests are on advancing the theory and methodology of simulation modeling, agent-directed simulation (to explore dynamics of socio-technical systems, organizations, and human/team behavior), and education in simulation modeling. Dr. Yilmaz is a member of ACM, IEEE Computer Society, Society for Computer Simulation International, and Upsilon Pi Epsilon. URL: http://www.eng.auburn.edu/~yilmaz     

Coordination and ligand exchange dynamics of solvated metal ions

Recent developments in computer speed and capacity have opened the access to highly accurate molecular dynamics simulations based on quantum mechanically calculated forces for the chemically relevant region around ions in solution (QM/MM formalism). This accuracy, although still extremely consuming (30-300 days of CP time per simulation), is needed for reliable structural details and ligand exchange rates. A large number of main group and transition metal ions have been investigated by this approach, giving very detailed insight into the properties of these ions in solution and allowing to classify the ions by various characteristics. Most first-row transition metal ions have a very stable first hexa-coordinated solvation shell, whose vibrational distortions, however, strongly influence the dynamics of the second shell. The dynamical Jahn-Teller effect - shown to be a femto- and picosecond phenomenon - can strongly influence ligand coordination and exchange dynamics. A large number of ions with very labile solvation shell such as most main group ions, but also transition metal ions, e.g. Ag(I) and Hg(II), can change their coordination within the picosecond scale, leading to an almost simultaneous presence of several species hardly accessible by present experimental techniques. Among these ions, the structure breakers are of particular interest, and it could be shown that there are two types of them, one with a large and very labile first coordination shell such as Cs(I), the other characterised by a small first but an unusually large second solvation shell such as Au(I). Investigations of metal ions coordinated to ammonia ligands have shown that coordination to hetero-atoms can accelerate the ligand exchange reaction rates by several orders of magnitude, e.g. for Cu(II) and Ni(II). Simulations of ions in aqueous ammonia gave a very detailed picture of the complexity of species almost simultaneously present and illustrate the enormous difficulties encountered when trying to fit X-ray or neutron diffraction data for such systems. In general, ligand exchange rates situated in the picosecond range are far below the NMR scale, and as femtosecond laser pulse spectroscopy could not be applied so far to ionic solutions, accurate simulations have become a very important tool to access structure and dynamics of solvated ions. A number of VIDEO clips supplied on the Web as supporting material illustrates the processes occurring in solutions of the metal ions.     

Dual microband electrodes: current distributions and diffusion layer ‘titrations’. Implications for electroanalytical measurements

The simulation of transport to double microband electrodes in generator–collector mode is reported focusing especially on the ‘titration curve’ approach to electroanalysis in which a titrant is electrogenerated from a redox active precursor on the generator electrode and reacts homogeneously with the target analyte. The current on the detector electrode reflects the amount of titrant ‘surviving’ passage between the two electrodes. The form of the titration curve – plots of detector current as a function of generator current – is shown to be highly sensitive to the electrode kinetics of the redox couple driven at the generator electrode. Accordingly the naïve use of such methodology for analysis without accompanying simulation and kinetic analysis is fraught with danger. Use of the conformal mapping approach in combination with the ADI method for investigation of the ‘titration’ current distributions at the double band system gives fast and precise simulation of this and similar problems. Convergence analysis is described which allows for the automatic selection of the simulation grid size so as to obtain a chosen accuracy (for example 1%) of the current for all experimentally meaningful values of the geometrical and physico-chemical parameters of the system to be investigated.     

Dose prediction and process optimization in a gamma sterilization facility using 3-D Monte Carlo code

A model of a gamma sterilizer was built using the ITS/ACCEPT Monte Carlo code and verified through dosimetry. Individual dosimetry measurements in homogeneous material were pooled to represent larger bodies that could be simulated in a reasonable time. With the assumptions and simplifications described, dose predictions were within 2–5% of dosimetry. The model was used to simulate product movement through the sterilizer and to predict information useful for process optimization and facility design.     

超临界NaCl水溶液的分子动力学模拟

采用分子动力学模拟的方法对超临界NaCl水溶液的微观结构进行了研究.模拟发现在所研究超临界条件下,密度的变化比温度的变化对超临界NaCl水溶液的微观结构影响更大.温度及密度对Cl－ H2O径向分布函数的影响比对Na＋ H2O径向分布函数的影响要大.超临界条件下,各gNa＋－Cl－在0.261 nm处出现峰值,表明Na＋、Cl－之间发生了离子的缔合.超临界条件下,随温度增加,缔合作用增强；随密度增加,缔合作用减弱.本文工作为建立可适用于超临界条件下的电解质热力学模型提供了依据.     

流体熵相关性质的Monte Carlo模拟新方法

提出了用于计算实际体系熵相关性质的Monte Carlo多级取样分子模拟方法. 应用这一方法, 对硬球流体的化学势及Helmholtz自由能进行了估算, 得到了满意的结果. 计算化学势时, 不存在通常试验粒子方法所遇到的高密度问题. 该方法特别适合规律性的系统研究, 较之普通模拟方法要有效得多. 模拟得到的硬球体系无限稀释组份的超额化学势与对比直径的关系, 在相变区域为一条双凹曲线; 无论是在相变区还是在单相区, Carnahan-Starling公式对这一关系的描述均有较大偏差.     

The Discrete Approximation of Continuous-state Nonlinear Branching Processes in Lévy Random Environments北大核心CSCD

In this paper, we establish the discrete approximation of continuous-state nonlinear branching processes in Lévy random environments by using tightness and convergence sequence in infinite dimensional product space via stochastic differential equations. Taking α-stable branching as an example, the conditions which are given to discretize continuous-state nonlinear branching processes in Lévy random environments are verified. © 2022 Chinese Academy of Sciences. All rights reserved.