Numerical Investigation on Two-dimensional Boundary Layer Flow with Transition

As a basic problem in many engineering applications, transition from laminar to turbulence still remains a difficult problem in computational fluid dynamics （CFD）. A numerical study of one transitional flow in two-dimensional is conducted by Reynolds averaged numerical simulation （RANS） in this paper. Turbulence model plays a significant role in the complex flows＇ simulation, and four advanced turbulence models are evaluated. Numerical solution of frictional resistance coefficient is compared with the measured one in the transitional zone, which indicates that Wilcox （2006） k-ω model with correction is the best candidate. Comparisons of numerical and analytical solutions for dimensionless velocity show that averaged streamwise dimensionless velocity profiles correct the shape rapidly in transitional region. Furthermore, turbulence quantities such as turbulence kinetic energy, eddy viscosity, and Reynolds stress are also studied, which are helpful to learn the transition＇s behavior.

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.

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.

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.

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.

Topology Optimization of Sound-Absorbing Materials for Two-Dimensional Acoustic Problems Using Isogeometric Boundary Element Method

In this work,an acoustic topology optimizationmethod for structural surface design covered by porous materials is proposed.The analysis of acoustic problems is performed using the isogeometric boundary elementmethod.Taking the element density of porousmaterials as the design variable,the volume of porousmaterials as the constraint,and the minimum sound pressure or maximum scattered sound power as the design goal,the topology optimization is carried out by solid isotropic material with penalization(SIMP)method.To get a limpid 0–1 distribution,a smoothing Heaviside-like function is proposed.To obtain the gradient value of the objective function,a sensitivity analysis method based on the adjoint variable method(AVM)is proposed.To find the optimal solution,the optimization problems are solved by the method of moving asymptotes(MMA)based on gradient information.Numerical examples verify the effectiveness of the proposed topology optimization method in the optimization process of two-dimensional acoustic problems.Furthermore,the optimal distribution of sound-absorbingmaterials is highly frequency-dependent and usually needs to be performed within a frequency band.

Two-Dimensional Riemann Problems:Transonic Shock Waves and Free Boundary Problems

We are concerned with global solutions of multidimensional(M-D)Riemann problems for nonlinear hyperbolic systems of conservation laws,focusing on their global configurations and structures.We present some recent developments in the rigorous analysis of two-dimensional(2-D)Riemann problems involving transonic shock waves through several prototypes of hyperbolic systems of conservation laws and discuss some further M-D Riemann problems and related problems for nonlinear partial differential equations.In particular,we present four different 2-D Riemann problems through these prototypes of hyperbolic systems and show how these Riemann problems can be reformulated/solved as free boundary problems with transonic shock waves as free boundaries for the corresponding nonlinear conservation laws of mixed elliptic-hyperbolic type and related nonlinear partial differential equations.

FINITE ELEMENT METHOD FOR SOLVING TWO-DIMENSIONAL DIFFUSION-REACTION EQUATIONS OF BOUNDARY LAYER TYPE IN POROUS CATALYST PELLET

In this paper,finite element method(FEM)is used to solve two-dimensional diffu-sion-reaction equations of boundary layer type.This kind of equations are usually too complicatedand diffcult to be solved by applying the traditional methods used in chemical engineering becauseof the steep gradients of concentration and temperature.But,these difficulties are easy to be over-comed when the FEM is used.The integraded steps of solving this kind of problems by the FEMare presented in this paper.By applying the FEM to the two actual examples,the conclusion can bereached that the FEM has the advantages of simplicity and good accuracy.

NUMERICAL STUDY FOR RESONANT TRIAD INTERACTION AND SECONDARY INSTABILITY IN TWO-DIMENSIONAL SUPERSONIC BOUNDARY LAYER

By direct numerical simulation, a supersonic boundary layer was studied to see whether the mechanism for the generation of sub-harmonic waves, similar to those for the incompressible flows, existed in the process of laminar-turbulent transition. The results showed that mechanisms of both resonant triad and secondary instability did exist. Discussions were made on whether these two mechanisms are really important in laminarturbulent transition.

ANALYSIS OF THE SECONDARY STABILITY OF COMPRESSIBLE BOUNDARY LAYER FLOWS OVER A TWO-DIMENSIONAL PLATE

This paper used the Floquet's three _dimensional linear stability theory in the analysis of two_dimensional compress ible boundary layer, a set of stability equations is constructed, the effect of three dimensional linear small perturbation on the two_dimensional compressible boundary layer transition is studied, and the effect of coming flow Ma number on growth and development of the subharmonics is calculated. It can be seen from t he calculations, the effect caused by the interaction of two_dimensional and thr ee_dimensional perturbation waves on the development of two_dimensional compress ible laminar boundary layer.

Static response of functionally graded multilayered two-dimensional quasicrystal plates with mixed boundary conditions

The unusual properties of quasicrystals(QCs)have attracted tremendous attention from researchers.In this paper,a semi-analytical solution is presented for the static response of a functionally graded(FG)multilayered two-dimensional(2 D)decagonal QC rectangular plate with mixed boundary conditions.Based on the elastic theory of FG 2 D QCs,the state-space method is used to derive the state equations composed of partial differential along the thickness direction.Besides,the Fourier series expansion and the differential quadrature technique are utilized to simulate the simply supported boundary conditions and the mixed boundary conditions,respectively.Then,the propagator matrix which connects the field variables at the upper interface to those at the lower interface of any homogeneous layer can be derived based on the state equations.Combined with the interface continuity condition,the static response can be obtained by imposing the sinusoidal load on the top surfaces of laminates.Finally,the numerical examples are presented to verify the effectiveness of this method,and the results are very useful for the design and understanding of the characterization of FG QC materials in their applications to multilayered systems.

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.

Inversion of two-dimensional tidal open boundary conditions of M_2 constituent in the Bohai and Yellow Seas

Two-dimensional tidal open boundary conditions of the M2 constituent in the Bohai and Yellow Seas(BYS) have been estimated by assimilating T/P altimeter data.During inversion,independent point(IP) strategy was used,in which several IPs on the open boundary is assumed,values at these IPs can be optimized with an adjoint method,and those at other grid points are determined by linearly interpolating the values at IPs.The reasonability and feasibility of the model are tested by ideal twin experiments.In the practical experiment(PE) after assimilation,the cost function may reach 1% or less of its initial value.Mean absolute errors in amplitude and phase can be less than 5 cm and 5°,respectively,and the obtained co-chart can show the character of the M2 constituent in the BYS.The results of the PE indicate that using only two IPs on the open boundary can yield better simulated 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.

Room temperature synthesis of two-dimensional multi layer magnets based on α-Co^(Ⅱ) layered hydroxides

Research on two-dimensional(2D) materials is one of the most active fields in materials science and nanotechnology. Among the members of the 2D family, layered hydroxides(LHs) represent an exceptional case of study due to their unparalleled chemical versatility which allows the modulation of their physicochemical properties at will. Nowadays, LHs based on earth-abundant metals are key materials in the areas of energy storage and conversion, hybrid materials or magnetism. α-Co hydroxides(Simonkolleite-like structures) are promising phases with tuneable electronic and magnetic properties by ligand modification. However, even in the simple case of α-Co^(Ⅱ) hydroxychlorides, the preparation of well-defined large 2D crystals is not straightforward, hindering the development of fundamental studies. Herein, we present the synthesis of 2D hexagonal crystals with outstanding sizethickness relationship(diameter > 5 μm and thickness of 20 ± 7 nm) by a simple homogeneous synthesis taking place at room temperature. In structural terms, no differences are observed between our layered materials and those obtained hydrothermally. However, dynamic susceptibility measurements alert about different arrangements of the magnetic sublattices, which have been rationalized with structural DFT calculations. This work provides an extremely easy bottom-up method to obtain high-quality 2D crystals based on α-CoIIhydroxides,paving the way for the development of fundamental studies and applications.