# Impact of Modeling Simplifications on Lightning Strike Simulation for Aeroengine

Authors: Yi-fan Qian; Zhi-feng Ye; Hai-bo Zhang
Journal: Mathematical Problems in Engineering
 (2019)

---

## Abstract

With the development of electromagnetic simulation software and affordable hardware, it is allowed for us to complete simulations for EMC purposes. However, simulation demands will be immense when simulations for models with complex structures, especially aircraft components, have to be solved. Hence, it is meaningful to investigate how to minimize the computational demands. One of the solutions to reduce the simulation expense is the simplification for the simulated model. But the simplified model should be guaranteed to provide credible simulation results which do not deviate from the original model apparently. Generally, the difference between the simulation results and experimental data is estimated, or if the experimental conditions are not achieved, the comparison between the simplified model and the original one has to be analyzed, at least. This paper explores the electromagnetic simulation of a turbofan engine encountering lightning strike. With the simplifications of different components on the turbofan engine, the influences on induced currents of engine controller cables are simulated and analyzed based on the transmission-line matrix method. A combining method of components removal and geometric structure simplification is proposed to simplify the whole engine model. Simplified components include compressor, combustion chamber, turbine, and nozzle. The effects of different simplification methods are quantified, and the rationality of the simplified model is verified by simulation analysis.


---

## Body

## 1. Introduction

A commercial aircraft suffers lightning strike once a year on average in the world. The lightning impact can be divided into two categories: direct and indirect. Direct impact causes physical damage to the structure due to the direct contact of lightning. Indirect effects mainly refer to the transient changes of induced currents in cables of aircraft caused by electromagnetic effect of lightning strikes. These transient changes are caused by electromagnetic fields generated by lightning currents, which diffuse along the surface of aircraft through components such as windows or composites. In the interior of aircraft, these transients can exceed the immunity limit of the equipment, cause malfunction or damage, and endanger flight safety [1]. In order to solve the problem of lightning protection, electromagnetic simulation technology is increasingly applied to aircraft designs. Considering the lightning protection problem in the development stage of aircraft, it can reduce the cost of development and shorten the development cycle. Therefore, relevant research institutions have devoted large amount of manpower and material resources in the electromagnetic simulation technology of aircraft [2].Aeroengine is the most significant part of aircraft. Once struck by lightning, if the shielding protection is weak, it will cause a devastating impact. In recent years, in order to investigate the electromagnetic protection problems often encountered by aeroengines in flight, various algorithms have been applied to such large-scale simulation. Da Silva and Bastos used the finite element method (FEM) to simulate and analyze large-scale electronic devices and simplified the simulation model [3]. Guiffaut and Reineix applied the finite-difference time domain method (FDTD) in the simulation of electromagnetic compatibility (EMC) for a small aircraft, and they analyzed the distribution of electromagnetic field for the aircraft under high-intensity radiation [4]. High-precision simulation can provide guidance for the design, certification, and life cycle prediction of aeroengine, thus improving flight safety. At the same time, with the increasing computing power of calculating machine, it is possible to simulate the electromagnetic effect for the whole aeroengine structure. However, the electromagnetic simulation of such large complex structures is a time-consuming process, and the calculation process is arduous. For the lightning strike simulation of aeroengine, the calculation time of the simulation can be even up to several weeks. It is definitely meaningful to reduce the cost of such simulations so that the development cycle can be shortened [5]. Recently, the development of electromagnetic simulation software and affordable hardware allows scholars to complete simulations for large-scale structures. But electromagnetic simulation demands are still enormous for a model with complex structures, especially aeroengine. Hence the need to minimize the computational demands is worthy of investigation. One of the feasible methods to reduce the computational expense is the simplification for the original model. But the simplified model should be guaranteed to provide credible results of simulation. Generally, the comparison of the simulation results with experimental data is estimated. Actually, if there are no experimental data, it is necessary to analyze the difference of simulations between the simplified model and the original one. Therefore, a simplification method is urgently needed to reduce the overall cost and the required time of simulation. In the past 10 years, a few scholars began to discuss the simplification methods of such kinds of simulation, gradually. The existing simplification methods mainly focus on the details and parts of the electromagnetic simulation model. Arnaud Christophe of Fiat Automobile Companys in Brazil has roughly classified the structural components of an automobile model, the components were removed and simulated one by one, and comparison between the effects of different simplified classifications on the electromagnetic field intensity in the vehicle was obtained [6], but the model is simple so the revivification is low. Reznicek and Raida of EVEKTOR Aircraft Manufacturing Company in Czech simulated the electromagnetic field of its small-sized turboprop aircraft, taking the electromagnetic field intensity in the cabin as the simulation target, some parts of the aircraft were removed and simplified. The simulation results after removal and simplification were compared, and the rationality of geometric simplification in electromagnetic simulation was verified [7]. However, effects of induced currents in cables caused by lightning strike were not studied. Engineer Guadalupe Gutierrez et al. of Airbus simplified the three components in the structure of a turboprop loaded by A400M transporter, i.e., the air particulate separator, the oil cooler, and the electric contact. The effects of simplification on the induced currents in the internal cables of turboprop were studied in detail to prove the rationality of simplification [8]. The peak values of the induced currents were taken as the judgment in this paper with the overall error of currents.This paper is organized as follows. Section2 provides the simplification method of aeroengine, Section 3 shows the modeling process of aeroengines and cables, simulation configurations are presented in Section 4, and comparison and analysis are shown in Section 5. Finally, simulation scheme and verification are presented in Section 6. In Section 7, we draw the conclusions briefly.


## 2. Simplification Method

The primary task of geometric simplification is to define the simulation target. For an aeroengine, the largest impact of indirect lightning strike effect is that the transient current in the cable exceeds the immunity of the aeroengine controller, which leads to the malfunction and causes potential safety hazards. Secondly, the establishment of simulation model is also very significant. The higher the accuracy of structure and cable layout, the higher the simulation revivification, and the more complex the model details are, thus, more detailed grouping and simplification of components can be carried out. Although simplification is an important way to reduce simulation cost, the rationality of simplification is also needed to be guaranteed. Each simplified simulation should be compared with the simulation results of the original model, and the data should be analyzed in detail. The simplification methods of components are classified (removal, detail processing, and cannot be removed), so as to guarantee the rationality for simplification. For the electromagnetic simulation of a whole aeroengine structure, whether the model is simplified reasonably or not will greatly affect the accuracy of the final results, otherwise the simulation will be meaningless. The flow chart is shown in Figure1.Figure 1
Flow chart of simplification.


## 3. Modeling

### 3.1. Geometric Structure

A certain turbofan engine is built as the research object. The model structure includes nacelle, inlet, compressor, combustor, part of fuel lines, turbine, and nozzle. The total length is about 4 meters, as shown in Figure2.Figure 2
3D model of turbofan engine.


### 3.2. Modeling of Cables

Four cables are established in the model according to the real position of the wiring harnesses in the turbofan engine. Different areas are covered by the cables imitating the connection with few major components of the aeroengine, as shown in Figure3.Figure 3
Diagram of internal cables.In order to show the effect of various types of cables, the harnesses include four types of cables which are most commonly used in aeroengines [9], namely shielded twisted pair, twisted pair, double-shielded twisted pair, and coaxial line, respectively. Table 1 lists the details of the cables.Table 1
Information of cables.
Cables
Areas
Types
Cross section
A
Compressor
Shielded twisted pair
B
Combustor
Twisted pair
C
Turbine
Double-shielded twisted pair
D
Nozzle
Coaxial line


### 3.3. Components Grouping

The influence of induced current in the cables for different components is the focus of this paper. Because the existence of a nacelle has a great impact on the electromagnetic shielding ability of the engine, it cannot be removed, and it is not considered when grouping. Therefore, several main components of the engine are divided into nine groups, regarding the whole engine model as one group, there are ten groups in total, and their English name and abbreviation are shown in Table2.Table 2
Component grouping of aeroengine.
Groups
Abbreviations
Aeroengine model
ALL
Fan
Fan
Intermediate-pressure compressor
IPC
High-pressure compressor
HPC
Fuel nozzle
FN
Combustor
Combustor
High-pressure turbine
HPT
Intermediate-pressure turbine
IPT
Low-pressure turbine
LPT
Nozzle
Nozzle


## 3.1. Geometric Structure

A certain turbofan engine is built as the research object. The model structure includes nacelle, inlet, compressor, combustor, part of fuel lines, turbine, and nozzle. The total length is about 4 meters, as shown in Figure2.Figure 2
3D model of turbofan engine.


## 3.2. Modeling of Cables

Four cables are established in the model according to the real position of the wiring harnesses in the turbofan engine. Different areas are covered by the cables imitating the connection with few major components of the aeroengine, as shown in Figure3.Figure 3
Diagram of internal cables.In order to show the effect of various types of cables, the harnesses include four types of cables which are most commonly used in aeroengines [9], namely shielded twisted pair, twisted pair, double-shielded twisted pair, and coaxial line, respectively. Table 1 lists the details of the cables.Table 1
Information of cables.
Cables
Areas
Types
Cross section
A
Compressor
Shielded twisted pair
B
Combustor
Twisted pair
C
Turbine
Double-shielded twisted pair
D
Nozzle
Coaxial line


## 3.3. Components Grouping

The influence of induced current in the cables for different components is the focus of this paper. Because the existence of a nacelle has a great impact on the electromagnetic shielding ability of the engine, it cannot be removed, and it is not considered when grouping. Therefore, several main components of the engine are divided into nine groups, regarding the whole engine model as one group, there are ten groups in total, and their English name and abbreviation are shown in Table2.Table 2
Component grouping of aeroengine.
Groups
Abbreviations
Aeroengine model
ALL
Fan
Fan
Intermediate-pressure compressor
IPC
High-pressure compressor
HPC
Fuel nozzle
FN
Combustor
Combustor
High-pressure turbine
HPT
Intermediate-pressure turbine
IPT
Low-pressure turbine
LPT
Nozzle
Nozzle


## 4. Simulation

In this paper, CST software based on the transmission-line matrix method (TLM) is used to simulate the electromagnetic effect of lightning strike on a turbofan engine. The simulation conditions are set according to the current waveform and lightning test method defined in SAE ARP 5412 [10] and SAE ARP 5416 standards [11].


### 4.1. Transmission-Line Matrix

The transmission-line matrix method (TLM) is used to solve the three-dimensional time-domain problem of Maxwell equation, in which the transmission line is a facility capable of transmitting electromagnetic energy between two points with minimum radiation. In the TLM method, firstly, the engineering problem is simulated through the continuous transmission line network, each node in the network is represented by a lumped element, whose physical parameters correspond to the actual problem, i.e., spatial discretization. The lumped element is simulated by each node of the transmission line, i.e., time discretization. One-dimensional transmission line is taken as an example, and its equivalent circuit is shown in Figure4.Figure 4
Lumped-parameter model of transmission line.According to Kirchhoff’s law, we have(1)∂2Vz,t∂z2=lc∂2Vz,t∂t2,∂2Iz,t∂z2=lc∂2Iz,t∂t2,where V and I stand for the voltage and current in the circuit, z is the position coordinate, t is the time coordinate, and l and c represent the distributed inductance and capacitance parameters, respectively.The two-dimensional TLM is composed of parallel-connected transmission line grids. Four pulses incident to a node from four branches, scattering first, and then incident to adjacent nodes, respectively. It can be expressed as(2)Vk+1i=CVk+1r,Vk+1r=SVki,where Vi is the excitation voltage matrix, Vr is the reflection voltage matrix, S stands for the impulse scattering matrix of the node, C is the connection matrix describing the network topology, and subscripts k and k+1 represent the discrete time interval of scattering. Thus, the scattering formula can be obtained as follows:(3)Vk+1ni=12∑m=14Vkmi−Vkmi,where r stands for scattering, i stands for incident, and m and n represent the port numbers, respectively.The three-dimensional TLM is based on a symmetrical condensation node (SCN) algorithm to simulate the electromagnetic propagation in the space unit. Each three-dimensional SCN has six branches, one branch is composed of two vertical transmission lines to simulate the propagation of electromagnetic field. Because the slot model can be regarded as the dual form of a dipole antenna, according to the Babinet theorem, the radiation field vector on the slot is the same as the equivalent dipole, but the electric and magnetic field vectors are interchangeable. For the narrow slot model, the field along the long edge of the slot forms an interface through a one-dimensional transmission line connected with a symmetrical condensation node, and it synchronizes the slot with the residual pulse and the time step of the three-dimensional TLM unit through the capacitance and inductance of the slot unit length, thus the global scattering model is established [12]. The formulation is as follows:(4)YL=ΔtLsΔl,Ys=2CsΔlΔt−YL,Cs=0.637εln0.563Δlω+εdω,Ls−1=πμ2ln1.591Δl/ω−1+μωd−1.In the formula,YL is the characteristic admittance of the transmission line, YS is the capacitive characteristic admittance, ε is the relative permittivity in free space, μ is the permeability, and Δl is the distance of electromagnetic wave transmission.The TLM can deal with complex nodes, simulate wave propagation in time domain, simplify calculation, and avoid possible problems such as convergence, stability, and pseudosolution. The TLM can also get the main- and high-order mode eigenvalues through the Fourier transform, which can accurately simulate the physical characteristics and behavior of waves and can explain the propagation characteristics of different states. Therefore, it is very suitable for the simulation of three-dimensional time-domain electromagnetic field [13].


### 4.2. Excitation Source

In this paper, the expression of lightning impulse excitation waveform of airborne electronic and electrical system isIt=I0e−αt−e−βt, I0=218810A, α=11354s−1, and β=647265s−1, which is given by SAE ARP 5412 and SAE ARP 5416 standards. Lightning strike point is set at the axle tip of the turbofan engine which is most likely to be struck by lightning [14], and the exit point of lightning strike signal is at the end of the nozzle, as shown in Figure 5. The boundary condition is set to “open” to simulate the infinite free space which means electromagnetic waves will be absorbed completely without any reflection when they hit these boundaries. The waveform of the excitation signal is shown in Figure 6. The peak current is 218.81kA, peak time is 6.4μs, and duration is 500μs; frequency is set at 0–50 MHz, and there is almost no current component above 50 MHz. It can represent the most extreme lightning strike situation and save time [15].Figure 5
Lightning strike point.Figure 6
Lightning signal waveform.The schematic diagram is shown in Figure7, in which the yellow part stands for the lightning signal, the resistance between twisted pairs is 100Ω, the matching resistance of coaxial cable is 50Ω, and the shielding layers are all grounded.Figure 7
Schematic diagram.


## 4.1. Transmission-Line Matrix

The transmission-line matrix method (TLM) is used to solve the three-dimensional time-domain problem of Maxwell equation, in which the transmission line is a facility capable of transmitting electromagnetic energy between two points with minimum radiation. In the TLM method, firstly, the engineering problem is simulated through the continuous transmission line network, each node in the network is represented by a lumped element, whose physical parameters correspond to the actual problem, i.e., spatial discretization. The lumped element is simulated by each node of the transmission line, i.e., time discretization. One-dimensional transmission line is taken as an example, and its equivalent circuit is shown in Figure4.Figure 4
Lumped-parameter model of transmission line.According to Kirchhoff’s law, we have(1)∂2Vz,t∂z2=lc∂2Vz,t∂t2,∂2Iz,t∂z2=lc∂2Iz,t∂t2,where V and I stand for the voltage and current in the circuit, z is the position coordinate, t is the time coordinate, and l and c represent the distributed inductance and capacitance parameters, respectively.The two-dimensional TLM is composed of parallel-connected transmission line grids. Four pulses incident to a node from four branches, scattering first, and then incident to adjacent nodes, respectively. It can be expressed as(2)Vk+1i=CVk+1r,Vk+1r=SVki,where Vi is the excitation voltage matrix, Vr is the reflection voltage matrix, S stands for the impulse scattering matrix of the node, C is the connection matrix describing the network topology, and subscripts k and k+1 represent the discrete time interval of scattering. Thus, the scattering formula can be obtained as follows:(3)Vk+1ni=12∑m=14Vkmi−Vkmi,where r stands for scattering, i stands for incident, and m and n represent the port numbers, respectively.The three-dimensional TLM is based on a symmetrical condensation node (SCN) algorithm to simulate the electromagnetic propagation in the space unit. Each three-dimensional SCN has six branches, one branch is composed of two vertical transmission lines to simulate the propagation of electromagnetic field. Because the slot model can be regarded as the dual form of a dipole antenna, according to the Babinet theorem, the radiation field vector on the slot is the same as the equivalent dipole, but the electric and magnetic field vectors are interchangeable. For the narrow slot model, the field along the long edge of the slot forms an interface through a one-dimensional transmission line connected with a symmetrical condensation node, and it synchronizes the slot with the residual pulse and the time step of the three-dimensional TLM unit through the capacitance and inductance of the slot unit length, thus the global scattering model is established [12]. The formulation is as follows:(4)YL=ΔtLsΔl,Ys=2CsΔlΔt−YL,Cs=0.637εln0.563Δlω+εdω,Ls−1=πμ2ln1.591Δl/ω−1+μωd−1.In the formula,YL is the characteristic admittance of the transmission line, YS is the capacitive characteristic admittance, ε is the relative permittivity in free space, μ is the permeability, and Δl is the distance of electromagnetic wave transmission.The TLM can deal with complex nodes, simulate wave propagation in time domain, simplify calculation, and avoid possible problems such as convergence, stability, and pseudosolution. The TLM can also get the main- and high-order mode eigenvalues through the Fourier transform, which can accurately simulate the physical characteristics and behavior of waves and can explain the propagation characteristics of different states. Therefore, it is very suitable for the simulation of three-dimensional time-domain electromagnetic field [13].


## 4.2. Excitation Source

In this paper, the expression of lightning impulse excitation waveform of airborne electronic and electrical system isIt=I0e−αt−e−βt, I0=218810A, α=11354s−1, and β=647265s−1, which is given by SAE ARP 5412 and SAE ARP 5416 standards. Lightning strike point is set at the axle tip of the turbofan engine which is most likely to be struck by lightning [14], and the exit point of lightning strike signal is at the end of the nozzle, as shown in Figure 5. The boundary condition is set to “open” to simulate the infinite free space which means electromagnetic waves will be absorbed completely without any reflection when they hit these boundaries. The waveform of the excitation signal is shown in Figure 6. The peak current is 218.81kA, peak time is 6.4μs, and duration is 500μs; frequency is set at 0–50 MHz, and there is almost no current component above 50 MHz. It can represent the most extreme lightning strike situation and save time [15].Figure 5
Lightning strike point.Figure 6
Lightning signal waveform.The schematic diagram is shown in Figure7, in which the yellow part stands for the lightning signal, the resistance between twisted pairs is 100Ω, the matching resistance of coaxial cable is 50Ω, and the shielding layers are all grounded.Figure 7
Schematic diagram.


## 5. Simulation and Analysis

### 5.1. Simulation

Each cable is equipped with a probe to measure the induced current, according to the groups divided from the components in Section3.2, and the simulation results of the models after removing each group are obtained. Figures 8–11 show the simulation results of four cables. The abscissa is time, the unit is μs, the longitudinal coordinate is induced current, and the unit is A. As can be seen from Figures 8–11, compared to the curves of the original model, the induced current corresponding to the removal of some components has changed considerably, such as fan (Fan) and intermediate-pressure compressor (IPC). The simulation results of the induced current corresponding to some simplified models are similar to those of the original model, such as nozzle (Nozzle) and low-pressure turbine (LPT).Figure 8
Comparison of simulation results on Cable (A).Figure 9
Comparison of simulation results on Cable (B).Figure 10
Comparison of simulation results on Cable (C).Figure 11
Comparison of simulation results on Cable (D).


### 5.2. Data Analysis

In the previous section, the induced current data in each cable of the simplified model has been obtained after removing different components. In order to obtain the quantified results, the simulation results of the original model are taken as the standard values. The data of each simplified model is compared with the data of the original one. Finally, the influence of the simplified components for the induced current in cables is obtained. Because the peak value of induction current of various types of cables differs greatly, some scholars have used decibel error to measure the error in the study of induction current errors. This analysis method is worth for reference [16]. Therefore, the decibel error in electronic instrument measurement and the root mean square error (RMSE) in mathematical statistics are used. In numerical calculation, the decibel error represents the difference between the accurate value and the approximation, that is, the error between the simulation result Ji of the original model and the simulation result Ii of the simplified model. For the induced current in a single cable, the expression is as follows:(5)ξi=20lgIiJidB,i=1,2,…,n,where n stands for the sampling time points.The RMSE is used to measure the overall situation of the error, and its expression is as follows:(6)γi,j=1n∑k=1nξi,j,k2,where γi,j represents the RSME of the decibel error in the cable j with the removal of group i in the simulation, which denotes the effect of i group on the induced current in cable j, n stands for the number of samples of the induced current, and ξi,j,k represents the decibel error of the induced current at the k sampling point on cable j and the corresponding time point for the original model after removing the i group of components.In this paper, the RMSE represents the overall variation of the induced current simulation curve corresponding to the removal of one group in the model and the induced current curve of the original one. The larger the RMSE is, the greater the influence of the component on the induced current of the target cable is. The RMSE is the manifestation of the comprehensive deviation of induction current, so it is the most reliable basis for judging the influence of components on the induced current [17].The RMSE of the induced current corresponding to the simplified model can be obtained by taking the simulated data into formulation (5) and (6). As shown in Table 3, γ represents the RMSE in unit dB.Table 3
Induced current error for each simplified model.
Δ
dB
Cable (A)
Cable (B)
Cable (C)
Cable (D)
γ
Fan
12.23
1.93
2.96
1.64
γ
IPC
4.51
2.17
2.09
1.60
γ
HPC
2.21
2.23
1.81
1.76
γ
FN
0.01
0.32
0.08
0.01
γ
Combustor
0.09
∞
3.81
1.38
γ
HPT
0.11
0.35
0.79
0.20
γ
IPT
0.17
1.15
1.85
0.99
γ
LPT
0.10
1.59
4.11
0.61
γ
Nozzle
0.03
1.51
1.92
∞From Table3, it can be seen that the error of fuel nozzle and high-pressure turbine are relatively small, which indicates that the two components have less influence on the induced current in the cable, and the error of fan and combustor is relatively large, which indicates that the two components have greater influence on the induced current. Guo [18] has studied the influence of small concave-convex structure in electromagnetic simulation, and the corresponding model simplification scheme was designed. Gutierrez et al. [19] have studied the evaluation criteria of the influence on aeroengine cable induced current. Combining with the investigation results of various scholars, this paper will directly judge the evaluation for the influence of components on induced current according to the combination of investigation results by previous scholars. Combining the above data, a simplification method for the simulation model can be obtained:(1)
Groups with RMSE less than or equal to 1 dB can be removed directly.(2)
Groups with RMSE greater than 1 dB and less than 5 dB can be adjusted structurally, holes with diameter less than 20 mm can be filled, and microstructures can be smoothed.(3)
Groups with RMSE greater than 5 dB can only be adjusted slightly. The model is fine-tuned to simplify the structures such as slots, bolts, bulges, and grooves, which are less than 10 mm in size, and basically maintain the original structure.


## 5.1. Simulation

Each cable is equipped with a probe to measure the induced current, according to the groups divided from the components in Section3.2, and the simulation results of the models after removing each group are obtained. Figures 8–11 show the simulation results of four cables. The abscissa is time, the unit is μs, the longitudinal coordinate is induced current, and the unit is A. As can be seen from Figures 8–11, compared to the curves of the original model, the induced current corresponding to the removal of some components has changed considerably, such as fan (Fan) and intermediate-pressure compressor (IPC). The simulation results of the induced current corresponding to some simplified models are similar to those of the original model, such as nozzle (Nozzle) and low-pressure turbine (LPT).Figure 8
Comparison of simulation results on Cable (A).Figure 9
Comparison of simulation results on Cable (B).Figure 10
Comparison of simulation results on Cable (C).Figure 11
Comparison of simulation results on Cable (D).


## 5.2. Data Analysis

In the previous section, the induced current data in each cable of the simplified model has been obtained after removing different components. In order to obtain the quantified results, the simulation results of the original model are taken as the standard values. The data of each simplified model is compared with the data of the original one. Finally, the influence of the simplified components for the induced current in cables is obtained. Because the peak value of induction current of various types of cables differs greatly, some scholars have used decibel error to measure the error in the study of induction current errors. This analysis method is worth for reference [16]. Therefore, the decibel error in electronic instrument measurement and the root mean square error (RMSE) in mathematical statistics are used. In numerical calculation, the decibel error represents the difference between the accurate value and the approximation, that is, the error between the simulation result Ji of the original model and the simulation result Ii of the simplified model. For the induced current in a single cable, the expression is as follows:(5)ξi=20lgIiJidB,i=1,2,…,n,where n stands for the sampling time points.The RMSE is used to measure the overall situation of the error, and its expression is as follows:(6)γi,j=1n∑k=1nξi,j,k2,where γi,j represents the RSME of the decibel error in the cable j with the removal of group i in the simulation, which denotes the effect of i group on the induced current in cable j, n stands for the number of samples of the induced current, and ξi,j,k represents the decibel error of the induced current at the k sampling point on cable j and the corresponding time point for the original model after removing the i group of components.In this paper, the RMSE represents the overall variation of the induced current simulation curve corresponding to the removal of one group in the model and the induced current curve of the original one. The larger the RMSE is, the greater the influence of the component on the induced current of the target cable is. The RMSE is the manifestation of the comprehensive deviation of induction current, so it is the most reliable basis for judging the influence of components on the induced current [17].The RMSE of the induced current corresponding to the simplified model can be obtained by taking the simulated data into formulation (5) and (6). As shown in Table 3, γ represents the RMSE in unit dB.Table 3
Induced current error for each simplified model.
Δ
dB
Cable (A)
Cable (B)
Cable (C)
Cable (D)
γ
Fan
12.23
1.93
2.96
1.64
γ
IPC
4.51
2.17
2.09
1.60
γ
HPC
2.21
2.23
1.81
1.76
γ
FN
0.01
0.32
0.08
0.01
γ
Combustor
0.09
∞
3.81
1.38
γ
HPT
0.11
0.35
0.79
0.20
γ
IPT
0.17
1.15
1.85
0.99
γ
LPT
0.10
1.59
4.11
0.61
γ
Nozzle
0.03
1.51
1.92
∞From Table3, it can be seen that the error of fuel nozzle and high-pressure turbine are relatively small, which indicates that the two components have less influence on the induced current in the cable, and the error of fan and combustor is relatively large, which indicates that the two components have greater influence on the induced current. Guo [18] has studied the influence of small concave-convex structure in electromagnetic simulation, and the corresponding model simplification scheme was designed. Gutierrez et al. [19] have studied the evaluation criteria of the influence on aeroengine cable induced current. Combining with the investigation results of various scholars, this paper will directly judge the evaluation for the influence of components on induced current according to the combination of investigation results by previous scholars. Combining the above data, a simplification method for the simulation model can be obtained:(1)
Groups with RMSE less than or equal to 1 dB can be removed directly.(2)
Groups with RMSE greater than 1 dB and less than 5 dB can be adjusted structurally, holes with diameter less than 20 mm can be filled, and microstructures can be smoothed.(3)
Groups with RMSE greater than 5 dB can only be adjusted slightly. The model is fine-tuned to simplify the structures such as slots, bolts, bulges, and grooves, which are less than 10 mm in size, and basically maintain the original structure.


## 6. Simplification Scheme and Simulation

According to the simplification method and the simulation data for lightning strike simulation in Table3, the aeroengine model is processed as follows:(1)
Because the RMSE of induced current for the removal of fuel nozzle (FN) or high-pressure turbine (HPT) is less than 1 dB, it shows that the effect of these two components on induced currents in cables can be neglected, and the corresponding components can be removed completely.(2)
The RMSE of intermediate-pressure compressor (IPC), high-pressure compressor (HPC), intermediate-pressure turbine (IPT), and low-pressure turbine (LPT) are all more than 1 dB and less than 5 dB. So, the slots with diameter less than 20 mm can be filled completely, and the size of the fine structure which is smaller than the mesh size can be flattened.(3)
The RMSE of fan, combustor, and nozzle are much larger than those of other components, so their main structures are all retained, and the fine structures with sizes less than 10 mm are smoothed and the holes with diameters less than 10 mm can be filled.According to the scheme above, the whole model of the turbofan engine in Section3.2 is simplified. The induced current curves obtained by the simulations of simplified models are compared with the corresponding results of the original one. Because the maximum value of the induced currents corresponding to different cables are quite different, the simulation results of four cables are normalized, as shown in Figure 12. It can be seen that the variation trends of the curves are the same and the magnitudes are similar. The error parameters are obtained by processing the simulated data as shown in Table 4, and the RMSE are very small. This means that the induced current of the simplified model changes quite little.Figure 12
Comparison of induced currents between the simplified model and original model.Table 4
Errors of the simplified model relative to the original model.
Δ
dB
Cable (A)
Cable (B)
Cable (C)
Cable (D)
γ
Original
−
simplified
0.75
0.30
1.09
0.98However, for the verification of EMC data in modern aeronautical industry, RMSE seems to be a humble standard to measure the data. In 2006, Duffy et al. formally proposed the Feature Selective Validation (FSV) method [20], and the principle of FSV is shown in Figure 13. Subsequently, this method shows a broad application prospect in various fields of computational electromagnetics (CEM) [21]. In recent years, especially with the introduction of the FSV method into the IEEE standard [22] by the working group on computer modeling and simulation verification standards for CEM, the application and expansion of the FSV method in practical engineering have become a hot spot of widespread concern.Figure 13
Assessment flow chart.The basic principle of the FSV method is to decompose the original data into DC, low-frequency, and high-frequency, and then extract the difference metrics, i.e., amplitude difference measure (ADM) and feature difference measure (FDM), and then form a global difference measure (GDM) [23]. The piecewise approach is shown in Table 5, and the “GRADE number” and “SPREAD number” are obtained based on the ADM and FDM. The GRADE number means how many agreement categories are required (starting at excellent) before 85% of the data. 85 percent of the data fall into the top quality categories, so lower the GRADE number, better the comparison will be. The SPREAD is like a typical standard deviation because it also determines how many categories are needed to include 85% of the data, but the category is started from the highest rated category (instead of the excellent category, as in GRADE number) [24].Table 5
Piecewise approach.
FSV value
Interpretation
Assessment
X
≤
0.1
1
+
10
∗
X
Excellent
0.1
<
X
≤
0.2
2
+
10
∗
X
−
0.099
Very good
0.2
<
X
≤
0.4
3
+
5
∗
X
−
0.199
Good
0.4
<
X
≤
0.8
4
+
2.5
∗
X
−
0.399
Fair
0.8
<
X
≤
1.6
5
+
1.25
∗
X
−
0.799
Poor
X
>
1.6
6
Very poorSo, the FSV method is applied to the non-normalized data of the simplification simulation results and the original results. The histogram distribution is shown in Figure14, and the assessments are shown in Table 6, from which it can be seen that the reliability of the simplified model is more strongly illustrated by the FSV, that is to say, in this lightning strike simulation for induced currents of cables in the turbofan engine, the simplification scheme is reasonable and the original model can be replaced.Figure 14
Histogram distribution of the GDM obtained from the FSV analysis of data.Table 6
FSV results for the comparison.
FSV
GDM
GRADE
SPREAD
A
0.281
3
4
B
0.349
3
2
C
0.123
2
3
D
0.208
3
3From Table7 it can be seen that the number of meshes, memory requirement, and computational time are reduced by 53.5%, 40.0%, and 69.8%, respectively. It is obvious that the simplified model not only reduces the requirement of hardware, but also reduces the computing time.Table 7
Comparison of computing requirements.
Model
Mesh quantity
Memory requirement (G)
Computational time required (h)
Original
193904740
120
158
Simplified
90090975
72
49


## 7. Conclusion

In this paper, the lightning strike simulation task on a turbofan engine is investigated. The induced currents of the aeroengine controller cables are taken as the simulation objective, and the aeroengine model is simplified by removing components and simplifying the structures. The transmission-line matrix method (TLM) is used to simulate and analyze the effects of different components on the induced currents of cables, and the simulation results are obtained. The error of the simplified simulation scheme is less than 1.3 dB compared with the original. The computing time and memory requirement are definitely reduced.Once the simulation target is unequivocal, some components of this kind of large-scale complex device will be redundant. Therefore, the complexity of the model is proportional to the effect of simplification. Further research on simplification of EMC simulation for an aeroengine which is more complex will be needed, which will make this technique more meaningful.



---

*Source: 5176560-2019-10-10.xml*