Volume-averaged Model Research Articles

Experimental and theoretical studies on 3D printed short and continuous carbon fiber hybrid reinforced composites

3D-printed carbon fiber composites hold significant potential in aerospace applications because of their lightweight, high strength and complex structure fabrication capabilities. However, additive manufacturing characteristics such as high porosity, anisotropy and continuous fiber content significantly affect the mechanical properties of printed parts. The aim of this study is to investigate the influence of hybrid printing of short and continuous carbon fibers on the mechanical properties and failure mechanisms of composites and to develop a model to predict elastic properties considering porosity. First, mechanical testing and characterization of short and continuous fiber reinforced polyamide (nylon) print filaments are performed. The results indicate that the elastic modulus and ultimate strength of continuous carbon fiber filaments reach 60 GPa and 1500 MPa, respectively, while the strength and elastic modulus of short carbon fiber filaments reach 60 MPa and 1 GPa, respectively. Then tensile specimens with different continuous fiber orientations and different continuous fiber placement sequences are printed and tested using a multi-material 3D printing technique. The specimens are characterized before and after testing using scanning electron microscopy and microscopy to assess the porosity and failure mechanisms of specimens with different configurations. The results show that the mechanical properties of the printed parts are much lower than those of the print filaments, which proves the serious negative impact of the material extrusion process on the mechanical properties of the printed structural parts. Finally, an analytical method for predicting the elastic behavior of printed composites is developed by introducing the porosity factor in the volume-averaged stiffness model. For continuous and short fiber filled composites with different fiber contents and printing orientations, the predicted results are in agreement with the experiments, and the prediction error is greatly reduced from 30 % to <5 %.

Evaluation of a Volume-Averaged Species Transport Model with Micro-Macro Coupling for Breakthrough Curve Prediction.

In porous water filters, the transport and entrapment of contaminants can be modeled as a classic mass transport problem, which employs the conventional convection-dispersion equation to predict the transport of species existing in trace amounts. Using the volume-averaging method (VAM), the upscaling has revealed two possible macroscopic equations for predicting contaminant concentrations in the filters. The first equation is the classical convection-dispersion equation, which incorporates a total dispersion tensor. The second equation involves an additional transport coefficient, identified as the adsorption-induced vector. In this study, the aforementioned equations were solved in 1D for column tests using 3D unit cells. The simulated breakthrough curves (BTCs), using the proposed micro-macro-coupling-based VAM model, are compared with the direct numerical simulation (DNS) results based on BCC-type unit cells arranged one-after-another in a daisy chain manner, as well as with three previously reported experimental works, in which the functionalized zeolite and zero-valent iron fillings were used as an adsorbent to remove phosphorous and arsenic from water, respectively. The disagreement of VAM BTC predictions with DNS and experimental results reveals the need for an alternative closure formulation in VAM. Detailed investigations reveal time constraint violations in all the three cases, suggesting this as the main cause of VAM's failure.

A multi-field couped model for controlling segregation in twin-roll casting

So-called center and edge segregation are formed during twin-roll casting, which cannot be effectively eliminated in subsequent rolling and heat treatment, thus significantly diminishing the quality of final product. Numerous studies have been conducted on such segregations so far, unfortunately, without fully considering conservation of mass, momentum, energy and solute throughout the process, as well as synergy of thermodynamics and kinetics for interfacial migration. Accordingly, a two-phase Eulerian-Eulerian volume-averaged model is herein established to predict the segregation upon twin-roll casting for Al–Si alloys, which, for the first time, addresses effects due to the interfacial migration and the solute diffusion on the solute transfer process, and so-called correlation between thermodynamics and kinetics (i.e., thermo-kinetic correlation). On this basis, the model accurately predicts the flow field upon twin-roll casting, where, the increased casting speed leads to the descended kissing point, thus gradually causing the decreased vortex area and the increased cooling rate. Whereas, the increased heat transfer coefficient leads to the ascended kissing point, thus gradually causing the increased vortex area and the decreased cooling rate. In the modeling, the increased thermodynamic driving force for interfacial migration is always accompanied by the decreased kinetic energy barrier, and vice versa. Furthermore, the model also accurately predicts the segregation upon twin-roll casting, where, the thermodynamic driving force increases with the casting speed, thus decreasing the central segregation, whereas, the heat transfer coefficient increases, thus enhancing the generalized stability and in turn, effectively suppressing the edge segregation. By employing a strategy involving high driving force in the early stage and high generalized stability in the late stage, it was determined that the anticipated center and edge segregation could be effectively suppressed.

Numerical modeling of wave interaction with a porous floating structure consisting of uniform spheres

Porous materials are increasingly being used in the design of floating structures in coastal and ocean engineering, but there is a lack of numerical tools that can aid in the design of a movable floating porous structure. To close this gap, the existing volume-averaged numerical model for flow interacting with a fixed porous body was extended to floating scenarios by (1) using the relative velocity in the porous friction force, (2) calculating rigid body motion using volume integral of porous body force, and (3) modifying a dynamic mesh algorithm for a mobile porous body. As a demonstration, the developed model was applied to a porous floating structure consisting of cubically packed uniform spheres. Two sets of model applications were involved. The first set considered three-dimensional flow around a fixed porous block placed beneath the free surface. The measured total force on the block under wave or steady flows was predicted accurately with an error less than 10%. The second set involved a two-dimensional wave interacting with a floating porous block representing a breakwater. For free-floating conditions, the model can accurately predict the dynamic response of the structure, including the time varying movement of its rigid body and the mean drift. For mooring-restrained conditions, the mooring force and wave transmission coefficients were also predicted well with an error less than 20%. The proposed numerical approach can be applied to other floating structures with a rigid volumetric porous body. Future research is also required to study the microscopic pore flows, upon which more detailed parameterization of the porous media can be derived.

Features of a fusible solid fuel gasification in a low-temperature gas generator satisfying the assumptions of the volume average model

A solid porous fuel gasification in the low-temperature gas generator for a high-speed flying vehicle is studied numerically in the framework of volume averaged model. It is shown that during the filtration of hot gas (∼ 1000-1200 °C) through the porous medium the significant part of the solid fuel has relatively high temperature (∼ 180-450 °C) for a long time (10–70 s). This suggests that melting of common materials such as polyethylene or polymethylmethacrylate considered in the literature as the fuel for the gas generator is possible, if the porous media satisfies assumptions of the volume averaged approach. Based on the results of a parametric study, possible ways to prevent the melting of the gasifier are discussed.

Volume-averaged modeling of multiphase solidification with equiaxed crystal sedimentation in a steel ingot

Volume-averaged modeling of multiphase solidification with equiaxed crystal sedimentation in a steel ingot

Fuel-Cell Performance and Stability during Liquid-Water Removal Cycles

Water management is crucial for achieving high-performance proton-exchange-membrane fuel cells (PEMFCs), as it helps keep the electrolyte hydrated while avoiding performance degradation and cell shutdown due to liquid-water accumulation. According to ex-situ measurements [1-3], accumulation of liquid water in gas-diffusion layers (GDLs) of PEMFCs is a transient process that can be accompanied with oscillations in capillary pressure and saturation. To understand how liquid-water accumulation and drainage impact PEMFC performance hysteresis and stability, a transient cell-level model that accounts for electrode structure and composition and is computationally efficient is needed. A number of volume-averaged models that describe the electrode structure through pore-size distribution have been developed in the past [4-7], but they are steady-state and thus cannot predict dynamic PEMFC performance. Existing transient models have also not been used to analyze cyclic liquid-water accumulation [8-10].In this work, a transient two-phase 2D PEMFC model is developed in the open-source fuel-cell modeling software OpenFCST [11] and applied to analyze PEMFC performance hysteresis and stability during liquid-water accumulation and drainage cycles. The model accounts for the electrode structure through a mixed-wettability pore-size-distribution framework and incorporates a novel dynamic boundary condition to describe the experimentally observed cyclic liquid-water accumulation. Results of the numerical simulations are compared to transient current-density and resistance data at two polarization scan rates and during voltage steps measured with an in-house single-channel cell at multiple operating conditions. This work demonstrates how high breakthrough pressure and rapid liquid-water removal from GDLs may cause highly unstable fuel-cell operation with strong hysteresis and oscillations in the polarization curve (such as those in the figure below) that closely resemble experimental reports [12]. Numerical simulations also reveal the existence of a scan rate that maximizes polarization hysteresis due to a match between the time scale of GDL flooding and of a complete voltage sweep. Fast-scan polarization sweeps are shown most suitable for detecting catalyst-layer flooding that depends on its wettability and occurs within single seconds in contrast to GDL flooding that takes hundreds of seconds. Overall, this work brings more attention to the transient analysis of fuel-cell performance under wet conditions.Figure: Polarization-curve oscillations caused by cyclic liquid-water removal from the cathode GDL. Similar fluctuations have been experimentally observed in [12]. Transient graphs show the dynamics of current density and cathode GDL saturation. Operating conditions are 60 °C, 90% RH, 1.5 atm; scan rate is 0.5 mV/s. References J. T. Gostick et al., J. Electrochem. Soc. 157.4 (2010) C. Quesnel et al., J. Phys. Chem. C 119.40 (2015) D. Ziegler, B.Sc. thesis, Hochschule Mannheim / University of Alberta (2020) A. Z. Weber et al., J. Electrochem. Soc. 151.10 (2004) M. Eikerling, J. Electrochem. Soc. 153.3 (2006) V. Mulone and K. Karan, Int. J. Hydrog. Energy 38.1 (2013) Zhou et al., J. Electrochem. Soc. 164.6 (2017) R. J. Balliet and J. Newman, J. Electrochem. Soc. 158.8 (2011) I. V. Zenyuk et al., J. Electrochem. Soc. 163.7 (2016) A. Goshtasbi et al., J. Electrochem. Soc. 166.7 (2019) M. Secanell et al. ECS Transactions 64.3 (2014) C. Ziegler and D. Gerteisen, J. Power Sources 188.1 (2009) Figure 1

Effect of the Cross-Section of a Porous Burner on the Combustion Stability Limit of Premixed Oxy-Methane Flames.

This study presents the idea of using a porous media burner to improve the oxy-methane combustion reaction rate and broaden the stability limit. Numerical studies on the premixed combustion of CH4/O2/CO2 in a two-layer porous medium burner using a two-dimensional symmetrical volume-average model with the skeleton mechanism based on OpenFOAM. The combustion characteristics of burners with variable cross-sectional (VC) and straight cylindrical (SC) structures were compared, including stable range, temperature field, thermal cycle efficiency, and CO emissions. It is confirmed that the cross-sectional structure is effective for broadening the stable range, and the broadening rate is more than 4 times. As more heat is lost from the outlet due to the increased cross-section, the solid-phase temperature of VC is lower than that of SC. As a result, the flame temperature of VC will be about 200 K lower than that of SC under the condition of insulated walls. It also leads to a reduction of about 2% in thermal cycle efficiency compared with SC. Meanwhile, it is shown that the VC structure slightly increases CO emissions at low thermal power but is smaller than that of SC at high power. When the wall heat loss of the burner is considered, the VC structure is still effective for improving the stable range. In addition, the VC structure slightly affects flame tilt and temperature uniformity at the burner outlet when the thermal power is large. It is confirmed that the VC structure is effective for increasing the power adjustment range and reducing pollutant emissions at high power.

Numerical evaluation of heat-triggered drug release via thermo-sensitive liposomes: A comparison between image-based vascularized tumor and volume-averaged porous media models

In this study, two numerical models for controlled drug delivery through thermo-sensitive liposomes under mild microwave hyperthermia are compared. The first model simulates blood flow directly in a realistic tumor vascular network. The second one, at the volume-averaged macro scale, employs the porous media model. Using the vascular-scale approach is helpful for these kind of applications since the porous media assumption does not allow to account for local vascular blood flow, which is a crucial variable to quantify outputs like local drug accumulations. In the present study, the vascularized tumor model is obtained from a photo-acoustic scan and used as computational domain. The vascular-scale model considers heat and drug transfer in vascular and intracellular spaces, employing energy equations for heat transfer and a three-compartment convection-diffusion model for drug diffusion. Blood flow is simulated using the momentum equation for non-Newtonian fluids. The volume-averaged model shows good agreement with the vascular-scale model even if provides less precise information about temperature field distribution and drug localization within the tumor tissue. Also, the volume-averaged model underestimates Fraction of Killed Cells by 3 % post-treatment, presenting a more uniform temperature field and drug distribution compared to the vascular-scale one. Despite these small limitations, the porous media model is reliable for obtaining average results, while the vascular-scale model is suitable when more accurate results are needed. This study's findings aid those who aim to predict hyperthermia-driven drug delivery, suggesting the vascular-scale approach for precision and accuracy and the volume-averaged approach for average outcomes.

Charging characteristics of finned thermal energy storage tube under variable rotation

Solar photovoltaic-thermal (PVT) systems effectively offset the drawbacks of intermittent solar and low photovoltaic conversion efficiency. Thermal energy storage (TES) tanks of PVT systems with high charging efficiency and consistent thermal safety might achieve efficient utilization of solar energy for building. A new variable rotational strategy has been proposed to optimize the charging characteristics for TES tubes, taking into consideration the non-uniformity of melting. A series of simulations based on the volume-averaged model are conducted to investigate the thermal energy storage property of TES tubes under variable rotary mechanism. Qualitative and quantitative comparisons are made between variable rotation (ω = 1.5–0.5, 1.5–1.0, 1.5–2.0 rad·s−1), constant rotation (ω = 1.5 rad·s−1), and stationary systems. The focus of the comparison is melting efficiency, temperature distribution, and natural convection. The results indicate that rotation effectively shortens charging time, with a 57.62% and 15.73% reduction with variable rotary mechanisms of 1.5–1.0 rad·s−1 when compared with stationary and constant rotating tubes. Meanwhile, in the final moment, the greatest medium-temperature (55–65 °C) proportion of 90.58% and less low-temperature (25–55 °C) and high-temperature (65–70 °C) paraffin occupation of 4.67% and 4.75% could be obtained, reflecting the completed latent heat storage and stable thermal safety. The optimal variable rotation achieves improvements of 47.84% and 106.73% in time-integral Grashof number (Gr) and heat storage rate, compared with traditional stationary tubes.

Hydrodynamic analysis of the combined structure of offshore monopile wind turbine foundation and aquaculture cage

Reliable prediction in hydrodynamics is crucial for the design of marine structures, especially for columnar structures such as offshore wind turbine foundations. A combined structure of fixed monopile wind turbine foundation and bottom-seated circular aquaculture cage is proposed for efficiently utilizing the offshore space and reducing the wave forces on the monopile wind turbine foundation. A further developed volume-averaged numerical model is adopted to simulate the pressure drop due to clean and fouled cages by introducing the porous media theory, in which the coefficients of source terms are determined by the least square method based on hydrodynamic force on cage from existing experimental data. The Stokes fifth-order wave theory is considered reliable in simulating the nonlinear waves on slopes and corresponding inline wave force by comparing with multiple existing experimental data. Then the effects of incident wave steepness, seabed bottom slopes and cage diameters on simulation results are analyzed. Results reveal the fouled cage can effectively reduce inline wave forces relative to the clean cage of the combined structure, accompanying a significant resonance effect of wave motion inside the cage.

Improving volume-averaged simulations of matrix-stabilized combustion through direct X-ray µCT characterization: Application to NH3/H2-air combustion

Porous media combustion (PMC) relies on internal heat recirculation in an open-cell ceramic foam matrix to enhance the flame speed of fuels with poor combustion properties. Volume-averaged simulations are often used to study the combustion performance and pollutant emissions of such systems. However, due to the varying complexity of matrix geometries found in practical burners, as well as the wide range of closure models for the constitutive relations of the solid phase, contradicting statements about the predictive accuracy of these volume-averaged models can be found in the literature. In this work, we propose an open-source modeling framework for accurate volume-averaged PMC simulations by using first-principles methods to determine effective properties used in closure models. This framework relies on adequately characterizing the topology of the solid matrix, using commonly available X-ray computed microtomography. With this approach, significant improvements in accuracy are reported compared to empirical models from the literature. The framework based on first-principle evaluations of constitutive relations is compared against experimental measurements conducted on an interface-stabilized burner operated with premixed NH3/H2-air. The model shows good agreement for exhaust gas composition and stability limits. The proposed simulation framework performs significantly better than state-of-the-art techniques that employ commonly used empirical correlations for effective matrix properties. Statement of SignificanceWe present a new open-source simulation framework for improved characterization of porous media combustion. By utilizing µCT techniques, accurate effective matrix properties can be determined from first-principle simulations. These effective properties are used in closure models for 1D volume-averaged reacting flow simulations using appropriate sub-models for heat recirculation. This modeling framework is able to reliably predict stability limits while conventional closure models yield erroneous trends. Assessment of the resulting modeling framework is performed using experiments with exhaust gas characterization performed on a NH3/H2-air porous media burner.

Numerical simulation of the interaction between nonlinear sloshing flow and side-mounted perforated baffles

This study investigates the effect of nonlinear strengths of sloshing on damping performance of side-mounted perforated baffles and discusses the influence of perforated baffles on resonance characteristics of sloshing flow. The model performance score is adopted to evaluate the nonlinear strength of water sloshing. A microscopic model and a volume-averaged macroscopic model are established to simulate the sloshing flow in a two-dimensional (2D) water tank with perforated baffles. In addition, a series of experimental tests are conducted to validate the numerical model. By analyzing the effects of nonlinear strength on damping performance of perforated baffles, a critical nonlinear strength of Reynolds number (Re) equal to 56 000 is obtained: For nonlinear strength smaller than critical nonlinear strength, the damping effect due to perforated baffles is mainly achieved by shifting the energy density to the excitation frequency; for nonlinear strength enhanced to the critical value, further increasing nonlinear strength shows slight influence on damping performance, and an additional increase in sloshing response occurred due to violent interaction between baffles and fluid flow. Importantly, the Reynolds number corresponding to the critical nonlinear strength is independent of the properties of perforated baffles, such as porosity and mounted form. Furthermore, the resonance effect of sloshing flow can be significantly weakened with an increasing blockage of perforated baffles, accompanying negative growth of maximum force amplitude on bulkhead.

Investigation of thermal hydraulic behavior of the High Temperature Test Facility's lower plenum via large eddy simulation

A high-fidelity computational fluid dynamics (CFD) analysis was performed using the Large Eddy Simulation (LES) model for the lower plenum of the High–Temperature Test Facility (HTTF), a ¼ scale test facility of the modular high temperature gas-cooled reactor (MHTGR) managed by Oregon State University. In most next–generation nuclear reactors, thermal stress due to thermal striping is one of the risks to be curiously considered. This is also true for HTGRs, especially since the exhaust helium gas temperature is high. In order to evaluate these risks and performance, organizations in the United States led by the OECD NEA are conducting a thermal hydraulic code benchmark for HTGR, and the test facility used for this benchmark is HTTF. HTTF can perform experiments in both normal and accident situations and provide high-quality experimental data. However, it is difficult to provide sufficient data for benchmarking through experiments, and there is a problem with the reliability of CFD analysis results based on Reynolds–averaged Navier–Stokes to analyze thermal hydraulic behavior without verification. To solve this problem, high-fidelity 3-D CFD analysis was performed using the LES model for HTTF. It was also verified that the LES model can properly simulate this jet mixing phenomenon via a unit cell test that provides experimental information. As a result of CFD analysis, the lower the dependency of the sub-grid scale model, the closer to the actual analysis result. In the case of unit cell test CFD analysis and HTTF CFD analysis, the volume-averaged sub-grid scale model dependency was calculated to be 13.0% and 9.16%, respectively. As a result of HTTF analysis, quantitative data of the fluid inside the HTTF lower plenum was provided in this paper. As a result of qualitative analysis, the temperature was highest at the center of the lower plenum, while the temperature fluctuation was highest near the edge of the lower plenum wall. The power spectral density of temperature was analyzed via fast Fourier transform (FFT) for specific points on the center and side of the lower plenum. FFT results did not reveal specific frequency-dominant temperature fluctuations in the center part. It was confirmed that the temperature power spectral density (PSD) at the top increased from the center to the wake. The vortex was visualized using the well-known scalar Q-criterion, and as a result, the closer to the outlet duct, the greater the influence of the mainstream, so that the inflow jet vortex was dissipated and mixed at the top of the lower plenum. Additionally, FFT analysis was performed on the support structure near the corner of the lower plenum with large temperature fluctuations, and as a result, it was confirmed that the temperature fluctuation of the flow did not have a significant effect near the corner wall. In addition, the vortices generated from the lower plenum to the outlet duct were identified in this paper. It is considered that the quantitative and qualitative results presented in this paper will serve as reference data for the benchmark.

A volume-averaged model for acoustic streaming induced by focused ultrasound in soft porous media.

Equations describing acoustic streaming in soft, porous media driven by focused ultrasound are derived based on the assumption that acoustic waves pass through the porous material as if it were homogeneous. From these equations, a model that predicts the time-averaged flow on the macroscopic scale, as well as the advective transport of the trace components, is created. The model is used to perform simulations for different shapes of the focused ultrasound beam. For a given shape, and using the paraxial approximation for the ultrasound, the acoustic streaming is found to be linearly proportional to the applied ultrasound intensity, to the permeability of the porous material and to the attenuation coefficient, and inversely proportional to the liquid viscosity. Results from simulations are compared to a simplified expression stating that the dimensionless volumetric liquid flux is equal to the dimensionless acoustic radiation force. This approximation for the acoustic streaming is found to be reasonable near the beam axis for focused ultrasound beam shapes that are long in the axial direction, compared to their width. Finally, a comparison is made between the model and experimental results on acoustic streaming in a gel, and good agreement is found.

Contribution of Ion Energy and Flux on High-Aspect Ratio SiO2 Etching Characteristics in a Dual-Frequency Capacitively Coupled Ar/C4F8 Plasma: Individual Ion Energy and Flux Controlled.

As the process complexity has been increased to overcome challenges in plasma etching, individual control of internal plasma parameters for process optimization has attracted attention. This study investigated the individual contribution of internal parameters, the ion energy and flux, on high-aspect ratio SiO2 etching characteristics for various trench widths in a dual-frequency capacitively coupled plasma system with Ar/C4F8 gases. We established an individual control window of ion flux and energy by adjusting dual-frequency power sources and measuring the electron density and self-bias voltage. We separately varied the ion flux and energy with the same ratio from the reference condition and found that the increase in ion energy shows higher etching rate enhancement than that in the ion flux with the same increase ratio in a 200 nm pattern width. Based on a volume-averaged plasma model analysis, the weak contribution of the ion flux results from the increase in heavy radicals, which is inevitably accompanied with the increase in the ion flux and forms a fluorocarbon film, preventing etching. At the 60 nm pattern width, the etching stops at the reference condition and it remains despite increasing ion energy, which implies the surface charging-induced etching stops. The etching, however, slightly increased with the increasing ion flux from the reference condition, revealing the surface charge removal accompanied with conducting fluorocarbon film formation by heavy radicals. In addition, the entrance width of an amorphous carbon layer (ACL) mask enlarges with increasing ion energy, whereas it relatively remains constant with that of ion energy. These findings can be utilized to optimize the SiO2 etching process in high-aspect ratio etching applications.

Modelling challenges of volume-averaged combustion in inert porous media

In porous burners, the flame thickness can be smaller than the pore size, resulting in sharp and locally-anchored flame fronts. The presence of such steep gradients at the pore level is a major hurdle for the derivation of volume-averaged models, particularly for the highly non-linear reaction rates. With the intent to address the difficulties associated with volume-averaging for porous media combustion, this work makes use of 3D Pore-Level Direct Numerical Simulations including conjugate heat transfer and complex chemistry in burners of finite length. These detailed 3D simulations are compared to their 1D volume-averaged counterpart, with effective properties estimated directly on the computational domains and of identical thermo-chemical scheme. Discrepancies in terms of burning rate, profiles, and a priori analysis on the reaction rates are discussed. Various pore sizes and geometries are considered. At the pore level, it is shown that preheating, wrinkling and wall quenching are the three main factors driving the global burning rate. Importantly, hydrodynamic dispersion is shown to have an indirect role on combustion processes. From the observations of combustion at pore scale, a new closure for reaction rates based on a flamelet assumption is proposed. It accounts for flame wrinkling and eliminates the unwanted effect of hydrodynamic dispersion on burning rate.

Comparison between 1D radial and 0D global models for low-pressure oxygen DC glow discharges

This work focuses on the comparison between a zero-dimensional (0D) global model (LoKI) and a one-dimensional (1D) radial fluid model for the positive column of oxygen DC glow discharges in a tube of 1 cm inner radius at pressures between 0.5 Torr and 10 Torr. The data used in the two models are the same, so that the difference between the models is reduced to dimensionality. A good agreement is found between the two models on the main discharge parameters (gas temperature, electron density, reduced electric field and dissociation fraction), with relative differences below 5%. The agreement on other species average number densities, charged and neutral, is slightly worse, with relative differences increasing with pressure from 11% at 0.5 Torr to 57% at 10 Torr. The success of the 0D global model in describing these plasmas through volume averaged quantities decreases with pressure, due to pressure-driven narrowing of radial profiles. Hence, in the studied conditions, we recommend the use of volume-averaged models only in the pressure range up to 10 Torr.

Prediction and Validation of Flow Properties in Porous Lattice Structures

Abstract High-porosity metal foams have been extensively studied as an attractive candidate for efficient and compact heat exchanger design. With the advancements in additive manufacturing, such foams can be manufactured with controlled topology to yield highly tailorable mechanical and transport properties. In this study, a lattice Boltzmann method (LBM)-based pore-scale model is implemented to simulate the fluid flow in additively manufactured (AM) metal foams with unit cell topologies of Cube, Face Diagonal (FD)-Cube, Tetrakaidecahedron (TKD), and Octet lattices. The pressure gradient versus average velocity profiles predicted by the LBM model were validated against in-house measurements on the AM lattice samples with the same unit cell topologies. Based on the simulation results, a novel hybrid model is proposed to accurately predict the volume averaged flow properties (permeability and inertial coefficients) of the four structures. Specifically, the linear LBM (neglecting inertial forces) is first implemented to obtain the intrinsic permeability, and then the standard LBM is applied to obtain the inertial coefficient. Convenient correlations for those flow properties as a function of porosity and fiber diameter are constructed. The effects of the AM print qualities on the flow properties are also discussed. The advantages of the hybrid model compared to the polynomial fitting approach for determining flow properties are discussed and compared quantitatively. The hybrid model and presented results are valuable for flow and thermal transport evaluation when designing new metal foams for specific applications and with different materials and topologies. The presented correlations based on pore-scale simulations can also be conveniently used in volume-averaged models to predict the macroscale flow behavior in such complex structures.

Design controllable TPMS structures for solar thermal applications: A pore-scale vs. volume-averaged modeling approach

The performance of solar energy conversion processes is largely dependent on the porous material which absorbs the irradiation, provides reaction sites, and facilitates mass transport. Porous media is the key energy conversion media in solar thermal/thermochemical conversion devices. The conventional way of fabricating porous structure shows limited flexibility in varying local morphology for better transport properties. The material structure optimization is hindered by the lack of methodology for precisely tuning material structure and lagged modeling framework to assess the structure-property relationship of designed materials. This structure-property relation can be further linked to the final conversion performance of receivers/reactors using such material. The triply periodic minimum surface (TPMS) structures are promising candidates due to their flexible lattice structure design, well-defined analytical expression, and easiness in anisotropic feature introduction. Moreover, the recent development of direct ceramic 3D printing technologies enables a unique route for the fabrication of gradient functional materials with complex structures. In this study, we focused on the quantification of TPMS structures’ transport properties with a case study to show their impact on the solar thermal conversion behavior. The identified effective properties were then used in a volume-averaged multi-physics model for the evaluation of an FKS structure under direct solar illuminations. This modeling framework offers a highfidelity numerical tool for the quantification of the link between porous material structure and their effective properties in solar thermal/thermochemical conversion applications.