SEARCH · Engineering Papers

Results for “Stress strain relations”

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 685 records · Page 38

All sources NASA NTRS DOE OSTI

NASA NTRS1984

Elevated temperature compressive steady state deformation and failure in the oxide dispersion strengthened alloy MA 6000E

The compressive flow strength-strain rate behavior of the oxide-dispersion-strengthened alloy MA 6000E has been studied in the temperature range 1144-1366 K, with strain rates ranging from 2.1 x 10 to the -5th to 2.1 x 10 to the -7th per s. It is found that the inherent strength of the alloy is essentially the same in all test directions and that the low strength observed in tensile tests results from the inability of grain boundaries to support high tensile stresses. The failure of MA 6000E under high-temperature, slow plastic flow conditions is shown to be the result of concentrated slip. Slow plastic deformation in MA 6000E can be described by a threshold stress model of creep where threshold stresses are calculated from relatively fast testing procedures and the effective stress exponent for creep is assumed to be 3.5.

Whittenberger, J. D.↗

DOE OSTI2023

Ion-beam radiation-induced Eshelby transformations: The mean and variance in hydrostatic and shear residual stresses

Ion beam plays a pivotal role in ion implantations and the fabrication of nanostructures. However, there lacks a quantitative model to describe the residual stresses associated with the ion-beam radiation. Radiation-induced residual stress/transformation strain have been mostly recognized in the hydrostatic sub strain space. Here, in this study, we use molecular dynamics (MD) simulations to show that the response of a material to irradiation is generally anisotropic that depends on the ion-beam direction, and should be described using tensorial quantities. We demonstrate that accelerator-based ion beam irradiation, combined with the intrinsic lattice anisotropy and externally induced anisotropy (such as anisotropic mechanical loadings), causes radiation-actuated shear transformation strains in addition to hydrostatic expansion. We map out these complex correlations for several materials. Radiation-induced defects are shown to be responsible for residual shear stresses in the manner of Eshelby inclusion transformation. We propose such tensorial response model should be considered for accurate nanoscale fabrication using ion-beam irradiation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2021

Estimating Internal Stress of an Alteration Layer Formed on Corroded Boroaluminosilicate Glass through Spectroscopic Ellipsometry Analysis

Aqueous corrosion of glass may result in the formation of an alteration layer in the glass surface of which chemical composition and network structure are different from those of the bulk glass. Since corrosion occurs far below the glass-transition temperature, the alteration layer cannot fully relax to the new structure with the lowest possible energy. Molecular dynamics simulations suggested that such a network will contain highly strained chemical bonds, which can be manifested as a stress in the alteration layer. Common techniques to measure stress in thin films or surface layers were found inadequate for thick monolithic glass samples corroded in water. Here, we explored the use of spectroscopic ellipsometry to test the presence of internal stress in the alteration layer formed by aqueous corrosion of glass. Furthermore, a procedure for analyses of spectroscopic ellipsometry data to determine birefringence in the alteration layer was developed. Findings with the established fitting procedure suggested that a stress builds up in the corroded surface layer of a boroaluminosilicate glass if there is a change in relative humidity, pH, or electrolyte concentration of the environment to which the glass surface is exposed. A similar process may occur in other types of glass, and it may affect the surface properties of corroded glass objects.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2026

A stress-based fracture model for reacting metal ejecta

The evolution of reacting metal ejecta continues to be a topic of interest at the forefront of metals in reactive and extreme environments. Ejecta are small particles formed when the surface of a metal undergoes Richtmyer–Meshkov instability from a strong shock. Experiments have shown that in the case where ejecta are in ambient conditions that induce a reaction, the ejecta behave irregularly. The ejecta temperature rises and then plateaus, and the acceleration profile shows unexpected jumps. These variations are assumed to be related to the exothermic heat release and particle mass loss caused by the reaction. To explain this phenomenon, efforts to model this in simulations have increased. While current models can capture many of these physical processes, they currently assign a constant reaction shell thickness with little physical reasoning. This work remedies this problem by assigning a dynamic physically informed shell thickness to the reacting particles, using solid analysis. The shell thickness of the particles impacts the rate of change of reacted mass in the system, as well as the rate at which the particles react. The model is based on a simple stress–strain relationship and gives a dynamic assignment for when the reacting particle should begin to fracture. We compare our model to the previous computational and simulation data to analyze the effects of different model parameters.

42 ENGINEERING↗

NASA NTRS2019

Through-Thickness Tensile Properties of Triaxial Braided Composites

The through-thickness strength of triaxial braided composites was measured using thin specimens bonded to metallic tabs. Despite providing reasonable engineering estimates of the bulk properties, some test specimen limitations were noted: For strong epoxies, significant shear and tensile stresses occur and induce failure between the first layer and bond-line; the strains over much of the exterior are ~10% lower than interior strains, which are representative of the nominal state; and failure often occurred near or along the adhesive bond line. These interference were partially overcome by adding a large radius to the gage section. Based on the finite element analysis and experimental measurements, the nominal elastic modulus and Poisson's ratio are E(sub 33) = 10.7 GPa and v(sub 31) = 0.08 for the three triaxial braided composites tested. Strengths vary from 14 to 40 MPa depending on the resin. Relative to the in-plane properties, the through thickness modulus, fracture strength, failure strain, and Poisson's ratio are 1/5th, 1/30th, 1/6th, and 1/4th. Interlaminar failure occurred with varying degrees of fiber tear out along with little damage to axial tows, resulting is a smooth undulating fracture surface. Post cure temperature excursions to 177 deg C (350 deg F) had little effect on elastic constants but reduced the strength E862 based composites severely. The strength of PR520 and 5208 based composites increased.

composite↗

NASA NTRSDate not supplied

Laser Forming of Sheet Metal for In-Space Manufacturing Applications

As the in-space economy matures, the migration of typically Earth-bound manufacturing methods, such as metal forming and joining, up to space will be required to meet the increasing demand for refueling and supply depots, habitats, laboratories, lunar surface structures, and sites for deep-space exploration shipbuilding. Only in this way can the magnitude of individual structures begin to exceed what can be launched, deployed, docked, or inflated and allow true economies of scale to be unlocked in space. However, typical terrestrial methods of forming and shaping sheet metal require massive machinery with the ability to deflect reaction forces to a large mass in a fixed reference frame (Earth). To both limit up-mass required for manufacture and to mitigate the issues of reaction forces in the relative reference frames of space, a non-contact, laser-based method of forming metal sheet is proposed. NASA Marshall Space Flight Center (MSFC), in cooperation with DARPA and the University of Florida (UF), have begun investigation of non-contact methods of laser forming sheet metal in thermal vacuum (TVAC) conditions to great success. To date, several grades of aluminum, stainless steel, and titanium coupons have been formed via laser-induced internal stresses in these metals at both ambient atmosphere and vacuum conditions. Temperatures in these tests ranged from -120 to 60 C, and though the total energy required to induce strain does vary with workpiece temperature and environmental pressure, the method remains predictable and controllable for the materials investigated. It is believed that this method will revolutionize manufacturing in-space as it not only continues to build upon the capabilities of the laser (sensing, marking, cutting, drilling, joining, and now forming) as an all-in-one tool, but also reduces the need to counteract forming forces as would be required in mechanical bending. In theory, all required stress to induce bending is generated through thermal gradient and strain within the workpiece, meaning that the need for reaction control and fuel expense on-orbit will be limited. Archimedes once stated, “Give me a lever long enough and a fulcrum on which to place it, and I shall move the world”. With a laser as our lever and our workpiece its own fulcrum, we can move the world’s manufacturing into space.

Laser↗

DOE OSTI2022

Characterization of the Fusion Zone in Laser Beam Welds on SS304L

Designing a welded joint requires consideration of the welding process, component/assembly structure and weld geometry, material composition, properties, and desired joint performance. The designer must also specify a set of acceptance criteria used to evaluate the welds. This study was conducted to support the design team responsible for specifying the spot weld profile by investigating the relationship between weld process parameters and the weld geometry and properties to further inform design optimization. The subjects of this study are two thin stainless steel 304L (SS304L) components joined by Laser Beam Welding (LBW). The components are the Housing, which contains flexible circuitry, and the Cover that protects the circuitry during handling, installation, and under service loads. Figure 1 illustrates coupon versions of the Cover and Housing used to simulate the interface between these components for weld testing. These coupon representations of assembly components were designed to match the material and relevant geometry related to the welded joint while other features not relevant to the welded joint are excluded.The objective of this study was to characterize the sensitivity of the weld geometry and mechanical properties to the input weld parameters to inform further Finite Element Analysis (FEA) efforts in balancing joint strength against assembly deformation caused by the welding process. Characterization was accomplished by varying laser welding parameters around the current design values and measuring weld geometry on the cross-sections. Nanoindentation with a spherical and flat-ended tip was used on polished section cuts of the weld to obtain direct estimates of the stress-strain response of the Fusion Zone (FZ) for comparison to the Base Metal (BM).

42 ENGINEERING↗

NASA NTRS1981

Thin film strain transducer

Previous attempts to develop an appropriate sensor for measuring the stress or strain of high altitude balloons during flight are reviewed as well as the various conditions that must be met by such a device. The design, development and calibration of a transducer which promises to satisfy the necessary design constraints are described. The thin film strain transducer has a low effective modulus so as not to interfere with the strain that would naturally occur in the balloon. In addition, the transducer has a high sensitivity to longitudinal strain (7.216 mV/V/unit strain) which is constant for all temperature from room temperature to -80 C and all strains from 5 percent compression to 10 percent tensile strain. At the same time, the sensor is relatively insensitive (0.27 percent) to transverse forces. The device has a standard 350 ohm impedance which is compatible with available bridge balance, amplification and telemetry instrumentation now available for balloon flight. Recommendations are included for improved coatings to provide passive thermal control as well as model, tethered and full scale flight testing.

Rand, J. L.↗

DOE OSTI2021

A mechanistic damage model for solid oxide fuel cell ceramic materials - Part I: Constitutive modeling

A multiscale mechanistic damage model for solid oxide fuel cell (SOFC) ceramics that combines micromechanics of stiffness reduction due to material porosity change and microcracking with a continuum damage mechanics (CDM) description for the evolution of microcracks up to fracture is developed in this work. The model also accounts for volumetric swelling of the anode due to redox cycling. Porosity change and swelling during redox combined with normal operating thermomechanical loads could increase the stress, strain and damage distributions in SOFC stacks leading them to failure. The proposed model involves three governing parameters: (i) the porosity in terms of the pore or void volume fraction, (ii) material swelling magnitude prescribed at a given loading step, and (iii) the damage variable that describes microcracking caused by thermomechanical loads and swelling. Pores and microcracks in the ceramics are modeled as randomly distributed ellipsoidal inclusions with negligible stiffness by an Eshelby-Mori-Tanaka formulation combined with an inclusion orientation distribution method. Microcracking damage evolution is described by a CDM formulation. Volumetric swelling is treated in a similar way to thermal expansion in the constitutive relations. After validation and numerical checks, the damage model is used to analyze a simplified SOFC positive-electrode/electrolyte/negative-electrode (PEN) structure subjected a redox cycle in addition to thermomechanical loads experienced during normal operation of a SOFC stack.

30 DIRECT ENERGY CONVERSION↗

NASA NTRS1989

Mesh refinement in finite element analysis by minimization of the stiffness matrix trace

Most finite element packages provide means to generate meshes automatically. However, the user is usually confronted with the problem of not knowing whether the mesh generated is appropriate for the problem at hand. Since the accuracy of the finite element results is mesh dependent, mesh selection forms a very important step in the analysis. Indeed, in accurate analyses, meshes need to be refined or rezoned until the solution converges to a value so that the error is below a predetermined tolerance. A-posteriori methods use error indicators, developed by using the theory of interpolation and approximation theory, for mesh refinements. Some use other criterions, such as strain energy density variation and stress contours for example, to obtain near optimal meshes. Although these methods are adaptive, they are expensive. Alternatively, a priori methods, until now available, use geometrical parameters, for example, element aspect ratio. Therefore, they are not adaptive by nature. An adaptive a-priori method is developed. The criterion is that the minimization of the trace of the stiffness matrix with respect to the nodal coordinates, leads to a minimization of the potential energy, and as a consequence provide a good starting mesh. In a few examples the method is shown to provide the optimal mesh. The method is also shown to be relatively simple and amenable to development of computer algorithms. When the procedure is used in conjunction with a-posteriori methods of grid refinement, it is shown that fewer refinement iterations and fewer degrees of freedom are required for convergence as opposed to when the procedure is not used. The mesh obtained is shown to have uniform distribution of stiffness among the nodes and elements which, as a consequence, leads to uniform error distribution. Thus the mesh obtained meets the optimality criterion of uniform error distribution.

Kittur, Madan G.↗

NASA NTRS1989

High temperature fatigue behavior of tungsten copper composites

The high temperature fatigue behavior of a 9 vol percent, tungsten fiber reinforced copper matrix composite was investigated. Load-controlled isothermal fatigue experiments at 260 and 560 C and thermomechanical fatigue (TMF) experiments, both in phase and out of phase between 260 and 560 C, were performed. The stress-strain response displayed considerable inelasticity under all conditions. Also, strain ratcheting was observed during all the fatigue experiments. For the isothermal fatigue and in-phase TMF tests, the ratcheting was always in a tensile direction, continuing until failure. The ratcheting during the out-of-phase TMF test shifted from a tensile direction to a compressive direction. This behavior was thought to be associated with the observed bulging and the extensive cracking of the out-of-phase specimen. For all cases, the fatigue lives were found to be controlled by damage to the copper matrix. Grain boundary cavitation was the dominant damage mechanism of the matrix. On a stress basis, TMF loading reduced lives substantially, relative to isothermal cycling. In-phase cycling resulted in the shortest lives, and isothermal fatigue at 260 C, the longest.

Verrilli, Michael J.↗

NASA NTRS1998

Constitutive properties of adult mammalian cardiac muscle cells

BACKGROUND: The purpose of this study was to determine whether changes in the constitutive properties of the cardiac muscle cell play a causative role in the development of diastolic dysfunction. METHODS AND RESULTS: Cardiocytes from normal and pressure-hypertrophied cats were embedded in an agarose gel, placed on a stretching device, and subjected to a change in stress (sigma), and resultant changes in cell strain (epsilon) were measured. These measurements were used to examine the passive elastic spring, viscous damping, and myofilament activation. The passive elastic spring was assessed in protocol A by increasing the sigma on the agarose gel at a constant rate to define the cardiocyte sigma-versus-epsilon relationship. Viscous damping was assessed in protocol B from the loop area between the cardiocyte sigma-versus-epsilon relationship during an increase and then a decrease in sigma. In both protocols, myofilament activation was minimized by a reduction in [Ca2+]i. Myofilament activation effects were assessed in protocol C by defining cardiocyte sigma versus epsilon during an increase in sigma with physiological [Ca2+]i. In protocol A, the cardiocyte sigma-versus-epsilon relationship was similar in normal and hypertrophied cells. In protocol B, the loop area was greater in hypertrophied than normal cardiocytes. In protocol C, the sigma-versus-epsilon relation in hypertrophied cardiocytes was shifted to the left compared with normal cells. CONCLUSIONS: Changes in viscous damping and myofilament activation in combination may cause pressure-hypertrophied cardiocytes to resist changes in shape during diastole and contribute to diastolic dysfunction.

NASA Discipline Cardiopulmonary↗

NASA NTRS1999

Crack Turning and Arrest Mechanisms for Integral Structure

In the course of several years of research efforts to predict crack turning and flapping in aircraft fuselage structures and other problems related to crack turning, the 2nd order maximum tangential stress theory has been identified as the theory most capable of predicting the observed test results. This theory requires knowledge of a material specific characteristic length, and also a computation of the stress intensity factors and the T-stress, or second order term in the asymptotic stress field in the vicinity of the crack tip. A characteristic length, r(sub c), is proposed for ductile materials pertaining to the onset of plastic instability, as opposed to the void spacing theories espoused by previous investigators. For the plane stress case, an approximate estimate of r(sub c), is obtained from the asymptotic field for strain hardening materials given by Hutchinson, Rice and Rosengren (HRR). A previous study using of high order finite element methods to calculate T-stresses by contour integrals resulted in extremely high accuracy values obtained for selected test specimen geometries, and a theoretical error estimation parameter was defined. In the present study, it is shown that a large portion of the error in finite element computations of both K and T are systematic, and can be corrected after the initial solution if the finite element implementation utilizes a similar crack tip discretization scheme for all problems. This scheme is applied for two-dimensional problems to a both a p-version finite element code, showing that sufficiently accurate values of both K(sub I) and T can be obtained with fairly low order elements if correction is used. T-stress correction coefficients are also developed for the singular crack tip rosette utilized in the adaptive mesh finite element code FRANC2D, and shown to reduce the error in the computed T-stress significantly. Stress intensity factor correction was not attempted for FRANC2D because it employs a highly accurate quarter-point scheme to obtain stress intensity factors.

Pettit, Richard↗

DOE OSTI2022

Comparative chemical genomic profiling across plant-based hydrolysate toxins reveals widespread antagonism in fitness contributions

Abstract The budding yeast Saccharomyces cerevisiae has been used extensively in fermentative industrial processes, including biofuel production from sustainable plant-based hydrolysates. Myriad toxins and stressors found in hydrolysates inhibit microbial metabolism and product formation. Overcoming these stresses requires mitigation strategies that include strain engineering. To identify shared and divergent mechanisms of toxicity and to implicate gene targets for genetic engineering, we used a chemical genomic approach to study fitness effects across a library of S. cerevisiae deletion mutants cultured anaerobically in dozens of individual compounds found in different types of hydrolysates. Relationships in chemical genomic profiles identified classes of toxins that provoked similar cellular responses, spanning inhibitor relationships that were not expected from chemical classification. Our results also revealed widespread antagonistic effects across inhibitors, such that the same gene deletions were beneficial for surviving some toxins but detrimental for others. This work presents a rich dataset relating gene function to chemical compounds, which both expands our understanding of plant-based hydrolysates and provides a useful resource to identify engineering targets.

09 BIOMASS FUELS↗

DOE OSTI2024

Optimal Geometry for Focused Ion Beam-Milled Samples for Direct-Pull Micro-Tensile Testing Performed In Situ in a Scanning Electron Microscope

A thorough procedure was developed to efficiently manufacture dogbone samples using focused ion beam (FIB) milling for micro-tensile testing. A Bruker PI 89 PicoIndenter, Billerica, MA, USA, was used as a case study, although the analysis and results are applicable to other micro-mechanical testing systems capable of mounting a standard, Ø12.7 mm × Ø3.2 mm pin, scanning electron microscopy (SEM) pin stub (Ted Pella, Redding, CA, USA). Nine dogbones were made from an Fe-45Cu alloy additively manufactured using powder-fed laser-directed energy deposition (DED-LB). Testing showed that fracture was confined to the gauge section for all dogbones and that the fracture mode, ductile vs. brittle, was entirely dependent on the grain orientation relative to the loading direction. The analysis showed that the measured plastic strain to failure can vary from >11% (optimal geometry) to <1% (non-optimal geometry) in micro-tensile testing of high-tensile-strength (>1 GPa) metallic materials. Subsequently, a finite element analysis (FEA) was conducted to identify the improved dogbone geometries. A total of ten thousand dogbone geometries were tested, and their dimensions were defined by a set of four adjustable parameters (corner radius, load surface angle, load surface length, and dogbone head length). The gauge width and gauge length were fixed to 4 µm and 10 µm, respectively. Three-dimensional surface plots of the stress concentration as a function of two parameters were used to identify the optimal ranges of parameter values. The addition of maximum width and length constraints, measuring 25 µm and 30 µm, respectively, allowed us to identify an optimal geometry at load surface angles of 30° and 45°. Their respective dimensions (corner radius, load surface length, and dogbone head length) are, in µm, 12, 6, and 7 and 10, 7, and 7. Testing these two optimal geometries with a range of gauge lengths from 4 to 20 µm showed that smaller gauge lengths only slightly reduced the detrimental stress concentration outside the gauge section. However, smaller gauge lengths will notably improve the FIB surface polishing step as tapering is reduced with smaller dogbone lengths.

Chemistry↗

NASA NTRS1993

A nonlinear theory for spinning anisotropic beams using restrained warping functions

A geometrically nonlinear theory is developed for spinning anisotropic beams having arbitrary cross sections. An assumed displacement field is developed using the standard 3D kinematics relations to describe the global beam behavior supplemented with an additional field that represents the local deformation within the cross section and warping out of the cross section plane. It is assumed that the magnitude of this additional field is directly proportional to the local stress resultants. In order to take into account the effects of boundary conditions, a restraining function is introduced. This function plays the role of reducing the amount of free warping deformation throughout the field due to the restraint of the cross section(s) at the end(s) of the beam, e.g., in the case of a cantilever beam. Using a developed ordering scheme, the nonlinear strains are calculated to the third order. The FEM is developed using the weak form variational formulation. Preliminary interesting numerical results have been obtained that indicate the role of the restraining function in the case of a cantilever beam with circular cross section. These results are for the cases of a tip displacement (static) and free vibration studies for both isotropic and anisotropic materials with varied fiber orientations.

Ie, C. A.↗

DOE OSTI2026

Characteristics of Fluid‐Solid Interaction Constitutive Models Within Poroelastodynamics at Higher Strain‐Rates and Large Deformations Implemented in 1D

The large deformation, mixed formulation, finite element (FE) modeling approach presented in Irwin et al. 2024 is extended herein to include improved constitutive models for representing dynamic solid-fluid interactions at higher strain rates (𝒪⁢(1⁢0 2 −1⁢0 3 )⁢s −1 ) and larger overpressure magnitudes (𝒪⁡(1⁢0 2 )⁢kPa) within a biphasic soft porous material using Theory of Porous Media (TPM) at finite strain. Specifically, these constitutive modeling improvements are the following: (i) a more physically robust constitutive model for pore fluid seepage velocity via inclusion of pore fluid viscous stress, and (ii) a modified deformation-dependent-permeability model and updated hyperelastic constitutive model better suited for handling larger volumetric compressions and extensions. The novelty of the present work is mainly the contribution (i): inclusion of pore fluid viscous stress at higher strain-rate and large deformations, which requires 𝐶 1 continuity in the weak formulation, accomplished by employing Hermite cubic interpolation functions within a mixed nonlinear poromechanical finite element formulation. In (ii), the model is updated to weakly enforce solid phase incompressibility, such that this assumption is not violated numerically, which provides improved numerical stability for achieving larger overpressure magnitudes on 𝒪⁡(1⁢0 2 ) kPa, which were not achievable with the previous Kozeny–Carman model in Irwin et al. 2024. Also in (ii), the volumetric part of the solid skeleton free energy function is modified to ensure proper bounds on the solid skeleton Jacobian of deformation 𝐽 s related to incompressibility of the solid phase. Uniaxial strain, unidirectional flow examples at higher strain rates (𝒪⁢(1⁢0 2 −1⁢0 3 )⁢s −1 ) and larger deformations (up to 0.2 (or 20%) nominal axial strain) demonstrate the improved physical representation—and numerical stability—of these constitutive model improvements.

42 ENGINEERING↗

NASA NTRS1996

Grazing-Activated Production of Dimethyl Sulfide (DMS) by two clones of Emiliania huxleyi

Emiliania huxleyi clones CCMP 370 and CCMP 373 produced similar amounts of dimethylsulfoniopropionate (DMSP) during axenic exponential growth, averaging 109 mM internal DMSP. Both clones had detectable DMSP lyase activity, as measured by production of dimethyl sulfide (DMS) during in vitro assays of crude cell preparations, but activities and conditions differed considerably between clones. Clone 373 had high activity; clone 370 had low activity and required chloride. For both strains, enzyme activity per cell was constant during exponential growth, but little DMS was produced by healthy cells. Rather, DMS production was activated when cells were subjected to physical or chemical stresses that caused cell lysis. We propose that DMSP lyase and DMSP are segregated within these cells and re-action only under conditions that result in cell stress or damage. Such activation occurs during microzooplankton grazing. When these clones were grazed by the dinoflagellate Oxyrrhis marina, DMS was produced; ungrazed cells, as well as those exposed to grazer exudates and associated bacteria, generated no DMS. Grazing of clone 373 produced much more DMS than grazing of clone 370, consistent with their relative in vitro DMSP lyase activities. DMS was only generated when cells were actually being grazed, indicating that ingested cells were responsible for the DMS formation. We suggest that even low levels of grazing can greatly accelerate DMS production.

Wolfe, Gordon V.↗