Move a Tropical Cyclone with 4D-Var and Vortex Dynamical Initialization in WRF Model

Previous studies showed that 4 D-Var technique used for data assimilation could be modified for weather control. This study demonstrates the ability of 4 D-Var to influence the future path of a tropical cyclone by calculating perturbations in WRF simulation. Given the background error covariance matrix, the initial field is improved by the vortex dynamic initialization technique. Our results show that 4 D-Var can be applied to control the trajectory of simulated tropical cyclones by producing "optimal" perturbations. In the numerical simulation experiment of Typhoon Mitag in 2019, after this kind of weather control similar to data assimilation, the tropical cyclone moved obviously,and the damaging wind over the coastline weakened. The prediction results after the initial field modified by 4 D-Var have a great change, and the position of the tropical cyclone moved about 0.5° southeastward after assimilation,which misses the southeast coast of China. Moreover, the damaging wind is also weakened. Since the 4 D-Var is premised on the assumption that the model is perfect and does not consider the model error, then the research plan to consider model error and introduce new methods is discussed in the paper.

Liutex-based vortex dynamics:A preliminary study

Vortex dynamics,with the possibility of efficient flow control,is explored in this study based on the new introduced vortex definition and identification system of Liutex.With the six core elements of vortex identification,including(1)absolute strength,(2)relative strength,(3)local rotational axis,(4)global rotational axis,(5)vortex core size and(6)vortex boundary,provided by the Liutex system,it is possible to numerically devise strategies,primarily by introducing additional source terms in Navier-Stokes equations,which we call Liutex force field model here,to control the vortex regions.Two methodologies of centripetal force model and counter-rotation force model are preliminarily investigated in a cavitating flow around two-dimensional Clark-Y hydrofoil.It is found that Liutex based models are capable of illustrating the vortex dynamics and possibly strengthening or weakening the vortices.

Vortex dynamics and entropy generation in separated transitional flow over a compressor blade at various incidence angles

The transition process within a Laminar Separation Bubble(LSB)that formed on a compressor blade surface was investigated using Large Eddy Simulations(LESs)at a Reynolds number of 1.5×10;and incidence angles of 0°,+3°,and+5°.The vortex dynamics in the separated shear layers were compared at various incidence angles and its effects on the loss generation were clarified through entropy analysis.Results showed that transition onset,which was accurately identified by the Linear Stability Theory(LST),was significantly promoted at the increased incidence angle.As such,the development of LSB was suppressed and the relative role of viscous instability played in the transition process was weakened.At the incidence angle of 0°,two-dimensional spanwise vortices detached from the blade surface and roiled up periodically,which were further stretched and eventually evolved into large-scale hairpin vortices.As time passed,the fully developed hairpin vortices broke down into small-scale eddies.Meanwhile,the flow near the wall reversely ejected into the outer separated shear layers and a sweeping process happened subsequently,forcing the separated shear layers to reattach and accelerating the generation of turbulent fluctuations.By comparison,the strength of vortex rolling-up was weakened at higher incidence angles,and the vortex pairing and breakdown of large-scale vortices were less pronounced.Therefore,the level of turbulent fluctuations that generated in the separated shear layers was reduced.Detailed entropy analysis showed that the turbulent dissipation effect related to the Reynolds shear stresses determined the largest amount of positive entropy generation,which declined to a lower level as the incidence angle increased from 0°to+5°.Correspondingly,the profile loss was reduced by 50.4%.

Cascade excitation of vortex motion and reentrant superconductivity in flexible Nb thin films

High quality Nb films were successfully prepared on both flexible polyimide(PI)and rigid Al2O3substrates and their transport properties were systematically studied at various applied currents,external magnetic fields,and sample orientations.It is found that a curved Nb/PI film exhibits quite different superconducting transition and vortex dynamics compared to the flat Nb/Al2O3film.For the curved Nb/PI film,smooth superconducting transitions were obtained at low currents,while unexpected cascade structures were revealed in theρ(T)curves at high currents.We attribute this phenomenon to the gradient distribution of vortex density together with a variation of superconductivity along the curved film.In addition,reentrant superconductivity was induced in the curved Nb/PI thin film by properly choosing the measurement conditions.We attribute this effect to the vortex pinning from both in-plane vortices and out-of-plane vortices.This work reveals the complex transport properties of curved superconducting thin films,providing important insights for further theoretical investigations and practical developments of flexible superconductors.

Numerical investigation of flow over a two-dimensional square cylinder with a synthetic jet generated by a bi-frequency signal

The flow around a square cylinder with a synthetic jet positioned at the rear surface is numerically investigated with the unsteady Reynolds-averaged Navier-Stokes(URANS)method.Instead of the typical sinusoidal wave,a bi-frequency signal is adopted to generate the synthetic jet.The bi-frequency signal consists of a basic sinusoidal wave and a high-frequency wave.Cases with various amplitudes of the high-frequency component are simulated.It is found that synthetic jets actuated by bi-frequency signals can realize better drag reduction with lower energy consumption when appropriate parameter sets are applied.A new quantity,i.e.,the actuation efficiency Ae,is used to evaluate the controlling efficiency.The actuation efficiency Ae reaches its maximum of 0.2668 when the amplitude of the superposed high-frequency signal is 7.5%of the basic signal.The vortex structures and frequency characteristics are subsequently analyzed to investigate the mechanism of the optimization of the bi-frequency signal.When the synthetic jet is actuated by a single-frequency signal with a characteristic velocity of 0.112 m/s,the wake is asymmetrical.The alternative deflection of vortex pairs and the peak at half of the excitation frequency in the power spectral density(PSD)function are detected.In the bi-frequency cases with the same characteristic velocity,the wake gradually turns to be symmetrical with the increase in the amplitude of the high-frequency component.Meanwhile,the deflection of the vortex pairs and the peak at half of the excitation frequency gradually disappear as well.

Finite superconducting square wire-network based on two-dimensional crystalline Mo_(2)C

Superconducting wire-networks are paradigms to study Cooper pairing issues,vortex dynamics and arrangements.Recently,emergent low-dimensional crystalline superconductors were reported in the minimal-disorder limit,providing novel platforms to reveal vortices-related physics.Study on superconducting loops with high-crystallinity is thus currently demanded.Here,we report fabrication and transport measurement of finite square-network based on two-dimensional crystalline superconductor Mo_(2)C.We observe oscillations in the resistance as a function of the magnetic flux through the loops.Resistance dips at both matching field and fractional fillings are revealed.Temperature and current evolutions are carried out in magnetoresistance to study vortex dynamics.The amplitude of oscillation is enhanced due to the interaction between thermally activated vortices and the currents induced in the loops.The driving current reduces the effective activation energy for vortex,giving rise to stronger vortex interaction.Moreover,by the thermally activated vortex creep model,we derive the effective potential barrier for vortex dissipation,which shows well-defined correspondence with structures in magnetoresistance.Our work shows that low-dimensional crystalline superconducting network based on Mo_(2)C possesses pronounced potential in studying the modulation of vortex arrangements and dynamics,paving the way for further investigations on crystalline superconducting network with various configurations.

Slat cove dynamics of multi-element airfoil at low Reynolds number

The flow around the slat cove of a two-dimensional 30P30N multi-element airfoil is investigated with time-resolved particle image velocimetry(TR-PIV)at low Reynolds number(Rec=2.41×10^(4)and 4.61×10^(4)).The effects of angle of attack(α=8°,12°,and 16°)on the mean flow characteristics and vortex dynamics are discussed.The size of the recirculation within the slat cove and the intensity of the shed vortices originating from the slat cusp shear layer are found to generally decrease as the angle of attack increases.The joint time-frequency analyses show that disturbances of different frequencies exist in the slat cusp shear layer and they trigger the different vortex shedding patterns of the slat cusp shear layer.The self-sustained oscillation within the slat cove,normally observed at high Reynolds number(Re_(c)~10^(6)),is proved to be responsible for the disturbances of different frequencies and the related vortex dynamics in the current study.

Probing the fractional topological charge of a vortex light beam by using dynamic angular double slits

Vortex beams with fractional topological charge(FTC) have many special characteristics and novel applications.However, one of the obstacles for their application is the difficulty of precisely determining the FTC of fractional vortex beams. We find that when a vortex beam with an FTC illuminates a dynamic angular double slit(ADS), the far-field interference patterns that include the information of the FTC of the beam at the angular bisector direction of the ADS vary periodically. Based on this property, a simple dynamic ADS device and data fitting method can be used to precisely measure the FTC of a vortex light beam with an error of less than 5%.

Liutex-based vortex control with implications for cavitation suppression

Viewed as sinews and muscles of fluid motion,coherent vortical structures with their interactions are key to understanding the flow dynamics.Based upon this observation,we explore the possibility of efficient flow control by directly manipulating vortices numerically inside the flow field based on the vortex definition and identification system of Liutex.The objective is twofold:(1)to study the vortex dynamics,for example,by observing the response of the flow to strengthening or weakening of certain vortices,and(2)to obtain efficient vortex-based control strategies which might lead us to practical applications.In the present numerical study,the manipulating of vortices is achieved by introducing additional source(force)terms to the Navier-Stokes equations,which hereafter will be collectively called Liutex force field model.Methodologies including controlling the rotation strength and centripetal force of particular vortices are detailed in a flow past a cylinder with different control purposes at Reynolds number of 200.Further examples are provided with a cavitating flow around two-dimensional Clark-Y hydrofoil,with particular interests on cavitation suppression.It is illustrated particular vortex with cavitation encircled could be effectively suppressed.

Individual influence of pitching and plunging motions on flow structures over an airfoil during dynamic stall

The individual influence of pitching and plunging motions on flow structures is studied experimentally by changing the phase lag between the geometrical angle of attack and the plunging angle of attack.Five phase lags are chosen as the experimental parameters,while the Strouhal number,the reduced frequency and the Reynolds number are fixed.During the motion of the airfoil,the leading edge vortex,the reattached vortex and the secondary vortex are observed in the flow field.The leading edge vortex is found to be the main flow structure through the proper orthogonal decomposition.The increase of phase lag results in the increase of the leading edge velocity,which strongly influences the leading edge shear layer and the leading edge vortex.The plunging motion contributes to the development of the leading edge shear layer,while the pitching motion is the key reason for instability of the leading edge shear layer.It is also found that a certain increase of phase lag,around 34.15°in this research,can increase the airfoil lift.

On the role of hemodynamics in predicting rupture of the abdominal aortic aneurysm

Hemodynamics plays a crucial role in the growth of an abdominal aortic aneurysm(AAA)and its possible rupture.Due to the serious consequences that arise from the aneurysm rupture,the ability to predict its evolution and the need for surgery are of primary importance in the medical field.Furthermore,the presence of intraluminal thrombus(ILT)strongly affects the evolution of the pathology.In this study,we analyzed the influence of hemodynamics on the growth and possible rupture of AAAs.Numerical investigations of pulsatile non-Newtonian blood flow were performed in six patient-specific AAAs reconstructed from diagnostic images,having different sizes and shapes,and with or without ILT.Wall shear stress and vorticity distribution in the bulge and their evolution during the cardiac cycle were analyzed.The results indicate that blood flow dynamics acts synergistically with atherosclerotic degeneration in the development of the disease.The high surface complexity and tortuosity of the aneurysms significantly affect the blood motion,and the presence of inflection in the aneurysm centerline has a noticeable effect on the vortex dynamics.Links between regions of slow recirculating flows,low values of time-averaged wall shear stress,high values of oscillatory shear index,and zones of ILT deposition were found.In the absence of ILT,possible thrombus accumulation areas and consequent aneurysm growth were identified.The findings of this study highlight the importance of hemodynamics in assessing the vulnerability of the aortic wall and underline the crucial role of patient-specific investigations in predicting the rupture of individual aneurysms.

Criteria of tracking vortex surfaces in turbulent-like flows

We propose criteria of tracking vortex surfaces in complex flows based on the vortex-surface field (VSF). The criteria characterize the accuracy and Lagrangian tracking performance of the numerical VSF solution, and determine the time period when the vortex surface tracking is satisfactory. Moreover, we develop a turbulent-like flow combining large-scale coherent structures in the Taylor-Green flow and small-scale turbulent structures in homogeneous isotropic turbulence (HIT). From tracking of vortex surfaces during the effective tracking period, we find that the imposed HIT disturbance significantly wrinkles vortex surfaces. Subsequently, the wrinkled vortex tube with large vorticity magnitude tends to be further twisted, contributing to energy cascade, while the wrinkling is mitigated in the region with small vorticity magnitude.

A 3D Numerical Method for Studying Vortex Formation Behind a Moving Plate

In this paper,we introduce a three-dimensional numerical method for computing the wake behind a flat plate advancing perpendicular to the flow.Our numerical method is inspired by the panel method of J.Katz and A.Plotkin[J.Katz and A.Plotkin,Low-speed Aerodynamics,2001]and the 2D vortex blob method of Krasny[R.Krasny,Lectures in Appl.Math.,28(1991),pp.385–402].The accuracy of the method will be demonstrated by comparing the 3D computation at the center section of a very high aspect ratio plate with the corresponding two-dimensional computation.Furthermore,we compare the numerical results obtained by our 3D numerical method with the corresponding experimental results obtained recently by Ringuette[M.J.Ringuette,Ph.D.Thesis,2004]in the towing tank.Our numerical results are shown to be in excellent agreement with the experimental results up to the so-called formation time.

圆和非圆湍流合成射流结构演化和质量输运的拉格朗日方法分析

旋涡卷吸特性对合成射流的流动控制应用具有重要意义.为了研究湍流合成射流涡环的结构演化和瞬态卷吸特性,采用时间解析层析粒子图像测速技术,测量了相同参数下,圆形和长宽比为3的矩形合成射流的三维流场.利用拉格朗日拟序结构方法,刻画了涡环的三维拓扑变形和质量输运过程,并量化了连续涡环的卷吸特性.特别是对矩形涡环,由于涡环的对流速度沿周向分布不均匀,导致传统的移动参考系方法无法有效计算卷吸.针对该问题,采用基于有限时间李雅普诺夫指数“脊线”辨识旋涡边界的方法,估计涡泡体积,从而计算涡环的卷吸和动能转换效率.结果表明,矩形涡环能够卷吸更多的周围流质,产生比圆形涡环更大的涡旋强度;此外,矩形涡环从孔口喷出的流质中获得了更多动能,以维持涡环的轴系转换,这些特征有利于增强射流掺混.上述工作通过发展旋涡动力学的数学描述,为量化流动控制中的旋涡卷吸特性提供了重要参考.

Large eddy simulations of a turbulent mixing layer periodically excited with fundamental and third harmonic frequency

A better understanding of the mixing behavior of excited turbulent mixing layers is critical to a number of aerospace applications.Previous studies of excited turbulent mixing layers focused on single frequency excitation or the excitation with fundamental and its second harmonic frequency.There is a lack of detailed studies on applying low and higher frequency excitation.In this study,we have performed large-eddy simulations of periodically excited turbulent mixing layers.The excitation consists of a fundamental frequency and its third harmonic.We have used phase-averaging to identify the vortex structure and strength in the mixing layer,and we have studied the vortex dynamics.Two different vortex paring mechanisms are observed depending on the phase shift between the two excitation frequencies.The influence of these two mechanisms on the mixing of a passive scalar is also studied.It is found that exciting the mixing layer with these low and high frequencies has initially an adverse influence on the mixing process;however,it improves the mixing further downstream of the splitter plate with the excitation using a phase shift ofΔφ=πshowing the best mixing performance.The present works shed lights on the fundamental vortex dynamics,and has great potential for aeronautical,automotive and combustion engineering applications.

Liutex-based centripetal force field model for improving the resistance and wake performances of JBC ship sailing in calm water

For complex aerodynamic and hydrodynamic problems,the analysis of vortex is very important.The Liutex method is an eigenvalue-based method which is local,accurate,and unique,which can give an accurate definition of vortex,so the control of vortex can be implemented and effectively guaranteed.Based on Liutex method,two methodologies of centripetal force model and counter-rotation force model were proposed to illustrate the vortex dynamics and possibly strengthen or weaken the vortices.In this paper,the Liutex-based centripetal force model is applied by adding a source term to the Navier-Stokes equations.In order to investigate the influence of the constructed Liutex force model on the 3-dimensional flow around a slow-fat ship,the calm-water drag calculation result of JBC ship is regarded as the initial flow field,and the new resistance and wake performances of the ship are obtained after applying the centripetal force model to the flow field with different strengths.Several views of the comparisons of the new steady flow fields are shown,and the parametric study results indicate that the Liutex-based centripetal force model can effectively change the resistance and wake performances of the JBC ship,which provides a new idea and theoretical basis for the comprehensive hydrodynamic performance optimization of the ship hull.

Parametric study of Liutex-based force field models

The motion of fluid consists of different scales of coherent vortical structures.These vortical structures determine the characteristics of fluid motion and are key to understand fluid dynamics.In this paper,we study the fine control method of vortical structures based on the Liutex force field model.This is achieved by constructing a source term using Liutex and directly add it to the Navier-Stokes equations.To investigate the influence of the constructed Liutex force model on vortical structures,a flow past a cylinder at Reynolds number of 100 is numerically studied with different source term magnitude and region.The drag and lift forces on the cylinder,as well as the flow field near and behind the cylinder are compared and analyzed.Results show that Liutex force model can effectively strengthen or weaken the vortical structures based on different purpose.

雷诺数和施密特数对激波-气泡相互作用中变密度混合的影响

惯性约束聚变、水动力空化侵蚀和微观尺度激波加速流动中的激波-气泡相互作用(SBI)构型处于高黏度(低雷诺数)或低扩散率(高施密特数)环境中,其混合性能尚不清楚.通过测定宽范围内的粘度和扩散系数,采用高分辨率数值模拟研究了雷诺数Re和施密特数Sc对SBI流动结构和混合行为的影响.次级斜压涡度(SBV)增长明显受雷诺数的影响,表现为斜压涡和黏性耗散之间的生长抑制涡度动力学平衡.混合速率受施密特数的影响,标度为χ^(∗)=β·Sc^(−α),其中系数β∼Re−0.2和标度指数α∼Re−0.385.混合标度表明,较高雷诺数下的SBV拉伸增强导致混合速率对施密特数的依赖性减弱,这可能为实际应用中的混合增强设计开辟新的途径.

Topology Optimization of the Caudal Fin of the Three-Dimensional Self-Propelled Swimming Fish

Based on the boundary vorticity-flux theory,topology optimization of the caudal fin of the three-dimensional self-propelled swimming fish is investigated by combining unsteady computational fluid dynamics with moving boundary and topology optimization algorithms in this study.The objective functional of topology optimization is the function of swimming efficiency,swimming speed and motion direction control.The optimal caudal fin,whose topology is different from that of the natural fish caudal fin,make the 3D bionic fish achieve higher swimming efficiency,faster swimming speed and better maneuverability.The boundary vorticity-flux on the body surface of the 3D fish before and after optimization reveals the mechanism of high performance swimming of the topology optimization bionic fish.The comparative analysis between the swimming performance of the 3D topology optimization bionic fish and the 3D lunate tail bionic fish is also carried out,and the wake structures of two types of bionic fish show the physical nature that the swimming performance of the 3D topology optimization bionic fish is significantly better than the 3D lunate tail bionic fish.