共查询到20条相似文献，搜索用时 46 毫秒
1.
In this work we study mixed finite element approximations of Richards’ equation for simulating variably saturated subsurface flow and simultaneous reactive solute transport. Whereas higher order schemes have proved their ability to approximate reliably reactive solute transport (cf., e.g. [Bause M, Knabner P. Numerical simulation of contaminant biodegradation by higher order methods and adaptive time stepping. Comput Visual Sci 7;2004:61–78]), the Raviart–Thomas mixed finite element method (RT0) with a first order accurate flux approximation is popular for computing the underlying water flow field (cf. [Bause M, Knabner P. Computation of variably saturated subsurface flow by adaptive mixed hybrid finite element methods. Adv Water Resour 27;2004:565–581, Farthing MW, Kees CE, Miller CT. Mixed finite element methods and higher order temporal approximations for variably saturated groundwater flow. Adv Water Resour 26;2003:373–394, Starke G. Least-squares mixed finite element solution of variably saturated subsurface flow problems. SIAM J Sci Comput 21;2000:1869–1885, Younes A, Mosé R, Ackerer P, Chavent G. A new formulation of the mixed finite element method for solving elliptic and parabolic PDE with triangular elements. J Comp Phys 149;1999:148–167, Woodward CS, Dawson CN. Analysis of expanded mixed finite element methods for a nonlinear parabolic equation modeling flow into variably saturated porous media. SIAM J Numer Anal 37;2000:701–724]). This combination might be non-optimal. Higher order techniques could increase the accuracy of the flow field calculation and thereby improve the prediction of the solute transport. Here, we analyse the application of the Brezzi-Douglas-Marini element (BDM1) with a second order accurate flux approximation to elliptic, parabolic and degenerate problems whose solutions lack the regularity that is assumed in optimal order error analyses. For the flow field calculation a superiority of the BDM1 approach to the RT0 one is observed, which however is less significant for the accompanying solute transport. 相似文献
2.
Groundwater models need to account for detailed but generally unknown spatial variability (heterogeneity) of the hydrogeologic model inputs. To address this problem we replace the large, m-dimensional stochastic vector β that reflects both small and large scales of heterogeneity in the inputs by a lumped or smoothed m-dimensional approximation γθ∗, where γ is an interpolation matrix and θ∗ is a stochastic vector of parameters. Vector θ∗ has small enough dimension to allow its estimation with the available data. The consequence of the replacement is that model function f(γθ∗) written in terms of the approximate inputs is in error with respect to the same model function written in terms of β, f(β), which is assumed to be nearly exact. The difference f(β) − f(γθ∗), termed model error, is spatially correlated, generates prediction biases, and causes standard confidence and prediction intervals to be too small. Model error is accounted for in the weighted nonlinear regression methodology developed to estimate θ∗ and assess model uncertainties by incorporating the second-moment matrix of the model errors into the weight matrix. Techniques developed by statisticians to analyze classical nonlinear regression methods are extended to analyze the revised method. The analysis develops analytical expressions for bias terms reflecting the interaction of model nonlinearity and model error, for correction factors needed to adjust the sizes of confidence and prediction intervals for this interaction, and for correction factors needed to adjust the sizes of confidence and prediction intervals for possible use of a diagonal weight matrix in place of the correct one. If terms expressing the degree of intrinsic nonlinearity for f(β) and f(γθ∗) are small, then most of the biases are small and the correction factors are reduced in magnitude. Biases, correction factors, and confidence and prediction intervals were obtained for a test problem for which model error is large to test robustness of the methodology. Numerical results conform with the theoretical analysis. 相似文献
3.
4.
We compare two methods for determining the upscaled water characteristics and saturation-dependent anisotropy in unsaturated hydraulic conductivity from a field-scale injection test. In both approaches an effective medium approximation is used to reduce a porous medium of M textures to an equivalent homogenous medium. The first approach is a phenomenological approach based on homogenization and assumes that moisture-based Richards’ equation can be treated like the convective–dispersive equation (CDE). The gravity term, dKz(θ)/d(θ), analogous to the vertical convective velocity in the CDE, is determined from the temporal evolution of the plume centroid along the vertical coordinate allowing calculation of an upscaled Kz(θ). As with the dispersion tensor in the CDE, the rate of change of the second spatial moment in 3D space is used to calculate the water diffusivity tensor, D(θ), from which an upscaled K(θ) is calculated. The second approach uses the combined parameter scale inverse technique (CPSIT). Parameter scaling is used first to reduce the number of parameters to be estimated by a factor M. Upscaled parameters are then optimized by inverse modeling to produce an upscaled K(θ) characterized by a pore tortuosity–connectivity tensor, L. Parameters for individual textures are finally determined from the optimized parameters by inverse scaling using scale factors determined a priori. Both methods produced upscaled K(θ) that showed evidence of saturation dependent anisotropy. Flow predictions with the STOMP simulator, parameterized with upscaled parameters, were compared with field observations. Predictions based on the homogenization method were able to capture the mean plume behavior but could not reproduce the asymmetry caused by heterogeneity and lateral spreading. The CPSIT method captured the effects of heterogeneity and anisotropy and reduced the mean squared residual by nearly 90% compared to local-scale and upscaled parameters from the homogenization method. The Pacific Northwest National Laboratory is operated for the US Department of Energy by Battelle under Contract DE-AC05-76RL01830. 相似文献
5.
Hamid Zafarani Assadollah NoorzadAnooshirvan Ansari Khosro Bargi 《Soil Dynamics and Earthquake Engineering》2009
Strong-motion data from eight significant well-documented earthquakes in Iran have been simulated using a stochastic modeling technique for finite faults proposed by Beresnev and Atkinson [Bull Seismol Soc Am 87 (1997) 67–84; Seism Res Lett 69 (1998) 27–32]. The database consists of 61 three-component records from eight earthquakes of magnitude ranging from M 6.3 to M 7.4, recorded at hypocentral distances up to 200 km. The model predictions are in good agreement with available Iranian strong-motion data as evidenced by near-zero average of differences between logarithms of the observed and predicted values for all frequencies. The strength factor, sfact, a quantity that controls the high-frequency radiation from the source is determined, on an event-by-event basis, by fitting simulated to observed response spectra. 相似文献
6.
7.
Pterois volitans and Pterois miles, two sub-species of lionfish, have become the first non-native, invasive marine fish established along the United States Atlantic coast and Caribbean. The route and timing of the invasion is poorly understood, however historical sightings and captures have been robustly documented since their introduction. Herein we analyze these records based on spatial location, dates of arrival, and prevailing physical factors at the capture sights. Using a cellular automata model, we examine the relationship between depth, salinity, temperature, and current, finding the latter as the most influential parameter for transport of lionfish to new areas. The model output is a synthetic validated reproduction of the lionfish invasion, upon which predictive simulations in other locations can be based. This predictive model is simple, highly adaptable, relies entirely on publicly available data, and is applicable to other species. 相似文献
8.
Gao Chen Jason J. Venkiteswaran Sherry L. Schiff William D. Taylor 《Aquatic Sciences - Research Across Boundaries》2014,76(3):313-329
A transient model, hereafter referred to as ROM-TM, was developed to quantify river ecosystem metabolic rates and reaeration rates from field observation of changes in dissolved O2 (DO) and the ratio of 18O to 16O in DO (δ18O-DO). ROM-TM applies an inverse modeling approach and is programmed using MATLAB. Parameters describing photosynthesis, ecosystem respiration, gas exchange, and isotopic fractionation, such as maximum photosynthetic rate (P m ), photosynthetic efficiency parameter (a), respiration rate at 20 °C (R 20 ), gas exchange coefficient (K), respiration isotopic fractionation factor (a R ), and photorespiration coefficient (β R ), can be abstracted by minimizing the sum of square errors between the fitted data and the observed field data. Then DO and δ18O-DO time series can be reconstructed using estimated parameters and input variables. Besides being capable of teasing apart metabolic processes and gas exchange to provide daily average estimates of metabolic parameters at the ecosystem scale, ROM-TM can be used to address issues related to light including light saturation phenomena at the ecosystem level, the effect of cloud cover on the metabolic balance, and photorespiration. Error and uncertainty analysis demonstrates that ROM-TM is stable and robust for the random errors of DO time series. The photosynthetic parameters P m and a are more sensitive than other parameters to lower-resolution time series data. 相似文献
9.
The variations in the density of the ionospheric F2 layer maximum (NmF2) under the action of the zonal plasma drift perpendicularly to the magnetic (B) and electric (E) fields in the direction geomagnetic west-geomagnetic east have been studied using the three-dimensional nonstationary theoretical model of electron and ion densities (N e and N i ) and temperatures (T e and T i ) in the low-latitude and midlatitude ionospheric F region and plasmasphere. The method of numerical calculations of N e , N i , T e , and T i , including the advantages of the Lagrangian and Eulerian methods, is used in the model. A dipole approximation of the geomagnetic field (B), taking into account the non-coincidence of the geographic and geomagnetic poles and differences between the positions of the Earth’s and geomagnetic dipole centers, is accepted in the calculations. The calculated NmF2 and altitudes of the F2 layer maximum (hmF2) have been compared with these quantities measured at 16 low-latitude ionospheric sounding stations during the geomagnetically quiet period October 11–12, 1958. This comparison made it possible to correct the input model parameters: the NRLMSISE-00 model [O], the meridional component of the neutral wind velocity according to the HWW90 model, and the meridional component of the equatorial plasma drift due to the electric field specified by the empirical model. It has been indicated that the effect of the zonal E × B plasma drift on NmF2 can be neglected under daytime conditions and changes in NmF2 and hmF2 under the action of this drift are insignificant under nighttime conditions north of 25° and south of ?26° geomagnetic latitude. The effect of the zonal E × B plasma drift on NmF2 and hmF2 is most substantial in the nightside ionosphere approximately from ?20° to 20° geomagnetic latitude, and the neglect of this drift results in an up to 2.4-fold underestimation of NmF2. The found dependence of the effect of the zonal E × B plasma drift on NmF2 and hmF2 on geomagnetic latitude is related to the longitudinal asymmetry of B, asymmetry of the neutral wind about the geomagnetic equator, and changes in the meridional E × B plasma drift at a change in geomagnetic longitude. 相似文献
10.
V. V. Vodyannikov G. I. Gordienko S. A. Nechaev O. I. Sokolova S. Yu. Khomutov A. F. Yakovets 《Geomagnetism and Aeronomy》2006,46(6):809-813
The time derivative (d H/dt) of the geomagnetic field horizontal component (H) for seven intervals of high geomagnetic activity in 2003–2005 has been calculated, based on the data of Alma-Ata, Novosibirsk, and Irkutsk observatories, in order to estimate the probability of appearance of geomagnetically induced current (GIC), the value of which is linearly dependent on d H/dt, in power lines on the territory of Kazakhstan. The distributions of the H and d H/dt directions have been constructed; in this case it was most interesting that these distributions were narrow and extended along the magnetic meridian for Alma-Ata and were wider angular for Novosibirsk and Irkutsk. It has been indicated that large H values, determining significant GIC values, took place at a sudden commencement of strong storms, which had a character of a pulsed disturbance of the geomagnetic field, and during large-amplitude geomagnetic field pulsations. The duration of the periods of large |d H/dt| values exceeding the threshold (30 nT/min), when GIC could cause unwanted consequences in power circuits, has been determined. 相似文献
11.
Sh. A. Mukhamediev 《Izvestiya Physics of the Solid Earth》2014,50(5):655-691
Rock masses contain ubiquitous multiscale heterogeneities, which (or whose boundaries) serve as the surfaces of discontinuity for some characteristics of the stress state, e.g., for the orientation of principal stress axes. Revealing the regularities that control these discontinuities is a key to understanding the processes taking place at the boundaries of the heterogeneities and for designing the correct procedures for reconstructing and theoretical modeling of tectonic stresses. In the present study, the local laws describing the refraction of the axes of extreme principal stresses T 1 (maximal tension in the deviatoric sense) and T 3 (maximal compression) of the Cauchy stress tensor at the transition over the elementary area n of discontinuity whose orientation is specified by the unit normal n are derived. It is assumed that on the area n of discontinuity, frictional contact takes place. No hypotheses are made on the constitutive equations, and a priori constraints are not posed on the orientation on the stress axes. Two domains, which adjoin area n on the opposite sides and are conventionally marked + and ?, are distinguished. In the case of the two-dimensional (2D) stress state, any principal stress axis on passing from domain ? to domain + remains in the same quadrant of the plane as the continuation of this axis in domain +. The sign and size of the refraction angle depend on the sign and amplitude of the jump of the normal stress, which is tangential to the surface of discontinuity. In the three-dimensional (3D) case, the refraction of axes T 1 and T 3 should be analyzed simultaneously. For each side, + and ?, the projections of the T 1 and T 3 axes on the generally oriented plane n form the shear sectors S + and S ?, which are determined unambiguously and to whose angular domains the possible directions p + and p ? of the shear stress vectors belong. In order for the extreme stress axes T 1 + ,T 3 + and T 1 ? , T 3 ? to be statically compatible on the generally oriented plane n, it is required that sectors S + and S ? had a nonempty intersection. The direction vectors p + and p ? are determined uniquely if, besides axes T 1 ? , T 3 ? and T 1 + , T 3 + , also the ratios of differential stresses R + and R ? (0 ≤ R ± ≤ 1) are known. This is equivalent to specifying the reduced stress tensors T R + and T R ? The necessary condition for tensors T R + and T R ? being statically compatible on plane n is the equality p + = p ?. In this paper, simple methods are suggested for solving the inverse problem of constructing the set of the orientations of the extreme stress axes from the known direction p of the shear stress vector on plane n and from the data on the shear sector. Based on these methods and using the necessary conditions of local equilibrium on plane n formulated above, all the possible orientations of axes T 1 + , T 3 + are determined if the projections of axes T 1 ? , T 3 ? axes on side — are given. The angle between the projections of axes T 1 + , T 1 ? and/or T 3 + , T 3 ? on the plane can attain 90°. Besides the general case, also the particular cases of the contact between the degenerate stress states and the special position of plane n relative to the principal stress axes are thoroughly examined. Generalization of the obtained results makes it possible to plot the local diagram of the orientations of axes T 1 + , T 3 + for a given sector S ?. This diagram is a so-called stress orientation sphere, which is subdivided into three pairs of areas (compression, tension, and compression-extension). The tension and compression zones cannot contain the poles of T 3 + and T 1 + axes, respectively. The compression-extension zones can contain the poles of either T 1 + or T 3 + axis but not both poles simultaneously. In the particular case when the shear stress vector has a unique direction p ? on side ?, the areas of compression-extension disappear and the diagram is reduced to a beach-ball plot, which visualizes the focal mechanism solution of an earthquake. If area n is a generally oriented plane and if the orientation of the pairs of the statically compatible axes T 1 ? , T 3 ? and T 1 + , T 3 + is specified, then, the stress values on side + are uniquely determined from the known stress values on side ?. From the value of differential stress ratio R ?, one can calculate the value of R +, and using the values of the principal stresses on side ?, determine the total stress tensor T + on side +. The obtained results are supported by the laboratory experiments and drilling data. In particular, these results disclose the drawbacks of some established notions and methods in which the possible refraction of the stress axes is unreasonably ignored or taken into account improperly. For example, it is generally misleading to associate the slip on the preexisting fault with the orientation of any particular trihedron of the principal stress axes. The reconstruction should address the potentially statically compatible principal stress axes, which are differently oriented on opposite sides of the fault plane. The fact that, based on the orientation of the intraplate principal stresses at the base of the lithosphere, one cannot make a conclusion on the active or passive influence of the mantle flows on the lithospheric plate motion is another example. The present relationships linking the stress values on the opposite sides of the fault plane on which the orientations of the principal stress axes are known demonstrate the incorrectness of the existing methods, in which the reduced stress tensors within the material domains are reconstructed without allowance for the dynamic interaction of these domains with their neighbors. In addition, using the obtained results, one can generalize the notion of the zone of dynamical control of a fault onto the case of the existence of discontinuities in this region and analyze the stress transfer across the system of the faults. 相似文献
12.
Using model simulations, the morphological picture (revealed earlier) of the disturbances in the F 2 region of the equatorial ionosphere under quiet geomagnetic conditions (Q-disturbances) is interpreted. It is shown that the observed variations in the velocity of the vertical E × B plasma drift, related to the zonal E y component of the electric field, are responsible for the formation of Q-disturbances. The plasma recombination at altitudes of the lower part of the F 2 region and the dependence of the rate of this process on heliogeophysical conditions compose the mechanism of Q-disturbance formation at night. The daytime positive Q-disturbances are caused exclusively by a decrease in the upward E × B drift, and this type of disturbances could be related to the known phenomenon of counter electrojet. Possible causes of formation of the daytime negative Q-disturbances are discussed. 相似文献
13.
On the basis of the model of the three-dimensional (3D) generalized Kadomtsev-Petviashvili equation for magnetic field h = B ~/B the formation, stability, and dynamics of 3D soliton-like structures, such as the beams of fast magnetosonic (FMS) waves generated in ionospheric and magnetospheric plasma at a low-frequency branch of oscillations when β = 4πnT/B 2 ? 1 and β > 1, are studied. The study takes into account the highest dispersion correction determined by values of the plasma parameters and the angle θ = (B, k), which plays a key role in the FMS beam propagation at those angles to the magnetic field that are close to π/2. The stability of multidimensional solutions is studied by an investigation of the Hamiltonian boundness under its deformations on the basis of solving of the corresponding variational problem. The evolution and dynamics of the 3D FMS wave beam are studied by the numerical integration of equations with the use of specially developed methods. The results can be interpreted in terms of the self-focusing phenomenon, as the formation of a stationary beam and the scattering and self-focusing of the solitary beam of FMS waves. These cases were studied with a detailed investigation of all evolutionary stages of the 3D FMS wave beams in the ionospheric and magnetospheric plasma. 相似文献
14.
The three-dimensional nonstationary theoretical model of the concentrations and temperatures of electrons and ions in the ionospheric F region and plasmasphere at low and middle latitudes is used to study variations in the concentration NmF2 and height hmF2 of the ionospheric F2 layer under the action of the plasma zonal drift in the direction geomagnetic west-geomagnetic east perpendicularly to the electric E and geomagnetic B fields. The calculated and measured values of NmF2 and hmF2 for 16 ionospheric sounding stations during the quiet geomagnetic period on March 28–29, 1964 at low solar activity are compared. This comparison made it possible to correct the input parameters of the model: [O] from the NRLMSISE-00 model and the meridional component of the neutral wind velocity from the HWW90 model. It is shown that the nighttime NmF2 values decrease up to twice at low solar activity in the low-latitude ionosphere, and the hmF2 values change by up to 16 km, if the plasma zonal E×B drift is not taken into account. Under the daytime conditions, the influence of the plasma zonal E×B drift on NmF2 can be neglected. 相似文献
15.
16.
V. I. Starostenko V. N. Shuman I. K. Pashkevich O. V. Legostaeva A. S. Savchenko 《Izvestiya Physics of the Solid Earth》2013,49(1):144-153
We give an overview of the methods designed for reconstructing close-to-harmonic functions from the magnetic field ΔT. The formula of Yu.P. Tafeev is refined. It is shown that this refined formula directly leads to the relation derived by V.M. Gordin and his colleagues that allows isolating the harmonic component in the function ΔT. V.N. Strakhov’s linearized representation of the function ΔT is immediately derived from the main approximate Tafeev formula for Q ΔT. The experience of using Strakhov’s ΔS function in the interpretation of the magnetic anomaly ΔT generated by the Krivoi Rog structure is described. It is noted that the problem of reconstructing the corresponding harmonic functions from the data of magnetic and gravity surveys has much in common. The specific features of measuring the magnetic field H and magnetic induction B in the material media are considered, and the physical interpretation of these fields is presented. 相似文献
17.
18.
The variations in the electron number density of the ionospheric F2 layer maximum (NmF2) under the action of the zonal plasma drift in the geomagnetic west-geomagnetic east direction perpendicularly to the electric (E) and geomagnetic (B) fields during a geomagnetically quiet period on December 7, 1989, at high solar activity have been studied based on a three-dimensional nonstationary theoretical model of electron number densities and temperatures in the ionospheric F region. Calculated and measured NmF2 values for 12 low-latitude ionospheric sounding stations have been compared. When the zonal E × B plasma drift is ignored, the NmF2 values become smaller by up to a factor of 3 under nighttime conditions in the low-latitude ionosphere. The average effect of the zonal E × B plasma drift on NmF2 in the low-latitude ionosphere is larger during winter nights than under summer nighttime conditions. 相似文献
19.
Correlation between abnormal trends in the spontaneous fields of tectonic plates and strong seismicities 
下载免费PDF全文
下载免费PDF全文 Tectonic activities, electrical structures, and electromagnetic environments are major factors that affect the stability of spontaneous fields. The method of correlating regional synchronization contrasts(CRSC) can determine the reliability of multi-site data trends or shortimpending anomalies. From 2008 to 2013, there were three strong earthquake cluster periods in the North–South seismic belt that lasted for 8–12 months. By applying the CRSC method to analyze the spontaneous field E_(SP) at 25 sites of the region in the past 6 years, it was discovered that for each strong earthquake cluster period, the E_(SP) strength of credible anomalous trends was present at minimum 30%of the stations. In the southern section of the Tan-Lu fault zone, the E_(SP) at four main geoelectric field stations showed significant anomalous trends after June 2015, which could be associated with the major earthquakes of the East China Sea waters(MS7.2) in November 2015 and Japan's Kyushu island(MS7.3) in April 2016. 相似文献
20.
N. A. Barkhatov L. I. Gromova A. E. Levitin S. E. Revunov 《Geomagnetism and Aeronomy》2008,48(6):713-718
The relation of the Kp index of geomagnetic activity to the solar wind electric field (E SW) and the projection of this field onto the geomagnetic dipole has been estimated. An analysis indicated that the southward component of the IMF vector (B z < 0) is the main geoeffective parameter, as was repeatedly indicated by many researchers. The presence of this component in any combinations of the interplanetary medium parameters is responsible for a high correlation between such combinations and geomagnetic activity referred to by the authors of different studies. Precisely this field component also plays the main role in the relation between the Kp index and the relative orientation of E SW and the Earth’ magnetic moment. 相似文献