Econometric analysis of microscopic simulation models

Microscopic simulation models are often evaluated based on visual inspection of the results. This paper presents formal econometric techniques to compare microscopic simulation (MS) models with real-life data. A related result is a methodology to compare different MS models with each other. For this purpose, possible parameters of interest, such as mean returns, or autocorrelation patterns, are classified and characterized. For each class of characteristics, the appropriate techniques are presented. We illustrate the methodology by comparing the MS model developed by He and Li [J. Econ. Dynam. Control, 2007, 31, 3396-3426, Quant. Finance, 2008, 8, 59-79] with actual data.

Simulation of coupled heat and moisture transfer in air-conditioned buildings

The simultaneous heat and moisture transfer in the building envelope has an important influence on the indoor environment and the overall performance of buildings. In this paper, a model for predicting whole building heat and moisture transfer was presented. Both heat and moisture transfer in the building envelope and indoor air were simultaneously considered; their interactions were modeled. The coupled model takes into account most of the main hygrothermal effects in buildings. The coupled system model was implemented in MATLAB-Simulink, and validated by using a series of published testing tools. The new program was applied to investigate the moisture transfer effect on indoor air humidity and building energy consumption under different climates. The results show that the use of more detailed simulation routines can result in improvements to the building's design for energy optimisation through the choice of proper hygroscopic materials, which would not be indicated by simpler calculation techniques.

Simulation of whole building coupled hygrothermal-airflow transfer in different climates

The coupled heat, air and moisture transfer between building envelopes and indoor air is complicated, and has a significant influence on the indoor environment and the energy performance of buildings. In the paper, a model for predicting coupled multi-zone hygrothermal-airflow transfer is presented. Both heat and moisture transfer in the building envelope and multi-zone indoor airflow are simultaneously considered; their interactions are modeled. The coupled system model is implemented into Matlab–Simulink, and is validated by using a series of testing tools and experiments. The new program is applied to investigate the moisture transfer effect on indoor air humidity and building energy consumption in different climates (hot-humid, temperate and hot-dry climates). The results show that not accounting for hygrothermal effects in modeling will result in overestimation of energy costs for hot and humid climate situations and possible over sizing of plant leading to inefficient operation.

Particle simulation study of electron heating by counter-streaming ion beams ahead of supernova remnant shocks

The growth and saturation of Buneman-type instabilities is examined with a particle-in-cell (PIC) simulation for parameters that are representative for the foreshock region of fast supernova remnant shocks. A dense ion beam and the electrons correspond to the upstream plasma and a fast ion beam to the shock-reflected ions. The purpose of the 2D simulation is to identify the nonlinear saturation mechanisms, the electron heating and potential secondary instabilities that arise from anisotropic electron heating and result in the growth of magnetic fields. We confirm that the instabilities between both ion beams and the electrons saturate by the formation of phase space holes by the beam-aligned modes. The slower oblique modes accelerate some electrons, but they cannot heat up the electrons significantly before they are trapped by the faster beam-aligned modes. Two circular electron velocity distributions develop, which are centred around the velocity of each ion beam. They develop due to the scattering of the electrons by the electrostatic wave potentials. The growth of magnetic fields is observed, but their amplitude remains low.

A Physically-Constrained Source Model For FDTD Acoustic Simulation

The Finite Difference Time Domain (FDTD) method is becoming increasingly popular for room acoustics simulation. Yet, the literature on grid excitation methods is relatively sparse, and source functions are traditionally implemented in a hard or additive form
using arbitrarily-shaped functions which do not necessarily obey the physical laws of sound generation. In this paper we formulate
a source function based on a small pulsating sphere model. A physically plausible method to inject a source signal into the grid
is derived from first principles, resulting in a source with a near-flat spectrum that does not scatter incoming waves. In the final
discrete-time formulation, the source signal is the result of passing a Gaussian pulse through a digital filter simulating the dynamics of the pulsating sphere, hence facilitating a physically correct means to design source functions that generate a prescribed sound field.

Multidimensional modelling of X-ray spectra for AGN accretion disc outflows - III. Application to a hydrodynamical simulation

We perform multidimensional radiative transfer simulations to compute spectra for a hydrodynamical simulation of a line-driven accretion disc wind from an active galactic nucleus. The synthetic spectra confirm expectations from parametrized models that a disc wind can imprint a wide variety of spectroscopic signatures including narrow absorption lines, broad emission lines and a Compton hump. The formation of these features is complex with contributions originating from many of the different structures present in the hydrodynamical simulation. In particular, spectral features are shaped both by gas in a successfully launched outflow and in complex flows where material is lifted out of the disc plane but ultimately falls back. We also confirm that the strong Fe Ka line can develop a weak, red-skewed line wing as a result of Compton scattering in the outflow. In addition, we demonstrate that X-ray radiation scattered and reprocessed in the flow has a pivotal part in both the spectrum formation and determining the ionization conditions in the wind. We find that scattered radiation is rather effective in ionizing gas which is shielded from direct irradiation from the central source. This effect likely makes the successful launching of a massive disc wind somewhat more challenging and should be considered in future wind simulations. © 2010 The Authors. Journal compilation © 2010 RAS.

Molecular dynamics simulation model for the quantitative assessment of tool wear during single point diamond turning of cubic silicon carbide

Silicon carbide (SiC) is a material of great technological interest for engineering applications concerning hostile environments where silicon-based components cannot work (beyond 623 K). Single point diamond turning (SPDT) has remained a superior and viable method to harness process efficiency and freeform shapes on this harder material. However, it is extremely difficult to machine this ceramic consistently in the ductile regime due to sudden and rapid tool wear. It thus becomes non trivial to develop an accurate understanding of tool wear mechanism during SPDT of SiC in order to identify measures to suppress wear to minimize operational cost.

In this paper, molecular dynamics (MD) simulation has been deployed with a realistic analytical bond order potential (ABOP) formalism based potential energy function to understand tool wear mechanism during single point diamond turning of SiC. The most significant result was obtained using the radial distribution function which suggests graphitization of diamond tool during the machining process. This phenomenon occurs due to the abrasive processes between these two ultra hard materials. The abrasive action results in locally high temperature which compounds with the massive cutting forces leading to sp3–sp2 order–disorder transition of diamond tool. This represents the root cause of tool wear during SPDT operation of cubic SiC. Further testing led to the development of a novel method for quantitative assessment of the progression of diamond tool wear from MD simulations.

A Hamiltonian approach to simulation of acoustic systems involving nonlinear interactions

Nonlinear interactions take place in most systems that arise in music acoustics, usually as a result of player-instrument coupling. Several time-stepping methods exist for the numerical simulation of such systems. These methods generally involve the discretization of the Newtonian description of the system. However, it is not always possible to prove the stability of the resulting algorithms, especially when dealing with systems where the underlying force is a non-analytic function of the phase space variables. On the other hand, if the discretization is carried out on the Hamiltonian description of the system, it is possible to prove the stability of the derived numerical schemes. This Hamiltonian approach is applied to a series of test models of single or multiple nonlinear collisions and the energetic properties of the derived schemes are discussed. After establishing that the schemes respect the principle of conservation of energy, a nonlinear single-reed model is formulated and coupled to a digital bore, in order to synthesize clarinet-like sounds.

Defining Simulation Intent

Defining Simulation Intent involves capturing high level modelling and idealisation decisions in order to create an efficient and fit-for-purpose analysis. These decisions are recorded as attributes of the decomposed design space.

An approach to defining Simulation Intent is described utilising three known technologies: Cellular Modelling, the subdivision of space into volumes of simulation significance (structures, gas paths, internal and external airflows etc.); Equivalencing, maintaining a consistent and coherent description
of the equivalent representations of the spatial cells in different analysis models; and Virtual Topology, which offers tools for partitioning and de-partitioning the model without disturbing the manufacturing oriented design geometry. The end result is a convenient framework to which high level analysis attributes can be applied, and from which detailed analysis models can be generated
with a high degree of controllability, repeatability and automation. There are multiple novel aspects to the approach, including its reusability, robustness to changes in model topology and the inherent links created between analysis models at different levels of fidelity and physics.

By utilising Simulation Intent, CAD modelling for simulation can be fully exploited and simulation work-flows can be more readily automated, reducing many repetitive manual tasks (e.g. the definition of appropriate coupling between elements of different types and the application of boundary conditions). The approach has been implemented and tested with practical examples, and
significant benefits are demonstrated.