Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Micro”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Enabling high rate capability, low internal resistance, and excellent cyclability for vanadium redox flow batteries utilizing ultrafast laser-structured graphite felt

The electrochemical performance of vanadium redox flow batteries (VRFBs) was enhanced via laser-patterned graphite felt (GF) electrodes. The laser-structured GF electrodes engineered via preparing a series of well-ordered microscopic channel structures with an average width of 200 μm, creating a three-dimensional carbon framework. The ultrafast laser patterning increased the porosity by 10%, as compared to pristine electrode. The analysis of the overpotential distribution using electrochemical impedance spectroscopy revealed that the electrode polarization involving charge transfer and diffusion resistance is strongly alleviated with the aid of carbon micro-perforation. The advanced design of laser-structured GF electrode also exhibits high rate capability, low internal resistance, and excellent cyclability. The improved performance was attributed to the synergistic effect involving (i) high electrochemically active surface area for rapid electrochemical reactions (i.e., V(II)/V(III) and V(IV)/V(V) redox couples) and (ii) excellent transport properties for facile electron/ion/species transport. The micro-scale channels created with the laser-ablation technique act as capillary structures and enabled homogeneous and rapid wetting of GF electrodes as formulated with the classical Washburn equation. The VRFB equipped with the laser-structured GF electrodes delivered a high discharge capacity with exceptional capacity retention upon extended cycling (>84.9%). In conclusion, accordingly, the design of laser-structured electrode paves the pathway towards finely tuning the electrode internal microstructure for improved cyclability, durability, and capacity retention in various electrochemical energy storage/conversion devices (e.g. all-vanadium redox flow batteries).

25 ENERGY STORAGE↗
DOE OSTI2025

Effect of cell compression on the performance and the structure of proton exchange membrane water electrolyzer (PEMWE) assembly

Here, in the field of water electrolysis, the proton exchange membrane water electrolyzer (PEMWE) is currently the most advanced technology for producing hydrogen without emitting CO 2 . Although PEMWE plants are already in operation, further research is needed to improve cell efficiency and reduce the use of rare materials, such as iridium oxide catalysts for the oxygen evolution reaction (OER). One of the main causes of performance loss in PEMWE is the relatively low electric conductivity of the porous transport layer (PTL) and of the anode catalyst layer, which results in ohmic losses and low catalyst utilization during high current density operation. The objective of this study is to investigate how optimization of the PTL and electrode interface can increase the cell performance. To this end, we tested different cell assemblies using fibrous and sintered PTLs, decreasing membrane thickness, reducing iridium loading, and inserting a microporous layer to increase contact surface area. Electrochemical characterization of each cell configuration was systematically performed at various compression levels as the pressure is a crucial parameter influencing the electrode/PTL contact area. In parallel, X-ray microcomputed tomography (micro-CT) was employed to investigate the effects of cell hydration and compression on the structure of PEMWE components. This study combining electrochemistry and micro-CT imaging presents how optimizing the electrode/PTL contact surface area, minimizes ohmic losses, and enables PEMWE operation with low iridium loading at high current densities.

Catalyst - PTL interface↗
DOE OSTI2021

Dynamic in situ imaging of methane hydrate formation in coal media

Fast 3D X-ray imaging has proved to provide crucial insights into multi-phase dynamic processes in various geomaterials. In particular, many in situ imaging experiments have been made to study gas-hydrate formation in porous sandy samples. Such imaging is challenging for the methane gas-hydrate formation in coal samples because of the coal micro-porosity structure and lower X-ray contrast. Here we present results of the first dynamic in situ micro-computed tomography experiment of methane-hydrate formation in coal samples. Synchrotron phase-contrast tomography techniques allowed to achieve necessary contrast levels to separate all the materials of interest (gas, water, coal, and gas hydrate) in reconstructed images with high spatial and temporal resolution. The imaging results are compared to the ones from a similar tomographic experiment with sand samples. Methane-hydrate formation is accompanied by the water movements caused by cryogenic water suction that happens in sequences of short fast movements with longer equilibrium states in between. Only one type of the hydrate formation was observed in coal (growth as shells on grain boundaries) as opposed to three types in sand (shells, growth into gas pockets, and inside water volumes). In particular, this leads to a slower hydrateformation speed in coal. For the coal sample, we also observed water extraction from grains, and interpreted it as competitive sorption of methane. Lastly, we visualized the dynamic behavior of this water extraction via microchannels inside the coal grains, and performed nano-tomography imaging of these channels for a better understanding of this phenomenon.

58 GEOSCIENCES↗
DOE OSTI2023

Assessing pore network heterogeneity across multiple scales to inform CO2 injection models

Geologic heterogeneity is a key feature that must be considered when translations of scaled data are performed. This paper presents the assessment of geologic heterogeneity using a multiscale workflow that includes image analysis-based methods coupled with well log analysis to provide data in which fractals and machine learning methods estimate the carbon dioxide (CO 2 ) storage resource potential of a reservoir. The heterogeneity of rock properties of the complex Bell Creek reservoir in Montana, USA, was explored at the pore scale (~nm to mm), core scale (~mm to m), and well scale (~cm to m). The data used in this study included advanced image analysis of micro-CT (computed tomography) images (pore scale), thin sections (pore scale), plugs and core images (core scale) and well logs (well scale). The micro-CT images were segmented using a U-net segmentation approach into objects of pores and grains. Further, the segmented images were reconstructed into subvolumes of different sizes. Physical properties (porosity and permeability) and fractal dimensions were calculated for the various subvolumes, and Lorenz coefficient (Lc) values, a single parameter to describe the degree of heterogeneity within a pay zone section, were calculated from thin-section images and well logs. Porosity and fractal dimension values were used to estimate the 188-µm threshold of representative elementary volume (REV) in this study. Both the Lc and fractal dimension values were found to be negatively correlated. When these two parameters are combined, it is possible to discern differences in the complex porous networks of the samples analyzed in this study.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

A remarkable CHF of 345W/cm 2 is achieved in a wicked-microchannel using HFE-7100

Boiling heat transfer of dielectric fluids is a promising cooling technique for thermal management of microelectronic systems. However, the critical heat flux achievable is generally low, because of the poor thermophysical properties of these fluids. Here, to address this dilemma, we propose a new cooling concept to substantially enhance global liquid supply during phase-change process using enhanced capillary-driven force. Additionally, dedicated vapor pathways are designed among a bank of micro-pillars to facilitate vapor removal. In this work, new wicks comprised of silicon micro-pinfin arrays are explored to significantly enhance the flow boiling heat transfer performance. To examine the functionalities of this wick, experiments on HFE-7100 were carried out with mass velocities varying from 247 to 3,465 kg/m 2 s. To explore the enhancement mechanisms and to analyze the capillary-assisted flow boiling process, visualization studies were conducted. The results indicate that sustainable evaporation induced by wick microstructures and efficient liquid supply are the enhancement mechanisms compared to parallel microchannels with solid walls. It is found that the overall heat transfer coefficient is substantially increased up to 75%. Remarkably, a high critical heat flux (CHF) of approximately 345 W/cm 2 is recorded at G = 3,465 kg/m 2 s at coolant inlet temperature of ~20 °C. Equally importantly, this noticeable enhancement of CHF value is associated with drastically decreased pressure drops compared to microchannels decorated with μ-pinfin fences.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Experimental investigation of multi-scale strain-field heterogeneity in rocks

Rocks have inherently heterogeneous micro-structures. The mechanical response of these heterogeneous micro-structures primarily governs the overall macroscopic behavior of intact rocks, as these features cause heterogeneity in the distribution of stress and strain within a rock volume. These attributes have generally been studied through numerical tools, but in this study, the stress-induced multi-scale spatial heterogeneity in the strain-field of three types of rocks was analyzed experimentally. The full-field non-contact strain measurement approach of 2D Digital Image Correlation (DIC) was employed for real-time evaluation of the progression of strain-field heterogeneity across the surface of the rock specimens when subjected to uniaxial compression. The results show that the strain-field heterogeneity is a characteristic of the inherent heterogeneity of the rock specimens, and that it amplifies with the magnitude of stress applied on the rock specimens. The changes in the strain-field heterogeneity with loading were consistent with established microcrack evolution processes in the rock specimens. The statistical fluctuation in the strain-field heterogeneity of the three rocks at varying gauge lengths (GL) was then analyzed for the determination of the Representative Volume Element (RVE) length-scale at increasing load levels. The evolution of the RVE length-scale with loading facilitates an objective way of determination of the stress at which the strain-field heterogeneity increases and begins to influence the overall global response of the material. The RVE length-scale variation with loading showed that at the RVE length-scale reaches the specimen-scale at load levels well below the peak strength.

2D-digital image correlation↗
DOE OSTI2020

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↗
DOE OSTI2021

Iron-rich microstructure records of high temperature multi-component silicate melt behavior in nuclear fallout

Above-ground nuclear explosions that interact with the surface of the earth entrain materials from the surrounding environment, influencing the resulting physical and chemical evolution of the fireball, which can affect the final chemical phase and mobility of hazardous radionuclides that are dispersed in the environment as fallout particles. The interaction of iron with a nuclear explosion is of specific interest due to the potential for iron to act as a redox buffer and because of the likelihood of significant masses of metals to be present in urban environments. We investigated fallout from a historic surface interacting nuclear explosion conducted on a steel tower and report the discovery of widespread and diverse iron-rich micro-structures preserved within the samples, including crystalline dendrites and micron-scale iron-rich spheres with liquid immiscibility textures. We assert these micro-structures reflect local redox conditions and cooling rates and can inform interpretation of high temperature events, enabling new insights into fireball condensation physics and chemistry when metals from the local environment (i.e. structural steel) are vaporized or entrained. Finally, these observations also significantly expand the availability of silicate immiscibility datasets applicable to rapidly quenched systems such as meteorite impact melt glass.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2025

Spatial mapping of the localized corrosion behavior of a magnesium alloy AZ31B tungsten inert gas weld

Sections of a magnesium alloy, AZ31B, joined with tungsten inert gas (TIG) welding, were examined with scanning electrochemical microscopy (SECM) and scanning Kelvin probe force microscopy (SKPFM) to investigate corrosion mechanisms by correlating observed corrosion behavior with weld-affected microstructural variations. Insight into the changing nature of the galvanic couples between weld zones and at localized microgalvanic sites were investigated using SECM and SKPFM to map both electrochemically active regions and Volta potential differences across the weld-affected zones. The formation of an Al-Zn solidification network in the fusion zone (FZ) at and near the TIG weld epicenter differs from the outer heat-affected zone (HAZ), where intermetallic particles (IMPs) are the notable secondary phase from the magnesium matrix. These microstructures were mapped with SKPFM before and after brief exposure to a salt solution, revealing micro-galvanic couples as the main driving force to corrosion initiation and propagation within each zone. The IMPs and Al-Zn solidification network act as strong cathodes and govern the corrosion processes. The galvanic coupling and evolution of the intrinsic corrosion behavior between the weld zones is explained by monitoring the hydrogen evolution reaction (HER) with SECM over time. Anodically induced cathodic activation is confirmed for this welded material, as micro-galvanic couples between microstructural features are found to transition over time to broad electrochemically active areas within the weld-affected zones, resulting in polarity reversal as time of exposure proceeds.

Anodically induced cathodic activation↗
DOE OSTI2023

Processing aluminum alloy with hybrid wire arc additive manufacturing and ultrasonic nanocrystalline surface modification to improve porosity, surface finish, and hardness

Here, in this paper, a novel hybrid wire arc additive manufacturing (WAAM) and ultrasonic nanocrystal surface modification (UNSM) on porosity manipulation and surface properties of aluminum 5356 alloys was studied. The goal is to improve the quality of the WAAM-built part by eliminating bigger pores and reducing its size, reducing surface roughness, and increasing surface hardness. The as-built WAAM and WAAM-UNSM-treated samples were quantitatively studied for porosity using an X-ray micro-computed tomography (μ-CT). The surface roughness was measured on the surface profile of the same samples before and after UNSM treatment. Followed by the Vickers micro-hardness tests to evaluate the hardness modified by the influence of the UNSM treatment. It was found that the bigger pores in the as-built WAAM samples were eliminated and the medium-sized pores were shrunk to almost half the size after the UNSM treatment. Further, the UNSM treatment showed a significant improvement in both surface roughness and hardness on the WAAM Al5356 samples. This experimental work demonstrates the critical advantages of hybrid WAAM-UNSM in improving the qualities of the WAAM processed parts.

Hardness↗
DOE OSTI2020

Quantitative and qualitative bone imaging: A review of synchrotron radiation microtomography analysis in bone research

All levels of the unique hierarchical structure of bone, consisting of collagen and hydroxyapatite crystals at the nanoscale to osteon/lamellae structures at the microscale, contribute to its characteristic toughness and material properties. Elements of bone's density and size contribute to bone quantity (or bone mass), whereas elements of bone's material composition, material properties, internal structure, and organization describe bone quality. Furthermore, bone quantity and quality can be degraded by factors such as aging, disease, treatments, and irradiation, compromising its ability to resist fracture and sustain loading. Accessing the morphology and architecture of bone at the microscale to quantify microstructural features and assess the degree of mineralization and path of crack propagation in bone provides crucial information on how these factors are influencing bone quantity and quality. Synchrotron radiation micro-computed tomography (SRμCT) was first used to assess bone structure at the end of the 1990's. One of the main advantages of the technique is that it enables accurate three-dimensional (3D), non-destructive quantification of structure while traditional histomorphometry on histological sections is inherantly destructive to the sample and two-dimensional (2D). Additionally, SRμCT uses monochromatic, high-flux X-ray beams to provide high-resolution and high-contrast imaging of bone samples. This allows the quantification of small microstructural features (e.g. osteocyte lacunae, canals, trabeculae, microcracks) and direct gray value compositional mapping (e.g. mineral quantification, cement lines) with greater speed and fidelity than lab-based micro-computed tomography. In this article, we review how SRμCT has been applied to bone research to elucidate the mechanisms by which bone aging, disease, and other factors affect bone fragility and resistance to fracture.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Pore-scale simulation of gas flow in microscopic permeable media with complex geometries

Lattice Boltzmann method (LBM) is an efficient tool to perform direct numerical simulation of gas flow in micro-size pores. Despite great advances in LBM and its use in describing gas flow in micro-size channels in slip and transitional flow regimes, exploring formulations that are able to capture flow behavior in domains with complex boundary geometries remains a challenging task. As a result, the impact of complexities in pore structures on gas flow may not be fully explored using LBM. In this study, we propose to use the combined bounce-back and Maxwellian diffusive reflection scheme to capture the slip velocity with a new set of slip coefficients. Optimal values of slip coefficients are determined using a gradient based method, where the results from the linearized Boltzmann equation is used as the reference solution. The proposed formulation is further validated against predictions from molecular dynamics simulation in the presence of complex geometries by introducing obstacles of different shapes in straight channels. LBM with the proposed boundary treatment is then used to investigate gas flow in a synthetic microscopic permeable medium. Here, results indicate that gas exhibits different flow configurations as Knudsen number varies, and apparent gas permeability appears to have up to a second-order dependency on the reciprocal mean pressure in slip and transitional flow regimes.

03 NATURAL GAS↗
DOE OSTI2021

Analysis of time-dependent strain heterogeneity in shales using X-ray microscopy and digital volume correlation

In this paper, we report the results of a 4-dimensional (3 spatial dimensions plus time) in-situ creep experiment on a milli-core shale plug from Barnett formation in Texas. An outstanding challenge in shale mechanics is connecting observable characteristics to predictable mechanical and hydraulic behavior. We employ novel imaging and analysis techniques to statistically link sample mineralogy and structure to observed deformation. In particular, we focus on mechanisms of porosity reduction during transient response to step changes in stress, analogous to the aseismic deformation that occurs after an episode of fracturing. Our main goal in this study is to take the first steps to visualize and capture the small-scale deformation mechanisms involved in creep of shales, using imaging techniques. To do so, we conducted a uniaxial creep experiment inside an Xradia Versa 520 micro-CT system for 12 h at room temperature under a constant uniaxial stress of 25 MPa. We employed a Digital Volume Correlation (DVC) method to measure micro-scale deforspmation of the sample by comparing the resulting volumetric images with a high-resolution scan of the sample obtained prior to the creep experiment. The results indicate that the time-dependent deformation in the compliant clay plus kerogen matrix constituents facilitates rigid grain rotation and compaction of intergranular porosity at the boundary between rigid grains and matrix.

58 GEOSCIENCES↗
DOE OSTI2022

Intense cyclic heating effects on thermo-fracture and thermal shock of solid tungsten and open-cell tungsten foam

Here we investigate here the effects of transient (cyclic) arc-jet plasma and laser heating on fracture behavior of W-foam and solid tungsten. The two key parameters that control the foam thermomechanical response are its density and mean cell size expressed in Pores Per Inch (PPI). Tungsten foam samples were fabricated with Chemical Vapor Deposition (CVD) with variety of PPI and relative density. These were tested under two types of qualitatively different conditions: (1) high-enthalpy arc-jet, and (2) high-power cyclic laser heating. None of the foam samples showed macroscopic through-thickness cracks. However, distributed micro-cracks were observed on ligaments and their triple junctions. Under the same loading conditions, W-foam and solid tungsten showed similar crack network pattern and characteristic length-scale. However, Crack Opening Displacement (COD) was twice as large in solid W as compared to W-foam. Foam samples that have been previously exposed to a low-pressure helium plasma showed significant changes in their surface forming nano-texture fuzz which was removed by subsequent testing in the arc-jet. Extensive fracture and re-crystallization were observed in the thin solid W disk that was fully-constrained from expansion. Thicker and fully-constrained solid W disks did not display recrystallization, grain growth, and extensive cracking. However, thicker disks that were free to expand showed some recrystallization and extensive through-thickness cracks due to less effective cooling and thus higher temperatures. Laser beam testing showed no visible damage formation at 0.19 GW/m 2 and 0.38 GW/m 2 for both low-density (23%) and high-density (43%) foams at low pulses (100-1000). Micro-cracks were observed after 10,000 pulses at 0.19 GW/m 2 in both foams, and in low-density foam after 100,000 at 0.38 GW/m 2 . The nature of thermomechanical damage in W-foam exposed to extreme power (GW/m 2 ) short-duration laser pulses was found to be qualitatively similar to that of high power (MW/m 2 ) long-duration arc-jet.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Temperature-Dependent Mechanical Anisotropy in Textured Zircaloy Cladding

This work evaluates deformation activity in textured zircaloy-4 cladding with the aim of elucidating the mechanisms that govern mechanical behavior across a temperature range anticipated in reactor service. In particular, the dependence of mechanical anisotropy on temperature and stress state (ratio of applied stresses) is examined via uniaxial and biaxial tensile experiments at room temperature and 400 °C. The mechanical behavior was interpreted based on micro-texture and dominant slip system activity. Thermal activation of basal systems is found as a key mechanism affecting mechanical behavior across the temperature range. The biaxial stress states produce resolved shear stresses contrasting to uniaxial stress states. As a result, mechanical testing and implications for reactor pellet-cladding mechanical interaction stress conditions are discussed.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Optimizing the hypre solver for manycore and GPU architectures

The solution of large-scale combustion problems with codes such as Uintah on modern computer architectures requires the use of multithreading and GPUs to achieve performance. Uintah uses a low-Mach number approximation that requires iteratively solving a large system of linear equations. The Hypre iterative solver has solved such systems in a scalable way for Uintah, but the use of OpenMP with Hypre leads to at least slowdown due to OpenMP overheads. The proposed solution uses the MPI Endpoints within Hypre, where each team of threads acts as a different MPI rank. This approach minimizes OpenMP synchronization overhead and performs as fast or (up to 1.44) faster than Hypre's MPI-only version, and allows the rest of Uintah to be optimized using OpenMP. The profiling of the GPU version of Hypre shows the bottleneck to be the launch overhead of thousands of micro-kernels. The GPU performance was improved by fusing these micro-kernels and was further optimized by using Cuda-aware MPI, resulting in an overall speedup of 1.16—1.44 compared to the baseline GPU implementation. The above optimization strategies were published in the International Conference on Computational Science 2020 [1]. This work extends the previously published research by carrying out the second phase of communication-centered optimizations in Hypre to improve its scalability on large-scale supercomputers. Additionally, this includes an efficient non-blocking inter-thread communication scheme, communication-reducing patch assignment, and expression of logical communication parallelism to a new version of the MPICH library that utilizes the underlying network parallelism [2]. The above optimizations avoid communication bottlenecks previously observed during strong scaling and improve performance by up to 2 on 256 nodes of Intel Knight's Landing processor.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Development of a multi-layer canopy model for E3SM Land Model with support for heterogeneous computing

The vertical structure of vegetation canopies creates micro-climates. However, the land components of most Earth System Models, including the Energy Exascale Earth System Model (E3SM), typically neglect vertical canopy structure by using a single layer big-leaf representation to simulate water, CO 2 , and energy exchanges between the land and the atmosphere. In this study, we developed a Multi-Layer Canopy Model for the E3SM Land Model to resolve the micro-climate created by vegetation canopies. The model developed in this study re-implements the CLM-ml_v1 to support heterogeneous computing architectures consisting of CPUs and GPUs and includes three additional optimization-based stomatal conductance models. The use of Portable, Extensible Toolkit for Scientific Computation provides a speedup of 25–50 times on a GPU relative to a CPU. The numerical implementation of the model was verified against CLM-ml_v1 for a month-long simulation using data from the Ameriflux US-University of Michigan Biological Station site. Model structural uncertainty was explored by performing control simulations for five stomatal conductance models that exclude and include the control of plant hydrodynamics (PHD) on photosynthesis. The bias in simulated sensible and latent heat fluxes was lower when PHD was accounted for in the model. Additionally, six idealized simulations were performed to study the impact of three environmental variables (i.e. air temperature, atmospheric CO 2 , and soil moisture) on canopy processes (i.e. net CO 2 assimilation, leaf temperature, and leaf water potential). Increasing air temperature reduced net CO 2 assimilation and increased air temperature. Net CO 2 assimilation increased at higher atmospheric CO 2 , while decreasing soil moisture resulted in lower leaf water potential.

54 ENVIRONMENTAL SCIENCES↗