VnV: A Self Documenting Testing Framework for In-situ Verification and Validation in High Performance Computing Applications
Explore the source record for details and available documents.
96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
SEARCH · Search NASA
Results for “performance verification”
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.
Validation of SAS4A/SASSYS-1 for predicting steady-state single-phase natural circulation
The validation of system analysis codes for nuclear reactor systems is required for the development and application of these computational tools. Designed as a comprehensive system analysis code for advanced nuclear reactors, SAS4A/SASSYS-1 requires validation of its physics model for capturing single-phase natural circulation behavior. To support the validation of SAS4A/SASSYS-1, high-precision experiments are performed capturing steady-state single-phase natural circulation on a scaled facility with comprehensive instrumentation. In this work, dedicated tests are performed quantifying the critical modeling parameters, and a single-phase natural circulation benchmark dataset is obtained with well-documented uncertainty and comprehensive facility description. The validation is then performed against the dataset examining the capability of SAS4A/SASSYS-1 in simulating steady-state single-phase natural circulation. The experimental facility is modeled in the candidate code. Solution verification is performed using Richardson-extrapolation-based estimators which quantify and restrict numerical errors from discretization. Input uncertainty provided by the benchmark dataset is forward propagated through the candidate code, quantifying the output uncertainty in a Monte Carlo approach. The composition of the output uncertainty is also quantified through a variance-based sensitivity analysis. With the uncertainty quantified for each individual condition, a detailed comparison between the simulation results and experimental data is performed covering the whole dataset. The results show consistent agreement for all primary parameters. The current validation activity provides a valuable benchmark dataset for the validation of system analysis codes in capturing single-phase natural circulation and demonstrates satisfactory prediction capability of SAS4A/SASSYS1 for steady-state single-phase natural circulation.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
Validation of SAS4A/SASSYS-1 for Steady-State Single-Phase Natural Circulation
The development of numerical software for engineering applications requires the validation of the code against experimental benchmark datasets. To support the validation of SAS4A/SASSYS-1 in simulating the single-phase natural circulation within nuclear reactor systems, high-precision steady-state experiments are performed to capture single-phase natural circulation phenomena on an existing scaled facility with comprehensive instrumentation. Forced convection tests are designed precisely capturing the facility’s critical thermal-hydraulic parameters required for one-dimensional modeling, and a comprehensive single-phase natural circulation dataset is obtained with well-documented experimental uncertainty and facility description. Analysis and discussion based on the finalized dataset confirm the dataset’s ability in capturing dominant phenomena and address important physical interpretation of parameters. The dataset reported in this project provides a valuable benchmark resource for the validation of system analysis codes under single-phase natural circulation inside nuclear reactors. Validation of SAS4A/SASSYS-1 is then performed against the obtained benchmark dataset to confirm the capability of its physics model in simulating the steady-state single-phase natural circulation. The experimental facility is modeled in the candidate code. Solution verification is performed using Richardson-extrapolation-based estimators to quantify and restrict numerical errors from discretization. Numerical uncertainty originating from finite maximum pseudo-transient time is also quantified and restricted. Input uncertainty provided by the benchmark dataset is forward propagated through the candidate code, directly quantifying the simulation output in a Monte Carlo approach. The composition of the output uncertainty is further determined by the estimators of Sobol’ indices through a variance-based sensitivity analysis. With uncertainty quantified for each individual condition, detailed comparison between the simulation results and experimental data is performed covering the whole dataset, which shows consistent agreement for all important quantities. The validation activity reported in this project demonstrates that SAS4A/SASSYS-1 can predict the primary parameters with satisfactorily accuracy under steady-state single-phase natural circulation.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
Computational Modeling of Molten Salt Infiltration and Oxidation in Nuclear Graphite
Graphite is utilized as a moderator and reflector in advanced nuclear reactor designs due to its high thermal conductivity, neutron moderation properties, and resistance to radiation damage. However, its longterm performance and reliability are challenged by degradation mechanisms such as molten salt infiltration in molten salt reactors (MSRs) and oxidation in gas-cooled reactors (GCRs). These mechanisms can compromise the structural integrity and operational lifetime of graphite components, necessitating a more detailed assessment of their physical behavior. This report focuses on the development of computational models for molten salt infiltration and oxidation of graphite to aid the design and performance analysis of graphite components. For molten salt infiltration, a computational framework is developed that couples incompressible Navier-Stokes and phase-field model to simulate the penetration of molten salt into graphite?s interconnected pore structure. Initial model verification is performed using two-phase flows in two dimensions, demonstrating the models ability to capture fundamental physical behavior and agree with analytical solution. This framework is then applied to a realistic IG110 nuclear graphite , where a computed tomography extracted pore geometry is used to analyse the infiltration behavior of FLiNaK molten salt. This model provides insights into how the microstructure and other relevant parameters influence the transport pathways of molten salt into graphite, potentially offering a means to rapidly evaluate a graphite grade?s resistance to infiltration. For oxidation, the report details pore-scale mass and heat transport models, describing the diffusion of gases, reaction kinetics, and thermal effects. Additionally, this report highlights inconsistencies in the existing volume-averaged macroscopic model, particularly in upscaling of reaction kinetics and flux terms, and surface to volume transformations. These inconsistencies suggest that current formulations may not accurately capture the experimentally observed graphite oxidation process, highlighting the need for improved model development. This work advances the development of physics-based computational models for graphite degradation, contributing to improved predictive models for next-generation nuclear reactor designs. Future efforts will focus on refining the infiltration model to address non-physical behaviors and enhance its robustness. Additionally, for oxidation, further studies will employ the principles of volume averaging to rigorously derive the upscaled equations, potentially in collaboration with subject matter experts.
Computational Modeling of Molten Salt Infiltration↗
The eROSITA Final Equatorial-Depth Survey (eFEDS)
Here, we investigate the physical properties – such as the stellar mass (M * ), star-formation rate, infrared (IR) luminosity (L IR ), X-ray luminosity (L X ), and hydrogen column density (N H ) – of mid-IR (MIR) galaxies and active galactic nuclei (AGN) at z < 4 in the 140 deg 2 field observed by eROSITA on SRG using the Performance-and-Verification-Phase program named the eROSITA Final Equatorial Depth Survey (eFEDS). By cross-matching the WISE 22 μm (W4)-detected sample and the eFEDS X-ray point-source catalog, we find that 692 extragalactic objects are detected by eROSITA. We have compiled a multiwavelength dataset extending from X-ray to far-IR wavelengths. We have also performed (i) an X-ray spectral analysis, (ii) spectral-energy-distribution fitting using X-CIGALE, (iii) 2D image-decomposition analysis using Subaru Hyper Suprime-Cam images, and (iv) optical spectral fitting with QSFit to investigate the AGN and host-galaxy properties. For 7088 WISE 22 μm objects that are undetected by eROSITA, we have performed an X-ray stacking analysis to examine the typical physical properties of these X-ray faint and probably obscured objects. We find that (i) 82% of the eFEDS–W4 sources are classified as X-ray AGN with log L X > 42 erg s –1 ; (ii) 67 and 24% of the objects have log(L IR /L ⊙ ) > 12 and 13, respectively; (iii) the relationship between L X and the 6 μm luminosity is consistent with that reported in previous works; and (iv) the relationship between the Eddington ratio and N H for the eFEDS–W4 sample and a comparison with a model prediction from a galaxy-merger simulation indicates that approximately 5.0% of the eFEDS–W4 sources in our sample are likely to be in an AGN-feedback phase, in which strong radiation pressure from the AGN blows out the surrounding material from the nuclear region. Thanks to the wide area coverage of eFEDS, we have been able to constrain the ranges of the physical properties of the WISE 22 μm-selected sample of AGNs at z < 4, providing a benchmark for forthcoming studies on a complete census of MIR galaxies selected from the full-depth eROSITA all-sky survey.
79 ASTRONOMY AND ASTROPHYSICS↗
A conservative, interface-resolved, compressible framework for the modeling and simulation of liquid/gas phase change
This paper presents a method for simulating evaporation in a compressible, interface-resolved framework appropriate for modeling problems of engineering interest. In order to achieve robustness and broad applicability, the method has been designed to discretely enforce consistent mass and thermal energy transport at the phase interface, to globally conserve mass, momentum, and energy, and to be capable of modeling compressible and incompressible systems. Verification is performed via the Sod-shock test, one-dimensional heat conduction, evaporation from a planar interface, and evaporation of three-dimensional droplets. Convergence with increasing mesh resolution is demonstrated in all tested configurations, and conservation is maintained near machine precision for a translating droplet. Conservation and accurate phase change rates are preserved at the low numerical resolutions commonly encountered in engineering calculations. Following verification, the method is validated by comparison to an empirical correlation for evaporating droplets in high temperature crossflow, and the presentation concludes with the simulation of an iso-octane spray at conditions representative of gasoline direct injection. In conclusion, successful verification, validation, and demonstrated practical utility suggest the method to be an accurate, efficient, and robust approach for the study of phase change in engineering systems.
97 MATHEMATICS AND COMPUTING↗
In anaerobic reactors the microbial community structure depends on feed type, with no “keystone” species tied to COD removal
Two-stage anaerobic digestion (AD) systems provide treatment for high strength wastewater with high stability and performance. Encapsulation technology can intensify AD to facilitate the separation of the solids retention time from the hydraulic retention time (HRT), offering lower HRTs, smaller reactors, and high effluent quality. To support successful deployment, however, the encapsulated community must contain all the needed microorganisms for successful treatment and be flexible enough to treat a variety of wastewaters. Here, a two-stage system was investigated in which microbial cultures were enriched on various high-strength wastewaters in suspended flow-through systems to determine how feed type influenced performance and microbial community structure. The hypothesis was that specific genera, or so-called “keystone species” would positively correlate to organic carbon degradation for a given feed, enabling construction of a well-functioning community for encapsulation. Results showed that the number of total bacteria (as 16S rRNA gene copies) did not correlate to soluble chemical oxygen demand (sCOD) removal, indicating that the community structure and/or members were important for good performance. Results also showed that feed type strongly influenced carbon removal and microbial community structure for 1st-stage fermenting communities, but not 2nd-stage methanogenic communities. In this study, the “core” community members were defined as organisms common to all of either the 1st- or 2nd-stage reactors irrespective of the feed they received and were present in at least 50% of the samples throughout the entire experiment. “Unique” community members were specific to a single feed, and hence, only present in either the 1st- or 2nd-stage reactors receiving that feed. In both 1st- and 2nd-stage communities, only one core genera and no unique genera were positively and significantly correlated to sCOD removal. Verification experiments performed with encapsulated communities showed that organisms identified in flow-through system and correlated with carbon degradation, though not significantly, seemed to be important for performance. Our results suggest that one cannot construct a community containing specific populations in lieu of enrichment. Nevertheless, a single diverse encapsulated anaerobic community should provide good (>80%) carbon removal when fed a variety of influents, if time is provided for enrichment after deployment.
54 ENVIRONMENTAL SCIENCES↗
Implementation of Isotopic Removal Capability in Griffin for Multi-Region MSR Depletion Analysis
This work reports the first step in such an analysis, the verification of isotope removal capability for multi-region MSR depletion in Griffin. Verification studies are performed against SERPENT for chemical species removal from the core to an off-gas system for MSR analysis. A brief review of the Bateman equations including the isotopic removal term is provided along with a description of the example problem used in the verification studies performed in Griffin versus SERPENT. Next, the results of the verification study are presented for isotopic depletion, isotopic depletion with removal, and the specific analysis case of varying the solubility of Iodine in a salt system and how this can affect material buildup in the off-gas system. Lastly, the results are analyzed and future work detailing the fully coupled MSR depletion analysis that incorporates the true multiphysics of the system is discussed.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
Verification, Validation, and Calibration Through a Causal Lens
While typical validation and verification approaches focus on identifying the associations between data elements using statistical and machine learning methods, the novel methods in this paper focus instead on identifying causal relationships between data elements. Statistical and machine-learning-based approaches are strictly data-driven, meaning that they provide quantitative comparison measures between data sets without explicitly considering the hypotheses behind them. This can lead to the erroneous conclusion that, if two data sets are close enough, the models that generated them are similar. In addition, when experimental and simulated data differ to an extent that fails to meet the acceptance criteria, calibration techniques are used to tweak simulation model parameters to reduce the gap between the two types of data. This produces the false expectation that a simulation model will match reality. The methods presented in this paper move away from these strictly data-driven methods for validation and calibration toward more robust, model-driven methods based on causal inference. Causal inference aims to identify the possible mechanisms that might have generated data. Thus, this analysis targets the prediction of the effects when one (or more) of the identified mechanisms are altered. There are many approaches to identify, quantify, and illustrate causal relationships. For the scope of this paper, directed graphs are employed as causal models. If the directed graph lacks cycles, it is known as a directed acyclic graph. A node in such a graph represents an observed data element while a directed edge connecting two nodes represents a causal relationship between two variables. The developed causal methods are designed to extract causal models from simulation models and experimental data. Causal models capture the causal relationships between data elements (e.g., simulated and experimental data). In this context, validation and verification are performed by comparing causal models. The proposed approach does not only inform system analysts on how a simulation model matches real-world data, but also identifies elements of the simulation model that should be revised when discrepancies between simulation and experimental data are observed. Through these causal methods, analysts can identify the portion of the model equation(s) that are behind an edge connecting two variables. Hence, once the structural differences between causal models have been determined, model calibration can occur by changing only those model parameters that impact the identified causal relationships.
97 MATHEMATICS AND COMPUTING↗
Skill of Mesoscale Models in Forecasting Springtime Macrophysical Cloud Properties at the Savannah River Site in the Southeastern US
Predicting boundary layer clouds is important for the accurate modeling of pollutant dispersion. Higher resolution mesoscale models would be expected to produce better forecasts of cloud properties that affect dispersion. Using ceilometer observations, we assess the skill of two operational mesoscale models (RAMS and WRF) to forecast cloud base altitude and cloud fraction at the Savannah River Site in the southeastern US during the springtime. Verifications were performed at small spatial and temporal scales necessary for dispersion modeling. Both models were unreliable with a 50% (RAMS) and a 46% (WRF) rate of predicting clouds observed by the ceilometer which led to low cloud fraction predictions. Results indicated that WRF better predicted daytime cloud bases from convection that occurred frequently later in the period and RAMS better predicted nighttime cloud bases. Using root mean squared error (RMSE) to score the forecast periods also highlighted this diurnal dichotomy, with WRF scores better during the day and RAMS scores better at night. Analysis of forecast errors revealed divergent model cloud base biases—WRF low and RAMS high. A hybrid solution which weighs more heavily the RAMS nighttime forecasts and WRF daytime forecasts will likely provide the best prediction of cloud properties for dispersion.
54 ENVIRONMENTAL SCIENCES↗
Progress with the 5D full-F continuum gyrokinetic code COGENT
COGENT is an Eulerian gyrokinetic code being developed for edge plasma modelling. The code is distinguished by the use of a high-order finite-volume (conservative) discretization combined with mapped multi-block grid technology. Our recent work is focused on the development of a 5D full-F COGENT version.Anumerical algorithm utilizing locally a field-aligned multi-block coordinate system is implemented to facilitate simulations of highly anisotropic microturbulence in the presence of a strong magnetic shear. In this approach, the toroidal direction is divided into blocks such that,within each block, the cells are field-aligned and a non-matching (non-conformal) grid interface is allowed at the block boundaries. In this paper we report on details of the numerical implementation and present preliminary results of verification studies performed for the case of the ion temperature gradient (ITG) instability in a sheared toroidal annulus geometry.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Developing Ultrahigh-Resolution E3SM Land Model for GPU Systems
Designing and refactoring complex scientific code, such as the E3SM land model (ELM), for new computing architectures is challenging. This paper presents design strategies and technical approaches to develop a data-oriented, GPU-ready ELM model using compiler directives (OpenACC/OpenMP). We first analyze the datatypes and processes in the original ELM code. Then we present design considerations for ultrahigh-resolution ELM (uELM) development for massive GPU systems. These techniques include the global data-oriented simulation workflow, domain partition, code porting and data copy, memory reduction, parallel loop restructure and flattening, and race condition detection. We implemented the first version of uELM using OpenACC targeting the NVidia GPUs in the Summit supercomputer at Oak Ridge National Laboratory. During the implementation, we developed a software tool (named SPEL) to facilitate code generation, verification, and performance tuning using these techniques. The first uELM implementation for Nvidia GPUs on Summit delivered promising results: 1) over 98% of the ELM code was automatically generated and tuned by scripts. Most ELM modules had better computational performances than the original ELM code for CPUs. The GPU-ready uELM is more scalable than the CPU code on fully-loaded Summit nodes. Example profiling results from several modules are also presented to illustrate the performance improvements and race condition detection. The lessons learned and toolkit developed in the study are also suitable for further uELM deployment using OpenMP on the first US exascale computer, Frontier, equipped with AMD CPUs and GPUs.
Schwartz, Peter↗
A methodology for corrosion testing ODS steels in liquid tin under reactor conditions in HFIR
Liquid metal–based divertor concepts have promising attributes for application in fusion reactors. Liquid tin is being considered as one of those liquid metals due to its excellent thermophysical properties. However, tin is extremely corrosive to steels at elevated temperatures, and the coupled effects of corrosion with neutron irradiation for application in fusion reactors have not been quantified. Researchers at Oak Ridge National Laboratory (ORNL) have designed and irradiated a series of capsules in the High Flux Isotope Reactor (HFIR) to help address gaps in the literature with respect to coupled corrosion and neutron irradiation effects. The capsules are filled with solid tin designed to melt due to the gamma heating in HFIR and allow interaction with the specimens. The capsules contain a SiC thermometer for post-irradiation temperature verification. Five capsules were irradiated in HFIR for 10.5 days, accumulating 0.93 dpa in the FeCrAl specimens. Post-irradiation temperature verification was performed on the SiC thermometers, and temperatures were observed to be higher than what was expected from thermal models. An explanation for the model underprediction is given herein. In conclusion, the capsules described within this article are the first to irradiate liquid metal within the HFIR flux trap and demonstrated the safety basis for continued liquid coolant studies in HFIR.
FeCrAl↗
Radiation hot spot formation on UF 6 30B cylinders
The accurate material accountancy of a country's enriched uranium hexafluoride (UF 6 ) stockpiles is imperative to the International Atomic Energy Agency (IAEA) for determining appropriate material safeguards. Uranium hexafluoride is typically stored in large cylinders with enrichment verification typically performed with a NaI gamma detector using the enrichment meter method. This approach for calculating enrichment assumes that the UF 6 is homogeneous inside the cylinder. Here, the measurements and tests presented in this manuscript show that the UF 6 in cylinders left outside in the elements can experience fractionation, leading to inhomogeneity, and subsequent development of radiation hot spots. These hot spot locations can produce errant enrichment measurements. In dark-colored cylinders that experience significant surface heating from the sun, UF 6 sublimates off the cylinder walls, leaving only the uranium daughter products in a horizontal patch (stripe) along the cylinder's midline. This daughter patch is hypothesized to be the source of the radiation hot spots. The patch forms during the summer and tends to decay away during the winter, when solar intensity and temperatures are reduced.
46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
Measuring very low radiation doses in PTFE for nuclear forensic enrichment reconstruction
Every country that has made nuclear weapons has used uranium enrichment to do so. Despite the centrality of this technology to international security, there is still no reliable physical marker of past enrichment in the open literature that can be used to perform forensic verification of historically produced weapons on gas centrifuges. We show that the extremely low radioactivity from uranium alpha emissions during enrichment leaves detectable and irreversible calorimetric signatures in the common enrichment gasket material PTFE, allowing for historical reconstruction of past enrichment activities at a sensitivity better than one weapon’s quantity of highly enriched uranium. Fast scanning calorimetry also enables the measurement of recrystallization enthalpies of sequentially microtomed slices, confirming the magnitude and the type of radiation exposure while also providing detection of tampering and a method for analyzing field samples useful for treaty verification. Furthermore, this work opens the door for common items to be turned into precise dosimeters to detect the past presence of radioactivity, nuclear materials, and related activities with high confidence.
98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
The eROSITA Final Equatorial-Depth Survey (eFEDS)
Theoretical models of the co-evolution of galaxies and active galactic nuclei (AGNs) ascribe an important role in the feedback process to a short, luminous, obscured, and dust-enshrouded phase during which the accretion rate of the supermassive black hole is expected to be at its maximum and the associated AGN-driven winds are also predicted to be maximally developed. Here, to test this scenario, we have isolated a textbook candidate from the eROSITA Final Equatorial-Depth Survey (eFEDS) obtained within the performance and verification program of the eROSITA telescope on board the Spectrum Röntgen Gamma mission. From an initial catalogue of 246 hard X-ray selected sources that are matched with the photometric and spectroscopic information available within the eROSITA and Hyper Suprime-Cam consortia, three candidates quasars in the feedback phase have been isolated applying a diagnostic proposed previously. Only one source (eFEDS J091157.4+014327) has a spectrum already available (from SDSS-DR16, z = 0.603) and it unambiguously shows abroad component (full width at half maximum ~1650 kms –1 ) in the [OIII]5007 line. The associated observed L [OIII] is ~2.6 × 10 42 erg s –1 , one to two orders of magnitude higher than that observed in local Seyfert galaxies and comparable to those observed in a sample of z ~ 0.5 type 1 quasars. From the multi-wavelength data available, we derive an Eddington ratio (L bol /L Edd ) of ~0.25 and a bolometric correction in the hard X-ray band of k bol ~ 10, which is lower than the corrections observed for objects at similar bolometric luminosity. These properties, along with the outflow, the high X-ray luminosity, the moderate X-ray obscuration (L X ~10 44.8 erg s –1 , N H ~2.7 × 10 22 cm –2 ), and the red optical colour, all match the prediction of quasars in the feedback phase from merger-driven models. Forecasting to the full eROSITA all-sky survey with its spectroscopic follow-up, we predict that by the end of 2024, we will have a sample of few hundred such objects at z= 0.5–2.
79 ASTRONOMY AND ASTROPHYSICS↗
Dynamic inflow and unsteady aerodynamics models for modal and stability analyses in OpenFAST
This work presents two aerodynamic models that are used to perform modal and stability analyses of wind turbines. Including aerodynamic effects in stability analyses is important to capture aerodynamic inertia, stiffness, and damping. The two models consist of a dynamic inflow model based on Øye’s model and an unsteady aerodynamics model based on a simplified Beddoes-Leishman dynamic stall model. Both models are written in a continuous-time state-space form, which is the form needed to perform linearizations within OpenFAST. In this article, we present the new models, perform some simple verifications against HAWC2 results, and then perform stability analyses. We investigate the contribution of the aerodynamic effects on the turbine modes, particularly the flapwise modes, which are the most affected. Because of the development presented in this work, Campbell diagrams accounting for dynamic stall and dynamic wake effects were generated for the first time using OpenFAST.
17 WIND ENERGY↗
Computational Evaluation of Thermal Response of Open-Cell Foam With Circular Pore (Computional Evaluation of Thermal Response of Open-Cell Foam with Circular Pore)
The evaluation of effective material properties in heterogeneous materials (e.g., composites or multicomponent structures) critically relevant to a wide spectrum of applications, including nuclear power, electronic packaging, flame retardants, hypersonics, and gas turbine power. The work described in this paper is centered around the numerical assessment of the thermal behavior of porous materials obtained from finite element thermal modeling and simulation. Here, two-dimensional, steady state analyses were performed on unit cells with centered, circular pores using a second order accurate Galerkin finite element method (FEM). The effective thermal conductivities of the porous systems were examined, encompassing a range of porosities from 4.9% to 60.1%. The geometries of the models were generated based on ordered circular pores for each modeled porosity level. The system response quantity (SRQ) under investigation was the dimensionless effective thermal conductivity across the unit cell. The dimensionless effective thermal conductivity was compared across all simulated cases, producing a trend between porosity and effective thermal conductivity. In the presented investigation, the method of manufactured solutions (MMS) was used to perform code verification, and the grid convergence index (GCI) was employed to estimate discretization uncertainty as solution verification. Code verification concluded an approximately second order accurate Galerkin FEM solver. It was found that the introduction of porosity to the unit cell material structure reduces effective thermal conductivity, as anticipated. Numerical results obtained in this study are compared to an analytical solution and to a sample of empirical data.
36 MATERIALS SCIENCE↗