Generalized nth-Order Perturbation Method Based on Loop Subdivision Surface Boundary Element Method for Three-Dimensional Broadband Structural Acoustic Uncertainty Analysis

In this paper,a generalized nth-order perturbation method based on the isogeometric boundary element method is proposed for the uncertainty analysis of broadband structural acoustic scattering problems.The Burton-Miller method is employed to solve the problem of non-unique solutions that may be encountered in the external acoustic field,and the nth-order discretization formulation of the boundary integral equation is derived.In addition,the computation of loop subdivision surfaces and the subdivision rules are introduced.In order to confirm the effectiveness of the algorithm,the computed results are contrasted and analyzed with the results under Monte Carlo simulations(MCs)through several numerical examples.

COMPACT FOUR-COMPONENT 2-D FDFD METHOD WITH EQUIVALENT SURFACE IMPEDANCE BOUNDARY CONDITION FOR MULTILAYER METAL-COATED WAVEGUIDE

A compact four-component two-dimensional （2-D） finite-difference frequency domain （FDFD） method with the equivalent surface impedance boundary condition is used to analyze the dispersion characteristics of multilayer metal-coated waveguides. According to the equivalent surface impedance boundary condition,the relationship between transverse field components on the boundary can be easily depicted. Once the eigen equation is solved,the propagation constant can be obtained as the eigen value for a given frequency. Results of the proposed method agaree well with those of high frequency structure simulator（HFSS）.

Using a Modified Soil-Plant-Atmosphere Scheme(MSPAS)to Simulate the Interaction between Land Surface Processes and Atmospheric Boundary Layer in Semi-Arid Regions

This paper uses a Modified Soil-Plant-Atmosphere Scheme (MSPAS) to study the interaction between land surface and atmospheric boundary layer processes. The scheme is composed of two main parts: atmospheric boundary layer processes and land surface processes. Compared with SiB and BATS, which are famous for their detailed parameterizations of physical variables, this simplified model is more convenient and saves much more computation time. Though simple, the feasibility of the model is well proved in this paper. The numerical simulation results from MSPAS show good agreement with reality. The scheme is used to obtain reasonable simulations for diurnal variations of heat balance, potential temperature of boundary layer, and wind field, and spatial distributions of temperature, specific humidity, vertical velocity, turbulence kinetic energy, and turbulence exchange coefficient over desert and oasis. In addition, MSPAS is used to simulate the interaction between desert and oasis at night, and again it obtains reasonable results. This indicates that MSPAS can be used to study the interaction between land surface processes and the atmospheric boundary layer over various underlying surfaces and can be extended for regional climate and numerical weather prediction study.

Geometrical design of blade's surface and boundary control of Navier-Stokes equations

In this article a new principle of geometric design for blade's surface of an impeller is provided.This is an optimal control problem for the boundary geometric shape of flow and the control variable is the surface of the blade.We give a minimal functional depending on the geometry of the blade's surface and such that the flow's loss achieves minimum.The existence of the solution of the optimal control problem is proved and the Euler-Lagrange equations for the surface of the blade are derived.In addition,under a new curvilinear coordinate system,the flow domain between the two blades becomes a fixed hexahedron,and the surface as a mapping from a bounded domain in R2 into R3,is explicitly appearing in the objective functional.The Navier-Stokes equations,which include the mapping in their coefficients,can be computed by using operator splitting algorithm.Furthermore,derivatives of the solution of Navier-Stokes equations with respect to the mapping satisfy linearized Navier-Stokes equations which can be solved by using operator splitting algorithms too.Hence,a conjugate gradient method can be used to solve the optimal control problem.

Experimental study on surface arc plasma actuation-based hypersonic boundary layer transition flow control

Effective control of hypersonic transition is essential.In order to avoid affecting the structural proflle of the aircraft,as well as reducing power consumption and electromagnetic interference,a low-frequency surface arc plasma disturbance experiment to promote hypersonic transition was carried out in theΦ0.25 m double-throat Ludwieg tube wind tunnel at Huazhong University of Science and Technology.Contacting printed circuit board sensors and non-contact focused laser differential interferometry testing technology were used in combination.Experimental results showed that the low-frequency surface arc plasma actuation had obvious stimulation effects on the second-mode unstable wave and could promote boundary layer transition by changing the spectral characteristics of the second-mode unstable wave.At the same time,the plasma actuation could promote energy exchange between the second-mode unstable wave and other unstable waves.Finally,the corresponding control mechanism is discussed.

Boundary scheme for lattice Boltzmann modeling of micro-scale gas flow in organic-rich pores considering surface diffusion

We propose a boundary scheme for addressing multi-mechanism flow in a porous medium in slip and early transition flow regimes, which is frequently encountered in shale gas reservoirs. Micro-gaseous flow in organic-rich shale involves a complex flow mechanism. A self-developed boundary scheme that combines the non-equilibrium extrapolation scheme and the combined diffusive reflection and bounce-back scheme(half-way DBB) to embed the Langmuir slip boundary into the single-relaxation-time lattice Boltzmann method(SRT-LBM) enables us to describe this process, namely, the coupling effect of micro-gaseous flow and surface diffusion in organic-rich nanoscale pores. The present LBM model comes with the careful consideration of the local Knudsen number, local pressure gradient, viscosity correction model, and regularization procedure to account for the rarefied gas flows in irregular pores. Its validity and accuracy are verified by several benchmarking cases, and the calculated results by this boundary scheme accord well with our analytical solutions.This boundary scheme shows a higher accuracy than the existing studies. Additionally, a subiteration strategy is presented to tackle the coupled micro-gaseous flow and surface diffusion, which necessitates the iteration process matching of these two mechanisms. The multi-mechanism flow in the self-developed irregular pores is also numerically investigated and analyzed over a wide range of parameters. The results indicate that the present model can effectively capture the coupling effect of micro-gaseous flow and surface diffusion in a tree-like porous medium.

Evaluation of Haney-Type Surface ThermalBoundary Conditions Using a CoupledAtmosphere and Ocean Model

A coupled atmosphere-ocean model developed at the Institute for Space Studies at NASA Goddard Space Flight Center (Russell et al., 1995) was used to verify the validity of Haney-type surface thermal boundary condition, which linearly connects net downward surface heat fluxQ to air/sea temperature difference ΔT by a relaxation coefficientk. The model was initiated from the National Centers for Environmental Prediction (NCEP) atmospheric observations for 1 December 1977, and from the National Ocean Data Center (NODC) global climatological mean December temperature and salinity fields at 1°x 1° resolution. The time step is 7.5 minutes. We integrated the model for 450 days and obtained a complete model-generated global data set of daily mean downward net surface fluxQ, surface air temperatureT A, and sea surface temperatureT O. Then, we calculated the cross-correlation coefficients (CCC) betweenQ and ΔT. The ensemble mean CCC fields show (a) no correlation betweenQ and ΔT in the equatiorial regions, and (b) evident correlation (CCC≥0.7) betweenQ and ΔT in the middle and high latitudes. Additionally, we did the variance analysis and found that whenk=120 W m?2K?1, the two standard deviations, σQ and σκδT , are quite close in the middle and high latitudes. These results agree quite well with a previous research (Chu et al., 1998) on analyzing the NCEP re-analyzed surface data, except that a smaller value ofk (80 W m?2K?1) was found in the previous study. Key words Air-sea coupled system - Ocean surface fluxes - Surface thermal boundary condition

Subdivision Surface-Based Isogeometric Boundary Element Method for Steady Heat Conduction Problems with Variable Coefficient

The present work couples isogeometric analysis(IGA)and boundary element methods(BEM)for three dimensional steady heat conduction problems with variable coefficients.The Computer-Aided Design(CAD)geometries are built by subdivision surfaces,and meantime the basis functions of subdivision surfaces are employed to discretize the boundary integral equations for heat conduction analysis.Moreover,the radial integration method is adopted to transform the additional domain integrals caused by variable coefficients to the boundary integrals.Several numerical examples are provided to demonstrate the correctness and advantages of the proposed algorithm in the integration of CAD and numerical analysis.

A boundary element method for the simulation of non-spherical bubbles and their interactions near a free surface

The basic principle and numerical technique for simulating two three-dimensional bubbles near a free surface are studied in detail by using boundary element method. The singularities of influence coefficient matrix are eliminated using coordinate transformation and so-called 4π rule. The solid angle for the open surface is treated in direct method based on its definition. Several kinds of configurations for the bubbles and free surface have been investigated. The pressure contours during the evolution of bubbles are obtained in our model and can better illuminate the mechanism underlying the motions of bubbles and free surface. The bubble dynamics and their interactions have close relation with the standoff distances, buoyancy parameters and initial sizes of bubbles. Completely different bubble shapes, free surface motions, jetting patterns and pressure distributions under different parameters can be observed in our model, as demon- strated in our calculation results.

The bounding-surfaces record of a barrier spit from Huangqihai Lake,North China:implications for coastal barrier boundary hierarchy

Ground-penetrating radar and trenching studies of a barrier spit on the north shore of Huangqihai Lake were made,that reveal important implications for the coastal washover barrier boundary hierarchy and interpretations of this depositional record.A four-fold hierarchy bounding-surface model,representing different levels of impact and genesis,is defined.Each level of the hierarchy is enclosed by a distinct kind of surface characterized by different ground-penetrating radar reflection features,sedimentary characteristics(color,grain size,sorting,rounding and sedimentary structures) and origin.We suggest that this hierarchical model can be applied to any coastal washover barrier deposits.

MOLECULAR DYNAMICS SIMULATION OF ENTROPY AND SURFACE TENSION FOR GRAIN BOUNDARY OF α-Fe

The grain boundary is an interface and the surface tension is one of its important thermodynamic properties. In this paper, the surface tension of the Σ9 grain boundary for α-Fe at various temperatures and pressures is calculated by means of Computer Molecular Dynamics (CMD). The results agree satisfactorily with the experimental data. It. is shown that the contribution of entropy to surface tension of grain boundary can be ignored.

Boundary layer flow on a moving surface in otherwise quiescent pseudo-plastic non-Newtonian fluids

A theoretical analysis for the boundary layer flow over a continuous moving surface in an otherwise quiescent pseudo-plastic non-Newtonian fluid medium was presented. The types of potential flows necessary for similar solutions to the boundary layer equations were determined and the solutions were numerically presented for different values of power law exponent.

Variable-coordinate forward modeling of irregular surface based on dual-variable grid

The mapping method is a forward-modeling method that transforms the irregular surface to horizontal by mapping the rectangular grid as curved; moreover, the wave field calculations move from the physical domain to the calculation domain. The mapping method deals with the irregular surface and the low-velocity layer underneath it using a fine grid. For the deeper high-velocity layers, the use of a fine grid causes local oversampling. In addition, when the irregular surface is transformed to horizontal, the flattened interface below the surface is transformed to curved, which produces inaccurate modeling results because of the presence of ladder-like burrs in the simulated seismic wave. Thus, we propose the mapping method based on the dual-variable finite-difference staggered grid. The proposed method uses different size grid spacings in different regions and locally variable time steps to match the size variability of grid spacings. Numerical examples suggest that the proposed method requires less memory storage capacity and improves the computational efficiency compared with forward modeling methods based on the conventional grid.

The Mechanism of Drag Reduction around Bodies of Revolution Using Bionic Non-Smooth Surfaces

Bionic non-smooth surfaces （BNSS） can reduce drag. Much attention has been paid to the mechanism of shear stress reduction by riblets. The mechanism of pressure force reduction by bionic non-smooth surfaces on bodies of revolution has not been well investigated. In this work CFD simulation has revealed the mechanism of drag reduction by BNSS, which may work in three ways. First, BNSS on bodies of revolution may lower the surface velocity of the medium, which prevents the sudden speed up of air on the cross section. So the bottom pressure of the model would not be disturbed sharply, resulting in less energy loss and drag reduction. Second, the magnitude of vorticity induced by the bionic model becomes smaller because, due to the sculpturing, the growth of tiny air bubbles is avoided. Thus the large moment of inertia induced by large air bubble is reduced. The reduction of the vorticity could reduce the dissipation of the eddy. So the pressure force could also be reduced. Third, the thickness of the momentum layer on the model becomes less which, according to the relationship between the drag coefficient and the momentum thickness, reduces drag.

Irregular surface seismic forward modeling by a body-fitted rotated–staggered-grid finite-difference method

Finite-difference(FD) methods are widely used in seismic forward modeling owing to their computational efficiency but are not readily applicable to irregular topographies. Thus, several FD methods based on the transformation to curvilinear coordinates using body-fitted grids have been proposed, e.g., stand staggered grid(SSG) with interpolation, nonstaggered grid, rotated staggered grid(RSG), and fully staggered. The FD based on the RSG is somewhat superior to others because it satisfies the spatial distribution of the wave equation without additional memory and computational requirements; furthermore, it is simpler to implement. We use the RSG FD method to transform the firstorder stress–velocity equation in the curvilinear coordinates system and introduce the highprecision adaptive, unilateral mimetic finite-difference(UMFD) method to process the freeboundary conditions of an irregular surface. The numerical results suggest that the precision of the solution is higher than that of the vacuum formalism. When the minimum wavelength is low, UMFD avoids the surface wave dispersion. We compare FD methods based on RSG, SEM, and nonstaggered grid and infer that all simulation results are consistent but the computational efficiency of the RSG FD method is higher than the rest.

Effects of Sea-Surface Waves and Ocean Spray on Air-Sea Momentum Fluxes

The effects of sea-surface waves and ocean spray on the marine atmospheric boundary layer （MABL） at different wind speeds and wave ages were investigated. An MABL model was developed that introduces a wave-induced component and spray force to the total surface stress. The theoretical model solution was determined assuming the eddy viscosity coefficient varied linearly with height above the sea surface. The wave-induced component was evaluated using a directional wave spectrum and growth rate. Spray force was described using interactions between ocean-spray droplets and wind-velocity shear. Wind profiles and sea-surface drag coefficients were calculated for low to high wind speeds for wind-generated sea at different wave ages to examine surface-wave and ocean-spray effects on MABL momentum distribution. The theoretical solutions were compared with model solutions neglecting wave-induced stress and/or spray stress. Surface waves strongly affected near-surface wind profiles and sea-surface drag coefficients at low to moderate wind speeds. Drag coefficients and near-surface wind speeds were lower for young than for old waves. At high wind speeds, ocean-spray droplets produced by wind-tearing breaking-wave crests affected the MABL strongly in comparison with surface waves, implying that wave age affects the MABL only negligibly. Low drag coefficients at high wind caused by ocean-spray production increased turbulent stress in the sea-spray generation layer, accelerating near-sea-surface wind. Comparing the analytical drag coefficient values with laboratory measurements and field observations indicated that surface waves and ocean spray significantly affect the MABL at different wind speeds and wave ages.

Laminar MHD natural convection of nanofluid containing gyrotactic microorganisms over vertical wavy surface saturated non-Darcian porous media

Magnetohydrodynamic （MHD） bioconvection of an incompressible electrically conducting nanofluid near a vertical wavy surface saturated porous medium containing both nanoparticle and gyrotactic microorganisms is investigated. The nanofluid is represented by a model that includes both Brownian motion and thermophoresis effects. A suitable set of non-dimensional variables are used to transform the governing boundary layer equations into a dimensionless form. The resulting nonlinear system is mapped to the vertical flat plate domain, and a non-similar solution is used to the obtained equations. The obtained non-similar system is then solved numerically using the fourth-order Runge-Kutta method. The influence of various physical parameters on the local Nusselt number, the local Sherwood number, the local density number of the motile microorganisms, the dimensionless velocity, the dimensionless temperature, and the rescaled density of motile microorganisms is studied. It is found that the local Nusselt number, the local Sherwood number, and the local density number of the motile microorganisms decrease by increasing either the Grashof number or the magnetic field parameter.

Some Characteristics of the Surface Boundary Layer of a Strong Cold Air Process over Southern China

In southern China, cold air is a common weather process during the winter season; it can cause strong wind, sharp temperature decreases, and even the snow or freezing rain events. However, the features of the atmospheric boundary layer during cold air passage are not clearly understood due to the lack of comprehensive observation data, especially regarding turbulence. In this study, four-layer gradient meteorological observation data and one-layer, 10-Hz ultrasonic anemometer-thermometer monitoring data from the northern side of Poyang Lake were employed to study the main features of the surface boundary layer during a strong cold-air passage over southern China. The results show that, with the passage of a cold air front, the wind speed exhibits low-frequency variations and that the wind systematically descends. During the strong wind period, the wind speed increases with height in the surface layer. Regular gust packets are superimposed on the basic strong wind flow. Before the passage of cold air, the wind gusts exhibit a coherent structure. The wind and turbulent momentum fluxes are small, although the gusty wind momentum flux is slightly larger than the turbulent momentum flux. However, during the invasion of cold air, both the gusty wind and turbulent momentum fluxes increase rapidly with wind speed, and the turbulent momentum flux is larger than the gusty wind momentum flux during the strong wind period. After the cold air invasion, this structure almost disappears.

ASYMPTOTIC BEHAVIOR FOR A CLASS OF ELLIPTICEQUIVALUED SURFACE BOUNDARY VALUE PROBLEM WITH DISCONTINUOUS INTERFACE CONDITIONS

Spontaneous potential well-logging is one of the important techniques in petroleum exploitation. A spontaneous potential satisfies an elliptic equivalued surface boundary value problem with discontinuous interface conditions. In practice, the measuring electrode is so small that we can simplify the corresponding equivalued surface to a point. In this paper, we give a positive answer to this approximation process:when the equivalued surface shrinks to a point, the solution of the original equivalued surface boundary value problem converges to the solution of the corresponding limit boundary value problem.

Large eddy simulation of turbulence in ocean surface boundary layer at Zhangzi Island offshore

This study uses a large eddy simulation （LES） model to investigate the turbulence processes in the ocean surface boundary layer at Zhangzi Island offshore. Field measurements at Zhangzi Island （39°N, 122°E） during July 2009 are used to drive the LES model. The LES results capture a clear diurnal cycle in the oceanic turbulence boundary layer. The process of the heat penetration and heat distribution characteristics are analyzed through the heat flux results from the LES and their differences between two diurnal cycles are discussed as well. Energy balance and other dynamics are investigated which show that the tide-induced shear production is the main source of the turbulence energy that balanced dissipation. Momentum flux near the surface shows better agreement with atmospheric data computed by the eddy correlation method than those computed by bulk formula.