SEARCH · Engineering Papers

Results for “Mixing”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

NASA NTRS1993

A bilinear failure criterion for mixed-mode delamination

Many different failure criteria have been suggested for mixed-mode delamination toughness, but few sets of mixed-mode data exist that are consistent over the full range of Mode 1 opening load to Mode 2 shear load range. The mixed-mode bending (MMB) test was used to measure the delamination toughness of a brittle epoxy composite, a state-of-the-art toughened epoxy composite, and a tough thermoplastic composite over the full mixed-mode range. To gain insight into the different failure responses of the different materials, the delamination fracture surfaces were also examined. An evaluation of several failure criteria that have been reported in the literature was performed, and the range of responses modeled by each criterion was analyzed. A bilinear failure criterion was introduced based on a change in the failure mechanism observed from the delamination surfaces. The different criteria were compared to the failure response of the three materials tested. The responses of the two epoxies were best modeled with the new bilinear failure criterion. The failure response of the tough thermoplastic composite could be modeled well with the bilinear criterion but could also be modeled with the more simple linear failure criterion. Since the materials differed in their mixed-mode failure response, mixed-mode delamination testing will be needed to characterize a composite material. This paper presents consistent sets of mixed-mode data, provides a critical evaluation of the mixed-mode failure criteria, and should provide general guidance for selecting an appropriate criterion for other materials.

Reeder, James R.↗

NASA NTRS1997

Observational and Model Studies of Large-Scale Mixing Processes in the Stratosphere

The following is the final technical report for grant NAGW-3442, 'Observational and Model Studies of Large-Scale Mixing Processes in the Stratosphere'. Research efforts in the first year concentrated on transport and mixing processes in the polar vortices. Three papers on mixing in the Antarctic were published. The first was a numerical modeling study of wavebreaking and mixing and their relationship to the period of observed stratospheric waves (Bowman). The second paper presented evidence from TOMS for wavebreaking in the Antarctic (Bowman and Mangus 1993). The third paper used Lagrangian trajectory calculations from analyzed winds to show that there is very little transport into the Antarctic polar vortex prior to the vortex breakdown (Bowman). Mixing is significantly greater at lower levels. This research helped to confirm theoretical arguments for vortex isolation and data from the Antarctic field experiments that were interpreted as indicating isolation. A Ph.D. student, Steve Dahlberg, used the trajectory approach to investigate mixing and transport in the Arctic. While the Arctic vortex is much more disturbed than the Antarctic, there still appears to be relatively little transport across the vortex boundary at 450 K prior to the vortex breakdown. The primary reason for the absence of an ozone hole in the Arctic is the earlier warming and breakdown of the vortex compared to the Antarctic, not replenishment of ozone by greater transport. Two papers describing these results have appeared (Dahlberg and Bowman; Dahlberg and Bowman). Steve Dahlberg completed his Ph.D. thesis (Dahlberg and Bowman) and is now teaching in the Physics Department at Concordia College. We also prepared an analysis of the QBO in SBUV ozone data (Hollandsworth et al.). A numerical study in collaboration with Dr. Ping Chen investigated mixing by barotropic instability, which is the probable origin of the 4-day wave in the upper stratosphere (Bowman and Chen). The important result from this paper is that even in the presence of growing, unstable waves, the mixing barriers around

Bowman, Kenneth P.↗

NASA NTRS2001

Storm-driven Mixing and Potential Impact on the Arctic Ocean

Observations of the ocean, atmosphere, and ice made by Ice-Ocean Environmental Buoys (IOEBs) indicate that mixing events reaching the depth of the halocline have occurred in various regions in the Arctic Ocean. Our analysis suggests that these mixing events were mechanically forced by intense storms moving across the buoy sites. In this study, we analyzed these mixing events in the context of storm developments that occurred in the Beaufort Sea and in the general area just north of Fram Strait, two areas with quite different hydrographic structures. The Beaufort Sea is strongly influenced by inflow of Pacific water through Bering Strait, while the area north of Fram Strait is directly affected by the inflow of warm and salty North Atlantic water. Our analyses of the basin-wide evolution of the surface pressure and geostrophic wind fields indicate that the characteristics of the storms could be very different. The buoy-observed mixing occurred only in the spring and winter seasons when the stratification was relatively weak. This indicates the importance of stratification, although the mixing itself was mechanically driven. We also analyze the distribution of storms, both the long-term climatology as well as the patterns for each year in the last two decades. The frequency of storms is also shown to be correlated- (but not strongly) to Arctic Oscillation indices. This study indicates that the formation of new ice that leads to brine rejection is unlikely the mechanism that results in the type of mixing that could overturn the halocline. On the other hand, synoptic-scale storms can force mixing deep enough to the halocline and thermocline layer. Despite a very stable stratification associated with the Arctic halocline, the warm subsurface thermocline water is not always insulated from the mixed layer.

Yang, Jiayan↗

NASA NTRS2003

Mixing and NO(x) Emission Calculations of Confined Reacting Jet Flows in a Cylindrical Duct

Rapid mixing of cold lateral jets with hot cross-stream flows in confined configurations is of practical interest in gas turbine combustors as it strongly affects combustor exit temperature quality, and gaseous emissions in for example rich-lean combustion. It is therefore important to further improve our fundamental understanding of the important processes of dilution jet mixing especially when the injected jet mass flow rate exceeds that of the cross-stream. The results reported in this report describe some of the main flow characteristics which develop in the mixing process in a cylindrical duct. A 3-dimensional tool has been used to predict the mixing flow field characteristics and NOx emission in a quench section of an RQL combustor, Eighteen configurations have been analyzed in a circular geometry in a fully reacting environment simulating the operating condition of an actual RQL gas turbine combustion liner. The evaluation matrix was constructed by varying three parameters: 1) jet-to-mainstream momentum-flux ratio (J), 2) orifice shape or orifice aspect ratio, and 3) slot slant angle. The results indicate that the mixing flow field significantly varies with the value of the jet penetration and subsequently, slanting elongated slots generally improve the mixing uniformity at high J conditions. Round orifices produce more uniform mixing and low NO(x) emissions at low J due to the strong and adequate jet penetration. No significant correlation was found between the NO(x) production rates and the mixing deviation parameters, however, strong correlation was found between NO(x) formation and jet penetration. In the computational results, most of the NO(x) formation occurred behind the orifice starting at the orifice wake region. Additional NO(x) is formed upstream of the orifice in certain configurations with high J conditions due to the upstream recirculation.

Holdeman, James D.↗

NASA NTRS2013

Evaluation of Aerosol Mixing State Classes in the GISS Modele-matrix Climate Model Using Single-particle Mass Spectrometry Measurements

Aerosol particles in the atmosphere are composed of multiple chemical species. The aerosol mixing state, which describes how chemical species are mixed at the single-particle level, provides critical information on microphysical characteristics that determine the interaction of aerosols with the climate system. The evaluation of mixing state has become the next challenge. This study uses aerosol time-of-flight mass spectrometry (ATOFMS) data and compares the results to those of the Goddard Institute for Space Studies modelE-MATRIX (Multiconfiguration Aerosol TRacker of mIXing state) model, a global climate model that includes a detailed aerosol microphysical scheme. We use data from field campaigns that examine a variety of air mass regimens (urban, rural, and maritime). At all locations, polluted areas in California (Riverside, La Jolla, and Long Beach), a remote location in the Sierra Nevada Mountains (Sugar Pine) and observations from Jeju (South Korea), the majority of aerosol species are internally mixed. Coarse aerosol particles, those above 1 micron, are typically aged, such as coated dust or reacted sea-salt particles. Particles below 1 micron contain large fractions of organic material, internally-mixed with sulfate and black carbon, and few external mixtures. We conclude that observations taken over multiple weeks characterize typical air mass types at a given location well; however, due to the instrumentation, we could not evaluate mass budgets. These results represent the first detailed comparison of single-particle mixing states in a global climate model with real-time single-particle mass spectrometry data, an important step in improving the representation of mixing state in global climate models.

climate models↗

NASA NTRS2018

Comparison of Mixing Characteristics for Several Fuel Injectors at Mach 8, 12, and 15 Hypervelocity Flow Conditions

CFD analysis is presented of the mixing characteristics and performance of three fuel injectors at flight Mach numbers of 8, 12, and 15. The Reynolds-averaged simulations (RAS) were carried out using the VULCAN-CFD solver. The high Mach number flow conditions match those of the experiments conducted as a part of the Enhanced Injection and Mixing Project (EIMP) at the NASA Langley Research Center. The EIMP aims to investigate scramjet fuel injection and mixing physics, improve the understanding of underlying physical processes, and develop enhancement strategies relevant to flight Mach numbers greater than 8. The injectors include a fuel placement device, a strut, and a fluidic vortical mixer, a ramp. These fuel injectors accomplish the necessary task of distributing and mixing fuel into the supersonic cross-flow, albeit via different strategies. For comparison, a flush-wall injector is also included. This type of injector generally represents the simplest method of introducing fuel into a scramjet combustor. The three injectors represent the baseline configurations of the EIMP experiments. The mixing parameters of interest, such as mixing efficiency and total pressure recovery, are computed from the RAS and compared for the three flight conditions and injector configurations. In addition to mixing efficiency and total pressure recovery, the combustion efficiency and thrust potential are also computed for the reacting simulations. Plotting the total pressure recovery and thrust potential as a function of mixing efficiency provides added insight into critical aspects of combustor performance as the flight condition and injector type are varied.

Shenoy, Rajiv R.↗

DOE OSTI2023

Fluid mixing apparatus using high- and low-pressure fluid streams

A fluid mixing apparatus includes mixing conduits that extend through a fluid plenum and that define injection holes therethrough. The fluid plenum, which surrounds a first wall defining a main passage fluidly coupled to a low-pressure fluid source, is surrounded itself by a second wall defining a high-pressure plenum fluidly coupled to a high-pressure fluid source. The mixing conduits fluidly couple the high-pressure plenum to the main passage, and the fluid from the fluid plenum is delivered with the high-pressure fluid to the main passage, where the fluids mix before being discharged from an outlet of the main passage. The fluid mixing apparatus may be used to mix one or more fuels with high- and low-pressure air in a gas turbine combustor. Alternately, the fluid mixing apparatus may mix a fluid with high- and low-pressure water streams.

Berry, Jonathan Dwight↗

DOE OSTI2024

Numerical Mixing Suppresses Submesoscale Baroclinic Instabilities Over Sloping Bathymetry

Abstract The impacts of spurious numerical salinity mixing on the larger‐scale flow and tracer fields are characterized using idealized simulations. The idealized model is motivated by realistic simulations of the Texas‐Louisiana shelf and features oscillatory near‐inertial wind forcing. can exceed the physical mixing from the turbulence closure in frontal zones and within the mixed layer. This suggests that simulated mixing processes in frontal zones are driven largely by . Near‐inertial alongshore wind stress amplitude is varied to identify a base case that maximizes the ratio of to in simulations with no prescribed horizontal mixing. We then test the sensitivity of the base case with three tracer advection schemes (MPDATA, U3HC4, and HSIMT) and conduct ensemble runs with perturbed bathymetry. Instability growth is evaluated using the volume‐integrated eddy kinetic energy and available potential energy . While all schemes have similar total mixing, the HSIMT simulations have over double the volume‐integrated and 20% less relative to other schemes, which suppresses the release of and reduces the by roughly 25%. This results in reduced isohaline variability and steeper isopycnals, evidence that enhanced suppresses instability growth. Differences in and between the MPDATA and U3HC4 simulations are marginal. However, the U3HC4 simulations have 25% more . Experiments with variable horizontal viscosity and diffusivity coefficients show that small amounts of prescribed horizontal mixing improve the representation of the ocean state for all advection schemes by reducing the and increasing the .

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Transient mixing behavior of a supercritical fluid injected into supercritical and subcritical environments

Transient, near-field mixing of a supercritical fluid injected into supercritical and subcritical environments is investigated experimentally. Fuel (propane) from a 100 μm diameter Bosch diesel injector is injected into a high-pressure test rig with a steady flow of heated, low-velocity nitrogen to control the ambient pressure and temperature. The fuel supply pressure is approximately 190 bar, allowing short injections lasting 3–4 ms to be performed in quick succession to acquire statistically significant datasets to represent the unsteady fuel–air mixing. Two test cases with the same ambient density and injector inlet pressure and temperature were chosen in an attempt to isolate effects of mixture phase behavior. Microscopic rainbow Schlieren deflectometry (RSD) is used for visualization and quantitative data acquisition. In the case of the supercritical ambient, the injector inlet and exit and ambient conditions are all supercritical. The RSD analysis reveals a mixing process similar to the single phase fluid mixing, though condensation is also observed locally, and can be explained by the adiabatic mixing theory presented in the paper. For the subcritical ambient case, the injector inlet is supercritical, but the injector exit and ambient conditions are both subcritical. The results show a rapid expansion in the near-field promoting mixing between propane and nitrogen. At downstream locations, significant condensation is also observed as predicted by the adiabatic mixing theory. As a result, the supercritical ambient case exhibited longer and narrower jet penetration compared to the subcritical ambient case.

02 PETROLEUM↗

DOE OSTI2025

Mixed-valence halide perovskites

Mixed-valence compounds—which feature an element in at least two different oxidation states—can display emergent optical and transport phenomena stemming from electron transfer between the different valences (intervalence charge-transfer; IVCT). As halide perovskites show promise as active materials in numerous optoelectronic devices, it is an opportune moment to incorporate and study the effects of mixed-valence in this versatile materials family, to access tunable electronic structures ranging from insulators to semiconductors to metals. Herein, we introduce the basic concepts of mixed-valence in molecules and discuss how these concepts may be extended to mixed-valence in extended solids. We then review the few studies of mixed valence in 3D and 2D halide perovskites and halide perovskites with mixed-valence impurities, ranging from studies in the early 1900s to the present day. In conclusion, through judicious choice of metal ion, its coordinating ligands and their geometry, and overall structural dimensionality, chemists can exert powerful synthetic control over electronic delocalization in mixed-valence perovskites, and we hope to see this intriguing materials class expand to encompass new compositions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Numerical study of multi-component flow and mixing in a scaled fission product venting system

Numerical simulation of isothermal multi-component gas flow was conducted to study fluid flow and mixing in a scaled low-pressure, low-temperature test facility for the fission product venting system (FPVS) of a gas cooled fast reactor (GFR). The geometry of FPVS test facility was an open loop including gradual expansion coupling and two 90°pipe elbows. First, the validation of large eddy simulation (LES) wall-adapting local eddy-viscosity (WALE) turbulent model was performed for transitional flow in a circular pipe with gradual expansion. The reactingFOAM solver in OpenFOAM v8 was employed. By introducing appropriate turbulence disturbance at the flow inlet, such as turbulence intensity, the numerical results showed good agreement in (1) the velocity profile downstream of the pipe expansion measured by particle image velocimetry (PIV) in this study, and (2) the reattachment length reported in the literature. Using this validated model, the numerical simulation of flow and component mixing was performed for the FPVS with an inlet flowrate corresponding to a Reynolds number of 2,400 to investigate the flow behavior and component mixing. Standard deviation of the mass fraction of component, referred to as absolute mixing index (AMI), was calculated to quantify the mixing of components, showing that the component mixing length determined by AMI is related to the reattachment length downstream of the expansion. Additionally, this study was able to identify the best location in the FPVS test facility to measure the velocity and concentration profies where multicomponent flow is well mixed.

15 GEOTHERMAL ENERGY↗

DOE OSTI2026

Models and Algorithms for Equilibrium Analysis of Mixed-Material Nucleic Acid Systems

Dynamic programming algorithms within the NUPACK software suite enable analysis of equilibrium base-pairing properties for complex and test tube ensembles containing arbitrary numbers of interacting nucleic acid strands. Currently, calculations are limited to single-material systems that are either all-RNA or all-DNA. Here, to enable analysis of mixed-material systems that are critical for modern applications in vitro, in situ, and in vivo, we develop physical models and dynamic programming algorithms that allow the material of the system to be specified at nucleotide resolution. Free energy parameter sets are constructed for both RNA/DNA and RNA/2'OMe-RNA mixed-material systems by combining available empirical mixed-material parameters with single-material parameter sets to enable treatment of the full complex and test tube ensembles. New dynamic programming recursions account for the material of each nucleotide throughout the recursive process. For a complex with N nucleotides, the mixed-material dynamic programming algorithms maintain the O(N 3 ) time complexity of the single-material algorithms, enabling efficient calculation of diverse physical quantities over complex and test tube ensembles (e.g., complex partition function, equilibrium complex concentrations, equilibrium base-pairing probabilities, minimum free energy secondary structure(s), and Boltzmann-sampled secondary structures) at a cost increase of roughly 2.0-3.5×. The results of existing single-material algorithms are exactly reproduced when applying the new mixed-material algorithms to single-material systems. Accuracy is significantly enhanced using mixed-material models and algorithms to predict RNA/DNA and RNA/2'OMe-RNA duplex melting temperatures from the experimental literature as well as RNA/DNA melt profiles from new experiments. In conclusion, mixed-material analyses can be performed online using the NUPACK web app (www.nupack.org) or locally using the NUPACK Python module.

2′OMe-RNA↗

DOE OSTI2026

Geologic stress modulates fluid mixing at fracture intersections

Fracture intersections are critical links that enable flow and transport in subsurface fracture networks, and their behavior strongly influences fluid mixing in a network. Although all subsurface fractures are subjected to geological stress, we lack a fundamental understanding of how fracture intersection geometry evolves under stress and how these changes influence fluid mixing. Here, we combine 3D printing, 3D X-ray tomographic imaging, and 3D pore-scale numerical simulations to reveal stress-induced changes in intersection geometry and their impact on mixing. Mixing is found to be strongly affected by partial closure of an intersection under stress. As an intersection closes, the void area for fluid flow and diffusion decreases leading to substantial deviations between conventional mixing models and full pore-scale modeling. To address this, we propose a modified mixing model that accounts for intersection deformation, which is essential for accurate modeling of solute transport and mixing through fracture networks.

15 GEOTHERMAL ENERGY↗

DOE OSTI2024

Assessment of the literature about Be-W mixed material layer formation in the fusion reactor environment

Abstract All plasma facing surfaces in a fusion reactor, whether initially pure or an alloy, will rapidly evolve into a mixed material due to plasma-induced erosion, migration and redeposition. Beryllium (Be) erosion from the main chamber, and its transport and deposition on to a tungsten (W) divertor results in the growth of mixed Be-W layers, which can evolve to form beryllides. These Be-W mixed materials exhibit generally less desirable properties than pure tungsten or pure beryllium, such as lower melting points. In order to better understand the parameter space for growth of these alloys, this paper reviews the literature on Be-W mixed material formation experiments—in magnetically confined fusion reactors, in linear plasma test stands, and during thin-film deposition—and on computational modeling of Be-W interactions, as well as briefly assesses the Be-W growth kinetics. We conclude that the following kinetic steps drive the material mixing: adsorption of the implanted/deposited ion on the metal surface; diffusion of the implanted/deposited ion from surface into the bulk, which is accelerated by defects; and loss of deposited material through erosion. Adsorption dominates (or prevents) material mixing in thin-film deposition experiments, whereas diffusion drives material mixing in plasma exposures due to the energetic ion implantation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Improving the representation of isopycnal mixing in E3SM (Final Report)

The DOE-supported Energy Exascale Earth System Model (E3SM) is a major attempt by the Department of Energy to develop a new Earth System Model using unstructured grids, which allow for the model to concentrate resolution where it is most needed and avoid numerical artifacts experienced by previous Earth System Models (ESMs) using regular grids in which the lateral resolution becomes extremely fine as the resolution shrinks at the poles. This effort required rewriting many of the algorithms previously developed for regular grids. One algorithm that was problematic in the first version of E3SM was the representation of mixing due to turbulent ocean eddies with scales smaller than the model grid. These mesoscale eddies are the primary way by which tracers are stirred along density surfaces in the ocean interior. In particular, the process of isopycnal mixing exchanges fresher, more oxygenated waters from polar regions with saltier, nutrient-rich and oxygen poor waters from tropical regions. Previous work in our group has shown that this mechanism is vitally important for bringing oxygen into poorly ventilated tropical regions (Gnanadesikan et al., 2012, 2013) and plays a significant role in determining the uptake of anthropogenic carbon dioxide (Gnanadesikan et al., 2015). However, in E3SMv1 this process was turned off, as turning it on caused the model to become unstable. Additionally, the rate of isopycnal mixing is directly proportional to a mixing coefficient $A_{Redi}$ whose value varies from less than 400 m 2 /s to 2000 m 2 /s across contemporary climate models. The main thrusts of this proposal therefore were to: 1.) Improve the numerics of isopycnal mixing in E3SM; 2.) Identify ways in which climate models are sensitive to the isopycnal mixing coefficient; and 3.) Develop better representations of the isopycnal mixing coefficient.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Assessment of the literature about Be-W mixed material layer formation in the fusion reactor environment

All plasma facing surfaces in a fusion reactor, whether initially pure or an alloy, will rapidly evolve into a mixed material due to plasma-induced erosion, migration and redeposition. Beryllium (Be) erosion from the main chamber, and its transport and deposition on to a tungsten (W) divertor results in the growth of mixed Be-W layers, which can evolve to form beryllides. These Be-W mixed materials exhibit generally less desirable properties than pure tungsten or pure beryllium, such as lower melting points. In order to better understand the parameter space for growth of these alloys, this paper reviews the literature on Be-W mixed material formation experiments—in magnetically confined fusion reactors, in linear plasma test stands, and during thin-film deposition—and on computational modeling of Be-W interactions, as well as briefly assesses the Be-W growth kinetics. We conclude that the following kinetic steps drive the material mixing: adsorption of the implanted/deposited ion on the metal surface; diffusion of the implanted/deposited ion from surface into the bulk, which is accelerated by defects; and loss of deposited material through erosion. Adsorption dominates (or prevents) material mixing in thin-film deposition experiments, whereas diffusion drives material mixing in plasma exposures due to the energetic ion implantation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Aligning Standards Communities for Omics Biodiversity Data: Sustainable Darwin Core-MIxS Interoperability

The standardization of data, encompassing both primary and contextual information (metadata), plays a pivotal role in facilitating data (re-)use, integration, and knowledge generation. However, the biodiversity and omics communities, converging on omics biodiversity data, have historically developed and adopted their own distinct standards, hindering effective (meta)data integration and collaboration. In response to this challenge, the Task Group (TG) for Sustainable DwC-MIxS Interoperability was established. Convening experts from the Biodiversity Information Standards (TDWG) and the Genomic Standards Consortium (GSC) alongside external stakeholders, the TG aimed to promote sustainable interoperability between the Minimum Information about any (x) Sequence (MIxS) and Darwin Core (DwC) specifications. To achieve this goal, the TG utilized the Simple Standard for Sharing Ontology Mappings (SSSOM) to create a comprehensive mapping of DwC keys to MIxS keys. This mapping, combined with the development of the MIxS-DwC extension, enables the incorporation of MIxS core terms into DwC-compliant metadata records, facilitating seamless data exchange between MIxS and DwC user communities. Through the implementation of this translation layer, data produced in either MIxS- or DwC-compliant formats can now be efficiently brokered, breaking down silos and fostering closer collaboration between the biodiversity and omics communities. To ensure its sustainability and lasting impact, TDWG and GSC have both signed a Memorandum of Understanding (MoU) on creating a continuous model to synchronize their standards. These achievements mark a significant step forward in enhancing data sharing and utilization across domains, thereby unlocking new opportunities for scientific discovery and advancement.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS2005

The Closure of the Ocean Mixed Layer Temperature Budget using Level-Coordinate Model Fields

Entrainment is an important element of the mixed layer mass, heat, and temperature budgets. Conventional procedures to estimate entrainment heat advection often do not permit the closure of heat and temperature budgets because of inaccuracies in its formulation. In this study a rigorous approach to evaluate the effect of entrainment using the output of a general circulation model (GCM) that does not have an explicit prognostic mixed layer model is described. The integral elements of the evaluation are 1) the rigorous estimates of the temperature difference between mixed layer water and entrained water at each horizontal grid point, 2) the formulation of the temperature difference such that the budget closes over a volume greater than one horizontal grid point, and 3) the apparent warming of the mixed layer during the mixed layer shoaling to account for the weak vertical temperature gradient within the mixed layer. This evaluation of entrainment heat advection is compared with the estimates by other commonly used ad hoc formulations by applying them in three regions: the north-central Pacific, the Kuroshio Extension, and the Nino-3 areas in the tropical Pacific. In all three areas the imbalance in the mixed layer temperature budget by the ad hoc estimates is significant, reaching a maximum of about 4 K yr(exp -1).

temperature gradient↗