Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Flow Conditions”

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 1,153 records · Page 64

Control methods for mitigating flow oscillations in a supercritical CO 2 recompression closed Brayton cycle

A dynamic model of a 10 MWe supercritical CO 2 (sCO 2 ) recompression closed Brayton power generation cycle is used to investigate control methods for mitigating oscillations in flow conditions and power demand load, especially as the load ramp rate is increased and the load setpoint is more closely tracked. The focus is on control of the main compressor inlet temperature (MCIT) since oscillations are prompted by the strong nonlinearities in sCO 2 fluid properties near the critical point and by changes in system inventory for efficient load management. In this study, the oscillatory behavior during load ramping is first shown and then two potential control solutions for substantially reducing oscillations are given: 1) the use of the main compressor inlet guide vanes (IGV) and 2) bypassing a portion of the sCO 2 around the cooler to maintain the MCIT. These control methods reduce the effect of cycle feedback from interactions between the MCIT and load. When implementing either control method, a load setpoint ramp of 7.5%/min was maintainable while achieving closer load setpoint tracking with limited oscillations. Finally, the impact of sensor noise on the control was examined since derivative action was used in the IGV and cooler sCO 2 bypass control. Finally, while some effect was evident, signal noise was not problematic to the control.

24 POWER TRANSMISSION AND DISTRIBUTION↗

On the prediction of critical heat flux using a physics-informed machine learning-aided framework

The critical heat flux (CHF) corresponding to the departure from nucleate boiling (DNB) crisis is essential to the design and safety of a two-phase flow boiling system. Despite the abundance of predictive tools available to the thermal engineering community, the path for an accurate, robust CHF model remains elusive due to lack of consensus on the DNB triggering mechanism. This work aims to apply a physics-informed machine learning (ML)-aided hybrid framework to achieve superior predictive capabilities. Such a hybrid approach takes advantage of existing understanding in the field of interest (i.e., domain knowledge) and uses ML to capture undiscovered information from the mismatch between the actual and domain knowledge-predicted target. A detailed case study is carried out with an extensive DNB-specific CHF database to demonstrate (1) the improved performance of the hybrid approach as compared to traditional domain knowledge-based models, and (2) the hybrid model’s superior generalization capabilities over standalone ML methods across a wide range of flow conditions. The hybrid framework could also readily extend its applicability domain and complexity on the fly, showing an elevated level of flexibility and robustness. Based on the case study conclusions, the window-type extrapolation mapping methodology is further proposed to better inform high-cost experimental work.

42 ENGINEERING↗

CFD modeling of turbulent air flow in self-heated gyroid TPMS structures: Thermal-hydraulic performance and validation

The application of mathematically derived geometries, such as triply periodic minimal surface (TPMS) lattices, has garnered significant interest across various fields, including the nuclear sector, due to their superior thermal-hydraulic characteristics for heat transfer compared to traditional plain or finned tubes. Here, this study validates a computational fluid dynamics (CFD) model, evaluates different turbulence models and CFD model settings, and performs uncertainty quantification to provide a comprehensive analysis. Despite extensive research on CFD modeling of TPMS lattices, such as gyroid and diamond geometries, there is a notable lack of publicly available literature providing comprehensive details on numerical analysis aspects, including convergence and methodological best practices. This study embarks on a benchmark analysis of a gyroid geometry to evaluate its thermal-hydraulic performance under turbulent flow conditions and scrutinize various CFD model configurations. The main contributions of this work include validating the CFD model, assessing and comparing different turbulence models, and enhancing pressure drop and temperature prediction capabilities. The results aim to support the development of methodologies needed to benchmark and enhance numerical analysis techniques for TPMS lattices. This work seeks to complement the existing body of knowledge, support the development of TPMS reactor concepts, and improve best practices for CFD modeling of TPMS lattices, ultimately advancing methodologies to support future applications in this domain.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Analysis of convection heat transfer on multiscale rough superhydrophobic and liquid infused surfaces

Multiscale rough superhydrophobic or slippery liquid infused porous surfaces have gained much interest in recent years for their improved transport phenomena properties. While there have been several studies on drag reduction and condensation on non-wetting surfaces, convection heat transfer that is important in many thermal and thermochemical applications has not been addressed systematically. Here, this article utilizes a fractal description of rough surface topographies to develop analytical models for the Nusselt number and the thermal hydraulic factor for fluid flow and heat transfer inside a cylinder with non-wetting surfaces. For air-infused superhydrophobic surfaces, the model considers the dynamic stability of the air/fluid interface in the asperities. Using the analytical formulations and the stability criteria, systematic studies are presented on the effects of the fractal surface parameters, cylinder radius and Reynolds number on the convective heat transfer characteristics, from which surface texture design maps are developed for maximizing the convection heat transfer. It is shown that multiscale non-wetting surfaces are most effective in the range of lower Reynolds number and small cylinder radius for achieving the best convective heat transfer and thermal hydraulic performance. Applying the models to actual non-wetting surface topographies fabricated using electrodeposition and chemical etching, it is shown that contrary to prevailing notion, superhydrophobicity, characterized by the highest contact angles, does not always lead to the maximum convective heat transfer performance, and that under certain fluid flow conditions, hydrophobic surfaces may offer a greater thermal performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Insights into mixing of non-isothermal multi-polymer melts for complex plastics recycling

Catalytic recycling or upcycling of plastics is often limited not by catalyst performance, but by transport, arising from highly viscous, non-Newtonian polymer melts. In this work, we develop a reactor-scale framework that integrates rheological measurements, constitutive modeling, computational fluid dynamics (CFD), and experiments to quantify mixing, heat transfer, and dispersion in surrogate hydrocarbon melts representing mixed plastics systems. Temperature- and shear rate-dependent viscosity of low-density polyethylene (LDPE) and high-density polyethylene (HDPE) is measured to create two surrogate polymers (PLD and PHD) that capture the dominant shear-thinning flow behavior while neglecting strong elastic effects, enabling tractable simulation of non-isothermal, polymer-melt mixing using a Carreau-Arrhenius generalized Newtonian framework. Three-dimensional CFD simulations are employed to evaluate impeller performance in PLD using mixing time, cavern volume, thermal uniformity, and interfacial area for regimes in which viscoelastic effects are not dominant. We show that magnetic stir bars commonly used in lab-scale studies produce large thermal gradients (~60 °C) and poor mixing, even under idealized power delivery and polymer flow conditions. In contrast, close-clearance anchor impellers achieve near-isothermal operation, reduce mixing times by up to 5×, and provide >90% active circulation volume. We further demonstrate that, at low pseudo-Deborah number (De*), motor power requirements can be predicted directly from shear rate-dependent rheology using the Carreau-Arrhenius framework, enabling rational selection of operating conditions. Extension to surrogate immiscible multi-polymer systems based on PLD and PHD shows that interfacial area is highly sensitive to operating conditions and impeller design, with coaxial anchor-turbine configurations enhancing dispersion by up to 4 × .

Close-clearance impellers↗

Rheological and pumping behaviour of 3D printable cementitious materials with varying aggregate content

3D printing with cementitious materials is an emerging technology and provides an active area of research. Despite this, there are only limited studies that examine the pumping behaviour of high yield stress fluids such as printable concrete. We present a study on the rheological behaviour of 3D printable cementitious materials with different aggregate-to-binder ratios (a/b) and analyse the pumping and printing characteristics with the measured rheological parameters. We observed that an increase in a/b resulted in a significant increase of plastic viscosity, and a nominal increase in yield stress and storage modulus. The variation in plastic viscosity and yield stress with increasing aggregate content was analysed with the Krieger-Dougherty and Chateau-Ovarlez-Trung model, respectively. Also, the rheological properties of the paste rich lubrication layer during pumping were determined with the help of a tribometer. The yield stress and viscous constant of the lubrication layer increase with the increase in a/b. The pumping pressure predictions were carried out by a suitable theoretical model by considering appropriate flow conditions for different mixtures. The variation in critical height of printed elements with an increase in a/b was correlated with the observed rheological parameters.

36 MATERIALS SCIENCE↗

Quantifying transport and electrocatalytic reaction processes in a gastight rotating cylinder electrode reactor via integration of Computational Fluid Dynamics modeling and experiments

Understanding the complexity of the multiple processes of mass, momentum, charge, and heat transport, and how these affect reaction kinetics at the electrode/electrolyte interface is one of the major challenges in the field of energy and catalysis. The rapid and rational scale-up of electrocatalytic systems to industrial scales require a detailed understanding of nonlinear transport-reaction processes, accessible only through the building of multi-physics models that capture with high fidelity the complexity of real-world devices. The gastight rotating cylinder electrode (RCE) reactor is a promising lab-scale tool that can decouple transport from intrinsic kinetics to generate data for first-principle models useful in the design of industrial, electrochemical reactors. Computational Fluid Dynamics (CFD) studies have previously been used to investigate the bulk flow in RCE reactors for simple corrosion and electroplating processes. However, the quantification of changes in local concentration within the viscous layer where catalysis takes place requires capturing the correct flow conditions inside the hydrodynamic boundary layer near the surface of the electrode. Further, this requires simulations with spatial resolution in the nm and μm scale and temporal resolutions between ms and s scales that are similar to the timescales for reactions on the electrode surface. In this study, experimental electrocatalysis is combined with CFD modeling to elucidate and parameterize the hydrodynamics in a gastight RCE reactor. CFD simulations of the electrochemical ferricyanide reduction reaction under mass transport limited conditions are used to evaluate the validity of the CFD model parameters by comparing calculated dimensionless mass transport descriptors to dimensionless correlations obtained experimentally. Justifications for assumptions and details of the simulation methods used in this study are presented to provide a detailed understanding of the effect that each model parameter has on the ability to accurately simulate electrocatalysis in RCE systems. The simulation methodology reported here is a first step towards the development of multi-scale models for the study of transport dependent electrocatalytic processes, such as the electrochemical transformation of CO 2 to fuels and chemicals.

42 ENGINEERING↗

Modeling and mitigating airborne pathogen risk factors in school buses

Computational fluid dynamics (CFD) models were developed to simulate the impact of different ventilation scenarios on airborne exposure risks in a 72-passenger school bus. Here, scenarios and factors that were investigated included a moving vs. stationary bus, impacts of a heating unit within the bus, and impacts of alternative ventilation scenarios with different combinations of openings (e.g., windows, door, emergency hatch). Results of the simulations showed that when the bus was stationary, use of the heater increased receptor concentrations unless there was another opening. When the bus was moving, simulations with at least two sets of openings separated from each other in the forward and aft directions produced a through-flow condition that reduced concentrations via dilution from outside air by a factor of ten or more. A single opening in a moving bus generally increased concentrations throughout the cabin due to increased mixing with minimal ventilation. The cumulative exposure risk (time-averaged concentrations) was found to be inversely correlated to the air exchange rate. Stationary and moving-bus scenarios that yielded above ~20 air changes per hour resulted in the lowest cumulative exposures. Recommendations from this study were implemented in new safety and operating procedures by the Albuquerque Public Schools Transportation Center.

60 APPLIED LIFE SCIENCES↗

A fast reduced model for a shell-and-tube based latent heat thermal energy storage heat exchanger and its application for cost optimal design by nonlinear programming

Numerical simulation of latent heat thermal energy storage (LHTES) systems plays a fundamental role in studying the physical process and guiding the engineering design. Discretization of the PDEs describing the nonlinear solidification/melting process of phase change materials (PCMs) leads to a large-scale complex dynamics system, where the system behavior depends on a set of parameters. In a design setting, repeated model evaluations are required over the set of parameters results in significant computational burden. In this paper, an explicit analytic solution was built for the propagation of the solidification front in a cylindrical coordinate. The analytic solution approach is further employed to develop a low computational reduced model (RM) as a module for a shell-and-tube based LHTES heat exchanger. The levelized Cost of Energy (LCOE) is used as a design metric and the RM model is used to apply system-level constraints in the nonlinear programming formulation that facilitates efficient global optimal design of the PCM properties, flow conditions and tube geometries. The use of LCOE as the design metric prevents over design of the heat transfer rate and also establishes a fair ground for evaluation of different thermal storage technologies and their integrated applications with other systems. Optimal results showed that a higher effectiveness results in a higher LCOE; the velocity of the HTF and the length of the channel are highly correlated with each other; both larger PCM latent energy and conductivity result in lower LCOE.

25 ENERGY STORAGE↗

Ligand-coordination effects on the selective hydrogenation of acetylene in single-site Pd-ligand supported catalysts

The selective hydrogenation of acetylene to ethylene is a critical step in the synthesis of polyethylenes. Achieving high conversion to ethylene without over-hydrogenation to ethane is a challenge that requires control of the transition metal site, which we achieve through a ligand-coordinated supported catalyst (LCSC) strategy. Using Pd catalysts coordinated to 1,10-phenanthroline-5,6-dione (PDO) ligands on CeO2 supports, we have discovered that the reaction selectivity depends strongly on the ligand:metal ratio with higher selectivity when more ligand is present in the catalyst. Catalyst structure was examined by extended X-ray absorption fine structure spectroscopy, transmission electron microscopy, and CO adsorption, which indicate single-atom character of the Pd. The ligand:metal ratio is determined by Xray photoelectron spectroscopy measurements and correlated with hydrogenation reactions under steady-state flow conditions to examine trends in hydrogenation activity and selectivity. Those trends can be better understood by density functional theory calculations that indicate hydrogen binding on the ligand to guide reaction selectivity toward the desired ethylene hydrogenation product. Here these results demonstrate the importance of considering the dynamic character of LCSCs and inform the design of future single-site heterogeneous catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fabrication and durability characterization of superhydrophobic and lubricant-infused surfaces

Hypothesis: Practical applications of non-wetting surfaces require good mechanical durability in the wet environments for which they are intended to be used. Durability of non-wetting surfaces is influenced by the surface features, interaction with the functionalization agent, and the lubricant properties that can be tuned independently to identify optimal combination. Experiments: In this study, superhydrophobic and lubricant-infused surfaces are fabricated on copper tubes using chemical etching and electrodeposition texturing techniques, six different functionalizing agents, and five different infused lubricants. Through 180 fabrication combinations and 102 durability tests, each parameter is systematically studied for contributions to initial non-wetting behavior and its durability in heated, wet environment, under high-energy water jet impingement, and under accelerated flow conditions. Findings: Among the adsorbing and curing functionalization agents investigated, n-Hexadecyl mercaptan that belongs to the sulfhydryl group and Sylgard-184, respectively, showed high durability in heated water immersion and under jet impingement tests. For lubricant-infused surfaces, lubricants with high surface tension demonstrated high durability in heated water immersion test, whereas durability in hydrodynamic conditions is closely correlated to lubricant viscosity. Results showed that a lubricant-infused surface will maintain its non-wetting properties in dropwise condensation conditions for approximately 1.5 years.

36 MATERIALS SCIENCE↗

Numerical simulation of vortex-induced vibration response of a single IEA 10-MW wind turbine blade

Three-dimensional simulation of vortex-induced vibration (VIV) of a single International Energy Agency (IEA) 10-MW reference wind turbine blade with a length of 97.325 m is performed using the ExaWind stack, an open-source suite of codes. This study aims to illustrate the spanwise VIV response characteristics and cross-validate the results with an existing commercial framework. Five near-body meshes and three time steps are selected for the convergence study. To improve computational efficiency, several VIV triggering methods are also compared to shorten the VIV development period. The ExaWind-based VIV simulation strategy for a single IEA 10-MW blade is determined. First, the modal shape is validated against published results. Then, spanwise VIV responses of four blade configurations under a fixed and varied incoming flow velocity are analyzed. Results show that the VIV response is dominated by the first edgewise (second overall) mode. Little first-mode contributions appear near the second-mode node, producing a pi phase jump, and a higher harmonics response occurs near the blade root. Rotational degrees of freedom are minor compared with translational motion. The response versus reduced velocity is analyzed, showing a two-branch behavior similar to that of VIV for a bluff cylinder. Across all tested cases, the dominant frequency remains locked to the natural frequency of the second mode with no observed desynchronization. A mild deviation is observed for the case of 90-degree pitch and 310-degree azimuth rotation near a reduced velocity of 6, which will be examined with additional cases in future work. These findings indicate that severe VIV responses can arise under specific configurations and flow conditions, thereby increasing the potential for VIV fatigue damage and requiring greater attention during operation.

17 WIND ENERGY↗

Wind power production from very large offshore wind farms

In this work, we provide the first quantitative assessment of power production and wake generation from offshore wind energy lease areas along the U.S. east coast. Deploying 15-MW wind turbines, with spacing equal to the European average, yields electricity production of 116 TWh/year or 3% of current national supply. However, power production is reduced by one-third due to wakes caused by upwind wind turbines and wind farms. Under some flow conditions whole wind-farm wakes can extend up to 90 km downwind of the largest lease areas, and the frequency-weighted average area with a 5% velocity deficit is 2.6 times the footprint of the lease areas. Simulations including maritime corridors demonstrate reduction in the wake effects leading to power-efficiency gains and may offer contingent benefits. First-order scaling rules are developed that describe how “wake shadows” from large offshore wind farms scale with prevailing meteorology and wind turbine installed densities.

15 MW↗

Ethanol-fueled metal supported solid oxide fuel cells with a high entropy alloy internal reforming catalyst

High-performance metal supported solid oxide fuel cells (MS-SOFC) with an integrated high entropy alloy (HEA) internal reforming catalyst (IRC) are demonstrated for transportation applications using ethanol and methanol as fuels. Addition of the HEA IRC dramatically improves cell performance and stability when using ethanol/water blend fuel. Absence of carbon deposition predicted by thermodynamic calculations is confirmed by Raman spectroscopy analysis of posttest anodes. Optimal catalyst processing (deposition technique, loading, firing temperature) and cell operation conditions (flow rates, temperature, fuel compositions) are explored. Infiltrated HEA reforming catalyst provides a highly porous structure and low catalyst loading (6 mg cm –2 ). The designed structure and catalysts achieve small mass transport resistances in the fuel electrode (26.2 s m –1 ) and oxygen electrode (41.6 s m –1 ). The best ethanol concentration (60:40 v% ethanol: water) provides 0.83 W cm –1 at 700 °C, without carbon deposition. The ethanol-fueled MS-SOFC is operated for 500 h, including five thermal cycles. As a result, cell evolution is similar to that reported previously for hydrogen fuel; nickel aggregation and chromia deposition were the major observed changes, and carbon formation can be avoided even after long-term operation.

30 DIRECT ENERGY CONVERSION↗

Enriching OpenStreetMap network data for transportation applications: Insights into the impact of urban congestion on accessibility

OpenStreetMap (OSM) data is a valuable open-source resource for various transportation, traffic, and planning applications. However, OSM network data lack operating traffic speed information, which is critical for transport planning and operations. Addressing this shortcoming, this study leverages commercial vendor data (to serve as ground truth) with exogenous, open-source variables characterizing local transport infrastructure, land use, and demographic information to predict average congested traffic speeds on OSM networks. Three machine-learning models were tested and estimated for OSM links with and without speed limit information in the Denver metropolitan region. Among these, XGBoost performed best, with mean absolute errors of 3.27 and 3.62 mph for links with and without speed limits, respectively. The developed models accurately predicted traffic speeds for different hours and days of the week compared to ground truth data. Using these predicted speeds, drive accessibility scores were computed for the Denver region for different time periods using the Mobility Energy Productivity (MEP) metric to understand the impact of congestion on energy-efficient accessibility. Results show that congestion-adjusted drive accessibility can be significantly lower compared to accessibility calculated using free flow speeds. Specifically, weekday evening hours saw a 42 % drop in accessibility due to reduced speeds, particularly around downtown Denver. Across the Denver metro region, approximately half as many opportunities and jobs are accessible in under 20 min by car during the evening peak period relative to free flow conditions. These findings underscore the importance of using congestion-adjusted operating speeds rather than speed limits in accessibility calculations, as reliance on speed limits can substantially overestimate energy-efficient drive accessibility in large, car-centric cities susceptible to significant congestion. In conclusion, the methodology presented here could further enrich OSM network data, making them useful for an even broader range of transportation applications.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Cost optimization of high recovery single stage gap membrane distillation

Estimating the cost of membrane distillation (MD) processes at scale is confounded by the vast array of possible process configurations, module configurations, and operating conditions and the extreme sensitivity of MD costs to these parameters. This work develops a comprehensive MD cost optimization framework for minimizing the levelized cost of water (LCOW) as a function of feed salinity and temperature. Our single stage gap MD model fully captures solution properties, heat transfer, and mass transfer relationships in the modules to determine the cost optimal design as a function of multiple decision variables, including process configuration (feed location and connectivity between the hot and cold-side), the gap type (air, permeate, conductive), equipment sizes (membrane, heater, chiller, and heat exchanger units), and operating conditions (flow rates, hot and cold side temperatures, and heating and cooling duties). We estimate that the minimum LCOW of conductive gap MD ranges between $\$$10 and $\$$16 per m 3 for feed salinities of 25 g/L to 200 g/L and water recoveries of 30%–75%. Further, we update several heuristic assumptions for cost optimal system operation, including: increased thermal conductivity of the gap decreases system costs; cost optimal MD operates at the highest allowable temperature; and the LCOW is a strong function of brine salinity. Finally, we find that improving membrane properties and lowering heating and cooling costs will provide only modest reductions in costs (<25%).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Finite element modeling of electropolishing of radio frequency dipole Nb crab cavity in hydrofluoric-sulfuric acid electrolyte

The superior performance of niobium superconducting radio frequency (SRF) cavities is achieved through meticulous surface treatment, notably via chemical electropolishing, ensuring exceptionally smooth surfaces. While this technique has been extensively employed for cylindrically symmetric structures, addressing more intricate geometries poses significant challenges in achieving uniform polishing and controlled material removal, especially when moving away from retractable cathodes. This presents a multifaceted electrochemical, thermal, and fluid dynamics conundrum. A prime example is the 197 MHz radio frequency dipole (RFD) crabbing cavity proposed for the Electron Ion Collider (EIC) project, exemplifying such complex structures. Our groundbreaking work integrates the localized oxide thickness variation, considering its impact on current distribution and Joule heating, within the framework of multi-physics modeling using the COMSOL® simulation suite. This approach was applied to a comprehensive model of the RFD cavity, allowing us to investigate optimal external cooling water flow conditions necessary for achieving desirable outcomes. In conclusion, this illustrates the increasing utility of such multi-physics codes to greatly aid in designing solutions to challenging processing requirements for increasingly complex accelerator cavities.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Improved departure from nucleate boiling prediction in rod bundles using a physics-informed machine learning-aided framework

The critical heat flux (CHF) corresponding to the departure from nucleate boiling (DNB) crisis is a regulatory limit for the licensing of pressurized water reactors (PWRs) worldwide. Despite the abundance of predictive tools available to the reactor thermal-hydraulics community, the path for an accurate CHF model remains elusive. This work approaches the prediction of DNB through a physics-informed machine learning-aided framework (PIMLAF) with the objective of achieving superior predictive capabilities for a rod bundle. In view of the limitations in existing macro-scale physics-driven tools, an improved mechanistic model is first proposed, leveraging key concepts in the liquid sublayer dryout and bubble crowding mechanisms. Furthermore, the proposed mechanistic model is able to predict DNB in different heater geometries for a broad range of flow conditions without the need for recalibration. This model is then incorporated as the physics-informed component of the hybrid framework PIMLAF, which takes advantage of established understanding in the field (i.e., domain knowledge [DK]) and uses machine learning (ML) to capture undiscovered information from the mismatch between the actual and DK-predicted output. Two bundle-related case studies using the PWR subchannel and bundle tests (PSBT) database are carried out to illustrate the PIMLAF’s improved performance over traditional approaches for both interpolation and extrapolation purposes. In light of the PIMLAF’s promising potential to reduce prediction error, reactor vendors are encouraged to leverage their in-house experimental efforts and apply the hybrid framework to potentially achieve margin reductions in the minimum DNB ratio (MDNBR) for the designs of interest.

42 ENGINEERING↗