Multifunctional ToF-SIMS: combinatorial mapping of gradient energy substrates

We present a simple method for chemical modification of chlorosilane self-assembled monolayers (SAMs) on Si surfaces by exposure to a gradient of UV-ozone radiation to create stable substrates with a range of contact angles (θ_H2O≈5–95°) and surface energies on a single substrate. These gradient energy substrates are developed to potentially generate libraries for combinatorial studies of thin film phenomenology, where a systematic variation of interfacial surface energy represents one of the significant parameters along one axis. The graded oxidation process presents a systematic variation of surface chemical composition. We have utilized contact angle measurements and time-of-flight secondary ion mass spectrometry (ToF-SIMS) to investigate this variation for a series of ions, among which are SiCH₃⁺, SiOH⁺ and COOH⁻. We show that the macroscopic measurements of surface free energy/contact angle correlate with the detailed analysis of surface chemistry (as assessed by ToF-SIMS) on these test substrates.     

Simulation of a turbulent premixed open V-shaped flame using contour advection with surgery

Despite its capability of high spatial resolution, simulation of turbulent flows with traditional Lagrangian (front tracking) scheme is often discouraged by numerical instability caused by clustering of marker nodes and topological changes of fronts. Contour advection surgery (CAS), being a robust front tracking scheme, can limit the growth of front complexity during simulation without jeopardizing accuracy or efficiency. This allows it to open up an advantage over traditional front-tracking schemes. It has already been demonstrated that CAS, with incorporation of the reaction sheet model, can accurately simulate the propagation and advection of a turbulent premixed V-shaped flame. In this study, it is further tested with 10 values of vortex circulation. A range of upstream turbulence levels of 1.8–19.8% was obtained. Results indicate that upstream turbulence increase the average flame length, flame zone area and the overall burning rate. Flame surface density Σ was also estimated. Maximum values of Σ obtained lie in the range 0.1–1.4 mm⁻¹. Skewness towards the burnt region was observed in all profiles of Σ. Similar to results from laboratory experiments, it was found that Σ values decreases with upstream turbulence. From this study, the ability of CAS to cope with intense turbulence is demonstrated and a better quantitative understanding on the scheme has also been acquired.     

Tight-binding method modelling of photonic crystal waveguides

We present here a tight-binding-like modelling of two-dimensional (2D) photonic crystals (PCs). Adopted from solid-state physics, the concept of generalized Wannier functions is used to construct a localized state basis that allows a parameter-free ab initio study of defects in PCs. We demonstrate here for a 2D triangular lattice of dielectric rods in air, the existence of this localized basis and the possibility to study large scale complex dielectric structures deviating from periodicity. Specific numerical simulations on a split waveguide embedded in this triangular lattice are performed, and they demonstrate the superiority of this method over plane wave based techniques.     

On a nonlinear multigrid algorithm with primal relaxation for the image total variation minimisation

Digital image restoration has drawn much attention in the recent years and a lot of research has been done on effective variational partial differential equation models and their theoretical studies. However there remains an urgent need to develop fast and robust iterative solvers, as the underlying problem sizes are large. This paper proposes a fast multigrid method using primal relaxations. The basic primal relaxation is known to get stuck at a ‘local’ non-stationary minimum of the solution, which is usually believed to be ‘non-smooth’. Our idea is to utilize coarse level corrections, overcoming the deadlock of a basic primal relaxation scheme. A further refinement is to allow non-regular coarse levels to correct the solution, which helps to improve the multilevel method. Numerical experiments on both 1D and 2D images are presented.     

3D molecular imaging SIMS

Thin monolayer and bilayer films of spin cast poly(methyl methacrylate) (PMMA), poly(2-hydroxyethyl methacrylate) (PHEMA), poly(lactic) acid (PLA) and PLA doped with several pharmaceuticals have been analyzed by dynamic SIMS using SF₅⁺ polyatomic primary ion bombardment. Each of these systems exhibited minimal primary beam-induced degradation under cluster ion bombardment allowing molecular depth profiles to be obtained through the film. By combing secondary ion imaging with depth profiling, three-dimensional molecular image depth profiles have been obtained from these systems. In another approach, bevel cross-sections are cut in the samples with the SF₅⁺ primary ion beam to produce a laterally magnified cross-section of the sample that does not contain the beam-induced damage that would be induced by conventional focussed ion beam (FIB) cross-sectioning. The bevel surface can then be examined using cluster SIMS imaging or other appropriate microanalysis technique.     

Elevation of the glass transition temperature in flexible‐chain semicrystalline polymers

In the idealized two‐phase model of a semicrystalline polymer, the amorphous intercrystalline layers are considered to have the same properties as the fully‐amorphous polymer. In reality, these thin intercrystalline layers can be substantially influenced by the presence of the crystals, as individual polymer molecules traverse both crystalline and amorphous phases. In polymers with rigid backbone units, such as poly(etheretherketone), PEEK, previous work has shown this coupling to be particularly severe; the glass transition temperature (T_g) can be elevated by tens of degrees celsius, with the magnitude of the elevation correlating directly with the thinness of the amorphous layer. However, this connection has not been explored for flexible‐chain polymers, such as those formed from vinyl‐type monomers. Here, we examine T_g in both isotactic polystyrene (iPS) and syndiotactic polystyrene (sPS), crystallized under conditions that produce a range of amorphous layer thicknesses. T_g is indeed shown to be elevated relative to fully‐amorphous iPS and sPS, by an amount that correlates with the thinness of the amorphous layer; the magnitude of the effect is severalfold less than that in PEEK, consistent with the minimum lengths of polymer chain required to make a fold in the different cases. © 2007 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 45: 1198–1204, 2007     

Polynomial-time approximation schemes for packing and piercing fat objects

We consider two problems: given a collection of n fat objects in a fixed dimension, (1) ( packing) find the maximum subcollection of pairwise disjoint objects, and (2) ( piercing) find the minimum point set that intersects every object. Recently, Erlebach, Jansen, and Seidel gave a polynomial-time approximation scheme (PTAS) for the packing problem, based on a shifted hierarchical subdivision method. Using shifted quadtrees, we describe a similar algorithm for packing but with a smaller time bound. Erlebach et al.'s algorithm requires polynomial space. We describe a different algorithm, based on geometric separators, that requires only linear space. This algorithm can also be applied to piercing, yielding the first PTAS for that problem.     

Numerical simulation on an interaction of a vortex street with an elliptical leading edge using an unstructured grid

An algorithm for a time accurate incompressible Navier–Stokes solver on an unstructured grid is presented. The algorithm uses a second order, three‐point, backward difference formula for the physical time marching. For each time step, a divergence free flow field is obtained based on an artificial compressibility method. An implicit method with a local time step is used to accelerate the convergence for the pseudotime iteration. To validate the code, an unsteady laminar flow over a circular cylinder at a Reynolds number of 200 is calculated. The results are compared with available experimental and numerical data and good agreements are achieved. Using the developed unsteady code, an interaction of a Karman vortex street with an elliptical leading edge is simulated. The incident Karman vortex street is generated by a circular cylinder located upstream. A clustering to the path of the vortices is achieved easily due to flexibility of an unstructured grid. Details of the interaction mechanism are analysed by investigating evolutions of vortices. Characteristics of the interactions are compared for large‐ and small‐scale vortex streets. Different patterns of the interaction are observed for those two vortex streets and the observation is in agreement with experiment. Copyright © 2004 John Wiley & Sons, Ltd.     

Fully adaptive multiresolution finite volume schemes for conservation laws

The use of multiresolution decompositions in the context of finite volume schemes for conservation laws was first proposed by A. Harten for the purpose of accelerating the evaluation of numerical fluxes through an adaptive computation. In this approach the solution is still represented at each time step on the finest grid, resulting in an inherent limitation of the potential gain in memory space and computational time. The present paper is concerned with the development and the numerical analysis of fully adaptive multiresolution schemes, in which the solution is represented and computed in a dynamically evolved adaptive grid. A crucial problem is then the accurate computation of the flux without the full knowledge of fine grid cell averages. Several solutions to this problem are proposed, analyzed, and compared in terms of accuracy and complexity.