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.

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.

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.

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.

Wavelet Multi-Resolution Interpolation Galerkin Method for Linear Singularly Perturbed Boundary Value Problems

In this study,a wavelet multi-resolution interpolation Galerkin method(WMIGM)is proposed to solve linear singularly perturbed boundary value problems.Unlike conventional wavelet schemes,the proposed algorithm can be readily extended to special node generation techniques,such as the Shishkin node.Such a wavelet method allows a high degree of local refinement of the nodal distribution to efficiently capture localized steep gradients.All the shape functions possess the Kronecker delta property,making the imposition of boundary conditions as easy as that in the finite element method.Four numerical examples are studied to demonstrate the validity and accuracy of the proposedwavelet method.The results showthat the use ofmodified Shishkin nodes can significantly reduce numerical oscillation near the boundary layer.Compared with many other methods,the proposed method possesses satisfactory accuracy and efficiency.The theoretical and numerical results demonstrate that the order of theε-uniform convergence of this wavelet method can reach 5.

Thermodynamics-directed bulk/grain-boundary engineering for superior electrochemical durability of Ni-rich cathode

Introducing high-valence Ta element is an essential strategy for addressing the structu ral deterioration of the Ni-rich LiNi_(1-x-y)Co_(x)Mn_(y)O_(2)(NCM)cathode,but the enlarged Li/Ni cation mixing leads to the inferior rate capability originating from the hindered Li~+migration.Note that the non-magnetic Ti~(4+)ion can suppress Li/Ni disorder by removing the magnetic frustration in the transition metal layer.However,it is still challenging to directionally design expected Ta/Ti dual-modification,resulting from the complexity of the elemental distribution and the uncertainty of in-situ formed coating compounds by introducing foreign elements.Herein,a LiTaO_3 grain boundary(GB)coating and bulk Ti-doping have been successfully achieved in LiNi_(0.834)Co_(0.11)Mn_(0.056)O_(2) cathode by thermodynamic guidance,in which the structural formation energy and interfacial binding energy are employed to predict the elemental diffusion discrepancy and thermodynamically stable coating compounds.Thanks to the coupling effect of strengthened structural/interfacial stability and improved Li~+diffusion kinetics by simultaneous bulk/GB engineering,the Ta/Ti-NCM cathode exhibits outstanding capacity retention,reaching 91.1%after 400 cycles at 1 C.This elaborate work contributes valuable insights into rational dual-modification engineering from a thermodynamic perspective for maximizing the electrochemical performances of NCM cathodes.

Boundary Conditions for Momentum and Vorticity at an Interface between Two Fluids

Boundary conditions for momentum and vorticity have been precisely derived, paying attention to the physical meaning of each mathematical expression of terms rigorously obtained from the basic equations: Navier-Stokes equation and the equation of vorticity transport. It has been shown first that a contribution of fluid molecules crossing over a conceptual surface moving with fluid velocity due to their fluctuating motion is essentially important to understanding transport phenomena of momentum and vorticity. A notion of surface layers, which are thin layers at both sides of an interface, has been introduced next to elucidate the transporting mechanism of momentum and vorticity from one phase to the other at an interface through which no fluid molecules are crossing over. A fact that a size of δV, in which reliable values of density, momentum, and velocity of fluid are respectively defined as a volume-averaged mass of fluid molecules, a volume-averaged momentum of fluid molecules and a mass-averaged velocity of fluid molecules, is not infinitesimal but finite has been one of the key factors leading to the boundary conditions for vorticity at an interface between two fluids. The most distinguished characteristics of the boundary conditions derived here are the zero-value conditions for a normal component of momentum flux and tangential components of vorticity flux, at an interface.

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.

Perfectly matched layer-absorbing boundary condition for finite-element time-domain modeling of elastic wave equations

The perfectly matched layer （PML） is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-element time-domain numerical modeling of elastic wave equation. However, the finite-element time-domain scheme is based on the second- order wave equation in displacement formulation. Thus, the first-order PML in velocity-stress formulation cannot be directly applied to this scheme. In this article, we derive the finite- element matrix equations of second-order PML in displacement formulation, and accomplish the implementation of PML in finite-element time-domain modeling of elastic wave equation. The PML has an approximate zero reflection coefficients for bulk and surface waves in the finite-element modeling of P-SV and SH wave propagation in the 2D homogeneous elastic media. The numerical experiments using a two-layer model with irregular topography validate the efficiency of PML in the modeling of seismic wave propagation in geological models with complex structures and heterogeneous media.

Numerical Simulation of Two-Dimensional Shock/Boundary-Layer Interaction between a Rocket and Booster

A two-dimensional Reynolds-averaged Navier-Stokes solver is applied to analyze the aerodynamic behavior of the Shock/Boundary-Layer interaction of rocket with a boosted The K-ε turbulence model and a finite volume method in a unstructured body-fitted curvilinear coordinates have been used. The results indicate that the separation and the reattachment occur in the Boundary-Layer of the main rocket because of the shock interaction. The shape of the booster nose effects the flow field obviously. In the case of the hemisphere booster nose the pressure has complicate distributions and the separation is very clear. The distance between the booster and main rocket has the evident effect on the flow field. If the distance is smaller the pressure coefficient is bigger the separation zone even the separation bubble occurs.

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.

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.