Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mesoscale modeling”

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 OSTI2024

Resource Assessment for Distributed Wind Energy: An Evaluation of Best-Practice Methods in the Continental US

Current wind resources within the United States (US) indicate a potential to profitably install nearly 1,400 gigawatts of distributed wind (DW) capacity. This amount is equivalent to over half of the United States’ current energy demand from electricity, making it enough to power millions of homes and businesses and replace countless fossil fuel-based generating plants. Despite the potential growth of DW in the US, deployments are presently hindered by a lack of confidence in resource estimation methods. One potential challenge is that smaller-scale turbines, with hub heights of 40 meters or less, are disproportionately impacted by obstacles such as buildings and vegetation. These obstacles may produce complex wake effects, best modeled with high-fidelity complex fluid dynamics (CFD) models that are too computationally expensive to use for routine siting and resource assessment. Thus, installers today make use of heuristics and simple equations to approximate the impact of obstacles while also leveraging long-term resource data from commercial or publicly available atmospheric models. This study evaluates these historical and commonly used methods alongside new lower-order obstacle models produced from CFD simulations and measurement-based bias correction. The preliminary results from this study show the importance of taking care in the choice and application of mesoscale atmospheric models and the significant value of bias correction using measurements from nearby meteorological towers. Detailed obstacle modeling provides only modest additional gains in performance and, in some cases, can add error, especially at sites where turbines have already been located to avoid obvious impact from upwind obstacles. These findings reinforce the importance of collecting in situ measurements and suggest that obstacle models may be better applied in practice to automated or computer-aided siting, rather than in economic wind resource assessments.

17 WIND ENERGY↗
DOE OSTI2022

Evaluation of the Fitch Wind-Farm Wake Parameterization with Large-Eddy Simulations of Wakes Using the Weather Research and Forecasting Model

Abstract Wind-farm parameterizations in weather models can be used to predict both the power output and farm effects on the flow; however, their correctness has not been thoroughly assessed. We evaluate the wind-farm parameterization of the Weather Research and Forecasting Model with large-eddy simulations (LES) of the wake performed with the same model. We study the impact on the velocity and turbulence kinetic energy (TKE) of inflow velocity, roughness, resolution, number of turbines (one or two), and inversion height and strength. We compare the mesoscale with the LES by spatially averaging the LES within areas correspondent to the mesoscale horizontal spacing: one covering the turbine area and two downwind. We find an excellent agreement of the velocity within the turbine area between the two types of simulations. However, within the same area, we find the largest TKE discrepancies because in mesoscale simulations, the turbine-added TKE has to be highest at the turbine position to be advected downwind. Within the downwind areas, differences between velocities increase as the wake recovers faster in the LES, whereas for the TKE both types of simulations show similar levels. From the various configurations, the impact of inversion height and strength is small for these heights and inversion levels. The highest impact for the one-turbine simulations appears under the low-speed case due to the higher thrust, whereas the impact of resolution is low for the large-eddy simulations but high for the mesoscale simulations. Our findings demonstrate that higher-fidelity simulations are needed to validate wind-farm parameterizations.

17 WIND ENERGY↗
DOE OSTI2022

Computing diffraction patterns of microstructures from phase-field simulations

The diffraction pattern of a material contains information not only on the crystal structures of its constituting phases, but also on its mesoscale spatial distributions of phases, grains, and ferroelastic, ferroelectric, and ferromagnetic domains. While diffraction patterns from experiments such as X-ray diffraction are presented in the reciprocal or Fourier space, mesoscale microstructure models such as the phase-field method naturally produce real-space images of spatial distribution of chemical composition, structural, and ferroic domains. Although one could rather readily compute the Fourier amplitudes of chemical and structural domain distributions generated by mesoscale simulations, they only contain information about the length scale and alignment of the real-space chemical and structure domains. Therefore, a direct comparison between diffraction experiments and mesoscale microstructure simulations is not possible. Here, we develop a theoretical approach to directly compute the crystal diffraction patterns of microstructures predicted by phase-field simulations. In particular, we consider five representative examples of microstructure patterns involving purely compositional domains, a single pair of tetragonal twin structures, multiple twin variants in a hexagonal system, ferroelectric polar vortices, and polycrystalline grains. The results are compared with previous experimental observations as well as X-ray diffraction experiments performed in the present study. Furthermore, the theoretical framework allows one to directly connect material microstructures and diffraction patterns predicted from phase-field simulations and the corresponding diffraction patterns from experiments, and thus providing guidance to experimental diffraction characterization and interpretation of microstructures.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Multilayer interface tracking model of pure tungsten oxidation

Here, we present a numerical model to predict oxide scale growth on tungsten surfaces under exposure to oxygen at high temperatures. The model captures the formation of four thermodynamically-compatible oxide sublayers, WO 2 , WO 2.72 , WO 2.9 , and WO 3 , on top of the metal substrate. Oxide layer growth is simulated by tracking the oxide/oxide and oxide/metal interfaces using a sharp-interface Stefan model coupled to diffusion kinetics. The model is parameterized using selected experimental measurements and electronic structure calculations of the diffusivities of all the oxide subphases involved. We simulate oxide growth at temperatures of 600°C and above, extracting the power law growth exponents in each case, which we find to deviate from classical parabolic growth in several cases. We conduct a comparison of the model predictions with an extensive experimental data set, with reasonable agreement at most temperatures. While many gaps in our understanding still exist, this work is a first attempt at embedding the thermodynamic and kinetic complexity of tungsten oxide growth into a comprehensive mesoscale kinetic model that attempts to capture the essential features of tungsten oxidation to fill existing knowledge gaps and guide and enhance future tungsten oxidation models.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Elucidating the Morphological Instability Mechanisms During Electrodeposition of Active Metals

The primary aim of this investigation is to elucidate the mechanisms by which a group of metals pertinent to battery technology exhibit unstable growth morphology during electrodeposition, which are crucial for the operational efficacy of batteries. The metals under consideration include alkali metals alongside zinc, all of which are characterized by high atomic mobility and reactivity with the electrolytes employed in battery systems. These materials represent promising electrode candidates for next-generation rechargeable batteries that are anticipated to deliver significantly enhanced energy density at a reduced cost. Nevertheless, the widespread occurrence of distinctly non-uniform morphologies, exemplified by the emergence of filamentous and moss-like structures, constitutes a significant impediment to their effective utilization, as it leads to accelerated capacity degradation and potential internal short circuits. This study addresses a critical question that has thus far garnered limited attention: Does the residual stress induced by the plating process play an instrumental role in inducing the morphological instability? Utilizing zinc as the primary model system, this project will be executed through three interconnected tasks: (1) the quantification of the relationship between plating stress and electroplating morphologies, (2) the three-dimensional imaging of the microstructure of the mossy structure at the nanoscale, and (3) the development of a mesoscale electroplating model to interpret empirical findings and generate additional predictive insights.

25 ENERGY STORAGE↗
DOE OSTI2022

Effect of Processing on Precipitation Kinetics in Nano Polycrystalline 7075 Aluminum Alloy

Fairmount Technologies (FT) seeks to optimize the thermomechanical processing of Aluminum Alloy 7075. The target yield strength is 20% higher than commercially available 7075-T6 temper sheet which will enable structural light weighting by a similar magnitude and reduce the carbon footprint of the transportation sector. To reduce processing cost and make this commercially viable, it is important to understand physical mechanisms driving microstructural changes. A mesoscale phase-field model seeded with experimentally measured microstructures and assessed thermodynamic and kinetic properties will be developed and exercised using HPC at PNNL. It will account for the effect of deformation, temperature, and chemistry on precipitation kinetics in nanocrystalline microstructures. Our objective is to understand this “stress-aging” coupling by performing parametric simulations with mesoscale models for a better understanding of GP zone and η' formation in nanocrystalline 7075 Al. This will help FT optimize processing schedules to achieve desired microstructures that increase both strength and ductility.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Linking Large-Scale Double-ITCZ Bias to Local-Scale Drizzling Bias in Climate Models

Tropical precipitation in climate models presents significant biases in both the large-scale pattern (i.e., double intertropical convergence zone bias) and local-scale characteristics (i.e., drizzling bias with too frequent drizzle/convection and reduced occurrences of no and heavy precipitation). By untangling the coupled system and analyzing the biases in precipitation, cloud, and radiation, this study shows that local-scale drizzling bias in atmospheric models can lead to large-scale double-ITCZ bias in coupled models by inducing convective-regime-dependent biases in precipitation and cloud radiative effects (CRE). The double-ITCZ bias consists of a hemispherically asymmetric component that arises from the asymmetric SST bias and a nearly symmetric component that exists in atmospheric models without the SST bias. By increasing light rain but reducing heavy rain, local-scale drizzling bias induces positive (negative) precipitation bias in the moderate (strong) convective regime, leading to the nearly symmetric wet bias in atmospheric models. By affecting the cloud profile, local-scale drizzling bias induces positive (negative) CRE bias in the stratocumulus (convective) regime in atmospheric models. Because the stratocumulus (convective) region is climatologically more pronounced in the southern (northern) tropics, the CRE bias is deemed to be hemispherically asymmetric and drives warm and wet (cold and dry) biases in the southern (northern) tropics when coupled to ocean. Our results suggest that correcting local-scale drizzling bias is critical for fixing large-scale double-ITCZ bias. Finally, the drizzling and double-ITCZ biases are not alleviated in models with mesoscale (0.25°–0.5°) or even storm-resolving (~3 km) resolution, implying that either large-eddy simulation or fundamental improvement in small-scale subgrid parameterizations is needed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

First-principles Modeling and Design of Solid-State Interfaces for the Protection and Use of Lithium Metal Anodes

Li-ion batteries are one of the most advanced energy storage technologies in use today. Li-ion batteries are used in a multitude of applications ranging from consumer electronics, medical devices, sensors and grid storage. However, improving the capacity and energy density delivered by current Li-ion technology requires advanced materials research into novel chemical systems. In this project we have focused particularly in the use of solid-state electrolytes with lithium metal electrodes. Research into all solid-state batteries (ASSB) with Li metal electrodes has significantly expanded in recent years, however most studies reported experimental findings, which left substantial room for theoretical and modeling work as a tool to understand and determine design principles allowing reliable and safe use of ASSBs with Li metal. Among the remaining obstacles preventing reliable use of ASSBs with a Li metal electrode, the stability of the interface between the solid electrolyte and Li metal, and the propagation/dendrite formation of Li metal and resulting mechanical degradation of the electrolyte are key phenomenon that are yet to be fully understood. In the current project we have addressed these two coupled phenomena using first principles calculations and mesoscale continuum modeling. We have obtained chemical and electrochemical stability windows for several solid electrolyte materials. Additionally, from mathematical and numerical modeling of Li protrusion and dendrite initiation during plating and stripping we have determined design criteria in terms of chemical, electrochemical, and mechanical properties and operating conditions for which stable deposition can occur. We also considered the effects of mixed electronic-ionic conduction in solid electrolytes, which has more recently been suggested as another important mechanism involved in ASSB failure. Throughout our work we have successfully addressed important questions necessary for the use of ASSB’s. We have determined guiding principles for materials properties and operating conditions necessary to operate ASSB’s. And have proposed novel solid electrolyte materials with predicted chemical stability an ionic conductivity. Although this represents significant progress in our understanding, open questions remain in order to fully develop reliable and safely operate ASSBs with Li metal. Future work, building on this project will require further experimental, theoretical and simulation efforts to address remaining questions.

25 ENERGY STORAGE↗
DOE OSTI2022

Scientific challenges to characterizing the wind resource in the marine atmospheric boundary layer

Abstract. With the increasing level of offshore wind energy investment, it is correspondingly important to be able to accurately characterize the wind resource in terms of energy potential as well as operating conditions affecting wind plant performance, maintenance, and lifespan. Accurate resource assessment at a particular site supports investment decisions. Following construction, accurate wind forecasts are needed to support efficient power markets and integration of wind power with the electrical grid. To optimize the design of wind turbines, it is necessary to accurately describe the environmental characteristics, such as precipitation and waves, that erode turbine surfaces and generate structural loads as a complicated response to the combined impact of shear, atmospheric turbulence, and wave stresses. Despite recent considerable progress both in improvements to numerical weather prediction models and in coupling these models to turbulent flows within wind plants, major challenges remain, especially in the offshore environment. Accurately simulating the interactions among winds, waves, wakes, and their structural interactions with offshore wind turbines requires accounting for spatial (and associated temporal) scales from O(1 m) to O(100 km). Computing capabilities for the foreseeable future will not be able to resolve all of these scales simultaneously, necessitating continuing improvement in subgrid-scale parameterizations within highly nonlinear models. In addition, observations to constrain and validate these models, especially in the rotor-swept area of turbines over the ocean, remains largely absent. Thus, gaining sufficient understanding of the physics of atmospheric flow within and around wind plants remains one of the grand challenges of wind energy, particularly in the offshore environment. This paper provides a review of prominent scientific challenges to characterizing the offshore wind resource using as examples phenomena that occur in the rapidly developing wind energy areas off the United States. Such phenomena include horizontal temperature gradients that lead to strong vertical stratification; consequent features such as low-level jets and internal boundary layers; highly nonstationary conditions, which occur with both extratropical storms (e.g., nor'easters) and tropical storms; air–sea interaction, including deformation of conventional wind profiles by the wave boundary layer; and precipitation with its contributions to leading-edge erosion of wind turbine blades. The paper also describes the current state of modeling and observations in the marine atmospheric boundary layer and provides specific recommendations for filling key current knowledge gaps.

17 WIND ENERGY↗
DOE OSTI2021

From classical thermodynamics to phase-field method

Phase-field method is a density-based computational method at the mesoscale for modeling and predicting the temporal microstructure and property evolution during materials processes. The focus of this article is on connecting the most common phase-field equations to the very basic first and second laws of classical thermodynamics through rudimentary irreversible thermodynamics. It briefly discusses the relations of the continuum phase-field equations to their counter parts at the microscopic and atomic levels. It attempts to clarify the contributions of long-range elastic, electrostatic, and magnetic interactions to domain structure evolution during structural, ferroelectric, and ferromagnetic phase transformations by separating order parameter changes due to the presence of quasi-static fields and those arising from phase transformations. A few examples are presented to demonstrate the possibility of employing the phase-field method to provide guidance to designing materials for optimum properties or discovering novel mesoscale phenomena or new materials functionalities. Here, the article ends with a brief perspective on a number of potential future directions on the development and applications of phase-field method beyond its traditional applications to structural alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2022

TRITIUM TRANSPORT MODELING AT THE PORE SCALE IN CERAMIC BEEDER MATERIALS USING TMAP8

Fusion reactors depend on the blanket material to breed and release tritium at the same rate or faster than it is consumed by the fusion reaction. Cellular ceramic breeders are dense materials that can maintain a high tritium breeding ratio while promoting tritium release thanks to highly connected pores. Assessing the tritium breeding capabilities of these materials requires a combination of extensive experimental and modeling efforts. In this work, we develop and calibrate a multiphysics model of tritium transport. This novel model accounts for ceramic and pore diffusion, trapping and detrapping, and several surface reactions at the pore surface. We perform a sensitivity analysis and calibrate the model by comparing its predictions against experimental measurements of deuterium absorption. The calibrated model is then used to model tritium absorption in samples with different pore microstructures to investigate the effect of pore interconnectivity on tritium absorption. The model is part of the development of the multiscale, multiphysics framework for tritium transport (i.e., the Tritium Migration Analysis Program [TMAP8]), which is itself built on top of the finite element multiphysics framework MOOSE. This study demonstrates some of TMAP8’s capabilities, and is the first step toward assessing the tritium breeding capabilities of ceramic breeder material designs.

08 HYDROGEN↗
DOE OSTI2026

Assessing the Impact of a Forest Canopy on Near-Surface Wind Statistics

Representing the forest canopy in atmospheric numerical models should improve simulated winds within and above the canopy up to a few hundred meters above the ground. Here, in this study, we implement a forest canopy parameterization into the Weather Research and Forecasting (WRF) Model in a large-eddy simulation (LES) mode by applying drag forces across multiple layers within the canopy height. We use unique observations from the Lidar Experiments for Assessing Flow over Forests (LEAFF) field campaign at the Wind River Experimental Forest (WREF) in the U.S. Pacific Northwest to evaluate model performance. In a 2-day case study, the canopy parameterization improved wind predictions both within and above the canopy, particularly during the daytime and at finer grid resolution. Without it, winds were frequently overpredicted above the canopy. Similarly, derived quantities such as the wind shear index also yielded estimates closer to observations with the canopy parameterization implemented. These findings suggest that representing the canopy using drag forces alone can improve simulated mean winds up to 200 m above the surface. Furthermore, second-order statistical moments of wind were more sensitive to canopy density than first-order moments, especially during the daytime. This increased sensitivity and the improved daytime performance in wind speed—evidenced by the lowest bias from observations (3% compared to 20% over diurnal cycle)—imply that winds above the canopy layer are strongly influenced by how well turbulence above the canopy is modeled. The results of this study can serve as a foundation for parameterizing forest canopy effects in coarser weather forecast models.

Energy - Wind↗
DOE OSTI2021

Evolution of Idealized Vortices in Monsoon-Like Shears: Application to Monsoon Depressions

Abstract This study examines processes fundamental to the development of South Asian monsoon depressions using an array of integrations of an idealized convection-permitting numerical model. In each integration, a wave of initially small amplitude is subjected to a different amount of vertical and meridional wind shear, with temperature and moisture fields constructed according to realistic constraints. Based on the evolution of this disturbance into monsoon depression–like vortices, two features of the background environment emerge as important: the low-level gradient of moist static energy (MSE) and the low-level meridional shear. As the low-level MSE gradient steepens, the disturbance becomes stronger and produces more rain. This strengthening results from the interaction of the vortex with latent heat release by convection that is in turn organized by positive MSE advection in the northerly flow west of the vortex. In this region of advection, moister air from the north ascends along upward-sloping isentropes, driving moist convection. The disturbance also becomes stronger with increasing meridional shear, which makes the environment more barotropically unstable. The absence of either of these two features of the background environment prevents substantial growth of the disturbance. Our results suggest that monsoon depression growth in South Asia is fostered by the coexistence of a strong low-level MSE gradient with strong meridional wind shear associated with the monsoon trough.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Lake Huron Enhances Snowfall Downwind of Lake Erie: A Modeling Study of the 2010 New Year’s Lake Effect Snowfall Event

In the Northeast and Great Lakes regions of the United States, the influence of multiple lakes on the overlying air can greatly affect lake-effect snowfall over downwind communities. To assess the impact of Lake Huron on snowfall downwind of Lake Erie, we simulated a lake-effect snow event which occurred from 1 to 6 January 2010 using the Regional Atmospheric Modeling System (RAMS). We found that the presence of Lake Huron enhances snowfall downwind of Lake Erie by almost 20% and leads to much heavier local snowfall totals than when Lake Huron is not present. This increase in snowfall is due to a lake-to-lake (L2L) convective band, as secondary circulations associated with lake-effect convection form over Lake Huron and persist overland between the lakes before reintensifying over Lake Erie. As these secondary circulations move over Lake Erie, lifting within the upward branch of the circulations accelerates the development of convection within the L2L band. Surface fluxes and convection over Lake Huron deepen the boundary layer, promoting deeper vertical development of this L2L band over Lake Erie. However, although this boundary layer modification strengthens the L2L band, we found that it actually reduced lake-effect snowfall produced by wind-parallel bands (WPB) in other parts of the lake. This indicates that boundary layer modification from upstream lakes may impact L2L bands differently than WPB which are not connected to upstream secondary circulations.

Lake effects↗
DOE OSTI2022

Assimilation of a Coordinated Fleet of Uncrewed Aircraft System Observations in Complex Terrain: Observing System Experiments

Uncrewed aircraft system (UAS) observations from the Lower Atmospheric Profiling Studies at Elevation–A Remotely-Piloted Aircraft Team Experiment (LAPSE-RATE) field campaign were assimilated into a high-resolution configuration of the Weather Research and Forecasting (WRF) Model. The impact of assimilating targeted UAS observations in addition to surface observations was compared to that obtained when assimilating surface observations alone using observing system experiments (OSEs) for a terrain-driven flow case and a convection initiation (CI) case observed within Colorado’s San Luis Valley (SLV). The assimilation of UAS observations in addition to surface observations results in a clear increase in skill for both flow regimes over that obtained when assimilating surface observations alone. For the terrain-driven flow case, the UAS observations improved the representation of thermal stratification across the northern SLV, which produced stronger upvalley flow over the eastern half of the SLV that better matched the observations. For the CI case, the UAS observations improved the representation of the pre-convective environment by reducing dry biases across the SLV and over the surrounding terrain. This led to earlier CI and more organized convection over the foothills that spilled outflows into the SLV, ultimately helping to increase low-level convergence and CI there. In addition, the importance of UAS capturing an outflow that originated over the Sangre de Cristo Mountains and triggered CI is discussed. These outflows and subsequent CI were not well captured in the simulation that assimilated surface observations alone. We report that observations obtained with a fleet of UAS are shown to notably improve high-resolution analyses and short-term predictions of two very different mesogamma-scale weather events.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Influence of the distribution of expansive sites in aggregates on microscopic damage caused by alkali-silica reaction: Insights into the mechanical origin of expansion

The origin of damage in concrete due to the alkali-silica reaction (ASR) is attributed to the expansion site in the aggregate. To investigate the cracking process of the aggregate during ASR and its consequences on concrete damage, the effect of the distribution of the expansive sites in the aggregate on ASR expansion and the crack patterns must be evaluated. Thus, in this study, a mesoscale discrete model was applied to ASR modeling to represent the propagation of cracking during ASR accurately. The distribution of the expansive sites in the aggregate was based on the gel pocket and reaction rim models, which are two ASR mechanisms reported in the literature. These two expansion models highlight the different crack patterns obtained based on the aggregate characteristics. Further, the expansion cracking processes determined based on the gel pocket and reaction rim models are consistent with the evolution of cracking with the expansion level.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Sensitivity of the strength and toughness of concrete to the properties of the interfacial transition zone

We report Civil infrastructure is made primarily of concrete structures or components and therefore understanding durability and fracture behavior of concrete is of utmost importance. Concrete contains an interfacial transition zone (ITZ), a porous region surrounding the aggregates, that is often considered to be the weakest region in the concrete. The ITZ is poorly characterized and property estimates for the ITZ differ considerably. In this simulation study, representative concrete mesostructures are produced by packing coarse aggregates with realistic geometries into a mortar matrix. A meshless numerical method, peridynamics, is utilized to simulate the mechanical response including fracture under uniaxial compression and tension. The sensitivity of the stiffness and fracture toughness of the samples to the ITZ properties is computed, showing strong relationships between the ITZ properties and the effective modulus and effective yield strength of the concrete. These results provides insight into the influence of the poorly characterized ITZ on the stiffness and strength of concrete. This work showcases the applicability of peridynamics to concrete systems, matching experimental strength and modulus values. Additionally, relationships between the ITZ’s mechanical properties and the overall concrete strength and stiffness are presented to enable future design decisions.

36 MATERIALS SCIENCE↗
DOE OSTI2021

On the computational solution of vector-density based continuum dislocation dynamics models: A comparison of two plastic distortion and stress update algorithms

Continuum dislocation dynamics models of mesoscale plasticity consist of dislocation transport-reaction equations coupled with crystal mechanics equations. The coupling between these two sets of equations is such that dislocation transport gives rise to the evolution of plastic distortion (strain), while the evolution of the latter fixes the stress from which the dislocation velocity field is found via a mobility law. Earlier solutions of these equations employed a staggered solution scheme for the two sets of equations in which the plastic distortion was updated via time integration of its rate, as found from Orowan's law. In this work, we show that such a direct time integration scheme can suffer from accumulation of numerical errors. We introduce an alternative scheme based on field dislocation mechanics that ensures consistency between the plastic distortion and the dislocation content in the crystal. The new scheme is based on calculating the compatible and incompatible parts of the plastic distortion separately, and the incompatible part is calculated from the current dislocation density field. Stress field and dislocation transport calculations were implemented within a finite element based discretization of the governing equations, with the crystal mechanics part solved by a conventional Galerkin method and the dislocation transport equations by the least squares method. A simple test was first performed to show the accuracy of the two schemes for updating the plastic distortion, which shows that the solution method based on field dislocation mechanics is more accurate. This method then was used to simulate an austenitic steel crystal under uniaxial loading and multiple slip conditions. By considering dislocation interactions caused by junctions, a hardening rate similar to discrete dislocation dynamics simulation results was obtained. Finally, the simulations show that dislocations exhibit some self-organized structures as the strain is increased.

36 MATERIALS SCIENCE↗