共查询到20条相似文献，搜索用时 93 毫秒
1.
The implementation of Monte Carlo simulations (MCSs) for the propagation of uncertainty in real-world seawater intrusion (SWI) numerical models often becomes computationally prohibitive due to the large number of deterministic solves needed to achieve an acceptable level of accuracy. Previous studies have mostly relied on parallelization and grid computing to decrease the computational time of MCSs. However, another approach which has received less attention in the literature is to decrease the number of deterministic simulations by using more efficient sampling strategies. Sampling efficiency is a measure of the optimality of a sampling strategy. A more efficient sampling strategy requires fewer simulations and less computational time to reach a certain level of accuracy. The efficiency of a sampling strategy is highly related to its space-filling characteristics.This paper illustrates that the use of optimized Latin hypercube sampling (OLHS) strategies instead of the widely employed simple random sampling (SRS) and Latin hypercube sampling (LHS) strategies, can significantly improve sampling efficiency and hence decrease the simulation time of MCSs. Nine OLHS strategies are evaluated including: improved Latin hypercube sampling (IHS); optimum Latin hypercube (OLH) sampling; genetic optimum Latin hypercube (GOLH) sampling; three sampling strategies based on the enhanced stochastic evolutionary (ESE) algorithm namely φp-ESE which employs the φp space-filling criterion, CLD-ESE which utilizes the centered L2-discrepancy (CLD) space-filling criterion, and SLD-ESE which uses the star L2-discrepancy (SLD) space-filling criterion; and three sampling strategies based on the simulated annealing (SA) algorithm namely φp-SA which employs the φp criterion, CLD-SA which uses the CLD criterion, and SLD-SA which utilizes the SLD criterion. The study applies SRS, LHS and the nine OLHS strategies to MCSs of two synthetic test cases of SWI. The two test cases are the Henry problem and a two-dimensional radial representation of SWI in a circular island. The comparison demonstrates that the CLD-ESE strategy is the most efficient among the evaluated strategies. This paper also demonstrates how the space-filling characteristics of different OLHS designs change with variations in the input arguments of their optimization algorithms. 相似文献
2.
A family of unconditionally stable direct integration algorithm with controllable numerical dissipations is proposed. The numerical properties of the new algorithms are controlled by three parameters α, β and γ. By the consistent and stability analysis, the proposed algorithms achieve the second-order accuracy and are unconditionally stable under the condition that α≥-0.5, β≤ 0.5 and γ≥-(1+α)/2. Compared with other unconditionally stable algorithms, such as Chang's algorithms and CR algorithm, the proposed algorithms are found to be superior in terms of the controllable numerical damping ratios. The unconditional stability and numerical damping ratios of the proposed algorithms are examined by three numerical examples. The results demonstrate that the proposed algorithms have a superior performance and can be used expediently in solving linear elastic dynamics problems. 相似文献
3.
At Delft Geotechnics the technique of ground-penetrating radar is in use for the detection of buried objects such as pipes. To enable us to give our ‘measurements in the field’ a more quantitative interpretation than can be deduced from these alone, a series of experiments has been started under well-defined conditions. A cylindrical vessel containing water simulates wet soil. Mounted horizontally above the water surface is a pulsed triangular half-wave dipole which is used as a transmitting antenna (TA). It has a carrier-frequency of about 160 MHz and a pulse repetition-frequency of about 50 kHz. A movable receiving dipole (‘probe’) in the water measures the transverse, mutually orthogonal Eφ,- and Eθ-components of the pulses as a function of probe-position (r, θ, φ) and of the height h of the TA above the water surface. When these pulses are Fourier-transformed, the transverse electric fields Eφ and Eθ at 200 MHz are obtained. The resulting field patterns are compared with computational results on the basis of the theory of the continuous wave, infinitesimal electric dipole (‘point dipole’). It can be concluded that:
- 1 Far-field conditions have not fully developed at a depth of about 2.50 m, the largest value of the radius r at which field patterns were measured, although it represents a distance of about 15 wavelengths.
- 2 The attenuation constant of the tapwater used, as deduced from E-field measurements for θ= 0, 2.50 m < r < 2.75 m, is slightly less than the value measured using a network analyser and air line combination, in agreement with (1).
- 3 E φ field patterns calculated using the value of the conductivity σ corresponding to the former value of the attenuation constant agree reasonably well with the measured patterns for r≤ 2.50 m and for θ < 20° at all antenna heights considered. Calculated Eφ patterns do not agree so well with the measured patterns when h is close to zero. With increasing height the agreement inproves.
- 4 In accordance with the theory of the point-dipole, the angular distribution of the radiation patterns of the TA becomes wider as the frequency decreases.
- 5 The normalized underwater pulse-spectra shift to lower frequencies with increasing r. This can be explained since the attenuation constant of the water rises with rising frequency.
4.
Mark Meadows Don Adams Rich Wright † Ali Tura Steve Cole David Lumley 《Geophysical Prospecting》2005,53(2):205-213
Rock physics analysis plays a vital role in time‐lapse seismic interpretation because it provides the link between changes in rock and fluid properties and the resulting seismic data response. In this case study of the Schiehallion Field, we discuss a number of issues that commonly arise in rock physics analyses for time‐lapse studies. We show that:
- 1 Logarithmic fits of dry bulk (Kdry) and shear (Gdry) moduli vs. effective pressure (Peff) are superior to polynomial fits.
- 2 2D surface fits of Kdry and Gdry over porosity (φ) and effective pressure using all the core data simultaneously are more useful and accurate than separate 1D fits over φ and Peff for each individual core.
- 3 One average set (facies) of Kdry(φ, Peff) and Gdry(φ, Peff) can be chosen to represent adequately the entire Schiehallion reservoir.
- 4 Saturated velocities and densities modelled by fluid substitution of Kdry(φ, Peff), Gdry(φ, Peff) and the dry bulk density ρdry(φ) compare favourably with well‐log velocities and densities.
- 5 P‐ and S‐wave impedance values resulting from fluid substitution of Kdry(φ, Peff), Gdry(φ, Peff) and ρdry(φ) show that the largest impedance changes occur for high porosities and low effective pressures.
- 6 Uncertainties in Kdry(φ, Peff) and Gdry(φ, Peff) derived for individual cores can be used to generate error surfaces for these moduli that represent bounds for quantifying uncertainties in seismic modelling or pressure–saturation inversion.
5.
The damping‐solvent extraction method for the analysis of unbounded visco‐elastic media is evaluated numerically in the frequency domain in order to investigate the influence of the computational parameters—domain size, amount of artificial damping, and mesh density—on the accuracy of results. An analytical estimate of this influence is presented, and specific questions regarding the influence of the parameters on the results are answered using the analytical estimate and numerical results for two classical problems: the rigid strip and rigid disc footings on a visco‐elastic half‐space with constant hysteretic material damping. As the domain size is increased, the results become more accurate only at lower frequencies, but are essentially unaffected at higher frequencies. Choosing the domain size to ensure that the static stiffness is computed accurately leads to an unnecessarily large domain for analysis at higher frequencies. The results improve by increasing artificial damping but at a slower rate as the total (material plus artificial) damping ratio ζt gets closer to 0.866. However, the results do not deteriorate significantly for the larger amounts of artificial damping, suggesting that ζt≈0.6 is appropriate; a larger value is not likely to influence the accuracy of results. Presented results do not support the earlier suggestion that similar accuracy can be achieved by a large bounded domain with small damping or by a small domain with larger damping. Copyright © 2002 John Wiley & Sons, Ltd. 相似文献
6.
Scattering attenuation in short wavelengths has long been interesting to geophysicists. Ultrasonic coda waves, observed as the tail portion of ultrasonic wavetrains in laboratory ultrasonic measurements, are important for such studies where ultrasonic waves interact with small-scale random heterogeneities on a scale of micrometers, but often ignored as noises because of the contamination of boundary reflections from the side ends of a sample core. Numerical simulations with accurate absorbing boundary can provide insight into the effect of boundary reflections on coda waves in laboratory experiments. The simulation of wave propagation in digital and heterogeneous porous cores really challenges numerical techniques by digital image of poroelastic properties, numerical dispersion at high frequency and strong heterogeneity, and accurate absorbing boundary schemes at grazing incidence. To overcome these difficulties, we present a staggered-grid high-order finite-difference (FD) method of Biot’s poroelastic equations, with an arbitrary even-order (2L) accuracy to simulate ultrasonic wave propagation in digital porous cores with strong heterogeneity. An unsplit convolutional perfectly matched layer (CPML) absorbing boundary, which improves conventional PML methods at grazing incidence with less memory and better computational efficiency, is employed in the simulation to investigate the influence of boundary reflections on ultrasonic coda waves. Numerical experiments with saturated poroelastic media demonstrate that the 2L FD scheme with the CPML for ultrasonic wave propagation significantly improves stability conditions at strong heterogeneity and absorbing performance at grazing incidence. The boundary reflections from the artificial boundary surrounding the digital core decay fast with the increase of CPML thicknesses, almost disappearing at the CPML thickness of 15 grids. Comparisons of the resulting ultrasonic coda Q sc values between the numerical and experimental ultrasonic S waveforms for a cylindrical rock sample demonstrate that the boundary reflection may contribute around one-third of the ultrasonic coda attenuation observed in laboratory experiments. 相似文献
7.
Boundary conditions can significantly affect a slope’s behavior under strong earthquakes. To evaluate the importance of boundary conditions for finite element (FE) simulations of a shake-table experiment on the slope response, a validated three-dimensional (3D) nonlinear FE model is presented, and the numerical and experimental results are compared. For that purpose, the robust graphical user-interface “SlopeSAR”, based on the open-source computational platform OpenSees, is employed, which simplifies the effort-intensive pre- and post-processing phases. The mesh resolution effect is also addressed. A parametric study is performed to evaluate the influence of boundary conditions on the FE model involving the boundary extent and three types of boundary conditions at the end faces. Generally, variations in the boundary extent produce inconsistent slope deformations. For the two end faces, fixing the y-direction displacement is not appropriate to simulate the shake-table experiment, in which the end walls are rigid and rough. In addition, the influence of the length of the 3D slope’s top face and the width of the slope play an important role in the difference between two types of boundary conditions at the end faces (fixing the y-direction displacement and fixing the (y, z) direction displacement). Overall, this study highlights that the assessment of a comparison between a simulation and an experimental result should be performed with due consideration to the effect of the boundary conditions. 相似文献
8.
In computational geodynamics, the Finite Element (FE) method is frequently used. The method is attractive as it easily allows employment of body-fitted deformable meshes and a true free surface boundary condition. However, when a Lagrangian mesh is used, remeshing becomes necessary at large strains to avoid numerical inaccuracies (or even wrong results) due to severely distorted elements. For this reason, the FE method is oftentimes combined with the particle-in-cell (PIC) method, where particles are introduced which track history variables and store constitutive information. This implies that the respective material properties have to be interpolated from the particles to the integration points of the finite elements. In numerical geodynamics, material parameters (in particular the viscosity) usually vary over a large range. This may be due to strongly temperature-dependent rheologies (which result in large but smooth viscosity variations) or material interfaces (which result in viscosity jumps). Here, we analyze the accuracy and convergence properties of velocity and pressure of the hybrid FE-PIC method in the presence of large viscosity variations. Standard interpolation schemes (arithmetic and harmonic) are compared to a more sophisticated interpolation scheme which is based on linear least squares interpolation for two types of elements ( \(Q_1P_0\) and \(Q_2P_{-1}\) ). In the case of a smooth viscosity field, the accuracy and convergence is significantly improved by the new interpolation scheme. In the presence of viscosity jumps, the order of accuracy is strongly decreased. 相似文献
9.
Olga-Joan Ktenidou Zafeiria Roumelioti Norman Abrahamson Fabrice Cotton Kyriazis Pitilakis Fabrice Hollender 《Bulletin of Earthquake Engineering》2018,16(6):2311-2336
Accelerometric data from the well-studied valley EUROSEISTEST are used to investigate ground motion uncertainty and variability. We define a simple local ground motion prediction equation (GMPE) and investigate changes in standard deviation (σ) and its components, the between-event variability (τ) and within-event variability (φ). Improving seismological metadata significantly reduces τ (30–50%), which in turn reduces the total σ. Improving site information reduces the systematic site-to-site variability, φ S2S (20–30%), in turn reducing φ, and ultimately, σ. Our values of standard deviations are lower than global values from literature, and closer to path-specific than site-specific values. However, our data have insufficient azimuthal coverage for single-path analysis. Certain stations have higher ground-motion variability, possibly due to topography, basin edge or downgoing wave effects. Sensitivity checks show that 3 recordings per event is a sufficient data selection criterion, however, one of the dataset’s advantages is the large number of recordings per station (9–90) that yields good site term estimates. We examine uncertainty components binning our data with magnitude from 0.01 to 2 s; at smaller magnitudes, τ decreases and φ SS increases, possibly due to κ and source-site trade-offs Finally, we investigate the alternative approach of computing φ SS using existing GMPEs instead of creating an ad hoc local GMPE. This is important where data are insufficient to create one, or when site-specific PSHA is performed. We show that global GMPEs may still capture φ SS , provided that: (1) the magnitude scaling errors are accommodated by the event terms; (2) there are no distance scaling errors (use of a regionally applicable model). Site terms (φ S2S ) computed by different global GMPEs (using different site-proxies) vary significantly, especially for hard-rock sites. This indicates that GMPEs may be poorly constrained where they are sometimes most needed, i.e., for hard rock. 相似文献
10.
B. J. Travis C. Anderson J. Baumgardner C. W. Gable B. H. Hager R. J. O'Connell 《地球物理与天体物理流体动力学》2013,107(3-4):137-160
Abstract A comparison is made between seven different numerical methods for calculating two-dimensional thermal convection in an infinite Prandtl number fluid. Among the seven methods are finite difference and finite element techniques that have been used to model thermal convection in the Earth's mantle. We evaluate the performance of each method using a suite of four benchmark problems, ranging from steady-state convection to intrinsically time-dependent convection with recurring thermal boundary layer instabilities. These results can be used to determine the accuracy of other computational methods, and to assist in the development of new ones. 相似文献
11.
A stochastic approach to nonlinear unconfined flow subject to multiple random fields 总被引：3，自引：3，他引：0
Liangsheng Shi Jinzhong Yang Dongxiao Zhang 《Stochastic Environmental Research and Risk Assessment (SERRA)》2009,23(6):823-835
In this study, the KLME approach, a moment-equation approach based on the Karhunen–Loeve decomposition developed by Zhang
and Lu (Comput Phys 194(2):773–794, 2004), is applied to unconfined flow with multiple random inputs. The log-transformed hydraulic conductivity F, the recharge R, the Dirichlet boundary condition H, and the Neumann boundary condition Q are assumed to be Gaussian random fields with known means and covariance functions. The F, R, H and Q are first decomposed into finite series in terms of Gaussian standard random variables by the Karhunen–Loeve expansion. The
hydraulic head h is then represented by a perturbation expansion, and each term in the perturbation expansion is written as the products of
unknown coefficients and Gaussian standard random variables obtained from the Karhunen–Loeve expansions. A series of deterministic
partial differential equations are derived from the stochastic partial differential equations. The resulting equations for
uncorrelated and perfectly correlated cases are developed. The equations can be solved sequentially from low to high order
by the finite element method. We examine the accuracy of the KLME approach for the groundwater flow subject to uncorrelated
or perfectly correlated random inputs and study the capability of the KLME method for predicting the head variance in the
presence of various spatially variable parameters. It is shown that the proposed numerical model gives accurate results at
a much smaller computational cost than the Monte Carlo simulation. 相似文献
12.
Existing activity indices (magnetic indices like AE, Kp, Dst or indices based on solar wind parameters) are poor predictors of the instantaneous magnetospheric configuration. We suggest a new activity index - the MT-index. It is defined as the invariant latitude of the isotropic boundary (IB) of \uparrow100 keV protons reduced to the midnight meridian. This IB is a low-altitude signature of the boundary between regions of adiabatic and chaotic regimes of particle motion in the tail current sheet which is controlled by the magnetic field in the equatorial near-Earth tail (at 5-10Re). We have investigated the local time and activity dependence of the IB latitude based on data from about 2000 orbits of NOAA spacecraft. By finding the formula to reduce the IB latitude to midnight meridian, we then evaluate the accuracy of the derived index. We compared the MT-index with the magnetic field measured simultaneously by geosynchronous GOES-2 spacecraft and showed that, unlike the traditional indices, the MT-index displays a good correlation (r\uparrow0.9) with the magnetic field inclination in the nightside portion of the geosynchronous orbit. It is, thus, a good measure to characterize quantitatively the tailward stretching of the tail magnetic field. Based on the measured MT value, a simple numerical procedure is suggested to choose the version of the T89 magnetospheric model. We conclude that the MT-index is the best known predictor of the instantaneous magnetic configuration in the near-Earth magnetotail. It may be available on a regular basis and can be implemented for scientific studies. 相似文献
13.
A modified multiple tension upward infiltration method to estimate the soil hydraulic properties 
下载免费PDF全文
下载免费PDF全文 Determination of saturated hydraulic conductivity, Ks, and the shape parameters α and n of the water retention curve, θ(h), is of paramount importance to characterize the water flow in the vadose zone. This work presents a modified upward infiltration method to estimate Ks, α and n from numerical inverse analysis of the measured cumulative upward infiltration (CUI) at multiple constant tension lower boundary conditions. Using the HYDRUS‐2D software, a theoretical analysis on a synthetic loam soil under different soil tensions (0, 0–10, 0–50 and 0–100 cm), with and without an overpressure step of 10 cm high from the top boundary condition at the end of the upward infiltration process, was performed to check the uniqueness and the accuracy of the solutions. Using a tension sorptivimeter device, the method was validated in a laboratory experiment on five different soils: a coarse and a fine sand, and a 1‐mm sieved loam, clay loam and silt‐gypseous soils. The estimated α and n parameters were compared to the corresponding values measured with the TDR‐pressure cell method. The theoretical analysis demonstrates that Ks and θ(h) can be simultaneously estimated from measured upward cumulative infiltration when high (>50 cm) soil tensions are initially applied at the lower boundary. Alternatively, satisfactory results can be also obtained when medium tensions (<50 cm) and the Ks calculated from the overpressure step at the end of the experiment are considered. A consistent relationship was found between the α (R2 = 0.86, p < 0.02) and n (R2 = 0.97, p < 0.001) values measured with the TDR‐pressure cell and the corresponding values estimated with the tension sorptivimeter. The error between the α (in logarithm scale) and n values estimated with the inverse analysis and the corresponding values measured with pressure chamber were 3.1 and 6.1%, respectively. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
14.
We introduce a simple correction to coastal heads for constant‐density groundwater flow models that contain a coastal boundary, based on previous analytical solutions for interface flow. The results demonstrate that accurate discharge to the sea in confined aquifers can be obtained by direct application of Darcy's law (for constant‐density flow) if the coastal heads are corrected to ((α + 1)/α)hs ? B/2α, in which hs is the mean sea level above the aquifer base, B is the aquifer thickness, and α is the density factor. For unconfined aquifers, the coastal head should be assigned the value . The accuracy of using these corrections is demonstrated by consistency between constant‐density Darcy's solution and variable‐density flow numerical simulations. The errors introduced by adopting two previous approaches (i.e., no correction and using the equivalent fresh water head at the middle position of the aquifer to represent the hydraulic head at the coastal boundary) are evaluated. Sensitivity analysis shows that errors in discharge to the sea could be larger than 100% for typical coastal aquifer parameter ranges. The location of observation wells relative to the toe is a key factor controlling the estimation error, as it determines the relative aquifer length of constant‐density flow relative to variable‐density flow. The coastal head correction method introduced in this study facilitates the rapid and accurate estimation of the fresh water flux from a given hydraulic head measurement and allows for an improved representation of the coastal boundary condition in regional constant‐density groundwater flow models. 相似文献
15.
Magnetic instability in a rapidly rotating cylindrical annulus with a finitely conducting inner core
C. J. Lamb 《地球物理与天体物理流体动力学》2013,107(2-4):227-250
Abstract We consider the stability of a toroidal magnetic field B = B(s*) (where (s*,φ,z*) are cylindrical polar coordinates) in a cylindrical annulus of conducting fluid with inner and outer radii si and s o rotating rapidly about its axis. The outer boundary is taken to be either insulating or perfectly conducting, and the effect of a finite magnetic diffusivity in the inner core is investigated. The ratio of magnetic diffusivity in the inner core to that of the outer core is denoted by ηη→0 corresponding to a perfectly conducting inner core and η→∞ to an insulating one. Comparisons with the results of Fearn (1983b, 1988) were made and a good match found in the limits η→0 and η→∞ with his perfectly conducting and insulating results, respectively. In addition a new mode of instability was found in the eta;→0 regime. Features of this new mode are low frequency (both the frequency and growth rate →0 as η→0) and penetration deep into the inner core. Typically it is unstable at lower magnetic field strengths than the previously known instabilities. 相似文献
16.
Benoit Cushman-Roisin 《地球物理与天体物理流体动力学》2013,107(1-2):61-91
Abstract A new non-linear model of mixing and convection based on a modelling of two buoyant interacting fluids is applied to penetrative convection in the upper ocean due to surface cooling. In view of simple algebra, the model is one-dimensional. Dissipation is included, but no mean shear is present. A non-similar analytical solution is found in the case of a well-mixed layer bounded below by a sharp thermocline treated as a boundary layer. This solution is valid if the Richardson number, R i , defined as the ratio of the total mixed-layer buoyancy to a characteristic rms vertical velocity, is much greater than unity. The model predicts a deepening rate proportional to R i ?3/4. The thermocline remains of constant thickness, and the ratio thermocline thickness to mixed-layer depth decreases as R i ?3/4 as the mixed layer deepens. If the surface flux is constant, the mixed-layer depth increases with time as t ½. The vertical structure throughout the mixed layer and thermocline is given by the analytical solution, and vertical profiles of mean temperature and vertical fluxes are plotted. Computed profiles and available laboratory data agree remarkably well. Moreover, the accuracy of the simple analytical results presented here is comparable to that of sophisticated turbulence numerical models. 相似文献
17.
A new formulation is proposed to model pounding between two adjacent structures, with natural periods T1 and T2 and damping ratios ζ1 and ζ2 under harmonic earthquake excitation, as non‐linear Hertzian impact between two single‐degree‐of‐freedom oscillators. For the case of rigid impacts, a special case of our analytical solution has been given by Davis (‘Pounding of buildings modelled by an impact oscillator’ Earthquake Engineering and Structural Dynamics, 1992; 21 :253–274) for an oscillator pounding on a stationary barrier. Our analytical predictions for rigid impacts agree qualitatively with our numerical simulations for non‐rigid impacts. When the difference in natural periods between the two oscillators increases, the impact velocity also increases drastically. The impact velocity spectrum is, however, relatively insensitive to the standoff distance. The maximum relative impact velocity of the coupled system can occur at an excitation period Tn* which is either between those of the two oscillators or less than both of them, depending on the ratios T1/T2 and ζ1/ζ2. Although the pounding force between two oscillators has been primarily modelled by the Hertz contact law, parametric studies show that the maximum relative impact velocity is not very sensitive to changes in the contact parameters. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
18.
David R. Fearn 《地球物理与天体物理流体动力学》2013,107(1-2):137-162
Abstract A cylindrical annulus containing a conducting fluid and rapidly rotating about its axis is a useful model for the Earth's core. With a shear flow U 0(s)∮, magnetic field B 0(s)∮, and temperature distribution T o(s) (where (s, ∮, z) are cylindrical polar coordinates), many important properties of the core can be modelled while a certain degree of mathematical simplicity is maintained. In the limit of rapid rotation and at geophysically interesting field strengths, the effects of viscous diffusion and fluid inertia are neglected. In this paper, the linear stability of the above basic state to instabilities driven by gradients of B 0 and U 0 is investigated. The global numerical results show both instabilities predicted by a local analysis due to Acheson (1972, 1973, 1984) as well as a new resistive magnetic instability. For the non-diffusive field gradient instability we looked at both monotonic fields [for which the local stability parameter Δ, defined in (1.4), is a constant] and non-monotonic fields (for which Δ is a function of s). For both cases we found excellent qualitative agreement between the numerical and local results but found the local criterion (1.6) for instability to be slightly too stringent. For the non-monotonic fields, instability is confined approximately to the region which is locally unstable. We also investigated the diffusive buoyancy catalysed instability for monotonic fields and found good quantitative agreement between the numerical results and the local condition (1.9). The new resistive instability was found for fields vanishing (or small) at the outer boundary and it is concentrated in the region of that boundary. The resistive boundary layer plays an important part in this instability so it is not of a type which could be predicted using a local stability analysis (which takes no account of the presence of boundaries). 相似文献
19.
Abstract The results of laboratory experiments and numerical model simulations are described in which the motion of a round, negatively-buoyant, turbulent jet discharged horizontally above a slope into a rotating homogeneous fluid has been investigated. For the laboratory study, flow visualisation data are presented to show the complex three-dimensional flow fields generated by the discharge. Analysis of the experimental data indicates that the spatial and temporal developments of the flow field are controlled primarily by the lateral and vertical discharge position of the jet (with respect to the bounding surfaces of the container of width W) and the specific momentum (M 0) and buoyancy (B 0) fluxes driving the jet. The flow is seen to be characterised by the formation of (i) a primary anticyclonic eddy (PCC) close to the source, (ii) an associated secondary cyclonic eddy (SCE) and (iii) a buoyancy-driven bottom boundary current along the right side boundary wall. For the parameter ranges studied, the size L p, s and spatial location x p, s of the PCC and SCE (and the nose velocity u N of the boundary current) are shown to be only weakly-dependent upon the value of the mixed parameter M 0Ω/B 0, where Ω is the background rotation rate. Both L p and x p are shown to scale with the separation distance y?/W of the right side wall (y = 0) from the source (y = y?), both L s and x s scale satisfactorily with the length scale l M (= M 0 3/4/B 0 ½) and u N is determined by the appropriate gravity current speed [(g']0 H]½ and the separation distance y?/W. Numerical model results show good qualitative agreement with the laboratory data with regard to the generation of the PCC, SCE and boundary current as characteristic features of the flow in question. In addition, extension of the numerical model to diagnose potential vorticity and plume thickness distributions for the laboratory cases allow the differences in momentum-and buoyancy-dominated flows to be clearly delineated. Specifically, the characteristic features of the SCE are shown to be strongly dependent upon the value of M 0Ω/B 0 for the buoyant jet flow; not least, the numerical model data are able to confirm the controlling role played by the boundary walls in the laboratory experiments. Quantitative agreement between the numerical and laboratory model data is fair; most significantly, the success of the former model in simulating the dominant flow features from the latter enables the reliable extension of the numerical model to be made to cases of direct oceanic interest. 相似文献
20.
Abstract In the presence of a magnetic field, convection may set in at a stationary or an oscillatory bifurcation, giving rise to branches of steady, standing wave and travelling wave solutions. Numerical experiments provide examples of nonlinear solutions with a variety of different spatiotemporal symmetries, which can be classified by establishing an appropriate group structure. For the idealized problem of two-dimensional convection in a stratified layer the system has left-right spatial symmetry and a continuous symmetry with respect to translations in time. For solutions of period P the latter can be reduced to Z 2 symmetry by sampling solutions at intervals of ½P. Then the fundamental steady solution has the spatiotemporal symmetry D 2 = Z 2 ? Z 2 and symmetry-breaking yields solutions with Z 2 symmetry corresponding to travelling waves, standing waves and pulsating waves. A further loss of symmetry leads to modulated waves. Interactions between the fundamental and its first harmonic are described by the group D 2h = D 2 ? Z 2 and its invariant subgroups, which describe solutions that are either steady or periodic in a uniformly moving frame. For a Boussinesq fluid in a layer with identical top and bottom boundary conditions there is also an up-down symmetry. With fixed lateral boundaries the spatiotemporal symmetries, again described by D 2h and its invariant subgroups, can be related to results obtained in numerical experiments and analysed by Nagata et al. (1990). With periodic boundary conditions, the full symmetry group, D 2h ?Z 2, is of order 16. Its invariant subgroups describe pure and mixed-mode solutions, which may be steady states, standing waves, travelling waves, pulsating waves or modulated waves. 相似文献