The evolution of the three-dimensional time-developing mixing layer was simulated numerically using the pseudo-spectral method. The initial perturbations consisted of the two-dimensional fundamental wave and the stre...The evolution of the three-dimensional time-developing mixing layer was simulated numerically using the pseudo-spectral method. The initial perturbations consisted of the two-dimensional fundamental wave and the streamwise-invariant three-dimensional disturbance. A comparison of the formations of the streamwise vortices with different amplitude functions for three-dimensional disturbances was made. In one case the results are similar to that of Rogers and Moser (1992), whereas a different way in which the quadrupole forms and sudden expansion of the rib were observed in another case. The simulation also confirms that stretching by the forming roller rather than Rayleigh centrifugal instability is responsible for the formation of the rib. Finally, numerical flow visualization results were presented. (Edited author abstract) 9 Refs.展开更多
Discontinuous Galerkin(DG) method is known to have several advantages for flow simulations,in particular,in fiexible accuracy management and adaptability to mesh refinement. In the present work,the DG method is deve...Discontinuous Galerkin(DG) method is known to have several advantages for flow simulations,in particular,in fiexible accuracy management and adaptability to mesh refinement. In the present work,the DG method is developed for numerical simulations of both temporally and spatially developing mixing layers. For the temporally developing mixing layer,both the instantaneous fiow field and time evolution of momentum thickness agree very well with the previous results. Shocklets are observed at higher convective Mach numbers and the vortex paring manner is changed for high compressibility. For the spatially developing mixing layer,large-scale coherent structures and self-similar behavior for mean profiles are investigated. The instantaneous fiow field for a three-dimensional compressible mixing layer is also reported,which shows the development of largescale coherent structures in the streamwise direction. All numerical results suggest that the DG method is effective in performing accurate numerical simulations for compressible shear fiows.展开更多
In this paper, several mathematical models for the large scale structures in some special kinds of mixing layers, which might be practically useful for enhancing the mixing, are proposed. First, the linear growth rate...In this paper, several mathematical models for the large scale structures in some special kinds of mixing layers, which might be practically useful for enhancing the mixing, are proposed. First, the linear growth rate of the large scale structures in the mixing layers was calculated. Then, using the much improved weakly non-linear theory, combined with the energy method, the non-linear evolution of large scale structures in two special mixing layer configurations is calculated. One of the mixing lavers has equal magnitudes of the upstream velocity vectors, while the angles between the velocity vectors and the trailing edge were pi /2 - phi and pi /2 + phi, respectively. The other mixing layer was generated by a splitter-plate with a 45-degree-sweep trailing edge.展开更多
Numerical simulations have been performed in time-developing plane mixing layers of the viscoelastic second-order fluids with pseudo-spectral method. Roll-up, pairing and merging of large eddies were examined at high ...Numerical simulations have been performed in time-developing plane mixing layers of the viscoelastic second-order fluids with pseudo-spectral method. Roll-up, pairing and merging of large eddies were examined at high Reynolds numbers and low Deborah numbers. The effect of viscoelastics on the evolution of the large coherent structure was shown by making a comparison between the second-order and Newtonian fluids at the same Reynolds numbers.展开更多
Semi-periodic structures namely inclined wavy structures (IWS) are experimentally observed in compressible mixing layers at two convective Mach numbers (Mc = 0.11 and 0.47). Flow structures are visualized by the l...Semi-periodic structures namely inclined wavy structures (IWS) are experimentally observed in compressible mixing layers at two convective Mach numbers (Mc = 0.11 and 0.47). Flow structures are visualized by the laserinduced planar laser Mie scattering (PLMS) technique. Two methods are developed to investigate the spatial distribu- tion and geometry of IWS: (1) the dominant mode extrac- tion (DME) method, to extract the dominant modes of IWS from the streamwise gray-level fluctuation, and (2) the phase tracking (PT) method, to identify the shape of IWS. The re- sults suggest that pressure perturbations account for the for- marion of IWS in the initial mixing region and the joint effect of dilatation and coherent vortices enhances IWS in the well- developed region. The large transverse (cross-flow) scale of the IWS and their relation to coherent vortices (CV) indicate that the disturbance originated from CV in the mixing center propagates far into the free streams. The DME and the PT method are shown to be the effective tools to study the geometrical features of wavy structures in compressible shear flows.展开更多
The method of nonlinear parabolized stability equations(PSE) is applied in the simulation of vortex structures in compressible mixing layer.The spatially-evolving unstable waves,which dominate the vortex structure,a...The method of nonlinear parabolized stability equations(PSE) is applied in the simulation of vortex structures in compressible mixing layer.The spatially-evolving unstable waves,which dominate the vortex structure,are investigated through spatial marching method.The instantaneous flow field is obtained by adding the harmonic waves to basic flow.The results show that T-S waves do not keep growing exponentially as the linear evolution,the energy transfer to high order harmonic modes,and that finally all harmonic modes get saturated due to nonlinear interaction.The mean flow distortion induced by the nonlinear interaction between the harmonic modes and their conjugate harmonic ones,makes great change of the average flow and increases the thickness of mixing layer. PSE methods can well capture the two- and three-dimensional large scale nonlinear vortex structures in mixing layers such as vortex roll-up,vortex pairing,and A vortex.展开更多
The infrared imaging windows of the hyper/supersonic optical dome are encountering severe aero-optical effects[AOEs],so a flow control device,the ramp vortex generator array[RVGA]is proposed based on the ramp vortex g...The infrared imaging windows of the hyper/supersonic optical dome are encountering severe aero-optical effects[AOEs],so a flow control device,the ramp vortex generator array[RVGA]is proposed based on the ramp vortex generator to inhibit the supersonic mixing layers’AOE,which is done by the nanotracer-based planar laser scattering technique and ray-tracing method.The experiments prove that under different pressure conditions,RVGA can reduce the mean and standard deviation of the root mean square of the optical path difference[OPDrms]and reduce the supersonic mixing layers’thickness and mixture a great deal.The AOE of the pressure-matched mixing layer is the weakest.Higher RVGA results in better optical performance.RVGA has the potential to be applied to supersonic film cooling to reduce aero-optical aberrations.展开更多
Surface waves have a considerable effect on vertical mixing in the upper ocean.In the past two decades,the vertical mixing induced through nonbreaking surface waves has been used in ocean and climate models to improve...Surface waves have a considerable effect on vertical mixing in the upper ocean.In the past two decades,the vertical mixing induced through nonbreaking surface waves has been used in ocean and climate models to improve the simulation of the upper ocean.Thus far,several nonbreaking wave-induced mixing parameterization schemes have been proposed;however,no quantitative comparison has been performed among them.In this paper,a one-dimensional ocean model was used to compare the performances of five schemes,including those of Qiao et al.(Q),Hu and Wang(HW),Huang and Qiao(HQ),Pleskachevsky et al.(P),and Ghantous and Babanin(GB).Similar to previous studies,all of these schemes can decrease the simulated sea surface temperature(SST),increase the subsurface temperature,and deepen the mixed layer,thereby alleviating the common thermal deviation problem of the ocean model for upper ocean simulation.Among these schemes,the HQ scheme exhibited the weakest wave-induced mixing effect,and the HW scheme exhibited the strongest effect;the other three schemes exhibited roughly the same effect.In particular,the Q and P schemes exhibited nearly the same effect.In the simulation based on observations from the Ocean Weather Station Papa,the HQ scheme exhibited the best performance,followed by the Q scheme.In the experiment with the HQ scheme,the root-mean-square deviation of the simulated SST from the observations was 0.43℃,and the mixed layer depth(MLD)was 2.0 m.As a contrast,the deviations of the SST and MLD reached 1.25℃ and 8.4 m,respectively,in the experiment without wave-induced mixing.展开更多
A recent study by Liu et al.(2020)suggested that due to the saturation of equatorially trapped planetary waves with different dynamical types,temporal periods,meridional and baroclinic modes,complex layer structures o...A recent study by Liu et al.(2020)suggested that due to the saturation of equatorially trapped planetary waves with different dynamical types,temporal periods,meridional and baroclinic modes,complex layer structures of vertical velocity shear and hence turbulent mixing could frequently occur in the thermocline of the eastern equatorial Pacific.We investigated the occurrence of the interior turbulent mixing as indicated by shear instabilities,above the Equatorial Undercurrent(EUC)core at three equatorial sites along 140°W,170°W,and 165°E,respectively,based mainly on data from the Tropical Atmosphere and Ocean(TAO)mooring array.We found that turbulent mixing bursts persisted in the thermocline of all three sites.Specifically,the interior turbulent mixing layers(ITMLs)could occur in probability of approximately 68%,53%,and 48%at the three sites,respectively.The overall occurrence probability shows obvious and similar biannual variations at 140°W and 170°W,which is higher in boreal from late summer to winter and lower in spring.Vertically,the ITMLs are primarily located above the EUC core and prevail in deeper(shallower)layers from late summer to winter(spring).Most ITMLs(70%)lasted for hours to 3 days,and a few of them(15%)for more than 7 days.The thicknesses of ITMLs were concentrated between 15 and 55 m.At 165°E,the vertical distribution of ITML occurrence probability was different from that at 140°W and 170°W,as it did not show a preference for depths;the durations of ITMLs are short(also from hours to several days)and their thicknesses were between 5 and 25 m.These properties,particularly the high occurrence probability,and short durations demonstrated the persistence of thermocline mixing in the western to eastern equatorial Pacific thermocline and confirmed the generation mechanism by persistent equatorial waves as well.展开更多
Three high order shock-capturing schemes are compared for large eddy simulations(LES)of temporally evolving mixing layers for different convective Mach numbers ranging from the quasi-incompressible regime to highly co...Three high order shock-capturing schemes are compared for large eddy simulations(LES)of temporally evolving mixing layers for different convective Mach numbers ranging from the quasi-incompressible regime to highly compressible supersonic regime.The considered high order schemes are fifth-order WENO(WENO5),seventh-order WENO(WENO7)and the associated eighth-order central spatial base scheme with the dissipative portion of WENO7 as a nonlinear post-processing filter step(WENO7fi).This high order nonlinear filter method of Yee&Sjogreen is designed for accurate and efficient simulations of shock-free compressible turbulence,turbulence with shocklets and turbulence with strong shocks with minimum tuning of scheme parameters.The LES results by WENO7fi using the same scheme parameter agree well with experimental results compiled by Barone et al.,and published direct numerical simulations(DNS)work of Rogers&Moser and Pantano&Sarkar,whereas results by WENO5 andWENO7 compare poorly with experimental data and DNS computations.展开更多
El Ni?o–Southern Oscillation(ENSO) exhibits a distinctive phase-locking characteristic, first expressed during its onset in boreal spring, developing during summer and autumn, reaching its peak towards winter, and de...El Ni?o–Southern Oscillation(ENSO) exhibits a distinctive phase-locking characteristic, first expressed during its onset in boreal spring, developing during summer and autumn, reaching its peak towards winter, and decaying over the next spring. Several studies have demonstrated that this feature arises as a result of seasonal variation in the growth rate of ENSO as expressed by the sea surface temperature(SST). The bias towards simulating the phase locking of ENSO by many state-of-the-art climate models is also attributed to the unrealistic depiction of the growth rate. In this study, the seasonal variation of SST growth rate in the Ni?o-3.4 region(5°S–5°N, 120°–170°W) is estimated in detail based on the mixed layer heat budget equation and recharge oscillator model during 1981–2020. It is suggested that the consideration of a variable mixed layer depth is essential to its diagnostic process. The estimated growth rate has a remarkable seasonal cycle with minimum rates occurring in spring and maximum rates evident in autumn. More specifically, the growth rate derived from the meridional advection(surface heat flux) is positive(negative) throughout the year. Vertical diffusion generally makes a negative contribution to the evolution of growth rate and the magnitude of vertical entrainment represents the smallest contributor. Analysis indicates that the zonal advective feedback is regulated by the meridional immigration of the intertropical convergence zone, which approaches its southernmost extent in February and progresses to its northernmost location in September, and dominates the seasonal variation of the SST growth rate.展开更多
As an economically critical pelagic migratory species,yellowfin tuna(Thunnus albacores,YFT)is very sensible to physical and environmental conditions,such as sea surface temperature(SST),ocean heat content(OHC),and the...As an economically critical pelagic migratory species,yellowfin tuna(Thunnus albacores,YFT)is very sensible to physical and environmental conditions,such as sea surface temperature(SST),ocean heat content(OHC),and the mixed layer depth(MLD).We investigated the impact of SST,OHC,and MLD on fluctuations of YFT catch in the western/eastern Indian Ocean using the long time series of 63-year environmental and YFT datasets.We found that the impact of SST on YFT was heavily overestimated in the past,and MLD plays a more critical role in the YFT catch fluctuation.When the MLD deepens(>34.8 m),SST was more influential in predicting the catches of YFT than OHC in the western Indian Ocean,and OHC was more critical to YFT than SST in the eastern Indian Ocean.However,when the MLD shallows(<34.8 m),MLD was more vital to predict the catch per unit effort(CPUE)of YFT than SST/OHC in the western.After 2000,there was an asynchronous pattern of YFT CPUE induced by higher frequency variations and ocean hiatus of SST/OHC signals in the western and eastern Indian Oceans basins.The impact of the subsurface hiatus may induce the decrease of YFT in the eastern Indian Ocean.The above findings clarified a non-stationary relationship between the environmental factors and catches of YFT and provided new insights into variations in YFT abundance.展开更多
Antarctic coastal polynyas are biological hotspots in the Southern Ocean that support the abundance of hightrophic-level predators and are important for carbon cycling in the high-latitude oceans.In this study,we exam...Antarctic coastal polynyas are biological hotspots in the Southern Ocean that support the abundance of hightrophic-level predators and are important for carbon cycling in the high-latitude oceans.In this study,we examined the interannual variation of summertime phytoplankton biomass in the Marguerite Bay polynya(MBP)in the western Antarctic Peninsula area,and linked such variability to the Southern Annular Mode(SAM)that dominated the southern hemisphere extratropical climate variability.Combining satellite data,atmosphere reanalysis products and numerical simulations,we found that the interannual variation of summer chlorophyll-a(Chl-a)concentration in the MBP is significantly and negatively correlated with the spring SAM index,and weakly correlated with the summer SAM index.The negative relation between summer Chl-a and spring SAM is due to weaker spring vertical mixing under a more positive SAM condition,which would inhibit the supply of iron from deep layers into the surface euphotic layer.The negative relation between spring mixing and spring SAM results from greater precipitation rate over the MBP region in positive SAM phase,which leads to lower salinity in the ocean surface layer.The coupled physical-biological mechanisms between SAM and phytoplankton biomass revealed in this study is important for us to predict the future variations of phytoplankton biomasses in Antarctic polynyas under climate change.展开更多
An obvious trend shift in the annual mean and winter mixed layer depth(MLD)in the Antarctic Circumpolar Current(ACC)region was detected during the 1960–2021 period.Shallowing trends stopped in mid-1980s,followed by a...An obvious trend shift in the annual mean and winter mixed layer depth(MLD)in the Antarctic Circumpolar Current(ACC)region was detected during the 1960–2021 period.Shallowing trends stopped in mid-1980s,followed by a period of weak trends.The MLD deepening trend difference between the two periods were mainly distributed in the western areas in the Drake Passage,the areas north to Victoria Land and Wilkes Land,and the central parts of the South Indian sector.The newly formed ocean current shear due to the meridional shift of the ACC flow axis between the two periods is the dominant driver for the MLD trends shift distributed in the western areas in the Drake Passage and the central parts of the South Indian sector.The saltier trends in the regions north to Victoria Land and Wilkes Land could be responsible for the strengthening mixing processes in this region.展开更多
Wavelet analysis is applied to the results obtained by the direct numerical simulation of a three-dimensional (3D) mixing layer in order to investigate coherent structures in dimension of scale. First, 3D orthonormal ...Wavelet analysis is applied to the results obtained by the direct numerical simulation of a three-dimensional (3D) mixing layer in order to investigate coherent structures in dimension of scale. First, 3D orthonormal wavelet bases are constructed, and the corresponding decomposition algorithm is developed. Then the Navier-Stokes equations are transformed into the wavelet space and the architecture for multi-scale analysis is established. From this architecture, the coarse field images in different scales are obtained and some local statistical quantities are calculated. The results show that, with the development of a mixing layer, the energy spectrum densities for different wavenumbers increase and the energy is transferred from the average flow to vortex structures in different scales. Due to the non-linear interactions between different scales, cascade processes of energy are very complex. Because vortices always roll and pair at special areas, for a definite scale, the energy is obtained from other scales at some areas while it is transferred to other scales at other areas. In addition, energy dissipation and transfer always occur where an intense interaction between vortices exists.展开更多
The seasonal variation of mixing layer depth (MLD) in the ocean is determined by a wind stress and a buoy- ance flux. A South China Sea (SCS) ocean data assimilation system is used to analyze the seasonal cycle of...The seasonal variation of mixing layer depth (MLD) in the ocean is determined by a wind stress and a buoy- ance flux. A South China Sea (SCS) ocean data assimilation system is used to analyze the seasonal cycle of its MLD. It is found that the variability of MLD in the SCS is shallow in summer and deep in winter, as is the case in general. Owing to local atmosphere forcing and ocean dynamics, the seasonal variability shows a regional characteristic in the SCS. In the northern SCS, the MLD is shallow in summer and deep in winter, affected coherently by the wind stress and the buoyance flux. The variation of MLD in the west is close to that in the central SCS, influenced by the advection of strong western boundary currents. The eastern SCS presents an annual cycle, which is deep in summer and shallow in winter, primarily impacted by a heat flux on the air-sea interface. So regional characteristic needs to be cared in the analysis about the MLD of SCS.展开更多
The equations for fiber suspensions in an evolving mixing layer were solved by the spectral method, and the trajectory and orientation of fibers were calculated based on the slender body theory. The calculated spatial...The equations for fiber suspensions in an evolving mixing layer were solved by the spectral method, and the trajectory and orientation of fibers were calculated based on the slender body theory. The calculated spatial and orientation distributions of fibers are consistent with the experimental ones that were performed in this paper. The relationship between the microstructure of fibers and additional stress was examined. The results show that the spatial and orientation distributions of fibers are heterogeneous because of the influence of coherent vortices in the flow, which leads to the heterogeneity of the additional stress. The degree of heterogeneity increases with the increasing of St number and fiber aspect ratio. The fibers in the flow make the momentum loss thickness of the mixing layer thicker and accelerate the vorticity dispersion.展开更多
Both the level 2.5 Mellor-Yamada turbulence closure scheme (MY) and K-profile parameterization (KPP) are popularly used by the ocean modeling community. The MY and the KPP are improved through including the non-br...Both the level 2.5 Mellor-Yamada turbulence closure scheme (MY) and K-profile parameterization (KPP) are popularly used by the ocean modeling community. The MY and the KPP are improved through including the non-breaking surface wave-induced vertical mixing (Bv), and the improved schemes were tested by using continuous data at the Papa ocean weather station (OWS) during 1961-1965. The numerical results showed that the Bv can make the temperature simulations fit much better with the continuous data from Papa Sta- tion. The two improved schemes overcame the shortcomings of predicting too shallow upper mixed layer depth and consequently overheated sea surface temperature during summertime, which are in fact com- mon problems for all turbulence closure models. Statistical analysis showed that the Bv effectively reduced the mean absolute error and root mean square error of the upper layer temperature and increased the corre- lation coefficient between simulation and the observation. Furthermore, the performance of vertical mixing induced by shear instability and the Bv is also compared. Both the temperature structure and its seasonal cycle significantly improved by including the Bv, regardless of whether shear instability was included or not, especially for the KPP mixing scheme, which suggested that Bv played a dominant role in the upper ocean where the mean current was relatively weak, such as at Papa Station. These results may provide a clue to improve ocean circulation models.展开更多
Climatology of the isothermal layer depth (ILD) and the mixed layer depth (MLD) has been produced from in-situ temperaturesalinity observations in the East China Sea (ECS) since 1925. The methods applied on the ...Climatology of the isothermal layer depth (ILD) and the mixed layer depth (MLD) has been produced from in-situ temperaturesalinity observations in the East China Sea (ECS) since 1925. The methods applied on the global are used to compute the ILD and the MLD in the ECS with a temperature criterion AT=0. 8 ℃ for the ILD, and a density criterion with a threshold △σθ corresponding to fixed △T=0. 8 ℃ for the MLD, respectively. With the derived climatology ILD and MLD, the monthly variations of the barrier layer (BL) and the compensation layer (CL) in the ECS are analyzed. The BL mainly exists in the shallow water region of the ECS during April-June with thickness larger than 15 m. From December to next March, the area along the shelf break from northeast of Taiwan Island to the northeast ECS is characterized by the CL. Two kinds of main temperature - salinity structures of the CL in this area are given.展开更多
The nano-particle-based planar laser scattering (NPLS) technique is used to measure the density distribution in the supersonic mixing layer of the convective Mach number 0.12, and the optical path difference (OPL)...The nano-particle-based planar laser scattering (NPLS) technique is used to measure the density distribution in the supersonic mixing layer of the convective Mach number 0.12, and the optical path difference (OPL), which is quite crucial for the study of aero-optics, is obtained by post processing. Based on the high spatiotemporal resolutions of the NPLS, the structure of the OPL is ana]ysed using wavelet methods. The coherent structures of the OPL are extracted using three methods, including the methods of thresholding the coefficients of the orthogonal wavelet transform and the wavelet packet transform, and preserving a number of wavelet packet coefficients with the largest amplitudes determined by the entropy dimension. Their performances are compared, and the method using the wavelet packet is the best. Based on the viewpoint of multifractals, we study the OPL by the wavelet transform maxima method (WTMM), and the result indicates that its scaling behaviour is evident.展开更多
基金The project supported by the Zhejiang Province Natural Science Special Fund for Youth Scientists' Cultivation.
文摘The evolution of the three-dimensional time-developing mixing layer was simulated numerically using the pseudo-spectral method. The initial perturbations consisted of the two-dimensional fundamental wave and the streamwise-invariant three-dimensional disturbance. A comparison of the formations of the streamwise vortices with different amplitude functions for three-dimensional disturbances was made. In one case the results are similar to that of Rogers and Moser (1992), whereas a different way in which the quadrupole forms and sudden expansion of the rib were observed in another case. The simulation also confirms that stretching by the forming roller rather than Rayleigh centrifugal instability is responsible for the formation of the rib. Finally, numerical flow visualization results were presented. (Edited author abstract) 9 Refs.
基金supported by the National Natural Science Foundation of China (90716008,10572004 and 10921202)MOST 973 Project (2009CB724100) and CSSA
文摘Discontinuous Galerkin(DG) method is known to have several advantages for flow simulations,in particular,in fiexible accuracy management and adaptability to mesh refinement. In the present work,the DG method is developed for numerical simulations of both temporally and spatially developing mixing layers. For the temporally developing mixing layer,both the instantaneous fiow field and time evolution of momentum thickness agree very well with the previous results. Shocklets are observed at higher convective Mach numbers and the vortex paring manner is changed for high compressibility. For the spatially developing mixing layer,large-scale coherent structures and self-similar behavior for mean profiles are investigated. The instantaneous fiow field for a three-dimensional compressible mixing layer is also reported,which shows the development of largescale coherent structures in the streamwise direction. All numerical results suggest that the DG method is effective in performing accurate numerical simulations for compressible shear fiows.
基金The project supported by the National Natural Science Foundation of China (19642001) and Deutsche Forschungsgemeinschaft (DFG)
文摘In this paper, several mathematical models for the large scale structures in some special kinds of mixing layers, which might be practically useful for enhancing the mixing, are proposed. First, the linear growth rate of the large scale structures in the mixing layers was calculated. Then, using the much improved weakly non-linear theory, combined with the energy method, the non-linear evolution of large scale structures in two special mixing layer configurations is calculated. One of the mixing lavers has equal magnitudes of the upstream velocity vectors, while the angles between the velocity vectors and the trailing edge were pi /2 - phi and pi /2 + phi, respectively. The other mixing layer was generated by a splitter-plate with a 45-degree-sweep trailing edge.
文摘Numerical simulations have been performed in time-developing plane mixing layers of the viscoelastic second-order fluids with pseudo-spectral method. Roll-up, pairing and merging of large eddies were examined at high Reynolds numbers and low Deborah numbers. The effect of viscoelastics on the evolution of the large coherent structure was shown by making a comparison between the second-order and Newtonian fluids at the same Reynolds numbers.
基金supported by National Nature Science Foundation of China(90716008,10572004,and 11172006)by MOST 973 Project(2009CB724100)
文摘Semi-periodic structures namely inclined wavy structures (IWS) are experimentally observed in compressible mixing layers at two convective Mach numbers (Mc = 0.11 and 0.47). Flow structures are visualized by the laserinduced planar laser Mie scattering (PLMS) technique. Two methods are developed to investigate the spatial distribu- tion and geometry of IWS: (1) the dominant mode extrac- tion (DME) method, to extract the dominant modes of IWS from the streamwise gray-level fluctuation, and (2) the phase tracking (PT) method, to identify the shape of IWS. The re- sults suggest that pressure perturbations account for the for- marion of IWS in the initial mixing region and the joint effect of dilatation and coherent vortices enhances IWS in the well- developed region. The large transverse (cross-flow) scale of the IWS and their relation to coherent vortices (CV) indicate that the disturbance originated from CV in the mixing center propagates far into the free streams. The DME and the PT method are shown to be the effective tools to study the geometrical features of wavy structures in compressible shear flows.
基金supported by the National Natural Sci- ence Foundation of China(11102198)
文摘The method of nonlinear parabolized stability equations(PSE) is applied in the simulation of vortex structures in compressible mixing layer.The spatially-evolving unstable waves,which dominate the vortex structure,are investigated through spatial marching method.The instantaneous flow field is obtained by adding the harmonic waves to basic flow.The results show that T-S waves do not keep growing exponentially as the linear evolution,the energy transfer to high order harmonic modes,and that finally all harmonic modes get saturated due to nonlinear interaction.The mean flow distortion induced by the nonlinear interaction between the harmonic modes and their conjugate harmonic ones,makes great change of the average flow and increases the thickness of mixing layer. PSE methods can well capture the two- and three-dimensional large scale nonlinear vortex structures in mixing layers such as vortex roll-up,vortex pairing,and A vortex.
基金supported by the National Natural Science Foundation of China(No.12102463)the National Defense Basic Scientific Research Program of China(No.2022-JCJQJJ-1123)the Natural Science Foundation of Hunan Province(No.2021JJ40652)。
文摘The infrared imaging windows of the hyper/supersonic optical dome are encountering severe aero-optical effects[AOEs],so a flow control device,the ramp vortex generator array[RVGA]is proposed based on the ramp vortex generator to inhibit the supersonic mixing layers’AOE,which is done by the nanotracer-based planar laser scattering technique and ray-tracing method.The experiments prove that under different pressure conditions,RVGA can reduce the mean and standard deviation of the root mean square of the optical path difference[OPDrms]and reduce the supersonic mixing layers’thickness and mixture a great deal.The AOE of the pressure-matched mixing layer is the weakest.Higher RVGA results in better optical performance.RVGA has the potential to be applied to supersonic film cooling to reduce aero-optical aberrations.
基金supported by the Laoshan Laboratory(No.LSKJ202201600)the National Key Research and Development Program of China(No.2022YFC2808304).
文摘Surface waves have a considerable effect on vertical mixing in the upper ocean.In the past two decades,the vertical mixing induced through nonbreaking surface waves has been used in ocean and climate models to improve the simulation of the upper ocean.Thus far,several nonbreaking wave-induced mixing parameterization schemes have been proposed;however,no quantitative comparison has been performed among them.In this paper,a one-dimensional ocean model was used to compare the performances of five schemes,including those of Qiao et al.(Q),Hu and Wang(HW),Huang and Qiao(HQ),Pleskachevsky et al.(P),and Ghantous and Babanin(GB).Similar to previous studies,all of these schemes can decrease the simulated sea surface temperature(SST),increase the subsurface temperature,and deepen the mixed layer,thereby alleviating the common thermal deviation problem of the ocean model for upper ocean simulation.Among these schemes,the HQ scheme exhibited the weakest wave-induced mixing effect,and the HW scheme exhibited the strongest effect;the other three schemes exhibited roughly the same effect.In particular,the Q and P schemes exhibited nearly the same effect.In the simulation based on observations from the Ocean Weather Station Papa,the HQ scheme exhibited the best performance,followed by the Q scheme.In the experiment with the HQ scheme,the root-mean-square deviation of the simulated SST from the observations was 0.43℃,and the mixed layer depth(MLD)was 2.0 m.As a contrast,the deviations of the SST and MLD reached 1.25℃ and 8.4 m,respectively,in the experiment without wave-induced mixing.
基金Supported by the National Natural Science Foundation of China(NSFC)(No.41730534)the Laoshan Laboratory Science and Technology Innovation Program(No.LSKJ 202202502)+1 种基金the NSFC(Nos.41976012,42090044)the Strategic Priority Research Program of Chinese Academy of Sciences(No.XDB42000000)。
文摘A recent study by Liu et al.(2020)suggested that due to the saturation of equatorially trapped planetary waves with different dynamical types,temporal periods,meridional and baroclinic modes,complex layer structures of vertical velocity shear and hence turbulent mixing could frequently occur in the thermocline of the eastern equatorial Pacific.We investigated the occurrence of the interior turbulent mixing as indicated by shear instabilities,above the Equatorial Undercurrent(EUC)core at three equatorial sites along 140°W,170°W,and 165°E,respectively,based mainly on data from the Tropical Atmosphere and Ocean(TAO)mooring array.We found that turbulent mixing bursts persisted in the thermocline of all three sites.Specifically,the interior turbulent mixing layers(ITMLs)could occur in probability of approximately 68%,53%,and 48%at the three sites,respectively.The overall occurrence probability shows obvious and similar biannual variations at 140°W and 170°W,which is higher in boreal from late summer to winter and lower in spring.Vertically,the ITMLs are primarily located above the EUC core and prevail in deeper(shallower)layers from late summer to winter(spring).Most ITMLs(70%)lasted for hours to 3 days,and a few of them(15%)for more than 7 days.The thicknesses of ITMLs were concentrated between 15 and 55 m.At 165°E,the vertical distribution of ITML occurrence probability was different from that at 140°W and 170°W,as it did not show a preference for depths;the durations of ITMLs are short(also from hours to several days)and their thicknesses were between 5 and 25 m.These properties,particularly the high occurrence probability,and short durations demonstrated the persistence of thermocline mixing in the western to eastern equatorial Pacific thermocline and confirmed the generation mechanism by persistent equatorial waves as well.
基金Department of Energy by Lawrence Livermore National Laboratory under Contract DE-AC52-07NA27344.
文摘Three high order shock-capturing schemes are compared for large eddy simulations(LES)of temporally evolving mixing layers for different convective Mach numbers ranging from the quasi-incompressible regime to highly compressible supersonic regime.The considered high order schemes are fifth-order WENO(WENO5),seventh-order WENO(WENO7)and the associated eighth-order central spatial base scheme with the dissipative portion of WENO7 as a nonlinear post-processing filter step(WENO7fi).This high order nonlinear filter method of Yee&Sjogreen is designed for accurate and efficient simulations of shock-free compressible turbulence,turbulence with shocklets and turbulence with strong shocks with minimum tuning of scheme parameters.The LES results by WENO7fi using the same scheme parameter agree well with experimental results compiled by Barone et al.,and published direct numerical simulations(DNS)work of Rogers&Moser and Pantano&Sarkar,whereas results by WENO5 andWENO7 compare poorly with experimental data and DNS computations.
基金supported by the National Natural Science Foundation of China (Grant No. 42192564)Guangdong Major Project of Basic and Applied Basic Research (Grant No. 2020B0301030004)the Ministry of Science and Technology of the People's Republic of China (Grant No.2020YFA0608802)。
文摘El Ni?o–Southern Oscillation(ENSO) exhibits a distinctive phase-locking characteristic, first expressed during its onset in boreal spring, developing during summer and autumn, reaching its peak towards winter, and decaying over the next spring. Several studies have demonstrated that this feature arises as a result of seasonal variation in the growth rate of ENSO as expressed by the sea surface temperature(SST). The bias towards simulating the phase locking of ENSO by many state-of-the-art climate models is also attributed to the unrealistic depiction of the growth rate. In this study, the seasonal variation of SST growth rate in the Ni?o-3.4 region(5°S–5°N, 120°–170°W) is estimated in detail based on the mixed layer heat budget equation and recharge oscillator model during 1981–2020. It is suggested that the consideration of a variable mixed layer depth is essential to its diagnostic process. The estimated growth rate has a remarkable seasonal cycle with minimum rates occurring in spring and maximum rates evident in autumn. More specifically, the growth rate derived from the meridional advection(surface heat flux) is positive(negative) throughout the year. Vertical diffusion generally makes a negative contribution to the evolution of growth rate and the magnitude of vertical entrainment represents the smallest contributor. Analysis indicates that the zonal advective feedback is regulated by the meridional immigration of the intertropical convergence zone, which approaches its southernmost extent in February and progresses to its northernmost location in September, and dominates the seasonal variation of the SST growth rate.
基金Supported by the National Natural Science Foundation of China(Nos.42090044,42376175,U2006211)the Marine S&T Fund of Laoshan Laboratory(Qingdao)(No.LSKJ202204302)。
文摘As an economically critical pelagic migratory species,yellowfin tuna(Thunnus albacores,YFT)is very sensible to physical and environmental conditions,such as sea surface temperature(SST),ocean heat content(OHC),and the mixed layer depth(MLD).We investigated the impact of SST,OHC,and MLD on fluctuations of YFT catch in the western/eastern Indian Ocean using the long time series of 63-year environmental and YFT datasets.We found that the impact of SST on YFT was heavily overestimated in the past,and MLD plays a more critical role in the YFT catch fluctuation.When the MLD deepens(>34.8 m),SST was more influential in predicting the catches of YFT than OHC in the western Indian Ocean,and OHC was more critical to YFT than SST in the eastern Indian Ocean.However,when the MLD shallows(<34.8 m),MLD was more vital to predict the catch per unit effort(CPUE)of YFT than SST/OHC in the western.After 2000,there was an asynchronous pattern of YFT CPUE induced by higher frequency variations and ocean hiatus of SST/OHC signals in the western and eastern Indian Oceans basins.The impact of the subsurface hiatus may induce the decrease of YFT in the eastern Indian Ocean.The above findings clarified a non-stationary relationship between the environmental factors and catches of YFT and provided new insights into variations in YFT abundance.
基金The Key Research&Development Program of the Ministry of Science and Technology of China under contract No.2022YFC2807601the National Natural Science Foundation of China under contract Nos 41941008 and 41876221+3 种基金the Fund of Shanghai Science and Technology Committee under contract Nos 20230711100 and 21QA1404300the Impact and Response of Antarctic Seas to Climate Change funded by the Chinese Arctic and Antarctic Administration under contract No.IRASCC 1-02-01Bthe National Key Research and Development Program of China under contract No.2019YFC1509102the Shanghai Pilot Program for Basic Research—Shanghai Jiao Tong University under contract No.21TQ1400201。
文摘Antarctic coastal polynyas are biological hotspots in the Southern Ocean that support the abundance of hightrophic-level predators and are important for carbon cycling in the high-latitude oceans.In this study,we examined the interannual variation of summertime phytoplankton biomass in the Marguerite Bay polynya(MBP)in the western Antarctic Peninsula area,and linked such variability to the Southern Annular Mode(SAM)that dominated the southern hemisphere extratropical climate variability.Combining satellite data,atmosphere reanalysis products and numerical simulations,we found that the interannual variation of summer chlorophyll-a(Chl-a)concentration in the MBP is significantly and negatively correlated with the spring SAM index,and weakly correlated with the summer SAM index.The negative relation between summer Chl-a and spring SAM is due to weaker spring vertical mixing under a more positive SAM condition,which would inhibit the supply of iron from deep layers into the surface euphotic layer.The negative relation between spring mixing and spring SAM results from greater precipitation rate over the MBP region in positive SAM phase,which leads to lower salinity in the ocean surface layer.The coupled physical-biological mechanisms between SAM and phytoplankton biomass revealed in this study is important for us to predict the future variations of phytoplankton biomasses in Antarctic polynyas under climate change.
基金The National Natural Science Foundation of China under contract No.41605052。
文摘An obvious trend shift in the annual mean and winter mixed layer depth(MLD)in the Antarctic Circumpolar Current(ACC)region was detected during the 1960–2021 period.Shallowing trends stopped in mid-1980s,followed by a period of weak trends.The MLD deepening trend difference between the two periods were mainly distributed in the western areas in the Drake Passage,the areas north to Victoria Land and Wilkes Land,and the central parts of the South Indian sector.The newly formed ocean current shear due to the meridional shift of the ACC flow axis between the two periods is the dominant driver for the MLD trends shift distributed in the western areas in the Drake Passage and the central parts of the South Indian sector.The saltier trends in the regions north to Victoria Land and Wilkes Land could be responsible for the strengthening mixing processes in this region.
基金The project supported by the Research Fund for the Doctoral Program of Higher Educationthe National Natural Science Foundation for Outstanding Youth of China (19925210)
文摘Wavelet analysis is applied to the results obtained by the direct numerical simulation of a three-dimensional (3D) mixing layer in order to investigate coherent structures in dimension of scale. First, 3D orthonormal wavelet bases are constructed, and the corresponding decomposition algorithm is developed. Then the Navier-Stokes equations are transformed into the wavelet space and the architecture for multi-scale analysis is established. From this architecture, the coarse field images in different scales are obtained and some local statistical quantities are calculated. The results show that, with the development of a mixing layer, the energy spectrum densities for different wavenumbers increase and the energy is transferred from the average flow to vortex structures in different scales. Due to the non-linear interactions between different scales, cascade processes of energy are very complex. Because vortices always roll and pair at special areas, for a definite scale, the energy is obtained from other scales at some areas while it is transferred to other scales at other areas. In addition, energy dissipation and transfer always occur where an intense interaction between vortices exists.
基金The National Basic Research Program of China under contract Nos 2011CB403505 and 2011CB403504the National NaturalScience Foundation of China under contract No.41206007+2 种基金the City University of Hong Kong Stritegic Research Grants under contract Nos 7002917 and 7002780the Knowledge Innovation Project for Distinguished Young Scholar of The Chinese Academy of Sciences under contract KZCX2-EWQN203the foundation for operational development of the National Marine Environment Forecasting Center under contract No.2013006
文摘The seasonal variation of mixing layer depth (MLD) in the ocean is determined by a wind stress and a buoy- ance flux. A South China Sea (SCS) ocean data assimilation system is used to analyze the seasonal cycle of its MLD. It is found that the variability of MLD in the SCS is shallow in summer and deep in winter, as is the case in general. Owing to local atmosphere forcing and ocean dynamics, the seasonal variability shows a regional characteristic in the SCS. In the northern SCS, the MLD is shallow in summer and deep in winter, affected coherently by the wind stress and the buoyance flux. The variation of MLD in the west is close to that in the central SCS, influenced by the advection of strong western boundary currents. The eastern SCS presents an annual cycle, which is deep in summer and shallow in winter, primarily impacted by a heat flux on the air-sea interface. So regional characteristic needs to be cared in the analysis about the MLD of SCS.
基金The project supported by the Doctoral Program of Higher Education in China(20030335001)
文摘The equations for fiber suspensions in an evolving mixing layer were solved by the spectral method, and the trajectory and orientation of fibers were calculated based on the slender body theory. The calculated spatial and orientation distributions of fibers are consistent with the experimental ones that were performed in this paper. The relationship between the microstructure of fibers and additional stress was examined. The results show that the spatial and orientation distributions of fibers are heterogeneous because of the influence of coherent vortices in the flow, which leads to the heterogeneity of the additional stress. The degree of heterogeneity increases with the increasing of St number and fiber aspect ratio. The fibers in the flow make the momentum loss thickness of the mixing layer thicker and accelerate the vorticity dispersion.
基金The National Basic Research Program("973"Program)of China under contract No.2010CB950303the Public Science and Technology Research Funds Projects of Ocean under contract No.201105019the National Key Technology Research and Development Program of China under contract No.2011BAC03B02
文摘Both the level 2.5 Mellor-Yamada turbulence closure scheme (MY) and K-profile parameterization (KPP) are popularly used by the ocean modeling community. The MY and the KPP are improved through including the non-breaking surface wave-induced vertical mixing (Bv), and the improved schemes were tested by using continuous data at the Papa ocean weather station (OWS) during 1961-1965. The numerical results showed that the Bv can make the temperature simulations fit much better with the continuous data from Papa Sta- tion. The two improved schemes overcame the shortcomings of predicting too shallow upper mixed layer depth and consequently overheated sea surface temperature during summertime, which are in fact com- mon problems for all turbulence closure models. Statistical analysis showed that the Bv effectively reduced the mean absolute error and root mean square error of the upper layer temperature and increased the corre- lation coefficient between simulation and the observation. Furthermore, the performance of vertical mixing induced by shear instability and the Bv is also compared. Both the temperature structure and its seasonal cycle significantly improved by including the Bv, regardless of whether shear instability was included or not, especially for the KPP mixing scheme, which suggested that Bv played a dominant role in the upper ocean where the mean current was relatively weak, such as at Papa Station. These results may provide a clue to improve ocean circulation models.
基金The National Natural Science Foundation of China under contract Nos40776018 and 40730842the National Basic Research Program of China under contract No.2007CB816002
文摘Climatology of the isothermal layer depth (ILD) and the mixed layer depth (MLD) has been produced from in-situ temperaturesalinity observations in the East China Sea (ECS) since 1925. The methods applied on the global are used to compute the ILD and the MLD in the ECS with a temperature criterion AT=0. 8 ℃ for the ILD, and a density criterion with a threshold △σθ corresponding to fixed △T=0. 8 ℃ for the MLD, respectively. With the derived climatology ILD and MLD, the monthly variations of the barrier layer (BL) and the compensation layer (CL) in the ECS are analyzed. The BL mainly exists in the shallow water region of the ECS during April-June with thickness larger than 15 m. From December to next March, the area along the shelf break from northeast of Taiwan Island to the northeast ECS is characterized by the CL. Two kinds of main temperature - salinity structures of the CL in this area are given.
基金Projected supported by the Innovation Research Foundations for Postgraduates of National University of Defense Technology and Hunan Provincethe National Natural Science Foundation of China (Grant No. 61008037)
文摘The nano-particle-based planar laser scattering (NPLS) technique is used to measure the density distribution in the supersonic mixing layer of the convective Mach number 0.12, and the optical path difference (OPL), which is quite crucial for the study of aero-optics, is obtained by post processing. Based on the high spatiotemporal resolutions of the NPLS, the structure of the OPL is ana]ysed using wavelet methods. The coherent structures of the OPL are extracted using three methods, including the methods of thresholding the coefficients of the orthogonal wavelet transform and the wavelet packet transform, and preserving a number of wavelet packet coefficients with the largest amplitudes determined by the entropy dimension. Their performances are compared, and the method using the wavelet packet is the best. Based on the viewpoint of multifractals, we study the OPL by the wavelet transform maxima method (WTMM), and the result indicates that its scaling behaviour is evident.