General-domain compressible Navier-Stokes solvers exhibiting quasi-unconditional stability and high-order accuracy in space and time

This thesis presents a new class of solvers for the subsonic compressible Navier-Stokes equations in general two- and three-dimensional spatial domains. The proposed methodology incorporates: 1) A novel linear-cost implicit solver based on use of higher-order backward differentiation formulae (BDF) and the alternating direction implicit approach (ADI); 2) A fast explicit solver; 3) Dispersionless spectral spatial discretizations; and 4) A domain decomposition strategy that negotiates the interactions between the implicit and explicit domains. In particular, the implicit methodology is quasi-unconditionally stable (it does not suffer from CFL constraints for adequately resolved flows), and it can deliver orders of time accuracy between two and six in the presence of general boundary conditions. In fact this thesis presents, for the first time in the literature, high-order time-convergence curves for Navier-Stokes solvers based on the ADI strategy---previous ADI solvers for the Navier-Stokes equations have not demonstrated orders of temporal accuracy higher than one. An extended discussion is presented in this thesis which places on a solid theoretical basis the observed quasi-unconditional stability of the methods of orders two through six. The performance of the proposed solvers is favorable. For example, a two-dimensional rough-surface configuration including boundary layer effects at Reynolds number equal to one million and Mach number 0.85 (with a well-resolved boundary layer, run up to a sufficiently long time that single vortices travel the entire spatial extent of the domain, and with spatial mesh sizes near the wall of the order of one hundred-thousandth the length of the domain) was successfully tackled in a relatively short (approximately thirty-hour) single-core run; for such discretizations an explicit solver would require truly prohibitive computing times. As demonstrated via a variety of numerical experiments in two- and three-dimensions, further, the proposed multi-domain parallel implicit-explicit implementations exhibit high-order convergence in space and time, useful stability properties, limited dispersion, and high parallel efficiency.

MEASURED AND PREDICTED LOSS GENERATION IN TRANSONIC TURBINE BLADING.

For a typical transonic turbine rotor blade, designed for use with coolant ejection, the trailing edge, or base loss is three to four times the profile boundary layer loss. The base region of such a profile is dominated by viscous effects and it seems essential to attack the problem of loss prediction by solving the compressible Navier-Stokes equations. However, such an approach is inevitably compromised by both numerical accuracy and turbulence modelling constraints. This paper describes a Navier-Stokes solver written for 2D blade-blade flows and employing a simple two-layer mixing length eddy viscosity model. Then, measured and predicted losses and base pressures are presented for two transonic rotor blades and attempts are made to assess the capabilities of the Navier-Stokes solver and to outline areas for future work.

Recent progress in high-bit-rate applications of four-wave-mixing in SOAs

Non-linearities in semiconductor optical amplifiers have been used to demonstrate a wide range of functions applicable to future optical networks such as wavelength conversion and optical switching. Four-wave-mixing effects in SOAs have been studied extensively in many laboratories with respect to the underlying physical processes and system applications. At BT Labs an optimization of SOAs for FWM has been achieved by altering the device active layer composition and by increasing the device length. We will review recent progress at BT Labs in dispersion compensation, wavelength conversion and demultiplexing at bit-rates of 40 Gbit/s using these devices.

Recent progress in high-bit-rate applications of four-wave-mixing in SOAs

Non-linearities in semiconductor optical amplifiers have been used to demonstrate a wide range of functions applicable to future optical networks such as wavelength conversion and optical switching. Four-wave-mixing effects in SOAs have been studied extensively in many laboratories with respect to the underlying physical processes and system applications. At BT Labs an optimisation of SOAs for FWM has been achieved by altering the device active layer composition and by increasing the device length. We will review recent progress at BT Labs in dispersion compensation, wavelength conversion and demultiplexing at bit-rates of 40Gbit/s using these devices.

Displacement and mixing ventilation driven by opposing wind and buoyancy

The effect of an opposing wind on the stratification and flow produced by a buoyant plume rising from a heat source on the floor of a ventilated enclosure is investigated. Ventilation openings located at high level on the windward side of the enclosure and at low level on the leeward side allow a wind-driven flow from high to low level, opposite to the buoyancy-driven flow. One of two stable steady flow regimes is established depending on a dimensionless parameter F that characterizes the relative magnitudes of the wind-driven and buoyancy-driven velocities within the enclosure, and on the time history of the flow. A third, unstable steady flow solution is identified. For small opposing winds (small F) a steady, two-layer stratification and displacement ventilation is established. Exterior fluid enters through the lower leeward openings and buoyant interior fluid leaves through the upper windward openings. As the wind speed increases, the opposing wind may cause a reversal in the flow direction. In this case, cool exterior fluid enters through the high windward openings and mixes the interior fluid, which exits through the leeward openings. There are now two possibilities. If the rate of heat input by the source exceeds the rate of heat loss through the leeward openings, the temperature of the interior increases and this flow reversal is only maintained temporarily. The buoyancy force increases with time, the flow reverts to its original direction, and steady two-layer displacement ventilation is re-established and maintained. In this regime, the increase in wind speed increases the depth and temperature of the warm upper layer, and reduces the ventilation flow rate. If, on the other hand, the heat loss exceeds the heat input, the interior cools and the buoyancy-driven flow decreases. The reversed flow is maintained, the stratification is destroyed and mixing ventilation occurs. Further increases in wind speed increase the ventilation rate and decrease the interior temperature. The transitions between the two ventilation flow patterns exhibit hysteresis. The change from displacement ventilation to mixing ventilation occurs at a higher F than the transition from mixing to displacement. Further, we find that the transition from mixing to displacement ventilation occurs at a fixed value of F, whereas the transition from displacement to mixing flow is dependent on the details of the time history of the flow and the geometry of the openings, and is not determined solely by the value of F. Theoretical models that predic t the steady stratification profiles and flow rates for the displacement and mixing ventilation, and the transitions between them, are presented and compared with measurements from laboratory experiments. The transition between these ventilation patterns completely changes the internal environment, and we discuss some of the implications for the natural ventilation of buildings. © 2004 Cambridge University Press.

Average corotation of line segments near a point and vortex identification

An easy-to-interpret kinematic quantity measuring the average corotation of material line segments near a point is introduced and applied to vortex identification. At a given point, the vector of average corotation of line segments is defined as the average of the instantaneous local rigid-body rotation over "all planar cross sections" passing through the examined point. The vortex-identification method based on average corotation is a one-parameter, region-type local method sensitive to the axial stretching rate as well as to the inner configuration of the velocity gradient tensor. The method is derived from a well-defined interpretation of the local flow kinematics to determine the "plane of swirling" and is also applicable to compressible and variable-density flows. Practical application to direct numerical simulation datasets includes a hairpin vortex of boundary-layer transition, the reconnection process of two Burgers vortices, a flow around an inclined flat plate, and a flow around a revolving insect wing. The results agree well with some popular local methods and perform better in regions of strong shearing. Copyright © 2013 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

Effect of In-mole-fraction in InGaAs overgrowth layer on self-assembled InAs/GaAs quantum dots

We have studied the optical and structural properties of InAs/GaAs QDs covered by InxGa1-xAs (0 less than or equal to x less than or equal to 0.3) layer using transmission electron microscopy, photoluminescence (PL) spectra and atomic force microscopy. We find that the strain reduces in the growth direction of InAs islands covered by InGaAs instead of GaAs layer. Significant redshift of PL peak energy and narrowing of PL linewidth are observed for the InAs QDs covered by 3 nm thick InGaAs layer. In addition, atomic force microscopy measurements indicate that the InGaAs islands will nucleate on top of InAs quantum dots, when 3 nm In0.3Ga0.7As overgrowth layer is deposited. This result can well explain the PL intensify degradation and linewidth increment of quantum dots with a higher In-mole-fraction InGaAs layer. The energy gap change of InAs QDs covered by InGaAs may be explained in terms of reducing strain, suppressing compositional mixing and increasing island height. (C) 2000 Elsevier Science B.V. All rights reserved.

A Planar InGaAs/InP Geiger Mode Avalanche Photodiode with Cascade Edge Breakdown Suppression

A Geiger mode planar InGaAs/InP avalanche photodiode （APD） with a cascade peripheral junction structure to suppress edge breakdowns is designed by finite-element analysis. The photodiode breakdown voltage is reduced to 54.3V by controlling the central junction depth, while the electric field distribution along the device central axis is controlled by adjusting doping level and thickness of the lnP field control layer. Using a cascade junction structure at the periphery of the active area, premature edge breakdowns are effectively suppressed. The simulations show that the quadra-cascade structure is a good trade-off between suppression performance and fabrication complexity, with a reduced peak electric field of 5.2 × 10~5 kV/cm and a maximum hole ionization integral of 1. 201. Work presented in this paper provides an effective way to design high performance photon counting InGaAs/InP avalanche photodiodes.

Weak epitaxy growth of metal-free phthalocyanine on p-sexiphenyl monolayer and double-layer films

Weak epitaxy growth (WEG) can afford high-mobility thin films of disk-like organic semiconductor of which mobility is up to the level of the corresponding single crystals. We investigated the WEG behavior and mechanism of planar phthalocyanine in the model system of metal-free phthalocyanine (H2Pc) grown on p-sexiphenyl (p-6P) ultrathin films (monolayers and double layers). Highly oriented H2Pc films with molecules standing up exhibited two kinds of different in-plane orientations, i.e., three sets of in-plane orientations and only one set of in-plane orientation, on p-6P monolayer and double-layer films, respectively.

Simulation of the ocean surface mixed layer under the wave breaking

A one-dimensional mixed-layer model, including a Mellor-Yamada level 2.5 turbulence closure scheme, was implemented to investigate the dynamical and thermal structures of the ocean surface mixed layer in the northern South China Sea. The turbulent kinetic energy released through wave breaking was incorporated into the model as a source of energy at the ocean surface, and the influence of the breaking waves on the mixed layer was studied. The numerical simulations show that the simulated SST is overestimated in summer without the breaking waves. However, the cooler SST is simulated when the effect of the breaking waves is considered, the corresponding discrepancy with the observed data decreases up to 20% and the MLD calculated averagely deepens 3.8 m. Owing to the wave-enhanced turbulence mixing in the summertime, the stratification at the bottom of the mixed layer was modified and the temperature gradient spread throughout the whole thermocline compared with the concentrated distribution without wave breaking.

Effect of Stokes drift on upper ocean mixing

Stokes drift is the main source of vertical vorticity in the ocean mixed layer. In the ways of Coriolis - Stokes forcing and Langmuir circulations, Stokes drift can substantially affect the whole mixed layer. A modified Mellor-Yamada 2.5 level turbulence closure model is used to parameterize its effect on upper ocean mixing conventionally. Results show that comparing surface heating with wave breaking, Stokes drift plays the most important role in the entire ocean mixed layer, especially in the subsurface layer. As expected, Stokes drift elevates both the dissipation rate and the turbulence energy in the upper ocean mixing. Also, influence of the surface heating, wave breaking and wind speed on Stokes drift is investigated respectively. Research shows that it is significant and important to assessing the Stokes drift into ocean mixed layer studying. The laboratory observations are supporting numerical experiments quantitatively.

Wave breaking on turbulent energy budget in the ocean surface mixed layer

As an important physical process at the air-sea interface, wave movement and breaking have a significant effect on the ocean surface mixed layer (OSML). When breaking waves occur at the ocean surface, turbulent kinetic energy (TKE) is input downwards, and a sublayer is formed near the surface and turbulence vertical mixing is intensively enhanced. A one-dimensional ocean model including the Mellor-Yamada level 2.5 turbulence closure equations was employed in our research on variations in turbulent energy budget within OSML. The influence of wave breaking could be introduced into the model by modifying an existing surface boundary condition of the TKE equation and specifying its input. The vertical diffusion and dissipation of TKE were effectively enhanced in the sublayer when wave breaking was considered. Turbulent energy dissipated in the sublayer was about 92.0% of the total depth-integrated dissipated TKE, which is twice higher than that of non-wave breaking. The shear production of TKE decreased by 3.5% because the mean flow fields tended to be uniform due to wave-enhanced turbulent mixing. As a result, a new local equilibrium between diffusion and dissipation of TKE was reached in the wave-enhanced layer. Below the sublayer, the local equilibrium between shear production and dissipation of TKE agreed with the conclusion drawn from the classical law-of-the-wall (Craig and Banner, 1994).