skip to main content

Explore Research Products in the NSF-PAR

Title: Measuring Basal Force Fluctuations of Debris Flows Using Seismic Recordings and Empirical Green's Functions
Abstract

We present a novel method for measuring the fluctuating basal normal and shear stresses of debris flows by using along‐channel seismic recordings. Our method couples a simple parameterization of a debris flow as a seismic source with direct measurements of seismic path effects using empirical Green's functions generated with a force hammer. We test this method using two large‐scale (8 and 10 m3) experimental flows at the U.S. Geological Survey debris‐flow flume that were recorded by dozens of three‐component seismic sensors. The seismically derived basal stress fluctuations compare well in amplitude and timing to independent force plate measurements within the valid frequency range (15–50 Hz). We show that although the high‐frequency seismic signals provide band‐limited forcing information, there are systematic relations between the fluctuating stresses and independently measured flow properties, especially mean basal shear stress and flow thickness. However, none of the relationships are simple, and since the flow properties also correlate with one another, we cannot isolate a single factor that relates in a simple way to the fluctuating forces. Nevertheless, our observations, most notably the gradually declining ratio of fluctuating to mean basal stresses during flow passage and the distinctive behavior of the coarse, unsaturated flow front, imply that flow style may be a primary control on the conversion of translational to vibrational kinetic energy. This conversion ultimately controls the radiation of high‐frequency seismic waves. Thus, flow style may provide the key to revealing the nature of the relationship between fluctuating forces and other flow properties.

  more » « less
NSF-PAR ID:
10449065
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Earth Surface
Volume:
125
Issue:
9
ISSN:
2169-9003
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ( , Journal of Non-Newtonian Fluid Mechanics)
    Poole, Robert J (Ed.)
    The FENE-P (Finitely-Extensible Nonlinear Elastic) dumbbell constitutive equation is widely used in simulations and stability analyses of free and wall-bounded viscoelastic shear flows due to its relative simplicity and accuracy in predicting macroscopic properties of dilute polymer solutions. The model contains three independent material parameters, which expressed in dimensionless form correspond to a Weissenberg number (Wi), i.e., the ratio of the dumbbell relaxation time scale to a characteristic flow time scale, a finite extensibility parameter (L), corresponding to the ratio of the fully extended dumbbell length to the root mean square end-to-end separation of the polymer chain under equilibrium conditions, and a solvent viscosity ratio. An exact solution for the rheological predictions of the FENE-P model in steady simple shear flow is available [Sureshkumar et al., Phys Fluids (1997)], but the resulting nonlinear and nested set of equations do not readily reveal the key shear-thinning physics that dominates at high Wi as a result of the finite extensibility of the polymer chain. In this note we review a simple way of evaluating the steady material functions characterizing the nonlinear evolution of the polymeric contributions to the shear stress and first normal stress difference as the shear rate increases, provide asymptotic expansions as a function of Wi , and show that it is in fact possible to construct universal master curves for these two material functions as well as the corresponding stress ratio. Steady shear flow experiments on three highly elastic dilute polymer solutions of different finite extensibilities also follow the identified master curves. The governing dimensionless parameter for these master curves is Wi/L and it is only in strong shear flows exceeding Wi/L > 1 that the effects of finite extensibility of the polymer chains dominate the evolution of polymeric stresses in the flow field. We suggest that reporting the magnitude of Wi/L when performing stability analyses or simulating shear-dominated flows with the FENE-P model will help clarify finite extensibility effects. 
    more » « less
  2. ( , Dissertation)
    The rising global trend to reduce dependence on fossil fuels has provided significant motivation toward the development of alternative energy conversion methods and new technologies to improve their efficiency. Recently, oscillating energy harvesters have shown promise as highly efficient and scalable turbines, which can be implemented in areas where traditional energy extraction and conversion are either unfeasible or cost prohibitive. Although such devices are quickly gaining popularity, there remain a number of hurdles in the understanding of their underlying fluid dynamics phenomena. The ability to achieve high efficiency power output from oscillating airfoil energy harvesters requires exploitation of the complexities of the event of dynamic stall. During dynamic stall, the oncoming flow separates at the leading edge of the airfoil to form leading ledge vortex (LEV) structures. While it is well known that LEVs play a significant role in aerodynamic force generation in unsteady animal flight (e.g. insects and birds), there is still a need to further understand their spatiotemporal evolution in order to design more effective energy harvesting enhancement mechanisms. In this work, we conduct extensive experimental investigations to shed-light on the flow physics of a heaving and pitching airfoil energy harvester operating at reduced frequencies of k = fc=U1 = 0.06-0.18, pitching amplitude of 0 = 75 and heaving amplitude of h0 = 0:6c. The experimental work involves the use of two-component particle image velocimetry (PIV) measurements conducted in a wind tunnel facility at Oregon State University. Velocity fields obtained from the PIV measurements are analyzed qualitatively and quantitatively to provide a description of the dynamics of LEVs and other flow structures that may be present during dynamic stall. Due to the difficulties of accurately measuring aerodynamic forces in highly unsteady flows in wind tunnels, a reduced-order model based on the vortex-impulse theory is proposed for estimating the aerodynamic loadings and power output using flow field data. The reduced-order model is shown to be dominated by two terms that have a clear physical interpretation: (i) the time rate of change of the impulse of vortical structures and (ii) the Kutta-Joukowski force which indirectly represents the history effect of vortex shedding in the far wake. Furthermore, the effects of a bio-inspired flow control mechanism based on deforming airfoil surfaces on the flow dynamics and energy harvesting performance are investigated. The results show that the aerodynamic loadings, and hence power output, are highly dependent on the formation, growth rate, trajectory and detachment of the LEV. It is shown that the energy harvesting efficiency increases with increasing reduced frequency, peaking at 25% when k = 0.14, agreeing very well with published numerical results. At this optimal reduced frequency, the time scales of the LEV evolution and airfoil kinematics are matched, resulting in highly correlated aerodynamic load generation and airfoil motion. When operating at k > 0:14, it is shown that the aerodynamic moment and airfoil pitching motion become negatively correlated and as a result, the energy harvesting performance is deteriorated. Furthermore, by using a deforming airfoil surface at the leading and trailing edges, the peak energy harvesting efficiency is shown to increase by approximately 17% and 25% relative to the rigid airfoil, respectively. The performance enhancement is associated with enhanced aerodynamic forces for both the deforming leading and trailing edges. In addition, The deforming trailing edge airfoil is shown to enhance the correlation between the aerodynamic moment and pitching motion at higher reduced frequencies, resulting in a peak efficiency at k = 0:18 as opposed to k = 0:14 for the rigid airfoil. 
    more » « less
  3. ; ; ( , Soil Dynamics and Earthquake Engineering)
    We present a new method to obtain dynamic body force at virtual interfaces to reconstruct shear wave motions induced by a source outside a truncated computational domain. Specifically, a partial differential equation (PDE)-constrained optimization method is used to minimize the misfit between measured motions at a limited number of sensors on the ground surface and their counterparts reconstructed from optimized forces. Numerical results show that the optimized forces accurately reconstruct the targeted ground motions in the surface and the interior of the domain. The proposed optimization framework yields a particular force vector among other valid solutions allowed by the domain reduction method (DRM). Per this optimized or inverted force vector, the reconstructed wave field is identical to its reference counterpart in the domain of interest but may differ in the exterior domain from the reference one. However, we remark that the inverted solution is valid and introduce a simple post-process that can modify the solution to achieve an alternative force vector corresponding to the reference wave field. We also study the desired sensor spacing to accurately reconstruct the wave responses for a given dominant frequency of interest. We remark that the presented method is omnidirectionally applicable in terms of the incident angle of an incoming wave and is effective for any given material heterogeneity and geometry of layering of a reduced domain. The presented inversion method requires information on the wave speeds and dimensions of only a reduced domain. Namely, it does not need any informa- tion on the geophysical profile of an enlarged domain or a seismic source profile outside a reduced domain. Thus, the computational cost of the method is compact even though it leads to the high-fidelity reconstruction of wave re- sponse in the reduced domain, allowing for studying and predicting ground and structural responses using real seismic measurements. 
    more » « less
  4. ; ; ; ; ; ; ; ; ( , Proceedings of the National Academy of Sciences)
    Debris flows are dense and fast-moving complex suspensions of soil and water that threaten lives and infrastructure. Assessing the hazard potential of debris flows requires predicting yield and flow behavior. Reported measurements of rheology for debris flow slurries are highly variable and sometimes contradictory due to heterogeneity in particle composition and volume fraction ( ϕ ) and also inconsistent measurement methods. Here we examine the composition and flow behavior of source materials that formed the postwildfire debris flows in Montecito, CA, in 2018, for a wide range of ϕ that encapsulates debris flow formation by overland flow. We find that shear viscosity and yield stress are controlled by the distance from jamming, Δ ϕ = ϕ m − ϕ , where the jamming fraction ϕ m is a material parameter that depends on grain size polydispersity and friction. By rescaling shear and viscous stresses to account for these effects, the data collapse onto a simple nondimensional flow curve indicative of a Bingham plastic (viscoplastic) fluid. Given the highly nonlinear dependence of rheology on Δ ϕ , our findings suggest that determining the jamming fraction for natural materials will significantly improve flow models for geophysical suspensions such as hyperconcentrated flows and debris flows. 
    more » « less
  5. ( , GSA Bulletin)
    Abstract Many low-angle normal faults (dip ≤30°) accommodate tens of kilometers of crustal extension, but their mechanics remain contentious. Most models for low-angle normal fault slip assume vertical maximum principal stress σ1, leading many authors to conclude that low-angle normal faults are poorly oriented in the stress field (≥60° from σ1) and weak (low friction). In contrast, models for low-angle normal fault formation in isotropic rocks typically assume Coulomb failure and require inclined σ1 (no misorientation). Here, a data-based, mechanical-tectonic model is presented for formation of the Whipple detachment fault, southeastern California. The model honors local and regional geologic and tectonic history and laboratory friction measurements. The Whipple detachment fault formed progressively in the brittle-plastic transition by linking of “minidetachments,” which are small-scale analogs (meters to kilometers in length) in the upper footwall. Minidetachments followed mylonitic anisotropy along planes of maximum shear stress (45° from the maximum principal stress), not Coulomb fractures. They evolved from mylonitic flow to cataclasis and frictional slip at 300–400 °C and ∼9.5 km depth, while fluid pressure fell from lithostatic to hydrostatic levels. Minidetachment friction was presumably high (0.6–0.85), based upon formation of quartzofeldspathic cataclasite and pseudotachylyte. Similar mechanics are inferred for both the minidetachments and the Whipple detachment fault, driven by high differential stress (∼150–160 MPa). A Mohr construction is presented with the fault dip as the main free parameter. Using “Byerlee friction” (0.6–0.85) on the minidetachments and the Whipple detachment fault, and internal friction (1.0–1.7) on newly formed Reidel shears, the initial fault dips are calculated at 16°–26°, with σ1 plunging ∼61°–71° northeast. Linked minidetachments probably were not well aligned, and slip on the evolving Whipple detachment fault probably contributed to fault smoothing, by off-fault fracturing and cataclasis, and to formation of the fault core and fractured damage zone. Stress rotation may have occurred only within the mylonitic shear zone, but asymmetric tectonic forces applied to the brittle crust probably caused gradual rotation of σ1 above it as a result of: (1) the upward force applied to the base of marginal North America by buoyant asthenosphere upwelling into an opening slab-free window and/or (2) basal, top-to-the-NE shear traction due to midcrustal mylonitic flow during tectonic exhumation of the Orocopia Schist. The mechanical-tectonic model probably applies directly to low-angle normal faults of the lower Colorado River extensional corridor, and aspects of the model (e.g., significance of anisotropy, stress rotation) likely apply to formation of other strong low-angle normal faults. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email