Total Cost of Ownership and Its Potential Consequences for the Development of the Hydrogen Fuel Cell Powered Vehicle Market in Poland

: Electromobility is a growing technology for land transport, constituting an important element of the concept of sustainable economic development. The article presents selected research results concerning one of the segments of this market-vehicles powered by hydrogen fuel cells. The subject of the research was to gain extensive knowledge on the economic factors inﬂuencing the future purchasing decisions of the demand side in relation to this category of vehicles. The research was based on a numerical experiment. For this purpose, a comparative analysis of purchase prices in relation to the TCO of the vehicle after 3–5 years of use was performed. The research included selected models that are powered by both conventional and alternative fuels. The use of this method will allow to assess the real costs associated with the hydrogen vehicle. The authors emphasize the important role of economic factors in the form of the TCO index for the development of this market. The experimental approach may be helpful in understanding the essence of economic relations that affect the development of the electro-mobility market and the market demand for hydrogen fuel cell-powered vehicles in Poland.


Introduction
Along with the deteriorating quality of life, there is a reorientation of activities towards sustainable and intelligent transport solutions [1][2][3][4][5][6][7]. One of such precedent, is a plan for electrification of means of transport through electromobility development programs, promoted by the governments of individual countries [8][9][10][11][12][13][14][15]. One of the key problems that should be solved as soon as possible is reduction of environmental pollution from the use of conventionally powered vehicles [16][17][18][19]. Introduction of Rechargeable Hybrid Cars (PHEV) was the first step towards clean transport [20,21]. Although they were intended to reduce emissions by 38-41% compared to conventional vehicles [22][23][24]. Research published by the International Clean Transport Council (ICCT) shows, that PHEV users rarely use the possibility of connecting their cars to the power grid. As the outcome, the actual level of CO2 emissions, was 2-4 times higher than the catalog data [25]. The next step in this process was the introduction to the market of EV/BEV vehicles, which are distinguished by zero emissions of carbon dioxide, hydrocarbons and nitrogen oxides.
A significant German consulting company -Berylls Strategy Advisors-published a report which indicates that the production and use of electric cars is very energy-intensive, On the other hand, the category of operating costs includes min. fuel-energy cost. Experts from the automotive market emphasize that it is challenging to estimate TCO elements, mainly due to the unknown as to fuel-energy prices in the long run. With the glance on the current market, the prices of raw materials depends on many factors beyond the control of the vehicle user [45][46][47]. Estimating the fuel-energy consumption, based on the analysis of the average mileage of the vehicle, is a much easier process. In the case of this cost category, observations of the market reality indicates that one should not rely on catalog data, but on reliable combustion tests of renowned research institutes.
The costs of using the vehicle also includes insurance costs. While the estimation of insurance costs in the first year of vehicle ownership is not a complex process, it may be much more difficult in the following years, while the estimation of insurance costs in the first year of vehicle ownership is not a complex process, it may be much more difficult in the following years.
Insurance costs are influenced not only by the history of our accidents, but also by the selected brand, model or engine version [48]. A slightly more straightforward process is to estimate the cost of insurance in the case of a leasing vehicle or long-term rental. In this case, the price of the policy is a component of the instalments paid and, as a rule, it is a fixed value usually defined for the entire duration of the contract.
Other cost categories are service and operation. Estimating the costs of compulsory inspections from a procedural point of view should not cause problems, because information on the price for the service and the frequency of inspections is available at any dealers of a given brand. However, a much more complicated process is to estimate the costs of possible repairs and the cost of spare parts.
The costs of usage also includes the costs of purchasing tires and registration. Although the purchase costs of tires are part of the service costs, nevertheless, due to their cyclical nature, it is easier to estimate their cost over a more extended period of time in a separate cost category. The same applies to vehicle registration costs.
At this stage, it is worth emphasizing that in the case of TCO, some researchers mention other costs of car ownership [49] such as parking fees, tolls for highways, car wash costs, etc. As a result, their TCO share is negligible and most market analysts do not include them in their TCO analysis.
Market observations shows that despite many advantages resulting from including the parameter in the form of total costs in the purchasing strategy, this tool has no practical application. This state of affairs is influenced by certain factors, resulting, inter alia, from the lack of proper knowledge in the use of this device, or an incorrect process approach in the area of meeting the purchasing needs. The analysis of purchasing decisions on the polish automotive market allows us to state that the most important of them includes:

1.
Lack of awareness of vehicle users about overhead costs beyond the purchase price apart from the purchase price. The vehicle purchase price parameter is still the basic selection criterion. Many potential users still decides to buy a vehicle, analyzing only the purchase price. This leads to the lack of awareness of the demand side that, in general, the potential economic benefits obtained from the low purchase price of the vehicle may be offset precisely by the additional costs associated with its use.

2.
Lack of a process approach in facing the purchasing needs on the demand side.
The new procurement process based on TCO analysis requires additional time, for example to collect relevant data and perform analyzes and forecasts. 3.
Lack of special knowledge and experience in conducting a self-contained TCO analysis. Most future vehicle owners are unaware that a complete TCO analysis requires a set of data that takes into account many cost parameters. 4.
No access to specialized expert programs. In up to date market's reality, a small group of automotive experts have access to the tools, which significantly limits the possibility of popularizing and using this instrument in optimizing purchasing processes.
The correct estimation of TCO is not only a complex process, but also requires interdisciplinary knowledge. It is worth emphasizing that this variable includes many factors that may not be fully verifiable at the initial stage of car purchase planning, and even very difficult to estimate. As in the case of other European automotive markets [50][51][52][53], the knowledge of the essence of the TCO parameter, i.e., the total cost of usage of a vehicle, is of fundamental economic base, both for private users and for owners of a commercial vehicle fleet in terms of the growing importance of vehicles such as polish electromobility market. In order to create appropriate mathematical algorithms, formulas and experimental data on reference combustion engines are necessary [54][55][56]. These data will make it possible to define the reference values of energy consumption and failure frequency of ICEV vehicles in relation to other propulsion systems [57][58][59][60][61].
This study extends the current research in the field of electromobility by showing the important role of economic factors identified as total costs of ownership in the development of the hydrogen fuel cell-powered vehicle market in Poland.
The results constitute the foundation for every region, country in the world, because they mainly refer to the nature of the timetable broken down into individual types of drive systems of motor vehicles. Thus, the potential user, operator or group of interested manufacturers can see the results of the specific proposed solutions and cost of ownership analysis. The consequences of the introduced consumer incentives may be the appropriate selection of the propulsion system for a given timetable in a given urban or non-urban agglomeration. Additionally, a strategy for planning the implementation of individual stages of hydrogen fuel implementation in a given territorial area can be adopted. Total cost of ownership with the adoption of specific conditions of use, operation and reliability phenomena will contribute to the precise determination of the strengths and weaknesses of the implementation of given drive system solutions or transformation of the transport system in a universal manner for any terrain.
To the best of our knowledge, this article is the first to combine TCO considerations with the development of the hydrogen cell vehicle market (including a comparative analysis of models on the market that run on both conventional and alternative fuels) and the first to address in the study case-domestic market. In short, this article brings new insights to the existing literature in the following areas: (I) electromobility and development strategy, (II) demand modeling for hydrogen fuel cell vehicles, passenger vehicle segment, (II) technology (fuel cells) and (IV) sustainable development and the domestic market.
The scientific novelty of the work is the expansion of the total cost of ownership model with utility values, including in particular the aspects of periodic and aperiodic faults. The introduction of these variables with the division into appropriate driving modes and types of drive system preferred by users significantly changes the general perception as to the justification for using a given solution. This model makes it possible to very precisely forecast the future economic effects of the introduced solution. The nature of faults and their frequency of occurrence is a very complex problem, especially in the assessment of different drives and different driving mode preferences in a given area. These data indicate significant differences in the results and analyzes in relation to the current models. Additionally, the introduced adjustments are based on actual user data for given groups, which significantly increases the credibility of the forecasts obtained.
The article is organized as follows. Section 2 provides a detailed description of the research approach. Section 3 describes the results of experimental studies, and contains their interpretations. Section 4 presents the final conclusions of the research, indicating their limitations, practical application and future directions of research in this field.

Methodology
The comparison of the energy consumption of HFEV vehicles, taking into account all energy losses resulting from the vehicle movement depending on the adopted characteristics of their displacement in relation to vehicles with other drives, is crucial in the analysis of the economic use of these vehicles. Analytical research should take into account the model of the location of the refueling station and the model of the power system and analyze their interaction. The adopted correlation of these parameters allows to determine the nature of the vehicle timetable, the potential distribution of the drive system load, which translates into the total efficiency of the vehicle and its effective range.
The adoption of the total cost of ownership (TCO) method allows to create a model of the costs of purchasing an HFEV and ICEV (Internal combustion engine vehicle) vehicle and the costs of services related to their operation and operation for a specific group of suppliers and covers the total cost of living related to the ownership of a given vehicle. This assessment is based on many variables that must be taken into account in order to obtain reliable results of the economic efficiency of the given drive system solutions. TCO modes are implemented using aggregated forecasts (periodic and stochastic models) and historical data to determine costs. A vehicle mileage of approximately 60,000 km was selected as the mileage variable. Consistent with the objective of operating cost modeling, the cost assumptions include the operator's critical costs for typical activities, i.e., approximated conditions resulting from experimentally accepted timetables. The distribution and nature of vehicle use is of particular importance for the mileage and intensification of its consumption and operating costs. The assessment also covers the frequency and cost of unforeseen malfunctions as well as periodic operating costs, e.g., replacement of necessary components and operating means. The purchase cost of the vehicle is provided by the manufacturer as well as the estimated fuel consumption and other handling costs. Generally, there are capital costs and costs of ownership which stands out.
A detailed breakdown of the Total Cost of Ownership (TCO) is shown in Figure 1.

Methodology
The comparison of the energy consumption of HFEV vehicles, taking into account all energy losses resulting from the vehicle movement depending on the adopted characteristics of their displacement in relation to vehicles with other drives, is crucial in the analysis of the economic use of these vehicles. Analytical research should take into account the model of the location of the refueling station and the model of the power system and analyze their interaction. The adopted correlation of these parameters allows to determine the nature of the vehicle timetable, the potential distribution of the drive system load, which translates into the total efficiency of the vehicle and its effective range.
The adoption of the total cost of ownership (TCO) method allows to create a model of the costs of purchasing an HFEV and ICEV (Internal combustion engine vehicle) vehicle and the costs of services related to their operation and operation for a specific group of suppliers and covers the total cost of living related to the ownership of a given vehicle. This assessment is based on many variables that must be taken into account in order to obtain reliable results of the economic efficiency of the given drive system solutions. TCO modes are implemented using aggregated forecasts (periodic and stochastic models) and historical data to determine costs. A vehicle mileage of approximately 60,000 km was selected as the mileage variable. Consistent with the objective of operating cost modeling, the cost assumptions include the operator's critical costs for typical activities, i.e., approximated conditions resulting from experimentally accepted timetables. The distribution and nature of vehicle use is of particular importance for the mileage and intensification of its consumption and operating costs. The assessment also covers the frequency and cost of unforeseen malfunctions as well as periodic operating costs, e.g., replacement of necessary components and operating means. The purchase cost of the vehicle is provided by the manufacturer as well as the estimated fuel consumption and other handling costs. Generally, there are capital costs and costs of ownership which stands out.
A detailed breakdown of the Total Cost of Ownership (TCO) is shown in Figure 1. Capital costs includes the cost of assets or purchases without VAT (in most cases, under applicable law, companies can claim a VAT refund on vehicles, excluding vehicles Capital costs includes the cost of assets or purchases without VAT (in most cases, under applicable law, companies can claim a VAT refund on vehicles, excluding vehicles intended for private purposes, not related to the company's operations). These costs determine the initial cost of purchasing vehicles as presented by the manufacturers for a given period of time. The statement assumes that the vehicle purchase costs constitute an one-off financial transaction without taking into account other forms of payment, e.g., loans or leasing installments, to which variable interest should be added (individual costs depending on the bank). It is assumed that in most cases total capital costs increases with inflation in such a way that their provisional value remains the same, without a discount rate. The exception to this assumption is the residual value. It depends significant on the nature of usage and the operating costs, e.g., replacement of essential consumables or aperiodic faults depending on the operating conditions. Fuel costs were calculated for each vehicle in PLN per kilometer (PLN / km) using the current prices of diesel, gasoline, LPG, H 2 and electricity. Timetables were adopted as an average fuel consumption, taking into account manufacturers data, and appropriate approximation measures were performed with a slight statistical error. All costs related to the maintenance, repair and operation of the vehicle were also taken into account. The costs of periodic and aperiodic repairs are an important parameter in the economic balance of vehicle use over a given distance and under given conditions. The latter group largely influences the overall costs of using the vehicle. It is associated with high costs of replacing components. As that, parts consumables are not expensive, e.g., lubricants, filters, fuel additives, apart from the costs of battery replacement in BEV vehicles (depending on the total capacity of the cells), so the cost of removing an aperiodic fault is very high. These are mainly unforeseen faults related to the electronic equipment of the vehicle, the drive and transmission system, and the power supply system. The typical operating costs of electric and hydrogen vehicles have been established based on experimental and statistical studies. In particular, unforeseen stochastic phenomena, i.e., faults, were taken into account. The remaining data was taken into account on the basis of offers from automotive, insurance, repair and spare parts manufacturers.
In order to determine the TCO of costs and benefits for users related to the purchase and use of ICEV, PHEV, BEV and HFEV powered vehicles, an appropriate modeling environment must be established that integrates all additive data. Some of these data are quite precisely given in the literature, they are mainly fixed costs (provided by manufacturers of vehicles), fuels, additional power sources, registration and insurance costs. One of the key factors influencing the final value of the benefits of adopting a given solution are costs related to: fuel price, fuel consumption, vehicle energy consumption (efficiency of all mechanical and electrical systems of the vehicle's equipment) and the frequency of periodic and aperiodic faults. The created mathematical model includes all components of the cost of ownership included in the literature, as well as additional costs that are often overlooked in such analyzes that concern unforeseen vehicle operating costs. They are related with the failure rate of individual vehicle design solutions. It should be noted that there is very little technical data on the cost and frequency of these faults for HFEV's in the entity's literature. Therefore, the data was obtained on the basis of an analysis of the frequency of the occurrence of particular faults on the basis of the data of repair plants that service these vehicles. The incidence of aperiodic faults was established on the basis of analytical tests and analysis of the structure of individual electrical and electronic solutions, energy storage, energy transmission and energy transfer for final mechanical work to the vehicle wheels. Such data can be obtained from a number of foreign websites of HFEV users expressing their opinion during the extended use of the vehicle. HFEV vehicles in Poland are a very minor group, and all of them will supplement the H 2 energy factor outside the country, due to the lack of refueling stations. The same applies to their periodic repairs and maintenance services. Therefore, most of the data was obtained from foreign databases of car repair shops and road tests. Based on simulation models, the baseline was enriched with projected data on the longer-term operation of these vehicles, based on the PHEV and BEV models. Searching for the optimal solution and credibility of the obtained results, requires the development of the model with many additional procedures. For the purposes of the research, data on fuel consumption was obtained for individual solutions and the same vehicle classes and their additional equipment, especially comfort and safety systems. Therefore, the main factor determining the energy consumption of a given solution is the efficiency of the energy conversion and storage system (not the aerodynamics of the vehicle). Data on fuel or electricity consumption and the potential costs of driving the specified distance for a given drive system was determined on the basis of various timetables. Data for urban and extra-urban distribution were adopted. It is difficult to define for such a wide group of vehicles the data specifying very detailed conditions of energy consumption. Therefore, the data was adopted based on the collected information from a wide group of vehicle users (200) for each drive system and a given vehicle class for two driving modes. These data constitute a reliable source and on the basis of the collected information, an approximating and prognostic model was made. The model includes mainly SUVs due to the very similar parameters of the body kit and total power of the system, which is similar to the energy consumption parameters related to the geometric structure of the vehicle and its power demand. Regardless of the driving mode, it can be expected that these parameters (for different drives), will experience a slight statistical error. To a large extent, this is also due to the fact that these vehicles generate a relatively large number of aperiodic faults in the case of ICEV and PHEV drives, which will be an additional advantage during tests. The cost of the battery is based on the current vehicle manufacturer's prices A vehicle's residual value represents its value in the used car market and is modeled on the vehicle's retail price at purchase and depreciating over the life of the vehicle. These prices were taken from market data. They are by far, the most reliable representation of this parameter, for each of the ICEV, PHEV, BEV and HCEV powertrains, the annual changes in these values will be significant. In particular, it results from the current technical condition, technological development of a given drive system and the need to perform its repairs and maintenance.
To a large extent, the residual value is also determined by the potential failure rate of a given vehicle and driveline for a specified model and brand, which is a common case of a sharp drop in price for ICEV, PHEV and BEV vehicles. On the basis of previous analyzes, it was established that the frequency of aperiodic failures is the highest for PHEV vehicles due to the dual drive train. On one hand, BEVs are less emergency, but on the other, the costs to be incurred in the case of battery replacement after a certain period of operation on average from 100,000-150,000 km are incomparable to the costs of occurrence and removal of periodic and aperiodic defects for ICEV and PHEV vehicles. However, for HFEVs, the cost of failure is quite high and is well known, they mainly result from the necessity to perform periodic maintenance activities, the nature and frequency of aperiodic faults including drive is a big unknown due to the small population of vehicles in use in the world and their relatively low mileage. Nevertheless, the basic criterion for the evaluation of this parameter was established and the data was brought into the model. The maintenance cost model for each vehicle type is based on a maintenance plan that includes periodic vehicle maintenance, 12-volt battery replacement, main battery replacement (BEV and PHEV), major mileage-related component replacement, and energy system maintenance costs (HFEV). The remaining operating costs for these body types, e.g., the cost of replacing tires after a certain mileage, depends on the body type and vehicle model, which in this case means the same parameters for the component. This procedure was performed on purpose, in order to limit the individual replacements of cyclical and non-cyclical components to those that indicate the nature of the work and design of a given drive system, drive train and system electricity management.
Insurance costs vary depending on the actual condition of the vehicle, insurance company, type of insurance and class of vehicle. This model represents the cost of insurance premiums with liability, comprehensive and collision insurance (limited to coverage in Poland). When determining the costs, the parameters influencing the insurance price, resulting from the visual condition of the vehicle and resulting defects during use, were omitted. In the case of Poland, the year, model, engine capacity, drive type and mileage counts. Insurance costs are modeled as a function of vehicle class and type. In order to model these costs within a vehicle class, the same or similar type of vehicle with a comparable purchase price and engine power was intentionally assumed with the least possible error. In Poland, the method of calculating the liability price in the case of BEV and HFEV vehicles is taken into account, prior to the engine capacity, but the engine power. For the group in question of vehicles, the engine power is close to each other regardless of the source of its supply and the method of energy conversion. When calculating the liability premium for BEV and HFEV vehicles, parameters that are equally important in the case of ICEV combustion vehicles and PHEV hybrid vehicles, such as the driver's age, driving experience, place of residence and, above all, the current discounts for damage-free driving are of great importance. Thus, there are practically no major differences in the calculation of the vehicle's compulsory insurance for different powertrains. This study estimates that insurance costs vary by 4.5 per year over the life of the vehicle. This is mainly due to the offers of insurance companies and the current market conditions.
The registration process in Poland is the same, regardless of the type of drive system. The only difference is the color of the manufactured license plates which are green for HFEV and BEV vehicles. Technical inspection is required and relevant documents required during the inspection. Therefore, it can be assumed that these costs are constant.

Analysis of the Number of Vehicles in Poland for Classic, Hybrid, Electric and Hydrogen Drives-A Comparative Trend of Demand for BEV and HFEV Vehicles
In order to determine the forecasted development trends in the demand for BEVs and HFEVs in Poland, an appropriate analysis was performed showing the number of registered vehicles in Poland, taking into account individual power sources and territorial distribution. Based on this analysis, it is possible to determine the development of a network of refueling stations and vehicle charging stations It can be assumed that the costs of electric and hydrogen fuel will drop with the growing interest in alternative energy sources for passenger cars. The availability of replenishment points will also increase, which will certainly translate into lower operating costs of these vehicles. The increase in the population of a given entity in the economic environment reduces the costs of indirect factors, including: the costs of replacing the main components (batteries), the costs of services (basic parts related to the propulsion system), transmission and control system, the system of intelligent monitoring of the state of charge of the batteries and their operating conditions, costs of spare parts and others related to the use of a specified drive. The number of vehicles in use in Poland, divided by voivodeships, is shown in

Chemical Fuel and Electricity Costs for SUVs with ICEV, PHEV, BEV and HFEV for Urban and Extra-Urban Mode
Fuel consumption for the Hyundai Tucson III for the ICEV 1.6 T-GDI177KM 130 kW with gasoline engine is about 9.4 L (average users report approximated based on 563 entries) and for the ICEV 2.0 CRDI 185 km 136 kW equipped with diesel engines stands for about 8.2 L (report average users approximated on the basis of 489 entries) in mixed mode.
The average charging cost is PLN 0.617 per 1 kWh, which includes the price of electricity and distribution services. The driving priceof a 100 km range, to a large extent depends on the nature of the engine load, dynamics and driving mode (urban or extra-urban). Electric energy required to travel a 100 km for PHEV and BEV vehicles (SUV class) ranges from 15.0-15.4 kWh/100 km-depending on the efficiency and frequency of the kinetic energy recuperation system from braking and the driving habits (no energy recovery braking in long-distance mode or increased energy consumption during acceleration in city mode and electricity consumption when stationary), 4 × 4 or single-axle drive, total vehicle weight and electronic equipment-multimedia and comfort.   The average cost of hydrogen fuel in the case of retail customers at Clean Energy Partnership stations, is approx. EUR 9.50 (approx. PLN 39.9/kg) per 1 kg of hydrogen.
There are yet no refueling stations in Poland. On average, HFEV SUVs consume 1.15 kg-H2/100 km in city mode and 0.95 kg-H2/100 km in extra-urban mode. The model assumes that the vehicle will be operated in a mixed mode in the ratio (40/60 city mode/extra-city   The average cost of hydrogen fuel in the case of retail customers at Clean Energy Partnership stations, is approx. EUR 9.50 (approx. PLN 39.9/kg) per 1 kg of hydrogen.
There are yet no refueling stations in Poland. On average, HFEV SUVs consume 1.15 kg-H2/100 km in city mode and 0.95 kg-H2/100 km in extra-urban mode. The model assumes

The Main TSO Module, Including Periodic and Aperiodic Operating Costs of Selected Drives in SUVs
In order for the TCO model to be reliable and to fully capture all phenomena, even those of a stochastic nature, it is very important to identify and calculate the necessary cost categories for a specific type of vehicle and drive systems. The cost of purchasing and owning a vehicle is associated with several different cost categories, in particular in the sub-category vehicle operating and repair costs.
The following TCO equation is a tailor-made approach for any passenger vehicle drive depending on the energy source: In this case, TCO is the total cost of ownership during given period of ownership, where Pp-x is the purchase price of a new vehicle with the selected drive (to reduce the error, identical comfort equipment is assumed in the group of tested vehicles), Rp-x is the resale price of the vehicle at the end of the ownership period, the value of changes in the nature of the set of values is not linear but depends on the type of drive system and operating costs incurred. This means that the age of the vehicle and the mileage are crucial to establish the mandatory inspection intervals, i.e., the operating distance, without incurring high periodic costs, i.e., replacing the main battery set in BEV and PHEV vehicles or replacing hydrogen tanks and energy transmission equipment in HFEV vehicles and the current state of wear of the vehicle, taking into account the design of the drive . The nature of the timetable is influential for energy consumption, therefore, coefficients have been introduced defining the nature of the timetable and the system load driving and energy, including kinetic energy recuperation systems for BEV and PHEV vehicles, which corresponds to: timetable coefficient for urban mode , timetable coefficient for extra-urban mode and the amount of energy recovered from the systems auxiliary . As the vehicle wear and tear progresses, the final price decreases in a non-linear way, individually for the vehicle, and the difference between Pp-x and Rpx is the vehicle depreciation value. Ec is the cost of energy equal to the cost of fuel or the cost of electricity per kilometer. Dv-ee is the total number of kilometers traveled during the Based on this analysis, it can be concluded that the total demand for PHEV, BEV and HFEV vehicles may increase significantly. This is evidenced by the overwhelming number of BEVs on PHEVs. It is mainly related to the lower number of operational faults up to the mileage not exceeding 150,000 km where the main power batteries are not required to be replaced, which speaks for the cost of the PHEV. Nevertheless, buyers more often declare the purchase of a BEV car, which is mainly related to it's electric range and low operating costs.

Chemical Fuel and Electricity Costs for SUVs with ICEV, PHEV, BEV and HFEV for Urban and Extra-Urban Mode
Fuel consumption for the Hyundai Tucson III for the ICEV 1.6 T-GDI177KM 130 kW with gasoline engine is about 9.4 L (average users report approximated based on 563 entries) and for the ICEV 2.0 CRDI 185 km 136 kW equipped with diesel engines stands for about 8.2 L (report average users approximated on the basis of 489 entries) in mixed mode.
The average charging cost is PLN 0.617 per 1 kWh, which includes the price of electricity and distribution services. The driving priceof a 100 km range, to a large extent depends on the nature of the engine load, dynamics and driving mode (urban or extraurban). Electric energy required to travel a 100 km for PHEV and BEV vehicles (SUV class) ranges from 15.0-15.4 kWh/100 km-depending on the efficiency and frequency of the kinetic energy recuperation system from braking and the driving habits (no energy recovery braking in long-distance mode or increased energy consumption during acceleration in city mode and electricity consumption when stationary), 4 × 4 or single-axle drive, total vehicle weight and electronic equipment-multimedia and comfort.
The average cost of hydrogen fuel in the case of retail customers at Clean Energy Partnership stations, is approx. EUR 9.50 (approx. PLN 39.9/kg) per 1 kg of hydrogen. There are yet no refueling stations in Poland. On average, HFEV SUVs consume 1.15 kg-H 2 /100 km in city mode and 0.95 kg-H 2 /100 km in extra-urban mode. The model assumes that the vehicle will be operated in a mixed mode in the ratio (40/60 city mode/extra-city mode), it is related to the mileage ratio in individual modes and not the driving time.

The Main TSO Module, Including Periodic and Aperiodic Operating Costs of Selected Drives in SUVs
In order for the TCO model to be reliable and to fully capture all phenomena, even those of a stochastic nature, it is very important to identify and calculate the necessary cost categories for a specific type of vehicle and drive systems. The cost of purchasing and owning a vehicle is associated with several different cost categories, in particular in the sub-category vehicle operating and repair costs.
The following TCO equation is a tailor-made approach for any passenger vehicle drive depending on the energy source: In this case, TCO is the total cost of ownership during given period of ownership, where P p-x is the purchase price of a new vehicle with the selected drive (to reduce the error, identical comfort equipment is assumed in the group of tested vehicles), R p-x is the resale price of the vehicle at the end of the ownership period, the value of changes in the nature of the set of values is not linear but depends on the type of drive system and operating costs incurred. This means that the age of the vehicle w t and the mileage w d are crucial to establish the mandatory inspection intervals, i.e., the operating distance, without incurring high periodic costs, i.e., replacing the main battery set in BEV and PHEV vehicles or replacing hydrogen tanks and energy transmission equipment in HFEV vehicles and the current state of wear of the vehicle, taking into account the design of the drive w z . The nature of the timetable is influential for energy consumption, therefore, coefficients have been introduced defining the nature of the timetable and the system load driving and energy, including kinetic energy recuperation systems for BEV and PHEV vehicles, which corresponds to: timetable coefficient for urban mode w d−tm , timetable coefficient for extra-urban mode w d−tmp and the amount of energy recovered from the systems auxiliary w d−rek . As the vehicle wear and tear progresses, the final price decreases in a non-linear way, individually for the vehicle, and the difference between P p-x and R p-x is the vehicle depreciation value. E c is the cost of energy equal to the cost of fuel or the cost of electricity per kilometer. D v-ee is the total number of kilometers traveled during the lifetime of the electric vehicle, D v-ce is the total number of kilometers traveled during the lifetime of the vehicle with an internal combustion engine, therefore E c (D v−ee + D v−ce ) is the total cost of the injected energy necessary to produce the work transferred to the wheels of the vehicle. Costs operational related to indirect transactions to purchase a vehicle is characterized by: m-is the monthly interest rate, P c -borrowed amount, n-number of monthly interest payments, P c -borrowed amount in the case of loans, (mP c /(1 − (1 + r) −N )) N-Pc means total interest cost. Additionally, it should be added: the cost of compulsory third party liability insurance to CIC and approximate costs related to maintenance and operation, including periodic replacement of components in accordance with the recommendations of the vehicle manufacturer P f and costs related to aperiodic faults A f . The costs of subsidies and taxes are ignored, they are not of great importance in Poland. In the current legal state, they do not apply and are identical for all vehicles. When considering the costs of Pf and A f , the analyzed range of vehicle mileage should be taken into account, eg from 0-100,000 km. The make, model and technology of the drive as well as the structure of the vehicle body also play a significant role in this system. Therefore, certain values must be entered coefficients based on actual analyzes of W Pf-x periodic and aperiodic W Ap-x faults within a given group of vehicles according to the opinion of services and users. This condition is also related to the timetable and operating conditions, including environmental factors.
The TCO model captures the direct costs related to the use and ownership of the vehicle depending on the type of a drive system. As a result, indirect costs, such as the cost of traveling to a charging or refueling site, will be neglected. However, it is possible to take them into account, provided that the location of the place of use and potential places for refueling are precisely presented. In the current state of development of the refueling network infrastructure, the results related to the costs of owning HFEV vehicles due to the long distance to refueling stations for these vehicles, especially outside Poland, would be biased in favor of other drives. It should be taken into account that this technology is still under development. In the case of BEV and PHEV vehicles, these costs may already be included in smaller comparative studies limited to this group of vehicles. This is justified because both versions of vehicles use the same charging network, so the distance necessary to top up electricity is identical and depends on the location of the vehicle and the charging network in a given area, e.g., within a given urban agglomeration. BEVs require long charging times due to the increased battery capacity, and charging is needed more often due to the shorter range compared to conventional ICEV fuel vehicles. The charging time of PHEV vehicles is significantly reduced compared to BEVs, but their range is incomparably shorter, which creates problems with the frequent need for replenishing energy. These vehicles are usually charged overnight at home or during the day at work, where the user is only required to perform the quick-connect procedure. In most cases, it is required to charge on site for 30-40 min, if provided charging station is available and fully operational, which leads to potential logistical problems. In case of ICEV vehicles, it requires approximately 5 min at the service station to fully charge.
Depreciation is often massive cost in a vehicle's TCO analysis, and therefore of great importance to new vehicle buyers. Moreover, it is a complex process determined by: individual characteristics of the vehicle, including energy consumption and comfort (color, equipment, shape, quality of multimedia equipment, comfort and safety), market perception of the brand, average fuel or electricity prices in a given year, maintenance costs resulting from periodic and aperiodic exploitation phenomena, government regulations and other benefits resulting from the general social trend, e.g., the desire to reduce harmful substances to the environment. As a consequence, the depreciation rates can vary considerably depending on the make, model and driveline, in particular in the case of ICEV and HFEV vehicles.

Algorithmic Sub-Module for Periodic and Aperiodic Phenomena
The influence of periodic and aperiodic phenomena depends on the design of the drive system, energy factors class and purpose of the vehicle are difficult to determine and require complex calculations. In principle, the majority of the data can be accepted on the basis of the information provided by the manufacturer. This mainly applies to periodic phenomena such as periodic repairs, replacement of components in a given mileage range and the current prices of consumables and original and spare parts (if such are available on the automotive market). In terms of costs and availability of spare parts on the automotive market for ICEV vehicles, they are very advantageous due to a number of production companies and fierce competition. Sometimes the cost of purchasing a spare part of similar quality to the original is several times lower. The prices of parts for PHEV, BEV and HFEV vehicles are very high due to the fact that there is a very minor group of companies producing parts for these vehicles. The number of companies undertaking the process of regenerating parts of significant value, e.g., drive system components, kinetic energy recuperation, energy storage, high-pressure tanks, components of the drive train is negligible. This is due to the low number of vehicles in the automotive market during operation and the slow pace of their growth, especially in Poland. In the case of HFEV vehicles Thus, the user of this vehicle is forced to purchase spare parts exclusively from the vehicle manufacturer and sourcing them from foreign warehouses, which is very costly and significantly increases the TCO. It is of particular importance for forecasting this state. In addition, due to the small number of vehicles of this type in use, there is little data on their potential defects that were not provided or foreseen by the manufacturer when designing the vehicles.
Therefore, special models were created to predict the reliability of ICEV, PHEV, BEV and HFEV vehicles, which allow to determine the trend in the frequency of occurrence and cost of aperiodic faults after driving significant distances, even over 200,000 km. They are based on manufacturer, service, user data and modeling of reliability systems monitoring. These methods are based on acute and fuzzy algebra, adaptive, neuro-fuzzy algebra and classical prognostic metadata. Currently, drive systems are made more and more towards unification and typification. The unification and typification tendencies include entire sets of sub-assemblies and assemblies of a complex system. The complexity and individuality of the system impose specific requirements for shaping their reliability in the design, production and operation stages, taking into account various external variables. Modeling the random process of electrical and mechanical failures solves the important the issue of replacement elements, which guarantees the maintenance of the reliability of complex vehicle drive systems and the precise determination of additional costs included in the TCO. All statistical components related to the exploitation system are subjected to the active and adaptive effect of the created model. The basis for data processing are advanced statistical procedures subject to numerical calculations. The model includes parametric and non-parametric, sharp, fuzzy and neural methods, leading to the creation of high accuracy of analyzes. The data distribution in the model that can be taken into account is: Normal, exponential, logarithmic, gamma, power, Maxwell, Weibull, Erlang, Rayleigh. These procedures are used in the statistical control and forecasting of aperiodic faults occurring in a given vehicle drive system and its other executive and measuring systems.
In all options of the main model there are random generators working in given configurations. Then, on the basis of probalistic sub-models, the simulation phase is moved to the simulation phase, where random generators, based on statistical input data, simulate the frequency and failure processes in the operation of a given vehicle up to a specific mileage. The essence of using the models in relation to other models commonly used in TCO analyzes are original algorithms and computational procedures at each level of the model structure related to the analysis and prediction of the reliability of given drive systems.
Selected fragments of the calculation model in the identification of costs resulting from aperiodic faults for random variables based on experimental and statistical data (overloaded functions used): Procedure for the linear algorithm: 1.
Identification of vehicle input parameters: where: kX, kY, kZ-model coefficients (trend algorithms: linear, quadratic, exponential, exponential, hyperbolic and power), kXt-time base discretization step for measurements and forecasting algorithm, ix-dimension of the table containing measurement data to build a forecast, brD a table with the ix dimension containing the measurement results.

2.
Forecast based on the data provided to the model: Procedure for the exponential algorithm: 1.

Identification of vehicle input parameters:
2.
Forecast based on the data presented on the model or the measurement results: Procedure for the power algorithm: 1.
Identification of vehicle input parameters: 2.
Forecast based on the data presented on the model or the measurement results: Procedure for the exponential algorithm: 1.
Identification of vehicle input parameters: 2.
Forecast based on the data presented on the model or the measurement results: In the version of the algorithm, based on experimental data, the equations from points 1 and 2 are used.
Analysis using a linear algorithm with disturbed parameters for individual components of the vehicle's equipment. Procedure: Analysis using an algorithm with the use of harmonic weight for individual components of the vehicle's equipment. Procedure: kU1 and kU2-harmonic weights used in the algorithm. Adaptive forecasting algorithm based on experimental results. Procedure: The remaining algorithms are used as needed depending on the complexity of the team and in determining its nature of the failure frequency distribution and aperiodic costs.

Results and Discussion
The technical data of the ICEV reference vehicles and analyzed PHEV, BEV and HFEV in terms of energy consumption per kilometer are presented in Table 1. Other data placed into the model presented above were withdrawn from a wide group of reports (over 5000) based on users data, vehicle manufacturers and repair service enterprises. Earlier research results and research work undertaken by other researchers indicated a greater advantage of BEV and PHEV drives over ICEVs. If many additional variables are introduced into the model, including periodic and aperiodic defects, it turns out that in the current state of the charging infrastructure and vehicle service and repair institutions, they are still uncompetitive. The introduction of these variables shows the weaknesses of these solutions and allows for a more efficient determination of what actions should be taken when implementing electromobility. In particular, the reliability of individual components for given drive systems should be assessed and the costs of their replacement reduced. This is still related to the cost of purchasing the battery base material.
The data on the failure frequency for the SUV class was also grasped from the data of other makes and models of vehicles in this class to increase the accuracy of the analyzes it was not limited to the group presented in Table 1.
In this analysis, it was assumed that their maximum power transmitted to the wheels is similar, so the energy demand is similar. The difference in the characteristics of the drive system of these vehicles consists in the differentiated efficiency of the general drive system, drive train and systems converting chemical energy into mechanical work transmitted to the wheels in HFEV and ICEV vehicles. The energy consumption-related information was obtained on the basis of reports from users and vehicle manufacturers with the introduction of statistical and forecasting algorithms based on the assumed timetables in two urban and extra-urban modes.
In Table 2, the correction coefficients of energy consumption in comparison with the reference ICEV petrol vehicle were calculated by simulation. For the model of fuel consumption from volumetric to energy units, net calorific values of 32.56 MJ/l for gasoline and 35.66 MJ/1 for diesel oil are assumed. A division into comfort classes was introduced in relation to the electrical equipment of the vehicle, which affects it's additive energy consumption. The division was made for two modes of urban and extra-urban driving, taking into account the data of users and experimental data-appropriate correction models were introduced. The total vehicle weights are identical and the drag coefficients are very similar. The total efficiency of all components was taken into account, depending on the design of the drive system and transmission and auxiliary components.
The input data for the main model, taking into account the types of users and vehicle kilometers and the correlation of driving modes, are presented in Table 3. The calculations included data on the nature of timetables based on the correlation of mileage (percentage) in urban and extra-urban mode and the annual mileage depending on the user's location, place of work, residence and potential locations of travel destinations. Input data: A1professionally active people-traveling short distances in the range from several to 100 km, A2-people traveling or commuting within an average distance of up to 100 km to the place of employment, A3-people traveling long distances with high travel frequency during the day, A4-people who work intensively, often moving between cities over very long distances. The energy consumption of the vehicle and its demand for energy differ significantly depending on the urban and extra-urban driving mode. BEV and PHEV electric vehicles and derivative powertrains are more competitive in urban driving mode, therefore different types of potential users have different energy consumption. These types are verified by the annual mileage, the percentage of urban and extra-urban mileage, and the share of electrically driven mileage. The service life of the vehicle is 15 years for all types of users. On the basis of the input data concerning the nature of the vehicle operation, the values of standard deviation and standard uncertainty for the urban mode (1) and for the extra-urban mode (2) were calculated for the given groups. They define the differentiation of the approximate value from the extreme values in the statistical data. The smaller the value of the coefficient, the more similar the data values for a given group, determining the nature of the timetable. The results show very similar timetables depending on a given type of user, which facilitates the analysis. The data is based on a group of 1000 users. Table 2. Coefficients correcting the energy consumption in comparison with ICEV vehicles (spark-ignition and compressionignition engines).   Tables 4 and 5 considers the total costs of ownership of ICEVs powered by gasoline without presuming the aperiodic phenomena and other operating circumstances that may increase operating costs and face these dependencies. It can be detect that the inclusion of this model offers vastly differences in TCO, which makes it possible to assess more reliably the actual operating costs of all types of vehicles. This model reflects the current availability status of the technical infrastructure in relation to a given drive system solution, in particular the refrigerant energy conversion system. The greatest changes in these values occur for groups A1 and A2, which proves the higher frequency of these phenomena and their removal costs in the urban mode. For the A3 and A4 flu, these costs also increase, but at a lesser level. Based on the results, it can be assumed that the impact of comfort equipment and energy consumption of individual equipment systems electrical, electronic and mechanical is also important for TCO. This is mainly due to the increased frequency of unforeseen failure of additional components and increased fuel consumption. In estimating TCO changes, in order to adopt the condition of introducing aperiodic phenomena, it is necessary to adopt very detailed data on the frequency and cost of faults occurring at a given mileage and in given operating conditions for a given vehicle, and to adopt appropriate calculation algorithms consistent with the previous arrangements for forecasting and approximating variables. Table 4. Total cost of ownership of the reference vehicle with ICEV petrol engine for 3 equipment classes depending on the types of users, excluding the sub-model related to aperiodic faults and other unforeseen circumstances of use (EUR2020/km).

Total Cost of Ownership of the Reference Vehicle SUV with an ICEV Petrol Engine-No Fault Model
Drive For other models, the total operating costs of the PHEV and BEV are significantly understated. This translates into a misinterpretation in relation to the profitability of their use.
An additional parameter TCO (DTCO) has to be entered in order to demonstrate the cost performance of alternative electric and hydrogen propulsion systems for the reference ICEV propulsion (TCO) (Tables 6 and 7) In 2020, the DTCO for the ICEV with a diesel engine is negative for groups A2, A3 and A4 due to the lower unit cost of fuel purchase, which compensates for the higher purchase costs of a diesel powered vehicle. Nevertheless, operational factors approximated over a longer period of operation, bearing in mind the sub-model of aperiodic phenomena, indicate similar DTCO values in relation to ICEV vehicles powered by gasoline. In the case of group A1, which has a lower annual mileage and a high frequency of urban driving, a negative DTCO can only be identified for diesel ICEVs. Despite the immense share of the total mileage in urban mode for group A1, but assuming a slight mileage and very favorable operating conditions for vehicles using electric drive (large share of kinetic energy recuperation from braking), regardless of the capacity of the main energy source used (MHEV, FHEV, PHEV, EREV or BEV), this does not compensate for losses related to the operation and purchase of the vehicle. This can be seen in particular when considering the sub-model of aperiodic phenomena. For frequent travelers on longer distances, groups A2, A3 and A4 using a partially propelled electric system can be cost-effective. The MHEV and FHEV cars show negative DTCO values compared to the reference ICEV petrol vehicle depending on the comfort class. For BEVs, they also show positive DTCO values regardless of the comfort class and the addition of a calculation model introducing additional aperiodic variables. After substituting this data in the model and mainly the high purchase cost of components resulting from unforeseen failures, their use will not it is already so cost effective. The higher annual mileage for mainly A2-A4 groups leads to a lower Total Cost of Ownership for all alternative drives than for A1 group, as the additional repair costs for intermittent and aperiodic faults are offset by lower energy consumption.    HFEV vehicles powered by hydrogen cells are currently the least competitive, as the consumption and fuel costs are similar to ICEVs, and the operating and repair costs are very high due to the limited availability of refueling stations, service and repair workshops and spare parts. Additional measures are also required to maintain high efficiency of the installation in order to achieve high energy efficiency, which is associated with periodic replacements and repairs of the power system. For these vehicles, they can prove very competitive in ensuring greater availability of these utilities and reducing fuel costs. And such a state may already take place in 2030, as long as the technology and infrastructure will continue to be systematically introduced for general use, while increasing the availability of vehicles with this power system. It seems, however, that a lot of time and changes in social perception of this system must pass power plant to increase the number of HFEV buyers.

Summary and Conclusions
The paper analyzes the impact of decarbonization, i.e., a systematic process of reducing carbon dioxide (CO2) emissions to the atmosphere by analyzing the impact of energy efficiency and indirect costs related to the implementation of HFEV vehicles, as an alternative to ICEV, PHEV and BEV for passenger cars.The results manifest that with the current indirect costs related to the production process, expansion of the refueling network infrastructure, the implementation and operation of HFEV vehicles is still a more expensive solution than for PHEV and BEV vehicles. It must be taken into account that even if the costs of H 2 distribution are omitted and the network of refueling stations and infrastructure is not available in larger urban agglomerations in Poland, especially in rural areas, the HFEV drive is still economically inefficient. It must be mentioned that currently there are no H 2 refueling stations in Poland, as well as service and repair workshops.
− The lack of an adequately developed network of EFEV service and repair vehicles significantly increases the costs of maintenance and repair activities, which negatively suggests the mass implementation of this type of drive system in Poland. At present, the operating costs are several times higher than in the case of ICEV, PHEV and BEV vehicles. However, they may be a good alternative to these propulsion systems in the future, but with a more advanced degree of expansion of the intermediate infrastructure. − The results also express that the continuous improvement of H 2 production systems provides increasing system adaptation flexibility, with the intensification of production with surplus RES, in line with the current H 2 production system efficiency and power system modernization methods. − The analysis display that the willingness to mass introduce HFEV vehicles in the case of passenger vehicles and public transport vehicles must be based on a sustainable system introducing infrastructure solutions for the expansion of the production system network, raw material distribution and refueling network, based on the growth of vehicles used on the Polish market. − The most important consideration for implementing these vehicles is the H 2 production economy and the energy efficiency of HFEV powertrains. If it is assumed that, irrespective of the design solution, the structure of the ICEV, PHEV, BEV and HFEV vehicle body remains very similar geometrically for a given model and brand, then the legitimacy of using a given solution depends on: Efficiency and operating costs of the energy storage system, efficiency and failure frequency of drive systems and transmission, overall efficiency of the drive unit and the purchase cost of a given fuel energy unit. − The nature of the timetable and the operating conditions of a given drive system also have a great impact on the legitimacy of using a given drive system, which largely determines the number of aperiodic faults, which have an economic share in the case of using these vehicles. − The development of H 2 fuel cell technology producing electricity for engines of passenger vehicles and heavy vehicles generates a real opportunity for the development of the energy market in Poland, but with the optimal combination of the gradual development of production and infrastructure with the costs of own contribution and potential economic benefits in a given area. It is assumed that the real development of this technology with sufficient network expansion may take place only in a dozen or so years (if it occurs). − Currently, in the production area in Poland, despite the production of over 1.0 million tons of hydrogen per year, there is no technology for producing hydrogen with the purity grade necessary for the safe and effective use of current vehicle drive systems. It is not until the end of 2021 that several refueling stations are planned, which significantly discourages potential customers from purchasing HFEV vehicles. − Only the co-optimization of many energy sectors of the economy in Poland, with reasonable planning of investments in the power system, can contribute to the effective use of synergy of these sectors and reduction of the costs of production potential in terms of reducing H 2 fuel prices. − The use of H 2 fuel in road transport to a large extent makes the fuel market independent of crude oil imports and the production of petroleum fuels, which will certainly translate into a reduction in CO2 emissions to the environment and will cause economic losses in the petroleum fuel sector. Therefore, the introduction of this technology requires a careful calculation of profits and losses in the global approach to the energy sector in Poland.