SEARCH · Search NASA

Results for “self similarity”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 361 records · Page 20

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Relaxation and darkening of excitonic complexes in electrostatically doped monolayer $\mathrm{WSe_2}$: Roles of exciton-electron and trion-electron interactions

Here we present photoluminescence measurements in monolayer $\mathrm{WSe_2}$, which point to the importance of the interaction between charged particles and excitonic complexes. The theoretical analysis highlights the key role played by exchange scattering, referring to cases wherein the particle composition of the complex changes after the interaction. For example, exchange scattering renders bright excitonic complexes dark in monolayer $\mathrm{WSe_2}$ on accounts of the unique valley-spin configuration in this material. In addition to the ultrafast energy relaxation of hot excitonic complexes following their interaction with electrons or holes, our analysis sheds light on several key features that are commonly seen in the photoluminescence of this monolayer semiconductor. In particular, we can understand why the photoluminescence intensity of the neutral bright exciton is strongest when the monolayer is hole doped rather than charge neutral or electron doped. Similarly, we can understand the reason for the dramatic increase of the photoluminescence intensity of negatively charged excitons (trions) as soon as electrons are added to the monolayer. To self-consistently explain the findings, we further study the photoluminescence spectra at different excitation energies and analyze the behavior of the elusive indirect exciton.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Density functional based methods for light harvesting molecules

The goal of the proposal was to computationally study the photoinduced charge transfer excited states of organic molecules. The systems of interest are the large organic supramolecules that have potential for photovoltaic applications. The simulation of the charge transfer excited states in large supramolecular assemblies require quantum mechanical methods that can be applied to systems with hundreds of atoms. Through this project we have implemented a perturbative delta-self-consistent field method that can easily describe the charge transfer excited states with regular density functionals without range separation. The cost of the calculations of an excited state is similar to that of the ground state. The method also allows for calculation of the energy gradients of the excited state. The developed method was applied to study a number of different supramolecular systems containing hundreds of atoms with excellent agreement with experiment. The absorption characteristics and the influence of density functional approximations were also examined. These led to several studies on the variety of organic electron donor and acceptor molecules. Apart from these developments through this project we also improved the code capabilities by implementing new meta-GGA functionals. This was a major step since the code is purely density functional and earlier the functional capabilities were limited only to generalized gradient approximations. A total of 23 peer-reviewed publications in top chemical physics journals resulted from this grant.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Charge Regulation Stabilizes the Formation of Ionic Liquid‐Based Amphiphilic Oligomer Droplet Interface Bilayers

Amphiphilic charged oligomers (oligodimethylsiloxane – methylimidazolium cation, ODMS-MIM (+) ), assemble into bilayers using the droplet interface bilayer (DIB) platform, possess similar size and functionality as phospholipid bilayers, but exhibit increased stability. The oligomer ionic headgroups (MIM (+) ) are covalently bound to monodisperse, short-chain (n = 13) hydrophobic tails (ODMS). These self-assemble as monolayer brushes at the oil–aqueous interface of water droplets that are influenced by both the charged cationic headgroups, and the nature of the covalently attached tails in the organic phase. Charge regulation (CR) stabilizes the formation of ordered, molecularly close-packed brush phases, which results in highly insulating, stable DIB membranes, with contributions from specific ion-pairing effects, Debye screening, and voltage-dependent electrocompressive stresses. In the oil phase, interactions between hexadecane, a good solvent for ODMS, and the hydrophobic tails result in extended waiting times for bilayer formation compared to phospholipid DIBs, for which hexadecane is a poor solvent. Close agreement between experimental values and predictions for two key parameters, the critical membrane thickness, h c , and maximal grafted headgroup density, Γ 0 , validate an electrostatic CR model consisting of adsorption and partial neutralization of counterions at a charged interface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Two-dimensional simulations of internal gravity waves in the radiation zones of intermediate-mass stars

Intermediate-mass main-sequence stars have large radiative envelopes overlying convective cores. This configuration allows internal gravity waves (IGWs) generated at the convective–radiative interface to propagate towards the stellar surface. The signatures of these waves can be observed in the photometric and spectroscopic data from stars. Here, we have studied the propagation of these IGWs using two-dimensional (2D) fully non-linear hydrodynamical simulations with realistic stellar reference states from the 1D stellar evolution code, Modules for Stellar Astrophysics (MESA). When a single wave is forced, we observe wave self-interaction. When two waves are forced, we observe non-linear interaction (i.e. triadic interaction) between these waves forming waves at different wavelengths and frequencies. When a spectrum of waves similar to that found in numerical simulations is forced, we find that the surface IGW frequency slope is consistent with recent observations. This power law is similar to that predicted by linear theory for the wave propagation, with small deviations that can be an effect of non-linearities. When the same generation spectrum is applied to 3 M ⊙ models at different stellar rotation and ages, the surface IGW spectrum slope is very similar to the generation spectrum slope.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Spontaneous Cyclogenesis without Radiative and Surface-Flux Feedbacks

Tropical cyclones (TCs) are among the most intense and feared storms in the world. What physical processes lead to cyclogenesis remains the most mysterious aspect of TC physics. Here, we study spontaneous TC genesis in rotating radiative–convective equilibrium using cloud-resolving simulations over an f plane with constant sea surface temperature. Previous studies proposed that spontaneous TC genesis requires either radiative or surface-flux feedbacks. To test this hypothesis, we perform mechanism-denial experiments, in which we switch off both feedback processes in numerical simulations. We find that TCs can self-emerge even without radiative and surface-flux feedbacks. Although these feedbacks accelerate the genesis and impact the size of the TCs, TCs in the experiments without them can reach similar intensities as those in the control experiment. We show that TC genesis is associated with an increase in the available potential energy (APE) and that convective heating dominates the APE production. Overall, our result suggests that spontaneous TC genesis may result from a cooperative interaction between convection and circulation and that radiative and surface-flux feedbacks accelerate the process. Furthermore, we find that increasing the planetary rotation favors spontaneous TC genesis.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

QCD–Gravity Double Copy in Regge Asymptotics: From $2\rightarrow n$ Amplitudes to Radiation in Shockwave Collisions

This paper discusses multi-particle production in QCD and in gravity at ultrarelativistic energies, their double-copy relations, and strong parallels in emergent shockwave dynamics. Dispersive techniques are applied to derive the BFKL equation for multi-gluon production in Regge asymptotics. Identical methods apply in gravity and are captured by a gravitational Lipatov equation. The building blocks in both cases are Lipatov vertices and reggeized propagators satisfying double-copy relations; in gravity, Weinberg’s soft theorem is recovered as a limit of the Lipatov framework. BFKL evolution in QCD generates wee parton states of maximal occupancy characterized by an emergent semi-hard saturation scale. Renormalization group equations in the Color Glass Condensate (CGC) EFT describe wee parton correlations and their rapidity evolution. A shockwave picture of deeply inelastic scattering and hadron–hadron collisions follows, with multi-particle production described by Cutkosky’s rules in strong time-dependent fields. Gluon radiation in the CGC EFT has a double copy in gravitational shockwave collisions, with a similar correspondence applicable between gluon and graviton shockwave propagators. Possible extensions of this semi-classical double copy are outlined for computing multi-particle production in gravitational shockwave collisions, self-force and tidal contributions, and classical and quantum noise in the focusing of geodesics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Structural and chemical changes from CO2 exposure to self-healing polymer cement composites for geothermal wellbores

Wellbore cement is subjected to a number of mechanical, thermal and chemical stress regimes over its lifetime. Therefore, next generation wellbore cement formulations need to be evaluated in conditions relevant to these environments. In this work, we investigate the mechanism of the alteration of a novel self-healing polymer-cement composite recently reported by our group after exposure to a CO2-rich environment by using synchrotron based X-ray Fluorescence (XRF) and X-ray absorption near edge structure (XANES) and scanning electron microscopy coupled with energy dispersive spectroscopy. Results showed that chemical alteration of the polymer-cement follows the rim carbonation mechanism, similar to conventional cement although carbonation takes place to a lesser extent in polymer-cements despite the higher porosity. Along with detailed mechanistic insights on carbonation in polymer-cement composite, the performance of these in CO2-rich environment is further studied using standard compressive strength analysis.

Elbakhshwan, Mohamed↗

DOE OSTI2024

Anion Selectivities in Zwitterion Grafted Nanopores: Effect of Zwitterion Architecture

The separation of ions of similar charge is a crucial challenge in many applications, from water treatment to precious metal recovery. Membranes with cross-linked zwitterionic amphiphilic copolymer (ZAC-X) selective layers, which feature self-assembled, zwitterion-lined nanodomains for permeation, offer unique permselectivity between monovalent anions (e.g. $\mathrm{Cl^-}$/$\mathrm{F^-}$). This has motivated studies on the mechanisms of transport and selectivity in this family of materials. In this study, we conducted molecular dynamics simulations of aqueous salt solutions within zwitterion-functionalized nanopores to elucidate the influence of dipole orientation of the ZI ligands on anion diffusivities, partitioning, and permeabilities. Here, our model compares systems with contrasting ZI organization: surface-cation-anion ($\mathrm{S\text{-}ZI^{+}\text{-}ZI^{-}}$, Motif A) and surface-anion-cation ($\mathrm{S\text{-}ZI^{-}\text{-}ZI^{+}}$, Motif B). Our results reveal that Motif A exhibits less pronounced ion pairing due to a spatial separation in the radial profiles of cations and anions. Motif B demonstrates prominent ion pairing for smaller anions owing to their overlap with cation distributions. Further, our potential of mean force profiles reveal that anion partitioning increases with anion size in both ligand motifs, whereas Motif B exhibits significantly higher partitioning selectivity towards larger anions compared to Motif A. Our results for ion diffusivities show that the self-diffusivities of both anions and cations are lower for Motif B compared to Motif A. Such trends in anion partitioning and diffusivities can be explained by differences in the interactions and steric hindrance experienced by the anionic species in Motifs A and B. Finally, our results for anion permselectivity, obtained by combining partitioning and diffusivity, indicate that partitioning trends dominate over diffusivity trends. Consequently, anion permeability increases with anion size, and ligand Motif B yields much higher permselectivity towards larger anions compared to ligand Motif A.

anions↗

DOE OSTI2022

Positrons for Our Future

JLab Users want beams of polarized (and unpolarized) positrons, with beam qualities and modes of operation similar to those provided in the electron beams at CEBAF today [Gra19]. The creation of positrons with a high degree of spin polarization is however very challenging – limited to beta decay, self-polarization in rings, or via e+/e- pair-production, a process requiring a source of polarized gamma rays. At JLab we have a unique source of polarized gamma rays – the bremsstrahlung radiation produced when the spin polarized electron beams of CEBAF interacts with matter – called Polarized Electrons for Polarized Positrons (or PEPPo) [Abb04]. This project sought to extend the PEPPo concept by exploring a possible layout and parameter space for a polarized positron injector at CEBAF based on this approach.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2026

Hierarchical Chiral Self-Assembly of Nanocylinders Composed of Sequence-Defined Mesogenic Dimers

Chiral ensembles can arise through supramolecular curvature that resolves geometric frustrations in the packing of bent, achiral molecular or colloidal building blocks. Here, we leverage orthogonal protection−deprotection click chemistry to create sequence-defined mesogenic heterodimers exhibiting emergent chirality. We compare the hierarchical self-assembly of the synthesized asymmetric, achiral heterodimers, which differ only in the position of a methyl substituent. Both dimers form chiral spherulites composed of nanocylinders. However, the detailed arrangement of nanocylinders depends on the position of the methyl substituent and the crystallization conditions. Despite the chemical similarity, in one dimer, two crystalline forms are optically active. They form conglomerates of dextrorotatory and levorotatory spherulites. The other dimer forms more highly anisotropic spherulites that mask circular birefringence arising from the misorientation of nanocylinders, while mapping of nanocylinder directors reveals a sense at the spherulite surface. We propose that differences in nanocylinder arrangements may arise from changes in nanocylinder curvature and dimensions dictated by the methyl substituent position, inducing chirality. These results demonstrate multiscale hierarchical assembly relevant to dense systems of tubular structures and highlight the role of sequence and molecular design in directing the bottom-up hierarchical self-assembly and chirality of mesogenic systems.

Alkyls↗

DOE OSTI2025

Langevin Dynamics modeling of gas-phase ion-ion recombination (Final Technical Report)

A self-consistent trajectory simulation approach to model MN reactions (Fig. 1) which incorporates the probability of electron transfer as a Monte Carlo operator (Fig. 2) was developed and published as Liu et al. J. Chem. Phys. 159, 114111 (2023). The electron transfer probability p ET estimated using the two-state Landau-Zener (LZ) theory was incorporated into classical trajectory simulations to elicit predictions of MN reaction cross-section σ (vacuum) or rate constant β (finite pressure). Electronic structure calculations with multireference configuration interaction (MRCI) and large correlation consistent basis sets were used to derive inputs to the LZ theory. The key advance of our trajectory simulation approach is the incorporation of electron transfer probability and the inclusion of the effect of ion-neutral interactions on MN using a Langevin representation of the effect of neutral gas on ions. For H + – H - and Li + – H(D) - pairs, our approach quantitatively agrees with measured speed-dependent cross-sections for up to ~10 5 m/s. For the ion pair Ne + – Cl - , our predictions of the MN rate constant at ~1 torr are a factor of ~2 – 3 higher than the experimentally measured value. Similarly, for Xe + – F - in the pressure range of ~20000 – 80000 Pa, our predictions of the MN rate constant are ~20% lower but are in excellent qualitative agreement with experimental data. The paradigm of using trajectory simulations to self-consistently model MN reactions is the basis for inclusion of additional non-classical, and static magnetic and electric field effects. Subsequent work, published as Roy et al. focused on modeling recombination rate constant for three ion pairs (rare gas Ar + cation and halide anions): Ar + – Cl - , Ar + – Br - , Ar + – I - , 2) considering spin-orbit couplings in the electronic structure calculations to obtain high-fidelity estimates of the electron transfer probability and incorporated within the classical trajectory simulations to elicit predictions. In addition to calculations of ion-ion recombination rate constants, a classical trajectory simulation technique (published as Roy et al. J. Chem. Phys. 162(9), 094104 (2023)) that uses quaternions to represent orientation of non-spherical particles (ions or aerosol particles) was developed to simulate the recombination of diatomic or more generally, polyatomic molecules. Finally, several other ion pairs such as Ne + – Cl - , Kr + – Cl - , were explored using the developed semi-classical trajectory simulations to understand various challenges in tackling electronic structure calculations. Using empirical approaches to parameterize the electron transfer radius, trajectory simulations were also used to probe the effect of ion number density on MN rate constant.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Excitation of fishbone-like mode in tokamaks due to bounce resonances of trapped thermal ions

When the drift kinetic effect of thermal ions is taken into account, a high-frequency fishbone-like mode (FLM) is found to be driven unstable by trapped thermal ions (TTIs) in tokamak plasmas, according to self-consistent magneto-hydrodynamic (MHD)-kinetic hybrid modeling utilizing the MARS-K code [Liu et al., Phys. Plasmas 15, 112503 (2008)] as well as an analytic theory. It is found that, similar to energetic particles, TTIs can also stabilize the internal kink mode, whereas the FLM is excited when the effective beta of TTIs exceeds a threshold value. The real frequency of the FLM is comparable to the bounce frequency of TTIs. The mode structure of the FLM can be significantly different from the conventional step-like function for the associated plasma radial displacement. This drift kinetic induced modification of the mode structure near the q = 1 surface is captured by non-perturbative MHD-kinetic hybrid computations with MARS-K. Furthermore, the FLM can only be triggered by TTIs at sufficiently high thermal temperatures. Both the FLM and the internal kink can be stabilized by sufficiently fast plasma toroidal rotation and parallel sound wave damping. These two conditions of high thermal temperature and (fast) flow stabilization, though making it challenging to observe the TTI-driven FLM in present day experiments, are favorable for the mode excitation in future reactor scale devices.

36 MATERIALS SCIENCE↗

DOE OSTI2022

The Influence of Shear on Deep Convection Initiation. Part II: Simulations

Abstract This study evaluates a hypothesis for the role of vertical wind shear in deep convection initiation (DCI) that was introduced in Part I by examining behavior of a series of numerical simulations. The hypothesis states, “Initial moist updrafts that exceed a width and shear threshold will ‘root’ within a progressively deeper steering current with time, increase their low-level cloud-relative flow and inflow, widen, and subsequently reduce their susceptibility to entrainment-driven dilution, evolving toward a quasi-steady self-sustaining state.” A theoretical model that embodied key elements of the hypothesis was developed in Part I, and the behavior of this model was explored within a multidimensional environmental parameter space. Remarkably similar behavior is evident in the simulations studied here to that of the theoretical model, both in terms of the temporal evolution of DCI and in the sensitivity of DCI to environmental parameters. Notably, both the simulations and theoretical model experience a bifurcation in outcomes, whereby nascent clouds that are narrower than a given initial radius R 0 threshold quickly decay and those above the R 0 threshold undergo DCI. An important assumption in the theoretical model, which states that the cloud-relative flow of the background environment V CR determines cloud radius R , is scrutinized in the simulations. It is shown that storm-induced inflow is small relative to V CR beyond a few kilometers from the updraft edge, and V CR therefore plays a predominant role in transporting conditionally unstable air to the updraft. Thus, the critical role of V CR in determining R is validated.

58 GEOSCIENCES↗

DOE OSTI2023

Neutral gas pressure dependence of ion–ion mutual neutralization rate constants using Landau–Zener theory coupled with trajectory simulations

In this computational study, we describe a self-consistent trajectory simulation approach to capture the effect of neutral gas pressure on ion–ion mutual neutralization (MN) reactions. The electron transfer probability estimated using Landau–Zener (LZ) transition state theory is incorporated into classical trajectory simulations to elicit predictions of MN cross sections in vacuum and rate constants at finite neutral gas pressures. Electronic structure calculations with multireference configuration interaction and large correlation consistent basis sets are used to derive inputs to the LZ theory. The key advance of our trajectory simulation approach is the inclusion of the effect of ion-neutral interactions on MN using a Langevin representation of the effect of background gas on ion transport. For H + – H – and Li + – H(D) – , our approach quantitatively agrees with measured speed-dependent cross sections for up to ~10 5 m/s. For the ion pair Ne + – Cl – , our predictions of the MN rate constant at ~1 Torr are a factor of ~2 to 3 higher than the experimentally measured value. Similarly, for Xe + – F – in the pressure range of ~20 000–80 000 Pa, our predictions of the MN rate constant are ~20% lower but are in excellent qualitative agreement with experimental data. Here, the paradigm of using trajectory simulations to self-consistently capture the effect of gas pressure on MN reactions advanced here provides avenues for the inclusion of additional nonclassical effects in future work.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Vortex bursting and associated twist dynamics on helical vortex tubes and vortex rings

The interaction of opposite-signed twist waves on vortex tubes can lead to vortex bursting, a process where the core expands into a double ring-like structure with strong swirling flows. Previous works have studied vortex bursting on rectilinear vortices by axially perturbing the initial core size to generate the twist waves, and observed largely axisymmetric bursting dynamics. In this work, we numerically study bursting on vortical structures with curved centrelines, analysing the interaction between the centreline dynamics, twist wave generation and propagation, and vortex bursting. We focus on axially perturbed helical vortex tubes with small radius-to-pitch ratios up to$0.0625$, as well as vortex rings with a large radius-to-core size ratio$10$, both at a circulation-based Reynolds number$5000$. The results show that though the initial twist wave propagation speeds are relatively unaffected by the curvature and torsion of the centreline, the bursting process is altered significantly compared with rectilinear vortices. The self-induced rotation of the centreline of the helical tube induces a non-axisymmetric distortion of the bursting structure, which rapidly breaks up the vortex core into small-scale helical structures. A similar destabilization of the bursting structure also occurs on vortex rings. The enstrophy increase and accelerated energy decay associated with bursting are predominantly determined by the twist wave strength, rather than the curvature and torsion of the centreline. Combined, our findings imply that bursting could play an important role in transferring and dissipating energy of vortical structures in wakes, and turbulent flows in general.

Mechanics↗

DOE OSTI2024

Formation of Disordered Cocontinuous Phases by Randomly Linked Star Copolymers

Cocontinuous polymeric nanostructures have garnered significant interest due to their ability to combine different properties of two separate polymer domains. Randomly linked copolymer networks have proven to be especially robust for formation of disordered cocontinuous phases across wide composition ranges (≈30 wt % or more). While theoretical treatments of microphase-separated networks have focused primarily on the role of random elastic forces imposed on the self-assembled nanostructures by virtue of the network architecture, experimental studies seeking to disentangle these contributions from other potential effects, such as dispersity in preferred interfacial curvatures, have been scarce. To provide insight into this matter, we here study the self-assembly of randomly linked star copolymers (RSCs), constructed by linking premade polymer arms of polystyrene (PS) and poly(d,l-lactide) (PLA) using 3, 4, and 6-functional connectors. This architecture yields similar distributions of preferred curvature as networks made using corresponding difunctional strands, but lacks the elastic forces imposed by a network architecture. Gravimetry and small-angle X-ray scattering, coupled with scanning electron microscopy, were performed to identify the percolation of PS/PLA RSCs. Remarkably, the 4-arm RSC system exhibited a disordered cocontinuous window of ≈25 wt %, indicating that dispersity in preferred curvature can in some cases be sufficient to robustly drive formation of this morphology. However, the other RSC architectures showed smaller cocontinuous ranges, which we interpret in terms of the influence of homopolymer stars in the 3-arm case and the narrower distribution of preferred interfacial curvatures in the 6-arm case. Finally, thin layers of interconnected porous PS were achieved by solution-processing, suggesting that RSCs have the potential to serve as a robust and easily processable cocontinuous polymeric nanomaterials in both bulk and membrane geometries.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Observing neutron alignment at high angular momenta in heavy deformed actinides

An outstanding puzzle in the structure of very heavy nuclei lies in the rotation response of deformed transuranic isotopes. While the rotation alignment of pairs of 𝑖 13/2 protons is ubiquitous and well understood, observing the analogous alignment of the 𝑗 15/2 neutrons has been elusive. In this work we present a new quantitative analysis approach: from their moments of inertia as a function of frequency, we compare the alignments of rotational bands with both “low” and high-𝑗 intrinsic structures in the 𝑁 = 151 isotones of Pu, Cm, Cf, and No, thereby distinguishing between the alignments of neutron versus proton pairs through “blocking” effects. In this way, the alignment of 𝑗 15/2 neutrons can be clearly identified for the first time. The protons and neutrons begin aligning at a very similar frequency. The contribution of 𝑗 15/2 neutrons to the total alignment is found to grow with increasing 𝑍 compared to 𝑖 13/2 protons. To explain these observations, it is necessary to adopt a self-consistent set of single-particle energies corresponding to the nuclear shape that minimizes the mass. The shape description needs to include high-order deformations up to at least 𝛽 6 .

A ≥ 220↗

DOE OSTI2021

TEAMER: Electrically Engaged Undulation (EEL) System

The Electrically Engaged UnduLation (EEL) system is a buoyancy-driven submersible device for powering oceanographic instruments. Physically, EEL is a slender body whose flexible spine is made up of energy units interconnected by uniaxial hinges. Each unit consists of a pair of piezoelectric elements that converts the bending stress into electrical current to a battery charging circuit. An outer plastic skin forms a seal against water and allows for flexibility at hinge locations. At the top is a bluff body with electronics that holds a ballast for buoyancy adjustment. The bluff body is also responsible for creating fluid instabilities in its wake. When gliding through the water (mode 2), the spine will flex in response to the alternating vortices that shed from the head. This "lock-in" phenomenon occurs when the frequency at which vortices shed resonates with the EEL natural frequency, during which the efficient gaits were found in species of sea snake, eels, and fish. For active propulsion, a single motor can be placed at the first segment and provide the oscillatory input for propulsion similar to a dolphin's kick. Such anguilliform swimming is both efficient and nearly silent compared to a spinning propeller. Ultimately, mimicking bio-locomotion provides a viable path to a drag-reduced, self-propelled energy harvesting system for ocean monitoring.

16 TIDAL AND WAVE POWER↗