Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flux optimization”

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 145 records · Page 8

Forecasting Lightning Threat Using WRF Proxy Fields

Objectives: Given that high-resolution WRF forecasts can capture the character of convective outbreaks, we seek to: 1. Create WRF forecasts of LTG threat (1-24 h), based on 2 proxy fields from explicitly simulated convection: - graupel flux near -15 C (captures LTG time variability) - vertically integrated ice (captures LTG threat area). 2. Calibrate each threat to yield accurate quantitative peak flash rate densities. 3. Also evaluate threats for areal coverage, time variability. 4. Blend threats to optimize results. 5. Examine sensitivity to model mesh, microphysics. Methods: 1. Use high-resolution 2-km WRF simulations to prognose convection for a diverse series of selected case studies. 2. Evaluate graupel fluxes; vertically integrated ice (VII). 3. Calibrate WRF LTG proxies using peak total LTG flash rate densities from NALMA; relationships look linear, with regression line passing through origin. 4. Truncate low threat values to make threat areal coverage match NALMA flash extent density obs. 5. Blend proxies to achieve optimal performance 6. Study CAPS 4-km ensembles to evaluate sensitivities.

McCaul, E. W., Jr.↗

Evaluating the design of an earth radiation budget instrument with system simulations. Part 2: Minimization of instantaneous sampling errors for CERES-I

Much of the new record of broadband earth radiation budget satellite measurements to be obtained during the late 1990s and early twenty-first century will come from the dual-radiometer Clouds and Earth's Radiant Energy System Instrument (CERES-I) flown aboard sun-synchronous polar orbiters. Simulation studies conducted in this work for an early afternoon satellite orbit indicate that spatial root-mean-square (rms) sampling errors of instantaneous CERES-I shortwave flux estimates will range from about 8.5 to 14.0 W/m on a 2.5 deg latitude and longitude grid resolution. Rms errors in longwave flux estimates are only about 20% as large and range from 1.5 to 3.5 W/sq m. These results are based on an optimal cross-track scanner design that includes 50% footprint overlap to eliminate gaps in the top-of-the-atmosphere coverage, and a 'smallest' footprint size to increase the ratio in the number of observations lying within to the number of observations lying on grid area boundaries. Total instantaneous measurement error also depends on the variability of anisotropic reflectance and emission patterns and on retrieval methods used to generate target area fluxes. Three retrieval procedures from both CERES-I scanners (cross-track and rotating azimuth plane) are used. (1) The baseline Earth Radiaton Budget Experiment (ERBE) procedure, which assumes that errors due to the use of mean angular dependence models (ADMs) in the radiance-to-flux inversion process nearly cancel when averaged over grid areas. (2) To estimate N, instantaneous ADMs are estimated from the multiangular, collocated observations of the two scanners. These observed models replace the mean models in computation of satellite flux estimates. (3) The scene flux approach, conducts separate target-area retrievals for each ERBE scene category and combines their results using area weighting by scene type. The ERBE retrieval performs best when the simulated radiance field departs from the ERBE mean models by less than 10%. For larger perturbations, both the scene flux and collocation methods produce less error than the ERBE retrieval. The scene flux technique is preferable, however, because it involves fewer restrictive assumptions.

Stowe, Larry↗

The Development of Erosion and Impact Resistant Turbine Airfoil Thermal Barrier Coatings

Thermal barrier coatings are used in gas turbine engines to protect engine hot-section components in the harsh combustion environments and extend component lifetimes. For thermal barrier coatings designed for turbine airfoil applications, further improved erosion and impact resistance are crucial for engine performance and durability. Advanced erosion resistant thermal barrier coatings are being developed, with a current emphasis on the toughness improvements using a combined rare earth- and transition metal-oxide doping approach. The performance of the doped thermal barrier coatings has been evaluated in burner rig and laser heat-flux rig simulated engine erosion and thermal gradient environments. The results have shown that the coating composition optimizations can effectively improve the erosion and impact resistance of the coating systems, while maintaining low thermal conductivity and cyclic durability. The erosion and impact damage mechanisms of the thermal barrier coatings will also be discussed.

Zhu, Dongming↗

Comparison Between DSMC and CFD for Hypersonic Planetary Entry Simulations

Hypersonic planetary entry flows span a wide range of Knudsen numbers between rarefied and continuum flows. While computational fluid dynamics (CFD) techniques cannot provide an accurate solution for flows in the rarefied regime, the direct simulation Monte Carlo (DSMC) method is capable of providing accurate solutions for flows in both in the rarefied and continuum regimes but becomes prohibitively expensive as the Knudsen number decreases. For the purpose of thermal protection systems (TPS) design and post-flight reconstruction, various selected points along an entry trajectory are often solved using hypersonic solvers. The quantities of interest that are obtained from that exercise are generally surface quantities, such as pressure, heat flux and enthalpy. Then, material response solvers are used to either design the heat shield to an optimal thickness based on a choice of material, or to provide in-depth heating profiles through the material at various select locations, and compare with flight instrumentation such as the ones that flew on NASA’s two most recent Mars missions, MSL and Mars2020. While most of the heating is generally experienced during the continuum part of the entry, the heating within the rarefied regime is significant for some atmospheres, and hence the flow solutions need to be computed using the DSMC method. Ensuring consistency between hypersonic CFD and the DSMC is crucial so that reliable surface quantities can be passed to material response solvers. Studies were performed to compare the two methods at various select locations, for both non-reacting argon flows as well reacting CO2/N2 flows. Preliminary conclusions show that, for non-reacting flows, the agreement between the two methods for surface heating is excellent (within expected uncertainties) for a freestream Knudsen number of 0.0006, and gets progressively worse as the Knudsen number increases to 0.06. Continuum breakdown analyses were performed and showed that, in general, the Gradient Length Local Knudsen number (KnGLL) from Boyd [1] and associated criteria (KnGLL > 0.05) can be accurately used to predict zones of breakdown in the flow, and associated errors for surface quantities. Furthermore, when studying reacting flow, our studies showed that while flow quantities are highly dependent on relaxation parameters and chemistry rates, it is possible to obtain a good agreement for surface heating, as long as the continuum breakdown in minimal (see Fig. 1). Finally, we studied the impact of the grid refinement on the surface quantities in DSMC.

DSMC↗

Comparison Between DSMC and CFD for Hypersonic Planetary Entry Simulations

Hypersonic planetary entry flows span a wide range of Knudsen numbers between rarefied and continuum flows. While computational fluid dynamics (CFD) techniques cannot provide an accurate solution for flows in the rarefied regime, the direct simulation Monte Carlo (DSMC) method is capable of providing accurate solutions for flows in both in the rarefied and continuum regimes but becomes prohibitively expensive as the Knudsen number decreases. For the purpose of thermal protection systems (TPS) design and post-flight reconstruction, various selected points along an entry trajectory are often solved using hypersonic solvers. The quantities of interest that are obtained from that exercise are generally surface quantities, such as pressure, heat flux and enthalpy. Then, material response solvers are used to either design the heat shield to an optimal thickness based on a choice of material, or to provide in-depth heating profiles through the material at various select locations, and compare with flight instrumentation such as the ones that flew on NASA’s two most recent Mars missions, MSL and Mars2020. While most of the heating is generally experienced during the continuum part of the entry, the heating within the rarefied regime is significant for some atmospheres, and hence the flow solutions need to be computed using the DSMC method. Ensuring consistency between hypersonic CFD and the DSMC is crucial so that reliable surface quantities can be passed to material response solvers. Studies were performed to compare the two methods at various select locations, for both non-reacting argon flows as well reacting CO2/N2 flows. Preliminary conclusions show that, for non-reacting flows, the agreement between the two methods for surface heating is excellent (within expected uncertainties) for a freestream Knudsen number of 0.0006, and gets progressively worse as the Knudsen number increases to 0.06. Continuum breakdown analyses were performed and showed that, in general, the Gradient Length Local Knudsen number (KnGLL) from Boyd and associated criterion (KnGLL > 0.05) seems conservative in predicting zones of breakdown in the flow, and associated errors for surface quantities. Updated criteria of KnGLL = 2.0 and 0.5 appear to be more appropriate for surface and flow quantities, respectively. Furthermore, when studying reacting flows, our studies showed that while flow quantities are highly dependent on relaxation parameters and chemistry rates, it is possible to obtain a good agreement for surface heating, as long as the continuum breakdown is minimal.

DSMC↗

Occultation Modeling for Radiation Obstruction Effects on Spacecraft Systems

A geometric occultation model has been developed to determine line-of-sight obstruction of radiation sources expected for different NASA space exploration mission designs. Example applications includes fidelity improvements for surface lighting conditions, radiation pressure, thermal and power subsystem modeling. The model makes use of geometric two dimensional shape primitives to most effectively model space vehicles. A set of these primitives is used to represent three dimensional obstructing objects as a two dimensional outline from the perspective of an observing point of interest. Radiation sources, such as the Sun or a Moon's albedo is represented as a collection of points, each of which is assigned a flux value to represent a section of the radiation source. Planetary bodies, such as a Martian moon, is represented as a collection of triangular facets which are distributed in spherical height fields for optimization. These design aspects and the overall model architecture will be presented. Specific uses to be presented includes a study of the lighting condition on Phobos for a possible future surface mission, and computing the incident flux on a spacecraft's solar panels and radiators from direct and reflected solar radiation subject to self-shadowing or shadowing by third bodies.

de Carufel, Guy↗

Transition-Edge Sensor Pixel Parameter Design of the Microcalorimeter Array for the X-Ray Integral Field Unit on Athena

The focal plane of the X-ray integral field unit (X-IFU) for ESA's Athena X-ray observatory will consist of approximately 4000 transition edge sensor (TES) x-ray microcalorimeters optimized for the energy range of 0.2 to 12 kiloelectronvolts. The instrument will provide unprecedented spectral resolution of approximately 2.5 electronvolts at energies of up to 7 kiloelectronvolts and will accommodate photon fluxes of 1 milliCrab (90 counts per second) for point source observations. The baseline configuration is a uniform large pixel array (LPA) of 4.28 arcseconds pixels that is read out using frequency domain multiplexing (FDM). However, an alternative configuration under study incorporates an 18 by × 18 small pixel array (SPA) of 2 arcseconds pixels in the central approximately 36 arcseconds region. This hybrid array configuration could be designed to accommodate higher fluxes of up to 10 milliCrabs (900 counts per second) or alternately for improved spectral performance (less than 1.5 electronvolts) at low count-rates. In this paper we report on the TES pixel designs that are being optimized to meet these proposed LPA and SPA configurations. In particular we describe details of how important TES parameters are chosen to meet the specific mission criteria such as energy resolution, count-rate and quantum efficiency, and highlight performance trade-offs between designs. The basis of the pixel parameter selection is discussed in the context of existing TES arrays that are being developed for solar and x-ray astronomy applications. We describe the latest results on DC biased diagnostic arrays as well as large format kilo-pixel arrays and discuss the technical challenges associated with integrating different array types on to a single detector die.

Smith, S. J.↗

Optimizing a Coupled Solution for the Orion MPCV Using DPLR and NEQAIR

The goal of this project is to study the convective and radiative heat flux predictions using computational fluid dynamics (CFD) on a smooth body representation of the NASA Orion Multi-Purpose Crew Vehicle (MPCV). The Orion MPCV is designed for long duration missions in space and will play a crucial role in transporting astronauts in the upcoming Artemis missions. With Orion's aim for cost-effective reusability comes a need to study its Thermal Protection System (TPS), which allows for safe entry into the atmospheres of the Earth, Moon, and Mars. In order to study the effects of the TPS on test flights and better understand the aerothermal environment, CFD is a necessary tool used to create simulations of entry vehicles and calculate the heat load and heat flux. However, solutions are often performed using an uncoupled approach between the flow and the radiation solver. This leads to an overestimation of the overall heating effects because they do not consider radiative cooling. Furthermore, the effects of radiation on the aft-body is an area of research with significant uncertainty. Through the use of codes developed at NASA Ames, DPLR and NEQAIR, a coupled solution can be attained for the Orion capsule which will more accurately solve for the effects of radiation in both the fore and aft-body regions. This yields more accurate design margins for the TPS, which produces a better estimate of the thickness required. By developing an accurate coupled solution of the TPS on the Orion MPCV, planners will have greater flexibility in carrying out future missions to the Moon and beyond.

Design↗

A graphite crystal polarimeter for stellar X-ray astronomy.

The first crystal X-ray polarimeter to be used for X-ray astronomy is described. Polarization is measured by modulation of the X rays diffracted at an average 45 deg glancing angle from large, curved graphite crystal panels as these rotate about an axis parallel to the incident X-ray flux. Arrangement of the crystal panels, the design of the detector, and the signal-processing circuitry were optimized to minimize systematic effects produced by off-axis pointing of the rocket and cosmic ray induced events. The in-flight performance of the instrument in relation to the observed background signal is discussed.

Weisskopf, M. C.↗

Comparison of MERRA-2 and ECCO-V4 Ocean Surface Heat Fluxes: Consequences of Different Forcing Feedbacks on Ocean Circulation and Implications for Climate Data Assimilation

Forcing ocean models with reanalysis data is a common practice in ocean modeling. As part of this practice, prescribed atmospheric state variables and interactive ocean SST (Sea Surface Temperature) are used to calculate fluxes between the ocean and the atmosphere. When forcing an ocean model with reanalysis fields, errors in the reanalysis data, errors in the ocean model and errors in the forcing formulation will generate a different solution compared to other ocean reanalysis solutions (which also have their own errors). As a first step towards a consistent coupled ocean-atmosphere reanalysis, we compare surface heat fluxes from a state-of-the-art atmospheric reanalysis, the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2), to heat fluxes from a state-of-the-art oceanic reanalysis, the Estimating the Circulation and Climate of the Ocean Version 4, Release 2 (ECCO-v4). Then, we investigate the errors associated with the MITgcm (Massachusetts Institute of Technology general circulation model) ocean model in its ECCO-v4 ocean reanalysis configuration (1992-2011) when it is forced with MERRA- 2 atmospheric reanalysis fields instead of with the ECCO-v4 adjoint optimized ERA-interim state variables. This is done by forcing ECCO-v4 ocean with and without feedbacks from MERRA-2 related to turbulent fluxes of heat and moisture and the outgoing long wave radiation. In addition, we introduce an intermediate forcing method that includes only the feedback from the interactive outgoing long wave radiation. The resulting ocean circulation is compared with ECCO-v4 reanalysis and in-situ observations. We show that, without feedbacks, imbalances in the energy and the hydrological cycles of MERRA-2 (which are directly related to the fact it was created without interactive ocean) result in considerable SST drifts and a large reduction in sea level. The bulk formulae and interactive outgoing long wave radiation, although providing air-sea feedbacks and reducing model-data misfit, strongly relax the ocean to observed SST and may result in unwanted features such as large change in the water budget. These features have implications in a desired forcing recipe to be used. The results strongly and unambiguously argue for next generation data assimilation climate studies to involve fully coupled systems.

MERRA-2↗

Optimizing Crystal Volume for Neutron Diffraction

Neutron diffraction is uniquely sensitive to hydrogen positions and protonation state. In that context structural information from neutron data is complementary to that provided through X-ray diffraction. However, there are practical obstacles to overcome in fully exploiting the potential of neutron diffraction, Le. low flux and weak scattering. Several approaches are available to overcome these obstacles and we have investigated the simplest: increasing the diffracting volume of the crystals. Volume is a quantifiable metric that is well suited for experiment design and optimization techniques. By using response surface methods we have optimized xylose isomerase crystal volume, enabling neutron diffraction while we determined the crystallization parameters with the minimum of experiments. Our results suggest a systematic means of enabling neutron diffraction studies for a larger number of samples that require information on hydrogen position and/or protonation state.

Snell, Edward H.↗

Comparative Wide Temperature Core Loss Characteristics of Two Candidate Ferrites for the NASA/TRW 1500 W PEBB Converter

High frequency core loss and magnetization properties of commercial type MN8CX and PC40, high resistivity, MnZn based, power ferrites are presented over the temperature range of -l50 C to 150 C, at selected values of peak flux density (B (sub p)). Most of the data is at 100 kHz, with some data extended to 200 and 300 kHz for the MN8CX. Plots of the specific Core loss against temperature exhibit the minimal characteristic of such ferrites. These plots show that the MN8CX is optimized for minimum loss at about 25 C, whereas the PC40 is optimized at about 80 C. At the points of minimum loss and for the same B (sub p), the MN8CX has roughly half the losses of the PC40 at the lower flux densities. This loss ratio continues down to cryogenic temperatures. However, above about 80 C the losses are practically equal. The lowest 100 kHz loss recorded, 50 mW/cm3 for the MNGCX at a B (sub p) of 0.1T, equals that of a very low loss, Co based, transverse magnetically annealed, amorphous ribbon material. Except possibly at lower B (sub p) or much higher frequencies, these ferrites are not competitive for low losses over a wide temperature range with certain specialty amorphous materials. Permeability is computed from a linear model, plots against temperature are presented and again compared to the specialty amorphous materials.

Niedra, Janis M.↗

Thermoelectric System Economics - The Apex: New Paradigms in Manufacturing and Interface Performance Relationships Driving System Cost Optimizations

Thermoelectric energy recovery systems require cost-optimized designs with high thermal and thermoelectric (TE) performance to surmount common commercialization barriers controlling the acceptance of thermoelectric generator (TEG) systems. Recent work has addressed this requirement with new analytic tools and paradigms that allow integrated TEG cost-performance analysis and optimization. New work herein has identified the existence of optimized TE element design (i.e., TE element length) for minimizing TEG cost and provides design guidelines associated with critical thermal and electrical contact resistance effects and TE manufacturing sensitivities. Once optimum hot-side heat flux criteria are satisfied, this work highlights the tradeoffs between interface resistance effects and manufacturing cost effects in determining optimum-cost TE element lengths, and then quantifies the increase in optimum-cost TE lengths as manufacturing costs become more sensitive to TE element lengths. Optimized TE element design criteria and key non-dimensional design parameters (i.e., [𝐶′′′LTE/𝐶′′], [TE/con], [con/(TE,aveLTE)]) driving the design are presented and explained showing the interdependencies between manufacturing cost parameters and thermal and electrical interface impacts on integrated TE cost-performance analysis and design. Comprehensive TEG system cost-performance relationships clearly demonstrate the two parameters that most greatly impact the TEG system cost are heat exchanger costs ($/(W/K)) and TE hot side heat flux. Reducing heat exhanger cost from $1/(W/K) to $0.5/(W/K) can decrease TEG system cost by $2/W to $4/W depending on TE hot side heat fluxes from 4-10 W/cm2. The impact of TE hot side heat flux on TEG system cost is larger; reducing TE system cost by factors of 2-4 as TE hot side heat flux increases from 4 W/cm2 to high levels around 18 W/cm2. The levels of thermal contact resistance and electrical contact resistance required to even approach the critical TEG system cost level of $1/W are characterized by con = 1.0 x 10-10 -m2 and r = 0.1, respectively. This not only lowers TEG costs, but it also makes the TEG system costs less sensitive to TE manufacturing dependences on TE element length. New design paradigms and relationships create holistic, integrated TE performance-cost models that enhance understanding of crucial interrelationships between component costs, TE design parameters and material properties, heat exchanger design parameters, interfacial heat flux, and now thermal and electrical interface effects in minimizing TEG system costs.

Hendricks, Terry J.↗

An engineering code to analyze hypersonic thermal management systems

Thermal loads on current and future aircraft are increasing and as a result are stressing the energy collection, control, and dissipation capabilities of current thermal management systems and technology. The thermal loads for hypersonic vehicles will be no exception. In fact, with their projected high heat loads and fluxes, hypersonic vehicles are a prime example of systems that will require thermal management systems (TMS) that have been optimized and integrated with the entire vehicle to the maximum extent possible during the initial design stages. This will not only be to meet operational requirements, but also to fulfill weight and performance constraints in order for the vehicle to takeoff and complete its mission successfully. To meet this challenge, the TMS can no longer be two or more entirely independent systems, nor can thermal management be an after thought in the design process, the typical pervasive approach in the past. Instead, a TMS that was integrated throughout the entire vehicle and subsequently optimized will be required. To accomplish this, a method that iteratively optimizes the TMS throughout the vehicle will not only be highly desirable, but advantageous in order to reduce the manhours normally required to conduct the necessary tradeoff studies and comparisons. A thermal management engineering computer code that is under development and being managed at Wright Laboratory, Wright-Patterson AFB, is discussed. The primary goal of the code is to aid in the development of a hypersonic vehicle TMS that has been optimized and integrated on a total vehicle basis.

Vangriethuysen, Valerie J.↗

Preconditioning for the Navier-Stokes equations with finite-rate chemistry

The extension of Van Leer's preconditioning procedure to generalized finite-rate chemistry is discussed. Application to viscous flow is begun with the proper preconditioning matrix for the one-dimensional Navier-Stokes equations. Eigenvalue stiffness is resolved and convergence-rate acceleration is demonstrated over the entire Mach-number range from nearly stagnant flow to hypersonic. Specific benefits are realized at the low and transonic flow speeds typical of complete propulsion-system simulations. The extended preconditioning matrix necessarily accounts for both thermal and chemical nonequilibrium. Numerical analysis reveals the possible theoretical improvements from using a preconditioner for all Mach number regimes. Numerical results confirm the expectations from the numerical analysis. Representative test cases include flows with previously troublesome embedded high-condition-number areas. Van Leer, Lee, and Roe recently developed an optimal, analytic preconditioning technique to reduce eigenvalue stiffness over the full Mach-number range. By multiplying the flux-balance residual with the preconditioning matrix, the acoustic wave speeds are scaled so that all waves propagate at the same rate, an essential property to eliminate inherent eigenvalue stiffness. This session discusses a synthesis of the thermochemical nonequilibrium flux-splitting developed by Grossman and Cinnella and the characteristic wave preconditioning of Van Leer into a powerful tool for implicitly solving two and three-dimensional flows with generalized finite-rate chemistry. For finite-rate chemistry, the state vector of unknowns is variable in length. Therefore, the preconditioning matrix extended to generalized finite-rate chemistry must accommodate a flexible system of moving waves. Fortunately, no new kind of wave appears in the system. The only existing waves are entropy and vorticity waves, which move with the fluid, and acoustic waves, which propagate in Mach number dependent directions. The nonequilibrium vibrational energies and species densities in the unknown state vector act strictly as convective waves. The essential concept for extending the preconditioning to generalized chemistry models is determining the differential variables which symmetrize the flux Jacobians. The extension is then straight-forward. This algorithm research effort will be released in a future version of the production level computational code coined the General Aerodynamic Simulation Program (GASP), developed by Walters, Slack, and McGrory.

Godfrey, Andrew G.↗

Efficient Implementations of the Quadrature-Free Discontinuous Galerkin Method

The efficiency of the quadrature-free form of the dis- continuous Galerkin method in two dimensions, and briefly in three dimensions, is examined. Most of the work for constant-coefficient, linear problems involves the volume and edge integrations, and the transformation of information from the volume to the edges. These operations can be viewed as matrix-vector multiplications. Many of the matrices are sparse as a result of symmetry, and blocking and specialized multiplication routines are used to account for the sparsity. By optimizing these operations, a 35% reduction in total CPU time is achieved. For nonlinear problems, the calculation of the flux becomes dominant because of the cost associated with polynomial products and inversion. This component of the work can be reduced by up to 75% when the products are approximated by truncating terms. Because the cost is high for nonlinear problems on general elements, it is suggested that simplified physics and the most efficient element types be used over most of the domain.

Lockard, David P.↗

Salt Pill Design and Fabrication for Adiabatic Demagnetization Refrigerators

The performance of an adiabatic demagnetization refrigerator (ADR) is critically dependent on the design and construction of the salt pills that produce cooling. In most cases, the primary goal is to obtain the largest cooling capacity at the low temperature end of the operating range. The realizable cooling capacity depends on a number of factors, including refrigerant mass, and how efficiently it absorbs heat from the various instrument loads. The design and optimization of "salt pills" for ADR systems depend not only on the mechanical, chemical and thermal properties of the refrigerant, but also on the range of heat fluxes that the salt pill must accommodate. Despite the fairly wide variety of refrigerants available, those used at very low temperature tend to be hydrated salts that require a dedicated thermal bus and must be hermetically sealed, while those used at higher temperature - greater than about 0.5 K - tend to be single-­‐ or poly-­‐crystals that have much simpler requirements for thermal and mechanical packaging. This paper presents a summary of strategies and techniques for designing, optimizing and fabricating salt pills for both low-­‐ and mid-­‐temperature applications.

sub-kelvin cooling↗

Enhancement of Nucleate Boiling Heat Flux on Macro/Micro-Structured Surfaces Cooled by Multiple Impinging Jets

An experimental investigation of nucleate boiling heat transfer from modified surfaces cooled by multiple in-line impinging circular jets is reported and found to agree with single jet results. A copper block is heated from the back by two electrical arcs, and cooled on the opposite side by three identical liquid jets of distilled water at subcoolings of 25 C 50 C and 77 C and Freon 113 at 24 C subcooling. Liquid flow rates are held constant at 5, 10, and 15 GPH for each of the three jets with jet velocities ranging from 1.4 m/s to 1 1.2 m/s and jet diameters from 0.95 mm to 2.2 mm. To increase the maximum heat flux (CHF) and heat removal rate, the boiling surface was modified by both macro and micro enhancements. Macro modification consists of machined radial grooves in the boiling surface arranged in an optimally designed pattern to allow better liquid distribution along the surface. These grooves also reduce splashing of liquid droplets, and provide 'channels' to sweep away bubbles. Micro modification was achieved by flame spraying metal powder on the boiling surface, creating a porous, sintered surface. With the addition of both micro and macro structured enhancements, maximum heat flux and nucleate boiling can be enhanced by more than 200%. Examination of each surface modification separately and together indicates that at lower superheats, the micro structure provides the enhanced heat transfer by providing more nucleation sites, while for higher superheats the macro structure allows better liquid distribution and bubble removal. A correlation is presented to account for liquid subcoolings and surface enhancements, in addition to the geometrical and fluid properties previously reported in the literature.

Kugler, Scott Lee↗