Multi-scale physics-informed neural networks for solving high Reynolds number boundary layer flows based on matched asymptotic expansions

Multi-scale system remains a classical scientific problem in fluid dynamics,biology,etc.In the present study,a scheme of multi-scale Physics-informed neural networks is proposed to solve the boundary layer flow at high Reynolds numbers without any data.The flow is divided into several regions with different scales based on Prandtl's boundary theory.Different regions are solved with governing equations in different scales.The method of matched asymptotic expansions is used to make the flow field continuously.A flow on a semi infinite flat plate at a high Reynolds number is considered a multi-scale problem because the boundary layer scale is much smaller than the outer flow scale.The results are compared with the reference numerical solutions,which show that the msPINNs can solve the multi-scale problem of the boundary layer in high Reynolds number flows.This scheme can be developed for more multi-scale problems in the future.

Analytical solutions of turbulent boundary layer beneath forward-leaning waves

As a typical nonlinear wave,forward-leaning waves can be frequently encountered in the near-shore areas,which can impact coastal sediment transport significantly.Hence,it is of significance to describe the characteristics of the boundary layer beneath forward-leaning waves accurately,especially for the turbulent boundary layer.In this work,the linearized turbulent boundary layer model with a linear turbulent viscosity coefficient is applied,and the novel expression of the near-bed orbital velocity that has been worked out by the authors for forward-leaning waves of arbitrary forward-leaning degrees is further used to specify the free stream boundary condition of the bottom boundary layer.Then,a variable transformation is found so as to make the equation of the turbulent boundary layer model be solved analytically through a modified Bessel function.Consequently,an explicit analytical solution of the turbulent boundary layer beneath forward-leaning waves is derived by means of variable separation and variable transformation.The analytical solutions of the velocity profile and bottom shear stress of the turbulent boundary layer beneath forward-leaning waves are verified by comparing the present analytical results with typical experimental data available in the previous literature.

Responses of the field-aligned currents in the plasma sheet boundary layer to a geomagnetic storm

Geomagnetic storms can result in large magnetic field disturbances and intense currents in the magnetosphere and even on the ground.As an important medium of momentum and energy transport among the solar wind,magnetosphere,and ionosphere,field-aligned currents(FACs)can also be strengthened in storm times.This study shows the responses of FACs in the plasma sheet boundary layer(PSBL)observed by the Magnetospheric Multiscale(MMS)spacecraft in different phases of a large storm that lasted from May 27,2017,to May 29,2017.Most of the FACs were carried by electrons,and several FACs in the storm time also contained sufficient ion FACs.The FAC magnitudes were larger in the storm than in the quiet period,and those in the main phase were the strongest.In this case,the direction of the FACs in the main phase showed no preference for tailward or earthward,whereas the direction of the FACs in the recovery phase was mostly tailward.The results suggest that the FACs in the PSBL are closely related to the storm and could be driven by activities in the tail region,where the energy transported from the solar wind to the magnetosphere is stored and released as the storm is evolving.Thus,the FACs are an important medium of energy transport between the tail and the ionosphere,and the PSBL is a significant magnetosphere–ionosphere coupling region in the nightside.

Tide-Induced Mixing in the Bottom Boundary Layer in the Western East China Sea

The East China Sea(ECS)boasts a vast continental shelf,where strong tidal motions play an important role in the substance transport and energy budget.In this study,the tide-induced mixing in the bottom boundary layer in the western ECS is analyzed based on records measured by moored acoustic Doppler current profilers from June to October 2014.Results show that the M_(2) tide is strong and shows a barotropic feature,whereas the O_(1) tide is much weaker.Based on the M_(2) tidal currents,the eddy viscosity in the bottom Ekman boundary layer is estimated with three schemes.The estimated eddy viscosity values vary within 10^(-4)–10^(-2)m^(2) s^(−1),reaching a maximum at approximately 5 m height from the bottom and decreasing exponentially with the height at all three stations.Moreover,the shear production of turbulent kinetic energy is calculated to quantify the mixing induced by different tidal constituents.The results show that the shear production of the M_(2) tide is much stronger than that of the O_(1) tide and shows a bottom intensified feature.

Different Turbulent Regimes and Vertical Turbulence Structures of the Urban Nocturnal Stable Boundary Layer

Turbulence in the nocturnal boundary layer(NBL)is still not well characterized,especially over complex underlying surfaces.Herein,gradient tower data and eddy covariance data collected by the Beijing 325-m tower were used to better understand the differentiating characteristics of turbulence regimes and vertical turbulence structure of urban the NBL.As for heights above the urban canopy layer(UCL),the relationship between turbulence velocity scale(VTKE)and wind speed(V)was consistent with the“HOckey-Stick”(HOST)theory proposed for a relatively flat area.Four regimes have been identified according to urban nocturnal stable boundary layer.Regime 1 occurs where local shear plays a leading role for weak turbulence under the constraint that the wind speed V<VT(threshold wind speed).Regime 2 is determined by the existence of strong turbulence that occurs when V>VT and is mainly driven by bulk shear.Regime 3 is identified by the existence of moderate turbulence when upside-down turbulence sporadic bursts occur in the presence of otherwise weak turbulence.Regime 4 is identified as buoyancy turbulence,when V>VT,and the turbulence regime is affected by a combination of local wind shear,bulk shear and buoyancy turbulence.The turbulence activities demonstrated a weak thermal stratification dependency in regime 1,for which within the UCL,the turbulence intensity was strongly affected by local wind shear when V<VT.This study further showed typical examples of different stable boundary layers and the variations between turbulence regimes by analyzing the evolution of wind vectors.Partly because of the influence of large-scale motions,the power spectral density of vertical velocity for upsidedown structure showed an increase at low frequencies.The upside-down structures were also characterized by the highest frequency of the stable stratifications in the higher layer.

Three-Dimensional Direct Numerical Simulation of Flow past A Near-Wall Circular Cylinder:Combined Effects of Gap Ratio and Boundary Layer Thickness on Flow Profiles and Pressure Distribution

Three-dimensional direct numerical simulations of the wake flow downstream of a near-wall circular cylinder at different gap ratios and boundary layer thicknesses are carried out by using the iterative immersed boundary method.The non-dimensional gap between the cylinder and the wall,G/D=0.2,0.6 and 1.0,the non-dimensional boundary layer thickness,δ/D=0.0,0.7 and 1.6,the Reynolds number,Re=350,and the aspect ratio of the cylinder,L/D=25are adopted.High-resolution visualizations of the complex vortex structures at differentδ/D and G/D are presented.The transition of the streamwise vortex mode,the combined effects ofδ/D and G/D on the flow statistics,the pressure and shear stress distribution and the hydrodynamic forces are analyzed.Results show that with decreasing G/D and increasingδ/D,the gap flow and its vortex-shedding are significantly weakened,together with an elongated wake and an enlarged low-velocity area near the wall,leading to the wake mode transition from the two-sided to one-sided vortex-shedding.Different relative positions of the cylinder regarding the boundary layer alter the flow features of the shear layers.With an increase inδ/D,the front stagnation point shifts to the upper surface,and the distance between the flow divergence point and the maximum pressure position increases.The mean drag coefficient and r.m.s.values of drag and lift coefficients decrease with a decrease in G/D and an increase inδ/D,while the mean lift coefficient increases with decreasing G/D but decreases with increasingδ/D.

Characteristics of atmospheric boundary layer structure and its influencing factors under different sea and land positions in Europe

This study identifies quantitatively the dominant contributions of meteorological factors on the development of the boundary layer heights(BLH)in the European region,based on 32 years(1990-2021)of radiosonde observations.The spatial variability of the BLH is further discussed by location,by classifying recording stations as inland,coastal,or bay.We find that the BLH in Europe varies considerably from day to night and with the seasons.Nighttime BLH is higher in winter and lower in summer,with the highest BLH recorded at coastal stations.Daytime BLH at coastal stations shows a bimodal structure with peaks in spring and autumn;at inland and bay stations,daytime BLH is lower in winter and higher in summer.The daily amplitudes of BLH at the inland and bay stations are stronger than those at coastal stations.Based on our multiple linear regression analysis and our decoupling analysis of temperature and specific humidity,we report that the development of the nighttime BLH at all types of stations is strongly dominated by the variations of surface wind speed(and,at coastal stations,wind directions).The main contributors to daytime BLH are the near-surface temperature variability at most coastal and inland stations,and,at most bay stations,the variation of the near-surface specific humidity.

Effect of local wall temperature on hypersonic boundary layer stability and transition

Wall temperature significantly affects stability and receptivity of the boundary layer. Changing the wall temperature locally may therefore be an effective laminar flow control technique. However, the situation is complicated when the wall temperature distribution is nonuniform, and researchers have experimentally found that local wall cooling may delay the onset of transition. We attempt to clarify the physical mechanisms whereby the local wall temperature affects the transition and the stability of a hypersonic boundary layer. A numerical investigation of the disturbance evolution in a Mach-6 sharp cone boundary layer with local wall heating or cooling is conducted. Direct numerical simulation(DNS) is performed for the single-frequency and broadband disturbance evolution caused by random forcing. We vary the local wall temperature and the location of heating/cooling, and then use the eNmethod to estimate the transition onset. Our results show that local wall cooling amplifies high-frequency unstable waves while stabilizing low-frequency unstable waves, with local heating amplifying all unstable waves locally. The disturbance amplitude and second-mode peak frequency obtained by DNS agree well with the previous experimental results. Local cooling/heating has a dual effect on the stability of the hypersonic boundary layer. For local cooling, while it effectively inhibits the growth of the low-frequency unstable waves that dominate the transition downstream, it also further destabilizes the downstream flow. In addition, while upstream cooling can delay the transition, excessive cooling may promote it;local heating always slightly promotes the transition.Finally, recommendations are given for practical engineering applications based on the present results.

Turbulent boundary layer control with DBD plasma actuators

The flat-plate turbulent boundary layer at Reτ=1140 is manipulated using a spanwise array of bidirectional dielectric barrier discharge(DBD)plasma actuators.Based on the features of no moving mechanical parts in the DBD plasma control technology and hot-wire anemometer velocity measurements,a novel convenient method of local drag reduction(DR)measurement is proposed by measuring the single-point velocity within the linear region of the viscous sublayer.We analyze the premise of using the method,and the maximum effective measurement range of-73.1%<DR<42.2%is obtained according to the experimental environment in this work.The local drag decreases downstream of the center of two adjacent upper electrodes and increases downstream of the upper electrodes.The magnitude of the local DR increases with increasing voltage and decreases as it moves away from the actuators.For the spanwise position in between,the streamwise distribution of the local DR is very dependent on the voltage.The variable-interval time-average detection results reveal that all bursting intensities are reduced compared to the baseline,and the amount of reduction is comparable to the absolute values of the local DR.Compared with previous results,we infer that the control mechanism is that many meandering streaks are combined together into single stabilized streaks.

NUMERICAL SIMULATION AND ANALYSIS OF HORSESHOE-SHAPED VORTEX IN NEAR-WALL REGION OF TURBULENT BOUNDARY LAYER

The low-Reynolds-number full developed turbulent flow in channels is simulated using large eddy simulation（LES）method with the preconditioned algorithm and the dynamic subgrid-scale model,with a given disturbance in inlet boundary,after a short development section.The inlet Reynolds number based on momentum thickness is 670.The computed results show good agreement with direct numerical simulation（DNS）,which include root mean square fluctuated velocity distribution and average velocity distribution.It is also found that the staggered phenomenon of the coherent structures is caused by sub-harmonic.The results clearly show the formation and evolution of horseshoe vortex in the turbulent boundary layer,including horseshoe vortex structure with a pair of streamwise vortexes and one-side leg of horseshoe vortex.Based on the results,the development of the horseshoe-shaped coherent structures is analyzed in turbulent boundary layer.

Boundary Layer Characteristics and Numerical Simulation Analysis of Winter Dense Fog in Nanjing

Based on the boundary layer data of winter dense fog in 2007 from Nanjing University of Information Science & Technology,the profile characteristics of temperature,wind direction,wind speed and humidity in a dense fog weather on December 13-14 in 2007 were analyzed,as well as their evolution laws in the formation and dispersion of fog,and the boundary layer characteristics of winter dense fog in Nanjing were revealed,while the development of fog was simulated by means of mesoscale numerical model.The results showed that the formation and dispersion of fog was greatly affected by inversion and humidity in the surface layer,and the wind direction in the surface layer also had effect on the formation and dispersion of advection fog.Mesoscale numerical model could preferably simulate the evolution of temperature,humidity,vertical speed in the development of fog,and the simulation of water vapor content in the fog could forecast the formation and dispersion of fog.

Test Research on the Characteristics of Pollution Boundary Layer in the River Outlet Area of Grand Liao River

By using the meteorological data in the pollution boundary layer which was observed in two ground observation sites:coast and land in the river outlet area of Grand Liao River during January-February in 2007,the daily change characteristics of pollute boundary layer in winter in the area were discussed. The results showed that the pollute boundary layer in the river outlet area of Grand Liao River was affected by the sea and land. In the certain weather condition,maybe the sea-land breeze appeared in the low altitude which was below 200 m in the coastal zone. The stability change in the different height in the coastal zone was more stable than in the land zone,and the wind field change in the area was mainly in 300 m low altitude. At night,the temperature inversion often appears in the area,and the thickness of temperature inversion layer is stably during 200-300 m. The thermal internal boundary layer penetrated deeply into the land about 10 km,and the height could reach 800 m. The atmospheric diffusion ability in the coastal area was weaker and stronger in the land area.

Large Eddy Simulation and Study of the Urban Boundary Layer

Based on a pseudo-spectral large eddy simulation (LES) model, an LES model with an anisotropy turbulent kinetic energy (TKE) closure model and an explicit multi-stage third-order Runge-Kutta scheme is established. The modeling and analysis show that the LES model can simulate the planetary boundary layer (PBL) with a uniform underlying surface under various stratifications very well. Then, similar to the description of a forest canopy, the drag term on momentum and the production term of TKE by subgrid city buildings are introduced into the LES equations to account for the area-averaged effect of the subgrid urban canopy elements and to simulate the meteorological fields of the urban boundary layer (UBL). Numerical experiments and comparison analysis show that: (1) the result from the LES of the UBL with a proposed formula for the drag coefficient is consistent and comparable with that from wind tunnel experiments and an urban subdomain scale model; (2) due to the effect of urban buildings, the wind velocity near the canopy is decreased, turbulence is intensified, TKE, variance, and momentum flux are increased, the momentum and heat flux at the top of the PBL are increased, and the development of the PBL is quickened; (3) the height of the roughness sublayer (RS) of the actual city buildings is the maximum building height (1.5-3 times the mean building height), and a constant flux layer (CFL) exists in the lower part of the UBL.

Selective enhancement of oblique waves caused by finite amplitude second mode in supersonic boundary layer

Nonlinear interactions of the two-dimensional （2D） second mode with oblique modes are studied numerically in a Mach 6.0 fiat-plate boundary layer, focusing on its selective enhancement effect on amplification of different oblique waves. Evolution of oblique modes with various frequencies and spanwise wavenumbers in the presence of 2D second mode is simulated successively, using a modified parabolized stability equation （PSE） method, which is able to simulate interaction of two modes with different frequen- cies efficiently. Numerical results show that oblique modes in a broad band of frequencies and spanwise wavenumbers can be enhanced by the finite amplitude 2D second mode instability wave. The enhancement effect is accomplished by interaction of the 2D second mode, the oblique mode, and a forced mode with difference frequency. Two types of oblique modes are found to be more amplified, i.e., oblique modes with frequency close to that of the 2D second mode and low-frequency first mode oblique waves. Each of them may correspond to one type of transition routes found in transition experiments. The spanwise wavenumber of the oblique wave preferred by the nonlinear interaction is also determined by numerical simulations.

DIRECT NUMERICAL SIMULATION OF VIBRATIONINDUCED T-S WAVE ON BOUNDARY LAYER WALL

The spatial growth of the disturbance in the boundary layer is directly numerically simulated, and the receptivity of the Blasius basic flow to the local two-dimensional （2-D） sustainable micro-vibration is investigated. Results show that the disturbance velocity presents the sine vibration features with the change of time, and the vibration period is identical to the vibration of the local wall. The disturbance velocity presents the fluctuation feature downstream, and the streamwise wave length approximates to the results from the Orr-Sommerfeld equation （OSE）. The growth rate from direct numerical simulation（DNS） is a little greater than that from OSE, and their trends are almost consistent. Under the condition of Re= 2 800, the disturbance amplitude gradually grows in the given computational region with the period T=30. However, it firstly increases and then decreases with the period T= 20. The disturbance harmonic of the former is obviously larger than that of the latter. The maximum streamwise and vertical disturbance velocities from DNS do not fully coincide with those from OSE at the vicinity of the local vibration wall, but coincide well with the former when they travel downstream. The 2-D disturbance induced by the local micro-vibration represents the form of Tollmien-Schlichting （T-S） wave on the boundary layer.

STABILITY OF 3-D DISTURBANCE WAVES ON BOUNDARY LAYER WITH DISTRIBUTED LOCAL ROUGHNESS

The numerical solution of the stable basic flow on a 3-D boundary layer is obtained by using local ejection, local suction, and combination of local ejection and suction to simulate the local rough wall. The evolution of 3-D disturbance T-S wave is studied in the spatial processes, and the effects of form and distribution structure of local roughness on the growth rate of the 3-D disturbance wave and the flow stability are discussed. Numerical results show that the growth of the disturbance wave and the form of vortices are accelerated by the 3-D local roughness. The modification of basic flow owing to the evolvement of the finite amplitude disturbance wave and the existence of spanwise velocity induced by the 3-D local roughness affects the stability of boundary layer. Propagation direction and phase of the disturbance wave shift obviously for the 3-D local roughness of the wall. The flow stability characteristics change if the form of the 2-D local roughness varies.

NONLINEAR EFFECT OF LOCAL ROUGH WALL ON LARGE EDDY STRUCTURE IN BOUNDARY LAYER

Based on the building of a theoretical model for the large eddy structure, the nonlinear effect of the local rough wall on the large eddy structure in the boundary layer is studied by direct numerical simulation. Numerical results show that factors of the local rough feature, the distributing structure and the intensity, etc. play an important role in the evolution of the large eddy structure in the boundary layer.

NONLINEAR EVOLUTION OF TWO-DIMENSIONAL T-S WAVE IN LOCALIZED ROUGH BOUNDARY LAYER

Numerical simulations are carried out to investigate the mechanism of the nonlinear evolution of two- dimensional （2-D） Tollmien-Schlichting （T-S） wave on a localized rough boundary layer. The three-dimensional （3-D） numerical solution of a base flow on a boundary layer is obtained for the localized rough wall with the local- ized ejection, the localized suction and the combination of ejection and suction. Based on numerical simulations, the processes of stable and the most instable nonlinear evolution of the 2-D disturbance T-S wave are studied. The effects of the form on the localized roughness, the intensity, and the distribution structure on the nonlinear evolution of 2-D T-S wave and the growth rate are discussed. Results show that the basic flow induced by the lo- calized rough wall is a key factor causing the fast growth of the disturbance wave. Due to the change of the aver- age flow profile and the existence of the spanwise velocity, the localized rough wall enhances the instability of the flow. Consequently, the instable region of the neutral curve is enlargened, and the maximnum growth rate of the 2-D T-S wave is increased. In the process of the nonlinear evolution of 2-D disturbance T-S wave, with the in- crease of the nonlinear interaction, the most instable 2-D disturbance wave triggers the appearance of the 3-D dis- turbance wave and the high-frequency harmonic wave. Its streamwise wave number and the frequency are the same as those of 2-D disturbance wave. The spanwise velocity can excite the growth of the 2-D disturbance wave, the instability of 2-D wave, the formation of the streamwise vortex, and the generation of 3-D disturbance wave. Simulation results agree well with experimental results.

Similarity solutions of Prandtl mixing length modelled two dimensional turbulent boundary layer equations

The exact similarity solutions of two dimensional laminar boundary layer were obtained by Blasius in 1908,however,for two dimensional turbulent boundary layers,no Blasius type similarity solutions(special exact solutions)have ever been found.In the light of Blasius’pioneer works,we extend Blasius similarity transformation to the two dimensional turbulent boundary layers,and for a special case of flow modelled by Prandtl mixing-length,we successfully transform the two dimensional turbulent boundary layers partial differential equations into a single ordinary differential equation.The ordinary differential equation is numerically solved and some useful quantities are produced.For numerical calculations,a complete Maple code is provided.

Development of a Single-Column Model in RegCM4 and Its Preliminary Application for Evaluating PBL Schemes in Simulating the Dry Convection Boundary Layer

A single-column model(SCM)is developed in the regional climate model RegCM4.The evolution of a dry convection boundary layer(DCBL)is used to evaluate this SCM.Moreover,four planetary boundary layer(PBL)schemes,namely the Holtslag-Boville scheme(HB),Yonsei University scheme(YSU),and two University of Washington schemes(UW01,Grenier-Bretherton-Mc Caa scheme and UW09,Bretherton-Park scheme),are compared by using the SCM approach.A large-eddy simulation(LES)of the DCBL is performed as a benchmark to examine how well a PBL parameterization scheme reproduces the LES results,and several diagnostic outputs are compared to evaluate the schemes.The results show that the SCM is properly constructed.In general,with the DCBL case,the YSU scheme performs best for reproducing the LES results,which include well-mixed features and vertical sensible heat fluxes;the simulated wind speed,turbulent kinetic energy,entrainment flux,and height of the entrainment zone are all underestimated in the UW09;the UW01 has all those biases of the UW09 but larger,and the simulated potential temperature is not well mixed;the HB is the least skillful scheme,by which the PBL height,entrainment flux,height of the entrainment zone,and the vertical gradients within the mixed layer are all overestimated,and a inversion layer near the top of the surface layer is wrongly simulated.Although more cases and further testing are required,these simulations show encouraging results towards the use of this SCM framework for evaluating the simulated physical processes by the RegCM4.