The enrichment and accumulation of natural gas hydrates depend on sufficient gas supply and effective migration pathways.The upward migration of deep thermogenic gases through fault systems is critical for seepage-typ...The enrichment and accumulation of natural gas hydrates depend on sufficient gas supply and effective migration pathways.The upward migration of deep thermogenic gases through fault systems is critical for seepage-type hydrate formation.This study aims to elucidate the developmental characteristics of Cenozoic fault systems in the eastern offshore area of Dongsha Island and their influence on natural gas hydrate formation.Utilizing high-resolution 3D seismic data,this study conducted a detailed structural interpretation and seismic attribute analysis to systematically investigate the spatial distribution,developmental stages,and dynamic mechanisms of the Cenozoic fault systems in this region.In addition,this study explored the role of these fault systems in facilitating the migration of deep thermogenic gases to the shallow strata.The study area is dominated by extensional and transtensional normal faults characterized by inherited development and relatively small fault displacements.The Cenozoic strata exhibit a tectonic framework of block-faulted uplift and subsidence with alternating highs and lows.Faults on either side of the central uplift dip in opposite directions and commonly exhibit parallel,step-like patterns.Differences in fault system attitudes were observed between the southern and northern parts of the study area.In the south,fault strikes remained consistent from deep to shallow levels,predominantly trending NE and NEE.In the north,fault strikes varied significantly with depth,transitioning from predominantly NEE in deeper strata to EW and NWW in shallower strata.The study identifies two distinct phases of Cenozoic fault activity:(1)66–10 Ma,a regional extensional tectonic regime controlled fault development,resulting in the formation of NEE-trending normal faults;(2)10–2.6 Ma,the Dongsha Movement influenced fault activity,during which EW-and NW–W-trending transtensional faults with dextral strike-slip characteristics developed in the Miocene strata of the northern region.The Cenozoic fault system played a significant positive role in facilitating the migration of deep thermogenic gas to shallow levels,thereby enabling the formation of natural gas hydrates.展开更多
The sampling method is an important numerical technique for solving reliability problems in engineering systems.However,the evaluation of the failure probability using classical sampling methods is time-consuming for ...The sampling method is an important numerical technique for solving reliability problems in engineering systems.However,the evaluation of the failure probability using classical sampling methods is time-consuming for complex engineering structure.To address this issue,this paper proposes a gradient optimization assisted bubble sampling method(GOBSM)to reduce the computational costs,which enhances the coverage range of bubbles,thereby improving the computational efficiency without sacrificing the accuracy.Furthermore,the bubble gradient iterative algorithm is developed to efficiently construct bubbles.Eight complex numerical examples are tested for assessing the failure probability,and the results demonstrate the performance of GOBSM.展开更多
To the Editor:Local ablation is crucial in the treatment of hepatocellular car-cinoma(HCC)and is widely regarded as a good alternative to surgical intervention for achieving local control of HCC.However,the most commo...To the Editor:Local ablation is crucial in the treatment of hepatocellular car-cinoma(HCC)and is widely regarded as a good alternative to surgical intervention for achieving local control of HCC.However,the most commonly used methods for treating HCC are radiofre-quency ablation(RFA)and microwave ablation(MWA),which are thermal ablation methods with limited applicability for tumors in special locations.For example,thermal damage to adjacent impor-tant structures(such as blood vessels,bile ducts and intestines)leads to serious complications.Moreover,thermal ablation of tu-mors near large blood vessels is associated with a higher incidence of incomplete eradication(i.e.,the heat-sink effect)[1,2].展开更多
The stability of gaseous detonation waves is crucial for the operation of detonation-based propulsion systems and the assessment of industrial explosion hazards.However,research on the stability of detonation waves in...The stability of gaseous detonation waves is crucial for the operation of detonation-based propulsion systems and the assessment of industrial explosion hazards.However,research on the stability of detonation waves in complex reactive systems that are composed of actual fuels and oxidants and can be described by numerous elementary chemical reactions,has not been fully carried out.To investigate the relationship between linear and nonlinear stabilities in gaseous detonation wave propagation for complex reactive systems,the linear stability analysis and the one-dimensionally nonlinear numerical simulations of H2/O2/Ar(argon)detonations based on the reactive Euler equations and detailed reaction mechanisms are carried out.The results show that in complex reactive systems characterized by elementary chemical reactions,the results of linear stability computation of detonation are consistent with those from one-dimensionally nonlinear oscillations of detonation wave.Utilizing these linear stability results,a neutral stability curve and a perturbation frequency transition curve in the phase plane of initial pressure versus inert gas(Ar)dilution ratio are derived,especially the new frequency transition curve clearly describes the transition of perturbations from low-frequency to high-frequency mode.One-dimensional nonlinear simulations show that near the perturbation frequency transition curve,the oscillations of the detonation wave can also transform between the lowfrequency,high-amplitude oscillation mode and the high-frequency,low-amplitude oscillation mode,with the oscillation frequency corresponding to the mode that exhibits the maximum growth rate identified in the linear stability analysis.This investigation into detonation stability in complex reactive gases offers guidance for selecting appropriate initial conditions and gas compositions in practical applications of detonation.展开更多
The effect of micro-bubbles on the turbulent boundary layer in the channel flow with Reynolds numbers (Re) ranging from 0.87 × 105 to 1.23 × 105 is experimentally studied by using particle image velocimetry ...The effect of micro-bubbles on the turbulent boundary layer in the channel flow with Reynolds numbers (Re) ranging from 0.87 × 105 to 1.23 × 105 is experimentally studied by using particle image velocimetry (PIV) measurements.The microbubbles are produced by water electrolysis.The velocity profiles,Reynolds stress and instantaneous structures of the boundary layer,with and without micro-bubbles,are measured and analyzed.The presence of micro-bubbles changes the streamwise mean velocity of the fluid and increases the wall shear stress.The results show that micro-bubbles have two effects,buoyancy and extrusion,which dominate the flow behavior of the mixed fluid in the turbulent boundary layer.The buoyancy effect leads to upward motion that drives the fluid motion in the same direction and,therefore,enhances the turbulence intense of the boundary layer.While for the extrusion effect,the presence of accumulated micro-bubbles pushes the flow structures in the turbulent boundary layer away from the near-wall region.The interaction between these two effects causes the vorticity structures and turbulence activity to be in the region far away from the wall.The buoyancy effect is dominant when the Re is relatively small,while the extrusion effect plays a more important role when Re rises.展开更多
An oblique detonation wave for a Mach 7 inlet flow over a long enough wedge of 30 turning angle is simulated numerically using Euler equation and one-step rection model.The fifth-order WENO scheme is adopted to captur...An oblique detonation wave for a Mach 7 inlet flow over a long enough wedge of 30 turning angle is simulated numerically using Euler equation and one-step rection model.The fifth-order WENO scheme is adopted to capture the shock wave.The numerical results show that with the compression of the wedge wall the detonation wave front structure is divided into three sections:the ZND model-like strcuture,single-sided triple point structure and dual-headed triple point strucuture.The first structure is the smooth straight,and the second has the characteristic of the triple points propagating dowanstream only with the same velocity,while the dual-headed triple point structure is very complicated.The detonation waves facing upstream and downstream propagate with different velocities,in which the periodic collisions of the triple points cause the oscillation of the detonation wave front.This oscillation process has temporal and spatial periodicity.In addition,the triple point trace are recorded to obtain different cell structures in three sections.展开更多
The three-dimensional interactions of a perturbed premixed flame interface with a planar incident shock wave and its reflected shock waves are numerically simulated by solving the compressible,reactive Navier-Stokes e...The three-dimensional interactions of a perturbed premixed flame interface with a planar incident shock wave and its reflected shock waves are numerically simulated by solving the compressible,reactive Navier-Stokes equations with the high-resolution scheme and a single-step chemical reaction.The effects of the initial incident shock wave strength (Mach number) and the initial perturbation pattern of interface on the interactions are investigated.The distinct properties of perturbation growth on the flame interface during the interactions are presented.Our results show that perturbation growth is mainly attributed to the flame stretching and propagation.The flame stretching is associated with the larger-scale vortical flow due to RichtmyerMeshkov instability while the flame propagation is due to the chemical reaction.The mixing properties of unburned/burned gases on both sides of the flame are quantitatively analyzed by using integral and statistical diagnostics.The results show that the large-scale flow due to the vortical motion always plays a dominating role during the reactive interaction process;however,the effect of chemistry becomes more important at the later stage of the interactions,especially for higher Mach number cases.The scalar dissipation due to the molecular diffusion is always small in the present study and can be negligible.展开更多
Gaseous detonation propagating in a toroidal chamber was numerically studied for hydrogen/oxygenitrogen mixtures. The numerical method used is based on the three-dimensional Euler equations with detailed finiterate ch...Gaseous detonation propagating in a toroidal chamber was numerically studied for hydrogen/oxygenitrogen mixtures. The numerical method used is based on the three-dimensional Euler equations with detailed finiterate chemistry. The results show that the calculated streak picture is in qualitative agreement with the picture recorded by a high speed streak camera from published literature. The three-dimensional flow field induced by a continuously rotating detonation was visualized and distinctive features of the rotating detonations were clearly depicted. Owing to the unconfined character of detonation wavelet, a deficit of detonation parameters was observed. Due to the effects of wall geometries, the strength of the outside detonation front is stronger than that of the inside portion. The detonation thus propagates with a constant circular velocity. Numerical simulation also shows three-dimensional rotating detonation structures, which display specific feature of the detonation- shock combined wave. Discrete burning gas pockets are formed due to instability of the discontinuity. It is believed that the present study could give an insight into the interest- ing properties of the continuously rotating detonation, and is thus beneficial to the design of continuous detonation propulsion systems.展开更多
Observations are presented from experiments and calculations where a laminar spherical CH4/air flame is perturbed successively by incident and reflected shock waves. The experiments are performed in a standard shock t...Observations are presented from experiments and calculations where a laminar spherical CH4/air flame is perturbed successively by incident and reflected shock waves. The experiments are performed in a standard shock tube arrangement, in which a high-speed shadowgraph imaging system is used to record evolutions of the flame. Numerical simulations are conducted by using second-order wave propagation algorithms, based on two-dimensional axisymmetric Navier-Stokes equations with detailed chemical reactions. Qualitative agreements are obtained between the experimental and numerical results. Under actions of incident shock waves, Richtmyer-Meshkov instability responsible for the flame deformation is induced in the flame, and the distoned flame takes a barrel shape. Then, under subsequent actions of the shock wave reflected from a planar wall, the flame takes an inclined non-symmetrical kidney shape in a symmetric cross section, which means a mushroom-like shape of the flame comes finally into being. The vorticity direction in the ring cap has been altered by the reflected shock's action, which makes the head of the mushroom-like flame extend quickly to the side wall.展开更多
The instabilities of a three-dimensional sinusoidally premixed flame induced by an incident shock wave with Mach = 1.7 and its reshock waves were studied by using the Navier-Stokes (NS) equations with a single-step ch...The instabilities of a three-dimensional sinusoidally premixed flame induced by an incident shock wave with Mach = 1.7 and its reshock waves were studied by using the Navier-Stokes (NS) equations with a single-step chemical reaction and a high resolution, 9th-order weighted essentially non-oscillatory scheme. The computational results were validated by the grid independence test and the experimental results in the literature. The computational results show that after the passage of incident shock wave the flame interface develops in symmetric structure accompanied by large-scale transverse vortex structures. After the interactions by successive reshock waves, the flame interface is gradually destabilized and broken up, and the large-scale vortex structures are gradually transformed into small-scale vortex structures. The small-scale vortices tend to be isotropic later. The results also reveal that the evolution of the flame interface is affected by both mixing process and chemical reaction. In order to identify the relationship between the mixing and the chemical reaction, a dimensionless parameter, , that is defined as the ratio of mixing time scale to chemical reaction time scale, is introduced. It is found that at each interaction stage the effect of chemical reaction is enhanced with time. The enhanced effect of chemical reaction at the interaction stage by incident shock wave is greater than that at the interaction stages by reshock waves. The result suggests that the parameter can reasonably character the features of flame interface development induced by the multiple shock waves.展开更多
Based on a weighted essentially nonoscillatory scheme, the multiple interactions of a flame interface with an incident shock wave and its reshock waves are numerically simulated by solving the compressible reactive Na...Based on a weighted essentially nonoscillatory scheme, the multiple interactions of a flame interface with an incident shock wave and its reshock waves are numerically simulated by solving the compressible reactive Navier–Stokes equations with a single-step Arrhenius chemical reaction. The two-dimensional sinusoidally perturbed premixed flames with different initial perturbed amplitudes are used to investigate the effect of the initial perturbation on the flame evolutions. The results show that the development of the flame interface is directly affected by the initial perturbed amplitudes before the passages of reshock waves, and the perturbation development is mainly controlled by the Richtmyer–Meshkov instability(RMI). After the successive impacts of multiple reshock waves, the chemical reaction accelerates the consumption of reactants and leads to a gradual disappearance of the initial perturbed information. The perturbation developments in frozen flows with the same initial interface as those in reactive flows are also demonstrated.Comparisons of results between the reactive and frozen flows show that a chemical reaction changes the perturbation pattern of the flame interface by decreasing the density gradient,thereby weakening the baroclinic torque in the flame mixing region, and therefore plays a dominant role after the passage of reshock waves.展开更多
The coronavirus disease 2019(COVID-19)pandemic caused by frequently mutating severe acute respiratory syndrome coronavirus 2(SARS-CoV-2)has had a worldwide impact.However,detailed data on the potential aerosol transmi...The coronavirus disease 2019(COVID-19)pandemic caused by frequently mutating severe acute respiratory syndrome coronavirus 2(SARS-CoV-2)has had a worldwide impact.However,detailed data on the potential aerosol transmission of SARS-CoV-2 in real-world and controlled laboratory settings remain sparse.During the COVID-19 pandemic in Shanghai,China in 2022,samples were collected in a Fangcang shelter hospital,a large-scale temporary hospital rapidly built by converting the existing National Exhibition and Convention Center(Shanghai)into a health care facility.Aerosol samples at different sites and intervals around patients and in public areas,surface samples,and pharyngeal swab samples from corresponding patients were included.Samples were tested for SARS-CoV-2 using real-time quantitative polymerase chain reaction(RT-qPCR)assays,followed by sequencing if the cycle threshold(Ct)value was<30.The positivity rate for SARS-CoV-2 in aerosol samples was high in contaminated zones(37.5%,104/277),especially around the bed(41.2%,68/165)and near ventilation inlets(45.2%,14/31).The prevalence of SARS-CoV-2 around the bed,public areas,and air inlets of exhaust vents fluctuated and was closely related to the positivity rate among patients at corresponding sampling sites.Some surface samples of different personal protective equipment from medical staff had high positivity rates.Sixty sequences of joined ORF1ab and spike genes obtained from sixty samples represented two main clusters of Omicron SARS-CoV-2.There was consistency in virus sequences from the same patient and their environment,and the detected virus sequences matched those of virus strains in circulation during the collection periods,which indicated a high likelihood of cross-contamination in the Fangcang shelter hospital.In summary,the results provide a quantitative and real landscape of the aerosol transmission of SARS-CoV-2 and a patient-centered view of contamination in large and enclosed spaces and offer a useful guide for taking targeted measures to avoid nosocomial infections during the management of SARS-CoV-2 or other respiratory virus diseases in a Fangcang shelter hospital.展开更多
This work is motivated by previous experimental and numerical studies which reveal that the hairpin vortex could be formed by the interaction between spanwise adjacent low-speed streaks. To prove that such an interact...This work is motivated by previous experimental and numerical studies which reveal that the hairpin vortex could be formed by the interaction between spanwise adjacent low-speed streaks. To prove that such an interaction mechanism is still applicable in the normal direction, two sinuous low-speed streaks with the same streamwise phase are set to be in the upper half and bottom half of a small size channel, respectively, and their evolution and interaction are investigated by direct numerical simulation. A new kind of hairpin-like vortical structure, distributed in the normal direction and straddled across both halves of the channel, is found during the cross-interaction process of the low-speed streaks. The influence of such a normal-distributed hairpin-like vortex (NHV) on the turbulent statistical regularity is also revealed. It is observed that the NHV can lead to a sudden surge of wall skin friction, but the value of the normal velocity as well as the streamwise and spanwise vorticity sharply decrease to zero in the center of the channel.展开更多
Standing of an Oblique Detonation Wave(ODW)on a wedge within combustor is the prerequisite of thrust generation for ODW engine which is regarded as a novel and conceptual propulsion device with hypersonic flight Mach ...Standing of an Oblique Detonation Wave(ODW)on a wedge within combustor is the prerequisite of thrust generation for ODW engine which is regarded as a novel and conceptual propulsion device with hypersonic flight Mach number.Usually a standing window of ODW is defined as the wedge angle ranged from the ODW detached angle from wedge(upper limit)to the angle that a Chapman-Jouguet(CJ)detonation occurs(lower limit).For pathological detonation cases,however,the CJ detonation cannot be achieved,and thus the lower limit of the standing window of ODW should be revisited.In present study,two types of reactions in hypersonic incoming flow that include the behavior of pathological detonation,that is,the single-step irreversible reaction with mole variation and the two-step irreversible reactions with exothermic process followed by endothermic process,have been used for studying standing behavior of ODW.The steady detonation polar analysis of ODW is carried out for both reaction systems.The results reveal that the reaction with more mole decrement and the reactions with stronger endothermic process show the pathological detonation feature and therefore modify the lower limit of standing window of ODW.Three equivalent parameters are proposed to quantitatively measure the standing window range of ODW from points of view of thermodynamics,Mach number of incoming flow and heat effect of reactions.It is found that the standing window of ODW is determined by the specific heat ratio,the overdrive degree of detonation and the endothermic level of the hypersonic incoming flow,regardless of whether the detonation is pathological or not.展开更多
Given the clean and inexhaustible solar energy from solar light,photocatalytic ammonia synthesis is extremely appealing.However,high electron-hole recombination rates and insufficient active sites severely limited N2 ...Given the clean and inexhaustible solar energy from solar light,photocatalytic ammonia synthesis is extremely appealing.However,high electron-hole recombination rates and insufficient active sites severely limited N2 photoreaction reduction.Herein,we designed and fabricated FeCu alloy nanoparticles anchored on carbon nitride nanosheets with excellent photocatalyt ic ammonia synthesis performance.As a coupler between Fe and carbon nitride,Cu promotes the separation of photogenerated charge carriers in carbon nitride under solar light irradiation,and renters the semiconductor a forceful electron donor for the Fe active sites.The accumulated electrons at Fe sites furtherly facilitated the adsorption and activation of the molecular nitrogen.Besides,the uniform dispersed FeCu alloy nanoparticles were on carbon nitride nanosheets enhanced the stability of photocataly tic nitrogen reduction reaction,making the artificial photocataly tic ammonia synthesis more sustainable for application.This work highlights that a direct electron transfer channel can be used to regulate the photochemical nitrogen fixation network.展开更多
The interactions of a spherical flame with an incident shock wave and its reflected shock wave in a confined space were investigated using the three-dimensional reactive Navier-Stokes equations, with emphasis placed o...The interactions of a spherical flame with an incident shock wave and its reflected shock wave in a confined space were investigated using the three-dimensional reactive Navier-Stokes equations, with emphasis placed on the effect of chemical reactivity of mixture on the flame distortion and detonation initiation after the passage of the reflected shock wave. It is shown that the spatio-temporal characteristics of detonation initiation depend highly on the chemi- cal reactivity of the mixture. When the chemical reactivity enhances, the flame can be severely distorted to form a reactive shock bifurcation structure with detonations initiating at different three-dimensional spatial locations. Moreover, the detonation initiation would occur earlier in a mixture of more enhanced reactivity. The results reveal that the detona- tions arise from hot spots in the unburned region which are initiated by the shock-detonation-transition mechanism.展开更多
The evolution of two spanwise-aligned low-speed streaks in a wall turbulent flow, triggered by the instability of the subharmonic varicose (SV) mode, is studied by a direct numerical simulation (DNS) method in a s...The evolution of two spanwise-aligned low-speed streaks in a wall turbulent flow, triggered by the instability of the subharmonic varicose (SV) mode, is studied by a direct numerical simulation (DNS) method in a small spatial-periodic channel. The results show that the SV low-speed streaks are self-sustained at the early stage, and then transform into subharmonic sinuous (SS) low-speed streaks. Initially, the streamwise vortex sheets are formed by shearing, and then evolve into zigzag vortex sheets due to the mutual induction. As the intensification of the SV low-speed streaks becomes prominent, the tilted streamwise vortex tubes and the V-like streamwise vortex tubes can be formed simultaneously by increasing +~. When the SV low-speed streaks break down, new zigzag streamwise vortices will be generated, thus giving birth to the next sustaining cycle of the SV low-speed streaks. When the second breakdown happens, new secondary V-like streamwise vortices instead of zigzag streamwise vortices will be generated. Because of the sweep motion of the fluid induced by the secondary V-like streamwise vortices, each decayed low-speed streak can be divided into two parts, and each part combines with the part of another streak, finally leading to the formation of SS low-speed streaks.展开更多
The growth of mixing zone on an interface induced by Richtmyer-Meshkov(RM)instability occurs frequently in natural phenomena and in engineering applications.Usually,the medium on which the RM instability happens is in...The growth of mixing zone on an interface induced by Richtmyer-Meshkov(RM)instability occurs frequently in natural phenomena and in engineering applications.Usually,the medium on which the RM instability happens is inhomogeneous,the effect of medium inhomogeneity on the growth of the mixing zone during the RM instability is still not clear.Therefore,it is necessary to investigate the RM instability in inhomogeneous medium.Based on a high-order computational scheme,the interactions of a density interface with an incident shock wave(ISW)in inhomogeneous medium are numerically simulated by solving the compressible Navier-Stokes equations.The effect of the inhomogeneity on the interface evolution after the passage of ISW through the interface is investigated.The results show that the interface morphology develops in a distinctive "spike-spike"structure in inhomogeneous medium.Particularly,the spike structure on the bottom of the interface is due to the reverse induction of RM instability by curved ISW or reflected shock wave.With the increase of inhomogeneity,the growth rate of the mixing zone width on interface increases,and the wave patterns caused by interaction between the shock wave and interface are more complex.Compared with RM instability in homogeneous medium,the inhomogeneous distribution of the density in medium further enhances the baroclinic effect and induces larger vorticity in flow field.Therefore,the interface is stretched much more significantly under the induction of enhanced vorticity in inhomogeneous medium.Based on above analyses,a model for predicting the growth of mixing zone width on the interface after the passage of ISW is proposed,in order to provide a useful method for evaluations of perturbation growth behavior during the RM instability in inhomogeneous medium.展开更多
Electrochemical CO2reduction to C2H4can provide a sustainable route to reduce globally accelerating CO2emissions and produce energy-rich chemical feedstocks.However,the poor selectivity in C2H4electr...Electrochemical CO2reduction to C2H4can provide a sustainable route to reduce globally accelerating CO2emissions and produce energy-rich chemical feedstocks.However,the poor selectivity in C2H4electrosynthesis limits its implementation in industrially interesting processes.Herein,we report a composite structured catalyst composed of Ag and Cu2O with different crystal faces to achieve highly efficient reduction of CO2to C2H4.The catalyst composed of Ag and octahedral Cu2O enclosed with(111)facet exhibits the best CO2electroreduction performance,with the Faradaic efficiency(FE)and partial current density reaching 66.8%and 17.8 mA cm2 for C2H4product at-1.2 VRHE in 0.5 M KHCO3,respectively.Physical characterization and electrochemical test analysis indicate that the high selectivity for C2H4product stems from the synergistic effect of crystal faces control engineering and tandem catalysis.Specifically,Ag can provide optimal availability of CO intermediate by suppressing hydrogen evolution;subsequently,C-C coupling is promoted on the intimate surface of Cu2O with facetdependent selectivity.The insights gained from this work may be beneficial for designing efficient multicomponent catalysts for improving the selectivity of electrochemical CO2reduction reaction to generate C2þproducts.展开更多
This study develops a direct optimal growth algorithm for three-dimensional transient growth analysis of perturbations in channel flows which are globally stable but locally unstable. Different from traditional non-mo...This study develops a direct optimal growth algorithm for three-dimensional transient growth analysis of perturbations in channel flows which are globally stable but locally unstable. Different from traditional non-modal methods based on the Orr- Somrnerfeld and Squire (OSS) equations that assume simple base flows, this algorithm can be applied to arbitrarily complex base flows. In the proposed algorithm, a re- orthogonalization Arnoldi method is used to improve orthogonality of the orthogonal basis of the Krylov subspace generated by solving the linearized forward and adjoint Navier-Stokes (N-S) equations. The linearized adjoint N-S equations with the specific boundary conditions for the channel are derived, and a new convergence criterion is pro- posed. The algorithm is then applied to a one-dimensional base flow (the plane Poiseuille flow) and a two-dimensional base flow (the plane Poiseuille flow with a low-speed streak) in a channel. For one-dimensional cases, the effects of the spanwise width of the chan- nel and the Reynolds number on the transient growth of perturbations are studied. For two-dimensional cases, the effect of strength of initial low-speed streak is discussed. The presence of the streak in the plane Poiseuille flow leads to a larger and quicker growth of the perturbations than that in the one-dimensional case. For both cases, the results show that an optimal flow field leading to the largest growth of perturbations is character- ized by high- and low-speed streaks and the corresponding streamwise vortical structures. The lift-up mechanism that induces the transient growth of perturbations is discussed. The performance of the re-orthogonalization Arnoldi technique in the algorithm for both one- and two-dimensional base flows is demonstrated, and the algorithm is validated by comparing the results with those obtained from the OSS equations method and the cross- check method.展开更多
基金funded by the Geological Survey Program of the China Geological Survey(DD20230402,DD2026032030)the National Natural Science Foundation of China(42072069)。
摘要The enrichment and accumulation of natural gas hydrates depend on sufficient gas supply and effective migration pathways.The upward migration of deep thermogenic gases through fault systems is critical for seepage-type hydrate formation.This study aims to elucidate the developmental characteristics of Cenozoic fault systems in the eastern offshore area of Dongsha Island and their influence on natural gas hydrate formation.Utilizing high-resolution 3D seismic data,this study conducted a detailed structural interpretation and seismic attribute analysis to systematically investigate the spatial distribution,developmental stages,and dynamic mechanisms of the Cenozoic fault systems in this region.In addition,this study explored the role of these fault systems in facilitating the migration of deep thermogenic gases to the shallow strata.The study area is dominated by extensional and transtensional normal faults characterized by inherited development and relatively small fault displacements.The Cenozoic strata exhibit a tectonic framework of block-faulted uplift and subsidence with alternating highs and lows.Faults on either side of the central uplift dip in opposite directions and commonly exhibit parallel,step-like patterns.Differences in fault system attitudes were observed between the southern and northern parts of the study area.In the south,fault strikes remained consistent from deep to shallow levels,predominantly trending NE and NEE.In the north,fault strikes varied significantly with depth,transitioning from predominantly NEE in deeper strata to EW and NWW in shallower strata.The study identifies two distinct phases of Cenozoic fault activity:(1)66–10 Ma,a regional extensional tectonic regime controlled fault development,resulting in the formation of NEE-trending normal faults;(2)10–2.6 Ma,the Dongsha Movement influenced fault activity,during which EW-and NW–W-trending transtensional faults with dextral strike-slip characteristics developed in the Miocene strata of the northern region.The Cenozoic fault system played a significant positive role in facilitating the migration of deep thermogenic gas to shallow levels,thereby enabling the formation of natural gas hydrates.
基金supports of the National Natural Science Foundation of China(Grant No.12372195)the Anhui Provincial Natural Science Foundation(Grant No.2408085J007)+1 种基金the Dreams Foundation of Jianghuai Advance Technology Center(Grant No.2023-ZM01 X013)BIM Engineering Center of Anhui Province(No.AHBIM2022KF02)are greatly appreciated.
摘要The sampling method is an important numerical technique for solving reliability problems in engineering systems.However,the evaluation of the failure probability using classical sampling methods is time-consuming for complex engineering structure.To address this issue,this paper proposes a gradient optimization assisted bubble sampling method(GOBSM)to reduce the computational costs,which enhances the coverage range of bubbles,thereby improving the computational efficiency without sacrificing the accuracy.Furthermore,the bubble gradient iterative algorithm is developed to efficiently construct bubbles.Eight complex numerical examples are tested for assessing the failure probability,and the results demonstrate the performance of GOBSM.
基金supported by grants from the Natural Science Foundation Key Project of Henan Province(232300421124)Funding for Scientific Research and Innovation Team of the First Affiliated Hospital of Zhengzhou University(QNCXTD2023002)the National Natural Science Foundation of China(82470654).
摘要To the Editor:Local ablation is crucial in the treatment of hepatocellular car-cinoma(HCC)and is widely regarded as a good alternative to surgical intervention for achieving local control of HCC.However,the most commonly used methods for treating HCC are radiofre-quency ablation(RFA)and microwave ablation(MWA),which are thermal ablation methods with limited applicability for tumors in special locations.For example,thermal damage to adjacent impor-tant structures(such as blood vessels,bile ducts and intestines)leads to serious complications.Moreover,thermal ablation of tu-mors near large blood vessels is associated with a higher incidence of incomplete eradication(i.e.,the heat-sink effect)[1,2].
基金supported by the National Natural Science Foundation of China(Grant No.12472335).
摘要The stability of gaseous detonation waves is crucial for the operation of detonation-based propulsion systems and the assessment of industrial explosion hazards.However,research on the stability of detonation waves in complex reactive systems that are composed of actual fuels and oxidants and can be described by numerous elementary chemical reactions,has not been fully carried out.To investigate the relationship between linear and nonlinear stabilities in gaseous detonation wave propagation for complex reactive systems,the linear stability analysis and the one-dimensionally nonlinear numerical simulations of H2/O2/Ar(argon)detonations based on the reactive Euler equations and detailed reaction mechanisms are carried out.The results show that in complex reactive systems characterized by elementary chemical reactions,the results of linear stability computation of detonation are consistent with those from one-dimensionally nonlinear oscillations of detonation wave.Utilizing these linear stability results,a neutral stability curve and a perturbation frequency transition curve in the phase plane of initial pressure versus inert gas(Ar)dilution ratio are derived,especially the new frequency transition curve clearly describes the transition of perturbations from low-frequency to high-frequency mode.One-dimensional nonlinear simulations show that near the perturbation frequency transition curve,the oscillations of the detonation wave can also transform between the lowfrequency,high-amplitude oscillation mode and the high-frequency,low-amplitude oscillation mode,with the oscillation frequency corresponding to the mode that exhibits the maximum growth rate identified in the linear stability analysis.This investigation into detonation stability in complex reactive gases offers guidance for selecting appropriate initial conditions and gas compositions in practical applications of detonation.
基金the National Natural Science Foundation of China(Grant 51609115)the Foundation of National Key Laboratory of Transient Physics(Grant9140C300206150C30143).
摘要The effect of micro-bubbles on the turbulent boundary layer in the channel flow with Reynolds numbers (Re) ranging from 0.87 × 105 to 1.23 × 105 is experimentally studied by using particle image velocimetry (PIV) measurements.The microbubbles are produced by water electrolysis.The velocity profiles,Reynolds stress and instantaneous structures of the boundary layer,with and without micro-bubbles,are measured and analyzed.The presence of micro-bubbles changes the streamwise mean velocity of the fluid and increases the wall shear stress.The results show that micro-bubbles have two effects,buoyancy and extrusion,which dominate the flow behavior of the mixed fluid in the turbulent boundary layer.The buoyancy effect leads to upward motion that drives the fluid motion in the same direction and,therefore,enhances the turbulence intense of the boundary layer.While for the extrusion effect,the presence of accumulated micro-bubbles pushes the flow structures in the turbulent boundary layer away from the near-wall region.The interaction between these two effects causes the vorticity structures and turbulence activity to be in the region far away from the wall.The buoyancy effect is dominant when the Re is relatively small,while the extrusion effect plays a more important role when Re rises.
基金supported by the National Natural Science Foundation of China (10872096)the Open Fund of State Key Laboratory of Explosion Science and Technology,Beijing Institute of Technology (KFJJ09-13)
摘要An oblique detonation wave for a Mach 7 inlet flow over a long enough wedge of 30 turning angle is simulated numerically using Euler equation and one-step rection model.The fifth-order WENO scheme is adopted to capture the shock wave.The numerical results show that with the compression of the wedge wall the detonation wave front structure is divided into three sections:the ZND model-like strcuture,single-sided triple point structure and dual-headed triple point strucuture.The first structure is the smooth straight,and the second has the characteristic of the triple points propagating dowanstream only with the same velocity,while the dual-headed triple point structure is very complicated.The detonation waves facing upstream and downstream propagate with different velocities,in which the periodic collisions of the triple points cause the oscillation of the detonation wave front.This oscillation process has temporal and spatial periodicity.In addition,the triple point trace are recorded to obtain different cell structures in three sections.
基金The work was supported by the National Natural Science Foundation of China(11372140).
摘要The three-dimensional interactions of a perturbed premixed flame interface with a planar incident shock wave and its reflected shock waves are numerically simulated by solving the compressible,reactive Navier-Stokes equations with the high-resolution scheme and a single-step chemical reaction.The effects of the initial incident shock wave strength (Mach number) and the initial perturbation pattern of interface on the interactions are investigated.The distinct properties of perturbation growth on the flame interface during the interactions are presented.Our results show that perturbation growth is mainly attributed to the flame stretching and propagation.The flame stretching is associated with the larger-scale vortical flow due to RichtmyerMeshkov instability while the flame propagation is due to the chemical reaction.The mixing properties of unburned/burned gases on both sides of the flame are quantitatively analyzed by using integral and statistical diagnostics.The results show that the large-scale flow due to the vortical motion always plays a dominating role during the reactive interaction process;however,the effect of chemistry becomes more important at the later stage of the interactions,especially for higher Mach number cases.The scalar dissipation due to the molecular diffusion is always small in the present study and can be negligible.
基金supported by the National Natural Science Foundation of China (10872096)the Open Fund of State Key Laboratory of Explosion Science and Technology, Beijing University of Science and Technology (KFJJ09-13)
摘要Gaseous detonation propagating in a toroidal chamber was numerically studied for hydrogen/oxygenitrogen mixtures. The numerical method used is based on the three-dimensional Euler equations with detailed finiterate chemistry. The results show that the calculated streak picture is in qualitative agreement with the picture recorded by a high speed streak camera from published literature. The three-dimensional flow field induced by a continuously rotating detonation was visualized and distinctive features of the rotating detonations were clearly depicted. Owing to the unconfined character of detonation wavelet, a deficit of detonation parameters was observed. Due to the effects of wall geometries, the strength of the outside detonation front is stronger than that of the inside portion. The detonation thus propagates with a constant circular velocity. Numerical simulation also shows three-dimensional rotating detonation structures, which display specific feature of the detonation- shock combined wave. Discrete burning gas pockets are formed due to instability of the discontinuity. It is believed that the present study could give an insight into the interest- ing properties of the continuously rotating detonation, and is thus beneficial to the design of continuous detonation propulsion systems.
基金supported by the National Natural Science Foundation of China (10472047)the Open Fund of State Key Laboratory of Explosion Science Technology, Beijing University of Science and Technology (KFJJ06-3)
摘要Observations are presented from experiments and calculations where a laminar spherical CH4/air flame is perturbed successively by incident and reflected shock waves. The experiments are performed in a standard shock tube arrangement, in which a high-speed shadowgraph imaging system is used to record evolutions of the flame. Numerical simulations are conducted by using second-order wave propagation algorithms, based on two-dimensional axisymmetric Navier-Stokes equations with detailed chemical reactions. Qualitative agreements are obtained between the experimental and numerical results. Under actions of incident shock waves, Richtmyer-Meshkov instability responsible for the flame deformation is induced in the flame, and the distoned flame takes a barrel shape. Then, under subsequent actions of the shock wave reflected from a planar wall, the flame takes an inclined non-symmetrical kidney shape in a symmetric cross section, which means a mushroom-like shape of the flame comes finally into being. The vorticity direction in the ring cap has been altered by the reflected shock's action, which makes the head of the mushroom-like flame extend quickly to the side wall.
基金supported by the National Natural Science Foundation of China (Grant 11372140)
摘要The instabilities of a three-dimensional sinusoidally premixed flame induced by an incident shock wave with Mach = 1.7 and its reshock waves were studied by using the Navier-Stokes (NS) equations with a single-step chemical reaction and a high resolution, 9th-order weighted essentially non-oscillatory scheme. The computational results were validated by the grid independence test and the experimental results in the literature. The computational results show that after the passage of incident shock wave the flame interface develops in symmetric structure accompanied by large-scale transverse vortex structures. After the interactions by successive reshock waves, the flame interface is gradually destabilized and broken up, and the large-scale vortex structures are gradually transformed into small-scale vortex structures. The small-scale vortices tend to be isotropic later. The results also reveal that the evolution of the flame interface is affected by both mixing process and chemical reaction. In order to identify the relationship between the mixing and the chemical reaction, a dimensionless parameter, , that is defined as the ratio of mixing time scale to chemical reaction time scale, is introduced. It is found that at each interaction stage the effect of chemical reaction is enhanced with time. The enhanced effect of chemical reaction at the interaction stage by incident shock wave is greater than that at the interaction stages by reshock waves. The result suggests that the parameter can reasonably character the features of flame interface development induced by the multiple shock waves.
基金supported by the National Natural Science Foundation of China(Grant 11372140)
摘要Based on a weighted essentially nonoscillatory scheme, the multiple interactions of a flame interface with an incident shock wave and its reshock waves are numerically simulated by solving the compressible reactive Navier–Stokes equations with a single-step Arrhenius chemical reaction. The two-dimensional sinusoidally perturbed premixed flames with different initial perturbed amplitudes are used to investigate the effect of the initial perturbation on the flame evolutions. The results show that the development of the flame interface is directly affected by the initial perturbed amplitudes before the passages of reshock waves, and the perturbation development is mainly controlled by the Richtmyer–Meshkov instability(RMI). After the successive impacts of multiple reshock waves, the chemical reaction accelerates the consumption of reactants and leads to a gradual disappearance of the initial perturbed information. The perturbation developments in frozen flows with the same initial interface as those in reactive flows are also demonstrated.Comparisons of results between the reactive and frozen flows show that a chemical reaction changes the perturbation pattern of the flame interface by decreasing the density gradient,thereby weakening the baroclinic torque in the flame mixing region, and therefore plays a dominant role after the passage of reshock waves.
基金supported by National Key Research and Development Program of China(2019YFC1200500)supported by Emergency Key Program of Guangzhou Laboratory(EKPG21-01)supported by National Natural Science Foundation of China(32070166)。
摘要The coronavirus disease 2019(COVID-19)pandemic caused by frequently mutating severe acute respiratory syndrome coronavirus 2(SARS-CoV-2)has had a worldwide impact.However,detailed data on the potential aerosol transmission of SARS-CoV-2 in real-world and controlled laboratory settings remain sparse.During the COVID-19 pandemic in Shanghai,China in 2022,samples were collected in a Fangcang shelter hospital,a large-scale temporary hospital rapidly built by converting the existing National Exhibition and Convention Center(Shanghai)into a health care facility.Aerosol samples at different sites and intervals around patients and in public areas,surface samples,and pharyngeal swab samples from corresponding patients were included.Samples were tested for SARS-CoV-2 using real-time quantitative polymerase chain reaction(RT-qPCR)assays,followed by sequencing if the cycle threshold(Ct)value was<30.The positivity rate for SARS-CoV-2 in aerosol samples was high in contaminated zones(37.5%,104/277),especially around the bed(41.2%,68/165)and near ventilation inlets(45.2%,14/31).The prevalence of SARS-CoV-2 around the bed,public areas,and air inlets of exhaust vents fluctuated and was closely related to the positivity rate among patients at corresponding sampling sites.Some surface samples of different personal protective equipment from medical staff had high positivity rates.Sixty sequences of joined ORF1ab and spike genes obtained from sixty samples represented two main clusters of Omicron SARS-CoV-2.There was consistency in virus sequences from the same patient and their environment,and the detected virus sequences matched those of virus strains in circulation during the collection periods,which indicated a high likelihood of cross-contamination in the Fangcang shelter hospital.In summary,the results provide a quantitative and real landscape of the aerosol transmission of SARS-CoV-2 and a patient-centered view of contamination in large and enclosed spaces and offer a useful guide for taking targeted measures to avoid nosocomial infections during the management of SARS-CoV-2 or other respiratory virus diseases in a Fangcang shelter hospital.
基金Project supported by the Natural Science Foundation of Zhejiang Province,China(Grant No.LQ18A020002)the National Natural Science Foundation of China(Grant No.41776105)the Foundation of the Educational Committee of Zhejiang Province,China(Grant No.Y201737053)
摘要This work is motivated by previous experimental and numerical studies which reveal that the hairpin vortex could be formed by the interaction between spanwise adjacent low-speed streaks. To prove that such an interaction mechanism is still applicable in the normal direction, two sinuous low-speed streaks with the same streamwise phase are set to be in the upper half and bottom half of a small size channel, respectively, and their evolution and interaction are investigated by direct numerical simulation. A new kind of hairpin-like vortical structure, distributed in the normal direction and straddled across both halves of the channel, is found during the cross-interaction process of the low-speed streaks. The influence of such a normal-distributed hairpin-like vortex (NHV) on the turbulent statistical regularity is also revealed. It is observed that the NHV can lead to a sudden surge of wall skin friction, but the value of the normal velocity as well as the streamwise and spanwise vorticity sharply decrease to zero in the center of the channel.
基金co-supported by the National Natural Science Foundation of China(No.11872213)the Open Fund of Key Laboratory of Transient Physics,NJUST(No.6142604180205)the Opening Project of State Key Laboratory of Explosion Science and Technology(Beijing Institute of Technology)(No.KFJJ17-11M)。
摘要Standing of an Oblique Detonation Wave(ODW)on a wedge within combustor is the prerequisite of thrust generation for ODW engine which is regarded as a novel and conceptual propulsion device with hypersonic flight Mach number.Usually a standing window of ODW is defined as the wedge angle ranged from the ODW detached angle from wedge(upper limit)to the angle that a Chapman-Jouguet(CJ)detonation occurs(lower limit).For pathological detonation cases,however,the CJ detonation cannot be achieved,and thus the lower limit of the standing window of ODW should be revisited.In present study,two types of reactions in hypersonic incoming flow that include the behavior of pathological detonation,that is,the single-step irreversible reaction with mole variation and the two-step irreversible reactions with exothermic process followed by endothermic process,have been used for studying standing behavior of ODW.The steady detonation polar analysis of ODW is carried out for both reaction systems.The results reveal that the reaction with more mole decrement and the reactions with stronger endothermic process show the pathological detonation feature and therefore modify the lower limit of standing window of ODW.Three equivalent parameters are proposed to quantitatively measure the standing window range of ODW from points of view of thermodynamics,Mach number of incoming flow and heat effect of reactions.It is found that the standing window of ODW is determined by the specific heat ratio,the overdrive degree of detonation and the endothermic level of the hypersonic incoming flow,regardless of whether the detonation is pathological or not.
基金financially supported by the National Natural Science Foundation of China(Nos.52002361 and 22109120)the Science and Technology Research Project of Jiangxi Provincial Education Department(Nos.GJJ2201045 and GJJ2201007)+2 种基金Zhejiang Provincial Natural Science Foundation of China(No.LQ21B030002)the Key Project for Science and Technology Cooperation of Jiangxi Province(No.20212BDH80005)the Project of the Science and Technology of Jingdezhen City(No.20202GYZD013-16)。
摘要Given the clean and inexhaustible solar energy from solar light,photocatalytic ammonia synthesis is extremely appealing.However,high electron-hole recombination rates and insufficient active sites severely limited N2 photoreaction reduction.Herein,we designed and fabricated FeCu alloy nanoparticles anchored on carbon nitride nanosheets with excellent photocatalyt ic ammonia synthesis performance.As a coupler between Fe and carbon nitride,Cu promotes the separation of photogenerated charge carriers in carbon nitride under solar light irradiation,and renters the semiconductor a forceful electron donor for the Fe active sites.The accumulated electrons at Fe sites furtherly facilitated the adsorption and activation of the molecular nitrogen.Besides,the uniform dispersed FeCu alloy nanoparticles were on carbon nitride nanosheets enhanced the stability of photocataly tic nitrogen reduction reaction,making the artificial photocataly tic ammonia synthesis more sustainable for application.This work highlights that a direct electron transfer channel can be used to regulate the photochemical nitrogen fixation network.
基金supported by the National Natural Science Foundation of China (10972107)Open Fund of State Key Laboratory of Explosion Science and Technology, Beijing Institute of Technology (KFJJ12-4Y)Jiangsu Innovation Program for Graduate Education (CXLX11 0271)
摘要The interactions of a spherical flame with an incident shock wave and its reflected shock wave in a confined space were investigated using the three-dimensional reactive Navier-Stokes equations, with emphasis placed on the effect of chemical reactivity of mixture on the flame distortion and detonation initiation after the passage of the reflected shock wave. It is shown that the spatio-temporal characteristics of detonation initiation depend highly on the chemi- cal reactivity of the mixture. When the chemical reactivity enhances, the flame can be severely distorted to form a reactive shock bifurcation structure with detonations initiating at different three-dimensional spatial locations. Moreover, the detonation initiation would occur earlier in a mixture of more enhanced reactivity. The results reveal that the detona- tions arise from hot spots in the unburned region which are initiated by the shock-detonation-transition mechanism.
基金supported by the National Natural Science Foundation of China(Nos.11372140 and11202102)the Innovation Project for College Graduates of Jiangsu Province(No.CXZZ13-0189)
摘要The evolution of two spanwise-aligned low-speed streaks in a wall turbulent flow, triggered by the instability of the subharmonic varicose (SV) mode, is studied by a direct numerical simulation (DNS) method in a small spatial-periodic channel. The results show that the SV low-speed streaks are self-sustained at the early stage, and then transform into subharmonic sinuous (SS) low-speed streaks. Initially, the streamwise vortex sheets are formed by shearing, and then evolve into zigzag vortex sheets due to the mutual induction. As the intensification of the SV low-speed streaks becomes prominent, the tilted streamwise vortex tubes and the V-like streamwise vortex tubes can be formed simultaneously by increasing +~. When the SV low-speed streaks break down, new zigzag streamwise vortices will be generated, thus giving birth to the next sustaining cycle of the SV low-speed streaks. When the second breakdown happens, new secondary V-like streamwise vortices instead of zigzag streamwise vortices will be generated. Because of the sweep motion of the fluid induced by the secondary V-like streamwise vortices, each decayed low-speed streak can be divided into two parts, and each part combines with the part of another streak, finally leading to the formation of SS low-speed streaks.
摘要The growth of mixing zone on an interface induced by Richtmyer-Meshkov(RM)instability occurs frequently in natural phenomena and in engineering applications.Usually,the medium on which the RM instability happens is inhomogeneous,the effect of medium inhomogeneity on the growth of the mixing zone during the RM instability is still not clear.Therefore,it is necessary to investigate the RM instability in inhomogeneous medium.Based on a high-order computational scheme,the interactions of a density interface with an incident shock wave(ISW)in inhomogeneous medium are numerically simulated by solving the compressible Navier-Stokes equations.The effect of the inhomogeneity on the interface evolution after the passage of ISW through the interface is investigated.The results show that the interface morphology develops in a distinctive "spike-spike"structure in inhomogeneous medium.Particularly,the spike structure on the bottom of the interface is due to the reverse induction of RM instability by curved ISW or reflected shock wave.With the increase of inhomogeneity,the growth rate of the mixing zone width on interface increases,and the wave patterns caused by interaction between the shock wave and interface are more complex.Compared with RM instability in homogeneous medium,the inhomogeneous distribution of the density in medium further enhances the baroclinic effect and induces larger vorticity in flow field.Therefore,the interface is stretched much more significantly under the induction of enhanced vorticity in inhomogeneous medium.Based on above analyses,a model for predicting the growth of mixing zone width on the interface after the passage of ISW is proposed,in order to provide a useful method for evaluations of perturbation growth behavior during the RM instability in inhomogeneous medium.
基金This work was supported by the University of Science and Technology Beijing.DG acknowledges the financial support from 111 Project(no.B170003)Foshan Science and Technology Innovation Project(no.2018IT100363).
摘要Electrochemical CO2reduction to C2H4can provide a sustainable route to reduce globally accelerating CO2emissions and produce energy-rich chemical feedstocks.However,the poor selectivity in C2H4electrosynthesis limits its implementation in industrially interesting processes.Herein,we report a composite structured catalyst composed of Ag and Cu2O with different crystal faces to achieve highly efficient reduction of CO2to C2H4.The catalyst composed of Ag and octahedral Cu2O enclosed with(111)facet exhibits the best CO2electroreduction performance,with the Faradaic efficiency(FE)and partial current density reaching 66.8%and 17.8 mA cm2 for C2H4product at-1.2 VRHE in 0.5 M KHCO3,respectively.Physical characterization and electrochemical test analysis indicate that the high selectivity for C2H4product stems from the synergistic effect of crystal faces control engineering and tandem catalysis.Specifically,Ag can provide optimal availability of CO intermediate by suppressing hydrogen evolution;subsequently,C-C coupling is promoted on the intimate surface of Cu2O with facetdependent selectivity.The insights gained from this work may be beneficial for designing efficient multicomponent catalysts for improving the selectivity of electrochemical CO2reduction reaction to generate C2þproducts.
基金supported by the National Natural Science Foundation of China(No.11372140)
摘要This study develops a direct optimal growth algorithm for three-dimensional transient growth analysis of perturbations in channel flows which are globally stable but locally unstable. Different from traditional non-modal methods based on the Orr- Somrnerfeld and Squire (OSS) equations that assume simple base flows, this algorithm can be applied to arbitrarily complex base flows. In the proposed algorithm, a re- orthogonalization Arnoldi method is used to improve orthogonality of the orthogonal basis of the Krylov subspace generated by solving the linearized forward and adjoint Navier-Stokes (N-S) equations. The linearized adjoint N-S equations with the specific boundary conditions for the channel are derived, and a new convergence criterion is pro- posed. The algorithm is then applied to a one-dimensional base flow (the plane Poiseuille flow) and a two-dimensional base flow (the plane Poiseuille flow with a low-speed streak) in a channel. For one-dimensional cases, the effects of the spanwise width of the chan- nel and the Reynolds number on the transient growth of perturbations are studied. For two-dimensional cases, the effect of strength of initial low-speed streak is discussed. The presence of the streak in the plane Poiseuille flow leads to a larger and quicker growth of the perturbations than that in the one-dimensional case. For both cases, the results show that an optimal flow field leading to the largest growth of perturbations is character- ized by high- and low-speed streaks and the corresponding streamwise vortical structures. The lift-up mechanism that induces the transient growth of perturbations is discussed. The performance of the re-orthogonalization Arnoldi technique in the algorithm for both one- and two-dimensional base flows is demonstrated, and the algorithm is validated by comparing the results with those obtained from the OSS equations method and the cross- check method.