Three-dimensional Computational Fluid Dynamics Model Research Articles

Thermal resistance capacity model for transient simulation of capillary-box heat exchangers

In this paper, a novel simplified hybrid model is developed to simulate the transient thermal behaviors of capillary-box heat exchangers (CBHEs) buried in the seabed, serving as the front-end heat exchangers in seawater source heat pump system (SWHPs). The thermal resistance and capacity (TRC) approach is applied for derivation of the governing equations inside and outside the heat exchangers. Also, an analytical solution is integrated to model heat transfer in the seabed base. The effects of seawater seepage on the thermal performance of CBHEs are taken into account in the present model. The state-space representation is used to solve the governing equations. The model is verified against experimental data, achieving very good agreement with the mean bias error (MBE) of 7.2 %. A comparison with a three-dimensional computational fluid dynamic (CFD) model indicates that the TRC model's maximum relative error and MBE are 0.7 % and 2.0 % lower than those of the CFD model. Additionally, the ratio of the time required by the CFD and TRC models for a 31-day run was 138. These results demonstrate that the TRC model is sufficiently accurate and fast in the thermal simulation of CBHEs. Furthermore, the thermal properties of CBHEs are examined using the present model. The model in this study provides practical implications for heat transfer analysis and design improvement of CBHEs utilized in SWHPs.

Human exposure to air contaminants under the far-UVC system operation in an office: effects of lamp position and ventilation condition

The far-UVC (222 nm) system has emerged as a solution for controlling airborne transmission, yet its effect on indoor air quality, particularly concerning positioning, remains understudied. In this study, we examined the impact of far-UVC lamp position on the disinfection and secondary contaminant formation in a small office. We employed a three-dimensional computational fluid dynamics (CFD) model to integrate UV intensity fields formed by different lamp positions (ceiling-mounted, wall-mounted, and stand-alone types) along with the air quality model. Our findings reveal that the ceiling-mounted type reduces human exposure to airborne pathogens by up to 80% compared to scenarios without far-UVC. For all the lamp positions, O3 concentration in the breathing zone increases by 4–6 ppb after one hour of operation. However, it should be noted that a high concentration zone (> 25 ppb) forms near the lamp when it is turned on. Moreover, ventilation plays a crucial role in determining human exposure to airborne pathogens and secondary contaminants. Increasing the ventilation rate from 0.7 h−1 to 4 h−1 reduces airborne pathogen and secondary contaminant concentrations by up to 90%. However, caution is warranted as higher ventilation rates can lead to elevated O3 indoors, especially under conditions of high outdoor O3 concentrations.

Modeling and simulation of novel bio-emulated flow field with improved branch channels on PEMFC performance

The channel geometry of the bipolar plate in the proton exchange membrane fuel cell (PEMFC) significantly influences its performance. In this study, a bio-emulated flow field (BEFF) pattern in the bipolar plate was used based on the external carotid artery design with its branches. A comprehensive investigation was conducted using a three-dimensional (3-D) multiphase computational fluid dynamics (CFD) model on bio-emulated flow fields with branch channels of zero (BE-0), six (BE-6), and ten (BE-10). The results obtained were compared with traditional parallel flow fields (TPFF) and traditional serpentine flow fields (TSFF). The BE-10 results showed reduced pressure drop, uniform distribution of reactants, and higher power output compared with other configurations. The peak power density of BE-10 surpassed that of TPFF by approximately 45.07% and exceeded BE-0 by 10.39%. The novel BEFF structure with an increase in branch channels showed that the reactants, temperature, and static pressure distribution are more uniform with better water management.

Computational study on the impact of gasoline-ethanol blending on autoignition and soot/NOx emissions under low-load gasoline compression ignition conditions

In the present work, computational fluid dynamics (CFD) simulations of a single-cylinder gasoline compression ignition (GCI) engine are performed to investigate the impact of gasoline-ethanol blending on autoignition, nitrogen oxide (NOx), and soot emissions under low-load conditions. In order to represent the test gasoline (RD5-87), a four-component toluene primary reference fuel (TPRF)+ethanol (ETPRF) surrogate (with 10% ethanol by volume; E10) is employed. A three-dimensional (3D) engine CFD model employing finite-rate chemistry with a skeletal kinetic mechanism (including NOx sub-mechanism), adaptive mesh refinement (AMR), and hybrid method of moments (HMOM) is adopted to capture the in-cylinder combustion phenomena and soot/NOx emissions. The engine CFD model is validated against experimental data for three gasoline-ethanol blends: E10, E30 and E100, with varying ethanol content by volume. Model validation is carried out for a broad range of start-of-injection (SOI) timings (−21, −27, −36, and −45 crank angle degrees (°CA) after top-dead-center (aTDC)) with respect to in-cylinder pressure, heat release rate, combustion phasing, NOx and soot emissions. For relatively later injection timings (−21 and −27 °CA aTDC), E30 yields higher amount of soot than E10; while the trend reverses for early injection cases (−36 and −45 °CA aTDC ). On the other hand, E100 yields the lowest amount of soot among all fuels irrespective of SOI timing. Further, E10 shows a non-monotonic trend in soot emissions with SOI timing: SOI-36>SOI-45>SOI-21>SOI-27, while soot emissions from E30 exhibit monotonic decrease with advancing SOI timing. NOx emissions from various fuels follow a trend of E10>E30>E100. On the other hand, NOx emissions increase as SOI timing is advanced for all fuels, with an anomaly for E10 and E100 where NOx decreases when SOI is advanced beyond −36 °CA aTDC. Detailed analysis of the numerical results is performed to investigate the soot/NOx emission trends and elucidate the impact of chemical composition and physical properties on autoignition and emissions characteristics.

The Impact of Flow-Thermal Characteristics in Ship-Board Solid Oxide Fuel Cells

Hydrogen is increasingly recognized as a clean and reliable energy vector for decarbonization in the future. In the marine sector, marine solid oxide fuel cells (SOFCs) that employ hydrogen as an energy source have already been developed. In this study, a multi-channel plate-anode-loaded SOFC was taken as the research object. A three-dimensional steady-state computational fluid dynamics (CFD) model for anode-supported SOFC was established, which is based on the mass conservation, energy conservation, momentum conservation, electrochemical reactions, and charge transport equations, including detailed geometric shapes, model boundary condition settings, and the numerical methods employed. The polarization curves calculated from the numerical simulation were compared with experimental results from the literature to verify the model’s accuracy. The curved model was applied by enlarging the flow channels or adding blocks. Numerical calculations were employed to obtain the current density, temperature distribution, and component concentration distribution under the operating conditions of the SOFC. Subsequently, the distribution patterns of various physical parameters during the SOFC operation were analyzed. Compared to the classical model, the temperature of the curved model was reduced by 1.3%, and the velocities of the cathode and anode were increased by 4.9% and 5.0%, respectively, with a 2.42% enhancement in performance. The findings of this study provide robust support for research into and the application of marine SOFCs, and offer they insights into how we may achieving “dual carbon” goals.

Hydrogen enrichment in methanol SI engine at varying injection timing during compression stroke

This study investigates the effect of hydrogen enrichment in a direct-injected methanol-fuelled SI engine, with combustion and emission performance analysed by varying the methanol injection timing from 150° to 60° CA bTDC. A three-dimensional computational fluid dynamic model of the hydrogen-enriched methanol engine was developed to predict the engine performance, in-cylinder, and emission characteristics. The numerical model was solved using three-dimensional Reynolds-Averaged Navier Stokes, the RNG k-epsilon model for turbulence, the O'Rourke and Amsden sub-model for heat transfer, the extended Zeldovich mechanism for nitric oxide emissions, and the Hiroyasu-NSC model for soot emissions. The results indicated that retarding the injection timing resulted in a decrease in the indicated specific CO and soot emissions along with a rise in the indicated specific NOx emissions. Furthermore, hydrogen enrichment with methanol enhanced the hydroxyl radical concentration, reduced the combustion duration, reduced also the indicated specific CO and soot emissions while increasing the indicated specific NOx emissions. The study indicates that hydrogen enrichment could extend the late injection timing limit of methanol by enhancing fuel-air mixing, which improves and controls the combustion process more effectively.

Three-dimensional CFD model and aerodynamic analysis of dragonfly’s corrugated wing during gliding

Abstract The aerodynamic performance of micro air vehicles (MAVs) in low-speed flight would be improved by mimicking the dragonfly’s wing with corrugated airfoil. Research on two-dimensional corrugated airfoils has revealed that the local vortex in the corrugated structure increases the flow speed and the lift-to-drag ratio in low Reynolds numbers. However, studies seldom focus on the effects of three-dimensional corrugated structures on aerodynamics. In this paper, the mechanism of high aerodynamic performance in a dragonfly’s wing is studied considering a three-dimensional corrugated structure. The high-fidelity dragonfly forewing model is established through reverse engineering. The computational fluid dynamics (CFD) simulation is performed for the gliding corrugated wing with angles of attack (AOA) of 0°-24° under three Reynolds numbers. The flow characteristics of corrugated wings and the aerodynamic gain from it are analyzed through the simulated streamlines and pressure contours, as well as the comparison between the corrugated wing, flat wing, and flat plate. The results show that the aerodynamic characteristics of corrugated wings are generally superior. Although the introduction of the corrugated structure increases some drag, it brings a higher lift and lift-to-drag ratio than the flat wings.

Three-dimensional CFD modeling for hollow fiber gas separation membrane modules based on a dual porous media model

Mathematical modeling is crucial to optimizing the design of the hollow fiber membrane (HFM) module. A Computational Fluid Dynamics (CFD) model enables efficient three-dimensional (3D) simulations of HFM modules, which integrates a dual porous media approach, significantly reducing computational costs. Validation against alternative simulations and experimental data confirms its prediction capability (<7% deviation). Comparative analyses demonstrate the superior reliability of the 3D approach over one-dimensional simulations. Detailed velocity, pressure, and concentration profiles for the shell and tube sides reveal optimization opportunities such as dead zones within the module. Aerometric studies show significant pressure drops (>30%) for smaller fiber diameters (<100μm) on the permeate side. The model's adaptability suggests broader applications in membrane gas separation processes with modifications. This research highlights the efficacy of CFD modeling in enhancing the design and optimization of HFM modules, highlighting the cost savings of the dual porous media approximation.

Aircraft Ducted Heat Exchanger Aerodynamic Shape and Thermal Optimization

Abstract Interest in aircraft electrification and hydrogen fuel cells is driving demand for efficient waste heat management systems. Ultimately, most of the heat must be rejected to the freestream air. Ducted heat exchangers, also called ducted radiators, are the most common and effective way to do this. Engineers manually design ducted heat exchangers by adjusting the duct's shape and heat exchanger's configuration to reduce drag and transfer sufficient heat. This manual approach misses potential performance improvements because engineers cannot simultaneously consider all of the complex interactions between the detailed duct shape, heat exchanger design, and operating conditions. To find these potential gains, we apply gradient-based optimization to a three-dimensional ducted heat exchanger computational fluid dynamics (CFD) model. The optimizer determines the duct shape, heat exchanger size, heater exchanger channel geometry, and coolant flowrate that minimize the ducted heat exchanger's power requirements while rejecting enough heat. Gradient-based optimization enables the use of nearly 100 shape design variables, creating a large design space and allowing fine-tuning of the optimal design. When applied to an arbitrary, poorly performing baseline, our method produces a nuanced and sophisticated ducted heat exchanger design with five times less cruise drag. Employing this method in the design of electric and fuel cell aircraft thermal management could uncover performance not achievable with manual design practices.

Hydrogen leakage and diffusion in the operational cabin of hydrogen tube bundle containers：A CFD study

The operational cabin of hydrogen tube bundle containers (HTBCs) is susceptible to vibration and fatigue loads during transportation, which may result in hydrogen leakage and deflagration accidents. In this paper, a three-dimensional (3D) computational fluid dynamics (CFD) model is developed to simulate the effects of leakage diameter, leakage pressure, leakage direction and position on hydrogen leakage and diffusion in the operational cabin of HTBCs. The results demonstrate that the formation of vortex recirculation zones leads to the accumulation of flammable gas cloud (FGC) in the operational cabin. The medium leakage diameter (1.5 mm) and higher leakage pressure (20 MPa, 35 MPa) cause the operational cabin to be filled with FGC within 5 s, requiring the operational cabin to be designed with no possible ignition source inside. In addition, when hydrogen leaks from the +Z direction, the sensor responds within 0.31 s, which is 1.7 s earlier than the response time in the -Z direction, indicating that the existing sensor layout cannot meet the requirements of fast response. When the leak position (L3) is close to the vent, the FGC volume proportion at 2 s is 19.3 % and 15.88 % lower than that of L1 and L2, respectively, indicating that the leak position close to the vent can effectively slow down the accumulation of FGC. The research results have implications for the safety design of operational cabin of HTBCs, the layout of hydrogen sensors and vents, and the emergency response measures for hydrogen leakage.

Analysis of ultra-high concentration solar cells integrated with a confined microjet impingement cooling

The high solar concentration levels result in an intense influx of energy onto the concentrated triple-junction solar cell, leading to elevated temperatures that can degrade the cell’s performance and overall lifespan. This necessitates the development of innovative and efficient cooling strategies to mitigate temperature-induced losses and ensure optimal energy conversion. This study introduces an innovative micro jet impingement cooling system designed for ultrahigh concentrated solar cells, addressing the critical challenge of temperature management under high solar concentration levels. Furthermore, structural modifications were proposed, including reducing ceramic layer thickness to decrease the thermal resistance and using thermal interface materials to improve the cell temperature uniformity. A three-dimensional computational fluid dynamics model was developed to evaluate the proposed system’s performance and compare it with the traditional jet impingement heat sink. The model was verified and validated using empirical data in the literature. In addition, the Response Surface Method was employed to provide a comprehensive understanding of the system’s behavior under various weather conditions and cooling fluid effects. The results showed that concentration ratio and direct normal irradiance significantly affect the system’s performance. Moreover, the new cooling system maintains temperatures below the recommended operating temperature. Finally, the environmental assessment indicated substantial energy payback time reductions at 20,000 Suns, which is significantly lower than at 1000 Suns. The CO2 emissions analysis shows a marked potential for emissions mitigation, particularly at higher CRs and longer sunshine hours. These insights offer a comprehensive overview of ultrahigh concentrated solar cell thermal system performance, highlighting its efficiency and environmental sustainability in solar energy applications.

Catalytic and non-catalytic reactions of methane pyrolysis for hydrogen production in screening and electromagnetic levitation reactors using computational fluid dynamics

Molten metal (MM)-based methane (CH4) pyrolysis is a promising technology for clean hydrogen production with a minimal carbon footprint. Electromagnetic levitation (EMLR) and screening reactors (SR) were used to investigate the intrinsic catalytic reaction kinetics of MM-based CH4 pyrolysis. Heterogeneous catalytic reactions in the SR and EMLR occurred on the surfaces of the MM crucible and droplet, respectively. However, a homogenous non-catalytic reaction may occur just above the surface because the CH4 gas is preheated by the high-temperature MM. This study presents three-dimensional (3D) computational fluid dynamics (CFD) model coupled with chemical reactions to assess the extent to which non-catalytic reactions intervene in CH4 pyrolysis in both the SR and EMLR, which is difficult to measure experimentally. The CFD model validated against experimental data revealed that the non-catalytic reaction accounted for 4.0 % of CH4 pyrolysis in the CH4-Ni0.27Bi0.73 SR at 1045 °C, whereas it represented 0.02 % in the CH4-Sn EMLR at 1021 °C. This result highlights that the EMLR is suitable for studying the intrinsic catalytic reaction kinetics of MM-based CH4 pyrolysis because the non-catalytic reaction occurs to a lesser extent.

Flow flied inspired by sieve plate structure of plant leaf veins for proton exchange membrane fuel cells

Bionic flow fields have been applied in bipolar plates for proton exchange membrane fuel cells (PEMFCs) to significantly improve cell performance. This study, inspired by the vascular structure of plant leaf veins, proposes a bionic flow field that incorporates a sieve plate structure to enhance the distribution of reactants. To evaluate this design, a three-dimensional two-phase isothermal computational fluid dynamics (CFD) model was developed. The effect of hole parameters on the cell performance, such as opening ratio, hole number, and hole arrangement on the bionic sieve plate (BSP), was investigated and compared with a parallel flow field (PFF). The results showed that the peak power density of the bionic sieve plate flow field (BSPFF) was 0.991 W/cm2, 25.4% higher than that of the PFF, when the ratio of openings was 22.3%, the number of openings was 4, the sieve angle was 22.5°, and the number of sieves was 5. This BSPFF design provides insight into the development of bipolar plates to improve the performance of PEMFCs.

CFD modeling using reactions kinetics for selective hydrogenation for acetylene in a fixed-bed reactor

Selective hydrogenation of acetylene (SHA) reactions is usually performed in fixed-bed reactors (FBR). The traditional SHA kinetics, when coupled with a three-dimensional computational fluid dynamics (CFD) model, requires improvement to accurately predict SHA reaction outcomes. The SHA microdynamics, when integrated with a three-dimensional CFD model, have not been comprehensively examined. In this paper, a mathematical model based on CFD was developed to simulate the reactive flow behavior of SHA in FBR using a novel OleMax100 catalyst. In this study, a 3D catalyst bed numerical simulation of FBR with coupled microdynamics description of conjugate heat transfer and surface catalytic reaction was carried out, and the accuracy of CFD at different Reynolds numbers (Re) was verified by experimental results, which showed a high degree of agreement between the model and experimental results. The effect of the ratio of tube to pellet diameter (D/d=N) on the fluid flow characteristics and SHA reaction in the bed was investigated, and the catalyst shape was considered. The simulation results show that an increase in N significantly improves the homogeneity of the flow field and the heat and mass transfer between the phases in the FBR, which leads to an increase in the conversion efficiency of acetylene. When N was increased from 4.00 to 6.67, the conversion was enhanced by 6.7 %. Increasing the Re value affects the reactivity, and the SHA reaction exists in a reaction zone where excessive residence time exacerbates the decrease in selectivity. The adopted and novel catalyst microkinetic models can provide some theoretical guidance for the optimal design of the SHA reaction and process improvement in FBR.

Effect of direct water injection timing and quantity on combustion, performance, and emission characteristics of a homogeneous charge compression ignition engine − A CFD analysis

Direct water injection in an HCCI engine is a promising strategy to extend the operating load limit. The advantage of water injection depends on the percentage of water evaporated inside the cylinder, which further depends on the time available for the evaporation and the interaction of water spray plumes with the in-cylinder air-flow field. Here, CFD investigations are done to find the effect of water injection timing (WIT) and quantity on the load limit and performance and emission characteristics of the HCCI engine. Various WITs are considered for the study, varying from 100° CA before TDC to 20° CA before TDC. The water-to-fuel ratio is varied from 0 to 1 in steps of 0.2. A three-dimensional computational fluid dynamics model is used for the study, and CFD models used are validated from the available experimental data in the literature. The results show that WITs before 650° CA and after 670° CA are inappropriate for the HCCI engine, as heat release rate and MRPR are high in these cases. Also, the evaporation rate of water increases when the water injection is done in the region where the direction of the air-flow field is opposite to the direction of water injection. At the WIT of 60°, the HRR and MRPR are minimal, and the water vapor distribution is good compared to other cases. Also, it is seen that a water-to-fuel ratio of 0.8 gives more scope for extending the upper load limit with a marginal penalty of the IMEP; hence, a WIT of 60° with the water-to-fuel ratio of 0.8 is considered optimum. With the optimum case, it is feasible to increase the ER from 0.74 to 0.95, resulting in an increase of IMEP from 3.23 bar to 4.92 bar, which is 52.3 % more than without a water injection case at the ER of 0.74. Also, it showed minimum NOx emissions in all the water injection cases. Hence, these results address the primary challenges associated with HCCI engines, demonstrating that direct water injection can substantially improve the combustion controllability and operational range of HCCI engines.

A novel ‘3D + digital twin + 3D’ upscaling strategy for predicting the detailed multi-physics distributions in a commercial-size proton exchange membrane fuel cell stack

With the rapid development of proton exchange membrane fuel cell (PEMFC) commercialization, a comprehensive knowledge of multi-physics fields in large-scale PEMFC stacks has become ever more critical. Although conventional three-dimensional computational fluid dynamic (CFD) models have achieved great success, the application in the commercial-size stack-scale simulation remains inapplicable due to enormous computational resource requirements. Herein, based on the latest 3D CFD model, multi-physics digital twin (DT) technology and 3D stack flow distribution prediction model, a novel multi-scale upscaling prediction model is proposed. The voltage, water and thermal management characteristics of a 164-cell PEMFC stack with an active electrode area of 292.5 cm2 are studied and analyzed in details. For the analysis of commercial-size PEMFC stacks, the most comprehensive multi-physics fields are covered in this paper to date. And the results suggest that by introducing the DT technology, the time requirement of the multi-physics field prediction for unit scale prediction can be reduced by hundreds of thousands of times with a maximum global relative deviation of 1% under 10 groups of random test conditions, giving a solution from the cell scale to stack scale performance prediction, design, heat and thermal management in the PEMFC research and application.

Aerodynamic discrepancies of high-speed trains meeting within two types noise barriers: considering modeling scale ratio

PurposeSustainable urban rail transit requires noise barriers. However, these barriers’ durability varies due to the differing aerodynamic impacts they experience. The purpose of this paper is to investigate the aerodynamic discrepancies of trains when they meet within two types of rectangular noise barriers: fully enclosed (FERNB) and semi-enclosed with vertical plates (SERNBVB). The research also considers the sensitivity of the scale ratio in these scenarios.Design/methodology/approachA 1:16 scaled moving model test analyzed spatiotemporal patterns and discrepancies in aerodynamic pressures during train meetings. Three-dimensional computational fluid dynamics models, with scale ratios of 1:1, 1:8 and 1:16, used the improved delayed detached eddy simulation turbulence model and slip grid technique. Comparing scale ratios on aerodynamic pressure discrepancies between the two types of noise barriers and revealing the flow field mechanism were done. The goal is to establish the relationship between aerodynamic pressure at scale and in full scale.FindingsThe aerodynamic pressure on SERNBVB is influenced by the train’s head and tail waves, whereas for FERNB, it is affected by pressure wave and head-tail waves. Notably, SERNBVB's aerodynamic pressure is more sensitive to changes in scale ratio. As the scale ratio decreases, the aerodynamic pressure on the noise barrier gradually increases.Originality/valueA train-meeting moving model test is conducted within the noise barrier. Comparison of aerodynamic discrepancies during train meets between two types of rectangular noise barriers and the relationship between the scale and the full scale are established considering the modeling scale ratio.

Three-dimensional simulation of a proton exchange membrane fuel cell with non-integrated metal foam as flow distributor

In proton exchange membrane fuel cells (PEMFCs), the presence of metal foam in the flow channels improves the polarity diagram and increases the cell generated power, but this leads to an increase in pressure drop and, as a result, an increase in parasitic power. In this study, two types of non-integrated arrangements of metal foam are introduced, which, in addition to having high generated power, have lower pressure drop, and subsequently, lower parasitic power, and will also result in a more uniform distribution of temperature and current density. In the first arrangement, which is called IPTG, an abbreviation of Increased Porosity Toward the gas diffusion layer (GDL), the channels are divided into four equal parts with variable porosities from 0.65 to 0.95, in such a way that the metal foam layer with the highest porosity (0.95) is located in the closest area to the GDL. In the second arrangement which is called Reduced Porosity Toward the GDL, the metal foam layer with the least porosity (0.65) is located in the closest area to the GDL. To make a better comparison between the present model and the conventional one, the net power parameter with subtracting parasitic power from generated electrical power of the fuel cell is calculated for all models. The numerical study of this paper is conducted by developing a three-dimensional computational fluid dynamics model by employing user-defined functions. The results show that the net power of the PEMFC with IPTG arrangement is 3.8 % higher than that of the one with integrated metal foam arrangement with a porosity of 0.95. Also, the higher the pores per inch (PPI) of metal foam used in the non-integrated metal foam model, the greater the net power of the fuel cell, so that, the net power of PEMFC with IPTG arrangement with PPI=80 is 3.2 % more than that of PPI=20.

Performance Analysis of Finned-Tube Heat Exchanger Charged with Phase Change Material for Space Cooling

The performance of a latent heat storage unit comprised of phase change material (PCM) enclosed in a finned-tube heat exchanger was evaluated experimentally and theoretically to determine its viability to condition a space during summer. The internal and external design conditions of a typical building were selected and analyzed to determine the type of PCM, and the phase change temperature required for space cooling. Subsequently, a PCM of Plus-ice A17 was selected and charged into a small-scale finned-tube heat exchanger. Extensive measurements were conducted on the PCM heat exchanger at different operating conditions. Meanwhile, a three-dimensional computational fluid dynamics model for the PCM heat exchanger was developed and validated with the experimental measurements and thus simulated. When the airflow velocity increases from 1.3 m/s to 6 m/s, the phase change periods decrease by 25% and 13% for the PCM charging and discharging processes respectively. When the PCM thermal conductivity increases from 1 W/(m·K) to 8 W/(m·K), the phase change periods reduce by 36.3% and 47.7% for the PCM charging and discharging processes respectively. In addition, for the same increased range of PCM thermal conductivity, the charging energy efficiency increases by 16.3%, and the discharging energy ratio drops by 7.1%.

Numerical investigation gaseous ammonia basic jet and mixing characteristics in the constant volume vessel

Ammonia is considered a promising zero-carbon fuel under the dual-carbon targets. Systematically understanding ammonia's basic jet and mixing characteristics can promote its practical application in the engine field. Thus, a three-dimensional computational fluid dynamics (CFD) model of constant volume vessel (CVV) is established in this paper, and the forecasting abilities of the Large Eddy Simulation (LES) and Reynolds-Averaged Navier-Stokes (RANS) models for the gaseous ammonia jet process are compared based on the experimental data. Then, the ammonia jet development rule and mixing characteristics under different ambient temperatures, ambient pressures, and fuel temperatures are investigated. The variation laws of key data such as penetration distance, jet area, ammonia mass concentration distribution, equivalence ratio, and turbulent kinetic energy are obtained and analyzed. The results show that during the injection process, the airflow on the left and right sides of the jet rotates in opposite directions. A higher temperature and lower pressure within the CVV lead to a larger flammable range for ammonia jets. Moreover, the distribution of ammonia mass concentration on the radial plane shows a dynamic development law with time, meaning that it first increases to a maximum value, then gradually decreases, and finally remains stable. The farther the radial plane is from the nozzle outlet position, the lower the ammonia mass concentration. In addition, the ammonia jet's turbulent kinetic energy has a large core region and symmetrical distribution along the central axis. At the beginning of the injection, the core length is small, gradually increasing with time and finally stabilizing.