On The Spot Damage Detection Methodology for Highway Bridges, July 2010

Vibration-based damage identification (VBDI) techniques have been developed in part to address the problems associated with an aging civil infrastructure. To assess the potential of VBDI as it applies to highway bridges in Iowa, three applications of VBDI techniques were considered in this study: numerical simulation, laboratory structures, and field structures. VBDI techniques were found to be highly capable of locating and quantifying damage in numerical simulations. These same techniques were found to be accurate in locating various types of damage in a laboratory setting with actual structures. Although there is the potential for these techniques to quantify damage in a laboratory setting, the ability of the methods to quantify low-level damage in the laboratory is not robust. When applying these techniques to an actual bridge, it was found that some traditional applications of VBDI methods are capable of describing the global behavior of the structure but are most likely not suited for the identification of typical damage scenarios found in civil infrastructure. Measurement noise, boundary conditions, complications due to substructures and multiple material types, and transducer sensitivity make it very difficult for present VBDI techniques to identify, much less quantify, highly localized damage (such as small cracks and minor changes in thickness). However, while investigating VBDI techniques in the field, it was found that if the frequency-domain response of the structure can be generated from operating traffic load, the structural response can be animated and used to develop a holistic view of the bridge’s response to various automobile loadings. By animating the response of a field bridge, concrete cracking (in the abutment and deck) was correlated with structural motion and problem frequencies (i.e., those that cause significant torsion or tension-compression at beam ends) were identified. Furthermore, a frequency-domain study of operational traffic was used to identify both common and extreme frequencies for a given structure and loading. Common traffic frequencies can be compared to problem frequencies so that cost-effective, preventative solutions (either structural or usage-based) can be developed for a wide range of IDOT bridges. Further work should (1) perfect the process of collecting high-quality operational frequency response data; (2) expand and simplify the process of correlating frequency response animations with damage; and (3) develop efficient, economical, preemptive solutions to common damage types.

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.

Revisiting the S-index for soil physical quality and its use in Brazil

The S-index was introduced in 2004 in a publication by A.R. Dexter. S was proposed as an indicator of soil physical quality. A critical value delimiting soils with rich and poor physical quality was proposed. At present, Brazil is world leader in citations of Dexter's publication. In this publication the S-theory is mathematically revisited and extended. It is shown that S is mathematically correlated to bulk density and total porosity. As an absolute indicator, the value of S alone has proven to be incapable of predicting soil physical quality. The critical value does not always hold under boundary conditions described in the literature. This is to be expected because S is a static parameter, therefore implicitly unable to describe dynamic processes. As a relative indicator of soil physical quality, the S-index has no additional value over bulk density or total porosity. Therefore, in the opinion of the author, the fact that bulk density or total porosity are much more easily determined than the water retention curve for obtaining S disqualifies S as an advantageous indicator of relative soil physical quality. Among the several equations available for the fitting of water retention curves, the Groenevelt-Grant equation is preferable for use with S since one of its parameters and S are linearly correlated. Since efforts in soil physics research have the purpose of describing dynamic processes, it is the author's opinion that these efforts should shift towards mechanistic soil physics as opposed to the search for empirical correlations like S which, at present, represents far more than its reasonable share of soil physics in Brazil.

Accurate and efficient simulation of multiphase flow in a heterogeneous reservoir by using error estimate and control in the multiscale finite-volume framework

The multiscale finite-volume (MSFV) method is designed to reduce the computational cost of elliptic and parabolic problems with highly heterogeneous anisotropic coefficients. The reduction is achieved by splitting the original global problem into a set of local problems (with approximate local boundary conditions) coupled by a coarse global problem. It has been shown recently that the numerical errors in MSFV results can be reduced systematically with an iterative procedure that provides a conservative velocity field after any iteration step. The iterative MSFV (i-MSFV) method can be obtained with an improved (smoothed) multiscale solution to enhance the localization conditions, with a Krylov subspace method [e.g., the generalized-minimal-residual (GMRES) algorithm] preconditioned by the MSFV system, or with a combination of both. In a multiphase-flow system, a balance between accuracy and computational efficiency should be achieved by finding a minimum number of i-MSFV iterations (on pressure), which is necessary to achieve the desired accuracy in the saturation solution. In this work, we extend the i-MSFV method to sequential implicit simulation of time-dependent problems. To control the error of the coupled saturation/pressure system, we analyze the transport error caused by an approximate velocity field. We then propose an error-control strategy on the basis of the residual of the pressure equation. At the beginning of simulation, the pressure solution is iterated until a specified accuracy is achieved. To minimize the number of iterations in a multiphase-flow problem, the solution at the previous timestep is used to improve the localization assumption at the current timestep. Additional iterations are used only when the residual becomes larger than a specified threshold value. Numerical results show that only a few iterations on average are necessary to improve the MSFV results significantly, even for very challenging problems. Therefore, the proposed adaptive strategy yields efficient and accurate simulation of multiphase flow in heterogeneous porous media.

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.

BEYOND THE “LEAST LIMITING WATER RANGE”: RETHINKING SOIL PHYSICS RESEARCH IN BRAZIL

As opposed to objective definitions in soil physics, the subjective term “soil physical quality” is increasingly found in publications in the soil physics area. A supposed indicator of soil physical quality that has been the focus of attention, especially in the Brazilian literature, is the Least Limiting Water Range (RLL), translated in Portuguese as "Intervalo Hídrico Ótimo" or IHO. In this paper the four limiting water contents that define RLLare discussed in the light of objectively determinable soil physical properties, pointing to inconsistencies in the RLLdefinition and calculation. It also discusses the interpretation of RLL as an indicator of crop productivity or soil physical quality, showing its inability to consider common phenological and pedological boundary conditions. It is shown that so-called “critical densities” found by the RLL through a commonly applied calculation method are questionable. Considering the availability of robust models for agronomy, ecology, hydrology, meteorology and other related areas, the attractiveness of RLL as an indicator to Brazilian soil physicists is not related to its (never proven) effectiveness, but rather to the simplicity with which it is dealt. Determining the respective limiting contents in a simplified manner, relegating the study or concern on the actual functioning of the system to a lower priority, goes against scientific construction and systemic understanding. This study suggests a realignment of the research in soil physics in Brazil with scientific precepts, towards mechanistic soil physics, to replace the currently predominant search for empirical correlations below the state of the art of soil physics.

Genome-wide analysis of Sphingomonas wittichii RW1 behaviour during inoculation and growth in contaminated sand.

The efficacy of inoculation of single pure bacterial cultures into complex microbiomes, for example, in order to achieve increased pollutant degradation rates in contaminated material (that is, bioaugmentation), has been frustrated by insufficient knowledge on the behaviour of the inoculated bacteria under the specific abiotic and biotic boundary conditions. Here we present a comprehensive analysis of genome-wide gene expression of the bacterium Sphingomonas wittichii RW1 in contaminated non-sterile sand, compared with regular suspended batch growth in liquid culture. RW1 is a well-known bacterium capable of mineralizing dibenzodioxins and dibenzofurans. We tested the reactions of the cells both during the immediate transition phase from liquid culture to sand with or without dibenzofuran, as well as during growth and stationary phase in sand. Cells during transition show stationary phase characteristics, evidence for stress and for nutrient scavenging, and adjust their primary metabolism if they were not precultured on the same contaminant as found in the soil. Cells growing and surviving in sand degrade dibenzofuran but display a very different transcriptome signature as in liquid or in liquid culture exposed to chemicals inducing drought stress, and we obtain evidence for numerous 'soil-specific' expressed genes. Studies focusing on inoculation efficacy should test behaviour under conditions as closely as possible mimicking the intended microbiome conditions.

Pore-scale modeling of two-phase flow instabilities in porous media

Les problèmes d'écoulements multiphasiques en média poreux sont d'un grand intérêt pour de nombreuses applications scientifiques et techniques ; comme la séquestration de C02, l'extraction de pétrole et la dépollution des aquifères. La complexité intrinsèque des systèmes multiphasiques et l'hétérogénéité des formations géologiques sur des échelles multiples représentent un challenge majeur pour comprendre et modéliser les déplacements immiscibles dans les milieux poreux. Les descriptions à l'échelle supérieure basées sur la généralisation de l'équation de Darcy sont largement utilisées, mais ces méthodes sont sujettes à limitations pour les écoulements présentant de l'hystérèse. Les avancées récentes en terme de performances computationnelles et le développement de méthodes précises pour caractériser l'espace interstitiel ainsi que la distribution des phases ont favorisé l'utilisation de modèles qui permettent une résolution fine à l'échelle du pore. Ces modèles offrent un aperçu des caractéristiques de l'écoulement qui ne peuvent pas être facilement observées en laboratoire et peuvent être utilisé pour expliquer la différence entre les processus physiques et les modèles à l'échelle macroscopique existants. L'objet premier de la thèse se porte sur la simulation numérique directe : les équations de Navier-Stokes sont résolues dans l'espace interstitiel et la méthode du volume de fluide (VOF) est employée pour suivre l'évolution de l'interface. Dans VOF, la distribution des phases est décrite par une fonction fluide pour l'ensemble du domaine et des conditions aux bords particulières permettent la prise en compte des propriétés de mouillage du milieu poreux. Dans la première partie de la thèse, nous simulons le drainage dans une cellule Hele-Shaw 2D avec des obstacles cylindriques. Nous montrons que l'approche proposée est applicable même pour des ratios de densité et de viscosité très importants et permet de modéliser la transition entre déplacement stable et digitation visqueuse. Nous intéressons ensuite à l'interprétation de la pression capillaire à l'échelle macroscopique. Nous montrons que les techniques basées sur la moyenne spatiale de la pression présentent plusieurs limitations et sont imprécises en présence d'effets visqueux et de piégeage. Au contraire, une définition basée sur l'énergie permet de séparer les contributions capillaires des effets visqueux. La seconde partie de la thèse est consacrée à l'investigation des effets d'inertie associés aux reconfigurations irréversibles du ménisque causé par l'interface des instabilités. Comme prototype pour ces phénomènes, nous étudions d'abord la dynamique d'un ménisque dans un pore angulaire. Nous montrons que, dans un réseau de pores cubiques, les sauts et reconfigurations sont si fréquents que les effets d'inertie mènent à différentes configurations des fluides. A cause de la non-linéarité du problème, la distribution des fluides influence le travail des forces de pression, qui, à son tour, provoque une chute de pression dans la loi de Darcy. Cela suggère que ces phénomènes devraient être pris en compte lorsque que l'on décrit l'écoulement multiphasique en média poreux à l'échelle macroscopique. La dernière partie de la thèse s'attache à démontrer la validité de notre approche par une comparaison avec des expériences en laboratoire : un drainage instable dans un milieu poreux quasi 2D (une cellule Hele-Shaw avec des obstacles cylindriques). Plusieurs simulations sont tournées sous différentes conditions aux bords et en utilisant différents modèles (modèle intégré 2D et modèle 3D) afin de comparer certaines quantités macroscopiques avec les observations au laboratoire correspondantes. Malgré le challenge de modéliser des déplacements instables, où, par définition, de petites perturbations peuvent grandir sans fin, notre approche numérique apporte de résultats satisfaisants pour tous les cas étudiés. - Problems involving multiphase flow in porous media are of great interest in many scientific and engineering applications including Carbon Capture and Storage, oil recovery and groundwater remediation. The intrinsic complexity of multiphase systems and the multi scale heterogeneity of geological formations represent the major challenges to understand and model immiscible displacement in porous media. Upscaled descriptions based on generalization of Darcy's law are widely used, but they are subject to several limitations for flow that exhibit hysteric and history- dependent behaviors. Recent advances in high performance computing and the development of accurate methods to characterize pore space and phase distribution have fostered the use of models that allow sub-pore resolution. These models provide an insight on flow characteristics that cannot be easily achieved by laboratory experiments and can be used to explain the gap between physical processes and existing macro-scale models. We focus on direct numerical simulations: we solve the Navier-Stokes equations for mass and momentum conservation in the pore space and employ the Volume Of Fluid (VOF) method to track the evolution of the interface. In the VOF the distribution of the phases is described by a fluid function (whole-domain formulation) and special boundary conditions account for the wetting properties of the porous medium. In the first part of this thesis we simulate drainage in a 2-D Hele-Shaw cell filled with cylindrical obstacles. We show that the proposed approach can handle very large density and viscosity ratios and it is able to model the transition from stable displacement to viscous fingering. We then focus on the interpretation of the macroscopic capillary pressure showing that pressure average techniques are subject to several limitations and they are not accurate in presence of viscous effects and trapping. On the contrary an energy-based definition allows separating viscous and capillary contributions. In the second part of the thesis we investigate inertia effects associated with abrupt and irreversible reconfigurations of the menisci caused by interface instabilities. As a prototype of these phenomena we first consider the dynamics of a meniscus in an angular pore. We show that in a network of cubic pores, jumps and reconfigurations are so frequent that inertia effects lead to different fluid configurations. Due to the non-linearity of the problem, the distribution of the fluids influences the work done by pressure forces, which is in turn related to the pressure drop in Darcy's law. This suggests that these phenomena should be taken into account when upscaling multiphase flow in porous media. The last part of the thesis is devoted to proving the accuracy of the numerical approach by validation with experiments of unstable primary drainage in a quasi-2D porous medium (i.e., Hele-Shaw cell filled with cylindrical obstacles). We perform simulations under different boundary conditions and using different models (2-D integrated and full 3-D) and we compare several macroscopic quantities with the corresponding experiment. Despite the intrinsic challenges of modeling unstable displacement, where by definition small perturbations can grow without bounds, the numerical method gives satisfactory results for all the cases studied.

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

Low-cost MR-compatible moving heart phantom.

For the development and evaluation of cardiac magnetic resonance (MR) imaging sequences and methodologies, the availability of a periodically moving phantom to model respiratory and cardiac motion would be of substantial benefit. Given the specific physical boundary conditions in an MR environment, the choice of materials and power source of such phantoms is heavily restricted. Sophisticated commercial solutions are available; however, they are often relatively costly and user-specific modifications may not easily be implemented. We therefore sought to construct a low-cost MR-compatible motion phantom that could be easily reproduced and had design flexibility. A commercially available K'NEX construction set (Hyper Space Training Tower, K'NEX Industries, Inc., Hatfield, PA) was used to construct a periodically moving phantom head. The phantom head performs a translation with a superimposed rotation, driven by a motor over a 2-m rigid rod. To synchronize the MR data acquisition with phantom motion (without introducing radiofrequency-related image artifacts), a fiberoptic control unit generates periodic trigger pulses synchronized to the phantom motion. Total material costs of the phantom are US$ < 200.00, and a total of 80 man-hours were required to design and construct the original phantom. With schematics of the present solution, the phantom reproduction may be achieved in approximately 15 man-hours. The presented MR-compatible periodically moving phantom can easily be reproduced, and user-specific modifications may be implemented. Such an approach allows a detailed investigation of motion-related phenomena in MR images.

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.

The Dutung-Thaktote extensional fault zone and nappe structures documented by illite crystallinity and clay-mineral paragenesis in the Tethys Himalaya between Spiti river and Tso Morari, NW India

Numerous measurements by XRD of the Scherrer width at half-peak height (001 reflection of illite), coupled with analyses of clay-size assemblages, provide evidence for strong variations in the conditions of low temperature metamorphism in the Tethyan Himalaya metasediments between the Spiti river and the Tso Morari. Three sectors can be distinguished along the Spiti river-Tso Morari transect. In the SW, the Takling and Parang La area is characterised by a metamorphism around anchizone-epizone boundary conditions. Further north, in the Dutung area, the metamorphic grade abruptly decreases to weak diagenesis, with the presence of mixed-layered clay phases. At the end of the profile towards the NE, a progressive metamorphic increase up to greenschist facies is recorded, marked by the appearance of biotite and chloritoid. The combination of these data with the structural. observations permits to propose that a nappe stack has been crosscut by the younger Dutung-Thaktote extensional fault zone (DTFZ). The change in metamorphism across this zone helps to assess the displacements which occurred during synorogenic extension. In the SW and NE parts of the studied transect, a burial of 12 km has been estimated, assuming a geothermal gradient of 25 degrees C/km. In the SW part, this burial is due to the juxtaposition of the Shikar Beh and Mata nappes and in the NE part, solely to burial beneath the Mata nappe. In the central part of the profile, the effect of the DTFZ is to bring down diagenetic sediments in-between the two aforesaid metamorphic zones. The offset along the Dutung-Thaktote normal faults is estimated at 16 km.