SEARCH · Engineering Papers

Results for “Conditional diffusion models”

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 289 records · Page 16

All sources NASA NTRS DOE OSTI

DOE OSTI2014

Derivation of a continuum model and the energy law for moving contact lines with insoluble surfactants

A continuous model is derived for the dynamics of two immiscible fluids with moving contact lines and insoluble surfactants based on thermodynamic principles. The continuum model consists of the Navier-Stokes equations for the dynamics of the two fluids and a convection-diffusion equation for the evolution of the surfactant on the fluid interface. The interface condition, the boundary condition for the slip velocity, and the condition for the dynamic contact angle are derived from the consideration of energy dissipations. Different types of energy dissipations, including the viscous dissipation, the dissipations on the solid wall and at the contact line, as well as the dissipation due to the diffusion of surfactant, are identified from the analysis. A finite element method is developed for the continuum model. Numerical experiments are performed to demonstrate the influence of surfactant on the contact line dynamics. The different types of energy dissipations are compared numerically.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Quantitative Prediction of Perovskite Stability Using Accelerated Testing and Machine Learning

PV technologies based on hybrid perovskites offer the potential for reducing solar cell costs, but they are particularly vulnerable to degradation by environmental factors such as moisture, oxygen, and illumination. Commercialization will require not only stable materials and device architectures but also accelerated testing protocols and models that can predict degradation from the accelerated testing data. Here, we report results from in situ photoluminescence (PL) and photoconductivity (PC) measurements during perovskite degradation with simultaneous optical transmittance (Tr) measurements or reflected dark field (DF) imaging. From PL, PC, Tr, and DF, we determine (respectively) the steady-state quasi-Fermi level splitting, the mean effective carrier diffusion length, the extent of conversion of perovskite to higher bandgap degradation products, and the extent of scattering from domains with different orientation or composition, all as a function of time during degradation. Simultaneous measurement of PL-PC-Tr or PL-PC-DF on perovskite absorbers in an environmental chamber over a wide range of humidity, oxygen, temperature, and illumination levels yields a rich data set. We use machine learning to develop a model that accurately predicts (within 10%) the time for the diffusion length to decrease to 85% of its initial value. The model takes the environmental conditions and the first few measurements from the PL-PC-Tr or PL-PC-DF experiment as input. Thus, the model provides a framework to interpret the results of accelerated testing of absorber materials. One of the dominant features in the model of degradation for CH 3 NH 3 PbI 3 is the initial rate at which transmittance increase. For devices with opaque contacts, the pixel-averaged rate of change of the intensity in dark-field images can be used in place of transmittance. Results are also presented on full PV devices with in-situ current-voltage (JV) measurements. Here, the data from simultaneous PL-DF-JV under environmental stresses reveal intimate connection between degradation and shunts and provide a framework to extend accelerated testing to devices. Further, the presentation reveals underlying universal behavior in degradation pathways over a broad range of environmental stresses and perovskite compositions.

14 SOLAR ENERGY↗

DOE OSTI2021

Thermal effects mediating the flow induced by laser-induced optical breakdown

Here, a short focused laser pulse can generate a hot plasma, which expands then contracts and can eject a hot jet, the speed and direction of which is sensitive to the details of the plasma kernel. The coupling of thermal and chemical nonequilibrium plasma mechanisms with this hydrodynamic development is assessed with detailed simulations of a two-temperature, three-species plasma model. Time scales for electron recombination, thermal relaxation, and diffusion are compared to that of the plasma expansion to anticipate conditions in which these mechanisms might affect the vorticity generation that leads to the ultimate flow pattern. The effect of these mechanisms are analyzed through comparison with corresponding inert-gas and equilibrium models. Thermal-nonequilibrium effects are found to be weak due to rapid relaxation of the heavy-particle and electronic temperatures. In contrast, chemical equilibration occurs at a rate comparable to the expansion and thereby enhance both it and subsequent hydrodynamic mechanisms as the energy stored in ion formation during the preceding breakdown is released by electron recombination. Thermal conduction, enhanced by high-temperature free electrons, weakens the ejection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2013

Removal of water and iodine by solid sorbents: adsorption isotherms and kinetics

Tritium and iodine-129 are two major radioactive elements that are present in off-gases from spent fuel reprocessing plants. Adsorption by solid sorbents is the state-of-the-art technique for removal of these species from off-gases. Modeling and simulating adsorption processes require accurate adsorption equilibrium and kinetic data to permit reasonable estimates of process parameters. We have developed a continuous flow single-pellet adsorption system to gather accurate adsorption equilibrium and kinetic data for adsorption of water by molecular sieve 3A and for adsorption of iodine by silver exchanged mordenite. In this paper, the design of the water and iodine adsorption experimental systems are briefly described and results of water adsorption experiments are presented and discussed. Water uptake curves are fitted with the linear-driving force (LDF) model and the shrinking-core model to determine kinetic parameters. It is shown that the kinetics of water adsorption on zeolite 3A under current experimental conditions is controlled by both the external film resistance and the macro-pore diffusion and can be predicted by both the LDF model and the shrinking-core model with the former one performing slightly better. Preliminary results from iodine adsorption experiments will be presented in the conference.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Designing Particle Morphologies for Materials with Solid Transport Limitations: A Case Study of Lithium and Manganese Rich Cathode Oxides

A lithium and manganese rich nickel-manganese-cobalt oxide (LMR-NMC) cathode is a promising candidate for next-generation batteries due to its high specific capacity, low cost, and low cobalt content. However, the material suffers from poor rate capability due to the diffusion limitations of lithium in the cathode particles. Understanding the material performance requires careful control of the morphology of the cathode particles, taking into account the primary and agglomerated diffusion pathways and the presence of pores, some of which could be closed from electrolyte infiltration. Here, in this study, we use a microstructure-based mathematical model combined with experimental data to understand the role of the complex cathode particle morphology in the rate performance of the material. Scanning electron microscopy images of cathodes made under different synthesis conditions, which results in different agglomerate morphologies, serve as the input into the mathematical model. The model is then compared to rate data to understand the controlling parameters. The presence of intra-agglomerate closed pores results in a large agglomerate diffusion length in comparison to the ideal condition, where the primary particles are agglomerated in an open and dispersed manner such that the entire interfacial area is available for electrochemical reaction. Smaller primary and agglomerate diffusion lengths result in better electrochemical performance. This points us toward designing the morphology of the cathode particles to compensate for the diffusion limitation of LMR-NMC while maximizing the density.

Tewari, Deepti↗

DOE OSTI2023

Using Homogenized Models to Explore the Effect of Fracture Densities on Weathering

Despite its importance, only a few researchers have incorporated the effects of fracturing into models of reactive transport for rock weathering. Here we explore 2D simulations that describe weathering under conditions of diffusive and advective transport within heterogeneous media consisting of rocky blocks and fractures. In our simulations, the Darcy velocities vary in space and time and depend on weathering processes within the rock matrix. We explore simulations with saturated and unsaturated flow for weathering bedrock that consists of blocks separated by inert or weathered material. The simulations show that a simplified homogenized model can approximate exact solutions for some of the simulated columns and hills and can allow exploration of coupling between flow and reaction in fractured rock. These hillslope simulations document that, even in the presence of 2D water flow, i) an increase in fracture density results in faster weathering advance rates; and ii) the water table locates deeper for a rock system that is weathered and fractured rather than weathered and unfractured. Some of these patterns have also been observed for natural systems. But these simulations also highlight how simplified models that do not use appropriate averaging of heterogeneities can be inaccurate in predicting weathering rate for natural systems. For example, if water flows both vertically and laterally through the vadose zone of a hill, then a prediction of the depth of regolith that is based on modeling strictly unidirectional downward infiltration will be unrealistically large. Likewise, if the fracture density observed near the land surface is used in a model to predict depth of weathering for a system where the fracture density decreases downward, the model will overestimate regolith depth. Learning how to develop accurately homogenized models could thus enable better conceptual models and predictions of weathering advance in natural systems.

58 GEOSCIENCES↗

DOE OSTI2021

Use of Physics to Improve Solar Forecast: Physics-Informed Persistence Models for Simultaneously Forecasting GHI, DNI, and DHI

Observation-based statistical models have been widely used in forecasting solar energy; however, existing models often lack a clear relation to physics and are limited largely to global horizontal irradiance (GHI) forecasts over relatively short time horizons (< 1 hour). Incorporating physics into observation-based models, increasing forecast time horizons and developing a model system for forecasting not only GHI but also direct normal irradiance (DNI) and diffuse horizontal irradiance (DHI) remain challenging, especially under cloudy conditions because of complex cloud-radiation interactions. This work attempts to address these challenges by developing a hierarchy of four new physics-informed persistence models that can be used to simultaneously forecast GHI, DNI and DHI. The decade-long measurements (1998 to 2014) at the U.S. Department of Energy's Atmospheric Radiation Measurement (ARM)'s Southern Great Plains (SGP) Central Facility site are used to evaluate the performance of the new models. Overall, the results show that the new physics-informed forecast models generally outperform the simple and smart persistence models, and improve the forecast accuracy at lead times from 1.25 hours up to 6 hours. Further analysis reveals that the forecast error is highly related to the error and temporal variability of the assumed cloud predictor. The best model for forecasting different radiative components can be explained by the relationship between solar irradiances and cloud properties.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Influences of Cloud Microphysics on the Components of Solar Irradiance in the WRF-Solar Model

An accurate forecast of Global Horizontal solar Irradiance (GHI) and Direct Normal Irradiance (DNI) in cloudy conditions remains a major challenge in the solar energy industry. This study focuses on the impact of cloud microphysics on GHI and its partition into DNI and Diffuse Horizontal Irradiance (DHI) using the Weather Research and Forecasting model specifically designed for solar radiation applications (WRF-Solar) and seven microphysical schemes. Three stratocumulus (Sc) and five shallow cumulus (Cu) cases are simulated and evaluated against measurements at the US Department of Energy’s Atmospheric Radiation Measurement (ARM) user facility, Southern Great Plains (SGP) site. Results show that different microphysical schemes lead to spreads in simulated solar irradiance components up to 75% and 350% from their ensemble means in the Cu and Sc cases, respectively. The Cu cases have smaller microphysical sensitivity due to a limited cloud fraction and smaller domain-averaged cloud water mixing ratio compared to Sc cases. Cloud properties also influence the partition of GHI into DNI and DHI, and the model simulates better GHI than DNI and DHI due to a non-physical error compensation between DNI and DHI. The microphysical schemes that produce more accurate liquid water paths and effective radii of cloud droplets have a better overall performance.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2019

Stochastic and coarse-grained two-dimensional modeling of directional particle movement

In this work, we study the evolution of interacting groups of pedestrians in two-dimensional geometries. We introduce a microscopic stochastic model that includes floor fields modeling the global flow of individual groups as well as local interaction rules. From this microscopic model we derive an analytically-tractable system of conservation laws that governs the evolution of the macroscopic pedestrian densities. Numerical simulations show good agreement between the system of conservation laws and the microscopic model, though the latter is slightly more diffusive. We conclude by deriving second-order corrections to the system of conservation laws.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1997

Numerical Study of Stratified Charge Combustion in Wave Rotors

A wave rotor may be used as a pressure-gain combustor effecting non-steady flow, and intermittent, confined combustion to enhance gas turbine engine performance. It will be more compact and probably lighter than an equivalent pressure-exchange wave rotor, yet will have similar thermodynamic and mechanical characteristics. Because the allowable turbine blade temperature limits overall fuel/air ratio to sub-flammable values, premixed stratification techniques are necessary to burn hydrocarbon fuels in small engines with compressor discharge temperature well below autoignition conditions. One-dimensional, unsteady numerical simulations of stratified-charge combustion are performed using an eddy-diffusivity turbulence model and a simple reaction model incorporating a flammability limit temperature. For good combustion efficiency, a stratification strategy is developed which concentrates fuel at the leading and trailing edges of the inlet port. Rotor and exhaust temperature profiles and performance predictions are presented at three representative operating conditions of the engine: full design load, 40% load, and idle. The results indicate that peak local gas temperatures will result in excessive temperatures within the rotor housing unless additional cooling methods are used. The rotor itself will have acceptable temperatures, but the pattern factor presented to the turbine may be of concern, depending on exhaust duct design and duct-rotor interaction.

Nalim, M. Razi↗

DOE OSTI2018

Mechanism-based modeling of solute strengthening: Application to thermal creep in Zr alloy

In this paper, a crystallographic thermal creep model is proposed for Zr alloys that accounts for the hardening contribution of solutes via their time-dependent pinning effect on dislocations. The core-diffusion model proposed by Soare and Curtin (2008a) is coupled with a recently proposed constitutive modeling framework (Wang et al., 2017, 2016) accounting for the heterogeneous distribution of internal stresses within grains. The Coble creep mechanism is also included. This model is, in turn, embedded in the effective medium crystallographic VPSC framework and used to predict creep strain evolution of polycrystals under different temperature and stress conditions. The simulation results reproduce the experimental creep data for Zircaloy-4 and the transition between the low (n~1), intermediate (n~4) and high (n~9) power law creep regimes. This is achieved through the dependence on local aging time of the solute-dislocation binding energy. The anomalies in strain rate sensitivity (SRS) are discussed in terms of core-diffusion effects on dislocation junction strength. The mechanism-based model captures the primary and secondary creep regimes results reported by Kombaiah and Murty (2015a, 2015b) for a comprehensive set of testing conditions covering the 500–600 °C interval, stresses spanning 14–156 MPa, and steady state creep rates varying between 1.5·10 -9 s -1 to 2·10 -3 s -1 . There are two major advantages to this model with respect to more empirical ones used as constitutive laws for describing thermal creep of cladding: 1) specific dependences on the nature of solutes and their concentrations are explicitly accounted for; 2) accident conditions in reactors, such as RIA and LOCA, usually take place in short times, and deformation takes place in the primary, not the steady-state creep stage. Finally, as a consequence, a model that accounts for the evolution with time of microstructure is more reliable for this kind of simulation.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Mixing Cell Model: A One-Dimensional Numerical Model for Assessment of Water Flow and Contaminant Transport in the Unsaturated Zone

This report describes the Mixing Cell Model code, a one dimensional model for water flow and solute transport in the unsaturated zone under steady state or transient flow conditions. The model is based on the principles and assumptions underlying mixing-cell model formulations. The unsaturated zone is discretized into a series of independent mixing cells. Each cell may have unique hydrologic, lithologic, and sorptive properties. Ordinary differential equations describe the material (water and solute) balance within each cell. Water-flow equations are derived from the continuity equation, assuming that unit gradient conditions exist at all times in each cell. Pressure gradients are considered implicitly through model discretization. Unsaturated hydraulic conductivity and moisture contents are determined by the material specific moisture-characteristic curves. Solute-transport processes include explicit treatment of advective processes, first order chain decay, and linear sorption reactions. Dispersion is addressed through implicit and explicit dispersion. Implicit dispersion is an inherent feature of all mixing-cell models and originates from the formulation of the problem in terms of mass balance around fully mixed volume elements. Expressions are provided that relate implicit dispersion to the physical dispersion of the system. Two FORTRAN codes were developed to solve the water flow and solute-transport equations: (1) the Mixing Cell Model for Flow (MCMF) solves transient water-flow problems and (2) the Mixing Cell Model for Transport (MCMT) solves the solute-transport problem. The transient water-flow problem is typically solved first by estimating the water flux through each cell in the model domain as a function of time using the MCMF code. These data are stored in either ASCII or binary files that are later read by the solute transport code (MCMT). Code output includes solute pore water concentrations, water and solute inventories in each cell and at each specified output time, and water and solute fluxes through each cell and specified output time. Computer run times for coupled transient water flow and solute transport were typically several seconds on a 2 GHz Intel Pentium IV desktop computer. The model was benchmarked against analytical solutions and finite element approximations to the partial differential equations (PDE) describing unsaturated flow and transport. Differences between the maximum solute flux estimated by the mixing cell model and the PDE models were typically less than two percent. This revision includes an option for a fixed concentration lower boundary condition for diffusive fluxes for versions 020321 and later.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2016

The Controversial Role of Inter-diffusion in Glass Alteration

Current kinetic models for nuclear waste glasses (e.g. GM2001, GRAAL) are based on a set of mechanisms that have been generally agreed upon within the international waste glass community. These mechanisms are: hydration of the glass, ion exchange reactions (the two processes are referred as inter-diffusion), hydrolysis of the silicate network, and condensation/precipitation of partly or completely hydrolyzed species that produces a porous and amorphous layer and crystalline phases on surface of the altered glass. Recently, a new idea with origins in the mineral dissolution community has been proposed that excludes inter-diffusion process as a potential rate-limiting mechanism. To understand how the so-called interfacial dissolution/precipitation model can change the current understanding of glass behavior, a key experiment used to account for this model was replicated to further revisit the interpretation. This experiment was performed at 50°C, with SON68 glass, in static mode, deionized water and S/V ratio of 10 m-1 for 6 months. It turn out that glass alters in an intermediate kinetic regime between the forward and the residual rate. According to previous and new solid characterizations, it is concluded that neither a simple inter-diffusion model nor the interfacial dissolution precipitation model can account for the observed elemental profiles within the alteration layer. More generally, far and close-to-saturation conditions must be distinguished and literature provides evidences that inter-diffusion takes place in slightly acidic conditions and far from saturation. However, closer to saturation, when a sufficiently dense layer is formed, a new approach is proposed requiring a full description of chemical reactions taking place within the alteration layer and involving water molecules as it is thought that water accessibility to the pristine glass is the rate-limiting process.

Gin, Stephane↗

DOE OSTI2020

Colloid-facilitated transport of 238 Pu, 233 U and 137 Cs through fractured chalk: Laboratory experiments, modelling, and implications for nuclear waste disposal

The influence of montmorillonite colloids on the mobility of 238 Pu, 233 U and 137 Cs through a chalk fracture was investigated to assess the transport potential for radioactive waste. Radioisotopes of each element, along with the conservative tracer tritium, were injected in the presence and absence of montmorillonite colloids into a naturally fractured chalk core. In parallel, batch experiments were conducted to obtain experimental sorption coefficients (K d , mL/g) for both montmorillonite colloids and the chalk fracture material. Breakthrough curves were modelled to determine diffusivity and sorption of each radionuclide to the chalk and the colloids under advective conditions. Uranium sorbed sparingly to chalk (log K d = 0.7 ± 0.2) in batch sorption experiments. 233 U(VI) breakthrough was controlled primarily by the matrix diffusion and sorption to chalk (15 and 25% recovery with and without colloids, respectively). Cesium, in contrast, sorbed strongly to both the montmorillonite colloids and chalk (batch log K d = 3.2 ± 0.01 and 3.9 ± 0.01, respectively). The high affinity to chalk and low colloid concentrations overwhelmed any colloidal Cs transport, resulting in very low 137 Cs breakthrough (1.1–5.5% mass recovery). Batch and fracture transport results, and the associated modelling revealed that Pu migrates both as Pu (IV) sorbed to montmorillonite colloids and as dissolved Pu(V) (7% recovery). Transport experiments revealed differences in Pu(IV) and Pu(V) transport behavior that could not be quantified in simple batch experiments but are critical to effectively predict transport behavior of redox-sensitive radionuclides. Finally, a brackish groundwater solution was injected after completion of the fracture flow experiments and resulted in remobilization and recovery of 2.2% of the total sorbed radionuclides which remained in the core from previous experiments. In general, our study demonstrates consistency in sorption behavior between batch and advective fracture transport. Finally, the results suggest that colloid-facilitated radionuclide transport will enhance radionuclide migration in fractured chalk for those radionuclides with exceedingly high affinity for colloids.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2017

Stochastic modeling of catalytic processes in nanoporous materials: Beyond mean-field approach

Transport and reaction in zeolites and other porous materials, such as mesoporous silica particles, has been a focus of interest in recent years. This is in part due to the possibility of anomalous transport effects (e.g. single-file diffusion) and its impact in the reaction yield in catalytic processes. Computational simulations are often used to study these complex nonequilibrium systems. Computer simulations using Molecular Dynamics (MD) techniques are prohibitive, so instead coarse grained one-dimensional models with the aid of Kinetic Monte Carlo (KMC) simulations are used. Both techniques can be computationally expensive, both time and resource wise. These coarse-grained systems can be exactly described by a set of coupled stochastic master equations, that describe the reaction-diffusion kinetics of the system. The equations can be written exactly, however, coupling between the equations and terms within the equations make it impossible to solve them exactly; approximations must be made. One of the most common methods to obtain approximate solutions is to use Mean Field (MF) theory. MF treatments yield reasonable results at high ratios of reaction rate k to hop rate h of the particles, but fail completely at low k=h due to the over-estimation of fluxes of particles within the pore. We develop a method to estimate fluxes and intrapore diffusivity in simple one- dimensional reaction-diffusion models at high and low k=h, where the pores are coupled to an equilibrated three-dimensional fluid. We thus successfully describe analytically these simple reaction-diffusion one-dimensional systems. Extensions to models considering behavior with long range steric interactions and wider pores require determination of multiple boundary conditions. We give a prescription to estimate the required parameters for these simulations. For one dimensional systems, if single-file diffusion is relaxed, additional parameters to describe particle exchange have to be introduced. We use Langevin Molecular Dynamics (MD) simulations to assess these parameters.

36 MATERIALS SCIENCE↗

DOE OSTI2014

Protection of rubbers against aging with the use of structural, diffusion and kinetic effects

Processes of rubber aging in the different kinetic modes and mathematical models of the processes including aging and destruction processes under extreme conditions have been studied. The aging process for rubbers as thermodynamically open nonlinear systems has also been considered in the paper. It has been revealed that the aging process can be controlled by the internal physical-chemical processes organization and by creation of external influences (by organization of thermodynamic forces and flows). According to the Onsager principle, in certain conditions of aging the conjugation of thermodynamic forces and flows is possible. The diffusion and structural aspects of aging in elastomeric rubbers and articles thereof have been considered. The composite antiaging systems of prolonged action based on the use of microparticles with a saturated solution of the antiagers migrating into an elastomeric matrix at the exploitation have been proposed.

36 MATERIALS SCIENCE↗

NASA NTRS1988

Fatigue life prediction modeling for turbine hot section materials

This paper presents a summary of the life prediction methods developed under the NASA Lewis Research Center's Hot Section Technology (HOST) program. A major objective of the fatigue and fracture efforts under the HOST program was to significantly improve the analytic life prediction tools used by the aeronautical gas turbine engine industry. This has been achieved in the areas of high-temperature thermal and mechanical fatigue of bare and coated high-temperature superalloys. Such technical improvements will eventually reduce life cycle costs. The cyclic crack initiation and propagation resistance of nominally isotropic polycrystalline alloys and highly anisotropic single crystal alloys has been addressed. A sizeable data base has been generated for three alloys [cast PWA 1455 (B-1900 + Hf), wrought Inconel 718, and cast single crystal PWA 1480] in bare and coated conditions. Two coatings systems, diffusion aluminide (PWA 273) and plasma sprayed MCrAIY overlay (PWA 286) were employed. Life prediction modeling efforts were devoted to creep-fatigue interaction, oxidation, coatings interactions, multiaxiality of stress-strain states, mean stress effects, cumulative damage, and thermomechanical fatigue. The fatigue crack initiation life models developed to date include the Cyclic Damage Accumulation (CDA) Model of Pratt & Whitney and the Total Strain Version of Strainrange Partitioning (TS-SRP) of NASA Lewis for nominally isotropic materials, and the Tensile Hysteretic Energy Model of Pratt & Whitney for anisotropic superalloys. The fatigue model being developed by the General Electric Company is based upon the concepts of Path-lndependent Integrals (PII) for describing cyclic crack growth under complex nonlinear response at the crack tip due to thermomechanical loading conditions. A micromechanistic oxidation crack extension model has been derived by researchers at Syracuse University. The models are described and discussed in the paper. Only limited verification has been achieved to-date as several of the technical programs are still in progress and the verification tasks are scheduled, quite naturally, near the conclusion of the program. To-date, efforts have concentrated on development of independent models for cyclic constitutive behavior, cyclic crack initiation, and cyclic crack propagation. The transition between crack initiation and crack propagation has not been thoroughly researched as yet, and the integration of these models into a unified life prediction method has not been addressed.

Metal fatigue↗

DOE OSTI2013

ENERGETIC PARTICLE DIFFUSION IN CRITICALLY BALANCED TURBULENCE

Observations and modeling suggest that the fluctuations in magnetized plasmas exhibit scale-dependent anisotropy, with more energy in the fluctuations perpendicular to the mean magnetic field than in the parallel fluctuations and the anisotropy increasing at smaller scales. The scale dependence of the anisotropy has not been studied in full-orbit simulations of particle transport in turbulent plasmas so far. In this paper, we construct a model of critically balanced turbulence, as suggested by Goldreich and Sridhar, and calculate energetic particle spatial diffusion coefficients using full-orbit simulations. The model uses an enveloped turbulence approach, where each two-dimensional wave mode with wavenumber k is packed into envelopes of length L following the critical balance condition, L{proportional_to}k {sup -2/3} , with the wave mode parameters changing between envelopes. Using full-orbit particle simulations, we find that both the parallel and perpendicular diffusion coefficients increase by a factor of two, compared to previous models with scale-independent anisotropy.

79 ASTRONOMY AND ASTROPHYSICS↗