Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic mass flow”

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

Fractionation of the Early Terrestrial Atmospheres: Dynamical Escape

Hydrodynamic escape may have played a significant role in the early fractionation of the atmospheres of the terrestrial planets. This possibility has been demonstrated in the last two decades by numerous models that show radial, transonic flow of hydrogen can occur in the presence of sufficient solar EUV Hydrodynamic escape may have played a significant role in the early fractionation of the atmospheres of the terrestrial planets. This possibility has been demonstrated in the last two decades by numerous models that show radial, transonic flow of hydrogen can occur in the presence of sufficient solar EUV flux, thought to exist in the first 500 My. The models show that the larger the solar flux the greater the mass of the fractionating species, which are accelerated to escape speeds by the hydrogen wind through drag processes. As the atmospheres evolve and the solar EUV flux wanes, the maximum mass of flowing gas constituents decreases until all gases become static. We show that fractionation can continue beyond this point when non-radial flow and dynamically enhanced Jeans escape are considered. For example, the early terrestrial atmospheres are thought to have had large hydrogen contents, resulting in exobase altitudes of a planetary radius or more. In this case, rotational speeds at the exobases of Earth and Mars would be large enough so that light constituents would "spin" off and fractionate, especially at equatorial latitudes. Also, in the presence of transonic flow of hydrogen only, non-radial expansion throws heavier gases to high altitudes in the exosphere, accompanied by strong bulk speeds at the exobase, which results in enhanced thermal escape fluxes and fractionation. flux, thought to exist in the first 500 My. The models show that the larger the solar flux the greater the mass of the fractionating species, which are accelerated to escape speeds by the hydrogen wind through drag processes. As the atmospheres evolve and the solar EUV flux wanes, the maximum mass of flowing gas constituents decreases until all gases become static. We show that fractionation can continue beyond this point when non-radial flow and dynamically enhanced Jeans escape are considered. For example, the early terrestrial atmospheres are thought to have had large hydrogen contents, resulting in exobase altitudes of a planetary radius or more. In this case, rotational speeds at the exobases of Earth and Mars would be large enough so that light constituents would "spin" off and fractionate, especially at equatorial latitudes. Also, in the presence of transonic flow of hydrogen only, non-radial expansion throws heavier gases to high altitudes in the exosphere, accompanied by strong bulk speeds at the exobase, which results in enhanced thermal escape fluxes and fractionation.

Hartle, Richard E.↗

Analytical Model For Fluid Dynamics In A Microgravity Environment

Report presents analytical approximation methodology for providing coupled fluid-flow, heat, and mass-transfer equations in microgravity environment. Experimental engineering estimates accurate to within factor of 2 made quickly and easily, eliminating need for time-consuming and costly numerical modeling. Any proposed experiment reviewed to see how it would perform in microgravity environment. Model applied in commercial setting for preliminary design of low-Grashoff/Rayleigh-number experiments.

Naumann, Robert J.↗

Modeling scramjet combustor flowfields with a grid adaptation scheme

The accurate description of flow features associated with the normal injection of fuel into supersonic primary flows is essential in the design of efficient engines for hypervelocity aerospace vehicles. The flow features in such injections are complex with multiple interactions between shocks and between shocks boundary layers. Numerical studies of perpendicular sonic N2 injection and mixing in a Mach 3.8 scramjet combustor environment are discussed. A dynamic grid adaptation procedure based on the equilibration of spring-mass system is employed to enhanced the description of the complicated flow features. Numerical results are compared with experimental measurements and indicate that the adaptation procedure enhances the capability of the modeling procedure to describe the flow features associated with scramjet combustor components.

Ramakrishnan, R.↗

Computational analysis of forebody tangential slot blowing on the high alpha research vehicle

Current and future fighter aircraft can maneuver in the high-angle-of-attack flight regime while flying at low subsonic and transonic freestream Mach numbers. However, at any flight speed, the ability of the vertical tails to generate yawing moment is limited in high-angle-of-attack flight. Thus, any system designed to provide the pilot with additional side force and yawing moment must work in both low subsonic and transonic flight. However, previous investigations of the effectiveness of forebody tangential slot blowing in generating the desired control forces and moments have been limited to the low subsonic freestream flow regime. In order to investigate the effectiveness of tangential slot blowing in transonic flight, a computational fluid dynamics analysis was carried out during the grant period. Computational solutions were obtained at three different freestream Mach numbers and at various jet mass flow ratios. All results were obtained using the isolated F/A-18 forebody grid geometry at 30.3 degrees angle of attack. One goal of the research was to determine the effect of freestream Mach number on the effectiveness of forebody tangential slot blowing in generating yawing moment. The second part of the research studied the force onset time lag associated with blowing. The time required for the yawing moment to reach a steady-state value from the onset of blowing may have an impact on the implementation of a pneumatic system on a flight vehicle.

Gee, Ken↗

Sensitivity Study of Dust-Laden Flows in the DLR GBK Facility using the DUST-US3D Framework

An unstructured integrated simulation framework DUST-US3D has been used to study particle-laden flows in the German Aerospace Center (DLR) GBK facility. Air/particle inflow conditions and the overall simulation methodology has been adjusted to ensure consonance with recent experimental work. The resulting particle dynamics are found to be extremely sensitive to the configuration of the particle injection process. However, owing to relatively low particle mass flow rates incorporating 2-way coupling (particles-gas exchange momentum/energy) or particle-particle collisions has no significant impact on predictions. Numerical estimates for particle velocity magnitude are consistently lower than those observed during experiments which could be indicative of particles coalescing into larger, possibly non-spherical aggregates.

Dusty Flows↗

MOSCATO Development and Integration in Fiscal Year 2025: Implementation of Multiphase, Multiphysics Modeling Capabilities for Molten Salt Systems

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy corrosion for molten salt systems. In FY25, significant developments were made to the code package, enhancing its capabilities for modeling all relevant phenomena within operating moltens salt reactors (MSRs). The developments and activities in FY25 included: 1. Implementation of Level-Set methods to enable modeling of single-bubble behavior in molten salts. In FY25, the Level-Set two-phase flow modeling implementation was improved to simulate single bubble behavior with molten salt media. The large density and viscosity ratios between typical gases and molten salt liquids present challenges for these types of numerical solvers. With enhancements to the pressure projection method, MOSCATO’s Level-Set solver was able to be successfully validated to experiments related to helium bubble rise in stagnant molten salt. The simulated bubble rising velocity showed reasonable good agreement with experimental measurements. The bubble shape and dynamics were also visually compared with experimental snapshots, demonstrating a good qualitative match. 2. Generation of mass transfer correlations for multiphase flow systems. To enable calculations of the tritium transport across the interface between gas bubbles and salt, we modeled high- Schmidt-number mass transfer around a sphere across a broad range of Reynolds numbers. The mesh near the sphere surface was highly refined to resolve steep concentration gradients caused by the low diffusion coefficient. Literature-based mass transfer correlations were compared with the numerical results, and modifications were proposed to improve agreement, particularly at higher Schmidt numbers. These mass transfer correlations were subsequently provided to other national laboratories to help enable high quality mass transfer simulations using lower-order solvers under development within the NEAMS program. 3. Preliminary implementation of a bubbly flow solver. To model bubbly flow in molten salt, we implemented a bubbly flow solver for void fractions less than 5%. To do so, an algebraic relative velocity model that assumes small bubbles with rapid momentum equilibration was added to MOSCATO to compute bubble velocities. Preliminary comparisons with experimental data showed reasonable agreement, and further improvements are underway. 4. Generation of mass transfer correlations for MSRE subchannel The Molten-Salt Reactor Experiment (MSRE) was a landmark historical project that demonstrated the feasibility of molten-salt reactor technology. The MSRE campaign also generated a significant body of experimental data and reports that continue to support molten-salt–related research. In this report, the MSRE core subchannel was used as the reference geometry for a mass transfer study performed with MOSCATO. The geometry and computational mesh were adapted from a previous study, providing adequate resolution for the relatively low Reynolds number in this case. Additional mesh refinement was applied to reach higher Schmidt numbers, enabling the derivation of a reliable mass-transfer correlation for the present scenario. 5. Simulations of oxygen ingressions into molten salt. In the previous fiscal year, we initiated a study to simulate oxygen ingression in stagnant salt. As oxygen enters the salt through its surface, it reacts with Ce 3+ to form solid CeO 2 and other reaction products. To more fully capture the complex diffusion-convection-reaction mechanisms, capabilities for modeling natural convection in the salt vessel were added. These were needed as the flow of the ingressed gas induced flow in the salt caused by surface shear and non-isothermal effects. With these updated physics in place, we were able to successfully reproduce the experimental results for the rate of change of CeCl 3 concentrations versus time. 6. Flow corrosion model validation. In FY24, MOSCATO’s corrosion model was validated against static corrosion experiments. In FY25, this work was extended to a flow corrosion experiment, where FLiNaK salt was driven by natural convection, with initial salt impurities to initiate corrosion. Despite uncertainties in parameters such as elemental diffusion coefficients in the alloy and unknown H + concentrations, the simulations achieved good agreement with experimental results, especially in predicting sample mass losses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Constraining Stellar Parameters and Atmospheric Dynamics of the Carbon AGB star V Oph

Molecules and dust produced by the atmospheres of cool evolved stars contribute to a significant amount of the total material found in the interstellar medium. To understand the mechanism behind the mass loss of these stars, it is of pivotal importance to investigate the structure and dynamics of their atmospheres. Our goal is to verify if the extended molecular and dust layers of the carbon-rich asymptotic giant branch (AGB) star V Oph, and their time variations, can be explained by dust-driven winds triggered by stellar pulsation alone, or if other mechanisms are in play. We model V Oph mid-infrared interferometric VLTI-MIDI data (8–13 μm), at phases 0.18, 0.49, and 0.65, together with literature photometric data, using the latest-generation self-consistent dynamic atmosphere models for carbon-rich stars: DARWIN. We determine the fundamental stellar parameters: T(eff) = 2600 K, L(bol) = 3585 L⊙, M = 1.5 M⊙, C/O = 1.35, Mass flow rate= 2.50 x 10(exp -6) M⊙/yr. We calculate the stellar photospheric radii at the three phases: 479, 494, 448 R⊙; and the dust radii: 780, 853, 787 R⊙. The dynamic models can fairly explain the observed N-band visibility and spectra, although there is some discrepancy between the data and the models, which is discussed in the text. We discuss the possible causes of the temporal variations of the outer atmosphere, deriving an estimate of the magnetic field strength, and computing upper limits for the Alfvén waves velocity. In addition, using period–luminosity sequences, and interferometric modeling, we suggest V Oph as a candidate to be reclassified as a semi-regular star.

Gioia Rau↗

Modelling the Centers of Galaxies

The key to studying central regions by means of nobody numerical experiments is to concentrate on the central few parsecs of a galaxy, replacing the remainder of the galaxy by a suitable boundary condition, rather after the manner in which stellar interiors can be studied without a detailed stellar atmosphere by replacing the atmosphere with a boundary condition. Replacements must be carefully designed because the long range gravitational force means that the core region is sensitive to mass outside that region and because particles can exchange between the outer galaxy and the core region. We use periodic boundary conditions, coupled with an iterative procedure to generate initial particle loads in isothermal equilibrium. Angular momentum conservation is ensured for problems including systematic rotation by a circular reflecting boundary and by integrating in a frame that rotates with the mean flow. Mass beyond the boundary contributes to the gravitational potential, but does not participate in the dynamics. A symplectic integration scheme has been developed for rotating coordinate systems. This combination works well, leading to robust configurations. Some preliminary results with this combination show that: (1) Rotating systems are extremely sensitive to non-axisymmetric external potentials, and (2) that a second core, orbiting near the main core (like the M31 second core system), shows extremely rapid orbital decay. The experimental setups will be discussed, along with preliminary results.

Smith, B. F.↗

High Fidelity Simulations of a Methane-Oxygen Rotating Detonation Rocket Engine

The rotating detonation combustor (RDC) is an important realization of pressure gain combustion for practical applications. The RDC system is characterized by a highly unsteady flow field, with multiple reflected pressure waves following the strong detonations and the existence of partially-burnt gases entrained within the detonation annulus during the injector recovery process. While experimental efforts have provided macroscopic properties of an RDC with a given fuel-oxidizer mixture and geometry, limited optical access and reduced lifespan of flow-field diagnostic equipment constrain understanding of the mechanisms that lend to wave stability, controllability, and sustainability. To this end, high-fidelity numerical simulations of the Air Force Research Laboratory (AFRL) methane-oxygen rotating detonation rocket engine (RDRE) are performed to provide detailed insight into the detonation physics and anomalous behavior within the combustor. The results highlight the merit of high-fidelity numerical studies in analyzing the combustion processes and injector dynamics of a complex combustion system. The simulations capture the general wave behavior observed within the experimental studies. Three different operating conditions, with varying global equivalence ratio and system mass flow rate, are considered: (1) baseline configuration, (2) high equivalence ratio, and (3) high mass flow rate. The high mass flow rate case results in increased detonation wave strength whereas the high equivalence ratio results in increased deflagrative heat release. The high equivalence ratio case results in chaotic behavior with two pairs of counter-rotating detonation waves, whereas a single pair of co-rotating waves is observed in the baseline and high mass flow rate cases. Further, parasitic combustion of the reactant mixture ahead of the detonation wave diminishes its strength and reduces propagation speed to approximately 50-60\% of the Chapman-Jouguet (CJ) condition. In this work, a detailed analysis of the detonation wave behavior, injector dynamics, and the statistics of combustion parameters are presented.

Prakash, Supraj↗

An integrated analytical aeropropulsive/aeroelastic model for the dynamic analysis of hypersonic vehicles

The development of an approach to the determination of the dynamic characteristics of hypersonic vehicles which is intentionally generic and basic is given. The approach involves a 2D hypersonic aerodynamic analysis utilizing Newtonian theory, coupled with a 1D aero/thermoanalysis of the flow in a scramjet-type propulsion system. In addition, the airframe is considered to be elastic, and the structural dynamics are characterized in terms of a simple lumped-mass model of the invacuo vibration modes. The vibration modes are coupled to the rigid-body modes through the aero/propulsive forces acting on the structure. The control effectors considered on a generic study configuration include aerodynamic pitch-control surfaces, as well as engine fuel flow and diffuser area ratio. The study configuration is shown to be highly statically unstable in pitch, and to exhibit strong airframe/engine/elastic coupling in the aeroelastic and attitude dynamics, as well as the engine responses.

Chavez, Frank R.↗

A dynamic grid adaptation procedure for mixing and reacting processes in scramjet combustors with injection

The accurate description of flow features associated with the normal injection of fuel into supersonic primary flows is essential in the design of efficient engines for hypervelocity aerospace vehicles. The flow features in such injections are complex with multiple interactions between shocks and between shocks and boundary layers. Numerical studies of perpendicular sonic N2 injection and mixing in a Mach-3.8 stream and of the mixing and reaction effects from normal injection of H2 in a Mach-2.5 scramjet combustor environment are discussed. A dynamic grid adaptation procedure based on equilibration of spring-mass systems is employed to enhance the description of the complicated flow features. Numerical results are compared with experimental measurements, and conclusions are drawn in regard to the capability of the modeling procedure to describe the flow features associated with scramjet combustor components.

Ramakrishnan, R.↗

Gas-Liquid Flow Modeling for Renewable Fuels Production

Aerobic/anaerobic and gas fermentation pathways have emerged as promising new technologies for the generation of renewable fuels/chemicals from biomass derived sugars, and mixtures of greenhouse/energy rich gas streams (CO2/CH4/H2/CO) via microbial action. Example pathways include sugars-to-ethanol conversion, biomethanation (CO2/H2 to CH4), biogas upgrading, CO fermentation and wet-waste conversion. Gas and liquid phase transport, mass-transfer, and mixing physics at large length scales can significantly affect microbial conversion rates, particularly when the microbial reaction requires a narrow set of conditions. These phenomena are difficult to study in small-scale bench-top reactors that are typically well-mixed. Predictive computational fluid dynamics (CFD) based simulations can therefore aid in the scale-up, design and optimization of these reactors. This work presents multiphase Euler-Euler CFD simulations of at-scale (~500 m3) bioreactors. Our mathematical model treats the gas and liquid as interpenetrating phases. This approach reduces the computational complexity of tracking individual gas bubbles that are several orders of magnitude smaller than reactor dimensions. We solve the Reynolds averaged Navier-Stokes (RANS) multiphase equations that account for phase and chemical species transport, interphase mass and momentum transfer and uses a phenomenological model for gas uptake by microbes. We use a customized solver derived from open-source CFD toolbox, OpenFOAM [1], to perform these simulations, which has been validated against small-scale reactors in our previous work [2]. There is currently a knowledge-gap regarding bubble-size distributions when using gas mixtures with vastly different properties, which can have a significant impact overall mass-transfer. For example, hydrogen bubbles are more buoyant compared to other relatively heavier gases (CO2/CH4/CO), resulting in a large distribution of residence times and bubble sizes. This work therefore develops a deeper understanding of bubble dynamics and interphase mass transfer in such heterogenous gas mixtures through well-resolved computational models. We use a population balance model (PBM) for bubble-size-distribution modeling that is validated against small-scale experiments in our solver with an uncertainty quantification study for bubble coalescence and break-up model parameters. Results pertaining to multiple simulations of gas-fermentation reactors are presented where gas mixtures with varying compositions of CO2/CH4/CO/H2 are imposed at the sparger boundaries. The spatio-temporal variations in bubble-size distribution and mass transfer coefficient are analyzed for varying superficial velocities and gas-compositions for varying sizes of bubble-column and airlift reactors. This work will also examine the performance of different reactor designs, viz. bubble column reactor, airlift reactor with an internal draft tube, and a stirred-tank reactor with Rushton impellers. Reactor mass-transfer coefficient, gas hold-up, and dissolved gas distribution are critically analyzed among reactors, and sensitivity studies pertaining to gas flow rates and reactor geometry will be presented. [1] Weller, H., Tabor, G., Jasak, H. and Fureby, C., A tensorial approach to computational continuum mechanics using object-oriented techniques, Computers in physics, 12, 6, 620--631, 1998. [2] Rahimi, M., Sitaraman, H., Humbird, D. and Stickel, J., Computational fluid dynamics study of full-scale aerobic bioreactors: Evaluation of gas-liquid mass transfer, oxygen uptake, and dynamic oxygen distribution, Chemical Engineering Research and Design, 139: 283-295.

BIOMASS FUELS↗

Investigation of Atmospheric Boundary-Layer Effects on Launch-Vehicle Ground Wind Loads

A launch vehicle ground-wind-loads program is underway at the NASA Langley Transonic Dynamics Tunnel. The objectives are to quantify key aerodynamic and structural characteristics that impact the occurrence of large wind-induced oscillations of a launch vehicle when exposed to ground winds prior to launch. Of particular interest is the dynamic response of a launch vehicle when a von Kármán vortex street forms in the wake of the vehicle resulting in quasiperiodic lift and drag forces. Vehicle response to these quasiperiodic forces can become quite large when the frequency of vortex shedding nears that of a lowly-damped structural mode thereby exciting a resonant response. Wind approaching the vehicle can be characterized by a varying speed with height and turbulence content. The combination of both the varying speed and turbulence content is referred to herein as the atmospheric boundary-layer. The importance of the atmospheric boundary-layer upon launch vehicle wind-induced oscillation response has long been questioned, and its effects are not well understood. Although there are several facilities around the world dedicated to replicating atmospheric boundary layers, the development of such a boundary layer in a wind tunnel capable of producing flight-representative Reynolds numbers for aeroelastically-scaled launch vehicle models has only recently been accomplished. The NASA Langley Transonic Dynamics Tunnel is capable of simulating flight-representative Reynolds numbers of launch vehicles on the pad and is uniquely capable of replicating many fluid-structure scaling parameters typical of aeroelastic tests. Recent test efforts successfully developed representative atmospheric boundary-layers for three launch sites in the Transonic Dynamics Tunnel, thereby allowing all known aerodynamic and fluid-structure coupling parameters to be simultaneously simulated for those sites. Dynamic aeroelastically-scaled models representative of typical large launch vehicles were constructed for testing. Aeroelastic scaling includes matching geometry, mode shapes, reduced frequencies, damping, running mass ratios, and running stiffness ratios. The models were tested in smooth uniform flow and then immersed in the atmospheric boundary-layer for comparison of these effects. Dynamic data were acquired measuring unsteady pressure, acceleration, and base bending moment. It was discovered that peak dynamic loads resulting from resonant wind-induced oscillation response are similar when acquired in either smooth uniform flow or an atmospheric boundary-layer. This indicates that resonant lock-in events are minimally impacted by representative turbulence and/or wind profile. Alternately, nonresonant wind-induced oscillation response events are stronger when acquired in an atmospheric boundary-layer. This indicates that a lowly-damped structural response will increase when exposed to an increased magnitude of random excitation, which is consistent with historical comparisons. Loads created by the resonant response events were substantially stronger than those from the nonresonant response events. Therefore, if testing is done to simply identify worst-case conditions and load magnitude, then smooth uniform flow is likely an adequate test technique. However, if nonresonant response loads are of primary interest, then atmospheric boundary-layer simulation is required.

Thomas G Ivanco↗

Exploring the Feasibility of INCONEL® ALLOY 740H® for Power Plant Headers: Integrating Machine Learning with Computational Fluid Dynamics (CFD)

This keynote presentation explores the behavior of headers—essential components of pipeline systems—using ANSYS simulation software and machine learning techniques. The study aims to predict the thermal and mechanical performance of headers under diverse conditions through both steady-state and transient simulations. We investigate critical parameters such as heat transfer coefficient, fluid velocity, and temperature to optimize header design. Conducted as part of a DOE project led by NCAT in collaboration with UNC Charlotte, this research encompasses multiple key topics. The initial section focuses on the behavior of header systems under steady-state conditions using ANSYS simulation. It underscores the importance of headers in industrial infrastructure, especially in the energy sector, and examines the implications of material selection and flow direction on heat transfer dynamics. Methodologically, we employ Computational Fluid Dynamics (CFD) analysis through ANSYS, detailing the development of models, material properties, geometry specifications, boundary conditions, and meshing strategies. Our simulations explore various operational parameters, including temperature and mass flow rates, crucial for predicting heat transfer coefficients and enhancing header design. Results from the study include parametric investigations into mesh sensitivity, viscosity model evaluations, and the effects of heat transfer locations, all validated against theoretical calculations. We conclude with insights on mesh optimization, the suitability of viscosity models, and recommendations for future research aimed at improving header system efficiency and sustainability in industrial applications.

20 FOSSIL-FUELED POWER PLANTS↗

Unsteady Phenomena During Operation of the SSME Fuel Flowmeter

This report describes a part of the analysis carried in support of the SSME (Space Shuttle Main Engine) Fuel Flowmeter redesign, addressing an intensely researched phenomenon known as "shifting" of the flowmeter constant value. It consists of a sudden change in the flowmeter indication, which occurs simultaneously with the onset of an oscillatory variation of the rotor speed. The change in the flowmeter indications does not correspond to a real change in the volumetric flow through the device. Several causes have been investigated in detail, in the past, without conclusive evidence towards a cause of this phenomenon. The present analysis addresses the flow physics through the flowmeter by assembling results from 3-D CFD (computational fluid dynamics) calculations, airfoil C(sub D)/C(sub L) performance curves and mass moment of inertia characteristics of the rotor into a synergistic calculation which simulates the unsteady regime of the flowmeter operation. The results show that the 4-bladed rotor interacts with the periodic flow pattern created behind the flow straightener upstream in a manner that generates a steady, periodic fluctuation in the rotor's speed. The amplitude of this fluctuation is significantly smaller than the 0.5% of mean speed threshold which constitutes a flight operational limit. When manufacturing errors occur, however, the fluctuations are amplified and can generate a significant apparent change in the flowmeter indication. Two types of possible fabrication errors-which can occur even for parts fabricated within the accepted tolerances for the blade airfoil-are presented, together with their effect on the flowmeter operation.

Marcu, Bogdan↗

Investigation of Gas Turbine Internal Cooling Using Supercritical CO 2 —Effect of Surface Roughness and Channel Aspect Ratio

In this paper, an experimental and numerical investigation of internal cooling channels with rib turbulators is presented with sCO 2 as the working fluid at process conditions (pressure-20.7 MPa and temperature up to 150 °C). The effect of channel aspect ratio up to 2:1 on thermal-hydraulic performance is explored in additively manufactured rectangular channels and square channels, both with and without 60 deg ribs on the top and bottom sides. The Wilson-plot method is employed to experimentally measure channel-averaged Nusselt number over a Reynolds number range up to 370,000. The friction factor is calculated from pressure drop and mass flow rate and additionally, the overall thermal performance factor (TPF) is reported. A companion computational fluid dynamics (CFD) simulation is performed for the rib turbulated cooling configurations reported in the experiments using the Reynolds average Navier–Stokes-based turbulence model. The objective of the numerical study is to gain insight into the local heat transfer augmentation in the ribbed channels as a result of varying the aspect ratio, channel configuration (square versus rectangular), operating conditions (Reynolds number) and the surface roughness, an inherent outcome of the additive manufacturing process. Surface roughness is simulated using sand grain roughness height (K S ) calculated from the experimental data, and a comparison is presented with the corresponding channel configuration with varying surface roughness heights starting from smooth surfaces (K S = 0). Experimental results indicate that the heat transfer augmentation is negligible in the rectangular channels with ribs on the long side compared to the square channel. However, it is enhanced by 60% in comparison to placing ribs on the shorter side. The TPF remains constant at around 1 for the entire range of Reynolds numbers consistent with prior work at the National Energy Technology Laboratory (NETL). Here, the simulation results highlight that increased surface roughness can have a favorable considerable influence on Nusselt number and overall thermal performance enhancement.

42 ENGINEERING↗