Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “comparable input ranges”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 343 records · Page 19

Full-Circle Resolver-to-Linear-Analog Converter

A circuit generates sinusoidal excitation signals for a shaft-angle resolver and, like the arctangent circuit described in the preceding article, generates an analog voltage proportional to the shaft angle. The disadvantages of the circuit described in the preceding article arise from the fact that it must be made from precise analog subcircuits, including a functional block capable of implementing some trigonometric identities; this circuitry tends to be expensive, sensitive to noise, and susceptible to errors caused by temperature-induced drifts and imprecise matching of gains and phases. These disadvantages are overcome by the design of the present circuit. The present circuit (see figure) includes an excitation circuit, which generates signals Ksin(Omega(t)) and Kcos(Omega(t)) [where K is an amplitude, Omega denotes 2(pi)x a carrier frequency (the design value of which is 10 kHz), and t denotes time]. These signals are applied to the excitation terminals of a shaft-angle resolver, causing the resolver to put out signals C sin(Omega(t)-Theta) and C cos(Omega(t)-Theta). The cosine excitation signal and the cosine resolver output signal are processed through inverting comparator circuits, which are configured to function as inverting squarers, to obtain logic-level or square-wave signals .-LL[cos(Omega(t)] and -LL[cos(Omega(t)-Theta)], respectively. These signals are fed as inputs to a block containing digital logic circuits that effectively measure the phase difference (which equals Theta between the two logic-level signals). The output of this block is a pulse-width-modulated signal, PWM(Theta), the time-averaged value of which ranges from 0 to 5 VDC as Theta ranges from .180 to +180deg. PWM(Theta) is fed to a block of amplifying and level-shifting circuitry, which converts the input PWM waveform to an output waveform that switches between precise reference voltage levels of +10 and -10 V. This waveform is processed by a two-pole, low-pass filter, which removes the carrier-frequency component. The final output signal is a DC potential, proportional to Theta that ranges continuously from -10 V at Theta = -180deg to +10 V at Theta = +180deg..

Alhorn, Dean C.↗

Predictions of the Azimuthal Variation of the Noise from Chevron Jets using an Acoustic Analogy

An existing acoustic analogy formulation is used to make predictions of the turbulent mixing noise from a series of chevron jets tested at NASA Glenn Research Center. Reynolds-averaged Navier-Stokes solutions are used as input to provide the mean flow and turbulence quantities needed to compute the Green’s function and parameterize the source model. The Green’s function for the non-axisymmetric mean flows generated by these nozzles is computed using a second-order finite volume method. A slightly modified version of an existing source model originally developed for round jets is used to assess its capability in these non-axisymmetric flows. Prediction results are compared with experimental data for a range of azimuthal orders of the nozzle geometry and flow speeds. The results show that the RANS-based acoustic analogy prediction method is able to capture the observed azimuthal variation of the sound field fairly well at relatively low jet exit velocities. At higher speeds, the prediction method still captures the observed trends in azimuthal directivity but not the absolute levels, particularly at higher frequencies. Areas for potential improvement of the method are suggested.

Jet Noise↗

Predictions of the Azimuthal Variation of the Noise from Chevron Jets using an Acoustic Analogy

An existing acoustic analogy formulation is used to make predictions of the turbulent mixing noise from a series of chevron jets tested at NASA Glenn Research Center. Reynolds-averaged Navier-Stokes solutions are used as input to provide the mean flow and turbulence quantities needed to compute the Green’s function and parameterize the source model. The Green’s function for the non-axisymmetric mean flows generated by these nozzles is computed using a second-order finite volume method. A slightly modified version of an existing source model originally developed for round jets is used to assess its capability in these non-axisymmetric flows. Prediction results are compared with experimental data for a range of azimuthal orders of the nozzle geometry and flow speeds. The results show that the RANS-based acoustic analogy prediction method is able to capture the observed azimuthal variation of the sound field fairly well at relatively low jet exit velocities. At higher speeds, the prediction method still captures the observed trends in azimuthal directivity but not the absolute levels, particularly at higher frequencies. Areas for potential improvement of the method are suggested.

Jet Noise↗

Performance of a Low-Power Cylindrical Hall Thruster

Recent mission studies have shown that a Hall thruster which operates at relatively constant thrust efficiency (45-55%) over a broad power range (300W - 3kW) is enabling for deep space science missions when compared with slate-of-the-art ion thrusters. While conventional (annular) Hall thrusters can operate at high thrust efficiency at kW power levels, it is difficult to construct one that operates over a broad power envelope down to 0 (100 W) while maintaining relatively high efficiency. In this note we report the measured performance (I(sub sp), thrust and efficiency) of a cylindrical Hall thruster operating at 0 (100 W) input power.

Polzin, Kurt A.↗

Tiltrotor Aeroacoustic Code (TRAC) Prediction Assessment and Initial Comparisons With TRAM Test Data

A prediction sensitivity assessment to inputs and blade modeling is presented for the TiltRotor Aeroacoustic Code (TRAC). For this study, the non-CFD prediction system option in TRAC is used. Here, the comprehensive rotorcraft code, CAMRAD.Mod 1, coupled with the high-resolution sectional loads code HIRES, predicts unsteady blade loads to be used in the noise prediction code WOPWOP. The sensitivity of the predicted blade motions, blade airloads, wake geometry, and acoustics is examined with respect to rotor rpm, blade twist and chord, and to blade dynamic modeling. To accomplish this assessment. an interim input-deck for the TRAM test model and an input-deck for a reference test model are utilized in both rigid and elastic modes. Both of these test models are regarded as near scale models of the V-22 proprotor (tiltrotor). With basic TRAC sensitivities established, initial TRAC predictions are compared to results of an extensive test of an isolated model proprotor. The test was that of the TiltRotor Aeroacoustic Model (TRAM) conducted in the Duits-Nederlandse Windtunnel (DNW). Predictions are compared to measured noise for the proprotor operating over an extensive range of conditions. The variation of predictions demonstrates the great care that must be taken in defining the blade motion. However, even with this variability, the predictions using the different blade modeling successfully capture (bracket) the levels and trends of the noise for conditions ranging from descent to ascent.

Burley, Casey L.↗

Tiltrotor Aeroacoustic Code (TRAC) Prediction Assessment and Initial Comparisons with Tram Test Data

A prediction sensitivity assessment to inputs and blade modeling is presented for the TiltRotor Aeroacoustic Code (TRAC). For this study, the non-CFD prediction system option in TRAC is used. Here, the comprehensive rotorcraft code, CAMRAD.Mod1, coupled with the high-resolution sectional loads code HIRES, predicts unsteady blade loads to be used in the noise prediction code WOPWOP. The sensitivity of the predicted blade motions, blade airloads, wake geometry, and acoustics is examined with respect to rotor rpm, blade twist and chord, and to blade dynamic modeling. To accomplish this assessment, an interim input-deck for the TRAM test model and an input-deck for a reference test model are utilized in both rigid and elastic modes. Both of these test models are regarded as near scale models of the V-22 proprotor (tiltrotor). With basic TRAC sensitivities established, initial TRAC predictions are compared to results of an extensive test of an isolated model proprotor. The test was that of the TiltRotor Aeroacoustic Model (TRAM) conducted in the Duits-Nederlandse Windtunnel (DNW). Predictions are compared to measured noise for the proprotor operating over an extensive range of conditions. The variation of predictions demonstrates the great care that must be taken in defining the blade motion. However, even with this variability, the predictions using the different blade modeling successfully capture (bracket) the levels and trends of the noise for conditions ranging from descent to ascent.

Burley, Casey L.↗

Flexible Generation of Kalman Filter Code

Domain-specific program synthesis can automatically generate high quality code in complex domains from succinct specifications, but the range of programs which can be generated by a given synthesis system is typically narrow. Obtaining code which falls outside this narrow scope necessitates either 1) extension of the code generator, which is usually very expensive, or 2) manual modification of the generated code, which is often difficult and which must be redone whenever changes are made to the program specification. In this paper, we describe adaptations and extensions of the AUTOFILTER Kalman filter synthesis system which greatly extend the range of programs which can be generated. Users augment the input specification with a specification of code fragments and how those fragments should interleave with or replace parts of the synthesized filter. This allows users to generate a much wider range of programs without their needing to modify the synthesis system or edit generated code. We demonstrate the usefulness of the approach by applying it to the synthesis of a complex state estimator which combines code from several Kalman filters with user-specified code. The work described in this paper allows the complex design decisions necessary for real-world applications to be reflected in the synthesized code. When executed on simulated input data, the generated state estimator was found to produce comparable estimates to those produced by a handcoded estimator

Richardson, Julian↗

Cordon screen: A cordon-based congestion pricing policy evaluation method for U.S. cities

Global trends toward urbanization will exacerbate traffic congestion, delays in economic productivity, and air pollution issues for growing cities. Traffic congestion pricing is one method available to help ameliorate these concerns. New York City is on the verge of implementing a cordon-based traffic congestion pricing policy around its central business district. For budget-constrained municipalities, evaluating implementation of such policy could be costly. This article proposes a sketch-planning methodology, called Cordon Screen, for major U.S. cities to evaluate the net income, traffic mitigation, and avoided pollution emissions from cordon-based traffic congestion pricing. This method relies on national datasets and limited user-specific data inputs, along with a range of user-selectable assumptions informed by academic literature to deliver order-of-magnitude results. The numerous limitations of this method are acceptable for preliminary policy evaluation to determine if greater financial investment to obtain more accurate results is justified. The Denver metropolitan area is used to demonstrate Cordon Screen capabilities, with mid-range assumption results suggesting the policy is most effective at generating net income and increasing vehicle speeds on major interstates. For Denver, the policy is comparably less effective at reducing air pollution and increasing speeds on minor roadways. Validation against early implementation results from the London cordon are acceptable. Still, users should discount revenue generation projections. Choice of cordon area may be the most difficult obstacle when using the Cordon Screen. With refinement, Cordon Screen could serve as a low-cost, open-source planning evaluation tool for growing and congested U.S. cities.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Comparison of predicted and measured temperatures of UH-60A helicopter transmission

The 2109-kW (2828-hp) UH-60A Black Hawk helicopter transmission was one of the transmissions used to obtain an experimental data base. Component improvements or new transmission concepts can thus be evaluated by comparison with the established data. Results of efficiency and vibration tests of the UH-60A have been reported previously. In this investigation the transmission was instrumented internally and tested over a range of operating conditions. The speed was varied from 50 to 100 percent of the full rated value, and the torque was varied from 10 to 100 percent of the full rated value. Temperatures of internal bearings and gears were measured. The computer program Planetsys was used to simulate the thermal performance of this transmission. The calculated temperatures were then compared with the corresponding measured values. The highest measured temperature was 405 K (270 F) on the outer race of the high-speed input shaft roller bearing, at the 100-percent power condition. In general, the predicted temperatures compared very well with the measured values, most of them being within 5 kelvin (9 F). Specifically, the temperatures predicted for the single-row spherical roller planetary bearing averaged only about 4 kelvin (8 F) lower than the highest measured value.

Coe, Harold H.↗

Dynamic Weakening (Extinction) of Simple Hydrocarbon-air Counterflow Diffusion Flames by Oscillatory Inflows

This study of laminar non-premixed HC-air flames used an Oscillatory-input Opposed Jet Burner (OOJB) system developed from a previously well-characterized 7.2-mm Pyrex-nozzle OJB system. Over 600 dynamic Flame Strength (FS) measurements were obtained on unanchored (free-floating) laminar Counterflow Diffusion Flames (CFDF's). Flames were stabilized using plug inflows having steady-plus-sinusoidal axial velocities of varied magnitude, frequency, f, up to 1600 Hz, and phase angle from 0 (most data) to 360 degrees. Dynamic FS is defined as the maximum average air input velocity (U(sub air), at nozzle exit) a CFDF can sustain before strain-induced extinction occurs due to prescribed oscillatory peak-to-peak velocity inputs superimposed on steady inputs. Initially, dynamic flame extinction data were obtained at low f, and were supported by 25-120 Hz Hot-Wire cold-flow velocity data at nozzle exits. Later, expanded extinction data were supported by 4-1600 Hz Probe Microphone (PM) pk/pk P data at nozzle exits. The PM data were first obtained without flows, and later with cold stagnating flows, which better represent speaker-diaphragm dynamics during runs. The PM approach enabled characterizations of Dynamic Flame Weakening (DFW) of CFDF's from 8 to 1600 Hz. DFW was defined as % decrease in FS per Pascal of pk/pk P oscillation, namely, DFW = - 100 d(U(sub air) / U(sub air),0Hz) / d(pkpk P). The linear normalization with respect to acoustic pressure magnitude (and steady state (SS) FS) led to a DFW unaffected by strong internal resonances. For the C2H4/N2-air system, from 8 to 20 Hz, DFW is constant at 8.52 plus or minus 0.20 (% weakening)/Pa. This reflects a quasi-steady flame response to an acoustically induced dU(sub air)/dP. Also, it is surprisingly independent of C2H4/N2 mole fraction due to normalization by SS FS. From 20 to approximately 150 Hz, the C2H4/N2 air-flames weakened progressively less, with an inflection at approximately 70 Hz, and became asymptotically insensitive (DFW approximately 0) at approximately 300 Hz, which continued to 1600 Hz. The DFW of CH4-air flames followed a similar pattern, but showed much greater weakening than C2H4/N2-air flames; i.e., the quasi-steady DFW (8 to approximately 15 Hz) was 44.3 %/Pa, or approximately 5x larger, even though the 0 Hz (SS) FS was only 3.0 x smaller. The quasi-steady DFW's of C3H8-air and C2H6-air were intermediate at 34.8 and 20.9 %Pa, respectively. The DFW profiles of all four fuels, at various frequencies, correlated well but non-linearly with respective SS FS's. Notably, the DFW profile for C3H8 air fell more rapidly in the range greater than 15 to 60 Hz, compared with the 1- and 2-carbon fuels. This may indicate a shift in chemical kinetics, and/or O2 transport to a flame that moved closer to the fuel-side. In conclusion, Dynamic Flame Weakening limits appear significant and unique for each fuel, and correlate closely, but non-linearly, with Steady-State Flame Strengths at any given frequency. For reasons unknown, the dynamic flames didn't weaken more at intermediate frequencies (e.g., at 20-50 Hz) than they did at low frequencies (less than 15 Hz), where quasi-steady weakening appears to dominate. Quasi-steady flame weakening ostensibly represents a transient input strain rate maximum that just exceeds the steady-state strain-rate-limited extinction limit for a few cycles. Clearly, further detailed mechanistic understanding is needed in the fall-off region.

Pellett, G.↗

Analysis and Control of a Three-Port Converter with High-Frequency AC Link for Hybrid Unmanned/Autonomous Aerial Vehicles

In this study, a three-port resonant converter with a high-frequency isolation link for hybrid unmanned/autonomous aerial vehicle (UAV/AAV) applications is analyzed. The rectified engine-generator set output supplies the first input ports while the other is connected to a battery energy storage system. Input ports consist of half-bridge high-frequency inverters, whereas output port utilizes a semi-bridgeless active rectifier. This rectifier is formed by replacing rectifier low-side diodes with active switching devices. A phase-shift duty-cycle control is implemented for switching signals to regulate the power share and flow between ports. The proposed control system implemented achieves zero voltage switching (ZVS), reduces switching losses, and minimum current flow optimization ensures regulation in a wide range of output power. Compared to the existing bidirectional pulse width modulation (PWM) converters, the proposed system decreases the number of power conversion stages and provides centralized control without needing communications between the ports. The proposed control technique operating in closed-loop control achieves the power flow between ports under various load conditions. Finally, theoretical analysis demonstrates the proposed converter’s viability and laboratory experimental validations at 500 W full power.

33 ADVANCED PROPULSION SYSTEMS↗

3D-Printed mmWave Quasi-Holographic Antenna for 2-D Beamforming

This paper presents the design and manufacturing of a novel 2D-scanning antenna that integrates a 3D-printed Rotman lens with a quasi-holographic leaky-wave antenna (HLWA). The proposed design achieved beam-scanning capabilities by leveraging the beamforming of the Rotman lens and the high-gain directive properties of the quasi-HLWA. The Rotman lens (RL) enables beam steering in the elevation plane by switching between input ports. The quasi-HLWA, designed using holographic principles, achieves frequency-controlled beam scanning in the azimuth plane. The entire antenna structure was fabricated using additive manufacturing with an Ink1092 substrate and silver ink for the conductive traces. This approach provides greater control over material placement and design freedom compared to traditional methods. A 25° transmission linear substrate taper was used to ensure good impedance matching between the Rotman lens and the quasi-HLWA, allowing greater gain while maintaining a good scanning range. The experimental results validate the 2D scanning capability of the proposed antenna. The antenna system provides coverage from −54° to 54° in the elevation θ plane and −28° to 28° in the azimuth plane ϕ , with a maximum measured gain of 21.3 dBi at 28 GHz with an average radiation efficiency η=60 %. The fabricated prototype was tested, and the performance was in good agreement with the simulated performance.

Rotman lens↗

Comparison of Radiative Forcing Calculations Due to Mineral Dust from a Transport Model, Satellite Measurements and an Assimilation System

This study uses information on mineral aerosol from a transport model to calculate global radiative forcing values. The transport model is driven by assimilated meteorology and outputs three-dimensional dust spatial information for various size ranges. The dust fields are input to an off-line radiative transfer calculation to obtain the direct radiative forcing due to the dust fields. During June, July and August of 1988 presence of dust 1) reduces the global net incoming radiation at the top of atmosphere (TOA) by 0.3 to 0.7 W/sq m and 2) reduces net incoming radiation at the earth's surface by 1.3 to 2.0 W/sq m. Over Africa our estimates of the reduction of radiation at the top of atmosphere compare well with TOA reductions derived from ERBE and TOMS satellite data. However, our heating rates are not consistent with analysis temperature increments produced by the assimilation system over regions of high aerosol loading. These increments are based on differences between temperature observations and temperatures from the assimilation general circulation model. One explanation is that the lower tropospheric temperatures retrieved by TOVS are being contaminated by mineral aerosol.

Weaver, Clark J.↗

Estimation of Aerosol Optical Depth at Different Wavelengths by Multiple Regression Method

This study aims to investigate and establish a suitable model that can help to estimate aerosol optical depth (AOD) in order to monitor aerosol variations especially during non-retrieval time. The relationship between actual ground measurements (such as air pollution index, visibility, relative humidity, temperature, and pressure) and AOD obtained with a CIMEL sun photometer was determined through a series of statistical procedures to produce an AOD prediction model with reasonable accuracy. The AOD prediction model calibrated for each wavelength has a set of coefficients. The model was validated using a set of statistical tests. The validated model was then employed to calculate AOD at different wavelengths. The results show that the proposed model successfully predicted AOD at each studied wavelength ranging from 340 nm to 1020 nm. To illustrate the application of the model, the aerosol size determined using measure AOD data for Penang was compared with that determined using the model. This was done by examining the curvature in the ln [AOD]-ln [wavelength] plot. Consistency was obtained when it was concluded that Penang was dominated by fine mode aerosol in 2012 and 2013 using both measured and predicted AOD data. These results indicate that the proposed AOD prediction model using routine measurements as input is a promising tool for the regular monitoring of aerosol variation during non-retrieval time.

Tan, Fuyi↗

Mathematical Morphological Filtering with a Self-Adaptive Reconstruction Technique and Application to Local Seismic Data

Recorded seismic data are generally contaminated by noise from different sources, which masks the signals of interest. In the seismology community, frequency filtering (FF) is the standard method for noise suppression. However, when the signal of interest and noise share the same frequency band, the latter cannot be filtered out without infringing on the former. We implemented a noise suppression approach based on the mathematical morphology theorem. The method involves compound operations of dilation and erosion using structuring elements of varying lengths and decomposes an input noisy waveform into several time functions with differing characteristics. Further, the filtered waveform is constructed from the time functions using a self-adaptive reconstruction technique. Application to a data set of >4700 local waveforms suggests that the implemented mathematical morphological filtering (MMF) approach is efficient for data with low signal-to-noise ratio (SNR) and significantly outperforms FF in that SNR range. For most of the dataset, FF, machine learning (ML) denoising, and continuous wavelet transform (CWT) thresholding result in higher SNR values compared with the MMF method. However, for ~42% of the waveforms, MMF outperforms FF, and the SNR gain achieved with MMF is as large as ~23 dB. Compared to ML denoising and CWT thresholding, this proportion drops to only ~10%–14%. Our results suggests that in an operational setting, MMF cannot replace the other noise suppression methods; however, signal detection can be improved if MMF is used to supplement them in some scenarios. MMF could help detect signals in problematic low-SNR data, which are currently being missed particularly when using FF alone.

58 GEOSCIENCES↗

Acoustic predictions using measured pressures from a model rotor in the DNW

A contemporary design, 4-bladed United Technologies model rotor with pressure-instrumented blades was tested in the Duits-Nederslandse Windtunnel. Simultaneous acoustic and pressure measurements were made for a wide range of operating conditions. Microphones were optimally positioned at a number of locations in the flow forward of the rotor to measure rotor thickness noise, high-speed impulsive noise (both in the rotor plane), and blade-vortex interaction noise (forward and 25 deg below the rotor plane). The blade surface pressure data are used as aerodynamic input to WOPWOP, which is a state-of-the-art rotor noise prediction program that predicts rotor thickness and loading noise. The predicted results using WOPWOP are compared to the measured noise levels for cases where either thickness noise, blade-vortex interaction noise, or high-speed impulsive noise is the dominant noise mechanism. The comparisons show regions of good agreement, as well as areas where further improvement is necessary.

Visintainer, Joseph A.↗

Study of the TRAC Airfoil Table Computational System

The report documents the study of the application of the TRAC airfoil table computational package (TRACFOIL) to the prediction of 2D airfoil force and moment data over a wide range of angle of attack and Mach number. The TRACFOIL generates the standard C-81 airfoil table for input into rotorcraft comprehensive codes such as CAM- RAD. The existing TRACFOIL computer package is successfully modified to run on Digital alpha workstations and on Cray-C90 supercomputers. A step-by-step instruction for using the package on both computer platforms is provided. Application of the newer version of TRACFOIL is made for two airfoil sections. The C-81 data obtained using the TRACFOIL method are compared with those of wind-tunnel data and results are presented.

Hu, Hong↗

Cabauw Experimental Results from the Project for Intercomparison of Land-Surface Parameterization Schemes

In the Project for Intercomparison of Land-Surface Parameterization Schemes phase 2a experiment, meteorological data for the year 1987 from Cabauw, the Netherlands, were used as inputs to 23 land-surface flux schemes designed for use in climate and weather models. Schemes were evaluated by comparing their outputs with long-term measurements of surface sensible heat fluxes into the atmosphere and the ground, and of upward longwave radiation and total net radiative fluxes, and also comparing them with latent heat fluxes derived from a surface energy balance. Tuning of schemes by use of the observed flux data was not permitted. On an annual basis, the predicted surface radiative temperature exhibits a range of 2 K across schemes, consistent with the range of about 10 W/m in predicted surface net radiation. Most modeled values of monthly net radiation differ from the observations by less than the estimated maximum monthly observational error (+/- 10 W/sq m). However, modeled radiative surface temperature appears to have a systematic positive bias in most schemes; this might be explained by an error in assumed emissivity and by models' neglect of canopy thermal heterogeneity. Annual means of sensible and latent heat fluxes, into which net radiation is partitioned, have ranges across schemes of 30 W/sq m and 25 W/sq m, respectively. Annual totals of evapotranspiration and runoff, into which the precipitation is partitioned, both have ranges of 315 mm. These ranges in annual heat and water fluxes were approximately halved upon exclusion of the three schemes that have no stomatal resistance under non-water-stressed conditions. Many schemes tend to underestimate latent heat flux and overestimate sensible heat flux in summer, with a reverse tendency in winter. For six schemes, root-mean-square deviations of predictions from monthly observations are less than the estimated upper bounds on observation errors (5 W/m for sensible beat flux and 10 W/m for latent heat flux). Actual runoff at the site is believed to be dominated by vertical drainage to ground water, but several schemes produced significant amounts of runoff as overland flow or interflow. There is a range across schemes of 184 mm (40% of total pore volume) in the simulated annual mean root-zone soil moisture. Unfortunately, no measurements of soil moisture were available for model evaluation. A theoretical analysis suggested that differences in boundary conditions used in various schemes are not sufficient to explain the large variance in soil moisture. However, many of the extreme values of soil moisture could be explained in terms of the particulars of experimental setup or excessive evapotranspiration.

Chen, Tian Hong↗