Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Heat transfer performance”

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 541 records · Page 30

Operating Hardware Impact on the Heat Transfer Properties of Windows

Despite significant advancements in fenestration technology in the last two decades, the thermal transmittance of fenestration products is still significantly higher than that of walls. This corresponds to 60% of the total energy loss of a modern building envelope through the windows. Hence, further development and improvements of fenestration products are necessary. Increasingly stringent codes and standards for fenestration stimulate industry to work on improved solutions. Thus, it is crucial that assessment techniques are able to account for innovations accurately. The thermal effects of non-continuous hardware in window frames are currently ignored by international rating procedures. A preliminary investigation conducted by our team showed significant performance degradation in two of the three out-opening casement profiles caused by the presence of operating hardware. Frames with the structure made of vinyl and fiberglass consist of many air cavities that are penetrated by operating hardware made of highly conductive materials. In these frames, in order to have an accurate assessment, it may be required to employ three-dimensional modeling due to the convective nature of heat transfer within the cavities. However, in this study, we demonstrate that the three-dimensional (3D) effects of non-continuous hardware can be approximated accurately with simpler two-dimensional (2D) simulations. We then develop a simplified model based on weighted average capable of replacing the time- and computation-intensive 3D simulations with 2D simulations and validate it against market available frames and their corresponding hardware. Validation results show that our approximation technique results in discrepancies lower than 0.05 W/(m2K), or 3% of the total thermal transmittance. Thus, we conclude that simplified 2D simulation models may be used for predicting hardware impact in window frames with reasonable accuracy. As windows and glazing structures are becoming ever better thermally insulated, it is becoming even more important to be able to model the impact of the operating hardware on the total thermal performance in order to design the best windows possible and not let the operating hardware ruin an otherwise well-proven design, which is hence addressed in this study.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Performance of Microreactor Test Article with Embedded Sensors During Testing in The Single Primary Heat Extraction and Removal Emulator

The nuclear industry is pursuing microreactors that can be factory assembled and deployed to remote regions for reliable power generation. One class of microreactors uses a monolithic metal core block coupled to heat pipes that use passive heat flow, increasing the surface area for heat transfer without requiring active coolant flow through the reactor core. Additional experimental testing is required to understand the heat rejection limitations of heat pipes and thermal stresses in the monolithic core block due to significant temperature gradients. This report describes the initial characterization and testing of a stainless-steel test article that was fabricated with embedded sensors to measure heat pipe performance limits, as well as spatially distributed temperatures and strains during electrically heated thermal testing to simulate nuclear heating. The electrically heated testing was performed in the Single Primary Heat Extraction and Removal Emulator facility located at Idaho National Laboratory. The ultimate goals of this work are to (1) accurately monitor temperature and strain distributions that result from differential thermal expansion in the test articles and (2) quantify the heat rejection limits of heat pipes as a function of operating temperature and working fluid during steady-state and transient operations. Initial tests focused on the feasibility of using advanced fiber-optic sensors and other sensor technologies to improve the understanding of temperature and strain distributions within the electrically heated experiments. However, these sensing capabilities could benefit the broader microreactor community if the sensors could be used to monitor component and system health during nuclear operations to inform a limited number of microreactor operators, ultimately reducing operation and maintenance costs and moving toward semiautonomous operation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Toward extreme high-temperature supercritical CO 2 power cycles: Leakage characterization of ceramic 3D-printed heat exchangers

We report future supercritical carbon dioxide (sCO 2 ) Brayton power cycles demand high-performance gas-to-gas heat exchangers (HXs) operating under extreme temperature and pressure conditions at which most existing superalloy materials fail to function safely. Ceramic HXs are deemed excellent candidates for advanced sCO 2 power plants as they can withstand high-temperature working environments. Particularly, ceramic 3D printing enables compact HX topologies employing complex and efficient heat transfer features. However, ceramic 3D-printed walls separating hot and cold flow streams are susceptible to a through-plane leakage inherent to a powder-based manufacturing process, including ceramic 3D printing. A potential leakage through ceramic separating walls poses a significant challenge in developing reliable ceramic 3D-printed HXs and could deteriorate thermal performance. In this study, various parameters, including feedstock slurry, 3D-printing direction, and post-processing conditions, are considered, for the first time, to characterize the argon gas leakage rate associated with alumina 3D-printed parts. Three 3D-printed ceramic structures of flat plates, curved tubes, and small-scale plate-and-frame HXs with various thicknesses are systematically studied to determine powder and ceramic 3D-printing conditions to eliminate the through-plane leakage. The results showed that an alumina 3D-printed plate with a thickness of 0.75 mm demonstrates a permeability of 6 x 10 -4 milli-darcy. An alumina 3D-printed tube with a wall thickness of 0.9 mm revealed a permeability of 9.6 x 10 -7 milli-darcy. Furthermore, leakage test results of functional 3D-printed modules showed a dependency on the 3D-printing direction. Particularly, alumina cell-scale HXs employing 1.5-mm-thick horizontal and vertical 3D-printed separating walls demonstrated impermeability and gas permeability of 7.2 x 10 -5 millidarcy, respectively. Insights gained from the present study facilitate the development of complex ceramic 3D-printed HXs and other balance of plant components for next-generation high-temperature high-pressure sCO 2 power cycles.

36 MATERIALS SCIENCE↗

Nonintrusive manufactured solutions for ablation

The study of heat transfer and ablation plays an important role in many problems of scientific and engineering interest. As with the computational simulation of any physical phenomenon, the first step toward establishing credibility in ablation simulations involves code verification. Code verification is typically performed using exact and manufactured solutions. However, manufactured solutions generally require the invasive introduction of an artificial forcing term within the source code such that the code solves a modified problem for which the solution is known. In this paper, we present a nonintrusive method for manufacturing solutions for a non-decomposing ablation code, which does not require the addition of a source term.

42 ENGINEERING↗

Oxygen enrichment combustion to reduce fossil energy consumption and emissions in hot rolling steel production

The reheating furnace operation in the hot mill is natural gas- and electricity-intensive. Oxygen enrichment combustion for reheating furnaces has been proposed to curb and replace natural gas use. In this study, heat transfer in steel slabs in the combustion environment of a push-type reheating furnace was simulated using a computational fluid dynamics (CFD) model. Two oxygen enrichment methods that optimized for performance were selected— a medium oxygen enrichment (MOE) case and an oxy-fuel (OF) case. A life cycle analysis (LCA) characterized the energy and emission profiles of an integrated iron and steel manufacturing process using the two oxygen enrichment cases for the hot mill. These conditions were evaluated for energy use and carbon intensity and compared with a baseline case. Results show that with oxygen enrichment, natural gas consumption can decrease by 19.6%–26.8%, total energy consumption (natural gas and electricity) can decrease by 15.1%–20.7% in the hot mill. Emissions of greenhouse gases can decrease by 11.1%–15.2% in the two optimized cases with 14%–27% reductions in regulated criteria pollutants (nitrogen oxides, carbon monoxide, particulate matter, volatile organic compounds, black carbon, organic carbon, and volatile organic carbons). There is a tradeoff between reducing natural gas consumption and increasing electricity demand from a life cycle perspective. Although the OF case resulted in higher energy- and emissions-related benefits, the MOE case showed the more desirable heat flux uniformity, which is key to maintaining product quality. The analysis suggests that oxygen enrichment in the reheating furnace process can have a significant impact on hot mill environmental performance and become a contributing factor in transitioning to low-carbon steel manufacturing.

36 MATERIALS SCIENCE↗

In situ embedment of type K sheathed thermocouples with directed energy deposition

Advanced nuclear reactor systems require new technologies for heat transfer and system monitoring. Additive manufacturing (AM) offers the design flexibility to allow in-situ sensor embedment through smart manufacturing for real-time monitoring, and performance of these systems. Here, this study focuses on experiments investigating the feasibility of in-situ sensor embedment using directed energy deposition (DED). Type K thermocouples are embedded into 316L stainless steel (SS) samples using two different configurations (e.g., exposed and embedded tips) and two designs (e.g., flush to substrate) within an AM base. Embedded sensor samples are analyzed via in-situ measurements and high-temperature performance validation tests at 350ºC and 900ºC. Temperature performance results at both temperature tests show good agreement with manufacturer specifications proving that these sensors could still capture accurate temperature readings after in-situ embedment during DED processing. An additional optimization experiment was conducted on the exposed tip configuration using a surrogate thermocouple to improve tolerances and the embedment process. Results improved tolerances, lower porosity, smaller gaps between the sensor and base, and better junction contact for the sensor. Although further optimization of this embedment strategy is necessary to improve the structural stability and tolerances within the component, this research strategy provides a proof-of-feasibility for DED embedment with commercial sheathed thermocouples. This research provides early impact on embedment of sensor for multiple materials and complex geometric components..

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Analytical investigation of the dew point water recovery system with potential application for drinking water production

The article presents the analytical investigation of the Dew Point Water Recovery (DPWR) system with potential application for drinking water production. Four different configurations of the DPWR system were proposed, differing in the number and type of dew point cooling exchangers used, and their operation was described using a validated convective heat transfer coefficient mathematical model. It was found that the three stage configu-ration of the DPWR system with a regenerative exchanger has the best water recovery performance, while the three stage configuration with a counter-flow exchanger has the best ratio of energy supplied to the system in relation to the amount of water produced. It was also found that the DPWR system can operate in a wide range of operating conditions with a specific electric energy consumption reaching below 0.9 kWh/m 3 and GOR in the range of 0.36-2.16. What distinguishes the DPWR system from other methods of drinking water production is the ability to work without an additional source of thermal energy apart from the energy of the air supplied to the system. Thanks to this property, the system has the potential to fill the technological gap of a low-energy and easily accessible method of drinking water production.

42 ENGINEERING↗

Temperature-dependent Battery Performance of a Na 3 V 2 (PO 4 ) 2 F 3 @MWCNT Cathode and In-situ Heat Generation on Cycling

Excellent structural stability, high operating voltage, and high capacity have made Na 3 V 2 (PO 4 ) 2 F 3 a promising cathode material for sodium-ion batteries. However, high-temperature battery performances and heat generation measurements have not been systematically reported yet. Carbon-coated Na 3 V 2 (PO 4 ) 2 F 3 @MWCNT (multi-walled carbon nanotube) samples are fabricated by a hydrothermal-assisted sol-gel method and the electrochemical performances are evaluated at three different temperatures (25, 45, and 55 °C). Furthermore, the well-crystallized Na 3 V 2 (PO 4 ) 2 F 3 @MWCNT samples exhibit good cycling stability at both low and high temperatures; they deliver an initial discharge capacity of 120–125 mAhg -1 at a 1 C rate with a retention of 53 % capacity after 1,400 cycles with 99 % columbic efficiency. The half-cell delivers a capacity of 100 mAhg -1 even at a high rate of 10 C at room temperature. Furthermore, the Na 3 V 2 (PO 4 ) 2 F 3 @MWCNT samples show good long-term durability; the capacity loss is an average of 0.05 % per cycle at a 1 C rate at 55 °C. Additionally, ionic diffusivity and charge transfer resistance are evaluated as functions of state of charge, and they explain the high electrochemical performance of the Na 3 V 2 (PO 4 ) 2 F 3 @MWCNT samples. In-situ heat generation measurements reveal reversible results upon cycling owing to the high structural stability of the material. Excellent electrochemical performances are also demonstrated in the full-cell configuration with hard carbon as well as antimony Sb/C anodes.

25 ENERGY STORAGE↗

Development and Validation of a Two-Phase Thermal-Hydraulic CFD Code NEK-2P

A project is underway to develop, verify and validate an advanced two-phase flow modeling capability for the highly-scalable, high-performance Computational Fluid Dynamics (CFD) code NEK5000. The goal of this work is to verify and validate the two-phase version of the NEK5000 code, named NEK-2P, to simulate the two-phase flow and heat transfer phenomena that occur in a Boiling Water Reactor (BWR) fuel bundle under various operating conditions. The NEK-2P two-phase flow models follow the approach used for the Extended Boiling Framework (EBF) previously developed at Argonne but include more fundamental physical models of boiling phenomena and advanced numerical algorithms for improved computational accuracy, robustness, and computational speed. The development of the NEK-2P two-phase solver and the implementation of the Extended Boiling Framework two-phase models were initially supported by Argonne National Laboratory (Argonne) through a Laboratory Directed Research and Development (LDRD) project during FY2014-2016. The development and validation of the two-phase models through analyses of selected two-phase boiling flow experiments was supported by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program in FY2017-2020. This report focuses on verification and validation of the water-steam boiling model NEK-2P Two-Phase, CFD code. The NEK-2P was validated with Nuclear Power Engineering Corporation (NUPEC) Pressurized Water Reactor (PWR) Sub-channel and Bundle Test (PSBT) void distribution benchmark. Three different simulations were performed and analyzed for various operating conditions such as wall-heat flux and sub-cooled inlet temperatures. Reasonably good agreement with measured data was obtained in predicting the measured void distributions. Simulations were performed for Virginia Tech. (VT) 3x3 rod bundle geometry with and without spacers. The preliminary results were presented for Simplified Spacer Grid (SSG). In addition, the implementation of interface reconstruction model was tested with one of the Becker benchmark Critical Heat Flux (CHF) experiments.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A Novel Concentrating Solar Weathering Apparatus for Experimental Validation of Multi-Modal Degradation Models

High-performance coatings for Concentrating Solar Power (CSP) receivers are subjected to remarkable environmental stressors during normal operations. Applied to the receiver tubes, these coatings serve to maximize the solar absorptivity of the receiver, transferring as much heat as possible from the solar collectors into the heat-transfer fluid (HTF). The lifecycle of these coatings is not well-defined, and the harsh operational conditions make them difficult to test. NREL has designed, built, and tested an apparatus to expose these samples to design levels of environmental stress and well beyond, into accelerated and destructive conditions. The chamber is actively cooled, monitored, and has the capability to supply humidification for cycling tests, allowing us to test multi-modal degradation and failure conditions at high temperature, high flux, and high humidity conditions. These conditions can catalyze high-temperature oxidation, mechanical degradation, and other modes of absorptivity loss seen in selective solar receiver coatings. The experimental data can feed lifecycle models for expensive and necessarily resilient materials, offering insights to aid maintenance schedules, technoeconomic analysis, and material industry performance benchmarks. This presentation will demonstrate the apparatus design and performance, as well as initial results for aging on a selective receiver coating.

14 SOLAR ENERGY↗

Design of a Combined Proportional Integral Derivative Controller to Regulate the Temperature Inside a High-Temperature Tubular Solar Reactor

Solar fuels are proven to be promising candidates for thermochemical energy storage. However, the transient nature of solar radiation is an obstacle to maintaining a stable operational temperature inside a solar reactor. To overcome this challenge, the temperature of a solar reactor can be regulated by controlling the incoming solar radiation or the feedstock flowrate inside the reactor. In this work, a combined proportional integral derivative (PID) controller is implemented to regulate the temperature inside a high-temperature tubular solar reactor with counter-current flowing gas/particles. The control model incorporates two control systems to regulate incoming solar radiation and gas flow simultaneously. The design of the controller is based on a reduced-order numerical model of a high-temperature tubular solar reactor that is vertically oriented with an upward gas flow and downward particle flow. The reactor receives heat circumferentially through its wall over a finite segment of its length. Formulation of the heat transfer model is presented by applying the energy balance for the reactor tube and considering heat and mass transfer inside. A set of governing differential equations are solved numerically by using the finite volume method to obtain reactor wall, particles, and gas temperatures along the reactor length with various boundary conditions. Simulation results are used to tune the PID controller parameters by utilizing the Ziegler–Nichols tuning method. Both the simulation results and the controller performance are visualized on the LABVIEW platform. The controller is challenged to track different temperature setpoints with different scenarios of transient solar radiation. Here, the performance of the PID controller was compared to experimental results obtained from an industrial PID controller embedded in a 7 kW electric furnace. Results show that the combined PID controller is successful in maintaining a stable temperature inside the reactor by regulating the incoming solar radiation and the flowrate via small steady-state error and reasonable settling time and overshoot.

14 SOLAR ENERGY↗

Modeling Molten Salt Spreading and Heat Transfer using MELTSPREAD – An Uncertainty Analysis

This report summarizes spreading and heat transfer results for molten salt pouring onto a flat stainless steel substrate that were calculated using MELTSPREAD to gain insight into bulk salt behavior during a molten salt spill accident. MELTSPREAD was developed at Argonne to model the one-dimensional flow and freezing behavior of molten corium and is being applied to model molten salts for the first time as part of this work. An uncertainty analysis of thermophysical properties and spill conditions was performed using MELTSPREAD to determine which properties and conditions have the greatest impact on model outcome. The model was run with and without the inclusion of decay heat. Eutectic FLiNaK was used as the salt composition because it is well-characterized and appropriate for model development. Small salt volumes were used for the initial model runs described in this report so that the model results can be compared to the results from ongoing experiments on molten salt spreading and heat transfer currently being conducted at Argonne at a benchtop scale. The spreading of a small volume of FLiNaK (50 mL) was found to be limited by the balance between the molten salt surface tension and gravity and not by freezing for the pour conditions and assumed mechanisms of heat transfer in the model. Heat transfer from the salt was highly inefficient to due to its low thermal conductivity, high heat capacity, and high heat of fusion and varying these properties had no effect on the spreading behavior. Changing the viscosity and density of the salt slightly affected the spreading behavior. Setting the salt surface tension to approximately zero to simulate the salt wetting the substrate had the greatest effect on the spreading behavior and this increased both the spreading velocity and overall spreading area. Including decay heat in the model had a significant effect on the long-term cooling behavior of the salt but had no effect on the overall spreading behavior. Future effort should involve modeling a larger volume of molten salt to represent a reactor case spill scenario because the spreading behavior of large volumes of salt may differ from the small volumes that were the focus of this document.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Development of CTF modeling of interfacial drag, wall shear, and interfacial heat transfer for bubbly and annular-mist flow regimes

This milestone report discusses improvements to CTF closure models for modeling of bubbly and annular mist two-phase flow conditions, which are important for accurate modeling and simulation for boiling water reactor (BWR) conditions. Models that were improved include the two-phase wall shear, wall boiling, and phase mass and heat transfer. A new flow regime map was implemented to improve on the prediction of transition to annular flow, which was generally largely over-predicted by the legacy CTF flow regime map. The annular mist interfacial drag models were also reviewed and new models were proposed for implementation into the code. In addition to implementing new models, an activity was performed to calibrate these new models to experimental data using a Markov Chain Monte Carlo technique, thus further improving their accuracy. Finally, new experiments were implemented into the CTF validation matrix which add more two-phase pressure drop and void measurements for the assessment of CTF for modeling of BWRs. The results show that the Chisholm model, which was implemented and then calibrated, results in significantly improved comparison against the Riso two-phase pressure drop data. Furthermore, improvements in the implicit heat transfer and wall boiling models resulted in good agreement with the void measurement data.

42 ENGINEERING↗

Nickel-based superalloy single crystals fabricated via electron beam melting

Additive manufacturing technologies have emerged as potential disruptive processes whose possible impacts range from supply chain logistics, prototyping, and novel materials synthesis. Numerous works illustrate the ability to control mi- crostructure in fusion based processes and a few recent authors have even produced single crystals. However, a number of open questions remain regarding the process window which enables printing of single crystals. Furthermore, it has been observed that these additively manufactured single crystals exhibit a preferred $\langle$011$\rangle$ secondary orientation. In this work we investigate fabrication conditions which enable printing of single crystals via electron beam melting. A space filling design of experiments is utilized to efficiently explore the fabrication space. Single crystals were success- fully obtained using both commercially available powders and custom melt alloys. Microstructures obtained via these exploratory experiments exhibited a continuum of columnar structures ranging from weakly textured polycrystals, near single crystal, and fully single crystalline material. Complex geometry experiments are performed to study the grain selection mechanism. We find that the grain selection mechanism is independent of the bulk scale geometry and must therefore by driven by local heat transfer and solidification dynamics. Furthermore, grain selection is shown to be driven by competing driving forces; one which prefers epitaxial growth and another which is driven by the imposed scan pattern geometry. Finally, a mechanism is proposed for the anomalous secondary orientation preference.

36 MATERIALS SCIENCE↗

Multi-stage load partitioning in additively manufactured Al6061+TiC nanocomposite characterized by in-situ neutron diffraction

Metal matrix composites (MMCs), combining metal matrix and ceramic particles, exhibit high mechanical strength and stiffness compared to conventional metallic materials. During fusion-based additive manufacturing (AM), MMCs undergo a rapid heating and cooling thermal history, which affects the interfacial bonding, thermal misfit stress and mechanical load transfer behavior between the metal matrix and particles, thus impacting the bulk mechanical performance. Here, we investigate the deformation dynamics and phase-specific load transfer in fusion-based AMed Al6061 +TiC MMCs through in-situ neutron diffraction. By calculating the phase-specific lattice strains, a multi-stage load transfer and deformation behavior during uniaxial compression is revealed: (1) elastic deformation in Stage I for both phases; (2) sudden stress rebalance between two phases in Stage II, evidenced by a sudden decrease in Al lattice strain and increase in TiC; (3) active load carrying by both phases in Stage III, with both phases experiencing an increase in lattice strains. Further, the deformation mechanism of each stage is deducted by correlating the evolutions of the interphase stresses and peaks’ broadening and intensity of the Al matrix. The local plastic deformation of the Al matrix near the phase interface is triggered and leads to stress rebalancing in Stage II. The global plastic deformation subsequently propagates throughout the Al matrix in Stage III, and the composite is further hardened through both dislocation multiplication and load transfer. The findings offer valuable insights into deformation and load-sharing behavior in AMed MMCs.

36 MATERIALS SCIENCE↗

Thermal performance of phase change materials with anisotropic carbon fiber inserts

In thermal energy storage systems, phase change materials (PCMs) are widely used for thermal energy management. Most PCMs have low thermal conductivities, which limits the heat transfer rate within PCM and thus makes the phase-changing process very slow. However, thermal conductivities of PCMs can be altered by inserting targeted additives. We hypothesize that the size and shape of these additive inserts play a key role in thermal management efficiency. To this end, the impacts of carbon fiber (CF) inserts on the phase change behavior and consequent heat transfer efficiencies of inorganic and organic PCMs were investigated using experiments and simulations. Long, anisotropic CFs with high thermal conductivities formed continuous fast heat flux tunnels inside PCMs to enhance the heat transfer. Such CFs could extend the phase change fronts from the limited container-shaped interface to the larger surface of numerous CF inserts inside the PCM. These special CF inserts work with a new heat transfer mechanism, different from conventional small additives or long isotropic CF inserts. The thermal energy release rate increased by 2.5 times with 1 wt.% anisotropic CF inserts in inorganic PCM. However, CF inserts in liquid organic PCM hindered the natural convection and compromised the improved heat conduction. In conclusion, the lab-scale multiphysics simulations support these experimental observations and indicate that CF inserts have potential to enhance heat transfer in inorganic PCMs, but they are less effective in organic PCMs.

25 ENERGY STORAGE↗

Understanding flow dynamics in membrane distillation: Effects of reactor design on polarization

Optimization and design of full-scale membrane distillation (MD) systems usually require Sherwood and Nusselt correlations that are developed from lab-scale systems. However, entrance effects in lab-scale systems can significantly impact heat, mass and momentum transfer in the reactor, therefore affect the accuracy of the developed experimental Sherwood and Nusselt correlations. Here, Computational Fluid Dynamics (CFD) simulations using OpenFOAM are performed to understand the effects of right-angled bends and inlet design on flow dynamics, temperature and concentration polarization in MD systems. Simulation results show that the presence of right-angled bends and inlets with sudden expansions lead to the formation of Dean vortices. Dean vortices enhance perpendicular mixing in MD systems and reduce both temperature and concentration polarization. Temperature and concentration polarization coefficients in MD systems with right-angled bends and inlets with sudden expansions vary significantly for the same volumetric flow rate. Our studies show that lab-scale systems with the same volumetric flow rate but different designs lead to significantly different Nusselt and Sherwood correlations. Here, this study demonstrates the importance of CFD-informed design of lab-scale systems to minimize entrance effects and suppress Dean vortices for consistent model development and calibration across multiple scales.

42 ENGINEERING↗

The discrete Green's function paradigm for two-way coupled Euler–Lagrange simulation

We outline a methodology for the simulation of two-way coupled particle-laden flows. The drag force that couples fluid and particle momentum depends on the undisturbed fluid velocity at the particle location, and this latter quantity requires modelling. We demonstrate that the undisturbed fluid velocity, in the low particle Reynolds number limit, can be related exactly to the discrete Green's function of the discrete Stokes equations. In addition to hydrodynamics, the method can be extended to other physics present in particle-laden flows such as heat transfer and electromagnetism. The discrete Green's functions for the Navier–Stokes equations are obtained at low particle Reynolds number in a two-plane channel geometry. We perform verification at different Reynolds numbers for a particle settling under gravity parallel to a plane wall, for different wall-normal separations. Compared with other point-particle schemes, the Stokesian discrete Green's function approach is the most robust at low particle Reynolds number, accurate at all wall-normal separations. To account for degradation in accuracy away from the wall at finite Reynolds number, we extend the present methodology to an Oseen-like discrete Green's function. The extended discrete Green's function method is found to be accurate within 6% at all wall-normal separations for particle Reynolds numbers up to 24. Furthermore, the discrete Green's function approach is well suited to dilute systems with significant mass loading and this is highlighted by comparison against other Euler–Lagrange as well as particle-resolved simulations of gas–solid turbulent channel flow. Strong particle–turbulence coupling is observed in the form of turbulence modification and turbophoresis suppression, and these observations are placed in context of the different methods.

42 ENGINEERING↗