Transactions of Nanjing University of Aeronautics & Astronautics Information for Contributors

Transactions of Nanjing University of Aeronautics&Astronautics(TNUAA)is a bimonthly journal facing international academic community.Emphasizing prompt and effective dissemination of key data and new scientific insights,TNUAA offers publication of new experimental and theoretical papers bearing on applications to all branches of aeronautics,astronautics and civil aviation.

Liquefaction susceptibility and deformation characteristics of saturated coral sandy soils subjected to cyclic loadings-a critical review

Coral sandy soils widely exist in coral island reefs and seashores in tropical and subtropical regions.Due to the unique marine depositional environment of coral sandy soils,the engineering characteristics and responses of these soils subjected to monotonic and cyclic loadings have been a subject of intense interest among the geotechnical and earthquake engineering communities.This paper critically reviews the progress of experimental investigations on the undrained behavior of coral sandy soils under monotonic and cyclic loadings over the last three decades.The focus of coverage includes the contractive-dilative behavior,the pattern of excess pore-water pressure(EPWP)generation and the liquefaction mechanism and liquefaction resistance,the small-strain shear modulus and strain-dependent shear modulus and damping,the cyclic softening feature,and the anisotropic characteristics of undrained responses of saturated coral sandy soils.In particular,the advances made in the past decades are reviewed from the following aspects:(1)the characterization of factors that impact the mechanism and patterns of EPWP build-up;(2)the identification of liquefaction triggering in terms of the apparent viscosity and the average flow coefficient;(3)the establishment of the invariable form of strain-based,stress-based,or energy-based EPWP ratio formulas and the unique relationship between the new proxy of liquefaction resistance and the number of cycles required to reach liquefaction;(4)the establishment of the invariable form of the predictive formulas of small strain modulus and strain-dependent shear modulus;and(5)the investigation on the effects of stress-induced anisotropy on liquefaction susceptibility and dynamic deformation characteristics.Insights gained through the critical review of these advances in the past decades offer a perspective for future research to further resolve the fundamental issues concerning the liquefaction mechanism and responses of coral sandy sites subjected to cyclic loadings associated with seismic events in marine environments.

Nontraditional energy-assisted mechanical machining of difficult-to-cut materials and components in aerospace community:a comparative analysis

The aerospace community widely uses difficult-to-cut materials,such as titanium alloys,high-temperature alloys,metal/ceramic/polymer matrix composites,hard and brittle materials,and geometrically complex components,such as thin-walled structures,microchannels,and complex surfaces.Mechanical machining is the main material removal process for the vast majority of aerospace components.However,many problems exist,including severe and rapid tool wear,low machining efficiency,and poor surface integrity.Nontraditional energy-assisted mechanical machining is a hybrid process that uses nontraditional energies(vibration,laser,electricity,etc)to improve the machinability of local materials and decrease the burden of mechanical machining.This provides a feasible and promising method to improve the material removal rate and surface quality,reduce process forces,and prolong tool life.However,systematic reviews of this technology are lacking with respect to the current research status and development direction.This paper reviews the recent progress in the nontraditional energy-assisted mechanical machining of difficult-to-cut materials and components in the aerospace community.In addition,this paper focuses on the processing principles,material responses under nontraditional energy,resultant forces and temperatures,material removal mechanisms,and applications of these processes,including vibration-,laser-,electric-,magnetic-,chemical-,advanced coolant-,and hybrid nontraditional energy-assisted mechanical machining.Finally,a comprehensive summary of the principles,advantages,and limitations of each hybrid process is provided,and future perspectives on forward design,device development,and sustainability of nontraditional energy-assisted mechanical machining processes are discussed.

Comparison of electrochemical behaviors of Ti-5Al-2Sn-4Zr-4Mo-2Cr-1Fe and Ti-6Al-4V titanium alloys in NaNO_(3) solution

The Ti-5Al-2Sn-4Zr-4Mo-2Cr-1Fe(β-CEZ)alloy is considered as a potential structural material in the aviation industry due to its outstanding strength and corrosion resistance.Electrochemical machining(ECM)is an efficient and low-cost technology for manufacturing theβ-CEZ alloy.In ECM,the machining parameter selection and tool design are based on the electrochemical dissolution behavior of the materials.In this study,the electrochemical dissolution behaviors of theβ-CEZ and Ti-6Al-4V(TC4)alloys in NaNO3solution are discussed.The open circuit potential(OCP),Tafel polarization,potentiodynamic polarization,electrochemical impedance spectroscopy(EIS),and current efficiency curves of theβ-CEZ and TC4 alloys are analyzed.The results show that,compared to the TC4 alloy,the passivation film structure is denser and the charge transfer resistance in the dissolution process is greater for theβ-CEZ alloy.Moreover,the dissolved surface morphology of the two titanium-based alloys under different current densities are analyzed.Under low current densities,theβ-CEZ alloy surface comprises dissolution pits and dissolved products,while the TC4 alloy surface comprises a porous honeycomb structure.Under high current densities,the surface waviness of both the alloys improves and the TC4 alloy surface is flatter and smoother than theβ-CEZ alloy surface.Finally,the electrochemical dissolution models ofβ-CEZ and TC4 alloys are proposed.

Triad-displaced ULAs configuration for non-circular sources with larger continuous virtual aperture and enhanced degrees of freedom

Non-uniform linear array(NULA)configurations are well renowned due to their structural ability for providing increased degrees of freedom(DOF)and wider array aperture than uniform linear arrays(ULAs).These characteristics play a significant role in improving the direction-of-arrival(DOA)estimation accuracy.However,most of the existing NULA geometries are primarily applicable to circular sources(CSs),while they limitedly improve the DOF and continuous virtual aperture for noncircular sources(NCSs).Toward this purpose,we present a triaddisplaced ULAs(Tdis-ULAs)configuration for NCS.The TdisULAs structure generally consists of three ULAs,which are appropriately placed.The proposed antenna array approach fully exploits the non-circular characteristics of the sources.Given the same number of elements,the Tdis-ULAs design achieves more DOF and larger hole-free co-array aperture than its sparse array competitors.Advantageously,the number of uniform DOF,optimal distribution of elements among the ULAs,and precise element positions are uniquely determined by the closed-form expressions.Moreover,the proposed array also produces a filled resulting co-array.Numerical simulations are conducted to show the performance advantages of the proposed Tdis-ULAs configuration over its counterpart designs.

Design and performance evaluation of a large field-of-view dual-particle time-encoded imager based on a depth-of-interaction detector

Time-encoded imaging is useful for identifying potential special nuclear materials and other radioactive sources at a distance.In this study,a large field-of-view time-encoded imager was developed for gamma-ray and neutron source hotspot imaging based on a depth-of-interaction(DOI)detector.The imager primarily consists of a DOI detector system and a rotary dual-layer cylindrical coded mask.An EJ276 plastic scintillator coupled with two SiPMs was designed as the DOI detector to increase the field of view and improve the imager performance.The difference in signal time at both ends and the log of the signal amplitude ratio were used to calculate the interaction position resolution.The position resolution of the DOI detector was calibrated using a collimated Cs-137 source,and the full width at half maximum of the reconstruction position of the Gaussian fitting curve was approximately 4.4 cm.The DOI detector can be arbitrarily divided into several units to independently reconstruct the source distribution images.The unit length was optimized via Am-Be source-location experiments.A multidetector filtering method is proposed for image denoising.This method can effectively reduce image noise caused by poor DOI detector position resolution.The vertical field of view of the imager was(-55°,55°)when the detector was placed in the center of the coded mask.A DT neutron source at 20 m standoff could be located within 2400 s with an angular resolution of 3.5°.

Assessment of Meteorological Threats to the Coordinated Search and Rescue of Unmanned/Manned Aircraft

The architecture and working principle of coordinated search and rescue system of unmanned/manned aircraft,which is composed of manned/unmanned aircraft and manned aircraft,were first introduced,and they can cooperate with each other to complete a search and rescue task.Secondly,a threat assessment method based on meteorological data was proposed,and potential meteorological threats,such as storms and rainfall,can be predicted by collecting and analyzing meteorological data.Finally,an experiment was carried out to evaluate the performance of the proposed method in different scenarios.The experimental results show that the coordinated search and rescue system of unmanned/manned aircraft can be used to effectively assess meteorological threats and provide accurate search and rescue guidance.

Improved Algorithm for Efficient Extraction of Relaxation Parameter Values from Wideband Permittivity of Baijiu

The complex permittivity of baijiu varies with frequency,and dielectric spectroscopy has been used to evaluate the quality.To simplify the analysis and reduce the number of the parameters,a dielectric relaxation model is often used to fit the permittivity data.However,existing fitting methods such as the least squares and particle swarm optimization methods are often computationally complex and require preset initial values.Therefore,a simpler calculation method of the relaxation parameters considering the geometric characteristics of the permittivity spectrum is proposed.It is based on the relationship between the Cole-Cole relaxation parameters and the Cole-Cole diagram,which is fitted by a geometric method.First,the concepts of the Cole-Cole parameters and the diagram are introduced,and then the process of obtaining the parameters from the complex permittivity measurement data is explained.Taking baijiu with 56%alcohol by volume(ABV)as an example,the fitting is better than the least squares method and similar to the particle swarm optimization.This method is then used for the parameter fitting of baijiu with ABV of 42-52%,and the average error is less than 1%,demonstrating its wider applicability.Finally,a prediction model is used for baijiu with 53%ABV,and the error is only 1.51%.Hence,the method can be applied to the measurement of ABV of baijiu.

"Three-in-One" Multi-Scale Structural Design of Carbon Fiber-Based Composites for Personal Electromagnetic Protection and Thermal Management

Wearable devices with efficient thermal management and electromagnetic interference(EMI) shielding are highly desirable for improving human comfort and safety. Herein, a multifunctional wearable carbon fibers(CF) @ polyaniline(PANI)/silver nanowires(Ag NWs) composites with a “branch-trunk” interlocked micro/nanostructure were achieved through "three-in-one" multi-scale design. The reasonable assembly of the three kinds of one-dimensional(1D) materials can fully exert their excellent properties i.e., the superior flexibility of CF, the robustness of PANI, and the splendid conductivity of Ag NWs. Consequently, the constructed flexible composite demonstrates enhanced mechanical properties with a tensile stress of 1.2 MPa, which was almost 6 times that of the original material. This is mainly attributed to the fact that the PNAI(branch) was firmly attached to the CF(trunk) through polydopamine(PDA), forming a robust interlocked structure. Meanwhile, the composite possesses excellent thermal insulation and heat preservation capacity owing to the synergistically low thermal conductivity and emissivity. More importantly, the conductive path of the composite established by the three 1D materials greatly improved its EMI shielding property and Joule heating performance at low applied voltage. This work paves the way for rational utilization of the intrinsic properties of 1D materials, as well as provides a promising strategy for designing wearable electromagnetic protection and thermal energy management devices.

Dynamic response of UHMWPE plates under combined shock and fragment loading

Ultra-high molecular weight polyethylene(UHMWPE)fiber composite has been extensively used to construct lightweight protective structures against ballistic impacts,yet little is known about its performance when subjected to combined blast and fragment impacts.Built upon a recently developed laboratory-scale experimental technique to generate simulated combined loading through the impact of a fragment-foam composite projectile launched from a light gas gun,the dynamic responses of fullyclamped UHMWPE plates subjected to combined loading were characterized experimentally,with corresponding deformation and failure modes compared with those measured with simulated blast loading alone.Subsequently,to explore the underlying physical mechanisms,three-dimensional(3D)numerical simulations with the method of finite elements(FE)were systematically carried out.Numerical predictions compared favorably well with experimental measurements,thus validating the feasibility of the established FE model.Relative to the case of blast loading alone,combined blast and fragment loading led to larger maximum deflections of clamped UHMWPE plates.The position of the FSP in the foam sabot affected significantly the performance of a UHMWPE target,either enhancing or decreasing its ballistic resistance.When the blast loading and fragment impact arrived simultaneously at the target,its ballistic resistance was superior to that achieved when subjected to fragment impact alone,and benefited from the accelerated movement of the target due to simultaneous blast loading.

GPU-accelerated three-dimensional reconstruction method of the Compton camera and its application in radionuclide imaging

A novel and fast three-dimensional reconstruction method for a Compton camera and its performance in radionuclide imaging is proposed and analyzed in this study. The conical surface sampling back-projection method with scattering angle correction(CSS-BP-SC) can quickly perform the back-projection process of the Compton cone and can be used to precompute the list-mode maximum likelihood expectation maximization(LM-MLEM). A dedicated parallel architecture was designed for the graphics processing unit acceleration of the back-projection and iteration stage of the CSS-BP-SC-based LM-MLEM. The imaging results of the two-point source Monte Carlo(MC) simulation demonstrate that by analyzing the full width at half maximum along the three coordinate axes, the CSS-BP-SC-based LM-MLEM can obtain imaging results comparable to those of the traditional reconstruction algorithm, that is, the simple back-projection-based LM-MLEM. The imaging results of the mouse phantom MC simulation and experiment demonstrate that the reconstruction results obtained by the proposed method sufficiently coincide with the set radioactivity distribution, and the speed increased by more than 664 times compared to the traditional reconstruction algorithm in the mouse phantom experiment. The proposed method will further advance the imaging applications of Compton cameras.

The growth and expansive applications of amorphous Ga_(2)O_(3)

As a promising ultra-wide bandgap semiconductor material,gallium oxide(Ga_(2)O_(3))is attracting extensive attention of researchers due to its feasible growth process,appropriate bandgap of 4.4 e V-5.3 e V allowing for deep-ultraviolet(deepUV)detection,good physical and chemical stability,high breakdown field strength and electron mobility,etc.Different from the strict processes for controllable crystalline Ga_(2)O_(3)(usually refer to as stable monoclinicβ-Ga_(2)O_(3)),amorphous Ga_(2)O_(3)(a-Ga_(2)O_(3))film can be prepared uniformly at low temperature on a large-area deposition substrate,suggesting great advantages such as low manufacturing cost and excellent flexibility,dispensing with high-temperature and high vacuum techniques.Thus,a-Ga_(2)O_(3)extremely facilitates important applications in various applied fields.Therefore,in this concise review,we summarize several major deposition methods for a-Ga_(2)O_(3)films,of which the characteristics are discussed.Additionally,potential methods to optimize the film properties are proposed by right of the inspiration from some recent studies.Subsequently,the applications of a-Ga_(2)O_(3)thin films,e.g.,in photodetectors,resistive random access memories(RRAMs)and gas sensors,are represented with a fruitful discussion of their structures and operating mechanisms.

Flexural wave bandgap properties of phononic crystal beams with interval parameters

Uncertainties are unavoidable in practical engineering,and phononic crystals are no exception.In this paper,the uncertainties are treated as the interval parameters,and an interval phononic crystal beam model is established.A perturbation-based interval finite element method(P-IFEM)and an affine-based interval finite element method(A-IFEM)are proposed to study the dynamic response of this interval phononic crystal beam,based on which an interval vibration transmission analysis can be easily implemented and the safe bandgap can be defined.Finally,two numerical examples are investigated to demonstrate the effectiveness and accuracy of the P-IFEM and A-IFEM.Results show that the safe bandgap range may even decrease by 10%compared with the deterministic bandgap without considering the uncertainties.

Oxidation mechanism of high-volume fraction SiCp/Al composite under laser irradiation and subsequent machining

Conventional mechanical machining of a composite material comprising an aluminum matrix reinforced with a high volume fraction of SiC particles(hereinafter referred to as an SiCp/Al composite)faces problems such as rapid tool wear,high specific cutting force,and poor surface integrity.Instead,a promising method for solving these problems is laser-induced oxidation-assisted milling(LOAM):under laser irradiation,the local workpiece material reacts with oxygen,thus forming loose and porous oxides that are easily removed.In the present work,the oxidation mechanism of SiCp/Al irradiated by a nanosecond pulsed laser is studied to better understand the laser-induced oxidation behavior and control the characteristics of the oxides,with laser irradiation experiments performed on a 65%SiCp/Al composite with various laser parameters and auxiliary gases(oxygen,nitrogen,and argon).With increasing laser pulse energy density,both the ablated groove depth and the width of the heat-affected zone increase.When oxygen is used as the auxiliary gas,an oxide layer composed of SiO_(2)and Al2O3 forms,and CO_(2)is produced and escapes from the material,thereby forming pores in the oxides.However,when nitrogen or argon is used as the auxiliary gas,a recast layer is produced that is relatively difficult to remove.Under laser irradiation,the sputtered material reacts with oxygen to form oxides on both sides of the ablated groove,and as the laser scanning path advances,the produced oxides accumulate to form an oxide layer.LOAM and conventional milling are compared using the same milling parameters,and LOAM is found to be better for reduced milling force and tool wear and improved machined surface quality.

Selecting suitable key supplier for core components during smart complex equipment central-private enterprises collaborative development process:from two different forms of evaluation information and matching perspective

With the development of central-private enterprises integration,selecting suitable key suppliers are able to provide core components for smart complex equipment.We consider selecting suitable key suppliers from matching perspective,for it not only satisfies natural development of smart complex equipment,it is also a good implementation of equipment project in central-private enterprises integration context.In in this paper,we carry out two parts of research,one is evaluation attributes based on comprehensive analysis,and the other is matching process between key suppliers and core components based on the matching attribute.In practical analysis process,we employ comprehensive evaluated analysis methods to acquire relevant attributes for the matching process that follows.In the analysis process,we adopt entropy-maximum deviation method(MDM)-decision-making trial and evaluation laboratory(DEMATEL)-technique for order preference by similarity to an ideal solution(TOPSIS)to obtain a comprehensive analysis.The entropy-MDM is applied to get weight value,DEMATEL is utilized to obtain internal relations,and TOPSIS is adopted to get ideal evaluated solution.We consider aggregating two types of evaluation information according to similarities of smart complex equipment based on the combination between geometric mean and arithmetic mean.Moreover,based on the aforementioned attributes and generalized power Heronian mean operator,we aggregate preference information to acquire relevant satisfaction degree,then combine the constructed matching model to get suitable key supplier.Through comprehensive analysis of selecting suitable suppliers,we know that two-sided matching and information aggregation can provide more research perspectives for smart complex equipment.Through analysis for relevant factors,we find that leading role and service level are also significant for the smart complex equipment development process.

Seismic response characteristics of nuclear island structure at generic soil and rock sites

The number of traditionally excellent coastal lithologic nuclear power plants is limited.It is a trend to develop nuclear power plants on soil sites in inland areas.Therefore,the seismic safety and adaptability of non-rock nuclear power plant(NPP)sites are the key concerns of nuclear safety researchers.Although the five site categories are clearly defined in the AP1000 design control documents,the effects of nuclear power five site conditions and soil nonlinearity on the seismic response characteristics of nuclear island buildings have not been systematically considered in previous related studies.In this study,targeting the AP1000 nuclear island structure as the research object,three-dimensional finite element models of a nuclear island structure at five types of sites(firm rock site(FR),soft rock site(SR),soft-to-medium soil site(SMS),upper bound soft-to-medium site(SMS-UB),and soft soil site(SS))are established.The partitioned analysis method of soil-structure interaction(PASSI)in the time-domain is used to investigate the effects of site hardness and nonlinearity on the acceleration,displacement,and acceleration response spectrum of the nuclear island structure under seismic excitation.The incremental equilibrium equation and explicit decoupling method are used to analyze the soil nonlinearity described by the Davidenkov model with simplified loading-reloading rules.The results show that,in the linear case,with the increase of site hardness,the peak ground acceleration(PGA)and the peak of acceleration response spectrum of the nuclear island structure increase except for the FR site,while the maximum displacement decreases.In nonlinear analysis,as the site hardness increases,the PGA,maximum displacement,and the peak of acceleration response spectrum of the nuclear island structure increase.The peak value of the acceleration response spectrum in the nonlinear case is greater than that in the linear case for FR,while smaller for SR and soil sites.The site nonlinearity reduces the peak values of the response spectrum for SR and soil sites much more as the site hardness decreases.The results of this study can provide a reference for design of nuclear island structures on soil and rock sites.

Design and analysis of longitudinal–flexural hybrid transducer for ultrasonic peen forming

Ultrasonic peen forming(UPF)is an emerging technology that exhibits great superiority in both its flexible operating modes and the deep residual stress that it produces compared with conventional plastic forming methods.Although ultrasonic transducers with longitudinal vibration have been widely studied,they have seldom been incorporated into UPF devices for machining in confined spaces.To meet the requirements of this type of machining,a sandwich-type piezoelectric transducer with coupled longitudinal-flexural vibrational modes is proposed.The basic structure of the transducer is designed to obtain large vibrational amplitudes in both modes.Experimental results obtained with a prototype device demonstrate the feasibility of the proposed transducer.The measured vibrational amplitude for the working face in the longitudinal vibrational mode is 1.0μm,and electrical matching increases this amplitude by 40%.The flexural vibration characteristics of the same prototype transducer are also tested and are found to be slightly smaller than those of longitudinal mode.The resultant working strokes of the UPF impact pins reach 1.7 mm and 1.2 mm in the longitudinal and flexural modes,respectively.The forming capability of the prototype has been evaluated via 15-min machining on standard 2024-T351 aluminum plates.After UPF,an improved surface morphology with lower surface roughness is obtained.The aluminum plate test piece has an apparent upper deformation with an arc height of 0.64 mm.The measured peak value of the compressive residual stress is around 250 MPa,appearing at a depth of 100μm.The proposed longitudinal-flexural hybrid transducer thus provides a high-performance tool for plate peen forming in confined spaces.

Peridynamic Study on Fracture Mode and Crack Propagation Path of a Plate with Multiple Cracks Subjected to Uniaxial Tension

How to simulate fracture mode and crack propagation path in a plate with multiple cracks is an attractive but difficult issue in fracture mechanics.Peridynamics is a recently developed nonlocal continuum formulation that can spontaneously predict the crack nucleation,branch and propagation in materials and structures through a meshfree discrete technique.In this paper,the peridynamic motion equation with boundary traction is improved by simplifying the boundary transfer functions.We calculate the critical cracking load and the fracture angles of the plate with multiple cracks under uniaxial tension.The results are consistent with those predicted by classical fracture mechanics.The fracture mode and crack propagation path are also determined.The calculation shows that the brittle fracture process of the plate with multiple cracks can be conveniently and correctly simulated by the peridynamic motion equation with boundary conditions.

Experimental Study and Failure Criterion Analysis of Rubber Fibre Reinforced Concrete under Biaxial Compression-Compression

In order to examine the biaxial compression-compression properties of rubber fibre reinforced concrete(RFRC),an experimental study on RFRC under different lateral compressive stresses was carried out by considering different rubber replacement rates and polypropylene fibre contents.The failure modes and mechanical property parameters of different RFRC working conditions were obtained from the experiment to explore the effects of rubber replacement rate and polypropylene fibre content on the biaxial compression-compression properties of RFRC.The following conclusions were drawn.Under the influence of lateral compressive stress,the biaxial compression-compression failure mode gradually developed from a columnar pattern to a flaky pattern,suggesting that the incorporation of rubber and polypropylene fibres into the concrete resulted in a significant change in the development of cracks.For different rubber replacement rates and polypropylene fibre contents,the vertical compressive stress exhibited the same developing trend under the influence of lateral compressive stress.Specifically,the lateral compressive stress imposed the minimum effect on the vertical compressive stress when the rubber replacement rate and polypropylene fibre content were 20%and 0.4%,respectively,and imposed the maximum effect when the rubber replacement rate and polypropylene fibre content were 20%and 0%,respectively.With the increase of rubber replacement rate,the vertical peak stress was significantly reduced,which implies that an appropriate amount of polypropylene fibres can increase the vertical peak stress to a certain extent.Then,the biaxial compression-compression mechanism of RFRC was analysed from the microscopic level by using scanning electron microscope(SEM).Meanwhile,based on Kupfer’s biaxial compression-compression failure criterion and the octahedral stress space,a biaxial compression-compression failure criterion for RFRC was proposed,which was proven to have good applicability.The research results of this study provide important theoretical basis for the engineering application and development of RFRC.

Experimental and numerical analyses of electrohydrodynamic force according to air pressure

Electrohydrodynamic(EHD) force produced by corona discharge is considered as a new thrust for solar-powered aircraft and stratosphere balloons in near space. However, its performance at low air pressures remains to be clarified.An experiment of measuring the EHD force at 0.02 atm–1.0 atm(1 atm = 1.01325×10^(5)Pa) is carried out with the wireto-cylinder geometric structure. The ion distribution is analyzed by using the drift-diffusion model with two-dimensional numerical simulation. The experimental result shows that the EHD force is not linearly related to the corona discharge current at low air pressures. Numerical simulation finds that the proportion of ions in the counter-direction electric field increases from approximately 0.36% to 30% when the pressure drops from 1.0 atm to 0.2 atm. As a result, the EHD force with a constant power supply drops faster than the previous theoretical prediction in the ground experiment, suggesting that the consideration of counter-direction EHD force is necessary to improve the ionic wind propulsion efficiency in near-space applications.