Templated i-PVD of metallic nanodot arrays

The results of numerical simulation of plasma-based, porous, template-assisted nanofabrication of Au nanodot arrays on highly-doped silicon taking into account typical electron density of low-temperature plasma of 1017-1018 m-3 and electron temperature of 2-5 eV are reported here. Three-dimensional microscopic topography of ion flux distribution over the outer and inner surfaces of the nanoporous template is obtained via numerical simulation of Au ion trajectories in the plasma sheath, in the close proximity of, and inside the nanopores. It is shown that, by manipulating the electron temperature, the cross-sheath potential drop, and by additionally altering the structure of the nanoporous template, one can control the ion fluxes within the nanopores, and eventually maximize the ion deposition onto the top surface of the developing crystalline Au nanodots (see top panel in the figure). In the same time, this procedure allows one to minimize amorphous deposits on the sidewalls that clutter and may eventually close the nanopores, thus disrupting the nanodot growth process, as it is shown in the bottom panel in the figure on the right.

Plasma-aided nanofabrication: "plasma-building block" approach

Unique features and benefits of the plasma-aided nanofabrication are considered by using the "plasma-building block" approach, which is based on plasma diagnostics and nanofilm characterization, cross-referenced by numerical simulation of generation and dynamics of building blocks in the gas phase, their interaction with nanostructured surfaces, and ab initio simulation of chemical structure of relevant nanoassemblies. The examples include carbon nanotip microemitter structures, semiconductor quantum dots and nanowires synthesized in the integrated plasma-aided nanofabrication facility.

Nanofabrication of single-crystalline flat-panel display microemitters : a plasma-building unit approach

This contribution is focused on plasma-enhanced chemical vapor deposition systems and their unique features that make them particularly attractive for nanofabrication of flat panel display microemitter arrays based on ordered patterns of single-crystalline carbon nanotip structures. The fundamentals of the plasma-based nanofabrication of carbon nanotips and some other important nanofilms and nanostructures are examined. Specific features, challenges, and potential benefits of using the plasma-based systems for relevant nanofabrication processes are analyzed within the framework of the "plasma-building unit" approach that builds up on extensive experimental data on plasma diagnostics and nanofilm/nanostructure characterization, and numerical simulation of the species composition in the ionized gas phase (multicomponent fluid models), ion dynamics and interaction with ordered carbon nanotip patterns, and ab initio computations of chemical structure of single crystalline carbon nanotips. This generic approach is also applicable for nanoscale assembly of various carbon nanostructures, semiconductor quantum dot structures, and nano-crystalline bioceramics. Special attention is paid to most efficient control strategies of the main plasma-generated building units both in the ionized gas phase and on nanostructured deposition surfaces. The issues of tailoring the reactive plasma environments and development of versatile plasma nanofabrication facilities are also discussed.

CFD simulations of flow and heat transfer through the porous interface of a metal foam heat exchanger

This paper offers numerical modelling of a waste heat recovery system. A thin layer of metal foam is attached to a cold plate to absorb heat from hot gases leaving the system. The heat transferred from the exhaust gas is then transferred to a cold liquid flowing in a secondary loop. Two different foam PPI (Pores Per Inch) values are examined over a range of fluid velocities. Numerical results are then compared to both experimental data and theoretical results available in the literature. Challenges in getting the simulation results to match those of the experiments are addressed and discussed in detail. In particular, interface boundary conditions specified between a porous layer and a fluid layer are investigated. While physically one expects much lower fluid velocity in the pores compared to that of free flow, capturing this sharp gradient at the interface can add to the difficulties of numerical simulation. The existing models in the literature are modified by considering the pressure gradient inside and outside the foam. Comparisons against the numerical modelling are presented. Finally, based on experimentally-validated numerical results, thermo-hydraulic performance of foam heat exchangers as waste heat recovery units is discussed with the main goal of reducing the excess pressure drop and maximising the amount of heat that can be recovered from the hot gas stream.

Including nonequilibrium interface kinetics in a continuum model for melting nanoscaled particles

The melting temperature of a nanoscaled particle is known to decrease as the curvature of the solid-melt interface increases. This relationship is most often modelled by a Gibbs--Thomson law, with the decrease in melting temperature proposed to be a product of the curvature of the solid-melt interface and the surface tension. Such a law must break down for sufficiently small particles, since the curvature becomes singular in the limit that the particle radius vanishes. Furthermore, the use of this law as a boundary condition for a Stefan-type continuum model is problematic because it leads to a physically unrealistic form of mathematical blow-up at a finite particle radius. By numerical simulation, we show that the inclusion of nonequilibrium interface kinetics in the Gibbs--Thomson law regularises the continuum model, so that the mathematical blow up is suppressed. As a result, the solution continues until complete melting, and the corresponding melting temperature remains finite for all time. The results of the adjusted model are consistent with experimental findings of abrupt melting of nanoscaled particles. This small-particle regime appears to be closely related to the problem of melting a superheated particle.

Deformation of a single red blood cell in a microvessel

Red blood cells (RBCs) are the most common type of cells in human blood and they exhibit different types of motions and deformed shapes in capillary flows. The behaviour of the RBCs should be studied in order to explain the RBC motion and deformation mechanism. This article presents a numerical simulation method for RBC deformation in microvessels. A two dimensional spring network model is used to represent the RBC membrane, where the elastic stretch/compression energy and the bending energy are considered with the constraint of constant RBC surface area. The forces acting on the RBC membrane are obtained from the principle of virtual work. The whole fluid domain is discretized into a finite number of particles using smoothed particle hydrodynamics concepts and the motions of all the particles are solved using Navier--Stokes equations. Minimum energy concepts are used to simulate the deformed shape of the RBC model. To verify the model, the motion of a single RBC is simulated in a Poiseuille flow and the characteristic parachute shape of the RBC is observed. Further simulations reveal that the RBC shows a tank treading motion when it flows in a linear shear flow.

Effects of strain artefacts arising from a pre-defined callus domain in models of bone healing mechanobiology

Iterative computational models have been used to investigate the regulation of bone fracture healing by local mechanical conditions. Although their predictions replicate some mechanical responses and histological features, they do not typically reproduce the predominantly radial hard callus growth pattern observed in larger mammals. We hypothesised that this discrepancy results from an artefact of the models’ initial geometry. Using axisymmetric finite element models, we demonstrated that pre-defining a field of soft tissue in which callus may develop introduces high deviatoric strains in the periosteal region adjacent to the fracture. These bone-inhibiting strains are not present when the initial soft tissue is confined to a thin periosteal layer. As observed in previous healing models, tissue differentiation algorithms regulated by deviatoric strain predicted hard callus forming remotely and growing towards the fracture. While dilatational strain regulation allowed early bone formation closer to the fracture, hard callus still formed initially over a broad area, rather than expanding over time. Modelling callus growth from a thin periosteal layer successfully predicted the initiation of hard callus growth close to the fracture site. However, these models were still susceptible to elevated deviatoric strains in the soft tissues at the edge of the hard callus. Our study highlights the importance of the initial soft tissue geometry used for finite element models of fracture healing. If this cannot be defined accurately, alternative mechanisms for the prediction of early callus development should be investigated.

A meshless method based on Point Interpolation Method (PIM) for the space fractional diffusion equation

This paper aims to develop a meshless approach based on the Point Interpolation Method (PIM) for numerical simulation of a space fractional diffusion equation. Two fully-discrete schemes for the one-dimensional space fractional diffusion equation are obtained by using the PIM and the strong-forms of the space diffusion equation. Numerical examples with different nodal distributions are studied to validate and investigate the accuracy and efficiency of the newly developed meshless approach.

Structural damage alarming and localization of cable-supported bridges using multi-novelty indices: A feasibility study

This paper presents a feasibility study on structural damage alarming and localization of long-span cable-supported bridges using multi-novelty indices formulated by monitoring-derived modal parameters. The proposed method which requires neither structural model nor damage model is applicable to structures of arbitrary complexity. With the intention to enhance the tolerance to measurement noise/uncertainty and the sensitivity to structural damage, an improved novelty index is formulated in terms of auto-associative neural networks (ANNs) where the output vector is designated to differ from the input vector while the training of the ANNs needs only the measured modal properties of the intact structure under in-service conditions. After validating the enhanced capability of the improved novelty index for structural damage alarming over the commonly configured novelty index, the performance of the improved novelty index for damage occurrence detection of large-scale bridges is examined through numerical simulation studies of the suspension Tsing Ma Bridge (TMB) and the cable-stayed Ting Kau Bridge (TKB) incurred with different types of structural damage. Then the improved novelty index is extended to formulate multi-novelty indices in terms of the measured modal frequencies and incomplete modeshape components for damage region identification. The capability of the formulated multi-novelty indices for damage region identification is also examined through numerical simulations of the TMB and TKB.

Structural optimization approach for specially shaped composite tank in spacecrafts

This study proposes an optimized approach of designing in which a model specially shaped composite tank for spacecrafts is built by applying finite element analysis. The composite layers are preliminarily designed by combining quasi-network design method with numerical simulation, which determines the ratio between the angle and the thickness of layers as the initial value of the optimized design. By adopting an adaptive simulated annealing algorithm, the angles and the numbers of layers at each angle are optimized to minimize the weight of structure. Based on this, the stacking sequence of composite layers is formulated according to the number of layers in the optimized structure by applying the enumeration method and combining the general design parameters. Numerical simulation is finally adopted to calculate the buckling limit of tanks in different designing methods. This study takes a composite tank with a cone-shaped cylinder body as example, in which ellipsoid head section and outer wall plate are selected as the object to validate this method. The result shows that the quasi-network design method can improve the design quality of composite material layer in tanks with complex preliminarily loading conditions. The adaptive simulated annealing algorithm can reduce the initial design weight by 30%, which effectively probes the global optimal solution and optimizes the weight of structure. It can be therefore proved that, this optimization method is capable of designing and optimizing specially shaped composite tanks with complex loading conditions.

Molecular orientational order of nitroxide radicals in liquid crystalline media

The orientational distribution of a set of stable nitroxide radicals in aligned liquid crystals 5CB (nematic) and 8CB (smectic A) was studied in detail by numerical simulation of EPR spectra. The order parameters up to the 10th rank were measured. The directions of the principal orientation axes of the radicals were determined. It was shown that the ordering of the probe molecules is controlled by their interaction with the matrix molecules more than the inherent geometry of the probes themselves. The rigid fused phenanthrene-based (A5) and 2-azaphenalene (A4) nitroxides as well as the rigid core elongated C11 and 5α-cholestane (CLS) nitroxides were found to be most sensitive to the orientation of the liquid crystal matrixes.

Modelling of the evaporation of a droplet suspended in a binary atmosphere

The process of spray drying is applied in a number of contexts. One such application is the production of a synthetic rock used for storage of nuclear waste. To establish a framework for a model of the spray drying process for this application, we here develop a model describing evaporation from droplets of pure water, such that the model may be extended to account for the presence of colloid within the droplet. We develop a spherically-symmetric model and formulate continuum equations describing mass, momentum, and energy balance in both the liquid and gas phases from first principles. We establish appropriate boundary conditions at the surface of the droplet, including a generalised Clapeyron equation that accurately describes the temperature at the surface of the droplet. To account for experiment design, we introduce a simplified platinum ball and wire model into the system using a thin wire problem. The resulting system of equations is transformed in order to simplify a finite volume solution scheme. The results from numerical simulation are compared with data collected for validation, and the sensitivity of the model to variations in key parameters, and to the use of Clausius–Clapeyron and generalised Clapeyron equations, is investigated. Good agreement is found between the model and experimental data, despite the simplicity of the platinum phase model.

On the sensitivity of elastic waves due to structural damages:Time-frequency based indexing method

Time-frequency analysis of various simulated and experimental signals due to elastic wave scattering from damage are performed using wavelet transform (WT) and Hilbert-Huang transform (HHT) and their performances are compared in context of quantifying the damages. Spectral finite element method is employed for numerical simulation of wave scattering. An analytical study is carried out to study the effects of higher-order damage parameters on the reflected wave from a damage. Based on this study, error bounds are computed for the signals in the spectral and also on the time-frequency domains. It is shown how such an error bound can provide all estimate of error in the modelling of wave propagation in structure with damage. Measures of damage based on WT and HHT is derived to quantify the damage information hidden in the signal. The aim of this study is to obtain detailed insights into the problem of (1) identifying localised damages (2) dispersion of multifrequency non-stationary signals after they interact with various types of damage and (3) quantifying the damages. Sensitivity analysis of the signal due to scattered wave based on time-frequency representation helps to correlate the variation of damage index measures with respect to the damage parameters like damage size and material degradation factors.

Analytical solution of unsteady three-dimensional MHD boundary layer flow and heat transfer due to impulsively stretched plane surface

The unsteady magnetohydrodynamic viscous flow and heat transfer of Newtonian fluids induced by an impulsively stretched plane surface in two lateral directions are studied by using an analytic technique, namely, the homotopy method. The analytic series solution presented here is highly accurate and uniformly valid for all time in the entire region. The effects of the stretching ratio and the magnetic field on the surface shear stresses and heat transfer are studied. The surface shear stresses in x- and y-directions and the surface heat transfer are enchanced by increasing stretching ratio for a fixed value of the magnetic parameter. For a fixed stretching ratio, the surface shear stresses increase with the magnetic parameter, but the heat transfer decreases. The Nusselt number takes longer time to reach the steady state than the skin friction coefficients. There is a smooth transition from the initial unsteady state to the steady state.

Smooth stabilisation of nonholonomic robots subject to disturbances

In this paper, we address the problem of stabilisation of robots subject to nonholonommic constraints and external disturbances using port-Hamiltonian theory and smooth time-invariant control laws. This should be contrasted with the commonly used switched or time-varying laws. We propose a control design that provides asymptotic stability of an manifold (also called relative equilibria)-due to the Brockett condition this is the only type of stabilisation possible using smooth time-invariant control laws. The equilibrium manifold can be shaped to certain extent to satisfy specific control objectives. The proposed control law also incorporates integral action, and thus the closed-loop system is robust to unknown constant disturbances. A key step in the proposed design is a change of coordinates not only in the momentum, but also in the position vector, which differs from coordinate transformations previously proposed in the literature for the control of nonholonomic systems. The theoretical properties of the control law are verified via numerical simulation based on a robotic ground vehicle model with differential traction wheels and non co-axial centre of mass and point of contact.