近海核电厂核素地下水释放通量的模型计算方法     

朱君  李婷  陈超  谢添  张艾明 《吉林大学学报(地球科学版)》2021,51(1):201-211

为了定量计算陵区近海核电站排水管线泄漏情景下核素通过地下水途径向海洋环境的释放通量，以某近海核电站为例进行研究。首先，应用GOCAD软件建立三维地形地质模型，刻画地层的分布、剥蚀以及倾向等特点；然后，运用地下水数值模拟软件FEFLOW精细刻画丘陵区地下水系统的补给、径流和排泄特征；最后，以不被吸附滞留的核素³H和被吸附滞留的核素⁹⁰Sr、¹³⁷Cs为对象，通过实验测定了⁹⁰Sr、¹³⁷Cs在不同岩土介质中的分配系数，模拟计算了排水管线连续渗漏60 a后³H、⁹⁰Sr、¹³⁷Cs在地下水中的放射性分布及释放。结果表明：³H迁移速度基本与地下水流速一致，地下水中的最大放射性浓度为0.285 0 Bq/L，第20 000天时向收纳水域的释放通量达到最大值，约526 Bq/d；⁹⁰Sr吸附性能相对较弱，最大迁移距离约80 m，地下水中的最大放射性浓度为0.032 1 Bq/L；¹³⁷Cs吸附能力较强，相当长的时间内被滞留在管线附近，地下水中最大放射性浓度分别为6.840×10^-3 Bq/L，释放通量为0 Bq/d。由弥散度的不确定分析可知，弥散度越大，地下水中³H的最大放射性浓度越小，向海洋环境的释放通量越多。  相似文献   

A note on the performance of the multiply-upstream semi-Lagrangian advection schemes for one-dimensional nonlinear momentum conservation equation     

Prof. Z. I. Janjić 《Meteorology and Atmospheric Physics》1995,55(1-2):1-16

Summary Two-time-level multiply-upstream semi-Lagrangian schemes were examined in the case of the self-advecting, one-dimensional nonlinear momentum conservation equation. The shock formation process was analyzed. It is pointed out that the shocks cannot be created in the truncated systems satisfying the Pudykiewicz, Benoit and Staniforth criterion.The numerical integrations were restricted to 12 h. It was shown that, at least in the sub-CFL range, increased complexity of the scheme can compensate reduced horizontal resolution. A considerable sensitivity of the schemes with respect to the time step was detected. In the super-CFL mode, several windows on various time scales were found within which the Pudykiewicz, Benoit and Staniforth criterion was satisfied. The time step of 1.44 times the maximum time step allowed by the CFL criterion was used in the semi-Lagrangian runs.The super-CFL, semi-Lagrangian solutions were diverging progressively from the sub-CFL ones as the forecasts advanced. This was also reflected in the energy spectra.Unacceptably large energy losses were encountered in the super-CFL, semi-Lagrangian runs. Most of these losses could be explained by the reduced mean wind speed, i.e., the amplitude of the zero wavenumber wave. At the same time, the energy content in the shorter waves increased. In a more complex model, such a situation would resemble a loss of zonal, and an increase of transient eddy kinetic energy.A trajectory error measure was defined as the maximum absolute value of the distance between the actual arriving point of the particle originating at the estimated departure point, and the grid point assumed to be the arrival point in the semi-Lagrangian procedure. In contrast to the sub-CFL regime, this measure could reach a considerable fraction of the grid distance in the computations with the super-CFL time steps.In the physical system considered, the trajectories are determined only by the velocities at the departure points. With the semi-Lagrangian schemes the distances traveled by the particles are estimated on the basis of the velocities at the points downstream with respect to the departure points. Thus, unless the solution is smooth (in space and time) on the scales of the extrapolation distances/times, the upstream extrapolation does not promise the convergence of the solution.With 16 Figures  相似文献   

The effect of an optically thin region on stability of an accretion disk around a supermassive black hole     

B. F. Liu  G. Z. Xie  K. F. Ji 《Astrophysics and Space Science》1995,225(2):299-314

We analyse conditions of the innermost portion of an accretion disk and establish a set of equations for this region. A stable innermost region may exist, which can probably explain the observed UV and X-ray spectra, avoiding the unstable emission. We then discuss the detailed radial structure of a disk around a black hole for typical AGN parameters and obtain different kinds of- relationships for different regions of a disk. On the basis of this, we discuss the stability. A new type of cycle is present, which we call a double S shaped cycle. In this cycle, the extent of accretion rate variability is much larger than that in dwarf nova cycles. This probably solves the problem of violent variability of AGN. In the meantime, the very high accretion rate at the hottest state in limit cycles in the unstable region may provide continuous injection of matter to the jet and power the relativistic motion of the jet.  相似文献   

On the internal structure and magnetic fields of     

C. Z. Zhang  K. Zhang 《Earth, Moon, and Planets》1995,69(3):237-247

First, a sequence of four-zone models for the interior of Venus is constructed under the assumption of hydrostatic equilibrium. While the equation of state for each zone is taken to be the Bullen's relation with its coefficients consistent with the PREM Earth model (Dziewonski and Anderson, 1981), the position of core-mantle boundary is determined by matching solutions of the Emden's equation in different regions. The results of hydrostatic models indicate the presence of a reasonably large molten iron core in Venus, broadly similar to the Earth. It is also found that the position of the core-mantle interface is nearly model-independent. Second, we focus on the question why Venus does not possess a significant global magnetic field and on what we can learn from this fact. Solutions of magnetohydrodynamic equations appropriate for the molten core of Venus are discussed. It is argued that, because the Elsasser number measuring the relative importance of Coriolis and Lorentz forces satisfies 1, equations for the problem of thermal convection in the Venusian fluid core must be nearly uncoupled with the dynamo equation. The existence of a global magnetic field, though small, then suggests that the size of the magnetic Reynolds numberR _m must beR _m =O(10), sustaining a dynamo action near its marginal state but not an active dynamo in the Venusian molten core. On the basis of asymptotic relations for finite amplitude convection, a useful constraint on important physical parameters for the liquid core of Venus is derived and discussed.  相似文献   

Kleť minor planet discovery frequency by months (comparison)     

J. Tichá  M. Tichý  Z. Moravec 《Earth, Moon, and Planets》1995,68(1-3):579-583

Abstrakt The Kle search programme for minor planets using 0,63m Maksutov camera of Kle Observatory was begun in 1977 as a long-term project continuing up to the present day. The discovery frequency according to seasons in the year shows multiple differences between particular months.The comparison with the number of observing nights and obtained search plates; weather conditions and other influences will be discussed. The Kle minor planet discovery frequency also will be compared with the data base of all numbered minor planets.  相似文献   

Fitting matrix-valued variogram models by simultaneous diagonalization (Part I: Theory)     

Tailiang Xie  Donald E. Myers 《Mathematical Geosciences》1995,27(7):867-875

Suppose that ¯(x₁),...,¯Z(x_n). are observations of vector-valued random function ¯(x). In the isotropic situation, the sample variogram γ^*(h) for a given lag h is $$\bar \gamma ^ * (h) = \frac{1}{{2N(h)}}\mathop \sum \limits_{s(h)} (\overline Z (x_1 ) - \overline Z (x_1 )) \overline {(Z} (x_1 ) - \overline Z (x_1 ))^T $$ where s(h) is a set of paired points with distance h and N(h) is the number of pairs in s(h).. For a selection of lags h₁, h₂, .... h_k such that N (h₁) > O. we obtain a ktuple of (semi) positive definite matrices $\bar \gamma ^ * (h_{ 1} ),. . . ., \bar \gamma ^ * (h_{ k} )$ . We want to determine an orthonormal matrix B which simultaneously diagonalizes the $\bar \gamma ^ * (h_{ 1} ),. . . ., \bar \gamma ^ * (h_{ k} )$ or nearly diagonalizes them in the sense that the sum of squares of offdiagonal elements is small compared to the sum of squares of diagonal elements. If such a B exists, we linearly transform $\overline Z (x)$ by $\overline Y (x) = B\overline Z (x)$ . Then, the resulting vector function $\overline Y (x)$ has less spatial correlation among its components than $\overline Z (x)$ does. The components of $\overline Y (x)$ with little contribution to the variogram structure may be dropped, and small crossvariograms fitted by straightlines. Variogram models obtained by this scheme preserve the negative definiteness property of variograms (in the matrix-valued function sense). A simplified analysis and computation in cokriging can be carried out. The principles of this scheme arc presented in this paper.  相似文献   

On transfer equation in a semi-infinite atmosphere with albedo ω>1     

S. R. Das Gupta  Z. Islam  S. K. Bishnu 《Astrophysics and Space Science》1980,70(1):191-209

In this note we derive an exact solution of transfer equation in a plane-parallel semiinfinite atmosphere with albedo >1, by the method of Laplace transform and Wiener-Hopf technique. The emergent intensityI(0, ) is obtained in terms of theH ⁰-functionH ⁰() (Das Gupta, 1978) for which some good approximations are given. Intensity at any depth is also obtained.I(0, )/I(0, 0) is plotted in graphs against [0,1], and shows a maximum which drops and shifts towards the origin as increases.  相似文献   

Observations and interpretation of the close binary system CG Cyg     

D. M. Z. Jassur 《Astrophysics and Space Science》1980,67(1):19-30

Photometric observations of the short-period (RS CVn-type) eclipsing binary system CG Cyg have been presented. Two sets of results, obtained from an analysis of theB, V andR light curves, represent ‘occultation’ and ‘transit’ solutions. The occultation solution is preferred as it gives a better fit to the colour curve. This hypothesis may also offer a more promising explanation of the observed peculiarities such as period changes and the light variation outside eclipses.  相似文献   

Transequatorial loops interconnecting McMath regions 12472 and 12474     

Švestka  Z.  Krieger  A. S.  Chase  R. C.  Howard  R. 《Solar physics》1977,52(1):69-90

We discuss the life-story of a transequatorial loop system which interconnected the newly born active region McMath 12474 with the old region 12472. The loop system was probably born through reconnection accomplished 1.5 to 5 days after the birth of 12474 and the loops were observed in soft X-rays for at least 1.5 days. Transient sharpenings of the interconnection and a striking brightening of the whole loop system for about 6 hr appear to be caused by magnetic field variations in the region 12474. A flare might have been related to the brightening, but only in an indirect way: the same emerging flux could have triggered the flare and at the same time strengthened the magnetic field at the foot-points of the loops. Electron temperature in the loop system, equal to 2.1 × 10⁶ K in its quiet phase, increased to 3.1 × 10⁶ K during the brightening. Electron density in the loop system was 1.3 × 10⁹ cm^–3 and it could be estimated to 7 × 10⁸ cm^–3 prior to the brightening. During the brightening the loops became twisted. There was no obvious effect whatsoever of the activity in 12474 upon the in erconnected old region. The final decay of the loop system reflected the decay of magnetic field in the region 12474.  相似文献   

Unusual coronal activity following the flare of 6 November 1980     

Z. Švestka  B. R. Dennis  M. Pick  A. Raoult  C. G. Rapley  R. T. Stewart  B. E. Woodgate 《Solar physics》1982,80(1):143-159

For almost 30 hr after the major (gamma-ray) two-ribbon flare on 6 November 1980, 03:30 UT, the Hard X-Ray Imaging Spectrometer (HXIS) aboard the SMM satellite imaged in > 3.5 keV X-rays a gigantic arch extending above the active region over the limb. Like a similar configuration on 22 May 1980, this arch formed the lowest part of a stationary post-flare radio noise storm recorded at metric wavelengths at Nançay and Culgoora. 6.5 hr after the flare a coronal region below the arch started quasi-periodic pulsations in X-ray brightness, observed by several SMM instruments. These brightness variations had no response in the chromosphere (H), very little in the transition layer (O v), but they clearly correlated with similar variations in brightness at 169 MHz. There were 13 pulses of this kind, with apparent periodicity of about 20 min, until another flare occurred in the active region at 15:00 UT. All the brightenings appeared within a localized area of about 30000 km² in the northern part of the active region, but they definitely did not occur all at the same place.The top of the X-ray arch, at an altitude of 155 000 km, was continuously and smoothly decaying, taking no part in the striking variations below it. Therefore, the area variable in brightness does not seem to be the footpoint of the arch, as we supposed for similar variations on 22 May. More likely, it is a separate region connected directly with the source of the radio storm; particles accelerated in the storm may be dumped into the low corona and cause the X-ray enhancements. The X-ray arch was enhanced by two orders of magnitude in 3.5–5.5 keV X-ray counts and the temperature increased from 7.3 × 10⁶ to 9 × 10⁶ K when the new two-ribbon flare occurred at 15:00 UT. Thus, it is possible that energy is brought into the arch via the upper parts of the reconnecting flare loops - a process that can continue for hours.  相似文献