Active cooling method coupling rotational flow cooling and regenerative cooling

Abstract

The invention relates to the technical field of engine thermal management, and particularly provides an active cooling method for coupling rotational flow cooling and regenerative cooling, which comprises the following steps: arranging a rectangular auxiliary inlet channel on the side wall of a rectangular main inlet channel of a regenerative cooling channel, and constructing a rotational flow regenerative cooling channel; when the rectangular main inlet channel performs regenerative cooling, coolant induced rotational flow is injected into the rectangular auxiliary inlet channel, and the rotational flow strength is regulated and controlled by adjusting the ratio of the flow of the coolant induced rotational flow to the total flow; and a cooling process simulation model is constructed, and the cooling effect of the rotational flow regeneration cooling channel is analyzed based on simulation calculation data. By means of the innovative rotational flow regeneration cooling channel design, controllable rotational flow is introduced into a traditional regeneration cooling channel, and the problems that a scramjet engine combustion chamber is insufficient in heat transfer in the extremely high-temperature environment and the fuel chemical heat sink utilization rate is low are solved. And fluid mixing is enhanced through rotational flow, a thermal boundary layer is destroyed, the fuel pyrolysis reaction is optimized, and efficient heat exchange is achieved.

Description

Technical Field

[0001] The present invention relates to the technical field of engine thermal management, and in particular to an active cooling method coupling swirl cooling and regenerative cooling. Background Art

[0002] Hypersonic vehicles face extreme thermal environments during flight. Scramjet combustion chambers, in particular, can reach temperatures exceeding 3000°C, far exceeding the tolerance limits of existing structural materials. Regenerative cooling, a highly efficient thermal management technology, utilizes endothermic hydrocarbon fuel (EHF) as both propellant and coolant, enabling thermal protection of the combustion chamber without significantly increasing mass.

[0003] However, conventional regenerative cooling systems suffer from insufficient heat transfer, leading to excessively high wall temperatures and low fuel chemical heat sink utilization, limiting their application in scramjets. Existing heat transfer enhancement technologies, such as turbulence generators and flow pattern optimization, while improving heat transfer performance to some extent, still suffer from drawbacks such as coking and blockage, and complex manufacturing. Therefore, an efficient and stable heat transfer enhancement method is urgently needed to overcome this "thermal bottleneck." Summary of the Invention

[0004] In view of this, the present invention proposes an active cooling method that couples swirl cooling and regenerative cooling to solve the problems existing in the above-mentioned prior art.

[0005] To achieve the above objectives, the present invention proposes an active cooling method that couples swirl cooling and regenerative cooling, comprising:

[0006] A rectangular secondary inlet channel is provided on the side wall of the rectangular main inlet channel of the regenerative cooling channel to construct a swirl regenerative cooling channel;

[0007] While the rectangular main inlet channel is performing regenerative cooling, a coolant is injected into the rectangular secondary inlet channel to induce a swirl flow, and the swirl intensity is controlled by adjusting the ratio of the flow rate of the coolant-induced swirl flow to the total flow rate;

[0008] A cooling process simulation model is constructed, simulation calculation data is obtained based on the cooling process simulation model, and the cooling effect of the swirl regeneration cooling channel is analyzed based on the simulation calculation data.

[0009] Furthermore, the rectangular main inlet channel and the rectangular secondary inlet channel are arranged asymmetrically, and the inlet flow direction of the rectangular secondary inlet channel points to the heating surface of the rectangular main inlet channel.

[0010] Furthermore, the size of the rectangular secondary inlet channel is smaller than half the size of the rectangular main inlet channel.

[0011] Furthermore, the ratio of the flow rate of the coolant-induced swirl flow to the total flow rate ranges from 0 to 100%.

[0012] Furthermore, the cooling process simulation model includes a fuel cracking reaction geometry model, a conventional regenerative cooling channel geometry model, and a swirl regenerative cooling channel geometry model;

[0013] The fuel cracking reaction geometric model is a two-dimensional rotationally symmetrical circular tube structure and does not contain wall solid domains; the traditional regenerative cooling channel geometric model and the swirl regenerative cooling channel geometric model are both three-dimensional rectangular cooling channel structures, with the side surfaces set as circulation boundaries and the solid domains on both sides taken as half the wall thickness;

[0014] All three geometric models are extended by 100 mm at the entrance and exit.

[0015] Furthermore, the process of obtaining simulation calculation data based on the cooling process simulation model includes:

[0016] Constructing high-density, medium-density, and low-density CFD calculation grids for each geometric model in the cooling process simulation model, and performing grid independence verification based on the CFD calculation grids;

[0017] Performing cracking reaction model verification based on the fuel cracking reaction geometric model and the corresponding CFD calculation grid to confirm whether the fuel cracking reaction model meets the accuracy requirements;

[0018] The fuel verified by the cracking reaction model is used to perform steady-state simulation calculations on the traditional regenerative cooling channel geometric model and the swirl regenerative cooling channel geometric model to calculate flow heat transfer data and thermal cracking data of the fuel.

[0019] Furthermore, after the mesh independence verification, a CFD calculation mesh with a near-wall y+ value less than 1.0 is selected for the subsequent verification process.

[0020] Furthermore, during the verification of the cracking reaction model, n-decane was used as fuel, a one-step total package reaction model was adopted to describe the thermal cracking reaction of n-decane, and the SST k-ω turbulence model was used; a pressure-based solver was used in conjunction with the SIMPLEC algorithm for calculation, a first-order upwind scheme was used for turbulent kinetic energy and turbulent dissipation rate, and a second-order upwind scheme was used for the remaining variables; and physical property calculation software was used to calculate the thermal properties of n-decane and its thermal cracking products.

[0021] Furthermore, during the steady-state simulation calculation, the ratio of the coolant-induced swirl flow to the total flow in the swirl regeneration cooling channel geometric model is set to 20% and 50% respectively, and the remaining boundary conditions are the same as those of the traditional regeneration cooling channel geometric model.

[0022] Compared with the prior art, the present invention has the following beneficial effects:

[0023] The present invention adopts innovative swirl / regeneration coupled cooling technology, which significantly improves the heat transfer efficiency of the regenerative cooling system and the utilization rate of the fuel chemical heat sink. The swirl destroys the thermal boundary layer and suppresses thermal stratification, reducing the average wall temperature by 67K, the maximum wall temperature by 25K, and the peak Nusselt number (Nu) by 60%; the swirl enhances the mixing in the near-wall area, and the total conversion rate of n-decane is increased from 10.7% to 12.9%, and the proportion of chemical heat sink is increased by 22.2%. The system accurately controls the swirl intensity through flow ratio regulation (20%-50%) to balance the local temperature rise and the overall heat transfer deterioration. Compared with traditional technologies, the present invention provides an efficient and reliable thermal protection solution for hypersonic aircraft to break through the "thermal bottleneck". BRIEF DESCRIPTION OF THE DRAWINGS

[0024] Various other advantages and benefits will become apparent to those skilled in the art upon reading the detailed description of the preferred embodiment below. The accompanying drawings are for illustration purposes only and are not to be considered as limiting the present invention. The same reference symbols are used throughout the drawings to represent the same components. In the drawings:

[0025] Figure 1 Schematic diagram of the structure of the swirl regeneration cooling channel in an embodiment of the present invention;

[0026] Figure 2 Schematic diagram of the process of the active cooling method of coupling swirl cooling and regenerative cooling in an embodiment of the present invention;

[0027] Figure 3 Schematic diagram of the geometric model structure of the fuel cracking reaction in an embodiment of the present invention;

[0028] Figure 4 This is the Nusselt number distribution diagram of the cyclone regeneration cooling channel described in an embodiment of the present invention, where (A) is the subcritical section, (B) is the supercritical section, and (C) is the thermal cracking section. DETAILED DESCRIPTION

[0029] The exemplary embodiments of the present disclosure will be described in more detail below with reference to the accompanying drawings. Although exemplary embodiments of the present disclosure are shown in the accompanying drawings, it should be understood that the present disclosure can be implemented in various forms and should not be limited by the embodiments set forth herein. On the contrary, these embodiments are provided to enable a more thorough understanding of the present disclosure and to fully convey the scope of the present disclosure to those skilled in the art. It should be noted that, unless there is a conflict, the embodiments of the present disclosure and the features described in the embodiments can be combined with each other. The present invention will be described in detail below with reference to the accompanying drawings and in combination with the embodiments.

[0030] The current heat exchange intensity of the regenerative cooling channel is far from meeting the application requirements, and heat transfer deterioration is prone to occur in the channel, seriously affecting the thermal safety of the aircraft. To this end, this embodiment develops a new regenerative / swirl coupled cooling method, which effectively improves the convective heat transfer intensity in the regenerative cooling channel (assessed according to the Nusselt number), and further improves the cracking rate and chemical heat sink of endothermic hydrocarbon fuels (the cracking rate is increased from 10.7% to 12.9%). In addition, the current regenerative cooling channels are mostly rectangular single-inlet channels. Under normal circumstances, swirl, a flow phenomenon that can enhance heat exchange, cannot appear in the channel. To this end, this embodiment develops a swirl / regenerative cooling channel structure, which successfully achieves the induction of swirl flow in the cooling channel.

[0031] like Figure 1 As shown, the swirl / regeneration cooling channel structure described in this embodiment includes a rectangular main inlet channel (Inlet-1) and a rectangular secondary inlet channel (Inlet-2). The two groups of channels are arranged asymmetrically, and the secondary inlet channel is located on the side wall of the main inlet channel. The size is 1 / 4 of the main inlet channel (1 / 4 in the current case, but less than 1 / 2 is sufficient in practice). The swirl is induced by injecting coolant into the lateral secondary inlet channel.

[0032] Based on the above channel structure, the active cooling method of coupling swirl cooling and regenerative cooling proposed in this embodiment is as follows: Figure 2 As shown, the steps include:

[0033] While the rectangular main inlet channel is performing regenerative cooling, coolant is injected into the rectangular secondary inlet channel to induce swirl, and the swirl intensity is regulated by adjusting the ratio of the flow rate of the coolant-induced swirl to the total flow rate; the flow ratio range is 0-100%, and the larger the flow ratio, the greater the swirl intensity.

[0034] A cooling process simulation model is constructed, simulation calculation data is obtained based on the cooling process simulation model, and the cooling effect of the swirl regeneration cooling channel is analyzed based on the simulation calculation data.

[0035] The cooling process simulation model includes: a fuel cracking reaction geometry model, a traditional regenerative cooling channel geometry model and a swirl / regenerative cooling channel geometry model;

[0036] Based on the fuel cracking reaction geometric model, traditional regenerative cooling channel geometric model, and swirl / regenerative cooling channel geometric model established in the above steps, three types of mass fluid mechanics CFD calculation grids with high density, medium density, and low density are established for each geometric model to perform grid independence verification; when establishing the CFD calculation grid, the grid selected after grid independence verification meets the following conditions: the y+ value near the wall is less than 1.0.

[0037] The fuel cracking reaction geometric model is used and the CFD calculation grid divided by the above steps is used to verify the fuel cracking reaction model. The verification includes: performing steady-state simulation calculations using boundary conditions consistent with the fuel cracking reaction model to obtain simulation results; the boundary conditions are the same as actual experimental conditions of the fuel cracking reaction, and the simulation results are compared with the experimental results to confirm whether the fuel cracking reaction model meets the predetermined accuracy requirements; in the verification of the fuel cracking reaction geometric model, n-decane is used as the fuel, a one-step total package reaction model is used to describe the thermal cracking reaction of n-decane, and an SSTk-ω turbulence model is used; a pressure-based solver is used in conjunction with the SIMPLEC algorithm for calculation, a first-order upwind scheme is used for turbulent kinetic energy and turbulent dissipation rate, and a second-order upwind scheme is used for other variables; and the thermal physical properties of n-decane and its thermal cracking products are calculated using the physical property calculation software SUPERTRAPP.

[0038] Since the thermal cracking reaction of n-decane needs to be considered during the simulation process, it is necessary to establish a reasonable cracking reaction model. In order to verify the established n-decane thermal cracking reaction model, it is necessary to establish a geometric model corresponding to the cracking experiment. In order to analyze the optimization effect of the swirl / regeneration cooling channel proposed in the present invention through numerical simulation, it is necessary to establish a swirl / regeneration cooling channel geometric model. In order to compare the cooling effect with that of a traditional regeneration cooling channel under the same operating conditions, it is necessary to establish a traditional regeneration cooling channel geometric model of the same size.

[0039] The present invention involves three geometries, namely the fuel cracking reaction geometry model required for n-decane cracking model verification (hereinafter referred to as channel A) (such as Figure 3 As shown in the figure), the traditional regenerative cooling channel geometric model (hereinafter referred to as channel B) and the swirl / regenerative cooling channel geometric model (hereinafter referred to as channel C). Since the cracking reaction model verification does not need to consider the influence of buoyancy and wall material, in order to reduce the resources required for simulation calculations, channel A is a two-dimensional rotational axisymmetric circular tube and does not contain wall solid domains. Channel A has a diameter of 0.5mm and a length of 375mm. Channel B and channel C are both three-dimensional rectangular cooling channels, and their channel cross-sections are all squares with a width of 2mm, a wall thickness of 1mm, and a length of 100mm. The second inlet of channel C is set at the inlet, the inlet flow direction points to the heating surface, and the side length of the second inlet is 1 / 4 of the first inlet. In particular, since the regenerative cooling channels are mostly arranged in parallel, in order to consider the influence of the fluids on both sides, the sides of channel B and channel C are set as circulation boundaries and the solid domains on both sides are taken as half the wall thickness. In order to enable the flow of fluid to be fully developed and optimize the convergence of the solution process, all three channels have been extended at the inlet and outlet, each extended by 100mm.

[0040] Divide the mesh and verify mesh independence:

[0041] Channel A in the present invention is a two-dimensional rotationally symmetric structure that only considers the fluid domain, and an encrypted grid is used at its flow edge. Channel B and Channel C both contain solid domains and fluid domains. The solid domain uses a relatively sparse grid, while the fluid domain uses a high-quality grid and an encrypted grid is used close to the solid wall. In particular, considering the computational stability of channel C, the grid is encrypted at the fluid intersection. In the grid independence verification, three sets of grids with different densities are established for each of the three channels. For channel A, the difference between the three sets of grids is the number of radial grid divisions, and the number of nodes is 15, 25 and 35 respectively. The rest of the grid parameters are consistent, including: the initial grid height of the fluid domain adjacent to the wall is 1μm, the growth rate along the vertical wall direction is 1.05, and the grid is evenly divided into 10,000 nodes along the axial direction. For channel B, the difference between the three sets of grids is the number of radial grid divisions in the fluid domain, and the number of nodes is 51×51, 61×61 and 71×71 respectively. All other mesh parameters remained the same, including: radial solid domain segmentation with 0.1 mm segments, an initial mesh height of 1 μm on the fluid wall, a growth rate of 1.05 along the perpendicular wall, and axial segmentation of 1 mm segments for both the solid and fluid domains. For channel C, the three mesh sets differed in the number of radial mesh divisions in the fluid domain: 51×51, 61×61, and 71×71 nodes, respectively. All other mesh parameters remained the same, including: radial solid domain segmentation with 0.1 mm segments, an initial mesh height of 1 μm on the fluid wall, a growth rate of 1.05 along the perpendicular wall, and axial segmentation of 1 mm segments for both the solid and fluid domains in the adiabatic section. The heating section after the second inlet was segmented with an initial mesh height of 10 μm, a growth rate of 1.05 along the wall, and 10,001 mesh nodes.

[0042] The grid independence of the divided grid is verified. When the grid quality does not have a significant impact on the simulation results (relative error is less than 5%), a medium-quality grid is used to carry out subsequent simulation work in order to take into account the computational efficiency.

[0043] Validation of the cracking reaction model:

[0044] The divided grid verifies the cracking reaction model. The present invention adopts a one-step total package reaction model to describe the cracking of n-decane. The following describes the two-step cracking model mechanism:

[0045] The cracking reaction equation of n-decane is as follows:

[0046]

[0047] The cracking reaction rate of n-decane can be expressed as:

[0048]

[0049] The chemical reaction rate constant kpy is described by the Arrhenius formula as follows:

[0050]

[0051] The pre-exponential factor Apy and activation energy Ea in formula (3) are: 1.6×1015s -1 and 2.63655×108 kcalmol -1 .

[0052] During the transcritical process, the thermophysical properties of n-decane undergo significant changes. Furthermore, since the cracking of n-decane produces a variety of reaction products, the composition of the fluid changes, further altering the overall thermophysical properties of the fluid. Accurate description of thermophysical properties is crucial for calculating flow and heat transfer characteristics. The present invention uses the property calculation software SUPERTRAPP to calculate the thermophysical properties of n-decane and its thermal cracking products. Specifically, the properties relevant to flow and heat transfer include density, specific heat at constant pressure, thermal conductivity, and viscosity.

[0053] This paper uses ANSYS Fluent as the platform for CFD simulation to carry out numerical simulation calculations. The following details the relevant settings for CFD simulation:

[0054] The SST k-ω turbulence model was selected to accurately simulate the turbulence in the cooling channels. To ensure that the initial mesh layer fell within the viscous substratum, the y+ value of the near-wall mesh was checked and controlled to be less than 1.0, ensuring a sufficiently accurate approximation of the flow and heat transfer of hydrocarbon fuels in the cooling channels.

[0055] Enable the component transport model, select the Volumetric option, and define the n-decane cracking reaction.

[0056] The solver uses a pressure-based solver and the SIMPLEC algorithm. Except for the turbulence kinetic energy and turbulence dissipation rate related to turbulence, which are calculated in the first-order upwind format, the rest use the second-order upwind format to obtain more accurate simulation results.

[0057] In the verification of the cracking reaction model, in order to compare with the experimental results, the same distributed temperature boundary as that in the cracking reaction model related experiments was used, and the other boundary conditions were consistent with the experiment, such as inlet flow rate, inlet temperature and pressure.

[0058] Perform steady-state simulation calculations:

[0059] The fuel used in the cracking reaction model is used to perform steady-state simulation calculations in the traditional regenerative cooling channel geometry model and the swirl / regenerative cooling channel geometry model to calculate the flow heat transfer data and the thermal cracking data of the fuel. In order to obtain the cooling effect of the swirl / regenerative cooling channel under different boundary conditions, the present invention sets two flow ratios, and the value of gravity acceleration is set to 9.8ms. -1 The heat flux is uniformly distributed, and its value is 0.8MWm -2 In the simulation calculation, the total mass flow rate is kept at 0.5gs -1 The flow ratio is defined as the ratio of the second inlet mass flow rate to the total mass flow rate, with the two flow ratios being 20% ​​and 50% respectively. To compare the cooling effect with that of the traditional regenerative cooling channel, the flow and heat transfer characteristics of n-decane in the traditional regenerative cooling channel were simulated under the same boundary conditions.

[0060] Post-processing and analysis of simulation results:

[0061] According to the above boundary conditions, steady-state simulation calculations are carried out. When the residual is less than 1.0×10 -5 When , the simulation results are considered to have converged. At this point, the simulation ends and enters the post-processing and result analysis stage.

[0062] Post-processing mainly includes analysis of streamlines, pressure drop, temperature distribution, maximum temperature analysis, Nusselt number distribution, degree of n-decane cracking, heat sink analysis, and flow field distribution.

[0063] Among them, the total cracking rate and chemical heat sink obtained after post-processing are shown in Table 1.

[0064] Table 1

[0065]

[0066] As can be seen from Table 1, compared with the traditional regenerative cooling channel, the total cracking rate and chemical heat sink ratio of the swirl / regenerative cooling channel proposed in the present invention are significantly improved, indicating that the swirl / regenerative cooling channel described in the present invention significantly improves the heat transfer efficiency and fuel chemical heat sink utilization rate of the regenerative cooling system, and further enhances the thermal cooling effect of the regenerative cooling system.

[0067] The Nusselt number axial distribution diagram obtained after post-processing is shown in the figure below. Figure 4 As shown, it can be seen that when the ratio of the flow rate of the coolant-induced swirl to the total flow rate is set to 50%, the peak value of the Nusselt number is effectively improved, reaching 14 in the thermal cracking section, indicating that the swirl / regeneration cooling channel proposed in the present invention significantly improves the heat transfer effect of the system.

[0068] In summary, this invention, through its innovative swirl / regenerative cooling channel design, introduces controlled swirl into conventional regenerative cooling channels, addressing the challenges of insufficient heat transfer and low fuel chemical heat sink utilization in scramjet engine combustion chambers operating under extremely high temperatures. This swirl enhances fluid mixing, disrupts the thermal boundary layer, and optimizes the fuel pyrolysis reaction, achieving efficient heat exchange.

[0069] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and not to limit it. Although the present invention has been described in detail with reference to the above embodiments, ordinary technicians in the field should understand that the specific implementation methods of the present invention can still be modified or replaced by equivalents. Any modification or equivalent replacement that does not depart from the spirit and scope of the present invention should be covered by the scope of protection of the claims of the present invention.

Claims

1. An active cooling method that couples swirl cooling and regenerative cooling, characterized in that: include: A rectangular secondary inlet channel is provided on the side wall of the rectangular main inlet channel of the regenerative cooling channel to construct a swirl regenerative cooling channel; While the rectangular main inlet channel is performing regenerative cooling, a coolant is injected into the rectangular secondary inlet channel to induce a swirl flow, and the swirl intensity is controlled by adjusting the ratio of the flow rate of the coolant-induced swirl flow to the total flow rate; A cooling process simulation model is constructed, simulation calculation data is obtained based on the cooling process simulation model, and the cooling effect of the swirl regeneration cooling channel is analyzed based on the simulation calculation data.

2. The active cooling method of coupled swirl cooling and regenerative cooling according to claim 1, characterized in that: The rectangular main inlet channel and the rectangular secondary inlet channel are arranged asymmetrically, and the inlet flow direction of the rectangular secondary inlet channel points to the heating surface of the rectangular main inlet channel.

3. The active cooling method of coupled swirl cooling and regenerative cooling according to claim 1, characterized in that: The size of the rectangular secondary inlet channel is less than half the size of the rectangular primary inlet channel.

4. The active cooling method of coupled swirl cooling and regenerative cooling according to claim 1, characterized in that: The ratio of the flow rate of the coolant-induced swirl flow to the total flow rate ranges from 0 to 100%.

5. The active cooling method of coupled swirl cooling and regenerative cooling according to claim 1, characterized in that: The cooling process simulation model includes a fuel cracking reaction geometry model, a traditional regenerative cooling channel geometry model and a swirl regenerative cooling channel geometry model; The fuel cracking reaction geometric model is a two-dimensional rotationally symmetrical circular tube structure and does not contain wall solid domains; the traditional regenerative cooling channel geometric model and the swirl regenerative cooling channel geometric model are both three-dimensional rectangular cooling channel structures, with the side surfaces set as circulation boundaries and the solid domains on both sides taken as half the wall thickness; All three geometric models are extended by 100 mm at the entrance and exit.

6. The active cooling method of coupled swirl cooling and regenerative cooling according to claim 5, characterized in that: The process of obtaining simulation calculation data based on the cooling process simulation model includes: Constructing high-density, medium-density, and low-density CFD calculation grids for each geometric model in the cooling process simulation model, and performing grid independence verification based on the CFD calculation grids; Performing cracking reaction model verification based on the fuel cracking reaction geometric model and the corresponding CFD calculation grid to confirm whether the fuel cracking reaction model meets the accuracy requirements; The fuel verified by the cracking reaction model is used to perform steady-state simulation calculations on the traditional regenerative cooling channel geometric model and the swirl regenerative cooling channel geometric model to calculate flow heat transfer data and thermal cracking data of the fuel.

7. The active cooling method of coupled swirl cooling and regenerative cooling according to claim 6, characterized in that: After the mesh independence verification, a CFD calculation mesh with a near-wall y+ value less than 1.0 is selected for the subsequent verification process.

8. The active cooling method of coupled swirl cooling and regenerative cooling according to claim 6, characterized in that: During the verification of the cracking reaction model, n-decane was used as fuel, a one-step total package reaction model was adopted to describe the thermal cracking reaction of n-decane, and the SST k-ω turbulence model was used. A pressure-based solver was used in conjunction with the SIMPLEC algorithm for calculations. A first-order upwind scheme was used for turbulent kinetic energy and turbulent dissipation rate, and a second-order upwind scheme was used for the remaining variables. Physical property calculation software was used to calculate the thermal properties of n-decane and its thermal cracking products.

9. The active cooling method of coupled swirl cooling and regenerative cooling according to claim 6, characterized in that: During the steady-state simulation calculation, the swirl regeneration cooling channel geometric model is set such that the flow rate of the coolant-induced swirl and the total flow rate ratio are 20% and 50% respectively, and the remaining boundary conditions are the same as those of the traditional regeneration cooling channel geometric model.

Images