Structure and energy transfer of unstable hot star winds     

A. Feldmeier 《Astrophysics and Space Science》1994,221(1-2):393-408

Due to the instability of the radiation line force, the winds of hot, luminous stars should show a pronounced time-dependence resulting from the nonlinear growth of initially small perturbations. Following the method of Owocki, Castor & Rybicki (1988), we describe the time-dependent wind structure obtained with an independently developed code. Under the central assumption ofisothermality, our results are in very good agreement with the ones by Owocki et al. We find that the response of the wind to periodic base perturbations remains largely periodic, at least up tor 2...3R _* , with no clear evidence of stochastic behaviour.In order to test the foregoing assumption of isothermality and to compute the X-ray emission from models of structured winds, we have also incorporated theenergy equation into our simulations. We encountered the numerical problem that all radiative cooling zones collapse because of the oscillatory thermal instability (cf. Langer et al. 1981). We present a method to hinder this collapse by changing the cooling function at low temperatures. The resulting wind showsresolved cooling zones; but, for a supergiant wind relatively close to the star (r 10R _* ), the macroscopic wind structure is very similar to isothermal calculations. Most of the hot material is caused by shell-shell collisions.  相似文献   

Estimation of covariance matrix from geochemical data with observations below detection limits     

Chang-Jo F. Chung 《Mathematical Geology》1993,25(7):851-865

Multivariate statistical analyses have been extensively applied to geochemical measurements to analyze and aid interpretation of the data. Estimation of the covariance matrix of multivariate observations is the first task in multivariate analysis. However, geochemical data for the rare elements, especially Ag, Au, and platinum-group elements, usually contain observations the below detection limits. In particular, Instrumental Neutron Activation Analysis (INAA) for the rare elements produces multilevel and possibly extremely high detection limits depending on the sample weight. Traditionally, in applying multivariate analysis to such incomplete data, the observations below detection limits are first substituted, for example, each observation below the detection limit is replaced by a certain percentage of that limit, and then the standard statistical computer packages or techniques are used to obtain the analysis of the data. If a number of samples with observations below detection limits is small, or the detection limits are relatively near zero, the results may be reasonable and most geological interpretations or conclusions are probably valid. In this paper, a new method is proposed to estimate the covariance matrix from a dataset containing observations below multilevel detection limits by using the marginal maximum likelihood estimation (MMLE) method. For each pair of variables, sayY andZ whose observations containing below detection limits, the proposed method consists of three steps: (i) for each variable separately obtaining the marginal MLE for the means and the variances, , , , and forY andZ: (ii) defining new variables by and and lettingA=C+D andB=C–D, and obtaining MLE for variances, and forA andB; (iii) estimating the correlation coefficient _YZ by and the covariance _YZ by . The procedure is illustrated by using a precious metal geochemical data set from the Fox River Sill, Manitoba, Canada.  相似文献   

新疆地区夏季平流层低层风场特征分析及数值模拟     

王晓婷  李森  史小康  黄江平  胡艳冰 《气象科技》2021,49(6):878-884

利用2017年6—8月的FNL再分析资料分析了新疆地区夏季平流层低层风场的时空演变特征，开展了平流层低层风场的高分辨率数值模拟和检验。结果表明：新疆地区夏季平流层纬向风随时间的变化存在经向差异，同时准零风层开始和结束的时间也存在经向差异，准零风层开始时间南部地区早于北部地区，而准零风层结束时间则相反。整个夏季新疆地区上空的准零风层处于70～40 hPa之间，其高度随时间呈先降低而后升高的变化趋势。准零风层数值模拟结果表明，模拟的准零风层参数与探空资料分析结果相比二者存在一致的变化趋势，准零风层起始高度的平均绝对误差为467 m，该高度对应的风速均方根误差为1.75 m/s。  相似文献   

参数的期望估计及其在形变分析中的应用   总被引：1，自引：0，他引：1  

黄杰 《华北地震科学》1997,15(2):25-31

利用“参数的期望估计”能准确定位、定量粗差及参数的期望估计不受粗差影响折特殊性质,寻找地壳变形区域和不变形区域,以确定拟稳点,进而进行拟稳变换,分析地菜变,将是一种有效方法。  相似文献   

Estimation for partially observed Markov processes     

M. E. Thompson  T. N. Kaseke 《Stochastic Hydrology and Hydraulics》1995,9(1):33-47

Many stochastic process models for environmental data sets assume a process of relatively simple structure which is in some sense partially observed. That is, there is an underlying process (X_n, n 0) or (X_t, t 0) for which the parameters are of interest and physically meaningful, and an observable process (Y_n, n 0) or (Y_t, t 0) which depends on the X process but not otherwise on those parameters. Examples are wide ranging: the Y process may be the X process with missing observations; the Y process may be the X process observed with a noise component; the X process might constitute a random environment for the Y process, as with hidden Markov models; the Y process might be a lower dimensional function or reduction of the X process. In principle, maximum likelihood estimation for the X process parameters can be carried out by some form of the EM algorithm applied to the Y process data. In the paper we review some current methods for exact and approximate maximum likelihood estimation. We illustrate some of the issues by considering how to estimate the parameters of a stochastic Nash cascade model for runoff. In the case of k reservoirs, the outputs of these reservoirs form a k dimensional vector Markov process, of which only the kth coordinate process is observed, usually at a discrete sample of time points.  相似文献   

The circumstellar shells of carbon miras     

Martin Groenewegen 《Astrophysics and Space Science》1995,224(1-2):321-324

The relation between mass loss rate and pulsation period in carbon Miras is discussed. The dust mass loss rate is very low (about 2 × 10^–10 M/yr) up to aboutP = 380 days, where there is a sudden increase. ForP > 400 days there is a linear relation between logM andP. The change in the mass loss rate near 380 days may be related to radiation pressure on dust becoming effective in driving the outflow.  相似文献   

短周期仪器记录在远震分析中的应用     

许康生 《地震地磁观测与研究》1992,13(2):15-19

地震泼由多种频率成份组成,在其传播过程中。由于地球介质的滤波作用,随着传播距离的增大。高频成份逐渐被衰减。周期越大的成份传播越远。因此,在地震观测上,不同频带的地震仪用于观测不同震中距范围的地震,使仪器较好地响应,就我国台网而言,短周期仪用于地方震、近震的监测,中长周期仪和长周期仪用于远震,极远震的监测,其中长周期仪偏重于记录极远震。但作者在实际分析工作中注意到,短周期仪的记录应用于远震,极远震的分析,显示出独特的作用,充分挖掘和利用短仪资料,将会更加丰富远震、极远震的震相资料。本文以高台地震台的资料为依据,从震中距和震级角度,与中长仪和长仪作比较,对短仪记录远震、极远震的震相及特征作了统计分析与初步探讨,并对新疆地区“影区”地震S波的记录、日本地区地震_PPcP等震相的出现提出了现象的存在,给出了相应的观测结果,以待步一步研究。  相似文献   

成都城市发展对热岛强度影响的主成份回归分析     

石凤云  唐佑民 《成都信息工程学院学报》1991,(Z1)

本文对成都市总人口、建成区面积等11个因子、作了主成份回归L-S估计和M-估计,讨论了成都城市发展对“热岛”强度的主要影响因子。结果表明,城区房屋建筑面积及总人口数是影响城市气候(气温)的主要因子,其次为城市人口总户数、建成面积等。 文中,对回归方程进行了拟合计算,回归效果比较满意(尤其是稳健回归)。  相似文献   

Comparison of structural and least-squares lines for estimating geologic relations   总被引：1，自引：0，他引：1  

Garnett P. Williams  Brent M. Troutman 《Mathematical Geology》1990,22(8):1027-1049

Two different goals in fitting straight lines to data are to estimate a true linear relation (physical law) and to predict values of the dependent variable with the smallest possible error. Regarding the first goal, a Monte Carlo study indicated that the structural-analysis (SA) method of fitting straight lines to data is superior to the ordinary least-squares (OLS) method for estimating true straight-line relations. Number of data points, slope and intercept of the true relation, and variances of the errors associated with the independent (X) and dependent (Y) variables influence the degree of agreement. For example, differences between the two line-fitting methods decrease as error in X becomes small relative to error in Y. Regarding the second goal—predicting the dependent variable—OLS is better than SA. Again, the difference diminishes as X takes on less error relative to Y. With respect to estimation of slope and intercept and prediction of Y, agreement between Monte Carlo results and large-sample theory was very good for sample sizes of 100, and fair to good for sample sizes of 20. The procedures and error measures are illustrated with two geologic examples.  相似文献   

Reservoir description using dynamic parameterisation selection with a combined stochastic and gradient search     

Lars Kristian Nielsen  Sam Subbey  Mike Christie  Trond Mannseth 《Computational Geosciences》2006,10(3):321-342

There is a correspondence between flow in a reservoir and large scale permeability trends. This correspondence can be derived by constraining reservoir models using observed production data. One of the challenges in deriving the permeability distribution of a field using production data involves determination of the scale of resolution of the permeability. The Adaptive Multiscale Estimation (AME) seeks to overcome the problems related to choosing the resolution of the permeability field by a dynamic parameterisation selection. The standard AME uses a gradient algorithm in solving several optimisation problems with increasing permeability resolution. This paper presents a hybrid algorithm which combines a gradient search and a stochastic algorithm to improve the robustness of the dynamic parameterisation selection. At low dimension, we use the stochastic algorithm to generate several optimised models. We use information from all these produced models to find new optimal refinements, and start out new optimisations with several unequally suggested parameterisations. At higher dimensions we change to a gradient-type optimiser, where the initial solution is chosen from the ensemble of models suggested by the stochastic algorithm. The selection is based on a predefined criterion. We demonstrate the robustness of the hybrid algorithm on sample synthetic cases, which most of them were considered insolvable using the standard AME algorithm.  相似文献