skip to main content

Title: Accuracy of slender body theory in approximating force exerted by thin fiber on viscous fluid

We consider the mapping properties of the integral operator arising in nonlocal slender body theory (SBT) for the model geometry of a straight, periodic filament. It is well known that the classical singular SBT integral operator suffers from high wavenumber instabilities, making it unsuitable for approximating theslender body inverse problem, where the fiber velocity is prescribed and the integral operator must be inverted to find the force density along the fiber. Regularizations of the integral operator must therefore be used instead. Here, we consider two regularization methods: spectral truncation and the‐regularization of Tornberg and Shelley (2004). We compare the mapping properties of these approximations to the underlying partial differential equation (PDE) solution, which for the inverse problem is simply the Stokes Dirichlet problem with data constrained to be constant on cross sections. For the straight‐but‐periodic fiber with constant radius, we explicitly calculate the spectrum of the operator mapping fiber velocity to force for both the PDE and the approximations. We prove that the spectrum of the original SBT operator agrees closely with the PDE operator at low wavenumbers but differs at high frequencies, allowing us to define a truncated approximation with a wavenumber cutoff. For both the truncated and‐regularized approximations, we obtain rigorous‐based convergence to the PDE solution as: A fiber velocity withregularity givesconvergence, while a fiber velocity with at leastregularity yieldsconvergence. Moreover, we determine the dependence of the‐regularized error estimate on the regularization parameter.

more » « less
Award ID(s):
Author(s) / Creator(s):
Publisher / Repository:
Date Published:
Journal Name:
Studies in Applied Mathematics
Page Range / eLocation ID:
p. 127-179
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Estimates of turbulence kinetic energy (TKE) dissipation rate (ε) are key in understanding how heat, gas, and other climate‐relevant properties are transferred across the air‐sea interface and mixed within the ocean. A relatively new method involving moored pulse‐coherent acoustic Doppler current profilers (ADCPs) allows for estimates ofεwith concurrent surface flux and wave measurements across an extensive length of time and range of conditions. Here, we present 9 months of moored estimates ofεat a fixed depth of 8.4 m at the Stratus mooring site (20°S, 85°W). We find that turbulence regimes are quantified similarly using the Obukhov length scaleand the newer Langmuir stability length scale, suggesting that ocean‐side friction velocityimplicitly captures the influence of Langmuir turbulence at this site. This is illustrated by a strong correlation between surface Stokes driftandthat is likely facilitated by the steady Southeast trade winds regime. In certain regimes,, whereis the von Kármán constant andis instrument depth, and surface buoyancy flux capture our estimates ofwell, collapsing data points near unity. We find that a newer Langmuir turbulence scaling, based onand, scalesεwell at times but is overall less consistent than. Monin‐Obukhov similarity theory (MOST) relationships from prior studies in a variety of aquatic and atmospheric settings largely agree with our data in conditions where convection and wind‐driven current shear are both significant sources of TKE, but diverge in other regimes.

    more » « less
  2. Abstract

    A detailed chemical kinetic model for oxidation of methylamine has been developed, based on theoretical work and a critical evaluation of data from the literature. The rate coefficients for the reactions of CHNH+ OCHNH/ CHNH + HO, CHNH+ HCH+ NH, CHNHCHNH, and CHNH + OCHNH + HOwere calculated from ab initio theory. The mechanism was validated against experimental results from batch reactors, flow reactors, shock tubes, and premixed flames. The model predicts satisfactorily explosion limits for CHNHand its oxidation in a flow reactor. However, oxidation in the presence of nitric oxide, which strongly promotes reaction at lower temperatures, is only described qualitatively. Furthermore, calculated flame speeds are higher than reported experimental values; the model does not capture the inhibiting effect of the NHgroup in CHNHcompared to CH. More work is desirable to confirm the products of the CHNH + NO reaction and to look into possible pathways to NHin methylamine oxidation.

    more » « less
  3. Abstract

    We present a statistical investigation of the effects of interplanetary magnetic field (IMF) on hemispheric asymmetry in auroral currents. Nearly 6 years of magnetic field measurements from Swarm A and C satellites are analyzed. Bootstrap resampling is used to remove the difference in the number of samples and IMF conditions between the local seasons and the hemispheres. Currents are stronger in Northern Hemisphere (NH) than Southern Hemisphere (SH) for IMF Bin NH (Bin SH) in most local seasons under both signs of IMF B. For Bin NH (Bin SH), the hemispheric difference in currents is small except in local winter when currents in NH are stronger than in SH. During Band Bin NH (Band Bin SH), the largest hemispheric asymmetry occurs in local winter and autumn, when the NH/SH ratio of field aligned current (FAC) is 1.180.09 in winter and 1.170.09 in autumn. During Band Bin NH (Band Bin SH), the largest asymmetry is observed in local autumn with NH/SH ratio of 1.160.07 for FAC. We also find an explicit Beffect on auroral currents in a given hemisphere: on average Bin NH and Bin SH causes larger currents than vice versa. The explicit Beffect on divergence‐free current during IMF Bis in very good agreement with the Beffect on the cross polar cap potential from the Super Dual Auroral Radar Network dynamic model except at SH equinox and NH summer.

    more » « less
  4. Abstract

    The potential for molecular hydrogen () generated via serpentinization to fuel subsurface microbial ecosystems independent from photosynthesis has prompted biogeochemical investigations of serpentinization‐influenced fluids. However, investigations typically sample via surface seeps or open‐borehole pumping, which can mix chemically distinct waters from different depths. Depth‐indiscriminate sampling methods could thus hinder understanding of the spatial controls on nutrient availability for microbial life. To resolve distinct groundwaters in a low‐temperature serpentinizing environment, we deployed packers (tools that seal against borehole walls during pumping) in two‐deep, peridotite‐hosted wells in the Samail Ophiolite, Oman. Isolation and pumping of discrete intervals as deep astobelow ground level revealed multiple aquifers that ranged in pH from 8 to 11. Chemical analyses and 16S rRNA gene sequencing of deep, highly reactedgroundwaters bearing up to,methane () andsulfate () revealed an ecosystem dominated by Bacteria affiliated with the class Thermodesulfovibrionia, a group of chemolithoheterotrophs supported byoxidation coupled toreduction. In shallower, oxidizedgroundwaters, aerobic and denitrifying heterotrophs were relatively more abundant. Highandof(up toand, respectively) indicated microbialoxidation, particularly inwaters with evidence of mixing withwaters. This study demonstrates the power of spatially resolving groundwaters to probe their distinct geochemical conditions and chemosynthetic communities. Such information will help improve predictions of where microbial activity in fractured rock ecosystems might occur, including beyond Earth.

    more » « less
  5. Abstract

    Letbe a graph,be an integer, and writefor the maximum number of edges in an‐vertex graph that is‐partite and has no subgraph isomorphic to. The functionhas been studied by many researchers. Findingis a special case of the Zarankiewicz problem. We prove an analog of the Kövári‐Sós‐Turán theorem for 3‐partite graphs by showingurn:x-wiley:10638539:media:jcd21654:jcd21654-math-0009for. Using Sidon sets constructed by Bose and Chowla, we prove that this upper bound is asymptotically best possible in the case thatandis odd, that is,for. In the cases ofand, we use a result of Allen, Keevash, Sudakov, and Verstraëte, to show that a similar upper bound holds for alland gives a better constant when. Finally, we point out an interesting connection between difference families from design theory and.

    more » « less