In this paper, the boundary layer equations (abbreviation BLE) for exterior flow around an obstacle are established using semi-geodesic coordinate system (S-coordinate) based on the curved two dimensional surface of t...In this paper, the boundary layer equations (abbreviation BLE) for exterior flow around an obstacle are established using semi-geodesic coordinate system (S-coordinate) based on the curved two dimensional surface of the obstacle. BLE are nonlinear partial differential equations on unknown normal viscous stress tensor and pressure on the obstacle and the existence of solution of BLE is proved. In addition a dimensional split method for dimensional three Navier-Stokes equations is established by applying several 2D-3C partial differential equations on two dimensional manifolds to approach 3D Navier-Stokes equations. The examples for the exterior flow around spheroid and ellipsoid are presents here.展开更多
Third order nonlinear ordinary differential equation, subject to appropriate boundary conditions, arising in fluid mechanics is solved exactly using more suggestive schemes- Dirichlet series and method of stretching v...Third order nonlinear ordinary differential equation, subject to appropriate boundary conditions, arising in fluid mechanics is solved exactly using more suggestive schemes- Dirichlet series and method of stretching variables. These methods have advantages over pure numerical methods in obtaining derived quantities accurately for various values of the parameters involved at a stretch and are valid in a much larger domain compared with classical numerical schemes.展开更多
This study evaluates the effectiveness of a new technique that transforms doma in integrals into boundary integrals that is applicable to the boundary element method.Si mulations were conducted in which two-dimensiona...This study evaluates the effectiveness of a new technique that transforms doma in integrals into boundary integrals that is applicable to the boundary element method.Si mulations were conducted in which two-dimensional surfaces were approximated by inter polation using radial basis functions with full and compact supports.Examples involving Poisson’s equation are presented using the boundary element method and the proposed te chnique with compact radial basis functions.The advantages and the disadvantages are e xamined through simulations.The effects of internal poles,the boundary mesh refinemen t and the value for the support of the radial basis functions on performance are assessed.展开更多
The meshless method is a new numerical technique presented in recent years.It uses the moving least square(MLS)approximation as a shape function.The smoothness of the MLS approximation is determined by that of the bas...The meshless method is a new numerical technique presented in recent years.It uses the moving least square(MLS)approximation as a shape function.The smoothness of the MLS approximation is determined by that of the basic function and of the weight function,and is mainly determined by that of the weight function.Therefore,the weight function greatly affects the accuracy of results obtained.Different kinds of weight functions,such as the spline function, the Gauss function and so on,are proposed recently by many researchers.In the present work,the features of various weight functions are illustrated through solving elasto-static problems using the local boundary integral equation method.The effect of various weight functions on the accuracy, convergence and stability of results obtained is also discussed.Examples show that the weight function proposed by Zhou Weiyuan and Gauss and the quartic spline weight function are better than the others if parameters c and α in Gauss and exponential weight functions are in the range of reasonable values,respectively,and the higher the smoothness of the weight function,the better the features of the solutions.展开更多
In the adjoint-state method, the forward-propagated source wavefield and the backward-propagated receiver wavefield must be available simultaneously either for seismic imaging in migration or for gradient calculation ...In the adjoint-state method, the forward-propagated source wavefield and the backward-propagated receiver wavefield must be available simultaneously either for seismic imaging in migration or for gradient calculation in inversion. A feasible way to avoid the excessive storage demand is to reconstruct the source wavefield backward in time by storing the entire history of the wavefield in perfectly matched layers. In this paper, we make full use of the elementwise global property of the Laplace operator of the spectral element method (SEM) and propose an efficient source wavefield reconstruction method at the cost of storing the wavefield history only at single boundary layer nodes. Numerical experiments indicate that the accuracy of the proposed method is identical to that of the conventional method and is independent of the order of the Lagrange polynomials, the element type, and the temporal discretization method. In contrast, the memory-saving ratios of the conventional method versus our method is at least N when using either quadrilateral or hexahedron elements, respectively, where N is the order of the Lagrange polynomials used in the SEM. A higher memorysaving ratio is achieved with triangular elements versus quadrilaterals. The new method is applied to reverse time migration by considering the Marmousi model as a benchmark. Numerical results demonstrate that the method is able to provide the same result as the conventional method but with about 1/25 times lower storage demand. With the proposed wavefield reconstruction method, the storage demand is dramatically reduced;therefore, in-core memory storage is feasible even for large-scale three-dimensional adjoint inversion problems.展开更多
The effects of streaks on boundary layer transition depend on the initial amplitude of T-S waves introduced to excite the transition. This problem was studied in a flat-plate boundary layer in water tunnel by using hy...The effects of streaks on boundary layer transition depend on the initial amplitude of T-S waves introduced to excite the transition. This problem was studied in a flat-plate boundary layer in water tunnel by using hydrogen bubble method. Three T-S wave initial amplitudes were tested. The results show that both narrow and wide-spacing streaks depress the transition excited by T-S waves with lower initial amplitude. However, when transition is excited by T-S waves of higher initial amplitude, the narrow-spacing streaks depress the transition, while the wide-spacing streaks promote the transition. Futrther the underlying mechanisms were also analyzed.展开更多
This paper presents a new method to estimate the height of the atmospheric boundary layer(ABL) by using COSMIC radio occultation bending angle(BA) data. Using the numerical differentiation method combined with the reg...This paper presents a new method to estimate the height of the atmospheric boundary layer(ABL) by using COSMIC radio occultation bending angle(BA) data. Using the numerical differentiation method combined with the regularization technique, the first derivative of BA profiles is retrieved, and the height at which the first derivative of BA has the global minimum is defined to be the ABL height. To reflect the reliability of estimated ABL heights, the sharpness parameter is introduced, according to the relative minimum of the BA derivative. Then, it is applied to four months of COSMIC BA data(January, April, July, and October in 2008), and the ABL heights estimated are compared with two kinds of ABL heights from COSMIC products and with the heights determined by the finite difference method upon the refractivity data. For sharp ABL tops(large sharpness parameters), there is little difference between the ABL heights determined by different methods, i.e.,the uncertainties are small; whereas, for non-sharp ABL tops(small sharpness parameters), big differences exist in the ABL heights obtained by different methods, which means large uncertainties for different methods. In addition, the new method can detect thin ABLs and provide a reference ABL height in the cases eliminated by other methods. Thus, the application of the numerical differentiation method combined with the regularization technique to COSMIC BA data is an appropriate choice and has further application value.展开更多
The density distribution of a supersonic turbulent boundary layer is measured with the nanoparticle-based planar laser scattering technique, and the temporal evolution of its optical path difference （OPD） in a short...The density distribution of a supersonic turbulent boundary layer is measured with the nanoparticle-based planar laser scattering technique, and the temporal evolution of its optical path difference （OPD） in a short time interval is characterized by proper orthogonal decomposition （POD）. Based on the advantage of POD in capturing the energy of a signal, a temporal evolution model is suggested for the POD coefficients of the OPD. In this model, the first few coefficients vary linearly with time, and the others are modeled by Gaussian statistics. As an application, this method is used to compute the shortexposure optical transfer function.展开更多
A boundary integral method with radial basis function approximation is proposed for numerically solving an important class of boundary value problems governed by a system of thermoelastostatic equations with variable ...A boundary integral method with radial basis function approximation is proposed for numerically solving an important class of boundary value problems governed by a system of thermoelastostatic equations with variable coe?cients. The equations describe the thermoelastic behaviors of nonhomogeneous anisotropic materials with properties that vary smoothly from point to point in space. No restriction is imposed on the spatial variations of the thermoelastic coe?cients as long as all the requirements of the laws of physics are satis?ed. To check the validity and accuracy of the proposed numerical method, some speci?c test problems with known solutions are solved.展开更多
In this article, the model of a non-Newtonian fluid (Thixotropic) flow past a vertical surface in the presence of exponential space and temperature dependent heat source in a thermally stratified medium is studied. It...In this article, the model of a non-Newtonian fluid (Thixotropic) flow past a vertical surface in the presence of exponential space and temperature dependent heat source in a thermally stratified medium is studied. It is assumed that free convection is induced by buoyancy and exponentially decaying internal heat source across the space. The dynamic viscosity is taken to be constant and thermal conductivity of this particular fluid model is assumed to vary linearly with temperature. Thermal stratification has been properly incorporated into the governing equation so that its effect can be revealed and properly reported. The governing partial differential equations describing the model are transformed and parameterized to a system of non-linear ordinary differential equation using similarity transformations. Approximate analytic solutions were obtained by adopting Optimal Homotopy Analysis Method (OHAM). The results show that for both cases of non-Newtonian parameters (Thixotropic) (K1=K2=0?& K1=K2=1.0), increasing stratification parameters, relate to decreasing in the heat energy entering into the fluid region and thus reducing the temperature of the Thixotropic fluid as it flows.展开更多
The present study reveals the effect of nonlinear thermal radiation and magnetic field on a boundary layer flow of a viscous fluid over a nonlinear stretching sheet with suction or an injection. Using suitable similar...The present study reveals the effect of nonlinear thermal radiation and magnetic field on a boundary layer flow of a viscous fluid over a nonlinear stretching sheet with suction or an injection. Using suitable similarity transformations, governing partial differential equations were reduced to higher order ordinary differential equations and further these are solved numerically using of Keller-Box method. Effect of flow controlling parameter on velocity, temperature and nanoparticle fluid concentration, local skin friction coefficient, local Nusselt number and local Sherwood numbers are discussed. It is found that the dimensionless velocity decreases and temperature, concentration are increased with the increasing of magnetic parameter. The temperature profile is an increasing function of thermal radiation when it is increasing.展开更多
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...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.展开更多
It is important to improve the speed of a ship,the friction resistance can be reduced by injection air at the bottom of a ship when the ship is running on the water.As the first part of the studying project,here numer...It is important to improve the speed of a ship,the friction resistance can be reduced by injection air at the bottom of a ship when the ship is running on the water.As the first part of the studying project,here numerical simulation study method,boundary condition and governing equations are presented.It is easy to study complicated problems from simple conditions,so the program concerning boundary layer condition is compiled to solve the problem.Here the spectral method is introduced,and the results are tested by Dorod’s results.展开更多
Non-singular fictitious boundary integral equations for orthotropic elastic plane problems were deduced according to boundary conditions by the techniques of singular-points-outside-domain. Then the unknown fictitious...Non-singular fictitious boundary integral equations for orthotropic elastic plane problems were deduced according to boundary conditions by the techniques of singular-points-outside-domain. Then the unknown fictitious load functions along the fictitious boundary were expressed in terms of basic spline functions, and the boundary-segment-least-squares method was proposed to eliminate the boundary residues obtained. By the above steps, numerical solutions to the integral equations can be achieved. Numerical examples are given to show the accuracy and efficiency of the proposed method.展开更多
A structure function approach is applied to estimate the turbulent kinetic energy （TKE） dissipation rate in the bottom boundary layer of the Pearl River Estuary （PRE）. Simultaneous measurements with an acoustic Do...A structure function approach is applied to estimate the turbulent kinetic energy （TKE） dissipation rate in the bottom boundary layer of the Pearl River Estuary （PRE）. Simultaneous measurements with an acoustic Doppler velocimeter （ADV） supplied independent data for the verification of the structure function method. The results show that, 1） the structure function approach is reliable and successfully applied method to estimate the TKE dissipation rate. The observed dissipation rates range between 8.3 ×10^-4 W/kg and 4.9× 10^-6 W/kg in YM01 and between 3.4×10^-4 W/kg and 4.8×10^-7 W/kg in YM03, respectively, while exhibiting a strong quarter-diurnal variation. 2） The balance between the shear production and viscous dissipation is better achieved in the straight river. This first-order balance is significantly broken in the estuary by non-shear production/dissipation due to wave-induced fluctuations.展开更多
In this paper,on the basis of Ref.[1],the author studies the boundary value problems of the second-order differential equations,the highest order derivatives of which contain the small parameters.The numerical example...In this paper,on the basis of Ref.[1],the author studies the boundary value problems of the second-order differential equations,the highest order derivatives of which contain the small parameters.The numerical examples show that the calculating process of this method is quite simple and its accuracy is even higher than that of the multiple scales method.展开更多
Electromagnetic scattering from a rough surface of layered medium is investigated, and the formulae of the scattering coefficients for different polarizations are derived using the small perturbation method. A rough s...Electromagnetic scattering from a rough surface of layered medium is investigated, and the formulae of the scattering coefficients for different polarizations are derived using the small perturbation method. A rough surface with exponential correlation function is presented for describing a rough soil surface of layered medium, the formula of its scattering coefficient is derived by considering the spectrum of the rough surface with exponential correlation function; the curves of the bistatic scattering coefficient of HH polarization with variation of the scattering angle are obtained by numerical calculation. The influence of the permittivity of layered medium, the mean layer thickness of intermediate medium, the roughness surface parameters and the frequency of the incident wave on the blstatic scattering coefficient is discussed. Numerical results show that the influence of the permittivity of layered medium, the mean layer thickness of intermediate medium, the rms and the correlation length of the rough surface, and the frequency of the incident wave on the bistatic scattering coefficient is very complex.展开更多
An iteration method similar to the thin-wing-expansion method for the compressible flow has been proposed to solve the boundary layer flow past a flat plate. Using such an iteration, the first step of which is Oseen’...An iteration method similar to the thin-wing-expansion method for the compressible flow has been proposed to solve the boundary layer flow past a flat plate. Using such an iteration, the first step of which is Oseen’s approximation, the boundary layer past a flat plate is studied. As proceeding from the first approximation to the second and third approximations, it is realized that our solution approaches to a well known Howarth’s bench mark one gradually. Hence, it is concluded that the usefulness of the present method has been confirmed.展开更多
The inverse problems of wave equation to recover unknown space-time dependent functions of wave speed and wave source are solved in this paper, without needing of initial conditions and no internal measurement of data...The inverse problems of wave equation to recover unknown space-time dependent functions of wave speed and wave source are solved in this paper, without needing of initial conditions and no internal measurement of data being required. After a homogenization technique, a sequence of spatial boundary functions at least the fourth-order polynomials are derived, which satisfy the homogeneous boundary conditions. The boundary functions and the zero element constitute a linear space, and then a new boundary functional is proved in the linear space, of which the energy is preserved for each dynamic energetic boundary function. The linear systems and iterative algorithms used to recover unknown wave speed and wave source functions with the dynamic energetic boundary functions as bases are developed, which converge fast at each time step. The input data are parsimonious, merely the measured boundary strains and the boundary values and slopes of unknown functions to be recovered. The accuracy and robustness of present methods are confirmed by comparing exact solutions with estimated results under large noises up to 20%.展开更多
In the design of the fatigue strength of dynamically loaded bearing in the equipmentssuch as internal combustion engines and roimg mun, the solution to the stress distribution on thebushing alloy layer is an important...In the design of the fatigue strength of dynamically loaded bearing in the equipmentssuch as internal combustion engines and roimg mun, the solution to the stress distribution on thebushing alloy layer is an important and difficult problem. In this paper, a new method has beenproposed by coupling BEM with etheticity method, The algorithm and its implementation were deseribed in details The calculation results verify that this up-dated method can provide us a moresimple and effective tool for solvingthe fatigue stress of the bushing alloy with tangible benefit oftime-saving and high computation accuraey. It may open a new vista in bearing fatigue strength design.展开更多
文摘In this paper, the boundary layer equations (abbreviation BLE) for exterior flow around an obstacle are established using semi-geodesic coordinate system (S-coordinate) based on the curved two dimensional surface of the obstacle. BLE are nonlinear partial differential equations on unknown normal viscous stress tensor and pressure on the obstacle and the existence of solution of BLE is proved. In addition a dimensional split method for dimensional three Navier-Stokes equations is established by applying several 2D-3C partial differential equations on two dimensional manifolds to approach 3D Navier-Stokes equations. The examples for the exterior flow around spheroid and ellipsoid are presents here.
文摘Third order nonlinear ordinary differential equation, subject to appropriate boundary conditions, arising in fluid mechanics is solved exactly using more suggestive schemes- Dirichlet series and method of stretching variables. These methods have advantages over pure numerical methods in obtaining derived quantities accurately for various values of the parameters involved at a stretch and are valid in a much larger domain compared with classical numerical schemes.
文摘This study evaluates the effectiveness of a new technique that transforms doma in integrals into boundary integrals that is applicable to the boundary element method.Si mulations were conducted in which two-dimensional surfaces were approximated by inter polation using radial basis functions with full and compact supports.Examples involving Poisson’s equation are presented using the boundary element method and the proposed te chnique with compact radial basis functions.The advantages and the disadvantages are e xamined through simulations.The effects of internal poles,the boundary mesh refinemen t and the value for the support of the radial basis functions on performance are assessed.
文摘The meshless method is a new numerical technique presented in recent years.It uses the moving least square(MLS)approximation as a shape function.The smoothness of the MLS approximation is determined by that of the basic function and of the weight function,and is mainly determined by that of the weight function.Therefore,the weight function greatly affects the accuracy of results obtained.Different kinds of weight functions,such as the spline function, the Gauss function and so on,are proposed recently by many researchers.In the present work,the features of various weight functions are illustrated through solving elasto-static problems using the local boundary integral equation method.The effect of various weight functions on the accuracy, convergence and stability of results obtained is also discussed.Examples show that the weight function proposed by Zhou Weiyuan and Gauss and the quartic spline weight function are better than the others if parameters c and α in Gauss and exponential weight functions are in the range of reasonable values,respectively,and the higher the smoothness of the weight function,the better the features of the solutions.
基金financial support for this work contributed by the National Key Research and Development Program of China (grant numbers 2016YFC0600101 and 2016YFC 0600201)the National Natural Science Foundation of China (grant numbers 41874065, 41604076, 41674102, 41674095, 41522401, 41574082, and 41774097)
文摘In the adjoint-state method, the forward-propagated source wavefield and the backward-propagated receiver wavefield must be available simultaneously either for seismic imaging in migration or for gradient calculation in inversion. A feasible way to avoid the excessive storage demand is to reconstruct the source wavefield backward in time by storing the entire history of the wavefield in perfectly matched layers. In this paper, we make full use of the elementwise global property of the Laplace operator of the spectral element method (SEM) and propose an efficient source wavefield reconstruction method at the cost of storing the wavefield history only at single boundary layer nodes. Numerical experiments indicate that the accuracy of the proposed method is identical to that of the conventional method and is independent of the order of the Lagrange polynomials, the element type, and the temporal discretization method. In contrast, the memory-saving ratios of the conventional method versus our method is at least N when using either quadrilateral or hexahedron elements, respectively, where N is the order of the Lagrange polynomials used in the SEM. A higher memorysaving ratio is achieved with triangular elements versus quadrilaterals. The new method is applied to reverse time migration by considering the Marmousi model as a benchmark. Numerical results demonstrate that the method is able to provide the same result as the conventional method but with about 1/25 times lower storage demand. With the proposed wavefield reconstruction method, the storage demand is dramatically reduced;therefore, in-core memory storage is feasible even for large-scale three-dimensional adjoint inversion problems.
基金supported by the National Natural Science Foundation of China (11072017)
文摘The effects of streaks on boundary layer transition depend on the initial amplitude of T-S waves introduced to excite the transition. This problem was studied in a flat-plate boundary layer in water tunnel by using hydrogen bubble method. Three T-S wave initial amplitudes were tested. The results show that both narrow and wide-spacing streaks depress the transition excited by T-S waves with lower initial amplitude. However, when transition is excited by T-S waves of higher initial amplitude, the narrow-spacing streaks depress the transition, while the wide-spacing streaks promote the transition. Futrther the underlying mechanisms were also analyzed.
基金supported by the National Natural Science Foundation of China (Grant No. 41475021)
文摘This paper presents a new method to estimate the height of the atmospheric boundary layer(ABL) by using COSMIC radio occultation bending angle(BA) data. Using the numerical differentiation method combined with the regularization technique, the first derivative of BA profiles is retrieved, and the height at which the first derivative of BA has the global minimum is defined to be the ABL height. To reflect the reliability of estimated ABL heights, the sharpness parameter is introduced, according to the relative minimum of the BA derivative. Then, it is applied to four months of COSMIC BA data(January, April, July, and October in 2008), and the ABL heights estimated are compared with two kinds of ABL heights from COSMIC products and with the heights determined by the finite difference method upon the refractivity data. For sharp ABL tops(large sharpness parameters), there is little difference between the ABL heights determined by different methods, i.e.,the uncertainties are small; whereas, for non-sharp ABL tops(small sharpness parameters), big differences exist in the ABL heights obtained by different methods, which means large uncertainties for different methods. In addition, the new method can detect thin ABLs and provide a reference ABL height in the cases eliminated by other methods. Thus, the application of the numerical differentiation method combined with the regularization technique to COSMIC BA data is an appropriate choice and has further application value.
基金Project supported by the Innovation Research Foundations for Postgraduates of National University of Defense Technology and Hunan Province, Chinathe National Natural Science Foundation of China (Grant No. 61008037)
文摘The density distribution of a supersonic turbulent boundary layer is measured with the nanoparticle-based planar laser scattering technique, and the temporal evolution of its optical path difference （OPD） in a short time interval is characterized by proper orthogonal decomposition （POD）. Based on the advantage of POD in capturing the energy of a signal, a temporal evolution model is suggested for the POD coefficients of the OPD. In this model, the first few coefficients vary linearly with time, and the others are modeled by Gaussian statistics. As an application, this method is used to compute the shortexposure optical transfer function.
文摘A boundary integral method with radial basis function approximation is proposed for numerically solving an important class of boundary value problems governed by a system of thermoelastostatic equations with variable coe?cients. The equations describe the thermoelastic behaviors of nonhomogeneous anisotropic materials with properties that vary smoothly from point to point in space. No restriction is imposed on the spatial variations of the thermoelastic coe?cients as long as all the requirements of the laws of physics are satis?ed. To check the validity and accuracy of the proposed numerical method, some speci?c test problems with known solutions are solved.
文摘In this article, the model of a non-Newtonian fluid (Thixotropic) flow past a vertical surface in the presence of exponential space and temperature dependent heat source in a thermally stratified medium is studied. It is assumed that free convection is induced by buoyancy and exponentially decaying internal heat source across the space. The dynamic viscosity is taken to be constant and thermal conductivity of this particular fluid model is assumed to vary linearly with temperature. Thermal stratification has been properly incorporated into the governing equation so that its effect can be revealed and properly reported. The governing partial differential equations describing the model are transformed and parameterized to a system of non-linear ordinary differential equation using similarity transformations. Approximate analytic solutions were obtained by adopting Optimal Homotopy Analysis Method (OHAM). The results show that for both cases of non-Newtonian parameters (Thixotropic) (K1=K2=0?& K1=K2=1.0), increasing stratification parameters, relate to decreasing in the heat energy entering into the fluid region and thus reducing the temperature of the Thixotropic fluid as it flows.
文摘The present study reveals the effect of nonlinear thermal radiation and magnetic field on a boundary layer flow of a viscous fluid over a nonlinear stretching sheet with suction or an injection. Using suitable similarity transformations, governing partial differential equations were reduced to higher order ordinary differential equations and further these are solved numerically using of Keller-Box method. Effect of flow controlling parameter on velocity, temperature and nanoparticle fluid concentration, local skin friction coefficient, local Nusselt number and local Sherwood numbers are discussed. It is found that the dimensionless velocity decreases and temperature, concentration are increased with the increasing of magnetic parameter. The temperature profile is an increasing function of thermal radiation when it is increasing.
基金Project financially supported by scientific research foundation coferring to Ph.D.
文摘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.
文摘It is important to improve the speed of a ship,the friction resistance can be reduced by injection air at the bottom of a ship when the ship is running on the water.As the first part of the studying project,here numerical simulation study method,boundary condition and governing equations are presented.It is easy to study complicated problems from simple conditions,so the program concerning boundary layer condition is compiled to solve the problem.Here the spectral method is introduced,and the results are tested by Dorod’s results.
文摘Non-singular fictitious boundary integral equations for orthotropic elastic plane problems were deduced according to boundary conditions by the techniques of singular-points-outside-domain. Then the unknown fictitious load functions along the fictitious boundary were expressed in terms of basic spline functions, and the boundary-segment-least-squares method was proposed to eliminate the boundary residues obtained. By the above steps, numerical solutions to the integral equations can be achieved. Numerical examples are given to show the accuracy and efficiency of the proposed method.
基金supported by the National Natural Science Foundation of China (Grant No. 41006050)the China Postdoctoral Science Foundation (Grant No. 20090460799)the Fundamental Research Funds for the Central Universities (Grant No. 11lgpy59)
文摘A structure function approach is applied to estimate the turbulent kinetic energy （TKE） dissipation rate in the bottom boundary layer of the Pearl River Estuary （PRE）. Simultaneous measurements with an acoustic Doppler velocimeter （ADV） supplied independent data for the verification of the structure function method. The results show that, 1） the structure function approach is reliable and successfully applied method to estimate the TKE dissipation rate. The observed dissipation rates range between 8.3 ×10^-4 W/kg and 4.9× 10^-6 W/kg in YM01 and between 3.4×10^-4 W/kg and 4.8×10^-7 W/kg in YM03, respectively, while exhibiting a strong quarter-diurnal variation. 2） The balance between the shear production and viscous dissipation is better achieved in the straight river. This first-order balance is significantly broken in the estuary by non-shear production/dissipation due to wave-induced fluctuations.
文摘In this paper,on the basis of Ref.[1],the author studies the boundary value problems of the second-order differential equations,the highest order derivatives of which contain the small parameters.The numerical examples show that the calculating process of this method is quite simple and its accuracy is even higher than that of the multiple scales method.
基金supported by the National Natural Science Foundation of China (Grant No 60571058)the Specialized Research Fund for the Doctoral Program of Higher Education Institutions of China (Grant No 20070701010)
文摘Electromagnetic scattering from a rough surface of layered medium is investigated, and the formulae of the scattering coefficients for different polarizations are derived using the small perturbation method. A rough surface with exponential correlation function is presented for describing a rough soil surface of layered medium, the formula of its scattering coefficient is derived by considering the spectrum of the rough surface with exponential correlation function; the curves of the bistatic scattering coefficient of HH polarization with variation of the scattering angle are obtained by numerical calculation. The influence of the permittivity of layered medium, the mean layer thickness of intermediate medium, the roughness surface parameters and the frequency of the incident wave on the blstatic scattering coefficient is discussed. Numerical results show that the influence of the permittivity of layered medium, the mean layer thickness of intermediate medium, the rms and the correlation length of the rough surface, and the frequency of the incident wave on the bistatic scattering coefficient is very complex.
文摘An iteration method similar to the thin-wing-expansion method for the compressible flow has been proposed to solve the boundary layer flow past a flat plate. Using such an iteration, the first step of which is Oseen’s approximation, the boundary layer past a flat plate is studied. As proceeding from the first approximation to the second and third approximations, it is realized that our solution approaches to a well known Howarth’s bench mark one gradually. Hence, it is concluded that the usefulness of the present method has been confirmed.
文摘The inverse problems of wave equation to recover unknown space-time dependent functions of wave speed and wave source are solved in this paper, without needing of initial conditions and no internal measurement of data being required. After a homogenization technique, a sequence of spatial boundary functions at least the fourth-order polynomials are derived, which satisfy the homogeneous boundary conditions. The boundary functions and the zero element constitute a linear space, and then a new boundary functional is proved in the linear space, of which the energy is preserved for each dynamic energetic boundary function. The linear systems and iterative algorithms used to recover unknown wave speed and wave source functions with the dynamic energetic boundary functions as bases are developed, which converge fast at each time step. The input data are parsimonious, merely the measured boundary strains and the boundary values and slopes of unknown functions to be recovered. The accuracy and robustness of present methods are confirmed by comparing exact solutions with estimated results under large noises up to 20%.
文摘In the design of the fatigue strength of dynamically loaded bearing in the equipmentssuch as internal combustion engines and roimg mun, the solution to the stress distribution on thebushing alloy layer is an important and difficult problem. In this paper, a new method has beenproposed by coupling BEM with etheticity method, The algorithm and its implementation were deseribed in details The calculation results verify that this up-dated method can provide us a moresimple and effective tool for solvingthe fatigue stress of the bushing alloy with tangible benefit oftime-saving and high computation accuraey. It may open a new vista in bearing fatigue strength design.
