Flutter suppression is an important measure to improve fatigue life and enhance the performance of aircraft in modern aircraft design. In order to design more effective controllers for flutter suppression with high efficiency, an efficient reduced-order framework for active/passive hybrid flutter suppression is proposed. The traditional CFD-based ROMs have been successfully applied to active flutter suppression with high accuracy and efficiency. But, when a structure modification is made such as in aeroelastic tailoring and aeroelastic structural optimisation, the structural model should be updated, and the expensive, time-consuming CFD-based ROMs have to be reconstructed; such a process is impractical for passive flutter suppression. To overcome the realistic challenge, an efficient reduced-order framework for active/passive hybrid flutter suppression is proposed by extending an efficient aeroelastic CFD-based POD/ROM which we have developed. The proposed framework is demonstrated and evaluated using an improved AGARD 445.6 wing model. The results show that the proposed framework can accurately predict the aeroelastic response for active/passive hybrid flutter suppression with high efficiency. It provides a powerful tool for active/passive hybrid flutter suppression, and therefore, is ideally suited to design more effective controllers, and may have the potential to reduce the overall cost of aircraft design.

A fast numerical method for unsteady aerodynamic calculation of 3D wing is established, which is suitable for the preliminary design. Based on the lifting-line method, the aerodynamic data of the 2D aerofoil obtained by the unsteady CFD simulation is used as the model input to solve the aerodynamic force of the 3D wing. Compared with the traditional steady lifting-line method, the augmented method adopts the unsteady Kutta-Jouowski (K-J) theorem to calculate the circulation and improve the accuracy of the method through the circulation correction. The pitching motion of 3D wing at different aspect ratio and reduction frequencies are studied. The results show that the aerodynamic forces obtained by the augmented lifting-line method have good agreement with the 3D unsteady CFD calculations. Compared with 3D CFD calculation, the calculation efficiency of the improved method is increased by more than 12 times. The improved method has extensive applicability and can be used to estimate the unsteady aerodynamic forces of 3D single or multiple wing configurations.

The simulations and assessment of transient performance of gas turbine engines during the conceptual and preliminary design stage may be conducted ignoring heat soakage and tip clearance variations due to lack of detailed geometrical and structural information. As a result, problems with transient performance stability may not be revealed correctly, and corresponding design iterations would be necessary and costly when those problems are revealed at a detailed design stage. To make an engine design more cost and time effective, it has become important to require better transient performance simulations during the conceptual and preliminary design stage considering all key impact factors such as fuel control schedule, rotor dynamics, inter-component volume effect as well as heat soakage and tip clearance variation effects. In this research, a novel transient performance simulation approach with generically simplified heat soakage and tip clearance models for major gas path components of gas turbine engines including compressors, turbines and combustors has been developed to support more realistic transient performance simulations of gas turbine engines at conceptual and preliminary design stages. Such heat soakage and tip clearance models only require thermodynamic design parameters as input, which is normally available during such design stages. The models have been implemented into in-house transient performance simulation software and applied to a model twin-spool turbojet engine to test their effectiveness. Comparisons between transient performance simulated with and without the heat soakage and tip clearance effects demonstrate that the results are promising. Although the introduced heat soakage and tip clearance models may not be as accurate as that using detailed component geometrical information, it is able to include the major heat soakage and tip clearance effects and make the transient performance simulations and analysis more realistic during conceptual and preliminary engine design stage.

Falling paper flutters and tumbles through air, whereas a paper airplane glides smoothly if its leading edge is appropriately weighted. We investigate this transformation from ‘plain paper’ to ‘paper plane’ through experiments, aerodynamic modelling and free flight simulations of thin plates with differing centre of mass (CoM) locations. Periodic modes such as fluttering, tumbling and bounding give way to steady gliding and then downward diving as the CoM is increasingly displaced towards one edge. To explain these observations, we formulate a quasi-steady aerodynamic model whose force and torque coefficients are informed by experimental measurements. The dependencies on angle of attack reflect the transition from attached to separated flow and a dynamic centre of pressure, effects that prove critical to reproducing the observed motions of paper planes in air and plates in water. Because the model successfully accounts for unsteady and steady flight modes, it may be usefully applied to further problems involving actuated motions, feedback control and interactions with ambient flows.

To investigate the downstream rim seal gas ingestion characteristics of a 1.5-stage turbine, the URANS equations were solved numerically using the SST turbulence model. The effects of different purge flow rates and the second vane on the ingestion characteristics of the aft cavity and the nonuniform fluctuations of the main gas path pressure are analysed. The results showed that the aft cavity is affected by the combined effects of the blade and the second vane, and the potential field at the leading edge of the second vane greatly influence the airflow variation in the aft cavity, which enhances the ingress of the mainstream into the wheel-space. The front purge flow weakens the egress between the suction side of the blade and the suction side of the second vane. The potential field at the leading edge of the second vane suppresses the nonuniform distribution of airflow in the aft cavity caused by the rotational effect of the blade.

The spatio-temporal variation of leaf chlorophyll content is an important crop phenotypic trait that is of great significance for evaluating crop productivity. This study used a soil-plant analysis development (SPAD) chlorophyll meter for non-destructive monitoring of leaf chlorophyll dynamics to characterize the patterns of spatio-temporal variation in the nutritional status of maize (Zea mays L.) leaves under three nitrogen treatments in two cultivars. The results showed that nitrogen levels could affect the maximum leaf SPAD reading (SPADmax) and the duration of high SPAD reading. A rational model was used to measure the changes in SPAD readings over time in single leaves. This model was suitable for predicting the dynamics of the nutrient status for each leaf position under different nitrogen treatments, and model parameter values were position dependent. SPADmax at each leaf decreased with the reduction of nitrogen supply. Leaves at different positions in both cultivars responded differently to higher nitrogen rates. Lower leaves (8th–10th positions) were more sensitive than the other leaves in response to nitrogen. Monitoring the SPAD reading dynamic of lower leaves could accurately characterize and assess the nitrogen supply in plants. The lower leaves in nitrogen-deficient plants had a shorter duration of high SPAD readings compared to nitrogen-sufficient plants; this physiological mechanism should be studied further. In summary, the spatio-temporal variation of plant nitrogen status in maize was analysed to determine critical leaf positions for potentially assisting in the identification of appropriate agronomic management practices, such as the adjustment of nitrogen rates in late fertilization.

The epidemic of tuberculosis has posed a serious burden in Qinghai province, it is necessary to clarify the epidemiological characteristics and spatial-temporal distribution of TB for future prevention and control measures. We used descriptive epidemiological methods and spatial statistical analysis including spatial correlation and spatial-temporal analysis in this study. Furthermore, we applied an exponential smoothing model for TB epidemiological trend forecasting. Of 43 859 TB cases, the sex ratio was 1.27:1 (M:F), and the average annual TB registered incidence was 70.00/100 000 of 2009–2019. More cases were reported in March and April, and the worst TB stricken regions were the prefectures of Golog and Yushu. High TB registered incidences were seen in males, farmers and herdsmen, Tibetans, or elderly people. 7132 cases were intractable, which were recurrent, drug resistant, or co-infected with other infections. Three likely cases clusters with significant high risk were found by spatial-temporal scan on data of 2009–2019. The exponential smoothing winters' additive model was selected as the best-fitting model to forecast monthly TB cases in the future. This research indicated that TB in Qinghai is still a serious threaten to the local residents' health. Multi-departmental collaboration and funds special for TB treatments and control are still needed, and the exponential smoothing model is promising which could be applied for forecasting of TB epidemic trend in this high-altitude province.

We report on experimental observation of non-laminar proton acceleration modulated by a strong magnetic field in laser irradiating micrometer aluminum targets. The results illustrate the coexistence of ring-like and filamentation structures. We implement the knife edge method into the radiochromic film detector to map the accelerated beams, measuring a source size of 30–110 μm for protons of more than 5 MeV. The diagnosis reveals that the ring-like profile originates from low-energy protons far off the axis whereas the filamentation is from the near-axis high-energy protons, exhibiting non-laminar features. Particle-in-cell simulations reproduced the experimental results, showing that the short-term magnetic turbulence via Weibel instability and the long-term quasi-static annular magnetic field by the streaming electric current account for the measured beam profile. Our work provides direct mapping of laser-driven proton sources in the space-energy domain and reveals the non-laminar beam evolution at featured time scales.

The present paper uses the detailed flow data produced by direct numerical simulation (DNS) of a three-dimensional, spatially developing plane free shear layer to assess several commonly used turbulence models in compressible flows. The free shear layer is generated by two parallel streams separated by a splitter plate, with a naturally developing inflow condition. The DNS is conducted using a high-order discontinuous spectral element method (DSEM) for various convective Mach numbers. The DNS results are employed to provide insights into turbulence modelling. The analyses show that with the knowledge of the Reynolds velocity fluctuations and averages, the considered strong Reynolds analogy models can accurately predict temperature fluctuations and Favre velocity averages, while the extended strong Reynolds analogy models can correctly estimate the Favre velocity fluctuations and the Favre shear stress. The pressure–dilatation correlation and dilatational dissipation models overestimate the corresponding DNS results, especially with high compressibility. The pressure–strain correlation models perform excellently for most pressure–strain correlation components, while the compressibility modification model gives poor predictions. The results of an a priori test for subgrid-scale (SGS) models are also reported. The scale similarity and gradient models, which are non-eddy viscosity models, can accurately reproduce SGS stresses in terms of structure and magnitude. The dynamic Smagorinsky model, an eddy viscosity model but based on the scale similarity concept, shows acceptable correlation coefficients between the DNS and modelled SGS stresses. Finally, the Smagorinsky model, a purely dissipative model, yields low correlation coefficients and unacceptable accumulated errors.

The extent of the reduction of maize (Zea mays L.) kernel moisture content through drying is closely related to field temperature (or accumulated temperature; AT) following maturation. In 2017 and 2018, we selected eight maize hybrids that are widely planted in Northeastern China to construct kernel drying prediction models for each hybrid based on kernel drying dynamics. In the traditional harvest scenario using the optimal sowing date (OSD), maize kernels underwent drying from 4th September to 5th October, with variation coefficients of 1.0–1.9. However, with a latest sowing date (LSD), drying occurred from 14th September to 31st October, with variation coefficients of 1.3–3.0. In the changed harvest scenario, the drying time of maize sown on the OSD condition was from 12th September to 9th November with variation coefficients of 1.3–3.0, while maize sown on the LSD had drying dates of 26th September to 28th October with variation coefficients of 1.5–3.6. In the future harvest scenario, the Fengken 139 (FK139) and Jingnongke 728 (JNK728) hybrids finished drying on 20th October and 8th November, respectively, when sown on the OSD and had variation coefficients of 2.7–2.8. Therefore, the maize kernel drying time was gradually delayed and was associated with an increased demand for AT ⩾ 0°C late in the growing season. Furthermore, we observed variation among different growing seasons likely due to differences in weather patterns, and that sowing dates impact variations in drying times to a greater extent than harvest scenarios.

Pneumatic launch systems for Unmanned Aerial Vehicles (UAVs), including mechanical and pneumatic systems, are complex and non-linear. They are subjected to system parameters during launch, which leads to difficulty in engineering research analysis. For example, the mismatch between the UAV parameters and the parameter design indices of the launch system as well as the unclear design indices of the launching speed and overload of UAVs have a great impact on launch safety. Considering this situation, some studies are presented in this paper. Taking the pneumatic launch system as a research object, a pneumatic launcher dynamic simulation model is built based on co-simulation considering the coupling characteristics of the mechanical structure and transmission system. Its accuracy was verified by laboratory test results. Based on this model, the paper shows the effects of the key parameters, including the mass of the UAV, cylinder volume, pressure and moment of inertia of the pulley block, on the performance of the dynamic characteristics of the launch process. Then, a method for matching the parameter characteristics between the UAV and launch system based on batch simulation is proposed. The set of matching parameter values of the UAV and launch system that satisfy the launch take-off safety criteria are calculated. Finally, the influence of the system parameters of the launch process on the launch performance was analysed in detail, and the design optimised. Meaningful conclusions were obtained. The analysis method and its results can provide a reference for engineering and theoretical research and development of pneumatic launch systems.

We investigate the European call option pricing problem under the fractional stochastic volatility model. The stochastic volatility model is driven by both fractional Brownian motion and standard Brownian motion. We obtain an analytical solution of the European option price via the Itô’s formula for fractional Brownian motion, Malliavin calculus, derivative replication and the fundamental solution method. Some numerical simulations are given to illustrate the impact of parameters on option prices, and the results of comparison with other models are presented.

We find a new refinement of Fine’s partition theorem on partitions into distinct parts with the minimum part odd. As a consequence, we obtain two companion partition identities. Both analytic and combinatorial proofs are provided.

The compressibility effects on energy exchange mechanisms in a three-dimensional, spatially developing plane free shear layer are investigated via data produced by direct numerical simulation. The compressible shear layer is simulated using a high-order discontinuous spectral element method for convective Mach numbers $M_c = 0.3$, 0.5 and 0.7. The energy exchange mechanisms in the flow are examined by analysing the budget terms of mean kinetic, internal and turbulent kinetic energy transport equations, in both transition and turbulent regions. The results show that turbulent production, turbulent viscous dissipation, mean viscous dissipation, pressure dilatation and enthalpic production are the main mechanisms responsible for energy exchange among different forms of energy. The effects of compressibility on energy transfer mechanisms are studied based on the analyses of those five budget terms. The primary budget terms evolve differently in the transition and turbulent regions and change significantly for varying compressibility. In the transition region, a double-peak variation becomes a single peak in the streamwise profile of the turbulent production as $M_c$ increases from 0.3 to 0.7, due to significant changes in the vortex pairing structures. The shear layer centre slightly shifts to the high-speed side due to the appearance of the velocity deficit. The velocity deficit presence distance (VDPD) becomes longer as compressibility increases. However, in the turbulent region, the cross-stream profiles of the main budget terms significantly shift to the low-speed side because of the asymmetric mass entrainment and shift even further as $M_c$ increases.

The aim of this study was to explore the frequency and distribution of gene mutations that are related to isoniazid (INH) and rifampin (RIF)-resistance in the strains of the multidrug-resistant tuberculosis (MDR-TB) Mycobacterium tuberculosis (M.tb) in Beijing, China. In this retrospective study, the genotypes of 173 MDR-TB strains were analysed by spoligotyping. The katG, inhA genes and the promoter region of inhA, in which genetic mutations confer INH resistance; and the rpoB gene, in which genetic mutations confer RIF resistance, were sequenced. The percentage of resistance-associated nucleotide alterations among the strains of different genotypes was also analysed. In total, 90.8% (157/173) of the MDR strains belonged to the Beijing genotype. Population characteristics were not significantly different among the strains of different genotypes. In total, 50.3% (87/173) strains had mutations at codon S315T of katG; 16.8% (29/173) of strains had mutations in the inhA promoter region; of them, 5.5% (15/173) had point mutations at −15 base (C→T) of the inhA promoter region. In total, 86.7% (150/173) strains had mutations at rpoB gene; of them, 40% (69/173) strains had mutations at codon S531L of rpoB. The frequency of mutations was not significantly higher in Beijing genotypic MDR strains than in non-Beijing genotypes. Beijing genotypic MDR-TB strains were spreading in Beijing and present a major challenge to TB control in this region. A high prevalence of katG Ser315Thr, inhA promoter region (−15C→T) and rpoB (S531L) mutations was observed. Molecular diagnostics based on gene mutations was a useful method for rapid detection of MDR-TB in Beijing, China.

A new near-infrared direct acceleration mechanism driven by Laguerre–Gaussian laser is proposed to stably accelerate and concentrate electron slice both in longitudinal and transversal directions in vacuum. Three-dimensional simulations show that a 2-μm circularly polarized ${\mathrm{LG}}_p^l$ (p = 0, l = 1, σz = −1) laser can directly manipulate attosecond electron slices in additional dimensions (angular directions) and give them annular structures and angular momentums. These annular vortex attosecond electron slices are expected to have some novel applications such as in the collimation of antiprotons in conventional linear accelerators, edge-enhancement electron imaging, structured X-ray generation, and analysis and manipulation of nanomaterials.

The present study was undertaken to evaluate the influence of rumen-protected folic acid (RPFA) on slaughter performance, visceral organ and gastrointestinal tract coefficients, and meat quality in lambs. Sixty-six lambs from 120 Hu ewes were selected based on body weight and maternal diets and then assigned to six groups using a randomised block experimental design in a 3 × 2 factorial arrangement. The first factor was folic acid (FA) as RPFA in the maternal diet (0 mg/kg (M0F), 16 mg/kg (M16F) or 32 mg/kg (M32F) on DM basis). The second factor was FA in the lambs’ diet from weaning until slaughter (0 mg/kg (OC) or 4·0 mg/kg (OF)). The results indicated that the addition of 16 mg/kg FA to the maternal diet increased pre-slaughter weight (PSW), dressing and meat percentage, the reticulum and omasum coefficients, length of the jejunum and ileum, tail fat and perirenal fat coefficient and a* value of the meat colour. The addition of RPFA to the lambs’ diet increased PSW, dressing and meat percentage, eye muscle area, abomasum weight, weight and length of the small intestine, but reduced the coefficients of tail fat. An M × O interaction was observed for the weights of heart, lungs, rumen and total stomach, weight and coefficient of omental fat and the girth rib value. Collectively, RPFA in the maternal and lambs’ diet improved slaughter performance and meat quality by stimulating the morphological development of the gastrointestinal tract and the distribution of fat in the body.