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 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.

Linear spatial instability analysis in 3D boundary layers using plane-marching 3D-LPSE

It is widely accepted that a robust and efficient method to compute the linear spatial amplified rate ought to be developed in three-dimensional （3D） boundary layers to predict the transition with the e^N method, especially when the boundary layer varies significantly in the spanwise direction. The 3D-linear parabolized stability equation （3D- LPSE） approach, a 3D extension of the two-dimensional LPSE （2D-LPSE）, is developed with a plane-marching procedure for investigating the instability of a 3D boundary layer with a significant spanwise variation. The method is suitable for a full Mach number region, and is validated by computing the unstable modes in 2D and 3D boundary layers, in both global and local instability problems. The predictions are in better agreement with the ones of the direct numerical simulation （DNS） rather than a 2D-eigenvalue problem （EVP） procedure. These results suggest that the plane-marching 3D-LPSE approach is a robust, efficient, and accurate choice for the local and global instability analysis in 2D and 3D boundary layers for all free-stream Mach numbers.

Linear stability theory with the equivalent spanwise wavenumber correction in 3D boundary layers

The prediction on small disturbance propagation in complex three-dimensional(3D) boundary layers is of great significance in transition prediction methodology, especially in the aircraft design. In this paper, the linear stability theory(LST) with the equivalent spanwise wavenumber correction(ESWC) is proposed in order to accurately predict the linear evolution of a disturbance in a kind of boundary layer flow with a vital variation in the spanwise direction. The LST with the ESWC takes not only the scale of the mean flow with the significant variation but also the wavenumber evolution of the disturbance itself. Compared with the conventional LST, the results obtained by the new method are in excellent agreement with those of the numerical simulations. The LST with the ESWC is an effective method on the prediction of the disturbance evolution in 3D boundary layers, which improves the prediction of the LST in the applications to complex 3D boundary layers greatly.

3-D Scattering of Obliquely Incident Plane P Waves by Alluvial Valley Embedded in Layered Half-Space

The indirect boundary element method(IBEM) was established to solve the problem of 3-D seismic responses of 2-D topographies,by calculating the free-field responses with the direct-stiffness method and simulating the scattering wave fields with the dynamic Green's functions of moving distributed loads.The proposed method yields accurate results,because the 3-D dynamic stiffness matrixes used are exact and the fictitious moving distributed loads can be acted directly on the interface between the alluvial valley and the layered half-space without singularity.The comparison with the published methods verifies the validity of the proposed method.And the numerical analyses are performed to give some beneficial conclusions.The study shows that 3-D scattering by an alluvial valley is essentially different from the 2-D case,and that the presence of soil layer affects not only the amplitude value of surface displacements but also the distribution of surface displacements.

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.

3-D seismic interpretation of stratigraphic and structural features in the Upper Jurassic to Lower Cretaceous sequence of the Gullfaks Field,Norwegian North Sea:A case study of reservoir development

The 3-D seismic dataset is a key tool to analyze and understand the mechanism of structural and stratigraphic hydrocarbon(HC)trapping in the subsurface.Conventionally used subsurface seismic characterization methods for fractures are based on the theory of effective anisotropy medium.The aim of this work is to improve the structural images with dense sampling of 3-D survey to evaluate structural and stratigraphic models for reservoir development to predict reservoir quality.The present study of the Gullfaks Field,located in the Norwegian North Sea Gullfaks sector,identifies the shallowest structural elements.The steepness of westward structural dip decreases eastward during the Upper Jurassic to Lower Cretaceous deposition.Reservoir sands consist of the Middle Jurassic deltaic deposits and Lower Jurassic fluvial channel and delta plain deposits.Sediment supply steadily prevails on sea-level rise and the succession displays a regressive trend indicated by a good continuous stacking pattern.The key factor for the development of reservoirs in the Gullfaks Field is fault transmissibility with spatially distributed pressure.The majority of mapped faults with sand-to-sand contacts are non-sealing,which provide restriction for the HC flow between the fault blocks.The traps for HC accumulation occur between the post-rift and syn-rift strata,i.e.antiform set by extensional system,unconformity trap at the top of syndeposition,and structural trap due to normal faults.Overall reservoir quality in the studied area is generally excellent with average 35%porosity and permeability in the Darcy range.Our findings are useful to better understand the development of siliciclastic reservoirs in similar geological settings worldwide.

Improved treatments for evaluating horizontal magnetic components through the 3-D FDM in E-polarization induction problems

To improve the accuracy of the numerical evaluation through the 3-D finite difference method, the surface boundary conditions are added to modify the old program. The author has tested the new program by making calculations for the model constructed by Wanamaker, et al (1984). The comparison between the numerical results obtained from this paper and those by Wannamaker, et al (1984) indicates that a pronounced improvement is realized in the evaluation of the horizontal magnetic components. Moreover, better calculations for the vertical magnetic components are also obtainable by using the new program.

氧化温度对空气中一次烧成SrTiO_(3)晶界层电容器的显微结构及电性能的影响

采用空气中的一次烧成工艺,研究了氧化温度对SrTiO_3.晶界层电容器试样的介电常数ε、绝缘电阻率p及介质损耗tgδ的影响;借助于TEM及STEM对试样进行了分析,观察到试样中有直线型和锯齿型两种不同类型的晶界及存在于多个晶粒交汇处的Ti_nO_(2n-1)第二晶相;提出了不同氧化温度下的晶界结构模型,较好地解释了试样介电常数随氧化温度的降低而单调增大的变化规律.

Exponential Convergence Theory of the Multipole and Local Expansions for the 3-D Laplace Equation in Layered Media

In this paper,we establish the exponential convergence theory for the multipole and local expansions,shifting and translation operators for the Green's function of 3-dimensional Laplace equation in layered media.An immediate application of the theory is to ensure the exponential convergence of the FMM which has been shown by the numerical results reported in[27].As the Green's function in layered media consists of free space and reaction field components and the theory for the free space components is well known,this paper will focus on the analysis for the reaction components.We first prove that the density functions in the integral representations of the reaction components are analytic and bounded in the right half complex wave number plane.Then,by using the Cagniard-de Hoop transform and contour deformations,estimates for the remainder terms of the truncated expansions are given,and,as a result,the exponential convergence for the expansions and translation operators is proven.

边界层方案与垂直混合对京津冀地区O_(3)模拟的影响研究

采用WRF-Chem模式中的3种边界层方案YSU、MYJ和ACM2对2019年6月京津冀及周边地区典型O_(3)污染月份开展模拟研究.详细对比了各方案对地面气象要素、NO_(2)和O_(3)浓度时空分布,以及温湿风要素和O_(3)浓度垂直分布的模拟效果.结果表明:3种方案对地面气象要素的时空分布和温湿风要素的垂直变化模拟较为合理.MYJ方案模拟地面气象要素整体效果最佳.各方案对边界层高度的日变化特征模拟较好,相关系数为0.58~0.69,但存在白天偏高、夜间偏低的现象,YSU方案相比效果最佳.3种边界层方案对NO2浓度模拟普遍高估,而O_(3)模拟结果则出现低估.白天模拟偏差较小而夜间偏差较显著.模拟最佳的是ACM2,其次为YSU和MYJ.3种方案均较好地模拟出了O_(3)的垂直分布特征,但整体低估了O_(3)浓度.对上午O_(3)垂直分布的模拟差异较下午更为明显.此外,基于YSU方案设置了3个敏感实验,通过调整化学模块所用的湍流扩散系数阈值,对比分析了垂直混合过程改变对O_(3)浓度模拟的影响,模拟的变化只反映由于边界层的垂直混合过程改变造成的污染差异,而不是由于热动力场的调整造成的变化.模拟结果表明3个方案均可改善区域上地面NO_(2)和O_(3)的模拟性能,尤其是对原3种边界层方案模拟O_(3)均明显低估的华北平原地区提升效果最显著,平均偏差降低了23.7%.在垂直方向上,湍流扩散系数阈值的调整增加了早间近地面模拟的O_(3)浓度,改善了模拟偏低的现象,但同时增大了高层O_(3)浓度的负偏差.敏感性方案显著改善了夜间的模拟,白天则并不明显.这些结果显示出湍流扩散系数对O_(3)垂直混合的重要影响.因此,改进湍流扩散系数的参数化对O_(3)模拟是必要的.

Scattering of seismic waves by three-dimensional large-scale hill topography simulated by a fast parallel IBEM

To solve seismic wave scattering by a large-scale three-dimensional(3-D) hill topography, a fast parallel indirect boundary element method(IBEM) is developed by proposing a new construction method for the wave field, modifying the generalized minimum residual(GMRES) algorithm and constructing an Open MP plus MPI parallel model. The validations of accuracy and efficiency show that this method can solve 3-D seismic response of a large-scale hill topography for broadband waves, and overcome the weakness of large storage and low efficiency of the traditional IBEM. Based on this new algorithm architecture, taking the broadband scattering of plane SV waves by a large-scale Gaussian-shaped hill of thousands-meters height as an example, the influence of several important parameters is investigated, including the incident frequency, the incident angle and the height-width and length-width ratio of the hill. The numerical results illustrate that the amplification effect on the ground motion by a near-hemispherical hill is more significant than the narrow hill. For low-frequency waves, the scattering effect of the higher hill is more pronounced, and there is only a single peak near the top of the hill. However, for high-frequency waves, rapid spatial variation of displacement amplitude appears on the hill surface.

Influence of Vertical Eddy Diffusivity Parameterization on Daily and Monthly Mean Concentrations of O3 and NOy

Two parameterization schemes for vertical eddy diffusivity were utilized to investigate their impacts on both the daily and monthly mean concentrations of ozone and NOy, which are the major fractions of the sum of all reactive nitrogen species, i.e., NOy=NO＋NO2＋HNO3＋PAN. Simulations indicate that great changes in the vertical diffusivity usually occur within the planetary boundary layer （PBL）. Daily and monthly mean concentrations of NOy are much more sensitive to changes in the vertical diffusivity than those of ozone and ozone and NOy levels only at or in （relatively） clean sites and areas, where long-range transport plays a crucial role, display roughly equivalent sensitivity. The results strongly suggest that a widely-accepted parameterization scheme be selected and the refinement of the model＇s vertical resolution in the PBL be required, even for regional and long-term studies, and ozone only being examined in an effort to judge the model＇s performance be unreliable, and NOy be included for model evaluations.

后续处理对SrTiO_(3)基晶界层电容器绝缘电阻的影响

采用二步法制备SrTiO_(3)晶界层电容器,并对其进行后续热、电和液氮处理,研究处理前后电容器电学性能的变化.实验结果表明,在50 V直流电压和200℃条件下对SrTiO_(3)晶界层电容器进行后续快速退火和液氮处理后,其介电常数和介电损耗在基本保持不变的情况下,其绝缘电阻值可得到大幅提升,从最初30 GΩ上升至200 GΩ.通过处理,最后可获得平均介电常数为30000,损耗为0.003,绝缘电阻(50 V测量)为200 GΩ的高性能SrTiO_(3)晶界层电容器.

FAST APPROXIMATE INVERSION OF MAGNETOTELLURIC(FAIMT) DATA FOR A HORIZONTALLY LAYERED EARTH WITH NEAR-SURFACE INHOMOGENEITIES

This paper discusses use of approximations and the Integral Mean Value Theorem to show that 6 coefficients approximately describe the distortions of near surface inhomogeneities on the MT field of a horizontally layered earth model. When these 6 coefficients are considered together with those of the magnetic field of a horizontally layered earth model,the analytic and approximate wave impedance equations can be derived for the MT response of a horizontally layered earth model with near-surface 2-D and 3-D inhomogeneities. These approximate wave impedance equations are used with inverted MT data for 2-D and 3-D forward modelling. Although these 6 coefficients cannot be determined before inversion,initial estimates can be used. The 6 coefficients and the asistivity and thickness of each layer of a horizontally layered earth can be obtained by using published inversion methods. The 6 coefficients give important informaion (depths and resistivities) on the near-surface inhomogenelties.The authors inverted 2-D and 3-D theoretical models for Fast Approximate Inversion of Magnetotelluric (FAIMT) data for a horizontally layered earth with near-surface inhomogeneities compares favorably with traditional invrsion methods, especially for inverting regional or basin structures. This method simplifies computation and gives a reasonable 1 -D geological model with fewer nonuniquenas problems.

Modeling of Large ELM Suppressions for High Confinement Plasma in D Ⅲ-D

The H-mode discharges with high edge pressure gradients are expected for the economic feasibility of future fusion reactors. However, the high edge pressure gradients easily produce ELM instability , which generally can expel large, heat and particle loading to the divertor targets. These ELMs limit the core plasma performance and reduce the lifetime of divertor target plates. The transports of heat and particles outward across plasma boundary are useful to control density and impurity profiles for achieving steady state, high performance plasmas. Consequently, any technique to eliminate or mitigate large fast ELM impulses must replace the transient heat and particle transports with another slow process. Such a technique is high priority for a burning plasma device such as ITER.

Analysis of 3-D Frictional Contact Mechanics Problems by a Boundary Element Method

The development of two boundary element algorithms for solving 3-D, frictional, and linear elastostatic contact problems is reported in this paper. The algorithms employ nonconforming discreti- zations for solving 3-D boundary element models, which provide much needed flexibility in the bound- ary element modeling for 3-D contact problems. These algorithms are implemented in a new 3-D boundary element code and verified using several examples. For the numerical examples studied, the results using the new boundary element algorithms match very well with the results using a commercial finite element code, and clearly demonstrate the feasibility of the new boundary element approach for 3-D contact analysis.

三维层状地基空沟主动隔振分析

基于薄层法在研究层状介质中波的传播问题的高效性、边界单元法处理无限域问题的精确性,结合二者的优点提出三维层状地基薄层法基本解答,建立了基于三维层状半空间薄层法位移基本解答的半解析动力边界元法。该方法可有效的分析多层场地的动力问题,解决土-结构动力相互作用问题。同时分别对粘弹性上软下硬地基及上硬下软地基的三维空沟主动隔振进行了详细的分析。结果表明,两种情况下采用空沟屏障隔振均可以取得一定的隔振效果;同时,地基分层参数对空沟隔振体系的隔振效果影响显著。

山丘地形的陆面过程及边界层特征的模拟

将模式NP 89的陆面过程参数化方法应用到北京大学的三维复杂地形中尺度数值模式中 ,得到了一个较理想的三维陆面过程及边界层模式 ,利用这个改进的三维模式对 2 0km× 2 0km范围的山丘地形的陆面过程及边界层特征进行了数值模拟。模拟结果表明 ,由于地形阻挡所造成山后的湍流较山前强 ,进而造成近地面温度梯度和感热支出小 ,最终造成山后的温度比山前的温度明显偏高 ;而且随着山高的增加 ,这种现象更加明显 ,即该模式对山丘地形条件下的陆面过程和大气边界层特征具有较强的模拟能力 ;模拟结果合理 ,对研究过山气流形成机制、起伏地形大气边界层物理特征和污染物的扩散具有理论和应用价值。

用三分量热线探头测量平板湍流边界层

本文用三分量热线探头,同时测量了零压梯度平板湍流边界层三个方向的脉动速度分量;用数字信号处理技术研究了平板湍流边界层的统计特性,并对测量结果进行了讨论。