Simulation of a collisionless planar electrostatic shock in a proton-electron plasma with a strong initial thermal pressure change

The localized deposition of the energy of a laser pulse, as it ablates a solid target, introduces high thermal pressure gradients in the plasma. The thermal expansion of this laser-heated plasma into the ambient medium (ionized residual gas) triggers the formation of non-linear structures in the collisionless plasma. Here an electron-proton plasma is modelled with a particle-in-cell simulation to reproduce aspects of this plasma expansion. A jump is introduced in the thermal pressure of the plasma, across which the otherwise spatially uniform temperature and density change by a factor of 100. The electrons from the hot plasma expand into the cold one and the charge imbalance drags a beam of cold electrons into the hot plasma. This double layer reduces the electron temperature gradient. The presence of the low-pressure plasma modifies the proton dynamics compared with the plasma expansion into a vacuum. The jump in the thermal pressure develops into a primary shock. The fast protons, which move from the hot into the cold plasma in the form of a beam, give rise to the formation of phase space holes in the electron and proton distributions. The proton phase space holes develop into a secondary shock that thermalizes the beam.

PIC simulation of a thermal anisotropy-driven Weibel instability in a circular rarefaction wave

The expansion of an initially unmagnetized planar rarefaction wave has recently been shown to trigger a thermal anisotropy-driven Weibel instability (TAWI), which can generate magnetic fields from noise levels. It is examined here whether the TAWI can also grow in a curved rarefaction wave. The expansion of an initially unmagnetized circular plasma cloud, which consists of protons and hot electrons, into a vacuum is modelled for this purpose with a two-dimensional particle-in-cell (PIC) simulation. It is shown that the momentum transfer from the electrons to the radially accelerating protons can indeed trigger a TAWI. Radial current channels form and the aperiodic growth of a magnetowave is observed, which has a magnetic field that is oriented orthogonal to the simulation plane. The induced electric field implies that the electron density gradient is no longer parallel to the electric field. Evidence is presented here that this electric field modification triggers a environments, which are needed to explain the electromagnetic emissions by astrophysical jets. It is outlined how this instability could be examined experimentally.second magnetic instability, which results in a rotational low-frequency magnetowave. The relevance of the TAWI is discussed for the growth of small-scale magnetic fields in astrophysical

Multiphysics Simulation of Thermal Plasma Produced by Lightning Strike on Aircraft Structures

A novel numerical technique is proposed to model thermal plasma of microseconds/milliseconds time-scale effect. Modelling thermal plasma due to lightning strike will allow the estimation of electric current density, plasma pressure, and heat flux at the surface of the aircraft structure. These input data can then be used for better estimation of the mechanical/thermal induced damage on the aircraft structures for better protection systems design. Thermal plasma generated during laser cutting, electric (laser) welding and other plasma processing techniques have been the focus of many researchers. Thermal plasma is a gaseous state that consists from a mixture of electrons, ions, and natural particles. Thermal plasma can be assumed to be in local thermodynamic equilibrium, which means the electrons and the heavy species have equal temperature. Different numerical techniques have been developed using a coupled Navier Stokes – Heat transfer – Electromagnetic equations based on the assumption that the thermal plasma is a single laminar gas flow. These previous efforts focused on generating thermal plasma of time-scale in the range of seconds. Lighting strike on aircraft structures generates thermal plasma of time-scale of milliseconds/microseconds, which makes the previous physics used not applicable. The difficulty comes from the Navier-Stokes equations as the fluid is simulated under shock load, this introducing significant changes in the density and temperature of the fluid.

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.

Neutral and Charged 1-Butyl-3-methylimidazolium Triflate Clusters: Equilibrium Concentration in the Vapor Phase and Thermal Properties of Nanometric Droplets

Ground state energy, structure, and harmonic vibrational modes of 1-butyl-3-methylimidazolium triflate ([bmim][Tf]) clusters have been computed using an all-atom empirical potential model. Neutral and charged species have been considered up to a size (30 [bmim][Tf] pairs) well into the nanometric range. Free energy computations and thermodynamic modeling have been used to predict the equilibrium composition of the vapor phase as a function of temperature and density. The results point to a nonnegligible concentration of very small charged species at pressures (P ~ 0.01 Pa) and temperatures (T 600 K) at the boundary of the stability range of [bmim][Tf]. Thermal properties of nanometric neutral droplets have been investigated in the 0 T 700 K range. A near-continuous transition between a liquidlike phase at high T and a solidlike phase at low T takes place at T ~ 190 K in close correspondence with the bulk glass point Tg ~ 200 K. Solidification is accompanied by a transition in the dielectric properties of the droplet, giving rise to a small permanent dipole embedded into the solid cluster. The simulation results highlight the molecular precursors of several macroscopic properties and phenomena and point to the close competition of Coulomb and dispersion forces as their common origin.

A hierarchical physical model-based approach to predictive control of a thermal power plant for efficient plant-wide disturbance rejection

Extending the work presented in Prasad et al. (IEEE Proceedings on Control Theory and Applications, 147, 523-37, 2000), this paper reports a hierarchical nonlinear physical model-based control strategy to account for the problems arising due to complex dynamics of drum level and governor valve, and demonstrates its effectiveness in plant-wide disturbance handling. The strategy incorporates a two-level control structure consisting of lower-level conventional PI regulators and a higher-level nonlinear physical model predictive controller (NPMPC) for mainly set-point manoeuvring. The lower-level PI loops help stabilise the unstable drum-boiler dynamics and allow faster governor valve action for power and grid-frequency regulation. The higher-level NPMPC provides an optimal load demand (or set-point) transition by effective handling of plant-wide interactions and system disturbances. The strategy has been tested in a simulation of a 200-MW oil-fired power plant at Ballylumford in Northern Ireland. A novel approach is devized to test the disturbance rejection capability in severe operating conditions. Low frequency disturbances were created by making random changes in radiation heat flow on the boiler-side, while condenser vacuum was fluctuating in a random fashion on the turbine side. In order to simulate high-frequency disturbances, pulse-type load disturbances were made to strike at instants which are not an integral multiple of the NPMPC sampling period. Impressive results have been obtained during both types of system disturbances and extremely high rates of load changes, right across the operating range, These results compared favourably with those from a conventional state-space generalized predictive control (GPC) method designed under similar conditions.

Plant-wide predictive control for a thermal power plant based on a physical plant model

A constrained non-linear, physical model-based, predictive control (NPMPC) strategy is developed for improved plant-wide control of a thermal power plant. The strategy makes use of successive linearisation and recursive state estimation using extended Kalman filtering to obtain a linear state-space model. The linear model and a quadratic programming routine are used to design a constrained long-range predictive controller One special feature is the careful selection of a specific set of plant model parameters for online estimation, to account for time-varying system characteristics resulting from major system disturbances and ageing. These parameters act as nonstationary stochastic states and help to provide sufficient degrees-of-freedom to obtain unbiased estimates of controlled outputs. A 14th order non-linear plant model, simulating the dominant characteristics of a 200 MW oil-fired pou er plant has been used to test the NPMPC algorithm. The control strategy gives impressive simulation results, during large system disturbances and extremely high rate of load changes, right across the operating range. These results compare favourably to those obtained with the state-space GPC method designed under similar conditions.

A new model based approach for the prediction and optimisation of thermal homogeneity in single screw extrusion

In polymer extrusion, delivery of a melt which is homogenous in composition and temperature is important for good product quality. However, the process is inherently prone to temperature fluctuations which are difficult to monitor and control via single point based conventional thermo- couples. In this work, the die melt temperature profile was monitored by a thermocouple mesh and the data obtained was used to generate a model to predict the die melt temperature profile. A novel nonlinear model was then proposed which was demonstrated to be in good agreement with training and unseen data. Furthermore, the proposed model was used to select optimum process settings to achieve the desired average melt temperature across the die while improving the temperature homogeneity. The simulation results indicate a reduction in melt temperature variations of up to 60%.

Autoclave Cure Simulation of Composite Structures Applying Implicit and Explicit FE

Simulation of the autoclave manufacturing technique of composites can yield a preliminary estimation of induced residual thermal stresses and deformations that affect component fatigue life, and required tolerances for assembly. In this paper, an approach is proposed to simulate the autoclave manufacturing technique for unidirectional composites. The proposed approach consists of three modules. The first module is a Thermo-chemical model to estimate the temperature and the degree of cure distributions in the composite part during the cure cycle. The second and third modules are a sequential stress analysis using FE-Implicit and FE-Explicit respectively. User-material subroutine is used to model the Viscoelastic properties of the material based on theory of micromechanics.

Numerical simulation of a patent technology for sealing of deep-sea oil wells using nonlinear finite element method

For over 50 years bridge plugs and cement have been used for well abandonment and work over and are still the material of choice. However the failures of cement abandonments using bridge plugs has been reported on many occasions, some of which have resulted in fatal consequences. A new patented product is designed to address the shortcomings associated with using bridge plugs and cement. The new developed tools use an alloy based on bismuth that is melted in situ using Thermite reaction. The tool uses the expansion properties of bismuth to seal the well. Testing the new technology in real field under more than 2 km deep sea water can be expensive. Virtual simulation of the new device under simulated thermal and mechanical environment can be achieved using nonlinear finite element method to validate the product and reduce cost. Experimental testing in the lab is performed to measure heat generated due to thermite reaction. Then, a sequential thermal mechanical explicit/implicit finite element solver is used to simulate the device under both testing lab and deep water conditions.

Atomistic investigation on the structure-property relationship during thermal spray nanoparticle impact

During thermal spraying, hot particles impact on a colder substrate. This interaction of crystalline copper nanoparticles and copper substrate is modelled, using MD simulation. The quantitative results of the impacts at different velocities and temperatures are evaluated using a newly defined flattening aspect ratio. This ratio between the maximum diameter after the impact and the height of the splat increases with increasing Reynolds numbers until a critical value is reached. At higher Reynolds numbers the flattening aspect ratio decreases again, as the kinetic energy of the particle leads to increasing substrate temperature and, therefore, decreases the substrate resistance. Thus, the particle penetrates into the substrate and deforms less.