Application of standard and refined heat balance integral methods to one-dimensional Stefan problems

The work in this paper concerns the study of conventional and refined heat balance integral methods for a number of phase change problems. These include standard test problems, both with one and two phase changes, which have exact solutions to enable us to test the accuracy of the approximate solutions. We also consider situations where no analytical solution is available and compare these to numerical solutions. It is popular to use a quadratic profile as an approximation of the temperature, but we show that a cubic profile, seldom considered in the literature, is far more accurate in most circumstances. In addition, the refined integral method can give greater improvement still and we develop a variation on this method which turns out to be optimal in some cases. We assess which integral method is better for various problems, showing that it is largely dependent on the specified boundary conditions.

Weak approximations. A Malliavin calculus approach

We introduce a variation of the proof for weak approximations that issuitable for studying the densities of stochastic processes which areevaluations of the flow generated by a stochastic differential equation on a random variable that maybe anticipating. Our main assumption is that the process and the initial random variable have to be smooth in the Malliavin sense. Furthermore if the inverse of the Malliavin covariance matrix associated with the process under consideration is sufficiently integrable then approximations fordensities and distributions can also be achieved. We apply theseideas to the case of stochastic differential equations with boundaryconditions and the composition of two diffusions.

Reply to Comment on 'Two-finger selection theory in the Saffman-Taylor problem

We clarify the meaning of the results of Phys. Rev. E 60, R5013 (1999). We discuss the use and implications of periodic boundary conditions, as opposed to rigid-wall ones. We briefly argue that the solutions of the paper above are physically relevant as part of a more general issue, namely the possible generalization to dynamics, of the microscopic solvability scenario of selection.

Monte Carlo simulation of kilovolt electron transport in solids

A Monte Carlo procedure to simulate the penetration and energy loss of low¿energy electron beams through solids is presented. Elastic collisions are described by using the method of partial waves for the screened Coulomb field of the nucleus. The atomic charge density is approximated by an analytical expression with parameters determined from the Dirac¿Hartree¿Fock¿Slater self¿consistent density obtained under Wigner¿Seitz boundary conditions in order to account for solid¿state effects; exchange effects are also accounted for by an energy¿dependent local correction. Elastic differential cross sections are then easily computed by combining the WKB and Born approximations to evaluate the phase shifts. Inelastic collisions are treated on the basis of a generalized oscillator strength model which gives inelastic mean free paths and stopping powers in good agreement with experimental data. This scattering model is accurate in the energy range from a few hundred eV up to about 50 keV. The reliability of the simulation method is analyzed by comparing simulation results and experimental data from backscattering and transmission measurements.

Bifurcations and chaos in single-roll natural convection with low Prandtl number

Convective flows of a small Prandtl number fluid contained in a two-dimensional cavity subject to a lateral thermal gradient are numerically studied by using different techniques. The aspect ratio (length to height) is kept at around 2. This value is found optimal to make the flow most unstable while keeping the basic single-roll structure. Two cases of thermal boundary conditions on the horizontal plates are considered: perfectly conducting and adiabatic. For increasing Rayleigh numbers we find a transition from steady flow to periodic oscillations through a supercritical Hopf bifurcation that maintains the centrosymmetry of the basic circulation. For a Rayleigh number of about ten times that of the Hopf bifurcation the system initiates a complex scenario of bifurcations. In the conductive case these include a quasiperiodic route to chaos. In the adiabatic one the dynamics is dominated by the interaction of two Neimark-Sacker bifurcations of the basic periodic solutions, leading to the stable coexistence of three incommensurate frequencies, and finally to chaos. In all cases, the complex time-dependent behavior does not break the basic, single-roll structure.

Properties of resonant activation phenomena

The phenomenon of resonant activation of a Brownian particle over a fluctuating barrier is revisited. We discuss the important distinctions between barriers that can fluctuate among up and down configurations, and barriers that are always up but that can fluctuate among different heights. A resonance as a function of the barrier fluctuation rate is found in both cases, but the nature and physical description of these resonances is quite distinct. The nature of the resonances, the physical basis for the resonant behavior, and the importance of boundary conditions are discussed in some detail. We obtain analytic expressions for the escape time over the barrier that explicitly capture the minima as a function of the barrier fluctuation rate, and show that our analytic results are in excellent agreement with numerical results.

Algebraic decay of velocity fluctuations near a wall

Computer simulations of the dynamics of a colloidal particle suspended in a fluid confined by an interface show that the asymptotic decay of the velocity correlation functions is algebraic. The exponents of the long-time tails depend on the direction of motion of the particle relative to the surface, as well as on the specific nature of the boundary conditions. In particular, we find that for the angular velocity correlation function, the decay in the presence of a slip surface is faster than the one corresponding to a stick one. An intuitive picture is introduced to explain the various long-time tails, and the simulations are compared with theoretical expressions where available.

Spectral analysis and limit behaviours in a spring-mass system

We consider a model for a damped spring-mass system that is a strongly damped wave equation with dynamic boundary conditions. In a previous paper we showed that for some values of the parameters of the model, the large time behaviour of the solutions is the same as for a classical spring-mass damper ODE. Here we use spectral analysis to show that for other values of the parameters, still of physical relevance and related to the effect of the spring inner viscosity, the limit behaviours are very different from that classical ODE

Microwave spectrometry for the Evaluation of the structural integrity of metallic stents

Purpose: To assess the feasibility of a method based on microwave spectrometry to detect structural distortions of metallic stents in open air conditions and envisage the prospects of this approach toward possible medical applicability for the evaluation of implanted stents. Methods: Microwave absorbance spectra between 2.0 and 18.0 GHz were acquired in open air for the characterization of a set of commercial stents using a specifically design setup. Rotating each sample over 360º, 2D absorbance diagrams were generated as a function of frequency and rotation angle. To check our approach for detecting changes in stent length (fracture) and diameter (recoil), two specific tests were performed in open air. Finally, with a few adjustments, this same system provides 2D absorbance diagrams of stents immersed in a water-based phantom, this time over a bandwidth ranging from 0.2 to 1.8 GHz. Results: The authors show that metallic stents exhibit characteristic resonant frequencies in their microwave absorbance spectra in open air which depend on their length and, as a result, may reflect the occurrence of structural distortions. These resonances can be understood considering that such devices behave like dipole antennas in terms of microwave scattering. From fracture tests, the authors infer that microwave spectrometry provides signs of presence of Type I to Type IV stent fractures and allows in particular a quantitative evaluation of Type III and Type IV fractures. Recoil tests show that microwave spectrometry seems able to provide some quantitative assessment of diametrical shrinkage, but only if it involves longitudinal shortening. Finally, the authors observe that the resonant frequencies of stents placed inside the phantom shift down with respect to the corresponding open air frequencies, as it should be expected considering the increase of dielectric permittivity from air to water. Conclusions: The evaluation of stent resonant frequencies provided by microwave spectrometry allows detection and some quantitative assessment of stent fracture and recoil in open air conditions. Resonances of stents immersed in water can be also detected and their characteristic frequencies are in good agreement with theoretical estimates. Although these are promising results, further verifica tion in a more relevant phantom is required in order to foresee the real potential of this approach.

Cosmic censorship in two-dimensional gravity

A weak version of the cosmic censorship hypothesis is implemented as a set of boundary conditions on exact semiclassical solutions of two-dimensional dilaton gravity. These boundary conditions reflect low-energy matter from the strong coupling region and they also serve to stabilize the vacuum of the theory against decay into negative energy states. Information about low-energy incoming matter can be recovered in the final state but at high energy black holes are formed and inevitably lead to information loss at the semiclassical level.

End point of Hawking radiation

The formation and semiclassical evaporation of two-dimensional black holes is studied in an exactly solvable model. Above a certain threshold energy flux, collapsing matter forms a singularity inside an apparent horizon. As the black hole evaporates the apparent horizon recedes and meets the singularity in a finite proper time. The singularity emerges naked, and future evolution of the geometry requires boundary conditions to be imposed there. There is a natural choice of boundary conditions which matches the evaporated black hole solution onto the linear dilaton vacuum. Below the threshold energy flux no horizon forms and boundary conditions can be imposed where infalling matter is reflected from a timelike boundary. All information is recovered at spatial infinity in this case.

Model of black hole evolution

We investigate a model where the quantum dynamics of black hole evaporation is determined by imposing a boundary on the apparent horizon with suitable boundary conditions. An unconventional scenario for the evolution emerges: only an insignificant fraction of energy of order (mG)-1 is radiated out; the outgoing wave carries a very small part of the quantum-mechanical information of the collapsed body, the bulk of the information remaining in the final stable black hole geometry.

Microwave spectrometry for the Evaluation of the structural integrity of metallic stents

Purpose: To assess the feasibility of a method based on microwave spectrometry to detect structural distortions of metallic stents in open air conditions and envisage the prospects of this approach toward possible medical applicability for the evaluation of implanted stents. Methods: Microwave absorbance spectra between 2.0 and 18.0 GHz were acquired in open air for the characterization of a set of commercial stents using a specifically design setup. Rotating each sample over 360º, 2D absorbance diagrams were generated as a function of frequency and rotation angle. To check our approach for detecting changes in stent length (fracture) and diameter (recoil), two specific tests were performed in open air. Finally, with a few adjustments, this same system provides 2D absorbance diagrams of stents immersed in a water-based phantom, this time over a bandwidth ranging from 0.2 to 1.8 GHz. Results: The authors show that metallic stents exhibit characteristic resonant frequencies in their microwave absorbance spectra in open air which depend on their length and, as a result, may reflect the occurrence of structural distortions. These resonances can be understood considering that such devices behave like dipole antennas in terms of microwave scattering. From fracture tests, the authors infer that microwave spectrometry provides signs of presence of Type I to Type IV stent fractures and allows in particular a quantitative evaluation of Type III and Type IV fractures. Recoil tests show that microwave spectrometry seems able to provide some quantitative assessment of diametrical shrinkage, but only if it involves longitudinal shortening. Finally, the authors observe that the resonant frequencies of stents placed inside the phantom shift down with respect to the corresponding open air frequencies, as it should be expected considering the increase of dielectric permittivity from air to water. Conclusions: The evaluation of stent resonant frequencies provided by microwave spectrometry allows detection and some quantitative assessment of stent fracture and recoil in open air conditions. Resonances of stents immersed in water can be also detected and their characteristic frequencies are in good agreement with theoretical estimates. Although these are promising results, further verifica tion in a more relevant phantom is required in order to foresee the real potential of this approach.

Modelo de comportamiento en sevicio para estructuras reticulares de hormigón. Parte I: Descripción del modelo y ejemplos de aplicación

A general method for instantaneous and time-dependent serviceability analysis of plane concrete frames is presented. The methodology is based in an extension of the classic matrix formulation for bars. The main aspects influencing the behaviour of the structural concrete are considered: cracking, creep, shrinkage or prestress losses. To simulate the effect of cracking a smeared model (developed in Part II) based on the modification of the tensile branch of the concrete stress-strain relationship is adopted. The general approach considered permits the application to different materials and constitutive laws. Sequential construction (sectional and structural), incorporation of reinforcement, consideration of the loads history; placing and removing shores, and restraining or releasing in boundary conditions are considered. Some examples are included to highlight the capabilities of the model

Leaving the Dark ages with AMIGA

We present an Analytic Model of Intergalactic-medium and GAlaxy (AMIGA) evolution since the dark ages. AMIGA is in the spirit of the popular semi-analytic models of galaxy formation, although it does not use halo merger trees but interpolates halo properties in grids that are progressively built. This strategy is less memory-demanding and allows one to start modeling at sufficiently high redshifts and low halo masses to have trivial boundary conditions. The number of free parameters is minimized by making a causal connection between physical processes usually treated as independent of each other, which leads to more reliable predictions. However, the strongest points of AMIGA are the following: (1) the inclusion of molecular cooling and metal-poor, population III (Pop III) stars with the most dramatic feedback and (2) accurate follow up of the temperature and volume filling factor of neutral, singly ionized, and doubly ionized regions, taking into account the distinct halo mass functions in those environments. We find the following general results. Massive Pop III stars determine the intergalactic medium metallicity and temperature, and the growth of spheroids and disks is self-regulated by that of massive black holes (MBHs) developed from the remnants of those stars. However, the properties of normal galaxies and active galactic nuclei appear to be quite insensitive to Pop III star properties due to the much higher yield of ordinary stars compared to Pop III stars and the dramatic growth of MBHs when normal galaxies begin to develop, which cause the memory loss of the initial conditions.