Capillary Waves at Liquid/Vapor Interfaces: A Molecular Dynamics Simulation

Evidence for capillary waves at a liquid/vapor interface are presented from extensive molecular dynamics simulations of a system containing up to 1.24 million Lennard-Jones particles. Careful measurements show that the total interfacial width depends logarithmically on $L_\parallel$, the length of the simulation cell parallel to the interface, as predicted theoretically. The strength of the divergence of the interfacial width on $L_\parallel$ depends inversely on the surface tension $\gamma$. This allows us to measure $\gamma$ two ways since $\gamma$ can also be obtained from the difference in the pressure parallel and perpendicular to the interface. These two independent measures of $\gamma$ agree provided that the interfacial order parameter profile is fit to an error function and not a hyperbolic tangent, as often assumed. We explore why these two common fitting functions give different results for $\gamma$

Shuffle operations on discrete paths

AbstractWe consider the shuffle operation on paths and study some parameters. In the case of square lattices, shuffling with a particular periodic word (of period 2) corresponding to paperfoldings reveals some characteristic properties: closed paths remain closed; the area and perimeter double; the center of gravity moves under a 45∘ rotation and a 2 zoom factor. We also observe invariance properties for the associated Dragon curves. Moreover, replacing square lattice paths by paths involving 2kπ/N-turns, we find analogous results using more general shuffles

Use of High-resolution WRF Simulations to Forecast Lightning Threat

Recent observational studies have confirmed the existence of a robust statistical relationship between lightning flash rates and the amount of large precipitating ice hydrometeors in storms. This relationship is exploited, in conjunction with the capabilities of recent forecast models such as WRF, to forecast the threat of lightning from convective storms using the output fields from the model forecasts. The simulated vertical flux of graupel at -15C is used in this study as a proxy for charge separation processes and their associated lightning risk. Six-h simulations are conducted for a number of case studies for which three-dimensional lightning validation data from the North Alabama Lightning Mapping Array are available. Experiments indicate that initialization of the WRF model on a 2 km grid using Eta boundary conditions, Doppler radar radial velocity and reflectivity fields, and METAR and ACARS data yield the most realistic simulations. An array of subjective and objective statistical metrics are employed to document the utility of the WRF forecasts. The simulation results are also compared to other more traditional means of forecasting convective storms, such as those based on inspection of the convective available potential energy field

Infrared Imaging of Capella with the IOTA Closure Phase Interferometer

We present infrared aperture synthesis maps produced with the upgraded IOTA interferometer. Michelson interferograms on the close binary system Capella (Alpha Aur) were obtained in the H-band between 2002 November 12 and 16 using the IONIC3 beam combiner. With baselines of 15m < B < 38m, we were able to determine the relative position of the binary components with milliarcsecond (mas) precision and to track their movement along the approx. 14 degree arc covered by our observation run. We briefly describe the algorithms used for visibility and closure phase estimation. Three different Hybrid Mapping and Bispectrum Fitting techniques were implemented within one software framework and used to reconstruct the source brightness distribution. By dividing our data into subsets, the system could be mapped at three epochs, revealing the motion of the stars. The precise position of the binary components was also determined with model fits, which in addition revealed I_Aa/I_Ab=1.49 +/- 0.10 and apparent stellar uniform-disk (UD) diameters of Theta_Aa=8.9 +/- 0.6 mas and Theta_Ab=5.8 +/- 0.8 mas. To improve the u, v-plane coverage, we compensated this orbital motion by applying a rotation-compensating coordinate transformation. The resulting model-independent map with a beam size of 5.4 x 2.6 mas allows the resolution of the stellar surfaces of the Capella giants themselves.Comment: Accepted by the Astronomical Journal (2005-03-21

A Model for the Elasticity of Compressed Emulsions

We present a new model to describe the unusual elastic properties of compressed emulsions. The response of a single droplet under compression is investigated numerically for different Wigner-Seitz cells. The response is softer than harmonic, and depends on the coordination number of the droplet. Using these results, we propose a new effective inter-droplet potential which is used to determine the elastic response of a monodisperse collection of disordered droplets as a function of volume fraction. Our results are in excellent agreement with recent experiments. This suggests that anharmonicity, together with disorder, are responsible for the quasi-linear increase of $G$ and $\Pi$ observed at $\varphi_c$.Comment: RevTeX with psfig-included figures and a galley macr