Direct numerical simulations of interface dynamics to link capillary pressure and total surface energy

The flow of two immiscible fluids through a porous medium depends on the complex interplay between gravity, capillarity, and viscous forces. The interaction between these forces and the geometry of the medium gives rise to a variety of complex flow regimes that are difficult to describe using continuum models. Although a number of pore-scale models have been employed, a careful investigation of the macroscopic effects of pore-scale processes requires methods based on conservation principles in order to reduce the number of modeling assumptions. In this work we perform direct numerical simulations of drainage by solving Navier-Stokes equations in the pore space and employing the Volume Of Fluid (VOF) method to track the evolution of the fluid-fluid interface. After demonstrating that the method is able to deal with large viscosity contrasts and model the transition from stable flow to viscous fingering, we focus on the macroscopic capillary pressure and we compare different definitions of this quantity under quasi-static and dynamic conditions. We show that the difference between the intrinsic phase-average pressures, which is commonly used as definition of Darcy-scale capillary pressure, is subject to several limitations and it is not accurate in presence of viscous effects or trapping. In contrast, a definition based on the variation of the total surface energy provides an accurate estimate of the macroscopic capillary pressure. This definition, which links the capillary pressure to its physical origin, allows a better separation of viscous effects and does not depend on the presence of trapped fluid clusters.

Design, Fabrication, and Operation of the Pavement Crack Monitor Gauge, HR-1054, 1990

A mechanical gauge was developed to monitor the movement of crack or joint openings in portland cement concrete structures, in general, and portland cement concrete pavements in particular. Designed to be inexpensive and simple to operate, this gauge is capable of recording maximum, minimum, and instantaneous crack or joint openings. Specific recommendations were made for recording minimum and maximum pavement temperature over the monitoring period. The report was written as a set of guidelines for design, fabrication, installation, and operation of the gauge as well as the temperature measuring device.

Effect of Density on Marshall Stability of Hot-Mix Asphaltic Concrete, MLR-82-01, 1982

The Iowa D.O.T. specifications do not require 100 percent of 50 blow Marshall density (generally 94%) for field compaction. However, stabilities are determined in the Laboratory on specimens compacted to 100 percent of Marshall density. The purpose of this study is to determine the stabilities of specimens compacted to various densities which are below 100 percent of the 50 blow Marshall density.

Investigation of Video Crack and Patch Surveying for Pavement Management, MLR-92-01, 1992

The Iowa DOT has been using the AASHTO Present Serviceability Index (PSI) rating procedure since 1968 to rate the condition of pavement sections. A ride factor and a cracking and patching factor make up the PSI value. Crack and patch surveys have been done by sending crews out to measure and record the distress. Advances in video equipment and computers make it practical to videotape roads and do the crack and patch measurements in the office. The objective of the study was to determine the feasibility of converting the crack and patch survey operation to a video recording system with manual post processing. The summary and conclusions are as follows: Video crack and patch surveying is a feasible alternative to the current crack and patch procedure. The cost per mile should be about 25 percent less than the current procedure. More importantly, the risk of accidents is reduced by getting the people and vehicles off the roadway and shoulder. Another benefit is the elimination of the negative public perceptions of the survey crew on the shoulder.

Rehabilitation of Concrete Pavements Utilizing Rubblization and Crack and Seat Methods (Phase II): Performance Evaluation of Rubblized Pavements in Iowa, TR-550, 2008

This Phase II follow-up study of IHRB Project TR-473 focused on the performance evaluation of rubblized pavements in Iowa. The primary objective of this study was to evaluate the structural condition of existing rubblized concrete pavements across Iowa through Falling Weight Deflectometer (FWD) tests, Dynamic Cone Penetrometer (DCP) tests, visual pavement distress surveys, etc. Through backcalculation of FWD deflection data using the Iowa State University's advanced layer moduli backcalculation program, the rubblized layer moduli were determined for various projects and compared with each other for correlating with the long-term pavement performance. The AASHTO structural layer coefficient for rubblized layer was also calculated using the rubblized layer moduli. To validate the mechanistic-empirical (M-E) hot mix asphalt (HMA) overlay thickness design procedure developed during the Phase I study, the actual HMA overlay thicknesses from the rubblization projects were compared with the predicted thicknesses obtained from the design software. The results of this study show that rubblization is a valid option to use in Iowa in the rehabilitation of portland cement concrete pavements provided the foundation is strong enough to support construction operations during the rubblization process. The M-E structural design methodology developed during Phase I can estimate the HMA overlay thickness reasonably well to achieve long-lasting performance of HMA pavements. The rehabilitation strategy is recommended for continued use in Iowa under those conditions conducive for rubblization.

Investigation of Electromagnetic Gauges for Determination of In-Place Density of HMA Pavements, TR-576, 2009

This study evaluated the use of electromagnetic gauges to determine the adjusted densities of HMA pavements. Field measurements were taken with two electromagnetic gauges, the Pavement Quality Indicator (PQI) 301 and the Pavetracker Plus 2701B. Seven projects were included in the study with 3 to 5 consecutive paving days. For each day/lot 20 randomly selected locations were tested along with seven core locations. The analysis of PaveTracker and PQI density consisted of determining which factors are statistically significant, and core density residuals and a regression analysis of core as a function of PaveTracker and PQI readings. The following key conclusions can be stated: 1. Core density, traffic and binder content were all found to be significant for both electromagnetic gauges studied, 2. Core density residuals are normally distributed and centered at zero for both electromagnetic gauges, 3. For PaveTracker readings, statistically one third of the lots do not have an intercept that is zero and two thirds of the lots do not rule out a scaler correction factor of zero, 4. For PQI readings, statistically the 95% confidence interval rules out the intercept being zero for all seven projects and six of the seven projects do not rule out the scaler correction factor being zero, 5. The PQI 301 gauge should not be used for quality control or quality assurance, and 6. The Pavetracker 2701B gauge can be used for quality control but not quality assurance. This study has found that with the limited sample size, the adjusted density equations for both electromagnetic gauges were determined to be inadequate. The PaveTracker Plus 2701B was determined to be better than the PQI 301. The PaveTracker 2701B could still be applicable for quality assurance if the number of core locations per day is reduced and supplemented with additional PaveTracker 2701B readings. Further research should be done to determine the minimum number of core locations to calibrate the gauges each day/lot and the number of additional PaveTracker 2701B readings required.

Crack Development in Ternary Mix Concrete Utilizing Various Saw Depths, TR-587, 2009

Early entry sawing applies sawing earlier and more shallowly than conventional sawing and is believed to increase sawing productivity and reduce the cost of the joint sawing operations. However, some early entry sawing joints (transverse joints) in Iowa were found to experience delayed cracking, sometimes up to 30 days. A concern is whether early entry sawing can lead to late-age random cracking. The present study investigated the effects of different sawing methods on random cracking in portland cement concrete (PCC) pavements. The approach was to assess the cracking potential at sawing joints by measuring the strain development of the concrete at the joints using concrete embedment strain gages. Ten joints were made with the early entry sawing method to a depth of 1.5 in., and two strain gages were installed in each of the joints. Another ten joints were made with the conventional sawing method, five of which were sawed to a depth of one-third of the pavement thickness (3.3 in.), and the other five of which were sawed to a depth of one-quarter of the pavement thickness (2.5 in.). One strain gage was installed in each joint made using conventional sawing. In total, 30 strain gages were installed in 20 joints. The results from the present study indicate that all 30 joints cracked within 25 days after paving, though most joints made using early entry sawing cracked later than the joints made using conventional sawing. No random cracking was observed in the early entry sawing test sections two months after construction. Additionally, it was found that the strain gages used were capable of monitoring the deformations at the joints. The joint crack times (or crack initiation time) measured by the strain gages were generally consistent with the visual observations.

Fly Ash Soil Stabilization for Non-Uniform Subgrade Soils, Volume II: Influence of Subgrade Non-Uniformity on PCC Pavement Performance, TR-461, 2005

To provide insight into subgrade non-uniformity and its effects on pavement performance, this study investigated the influence of non-uniform subgrade support on pavement responses (stress and deflection) that affect pavement performance. Several reconstructed PCC pavement projects in Iowa were studied to document and evaluate the influence of subgrade/subbase non-uniformity on pavement performance. In situ field tests were performed at 12 sites to determine the subgrade/subbase engineering properties and develop a database of engineering parameter values for statistical and numerical analysis. Results of stiffness, moisture and density, strength, and soil classification were used to determine the spatial variability of a given property. Natural subgrade soils, fly ash-stabilized subgrade, reclaimed hydrated fly ash subbase, and granular subbase were studied. The influence of the spatial variability of subgrade/subbase on pavement performance was then evaluated by modeling the elastic properties of the pavement and subgrade using the ISLAB2000 finite element analysis program. A major conclusion from this study is that non-uniform subgrade/subbase stiffness increases localized deflections and causes principal stress concentrations in the pavement, which can lead to fatigue cracking and other types of pavement distresses. Field data show that hydrated fly ash, self-cementing fly ash-stabilized subgrade, and granular subbases exhibit lower variability than natural subgrade soils. Pavement life should be increased through the use of more uniform subgrade support. Subgrade/subbase construction in the future should consider uniformity as a key to long-term pavement performance.

Higher glandular trichome density in tomato leaflets and repellence to spider mites

The objective of this work was to evaluate the feasibility of selection for higher glandular trichome densities, as an indirect criterion of selection for increasing repellence to spider mites Tetranychus urticae, in tomato populations derived from an interspecific cross between Lycopersicon esculentum x L. hirsutum var. glabratum PI 134417. Trichome densities were evaluated in 19 genotypes, including 12 from advanced backcross populations, derived from the original cross L. esculentum x L. hirsutum var. glabratum PI 134417. Counts were made both on the adaxial and abaxial leaf surfaces, and trichomes were classified into glandular types IV and VI, other glandular types (types I+VII), and nonglandular types. Mite repellence was measured by distances walked by mites onto the tomato leaf surface after 20, 40 and 60 min. Spider mite repellence biotests indicated that higher densities of glandular trichomes (especially type VI) decreased the distances walked by the mites onto the tomato leaf surface. Selection of plants with higher densities of glandular trichomes can be an efficient criterion to obtain tomato genotypes with higher resistance (repellence) to spider mites.

Bone mineral density in breast cancer patients treated with adjuvant letrozole, tamoxifen, or sequences of letrozole and tamoxifen in the BIG 1-98 study (SAKK 21/07).

BACKGROUND: The risk of osteoporosis and fracture influences the selection of adjuvant endocrine therapy. We analyzed bone mineral density (BMD) in Swiss patients of the Breast International Group (BIG) 1-98 trial [treatment arms: A, tamoxifen (T) for 5 years; B, letrozole (L) for 5 years; C, 2 years of T followed by 3 years of L; D, 2 years of L followed by 3 years of T]. PATIENTS AND METHODS: Dual-energy X-ray absorptiometry (DXA) results were retrospectively collected. Patients without DXA served as control group. Repeated measures models using covariance structures allowing for different times between DXA were used to estimate changes in BMD. Prospectively defined covariates were considered as fixed effects in the multivariable models. RESULTS: Two hundred and sixty-one of 546 patients had one or more DXA with 577 lumbar and 550 hip measurements. Weight, height, prior hormone replacement therapy, and hysterectomy were positively correlated with BMD; the correlation was negative for letrozole arms (B/C/D versus A), known osteoporosis, time on trial, age, chemotherapy, and smoking. Treatment did not influence the occurrence of osteoporosis (T score < -2.5 standard deviation). CONCLUSIONS: All aromatase inhibitor regimens reduced BMD. The sequential schedules were as detrimental for bone density as L monotherapy.

Effect of once-yearly zoledronic acid five milligrams on fracture risk and change in femoral neck bone mineral density.

CONTEXT: In the Health Outcomes and Reduced Incidence with Zoledronic Acid Once Yearly - Pivotal Fracture Trial (HORIZON-PFT), zoledronic acid (ZOL) 5 mg significantly reduced fracture risk. OBJECTIVE: The aim of the study was to identify factors associated with greater efficacy during ZOL 5 mg treatment. DESIGN, SETTING, AND PATIENTS: We conducted a subgroup analysis (preplanned and post hoc) of a multicenter, double-blind, placebo-controlled, 36-month trial in 7765 women with postmenopausal osteoporosis. Intervention: A single infusion of ZOL 5 mg or placebo was administered at baseline, 12, and 24 months. MAIN OUTCOME MEASURES: Primary endpoints were new vertebral fracture and hip fracture. Secondary endpoints were nonvertebral fracture and change in femoral neck bone mineral density (BMD). Baseline risk factor subgroups were age, BMD T-score and vertebral fracture status, total hip BMD, race, weight, geographical region, smoking, height loss, history of falls, physical activity, prior bisphosphonates, creatinine clearance, body mass index, and concomitant osteoporosis medications. RESULTS: Greater ZOL induced effects on vertebral fracture risk were seen with younger age (treatment-by-subgroup interaction, P = 0.05), normal creatinine clearance (P = 0.04), and body mass index >or= 25 kg/m(2) (P = 0.02). There were no significant treatment-factor interactions for hip or nonvertebral fracture or for change in BMD. CONCLUSIONS: ZOL appeared more effective in preventing vertebral fracture in younger women, overweight/obese women, and women with normal renal function. ZOL had similar effects irrespective of fracture risk factors or femoral neck BMD.

Models for vertical wood density of Scots pine, Norway spruce and birch stems, and their application to determine average wood density

Abstract

Reflective Crack Mitigation Guide for Flexible Pavements, TR-641, 2015

Reflective cracks form in pavements when hot-mix asphalt (HMA) overlays are placed over jointed and/or severely cracked rigid and flexible pavements. In the first part of the research, survival analysis was conducted to identify the most appropriate rehabilitation method for composite pavements and to evaluate the influence of different factors on reflective crack development. Four rehabilitation methods, including mill and fill, overlay, heater scarification (SCR), and rubblization, were analyzed using three performance indicators: reflective cracking, international roughness index (IRI), and pavement condition index (PCI). It was found that rubblization can significantly retard reflective cracking development compared to the other three methods. No significant difference for PCI was seen among the four rehabilitation methods. Heater scarification showed the lowest survival probability for both reflective cracking and IRI, while an overlay resulted in the poorest overall pavement condition based on PCI. In addition, traffic level was found not to be a significant factor for reflective cracking development. An increase in overlay thickness can significantly delay the propagation of reflective cracking for all four treatments. Soil types in rubblization pavement sites were assessed, and no close relationship was found between rubblized pavement performance and subgrade soil condition. In the second part of the research, the study objective was to evaluate the modulus and performance of four reflective cracking treatments: full rubblization, modified rubblization, crack and seat, and rock interlayer. A total of 16 pavement sites were tested by the surface wave method (SWM), and in the first four sites both falling weight deflectometer (FWD) and SWM were conducted for a preliminary analysis. The SWM gave close concrete layer moduli compared to the FWD moduli on a conventional composite pavement. However, the SWM provided higher moduli for the rubblized concrete layer. After the preliminary analysis, another 12 pavement sites were tested by the SWM. The results showed that the crack and seat method provided the highest moduli, followed by the modified rubblization method. The full rubblization and the rock interlayer methods gave similar, but lower, moduli. Pavement performance surveys were also conducted during the field study. In general, none of the pavement sites had rutting problems. The conventional composite pavement site had the largest amount of reflective cracking. A moderate amount of reflective cracking was observed for the two pavement sites with full rubblization. Pavements with the rock interlayer and modified rubblization treatments had much less reflective cracking. It is recommended that use of the modified rubblization and rock interlayer treatments for reflective cracking mitigation are best.

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.