浅水方程数值计算方法的研究   总被引：1，自引：5，他引：1       下载免费PDF全文

张修忠  王光谦  金生 《水科学进展》2003,14(4):417-423

求解以水位为变量的连续方程,并根据Navier-Stokes方程压力修正算法的基本思想,建立了浅水方程的水位修正算法,放宽了对离散时间步长的限制.通过对离散方程系数矩阵的重新构造,建立了高分辨率有限元格式,该格式既具有较高的离散精度又避免了数值解的伪振荡.对动量方程的阻力项做负坡线性化处理,提高了露滩计算的稳定性.数值模拟结果与解析解吻合良好,表明所建立的数值计算方法是正确的和可靠的.  相似文献   

Scaling effects on modeled surface energy-balance components using the NOAH-OSU land surface model   总被引：3，自引：0，他引：3  

V. Sridhar  Ronald L. Elliott  Fei Chen 《Journal of Hydrology》2003,280(1-4):105-123

As surface exchange processes are highly non-linear and heterogeneous in space and time, it is important to know the appropriate scale for the reasonable prediction of these exchange processes. For example, the explicit representation of surface variability has been vital in predicting mesoscale weather events such as late-afternoon thunderstorms initiated by latent heat exchanges in mid-latitude regions of the continental United States. This study was undertaken to examine the effects of different spatial scales of input data on modeled fluxes, so as to better understand the resolution needed for accurate modeling. A statistical procedure was followed to select two cells from the Southern Great Plains 1997 hydrology experiment region, each 20 km×20 km, representing the most homogeneous and the most heterogeneous surface conditions (based on soil and vegetation) within the study region. The NOAH-OSU (Oregon State University) Land Surface Model (LSM) was employed to estimate surface energy fluxes. Three scales of study (200 m, 2 and 20 km) were considered in order to investigate the impacts of the aggregation of input data, especially soil and vegetation inputs, on the model output. Model results of net radiation and latent, sensible and ground heat fluxes were compared for the three scales. For the heterogeneous area, the model output at the 20-km resolution showed some differences when compared with the 200-m and 2-km resolutions. This was more pronounced in latent heat (12% decrease), sensible heat (22% increase), and ground heat flux (44% increase) estimation than in net radiation. The scaling effects were much less for the relatively homogeneous land area with 5% increase in sensible heat and 4% decrease in ground heat flux estimation. All of the model outputs for the 2- and 20-km resolutions were in close agreement. The results suggested that, for this study region, soils and vegetation input resolution of about 2 km should be chosen for realistic modeling of surface exchange processes. This resolution was sufficient to capture the effects of sub-grid scale heterogeneity, while avoiding the data and computational difficulties associated with higher spatial resolutions.  相似文献   

Efficiency of base isolation systems in structural seismic protection and energetic assessment     

Giuseppe Carlo Marano  Rita Greco 《地震工程与结构动力学》2003,32(10):1505-1531

This paper concerns the seismic response of structures isolated at the base by means of High Damping Rubber Bearings (HDRB). The analysis is performed by using a stochastic approach, and a Gaussian zero mean filtered non‐stationary stochastic process is used in order to model the seismic acceleration acting at the base of the structure. More precisely, the generalized Kanai–Tajimi model is adopted to describe the non‐stationary amplitude and frequency characteristics of the seismic motion. The hysteretic differential Bouc–Wen model (BWM) is adopted in order to take into account the non‐linear constitutive behaviour both of the base isolation device and of the structure. Moreover, the stochastic linearization method in the time domain is adopted to estimate the statistical moments of the non‐linear system response in the state space. The non‐linear differential equation of the response covariance matrix is then solved by using an iterative procedure which updates the coefficients of the equivalent linear system at each step and searches for the solution of the response covariance matrix equation. After the system response variance is estimated, a sensitivity analysis is carried out. The final aim of the research is to assess the real capacity of base isolation devices in order to protect the structures from seismic actions, by avoiding a non‐linear response, with associated large plastic displacements and, therefore, by limiting related damage phenomena in structural and non‐structural elements. In order to attain this objective the stochastic response of a non‐linear n‐dof shear‐type base‐isolated building is analysed; the constitutive law both of the structure and of the base devices is described, as previously reported, by adopting the BWM and by using appropriate parameters for this model, able to suitably characterize an ordinary building and the base isolators considered in the study. The protection level offered to the structure by the base isolators is then assessed by evaluating the reduction both of the displacement response and the hysteretic dissipated energy. Copyright © 2003 John Wiley & Sons, Ltd.  相似文献   

两城水源地喀斯特水资源评价     

王延恩  李念国  韩汝钊 《地下水》2003,25(2):103-105

两城水源地具有较丰富的喀斯特水，根据长期观测资料，充分考虑水资源的各种补给和开采条件，应用水均衡原理，对其水资源量进行定量评价，获得允许可开采量，为工农业供水提供可靠依据。  相似文献   

A quantitative study of the energy release in the aftershocks of the Bhuj earthquake, 2001, India, using Lg phase     

G. Jayachandran  M.M. Abdul RazakA.G.V. Prasad  E. Unnikrishnan 《Physics of the Earth and Planetary Interiors》2003,138(2):131-145

The devastating earthquake on 26 January 2001 at Bhuj, India, resulted in large-scale death and destruction of properties of several million US dollars. The moment magnitude of the earthquake was 7.7 and its maximum focal intensity exceeded X in MM scale. The rate of aftershocks of this earthquake, recorded at Gauribidanur seismic array station (GBA), shows a monotonic decay with time superposed with oscillations. For the Indian continent the Lg phase is a prominent arrival at regional distances. The estimate of Lg amplitude is obtained by optimally fitting the Lg wave train to a exponential decay curve. The logarithm of these amplitudes and logarithm of root mean square (rms) value of actual amplitudes of the Lg are calibrated with USGS m_b to create a local m_bLg magnitude scale. The energy released from these aftershocks is calculated from the rms value of Lg phase. The plot of cumulative energy release with time follows the power law of the form t^p, superposed with oscillations. The exponent of the power law, p, is estimated both by a time-window scanning method and by an interpolation method. The value of p is 0.434 for time-window scanning method and 0.432 for the interpolation method. The predominant periods found in the oscillatory part of the cumulative energy, obtained by differencing the observed from the power law fit, are 10.6, 7.9, 5.4, 4.6 and 3.5 h for time-window scanning method. The corresponding periods for interpolation method are 13.4, 11.5, 7.4, 4.2, 3.5, 2.6 and 2.4 h.  相似文献   

Optimal distribution of linear control intensity over the frequency range     

Ottavia Corbi 《地震工程与结构动力学》2003,32(2):173-185

In this paper, a control strategy for structural systems is proposed and developed in the frequency domain. The algorithm is substantially based on a linear derivative feedback and a convolution of the control parameter, whose distribution in the frequency field is chosen in such a manner as to comply with the requirements of an ad hoc formulated constrained optimum problem, with the response data monitored until the instant of control action application. Some numerical testing is carried out by referring to given recorded accelerograms, showing a good performance of the proposed approach. Copyright © 2002 John Wiley & Sons, Ltd.  相似文献   

A procedure for evaluating seismic energy demand of framed structures     

Chung‐Che Chou  Chia‐Ming Uang 《地震工程与结构动力学》2003,32(2):229-244

Energy serves as an alternative index to response quantities like force or displacement to include the duration‐related seismic damage effect. A procedure to evaluate the absorbed energy in a multistorey frame from energy spectra was developed. For low‐ to medium‐rise frames, it required a static pushover analysis of the structure to determine the modal yield force and ductility factor of an equivalent single‐degree‐of‐freedom system for the first two modes. The energy spectra were then used to determine the energy contribution of each mode. A procedure was also developed to distribute the energy along the frame height based on energy shapes. This study showed that the second‐mode response in some cases needs to be considered to reflect the energy (or damage) concentration in the upper floors. Copyright © 2002 John Wiley & Sons, Ltd.  相似文献   

共线条件方程式教学中的几个问题     

冯文灏 《测绘信息与工程》2003,28(2):34-35

指出了共线条件方程式教学中应注意的一些问题：共线条件方程式是联立的两个平面方程式，存在双主距（fx,fy)时的几何概念，以及它的变换式与直接线性变换关系式的异同点。  相似文献   

A simple test and predictive models for assessing swell potential of Ankara (Turkey) Clay     

Zeynal Abiddin Erguler  Resat Ulusay   《Engineering Geology》2003,67(3-4):331-352

During the past three decades, damage due to swelling action of Ankara Clay has been observed more clearly in some parts of Ankara where rapid expansion of the city led to the construction of various kinds of structures. In this study, a comprehensive research program has been conducted (i) to investigate the effect of remoulding and desiccation on the swelling behavior of Ankara Clay and its swelling anisotropy, (ii) to estimate depth of active zone, (iii) to develop a simple technique in determining the magnitude of swelling based on water content of the soaked specimen after 24 and 72 h (w_max24,72), and (iv) to produce predictive models which could be used to estimate the swelling potential of Ankara Clay from its mineralogical and simply measured engineering characteristics. A laboratory testing program was carried out using both undisturbed, and remoulded and desiccated samples selected from 20 different locations. Montmorillonite was identified as being the main clay mineral present. Based on the moisture content variation with depth, the active expansive zone was considered to be about 2 m. The test results suggest that swelling pressure of the clay considerably decreases and/or dies out when the water content is greater than 30%. It is also noted that the measured lateral swelling is significantly in excess of the vertical equivalent indicating an anisotropy with respect to swelling. Statistical assessments indicate that a new parameter (w_max24,72) from a simple test, suggested in this study, is a very strong parameter for predicting swelling parameters of Ankara Clay. Based on 60 empirical predictive equations with coefficients of correlation between 0.96 and 0.66 from multiple regression analyses, w_max24,72, methylene blue value, liquid limit, dry unit weight and smectite content are the most important index and mineralogical properties to predict the swelling parameters of Ankara clay with small deviations from the measured values.  相似文献   

Simulation of marine boundary layer characteristics using a 1-D PBL model over the Bay of Bengal during BOBMEX-99     

N. V. Sam  U. C. Mohanty  A. N. V. Satyanarayana 《Journal of Earth System Science》2003,112(2):185-204

The characteristic features of the marine boundary layer (MBL) over the Bay of Bengal during the southwest monsoon and the factors influencing it are investigated. The Bay of Bengal and Monsoon Experiment (BOBMEX) carried out during July–August 1999 is the first observational experiment under the Indian Climate Research Programme (ICRP). A very high-resolution data in the vertical was obtained during this experiment, which was used to study the MBL characteristics off the east coast of India in the north and south Bay of Bengal. Spells of active and suppressed convection over the Bay were observed, of which, three representative convective episodes were considered for the study. For this purpose a one-dimensional multi-level PBL model with a TKE-ε closure scheme was used. The soundings, viz., the vertical profiles of temperature, humidity, zonal and meridional component of wind, obtained onboard ORV Sagar Kanya and from coastal stations along the east coast are used for the study. The temporal evolution of turbulent kinetic energy, marine boundary layer height (MBLH), sensible and latent heat fluxes and drag coefficient of momentum are simulated for different epochs of monsoon and monsoon depressions during BOBMEX-99.The model also generates the vertical profiles of potential temperature, specific humidity, zonal and meridional wind. These simulated values compared reasonably well with the observations available from BOBMEX.  相似文献