Semi-flexible filaments interacting with molecular motors and immersed in rheologically complex and viscoelastic media constitute a common motif in biology. Synthetic mimics of filament-motor systems also feature active or field-activated filaments. A feature common to these active assemblies is the spontaneous emergence of stable oscillations as a collective dynamic response. In nature, the frequency of these emergent oscillations is seen to depend strongly on the viscoelastic characteristics of the ambient medium. Motivated by these observations, we study the instabilities and dynamics of a minimal filament-motor system immersed in model viscoelastic fluids. Using a combination of linear stability analysis and full non-linear numerical solutions, we identify steady states, test the linear stability of these states, derive analytical stability boundaries, and investigate emergent oscillatory solutions. We show that the interplay between motor activity, filament and motor elasticity, and fluid viscoelasticity allows for stable oscillations or limit cycles to bifurcate from steady states. When the ambient fluid is Newtonian, frequencies are controlled by motor kinetics at low viscosities, but decay monotonically with viscosity at high viscosities. In viscoelastic fluids that have the same viscosity as the Newtonian fluid, but additionally allow for elastic energy storage, emergent limit cycles are associated with higher frequencies. The increase in frequency depends on the competition between fluid relaxation time-scales and time-scales associated with motor binding and unbinding. Our results suggest that both the stability and oscillatory properties of active systems may be controlled by tailoring the rheological properties and relaxation times of ambient fluidic environments. 
                        more » 
                        « less   
                    
                            
                            Identification of Characteristic Spatial Scales to Improve the Performance of Analytical Spectral Solutions to the Groundwater Flow Equation
                        
                    
    
            Abstract Analytical solutions for the three‐dimensional groundwater flow equation have been widely used to gain insight about subsurface flow structure and as an alternative to computationally expensive numerical models. Of particular interest are solutions that decompose prescribed hydraulic head boundaries (e.g., Dirichlet boundary condition) into a collection of harmonic functions. Previous studies estimate the frequencies and amplitudes of these harmonics with a least‐square approach where the amplitudes are fitted given a pre‐assigned set of frequencies. In these studies, an ad hoc and structured discretization of the frequency domain is typically used, excluding dominant frequencies while assigning importance to spurious frequencies, with significant consequences for estimating the fluxes and residence times. This study demonstrates the advantages of using a pre‐assigned frequency spectrum that targets the dominant frequencies based on rigorous statistical analysis with predefined significance levels. The new approach is tested for three hydrologic conceptualizations: (a) a synthetic periodic basin, (b) synthetic bedforms, and (c) a natural mountainous watershed. The performance of the frequency spectrum selection is compared with exact analytical or approximate numerical solutions. We found that the new approach better describes the fluxes and residence times for Dirichlet boundaries with well‐defined characteristics spatial scales (e.g., periodic basins and bedforms). For more complex scenarios, such as natural mountainous watersheds, both pre‐assigned frequency spectrums present similar performance. The spectral solutions presented here can play a central role in developing reduced‐complexity models for assessing regional water and solute fluxes within mountain watersheds and hyporheic zones. 
        more » 
        « less   
        
    
    
                            - PAR ID:
- 10370856
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Water Resources Research
- Volume:
- 57
- Issue:
- 12
- ISSN:
- 0043-1397
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
- 
            
- 
            Analytical formulas are obtained for frequency-dependent effective elastic modulus and effective mass density for a periodic layered structure. The proposed homogenization procedure is valid at sufficiently high frequencies well above the lowest band gap in the acoustic spectrum of the structure. It is shown that frequency-dependent effective parameters may take negative values either in different regions of frequencies or in the same quite narrow region. This property demonstrates that 1D elastic structure may behave in the limit of small Bloch wave vectors as a double-negative acoustic metamaterial.more » « less
- 
            Abstract A transient two-dimensional acoustic boundary element solver is coupled to a potential flow boundary element solver via Powell's acoustic analogy to determine the acoustic emission of isolated hydrofoils performing biologically-inspired motions. The flow-acoustic boundary element framework is validated against experimental and asymptotic solutions for the noise produced by canonical vortex-body interactions. The numerical framework then characterizes the noise production of an oscillating foil, which is a simple representation of a fish caudal fin. A rigid NACA 0012 hydrofoil is subjected to combined heaving and pitching motions for Strouhal numbers ($0.03 < St < 1$) based on peak-to-peak amplitudes and chord-based reduced frequencies ($0.125 < f^* < 1$) that span the parameter space of many swimming fish species. A dipolar acoustic directivity is found for all motions, frequencies, and amplitudes considered, and the peak noise level increases with both the reduced frequency and the Strouhal number. A combined heaving and pitching motion produces less noise than either a purely pitching or purely heaving foil at a fixed reduced frequency and amplitude of motion. Correlations of the lift and power coefficients with the peak root-mean-square acoustic pressure levels are determined, which could be utilized to develop long-range, quiet swimmers.more » « less
- 
            We study the electrothermal actuation of nanomechanical motion using a combination of numerical simulations and analytical solutions. The nanoelectrothermal actuator structure is a u-shaped gold nanoresistor that is patterned on the anchor of a doubly clamped nanomechanical beam or a microcantilever resonator. This design has been used in recent experiments successfully. In our finite-element analysis (FEA) based model, our input is an ac current; we first calculate the temperature oscillations due to Joule heating using Ohm’s law and the heat equation; we then determine the thermally induced bending moment and the displacement profile of the beam by coupling the temperature field to Euler–Bernoulli beam theory with tension. Our model efficiently combines transient and frequency-domain analyses: we compute the temperature field using a transient approach and then impose this temperature field as a harmonic perturbation for determining the mechanical response in the frequency domain. This unique modeling method offers lower computational complexity and improved accuracy and is faster than a fully transient FEA approach. Our dynamical model computes the temperature and displacement fields in the time domain over a broad range of actuation frequencies and amplitudes. We validate the numerical results by directly comparing them with experimentally measured displacement amplitudes of nano-electro-mechanical system beams around their eigenmodes in vacuum. Our model predicts a thermal time constant of 1.9 ns in vacuum for our particular structures, indicating that electrothermal actuation is efficient up to ∼80 MHz. We also investigate the thermal response of the actuator when immersed in a variety of fluids.more » « less
- 
            Abstract Hydrologic modeling has been a useful approach for analyzing water partitioning in catchment systems. It will play an essential role in studying the responses of watersheds under projected climate changes. Numerous studies have shown it is critical to include subsurface heterogeneity in the hydrologic modeling to correctly simulate various water fluxes and processes in the hydrologic system. In this study, we test the idea of incorporating geophysics‐obtained subsurface critical zone (CZ) structures in the hydrologic modeling of a mountainous headwater catchment. The CZ structure is extracted from a three‐dimensional seismic velocity model developed from a series of two‐dimensional velocity sections inverted from seismic travel time measurements. Comparing different subsurface models shows that geophysics‐informed hydrologic modeling better fits the field observations, including streamflow discharge and soil moisture measurements. The results also show that this new hydrologic modeling approach could quantify many key hydrologic fluxes in the catchment, including streamflow, deep infiltration, and subsurface water storage. Estimations of these fluxes from numerical simulations generally have low uncertainties and are consistent with estimations from other methods. In particular, it is straightforward to calculate many hydraulic fluxes or states that may not be measured directly in the field or separated from field observations. Examples include quickflow/subsurface lateral flow, soil/rock moisture, and deep infiltration. Thus, this study provides a useful approach for studying the hydraulic fluxes and processes in the deep subsurface (e.g., weathered bedrock), which needs to be better represented in many earth system models.more » « less
 An official website of the United States government
An official website of the United States government 
				
			 
					 
					
