Analysis of a parallelized nonlinear elliptic boundary value problem solver with application to reacting flows

A parallelized finite difference code based on the Newton method for systems of nonlinear elliptic boundary value problems in two dimensions is analyzed in terms of computational complexity and parallel efficiency. An approximate cost function depending on 15 dimensionless parameters is derived for algorithms based on stripwise and boxwise decompositions of the domain and a one-to-one assignment of the strip or box subdomains to processors. The sensitivity of the cost functions to the parameters is explored in regions of parameter space corresponding to model small-order systems with inexpensive function evaluations and also a coupled system of nineteen equations with very expensive function evaluations. The algorithm was implemented on the Intel Hypercube, and some experimental results for the model problems with stripwise decompositions are presented and compared with the theory. In the context of computational combustion problems, multiprocessors of either message-passing or shared-memory type may be employed with stripwise decompositions to realize speedup of O(n), where n is mesh resolution in one direction, for reasonable n

Analysis of the symbiotic star AG Pegasi

High and low dispersion IUE data are analyzed in conjunction with coincident ground based spectrophotometric scans and supplementary infrared photometry of the symbiotic object AG Pegasi. The IUE observations yield an improved value of E(B-V) = 0.12. The two stellar components are easily recognized in the spectra. The cool component may be an M1.7 III star and the hot component appears to have T (sub eff) of approximately 30000 K. The emission lines observed in the ultraviolet indicate two or three distince emitting regions. Nebular component ultraviolet intercombination lines suggest an electron density of several times 10 billion/cu cm

Domain-decomposed preconditionings for transport operators

The performance was tested of five different interface preconditionings for domain decomposed convection diffusion problems, including a novel one known as the spectral probe, while varying mesh parameters, Reynolds number, ratio of subdomain diffusion coefficients, and domain aspect ratio. The preconditioners are representative of the range of practically computable possibilities that have appeared in the domain decomposition literature for the treatment of nonoverlapping subdomains. It is shown that through a large number of numerical examples that no single preconditioner can be considered uniformly superior or uniformly inferior to the rest, but that knowledge of particulars, including the shape and strength of the convection, is important in selecting among them in a given problem

Challenges in imaging and predictive modeling of rhizosphere processes

Background Plant-soil interaction is central to human food production and ecosystem function. Thus, it is essential to not only understand, but also to develop predictive mathematical models which can be used to assess how climate and soil management practices will affect these interactions. Scope In this paper we review the current developments in structural and chemical imaging of rhizosphere processes within the context of multiscale mathematical image based modeling. We outline areas that need more research and areas which would benefit from more detailed understanding. Conclusions We conclude that the combination of structural and chemical imaging with modeling is an incredibly powerful tool which is fundamental for understanding how plant roots interact with soil. We emphasize the need for more researchers to be attracted to this area that is so fertile for future discoveries. Finally, model building must go hand in hand with experiments. In particular, there is a real need to integrate rhizosphere structural and chemical imaging with modeling for better understanding of the rhizosphere processes leading to models which explicitly account for pore scale processes

Analytical and experimental studies of shock interference heating in hypersonic flows

An analytical and experimental study is presented of the aerodynamic heating resulting from six types of shock interference patterns encountered in high speed flow. Centerline measurements of pressure and heat transfer distributions on basic bodies were obtained in four wind tunnels for Mach numbers from 6 to 20, specific heat ratios from 1.27 to 1.67, and free stream Reynolds numbers from 3 million to 25.6 million per meter. Peak heating and peak pressures up to 17 and 7.5 times stagnation values, respectively, were measured. In general, results obtained from semiempirical methods developed for each of the six types of interference agreed with the experimental peaks

Analysis of high excitation planetary nebulae

Combination of extensive ground-based spectroscopic observation of high excitation planetary with IUE data permit determination not only of improved diagnostics but also better abundances for elements such as C and N that are well represented in the ultraviolet spectra and also C, Ar and metals Na, Ca and K whose lines appear in the wavelength 3200-8100 A region

Electron-Acoustic Phonon Energy Loss Rate in Multi-Component Electron Systems with Symmetric and Asymmetric Coupling Constants

We consider electron-phonon (\textit{e-ph}) energy loss rate in 3D and 2D multi-component electron systems in semiconductors. We allow general asymmetry in the \textit{e-ph} coupling constants (matrix elements), i.e., we allow that the coupling depends on the electron sub-system index. We derive a multi-component \textit{e-ph}power loss formula, which takes into account the asymmetric coupling and links the total \textit{e-ph} energy loss rate to the density response matrix of the total electron system. We write the density response matrix within mean field approximation, which leads to coexistence of\ symmetric energy loss rate $F_{S}(T)$ and asymmetric energy loss rate $F_{A}(T)$ with total energy loss rate $F(T)=F_{S}(T)+F_{A}(T)$ at temperature $T$. The symmetric component F_{S}(T) $is equivalent to the conventional single-sub-system energy loss rate in the literature, and in the Bloch-Gr\"{u}neisen limit we reproduce a set of well-known power laws$F_{S}(T)\propto T^{n_{S}}$, where the prefactor and power$n_{S}$depend on electron system dimensionality and electron mean free path. For$F_{A}(T)$we produce a new set of power laws F_{A}(T)\propto T^{n_{A}}$. Screening strongly reduces the symmetric coupling, but the asymmetric coupling is unscreened, provided that the inter-sub-system Coulomb interactions are strong. The lack of screening enhances $F_{A}(T)$ and the total energy loss rate $F(T)$. Especially, in the strong screening limit we find $F_{A}(T)\gg F_{S}(T)$. A canonical example of strongly asymmetric \textit{e-ph} matrix elements is the deformation potential coupling in many-valley semiconductors.Comment: v2: Typos corrected. Some notations changed. Section III.C is embedded in Section III.B. Paper accepted to PR

Instantaneous Normal Mode analysis of liquid HF

We present an Instantaneous Normal Modes analysis of liquid HF aimed to clarify the origin of peculiar dynamical properties which are supposed to stem from the arrangement of molecules in linear hydrogen-bonded network. The present study shows that this approach is an unique tool for the understanding of the spectral features revealed in the analysis of both single molecule and collective quantities. For the system under investigation we demonstrate the relevance of hydrogen-bonding ``stretching'' and fast librational motion in the interpretation of these features.Comment: REVTeX, 7 pages, 5 eps figures included. Minor changes in the text and in a figure. Accepted for publication in Phys. Rev. Let

Configurational entropy of hard spheres

We numerically calculate the configurational entropy S_conf of a binary mixture of hard spheres, by using a perturbed Hamiltonian method trapping the system inside a given state, which requires less assumptions than the previous methods [R.J. Speedy, Mol. Phys. 95, 169 (1998)]. We find that S_conf is a decreasing function of packing fraction f and extrapolates to zero at the Kauzmann packing fraction f_K = 0.62, suggesting the possibility of an ideal glass-transition for hard spheres system. Finally, the Adam-Gibbs relation is found to hold.Comment: 10 pages, 6 figure