Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “compositeness”

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 973 records · Page 54

Magnetism in metastable and annealed compositionally complex alloys

Compositionally complex materials (CCMs) present a potential paradigm shift in the design of magnetic materials. These alloys exhibit long-range structural order coupled with limited or no chemical order. As a result, extreme local environments exist with a large variations in the magnetic energy terms, which can manifest large changes in the magnetic behavior. In the current work, the magnetic properties of (Cr, Mn, Fe, Ni) alloys are presented. These materials were prepared by room-temperature combinatorial sputtering, resulting in a range of compositions with a single bcc structural phase and no chemical ordering. The combinatorial growth technique allows CCMs to be prepared outside of their thermodynamically stable phase, enabling the exploration of otherwise inaccessible order. The mixed ferromagnetic and antiferromagnetic interactions in these alloys causes frustrated magnetic behavior, which results in an extremely low coercivity (<1mT), which increases rapidly at 50 K. At low temperatures, the coercivity achieves values of nearly 500 mT, which is comparable to some high-anisotropy magnetic materials. Further, commensurate with the divergent coercivity is an atypical drop in the temperature dependent magnetization. These effects are explained by a mixed magnetic phase model, consisting of ferro-, antiferro-, and frustrated magnetic regions, and are rationalized by simulations. A machine-learning algorithm is employed to visualize the parameter space and inform the development of subsequent compositions. Annealing the samples at 600 °C orders the sample, more-than doubling the Curie temperature and increasing the saturation magnetization by as much as 5×. Simultaneously, the large coercivities are suppressed, resulting in magnetic behavior that is largely temperature independent over a range of 350 K. The ability to transform from a hard magnet to a soft magnet over a narrow temperature range makes these materials promising for heat-assisted recording technologies.

36 MATERIALS SCIENCE↗

Interface-independent sound speed and thermal conductivity of atomic-layer-deposition-grown amorphous AlN / Al 2 O 3 multilayers with varying oxygen composition

Dielectric amorphous multilayers (AMLs) play a critical role in a wide array of technologies such as optical coatings, nanoelectronics, energy harvesting, and recovery devices. However, despite their wide applications, a robust understanding of the effect of the interplay between chemical and structural disorder on the thermal properties of AMLs is still lacking. Therefore, in this paper, we experimentally and numerically investigate the effects of composition and interface density on the sound speed and thermal conductivity of a series of amorphous aluminum nitride and aluminum oxide multilayers grown via plasma-enhanced atomic layer deposition. To systematically change the composition, the oxygen content of the AMLs is proportionally varied with interface density during growth. We find that the longitudinal sound speed of these AMLs is dictated by the oxygen content instead of the number of interfaces. The thermal conductivity, in contrast, is dictated by both interface density and oxygen content. The interfaces act to decrease the thermal conductivity, whereas the oxygen content increases the thermal conductivity. Due to the competing influence of the interfaces and oxygen content, the thermal conductivity of the AMLs remains nearly constant as a function of interface density. Furthermore, our study provides crucial insights into the effect of the interplay of composition and interfaces on the sound speed and thermal conductivity of AMLs.

36 MATERIALS SCIENCE↗

Extreme sensitivity of the magnetic ground state to halide composition in FeCl 3 – x Br x

Mixed halide chemistry has recently been utilized to tune the intrinsic magnetic properties of transition-metal halides—one of the largest families of magnetic van der Waals materials. Prior studies have shown that the strength of exchange interactions, hence the critical temperature, can be tuned smoothly with halide composition for a given ground state. Here we show that the ground state itself can be altered by a small change of halide composition in FeCl 3–$x$ Br $x$ . Specifically, we find a threefold jump in the Néel temperature and a sign change in the Weiss temperature at $x$=0.08 corresponding to only 3% bromine doping. Using neutron scattering, we reveal a change of the ground state from spiral order in FeCl 3 to $A$-type antiferromagnetic order in FeBr 3 . From first-principles calculations, we show that a delicate balance between nearest and next-nearest neighbor interactions is responsible for such a transition. In conclusion, these results demonstrate how varying the halide composition can tune the competing interactions and change the ground state of a spiral spin liquid system.

2-dimensional systems↗

Tuning the spin dynamics and magnetic phase transitions of the Cantor alloy via composition and sample processing protocols: A muon spin relaxation study

CrMnFeCoNi, also called the Cantor alloy, is a well-known high-entropy alloy whose magnetic properties have recently become a focus of attention. Here, we present a detailed muon spin relaxation study of the influence of chemical composition and sample processing protocols on the magnetic phase transitions and spin dynamics of several different Cantor alloy samples. Specific samples studied include a pristine equiatomic sample, samples with deficient and excess Mn content, and equiatomic samples magnetized in a field of 9 T or plastically deformed in pressures up to 0.5 GPa. The results confirm the sensitive dependence of the transition temperature on composition and demonstrate that post-synthesis pressure treatments cause the transition to become significantly less homogeneous throughout the sample volume. In addition, we observe critical spin dynamics in the vicinity of the transition in all samples, reminiscent of canonical spin glasses and magnetic materials with ideal continuous phase transitions. Application of an external magnetic field suppresses the critical dynamics in the Mn-deficient sample, while the equiatomic and Mn-rich samples show more robust critical dynamics. The spin-flip thermal activation energy in the paramagnetic phase increases with Mn content, ranging from 3.1⁢(3) × 10 -21 J for 0% Mn to 1.2⁢(2) × 10 -20 J for 30% Mn content. These results shed light on critical magnetic behavior in environments of extreme chemical disorder and demonstrate the tunability of spin dynamics in the Cantor alloy via chemical composition and sample processing.

36 MATERIALS SCIENCE↗

Nepheline crystallization and the residual glass composition: Understanding waste glass durability

The vast majority of High-Level Waste (HLW) originating from defense nuclear programs is sequestered and immobilized in borosilicate glass. Borosilicate glass is universally accepted for immobilizing HLW, but its efficiency has limitations based on the compositional makeup of the waste stream. The chemical durability of the glass is the most important factor in determining the longevity and usefulness of the final glass waste form. The primary detriment to this durability in glasses containing high levels of aluminum is nepheline (NaAlSiO 4 ) crystallization, as it is generally accompanied by a measurable decrease in the glass's chemical durability. This work seeks to understand nepheline crystallization, within the context of thermal history, and to elucidate the influence of compositional shifts in the residual glass (after crystallization) on the measured durability. The results presented within show a distinct deviation in leaching behavior as a function of structural makeup (calculated Q units). This understanding will provide practical information required for broadening glass compositional regions needed to more efficiently vitrify HLW.

36 MATERIALS SCIENCE↗

Phase transformation and electrochemical charge storage properties of vanadium oxide/carbon composite electrodes synthesized via integration with dopamine

Abstract Chemically preintercalated dopamine (DOPA) molecules were used as both a reducing agent and a carbon precursor to prepare δ‐V 2 O 5 · n H 2 O/C, H 2 V 3 O 8 /C, VO 2 (B)/C, and V 2 O 3 /C nanocomposites via hydrothermal treatment or hydrothermal treatment followed by annealing under Ar flow. We found that the phase composition and morphology of the produced composites are influenced by the DOPA:V 2 O 5 ratio used to synthesize (DOPA) x V 2 O 5 precursors through DOPA diffusion into the interlayer region of the δ‐V 2 O 5 · n H 2 O framework. The increase of DOPA concentration in the reaction mixture led to a more pronounced reduction of vanadium and a higher fraction of carbon in the composites’ structure, as evidenced by X‐ray photoelectron spectroscopy and Raman spectroscopy measurements. The electrochemical charge storage properties of the synthesized nanocomposites were evaluated in Li‐ion cells with nonaqueous electrolytes. δ‐V 2 O 5 · n H 2 O/C, H 2 V 3 O 8 /C, VO 2 (B)/C, and V 2 O 3 /C electrodes delivered high initial capacities of 214, 252, 279, and 637 mAh g –1 , respectively. The insights provided by this investigation open up the possibility of creating new nanocomposite oxide/carbon electrodes for a variety of applications, such as energy storage, sensing, and electrochromic devices.

Materials Science↗

Insight of BaCe 0.5 Fe 0.5 O 3– δ twin perovskite oxide composite for solid oxide electrochemical cells

One-pot synthesized twin perovskite oxide composite of BaCe 0.5 Fe 0.5 O 3–δ (BCF), comprising cubic and orthorhombic perovskite phases, shows triple-conducting properties for promising solid oxide electrochemical cells. Phase composition evolution of BCF under various conditions was systematically investigated, revealing that the cubic perovskite phase could be fully/partially reduced into the orthorhombic phase under certain conditions. The reduction happened between the two phases at the interface, leading to the microstructure change. As a result, the corresponding apparent conducting properties also changed due to the difference between predominant conduction properties for each phase. Based on the revealed phase composition, microstructure, and electrochemical properties changes, a deep understanding of BCF's application in different conditions (oxidizing atmospheres, reducing/oxidizing gradients, cathodic conditions, and anodic conditions) was achieved. Triple-conducting property (H + /O 2– /e – ), fast open-circuit voltage response (~16–~470 mV) for gradients change, and improved single-cell performance (~31% lower polarization resistance at 600°C) were comprehensively demonstrated. Besides, the performance was analyzed under anodic conditions, which showed that the microstructure and phase change significantly affected the anodic behavior.

08 HYDROGEN↗

Compositional evidence for chondrule origins of low‐Ca pyroxenes in comet Wild 2 and a giant cluster IDP

Abstract A literature compilation of 1136 low‐Ca pyroxene compositions from chondrules from 12 primitive type 2–3 carbonaceous, ordinary and enstatite chondrite groups define unique regions on an Al 2 O 3 and Cr 2 O 3 diagram when compared to low‐Ca pyroxenes from equilibrated type 4‐6 chondrites. Measured compositions of 100 low‐Ca pyroxenes from comet Wild 2 and a giant cluster IDP of probable cometary origin are similar to each other and fall in the type 2–3 chondrite chondrule region suggesting that most of the pyroxenes likely formed in the solar nebula like conventional chondrules. The data imply that most low Ca‐pyroxenes from comet Wild 2 and the giant cluster IDP formed from igneous crystallization processes and did not experience significant thermal metamorphism, indicating that the low‐Ca pyroxenes were unlikely incorporated into large parent bodies prior to accretion in their respective comet bodies. An intriguing group of nine low‐Ca pyroxenes from comet Wild 2 with low Cr and Al that fall where type 4–6 chondrites are located are interpreted as products of condensation. The compositional data combined with previously measured oxygen isotopes on 17 low‐Ca pyroxenes support earlier conclusions that comet samples have links with carbonaceous, ordinary, and possibly enstatite chondrite groups. Our results provide additional evidence that comets accreted materials from multiple chondrule reservoirs throughout the solar nebula.

Geochemistry & Geophysics↗

Three-Dimensional Permeability of Thick-Section Glass Fabric Reinforced Polymer Composite by Vacuum-Assisted Resin Infusion Molding

Abstract Determination of permeability of thick-section glass fabric preforms with fabric layers of different architectures is critical for manufacturing large, thick composite structures with complex geometry, such as wind turbine blades. The thick-section reinforcement permeability is inherently three-dimensional and needs to be obtained for accurate composite processing modeling and analysis. Numerical simulation of the liquid stage of vacuum-assisted resin infusion molding (VARIM) is important to advance the composite manufacturing process and reduce processing-induced defects. In this research, the 3D permeability of thick-section E-glass fabric reinforcement preforms is determined, and the results are validated by a comparison between flow front progressions from experiments and from numerical simulations using ansys fluent software. The orientation of the principal permeability axes were unknown prior to experiments. The approach used in this research differs from those in literature in that the through-thickness permeability is determined as a function of flow front positions along the principal axes and the in-plane permeabilities and is not dependent on the inlet radius. The approach was tested on reinforcements with fabric architectures which vary through-the-thickness direction, such as those in a spar cap of a wind turbine blade. The computational simulations of the flow-front progression through-the-thickness were consistent with experimental observations.

Engineering↗

A Subspace-Inclusive Sampling Method for the Computational Design of Compositionally Graded Alloys

Abstract Compositionally graded alloys, a subclass of functionally graded materials (FGMs), utilize localized variations in composition with a single metal part to achieve higher performance than traditional single material parts. In previous work [Kirk, T., Galvan, E., Malak, R., and Arroyave, R., 2018, “Computational Design of Gradient Paths in Additively Manufactured Functionally Graded Materials,” J. Mech. Des., 140, p. 111410. 10.1115/1.4040816], the authors presented a computational design methodology that avoids common issues which limit a gradient alloy’s feasibility, such as deleterious phases, and optimizes for performance objectives. However, the previous methodology only samples the interior of a composition space, meaning designed gradients must include all elements in the space throughout the gradient. Because even small amounts of additional alloying elements can introduce new deleterious phases, this characteristic often neglects potentially simpler solutions to otherwise unsolvable problems and, consequently, discourages the addition of new elements to the state space. The present work improves upon the previous methodology by introducing a sampling method that includes subspaces with fewer elements in the design search. The new method samples within an artificially expanded form of the state space and projects samples outside the true region to the nearest true subspace. This method is evaluated first by observing the sample distribution in each subspace of a 3D, 4D, and 5D state space. Next, a parametric study in a synthetic 3D problem compares the performance of the new sampling scheme to the previous methodology. Lastly, the updated methodology is applied to design a gradient from stainless steel to equiatomic NiTi that has practical uses such as embedded shape memory actuation and for which the previous methodology fails to find a feasible path.

Engineering↗

A FeCrAl-Al2O3 Composite Produced via Laser Powder Bed Fusion of a Mixed Powder for Porous Catalyst Scaffolds

This study proposes a novel approach for synthesizing and etching bicontinuous FeCrAl-Al2O3 composites as a means for replacing FeCrAl foams as catalyst scaffolds in bio-driven alcohol reactors for jet-fuel production. Conventional FeCrAl foams suffer from poor availability and consequent high costs. New additive manufacturing techniques provide an opportunity to produce tailored foams in reasonable times and at acceptable costs. This research aimed to generate a porous FeCrAl structure by etching a bicontinuous FeCrAl-Al2O3 composite produced by laser powder bed fusion of amalgamated FeCrAl and Al2O3 powders. The composite powder for laser powder bed fusion is created by ball-milling FeCrAl and Al2O3 powders. This research focuses on achieving a bi-continuous FeCrAl-Al2O3 structure, essential for the selective removal of the ceramic phase. The influence of laser processing parameters on the microstructure was examined across a range of laser powers (60-120 W) and scan speeds (100-400 mm/s), showing that higher powers and speeds produce finer metal struts. A bi-continuous microstructure was consistently obtained, marking a key achievement. The Al2O3 removal process involved a two-step etching method using hydrochloric and phosphoric acids, tested across various etching times. The alumina phase was reduced from 36 vol% to 17 vol% (corresponding to an increase in porosity from 24 vol% to 43 vol%), showing the potential for use as a porous catalyst scaffold. This research demonstrates the potential for using additive manufacturing to produce porous FeCrAl structures capable of replacing hard-to-source FeCrAl foams.

Son, Kwangtae↗

Analytical Modeling of Metal Foam Composite Phase Change Materials (PCM) in Thermal Energy Storage Using Asymptotic Analysis

The use of phase change materials (PCMs) for thermal energy storage can release or absorb a significant amount of latent heat during the freezing or melting process, offering a higher energy storage density. One of the main drawbacks of PCMs is their low thermal conductivity, resulting in poor thermal performance. Recent research has attempted to enhance heat transfer and increase the thermal conductivity of PCMs, including the use of metal foams. However, modeling the metal foam composite PCM using conventional methods is computationally expensive. This paper proposes an asymptotic solution for a Stefan-like problem subject to a convective boundary for outward solidification in a hollow cylinder, capable of predicting the freeze-melt cycle of the metal foam composite PCM. Specifically, three temporal regimes and four spatial layers are considered in the asymptotic analysis for each phase change process. The thermal conductivity is calculated by a theoretical three-dimensional tetrakaidecahedron model, while other thermophysical properties are obtained using the method of volume averaging. The results are verified with numerical data and validated against experimental data in the literature. The presented analytical modeling framework could have the potential to be applied to other types of composite PCMs with considerably lower computational costs compared with conventional methods.

analytical model↗

Magnetic phase diagram mapping in Fe 1- x Rh x composition-spread thin films

We have fabricated high-quality polycrystalline Fe 1- x Rh x composition-spread thin films by cosputtering Fe and Rh, and investigated their structural and magnetic transformations as a systematic function of composition. With increasing Rh concentration, Fe 1-x Rh x thin film undergoes from an α' phase to a disordered γ phase and also shows a magnetic transition from a ferromagnetic phase to a paramagnetic phase. Vibrating-sample magnetometry and x-ray magnetic circular dichroism measurements show an antiferromagnetic-ferromagnetic transition in the range of 0.52 < x < 0.58 in the Fe 1- x Rh x composition gradient at room temperature. Based on our structural and magnetic property mapping, we construct a thin-film phase diagram of Fe 1- x Rh x . Compared to reported results in bulk alloys, the antiferromagnetic-ferromagnetic transition in the Fe 1- x Rh x thin films was found to occur at slightly higher Rh concentrations, while the boundary between the pure γ phase and the α'/ γ mixed phase region is shifted to the lower concentration Rh.

36 MATERIALS SCIENCE↗

Composite entanglement topology and extensional rheology of symmetric ring-linear polymer blends

Extensive molecular simulations are applied to characterize the equilibrium dynamics, entanglement topology, and nonlinear extensional rheology of symmetric ring-linear polymer blends with systematically varied ring fraction ϕR. Chains with degree of entanglement Z≈14 are mixed to produce 10 well-entangled systems with ϕR varying from neat linear to neat ring melts. Z is large enough that except for very large ϕR, the rings are threaded by multiple linear chains in equilibrium. Primitive path analysis is used to visualize and quantify the structure of the composite ring-linear entanglement network. We measure the quantity of ring-linear threading and linear-linear entanglement as a function of ϕR and identify with simple arguments a ring fraction ϕR≈0.4 where the topological constraints of the composite entanglement network are maximized. These topological analyses are used to rationalize the ϕR-dependence of ring and linear chain dynamics, conformations, and rheology. Simulations of startup uniaxial elongation flows demonstrate the extensional stress overshoot observed in recent filament stretching experiments and characterize how it depends on the blend composition and entanglement topology. The overshoot is driven by an overstretching and recoil of ring polymers due to the convective unthreading of rings from linear chains.

O’Connor, Thomas C.↗

Multidecadal shifts in forest plant diversity and community composition across glacial landforms in northern lower Michigan, USA

Understanding how plant community assemblage is affected by spatial and temporal patterns is crucial to understanding forest ecosystem responses to disturbance, including future climate change. In this article, we tracked how diversity and composition are distributed through space and time in a midsuccessional mixed hardwood forest in northern lower Michigan, United States. This region’s geographically and abiotically distinct glacial landforms influence both the spatial and temporal dynamics of its forest communities. Vegetation sampling plots (n = 87) were established at the University of Michigan Biological Station in 1990 and resampled in 2015. Vegetation in the overstory, sapling, and groundcover layers was censused. Abiotic variables, including elevation, pH, and soil nutrients, were measured in a subset of plots (n = 40). There were strong differences in diversity and community composition among glacial landforms, with the moraine having a 31% greater species richness in the groundcover layer compared with those of the other glacial landforms. Surprisingly, plant communities across all three vegetation layers showed little change over the 25-year period, and we found no evidence of differences in successional rates among glacial landforms. Our findings indicate that glacial landforms have a large influence on the production and maintenance of local plant diversity and community composition in this area and suggest that successional dynamics may manifest themselves over much longer time periods in these northern biomes.

Forestry↗

Computational Modeling of Heterogeneity of Stress, Charge, and Cyclic Damage in Composite Electrodes of Li-Ion Batteries

Charge heterogeneity is a prevalent feature in many electrochemical systems. In a commercial cathode of Li-ion batteries, the composite is hierarchically structured across multiple length scales including the sub-micron single-crystal primary-particle domains up to the macroscopic particle ensembles. The redox kinetics of charge transfer and mass transport strongly couples with mechanical stresses. This interplay catalyzes substantial heterogeneity in the charge (re)distribution, stresses, and mechanical damage in the composite electrode during charging and discharging. We assess the heterogeneous electrochemistry and mechanics in a LiNi x Mn y Co z O 2 (NMC) cathode using a fully coupled electro-chemo-mechanics model at the cell level. A microstructure-resolved model is constructed based on the synchrotron X-ray tomography data. We calculate the stress field in the composite and then quantitatively evaluate the kinetics of surface charge transfer and Li transport biased by mechanical stresses. We further model the cyclic behavior of the cell. The repetitive deformation of the active particles and the weakening of the interfacial strength cause gradual increase of the interfacial debonding. The mechanical damage impedes electron transfer, incurs more charge heterogeneity, and results in the capacity degradation in batteries over cycles.

25 ENERGY STORAGE↗

Improved Fuel Cell Chemical Durability of an Heteropoly Acid Functionalized Perfluorinated Terpolymer-Perfluorosulfonic Acid Composite Membrane

Commercial proton exchange membrane heavy-duty fuel cell vehicles will require a five-fold increase in durability compared to current state-of-the art light-duty fuel cell vehicles. We describe a new composite membrane that incorporates silicotungstic heteroply acid (HPA), α -K 8 SiW 11 O 40 ▪13H 2 O, a radical decomposition catalyst and when acid-exchanged can potentially conduct protons. The HPA was covalently bound to a terpolymer of tetrafluoroethylene, vinylidene fluoride, and sulfonyl fluoride containing monomer (1,1,2,2,3,3,4,4-octafluoro-4-((1,2,2-trifluorovinyl)oxy)butane-1-sulfonyl fluoride) by dehydrofluorination followed by addition of diethyl (4-hydroxyphenyl) phosphonate, giving a perfluorosulfonic acid-vinylidene fluoride-heteropoly acid (PFSA-VDF-HPA). A composite membrane was fabricated using a blend of the PFSA-VDF-HPA and the 800EW 3M perfluoro sulfonic acid polymer. The bottom liner-side of the membrane tended to have a higher proportion of HPA moieties compared to the air-side as gravity caused the higher mass density PFSA-VDF-HPA to settle. The composite membrane was shown to have less swelling, more hydrophobic properties, and higher crystallinity than the pure PFSA membrane. The proton conductivity of the membrane was 0.130 ± 0.03 S cm −1 at 80 °C and 95% RH. Impressively, when the membrane with HPA-rich side was facing the anode, the membrane survived more than 800 h under accelerated stress test conditions of open-circuit voltage, 90 °C and 30% RH.

08 HYDROGEN↗

Modeling Single and Two-Phase Transport in Thin Porous Layers Using a Composite Continuum-Pore Network Formulation

In this work, a composite continuum-pore network formulation is presented to model single and two-phase transport thin porous layers, such as gas diffusion layers in polymer electrolyte fuel cells (PEFCs) and active electrodes in redox flow batteries (RFBs). The formulation can be integrated into CFD codes, thus combining the ease of implementation of continuum-based modeling and the computational power of pore-network modeling. The composite model includes a control volume (CV) mesh at the layer scale, which embeds a cubic pore network [1,2]. The pore-network model is used to determine analytically local anisotropic effective transport properties (local effective diffusivity and permeability), which are mapped onto the CV mesh to simulate transport in the porous transport layer [3,4]. Good agreement is found between the predicted global effective transport properties (global effective diffusivity and permeability) under dry and wet conditions and previous experimental data reported in the literature for Toray TGP-H series carbon paper. Water saturation distributions are also compared with results obtained using X-ray computed tomography [4]. [1] P.A. García-Salaberri, I.V. Zenyuk, J.T. Gostick, A.Z. Weber, Modeling Gas Diffusion Layer in Polymer Electrolyte Fuel Cells Using a Continuum-Based Pore-Network Formulation, ECS Trans. 97 (2020) 615. [2] P.A. García-Salaberri, Modeling diffusion and convection in thin porous transport layers using a composite continuum-network model: Application to gas diffusion layers in polymer electrolyte fuel cells, Int J. Heat Mass Transf. (2020), submitted. [3] P.A. García-Salaberri, J.T. Gostick, G. Hwang, A.Z. Weber, M. Vera, Effective diffusivity in partially-saturated carbon-fiber gas diffusion layers: Effect of local saturation and application to macroscopic continuum models, J. Power Sources 296 (2015) 440–453. [4] P.A. García-Salaberri, G. Hwang, M. Vera, A.Z. Weber, J.T. Gostick, Effective diffusivity in partially-saturated carbon-fiber gas diffusion layers: Effect of through-plane saturation distribution, Int. J. Heat Mass Transf. 86 (2015) 319–333.

Garcia-Salaberri, Pablo↗