Numerical studies of identification in nonlinear distributed parameter systems

An abstract approximation framework and convergence theory for the identification of first and second order nonlinear distributed parameter systems developed previously by the authors and reported on in detail elsewhere are summarized and discussed. The theory is based upon results for systems whose dynamics can be described by monotone operators in Hilbert space and an abstract approximation theorem for the resulting nonlinear evolution system. The application of the theory together with numerical evidence demonstrating the feasibility of the general approach are discussed in the context of the identification of a first order quasi-linear parabolic model for one dimensional heat conduction/mass transport and the identification of a nonlinear dissipation mechanism (i.e., damping) in a second order one dimensional wave equation. Computational and implementational considerations, in particular, with regard to supercomputing, are addressed

Energy and momentum of cylindrical gravitational waves. II

Recently Nathan Rosen and the present author obtained the energy and momentum densities of cylindrical gravitational waves in Einstein's prescription and found them to be finite and reasonable. In the present paper we calculate the same in prescriptions of Tolman as well as Landau and Lifshitz and discuss the results.Comment: 8 pages, LaTex, To appear in Pramana- J. Physic

ORFEUS II and IUE Spectroscopy of EX Hydrae

Using ORFEUS-SPAS II FUV spectra, IUE UV spectra, and archival EUVE deep survey photometry, we present a detailed picture of the behavior of the magnetic cataclysmic variable EX Hydrae. Like HUT spectra of this source, the FUV and UV spectra reveal broad emission lines of He II, C II-IV, N III and V, O VI, Si III-IV, and Al III superposed on a continuum which is blue in the UV and nearly flat in the FUV. Like ORFEUS spectra of AM Her, the O VI doublet is resolved into broad and narrow emission components. Consistent with its behavior in the optical, the FUV and UV continuum flux densities, the FUV and UV broad emission line fluxes, and the radial velocity of the O VI broad emission component all vary on the spin phase of the white dwarf, with the maximum of the FUV and UV continuum and broad emission line flux light curves coincident with maximum blueshift of the broad O VI emission component. On the binary phase, the broad dip in the EUV light curve is accompanied by strong eclipses of the UV emission lines and by variations in both the flux and radial velocity of the O VI narrow emission component. The available data are consistent with the accretion funnel being the source of the FUV and UV continuum and the O VI broad emission component, and the white dwarf being the source of the O VI narrow emission component.Comment: 21 pages, 10 Postscript figures; LaTeX format, uses aaspp4.sty; table2.tex included separately because it must be printed sideways - see instructions in the file; accepted on 1999 Feb 20 for publication in The Astrophysical Journa

Dynamics and Structure of Three-Dimensional Trans-Alfvenic Jets. II. The Effect of Density and Winds

Two three-dimensional magnetohydrodynamical simulations of strongly magnetized conical jets, one with a poloidal and one with a helical magnetic field, have been performed. In the poloidal simulation a significant sheath (wind) of magnetized moving material developed and partially stabilized the jet to helical twisting. The fundamental pinch mode was not similarly affected and emission knots developed in the poloidal simulation. Thus, astrophysical jets surrounded by outflowing winds could develop knotty structures along a straight jet triggered by pinching. Where helical twisting dominated the dynamics, magnetic field orientation along the line-of-sight could be organized by the toroidal flow field accompanying helical twisting. On astrophysical jets such structure could lead to a reversal of the direction of Faraday rotation in adjacent zones along a jet. Theoretical analysis showed that the different dynamical behavior of the two simulations could be entirely understood as a result of dependence on the velocity shear between jet and wind which must exceed a surface Alfven speed before the jet becomes unstable to helical and higher order modes of jet distortion.Comment: 25 pages, 15 figures, in press Astrophysical Journal (September

Geodetic, teleseismic, and strong motion constraints on slip from recent southern Peru subduction zone earthquakes

We use seismic and geodetic data both jointly and separately to constrain coseismic slip from the 12 November 1996 M_w 7.7 and 23 June 2001 M_w 8.5 southern Peru subduction zone earthquakes, as well as two large aftershocks following the 2001 earthquake on 26 June and 7 July 2001. We use all available data in our inversions: GPS, interferometric synthetic aperture radar (InSAR) from the ERS-1, ERS-2, JERS, and RADARSAT-1 satellites, and seismic data from teleseismic and strong motion stations. Our two-dimensional slip models derived from only teleseismic body waves from South American subduction zone earthquakes with M_w > 7.5 do not reliably predict available geodetic data. In particular, we find significant differences in the distribution of slip for the 2001 earthquake from models that use only seismic (teleseismic and two strong motion stations) or geodetic (InSAR and GPS) data. The differences might be related to postseismic deformation or, more likely, the different sensitivities of the teleseismic and geodetic data to coseismic rupture properties. The earthquakes studied here follow the pattern of earthquake directivity along the coast of western South America, north of 5°S, earthquakes rupture to the north; south of about 12°S, directivity is southerly; and in between, earthquakes are bilateral. The predicted deformation at the Arequipa GPS station from the seismic-only slip model for the 7 July 2001 aftershock is not consistent with significant preseismic motion

Phenomenology of loop quantum cosmology

After introducing the basic ingredients of Loop Quantum Cosmology, I will briefly discuss some of its phenomenological aspects. Those can give some useful insight about the full Loop Quantum Gravity theory and provide an answer to some long-standing questions in early universe cosmology.Comment: 16 pages, 3 figures; Invited talk in the First Mediterranean Conference on Classical and Quantum Gravity (Crete, Greece

A Quantitative Non-radial Oscillation Model for the Subpulses in PSR B0943+10

In this paper, we analyze time series measurements of PSR B0943+10 and fit them with a non-radial oscillation model. The model we apply was first developed for total intensity measurements in an earlier paper, and expanded to encompass linear polarization in a companion paper to this one. We use PSR B0943+10 for the initial tests of our model because it has a simple geometry, it has been exhaustively studied in the literature, and its behavior is well-documented. As prelude to quantitative fitting, we have reanalyzed previously published archival data of PSR B0943+10 and uncovered subtle but significant behavior that is difficult to explain in the framework of the drifting spark model. Our fits of a non-radial oscillation model are able to successfully reproduce the observed behavior in this pulsar.Comment: 45 pages, 16 figures, accepted Ap

FUSE Observations of the Dwarf Nova SW UMa During Quiescence

We present spectroscopic observations of the short-period cataclysmic variable SW Ursa Majoris, obtained by the Far Ultraviolet Spectroscopic Explorer (FUSE) satellite while the system was in quiescence. The data include the resonance lines of O VI at 1031.91 and 1037.61 A. These lines are present in emission, and they exhibit both narrow (~ 150 km/s) and broad (~ 2000 km/s) components. The narrow O VI emission lines exhibit unusual double-peaked and redshifted profiles. We attribute the source of this emission to a cooling flow onto the surface of the white dwarf primary. The broad O VI emission most likely originates in a thin, photoionized surface layer on the accretion disk. We searched for emission from H_2 at 1050 and 1100 A, motivated by the expectation that the bulk of the quiescent accretion disk is in the form of cool, molecular gas. If H_2 is present, then our limits on the fluxes of the H_2 lines are consistent with the presence of a surface layer of atomic H that shields the interior of the disk. These results may indicate that accretion operates primarily in the surface layers of the disk in SW UMa. We also investigate the far-UV continuum of SW UMa and place an upper limit of 15,000 K on the effective temperature of the white dwarf.Comment: 21 Pages, 3 figures, to be published in Ap

Faster Family-wise Error Control for Neuroimaging with a Parametric Bootstrap

In neuroimaging, hundreds to hundreds of thousands of tests are performed across a set of brain regions or all locations in an image. Recent studies have shown that the most common family-wise error (FWE) controlling procedures in imaging, which rely on classical mathematical inequalities or Gaussian random field theory, yield FWE rates that are far from the nominal level. Depending on the approach used, the FWER can be exceedingly small or grossly inflated. Given the widespread use of neuroimaging as a tool for understanding neurological and psychiatric disorders, it is imperative that reliable multiple testing procedures are available. To our knowledge, only permutation joint testing procedures have been shown to reliably control the FWER at the nominal level. However, these procedures are computationally intensive due to the increasingly available large sample sizes and dimensionality of the images, and analyses can take days to complete. Here, we develop a parametric bootstrap joint testing procedure. The parametric bootstrap procedure works directly with the test statistics, which leads to much faster estimation of adjusted \emph{p}-values than resampling-based procedures while reliably controlling the FWER in sample sizes available in many neuroimaging studies. We demonstrate that the procedure controls the FWER in finite samples using simulations, and present region- and voxel-wise analyses to test for sex differences in developmental trajectories of cerebral blood flow

Quantification of the Thermal Hazard from Metallic and Organic Dust Flash Fires

PresentationAn in-house constructed nominally 20-L flash fire apparatus was used to evaluate and compare the flash fires fueled by an organic dust (non-dairy coffee creamer) a metal dust (aluminum) and a flammable gas (methane). Dispersion was achieved using an injection system similar to the injection systems found in standard Siwek 20-L combustion chambers and a 10-J spark igniter was used to ignite the fuels. A heat flux gauge, thermocouples, an Infrared video camera, and an HD video camera were all used to evaluate the severity of the flash fires. Multiple concentrations of dusts and a stoichiometric methane mixture were tested. All measurement methods showed reasonable agreement when ranking the severity of different deflagrations, but thermocouple and heat flux gauge measurements were sensitive to the position of the flame, leading to some inconsistency. IR video measurements provided fireball dimensions and growth rates, and relative temperatures for dust-fueled deflagration, but were unable to accurately assess the high-turbulence premixed methane flash fires due to the high burning velocities and 30 frames per second limitation. The IR camera was also limited to 1200o C, which is inadequate for the temperature of some metal dust deflagrations, including aluminum and requires adjustment of the material emissivity, which would require additional analysis and testing. Measurement strategies for a next-generation flash fire testing apparatus are proposed based on the results of this study