Spectral estimates for saddle point matrices arising in weak constraint four‐dimensional variational data assimilation

We consider the large sparse symmetric linear systems of equations that arise in the solution of weak constraint four‐dimensional variational data assimilation, a method of high interest for numerical weather prediction. These systems can be written as saddle point systems with a 3 × 3 block structure but block eliminations can be performed to reduce them to saddle point systems with a 2 × 2 block structure, or further to symmetric positive definite systems. In this article, we analyse how sensitive the spectra of these matrices are to the number of observations of the underlying dynamical system. We also obtain bounds on the eigenvalues of the matrices. Numerical experiments are used to confirm the theoretical analysis and bounds.     

Smart Geodesic Tracing in GRworkbench

We have developed a new tool for numerical work in General Relativity: GRworkbench. We discuss how GRworkbench's implementation of a numerically-amenable analogue to Differential Geometry facilitates the development of robust and chart-independent numerical algorithms. We consider, as an example, geodesic tracing on two charts covering the exterior Schwarzschild space-time.     

Bringing the balance function into focus

Charge balance functions provide important insight concerning hadronization and charge transport in heavy ion collisions at RHIC. Unfortunately, this observable is clouded by several effects. One of these effects is resonance production. Here, we present a simple thermal model where resonance production is included via a Monte Carlo technique based on the canonical ensemble. Resonance production is not found to have a large effect on the width.     

Ultrafast creation and annihilation of space-charge domains in a semiconductor superlattice observed by use of Terahertz fields

We report an experimental study indicating ultrafast creation and annihilation of space-charge domains in a semiconductor superlattice under the action of a THz field. Our experiment was performed for an InGaAs/InAlAs superlattice with the conduction electrons undergoing miniband transport. We applied to a superlattice a dc bias that was slightly smaller than a critical bias necessary for the formation of space-charge domains caused by a static negative differential conductivity. Additionally subjecting the superlattice to a strong THz field, resulted in a dc transport governed by the formation of domains if the frequency of the field was smaller than an upper frequency limit (~3 THz). From this frequency limit for the creation and annihilation of domains we determined the characteristic time of the domain buildup. Our analysis shows that the buildup time of domains in a wide miniband and heavily doped superlattice is limited by the relaxation time due to scattering of the miniband electrons at polar optic phonons. Our results are of importance for both an understanding of ultrafast dynamics of pattern formation in nanostructures and the development of THz electronic devices.Received: 25 March 2004, Published online: 23 July 2004PACS: 72.20.Ht High-field and nonlinear effects - 72.30. + q High-frequency effects; plasma effects - 73.21.Cd SuperlatticesK.N. Alekseev: Permanent address: Department of Physical Sciences, P.O. Box 3000, University of Oulu FIN-90014, Finland.     

The effect of nanoparticle shape on polymer‐nanocomposite rheology and tensile strength

Nanoparticles can influence the properties of polymer materials by a variety of mechanisms. With fullerene, carbon nanotube, and clay or graphene sheet nanocomposites in mind, we investigate how particle shape influences the melt shear viscosity η and the tensile strength τ, which we determine via molecular dynamics simulations. Our simulations of compact (icosahedral), tube or rod‐like, and sheet‐like model nanoparticles, all at a volume fraction ? ≈ 0.05, indicate an order of magnitude increase in the viscosity η relative to the pure melt. This finding evidently can not be explained by continuum hydrodynamics and we provide evidence that the η increase in our model nanocomposites has its origin in chain bridging between the nanoparticles. We find that this increase is the largest for the rod‐like nanoparticles and least for the sheet‐like nanoparticles. Curiously, the enhancements of η and τ exhibit opposite trends with increasing chain length N and with particle shape anisotropy. Evidently, the concept of bridging chains alone cannot account for the increase in τ and we suggest that the deformability or flexibility of the sheet nanoparticles contributes to nanocomposite strength and toughness by reducing the relative value of the Poisson ratio of the composite. The molecular dynamics simulations in the present work focus on the reference case where the modification of the melt structure associated with glass‐formation and entanglement interactions should not be an issue. Since many applications require good particle dispersion, we also focus on the case where the polymer‐particle interactions favor nanoparticle dispersion. Our simulations point to a substantial contribution of nanoparticle shape to both mechanical and processing properties of polymer nanocomposites. © 2007 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 45: 1882–1897, 2007