Wideband finite difference time domain implementation of surface impedance boundary conditions for good conductors

Surface impedance boundary conditions are used to reduce the solution volume during the analysis of scattering from lossy dielectric objects. In a finite difference solution, they also can be used to avoid using small cells, made necessary by shorter wavelengths in conducting media, throughout the solution volume. A one dimensional implementation is presented for a surface impedance boundary condition for good conductors in the Finite Difference Time Domain (FDTD) technique. In order to illustrate the FDTD surface impedance boundary condition, a planar air-lossy dielectric interface is considered

Finite difference time domain implementation of surface impedance boundary conditions

Surface impedance boundary conditions are employed to reduce the solution volume during the analysis of scattering from lossy dielectric objects. In a finite difference solution, they also can be utilized to avoid using small cells, made necessary by shorter wavelengths in conducting media throughout the solution volume. The standard approach is to approximate the surface impedance over a very small bandwidth by its value at the center frequency, and then use that result in the boundary condition. Two implementations of the surface impedance boundary condition are presented. One implementation is a constant surface impedance boundary condition and the other is a dispersive surface impedance boundary condition that is applicable over a very large frequency bandwidth and over a large range of conductivities. Frequency domain results are presented in one dimension for two conductivity values and are compared with exact results. Scattering width results from an infinite square cylinder are presented as a 2-D demonstration. Extensions to 3-D should be straightforward

Multiwavelength Mass Comparisons of the z~0.3 CNOC Cluster Sample

Results are presented from a detailed analysis of optical and X-ray observations of moderate-redshift galaxy clusters from the Canadian Network for Observational Cosmology (CNOC) subsample of the EMSS. The combination of extensive optical and deep X-ray observations of these clusters make them ideal candidates for multiwavelength mass comparison studies. X-ray surface brightness profiles of 14 clusters with 0.17<z<0.55 are constructed from Chandra observations and fit to single and double beta-models. Spatially resolved temperature analysis is performed, indicating that five of the clusters in this sample exhibit temperature gradients within their inner 60-200 kpc. Integrated spectra extracted within R_2500 provide temperature, abundance, and luminosity information. Under assumptions of hydrostatic equilibrium and spherical symmetry, we derive gas and total masses within R_2500 and R_200. We find an average gas mass fraction within R_200 of 0.136 +/- 0.004, resulting in Omega_m=0.28 +/- 0.01 (formal error). We also derive dynamical masses for these clusters to R_200. We find no systematic bias between X-ray and dynamical methods across the sample, with an average M(dyn)/M(X-ray) = 0.97 +/- 0.05. We also compare X-ray masses to weak lensing mass estimates of a subset of our sample, resulting in a weighted average of M(lens)/M(X-ray) of 0.99 +/- 0.07. We investigate X-ray scaling relationships and find powerlaw slopes which are slightly steeper than the predictions of self-similar models, with an E(z)^(-1) Lx-Tx slope of 2.4 +/- 0.2 and an E(z) M_2500-Tx slope of 1.7 +/- 0.1. Relationships between red-sequence optical richness (B_gc,red) and global cluster X-ray properties (Tx, Lx and M_2500) are also examined and fitted.Comment: Astrophysical Journal, 48 pages, 11 figures, LaTeX. Added correction to surface brightness normalization of MS1512.4+3647, corrections to sample gas mass fractions and calculated value of Omega_m. Figure resolution has been reduced to comply with astro-ph upload requirement

A high order compact scheme for hypersonic aerothermodynamics

A novel high order compact scheme for solving the compressible Navier-Stokes equations has been developed. The scheme is an extension of a method originally proposed for solving the Euler equations, and combines several techniques for the solution of compressible flowfields, such as upwinding, limiting and flux vector splitting, with the excellent properties of high order compact schemes. Extending the method to the Navier-Stokes equations is achieved via a Kinetic Flux Vector Splitting technique, which represents an unusual and attractive way to include viscous effects. This approach offers a more accurate and less computationally expensive technique than discretizations based on more conventional operator splitting. The Euler solver has been validated against several inviscid test cases, and results for several viscous test cases are also presented. The results confirm that the method is stable, accurate and has excellent shock-capturing capabilities for both viscous and inviscid flows

Chemical aging of m-xylene secondary organic aerosol: laboratory chamber study

Secondary organic aerosol (SOA) can reside in the atmosphere for a week or more. While its initial formation from the gas-phase oxidation of volatile organic compounds tends to take place in the first few hours after emission, SOA can continue to evolve chemically over its atmospheric lifetime. Simulating this chemical aging over an extended time in the laboratory has proven to be challenging. We present here a procedure for studying SOA aging in laboratory chambers that is applied to achieve 36 h of oxidation. The formation and evolution of SOA from the photooxidation of m-xylene under low-NO_x conditions and in the presence of either neutral or acidic seed particles is studied. In SOA aging, increasing molecular functionalization leads to less volatile products and an increase in SOA mass, whereas gas- or particle-phase fragmentation chemistry results in more volatile products and a loss of SOA. The challenge is to discern from measured chamber variables the extent to which these processes are important for a given SOA system. In the experiments conducted, m-xylene SOA mass, calculated under the assumption of size-invariant particle composition, increased over the initial 12–13 h of photooxidation and decreased beyond that time, suggesting the existence of fragmentation chemistry. The oxidation of the SOA, as manifested in the O:C elemental ratio and fraction of organic ion detected at m/z 44 measured by the Aerodyne aerosol mass spectrometer, increased continuously starting after 5 h of irradiation until the 36 h termination. This behavior is consistent with an initial period in which, as the mass of SOA increases, products of higher volatility partition to the aerosol phase, followed by an aging period in which gas- and particle-phase reaction products become increasingly more oxidized. When irradiation is stopped 12.4 h into one experiment, and OH generation ceases, minimal loss of SOA is observed, indicating that the loss of SOA is either light- or OH-induced. Chemical ionization mass spectrometry measurements of low-volatility m-xylene oxidation products exhibit behavior indicative of continuous photooxidation chemistry. A condensed chemical mechanism of m-xylene oxidation under low-NO_x conditions is capable of reproducing the general behavior of gas-phase evolution observed here. Moreover, order of magnitude analysis of the mechanism suggests that gas-phase OH reaction of low volatility SOA precursors is the dominant pathway of aging in the m-xylene system although OH reaction with particle surfaces cannot be ruled out. Finally, the effect of size-dependent particle composition and size-dependent particle wall loss rates on different particle wall loss correction methods is discussed