Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “VISCOSITY”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,711 records · Page 95

SPREADING LAYERS IN ACCRETING OBJECTS: ROLE OF ACOUSTIC WAVES FOR ANGULAR MOMENTUM TRANSPORT, MIXING, AND THERMODYNAMICS

Disk accretion at a high rate onto a white dwarf (WD) or a neutron star has been suggested to result in the formation of a spreading layer (SL)—a belt-like structure on the object's surface, in which the accreted matter steadily spreads in the poleward (meridional) direction while spinning down. To assess its basic characteristics, we perform two-dimensional hydrodynamic simulations of supersonic SLs in the relevant morphology with a simple prescription for cooling. We demonstrate that supersonic shear naturally present at the base of the SL inevitably drives sonic instability that gives rise to large-scale acoustic modes governing the evolution of the SL. These modes dominate the transport of momentum and energy, which is intrinsically global and cannot be characterized via some form of local effective viscosity (e.g., α-viscosity). The global nature of the wave-driven transport should have important implications for triggering Type I X-ray bursts in low-mass X-ray binaries. The nonlinear evolution of waves into a system of shocks drives effective rearrangement (sensitively depending on thermodynamical properties of the flow) and deceleration of the SL, which ultimately becomes transonic and susceptible to regular Kelvin–Helmholtz instability. We interpret this evolution in terms of the global structure of the SL and suggest that mixing of the SL material with the underlying stellar fluid should become effective only at intermediate latitudes on the accreting object's surface, where the flow has decelerated appreciably. In the near-equatorial regions the transport is dominated by acoustic waves and mixing is less efficient. We speculate that this latitudinal nonuniformity of mixing in accreting WDs may be linked to the observed bipolar morphology of classical nova ejecta.

79 ASTRONOMY AND ASTROPHYSICS↗

SIMULATIONS OF VISCOUS ACCRETION FLOW AROUND BLACK HOLES IN A TWO-DIMENSIONAL CYLINDRICAL GEOMETRY

We simulate shock-free and shocked viscous accretion flows onto a black hole in a two-dimensional cylindrical geometry, where initial conditions were chosen from analytical solutions. The simulation code used the Lagrangian total variation diminishing plus remap routine, which enabled us to attain high accuracy in capturing shocks and to handle the angular momentum distribution correctly. The inviscid shock-free accretion disk solution produced a thick disk structure, while the viscous shock-free solution attained a Bondi-like structure, but in either case, no jet activity nor any quasi-periodic oscillation (QPO)-like activity developed. The steady-state shocked solution in the inviscid as well as in the viscous regime matched theoretical predictions well. However, increasing viscosity renders the accretion shock unstable. Large-amplitude shock oscillation is accompanied by intermittent, transient inner multiple shocks. This oscillation of the inner part of the disk is interpreted as the source of QPO in hard X-rays observed in micro-quasars. Strong shock oscillation induces strong episodic jet emission. The jets also show the existence of shocks, which are produced as one shell hits the preceding one. The periodicities of the jets and shock oscillation are similar; the jets for the higher viscosity parameter appear to be stronger and faster.

79 ASTRONOMY AND ASTROPHYSICS↗

DEEP CHANDRA OBSERVATIONS OF NGC 1404: CLUSTER PLASMA PHYSICS REVEALED BY AN INFALLING EARLY-TYPE GALAXY

The intracluster medium (ICM), as a magnetized and highly ionized fluid, provides an ideal laboratory to study plasma physics under extreme conditions that cannot be achieved on Earth. NGC 1404 is a bright elliptical galaxy that is being gas stripped as it falls through the ICM of the Fornax Cluster. We use the new Chandra X-ray observations of NGC 1404 to study ICM microphysics. The interstellar medium of NGC 1404 is characterized by a sharp leading edge, 8 kpc from the Galaxy center, and a short downstream gaseous tail. Contact discontinuities are resolved on unprecedented spatial scales (0.″5 = 45 pc) due to the combination of the proximity of NGC 1404, the superb spatial resolution of Chandra , and the very deep (670 ks) exposure. At the leading edge, we observe sub-kiloparsec-scale eddies generated by Kelvin–Helmholtz instability (KHI) and put an upper limit of 5% Spitzer on the isotropic viscosity of the hot cluster plasma. We also observe mixing between the hot cluster gas and the cooler galaxy gas in the downstream stripped tail, which provides further evidence of a low viscosity plasma. The assumed ordered magnetic fields in the ICM ought to be smaller than 5 μ G to allow KHI to develop. The lack of an evident magnetic draping layer just outside the contact edge is consistent with such an upper limit.

79 ASTRONOMY AND ASTROPHYSICS↗

Metal-silicate Partitioning and Its Role in Core Formation and Composition on Super-Earths

We use a thermodynamic framework for silicate-metal partitioning to determine the possible compositions of metallic cores on super-Earths. We compare results using literature values of the partition coefficients of Si and Ni, as well as new partition coefficients calculated using results from laser shock-induced melting of powdered metal-dunite targets at pressures up to 276 GPa, which approaches those found within the deep mantles of super-Earths. We find that larger planets may have little to no light elements in their cores because the Si partition coefficient decreases at high pressures. The planet mass at which this occurs will depend on the metal-silicate equilibration depth. We also extrapolate the equations of state (EOS) of FeO and FeSi alloys to high pressures, and present mass–radius diagrams using self-consistent planet compositions assuming equilibrated mantles and cores. We confirm the results of previous studies that the distribution of elements between mantle and core will not be detectable from mass and radius measurements alone. While observations may be insensitive to interior structure, further modeling is sensitive to compositionally dependent properties, such as mantle viscosity and core freeze-out properties. We therefore emphasize the need for additional high pressure measurements of partitioning as well as EOSs, and highlight the utility of the Sandia Z-facilities for this type of work.

79 ASTRONOMY AND ASTROPHYSICS↗

Turbulent Transport by Diffusive Stratified Shear Flows: From Local to Global Models. I. Numerical Simulations of a Stratified Plane Couette Flow

Shear-induced turbulence could play a significant role in mixing momentum and chemical species in stellar radiation zones, as discussed by Zahn. In this paper we analyze the results of direct numerical simulations of stratified plane Couette flows, in the limit of rapid thermal diffusion, to measure the turbulent viscosity and the turbulent diffusivity of a passive tracer as a function of the local shear and the local stratification. We find that the stability criterion proposed by Zahn, namely that the product of the gradient Richardson number and the Prandtl number must be smaller than a critical values (JPr){sub c} for instability, adequately accounts for the transition to turbulence in the flow, with (JPr){sub c}≃0.007. This result recovers and confirms the prior findings of Prat et al. Zahn’s model for the turbulent diffusivity and viscosity, namely that the mixing coefficient should be proportional to the ratio of the thermal diffusivity to the gradient Richardson number, does not satisfactorily match our numerical data. It fails (as expected) in the limit of large stratification where the Richardson number exceeds the aforementioned threshold for instability, but it also fails in the limit of low stratification where the turbulent eddy scale becomes limited by the computational domain size. We propose a revised model for turbulent mixing by diffusive stratified shear instabilities that properly accounts for both limits, fits our data satisfactorily, and recovers Zahn’s model in the limit of large Reynolds numbers.

79 ASTRONOMY AND ASTROPHYSICS↗

Radial Transport and Meridional Circulation in Accretion Disks

Radial transport of particles, elements and fluid driven by internal stresses in three-dimensional (3D) astrophysical accretion disks is an important phenomenon, potentially relevant for the outward dust transport in protoplanetary disks, origin of the refractory particles in comets, isotopic equilibration in the Earth–Moon system, etc. To gain better insight into these processes, we explore the dependence of meridional circulation in 3D disks with shear viscosity on their thermal stratification, and demonstrate a strong effect of the latter on the radial flow. Previous locally isothermal studies have normally found a pattern of the radial outflow near the midplane, switching to inflow higher up. Here we show, both analytically and numerically, that a flow that is inward at all altitudes is possible in disks with entropy and temperature steeply increasing with height. Such thermodynamic conditions may be typical in the optically thin, viscously heated accretion disks. Disks in which these conditions do not hold should feature radial outflow near the midplane, as long as their internal stress is provided by the shear viscosity. Our results can also be used for designing hydrodynamical disk simulations with a prescribed pattern of the meridional circulation.

79 ASTRONOMY AND ASTROPHYSICS↗

Save the Planet, Feed the Star: How Super-Earths Survive Migration and Drive Disk Accretion

Two longstanding problems in planet formation include (1) understanding how planets survive migration, and (2) articulating the process by which protoplanetary disks disperse—and in particular how they accrete onto their central stars. We can go a long way toward solving both problems if the disk gas surrounding planets has no intrinsic diffusivity (“viscosity”). In inviscid, laminar disks, a planet readily repels gas away from its orbit. On short timescales, zero viscosity gas accumulates inside a planet’s orbit to slow Type I migration by orders of magnitude. On longer timescales, multiple super-Earths (distributed between, say, ∼0.1–10 au) can torque inviscid gas out of interplanetary space, either inward to feed their stars, or outward to be blown away in a wind. We explore this picture with 2D hydrodynamics simulations of Earths and super-Earths embedded in inviscid disks, confirming their slow/stalled migration even under gas-rich conditions, and showing that disk transport rates range up to ∼10{sup −7} M{sub ⊙} yr{sup −1} and scale as M-dot ∝ΣM{sub p}{sup 3/2}, where Σ is the disk surface density and M {sub p} is the planet mass. Gas initially sandwiched between two planets is torqued past both into the inner and outer disks. In sum, sufficiently compact systems of super-Earths can clear their natal disk gas in a dispersal history that may be complicated and non-steady but which conceivably leads over Myr timescales to large gas depletions similar to those characterizing transition disks.

79 ASTRONOMY AND ASTROPHYSICS↗

Dust Density Distribution and Imaging Analysis of Different Ice Lines in Protoplanetary Disks

Recent high angular resolution observations of protoplanetary disks at different wavelengths have revealed several kinds of structures, including multiple bright and dark rings. Embedded planets are the most used explanation for such structures, but there are alternative models capable of shaping the dust in rings as it has been observed. We assume a disk around a Herbig star and investigate the effect that ice lines have on the dust evolution, following the growth, fragmentation, and dynamics of multiple dust size particles, covering from 1 μ m to 2 m sized objects. We use simplified prescriptions of the fragmentation velocity threshold, which is assumed to change radially at the location of one, two, or three ice lines. We assume changes at the radial location of main volatiles, specifically H{sub 2}O, CO{sub 2}, and NH{sub 3}. Radiative transfer calculations are done using the resulting dust density distributions in order to compare with current multiwavelength observations. We find that the structures in the dust density profiles and radial intensities at different wavelengths strongly depend on the disk viscosity. A clear gap of emission can be formed between ice lines and be surrounded by ring-like structures, in particular between the H{sub 2}O and CO{sub 2} (or CO). The gaps are expected to be shallower and narrower at millimeter emission than at near-infrared, opposite to model predictions of particle trapping. In our models, the total gas surface density is not expected to show strong variations, in contrast to other gap-forming scenarios such as embedded giant planets or radial variations of the disk viscosity.

79 ASTRONOMY AND ASTROPHYSICS↗

Physics and Chemistry of Te and HgTe-based Ternary Melts

Historically, the theoretical treatment of the liquid phases has always been more difficult and complicated than that for solid and gas phases. A liquid has no lattice structure as crystalline solids and the atoms/molecules in the liquid can migrate through it relatively rapidly. On the other hand, it is also interacting with many other atoms/molecules so that the simplifications of the kinetic theory of gases cannot be employed. For more complicated liquids, such as the liquids of high ionicity and those containing hydrogen bonds and electric dipoles, the understanding is far from complete. At the same time, accurate information on the physics and chemistry of semiconductor melts is needed for the quantitative descriptions of the process of crystal growth from melt. The pre-crystallization phenomena in the liquid phase are critical because the properties of the grown crystals depend on the state and structure of the melt as well as the thermal history of the melt during solidification process. However, the data on the liquid phase, such as thermophysical properties of semiconductor melts are scarce, especially for the HgTe-based II-VI ternary compound semiconductors because of their high vapor pressure and extreme toxicity. Analysis of the thermophysical properties of the melt can provide information about structural transitions of the melt during the solidification process. From a broader point of view, the structure of liquids is much more complicated than the crystalline solids, especially the relaxation behavior through different thermal histories. The theory of hetero-phase fluctuations of liquids is applicable to any many-body systems including condensed-matter physics, field theory, physics of nuclear-matter, cosmology, biology and even sociology. This book summarizes the physics and chemistry from the experimental measurements and theoretical analyses of phase diagram, thermodynamic properties, density, thermal conductivity, viscosity, and electrical conductivity on the binary, pseudo-binary and ternary melts of the most advanced IR-detector material systems of HgCdTe and HgZnTe as well as the analyses of these results. The main objectives of this study are: (1) to provide the phase diagrams and thermodynamic properties of Hg-Cd-Te and Hg-Zn-Te systems through quantitatively fitting the experimental data by assuming an associated solution model for the liquid phase, (2) to experimentally measure the thermophysical properties of the Hg-Cd-Te and Hg-Zn-Te melts, including density, viscosity, electrical conductivity and thermal conductivity as functions of temperature and composition and (3) to enhance the fundamental knowledge of hetero-phase fluctuations and relaxation phenomena in the melts and extend our understanding of the solidification process in order to interpret the experimental results of crystal growth so as to improve the melt growth processes of the compound semiconductor. The physics and chemistry of Te and HgTe-based ternary melts were explored through the studies of the structural transformation during melting, the supercooling during solidification, the relaxation phenomena after rapid cooling of the melts and the metal-semiconductor transition in the melts through the analyses of electrical conductivity and Lorenz number. An in-depth study on the thermophysical properties and their time-dependent structural dynamic processes taking place in the vicinity of the solid-liquid phase transition of the narrow homogeneity range HgTe-based ternary semiconductors as well as the analysis of the homogenization process in the melt will also be presented.

Thermophysical properties↗

Understanding turbulence in compressing plasmas and its exploitation or prevention.

Unprecedented densities and temperatures are now achieved in compressions of plasma, by lasers and by pulsed power, in major experimental facilities. These compressions, carried out at the largest scale at the National Ignition Facility and at the Z Pulsed Power Facility, have important applications, including fusion, X-ray production, and materials research. Several experimental and simulation results suggest that the plasma in some of these compressions is turbulent. In fact, measurements suggest that in certain laboratory plasma compressions the turbulent energy is a dominant energy component. Similarly, turbulence is dominant in some compressing astrophysical plasmas, such as in molecular clouds. Turbulence need not be dominant to be important; even small quantities could greatly influence experiments that are sensitive to mixing of non-fuel into fuel, such as compressions seeking fusion ignition. Despite its important role in major settings, bulk plasma turbulence under compression is insufficiently understood to answer or even to pose some of the most fundamental questions about it. This thesis both identifies and answers key questions in compressing turbulent motion, while providing a description of the behavior of three-dimensional, isotropic, compressions of homogeneous turbulence with a plasma viscosity. This description includes a simple, but successful, new model for the turbulent energy of plasma undergoing compression. The unique features of compressing turbulence with a plasma viscosity are shown, including the sensitivity of the turbulence to plasma ionization, and a sudden viscous dissipation'' effect which rapidly converts plasma turbulent energy into thermal energy. This thesis then examines turbulence in both laboratory compression experiments and molecular clouds. It importantly shows: the possibility of exploiting turbulence to make fusion or X-ray production more efficient; conditions under which hot-spot turbulence can be prevented; and a lower bound on the growth of turbulence in molecular clouds. This bound raises questions about the level of dissipation in existing molecular cloud models. Finally, the observations originally motivating the thesis, Z-pinch measurements suggesting dominant turbulent energy, are reexamined by self-consistently accounting for the impact of the turbulence on the spectroscopic analysis. This is found to strengthen the evidence that the multiple observations describe a highly turbulent plasma state.

fast ignition↗

Sedimentation behavior of noble metal particles in simulated high-level waste borosilicate glasses

Solubility of noble metal elements (NME) in the melted borosilicate glass is much smaller than its normal concentration of the high level liquid waste. Thus most of NME show small particles in the melted glass and tend to sediment in the bottom region of the vitrification melter due to their higher density than that of glass. Experiments of the sedimentation of NME particles in the melted glass were carried out under static condition. Three conditions of initial NME concentration (1.1, 3.0, 6.1 wt % with an equivalent for each oxide) in the simulated glass were set and held at 1100 C. degrees up to 2880 hours. The specimen with 1.1 wt % initial NME concentration indicated zone settling, and the settling rate of the interface is constant: 2.4 mm/h. This sedimentation behavior is the type of rapid settling. Following the rapid settling, the settling rate goes gradually slower; this is the type of compressive settling. The specimens with 3.0 wt % and 6.1 wt % initial NME concentration showed compression settling from the beginning. From the settling curve of the interface, the maximum concentration of NME in sediment was estimated to be around 23- 26 wt %. Growth of NME particles was observed by holding at 1100 C. degrees for up to 2880 hours. The viscosity becomes higher as NME concentration increases and the dependence on shear rate becomes simultaneously stronger. The effect of the particle growth to viscosity appears to be not significant.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Demonstration of Mixing and Transferring Settling Cohesive Slurry Simulants in the AY-102 Tank - 12323

In support of Hanford's feed delivery of high level waste (HLW) to the Waste Treatment and Immobilization Plant (WTP), pilot-scale testing and demonstrations with simulants containing cohesive particles were performed as a joint collaboration between Savannah River National Laboratory (SRNL) and the Pacific Northwest National Laboratory (PNNL) staff. The objective of the demonstrations was to determine the impact that cohesive particle interactions in the simulants, and the resulting non- Newtonian rheology, have on tank mixing and batch transfer of large and dense seed particles. The work addressed the impacts cohesive simulants have on mixing and batch transfer performance in a pilot-scale system. Kaolin slurries with a range of wt% concentrations to vary the Bingham yield stress were used in all the non-Newtonian simulants. To study the effects of just increasing the liquid viscosity (no yield stress) on mixing and batch transfers, a glycerol/water mixture was used. Stainless steel 100 micron particles were used as seed particles due to their density and their contrasting color to the kaolin and glycerol. Testing results show that water always transfers less seed particles, and is conservative when compared to fluids with a higher yield stress and/or higher viscosity at the same mixing/transfer parameters. The impact of non-Newtonian fluid properties depends on the magnitude of the yield stress. A higher yield stress in the carrier fluid resulted in more seed particles being transferred to the RTs. A dimensional analysis highlighting the role of a yield stress (due to cohesive particle interactions) defined four regions of behavior and indicates how the results obtained in this study can be applied to the full-scale mixing behavior of a high level waste tank. The analysis indicates that the regions of behavior for full-scale mixing have been adequately represented by the current small-scale tests. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Development of Cement Solidification Technique for Sodium Borate Waste Produced in PWR Plants - 16061

A cement solidification process for treating sodium borate waste produced in pressurized water reactor (PWR) plants was studied. It is known that the setting and hardening reaction of cement is retarded by the presence of boron compounds. This technique was developed to obtain high uniaxial compressive strength more than 10 MPa which can satisfy various countries disposal regulations. As examples of disposal regulations, France requires 8 MPa for homogeneous solidified products. The sodium/ boron (Na/B) ratio was controlled to about 1 by adding NaOH to simulated liquid waste for the purpose of improving the workablility and mechanical strength. To obtain high volume reduction, simulated liquid waste was dehydrated and powdered by using a wiped film evaporator. Ordinary portland cement (OPC) was used for the solidification. The viscosity of cement mixture was controlled by the water/cement ratio. The resulting viscosity became less than 25 dPa.s, and it was possible to continuously mix the mixture. Solidified cement with containing 73 kg/m{sup 3} of boron had a uniaxial compressive strength of more than 10 MPa. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A Study on Vitrification for LLW Radioactive Wastes - 16235

Ash and ion exchange resin wastes to be vitrified were considered among various radioactive wastes in the Japanese nuclear industry. The appropriate glass composition for the waste is calculated using the computer program. The glass composition is systematically studied using a ternary diagram, considering viscosity, electrical conductivity, and durability. The optimized glass composition based on calculated parameters of vitrification process and the convenience of mixing wastes is selected for crucible experiment. Since differences between calculated values and measured values may occur upon evaluation of glass formulations using glass formulation computer simulation program, the calculated glass formulations should be verified through experiments on a laboratory scale for more meticulous development of glass formulations. For the follow-up development of glass formulations, several glass formulations are planned to be developed. The candidate glasses for treatment of ash and ion exchange resin were melted and fabricated using reagents, respectively. The chemicals were weighed out to produce the desired batch compositions and mixed thoroughly and then put into clay crucibles for melting. Next, the crucibles were introduced into an electric furnace at 1,150 deg. C in the ambient atmosphere and pressure with occasional stirring by a quartz rod. A hold time of 1 hour was introduced to promote homogenization. Finally, the melted material was poured on a graphite plate to be rapidly cooled. The viscosities of all fabricated glasses were in the desired range. After rapid cooling, the micro-structure of the glasses was analyzed using the SEM. They exhibited amorphous structure with good homogeneity. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A Semi-Empirical Correlation for Counter-Current Flow Limitation in a Vertical Pipe

Counter-current two-phase flows may take place in the steam generator of a pressurized water reactor (PWR) during a reflux-cooling mode. Therefore developing reliable correlations of counter-current flow limitation (CCFL) is of great importance in safety analyses. Wallis derived a CCFL correlation, which is expressed in terms of the so-called Wallis parameters, which are the dimensionless gas and liquid volumetric fluxes. Although effects of the fluid properties are not sufficiently accounted for in this correlation, he introduced the inverse liquid viscosity number into it to account for the viscous effect on CCFL . The functional form of the factor for the viscous effect was however empirically determined without physical grounds. Zapke and Kroeger proposed an empirical CCFL correlation by introducing the Ohnesorge number into the original Wallis correlation to account for the effects of the liquid viscosity and surface tension on CCFL. The effects of the fluid properties have been thus added in the original Wallis correlation in an empirical manner. In this study, a semi-empirical correlation for CCFL in a vertical pipe is derived from one-dimensional equations for a steady counter-current annular two-phase flow in a vertical pipe. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Experimental study on nonmonotonicity of capillary desaturation curves in a 2-D pore-network

Immiscible displacement in a porous medium is important in many applications such as soil remediation and enhanced oil recovery. When gravitational forces are negligible, two-phase immiscible displacement at the pore level is controlled by capillary and viscous forces whose relative importance is quantified through the dimensionless capillary number Ca and the viscosity ratio M between liquid phases. Depending on the values of Ca and M, capillary fingering, viscous fingering, or stable displacement may be observed resulting in a variety of patterns affecting the phase entrapment. The Capillary Desaturation Curve (CDC), which represents the relationship between the residual oils saturation and Ca, is an important relation to describe the phase entrapment at a given Ca. In the present study, we investigate the CDC as influenced by the viscosity ratio. A comprehensive series of experiments using a high-resolution microscope and state-of-the-art micromodels were conducted. The CDCs were calculated and the effects of Ca and M on phase entrapments were quantified. The results show that CDCs are not necessarily monotonic for all M.

capillary desaturation curve↗

Statistical tests of simple earthquake cycle models: Testing Earthquake Cycle Models

A central goal of observing and modeling the earthquake cycle is to forecast when a particular fault may generate an earthquake: a fault late in its earthquake cycle may be more likely to generate an earthquake than a fault early in its earthquake cycle. Models that can explain geodetic observations throughout the entire earthquake cycle may be required to gain a more complete understanding of relevant physics and phenomenology. Previous efforts to develop unified earthquake models for strike-slip faults have largely focused on explaining both preseismic and postseismic geodetic observations available across a few faults in California, Turkey, and Tibet. An alternative approach leverages the global distribution of geodetic and geologic slip rate estimates on strike-slip faults worldwide. Here we use the Kolmogorov-Smirnov test for similarity of distributions to infer, in a statistically rigorous manner, viscoelastic earthquake cycle models that are inconsistent with 15 sets of observations across major strike-slip faults. We reject a large subset of two-layer models incorporating Burgers rheologies at a significance level of α = 0.05 (those with long-term Maxwell viscosities η M <~ 4.0 × 10 19 Pa s and η M >~ 4.6 × 10 20 Pa s) but cannot reject models on the basis of transient Kelvin viscosity η K . Finally, we examine the implications of these results for the predicted earthquake cycle timing of the 15 faults considered and compare these predictions to the geologic and historical record.

58 GEOSCIENCES↗

The Potential for Metamorphic Thermal Pulses to Develop During Compaction-Driven Fluid Flow

Compaction-driven fluid flow below the brittle-ductile transition may be a means of transporting fluids during metamorphism. In particular, when a decompaction weakening mechanism is introduced to account for the rock viscosity reduction due to fluid overpressures, channeling instabilities evolve into high-porosity/permeability fluid conduits that focus mass and energy transfer. In this study, we consider a crustal rheology that accounts simultaneously for upward-increasing viscosity and decompaction weakening to examine the nucleation and evolution of fluid channelization in two dimensions (2-D). The model shows that plume-shaped flow patterns can develop on time scales as short as 104 years, during which the plume tails act as fluid conduits and the plume heads act as fluid dispersion zones near the brittle-ductile transition. Collection of fluids into conduits is accomplished by a basal fluid catchment zone characterized by strong lateral fluid pressure gradients but low porosity/permeability. Relatively narrow ranges of viscous activation energy (~100 kJ mo -1 ) and decompaction weakening factor (~10 -4 ) are constrained if the fluid conduits are of kilometer scale in width. Significant thermal excursions (~65 °C) can be induced if a high flow rate, potentially from rapid intermittent dehydration, is realized within channels. Moreover, if the focused fluids emanate from external anomalously hot sources (e.g., magma intrusion), thermal pulses (>100°C), and steep lateral temperature gradients (>50°C km -1 ) can be generated. Given the focusing efficiency estimated from our 2-D compaction model, simple 3-D modeling further shows that tubular conduits have the potential to cause thermal pulses >200°C within 10 4 years

54 ENVIRONMENTAL SCIENCES↗