Transonic aeroelastic simulation for instability searches and uncertainty analysis

In this paper the use of eigenvalue stability analysis of very large dimension aeroelastic numerical models arising from the exploitation of computational fluid dynamics is reviewed. A formulation based on a block reduction of the system Jacobian proves powerful to allow various numerical algorithms to be exploited, including frequency domain solvers, reconstruction of a term describing the fluid–structure interaction from the sparse data which incurs the main computational cost, and sampling to place the expensive samples where they are most needed. The stability formulation also allows non-deterministic analysis to be carried out very efficiently through the use of an approximate Newton solver. Finally, the system eigenvectors are exploited to produce nonlinear and parameterised reduced order models for computing limit cycle responses. The performance of the methods is illustrated with results from a number of academic and large dimension aircraft test cases.

Mass media and polarisation processes in the bounded confidence model of opinion dynamics

This paper presents a social simulation in which we add an additional layer of mass media communication to the social network 'bounded confidence' model of Deffuant et al (2000). A population of agents on a lattice with continuous opinions and bounded confidence adjust their opinions on the basis of binary social network interactions between neighbours or communication with a fixed opinion. There are two mechanisms for interaction. 'Social interaction' occurs between neighbours on a lattice and 'mass communication' adjusts opinions based on an agent interacting with a fixed opinion. Two new variables are added, polarisation: the degree to which two mass media opinions differ, and broadcast ratio: the number of social interactions for each mass media communication. Four dynamical regimes are observed, fragmented, double extreme convergence, a state of persistent opinion exchange leading to single extreme convergence and a disordered state. Double extreme convergence is found where agents are less willing to change opinion and mass media communications are common or where there is moderate willingness to change opinion and a high frequency of mass media communications. Single extreme convergence is found where there is moderate willingness to change opinion and a lower frequency of mass media communication. A period of persistent opinion exchange precedes single extreme convergence, it is characterized by the formation of two opposing groups of opinion separated by a gradient of opinion exchange. With even very low frequencies of mass media communications this results in a move to central opinions followed by a global drift to one extreme as one of the opposing groups of opinion dominates. A similar pattern of findings is observed for Neumann and Moore neighbourhoods.

Fully kinetic simulation of ion acoustic and dust-ion acoustic waves

A series of numerical simulations is presented, based on a recurrence-free Vlasov kinetic model using kinetic phase point trajectories. All plasma components are modeled kinetically via a Vlasov evolution equation, then coupled through Poisson’s equation. The dynamics of ion acoustic waves in an electron-ion and in a dusty (electron-ion-dust) plasma configuration are investigated, focusing on wave decay due to Landau damping and, in particular, on the parametric dependence of the damping rate on the dust concentration and on the electron-to-ion temperature ratio. In the absence of dust, the occurrence of damping was observed, as expected, and its dependence to the relative magnitude of the electron vs ion temperature(s) was investigated. When present, the dust component influences the charge balance, enabling dust-ion acoustic waves to survive Landau damping even in the extreme regime where Te???Ti. The Landau damping rate is shown to be minimized for a strong dust concentration or/and for a high value of the electron-to-ion temperature ratio. Our results confirm earlier theoretical considerations and contribute to the interpretation of experimental observations of dust-ion acoustic wave characteristics.

Unsteady Gas Dynamic Engine Simulation

This paper provides an overview of the basic theory underlying 1D unsteady gas dynamics, the computational method developed at Queen’s University Belfast (QUB), the use of CFD as an alternative and some experimental results that demonstrate the techniques used to develop the mathematical models.

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.

Conformational Dynamics of the Human Propeller Telomeric DNA Quadruplex on a Microsecond Time Scale

The human telomeric DNA sequence with four repeats can fold into a parallel-stranded propeller-type topology. NMR structures solved under molecular crowding experiments correlate with the crystal structures found with crystal-packing interactions that are effectively equivalent to molecular crowding. This topology has been used for rationalization of ligand design and occurs experimentally in a number of complexes with a diversity of ligands, at least in the crystalline state. While G-quartet stems have been well characterised, the interactions of the TTA loop with the G-quartets are much less defined. To better understand the conformational variability and structural dynamics of the propeller-type topology, we performed molecular dynamics simulations in explicit solvent up to 1.5 µs. The analysis provides a detailed atomistic account of the dynamic nature of the TTA loops highlighting their interactions with the G-quartets including formation of an A:A base pair, triad, pentad and hexad. The results present a threshold in quadruplex simulations, with regards to understanding the flexible nature of the sugar-phosphate backbone in formation of unusual architecture within the topology. Furthermore, this study stresses the importance of simulation time in sampling conformational space for this topology.

Simulation of relativistically colliding laser-generated electron flows

The plasma dynamics resulting from the simultaneous impact, of two equal, ultra-intense laser pulses, in two spatially separated spots, onto a dense target is studied via particle-in-cell simulations. The simulations show that electrons accelerated to relativistic speeds cross the target and exit at its rear surface. Most energetic electrons are bound to the rear surface by the ambipolar electric field and expand along it. Their current is closed by a return current in the target, and this current configuration generates strong surface magnetic fields. The two electron sheaths collide at the midplane between the laser impact points. The magnetic repulsion between the counter-streaming electron beams separates them along the surface normal direction, before they can thermalize through other beam instabilities. This magnetic repulsion is also the driving mechanism for the beam-Weibel (filamentation) instability, which is thought to be responsible for magnetic field growth close to the internal shocks of gamma-ray burst jets. The relative strength of this repulsion compared to the competing electrostatic interactions, which is evidenced by the simulations, suggests that the filamentation instability can be examined in an experimental setting. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4768426]

Synthetic X-ray spectra for simulations of the dynamics of an accretion flow irradiated by a quasar

Ultraviolet and X-ray observations show evidence of outflowing gas around many active galactic nuclei. It has been proposed that some of these outflows are driven off gas infalling towards the central supermassive black hole. We perform radiative transfer calculations to compute the gas ionization state and the emergent X-ray spectra for both two- and three-dimensional (3D) hydrodynamical simulations of this outflow-from-inflow scenario. By comparison with observations, our results can be used to test the theoretical models and guide future numerical simulations. We predict both absorption and emission features, most of which are formed in a polar funnel of relatively dense (10 -10 g cm ) outflowing gas. This outflow causes strong absorption for observer orientation angles of ?35°. Particularly in 3D, the strength of this absorption varies significantly for different lines of sight owing to the fragmentary structure of the gas flow. Although infalling material occupies a large fraction of the simulation volume, we do not find that it imprints strong absorption features in the X-ray spectra since the ionization state is predicted to be very high. Thus, an absence of observed inflow absorption features does not exclude the models. The main spectroscopic consequence of the infalling gas is a Compton-scattered continuum component that partially re-fills the absorption features caused by the outflowing polar funnel. Fluorescence and scattering in the outflow are predicted to give rise to several emission features including a multicomponent Fe Ka emission complex for all observer orientations. For the hydrodynamical simulations considered, we predict both ionization states and column densities for the outflowing gas that are too high to be quantitatively consistent with well-observed X-ray absorption systems. Nevertheless, our results are qualitatively encouraging and further exploration of the model parameter space is warranted. Higher resolution hydrodynamic simulations are needed to determine whether the outflows fragment on scales unresolved in our current study, which may yield the denser lower ionization material that could reconcile the models and the observations. © 2012 The Authors Monthly Notices of the Royal Astronomical Society © 2012 RAS.

A self-learning algorithm for biased molecular dynamics

A new self-learning algorithm for accelerated dynamics, reconnaissance metadynamics, is proposed that is able to work with a very large number of collective coordinates. Acceleration of the dynamics is achieved by constructing a bias potential in terms of a patchwork of one-dimensional, locally valid collective coordinates. These collective coordinates are obtained from trajectory analyses so that they adapt to any new features encountered during the simulation. We show how this methodology can be used to enhance sampling in real chemical systems citing examples both from the physics of clusters and from the biological sciences.

High pressure CO2 absorption studies on imidazolium-based ionic liquids:Experimental and simulation approaches

A combined experimental-computational study on the CO absorption on 1-butyl-3-methylimidazolium hexafluophosphate, 1-ethyl-3-methylimidazolium bis[trifluoromethylsulfonyl]imide, and 1-butyl-3-methylimidazolium bis[trifluoromethylsulfonyl]imide ionic liquids is reported. The reported results allowed to infer a detailed nanoscopic vision of the absorption phenomena as a function of pressure and temperature. Absorption isotherms were measured at 318 and 338K for pressures up to 20MPa for ultrapure samples using a state-of-the-art magnetic suspension densimeter, for which measurement procedures are developed. A remarkable swelling effect upon CO absorption was observed for pressures higher than 10MPa, which was corrected using a method based on experimental volumetric data. The experimental data reported in this work are in good agreement with available literature isotherms. Soave-Redlich-Kwong and Peng-Robinson equations of state coupled with bi-parametric van der Waals mixing rule were used for successful correlations of experimental high pressure absorption data. Molecular dynamics results allowed to infer structural, energetic and dynamic properties of the studied CO+ionic liquids mixed fluids, showing the relevant role of the strength of anion-cation interactions on fluid volumetric properties and CO absorption. © 2012 Elsevier B.V.

Parameter estimation and equation formulation in Business Dynamics

System Dynamics enables modelling and simulation of highly non-linear feedback systems to predict future system behaviour. Parameter estimation and equation formulation are techniques in System Dynamics, used to retrieve the values of parameters or the equations for ?ows and/or variables. These techniques are crucial for the annotations and thereafter the simulation. This paper critically examines existing and well established approaches in parameter estimation and equation formulation along with their limitations, identifying performance gaps as well as providing directions for potential future research.

Dynamics of dark hollow Gaussian laser pulses in relativistic plasma

Optical beams with null central intensity have potential applications in the field of atom optics. The spatial and temporal evolution of a central shadow dark hollow Gaussian (DHG) relativistic laser pulse propagating in a plasma is studied in this article for first principles. A nonlinear Schrodinger-type equation is obtained for the beam spot profile and then solved numerically to investigate the pulse propagation characteristics. As series of numerical simulations are employed to trace the profile of the focused and compressed DHG laser pulse as it propagates through the plasma. The theoretical and simulation results predict that higher-order DHG pulses show smaller divergence as they propagate and, thus, lead to enhanced energy transport. © 2013 American Physical Society.

Calibrated CFD simulation to evaluate thermal comfort in a highly-glazed naturally ventilated room

Natural ventilation is a sustainable solution to maintaining healthy and comfortable environmental conditions in buildings. However, the effective design, construction and operation of naturally ventilated buildings require a good understanding of complex airflow patterns caused by the buoyancy and wind effects.The work presented in this article employed a 3D computational fluid dynamics (CFD) analysis in order to investigate environmental conditions and thermal comfort of the occupants of a highly-glazed naturally ventilated meeting room. This analysis was facilitated by the real-time field measurements performed in an operating building, and previously developed formal calibration methodology for reliable CFD models of indoor environments. Since, creating an accurate CFD model of an occupied space in a real-life scenario requires a high level of CFD expertise, trusted experimental data and an ability to interpret model input parameters; the calibration methodology guided towards a robust and reliable CFD model of the indoor environment. This calibrated CFD model was then used to investigate indoor environmental conditions and to evaluate thermal comfort indices for the occupants of the room. Thermal comfort expresses occupants' satisfaction with thermal environment in buildings by defining the range of indoor thermal environmental conditions acceptable to a majority of occupants. In this study, the thermal comfort analysis, supported by both field measurements and CFD simulation results, confirmed a satisfactory and optimal room operation in terms of thermal environment for the investigated real-life scenario. © 2013 Elsevier Ltd.

Applications of the generalized Langevin equation: Towards a realistic description of the baths

The generalized Langevin equation (GLE) method, as developed previously [L. Stella et al., Phys. Rev. B 89, 134303 (2014)], is used to calculate the dissipative dynamics of systems described at the atomic level. The GLE scheme goes beyond the commonly used bilinear coupling between the central system and the bath, and permits us to have a realistic description of both the dissipative central system and its surrounding bath. We show how to obtain the vibrational properties of a realistic bath and how to convey such properties into an extended Langevin dynamics by the use of the mapping of the bath vibrational properties onto a set of auxiliary variables. Our calculations for a model of a Lennard-Jones solid show that our GLE scheme provides a stable dynamics, with the dissipative/relaxation processes properly described. The total kinetic energy of the central system always thermalizes toward the expected bath temperature, with appropriate fluctuation around the mean value. More importantly, we obtain a velocity distribution for the individual atoms in the central system which follows the expected canonical distribution at the corresponding temperature. This confirms that both our GLE scheme and our mapping procedure onto an extended Langevin dynamics provide the correct thermostat. We also examined the velocity autocorrelation functions and compare our results with more conventional Langevin dynamics.

Ecological dynamics of extinct species in empty habitat networks. 2. The role of host plant dynamics

This paper explores the relative effects of host plant dynamics and butterfly-related parameters on butterfly persistence. It considers an empty habitat network where a rare butterfly (Cupido minimus) became extinct in 1939 in part of its historical range in north Wales, UK. Surviving populations of the butterfly in southern Britain were visited to assess use of its host plant (Anthyllis vulneraria) in order to calibrate habitat suitability and carrying capacity in the empty network in north Wales. These data were used to deduce that only a portion ( similar to 19%) of the host plant network from north Wales was likely to be highly suitable for oviposition. Nonetheless, roughly 65,460 eggs (3273 adult equivalents) could be expected to be laid in north Wales, were the empty network to be populated at the same levels as observed on comparable plants in surviving populations elsewhere. Simulated metapopulations of C. minimus in the empty network revealed that time to extinction and patch occupancy were significantly influenced by carrying capacity, butterfly mean dispersal distance and environmental stochasticity, although for most reasonable parameter values, the model system persisted. Simulation outputs differed greatly when host plant dynamics was incorporated into the modelled butterfly dynamics. Cupido minimus usually went extinct when host plant were at low densities. In these simulations host plant dynamics appeared to be the most important determinant of the butterfly's regional extirpation. Modelling the outcome of a reintroduction programme to C. minimus variation at high quality locations, revealed that 65% of systems survived at least 100 years. Given the current amount of resources of the north Wales landscape, the persistence of C. minimus under a realistic reintroduction programme has a good chance of being successful, if carried out in conjunction with a host plant management programme.