Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Constrained control”

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,315 records · Page 73

Modeling the Strahl in the Solar Wind Using CLUSTER/PEACE Electron Distribution Functions

We present a study of kinetic properties of the strahl electron velocity distribution functions (VDF's) in the solar wind. These are used to investigate the pitch-angle scattering and stability of the population to interactions with electromagnetic (whistler) fluctuations. The study is based on high time resolution data from the Cluster/PEACE electron spectrometer. Our study focuses on the mechanisms that control and regulate the pitch-angle and stability of strahl electrons in the solar wind; mechanisms that are not yet well understood. Various parameters are investigated such as the strahl-electron density, temperature anisotropy, and electron heat-flux. The goal is to check whether the strahl electrons are constrained by some instability (e.g., the whistler instability), or are maintained by other types of processes. The electron heat-flux and temperature anisotropy are determined by modeling the 3D-VDFs to a spectral spherical harmonic model from which the moments are obtained directly from the spectral coefficients.

Vinas, Adolfo F.↗

Exploration Systems Development (ESD) Approach to Enterprise Risk Management

The National Aeronautics and Space Administration (NASA) Exploration Systems Development (ESD) Division has implemented an innovative approach to Enterprise Risk Management under a unique governance structure and streamlined integration model. ESD's mission is to design and build the capability to extend human existence to deep space. The Enterprise consists of three Programs: Space Launch System (SLS), Orion, and Ground Systems Development and Operations (GSDO). The SLS is a rocket and launch system that will be capable of powering humans, habitats, and support systems to deep space. Orion will be the first spacecraft in history capable of taking humans to multiple destinations within deep space. GSDO is modernizing Kennedy's spaceport to launch spacecraft built and designed by both NASA and private industry. ESD's approach to Enterprise Risk Management is commensurate with affordability and a streamlined management philosophy. ESD Enterprise Risk Management leverages off of the primary mechanisms for integration within the Enterprise. The Enterprise integration approach emphasizes delegation of authority to manage and execute the majority of cross-program activities and products to the individual Programs, while maintaining the overall responsibility for all cross-program activities at the Division. The intent of the ESD Enterprise Risk Management approach is to improve risk communication, to avoid replication and/or contradictory strategies, and to minimize overhead process burden. This is accomplished by the facilitation and integration of risk information within ESD. The ESD Division risks, Orion risks, SLS risks, and GSDO risks are owned and managed by the applicable Program. When the Programs have shared risks with multiple consequences, they are jointly owned and managed. When a risk is associated with the integrated system that involves more than one Program in condition, consequence, or mitigation plan, it is considered an Exploration Systems Integration (ESI) Risk. An ESI risk may require visibility and risk handling by multiple organizations. The Integrated Risk Working Group (IRWG) is a small team of Risk experts that are responsible for collaborating and communicating best practices. In addition, the forum facilitates proper integration of risks across the Enterprise. The IRWG uses a Continuous Risk Management approach for facilitating the identification, analysis, planning, tracking, and controlling of ESI Risks. The ESD Division, Programs, and Integrated Task Teams identify ESI Risks. The IRWG maintains a set of metrics for understanding Enterprise Risk process and the overall Risk Posture. The team is also actively involved in the modeling of risk for Enterprise Performance Management. With the Enterprise being constrained in Schedule and Budget, and with significant technical complexity, the appropriate use of Risk Management techniques is crucial to the success of the Enterprise. The IRWG achieves this through the modified approach, providing a forum for collaboration on risks that cross boundaries between the separate entities.

Bauder, Stephen P.↗

Cislunar Near Rectilinear Halo Orbit for Human Space Exploration

In order to conduct sustained human exploration beyond Low Earth Orbit (LEO), spacecraft systems are designed to operate in a series of missions of increasing complexity. Regardless of the destination, Moon, Mars, asteroids or beyond, there is a substantial set of common objectives that must be met. Many orbit characterization studies have endeavored to evaluate the potential locations in cislunar space that are favorable for meeting common human exploration objectives in a stepwise approach. Multiple studies, by both NASA and other international space agencies, have indicated that Earth-­‐moon libration point orbits are attractive candidates for staging operations in the proving ground and beyond. In particular, the Near Rectilinear Orbit (NRO) has been demonstrated to meet multi-­‐mission and multi-­‐destination architectural constraints. However, a human mission to a selected NRO presents a variety of new challenges for mission planning. While a growing number of robotic missions have completed successful operations to various specific libration point orbits, human missions have never been conducted to orbits of this class. Human missions have unique challenges that differ significantly from robotic missions, including a lower tolerance for mission risk and additional operational constraints that are associated only with human spacecraft. In addition, neither robotic nor human missions have been operated in the NRO regime specifically, and NROs exhibit dynamical characteristics that can differ significantly as compared to other halo orbits. Finally, multi-­‐body orbits, such as libration point orbits, are identified to exist in a simplified orbit model known as the Circular Restricted Three Body Problem (CRTBP) and must then be re-­‐solved in the full ephemeris model. As a result, the behavior of multi-­‐body orbits cannot be effectively characterized within the classical two-­‐body orbit dynamics framework more familiar to the human spaceflight community. In fact, a given NRO is not identified by a set of Keplerian orbit parameters, and a valid epoch specific state vector must be first obtained from a multi-body dynamical model. In this paper, the significant performance and operational challenges of conducting human missions to the NRO are evaluated. First, a systematic process for generating full ephemeris based ballistic NROs of various families is outlined to demonstrate the relative ease in which a multi-­‐revolution orbit can be found for any epoch and for various orbit geometries. In the Earth-­‐Moon system, NROs, which are halo orbits with close passage over a lunar pole, can exist with respect to libration point 1 (L1) or libration point 2 (L2) and are either from a North or South family orbit class with respect to the ecliptic. Second, the ability to maintain the orbit over the lifetime of a habitat mission by applying a reliable station-keeping strategy is investigated. The NRO, while similar to the quasi-­‐halo orbits that the Artemis mission flew, requires an updated station keeping strategy. This is due to several dynamical differences such as the increased relative stability of the NRO compared to other halo orbits and the close passage over the lunar surface as shown in Figure 1. Multiple station-keeping strategies are being investigated to ensure a human spacecraft remains on a predictable path. As the NRO is not described in simple two-­‐body parameters, analysis must determine the best strategy for targeting a reference NRO as well as how closely a future state should be constrained. In addition, costs will be minimized by determining maneuver directionality based on an identified pattern in the optimal station-keeping solutions or an analytically derived relationship. The candidate station-keeping algorithm must be stable and robust to environmental and vehicle uncertainties as well to navigation estimation and flight control execution errors. To that end, navigation accuracies, the impact on the station-keeping execution errors as well as other vehicle uncertainties need to be assessed. Starting with Orion, current navigation accuracies are evaluated and then navigation requirements are derived assuming a desired station-keeping propellant budget. Third, the performance requirements to and from the NRO are evaluated. Important parameters for developing expected propellant costs include epoch of operation, size and type of NRO, Earth departure and return constraints, as well as abort or early-­‐return capability. Finally, rendezvous and proximity operations are vital aspects of multi-­‐mission human exploration endeavors. The ability to conduct rendezvous and the associated propellant costs are assessed as well as the impacts of various profile assumptions including the location within the NRO the rendezvous is performed. The results of these studies will influence plans for international cooperation on both nearer term proving ground missions and beyond.

Whitley, Ryan↗

NetFlow Dynamics

NetFlow Dynamics is a web-accessible analysis environment for simulating dynamic flows of materials on model networks. Performing a simulation requires both the NetFlow Dynamics application and a network model which is a description of the structure of the nodes and edges of a network including the flow capacity of each edge and the storage capacity of each node, and the sources and sinks of the material flowing on the network. NetFlow Dynamics consists of databases for storing network models, algorithms to calculate flows on networks, and a GIS-based graphical interface for performing simulations and viewing simulation results. Simulated flows are dynamic in the sense that flows on each edge of the network and inventories at each node change with time and can be out of equilibrium with boundary conditions. Any number of network models could be simulated using Net Flow Dynamics. To date, the models simulated have been models of petroleum infrastructure. The main model has been the National Transportation Fuels Model (NTFM), a network of U.S. oil fields, transmission pipelines, rail lines, refineries, tank farms, and distribution terminals. NetFlow Dynamics supports two different flow algorithms, the Gradient Flow algorithm and the Inventory Control algorithm, that were developed specifically for the NetFlow Dynamics application. The intent is to add additional algorithms in the future as needed. The ability to select from multiple algorithms is desirable because a single algorithm never covers all analysis needs. The current algorithms use a demand-driven capacity-constrained formulation which means that the algorithms strive to use all available capacity and stored inventory to meet desired flows to sinks, subject to the capacity constraints of each network component. The current flow algorithms are best suited for problems in which a material flows on a capacity-constrained network representing a supply chain in which the material supplied can be stored at each node of the network. In the petroleum models, the flowing materials are crude oil and refined products that can be stored at tank farms, refineries, or terminals (i.e. the nodes of the network). Examples of other network models that could be simulated are currency flowing in a financial network, agricultural products moving to market, or natural gas flowing on a pipeline network.

Corbet Jr., Thomas F↗

Multichannel Processing of Vibrational Response Measurements: An Application of Constrained Subspace Identification

Dynamic structural systems occur in a wide variety of applications such as autonomous ve-hicles employed militarily for battlefield management or non-military tasks like tracking awildfire in the hills of a drought-ridden area. Many such systems operating in the air or atsea employ a variety of subsystems ranging from airfoils for flight or machinery for propulsionor reactor components as simple as pipes carrying cooling for control.

42 ENGINEERING↗

PDE-constrained high-order mesh optimization

Here, we present a novel framework for PDE-constrained r-adaptivity of high-order meshes. The proposed method formulates mesh movement as an optimization problem, with an objective function defined as a convex combination of a mesh quality metric and a measure of the accuracy of the PDE solution obtained via finite element discretization. The proposed formulation achieves optimized, well-defined high-order meshes by integrating mesh quality control, PDE solution accuracy, and robust gradient regularization. We adopt the Target-Matrix Optimization Paradigm to control geometric properties across the mesh, independent of the PDE of interest. To incorporate the accuracy of the PDE solution, we introduce error measures that control the finite element discretization error. The implicit dependence of these error measures on the mesh nodal positions is accurately captured by adjoint sensitivity analysis. Additionally, a convolution-based gradient regularization strategy is used to ensure stable and effective adaptation of high-order meshes. We demonstrate that the proposed framework can improve mesh quality and reduce the error by up to 10 times for the solution of Poisson and linear elasto-static problems. The approach is general with respect to the dimensionality, the order of the mesh, the types of mesh elements, and can be applied to any PDE that admits well-defined adjoint operators.

Computer science↗

Quantifying Variability and Controls of Riverine Dissolved Organic Carbon Exported to Arctic Coastal Margins of North America (Final Report)

This project involved implementation and application of a coupled permafrost hydrology and dissolved organic carbon process models to investigate how spatial and seasonal variations in terrestrial hydrology and soil freeze/thaw dynamics influence the mobilization, loading, and export of organic carbon to the stream network for selected arctic basins across northern Alaska and northwest Canada. The model simulations were constrained by detailed observations of in- stream chemistry, soil active layer profile moisture and temperature dynamics, streamflow, soil carbon inventories and satellite microwave remote sensing based assessments of surface soil freeze-thaw dynamics. We developed and applied the numerical modeling and data analysis, incorporating observed data for calibration and validation, in order to investigate the terrestrial hydrology, permafrost dynamics, and associated DOC production and loading to rivers across a region encompassing watersheds draining to the coast. The project produced six publications and three datasets archived in public repositories.

54 ENVIRONMENTAL SCIENCES↗

Tidal Control of Jet Eruptions Observed by Cassini ISS

Observations by Cassini's Imaging Science Subsystem (ISS) of Enceladus' south polar region at high phase angles has revealed jets of material venting into space. Observations by Cassini's Composite Infrared Spectrometer (CIRS) have also shown that the south polar region is anomalously warm with hotspots associated with geological features called the Tiger Stripes. The Tiger Stripes are large rifts near the south pole of Enceladus, which are typically about 130 km in length, 2 km wide, with a trough 500 m deep, and are l1anked on each side by 100m tall ridges. Preliminary triangulation of jets as viewed at different times and with different viewing geometries in Cassini ISS images taken between 2005 and 2007 have constrained the locations of eight major eruptions of material and found all of them associated with the south polar fractures unofficially the 'Tiger Stripes', and found four of them coincident with the hotspots reported in 2006 by CIRS. While published ISS observations of jet activity suggest that individual eruption sites stay active on the timescale of years, any shorter temporal variability (on timescales of an orbital period, or 1.3 Earth days, for example) is more difficult to establish because of the spotty temporal coverage and the difficulty of visually isolating one jet from the forest of many seen in a typical image. Consequently, it is not known whether individual jets are continuously active, randomly active, or if they erupt on a predictable, periodic schedule. One mechanism that may control the timing of eruptions is diurnal tidal stress, which oscillates between compression/tension as well as right and left lateral shear at any given location throughout Enceladus' orbit and may allow the cracks to open and close regularly. We examine the stresses on the Tiger Stripe regions to see how well diurnal tidal stress caused by Enceladus' orbital eccentricity may possibly correlate with and thus control the observed eruptions. We then identify possible mechanisms by which tidal stress can provide access to the surface for volatile material and implications for observed jet activity.

Hurford, T. A.↗

Oxide scale exfoliation and regrowth in TP347H superheater tubes: Oxide scale exfoliation and regrowth

This paper provides an introduction to a comprehensive model being developed to predict and control oxide scale exfoliation from the steam-side of superheater and reheater tubes in steam boilers. The model deals with the main phenomena involved in scale growth and failure in steam, and incorporates major variables related to boiler design and operation. The considerations used to calculate oxide growth under the specific constrains of small diameter tubes carrying high-pressure steam and operating with large temperature gradients under temperature and pressure cycling conditions, as well as the evolution of stresses and strains in the scales, are indicated but only a cursory description is given of the details of the analytical treatments. An example is presented of calculations made with the model to predict the extent of blockage expected in a single superheater loop as a function of time and outlet steam temperature under several realistic service conditions. The results suggest that problems due to scale exfoliation would be expected early in the operating life of superheater tubes made from austenitic steel TP347H.

Sabau, A. S.↗

THE SDSS CO-ADD: COSMIC SHEAR MEASUREMENT

Stripe 82 in the Sloan Digital Sky Survey was observed multiple times, allowing deeper images to be constructed by co-adding the data. Here, we analyze the ellipticities of background galaxies in this 275 deg{sup 2} region, searching for evidence of distortions due to cosmic shear. We do so using measurements of both the shear-shear correlation function and power spectrum, with the latter determined using both ''quadratic'' and ''pseudo'' estimation techniques. We show how we verified these methods using mock catalogs. We also describe our methods for modeling and correcting for the effects of the point-spread function (PSF) in our shape measurements, and we also describe our prescription for estimating photometric redshifts (photo-z's) for our galaxy sample. In particular, we assess the impact of potential systematic effects due to the PSF and to photo-z's, and show that these are under control in our analysis. We find consistent correlation function and power spectrum results, where the E-mode cosmic shear signal is detected in both real and Fourier space with >5{sigma} significance on degree scales, while the B-mode is consistent with zero as expected. The amplitude of the signal constrains the combination of the matter density {Omega}{sub m} and fluctuation amplitude {sigma}{sub 8} to be {Omega}{sup 0.7}{sub m}{sigma}{sub 8} = 0.252{sup +0.032}{sub -{sub 0.052}}.

79 ASTRONOMY AND ASTROPHYSICS↗

Spatiotemporally Adaptive Compression for Scientific Dataset with Feature Preservation – A Case Study on Simulation Data with Extreme Climate Events Analysis

Scientific discoveries are increasingly constrained by limited storage space and I/O capacities. For time-series simulations and experiments, their data often need to be decimated over timesteps to accommodate storage and I/O limitations. In this paper, we propose a technique that addresses storage costs while improving post-analysis accuracy through spatiotemporal adaptive, error-controlled lossy compression. We investigate the trade-off between data precision and temporal output rates, revealing that reducing data precision and increasing timestep frequency lead to more accurate analysis outcomes. Additionally, we integrate spatiotemporal feature detection with data compression and demonstrate that performing adaptive error-bounded compression in higher dimensional space enables greater compression ratios, leveraging the error propagation theory of a transformation-based compressor. To evaluate our approach, we conduct experiments using the well-known E3SM climate simulation code and apply our method to compress variables used for cyclone tracking. Our results show a significant reduction in storage size while enhancing the quality of cyclone tracking analysis, both quantitatively and qualitatively, in comparison to the prevalent timestep decimation approach. Compared to three state-of-the-art lossy compressors lacking feature preservation capabilities, our adaptive compression framework improves perfectly matched cases in TC tracking by 26.4-51.3% at medium compression ratios and by 77.3-571.1% at large compression ratios, with a merely 5–11% computational overhead.

Gong, Qian↗

Multirate digital control system design

Methods for multirate digital control system design are discussed. A simple method for sampling rate selection based on control bandwidths is proposed. Methods for generating a discrete-time state model of a sampled-data plant and a discrete-time equivalent to an analog cost function for a sampled-data plant are described. The succesive loop closures and linear quadratic Gaussian synthesis methods are reviewed, and a constrained optimization synthesis method is introduced. The proposed sampling rate selection, discretization, and synthesis methods are applied to two example design problems. Multirate and single-rate compensators synthesized by the different methods are compared, based on closed-loop responses, with compensators having the same real-time computation load.

Berg, Martin C.↗

Options for Robust Airfoil Optimization under Uncertainty

A robust optimization method is developed to overcome point-optimization at the sampled design points. This method combines the best features from several preliminary methods proposed by the authors and their colleagues. The robust airfoil shape optimization is a direct method for drag reduction over a given range of operating conditions and has three advantages: (1) it prevents severe degradation in the off-design performance by using a smart descent direction in each optimization iteration, (2) it uses a large number of spline control points as design variables yet the resulting airfoil shape does not need to be smoothed, and (3) it allows the user to make a tradeoff between the level of optimization and the amount of computing time consumed. For illustration purposes, the robust optimization method is used to solve a lift-constrained drag minimization problem for a two-dimensional (2-D) airfoil in Euler flow with 20 geometric design variables.

Padula, Sharon L.↗

Modeling the Effect of Evaporative Flux and Colloidal Aggregation on Crud Deposition

During the last few years lot of research has been conducted to investigate and predict the Crud deposition in PWRs. One of the common findings of various experimental and computational studies is that thickness or amount of crud deposits is directly related to evaporative flux maps. Evaporation driven deposition profile on the surfaces takes different forms, i.e. from ring to disk shapes during initial stages and in the form of porous to impervious continuous layers covering the hot fuel surface at the later stages. The amount and thickness of these deposits although are directly related to evaporative flux, the relationship is much complex and can be dependent upon several other factors. A popular problem relating evaporative flux to particulate deposition has been well studied in the form of drying of sessile colloidal drops. It was proposed that, this phenomenon is due to a geometrical constraint: the free surface, constrained by a pinned contact line squeezes the fluid outwards to compensate for evaporative losses and thereby transporting particles to the contact line. Later on a quantitative relationship between evaporative flux and deposition rate was proposed and validated by the same group, which showed power-law relationship. Deposition studies conducted with controlled concentration of particles with subcooled nucleate boiling also showed a power-law relationship rather than a simpler linear relationship whereas physics of evaporation would tend to concentrate particles at linear rate. In a typical subcooled nucleate boiling system, bubble growth and departure time is of the order 1-2 ms, measured experimentally. Diffusion time constant for dilute spherical metallic colloids such as crud particulates, generally present as particulate in process piping, range from 0.0025-2.5 sec in water for a size range 0.1-1.0 μm. Colloidal behavioral characteristics during the growth of bubbles can be computed with escape time constant at collision, which is dependent on effective Hamaker constant. For minimum effective Hamaker constant of Hematite is -2.5 kT and diffusion time constant τ = 0.0025 sec, escape time constant is computed as Kramer rule, τexp[A/kT] magnitude is 0.03 sec or 30 ms. This clearly indicates that evaporation will cause particles to come together and they cannot escape at a fast enough rate which would lead them to form agglomerates. Therefore both effects of evaporative flux in conjunction with the agglomeration kinetics must be considered for modeling the crud deposition if the evaporative flux is known. In this paper, we propose a phenomenological model of colloidal transport and aggregation with evaporative flux. The dynamic of particles interactions while drying of sessile colloidal drops is developed by coupling general transport equation with Smoluchowshki Coagulation equation. Crud deposition on the nuclear fuel has been a serious problem with severe consequences in nuclear power plants. In several studies it has been consistently reported that crud deposition is directly related to evaporative flux on the nuclear fuel surfaces. However, it has not been established how to quantify the deposition rates of crud type particulates if the evaporative flux is known. This complex multiphysics problem is similar to coffee-ring effect, where the governing physics is understood more clearly. In this work, a simplified geometry is used to model the effects of evaporative flux acting at a boundary on the colloidal concentration. The effect of this condition on change in concentration is shown numerically. Further, agglomeration is suspected to play a significant role in high evaporative flux such as bubble nucleation. The coagulation kinetics is modeled using uniform particle distribution and the effective total concentration of colloids including original particles and agglomerates is studied. Once this model is validated with a simple separate effect experiments, it can be integrated to large scale PWR models of evaluating evaporative flux and chemistry models developed by Electrical Power Research Institute (EPRI). (authors)

42 ENGINEERING↗

Four-dimensional guidance algorithms for aircraft in an air traffic control environment

Theoretical development and computer implementation of three guidance algorithms are presented. From a small set of input parameters the algorithms generate the ground track, altitude profile, and speed profile required to implement an experimental 4-D guidance system. Given a sequence of waypoints that define a nominal flight path, the first algorithm generates a realistic, flyable ground track consisting of a sequence of straight line segments and circular arcs. Each circular turn is constrained by the minimum turning radius of the aircraft. The ground track and the specified waypoint altitudes are used as inputs to the second algorithm which generates the altitude profile. The altitude profile consists of piecewise constant flight path angle segments, each segment lying within specified upper and lower bounds. The third algorithm generates a feasible speed profile subject to constraints on the rate of change in speed, permissible speed ranges, and effects of wind. Flight path parameters are then combined into a chronological sequence to form the 4-D guidance vectors. These vectors can be used to drive the autopilot/autothrottle of the aircraft so that a 4-D flight path could be tracked completely automatically; or these vectors may be used to drive the flight director and other cockpit displays, thereby enabling the pilot to track a 4-D flight path manually.

Pecsvaradi, T.↗

A dynamic volumetric heat source model for laser additive manufacturing

Melt pool scale models of laser powder bed fusion (LPBF) offer insights into the process-structure-property relationships in additive manufacturing (AM). These models often neglect physical phenomena such as vapor cavity formation and fluid mechanics to reduce computational demands. Instead, volumetric heat source models are used to represent the effects that these phenomena have on the predicted melt pool dimensions. Generally, the dimensions and effective absorption of the volumetric heat source are calibrated to reproduce melt pool dimensions observed in metallographic cross sections taken from single-track experiments on bare plate. However, the transient nature of LPBF often deviates the melt pool dimensions from the assumed steady-state conditions of single-track experiments, motivating the need for a volumetric heat source model that more generally considers the dynamic relationship between melt pool shape and laser-material interactions. Here, we introduce a two-parameter volumetric heat source model that integrates several existing models into a generalized mathematical expression, providing independent control over the radial heat distribution via the parameter k and the volumetric shape of the heat source via the parameter m. This parameterization enables the calibration of melt pool shape predictions through simultaneous adjustment of these parameters, while keeping the radial heat source dimensions consistent with the experimental spot size (D4σ) and constraining the heat source depth and absorption to physically derived expressions for cavities. Consequently, the proposed volumetric heat source model adapts to changes in the local melt pool conditions due to scanning strategy and part geometry by dynamically adjusting the heat source depth and absorption. We demonstrate the capabilities of the proposed model through comparisons with a collection of experiments from the Additive Manufacturing Benchmark (AMBench).

36 MATERIALS SCIENCE↗

Friction-Stability-Permeability Evolution of a Fracture in Granite

The evolution of frictional strength, stability, and fracture permeability is intimately linked both to the seismic cycle and to the impact of hydraulic stimulation for fractured reservoirs. However, despite this importance, the poromechanical relationships between fault permeability and strength remain unclear. The present study explores this relationship via laboratory experiments for concurrent shear-flow on smooth fractures of Westerly granite. The novelty of these experiments is that the shear velocity is precisely controlled during the measurement of fracture permeability. Results indicate permeability enhancement during velocity-weakening (potentially unstable) frictional slip. To decipher key processes contributing to this response, we evaluate the state of contacting asperities and of fracture surface asperities via digital rock fracture modeling of statistically equivalent surfaces. We propose two plausible mechanisms constraining the relationship between friction and permeability evolution—one based on changes in asperity contact distribution and one on shear-induced dilation triggered by changes in fault slip velocity. Finally, these mechanisms should be taken into account in interpreting field observation such as the abrupt permeability increase of natural faults at the onset of seismic slip.

54 ENVIRONMENTAL SCIENCES↗