# A Computational Design Analysis of UAV’s Rotor Blade in Low-Temperature Conditions for the Defence Applications

Authors: Sreenadh Chevula; Sankeerth Chillamcharal; Satya Prasad Maddula
Journal: International Journal of Aerospace Engineering
 (2021)

---

## Abstract

This paper discusses about the critical situations faced by the Defence operations with drones in the area of Siachen Glacier in the Himalayas. The reasons for the structural failures in drone’s rotor blades and the low-performance efficiency of the drones at low-temperature conditions are highlighted. A possible solution to the above-mentioned problems has been addressed by introducing a new boundary design in the rotor blades and composite materials. The results which are shown in this paper are obtained by the computational analysis facility located at the Department of Aerospace Engineering, School of Technology, GITAM (Deemed to be University), Hyderabad. By mimicking the Siachen Glacier atmosphere conditions, the proposed rotor blade design has been analysed in CFD.


---

## Body

## 1. Introduction

Unmanned air vehicles (UAV) and drones play a major role in various Defence operations, namely in remote sensing, surveillance, and data collection in the critical areas. In such operations, the lifetime of the drones plays a major role. It is a fact that various parameters of the drone’s operational conditions show huge influence on the performance and structural integrity of the drones. In terms of structural components, the main component of the drones and most of the UAV are rotor blades, and in the design of these blades, the condition of the operational temperature is one of the major parameter. As a study of temperature influence and to overcome the structural failures in the rotor blades, in this paper, the case study of Siachen Glacier is under the India Army control [1], located at 35.421226° N, 77.109540° E. The glacier is nothing but moving location of densely formed ice which has been classified by their thermal characteristics, morphology, and behaviour [2]. The selected location map of the Siachen Glacier can be found in Figure 1 [3].Figure 1
Map location of Siachen Glacier [3].The temperature conditions (see Figure2) [4, 5] for the proposed work in this paper have been identified from the Siachen Glacier with latitude 35,500000 | 35, 30.000000 N | N35 30 00; longitude 77.000000 | 77 0.000000E | E0 77 00 00; and field elevation 22,000 ft/6706 m MSL. In the selected area, the winter snow fall will approximately 35 ft from the surface and the temperature can drop up to -50°C (-58°F).Figure 2
Location of Siachen Glacier [4, 5].Considering the applications of UAV and drones in glacier conditions, world widely, many studies had been conducted, such as analysing the snow depth [6–9]. A low-cost UAV has been used to test the SFM performance in topographic and lighting conditions of different types of snows which has been discussed in Cimoli et al. [6]. A comparison of photogrammetric maps collected by the UAV has been performed by Avanzi et al. [10]. A flow velocity over glacier land has been conducted by Jones C et al. [11]. Another studies relevant to identifying and analysing changes of surface features such as ice cliffs have been conducted by Steiner, J.F. et. al., Buri, P et al., and Brun, et al. [12–14],Beyond various studies conducted by worldwide researchers in terms of durability and structural integrity of the UAV’s and drones in glacier conditions, this paper focuses on the computational design analysis of rotor blades (further mentioned as propellers in this paper). An existing propeller blade design (Design-1) and proposed new designs (Design-2) and (Design-3) are created in Catia V5 (geometry information can be found in the following sections), and computational analysis and results are obtained by using ANSYS 16.0. A comparison between the results obtained from the computational analysis of Design-1 (existing rotor blade model), Design-2, and Design-3 (Proposed model) is shown in “Obtained Results.” Finally, conclusions has been plotted in “Conclusion.”


## 2. Designing the Solution

Several methods are used to calculate the performance parameters of a propeller, namely, momentum theories [15], blade element theories [16], hybrid blade element momentum theories [17, 18], and lifting line theories [19]. Without computational analysis, using some specified mathematical models will give the values of few propeller performance parameters. But in terms of high accuracy and optimization of the procedures, some new methods have been introduced which can be obtained only by using computational analysis. It is a well-known fact that a propeller acts like a rotating wing, and its cross-section is an airfoil. In terms of optimizing, a propeller performance optimizes the airfoil shapes used along the blade span, and an optimization of an airfoil shape can be challenging, if the shape is described using individual points. A common approach to describing the number of parameters used to describe an airfoil is to parameterize the shape. A study can be found with respective to the common approach can be found in Kulfan et al. [20]. This approach has been adopted in this paper and aimed at studying the propeller optimization with a coupled electric motor, and also, a hybrid blade element momentum theory has been proposed to estimate the propeller performance analysis.The lifting surface on a propeller is called a blade, and a propeller can have any number of blades. Most propellers have two to four blades. Any given point along a blade the cross-section has all the same characteristics as an airfoil such as leading and trailing edges, mean camber line, chord line, and thickness, where the blades connect is called the hub which is directly attached either to an engine or to a transmission. The root is the area between the hub and the blade, and the tip is end of the blade opposite the hub, as shown in Figure3 [21, 22] and Figure 4 [21, 22]. The propeller performance parameters, namely, the blade angle, β, is the “resultant angle between the free stream and rotational velocity components,” the effective pitch, pe, is “the distance a propeller advances in one rotation,” the geometric pitch, ge, is “the theoretical distance an element of a propeller blade would travel in one rotation”(ge is not constant along the length of blade), and the advance ratio, J, is “the ratio between the distance the propeller moves forward through one rotation and the blade diameter.”
(1)J=VnD,where n is in rotations per second, V is the velocity, and D is the diameter of the blade); the aspect ratio (AR) is the tip radius divided by the maximum blade width. A spinner is a conical or parabolic shaped fairing that is mounted over the canter of the centre of the propeller where it is connected to the hub, the thrust or driving face is the blade face is the lower surface of the propeller airfoil, and the blade back is the upper surface of the propeller. Most of the parameters can be visually found in Figure 5.Figure 3
Cross-section of the blade with velocity diagram [21, 22].Figure 4
Basic nomenclature of propellers [24].Figure 5
Cross-sections of the propeller blade [21, 22].Based on the operational areas of the drones or UAV, propellers are classified in four ways, namely, tractor, pusher propellers (based on thrust production), and fixed and variable pitched propellers (based on the pitch geometry). A variable pitch propeller’s blades can be adjusted either on the ground or during fight to allow the propeller to operate at maximum performance throughout its operation range.As mentioned above, to obtain the propeller performance parameters values with more accuracy, several computational methods are used in the design of propellers, which will utilizes some nondimensional numbers, namely, the thrust coefficientCt, power coefficient, Cq, and the efficiency η.
(2)Ct=Ftρn2D4,(3)Cq=Qρn3D5,(4)η=CtJCP,where n is the rotation speed, D is the diameter of the propeller, ρ is the density of air, and J is the advance ratio.


## 3. Blade Element Rotor Theory

In the propeller performance parameters analysis, the blade element theory is one of the simplest method. A relatively simple method of predicting the more detailed performance of a helicopter rotor using blade element theory can be found in ref [23]. From Figures 6 and 7, the difference in angle between thrust and lift directions is defined as
(5)φ=θ−α,where θ is the geometric pitch angle, α is the angle of attack with respective to the flow, and φ is the difference between the angles.Figure 6
Blade element subdivision [23].Figure 7
Resultant of force vectors and flow vectors [23].V0 is the axial flow at propeller disk, V1 is the section local flow velocity vector, summation of vectors V0 and V2, and V2 is the angular flow velocity vector.The elemental thrust and torque of this blade element can thus be written as(6)ΔT=ΔLcosφ−ΔDsinφ,ΔQr=ΔDcosφ+ΔLsinφ.Case 1.
For a single blade.
Substituting section data (CL and CD for the given α) then
(7)ΔL=CL12ρVI2c.dr,ΔD=CD12ρVI2c.dr,where ρ is the air density, and c is the blade chord so that the lift producing area of the blade element is c.dr.Case 2.
If the number of propeller blades isB then the area generating the lift will be the element area multiplied by the number of blades B covering the circuit, and the lift coefficient can be assumed to be linear.
(8)ΔT=12ρVI2cCLcosφ−CDsinφB.dr,(9)ΔT=12ρVI2cCDcosφ−CLsinφB.dr.The overall propeller thrust and torque can be found by the summation of the radial blade element values, i.e.,(10)T=ΣΔTforallelementsandQ=ΣΔQforallelements.The relation between nondimensional thrust (Cp) and torque coefficients (Ct) and the efficiency of the propeller (η) and advance ratio J can be found in Equations (1), (2), (3), and (4). The blade element theory is an iterative solution procedure and starts with an initial guess of induced flow component VI. With the guess of VI, the flow angle on the blade will be estimated, and based on the use of blade section properties, the element thrust and torque will be estimated. The approximated values obtained from the guessed values and the force balance can be used to give improved estimates of the induced velocity VI. Until the value of VI got converged to within a specified tolerance, the iteration process will be continued. (It should also be noted that convergence for a nonlinear system of equations is not guaranteed.) In this paper, linear airfoil section properties have been used to design propellers. So to obtain required results, a convergence enhancing technique (Crank-Nicholson underrelaxation) has been proposed to use.


## 4. Computational Analysis of the Propellers

Various experimental and computational analysis are available to obtain the results of the propeller analysis. Firstly, in an experimental method, the propeller blades will be tested in both static and advancing flow conditions in a wind tunnel setup. Secondly, in a computational analysis adopts three-dimensional computational fluid dynamics (CFD) simulation, utilizing the Reynolds-average Navier–Stokes (RANS) equation. CFD tools are the useful tools for propeller design and analysis. In this paper, a flow simulation (under low-temperature conditions) on propeller blades is proposed to conduct by using ANSYS 16.0. The geometry conditions of propellers Design-1 (Figure8), Design-2 (Figure 9), and Design-3 (Figure 10) can be found in Table 1. In general, the standard geometry of the propellers has a 0.254 m diameter. So in this paper, for computational design, approximate values of 0.254 m in diameter propellers have been designed in Catia V5.Figure 8
Design-1.Figure 9
Design-2.Figure 10
Design-3.Table 1
The geometry design of propellers.
NameDiameter in mPitchAmplitudeFrequencyDesign-10.2540.1778N/AN/ADesign-20.2540.1778N/AN/ADesign-30.2540.17781c20wIn the above given three designs, the pitch of 0.1778 m provides a pitch-to-diameter ratio of 0.7, which is common for an off-shelf propeller type. The value of the Reynolds number for the testing flow condition is proposed to set for 50804 which can be obtained by a rotational speed of the propeller 3008 RPM and 75% of blade station from the chord.


## 5. Computational Setup

The computational predictions have been obtained from the ANSYS 16.0 CFD solver. The approach of the multiple reference frame model (MRF) has been incorporated to obtain the predicted flow around the propellers which can be shown in the following figures (see Figures11–13). The computational domain has been split into a global, rotating region, and local domains. The rotational domain can be found where the propeller blade and hub enclosed without openings (see Figures 11–13).Figure 11
Rotation domain of Design-1.Figure 12
Rotation domain of Design-2.Figure 13
Rotation domain of Design-3.Considering the full development (without any unexpected interruptions) of the upstream and downstream flow simulations and results around the propellers, a good configuration of the flow domain plays a major role in the computational setup. In this line, for the proposed computational setup of this paper, the stationary regions inlet and outlet boundaries are located at 4D from the origin of the propellers, and the enclosure of the rotating domain is set to be 1.1D and 0.4D.Grid generation is another important factor to obtain good and accurate results in computational analysis. The rate of convergence, computational time, and performance depends on the quality of the computational grid generation. In this paper, the mesh grid has been generated using Ansys FLUENT 16.0 (see Figures14–17). The cell size of the mesh kept smaller size or range around the propeller blade in rotational domain. From the rotational domain, a gradual increment in cell size towards the stationary region has been implemented in the grid generation with enough grid refinements across the interface.Figure 14
Design-1 (generated mesh). Mesh type: unstructured tetrahedral mesh. Nodes: 120561. Elements: 81184.Figure 15
Design-2 (generated mesh). Mesh type: unstructured tetrahedral mesh. Nodes: 759610. Elements: 482119.Figure 16
Design-3 (generated mesh). Mesh type: unstructured tetrahedral mesh. Nodes: 120561. Elements: 81184.Figure 17
Design-3 (internal view of mesh). Mesh type: unstructured tetrahedral mesh. Nodes: 120561. Elements: 81184.In the computational analysis, especially in the grid generation, most of the researchers will select the unstructured tetrahedral grids because of its compatibility to solve complex geometries by discretizing. This paper also selected unstructured tetrahedral grids in the grid generation. Table2 shows the details of the grid generation of this work.Table 2
The parameters of the grid generation.
Use advance size functionCurvatureRelevance centreCoarseCurvature normal angle70.9Min size0.589180Max face size50.9180Max size117.840Growth rate1.850Minimum edge length1.507e−002


## 6. Obtained Results

Defence operational conditions for UAV such as operational altitude, geographical areas, and mission operational profiles can be found in Figure18.Figure 18
Operational mission profile.Operational Area Siachen Glacier is known as the world’s highest battlefield, and soldiers have been deployed at an elevation of up to 6000 m (19685 feet) from sea level. In this paper, based on the operational temperature with respective change in the altitudes (with gradual increment of 250) (see Figure18), the vertical drone operational area has been divided in to five zones. Complete information about the five zones, temperature, and required RPM of the Drone propellers can be found in Table 3.Table 3
Operational mission profile.
∗Altitudes (m)∗TemperatureRPM droneVelocity (m/s)Zone-10050008Zone-2250-25°C55008Zone-3500-50°C60008Zone-4750-75°C65008Zone-51000-100°C70008∗The altitude and temperatures are taken from above 6000 m from the sea level.


### 6.1. Pressure Contour for Propeller Design-1 (Zone-1 to Zone-5), Respectively

For the mimicking condition of Zone-1 with 0°C and velocity of 8 m/s obtained by 5000 RPM, the results of the pressure contour influence vary from the hub to the tailing edge of the rotor (−3.651e+003 to 8.385e+002Pa).For the mimicking condition of Zone-2 with -25°C and velocity of 8 m/s obtained by5500 RPM, the results of the pressure contour influence vary from the hub to the tailing edge of the rotor (−3.48e+003 to 9.796e0+002Pa).For the mimicking condition of Zone-3 with -50°C and velocity of 8 m/s obtained by6000 RPM, the results of the pressure contour influence vary from the hub to the tailing edge of the rotor (−4.34e+003 to 8.796e0+002Pa).For the mimicking condition of Zone-4 with -75°C and velocity of 8 m/s obtained by6500 RPM, the results of the pressure contour influence vary from the hub to the tailing edge of the rotor (−4.17e+003 to 1.129e0+002Pa).For the mimicking condition of Zone-5 with -100°C and velocity of 8 m/s obtained by7000 RPM, the results of the pressure contour influence vary from the hub to the tailing edge of the rotor (−5.33e+003 to 1.30e0+003Pa).Note: the obtained pressure contour for the propeller Design-2 (Zone-1 to Zone-5) is almost similar to Design-1. So those results are not shown in this paper.


### 6.2. Pressure Contour for Propeller Design-3 (Zone-1 to Zone-5), Respectively

For the mimicking condition of Zone-1 with -0°C and velocity of 8 m/s obtained by5000 RPM, the results of the pressure contour influence vary from the hub to the tailing edge of the rotor (−2.42e+003 to 1.434e0+003Pa).For the mimicking condition of Zone-2 with -25°C and velocity of 8 m/s obtained by5500 RPM, the results of the pressure contour influence vary from the hub to the tailing edge of the rotor (−2.74e+003 to 1.722+003Pa).For the mimicking condition of Zone-3 with -50°C and velocity of 8 m/s obtained by6000 RPM, the results of the pressure contour influence vary from the hub to the tailing edge of the rotor (−3.09e+003 to 2.03e0+003Pa).For the mimicking condition of Zone-4 with -75°C and velocity of 8 m/s obtained by6500 RPM, the results of the pressure contour influence vary from the hub to the tailing edge of the rotor (−3.65e+003 to 2.39e0+003Pa).For the mimicking condition of Zone-5 with -100°C and velocity of 8 m/s obtained by7000 RPM, the results of the pressure contour influence vary from the hub to the tailing edge of the rotor (−4.24e+003 to 2.78e0+003Pa).With the gradual increment in the RPM obtained the required velocity of 8 m/s with respect to the increments in the zone’s altitude (Zone-1 to Zone-5). The influence of the pressure distribution from the leading edge to the tailing edge in Design-1 and Design-3 has a drastic change, and finally, the proposed Design-3 performed smooth pressure distribution profiles (Figures19–22) in all operational conditions which are simulated for Siachen Glacier geographical location.Figure 19
Pressure contour at 0°C temperature.Figure 20
Pressure contour at -25°C temperature.Figure 21
Pressure contour at -50°C temperature.Figure 22
Pressure contour at -75°C temperature.


### 6.3. Results with Respect to the Influencing Parameters

The obtained results with respect to the drone’s performance influencing the parameters of Design-1, Design-2, and Design-3 are summarized (see Equations (1)–(4) and Equation (9)) in the tables and the graphical representation.From Figure23, we can find that Design-3 and Design-2 have a similar performance at temperature -25°C. But in other temperature conditions, low performance has been identified due to its geometrical design conditions. Finally, the proposed Design-3 results the better performance in the comparison of existing Design-1 and Design-2 with respect to the all operational conditions which are simulated for Siachen Glacier geographical location.Figure 23
Pressure contour at -100°C temperature.


## 6.1. Pressure Contour for Propeller Design-1 (Zone-1 to Zone-5), Respectively

For the mimicking condition of Zone-1 with 0°C and velocity of 8 m/s obtained by 5000 RPM, the results of the pressure contour influence vary from the hub to the tailing edge of the rotor (−3.651e+003 to 8.385e+002Pa).For the mimicking condition of Zone-2 with -25°C and velocity of 8 m/s obtained by5500 RPM, the results of the pressure contour influence vary from the hub to the tailing edge of the rotor (−3.48e+003 to 9.796e0+002Pa).For the mimicking condition of Zone-3 with -50°C and velocity of 8 m/s obtained by6000 RPM, the results of the pressure contour influence vary from the hub to the tailing edge of the rotor (−4.34e+003 to 8.796e0+002Pa).For the mimicking condition of Zone-4 with -75°C and velocity of 8 m/s obtained by6500 RPM, the results of the pressure contour influence vary from the hub to the tailing edge of the rotor (−4.17e+003 to 1.129e0+002Pa).For the mimicking condition of Zone-5 with -100°C and velocity of 8 m/s obtained by7000 RPM, the results of the pressure contour influence vary from the hub to the tailing edge of the rotor (−5.33e+003 to 1.30e0+003Pa).Note: the obtained pressure contour for the propeller Design-2 (Zone-1 to Zone-5) is almost similar to Design-1. So those results are not shown in this paper.


## 6.2. Pressure Contour for Propeller Design-3 (Zone-1 to Zone-5), Respectively

For the mimicking condition of Zone-1 with -0°C and velocity of 8 m/s obtained by5000 RPM, the results of the pressure contour influence vary from the hub to the tailing edge of the rotor (−2.42e+003 to 1.434e0+003Pa).For the mimicking condition of Zone-2 with -25°C and velocity of 8 m/s obtained by5500 RPM, the results of the pressure contour influence vary from the hub to the tailing edge of the rotor (−2.74e+003 to 1.722+003Pa).For the mimicking condition of Zone-3 with -50°C and velocity of 8 m/s obtained by6000 RPM, the results of the pressure contour influence vary from the hub to the tailing edge of the rotor (−3.09e+003 to 2.03e0+003Pa).For the mimicking condition of Zone-4 with -75°C and velocity of 8 m/s obtained by6500 RPM, the results of the pressure contour influence vary from the hub to the tailing edge of the rotor (−3.65e+003 to 2.39e0+003Pa).For the mimicking condition of Zone-5 with -100°C and velocity of 8 m/s obtained by7000 RPM, the results of the pressure contour influence vary from the hub to the tailing edge of the rotor (−4.24e+003 to 2.78e0+003Pa).With the gradual increment in the RPM obtained the required velocity of 8 m/s with respect to the increments in the zone’s altitude (Zone-1 to Zone-5). The influence of the pressure distribution from the leading edge to the tailing edge in Design-1 and Design-3 has a drastic change, and finally, the proposed Design-3 performed smooth pressure distribution profiles (Figures19–22) in all operational conditions which are simulated for Siachen Glacier geographical location.Figure 19
Pressure contour at 0°C temperature.Figure 20
Pressure contour at -25°C temperature.Figure 21
Pressure contour at -50°C temperature.Figure 22
Pressure contour at -75°C temperature.


## 6.3. Results with Respect to the Influencing Parameters

The obtained results with respect to the drone’s performance influencing the parameters of Design-1, Design-2, and Design-3 are summarized (see Equations (1)–(4) and Equation (9)) in the tables and the graphical representation.From Figure23, we can find that Design-3 and Design-2 have a similar performance at temperature -25°C. But in other temperature conditions, low performance has been identified due to its geometrical design conditions. Finally, the proposed Design-3 results the better performance in the comparison of existing Design-1 and Design-2 with respect to the all operational conditions which are simulated for Siachen Glacier geographical location.Figure 23
Pressure contour at -100°C temperature.


## 7. Conclusion

The aim of the paper is to perform a computational analysis by Ansys 16.0 of newly designed rotor blades at very low temperate conditions (which can fit the atmospheric conditions of Siachen Glacier, in India) which has been achieved. From the obtained pressure contours of Design-1 (Figures19–23; see Section 6.1) and Design-3 (Figures 24–28, see Section 6.2), it can be clearly found that Design-3 is the best solution for very low-temperature conditions with high-quality pressure distribution around the propeller surface in high altitudes at low-temperature conditions. From the obtained values of the performance parameters (J, T, Q, Ct, Cq, and η from Tables 4–6, and Figures 29–31), it can also be found that the values of the thrust coefficient Ct and the efficiency η of Design-3 result in reliable and more durable values compared to Design-1 and Design-2. Table 7 and Figure 32 show the comparison of the obtained values of the Design-1, Design-2, and Design-3 efficiencies, and as a result of comparison, we can again find that Design-3 is the best design propeller for very low-temperature operational conditions such as in Siachen Glacier, in India. As a continuation of this work, a static structural analysis of Design-3 has been proposed which is going to be presented in the upcoming paper.Figure 24
Pressure contour at 0°C temperature.Figure 25
Pressure contour at -25°C temperature.Figure 26
Pressure contour at -50°C temperature.Figure 27
Pressure contour at -75°C temperature.Figure 28
Pressure contour at -100°C temperature.Figure 29
Performance influencing parameters of Design-1.Table 4
Performance influencing parameters of Design-1.
Zone-1Zone-2Zone-3Zone-4Zone-5Temperature0°C-25°C-50°C-75°C-100°CDiameter (m)0.2540.2540.2540.2540.254Velocity (m/s)88888RPM50005500600065007000J0.000406400.000369460.000338670.000312620.00029029T5.67923656.46525556.88448966.90860278.8464213Q0.06892350.07920340.104893980.105339870.13501404Ct4.55e-054.28e-053.83e-053.27e-053.61e-05Cp5.52e-072.30e-062.30e-061.97e-062.17e-06η0.005329672.05e-130.0008985620.0008985620.000768898Table 5
Performance influencing parameters of Design-2.
Zone-1Zone-2Zone-3Zone-4Zone-5Temperature0°C-25°C-50°C-75°C-100°CDiameter (m)0.2540.2540.2540.2540.254Velocity (m/s)88888RPM50005500600065007000J0.000406400.000369460.000338670.000312620.00029029T6040768324.44982355.001458765.9845727.5842365Q0.087562570.10586940.05866680.01358790.03658415Ct5013e-052.95e-052.78e-052.84e-053.10e-05Cq7.01e-077.67e-101.28e-062.53e-075.86e-07η0.044730002.57e-010.0011675080.0055676770.0002433480Table 6
Performance influencing parameters of Design-3.
Zone-1Zone-2Zone-3Zone-4Zone-5Temperature0-25-50-75-100Diameter (m)0.2540.2540.2540.2540.254Velocity (m/s)88888RPM50005500600065007000J0.000406400.000369460.000338670.000312620.00029029T2.45520963.27122764.21980215.13332246.0548878Q0.74523160.09334910.110583570.091355520.15078502Ct1.97e-052.17e-052.35e-052.43e-052.47-05Cq1.57e-105.64e-105.14e-104.54e-103.98e-10η8.0790266322.23e+1002.4624873182.6676945942.872901871Figure 30
Performance influencing parameters of Design-2.Figure 31
Performance influencing parameters of Design-3.Figure 32
Operational zones vs. propeller efficiency.Table 7
Operational zones vs. propeller efficiency.
TemperatureDesign-1Design-2Design-3Zone-100.005329630.004734628.07902663Zone-2-252.05e-132.257280042.25728004Zone-3-500.000898560.001167512.46248732Zone-4-750.000898560.005567682.66769459Zone-5-1000.00076890.002433482.87290187



---

*Source: 8843453-2021-01-06.xml*