Search NASA⌕ Search

SEARCH · Search NASA

Results for “simulation and field testing”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

PROTOCALC, a W -band Polarized Calibrator for Cosmic Microwave Background Telescopes: Application to Simons Observatory and CLASS

Current- and next-generation cosmic microwave background (CMB) experiments will measure polarization anisotropies with unprecedented sensitivities. The need for high precision in these measurements underscores the importance of gaining a comprehensive understanding of instrument properties, with a particular emphasis on the study of the beam properties, and especially their polarization characteristics and the measurement of the polarization angle. In this context, a major challenge lies in the scarcity of millimeter polarized astrophysical sources with sufficient brightness and calibration knowledge to meet the stringent accuracy requirements of future CMB missions. This led to the development of a drone-borne calibration source designed for the frequency band centered on approximately 90 GHz, matching a commonly used channel in ground-based CMB measurements. The Prototype Calibrator for Cosmology, PROTOCALC, has undergone thorough in-lab testing, and its properties have been subsequently modeled through simulation software integrated into the standard Simons Observatory analysis pipeline. Moreover, the PROTOCALC system has been tested in the field, having been deployed twice on calibration campaigns with CMB telescopes in the Atacama Desert. The data collected constrain the roll angle of the source with a statistical accuracy of 0$^°_•$045.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Design study of a combined interferometer and polarimeter for a high-field, compact tokamak

This article is the first design study of a combined interferometer and polarimeter on a compact, high-field, high-density, net-energy tokamak. Recent advances in superconducting technology have made possible designs for compact, high magnetic field fusion power plants, such as ARC [Sorbom et al., Fusion Eng. Des. 100, 378 (2015)], and experiments, such as SPARC [Greenwald et al., PSFC Report No. RR-18-2 (2018)]. These new designs create both challenges and opportunities for plasma diagnostics. The diagnostic proposed in this work, called InterPol, takes advantage of unique opportunities provided by high magnetic field and density to measure both line-averaged density and poloidal magnetic field with a single set of CO 2 and quantum cascade lasers. These measurements will be used for fast density feedback control, constraint of density and safety factor profiles, and density fluctuation measurements. Synthetic diagnostic testing using a model machine geometry, called MQ1 (Mission Q ≥ 1), and profiles simulated with Tokamak Simulation Code indicate that InterPol will be able to measure steady state density and poloidal magnetic field, as well as fluctuations caused by toroidal Alfven eigenmodes and other phenomena on a high-field compact tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Data for Impact of Vertical and Seasonal Variation in Leaf Traits on Simulating Soybean Canopy Photosynthesis via 1D and 3D Modeling

Accurate modeling of photosynthesis is crucial for predicting crop productivity and quantifying the carbon cycle in agroecosystems. Leaf traits are essential inputs for modeling canopy photosynthesis. Yet, many existing models still use fixed plant functional type (PTF)-based values to parameterize leaf traits under a big-leaf or two-big-leaf assumption, neglecting their vertical profiles and seasonal changes. This simplification may introduce significant uncertainties in estimating gross primary productivity (GPP). In this study, we simulated soybean GPP and tested the effects of vertical and seasonal variation in three key leaf photosynthetic traits: the maximum carboxylation rate at 25 °C (Vcmax25), leaf chlorophyll content (LCC), and leaf mass per area (LMA) in the 1D-SCOPE and 3D-Helios models. Weekly field measurements were conducted during the growing season of 2024 to support the simulation. We designed ten leaf trait parameterization schemes by incorporating different combinations of vertical profiles and seasonal changes, while assuming homogeneous canopy architecture in both models. Our results revealed that Vcmax25 vertical and seasonal variation had the strongest influence on simulated GPP in both 1D and 3D models, while LCC and LMA effects were minimal. Particularly, the scheme with an empirically parameterized Vcmax25 profile achieved comparable performance to the scheme with the measured Vcmax25 profile. Both 1D-SCOPE and 3D-Helios accurately modeled GPP (SCOPE: R2 = 0.87, Bias = 0.55 µmol m⁻² s⁻¹; Helios: R2 = 0.9, Bias = 0.22 µmol m⁻² s⁻¹) under the most complex scheme, and their responses to vertical and seasonal variation in leaf traits were consistent, demonstrating the robustness of our findings. Based on our findings, we propose a scalable framework for parameterizing leaf traits to improve GPP simulations. This study contributes to improving the representation of leaf trait dynamics in canopy-level photosynthesis models, potentially enhancing our ability to predict crop productivity and understand agroecosystem carbon dynamics.

Photosynthesis↗
DOE OSTI2020

Scaling Law Describes the Spin-Glass Response in Theory, Experiments, and Simulations

The correlation length ξ, a key quantity in glassy dynamics, can now be precisely measured for spin glasses both in experiments and in simulations. However, known analysis methods lead to discrepancies either for large external fields or close to the glass temperature. We solve this problem by introducing a scaling law that takes into account both the magnetic field and the time-dependent spin-glass correlation length. Here, the scaling law is successfully tested against experimental measurements in a CuMn single crystal and against large-scale simulations on the Janus II dedicated computer.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

iMaNGA: mock MaNGA galaxies based on IllustrisTNG and MaStar SSPs – II. The catalogue

ABSTRACT Strengthening the synergy between simulations and observations is essential to test galaxy formation and evolution theories. To achieve this goal, in the first paper of this series, we presented a method to generate mock SDSS-IV/Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) integral-field spectroscopic galaxy observations from cosmological simulations. In this second paper, we build the iMaNGA catalogue consisting of ∼1000 unique galaxies from the TNG50 cosmological simulations, selected to mimic the SDSS-IV/MaNGA-Primary sample selection. Here, we present and discuss the iMaNGA sample and its comparison to the MaNGA Primary catalogue. The iMaNGA sample well recovers the MaNGA-Primary sample in terms of stellar mass versus angular size relation and spatial resolution. The Sérsic index versus angular size relation, instead, is not reproduced well by the simulations, mostly because of a paucity of high-mass elliptical galaxies in TNG50. We also investigate our ability to recover the galaxy kinematics and stellar population properties with full-spectral fitting. We demonstrate that ‘intrinsic’ and ‘recovered’ stellar kinematics, stellar ages, and metallicities are consistent, with residuals compatible with zero within 1σ. Also ‘intrinsic’ and ‘recovered’ star formation histories display a great resemblance. We conclude that our mock generation and spectral fitting processes do not distort the ‘intrinsic’ galaxy properties. Therefore, in the third paper of this series, we can meaningfully test the cosmological simulations, comparing the stellar population properties and kinematics of the iMaNGA mock galaxies and the MaNGA observational results.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Conservative discontinuous Galerkin interpolation: Sheared boundary conditions

Local studies of accretion disks and laboratory magnetized plasmas employ analytical coordinate mappings that introduce sheared boundary conditions (BCs). Here, we present a discontinuous Galerkin (DG) algorithm to apply such BCs based on projections and quadrature-free integration. The procedure is high-order accurate, preserves moments exactly and works in multiple dimensions. Tests of increasing complexity are provided, beginning with translations of one and two dimensional fields, followed by 3D and 5D simulations with sheared (twist-shift) BCs. Results show that the algorithm is (p + 1)-order accurate in the DG representation and (p + 2) -order accurate in the cell averages, with p being the order of the polynomial basis. Quantification of the algorithm's hyperdiffusion and discussion of aliasing errors are given. This technique enables conservative local simulations of plasma turbulence with DG, not possible until now.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Insights Into Seismicity Associated With Flexibly Operating Enhanced Geothermal System From Real‐Time Distributed Acoustic Sensing

Enhanced Geothermal Systems (EGS) have the capacity to broaden the accessible resource pool for geothermal power generation. Traditionally viewed as a “baseload” resource, their flexible operation might also enable dispatchable load‐following generation and long‐term energy storage, aligning them with the evolving landscape of decarbonized electricity systems. However, increasing permeability and extracting energy during EGS operations can induce microseismic events; for many prior EGS efforts, some associated seismicity has been observed. While energetically beneficial, the flexibility of EGS operations prompts our inquiry into whether new types of operations will yield previously unseen seismicity patterns. We demonstrate the use of distributed acoustic sensing (DAS) with real‐time edge computing to monitor seismicity during a pilot test of a cyclically operated EGS facility at the Blue Mountain geothermal field. Our focus lies in uncovering seismicity insights from the real‐time microseismic catalog, particularly during load‐following dispatchability tests simulating flexible EGS operation. Here, we find that variations in pore pressure consistently correlate with seismicity, and that controlling pressure cycles during flexible operations appears to constrain microseismic activity during subsequent cycles. The spatio‐temporal evolution of microseismic clouds recorded during cyclic injection cycles fits diffusive models over our available observation period. Additionally, seismicity elevation lags behind pore pressure increases, likely due to pressure diffusion to the fracture system boundary. Through real‐time monitoring, we offer novel insights into seismicity associated with flexibly operating EGS. Our findings suggest that leveraging DAS and edge computing can inform EGS operations and help mitigate induced seismicity.

Chamarczuk, Michal [Rice Univ., Houston, TX (Unite↗
DOE OSTI2022

Characterization of a 1D-imaging high-energy x-ray backlighter driven by the National Ignition Facility Advanced Radiographic Capability laser

Plastic deformation of samples compressed to Mbar pressures at high strain rates at the National Ignition Facility (NIF) forms the basis of ongoing material strength experiments in conditions relevant to meteor impacts, geophysics, armor development, and inertial confinement fusion. Hard x-ray radiography is the primary means of measuring the evolution of these samples, typically employing a slit-collimated high-Z microdot driven by the NIF laser to generate >40 keV x rays [E. Gumbrell et al., Rev. Sci. Instrum. 89, 10G118 (2018) and C. M. Huntington et al., Rev. Sci. Instrum. 89, 10G121 (2018)]. Alternatively, a dysprosium “micro-flag” target driven by the Advanced Radiographic Capability laser (~2 kJ, 10 ps) can deliver significantly higher spatiotemporal resolution [M. P. Hill et al., Rev. Sci. Instrum. 92, 033535 (2021)], especially in high-opacity samples. Here, initial experiments revealed problematic brightness and spectral gradients from this source, but by radiographing a set of diamond-turned, 105 µm-thick Pb test objects and supported by simulations using the 3D Monte Carlo code GEANT4, these geometry-dependent gradients across the field of view are quantified and mitigation strategies are assessed. In addition to significantly enhancing the modulation transfer function compared to the existing system, image stacking from multiple layers of image plate is shown to almost double the signal to noise ratio that will reduce uncertainties in future dynamic strength experiments.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Magnetized ICF implosions: non-axial magnetic field topologies

This paper explores four different magnetic field topologies for application to spherical inertial confinement fusion implosions: axial, mirror, cusp and closed field lines. A mirror field is found to enhance the impact of magnetization over an axial field; this is because the mirror field more closely follows the hot-spot surface. A cusp field, while simple to generate, is not found to have any benefits over the tried-and-tested axial field. Closed field lines are found to be of the greatest benefit to hot-spot performance, with the simulated design undergoing a 2× increase in ion temperature before α-heating is considered. The plasma properties of the simulation with closed field lines are radically different from the unmagnetized counterpart, with electron temperatures in excess of 100 keV, suggesting that a fundamental redesign of the capsule implosion is possible if this method is pursued.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Optimizing Bifacial PV Systems (Q3 FY2020 Project Report)

This project has four main technical objectives. Develop and improve bifacial performance models by adding the capability to evaluate electrical behavior and performance of bifacial modules and arrays under realistic field conditions including irradiance variability caused by racking, module frame, and position in the array. Instrument and monitor performance of fielded bifacial systems to validate performance models and to measure, analyze and publish on bifacial energy gain. These should include both research and commercial bifacial systems and cover a variety of deployment applications. Evaluate optimal bifacial system designs using simulations leveraging high performance computing, and also using full sized and miniaturized experimental field deployments. Establish and contribute to international test standards for bifacial system performance, testing, and safety, and work with the community to establish installation and siting best practices.

14 SOLAR ENERGY↗
DOE OSTI2026

Photon (Non)Conservation in the Reduced Speed of Light Approximation and How to (Almost) Fix It

The "Reduced Speed of Light" (RSL) approximation is commonly used to speed up radiative transfer calculations in cosmological simulations. However, it has been shown previously that the RSL approximation leads to photon non-conservation when the radiation field is rapidly evolving in time. I show that these missing photons can be counted exactly for some numerical schemes. Adding them back into a simulation, however, is a much harder task. I show one example of such a scheme, which achieves sub-percent accuracy on simple tests. Unfortunately, the scheme performs much worse on semi-realistic simulations of cosmic reionization, leading to a faster overlap and significant errors in the point-wise comparison of the RSL radiation field with the reference simulation that maintains the full speed of light for the radiative transfer.

Gnedin, Nickolay Y. [Fermilab; Chicago U., KICP; C↗
DOE OSTI2024

Dark energy survey year 3 results: likelihood-free, simulation-based w CDM inference with neural compression of weak-lensing map statistics

We present simulation-based cosmological wcold dark matter (wCDM) inference using dark energy survey year 3 weak-lensing maps, via neural data compression of weak-lensing map summary statistics: power spectra, peak counts, and direct map-level compression/inference with convolutional neural networks (CNN). Using simulation-based inference, also known as likelihood-free or implicit inference, we use forward-modelled mock data to estimate posterior probability distributions of unknown parameters. This approach allows all statistical assumptions and uncertainties to be propagated through the forward-modelled mock data; these include sky masks, non-Gaussian shape noise, shape measurement bias, source galaxy clustering, photometric redshift uncertainty, intrinsic galaxy alignments, non-Gaussian density fields, neutrinos, and non-linear summary statistics. We include a series of tests to validate our inference results. This paper also describes the Gower Street simulation suite: 791 full-sky pkdgrav3 dark matter simulations, with cosmological model parameters sampled with a mixed active-learning strategy, from which we construct over 3000 mock dark energy survey lensing data sets. For wCDM inference, for which we allow –1 < w < –$\frac{1}{3}$⁠, our most constraining result uses power spectra combined with map-level (CNN) inference. Using gravitational lensing data only, this map-level combination gives Ω m = 0.283$^{+0.020}_{–0.027}$⁠, S 8 = 0.804$^{+0.025}_{–0.017⁠}$, and w < –0.80 (with a 68 per cent credible interval); compared to the power spectrum inference, this is more than a factor of two improvement in dark energy parameter (Ω⁠ DE , w⁠) precision.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

GPU-accelerated kinetic Hall thruster simulations in WarpX

Abstract Two-dimensional (axial-azimuthal) simulations of a Hall thruster are performed using the open-source particle-in-cell code WarpX. The simulation conditions are chosen to match those of the axial-azimuthal benchmark first reported by Charoy et al. in 2019. A range of numerical and solver parameters is investigated in order to find those which yield the best performance. It is found that WarpX completes the benchmark case in 3.8 days on an Nvidia V100 GPU, and in as low as 1.5 days on a more recent Nvidia H100 GPU. Of the numerical parameters investigated, it is determined that the field-solve tolerance and particle resampling thresholds have the largest effect on the simulation wall time and that particle resampling may artificially widen electron velocity distribution functions, leading to unphysical heating. A semi-implicit scheme for the electrostatic field solve is tested and is found to produce results consistent to within 10% of the benchmark in less than twelve hours. The scaling properties of the electrostatic solver to multiple GPUs are also assessed on a uniform plasma test problem. The results of this work are discussed in the context of advancements in GPU hardware and the suitability of kinetic Hall thruster simulations for engineering applications.

Marks, Thomas A.↗
DOE OSTI2021

Effects of residue removal and tillage on greenhouse gas emissions in continuous corn systems as simulated with RZWQM2

Agricultural production is a major source of carbon dioxide (CO 2 ) and nitrous oxide (N 2 O) globally. The effects of conservation practices on soil CO 2 and N 2 O emissions remain a high degree of uncertainty. In this study, soil CO 2 and N 2 O emissions under different residue and tillage practices in an irrigated, continuous corn system, were investigated using the Root Zone Water Quality Model (RZWQM2). Combinations of no/high stover removal (NR and HR, respectively) and no-till/conventional tillage (NT and CT, respectively) field experiments were tested over the four crop-years (Apr. 2011–Apr. 2015). The model was calibrated using the NRCT, and validated with other treatments. The simulation results showed that soil volumetric water content (VWC) in the NR treatments (i.e., NRCT and NRNT) was 1.3%–1.9% higher than that in the HR treatments (i.e., HRCT and HRNT) averaged across the four years. A higher amount of CO 2 and N 2 O emissions were simulated in the NRCT across the four years (annual average: 7034 kg C/ha/yr for CO 2 and 3.8 kg N/ha/yr for N 2 O), and lower emissions were in the HRNT (annual average: 6329 kg C/ha/yr and 3.7 kg N/ha/yr for N 2 O). A long-term simulation (2001–2015) suggested that the CO 2 and N 2 O emissions were closely correlated with the stover removal degree (SRD), tillage, VWC, soil temperature (ST), years in management (Y), and fertilizer application. Stover and tillage practices had cumulative effects on CO 2 emissions. The simulated annual CO 2 emissions in 1st year from NRCT, NRNT, and HRCT were 7.8%, 0.0%, and 7.7% higher than that from HRNT, respectively; then the emissions in 15th year were 63.6%, 47.7%, and 29.1% higher, respectively. Meanwhile, there were no cumulative effects on N 2 O emissions. Furthermore, the results also demonstrated that the RZWQM2 is a promising tool for evaluating the long-term effects of CO 2 and N 2 O emissions on different conservation practices.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Microstructure Experiments-Enabled MARMOT Simulations of SiC/SiC-based Accident Tolerant Nuclear Fuel System

We have undertaken an experimental-computational project that addresses a few key technology gaps associated with the use of SiC/SiC composites for light water reactor fuel cladding. Two principal endeavors in this project are to assess the irradiation-induced microstructural changes and swelling in SiC/SiC composites, and to characterize and model the porous oxide surface layer (when exposed to steam) along SiC recession that is detrimental to the clad integrity under accident conditions. The project tasks are: (i) ion irradiation and characterization (University of Tennessee, Knoxville), (ii) mechanical tests and analysis (University of South Carolina), (iii) steam exposure tests (Oak Ridge National Laboratory), (iv) electron microscopy and spectroscopic characterization (NC State University), (v) x-ray microscopy and reconstruction (University of South Carolina and NC State University), (vi) phase field modeling and simulations (NC State University and Idaho National Laboratory). The effects of 10 MeV Au ion irradiation at 350°C on the microstructure evolution in SiC/SiC composites are investigated at doses up to 400 displacements per atom (dpa) at the University of Tennessee, Knoxville. Atomic force microscopy and optical profilometry reveal irradiation induced axial and radial shrinkage of the fibers for doses greater than 10 dpa. Based on detailed electron microscopy characterization, the primary cause of the fiber shrinkage is attributed to irradiation-induced loss of carbon packets. Additionally, the multilayer PyC interface is observed to portray high resistance to irradiation damage. The mechanical response and failure mechanisms of un-irradiated samples is also assessed through loading tests and X-ray imaging at the University of South Carolina. Steam exposure tests are performed at the Oak Ridge National Laboratory in a facility that represents a reactor pressure vessel under a loss-of-coolant-accident scenario. The samples that are analyzed methodically in this report are tested for 32/31 hours at 1200°C with a velocity of 0.25 cm/s for the pressures: 0.1 MPa, 0.45 MPa, 0.92 MPa and 1.38 MPa. Scanning/transmission electron microscopy analysis conducted at the NC State University (NCSU) shows that the oxide layer thickness increases with the steam pressure. While the oxide layer is crystalline (α- cristobalite) for the pressures 0.45 MPa, 0.92 MPa and 1.38 MPa, the layer is amorphous at 0.1 MPa. Results from Raman spectroscopy have also confirmed the formation of α-cristobalite phase of SiO₂. The abrupt increase in the integrated Raman intensity ratio between steam pressures 0.1 and 0.45 MPa suggests the onset of accelerated crystallization. Non-destructive three dimensional X-ray microscopy/tomography (XCT) and computational image processing techniques are employed by the University of South Carolina to probe the porous oxidation features on SiC samples at varying pressures. Interestingly, most pores are observed to be located away from the surface as well as the oxide-SiC interface. The average oxide layer thickness assessed from the XCT analysis is seen to be in excellent agreement with the values determined through scanning electron microscopy analysis for the highest pressures where the oxide layer is relatively more uniform. A phase-field model developed by the NC State University and Idaho National Laboratory for simulating oxidation of SiC by steam captures the paralinear kinetics of SiC oxidation with a high degree of fidelity. Results from quasi-one-dimensional and two dimensional simulations show that the pores with oxidizing species lead to a higher volatilization rate. These simulations indicate that the enhanced apparent volatilization rates at higher pressures observed in experiments can be rationalized by the increased volatilization from the pores. The project team has also successfully developed the capability to import images from experiments for realistic evolution of the oxidizing microstructure.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Activation Energy for Solder Bond Degradation: Thermal Cycling of Field-aged Modules at Multiple Upper Temperatures

The reliability of solder joints in the solar cell metallization-interconnect system influences the lifetime of photovoltaic modules. Two field-aged modules-one with Sn 62 Pb 36 Ag 2 solder at the solder joints (Solarex MSX 60), and the other with the standard Sn 60 Pb 40 solder (Siemens M55)-were subjected to a modified thermal cycling (TC) test of IEC 61215. Three sections in each module were maintained at 85 degrees C, 95 degrees C, and 105 degrees C during the 15-minute high temperature dwell time. Current equivalent to the module short-circuit current was injected through the module when the chamber temperature was above 25 degrees C to simulate regular field operation. This novel approach aims to induce thermomechanical fatigue (TMF) at the solder joints and intermetallic compound (IMC) formation at the metal/solder interfaces. The activation energy (E a ) for solder bond degradation was calculated based on the series resistance (R s ) increase in TC testing rather than power drop to avoid the effect of confounding variables. Module-level R s increase in MSX 60 module after TC800 cycles and in M55 module after TC400 cycles was 1.22% and 183.7%, respectively. The E a determined for the module with 2wt% Ag is 0.24 eV and for that with the standard solder is 0.27 eV. The solder bond degradation seems to have been driven by TMF rather than IMC formation for both modules.

14 SOLAR ENERGY↗
DOE OSTI2020

Dynamics and Control of Microgrids as a Resiliency Source

Microgrids can be operated in a utility-connected mode with the neighboring distribution system or in an islanded mode in separation with the distribution system. As the major disasters occur more and more often, islanded microgrid is proposed as a possible solution to serve essential loads until the utility service is restored. This paper evaluates the feasibility of using microgrids as a resiliency resource by a field test. A control scheme, integrating droop control and feedback control, is applied on a modified IEEE 13-node test feeders. The simulation results validate the proposed control scheme and the transient stability

Lee, Lung-An↗
DOE OSTI2025

Theoretical Investigations of Anharmonic Effects and Phonon Transport in the Cubic Phase of Crystalline Perovskite CsPbCl 3

The role of anharmonic effects on lattice dynamics and thermal transport has been investigated in the cubic phase of the CsPbCl 3 perovskite. Limitations of the harmonic approximation, which lead to phonon instabilities in the Brillouin zone, were addressed based on self-consistent phonon theory in combination with first-principles calculations and by incorporation of frequency renormalization effects from bubble and loop diagrams. This theoretical approach demonstrates significant improvements compared to results obtained using a self-consistent phonon dynamic matrix. Both the cubic-to-tetragonal phase transition temperature as well as the predicted lattice thermal conductivity values show good agreement with selected sets of available experimental data. In conclusion, the accuracy of the predicted thermal conductivity data is also tested against systems significantly larger than those allowed by quantum calculations, by using molecular dynamic simulations with machine-learning force fields.

Lattices↗