首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 62 毫秒
1.
A new computational method for the calculation of shallow water flows with moving physical boundaries is presented. The procedure can cope with shallow water problems having arbitrarily complex geometries and moving boundary elements. Although the method provides a fully boundary-fitted capability, no mesh generation is required in the conventional sense. Solid regions are simply cut out of a background Cartesian mesh with their boundaries represented by different types of cut cell. Moving boundaries are accommodated by up-dating the local cut cell information on a stationary background mesh as the boundaries move. No large-scale re-meshing is required. For the flow calculations, a multi-dimensional high resolution upwind finite volume scheme is used in conjunction with an efficient approximate Riemann solver at flow interfaces, and an exact Riemann solution for a moving piston at moving boundary elements. The method is validated for test problems that include a ship's hull moving at supercritical velocity and two hypothetical landslide events where material plunges laterally into a quiescent shallow lake and a fiord.  相似文献   

2.
A new Lagrangian particle model based on smoothed particle hydrodynamics (SPH) is developed and used to simulate Darcy scale flow and transport in porous media. The method has excellent conservation properties and treats advection exactly. The Lagrangian method is used in stochastic analysis of miscible density-driven fluid flows. Results show that heterogeneity significantly increases dispersion and slows development of Rayleigh–Taylor instability. The presented numerical examples illustrate the advantages of Lagrangian methods for stochastic transport simulations.  相似文献   

3.
A forward particle tracking Eulerian Lagrangian localized adjoint method (ELLAM) is applied to the multicomponent reactive transport problem using a split operator approach. Two split operator algorithms are compared, the Strang algorithm and the sequential non-iterative algorithm (SNIA). The reaction equations are integrated using a coupled predictor corrector algorithm with adaptive time stepping. Reaction time steps are adjusted at the inflow boundary to reflect the actual time of transport inside the solution domain.Results show that split operator ELLAM formulations are competitive with direct or fully coupled ELLAM solutions for reactive transport problems. The SNIA algorithm is more accurate than the Strang splitting algorithm when large time steps are used. The reaction algorithm employed dominates computational effort in runs with large time step sizes. To illustrate the use of the method in practical problems, the model is fitted to aerobic aniline degradation data from laboratory scale column experiments. Model inversion is achieved using non-linear regression with a shuffled complex evolution optimization algorithm and parameter uncertainty is assessed using a Bayesian uncertainty analysis procedure.  相似文献   

4.
The dispersion and deposition of particulate organic matter from a fish cage located in an idealized curved channel with a 90° bend are studied for different horizontal grid resolutions. The model system consists of a three-dimensional, random-walk particle tracking model coupled to a terrain-following ocean model. The particle tracking model is a Lagrangian particle tracking simulator which uses the local flow field, simulated by the ocean model, for advection of the particles and random walk to simulate the turbulent diffusion. The sinking of particles is modeled by imposing an individual particle settling velocity. As the homogeneous water flows through the bend in the channel, the results show that a cross-channel secondary circulation is developed. The motion of this flow is similar to a helical motion where the water in the upper layers moves towards the outer bank and towards the inner bank in the lower layers. The intensity of the secondary circulation will depend on the viscosity scheme and increases as the horizontal grid resolution decreases which significantly affects the distribution of the particles on the seabed. The presence of the secondary circulation leads to that most of the particles that settle, settle close to the inner bank of the channel.  相似文献   

5.
  相似文献   

6.
Numerical solution of large-scale ground water flow and transport problems is often constrained by the convergence behavior of the iterative solvers used to solve the resulting systems of equations. We demonstrate the ability of an algebraic multigrid algorithm (AMG) to efficiently solve the large, sparse systems of equations that result from computational models of ground water flow and transport in large and complex domains. Unlike geometric multigrid methods, this algorithm is applicable to problems in complex flow geometries, such as those encountered in pore-scale modeling of two-phase flow and transport. We integrated AMG into MODFLOW 2000 to compare two- and three-dimensional flow simulations using AMG to simulations using PCG2, a preconditioned conjugate gradient solver that uses the modified incomplete Cholesky preconditioner and is included with MODFLOW 2000. CPU times required for convergence with AMG were up to 140 times faster than those for PCG2. The cost of this increased speed was up to a nine-fold increase in required random access memory (RAM) for the three-dimensional problems and up to a four-fold increase in required RAM for the two-dimensional problems. We also compared two-dimensional numerical simulations of steady-state transport using AMG and the generalized minimum residual method with an incomplete LU-decomposition preconditioner. For these transport simulations, AMG yielded increased speeds of up to 17 times with only a 20% increase in required RAM. The ability of AMG to solve flow and transport problems in large, complex flow systems and its ready availability make it an ideal solver for use in both field-scale and pore-scale modeling.  相似文献   

7.
The Karhunen-Loeve (KL) decomposition and the polynomial chaos (PC) expansion are elegant and efficient tools for uncertainty propagation in porous media. Over recent years, KL/PC-based frameworks have successfully been applied in several contributions for the flow problem in the subsurface context. It was also shown, however, that the accurate solution of the transport problem with KL/PC techniques is more challenging. We propose a framework that utilizes KL/PC in combination with sparse Smolyak quadrature for the flow problem only. In a subsequent step, a Lagrangian sampling technique is used for transport. The flow field samples are calculated based on a PC expansion derived from the solutions at relatively few quadrature points. To increase the computational efficiency of the PC-based flow field sampling, a new reduction method is applied. For advection dominated transport scenarios, where a Lagrangian approach is applicable, the proposed PC/Monte Carlo method (PCMCM) is very efficient and avoids accuracy problems that arise when applying KL/PC techniques to both flow and transport. The applicability of PCMCM is demonstrated for transport simulations in multivariate Gaussian log-conductivity fields that are unconditional and conditional on conductivity measurements.  相似文献   

8.
We present a new numerical approach for simulating geomorphic and stratigraphic processes that combines open‐channel flow with non‐uniform sediment transport law and semi‐empirical diffusive mass wasting. It is designed to facilitate modelling of surface processes across multiple space‐ and time‐scales, and under a variety of environmental and tectonic conditions. The physics of open‐channel flow is primarily based on an adapted Lagrangian formulation of shallow‐water equations. The interaction between flow and surface geology is performed by a non‐uniform total‐load sediment transport law. Additional hillslope processes are simulated using a semi‐empirical method based on a diffusion approach. In the implementation, the resolution of flow dynamics is made on a triangulated grid automatically mapped and adaptively remeshed over a regular orthogonal stratigraphic mesh. These new methods reduce computational time while preserving stability and accuracy of the physical solutions. In order to illustrate the potential of this method for landscape and sedimentary system modelling, we present a set of three generic experiments focusing on assessing the influence of contrasting erodibilities on the evolution of an active bedrock landscape. The modelled ridges morphometrics satisfy established relationships for drainage network geometry and slope distribution, and provide quantitative information on the relative impact of hillslope and channel processes, sediment discharge and alluviation. Our results suggest that contrasting erodibility can stimulate autogenic changes in erosion rate and influence the landscape morphology and preservation. This approach offers new opportunities to investigate joint landscape and sedimentary systems response to external perturbations. The possibility to define and track a large number of materials makes the implementation highly suited to model source‐to‐sink problems where material dispersion is the key question that needs to be addressed, such as natural resources exploration and basin analysis. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

9.
We present a new algorithm for the inversion of full-waveform ground-penetrating radar (GPR) data. It is designed to tame the non-linearity issue that afflicts inverse scattering problems, especially in high contrast media. We first investigate the limitations of current full-waveform time-domain inversion schemes for GPR data and then introduce a much-improved approach based on a combined frequency-time-domain analysis. We show by means of several synthetic tests and theoretical considerations that local minima trapping (common in full bandwidth time-domain inversion) can be avoided by starting the inversion with only the low frequency content of the data. Resolution associated with the high frequencies can then be achieved by progressively expanding to wider bandwidths as the iterations proceed. Although based on a frequency analysis of the data, the new method is entirely implemented by means of a time-domain forward solver, thus combining the benefits of both frequency-domain (low frequency inversion conveys stability and avoids convergence to a local minimum; whereas high frequency inversion conveys resolution) and time-domain methods (simplicity of interpretation and recognition of events; ready availability of FDTD simulation tools).  相似文献   

10.
A Lagrangian particle tracking scheme has been extended to simulate advective transport through coupled surface and subsurface flows. This extended scheme assumes a two-dimensional flow field for the overland domain and a three-dimensional flow field for the subsurface domain. Moreover it is assumed that the flow fields are simulated by a cell centered finite difference method. Pathlines in both the subsurface and the overland domain are simulated by classical particle tracking methods. Exchange of particles between the two domains is simulated by newly-developed algorithms presented in this study. Different algorithms are used depending on the direction of the exchange across the interface between the two domains. In the subsurface domain knowledge about a particle’s pathline is enough to detect a transfer to the surface domain and the solution is straightforward. However, in the two-dimensional overland domain pathlines are parallel to the land surface. Therefore the velocity field in the overland domain alone cannot be used to detect a transfer to the subsurface. We propose a relatively simple algorithm to estimate the probability of transfer to the subsurface domain. It is shown that this algorithm can also be used to handle the transfer from the overland domain to the atmosphere domain. The algorithm to estimate the transfer probabilities is based on the mass balance of water on a streamtube aligned with the particle’s pathline. This newly developed technique ensures that transit time distributions can be simulated accurately. These new relationships are implemented in an existing particle tracking code and are verified using analytical solutions for transit times.  相似文献   

11.
A computational modeling analysis of the flow and sediment transport, and deposition in meandering-river models was performed. The Reynolds stress transport model of the FLUENTTM code was used for evaluating the river flow characteristics, including the mean velocity field and the Reynolds stress components. The simulation results were compared with the available experimental data of the river model and discussed. The Lagrangian tracking of individual particles was performed, and the transport and deposition of particles of various sizes in the meandering river were analyzed. Particular attention was given to the sedimentation patterns of different size particles in the river-bend model. The flow patterns in a physical river were also studied. A Froude number based scale ratio of 1:100 was used, and the flow patterns in the physical and river models are compared. The result shows that the mean-flow quantities exhibit dynamic similarity, but the turbulence parameters of the physical river are different from the model. More strikingly, the particle sedimentation features in the physical and river models do not obey the expected similarity scaling.  相似文献   

12.
The effect of mesh type on the accuracy and computational demands of a two-dimensional Godunov-type flood inundation model is critically examined. Cartesian grids, constrained and unconstrained triangular grids, constrained quadrilateral grids, and mixed meshes are considered, with and without local time stepping (LTS), to determine the approach that maximizes computational efficiency defined as accuracy relative to computational effort. A mixed-mesh numerical scheme is introduced so all grids are processed by the same solver. Analysis focuses on a wide range of dam-break type test cases, where Godunov-type flood models have proven very successful. Results show that different mesh types excel under different circumstances. Cartesian grids are 2–3 times more efficient with relatively simple terrain features such as rectilinear channels that call for a uniform grid resolution, while unstructured grids are about twice as efficient in complex domains with irregular terrain features that call for localized refinements. The superior efficiency of locally refined, unstructured grids in complex terrain is attributable to LTS; the locally refined unstructured grid becomes less efficient using global time stepping. These results point to mesh-type tradeoffs that should be considered in flood modeling applications. A mixed mesh model formulation with LTS is recommended as a general purpose solver because the mesh type can be adapted to maximize computational efficiency.  相似文献   

13.
Smoothed particle hydrodynamics(SPH) is a mesh-free adaptive Lagrangian particle method with attractive features for dealing with the free surface flow.This paper applies the SPH method to simulate the large-amplitude lateral sloshing both with and without a floating body,and the vertical parametrically-excited sloshing in a two-dimensional tank.The numerical results show that the SPH approach has an obvious advantage over conventional mesh-based methods in handling nonlinear sloshing problems such as violent fluid-solid interaction,and flow separation and wave-breaking on the free fluid surface.The SPH method provides a new alternative and an effective way to solve these special strong nonlinear sloshing problems.  相似文献   

14.
Modern ground water characterization and remediation projects routinely require calibration and inverse analysis of large three-dimensional numerical models of complex hydrogeological systems. Hydrogeologic complexity can be prompted by various aquifer characteristics including complicated spatial hydrostratigraphy and aquifer recharge from infiltration through an unsaturated zone. To keep the numerical models computationally efficient, compromises are frequently made in the model development, particularly, about resolution of the computational grid and numerical representation of the governing flow equation. The compromise is required so that the model can be used in calibration, parameter estimation, performance assessment, and analysis of sensitivity and uncertainty in model predictions. However, grid properties and resolution as well as applied computational schemes can have large effects on forward-model predictions and on inverse parameter estimates. We investigate these effects for a series of one- and two-dimensional synthetic cases representing saturated and variably saturated flow problems. We show that "conformable" grids, despite neglecting terms in the numerical formulation, can lead to accurate solutions of problems with complex hydrostratigraphy. Our analysis also demonstrates that, despite slower computer run times and higher memory requirements for a given problem size, the control volume finite-element method showed an advantage over finite-difference techniques in accuracy of parameter estimation for a given grid resolution for most of the test problems.  相似文献   

15.
Due to complex dynamics inherent in the physical models, numerical formulation of subsurface and overland flow coupling can be challenging to solve. ParFlow is a subsurface flow code that utilizes a structured grid discretization in order to benefit from fast and efficient structured solvers. Implicit coupling between subsurface and overland flow modes in ParFlow is obtained by prescribing an overland boundary condition at the top surface of the computational domain. This form of implicit coupling leads to the activation and deactivation of the overland boundary condition during simulations where ponding or drying events occur. This results in a discontinuity in the discrete system that can be challenging to resolve. Furthermore, the coupling relies on unstructured connectivities between the subsurface and surface components of the discrete system, which makes it challenging to use structured solvers to effectively capture the dynamics of the coupled flow. We present a formulation of the discretized algebraic system that enables the use of an analytic form of the Jacobian for the Newton–Krylov solver, while preserving the structured properties of the discretization. An effective multigrid preconditioner is extracted from the analytic Jacobian and used to precondition the Jacobian linear system solver. We compare the performance of the new solver against one that uses a finite difference approximation to the Jacobian within the Newton–Krylov approach, previously used in the literature. Numerical results explores the effectiveness of using the analytic Jacobian for the Newton–Krylov solver, and highlights the performance of the new preconditioner and its cost. The results indicate that the new solver is robust and generally outperforms the solver that is based on the finite difference approximation to the Jacobian, for problems where the overland boundary condition is activated and deactivated during the simulation. A parallel weak scaling study highlights the efficiency of the new solver.  相似文献   

16.
The tightly coupled, strongly nonlinear nature of non-isothermal multi-phase flow in porous media poses a tough challenge for numerical simulation. This trait is even more pronounced, if miscibility is also considered. A primary reason why inclusion of miscibility tends to be problematic are the difficulties stemming from phase transitions: on the one hand, phase transitions need to be included since the presence or absence of fluid phases has a major impact on the flow behavior; on the other hand, convergence of the nonlinear solver may be severely affected if they are not handled robustly.In this work, we present a mathematically sound approach to include phase transitions in the nonlinear system of equations: first, the transition conditions are formulated as a set of local inequality constraints, which are then directly integrated into the nonlinear solver using a nonlinear complementarity function. Under this scheme, Newton-Raphson solvers exhibit considerably more robust convergence behaviour compared to some previous approaches, which is then illustrated by several numerical examples.  相似文献   

17.
Water pollution by industrial and agricultural waste is an increasingly major public health issue. It is therefore important for water engineers and managers to be able to predict accurately the local behaviour of water-borne pollutants. This paper describes the novel and efficient coupling of dynamically adaptive hierarchical grids with standard solvers of the advection–diffusion equation. Adaptive quadtree grids are able to focus on regions of interest such as pollutant fronts, while retaining economy in the total number of grid elements through selective grid refinement. Advection is treated using Lagrangian particle tracking. Diffusion is solved separately using two grid-based methods; one is by explicit finite differences, the other a diffusion-velocity approach. Results are given in two dimensions for pure diffusion of an initially Gaussian plume, advection–diffusion of the Gaussian plume in the rotating flow field of a forced vortex, and the transport of species in a rectangular channel with side wall boundary layers. Close agreement is achieved with analytical solutions of the advection–diffusion equation and simulations from a Lagrangian random walk model. An application to Sepetiba Bay, Brazil is included to demonstrate the method with complex flows and topography.  相似文献   

18.
Numerical modelling of floating debris in the world's oceans   总被引：3，自引：0，他引：3  
A global ocean circulation model is coupled to a Lagrangian particle tracking model to simulate 30 years of input, transport and accumulation of floating debris in the world ocean. Using both terrestrial and maritime inputs, the modelling results clearly show the formation of five accumulation zones in the subtropical latitudes of the major ocean basins. The relative size and concentration of each clearly illustrate the dominance of the accumulation zones in the northern hemisphere, while smaller seas surrounded by densely populated areas are also shown to have a high concentration of floating debris. We also determine the relative contribution of different source regions to the total amount of material in a particular accumulation zone. This study provides a framework for describing the transport, distribution and accumulation of floating marine debris and can be continuously updated and adapted to assess scenarios reflecting changes in the production and disposal of plastic worldwide.  相似文献   

19.
The contaminant transport equation is solved in three dimensions using the Eulerian–Lagrangian Localized Adjoint Method (ELLAM). Trilinear and finite volume test functions defined by the characteristics of the governing equation are employed and compared. Integrations are simplified by forward tracking of integration points along the characteristics. The resulting equations are solved using a preconditioned conjugate gradient method. The algorithm is coupled to a block-centered finite difference approximation of the groundwater flow equation similar to that used in the popular MODFLOW code. The ELLAM is tested by comparison with 1D and 3D analytic solutions. The method is then applied with random, spatially correlated hydraulic conductivities in a simulation of a tracer experiment performed on Cape Cod, Massachusetts. The linear test function ELLAM was found to perform better than the finite volume ELLAM. Both ELLAM formulations were found to be robust, computationally efficient and relatively straightforward to implement. When compared to traditional particle tracking and characteristics codes commonly used with MODFLOW, the ELLAM retains the computational advantages of traditional characteristic methods with the added advantage of good mass conservation.  相似文献   

20.
Recent advancements in analytical solutions to quantify water and solute travel time distributions (TTDs) and the related StorAge Selection (SAS) functions synthesize catchment complexity into a simplified, lumped representation. Although these analytical approaches are efficient in application, they require rarely available long‐term and high‐frequency hydrochemical data for parameter estimation. Alternatively, integrated hydrologic models coupled to Lagrangian particle‐tracking approaches can directly simulate age under different catchment geometries and complexity, but at a greater computational expense. Here, we bridge the two approaches, using a physically based model to explore the uncertainty in the estimation of the SAS function shape. In particular, we study the influence of subsurface heterogeneity, interactions between distinct flow domains (i.e., the vadose zone and saturated groundwater), diversity of flow pathways, and recharge rate on the shape of TTDs and the SAS functions. We use an integrated hydrology model, ParFlow, linked with a particle‐tracking model, SLIM, to compute transient residence times (or ages) at every cell in the domain, facilitating a direct characterization of the SAS function. Steady‐state results reveal that the SAS function shape shows a wide range of variation with respect to the variability in the structure of subsurface heterogeneity. Ensembles of spatially correlated realizations of hydraulic conductivity indicate that the SAS functions in the saturated groundwater have an overall weak tendency toward sampling younger ages, whereas the vadose zone gives a strong preference for older ages. We further show that the influence of recharge rate on the TTD is tightly dependent on the variability of subsurface hydraulic conductivity.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号