Efficient calculation of the green function for acoustic propagation above a homogeneous impedance plane

This paper is concerned with the problem of propagation from a monofrequency coherent line source above a plane of homogeneous surface impedance. The solution of this problem occurs in the kernel of certain boundary integral equation formulations of acoustic propagation above an impedance boundary, and the discussion of the paper is motivated by this application. The paper starts by deriving representations, as Laplace-type integrals, of the solution and its first partial derivatives. The evaluation of these integral representations by Gauss-Laguerre quadrature is discussed, and theoretical bounds on the truncation error are obtained. Specific approximations are proposed which are shown to be accurate except in the very near field, for all angles of incidence and a wide range of values of surface impedance. The paper finishes with derivations of partial results and analogous Laplace-type integral representations for the case of a point source.

Ice internal friction: Standard theoretical perspectives on friction codified, adapted for the unusual rheology of ice, and unified

Sea ice contains flaws including frictional contacts. We aim to describe quantitatively the mechanics of those contacts, providing local physics for geophysical models. With a focus on the internal friction of ice, we review standard micro-mechanical models of friction. The solid's deformation under normal load may be ductile or elastic. The shear failure of the contact may be by ductile flow, brittle fracture, or melting and hydrodynamic lubrication. Combinations of these give a total of six rheological models. When the material under study is ice, several of the rheological parameters in the standard models are not constant, but depend on the temperature of the bulk, on the normal stress under which samples are pressed together, or on the sliding velocity and acceleration. This has the effect of making the shear stress required for sliding dependent on sliding velocity, acceleration, and temperature. In some cases, it also perturbs the exponent in the normal-stress dependence of that shear stress away from the value that applies to most materials. We unify the models by a principle of maximum displacement for normal deformation, and of minimum stress for shear failure, reducing the controversy over the mechanism of internal friction in ice to the choice of values of four parameters in a single model. The four parameters represent, for a typical asperity contact, the sliding distance required to expel melt-water, the sliding distance required to break contact, the normal strain in the asperity, and the thickness of any ductile shear zone.

Impact of variable atmospheric and oceanic form drag on simulations of Arctic sea ice

Over Arctic sea ice, pressure ridges and floe andmelt pond edges all introduce discrete obstructions to the flow of air or water past the ice and are a source of form drag. In current climate models form drag is only accounted for by tuning the air–ice and ice–ocean drag coefficients, that is, by effectively altering the roughness length in a surface drag parameterization. The existing approach of the skin drag parameter tuning is poorly constrained by observations and fails to describe correctly the physics associated with the air–ice and ocean–ice drag. Here, the authors combine recent theoretical developments to deduce the total neutral form drag coefficients from properties of the ice cover such as ice concentration, vertical extent and area of the ridges, freeboard and floe draft, and the size of floes and melt ponds. The drag coefficients are incorporated into the Los Alamos Sea Ice Model (CICE) and show the influence of the new drag parameterization on the motion and state of the ice cover, with the most noticeable being a depletion of sea ice over the west boundary of the Arctic Ocean and over the Beaufort Sea. The new parameterization allows the drag coefficients to be coupled to the sea ice state and therefore to evolve spatially and temporally. It is found that the range of values predicted for the drag coefficients agree with the range of values measured in several regions of the Arctic. Finally, the implications of the new form drag formulation for the spinup or spindown of the Arctic Ocean are discussed.

Occurrence probability, width and number of steps of cusp precipitation for fully pulsed reconnection at the dayside magnetopause

We discuss the characteristics of magnetosheath plasma precipitation in the “cusp” ionosphere for when the reconnection at the dayside magnetopause takes place only in a series of pulses. It is shown that even in this special case, the low-altitude cusp precipitation is continuous, unless the intervals between the pulses are longer than observed intervals between magnetopause flux transfer event (FTE) signatures. We use FTE observation statistics to predict, for this case of entirely pulsed reconnection, the occurrence frequency, the distribution of latitudinal widths, and the number of ion dispersion steps of the cusp precipitation for a variety of locations of the reconnection site and a range of values of the local de-Hoffman Teller velocity. It is found that the cusp occurrence frequency is comparable with observed values for virtually all possible locations of the reconnection site. The distribution of cusp width is also comparable with observations and is shown to be largely dependent on the distribution of the mean reconnection rate, but pulsing the reconnection does very slightly increase the width of that distribution compared with the steady state case. We conclude that neither cusp occurrence probability nor width can be used to evaluate the relative occurrence of reconnection behaviors that are entirely pulsed, pulsed but continuous and quasi-steady. We show that the best test of the relative frequency of these three types of reconnection is to survey the distribution of steps in the cusp ion dispersion characteristics.

Analysis of incoherent scatter radar data from non-thermal F-region plasma

A procedure is presented for fitting incoherent scatter radar data from non-thermal F-region ionospheric plasma, using theoretical spectra previously predicted. It is found that values of the shape distortion factor D∗, associated with deviations of the ion velocity distribution from a Maxwellian distribution, and ion temperatures can be deduced (the results being independent of the path of iteration) if the angle between the line-of-sight and the geomagnetic field is larger than about 15–20°. The procedure can be used with one or both of two sets of assumptions. These concern the validity of the adopted model for the line-of-sight ion velocity distribution in the one case or for the full three-dimensional ion velocity distribution function in the other. The distribution function employed was developed to describe the line-of-sight velocity distribution for large aspect angles, but both experimental data and Monte Carlo simulations indicate that the form of the field-perpendicular distribution can also describe the distribution at more general aspect angles. The assumption of this form for the line-of-sight velocity distribution at a general aspect angle enables rigorous derivation of values of the one-dimensional, line-of-sight ion temperature. With some additional assumptions (principally that the field-parallel distribution is always Maxwellian and there is a simple relationship between the ion temperature anisotropy and the distortion of the field-perpendicular distribution from a Maxwellian), fits to data for large aspect angles enable determination of line-of-sight temperatures at all aspect angles and hence, of the average ion temperature and the ion temperature anisotropy. For small aspect angles, the analysis is restricted to the determination of the line-of-sight ion temperature because the theoretical spectrum is insensitive to non-thermal effects when the plasma is viewed along directions almost parallel to the magnetic field. This limitation is expected to apply to any realistic model of the ion velocity distribution function and its consequences are discussed. Fit strategies which allow for mixed ion composition are also considered. Examples of fits to data from various EISCAT observing programmes are presented.

Deconstructing the climate change response of the Northern Hemisphere wintertime storm tracks

There are large uncertainties in the circulation response of the atmosphere to climate change. One manifestation of this is the substantial spread in projections for the extratropical storm tracks made by different state-of-the-art climate models. In this study we perform a series of sensitivity experiments, with the atmosphere component of a single climate model, in order to identify the causes of the differences between storm track responses in different models. In particular, the Northern Hemisphere wintertime storm tracks in the CMIP3 multi-model ensemble are considered. A number of potential physical drivers of storm track change are identified and their influence on the storm tracks is assessed. The experimental design aims to perturb the different physical drivers independently, by magnitudes representative of the range of values present in the CMIP3 model runs, and this is achieved via perturbations to the sea surface temperature and the sea-ice concentration forcing fields. We ask the question: can the spread of projections for the extratropical storm tracks present in the CMIP3 models be accounted for in a simple way by any of the identified drivers? The results suggest that, whilst the changes in the upper-tropospheric equator-to-pole temperature difference have an influence on the storm track response to climate change, the large spread of projections for the extratropical storm track present in the northern North Atlantic in particular is more strongly associated with changes in the lower-tropospheric equator-to-pole temperature difference.

Effects of torsional disorder on poly-para-phenylene

We study the effect of thermal disorder on the electronic structure of one-dimensional poly-para-phenylene (PPP). In a real chain the torsion angles between rings are bound to be distributed over a range of values, which depend on temperature, and thus the chain is intrinsically disordered. In this study we simulated this kind of thermally induced off-diagonal disorder through the simple Huckel method. We base our Hamiltonian on ab initio results for the effect of temperature on torsion angles, and the effect of torsion angles on the energy gap. We analyze the electronic structure of 200-monomer-long chains focusing on the density of states, and the associated localization character (measured by the inverse participation ratio). Our results contrast with the usually assumed Gaussian-shaped density of localized states for disordered systems. (C) 2009 Elsevier B.V. All rights reserved.

A quasispecies approach to the evolution of sexual replication in unicellular organisms

This study develops a simplified model describing the evolutionary dynamics of a population composed of obligate sexually and asexually reproducing, unicellular organisms. The model assumes that the organisms have diploid genomes consisting of two chromosomes, and that the sexual organisms replicate by first dividing into haploid intermediates, which then combine with other haploids, followed by the normal mitotic division of the resulting diploid into two new daughter cells. We assume that the fitness landscape of the diploids is analogous to the single-fitness-peak approach often used in single-chromosome studies. That is, we assume a master chromosome that becomes defective with just one point mutation. The diploid fitness then depends on whether the genome has zero, one, or two copies of the master chromosome. We also assume that only pairs of haploids with a master chromosome are capable of combining so as to produce sexual diploid cells, and that this process is described by second-order kinetics. We find that, in a range of intermediate values of the replication fidelity, sexually reproducing cells can outcompete asexual ones, provided the initial abundance of sexual cells is above some threshold value. The range of values where sexual reproduction outcompetes asexual reproduction increases with decreasing replication rate and increasing population density. We critically evaluate a common approach, based on a group selection perspective, used to study the competition between populations and show its flaws in addressing the evolution of sex problem.

DIAGNOSTIC TECHNIQUES OF LOCAL INFLUENCE IN SPATIAL ANALYSIS OF SOYBEAN YIELD

Modeling of spatial dependence structure, concerning geoestatistics approach, is an indispensable tool for fixing parameters that define this structure, applied on interpolation of values in places that are not sampled, by kriging techniques. However, the estimation of parameters can be greatly affected by the presence of atypical observations on sampled data. Thus, this trial aimed at using diagnostics techniques of local influence in spatial linear Gaussians models, applied at geoestatistics in order to evaluate sensitivity of maximum likelihood estimators and restrict maximum likelihood to small perturbations in these data. So, studies with simulated and experimental data were performed. Those results, obtained from the study of real data, allowed us to conclude that the presence of atypical values among the sampled data can have a strong influence on thematic maps, changing, therefore, the spatial dependence. The application of diagnostics techniques of local influence should be part of any geoestatistic analysis, ensuring that the information from thematic maps has better quality and can be used with greater security by farmers.

Influences of infrastructure and attitudes to health on value of travel time savings in bicycle journeys

In this paper we investigate how attitudes to health and exercise in connection with cycling influence the estimation of values of travel time savings in different kinds of bicycle environments (mixed traffic, bicycle lane in the road way, bicycle path next to the road, and bicycle path not in connection with the road). The results, based on two Swedish stated choice studies, suggest that the values of travel time savings are lower when cycling in better conditions. Surprisingly, the respondents do not consider cycling on a path next to the road worse than cycling on a path not in connection to the road, indicating that they do not take traffic noise and air pollution into account in their decision to cycle. No difference can be found between cycling on a road way (mixed traffic) and cycling in a bicycle lane in the road way. The results also indicate that respondents that include health aspects in their choice to cycle have lower value of travel time savings for cycling than respondents that state that health aspects are of less importance, at least when cycling on a bicycle path. The appraisals of travel time savings regarding cycling also differ a lot depending on the respondents’ alternative travel mode. The individuals who stated that they will take the car if they do not cycle have a much higher valuation of travel time savings than the persons stating public transport as the main alternative to cycling.

A simple analysis of the U.S. emission permits auctions

We examine a stylized version of EPA auctions when agents know the list of values of sellers and buyers. Sellers and buyers behave strategically. We show that there are two types of equilibria: inefficient equilibria where no goods are traded and efficient equilibria where alI exchange occurs at a uniform price. We also provide examples of the EPA auction game under incomplete information when the uniform price equilibrium holds and when it does not hold. When the uniform price equilibrium holds, sellers shade their bids up and buyers shade their bids down. In the example where the uniform price equilibrium does not hold, both buyers and sellers shade their bids down in an equilibrium.

Effect of haematocrit on pump and kidney dysfunction after myocardial revascularization

Objective - Kidney dysfunction is a common complication after cardiac surgery. It occurs in 7 to 31% of the patients. The lowest haematocrit after cardiopulmonary bypass surgery (LHCT) has been identified as a risk factor for kidney dysfunction after cardiac surgery. The aim of this study is to determine whether different levels of haematocrit during cardiopulmonary bypass surgery are related to kidney dysfunction.Methods and results-A prospective study was conducted on consecutive adult patients undergoing myocardial revascularization. Preoperative renal function was assessed by baseline serum creatinine level (CrPre). Peak postoperative creatinine (CrPost) was defined as the highest daily in-hospital postoperative value. Peak fractional change in creatinine (% Delta Cr) was defined as the difference between the CrPre and CrePost represented as a percentage of the preoperative value. The LHTC was defined as the lowest recorded haematocrit prior to weaning from the initial pump run. A category variable was created for haematocrit based on the distribution of values. The category variable had the following cut-off points: less than 23%, 23.1 to 28% and greater than 28.1 %. Lowest haematocrit (26.62 +/- 4.15%), CPB (74.71 +/- 24.90 min), CrPre (1.23 +/- 0.37 mg/dl) and highest CrPost (1.52 +/- 0.47 mg/dl) data varied in near-normal fashion. Statistical significance has been observed in the < 23% lowest haematocrit group (CrIPOD and Cr5POD; P = 0.006) and the 23.1 28% lowest haematocrit level group (CrPre and Cr2POD; P = 0.047). CrPre and Cr5POD did not differ between groups (P > 0.05). The multiple linear regression model confirmed that the determinants for higher %Delta Cr were age, body surface area and preoperative serum creatinine level.Conclusion - The LHTC was not identified as a risk factor for kidney dysfunction after myocardial revascularization.

Histomorphometric evaluation of periodontal compression and tension sides during orthodontic tooth movement in rats

OBJECTIVE: The purpose of this study was to evaluate the thickness of the periodontal ligament of rat molars during orthodontic tooth movement (OTM). METHODS: Thirty Wistar rats were divided into three groups of 10 animals each: GI, GII and GIII and the mice were euthanized at 7, 14 and 21 days, respectively. Experimental subjects were compared to their respective controls by the Mann-Whitney test. Comparison of values between compression and tension sides were performed during the same and different time periods through Analysis of Variance (ANOVA), Kruskal-Wallis test and, subsequently, Tukey's test. RESULTS: Groups GI and GII showed decreased PDL size in the apical regions of the mesiobuccal root and in the cervical region of the distobuccal root. There was also an increased PDL in the cervical regions of the mesiobuccal root, apical region of the distobuccal root and middle region of both roots. CONCLUSION: The reduction and increase in PDL size were seen in the same root, which characterizes tooth inclination. The apical, middle and cervical regions were compared with one another in each time period and at three times: 7, 14 and 21 days. They were also compared in each region, confirming a tipping movement in GI and GII and a gradual decreased intensity between GI to GII, reaching normal dimension in GIII.

The Origin of Nonmonotonic Complex Behavior and the Effects of Nonnative Interactions on the Diffusive Properties of Protein Folding

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Dynamics of Bose-Einstein condensates in cigar-shaped traps

The Gross-Pitaevskii equation for a Bose-Einstein condensate confined in an elongated cigar-shaped trap is reduced to an effective system of nonlinear equations depending on only one space coordinate along the trap axis. The radial distribution of the condensate density and its radial velocity are approximated by Gaussian functions with real and imaginary exponents, respectively, with parameters depending on the axial coordinate and time. The effective one-dimensional system is applied to a description of the ground state of the condensate, to dark and bright solitons, to the sound and radial compression waves propagating in a dense condensate, and to weakly nonlinear waves in repulsive condensate. In the low-density limit our results reproduce the known formulas. In the high-density case our description of solitons goes beyond the standard approach based on the nonlinear Schrodinger equation. The dispersion relations for the sound and radial compression waves are obtained in a wide region of values of the condensate density. The Korteweg-de Vries equation for weakly nonlinear waves is derived and the existence of bright solitons on a constant background is predicted for a dense enough condensate with a repulsive interaction between the atoms.