SBLOCK: A framework for efficient stencil-based PDE solvers on multi-core platforms

We present a new software framework for the implementation of applications that use stencil computations on block-structured grids to solve partial differential equations. A key feature of the framework is the extensive use of automatic source code generation which is used to achieve high performance on a range of leading multi-core processors. Results are presented for a simple model stencil running on Intel and AMD CPUs as well as the NVIDIA GT200 GPU. The generality of the framework is demonstrated through the implementation of a complete application consisting of many different stencil computations, taken from the field of computational fluid dynamics. © 2010 IEEE.

Global linear and nonlinear stability of viscous confined plane wakes with co-flow

The global stability of confined uniform density wakes is studied numerically, using two-dimensional linear global modes and nonlinear direct numerical simulations. The wake inflow velocity is varied between different amounts of co-flow (base bleed). In accordance with previous studies, we find that the frequencies of both the most unstable linear and the saturated nonlinear global mode increase with confinement. For wake Reynolds number Re = 100 we find the confinement to be stabilising, decreasing the growth rate of the linear and the saturation amplitude of the nonlinear modes. The dampening effect is connected to the streamwise development of the base flow, and decreases for more parallel flows at higher Re. The linear analysis reveals that the critical wake velocities are almost identical for unconfined and confined wakes at Re ≈ 400. Further, the results are compared with literature data for an inviscid parallel wake. The confined wake is found to be more stable than its inviscid counterpart, whereas the unconfined wake is more unstable than the inviscid wake. The main reason for both is the base flow development. A detailed comparison of the linear and nonlinear results reveals that the most unstable linear global mode gives in all cases an excellent prediction of the initial nonlinear behaviour and therefore the stability boundary. However, the nonlinear saturated state is different, mainly for higher Re. For Re = 100, the saturated frequency differs less than 5% from the linear frequency, and trends regarding confinement observed in the linear analysis are confirmed.

Characterization of dynamic nonlinear absorption of carbon nanotube saturable absorber

Dynamic nonlinear absorption of composite-type single-wall carbon nanotube saturable absorbers is characterized using both femtosecond and picosecond pump pulses. Results are compared with numerical simulations based on two commonly used saturable absorber models. © 2010 Optical Society of America.

The nonlinear heat release response of stratified lean-premixed flames to acoustic velocity oscillations

This paper analyzes the forced response of swirl-stabilized lean-premixed flames to high-amplitude acoustic forcing in a laboratory-scale stratified burner operated with CH4 and air at atmospheric pressure. The double-swirler, double-channel annular burner was specially designed to generate high-amplitude acoustic velocity oscillations and a radial equivalence ratio gradient at the inlet of the combustion chamber. Temporal oscillations of equivalence ratio along the axial direction are dissipated over a long distance, and therefore the effects of time-varying fuel/air ratio on the response are not considered in the present investigation. Simultaneous measurements of inlet velocity and heat release rate oscillations were made using a constant temperature anemometer and photomultiplier tubes with narrow-band OH*/CH* interference filters. Time-averaged and phase-synchronized CH* chemiluminescence intensities were measured using an intensified CCD camera. The measurements show that flame stabilization mechanisms vary depending on equivalence ratio gradients for a constant global equivalence ratio (φg=0.60). Under uniformly premixed conditions, an enveloped M-shaped flame is observed. In contrast, under stratified conditions, a dihedral V-flame and a toroidal detached flame develop in the outer stream and inner stream fuel enrichment cases, respectively. The modification of the stabilization mechanism has a significant impact on the nonlinear response of stratified flames to high-amplitude acoustic forcing (u'/U∼0.45 and f=60, 160Hz). Outer stream enrichment tends to improve the flame's stiffness with respect to incident acoustic/vortical disturbances, whereas inner stream stratification tends to enhance the nonlinear flame dynamics, as manifested by the complex interaction between the swirl flame and large-scale coherent vortices with different length scales and shedding points. It was found that the behavior of the measured flame describing functions (FDF), which depend on radial fuel stratification, are well correlated with previous measurements of the intensity of self-excited combustion instabilities in the stratified swirl burner. The results presented in this paper provide insight into the impact of nonuniform reactant stoichiometry on combustion instabilities, its effect on flame location and the interaction with unsteady flow structures. © 2011 The Combustion Institute.

Optimal control with stochastic PDE constraints and uncertain controls

The optimal control of problems that are constrained by partial differential equations with uncertainties and with uncertain controls is addressed. The Lagrangian that defines the problem is postulated in terms of stochastic functions, with the control function possibly decomposed into an unknown deterministic component and a known zero-mean stochastic component. The extra freedom provided by the stochastic dimension in defining cost functionals is explored, demonstrating the scope for controlling statistical aspects of the system response. One-shot stochastic finite element methods are used to find approximate solutions to control problems. It is shown that applying the stochastic collocation finite element method to the formulated problem leads to a coupling between stochastic collocation points when a deterministic optimal control is considered or when moments are included in the cost functional, thereby forgoing the primary advantage of the collocation method over the stochastic Galerkin method for the considered problem. The application of the presented methods is demonstrated through a number of numerical examples. The presented framework is sufficiently general to also consider a class of inverse problems, and numerical examples of this type are also presented. © 2011 Elsevier B.V.