Search NASA⌕ Search

SEARCH · Search NASA

Results for “Turbulent Heating”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Preface for frontiers of magnetic reconnection research in heliophysical, astrophysical, and laboratory plasmas

Magnetic reconnection—the topological rearrangement of magnetic field—underlies many explosive phenomena across a wide range of natural and laboratory plasmas.3 It plays a pivotal role in electron and ion heating, particle acceleration to high energies, energy transport, and self-organization. Reconnection can have a complex relationship with turbulence at both large and small scales, leading to various effects that are only beginning to be understood. In heliophysics, magnetic reconnection plays a key role in solar flares, coronal mass ejections, coronal heating, solar wind dissipation, the interaction of interplanetary plasma with magnetospheres, dynamics of planetary magnetospheres such as magnetic substorms, and the heliospheric boundary with the interstellar medium. Here, the magnetic reconnection is integral to the solar and planetary dynamo processes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Improving a Multilevel Turbulence Closure Model for a Shallow Lake in Comparison With Other 1-D Models

Lakes differ from lands in water availability, heat capacity, albedo, and roughness, which affect local surface-atmospheric interactions. This study modified a multilevel upper ocean model (UOM) for lake applications and evaluated its performance in Lake Taihu (China) with comprehensive measurements against three popular one-dimensional (1-D) lake models. These models were based on different concepts, including the self-similarity (FLake), the wind-driven eddy diffusion (LISSS), the k-ε turbulence closure (SIMSTRAT), and a simplified turbulence closure (UOM). The surface flux scheme in these models was unified to exclude the discrepancies in representing air-lake exchanges. All models in their default formulations presented obvious cold water temperature biases and largely underestimated the lake surface temperature (LST) diurnal range. For each model, these deficiencies were significantly reduced by incorporating new physics schemes or calibrated tunable parameters based on systematic sensitivity tests. The primary modifications for UOM included (1) a new scheme of decreased surface roughness lengths to better characterize the shallow lake, (2) a solar radiation penetration scheme with increased light extinction coefficient and surface absorption fraction to account for the high water turbidity, and (3) turbulent Prandtl number increased by a factor of 20 to reduce the turbulent vertical mixing. All other models were improved in these three aspects (roughness, extinction, and mixing) within their original formulations. Given these improvements, UOM showed superior performance to other models in capturing LST diurnal cycle and daily to seasonal variations, as well as summer-autumn vertical stratification changes. The new UOM is well suited for application in shallow lakes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Ducted fuel injection

Various technologies presented herein relate to heating a fuel and enhancing fuel and charge-gas mixing inside a combustion chamber to enable minimal, or no, generation of soot and/or other undesired emissions during ignition and subsequent combustion of the locally premixed mixtures. To enable sufficient mixing of the fuel and charge-gas, a jet of fuel can be directed to pass through a bore of a duct causing charge-gas to be drawn into the bore creating turbulence to mix the fuel and the drawn charge-gas. The duct can be heated to provision heating of the fuel and charge-gas mixture. The duct can be located proximate to an opening in a tip of a fuel injector. An ignition assist component can be located downstream of the duct to facilitate ignition of the fuel/charge-gas mixture.

42 ENGINEERING↗
DOE OSTI2021

Scattering of radio frequency waves by randomly modulated density interfaces in the edge of fusion plasmas

In the scrape-off layer and the edge region of a tokamak, the plasma is strongly turbulent and scatters the radio frequency (RF) electromagnetic waves that propagate through this region. It is important to know, whether used for diagnostics or for heating and current drive, the spectral properties of these scattered RF waves. The spectral changes influences the interpretation of the diagnostic-data obtained and the current and heating profiles. A full-wave, 3D electromagnetic code ScaRF (see Papadopoulos et al. 2019) has been developed for studying the RF wave propagation through turbulent plasma. ScaRF is a finite-difference frequency-domain (FDFD) method for solving Maxwell's equations. The magnetized plasma is defined through the cold plasma, anisotropic permittivity tensor. As a result, ScaRF can be used to study the scattering of any cold plasma RF wave. It can be for the study of scattering of electron cyclotron waves in ITER-type and medium-sized tokamaks such as TCV, ASDEX-U, DIII-D. For the case of medium-sized tokamaks, there's experimental evidence that drift waves and rippling modes are present in the edge region (see Ritz et al. 1984). Hence, we study the scattering of RF waves by periodic density interfaces (plasma gratings) in the form of a superposition of spatial modes with varying periodicity and random amplitudes (see Papadopoulos et al. 2019). The power reflection coefficient (a random variable) is calculated for different realizations of the density interface. In this work, the uncertainty of the power reflection coefficient is rigorously quantified by use of the Polynomial Chaos Expansion (see Xiu & Karniadakis 2002) method in conjunction with the Smolyak sparse grid integration (see Papadopoulos et al. 2018) (PCE-SG). The PCE-SG method is proven accurate and much more efficient (roughly 2-orders of magnitude shorter execution time) compared to alternative methods such as the Monte Carlo (MC) approach.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

In situ recording of Mars soundscape

Abstract Before the Perseverance rover landing, the acoustic environment of Mars was unknown. Models predicted that: (1) atmospheric turbulence changes at centimetre scales or smaller at the point where molecular viscosity converts kinetic energy into heat 1 , (2) the speed of sound varies at the surface with frequency 2,3 and (3) high-frequency waves are strongly attenuated with distance in CO 2 (refs. 2–4 ). However, theoretical models were uncertain because of a lack of experimental data at low pressure and the difficulty to characterize turbulence or attenuation in a closed environment. Here, using Perseverance microphone recordings, we present the first characterization of the acoustic environment on Mars and pressure fluctuations in the audible range and beyond, from 20 Hz to 50 kHz. We find that atmospheric sounds extend measurements of pressure variations down to 1,000 times smaller scales than ever observed before, showing a dissipative regime extending over five orders of magnitude in energy. Using point sources of sound (Ingenuity rotorcraft, laser-induced sparks), we highlight two distinct values for the speed of sound that are about 10 m s −1 apart below and above 240 Hz, a unique characteristic of low-pressure CO 2 -dominated atmosphere. We also provide the acoustic attenuation with distance above 2 kHz, allowing us to explain the large contribution of the CO 2 vibrational relaxation in the audible range. These results establish a ground truth for the modelling of acoustic processes, which is critical for studies in atmospheres such as those of Mars and Venus.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Multi-scale Interaction Mechanism for Edge-Localized-Mode Suppression in the Tokamak Edge

A central challenge in fusion energy is reconciling the high-confinement mode required for reactor performance with the intense intermittent relaxation events it produces, known as edge-localized modes. These instabilities arise in the steep pressure pedestal at the plasma edge when magnetohydrodynamic thresholds are crossed, inflicting damaging heat loads on reactor components. Here, we show that multiscale interactions between microscopic turbulence and macroscopic magnetohydrodynamic modes provide encouraging prospects for self-organized edge-localized modes regulation. Using direct quantitative measurements of multiscale modes, eddy dynamics, and turbulent flux, we show that small-scale electron drift wave turbulence actively scatters the large-scale peeling-ballooning modes. This scattering decorrelates the pressure and velocity fields of the instability, so arresting its growth. Our modeling and theoretical analysis confirm this suppression mechanism is effective even when conventional linear stability thresholds are exceeded. This work establishes a nonlinear principle for edge-localized modes stability, revealing how ambient micro-turbulence can be leveraged to maintain a macro-stable, high-performance pedestal for future fusion reactors.

experimental nuclear physics↗
DOE OSTI2020

NUMERICAL INVESTIGATION OF FUEL PROPERTY EFFECTS ON MIXED-MODE COMBUSTION IN A SPARK-IGNITION ENGINE

In the present study, mixed-mode combustion of an E30 fuel in a direct-injection spark-ignition engine is numerically investigated at a fuel-lean operating condition using multidimensional computational fluid dynamics (CFD). A fuel surrogate matching Research Octane Number (RON) and Motor Octane Number (MON) of E30 is first developed using neural network based non-linear regression model. To enable efficient 3D engine simulations, a 164-species skeletal reaction mechanism incorporating NO, chemistry is reducedfrom a detailed chemical kinetic model. A hybrid approach that incorporates the G-equation model for tracking turbulentflame front, and the multi-zone well-stirred reactor model for predicting auto-ignition in the end gas, is employed to account for turbulent combustion interactions in the engine cylinder. Predicted in-cylinder pressure and heat release rate traces agree well with experimental measurements. The proposed modelling approach also captures moderated cyclic variability. Two different types of combustion cycles, corresponding to purely deflagrative and mixed-mode combustion, are observed. In contrast to the purely deflagrative cycles, mixed-mode combustion cycles feature early flame propagation followed by end-gas auto-ignition, leading to two distinctive peaks in heat release rate traces. The positive correlation between mixed-mode combustion cycles and early flame propagation is well captured by simulations. With the validated numerical setup, effects of NOx chemistry on mixed-mode combustion predictions are investigated. NOx chemistry is found to promote auto-ignition through residual gas recirculation, while the deflagrative flame propagation phase remains largely unaffected. Local sensitivity analysis is then performed to understand effects of physical and chemical properties of the fuel, i.e., heat of evaporation (HoV) and laminar flame speed (S-L). An increased HoV tends to suppress end-gas auto-ignition due to increased vaporization cooling, while the impact of HoV on flame propagation is insignificant. In contrast, an increased SL is found to significantly promote both flame propagation and auto-ignition. The promoting effect of SL on auto-ignition is not a direct chemical effect; it is rather caused by an advancement of the combustion phasing, which increases compression heating of the end gas.

detailed chemical kinetics↗
DOE OSTI2022

Simulation and Reynolds-averaged Navier-Stokes modeling of a three-component Rayleigh-Taylor mixing problem with thermonuclear burn

Rayleigh-Taylor mixing in the presence of a third component with intermediate density is investigated through three-dimensional large-eddy simulation (LES) with a high-order compact finite-difference code. Two configurations are considered: (1) a symmetric configuration in which the Atwood number between the heavy and intermediate components matches the Atwood number between the intermediate and light components and (2) an asymmetric configuration in which the Atwood number between the heavy and intermediate components is an order of magnitude greater than the Atwood number between the intermediate and light components. Mass fraction covariances are extracted, and proposed Reynolds-averaged Navier-Stokes (RANS) closures for density-specific-volume and density-mass-fraction covariances are evaluated in an a priori fashion. Additionally, a multicomponent extension of the k - Φ - L - a - V RANS model [Morgan, Phys. Rev. E 104, 015107 (2021)] is presented which includes model equations for the upper-triangular elements of the mass fraction covariance matrix. This model, referred to as the k - Φ - L - a - C model, is compared against results from LES and against other RANS models. Profiles of average mass fraction, mass-fraction covariance, and density-specific-volume covariance obtained with the k - Φ - L - a - C model are found to agree well with LES data. Finally, the impact of three-component turbulent mixing on average reaction rate is investigated in both premixed and nonpremixed cases by heating the mixing layer and allowing it to undergo thermonuclear (TN) burn. A closure model for average reaction rate is proposed for use with the k - Φ - L - a - C model, and when this model is applied, improved agreement is obtained between LES and RANS in total TN neutron production.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Overview of Wendelstein 7-X high-performance operation

The Wendelstein 7-X (W7-X) stellarator has completed two consecutive experimental campaigns OP 2.2 (Sep.-Dec. 2024) and OP 2.3 (Feb.-May 2025) under a new operational strategy enabling more than one year of uninterrupted device availability. This approach, supported by exceptionally high subsystem reliability, allowed sustained high-efficiency plasma operations with up to 80–100 discharges per day across a broad range of magnetic configurations. Several key technical upgrades-most notably the first operation of a 1.5 MW class steady-state gyrotron, a new steady-state pellet injector, and advanced real-time feedback control systems significantly enhanced heating, fueling, and plasma control capabilities. Together, these improvements enabled major advances in long-pulse performance, high-β operation, and confinement optimization. Long-pulse discharges achieved 1.8 GJ of injected energy under fully detached divertor conditions, while reduced-field scenarios facilitated record volume-averaged β values approaching 3%. High-performance plasmas with centrally peaked density profiles, created via neutral beam injection (NBI) or sustained pellet fueling, demonstrated strongly reduced turbulent transport and stellarator-record fusion triple products. Complementary studies of power exhaust and divertor heat loads revealed the role of scrape-off-layer drift physics in shaping strike-line patterns under attached conditions. Together, the results from OP 2.2 and OP 2.3 significantly expand the operational space of W7-X and strengthen its role as a leading platform for steady-state stellarator research and reactor-relevant plasma scenarios.

Wendelstein 7-X↗
DOE OSTI2025

Dual-band Ion-scale Waves in the Solar Wind: Statistical Properties and Comparison with Single-band Ion-scale Waves

We present a comprehensive statistical analysis of ion-scale waves including dual-band ion-scale waves (DBIWs) observed by the Parker Solar Probe in the solar wind. DBIWs are characterized by the simultaneous occurrence of distinct left- and right-handed polarized wave packets at higher and lower frequencies, respectively, in the spacecraft frame, implying the presence of bidirectionally propagating L-mode waves. We identify and compare DBIWs to single-band ion-scale waves (SBIWs) and investigate their spatial distributions, background magnetic field conditions, and associated inertial range turbulence characteristics. We find that DBIWs exhibit a strong preferential occurrence at heliocentric distances between 0.2 and 0.25 au (equivalently 43 and 54 R S ) in contrast to SBIWs, the occurrence rate of which peaks closer to the Sun and decreases with the heliocentric distance. DBIWs are also preferentially observed under conditions of strong Alfvénic turbulence, characterized by low magnetic compressibility and perpendicular-dominant fluctuations. These findings support a scenario in which stochastic ion heating via sufficiently large perpendicular Alfvénic fluctuation generates the ion temperature anisotropies required to excite bidirectional ion cyclotron waves. This work advances our understanding of kinetic wave generation in the solar wind and its connection to large-scale turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Increasing Main Cooler Thermal Performance for sCO2 Power Cycles

Increasing sCO2 power cycle cooler performance will increase the efficiency and lower the cost of electricity for sCO2 power cycles. This work considered the heat transfer performance of two counter-flow, shell-and-tube heat exchangers using local measurements. The heat exchangers were constructed using commercial tubing. One heat exchanger was constructed with a conventional, smooth inner tube (7 mm nominal inner diameter). The other was constructed with an additively manufactured tube with a square cross section and angled rib turbulators. For the present study, the measured Nusselt number was 30% less than the Dittus-Boelter correlation. The decrease was attributed to heat transfer degradation. The cooler performance was also reported as heat exchanger effectiveness. A factor of two improvement in Nusselt number (94%) relative to smooth conventional tubes was observed, indicating that heat transfer enhancement features (such as the ribs) are a suitable technology to overcome heat transfer degradation in sCO2 coolers. Trends in heat exchanger effectiveness indicated that the ribs were most beneficial when the water flow rates were at the low end of the range considered: 0.016 – 0.126 kg/s. Using results from prior system modeling, it is anticipated that adding heat transfer enhancement features in a cycle cooler will increase the sCO2 cycle efficiency by 0.34 percentage points.

Searle, Matthew↗
DOE OSTI2021

Summary of 2020 edge-plasma modeling for FNSFand plans for 2021

This year, following the study of liquid-wall options, our edge-plasma modeling has returned to analyzing issues associated with a solid tungsten-wall, Fusion Nuclear Science Facility (FNSF)-size device, primarily using the 2D UEDGE plasma/neutral transport code. The focus is to predict plasma/neutral heat and particle fluxes to various wall components and determine the concentration of injected and wall-sputtered impurities that reach the core plasma region. Our models have been advanced in a number of areas compared to the initial modeling done in the 2017 timeframe: this year we include each individual charge-state of impurity ions instead of assuming a fixed concentration of impurities relative to the deuterium-tritium (DT) fuel plasma, including the interaction between inner and outer divertor legs. Here edge fueling of the impurity and pumping in the private-flux region is incorporated. Increasingly detailed atomic-physics rate files for impurity radiation have been compared, as well as the role of molecules. The impact of spatial variation of the turbulence-driven radial plasma transport is considered to contribute to the understanding of the heat-flux-limit effects of choosing a double-null or single-null divertor configuration. Detailed plasma/neutral wall fluxes, densities, and temperatures have been provided to the ERO group for their more detailed studies of wall sputtering and near-surface impurity transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Production and Persistence of Extreme Two-temperature Plasmas in Radiative Relativistic Turbulence

Turbulence is a predominant process for energizing electrons and ions in collisionless astrophysical plasmas, and thus is responsible for shaping their radiative signatures (luminosity, spectra, and variability). To better understand the kinetic properties of a collisionless radiative plasma subject to externally driven turbulence, we investigate particle-in-cell simulations of relativistic plasma turbulence with external inverse Compton cooling acting on the electrons. We find that ions continuously heat up while electrons gradually cool down (due to the net effect of radiation), and hence the ion-to-electron temperature ratio T i /T e grows in time. Furthermore, we show that Ti/Te is limited only by the size and duration of the simulations (reaching ${T}_{i}/{T}_{e}\sim {10}^{3}$), indicating that there are no efficient collisionless mechanisms of electron–ion thermal coupling. This result has implications for models of radiatively inefficient accretion flows, such as observed in the Galactic center and in M87, for which so-called two-temperature plasmas with ${T}_{i}/{T}_{e}\gg 1$ have been invoked to explain their low luminosity. Additionally, we find that electrons acquire a quasi-thermal distribution (dictated by the competition of turbulent particle energization and radiative cooling), while ions undergo efficient nonthermal acceleration (acquiring a harder distribution than in equivalent nonradiative simulations). There is a modest nonthermal population of high-energy electrons that are beamed intermittently in space, time, and direction; these beamed electrons may explain rapid flares in certain high-energy astrophysical systems (e.g., in the Galactic center). These numerical results demonstrate that extreme two-temperature plasmas can be produced and maintained by relativistic radiative turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Scattering of radiofrequency waves by randomly modulated density interfaces in the edge of fusion plasmas

In the scrape-off layer and the edge region of a tokamak, the plasma is strongly turbulent and scatters the radiofrequency (RF) electromagnetic waves that propagate through this region. It is important to know the spectral properties of these scattered RF waves, whether used for diagnostics or for heating and current drive. The spectral changes influence the interpretation of the obtained diagnostic data, and the current and heating profiles. A full-wave, three-dimensional (3-D) electromagnetic code ScaRF has been developed for studying the RF wave propagation through turbulent plasma. ScaRF is a finite-difference frequency-domain (FDFD) method used for solving Maxwell's equations. The magnetized plasma is defined through the cold plasma by the anisotropic permittivity tensor. As a result, ScaRF can be used to study the scattering of any cold plasma RF wave. It can also be used for the study of the scattering of electron cyclotron waves in ITER-type and medium-sized tokamaks such as TCV, ASDEX-U and DIII-D. For the case of medium-sized tokamaks, there is experimental evidence that drift waves and rippling modes are present in the edge region. Hence, we have studied the scattering of RF waves by periodic density interfaces (plasma gratings) in the form of a superposition of spatial modes with varying periodicity and random amplitudes. The power reflection coefficient (a random variable) is calculated for different realizations of the density interface. In this work, the uncertainty of the power reflection coefficient is rigorously quantified by use of the Polynomial Chaos Expansion method in conjunction with the Smolyak sparse-grid integration, which is known as the PCE-SG method. The PCE-SG method is proven to be accurate and more efficient compared with alternative methods such as the Monte Carlo (MC) approach.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Validation Assessment of Turbulent Reacting Flow Model Using The Area-Validation Metric on Medium-Scale Methanol Pool Fire Results

Accident analysis and ensuring power plant safety are pivotal in the nuclear energy sector. Significant strides have been achieved over the past few decades regarding fire protection and safety, primarily centered on design and regulatory compliance. Yet, after the Fukushima accident a decade ago, the imperative to enhance measures against fire, internal flooding, and power loss has intensified. Hence, a comprehensive, multilayered protection strategy against severe accidents is needed. Consequently, gaining a deeper insight into pool fires and their behavior through extensive validated data can greatly aid in improving these measures using advanced validation techniques. A model validation study was performed at Sandia National Laboratories in which a 30-cm diameter methanol pool fire was modeled using the SIERRA/Fuego turbulent reacting flow code. This validation study used a standard validation experiment to compare model results against, and conclusions have been published. The fire was modeled with a Large Eddy Simulation (LES) turbulence model with subgrid turbulent kinetic energy closure. Combustion was modeled using a strained laminar flamelet library approach. Radiative heat transfer was accounted for with a model utilizing the gray-gas approximation. In the present study, additional validation analysis is performed using the area validation metric (AVM). These activities are done on multiple datasets involving different variables and temporal/spatial ranges and intervals. In conclusion, the results provide insight into the use of the area validation metric on such temporally varying datasets and the importance of physics-aware use of the metric for proper analysis.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2021

Assessment of Turbulent Combustion Models for Simulating Prechamber Ignition in a Natural Gas Engine

Lean combustion is a promising strategy to increase thermal efficiency in an internal combustion engine, by exploiting a favorable specific heat ratio of the fresh mixture while simultaneously suppressing the heat losses to the cylinder wall. However, unstable ignition and slow flame propagation at fuel-lean conditions lead to large cycle-to-cycle variability and limit the high-efficiency engine operating range. Prechamber ignition is considered an effective concept to extend the lean operating limit, by providing spatially distributed ignition with multiple turbulent flame-jets and enabling a faster combustion rate compared to the conventional spark ignition approach. From a numerical modeling standpoint to date science base and available simulation tools are inadequate to properly understand and predict the combustion processes in prechamber ignited engines. In this paper, conceptually different Reynolds-averaged Navier–Stokes (RANS) combustion models widely adopted in the engine modeling community are used to simulate the ignition and combustion processes in a medium-duty natural gas engine with a prechamber spark-ignition system. Here, a flamelet-based turbulent combustion model, i.e., G-equation, and a multizone well-stirred reactor model are employed for this modeling study. Simulation results are compared with experimental data in terms of in-cylinder pressure and heat release rate. Finally, the analysis of the performance of the two models is carried out to highlight the strengths and limitations of the two evaluated approaches.

42 ENGINEERING↗
DOE OSTI2021

Application of Scalar Filtered Density Function to Turbulent Flows Under Supercritical Condition

Abstract The scalar filtered density function (FDF) methodology is extended and employed for large eddy simulation (LES) of turbulent flows under supercritical condition. To describe real fluid behavior, the extended methodology incorporates the generalized heat and mass diffusion models along with real fluid thermodynamic relations which are derived using the cubic Peng–Robinson equation of state. These models are implemented within the stochastic differential equations comprising the scalar FDF transport. Simulations are conducted of a temporally developing mixing layer under supercritical condition and the results are assessed by comparing with data generated by direct numerical simulation (DNS) of the same layer. The consistency of the proposed FDF methodology is assessed. The LES-FDF predictions are shown to agree favorably with the DNS data and exhibit several key features pertaining to supercritical turbulent flows.

Energy & Fuels↗
DOE OSTI2022

A Hybrid Bulk Algorithm to Predict Turbulent Fluxes over Dry and Wet Bare Soils

Measurements made in the Columbia River basin (Oregon) in an area of irregular terrain during the second Wind Forecast Improvement Project (WFIP2) field campaign are used to develop an optimized hybrid bulk algorithm to predict the surface turbulent fluxes from readily measured or modeled quantities over dry and wet bare or lightly vegetated soil surfaces. The hybrid (synthetic) algorithm combines (i) an aerodynamic method for turbulent flow, which is based on the transfer coefficients (drag coefficient and Stanton number), roughness lengths, and Monin–Obukhov similarity; and (ii) a modified Priestley–Taylor (P-T) algorithm with physically based ecophysiological constraints, which is essentially based on the surface energy budget (SEB) equation. Soil heat flux in the latter case was estimated from measurements of soil temperature and soil moisture. In the framework of the hybrid algorithm, bulk estimates of the momentum flux and the sensible heat flux are derived from a traditional aerodynamic approach, whereas the latent heat flux (or moisture flux) is evaluated from a modified P-T model. Direct measurements of the surface fluxes (turbulent and radiative) and other ancillary atmospheric/soil parameters made during WFIP2 for different soil conditions (dry and wet) are used to optimize and tune the hybrid bulk algorithm. The bulk flux estimates are validated against the measured eddy-covariance fluxes. We also discuss the SEB closure over dry and wet surfaces at various time scales based on the modeled and measured fluxes. Although this bulk flux algorithm is optimized for the data collected during the WFIP2, a hybrid approach can be used for similar flux-tower sites and field campaigns.

54 ENVIRONMENTAL SCIENCES↗