Search NASA⌕ Search

SEARCH · Search NASA

Results for “STRESS CALCULATION”

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.

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

At least 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Discovery of multimechanisms of screw dislocation interaction in bcc iron from open-ended saddle point searches

Dislocation motion and interactions determine mechanical properties in body-centered cubic (bcc) metallic materials. However, studying mechanisms for the screw dislocation interaction is fundamentally challenging since many underlying processes involve mesotimescales and atomistic resolution, currently inaccessible by either experimental techniques or continuum theoretical methods. Here, in this paper, we develop a computational capability based on self-evolving atomistic kinetic Monte Carlo (SEAKMC) to sample the critical events and saddle point energies related to screw dislocations and their junctions. The method is first validated by calculating the stress dependence of Peierls barriers and formation energies of kink pairs and cross-slip kink pairs on a single screw dislocation in bcc iron. Then the method is applied to a binary junction of a pair of intersecting screw dislocations, the structure of which is crucial for low-temperature plastic deformation. We identify three important mechanisms: coplanar cross-slipping, jog-pinning, and a previously unknown unzipping mechanism during the evolution of the binary junction. The mechanisms are then further validated using classical molecular dynamics simulations. The computational capability developed in this paper provides an effective tool to evaluate screw dislocation related thermally activated events in complex stress conditions. The mechanisms discovered in this paper provide critical insights into temperature dependence of the anomalous slip, a breakdown of the Schmidt law, during the plastic deformation in bcc iron and can be generalized to other bcc metals.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Understanding plasticity in irradiated alloys through TEM in situ compression pillar tests

Depth-sensing transmission electron microscopic (TEM) in situ mechanical testing has become widely utilized for understanding deformation in irradiated materials. Until now, compression pillars have primarily been used to study the elastic properties and yield of irradiated materials. In this study, we utilize TEM in situ compression pillars to investigate plastic deformation in two ion-irradiated alloys: Fe–9% Cr oxide dispersion strengthened (ODS) alloy and nanocrystalline Cu–24% Ta. We develop an algorithm to automate the extraction of instantaneous pillar dimensions from TEM videos, which we use to calculate true stress–strain curves and strain hardening exponents. True stress–strain curves reveal intermitted plastic flow in all specimen conditions. In the Fe–9% Cr ODS, intermitted plastic flow is linked to strain bursts observed in TEM videos. Low strain hardening or strain softening is observed in all specimen conditions. TEM videos link the strain softening in irradiated Fe–9% Cr ODS to dislocation cross-slip, and in Cu–24% Ta to grain boundary sliding.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Stress from long-range interactions in particulate system

In particulate systems, the ensemble average accounting for effects from all particles is expressed as the expected value related to the nearest particle probability. Using this expression, a stress is defined even for long-range particle interactions without the divergence difficulty because of the rapid far field decay of the probability. Similar to the potential part of the virial stress in molecular systems, this stress represents particle-particle interactions through the surrounding field. As examples, the stresses are calculated for systems of charged particles and disperse multiphase flows.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2012

Thermal Analysis to Calculate the Vessel Temperature and Stress in Alcator C-Mod Due to the Divertor Upgrade

Alcator C-Mod is planning an upgrade to its outer divertor. The upgrade is intended to correct the existing outer divertor alignment with the plasma, and to operate at elevated temperatures. Higher temperature operation will allow study of edge physics behavior at reactor relevant temperatures. The outer divertor and tiles will be capable of operating at 600oC. Longer pulse length, together with the plasma and RF heat of 9MW, and the inclusion of heater elements within the outer divertor produces radiative energy which makes the sustained operation much more difficult than before. An ANSYS model based on ref. 1 was built for the global thermal analysis of C-Mod. It models the radiative surfaces inside the vessel and between the components, and also includes plasma energy deposition. Different geometries have been simulated and compared. Results show that steady state operation with the divertor at 600oC is possible with no damage to major vessel internal components. The differential temperature between inner divertor structure, or "girdle" and inner vessel wall is ~70oC. This differential temperature is limited by the capacity of the studs that hold the inner divertor backing plates to the vessel wall. At a 70oC temperature differential the stress on the studs is within allowable limits. The thermal model was then used for a stress pass to quantify vessel shell stresses where thermal gradients are significant.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Advancements in 3D Structural Analysis of Geothermal Systems

Robust geothermal activity in the Great Basin, USA is a product of both anomalously high regional heat flow and active fault-controlled extension. Elevated permeability associated with some fault systems provides pathways for circulation of geothermal fluids. Constraining the local-scale 3D geometry of these structures and their roles as fluid flow conduits is crucial in order to mitigate both the costs and risks of geothermal exploration and to identify blind (no surface expression) geothermal resources. Ongoing studies have indicated that much of the robust geothermal activity in the Great Basin is associated with high density faulting at structurally complex fault intersection/interaction areas, such as accommodation/transfer zones between discrete fault systems, step-overs or relay ramps in fault systems, intersection zones between faults with different strikes or different senses of slip, and horse-tailing fault terminations. These conceptualized models are crucial for locating and characterizing geothermal systems in a regional context. At the local scale, however, pinpointing drilling targets and characterizing resource potential within known or probable geothermal areas requires precise 3D characterization of the system. Employing a variety of surface and subsurface data sets, we have conducted detailed 3D geologic analyses of two Great Basin geothermal systems. Using EarthVision (Dynamic Graphics Inc., Alameda, CA) we constructed 3D geologic models of both the actively producing Brady’s geothermal system and a ‘greenfield’ geothermal prospect at Astor Pass, NV. These 3D models allow spatial comparison of disparate data sets in 3D and are the basis for quantitative structural analyses that can aid geothermal resource assessment and be used to pinpoint discrete drilling targets. The relatively abundant data set at Brady’s, ~80 km NE of Reno, NV, includes 24 wells with lithologies interpreted from careful analysis of cuttings and core, a 1:24,000 scale detailed geologic map and cross-sections, 2D seismic reflection profiles and other geophysical data, and downhole temperature data. The 3D geologic model based on these data consists of 61 fault planes, 25 distinct stratigraphic units, and 2 intrusive bodies. Geothermal fluids are produced from a left step-over/relay ramp within the Brady’s Fault Zone (BFZ). Under local stress conditions, fault segments that strike NNE-to-NE are most likely to slip and/or dilate, and therefore transmit geothermal fluids. The 3D model defines the locations of discrete fault intersections within the BFZ and indicates that the densest zones of structurally controlled fracture permeability are ~10-to-10s of meters in diameter and plunge ~55° NW-NNW beneath the heart of the BFZ step over. The locations of high intersection density, high fault slip and dilation tendency, high subsurface temperature, and lithologies known to support high fracture permeability are combined to produce 3D ‘fairway’ maps useful in both assessments of geothermal resource potential and for defining drilling targets. Astor Pass is located on the Pyramid Lake Paiute Reservation, ~80 km north of Reno, NV. It is a prospective ‘greenfield’ geothermal area, and thus subsurface data are relatively sparse. Available data include: two relatively deep wells (~1400 m) and one shallower well (~500 m) with lithologies interpreted from drill cuttings, several 2D seismic reflection profiles, a 1:24,000 scale geologic map and cross-section, a shallow temperature survey, and downhole temperature data. 3D modeling based on these data has defined 19 distinct fault planes and 16 stratigraphic units. Based on the stress field calculated from borehole breakouts, drilling induced tensile cracks and petal-centerline cracks in the two relatively deep wells, 3D slip and dilation tendency analysis indicates that northerly striking fault segments are most likely to slip and/or dilate, and therefore transmit geothermal fluids. Analysis of fault intersection density indicates that the highest density of structurally controlled permeability within the field lies in a narrow (10-to-10s of m) zone plunging moderately (~35°) to the NNW beneath Pleistocene tufa deposits. This zone of increased fracture density, which we interpret as the primary upflow zone, is controlled by the intersection of N-to-NNW striking normal faults and a WNW striking dextral fault zone and represents the most promising target for future drilling. Construction of a 3D geologic model involves integration of a variety of data into an internally consistent framework. A robust model allows for spatial comparison between the various types of data (structural, stratigraphic, geophysical, temperature, etc.) that are commonly used independently to site geothermal wells. Furthermore, highly detailed 3D geologic models provide the basis for additional quantitative analysis, including 3D fault slip and dilation tendency analysis and the precise location of structurally controlled permeability pathways. These analyses provide detailed information relating to the internal dynamics of geothermal systems and can mitigate the costs and risks of geothermal exploration and development by contributing to better well targeting and more accurate evaluations of resource potential.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

NuScale Pressure and Temperature Limits Methodology Using Finite Element Analysis

Per 10 CFR 50 Appendix G, the pressure-temperature (P-T) limits curves and minimum temperature must be established to provide adequate margins for ferritic pressure-retaining components of the reactor coolant pressure boundary; this is to protect against brittle failure during any normal operating conditions, including anticipated operational occurrences and system hydrostatic tests, to which the pressure boundary may be subjected over its service lifetime. Specifically, ASME Code Section XI Appendix G procedures must be used for P-T limits calculation considering the pressure and temperature at various operating transient conditions. However, the elastic fracture mechanics solutions in Section XI are only suitable for cylindrical reactor pressure vessel (RPV) beltline without geometric discontinuities. Hence, these solutions are not suitable for postulated flaws near the core support blocks attached to the NuScale RPV inside surface, which is part of the beltline. As a result, NuScale has used finite element analysis to calculate thermal stress and stress intensity factor for the postulated flaws. The solutions using finite-element analysis have been validated using the formulations for straight cylinders for both axial and circumferential flaws. In addition, special considerations are given to the RPV beltline nil-ductility transition temperature (RTNDT) due to neutron irradiation occurring at lower temperatures than conventional plants. The paper summarizes the methodology and finite-element models used to develop P-T limits curves for NuScale RPV at the end of its 60-year design life.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2017

Nondestructive Evaluation of Alkali-Silica Reaction in High-Strength Concrete for Aging Structures Sustainability

Concrete is a multiphase geo-material, which is utilized extensively in passive supporting and shielding structures at nuclear facilities. Chemo-mechanical degradation of concrete in old structures is a challenge to sustainability of aging commercial light water reactors. Alkali-silica reaction (ASR) is a common slowly evolving chemo-mechanical damage mechanism of concrete, which could take tens of years to reach the threshold stage when visible cracks appear in a concrete structure. In ASR, concrete constituents react with each other, namely the alkalis in the pore solution of hydrated cement paste (mostly K{sup +}OH{sup -} and Na{sup +}OH{sup -}) react with non-crystalline silica (SiO{sub 2}) of the aggregates to form highly hydrophilic gel. At sufficient levels of ambient relative humidity (typically above 85%), imbibition of water causes expansion of the ASR gel, which results in cracking of constrained reinforced concrete and ultimate structural weakening below safety margin levels. Stability and functional performance of concrete structures potentially suffering from chemo-mechanical deterioration due to ASR is evaluated in a multi-stage process involving material evaluation and structural engineering analysis. Existing methods used by industry for concrete material evaluation are based on extrusion of cores from multiple locations in a structure, and destructive laboratory tests on cores. Because each destructive test requires separate core extrusion, the total number of cores taken from the structure at once could be numbering in hundreds. Laboratory tests include petrography (high resolution optical microscopy of core sections) to develop damage rating index (DRI) score of concrete cores. Another test may involve chemical pore solution extraction from cores, and solution speciation analysis using inductively coupled plasma and mass spectroscopy (ICP/MS). Other tests may consist of mechanical evaluation of cores strength and fracture toughness. The DRI score of concrete material as well as other information obtained from the laboratory tests is subsequently used in structural engineering calculations, where stresses in the entire structure are modeled to evaluate its safety and mechanical stability. Extrusion of a large number of cores is highly undesirable for evaluation of nuclear facility concrete structures. The amount of destructive material testing could be minimized by developing methods for in-situ NDE of ASR to quantify the amount of deterioration that is occurring in the structure, and to identify areas in the structure where more deterioration may be occurring.

36 MATERIALS SCIENCE↗
DOE OSTI2022

ELASTOMERIC MICROVASCULAR SELF-HEALING MATERIALS

Damaged elastomeric diaphragms within pneumatic controllers used in the oil and gas industry lead to an unintended release of methane. Self-healing microvascular materials capable of healing various types of damage have been fabricated. These microvascular materials are designed to replace currently available commercial diaphragms found in pneumatic controllers and provide a solution to reduce unintended methane leaks. Poly(dimethylsiloxane) (PDMS) was used as the main matrix material, with additional testing conducted on polyurethane and flexibilised epoxy materials. Microvascular networks were implemented into the elastomeric membranes to act as ves- sels to deliver healing agents to the damaged areas. The complex, interdigitated channel networks were created using a 3D printed custom compounded filament composed of polylactic acid (PLA) and tin(II) oxalate. Two-part liquid casting polymers were poured around the channels into a mould to create samples with the scaffolding of the microvascular channels intact. Hollow microvascular networks were created by placing the samples within a vacuum oven at 250◦C which causes the thermal depolymerisation of PLA into its gaseous monomers. The manufactured materials were placed in test stands and pressurized using nitrogen gas to determine their healing and mechanical behaviour. Self-healing behaviour was demonstrated using PDMS matrix materials and healing chemistries within a pressure test cell. This test cell was designed to detect any damage to a sample by record- ing an outlet pressure. Damages including puncture holes, diagonal cuts and star-shaped central iv cuts were applied to the samples and a reduction in the outlet pressure was recorded for all sam- ples. A variation in channel spacing and diameter was studied to determine the optimal design of the microvascular network for self-healing performance. The mechanical performance of elastomeric membrane materials within a diaphragmatic pressurised stand was observed. Deformations and strains around the channels and channel inter- sections at the midplane of the materials were recorded using digital image correlation. Normal stresses were calculated using Hooke’s law and the material properties of the PDMS matrix. A reduction in leak rate, implying the success of self-healing, was recorded within manu- factured diaphragm samples tested within a commercial valve. PDMS is the favourable material when compared to polyurethane and flexibilised epoxy for creating

03 NATURAL GAS↗
DOE OSTI2019

The role of diffusion-driven pure climb creep on the rheology of bridgmanite under lower mantle conditions

The viscosity of Earth’s lower mantle is poorly constrained due to the lack of knowledge on some fundamental variables that affect the deformation behaviour of its main mineral phases. This study focuses on bridgmanite, the main lower mantle constituent, and assesses its rheology by developing an approach based on mineral physics. Following and revising the recent advances in this field, pure climb creep controlled by diffusion is identified as the key mechanism driving deformation in bridgmanite. The strain rates of this phase under lower mantle pressures, temperatures and stresses are thus calculated by constraining diffusion and implementing a creep theoretical model. The viscosity of MgSiO3 bridgmanite resulting from pure climb creep is consequently evaluated and compared with the viscosity profiles available from the literature. We show that the inferred variability of viscosity in these profiles can be fully accounted for with the chosen variables of our calculation, e.g., diffusion coefficients, vacancy concentrations and applied stresses. A refinement of these variables is advocated in order to further constrain viscosity and match the observables.

58 GEOSCIENCES↗
DOE OSTI2015

RECENT PROGRESS OF CRACK BRIDGING MODELING OF DUCTILE-PHASE-TOUGHENED W-CU COMPOSITES

A crack bridging model using calculated Cu stress-strain curves has been developed to study the toughening of W-Cu composites. A strengthening factor and necking parameters have been added to the model for the ductile-phase bridges to incorporate constraint effects at small bridge sizes. Parametric studies are performed to investigate the effect of these parameters. The calculated maximum applied stress intensity, aKmax, to induce a 1-mm stable crack is compared to the experimental stress intensity at peak load, Kpeak. Without bridge necking, increasing the strengthening factor improves the agreement between aKmax and Kpeak when plotted vs. logarithm of the displacement rate. Improvement can also be achieved by allowing necking with a larger failure strain. While the slope is better matched with this latter approach, the calculated value of aKmax is significantly larger than Kpeak.

Setyawan, Wahyu↗
DOE OSTI2017

Design Guidelines for Predicting Stress in Cemented Doublets Undergoing Temperature Change

This work explores quick predictive methods for calculating potentially risky stresses in cemented doublets underdoing temperature change that agree well with finite element analysis. We also provide guidelines for avoiding stress concentrations. Cemented doublets are ubiquitous in the optics world, but it is not always clear if they will perform as intended in a specific thermal environment due to excessive stress. The stress in the optic can cause debonding or even fracture. Yoder addressed this problem using analytical results from Chen and Nelson but also expressed some skepticism as to the accuracy and usefulness of those results. Closed-form analytical solutions are extremely useful when doing trade-off studies to see if a particular geometry or glass choice will cause stress problems. The problem of bonded materials undergoing temperature change is of great interest in the electronics industry, and improved analytical solutions are now available that address the issues raised by Yoder, and they have been validated using finite element analysis (FEA). But the newer solutions still have limitations. The analytical solutions presented here cannot predict stresses due to stress concentrations or singularities at the edge of the bondline—this still has to be done using FEA—although some advanced methods are available. Also, the existing analytical methods assume a flat interface between adjoining elements; they do not account for the effect of curvature. We find that most doublet failures start at the edges and involve debonding of the glass. Glass fracture is rare for doublets, accounting for less than 10% of the failures. This paper provides the reader with a general introduction to a large and dense body of literature in order to demystify the subject of temperature-induced glass bonding failures.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Characterization of microstructure and residual stress following the friction stir welding of dissimilar aluminum alloys

Friction-stir-welding (FSW) has gained importance as an effective way to join dissimilar materials due to its solid-state nature. However, FSW produces substantial residual stresses in aluminum welding, often deteriorating the component's structural performance and reliability. In the current study, the d-spacing (inter-atomic distance) of several locations on a battery tray and stress-free reference samples are measured with neutron diffraction to calculate the residual stresses in a lap 6061 to A365 dissimilar friction-stir-weld. Furthermore, the present study connects the measured residual stress to microstructure properties through electron backscatter diffraction (EBSD) and optical microscopy analyses. The EBSD results indicate that the FSW operation affects the grains’ orientation, thereby producing location-specific directionally varying elasticity across the weldment. Further, this varied elasticity creates unique d-spacing profiles even in the stress-free material. Additionally, the resulting stressed and stress-free d-spacing profiles show dependencies on the FSW tool geometry, traverse direction, and rotational direction. Lastly, this study suggests that the final component residual stress can be influenced by the sequence of multi-pass welds contributing to a rise in temperature in the presence of a completed weld relaxing its residual stress.

36 MATERIALS SCIENCE↗
DOE OSTI2012

Cryogenic system design of 11 GEV/C super high momentum spectrometer superconducting magnets at Jefferson Lab

The design of the cryogenic system for the 11 GeV/c Super High Momentum Spectrometer (SHMS) is presented. A description of the cryogenic control reservoir and the cryogenic transfer line is given. Details of the cryogenic control reservoirs, cryogenic transfer lines, and pressure piping are summarized. Code compliance is ensured through following the requirements of the ASME Pressure Vessel Code and Pressure Piping Code. An elastic-plastic-analysis-based combined safety factor approach is proposed to meet the low stress requirement of ASME 2007 Section VIII, Division 2 so that Charpy V-notch (CVN) impact testing can be avoided through analysis. Material toughness requirements in ASME 2007 Section VIII, Division 2 are adopted as CVN impact testing rules of stainless steel 304 piping at 4.2 K and 77 K. A formula-based combined safety factor approach for pressure piping is also proposed to check whether the impact testing can be avoided due to low stress. Analysis and calculation have shown that no CVN impact testing of base metal and heat affected zones is required for the helium reservoir, nitrogen reservoir, and their relevant piping. Total heat loads to liquid helium and liquid nitrogen are studied also. The total heat load to LHe for SHMS is estimated to be 137 W, and the total load to LN2 is calculated to be 420 W.

36 MATERIALS SCIENCE↗
DOE OSTI2019

A higher-order Lagrangian discontinuous Galerkin hydrodynamic method for solid dynamics

We present a new multidimensional high-order Lagrangian discontinuous Galerkin (DG) hydrodynamic method that supports hypoelastic and hyperelastic strength models for simulating solid dynamics with higher-order elements. We also present new one-dimensional test problems that have an analytic solution corresponding to a hyperelastic–plastic wave. A modal DG approach is used to evolve fields relevant to conservation laws. These fields are approximated high-order Taylor series polynomials. The stress fields are represented using nodal quantities. The constitutive models used to calculate the deviatoric stress are either a hypoelastic–plastic, infinitesimal strain hyperelastic–plastic, or finite strain hyperelastic–plastic model. These constitutive models require new methods for calculating high-order polynomials for the velocity gradient and deformation gradient in an element. The plasticity associated with the strength model is determined using a radial return method to with a J 2 yield criterion and perfect plasticity. The temporal evolution of the governing equations is achieved with the total variation diminishing Runge–Kutta (TVD RK) time integration method. A diverse suite of 1D and 2D test problems are calculated. Here, the new 1D piston test problems, which have analytic solutions for each elastic–plastic model, are presented and calculated to demonstrate the stability and formal accuracy of the various models with the new Lagrangian DG method. 2D test problems are calculated to demonstrate the stability and robustness of the new Lagrangian DG method on multidimensional problems with high-order elements, which have faces that can bend.

42 ENGINEERING↗
DOE OSTI2014

On the calculation of internal forces in mechanically stressed polyatomic molecules

We discuss how to define and to compute internal forces in a molecule subjected to mechanical stress. Because of the inherently many-body character of intramolecular interactions, internal forces cannot be uniquely defined without specifying a set of internal coordinates used to describe the molecular structure. When such a set is comprised of 3N − 6 interactomic distances (N being the number of atoms) and includes the bond lengths of interest, we show that the associated forces, while satisfying the equation F = ∂V/∂R (where R is the bond length, F is the internal force in this bond, and V is the potential energy of the molecule), can be determined from the molecular geometry alone. We illustrate these ideas using several toy models ranging from small molecules to a graphene sheet and show that the magnitude of the internal force in a bond is not necessarily a good predictor of its strength in response to mechanical loading. At the same time, analysis of internal forces reveals interesting phenomena such as the force multiplication effect, where weak external forces may, e.g., be used to break strong bonds, and offers insight into the catch-bond phenomenon where chemical reactivity is suppressed through application of a force.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Pulse laser induced graphite-to-diamond phase transition: the role of quantum electronic stress

First-principles calculations show that the pulse laser induced graphite-to-diamond phase transition is related to the lattice stress generated by the excited carriers, termed as “quantum electronic stress (QES)”. We found that the excited carriers in graphite generate a large anisotropic QES that increases linearly with the increasing carrier density. Using the QES as a guiding parameter, structural relaxation spontaneously transforms the graphite phase into the diamond phase, as the QES is reduced and minimized. Our results suggest that the concept of QES can be generally applied as a good measure to characterize the pulse laser induced phase transitions, in analogy to pressure induced phase transitions.

Wang, ZhengFei↗
DOE OSTI2018

Dynamic Recrystallization Behavior of AISI 422 Stainless Steel During Hot Deformation Processes

In this work, hot compression tests were performed to investigate the dynamic recrystallization (DRX) process of a martensitic stainless steel (AISI 422) at temperatures of 950, 1000, 1050, 1100 and 1150 °C and strain rates of 0.01, 0.1 and 1 s{sup −1}. The dependency of strain-hardening rate on flow stress was used to estimate the critical stress for the onset of DRX. Accordingly, the critical stress to peak stress ratio was calculated as 0.84. Moreover, the effect of true strain was examined by fitting stress values to an Arrhenius type constitutive equation, and then considering material constants as a function of strain by using a third-order polynomial equation. Finally, two constitutive models were used to investigate the competency of the strain-dependent constitutive equations to predict the flow stress curves of the studied steel. It was concluded that one model offers better precision on the flow stress values after the peak stress, while the other model gives more accurate results before the peak stress.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Nonlinear parallel momentum transport in strong electrostatic turbulence

Most existing theoretical studies of momentum transport focus on calculating the Reynolds stress based on quasilinear theory, without considering the nonlinear momentum flux-〈v{sup ~}{sub r}n{sup ~}u{sup ~}{sub ∥}〉. However, a recent experiment on TORPEX found that the nonlinear toroidal momentum flux induced by blobs makes a significant contribution as compared to the Reynolds stress [Labit et al., Phys. Plasmas 18, 032308 (2011)]. In this work, the nonlinear parallel momentum flux in strong electrostatic turbulence is calculated by using a three dimensional Hasegawa-Mima equation, which is relevant for tokamak edge turbulence. It is shown that the nonlinear diffusivity is smaller than the quasilinear diffusivity from Reynolds stress. However, the leading order nonlinear residual stress can be comparable to the quasilinear residual stress, and so may be important to intrinsic rotation in tokamak edge plasmas. A key difference from the quasilinear residual stress is that parallel fluctuation spectrum asymmetry is not required for nonlinear residual stress.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗