Direct numerical simulations of Mach 6 hypersonic flow over a 34°compression corner subject to steady jet are conducted.Distributions of skin friction coefficient,wall pressure,mean velocity and temperature,bound...Direct numerical simulations of Mach 6 hypersonic flow over a 34°compression corner subject to steady jet are conducted.Distributions of skin friction coefficient,wall pressure,mean velocity and temperature,boundary layer thickness and Stanton number demonstrate that the flow changes dramatically in the shock wave/turbulent boundary layer interaction area.It is found that the steady jet has no effect on suppressing flow separation unexpectedly,but increases its spatial scale instead.Instantaneous flow structures show that the turbulence amplification can be observed after the application of flow control,and abundant Görtler-like vorticities appear,but the strength of the main shock decreases.Analyzing the wall fluctuating pressure signals using weighted power spectral density,we found an interesting thing.That is,although the low-frequency oscillation phenomenon induced by separation shock is suppressed by the steady jet,wall fluctuating pressure beneath the jet shock is oscillating at a frequency lower than 0.1u∞/δref.Results of coherent and intermittency factor reveal that it is related to the backand-forth movement of the jet shock itself.展开更多
Unsteady aerodynamic characteristics at high angles of attack are of great importance to the design and development of advanced fighter aircraft, which are characterized by post-stall maneuverability with multiple Deg...Unsteady aerodynamic characteristics at high angles of attack are of great importance to the design and development of advanced fighter aircraft, which are characterized by post-stall maneuverability with multiple Degrees-of-Freedom(multi-DOF) and complex flow field structure.In this paper, a special kind of cable-driven parallel mechanism is firstly utilized as a new suspension method to conduct unsteady dynamic wind tunnel tests at high angles of attack, thereby providing experimental aerodynamic data. These tests include a wide range of multi-DOF coupled oscillatory motions with various amplitudes and frequencies. Then, for aerodynamic modeling and analysis, a novel data-driven Feature-Level Attention Recurrent neural network(FLAR) is proposed. This model incorporates a specially designed feature-level attention module that focuses on the state variables affecting the aerodynamic coefficients, thereby enhancing the physical interpretability of the aerodynamic model. Subsequently, spin maneuver simulations, using a mathematical model as the baseline, are conducted to validate the effectiveness of the FLAR. Finally, the results on wind tunnel data reveal that the FLAR accurately predicts aerodynamic coefficients, and observations through the visualization of attention scores identify the key state variables that affect the aerodynamic coefficients. It is concluded that the proposed FLAR enhances the interpretability of the aerodynamic model while achieving good prediction accuracy and generalization capability for multi-DOF coupling motion at high angles of attack.展开更多
To assess the aerodynamic performance and vibration characteristics of rotor blades during rotation,a study of unsteady blade surface forces is conducted in a low-speed axial flow compressor under a rotating coordinat...To assess the aerodynamic performance and vibration characteristics of rotor blades during rotation,a study of unsteady blade surface forces is conducted in a low-speed axial flow compressor under a rotating coordinate system.The capture,modulation,and acquisition of unsteady blade surface forces are achieved by using pressure sensors and strain gauges attached to the rotor blades,in conjunction with a wireless telemetry system.Based on the measurement reliability verification,this approach allows for the determination of the static pressure distribution on rotor blade surfaces,enabling the quantitative description of loadability at different spanwise positions along the blade chord.Effects caused by the factors such as Tip Leakage Flow(TLF)and flow separation can be perceived and reflected in the trends of static pressure on the blade surfaces.Simultaneously,the dynamic characteristics of unsteady pressure and stress on the blade surfaces are analyzed.The results indicate that only the pressure signals measured at the mid-chord of the blade tip can distinctly detect the unsteady frequency of TLF due to the oscillation of the low-pressure spot on the pressure surface.Subsequently,with the help of one-dimensional continuous wavelet analysis method,it can be inferred that as the compressor enters stall,the sensors are capable of capturing stall cell frequency under a rotating coordinate system.Furthermore,the stress at the blade root is higher than that at the blade tip,and the frequency band of the vibration can also be measured by the pressure sensors fixed on the casing wall in a stationary frame.While the compressor stalls,the stress at the blade root can be higher,which can provide valuable guidance for monitoring the lifecycle of compressor blades.展开更多
Propeller design is a highly intricate and interdisciplinary task that necessitates careful trade-offs between radiated noise levels and aerodynamic efficiency.To achieve efficient trade-off designs,an enhanced on-the...Propeller design is a highly intricate and interdisciplinary task that necessitates careful trade-offs between radiated noise levels and aerodynamic efficiency.To achieve efficient trade-off designs,an enhanced on-the-fly unsteady adjoint-based aerodynamic and aeroacoustic optimization methodology is developed,which maintains the fidelity of the Navier-Stokes solution for unsteady flow and of the moving-medium Ffowcs Williams-Hawkings(FW-H)formulation for capturing tonal noise.Furthermore,this on-the-fly approach enables a unified architecture for discreteadjoint sensitivity analysis encompassing both aerodynamics and aeroacoustics,facilitating effective multi-objective weighted optimizations.Subsequently,this proposed methodology is applied to perform trade-off optimizations between aerodynamics and aeroacoustics for a propeller by employing varying weighting factors to comprehend their influence on optimal configurations.The results demonstrate a positive correlation between efficiency and noise sensitivities,and thus indicate an inherent synchronicity where pursing noise reduction through purely aeroacoustic optimization inevitably entails sacrificing aerodynamic efficiency.However,by effectively incorporating appropriate weighting factors(recommended to range from 0.25 to 0.5)into the multi-objective function combined with both aerodynamics and aeroacoustics,it becomes feasible to achieve efficiency enhancement and noise reduction simultaneously.Key findings show that reducing blade planform size and equipping“rotated-S”shaped airfoil profiles in the tip region can effectively restrain noise levels while maintaining aerodynamic performance.展开更多
Obtaining unsteady hydrodynamic performance is of great significance for seaplane design.Common methods for obtaining unsteady hydrodynamic performance data include tank test and Computational Fluid Dynamics(CFD)numer...Obtaining unsteady hydrodynamic performance is of great significance for seaplane design.Common methods for obtaining unsteady hydrodynamic performance data include tank test and Computational Fluid Dynamics(CFD)numerical simulation,which are costly and time-consuming.Therefore,it is necessary to obtain unsteady hydrodynamic performance in a low-cost and high-precision manner.Due to the strong nonlinearity,complex data distribution,and temporal characteristics of unsteady hydrodynamic performance,the prediction of it is challenging.This paper proposes a Temporal Convolutional Diffusion Model(TCDM)for predicting the unsteady hydrodynamic performance of seaplanes given design parameters.Under the framework of a classifier-free guided diffusion model,TCDM learns the distribution patterns of unsteady hydrodynamic performance data with the designed denoising module based on temporal convolutional network and captures the temporal features of unsteady hydrodynamic performance data.Using CFD simulation data,the proposed method is compared with the alternative methods to demonstrate its accuracy and generalization.This paper provides a method that enables the rapid and accurate prediction of unsteady hydrodynamic performance data,expecting to shorten the design cycle of seaplanes.展开更多
Efficient and accurate simulation of unsteady flow presents a significant challenge that needs to be overcome in computational fluid dynamics.Temporal discretization method plays a crucial role in the simulation of un...Efficient and accurate simulation of unsteady flow presents a significant challenge that needs to be overcome in computational fluid dynamics.Temporal discretization method plays a crucial role in the simulation of unsteady flows.To enhance computational efficiency,we propose the Implicit-Explicit Two-Step Runge-Kutta(IMEX-TSRK)time-stepping discretization methods for unsteady flows,and develop a novel adaptive algorithm that correctly partitions spatial regions to apply implicit or explicit methods.The novel adaptive IMEX-TSRK schemes effectively handle the numerical stiffness of the small grid size and improve computational efficiency.Compared to implicit and explicit Runge-Kutta(RK)schemes,the IMEX-TSRK methods achieve the same order of accuracy with fewer first derivative calculations.Numerical case tests demonstrate that the IMEX-TSRK methods maintain numerical stability while enhancing computational efficiency.Specifically,in high Reynolds number flows,the computational efficiency of the IMEX-TSRK methods surpasses that of explicit RK schemes by more than one order of magnitude,and that of implicit RK schemes several times over.展开更多
This study presents a two-dimensional,transient model to simulate the flow and thermal behavior of CO_(2) within a fracturing wellbore.The model accounts for high-velocity flow within the tubing and radial heat exchan...This study presents a two-dimensional,transient model to simulate the flow and thermal behavior of CO_(2) within a fracturing wellbore.The model accounts for high-velocity flow within the tubing and radial heat exchange between the wellbore and surrounding formation.It captures the temporal evolution of temperature,pressure,flow velocity,and fluid density,enabling detailed analysis of phase transitions along different tubing sections.The influence of key operational and geological parameters,including wellhead pressure,injection velocity,inlet temperature,and formation temperature gradient,on the wellbore’s thermal and pressure fields is systematically investigated.Results indicate that due to intense convective transport by the high-speed CO_(2) flow,the temperature and velocity within the tubing are primarily governed by the inlet temperature and injection velocity,with relatively minor influence from radial heat transfer with the formation.The pressure,flow velocity,and density of CO_(2) within the tubing are strongly dependent on wellhead conditions.Frictional losses and well depth contribute to pressure variations,particularly in the horizontal section of the wellbore,where a noticeable pressurization effect increases the fluid density.During injection,liquid CO_(2) initially undergoes a rapid transition to a supercritical state,with the depth at which this phase change occurs stabilizing as injection progresses.展开更多
A mixed algorithm of central and upwind difference scheme for the solution of steady/unsteady incompressible Navier-Stokes equations is presented. The algorithm is based on the method of artificial compressibility and...A mixed algorithm of central and upwind difference scheme for the solution of steady/unsteady incompressible Navier-Stokes equations is presented. The algorithm is based on the method of artificial compressibility and uses a third-order flux-difference splitting technique for the convective terms and the second-order central difference for the viscous terms. The numerical flux of semi-discrete equations is computed by using the Roe approximation. Time accuracy is obtained in the numerical solutions by subiterating the equations in pseudotime for each physical time step. The algebraic turbulence model of Baldwin-Lomax is ulsed in this work. As examples, the solutions of flow through two dimensional flat, airfoil, prolate spheroid and cerebral aneurysm are computed and the results are compared with experimental data. The results show that the coefficient of pressure and skin friction are agreement with experimental data, the largest discrepancy occur in the separation region where the lagebraic turbulence model of Baldwin-Lomax could not exactly predict the flow.展开更多
The unsteady 3D flow fields in a single-stage transonic compressor under designed conditions are simulated numerically to investigate the effects of the curved rotors on the stage performance and the aerodynamic inter...The unsteady 3D flow fields in a single-stage transonic compressor under designed conditions are simulated numerically to investigate the effects of the curved rotors on the stage performance and the aerodynamic interaction between the blade rows. The results show that, compared to the compressor with unurved rotors, the compressor under scrutiny acquires remarkable increases in efficiency with significantly reduced amplitudes of the time-dependent fluctuation. The amplitude of the pressure fluctuation around the stator leading edge decreases at both endwalls, but increases at the mid-span in the curved rotors. The pressure fluctuation near the stator leading edge, therefore, becomes more uniform in the radial direction of this compressor. Except for the leading edge area, the pressure fluctuatinn amplitude declines remarkably in the tip region of stator surface downstream of the curved rotor, but hardly changes in the middle and at the hub.展开更多
Gridless method is developed for unsteady viscous flows involving moving boundaries. The point distri- bution of gridless method is implemented in an isotropic or anisotropic way according to the features of viscous f...Gridless method is developed for unsteady viscous flows involving moving boundaries. The point distri- bution of gridless method is implemented in an isotropic or anisotropic way according to the features of viscous flows. In the area far away from the body, the traditional cloud of isotropic points is used, while in the adjacent area, the cloud of anisotropic points is distributed. In this way, the point spacing normal to the wall can be small enough for simulating the boundary layer, and meanwhile, the total number of points in the computational do- main can be controlled due to large spacing in other tangential direction through the anisotropic way. A fast mov- ing technique of clouds of points at each time-step is presented based on the attenuation law of disturbed motion for unsteady flows involving moving boundaries. In the mentioned cloud of points, a uniform weighted least- square curve fit method is utilized to discretize the spatial derivatives of the Navier-Stokes equations. The pro- posed gridless method, coupled with a dual time-stepping method and the Spalart-Allmaras turbulence model, is implemented for the Navier-Stokes equations. The computational results of unsteady viscous flows around a NLR7301 airfoil with an oscillating flap and a pitching NACA0012 airfoil are presented in a good agreement with the available experimental data.展开更多
The impulse and self starting characteristics of a mixed-compression hypersonic inlet designed at Mach number of 6.5 are studied by applying the unsteady computational fluid dynamics （CFD） method. The full Navier–S...The impulse and self starting characteristics of a mixed-compression hypersonic inlet designed at Mach number of 6.5 are studied by applying the unsteady computational fluid dynamics （CFD） method. The full Navier–Stokes equations are solved with the assumption of viscous perfect gas model, and the shear-stress transport （SST） k–x two-equation Reynolds averaged Navier– Stokes （RANS） model is used for turbulence modeling. Results indicate that during impulse starting, the flow field is divided into three zones with different aerodynamic parameters by primary shock and upstream-facing shock. The separation bubble on the shoulder of ramp undergoes a generating, growing, swallowing and disappearing process in sequence. But a separation bubble at the entrance of inlet exists until the freestream velocity is accelerated to the starting Mach number during self starting. The mass flux distribution of flow field is non-uniform because of the interaction between shock and boundary layer, so that the mass flow rate at throat is unsteady during impulse starting. The duration of impulse starting process increases almost linearly with the decrease of freestream Mach number but rises abruptly when the freestream Mach number approaches the starting Mach number. The accelerating performance of booster almost has no influence on the self starting ability of hypersonic inlet.展开更多
The pressure fluctuation in the flow passage of both impeller and casing is addressed on design condition. The initial conditions for the unsteady turbulent simulation are resulted from the steady calculations, and th...The pressure fluctuation in the flow passage of both impeller and casing is addressed on design condition. The initial conditions for the unsteady turbulent simulation are resulted from the steady calculations, and the three dimensional unsteady turbulent simulation concerning the rotor-stator interaction is executed by a Navier-Stoke solver embedded with k -ε turbulence model and with appropriate moving interface boundary conditions. Detecting points are distributed in the flow passage in different radial and circumferential positions to capture the static pressure fluctuation character for one cycle of the impeller. The time-domain spectrums show that the static pressure curves are periodic and have five peaks and five valleys. With the radius increasing, the pressure fluctuation peak-to-peak values in the impeller are increasing, and reach the maximum value on the interface. In the casing flow passage, those values are about 7% of local static pressure except some ones near the tongue. The values become decreasingly in the diffuser pipe. The frequency spectrums transformed by fast Fourier transform （FFT） show that the dominant frequency is approximate with the blade passing frequency, and the pressure fluctuations in impeller passage have high frequency content while those in casing ones have no such information.展开更多
In order to calculate the unsteady aerodynamic characteristics of a tilt-rotor in a conver- sion mode, a virtual blade model （VBM） and an real blade model （RBM） are established respec- tively. A new multi-layer mo...In order to calculate the unsteady aerodynamic characteristics of a tilt-rotor in a conver- sion mode, a virtual blade model （VBM） and an real blade model （RBM） are established respec- tively. A new multi-layer moving-embedded grid technique is proposed to reduce the numerical dissipation of the tilt-rotor wake in a conversion mode. In this method, a grid system generated abound the rotor accounts for rigid blade motions, and a new searching scheme named adaptive inverse map （AIM） is established to search corresponding donor elements in the present moving- embedded grid system to translate information among the different computational zones. A dual-time method is employed to fulfill unsteady calculations on the flowfield of the tilt-rotor, and a second-order centered difference scheme considering artificial viscosity is used to calculate the flux. In order to improve the computing efficiency, the single program multiple data （SPMD） model parallel acceleration technology is adopted, according to the characteristic of the current grid system. The lift and drag coefficients of an NACA0012 airfoil, the dynamic pressure distributions below a typical rotor plane, and the sectional pressure distributions on a three-bladed Branum- Tung tilt-rotor in hover flight are calculated respectively, and the present VBM and RBM are val- idated by comparing the calculated results with available experimental data. Then, unsteady aero- dynamic forces and flowfields of an XV-15 tilt-rotor in different modes, such as a fixed conversion mode at different tilt angles （15°, 30°, 60°） and a whole conversion mode which converses from 0° to 90°, are numerically simulated by the VBM and RBM respectively. By analyses and comparisons on the simulated results of unsteady aerodynamic forces of the tilt-rotor in different modes, some meaningful conclusions about distorted blade-tip vortex distribution and unsteady aerodynamic force variation in a conversion mode are obtained, and these investigation results could provide a good foundation for tilt-rotor aircraft design in the future.展开更多
A full-span free-wake method is coupled with an unsteady panel method to accurately predict the unsteady aerodynamics of helicopter rotor blades in hover and forward flight. The unsteady potential-based panel method i...A full-span free-wake method is coupled with an unsteady panel method to accurately predict the unsteady aerodynamics of helicopter rotor blades in hover and forward flight. The unsteady potential-based panel method is used to consider aerodynamics of finite thickness multi-bladed rotors, and the full-span free-wake method is applied to simulating dynamics of rotor wake. These methods are tightly coupled through trailing-edge Kutta condition and by converting doublet-wake panels to full-span vortex filaments. A velocity-field integration technique is also adopted to overcome singularity problem during the interaction between the rotor wake and blades. Helicopter rotors including Caradonna–Tung, UH-60A, and AH-1G rotors, are simulated in hover and forward flight to validate the accuracy of this approach. The predicted aerodynamic loads of rotor blades agree well with available measured data and computational fluid dynamics （CFD） results, and the unsteady dynamics of rotor wake is also well simulated. Compared to CFD, the present method obtains accurate results more efficiently and is suitable to rotorcraft aeroelastic analysis.展开更多
The fractional calculus approach in the constitutive relationship model of a generalized second grade fluid is introduced.Exact analytical solutions are obtained for a class of unsteady flows for the generalized secon...The fractional calculus approach in the constitutive relationship model of a generalized second grade fluid is introduced.Exact analytical solutions are obtained for a class of unsteady flows for the generalized second grade fluid with the fractional derivative model between two parallel plates by using the Laplace transform and Fourier transform for fractional calculus.The unsteady flows are generated by the impulsive motion or periodic oscillation of one of the plates.In addition,the solutions of the shear stresses at the plates are also determined.展开更多
Unsteady supersonic base flows around three afterbodies, cylindrical （Cy）, boattailed （BT） and three-step （MS）, are investigated in this paper. Reynolds-averaged Navier-Stokes （RANS） and two RANS/LES （large-e...Unsteady supersonic base flows around three afterbodies, cylindrical （Cy）, boattailed （BT） and three-step （MS）, are investigated in this paper. Reynolds-averaged Navier-Stokes （RANS） and two RANS/LES （large-eddy simulation） hybrid methods, detached eddy simulation （DES） and delayed-DES （DDES）, are used to predict the base flow characteristics around the baseline Cy afterbody. All the RANS and hybrid methods are based on the two-equation SST （shear-stress transport） model with compressible corrections （CC）. According to the comparison of measurements, both DES and DDES can produce more satisfactory results than RANS. RANS can only present the ＂stable＂ flow pat- terns, while the hybrid methods can demonstrate unsteady flow structures. DDES and DES results are little different from one another although the latter exhibits better agreement with the experiment. DES is taken to investigate the 5° BT and three-step afterbodies. The mean flow data and the instantaneous turbulent coherent structures are compared against available measurements.展开更多
Abstract Accurate aerodynamic models are the basis of flight simulation and control law design. Mathematically modeling unsteady aerodynamics at high angles of attack bears great difficulties in model structure determ...Abstract Accurate aerodynamic models are the basis of flight simulation and control law design. Mathematically modeling unsteady aerodynamics at high angles of attack bears great difficulties in model structure determination and parameter estimation due to little understanding of the flow mechanism. Support vector machines （SVMs） based on statistical learning theory provide a novel tool for nonlinear system modeling. The work presented here examines the feasibility of applying SVMs to high angle.-of-attack unsteady aerodynamic modeling field. Mainly, after a review of SVMs, several issues associated with unsteady aerodynamic modeling by use of SVMs are discussed in detail, such as sele, ction of input variables, selection of output variables and determination of SVM parameters. The least squares SVM （LS-SVM） models are set up from certain dynamic wind tunnel test data of a delta wing and an aircraft configuration, and then used to predict the aerodynamic responses in other tests. The predictions are in good agreement with the test data, which indicates the satisfving learning and generalization performance of LS-SVMs.展开更多
A type of flow unsteadiness with low frequencies and large amplitude was investigated experimentally for vortex wakes around an ogive-tangent cylinder. The experiments were carried out at angles of attack of 60–80 an...A type of flow unsteadiness with low frequencies and large amplitude was investigated experimentally for vortex wakes around an ogive-tangent cylinder. The experiments were carried out at angles of attack of 60–80 and subcritical Reynolds numbers of 0.6–1.8×105. The reduced frequencies of the unsteadiness are between 0.038 and 0.072, much less than the frequency of Karman vortex shedding. The unsteady flow induces large fluctuations of sectional side forces. The results of pressure measurements and particle image velocimetry indicate that the flow unsteadiness comes from periodic oscillation of the vortex wakes over the slender body. The time-averaged vortex patterns over the slender body are asymmetric, whose orientation is dependent on azimuthal locations of tip perturbations. Therefore, the vortex oscillation is a type of unsteady oscillation around a time-averaged asymmetric vortex structure.展开更多
基金supported by the National Natural Science Foundation of China(Grant Nos.12202488 and 12002377)the Natural Science Program of National University of Defense Technology(Grant No.ZK22-30)Independent Cultivation Project for Young Talents in College of Aerospace Science and Engineering.
文摘Direct numerical simulations of Mach 6 hypersonic flow over a 34°compression corner subject to steady jet are conducted.Distributions of skin friction coefficient,wall pressure,mean velocity and temperature,boundary layer thickness and Stanton number demonstrate that the flow changes dramatically in the shock wave/turbulent boundary layer interaction area.It is found that the steady jet has no effect on suppressing flow separation unexpectedly,but increases its spatial scale instead.Instantaneous flow structures show that the turbulence amplification can be observed after the application of flow control,and abundant Görtler-like vorticities appear,but the strength of the main shock decreases.Analyzing the wall fluctuating pressure signals using weighted power spectral density,we found an interesting thing.That is,although the low-frequency oscillation phenomenon induced by separation shock is suppressed by the steady jet,wall fluctuating pressure beneath the jet shock is oscillating at a frequency lower than 0.1u∞/δref.Results of coherent and intermittency factor reveal that it is related to the backand-forth movement of the jet shock itself.
基金supported by the National Natural Science Foundation of China(Nos.12172315,12072304,11702232)the Fujian Provincial Natural Science Foundation,China(No.2021J01050)the Aeronautical Science Foundation of China(No.20220013068002).
文摘Unsteady aerodynamic characteristics at high angles of attack are of great importance to the design and development of advanced fighter aircraft, which are characterized by post-stall maneuverability with multiple Degrees-of-Freedom(multi-DOF) and complex flow field structure.In this paper, a special kind of cable-driven parallel mechanism is firstly utilized as a new suspension method to conduct unsteady dynamic wind tunnel tests at high angles of attack, thereby providing experimental aerodynamic data. These tests include a wide range of multi-DOF coupled oscillatory motions with various amplitudes and frequencies. Then, for aerodynamic modeling and analysis, a novel data-driven Feature-Level Attention Recurrent neural network(FLAR) is proposed. This model incorporates a specially designed feature-level attention module that focuses on the state variables affecting the aerodynamic coefficients, thereby enhancing the physical interpretability of the aerodynamic model. Subsequently, spin maneuver simulations, using a mathematical model as the baseline, are conducted to validate the effectiveness of the FLAR. Finally, the results on wind tunnel data reveal that the FLAR accurately predicts aerodynamic coefficients, and observations through the visualization of attention scores identify the key state variables that affect the aerodynamic coefficients. It is concluded that the proposed FLAR enhances the interpretability of the aerodynamic model while achieving good prediction accuracy and generalization capability for multi-DOF coupling motion at high angles of attack.
基金funded by the National Natural Science Foundation of China(Nos.U24A20138 and No.52376039)the Beijing Natural Science Foundation,China(JQ24017)+1 种基金the National Science and Technology Major Project of China(Nos.J2019-II-0005-0025 and Y2022-II-0002-0005)the Special Fund for the Member of Youth Innovation Promotion Association of Chinese Academy of Sciences,China(No.2018173).
文摘To assess the aerodynamic performance and vibration characteristics of rotor blades during rotation,a study of unsteady blade surface forces is conducted in a low-speed axial flow compressor under a rotating coordinate system.The capture,modulation,and acquisition of unsteady blade surface forces are achieved by using pressure sensors and strain gauges attached to the rotor blades,in conjunction with a wireless telemetry system.Based on the measurement reliability verification,this approach allows for the determination of the static pressure distribution on rotor blade surfaces,enabling the quantitative description of loadability at different spanwise positions along the blade chord.Effects caused by the factors such as Tip Leakage Flow(TLF)and flow separation can be perceived and reflected in the trends of static pressure on the blade surfaces.Simultaneously,the dynamic characteristics of unsteady pressure and stress on the blade surfaces are analyzed.The results indicate that only the pressure signals measured at the mid-chord of the blade tip can distinctly detect the unsteady frequency of TLF due to the oscillation of the low-pressure spot on the pressure surface.Subsequently,with the help of one-dimensional continuous wavelet analysis method,it can be inferred that as the compressor enters stall,the sensors are capable of capturing stall cell frequency under a rotating coordinate system.Furthermore,the stress at the blade root is higher than that at the blade tip,and the frequency band of the vibration can also be measured by the pressure sensors fixed on the casing wall in a stationary frame.While the compressor stalls,the stress at the blade root can be higher,which can provide valuable guidance for monitoring the lifecycle of compressor blades.
基金supported by the National Science and Technology Major Project,China(No.Y2019-I-0018-0017)the National Natural Science Foundation of China(No.11602200)+1 种基金Hunan Innovative Province Construction Special Fund,China(No.2021GK1020)the Priority Academic Program Development of Jiangsu Higher Education Institutions,China。
文摘Propeller design is a highly intricate and interdisciplinary task that necessitates careful trade-offs between radiated noise levels and aerodynamic efficiency.To achieve efficient trade-off designs,an enhanced on-the-fly unsteady adjoint-based aerodynamic and aeroacoustic optimization methodology is developed,which maintains the fidelity of the Navier-Stokes solution for unsteady flow and of the moving-medium Ffowcs Williams-Hawkings(FW-H)formulation for capturing tonal noise.Furthermore,this on-the-fly approach enables a unified architecture for discreteadjoint sensitivity analysis encompassing both aerodynamics and aeroacoustics,facilitating effective multi-objective weighted optimizations.Subsequently,this proposed methodology is applied to perform trade-off optimizations between aerodynamics and aeroacoustics for a propeller by employing varying weighting factors to comprehend their influence on optimal configurations.The results demonstrate a positive correlation between efficiency and noise sensitivities,and thus indicate an inherent synchronicity where pursing noise reduction through purely aeroacoustic optimization inevitably entails sacrificing aerodynamic efficiency.However,by effectively incorporating appropriate weighting factors(recommended to range from 0.25 to 0.5)into the multi-objective function combined with both aerodynamics and aeroacoustics,it becomes feasible to achieve efficiency enhancement and noise reduction simultaneously.Key findings show that reducing blade planform size and equipping“rotated-S”shaped airfoil profiles in the tip region can effectively restrain noise levels while maintaining aerodynamic performance.
基金supported by the Aeronautical Science Foundation of China(Nos.2018ZA52002,2019ZA052011)the National Natural Science Foundation of China(No.12472236).
文摘Obtaining unsteady hydrodynamic performance is of great significance for seaplane design.Common methods for obtaining unsteady hydrodynamic performance data include tank test and Computational Fluid Dynamics(CFD)numerical simulation,which are costly and time-consuming.Therefore,it is necessary to obtain unsteady hydrodynamic performance in a low-cost and high-precision manner.Due to the strong nonlinearity,complex data distribution,and temporal characteristics of unsteady hydrodynamic performance,the prediction of it is challenging.This paper proposes a Temporal Convolutional Diffusion Model(TCDM)for predicting the unsteady hydrodynamic performance of seaplanes given design parameters.Under the framework of a classifier-free guided diffusion model,TCDM learns the distribution patterns of unsteady hydrodynamic performance data with the designed denoising module based on temporal convolutional network and captures the temporal features of unsteady hydrodynamic performance data.Using CFD simulation data,the proposed method is compared with the alternative methods to demonstrate its accuracy and generalization.This paper provides a method that enables the rapid and accurate prediction of unsteady hydrodynamic performance data,expecting to shorten the design cycle of seaplanes.
基金supported by the National Natural Science Foundation of China(No.92252201)the Fundamental Research Funds for the Central Universitiesthe Academic Excellence Foundation of Beihang University(BUAA)for PhD Students。
文摘Efficient and accurate simulation of unsteady flow presents a significant challenge that needs to be overcome in computational fluid dynamics.Temporal discretization method plays a crucial role in the simulation of unsteady flows.To enhance computational efficiency,we propose the Implicit-Explicit Two-Step Runge-Kutta(IMEX-TSRK)time-stepping discretization methods for unsteady flows,and develop a novel adaptive algorithm that correctly partitions spatial regions to apply implicit or explicit methods.The novel adaptive IMEX-TSRK schemes effectively handle the numerical stiffness of the small grid size and improve computational efficiency.Compared to implicit and explicit Runge-Kutta(RK)schemes,the IMEX-TSRK methods achieve the same order of accuracy with fewer first derivative calculations.Numerical case tests demonstrate that the IMEX-TSRK methods maintain numerical stability while enhancing computational efficiency.Specifically,in high Reynolds number flows,the computational efficiency of the IMEX-TSRK methods surpasses that of explicit RK schemes by more than one order of magnitude,and that of implicit RK schemes several times over.
基金funded by National Natural Science Foundation of China(Mechanisms of proppant-carrying transport by magnetic cross-linked microparticle grids and their degradation patterns in CO_(2) fractured cracks).
文摘This study presents a two-dimensional,transient model to simulate the flow and thermal behavior of CO_(2) within a fracturing wellbore.The model accounts for high-velocity flow within the tubing and radial heat exchange between the wellbore and surrounding formation.It captures the temporal evolution of temperature,pressure,flow velocity,and fluid density,enabling detailed analysis of phase transitions along different tubing sections.The influence of key operational and geological parameters,including wellhead pressure,injection velocity,inlet temperature,and formation temperature gradient,on the wellbore’s thermal and pressure fields is systematically investigated.Results indicate that due to intense convective transport by the high-speed CO_(2) flow,the temperature and velocity within the tubing are primarily governed by the inlet temperature and injection velocity,with relatively minor influence from radial heat transfer with the formation.The pressure,flow velocity,and density of CO_(2) within the tubing are strongly dependent on wellhead conditions.Frictional losses and well depth contribute to pressure variations,particularly in the horizontal section of the wellbore,where a noticeable pressurization effect increases the fluid density.During injection,liquid CO_(2) initially undergoes a rapid transition to a supercritical state,with the depth at which this phase change occurs stabilizing as injection progresses.
文摘A mixed algorithm of central and upwind difference scheme for the solution of steady/unsteady incompressible Navier-Stokes equations is presented. The algorithm is based on the method of artificial compressibility and uses a third-order flux-difference splitting technique for the convective terms and the second-order central difference for the viscous terms. The numerical flux of semi-discrete equations is computed by using the Roe approximation. Time accuracy is obtained in the numerical solutions by subiterating the equations in pseudotime for each physical time step. The algebraic turbulence model of Baldwin-Lomax is ulsed in this work. As examples, the solutions of flow through two dimensional flat, airfoil, prolate spheroid and cerebral aneurysm are computed and the results are compared with experimental data. The results show that the coefficient of pressure and skin friction are agreement with experimental data, the largest discrepancy occur in the separation region where the lagebraic turbulence model of Baldwin-Lomax could not exactly predict the flow.
基金National Natural Science Foundation of China (506460210) Chinese Specialized Research Fund for the Doctoral Program of Higher Education (20060213007)Development Program for Outstanding Young Teachers in Harbin Institute of Technology (HITQNJS.2006.046)
文摘The unsteady 3D flow fields in a single-stage transonic compressor under designed conditions are simulated numerically to investigate the effects of the curved rotors on the stage performance and the aerodynamic interaction between the blade rows. The results show that, compared to the compressor with unurved rotors, the compressor under scrutiny acquires remarkable increases in efficiency with significantly reduced amplitudes of the time-dependent fluctuation. The amplitude of the pressure fluctuation around the stator leading edge decreases at both endwalls, but increases at the mid-span in the curved rotors. The pressure fluctuation near the stator leading edge, therefore, becomes more uniform in the radial direction of this compressor. Except for the leading edge area, the pressure fluctuatinn amplitude declines remarkably in the tip region of stator surface downstream of the curved rotor, but hardly changes in the middle and at the hub.
基金Supported by the National Natural Science Foundation of China(10372043,11172134)the Fundingof Jiangsu Innovation Program for Graduate Education(CXZZ11-0192)~~
文摘Gridless method is developed for unsteady viscous flows involving moving boundaries. The point distri- bution of gridless method is implemented in an isotropic or anisotropic way according to the features of viscous flows. In the area far away from the body, the traditional cloud of isotropic points is used, while in the adjacent area, the cloud of anisotropic points is distributed. In this way, the point spacing normal to the wall can be small enough for simulating the boundary layer, and meanwhile, the total number of points in the computational do- main can be controlled due to large spacing in other tangential direction through the anisotropic way. A fast mov- ing technique of clouds of points at each time-step is presented based on the attenuation law of disturbed motion for unsteady flows involving moving boundaries. In the mentioned cloud of points, a uniform weighted least- square curve fit method is utilized to discretize the spatial derivatives of the Navier-Stokes equations. The pro- posed gridless method, coupled with a dual time-stepping method and the Spalart-Allmaras turbulence model, is implemented for the Navier-Stokes equations. The computational results of unsteady viscous flows around a NLR7301 airfoil with an oscillating flap and a pitching NACA0012 airfoil are presented in a good agreement with the available experimental data.
文摘The impulse and self starting characteristics of a mixed-compression hypersonic inlet designed at Mach number of 6.5 are studied by applying the unsteady computational fluid dynamics （CFD） method. The full Navier–Stokes equations are solved with the assumption of viscous perfect gas model, and the shear-stress transport （SST） k–x two-equation Reynolds averaged Navier– Stokes （RANS） model is used for turbulence modeling. Results indicate that during impulse starting, the flow field is divided into three zones with different aerodynamic parameters by primary shock and upstream-facing shock. The separation bubble on the shoulder of ramp undergoes a generating, growing, swallowing and disappearing process in sequence. But a separation bubble at the entrance of inlet exists until the freestream velocity is accelerated to the starting Mach number during self starting. The mass flux distribution of flow field is non-uniform because of the interaction between shock and boundary layer, so that the mass flow rate at throat is unsteady during impulse starting. The duration of impulse starting process increases almost linearly with the decrease of freestream Mach number but rises abruptly when the freestream Mach number approaches the starting Mach number. The accelerating performance of booster almost has no influence on the self starting ability of hypersonic inlet.
基金supported by National Outstanding Young Scientists Funds of China (Grand No.50825902)
文摘The pressure fluctuation in the flow passage of both impeller and casing is addressed on design condition. The initial conditions for the unsteady turbulent simulation are resulted from the steady calculations, and the three dimensional unsteady turbulent simulation concerning the rotor-stator interaction is executed by a Navier-Stoke solver embedded with k -ε turbulence model and with appropriate moving interface boundary conditions. Detecting points are distributed in the flow passage in different radial and circumferential positions to capture the static pressure fluctuation character for one cycle of the impeller. The time-domain spectrums show that the static pressure curves are periodic and have five peaks and five valleys. With the radius increasing, the pressure fluctuation peak-to-peak values in the impeller are increasing, and reach the maximum value on the interface. In the casing flow passage, those values are about 7% of local static pressure except some ones near the tongue. The values become decreasingly in the diffuser pipe. The frequency spectrums transformed by fast Fourier transform （FFT） show that the dominant frequency is approximate with the blade passing frequency, and the pressure fluctuations in impeller passage have high frequency content while those in casing ones have no such information.
基金supported by the National Natural Science Foundation of China(No.11272150)
文摘In order to calculate the unsteady aerodynamic characteristics of a tilt-rotor in a conver- sion mode, a virtual blade model （VBM） and an real blade model （RBM） are established respec- tively. A new multi-layer moving-embedded grid technique is proposed to reduce the numerical dissipation of the tilt-rotor wake in a conversion mode. In this method, a grid system generated abound the rotor accounts for rigid blade motions, and a new searching scheme named adaptive inverse map （AIM） is established to search corresponding donor elements in the present moving- embedded grid system to translate information among the different computational zones. A dual-time method is employed to fulfill unsteady calculations on the flowfield of the tilt-rotor, and a second-order centered difference scheme considering artificial viscosity is used to calculate the flux. In order to improve the computing efficiency, the single program multiple data （SPMD） model parallel acceleration technology is adopted, according to the characteristic of the current grid system. The lift and drag coefficients of an NACA0012 airfoil, the dynamic pressure distributions below a typical rotor plane, and the sectional pressure distributions on a three-bladed Branum- Tung tilt-rotor in hover flight are calculated respectively, and the present VBM and RBM are val- idated by comparing the calculated results with available experimental data. Then, unsteady aero- dynamic forces and flowfields of an XV-15 tilt-rotor in different modes, such as a fixed conversion mode at different tilt angles （15°, 30°, 60°） and a whole conversion mode which converses from 0° to 90°, are numerically simulated by the VBM and RBM respectively. By analyses and comparisons on the simulated results of unsteady aerodynamic forces of the tilt-rotor in different modes, some meaningful conclusions about distorted blade-tip vortex distribution and unsteady aerodynamic force variation in a conversion mode are obtained, and these investigation results could provide a good foundation for tilt-rotor aircraft design in the future.
文摘A full-span free-wake method is coupled with an unsteady panel method to accurately predict the unsteady aerodynamics of helicopter rotor blades in hover and forward flight. The unsteady potential-based panel method is used to consider aerodynamics of finite thickness multi-bladed rotors, and the full-span free-wake method is applied to simulating dynamics of rotor wake. These methods are tightly coupled through trailing-edge Kutta condition and by converting doublet-wake panels to full-span vortex filaments. A velocity-field integration technique is also adopted to overcome singularity problem during the interaction between the rotor wake and blades. Helicopter rotors including Caradonna–Tung, UH-60A, and AH-1G rotors, are simulated in hover and forward flight to validate the accuracy of this approach. The predicted aerodynamic loads of rotor blades agree well with available measured data and computational fluid dynamics （CFD） results, and the unsteady dynamics of rotor wake is also well simulated. Compared to CFD, the present method obtains accurate results more efficiently and is suitable to rotorcraft aeroelastic analysis.
基金The project supported by the National Natural Science Foundation of China (10372007,10002003) and CNPC Innovation Fund
文摘The fractional calculus approach in the constitutive relationship model of a generalized second grade fluid is introduced.Exact analytical solutions are obtained for a class of unsteady flows for the generalized second grade fluid with the fractional derivative model between two parallel plates by using the Laplace transform and Fourier transform for fractional calculus.The unsteady flows are generated by the impulsive motion or periodic oscillation of one of the plates.In addition,the solutions of the shear stresses at the plates are also determined.
基金supported by the National Natural Science Foundation of China (10502030 and 90505005)Innovation and Support Foundation of Chinese Astronautics
文摘Unsteady supersonic base flows around three afterbodies, cylindrical （Cy）, boattailed （BT） and three-step （MS）, are investigated in this paper. Reynolds-averaged Navier-Stokes （RANS） and two RANS/LES （large-eddy simulation） hybrid methods, detached eddy simulation （DES） and delayed-DES （DDES）, are used to predict the base flow characteristics around the baseline Cy afterbody. All the RANS and hybrid methods are based on the two-equation SST （shear-stress transport） model with compressible corrections （CC）. According to the comparison of measurements, both DES and DDES can produce more satisfactory results than RANS. RANS can only present the ＂stable＂ flow pat- terns, while the hybrid methods can demonstrate unsteady flow structures. DDES and DES results are little different from one another although the latter exhibits better agreement with the experiment. DES is taken to investigate the 5° BT and three-step afterbodies. The mean flow data and the instantaneous turbulent coherent structures are compared against available measurements.
文摘Abstract Accurate aerodynamic models are the basis of flight simulation and control law design. Mathematically modeling unsteady aerodynamics at high angles of attack bears great difficulties in model structure determination and parameter estimation due to little understanding of the flow mechanism. Support vector machines （SVMs） based on statistical learning theory provide a novel tool for nonlinear system modeling. The work presented here examines the feasibility of applying SVMs to high angle.-of-attack unsteady aerodynamic modeling field. Mainly, after a review of SVMs, several issues associated with unsteady aerodynamic modeling by use of SVMs are discussed in detail, such as sele, ction of input variables, selection of output variables and determination of SVM parameters. The least squares SVM （LS-SVM） models are set up from certain dynamic wind tunnel test data of a delta wing and an aircraft configuration, and then used to predict the aerodynamic responses in other tests. The predictions are in good agreement with the test data, which indicates the satisfving learning and generalization performance of LS-SVMs.
基金supported by the National Natural Science Foundation of China(No.11272033)
文摘A type of flow unsteadiness with low frequencies and large amplitude was investigated experimentally for vortex wakes around an ogive-tangent cylinder. The experiments were carried out at angles of attack of 60–80 and subcritical Reynolds numbers of 0.6–1.8×105. The reduced frequencies of the unsteadiness are between 0.038 and 0.072, much less than the frequency of Karman vortex shedding. The unsteady flow induces large fluctuations of sectional side forces. The results of pressure measurements and particle image velocimetry indicate that the flow unsteadiness comes from periodic oscillation of the vortex wakes over the slender body. The time-averaged vortex patterns over the slender body are asymmetric, whose orientation is dependent on azimuthal locations of tip perturbations. Therefore, the vortex oscillation is a type of unsteady oscillation around a time-averaged asymmetric vortex structure.
