Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “pore pressure”

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 19 records

Effect of pore pressure on the velocity of compressional waves in low-porosity rocks.

The velocity V sub p of compressional waves has been measured in rock samples of low porosity to confining pressures P sub c of 2 kb for a number of different constant pore pressures P sub p. An effective pressure defined by P sub e = P sub c-nP sub p, n less than or equal to 1, is found to be the determining factor in the behavior of V sub p rather than an effective pressure defined simply by the differential pressure Delta P = P sub c-P sub p. As pore pressure increases at constant effective pressure, the value of n increases and approaches 1, but as effective pressure increases at constant pore pressure, the value of n decreases. These observations are consistent with Biot's theory of the propagation of elastic waves in a fluid-saturated porous solid.

Todd, T.↗

Explicit Pore Pressure Material Model in Carbon-Cloth Phenolic

An explicit material model that uses predicted pressure in the pores of a carbon-cloth phenolic (CCP) composite has been developed. This model is intended to be used within a finite-element model to predict phenomena specific to CCP components of solid-fuel-rocket nozzles subjected to high operating temperatures and to mechanical stresses that can be great enough to cause structural failures. Phenomena that can be predicted with the help of this model include failures of specimens in restrained-thermal-growth (RTG) tests, pocketing erosion, and ply lifting

Gutierrez-Lemini, Danton↗

Ice in Channels and Ice-Rock Mixtures in Valleys on Mars: Did They Slide on Deformable Rubble Like Antarctic Ice Streams?

Recent studies of ice streams in Antarctica reveal a mechanism of basal motion that may apply to channels and valleys on Mars. The mechanism is sliding of the ice on deformable water-saturated till under high pore pressures. It has been suggested by Lucchitta that ice was present in outflow channels on Mars and gave them their distinctive morphology. This ice may have slid like Antarctic ice streams but on rubbly weathering products rather than till. However, to generate water under high pore pressures, elevated heatflow is needed to melt the base of the ice. Either volcanism or higher heatflow more than 2 b.y. ago could have raised the basal temperature. Regarding valley networks, higher heatflow 3 b.y. ago could have allowed sliding of ice-saturated overburden at a few hundred meters depth. If the original, pristine valleys were somewhat deeper than they are now, they could have formed by the same mechanism. Recent sounding of the seafloor in front of the Ross Ice Shelf in Antarctica reveals large persistent patterns of longitudinal megaflutes and drumlinoid forms, which bear remarkable resemblance to longitudinal grooves and highly elongated streamlined islands found on the floors of martian outflow channels. The flutes are interpreted to have formed at the base of ice streams during the last glacial advance. Additional similarities of Antarctic ice streams with martian outflow channels are apparent. Antarctic ice streams are 30 to 80 km wide and hundreds of kilometers long. Martian outflow channels have similar dimensions. Ice stream beds are below sea level. Carr determined that most common floor elevations of martian outflow channels lie below martian datum, which may have been close to or below past martian sea levels. The Antarctic ice stream bed gradient is flat and locally may go uphill, and surface slopes are exceptionally. Martian channels also have floor gradients that are shallow or go uphill locally and have low surface gradients. The depth to the bed in ice streams is 1 to 1.5 km. At bankful stage, the depth of the fluid in outflow channels was 1 to 2 km, according to the height of bordering scarps. The similarity between Antarctic ice streams and martian outflow channels suggests that ice may have flowed through and shaped the outflow channels, and that perhaps the mechanism of motion of Antarctic ice streams also operated in outflow channels. In addition, sliding on deformable rubble may explain the formation of small valley networks. The large Siple Coast Antarctic ice streams are thought to slide over longitudinally grooved, deforming till, where much of the movement is within the till. The till is saturated with water at high pore pressures that nearly supports all of the weight of the ice. The small differential between overburden pressure and pore pressure at the bed is more important than the volume of water, but water needs to be supplied to the till interface. For pore pressures to remain high, the ice streams have to act as a seal that blocks the flow of water through them, and the rock underneath has to be of low permeability to prevent the water from draining away.

Lucchitta, B. K.↗

Thermostructural Analysis of Carbon Cloth Phenolic Material Tested at the Laser Hardened Material Evaluation Laboratory

RSRM nozzle liner components have been analyzed and tested to explore the occurrence of anomalous material performance known as pocketing erosion. Primary physical factors that contribute to pocketing seem to include the geometric permeability, which governs pore pressure magnitudes and hence load, and carbon fiber high temperature tensile strength, which defines a material limiting capability. The study reports on the results of a coupled thermostructural finite element analysis of Carbon Cloth Phenolic (CCP) material tested at the Laser Hardened Material Evaluation Laboratory (the LHMEL facility). Modeled test configurations will be limited to the special case of where temperature gradients are oriented perpendicular to the composite material ply angle. Analyses were conducted using a transient, one-dimensional flow/thermal finite element code that models pore pressure and temperature distributions and in an explicitly coupled formulation, passes this information to a 2-dimensional finite element structural model for determination of the stress/deformation behavior of the orthotropic fiber/matrix CCP. Pore pressures are generated by thermal decomposition of the phenolic resin which evolve as a multi-component gas phase which is partially trapped in the porous microstructure of the composite. The nature of resultant pressures are described by using the Darcy relationships which have been modified to permit a multi-specie mass and momentum balance including water vapor condensation. Solution to the conjugate flow/thermal equations were performed using the SINDA code. Of particular importance to this problem was the implementation of a char and deformation state dependent (geometric) permeability as describing a first order interaction between the flow/thermal and structural models. Material property models are used to characterize the solid phase mechanical stiffness and failure. Structural calculations were performed using the ABAQUS code. Iterations were made between the two codes involving the dependent variables temperature, pressure and across-ply strain level. Model results comparisons are made for three different surface heat rates and dependent variable sensitivities discussed for the various cases.

Clayton, J. Louie↗

Extensional reactivation of an abandoned thrust - A bound on shallowing in the brittle regime

Shallow dip angles (not greater than 45 deg) suggested by field observations of continental extensional faults are not predicted by classical isotropic Mohr-Coulomb-Anderson theory. Earthquake data indicate that normal faults exist in the upper crust with dip angles commonly as shallow as 30 deg. One explanation of the apparent conflict between theory and structural/seismic observations is that intrinsically weak, shallow-dipping pre-existing faults are preferentially reactivated. An analytical treatment clarifies the roles of geometry, intact/preexisting fault strengths and fluid pore pressures. Frictional strength ratios of 3 or greater could account for extremely shallow normal faults (dips 10-20 deg) without consideration of pore pressures in excess of the least principal stress or of principal stress systems rotated away from the gravity vector. Moderate reduction in friction (3/4) with respect to wall rock can reduce the dip to 30 deg and can account for shallow normal-slip earthquakes.

Ivins, E. R.↗

The effects of pressure, temperature, and pore water on velocities in Westerly granite

A description is presented of an experimental assembly which has been developed to conduct concurrent measurements of compressional and shear wave velocities in rocks at high temperatures and confining pressures and with independent control of the pore pressure. The apparatus was used in studies of the joint effects of temperature, external confining pressure, and internal pore water on sonic velocities in Westerly granite. It was found that at a given temperature, confining pressure has a larger accelerating effect on compressional waves in dry rock, whereas at a given confining pressure, temperature has a larger retarding effect on shear waves.

Spencer, J. W., Jr.↗

Elasticity of water-saturated rocks as a function of temperature and pressure.

Compressional and shear wave velocities of water-saturated rocks were measured as a function of both pressure and temperature near the melting point of ice to confining pressure of 2 kb. The pore pressure was kept at about 1 bar before the water froze. The presence of a liquid phase (rather than ice) in microcracks of about 0.3% porosity affected the compressional wave velocity by about 5% and the shear wave velocity by about 10%. The calculated effective bulk modulus of the rocks changes rapidly over a narrow range of temperature near the melting point of ice, but the effective shear modulus changes gradually over a wider range of temperature. This phenomenon, termed elastic anomaly, is attributed to the existence of liquid on the boundary between rock and ice due to local stresses and anomalous melting of ice under pressure.

Takeuchi, S.↗

Thermo-Poro-Mechanical Modeling of RTV Intumescence

Room temperature vulcanizing (RTV) silicone is a high-temperature adhesive used as a gap-filler between heatshield tiles in numerous entry missions. Its propensity to intumesce, or swell upon exposure to heat, is a well-known effect that needs to be carefully quantified during design. At tile interfaces of charring ablators, intumescence, combined with differential recession, could cause the gap filler to protrude past the ablator outer mold line, forming a “fence”. Fencing can in turn cause transition to turbulence of the flow wetting the heat shield, leading to augmented surface heating. Recent experiments conducted at the Plasmatron X facility, the high enthalpy wind tunnel of the Center for Hypersonics and Entry Systems Studies, have shown prominent fencing of RTV gap fillers in PICA, under both nitrogen and air plasmas. Similar observations are well known in the arcjet literature. Further experiments under controlled environment, performed using in situ X-ray micro-computed tomography (micro-CT) at the Advanced Light Source (ALS), have shown heating rate-dependent swelling and shrinkage of RTV during pyrolysis. To simulate RTV intumescence, a novel model was introduced in the Porous Materials Analysis Toolbox based on OpenFOAM, PATO, to account for pore-pressure buildup within both closed- and open-pores. The governing equation for the thermo-poro-mechanical response were developed, assuming linear elasticity for the charring silicone. A new multi-pyrolysis model that tracks non-monotonic advancement of material properties with pyrolysis was proposed. This model addresses the limitations of state-of-the-art ablator models to capture the different stages of thermal degradation and coupled thermomechanics. Swelling of RTV was simulated using the new thermo-poro-mechanical model and compared against in situ micro-CT data. Results showed good agreement in intumescence height and temperature profiles at all heating rates, indicating that the key factor contributing to RTV swelling is the internal pressure build-up within closed- and open-pores. As RTV is cured into a soft (rubbery) compound with low-porosity and permeability, initial temperature increase and pyrolysis gas production cause a significant increase of internal pressure, causing a pronounced volume growth. As thermal degradation progresses, rigidization of the silicone occurs due to char hardening which counteract volume shrinkage after gas pressure relief. Overall, our model shows that accounting for changes in properties such porosity, permeability and key thermomechanical coefficients is crucial for capturing the RTV volume change during ablation and enable a predictive capability for heatshield tile interface response. A plan for future calibration of thermomechanical properties evolution during degradation is discussed, as a key next step to close the new model.

RTV↗

Thermo-Poro-Mechanical Modeling of RTV Intumescence

Room temperature vulcanizing (RTV) silicone is a high-temperature adhesive used as a gap-filler between heatshield tiles in numerous entry missions. Its propensity to intumesce, or swell upon exposure to heat, is a well-known effect that needs to be carefully quantified during design. At tile interfaces of charring ablators, intumescence, combined with differential recession, could cause the gap filler to protrude past the ablator outer mold line, forming a “fence”. Fencing can in turn cause transition to turbulence of the flow wetting the heat shield, leading to augmented surface heating. Recent experiments conducted at the Plasmatron X facility, the high enthalpy wind tunnel of the Center for Hypersonics and Entry Systems Studies, have shown prominent fencing of RTV gap fillers in PICA, under both nitrogen and air plasmas. Similar observations are well known in the arcjet literature. Further experiments under controlled environment, performed using in situ X-ray micro-computed tomography (micro-CT) at the Advanced Light Source (ALS), have shown heating rate-dependent swelling and shrinkage of RTV during pyrolysis. To simulate RTV intumescence, a novel model was introduced in the Porous Materials Analysis Toolbox based on OpenFOAM, PATO, to account for pore-pressure buildup within both closed- and open-pores. The governing equation for the thermo-poro-mechanical response were developed, assuming linear elasticity for the charring silicone. A new multi-pyrolysis model that tracks non-monotonic advancement of material properties with pyrolysis was proposed. This model addresses the limitations of state-of-the-art ablator models to capture the different stages of thermal degradation and coupled thermomechanics. Swelling of RTV was simulated using the new thermo-poro-mechanical model and compared against in situ micro-CT data. Results showed good agreement in intumescence height and temperature profiles at all heating rates, indicating that the key factor contributing to RTV swelling is the internal pressure build-up within closed- and open-pores. As RTV is cured into a soft (rubbery) compound with low-porosity and permeability, initial temperature increase and pyrolysis gas production cause a significant increase of internal pressure, causing a pronounced volume growth. As thermal degradation progresses, rigidization of the silicone occurs due to char hardening which counteract volume shrinkage after gas pressure relief. Overall, our model shows that accounting for changes in properties such porosity, permeability and key thermomechanical coefficients is crucial for capturing the RTV volume change during ablation and enable a predictive capability for heatshield tile interface response. A plan for future calibration of thermomechanical properties evolution during degradation is discussed, as a key next step to close the new model.

silicone intumescence↗

Finite Element Analysis of Poroelastic Composites Undergoing Thermal and Gas Diffusion

A theory for time-dependent thermal and gas diffusion in mechanically time-rate-independent anisotropic poroelastic composites has been developed. This theory advances previous work by the latter two authors by providing for critical transverse shear through a three-dimensional axisymmetric formulation and using it in a new hypothesis for determining the Biot fluid pressure-solid stress coupling factor. The derived governing equations couple material deformation with temperature and internal pore pressure and more strongly couple gas diffusion and heat transfer than the previous theory. Hence the theory accounts for the interactions between conductive heat transfer in the porous body and convective heat carried by the mass flux through the pores. The Bubnov Galerkin finite element method is applied to the governing equations to transform them into a semidiscrete finite element system. A numerical procedure is developed to solve the coupled equations in the space and time domains. The method is used to simulate two high temperature tests involving thermal-chemical decomposition of carbon-phenolic composites. In comparison with measured data, the results are accurate. Moreover unlike previous work, for a single set of poroelastic parameters, they are consistent with two measurements in a restrained thermal growth test.

Salamon, N. J.↗

Hot Isostatic Pressing for Reduction of Defects in Additively Manufactured Inconel-718 Preliminary Simulation Approach and Results

Selective laser melting (SLM) of Inconel-718 (IN718) is a promising advanced manufacturing technology; however, challenges exist for fracture-critical parts where pores induced during the manufacturing process reduce fatigue life of the material. Hot isostatic pressing (HIP) is a common post processing technique used to reduce the part porosity by exposing parts to high temperature and pressures. Pore closure depends on the pore shape, trapped gas, and HIP parameters, and the process can only be optimized by understanding the effect of these variables. An analysis capability is described to better understand the conditions necessary for pore closure, and the effect of several variables on porosity are explored and compared to a baseline case. The importance of valid material response models is identified as a significant factor in model accuracy; especially the creep material response model which, in the baseline case, predicts that creep strain accounts for 87.7% of the total strain at the conclusion of the HIP cycle. Lessons learned are described, in addition to future work including model validation.

Joshua M Fody↗

The seismic response of an aquifer to the propagation of an impact generated shock wave: A possible trigger of the Martian outflow channels?

Aquifer dilation from shock waves produced by the 8.4 magnitude Alaskan earthquake of 1964 led to water and sediment ejection from the ground up to 400 km away from the earthquake's epicenter. Groundwater disturbances were observed as far away as Perry, Florida (approximately 5500 km), where well water fluctuations with an amplitude of as much as 2.3 m were reported. The martian cratering record provides evidence that the planet has experienced numerous seismic events of a similar, and often much greater, magnitude. Given this fact, and the photogeologic evidence for abundant water in the early crust, the response of a basalt aquifer to the propagation of compressional waves (P-waves) produced by impacts in the 33-1000 km diameter size range were investigated. The resulting one-dimensional changes in effective stress and pore pressure were calculated - as a function of both distance and time - based on the following assumptions: (1) that all of the seismic energy radiated by an impact is transmitted as a single compressional wave; (2) that both the host rock and groundwater are compressible; and (3) that there is no net flow between the water-filled pores.

Leyva, Ivett A.↗

Formation of Martian flood features by release of water from confined aquifers

It is proposed that the rapid release of water under great pressure from deeply buried aquifers is responsible for the formation of the Martian channels suggestive of catastrophic flooding (outflow channels). Fine channels in the Martian surface suggest the presence of surface water early in the history of the planet, which would have entered the ground water system through the porous near-surface rocks. Subsequent global cooling would have trapped the ground water under a thick permafrost layer and formed a system of confined aquifers. High pore pressures within the aquifers are considered to have triggered the breakout of water from the aquifers at rates of from 10 to the 5th to 10 to the 7th cu m/sec, which would be prevented from reentering the ground water system by the layer of permafrost. Outflow from the aquifer is also considered to have caused the undermining of adjacent areas and the collapse of the surface to form areas of chaos, often associated with channels.

Carr, M. H.↗

Porous Media Approach for Modeling Closed Cell Foam

In order to minimize boil off of the liquid oxygen and liquid hydrogen and to prevent the formation of ice on its exterior surface, the Space Shuttle External Tank (ET) is insulated using various low-density, closed-cell polymeric foams. Improved analysis methods for these foam materials are needed to predict the foam structural response and to help identify the foam fracture behavior in order to help minimize foam shedding occurrences. This presentation describes a continuum based approach to modeling the foam thermo-mechanical behavior that accounts for the cellular nature of the material and explicitly addresses the effect of the internal cell gas pressure. A porous media approach is implemented in a finite element frame work to model the mechanical behavior of the closed cell foam. The ABAQUS general purpose finite element program is used to simulate the continuum behavior of the foam. The soil mechanics element is implemented to account for the cell internal pressure and its effect on the stress and strain fields. The pressure variation inside the closed cells is calculated using the ideal gas laws. The soil mechanics element is compatible with an orthotropic materials model to capture the different behavior between the rise and in-plane directions of the foam. The porous media approach is applied to model the foam thermal strain and calculate the foam effective coefficient of thermal expansion. The calculated foam coefficients of thermal expansion were able to simulate the measured thermal strain during heat up from cryogenic temperature to room temperature in vacuum. The porous media approach was applied to an insulated substrate with one inch foam and compared to a simple elastic solution without pore pressure. The porous media approach is also applied to model the foam mechanical behavior during subscale laboratory experiments. In this test, a foam layer sprayed on a metal substrate is subjected to a temperature variation while the metal substrate is stretched to simulate the structural response of the tank during operation. The thermal expansion mismatch between the foam and the metal substrate and the thermal gradient in the foam layer causes high tensile stresses near the metal/foam interface that can lead to delamination.

Ghosn, Louis J.↗

The measurement of long period and secular deformation with deep borehole tiltmeters

Two clusters of instruments were emplaced in fractured bedrock in eastern Massachusetts in 1970 and 1975. The intrasite agreement at tidal periods was about two percent, but there is no agreement at longer periods. A strong temperature-induced annual component ranging from 3 to 15 urads was present on instruments installed at depths of 15-20m; it was not apparent on those at 100-120m. One instrument, in continuous operation for three years at 100m, showed a net drift of 0.3 urads down to the SW, with a maximum departure of 2.0 urads from the trend. Pore pressure variations, material corrosion and creep, and local movements are apparently the limiting factors to long-term measurements.

Cabaniss, G. H.↗

Gravity-driven deformation of the crust on Venus

Intense horizontal deformation features are observed on the surface of Venus despite the apparent absence of terrestrial-style plate tectonics. High surface temperature, low erosion rate, and regions of elevated terrain suggest that gravity spreading may be responsible for some of these features. Gravity spreading on the Earth is known to cause extensive deformation in certain thin-skinned tectonic regions, where the uppermost layers become detached along a decollement. While the high pore pressures and soft sedimentary rocks frequently observed along decollements on the Earth are not plausible for Venus, temperature induced regions of ductile deformation may exist. Simple viscous models of gravity spreading above a ductile region of the crust are developed to examine the effects of rheology, thermal gradient, surface slope and crustal thickness on deformation. For a reasonable range of parameters, the models show that surface-slope-driven ductile deformation at relatively shallow depths in the crust is possible and may be capable of detaching the uppermost rigid portion of the crust, resulting in horizontal translation.

Smrekar, Suzanne↗

Mars low albedo regions: Possible map of near-surface

A freeze/thaw desorption mechanism is proposed in certain low albedo areas which could be the factor that instigated dust storms. It is widely accepted that the bulk of the episodic gas evolution (not necessarily the oxygen release) experienced during the humidification process in the Viking Gas Exhange Experiment (GEX) was caused by a familiar process in which more polar H2O molecules replace large quantities of other preadsorbed gas molecules on adsorption sites. The author suggests that a similar process could produce high pore pressures in soil that could disrupt the soil and eject dust at high velocity. The author also argued that association of sites of dust storms initiated with high thermal inertial areas may simply reflect repeated dust depletion.

Huguenin, R. L.↗