Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spectral 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 649 records · Page 36

Machine learning inversion of interatomic force constants from single-crystal inelastic neutron scattering

Atomic vibrations govern many macroscopic properties of materials, but experiments to comprehensively probe them remain challenging. Inelastic neutron scattering (INS) is a powerful technique to map phonon dispersions in crystals, especially when leveraging modern time-of-flight (ToF) spectrometers with large detectors. However, efficiently and robustly extracting interatomic force constants (FCs) parameterizing phonon dynamics from experimental spectra remains a bottleneck due to the complexity and high dimensionality of ToF INS datasets. Here, we present a machine learning approach for the direct inversion of FCs from single-crystal INS measurements. The framework leverages synthetic training data generated using universal machine-learned force fields and an efficient physics-based forward model. We benchmark two neural architectures–one emphasizing structured latent representation learning and the other direct, supervised spectral regression–across simulated datasets for two materials under idealized and noisy conditions. The latent-representation model is subsequently applied to experimental single-crystal INS data on germanium. The model is shown to reproduce FCs derived from both first-principles simulations and from iterative optimization, and furthermore achieves reliable inference even from sparse, single-orientation measurements representing short data acquisitions. Analysis of the learned latent space reveals semantically continuous and physically interpretable encodings that support strong cross-domain generalization. By bridging theoretical and experimental domains, we establish a path toward rapid inversion of experimental spectra and data-driven interpretation of temperature-dependent lattice dynamics.

42 ENGINEERING↗

Freestream turbulence effects on unsteady wind turbine loads and wakes: An IDDES study

We investigate numerically the effects of freestream turbulence on the unsteady aerodynamics and wakes of the National Renewable Energy Laboratory Phase VI wind turbine rotor for increasing wind speed. Turbulence is modeled using the Improved Delayed Detached-Eddy Simulation (IDDES) method. As a first step, a detailed mesh resolution study is conducted with the decaying freestream turbulence model at turbulence intensity of 0.5%. Our blade-resolved IDDES simulations show that grid-independent average torque and thrust results can be achieved with relatively coarse meshes, whereas dramatically higher mesh resolution is required for grid-independent results for power spectral densities of thrust force, especially in the deep-stall regime. Comparing the loads with the Shear-Stress Transport model demonstrates the superiority of IDDES in predicting massively separated flows. The aerodynamic performance and wake predictions with the decaying freestream turbulence model are compared with the synthetic freestream turbulence model. Both models predict nearly the same loads, spectral energy content, and wake characteristics. The properties of both the near- and far-wake regions are then examined. Furthermore, we show that separated boundary layers accelerate turbulent mixing and entrainment of the external flow, which results in faster wake recovery. The effect of increasing turbulence intensity to 6% is investigated using the synthetic freestream turbulence model. In contrast with the fully attached boundary layer, higher freestream turbulence in deep stall does not significantly affect the loads and vortex-shedding characteristics. However, the turbulent mixing in the wake is enhanced, which further hastens the recovery of the self-similar velocity profile. In general, increasing the wind speed at high turbulence intensity shifts the recovery farther upstream and increases the wake width.

17 WIND ENERGY↗

Accurate Prediction of Algal Biomass Lipid, Protein, and Carbohydrate Composition with Machine Learning Regression Modelling of Near-IR Spectra

During large scale algal biomass cultivation, it is difficult to reliably control relative composition to target levels. Rapid determination of chemical composition is feasible by using near infrared (NIR) spectral data. We sought to build and improve on reliable high-throughput screening prediction method based on partial least squares regression (PLSR) by the application of artificial neural networks (ANN) and associated optimization strategies. The algal biomass sample set was designed and created in an iterative process of culturing in physiologically diverse conditions at the GAI field site, followed by compositional analyses at NREL. The workflow allowed us to identify gaps in compositional space for informing the subsequent cultivation and sampling efforts and generated a high quality set of 210 unique samples with chemical analysis results, spectral scanning data, and cultivation metadata. We observed a significant improvement in the performance of carbohydrate content predictions using an optimized ANN model compared to PLSR, with > 16% reduction in mean absolute percent error (MAPE) when tested on the same set of reserved data. The optimized ANN models for FAME and protein prediction performed exceptionally well with 5.99% and 5.09% MAPE, respectively. Application of these methods to detection and quantification of minor biomass constituents that are relevant to certain product streams has shown positive preliminary results, opening the possibility for extensions to the outputs of this powerful data type. All models are accompanied by prediction uncertainties and unsupervised spectral outlier detection to alert an operator to unreliable spectral data. These tools can be deployed for rapid determination of algal culture status, and cultivation and biomass quality improvement.

algal biofuels↗

Active- and transfer-learning applied to microscale-macroscale coupling to simulate viscoelastic flows

Active- and transfer-learning are applied to microscale dynamics of polymer flows for the multiscale discovery of effective constitutive approximations required in viscoelastic flow simulation. The result is macroscopic rheology directly connected to a microstructural model. Micro and macroscale simulations are adaptively coupled by means of Gaussian process regression (GPR) to run the expensive microscale computations only as necessary. This multiscale method is demonstrated with flows of a polymer solution as a model system. At the microscale level dissipative particle dynamics (DPD) is employed to model the fluid as a suspension of bead-spring micro-structures subjected to steady shear flow. The results yield the non-Newtonian viscosity and the first normal stress difference at strain rates as training data used in a GPR model. DPD parameters are calibrated with respect to experimental data for a real polymer solution. Compliance with these data requires adjustment of the DPD model's cutoff radius, which then becomes a function of the second invariant of the strain rate tensor. The FENE-P model is chosen for the macroscale description using the spectral element method (SEM) to simulate channel flow and flow past a circular cylinder. The DPD results at the lowest possible shear strain rate yield an estimate of the zero-shear rate viscosity, which allows the initiation of the macroscale flow by SEM as a Newtonian fluid. The resulting strain-rate field is surveyed to determine additional shear strain rate sampling points for the DPD system. This new information allows an initial fitting of parameters of the constitutive equation followed by new SEM simulations at the macroscale. Additionally, guided by active-learning GPR to select new sampling points, this process continues until convergence is achieved. The effectiveness of this new simulation paradigm for viscoelastic flows is tested with different macroscale operating conditions. The effective closure learned in the channel simulation is then transferred directly to the flow past a circular cylinder at low Reynolds number, where the results show that only two additional DPD simulations are required to achieve a satisfactory constitutive model. With an increase of the Reynolds number, the active-learning scheme automatically detects the inaccuracy of the learned constitutive model, and initiates additional DPD simulations for the extra data needed to once again close the microscale-macroscale coupled system. This new paradigm of active- and transfer-learning for multiscale modeling is readily applicable to other microscale-macroscale coupled simulations of complex fluids and other materials. Furthermore, the coupling between microscale and macroscale solvers can be seamlessly implemented with our open source multiscale universal interface (MUI) library.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Experimental methods to test photoionized plasma models in emission

The expanding foil photoionized plasma platform was developed on the Z-machine at Sandia National Laboratories to produce terrestrial photoionized plasmas at some of the same conditions as those observed in astrophysical black hole accretion disks. The specific scientific goals of the platform are motivated by long-standing astrophysical puzzles related to accreting black hole systems. Ultimately, the experimental goal is to inform the supersolar Fe abundance problem by testing photoionized plasma models in emission with at-parameter laboratory data. High quality high resolution emission data with absolute intensity calibration from a laboratory photoionized plasma have never been collected prior to the experiments described here, providing broader motivation. A robust test of the model predictions requires the successful production of a laboratory photoionized plasma along with measurements of the supporting model inputs including the radiation drive spectrum used to create the photoionized plasma, the electron density, and the temperature. The measured absolute spectral radiance in emission collected from the independently diagnosed plasma can be compared against model calculations performed using the measured inputs. In this paper, we focus on the experimental platform and the model inputs, specifically detailing the methodology used to measure the plasma conditions. The moderate electron temperature of ∼41 eV ± 15 eV and electron density of ∼1e19 ± 1.6e18 e−/cm3 and the ∼3e12 W/cm2 x-ray irradiance at the sample confirm that photoionization dominates the plasma ionization and spectra.

Cho, P. B. (ORCID:0000000291632905)↗

Recent advances and applications of deep learning methods in materials science

Deep learning (DL) is one of the fastest-growing topics in materials data science, with rapidly emerging applications spanning atomistic, image-based, spectral, and textual data modalities. DL allows analysis of unstructured data and automated identification of features. The recent development of large materials databases has fueled the application of DL methods in atomistic prediction in particular. In contrast, advances in image and spectral data have largely leveraged synthetic data enabled by high-quality forward models as well as by generative unsupervised DL methods. In this article, we present a high-level overview of deep learning methods followed by a detailed discussion of recent developments of deep learning in atomistic simulation, materials imaging, spectral analysis, and natural language processing. For each modality we discuss applications involving both theoretical and experimental data, typical modeling approaches with their strengths and limitations, and relevant publicly available software and datasets. We conclude the review with a discussion of recent cross-cutting work related to uncertainty quantification in this field and a brief perspective on limitations, challenges, and potential growth areas for DL methods in materials science.

36 MATERIALS SCIENCE↗

Comparison of Numerical Methods for Modeling the Wave Field Effects Generated by Individual Wave Energy Converters and Multiple Converter Wave Farms

This numerical study compares the wave field generated by the spectral wave action balance code, SNL-SWAN, to the linear-wave boundary-element method (BEM) code, WAMIT. The objective of this study is to assess the performance of SNL-SWAN for modeling wave field effects produced by individual wave energy converters (WECs) and wave farms comprising multiple WECs by comparing results from SNL-SWAN with those produced by the BEM code WAMIT. BEM codes better model the physics of wave-body interactions and thus simulate a more accurate near-field wave field than spectral codes. In SNL-SWAN, the wave field’s energy extraction is modeled parametrically based on the WEC’s power curve. The comparison between SNL-SWAN and WAMIT is made over a range of incident wave conditions, including short-, medium-, and long-wavelength waves with various amounts of directional spreading, and for three WEC archetypes: a point absorber (PA), a pitching flap (PF) terminator, and a hinged raft (HR) attenuator. Individual WECs and wave farms of five WECs in various configuration were studied with qualitative comparisons made of wave height and spectra at specific locations, and quantitative comparisons of the wave fields over circular arcs around the WECs as a function of radial distance. Results from this numerical study demonstrate that in the near-field, the difference between SNL-SWAN and WAMIT is relatively large (between 20% and 50%), but in the far-field from the array the differences are minimal (between 1% and 5%). The resultant wave field generated by the two different numerical approaches is highly dependent on parameters such as: directional wave spreading, wave reflection or scattering, and the WEC’s power curve.

16 TIDAL AND WAVE POWER↗

Theorem on extensive spectral degeneracy for systems with rigid higher symmetries in general dimensions

We establish, in the spirit of the Lieb-Schultz-Mattis theorem, lower bounds on the spectral degeneracy of quantum systems with higher (gaugelike) symmetries with rather generic physical boundary conditions in an arbitrary number of spatial dimensions. Contrary to applying twists or equivalent adiabatic operations, we exploit the effects of modified boundary conditions. We find that when a general choice of boundary geometry is immaterial in determining spectral degeneracies while approaching the thermodynamic limit, systems that exhibit rigid noncommuting gaugelike symmetries, such as the orbital compass model, must have an exponential (in the size of the boundary) degeneracy of each of their spectral levels.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Neural Network Approaches for Mobile Spectroscopic Gamma-Ray Source Detection

Artificial neural networks (ANNs) for performing spectroscopic gamma-ray source identification have been previously introduced, primarily for applications in controlled laboratory settings. To understand the utility of these methods in scenarios and environments more relevant to nuclear safety and security, this work examines the use of ANNs for mobile detection, which involves highly variable gamma-ray background, low signal-to-noise ratio measurements, and low false alarm rates. Simulated data from a 2” × 4” × 16” NaI(Tl) detector are used in this work for demonstrating these concepts, and the minimum detectable activity (MDA) is used as a performance metric in assessing model performance.In addition to examining simultaneous detection and identification, binary spectral anomaly detection using autoencoders is introduced in this work, and benchmarked using detection methods based on Non-negative Matrix Factorization (NMF) and Principal Component Analysis (PCA). On average, the autoencoder provides a 12% and 23% improvement over NMF- and PCA-based detection methods, respectively. Additionally, source identification using ANNs is extended to leverage temporal dynamics by means of recurrent neural networks, and these time-dependent models outperform their time-independent counterparts by 17% for the analysis examined here. The paper concludes with a discussion on tradeoffs between the ANN-based approaches and the benchmark methods examined here.

Bilton, Kyle J. (ORCID:0000000184553689)↗

Fully compressible simulations of waves and core convection in main-sequence stars

Context: Recent, nonlinear simulations of wave generation and propagation in full-star models have been carried out in the anelastic approximation using spectral methods. Although it makes long time steps possible, this approach excludes the physics of sound waves completely and requires rather high artificial viscosity and thermal diffusivity for numerical stability. A direct comparison with observations is thus limited. Aims: In this study, we explore the capabilities of our compressible multidimensional Seven-League Hydro (SLH) code to simulate stellar oscillations. Methods: We compare some fundamental properties of internal gravity and pressure waves in 2D SLH simulations to linear wave theory using two test cases: (1) an interval gravity wave packet in the Boussinesq limit and (2) a realistic 3 M ⊙ stellar model with a convective core and a radiative envelope. Oscillation properties of the stellar model are also discussed in the context of observations. Results: Our tests show that specialized low-Mach techniques are necessary when simulating oscillations in stellar interiors. Basic properties of internal gravity and pressure waves in our simulations are in good agreement with linear wave theory. As compared to anelastic simulations of the same stellar model, we can follow internal gravity waves of much lower frequencies. The temporal frequency spectra of velocity and temperature are flat and compatible with the observed spectra of massive stars. Conclusion: The low-Mach compressible approach to hydrodynamical simulations of stellar oscillations is promising. Our simulations are less dissipative and require less luminosity boosting than comparable spectral simulations. The fully-compressible approach allows for the coupling of gravity and pressure waves in the outer convective envelopes of evolved stars to be studied in the future.

79 ASTRONOMY AND ASTROPHYSICS↗

Fast and Thermal Neutron Spectrum Dosimetry Measurements in the Advanced Test Reactor large-B and small-I Positions Following the Sixth Core Internals Change-out [Slides]

The Advanced Test Reactor (ATR) has a wide variety of irradiation positions that have had experiments that were developed by users from around the world. Most experiment irradiations rely on thoroughly benchmarked numerical models. However, some irradiation positions in ATR are not as well- characterized and are complicated by spectral perturbations from control cylinder orientation. The “small I” and “large-B” irradiation positions are located nearby control cylinders and suffer from these flux perturbations from control cylinder orientations that change during an irradiation cycle to maintain the desired core power distribution. Models predicted that these types of position would exhibit both spectral shifts and amplitude changes in neutron flux, but few measurements have been conducted to benchmark these predictions. Recently, requalification testing was performed to confirm the operational readiness of the ATR following the completion of the Core Internals Change-out (CIC). These tests provided a unique opportunity to validate the analytical methods used to simulate the ATR because nearly all components in the reactor were in a clean as-built state, significantly reducing modelling uncertainties. One subset of the testing included characterization of the fast and thermal neutron flux in the “small-I” and “large-B” positions using silver, cobalt, and nickel neutron dosimetry. In contrast to typical irradiation cycles, the control cylinders were held in position during the post-CIC nuclear testing. The specific activity of these dosimeter wires was measured following two low-power tests with different control cylinder positions.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

A317_Fast and Thermal Neutron Spectrum Dosimetry Measurements in the Advanced Test Reactor large-B and small-I Positions Following the Sixth Core Internals Change-out

The Advanced Test Reactor (ATR) has a wide variety of irradiation positions that have had experiments that were developed by users from around the world. Most experiment irradiations rely on thoroughly benchmarked numerical models. However, some irradiation positions in ATR are not as well- characterized and are complicated by spectral perturbations from control cylinder orientation. The “small I” and “large-B” irradiation positions are located nearby control cylinders and suffer from these flux perturbations from control cylinder orientations that change during an irradiation cycle to maintain the desired core power distribution. Models predicted that these types of position would exhibit both spectral shifts and amplitude changes in neutron flux, but few measurements have been conducted to benchmark these predictions. Recently, requalification testing was performed to confirm the operational readiness of the ATR following the completion of the Core Internals Change-out (CIC). These tests provided a unique opportunity to validate the analytical methods used to simulate the ATR because nearly all components in the reactor were in a clean as-built state, significantly reducing modelling uncertainties. One subset of the testing included characterization of the fast and thermal neutron flux in the “small-I” and “large-B” positions using silver, cobalt, and nickel neutron dosimetry. In contrast to typical irradiation cycles, the control cylinders were held in position during the post-CIC nuclear testing. The specific activity of these dosimeter wires was measured following two low-power tests with different control cylinder positions.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

A structured framework for predicting sustainable aviation fuel properties using liquid-phase FTIR and machine learning

Sustainable aviation fuels have the potential to improve efficiency, reduce emissions, and enhance energy security. To help identify viable sustainable aviation fuels and accelerate research, machine learning models have been developed to predict relevant physicochemical properties. However, many models have limited applicability, leverage data from complex analytical techniques with confined spectral ranges, or use feature decomposition methods that offer limited interpretability. Using liquid-phase Fourier Transform Infrared (FTIR) spectra, this study presents a structured method for creating accurate and interpretable property prediction models for neat molecules, aviation fuels, and blends. Liquid FTIR spectra can be collected quickly and consistently, offering high reliability, sensitivity, and component specificity using less than 2 ml of sample. The method first decomposes FTIR spectra into fundamental building blocks using non-negative matrix factorization (NMF) to enable scientific analysis of FTIR spectra attributes and fuel properties. The NMF features are then used to create five ensemble models for predicting final boiling point, flash point, freezing point, density at 15°C, and kinematic viscosity at -20°C. All models were trained using experimental property data from neat molecules, aviation fuels, and blends. The models accurately predict key properties across a broad range of neat molecules and representative fuels and blends, while enabling interpretation of relationships between compositional elements, such as functional groups or chemical classes, and their resulting properties. This demonstrates strong potential to support sustainable aviation fuel research and development. The models and data are available on an interactive web tool.

Fourier transform infrared spectroscopy↗

Frequency-domain modeling of floating wind arrays with shared mooring lines

A frequency-domain model for floating wind turbine dynamics has been extended to model floating wind farms with couplings from shared mooring systems. The model, called RAFT, could previously calculate the mean offsets and wave-induced response spectra for single floating wind turbines. Now, the model supports multiple floating wind turbines, each with their own properties and responses, along with mooring lines that run directly between floating wind turbines in the array, meaning that shared mooring lines or fully suspended dynamic power cables can be included. This capability is achieved by setting up an array-level solution of the system mean offsets and assembling the full system matrices for solving the dynamic response. The quasi-static mooring model MoorPy is used to linearize the mooring system properties. To compute the floating wind turbine relative motions, phase offsets are applied to each turbine's response as a function of wave frequency based on the wavelength and relative positions in the array. These differential motions are then applied to mooring system tension Jacobians to compute the tension loads in the shared mooring lines. Overall, the capability provides a frequency-domain analog to the modeling capabilities of the floating support structure in FAST.Farm. Mean offsets and power spectral density plots of responses are compared between RAFT and FAST.Farm to verify the implementation. The results indicate good agreement within the expectations of a frequency-domain modeling approach and suggest correct implementation of the shared mooring aspects. Additionally, a unique comb-like frequency response in the shared mooring line tensions has been observed. This phenomenon has a clear physical basis and may be an important design consideration for future shared mooring systems.

17 WIND ENERGY↗

In Situ Investigation of Upper Airway Occlusion in Sleep Disordered Breathing Using Ultrasonic Transducer Arrays

This work presents a novel application of ultrasound for the real-time, non-invasive investigation of occlusion of the upper airway during events of obstructive sleep apnea/hypopnea syndrome. It is hypothesized that ultrasonic pulses applied to the neck during apneic events produce spectral and temporal features that can detect apnea occurrence. Theoretical models of ultrasound propagation and an in vitro test were conducted to test this hypothesis in both transmission and reflection modes. Complete specifications and technical details of the system design and fabrication, which is mounted on each subject’s neck, are presented, including the methodology. Nine patients (seven male and two female, mean age of 42 years, with a range of 25 to 56 years, and body mass index 37.6 ± 6.6 kg/m2) were recruited for a full night study, which included simultaneous nocturnal polysomnography for the validation of the results. Nine temporal features and four spectral features were extracted from the envelope of the received pulse waveform. These were used to compute 26 metrics to quantify the changes in the ultrasonic waveforms between normal breathing and apneic events. The statistical analysis of the collected ultrasonic data showed that at least two or more of the proposed features could detect apneic events in all subjects. The findings establish the feasibility of the proposed method as a cost-effective and non-invasive OSAHS screening tool.

60 APPLIED LIFE SCIENCES↗

Hard X-Ray to Radio Multiwavelength SED Analysis of Local U/LIRGs in the GOALS Sample with a Self-consistent AGN Model including a Polar-dust Component

We conduct hard X-ray to radio multiwavelength spectral energy distribution (SED) decomposition for 57 local luminous and ultraluminous infrared galaxies observed with the Nuclear Spectroscopic Telescope Array and/or Swift/Burst Alert Telescope in the GOALS sample. We modify the latest SED-fitting code X-CIGALE by implementing the infrared (IR) CLUMPY model, allowing us to conduct the multiwavelength study with the X-ray torus model XCLUMPY self-consistently. Adopting the torus parameters obtained by the X-ray fitting, we estimate the properties of the host galaxies, active galactic nucleus (AGN) tori, and polar dust. The star formation rates (SFRs) become larger with merger stage and most of them are above the main sequence. The SFRs are correlated with radio luminosity, indicating starburst emission is dominant in the radio band. Although polar-dust extinction is much smaller than torus extinction, the UV-to-IR (mainly IR) polar dust luminosities are ~2 times larger than the torus ones. The polar-dust temperature decreases while the physical size, estimated by the temperature and dust sublimation radius, increases with AGN luminosity from a few tens of parsec (early mergers) to kiloparsec scales (late mergers), where the polar dust likely comes from expanding (i.e., evolving) dusty outflows. A comparison between the SFRs and intrinsic AGN luminosities suggests that starbursts occur first and AGNs arise later, and overall their growth rates follow the simultaneous coevolution local galaxy–SMBH mass relation. We confirm the coexistence of intense starbursts, AGNs, and large-scale outflows in late mergers, supporting a standard AGN feedback scenario.

79 ASTRONOMY AND ASTROPHYSICS↗

The coupled SYK model at finite temperature

Sachdev-Ye-Kitaev (SYK) model, which describes N randomly interacting Majorana fermions in 0+1 dimension, is found to be an solvable UV-complete toy model for holographic duality in nearly AdS 2 dilaton gravity. Ref. [1] proposed a modified model by coupling two identical SYK models, which at low-energy limit is dual to a global AdS 2 geometry. This geometry is an “eternal wormhole" because the two boundaries are causally connected. Increasing the temperature drives a Hawking-Page like transition from the eternal wormhole geometry to two disconnected black holes with coupled matter field. To gain more understanding of the coupled SYK model, in this work, we study the finite temperature spectral function of this system by numerical solving the Schwinger-Dyson equation in real-time. We find in the low-temperature phase the system is well described by weakly interacting fermions with renormalized single-particle gap, while in the high temperature phase the system is strongly interacting and the single-particle peaks merge. We also study the q dependence of the spectral function.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A FFT-based mesoscale continuum dislocation mechanics with defect energy: Applications to composites and polycrystals

A crystal plasticity elastoviscoplastic FFT (fast Fourier transform) formulation with a mesoscale continuum field dislocation mechanics model is presented, which incorporates a defect energy density that depends on GND densities and an associated material length scale. This allows to thermodynamically derive internal length scale dependent intra-crystalline backstress and Peach–Koehler force acting on GND densities. The model considers GND density evolution through a filtered numerical spectral approach, which is coupled with stress equilibrium through the elastoviscoplastic FFT algorithm. The discrete Fourier transform (DFT) method together with finite difference (FD) schemes is applied to solve both the backstress tensor and the Fourier–Green operator. Numerical results are first reported for two-phase laminate composites with plastic single crystal channels and elastic precipitates for shear loadings. Channel size effects are simulated and analyzed on the overall and local hardening behaviors during monotonous loadings. In addition, the evolutions of GND densities and the role of their associated backstress on size effects are examined during reversible shear loading. In a second part, the role of the defect energy internal length scale on polycrystal’s hardening during tension–compression is discussed. The results are compared to those obtained using FFT-based continuum field dislocation mechanics without defect energy.

36 MATERIALS SCIENCE↗