Impinging jet simulation of stationary downburst flow over topography

A non-translating, long duration thunderstorm downburst has been simulated experimentally and numerically by modelling a spatially stationary steady flow impinging air jet. Velocity profiles were shown to compare well with an upper-bound of velocity measurements reported for full-scale microbursts. Velocity speed-up over a range of topographic features in simulated downburst flow was also tested with comparisons made to previous work in a similar flow, and also boundary layer wind tunnel experiments. It was found that the amplification measured above the crest of topographic features in simulated downburst flow was up to 35% less than that observed in boundary layer flow for all shapes tested. From the computational standpoint we conclude that the Shear Stress Transport (SST) model performs the best from amongst a range of eddy-viscosity and second moment closures tested for modelling the impinging jet flow.

Pulsed wall jet simulation of a stationary thunderstorm downburst, Part A: Physical structure and flow field characterization

A pulsed wall jet has been used to simulate the gust front of a thunderstorm downburst. Flow visualization, wind speed and surface pressure measurements were obtained. The characteristics of the hypothesized ring vortex of a full-scale downburst were reproduced at a scale estimated to be 1:3000.

Numerical simulation of idealised three-dimensional downburst wind fields

Convective downburst wind storms generate the peak annual gust wind speed for many parts of the non-cyclonic world at return periods of importance for ultimate limit state design. Despite this there is little clear understanding of how to appropriately design for these wind events given their significant dissimilarities to boundary layer winds upon which most design is based. To enhance the understanding of wind fields associated with these storms a three-dimensional numerical model was developed to simulate a multitude of idealised downburst scenarios and to investigate their near-ground wind characteristics. Stationary and translating downdraft wind events in still and sheared environments were simulated with baseline results showing good agreement with previous numerical work and full-scale observational data. Significant differences are shown in the normalised peak wind speed velocity profiles depending on the environmental wind conditions in the vicinity of the simulated event. When integrated over the height of mid- to high rise structures, all simulated profiles are shown to produce wind loads smaller than an equivalent 10 m height matched open terrain boundary layer profile. This suggests that for these structures the current design approach is conservative from an ultimate loading standpoint. Investigating the influence of topography on the structure of the simulated near-ground downburst wind fields, it is shown that these features amplify wind speeds in a manner similar to that expected for boundary layer winds, but the extent of amplification is reduced. The level of reduction is shown to be dependent on the depth of the simulated downburst outflow.

Computer simulation of scaffold degradation

Scaffolds are porous biocompatible materials with suitable microarchitectures that are designed to allow for cell adhesion, growth and proliferation. They are used in combination with cells in regenerative medicine to promote tissue regeneration by means of a controlled deposition of natural extracellular matrix by the hosted cells therein. This healing process is in many cases accompanied by scaffold degradation up to its total disappearance when the scaffold is made of a biodegradable material. This work presents a computational model that simulates the degradation of scaffolds. The model works with three-dimensional microstructures, which have been previously discretised into small cubic homogeneous elements, called voxels. The model simulates the evolution of the degradation of the scaffold using a Monte Carlo algorithm, which takes into account the curvature of the surface of the fibres. The simulation results obtained in this study are in good agreement with empirical degradation measurements performed by mass loss on scaffolds after exposure to an etching alkaline solution.

Simulation of plant cell shrinkage during drying : A SPH–DEM approach

Plant based dried food products are popular commodities in global market where much research is focused to improve the products and processing techniques. In this regard, numerical modelling is highly applicable and in this work, a coupled meshfree particle-based two-dimensional (2-D) model was developed to simulate micro-scale deformations of plant cells during drying. Smoothed Particle Hydrodynamics (SPH) was used to model the viscous cell protoplasm (cell fluid) by approximating it to an incompressible Newtonian fluid. The visco-elastic characteristic of the cell wall was approximated to a Neo-Hookean solid material augmented with a viscous term and modelled with a Discrete Element Method (DEM). Compared to a previous work [H. C. P. Karunasena, W. Senadeera, Y. T. Gu and R. J. Brown, Appl. Math. Model., 2014], this study proposes three model improvements: linearly decreasing positive cell turgor pressure during drying, cell wall contraction forces and cell wall drying. The improvements made the model more comparable with experimental findings on dried cell morphology and geometric properties such as cell area, diameter, perimeter, roundness, elongation and compactness. This single cell model could be used as a building block for advanced tissue models which are highly applicable for product and process optimizations in Food Engineering.

STI-GMaS : an open-source environment for simulation of sexually-transmitted infections

Background Sexually-transmitted pathogens often have severe reproductive health implications if treatment is delayed or absent, especially in females. The complex processes of disease progression, namely replication and ascension of the infection through the genital tract, span both extracellular and intracellular physiological scales, and in females can vary over the distinct phases of the menstrual cycle. The complexity of these processes, coupled with the common impossibility of obtaining comprehensive and sequential clinical data from individual human patients, makes mathematical and computational modelling valuable tools in developing our understanding of the infection, with a view to identifying new interventions. While many within-host models of sexually-transmitted infections (STIs) are available in existing literature, these models are difficult to deploy in clinical/experimental settings since simulations often require complex computational approaches. Results We present STI-GMaS (Sexually-Transmitted Infections – Graphical Modelling and Simulation), an environment for simulation of STI models, with a view to stimulating the uptake of these models within the laboratory or clinic. The software currently focuses upon the representative case-study of Chlamydia trachomatis, the most common sexually-transmitted bacterial pathogen of humans. Here, we demonstrate the use of a hybrid PDE–cellular automata model for simulation of a hypothetical Chlamydia vaccination, demonstrating the effect of a vaccine-induced antibody in preventing the infection from ascending to above the cervix. This example illustrates the ease with which existing models can be adapted to describe new studies, and its careful parameterisation within STI-GMaS facilitates future tuning to experimental data as they arise. Conclusions STI-GMaS represents the first software designed explicitly for in-silico simulation of STI models by non-theoreticians, thus presenting a novel route to bridging the gap between computational and clinical/experimental disciplines. With the propensity for model reuse and extension, there is much scope within STI-GMaS to allow clinical and experimental studies to inform model inputs and drive future model development. Many of the modelling paradigms and software design principles deployed to date transfer readily to other STIs, both bacterial and viral; forthcoming releases of STI-GMaS will extend the software to incorporate a more diverse range of infections.

Idealised simulation of downburst winds

Convectively driven downburst winds pose a threat to structures and communities in many regions of Australia not subject to tropical cyclones. Extreme value analysis shows that for return periods of interest to engineering design these events produce higher gust wind speeds than synoptic scale windstorms. Despite this, comparatively little is known of the near ground wind structure of these potentially hazardous windstorms. With this in mind, a series of idealised three-dimensional numerical simulations were undertaken to investigate convective storm wind fields. A dry, non-hydrostatic, sub-cloud model with parameterisation of the microphysics was used. Simulations were run with a uniform 20 m horizontal grid resolution and a variable vertical resolution increasing from 1 m. A systematic grid resolution study showed further refinement did not alter the morphological structure of the outflow. Simulations were performed for stationary downbursts in a quiescent air field, stationary downbursts embedded within environmental boundary layer winds, and also translating downbursts embedded within environmental boundary layer winds.

Simulation of carbon arc discharge for the synthesis of nanotubes

Arc discharge ablation with a catalyst-filled carbon anode in a helium background was used for the synthesis of graphene and carbon nanotubes. In this paper, we present the results of the numerical simulation of the distribution of various plasma parameters in discharge, as well as the distribution of carbon flux on the nanotube surface, for the typical discharge with an arc current of 60 A and a background gas pressure of 68 kPa.

Diagnostics and two-dimensional simulation of low-frequency inductively coupled plasmas with neutral gas heating and electron heat fluxes

This article presents the results on the diagnostics and numerical modeling of low-frequency (∼460 KHz) inductively coupled plasmas generated in a cylindrical metal chamber by an external flat spiral coil. Experimental data on the electron number densities and temperatures, electron energy distribution functions, and optical emission intensities of the abundant plasma species in low/intermediate pressure argon discharges are included. The spatial profiles of the plasma density, electron temperature, and excited argon species are computed, for different rf powers and working gas pressures, using the two-dimensional fluid approach. The model allows one to achieve a reasonable agreement between the computed and experimental data. The effect of the neutral gas temperature on the plasma parameters is also investigated. It is shown that neutral gas heating (at rf powers≥0.55kW) is one of the key factors that control the electron number density and temperature. The dependence of the average rf power loss, per electron-ion pair created, on the working gas pressure shows that the electron heat flux to the walls appears to be a critical factor in the total power loss in the discharge.

Tailoring the composition of self-assembled Si1−xCx quantum dots : simulation of plasma/ion-related controls

Precise control of composition and internal structure is essential for a variety of novel technological applications which require highly tailored binary quantum dots (QDs) with predictable optoelectronic and mechanical properties. The delicate balancing act between incoming flux and substrate temperature required for the growth of compositionally graded (Si1-xC x; x varies throughout the internal structure), core-multishell (discrete shells of Si and C or combinations thereof) and selected composition (x set) QDs on low-temperature plasma/ion-flux-exposed Si(100) surfaces is investigated via a hybrid numerical simulation. Incident Si and C ions lead to localized substrate heating and a reduction in surface diffusion activation energy. It is shown that by incorporating ions in the influx, a steady-state composition is reached more quickly (for selected composition QDs) and the composition gradient of a Si1-xCx QD may be fine tuned; additionally (with other deposition conditions remaining the same), larger QDs are obtained on average. It is suggested that ionizing a portion of the influx is another way to control the average size of the QDs, and ultimately, their internal structure. Advantages that can be gained by utilizing plasma/ion-related controls to facilitate the growth of highly tailored, compositionally controlled quantum dots are discussed as well.

Nanosphere monolayer-templated, ion-assisted nanofeature etching in dielectric materials : a numerical simulation of nanoscale ion flux topography

The results of numerical simulations of nanometer precision distributions of microscopic ion fluxes in ion-assisted etching of nanoscale features on the surfaces of dielectric materials using a self-assembled monolayer of spherical nanoparticles as a mask are presented. It is shown that the ion fluxes to the substrate and nanosphere surfaces can be effectively controlled by the plasma parameters and the external bias applied to the substrate. By proper adjustment of these parameters, the ion flux can be focused onto the areas uncovered by the nanospheres. Under certain conditions, the ion flux distributions feature sophisticated hexagonal patterns, which may lead to very different nanofeature etching profiles. The results presented are generic and suggest viable ways to overcome some of the limitations of the existing plasma-assisted nanolithography.

Multi-scale hybrid numerical simulation of the growth of high-aspect-ratio nanostructures

Recently, a variety high-aspect-ratio nanostructures have been grown and profiled for various applications ranging from field emission transistors to gene/drug delivery devices. However, fabricating and processing arrays of these structures and determining how changing certain physical parameters affects the final outcome is quite challenging. We have developed several modules that can be used to simulate the processes of various physical vapour deposition systems from precursor interaction in the gas phase to gas-surface interactions and surface processes. In this paper, multi-scale hybrid numerical simulations are used to study how low-temperature non-equilibrium plasmas can be employed in the processing of high-aspect-ratio structures such that the resulting nanostructures have properties suitable for their eventual device application. We show that whilst using plasma techniques is beneficial in many nanofabrication processes, it is especially useful in making dense arrays of high-aspect-ratio nanostructures.

Simulation of gas-phase nanoparticle dynamics in the plasma-enhanced chemical vapor deposition of carbon nanostructures

The results of 1D simulation of nanoparticle dynamics in the areas adjacent to nanostructured carbon-based films exposed to chemically active complex plasma of CH4 + H2 + Ar gas mixtures are presented. The nanoparticle-loaded near-substrate (including sheath and presheath) areas of a low-frequency (0.5 MHz) inductively coupled plasma facility for the PECVD growth of the ordered carbon-based nanotip structures are considered. The conditions allowing one to predict the size of particles that can pass through the plasma sheath and softly land onto the surface are formulated. The possibility of soft nano-cluster deposition without any additional acceleration common for some existing nano-cluster deposition schemes is demonstrated. The effect of the substrate heating power and the average atomic mass of neutral species is studied numerically and verified experimentally.