首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 390 毫秒
1.
In our previous works, based on numerical models, it was shown that under certain conditions a hot material can rise in portions in the tails of thermal mantle plumes. The spectrum of these pulsations can correspond to the observed spectra of catastrophic hotspot eruptions. Since most of the existing numerical models of thermal convection for the mantle of the present Earth do not reveal these pulsations, in this work, we analyze the physical cause and initiation conditions of pulsations of thermal plumes. The results of a numerical solution of the thermal convection equations for a material with varying parameters in the extended Boussinesq approximation are presented. It is shown how the structure of the convection is transformed with the increase of convection intensity. At the Rayleigh numbers Ra > 106, convection becomes unsteady, and the configuration of the ascending and descending flows changes. The new flow emerging at the mantle bottom acquires a mushroom shape with a head and a tail. After the rise of the plume’s head to the surface, the tail remains in the mantle in the form of a quasi-stationary hot steam. It turns out that at Ra ~ 5 × 107, the thermal mantle plume becomes pulsating and its tail is in fact a heated channel through which the hot material rises in successive portions. At the Rayleigh numbers Ra > 5 × 108, the tail of the thermal plume breaks and the plume becomes a regular conveyor of separate ascending portions of the hot material, which are referred to as thermals. Thus, thermal convection with pulsating plumes takes place at the transitional stage from the regime of quasi-stationary plumes to the regime of thermals.  相似文献   

2.
The similarity solutions of the equations of kinematic flow on a cone-shaped surface, provide a means for checking the numerical models used for converging and diverging flows. These models can be further improved if the initial excess rainfall curve can be described as a power of time. Illustration of this improvement and assessment of numerical models are given for a constant excess rainfall and for excess rainfall proportional to the inverse of the square root of time.  相似文献   

3.
In this paper we consider one-dimensional capillary redistribution of two immiscible and incompressible fluids in a porous medium with a single discontinuity. We study a special time-dependent solution, a similarity solution, which is found when the initial saturation is discontinuous at the same point as the permeability and porosity, and is constant elsewhere. The similarity solution can be used to validate numerical algorithms describing two-phase flow in porous media with discontinuous heterogeneities. We discuss the construction of the similarity solution, in which we pay special attention to the interface conditions at the discontinuity, both for media with positive and zero entry pressure. Moreover, we discuss some qualitative properties of the solution, and outline a numerical procedure to determine its graph. Examples are given for the Brooks-Corey and Van Genuchten model. We also consider similarity solutions for unsaturated water flow, which is a limit case of two-phase flow for negligible nonwetting phase viscosity.  相似文献   

4.
Presently, the study of the mantle flow structure is mainly based on numerical modeling. The most important stage of the development of a computer program is its testing. For this purpose, results of various test models of convection flows with a given set of parameters are compared. The solution of the Stokes equation, involving the derivative of viscous stresses, is most difficult. Exact analytical solutions of the Stokes equation are obtained in this work for various cases of special loads. These solutions can be used as benchmarks for testing programs of numerical calculation of viscous flows in both geophysics and engineering. The advantage of this testing technique is the exceptional simplicity of the solution form, the admissibility of any spatial viscosity variations, and the fact that solutions can be compared not for a narrow set of the solution parameters but for any distributions of velocities, viscous stresses, and pressures at all points of the space.  相似文献   

5.
A strong Coriolis force in a rapidly rotating planet and star not only enforces two-dimensionality of fluid motion driven by thermal instabilities but also generates strong differential rotation even in the vicinity of the onset of instabilities. We derive an asymptotic solution describing convection-driven differential rotation in rotating, self-gravitating Boussinesq fluid spheres with the no-slip boundary condition, taking into account full spherical curvature and being valid for asymptotically small Ekman numbers. For the purpose of validating the asymptotic solution, the corresponding numerical analysis valid for large or small Ekman numbers is also carried out, showing a satisfactory agreement between the asymptotic and numerical solutions.  相似文献   

6.
In this study, a new method for conversion of solid finite element solution to beam finite element solution is developed based on the meta-modeling theory which constructs a model consistent with continuum mechanics. The proposed method is rigorous and efficient compared to a typical conversion method which merely computes surface integration of solid element nodal stresses to obtain cross-sectional forces. The meta-modeling theory ensures the rigorousness of proposed method by defining a proper distance between beam element and solid element solutions in a function space of continuum mechanics. Results of numerical verification test that is conducted with a simple cantilever beam are used to find the proper distance function for this conversion. Time history analysis of the main tunnel structure of a real ramp tunnel is considered as a numerical example for the proposed conversion method. It is shown that cross-sectional forces are readily computed for solid element solution of the main tunnel structure when it is converted to a beam element solution using the proposed method. Further, envelopes of resultant forces which are of primary importance for the purpose of design, are developed for a given ground motion at the end.  相似文献   

7.
任意圆弧形凸起地形对平面SH波的散射   总被引：22，自引：10，他引：22  
本文采用波函数展开方法提出弹性半空间表面一任意圆弧形凸起边界对平面ＳＨ波二维散射的封闭级数解答，利用引入的辅助函数和推广的外域型Ｇｒａｆ加法公式将解答归结为一代无穷代数方程组的求解，解答的数值结果可由无穷方程的截断计算得出，文中从级数项数增加时计算结果的收敛以及边界条件的收敛满足两方面检验了截断计算的精度，并指出了位移解答的收敛速度与主要参数之间的关系，同时，对凸起角点的应力奇异性及其对解答精度的  相似文献   

8.
《Advances in water resources》2003,26(11):1189-1198
A two-dimensional finite element based overland flow model was developed and used to study the accuracy and stability of three numerical schemes and watershed parameter aggregation error. The conventional consistent finite element scheme results in oscillations for certain time step ranges. The lumped and the upwind finite element schemes are tested as alternatives to the consistent scheme. The upwind scheme did not improve on the stability or the accuracy of the solution, while the lumped scheme provided stable and accurate solutions for time steps twice the size of time steps needed for the consistent scheme. A new accuracy based dynamic time step estimate for the two-dimensional overland flow kinematic wave solution is developed for the lumped scheme. The newly developed dynamic time step estimates are functions of the mesh size, and time of concentration of the watershed hydrograph. Due to lack of analytical solutions, the time step was developed by comparing numerical solutions of various levels of discretization to a reference solution using a very fine mesh and a very small time step. The time step criteria were tested on a different set of problems and proved to be adequate for accurate and stable solutions. A sensitivity analysis for the watershed slope, Manning’s roughness coefficient and excess rainfall rate was conducted in order to test the effect of parameter aggregation on the stability and accuracy of the solution. The results of this analysis show that aggregation of the slope data resulted in the highest error. The roughness coefficient had a smaller effect on the solution while the rainfall intensity did not show any significant effect on the flow rate solution for the range of rainfall intensity used. This work pioneers the challenge of providing guidelines for accurate and stable numerical solutions of the two-dimensional kinematic wave equations for overland flow.  相似文献   

9.
An ensemble of cold thermals in a convective layer of a fresh or salt water body is considered. The probability distribution density of the vertical velocities of thermals is assumed to satisfy the K-form of Fokker?Planck equation. The K-form of the Fokker?Planck equation is shown to have a steady-state solution in the form of a generalized Maxwell distribution. Natural and numerical experiments were carried out to demonstrate the existence of such distribution. Lyapunov functional in the form of Kullback relationship was used to establish the nonlinear instability of this distribution with respect to small perturbations in the Hilbert space norm.  相似文献   

10.
波动问题有限元离散后会引起数值误差, 数值频散的本质就是数值误差传播引起的非物理解. 数值频散不仅没有实际意义, 而且还会影响对真实波动现象的认识. 为厘清有限元三角网格中波动数值频散的影响因素, 本文推导了集中质量矩阵和一致质量矩阵的频散函数, 同时给出了组合质量矩阵的频散函数, 并对不同质量矩阵的数值频散进行了对比研究. 理论分析和数值计算结果表明: 有限元三角网格中波动的数值频散受网格布局、 波传播方向、 单元网格纵横比以及质量矩阵的影响; 一致质量矩阵的数值频散比集中质量矩阵更易受到波传播方向的影响; 不合理的三角网格单元会对数值相速度(数值频散)产生不良影响; 正三角网格中波动的数值频散几乎不受波传播方向的影响; 一致质量矩阵与集中质量矩阵的线性组合能够有效地压制数值频散.   相似文献   

11.
This paper presents a method for constructing polynomial-based approximate solutions to the Boussinesq equation with cylindrical symmetry. This equation models water injection at a single well in an unconfined aquifer; as a sample problem we examine recharge of an initially empty aquifer. For certain injection regimes it is possible to introduce similarity variables, reducing the original problem to a boundary-value problem for an ordinary differential equation. The approximate solutions introduced here incorporate both a singular part to model the behavior near the well and a polynomial part to model the behavior in the far field. Although the nonlinearity of the problem prevents decoupling of the singular and polynomial parts, the paper presents an approach for calculating the solution based on its spatial moments. This approach yields closed-form expressions for the position of the wetting front and for the form of the phreatic surface. Comparison with a highly accurate numerical solution verifies the accuracy of the newly derived approximate solutions.  相似文献   

12.
本文就非弹性地震反应谱的计算问题（运动方程的形式，谱参量的选取，方程解算中恢复力模型上拐折点的处理和时间积分方法），给出了几点理论考虑，介绍了一个改进的算法。这一算法具有精度好和效率高的优点，这对做好结构抗震设计和其它地震工程问题中需要处理大量强震动数据的分析研究工作是有利的。文中还就双线型恢复力模型具体计算了2个例子（一个为屈服强化情形，另一个为理想弹塑性情形）。最后联系计算工作和计算结果，对非弹性反应谱计算中的其它有关事项做了补充讨论和说明。  相似文献   

13.
With the increasing emphasis of performance‐based earthquake engineering in the engineering community, several investigations have been presented outlining simplified approaches suitable for performance‐based seismic design (PBSD). Central to most of these PBSD approaches is the use of closed‐form analytical solutions to the probabilistic integral equations representing the rate of exceedance of key performance measures. Situations where such closed‐form solutions are not appropriate primarily relate to the problem of extrapolation outside of the region in which parameters of the closed‐form solution are fit. This study presents a critical review of the closed‐form solution for the annual rate of structural collapse. The closed‐form solution requires the assumptions of lognormality of the collapse fragility and power model form of the ground motion hazard, of which the latter is more significant regarding the error of the closed‐form solution. Via a parametric study, the key variables contributing to the error between the closed‐form solution and solution via numerical integration are illustrated. As these key variables cannot be easily measured, it casts doubt on the use of such closed‐form solutions in future PBSD, especially considering the simple and efficient nature of using direct numerical integration to obtain the solution. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

14.
A closed-form solution of two-dimensional scattering of plane SH waves by a cylindrical hill of circular-arc cross-section in a half-space is presented using the wave functions expansion method. The solution is reduced to solving a set of infinite linear algebraic equations, using the auxiliary functions and the exterior region form of the Graf's addition theorem. Numerical solutions are obtained by truncation of the infinite equations and their accuracies are demonstrated by convergence of the numerical results and by the extent to which the numerical results fit the exact boundary conditions with increasing the truncation order. The numerical results for some typical cases are then presented for checking accuracies of various numerical methods. The effects of the height-to-width ratio of the hill on surface ground motion are finally illustrated.  相似文献   

15.
Nonlinear behavior of liquid sloshing inside a partially filled rectangular tank is investigated. The nonlinearity in the numerical modeling of the liquid sloshing originates from the nonlinear terms of the governing equations of the fluid flow and the liquid free surface motion as a not known boundary condition. The numerical simulations are performed for both linear and nonlinear conditions. The computed results using linear conditions are compared with readily available exact solution. In order to verify the results of the nonlinear numerical solution, a series of the shaking table tests on rectangular tank were conducted. Having verified linear and nonlinear numerical models, they are used for computation of near wall sloshing height at a series of real scale tanks (with various dimensions) under the both harmonic and earthquake base excitation. Finally, the nonlinear effects on liquid sloshing modeling are discussed and the practical limitations of the linear solution in evaluating the response of seismically excited liquids are also addressed.  相似文献   

16.
Summary A unified analysis has been made to obtain all possible similarity solutions of the steady and unsteady, forced flow, inside a boundary layer along a flat plate. Though previously, attempts were made to obtain similarity solutions of a steady boundary layer flow neglecting viscous dissipation term in the energy conservation equation but the treatments were not complete. Here we have taken account of the viscous dissipation term. In the steady case it has been shown that for a similarity solution of both velocity and temperature, there should be a relation between the undisturbed flow outside the boundary layer and the temperature of the plate. It has been shown that the similarity solution exists in the unsteady case if we neglet the viscous dissipation term in the energy equation.  相似文献   

17.
After 50 years of Prabhu’s paper on the exact solution of the stochastic reservoir equation for the important class of gamma inflow distributions with an integral shape parameter, a detailed implementation of the exact solution is still lacking, despite its potential usefulness from both theoretical and practical points of view. This paper explores some properties of Prabhu’s exact solution and investigates the numerical difficulties associated with its implementation. The solution is also extended to derive the distributions of deficit, spillage, yield, and actual release from the reservoir. Explicit analytical solutions for three relatively simple cases are given in detail as examples and comparisons with approximate numerical solutions are made, which reveal some shortcomings of approximate methods. The implementation of the solution in the general case reveals some numerical problems associated with large values of the shape parameter of the inflow distribution and large ratios of reservoir size to draft, mainly due to accumulation of round-off errors. A Matlab program has been developed to calculate emptying and filling probabilities over a wide range of reservoir parameters using extended precision. Comparison of Prabhu’s solution with the numerical solution of the reservoir integral equation highlights possible problems with the numerical solution, which may produce inaccurate or even invalid results for large reservoirs, large drift, and large skewness of the inflow distribution. A comparison between gamma and lognormal distributions as models of skew revealed that as the reservoir size, drift, and skewness increase, the probability of emptying of the reservoir becomes smaller for the case of gamma inflow than in the case of lognormal flow having the same skewness coefficient.  相似文献   

18.
Analytical solutions for the water flow and solute transport equations in the unsaturated zone are presented. We use the Broadbridge and White nonlinear model to solve the Richards’ equation for vertical flow under a constant infiltration rate. Then we extend the water flow solution and develop an exact parametric solution for the advection-dispersion equation. The method of characteristics is adopted to determine the location of a solute front in the unsaturated zone. The dispersion component is incorporated into the final solution using a singular perturbation method. The formulation of the analytical solutions is simple, and a complete solution is generated without resorting to computationally demanding numerical schemes. Indeed, the simple analytical solutions can be used as tools to verify the accuracy of numerical models of water flow and solute transport. Comparison with a finite-element numerical solution indicates that a good match for the predicted water content is achieved when the mesh grid is one-fourth the capillary length scale of the porous medium. However, when numerically solving the solute transport equation at this level of discretization, numerical dispersion and spatial oscillations were significant.  相似文献   

19.
The development of a displacement finite element formulation and its application to convective transport problems is presented. The formulation is based on the introduction of a generalized quantity defined as transport displacement. The governing equation is expressed in terms of this quantity and by using generalized coordinates a variational form of the governing equation is obtained. This equation may be solved by any numerical method, though it is of particular interest for application of the finite element method. Two finite element models are derived for the solution of convection-diffusion boundary value problems. The performance of the two element models is discussed and numerical results are given for different cases of convection and diffusion with two types of boundary conditions. The numerical results obtained show not only the efficiency of the numerical models in handling pure convection, pure diffusion and mixed convection-diffusion problems, but also good stability and accuracy. The applications of the developed numerical models are not limited to diffusion-convection problems but can also be applied to other types of problems such as mass transfer, hydrodynamics and wave propagation.  相似文献   

20.
A boundary element formulation for the dynamic analysis of axially and laterally loaded single piles and pile groups is presented. The piles are represented by compressible beam-column elements and the soil as a hysteretic elastic half-space. The governing equations of motion for the pile domain have been solved exactly for distributed periodic loading intensities. These solutions are then coupled with a numerical solution for the motion of the soil domain by satisfying equilibrium and compatibility at the pile-soil interface. The results obtained from the analysis compare favourably with those from alternative analyses, e.g. finite element, but at greatly reduced computational costs.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号