scholarly journals Hydrodynamic Simulation of an Irregularly Meandering Gravel-Bed River: Comparison of MIKE 21 FM and Delft3D Flow models

2018 ◽  
Vol 40 ◽  
pp. 02004 ◽  
Author(s):  
Parna Parsapour-moghaddam ◽  
Colin D. Rennie ◽  
Jonathan Slaney
Keyword(s):  
Flow Model ◽  
Numerical Models ◽  
Computational Cost ◽  
Grid Cell ◽  
Flood Management ◽  
Water Levels ◽  
Test Analysis ◽  
Flow Models ◽  
Gravel Bed ◽  
Mike 21

This study aims at hydrodynamic modelling of Bow River, which passes through the City of Calgary, Canada. Bow River has a mobile gravel bed. Erosion and deposition processes were exacerbated by a catastrophic flood in 2013. Channel banks were eroded at various locations, and large gravel bars formed, which could lead to water level changes and accordingly more flooding. This study investigates the performance of Delft3D-Flow and MIKE 21 FM to simulate the hydrodynamics of the river during the 2013 flood. MIKE 21FM employs unstructured triangular mesh while Delft3D-Flow model uses curvilinear structured grids. Performance of each model was evaluated by the available historical water levels. The results of this study demonstrated that, with approximately the same averaged grid resolution, MIKE 21 FM resulted in more accurate results with a higher computational cost compared to the Delft3DFlow model. It was shown that Delft3D-Flow model may require higher grid cell resolution to result in comparably same depth-averaged velocities throughout the study area. However, considering the balance between the computational cost and the accuracy of the results, both models were capable to adequately replicate the hydrodynamics of the river during the 2013 flood. Results of statistical KS and ANOVA test analysis showed that the model predictions were sensitive to the horizontal eddy viscosity and the Manning roughness. This confirms the necessity of an appropriate calibration of the generated numerical models. The findings of this study shed light on the Bow River flood modelling, which can guide flood management.

scholarly journals Coupling ice flow models of varying orders of complexity with the Tiling method

2012 ◽  
Vol 58 (210) ◽  
pp. 776-786 ◽  
Author(s):  
Helene Seroussi ◽  
Hachmi Ben Dhia ◽  
Mathieu Morlighem ◽  
Eric Larour ◽  
Eric Rignot ◽  
...  
Keyword(s):  
Numerical Models ◽  
Computational Cost ◽  
Higher Order ◽  
Ice Flow ◽  
Flow Models ◽  
Continental Scale ◽  
Critical Areas ◽  
Warming Climate ◽  
Computationally Intensive ◽  
A New Technique

AbstractIce flow numerical models are essential for predicting the evolution of ice sheets in a warming climate. Recent research emphasizes the need for higher-order and even full-Stokes flow models, instead of the traditional shallow-ice approximation, whose assumptions are not valid in certain critical areas. These higher-order models are, however, computationally intensive and difficult to use at the continental scale. Here we present a new technique, the Tiling method, to couple ice flow models of varying orders of complexity. The goal of the method is to limit the spatial extent of where higherorder models are applied to reduce the computational cost, while maintaining the model precision. We apply this method on synthetic geometries to demonstrate its practical use. We first use a geometry for which all models yield the same results to check the consistency of the method. Then we apply our method to a geometry for which a full-Stokes model is required in the vicinity of the ice front. Our results show that the hybrid models present significant improvements over mono-model approaches and reduce computational times.

scholarly journals Ensemble-based method with combined fractional flow model for waterflooding optimization

2021 ◽  
Vol 76 ◽  
pp. 7
Author(s):  
Dilayne Santos Oliveira ◽  
Bernardo Horowitz ◽  
Juan Alberto Rojas Tueros
Keyword(s):  
Optimization Problem ◽  
Flow Model ◽  
Computational Cost ◽  
Fractional Flow ◽  
University Of Texas ◽  
Flow Models ◽  
Proxy Models ◽  
The Difference ◽  
The University ◽  
Low Computational Cost

Proxy models are widely used to estimate parameters such as interwell connectivity in the development and management of petroleum fields due to their low computational cost and not require prior knowledge of reservoir properties. In this work, we propose a proxy model to determine both oil and water production to maximize reservoir profitability. The approach uses production history and the Capacitance and Resistance Model based on Producer wells (CRMP), together with the combination of two fractional flow models, Koval [Cao (2014) Development of a Two-phase Flow Coupled Capacitance Resistance Model. PhD Dissertation, The University of Texas at Austin, USA] and Gentil [(2005) The use of Multilinear Regression Models in patterned waterfloods: physical meaning of the regression coefficient. Master’s Thesis, The University of Texas at Austin, USA]. The proposed combined fractional flow model is called Kogen. The combined fractional flow model can be formulated as a constrained nonlinear function fitting. The objective function to be minimized is a measure of the difference between calculated and observed Water cut (Wcut) values or Net Present Values (NPV). The constraint limits the difference in water cuts of the Koval and Gentil models at the time of transition between the two. The problem can be solved using the Sequential Quadratic Programming (SQP) algorithm. The parameters of the CRMP model are the connectivity between wells, time constant and productivity index. These parameters can be found using a Nonlinear Least Squares (NLS) algorithm. With these parameters, it is possible to predict the liquid rate of the wells. The Koval and Gentil models are used to calculate the Wcut in each producer well over the concession period which in turn allows to determine the accumulated oil and water productions. To verify the quality of Kogen model to forecast oil and water productions, we formulated an optimization problem to maximize the reservoir profitability where the objective function is the NPV. The design variables are the injector and producer well controls (liquid rate or bottom hole pressure). In this work the optimization problem is solved using a gradient-based method, SQP. Gradients are approximated using an ensemble-based method. To validate the proposed workflow, we used two realistic reservoirs models, Brush Canyon Outcrop and Brugge field. The results are shown into three stages. In the first stage, we analyze the ensemble size for the gradient computation. Second, we compare the solutions obtained with the three fractional flow models (Koval, Gentil and Kogen) with results achieved directly from the simulator. Third, we use the solutions calculated with the proxy models as starting points for a new high-fidelity optimization process, using exclusively the simulator to calculate the functions involved. This study shows that the proposed combined model, Kogen, consistently generated more accurate results. Also, CRMP/Kogen proxy model has demonstrated its applicability, especially when the available data for model construction is limited, always producing satisfactory results for production forecasting with low computational cost. In addition, it generates a good warm start for high fidelity optimization processes, decreasing the number of simulations by approximately 65%.

Mathematical modelling of salt purification by recrystallization. A cross-current flow model and its comparison with the countercurrent arrangement

1986 ◽  
Vol 51 (11) ◽  
pp. 2489-2501
Author(s):  
Benitto Mayrhofer ◽  
Jana Mayrhoferová ◽  
Lubomír Neužil ◽  
Jaroslav Nývlt
Keyword(s):  
Flow Model ◽  
Current Flow ◽  
Countercurrent Flow ◽  
Flow Models ◽  
Multi Stage ◽  
Stage Crystallization ◽  
The Cross ◽  
Current Flows ◽  
Set Up ◽  
Cross Current

A model is derived for a multi-stage crystallization with cross-current flows of the solution and the crystals being purified. The purity of the product is compared with that achieved in the countercurrent arrangement. A suitable function has been set up which allows the cross-current and countercurrent flow models to be compared and reduces substantially the labour of computation for the countercurrent arrangement. Using the recrystallization of KAl(SO4)2.12 H2O as an example, it is shown that, when the cross-current and countercurrent processes are operated at the same output, the countercurrent arrangement is more advantageous because its solvent consumption is lower.

Study of the Response of PW Traffic Flow Model to a Bottleneck on a Circular Road and its Suitability for Heterogeneous Traffic Conditions

2021 ◽  
Vol 9 (4) ◽  
pp. 460-463
Keyword(s):  
Fluid Flow ◽  
Traffic Flow ◽  
Flow Model ◽  
Equation Model ◽  
Heterogeneous Traffic ◽  
Flow Conditions ◽  
Flow Models ◽  
Traffic Conditions ◽  
Model Response ◽  
Traffic Flow Models

The traffic flow conditions in developing countries are predominantly heterogeneous. The early developed traffic flow models have been derived from fluid flow to capture the behavior of the traffic. The very first two-equation model derived from fluid flow is known as the Payne-Whitham or PW Model. Along with the traffic flow, this model also captures the traffic acceleration. However, the PW model adopts a constant driver behavior which cannot be ignored, especially in the situation of heterogeneous traffic.This research focuses on testing the PW model and its suitability for heterogeneous traffic conditions by observing the model response to a bottleneck on a circular road. The PW model is mathematically approximated using the Roe Decomposition and then the performance of the model is observed using simulations.

scholarly journals A system of conservative regridding for ice–atmosphere coupling in a General Circulation Model (GCM)

2014 ◽  
Vol 7 (3) ◽  
pp. 883-907 ◽  
Author(s):  
R. Fischer ◽  
S. Nowicki ◽  
M. Kelley ◽  
G. A. Schmidt
Keyword(s):  
General Circulation Model ◽  
General Circulation ◽  
Circulation Model ◽  
Grid Cell ◽  
Surface Model ◽  
Ice Flow ◽  
Surface Mass ◽  
Flow Models ◽  
Ice Surface ◽  
Ice Model

Abstract. The method of elevation classes, in which the ice surface model is run at multiple elevations within each grid cell, has proven to be a useful way for a low-resolution atmosphere inside a general circulation model (GCM) to produce high-resolution downscaled surface mass balance fields for use in one-way studies coupling atmospheres and ice flow models. Past uses of elevation classes have failed to conserve mass and energy because the transformation used to regrid to the atmosphere was inconsistent with the transformation used to downscale to the ice model. This would cause problems for two-way coupling. A strategy that resolves this conservation issue has been designed and is presented here. The approach identifies three grids between which data must be regridded and five transformations between those grids required by a typical coupled atmosphere–ice flow model. This paper develops a theoretical framework for the problem and shows how each of these transformations may be achieved in a consistent, conservative manner. These transformations are implemented in Glint2, a library used to couple atmosphere models with ice models. Source code and documentation are available for download. Confounding real-world issues are discussed, including the use of projections for ice modeling, how to handle dynamically changing ice geometry, and modifications required for finite element ice models.

A Study of the Behaviour of a Single-Bladed Waste-Water Pump

2006 ◽  
Vol 220 (2) ◽  
pp. 79-87 ◽  
Author(s):  
J Keays ◽  
C Meskell
Keyword(s):  
Waste Water ◽  
Power Consumption ◽  
Turbulence Model ◽  
Numerical Models ◽  
Sewage Treatment ◽  
Computational Cost ◽  
Primary Source ◽  
Pressure Head ◽  
Primary Objective ◽  
Pump Efficiency

A single-vaned centrifugal pump, typical of the kind employed in waste-water applications (e.g. sewage treatment), has been investigated numerically. The primary objective was to identify a modelling approach that was accurate, but at an acceptable computational cost. A test program has been executed to provide data to validate the numerical models. The global performance of the pump was assessed in terms of the pressure head, the mass flowrate, the power consumption, and the pump efficiency. In addition, time-resolved surface-pressure measurements were made at the volute wall. Five combinations of three modelling approximations (two or 3D; k-ε or Reynolds stress model turbulence model; unsteady or quasi-steady) were investigated and compared with the experimental results. It was found that the choice of turbulence model did not have a significant effect on the predictions. In all cases, the head-discharge curve was well predicted. However, it was found that only the quasi-steady models could capture the trend of the power consumption curve, and hence that of the efficiency. Discrepancies in the magnitude of the power consumption can be accounted for by the lack of losses such as leakage in the numerical models. Qualitative analysis of the numerical results identifies the trailing edge of the impeller as the primary source of power loss, with the flow in the region of the cut water also contributing significantly to the poor overall efficiency of the design.

scholarly journals Assimilation of Lightning Data Using a Nudging Method Involving Low-Level Warming

2014 ◽  
Vol 142 (12) ◽  
pp. 4850-4871 ◽  
Author(s):  
Max R. Marchand ◽  
Henry E. Fuelberg
Keyword(s):  
Numerical Models ◽  
Deep Convection ◽  
Computational Cost ◽  
Wrf Model ◽  
Stage Iv ◽  
Simulated Precipitation ◽  
Synoptic Forcing ◽  
Total Lightning ◽  
Mu Method ◽  
Short Time

Abstract This study presents a new method for assimilating lightning data into numerical models that is suitable at convection-permitting scales. The authors utilized data from the Earth Networks Total Lightning Network at 9-km grid spacing to mimic the resolution of the Geostationary Lightning Mapper (GLM) that will be on the Geostationary Operational Environmental Satellite-R (GOES-R). The assimilation procedure utilizes the numerical Weather Research and Forecasting (WRF) Model. The method (denoted MU) warms the most unstable low levels of the atmosphere at locations where lightning was observed but deep convection was not simulated based on the absence of graupel. Simulation results are compared with those from a control simulation and a simulation employing the lightning assimilation method developed by Fierro et al. (denoted FO) that increases water vapor according to a nudging function that depends on the observed flash rate and simulated graupel mixing ratio. Results are presented for three severe storm days during 2011 and compared with hourly NCEP stage-IV precipitation observations. Compared to control simulations, both the MU and FO assimilation methods produce improved simulated precipitation fields during the assimilation period and a short time afterward based on subjective comparisons and objective statistical scores (~0.1, or 50%, improvement of equitable threat scores). The MU generally performs better at simulating isolated thunderstorms and other weakly forced deep convection, while FO performs better for the case having strong synoptic forcing. Results show that the newly developed MU method is a viable alternative to the FO method, exhibiting utility in producing thunderstorms where observed, and providing improved analyses at low computational cost.

Quantitative investigations of the hydraulic connection between a large reservoir and a buried valley aquifer in southern Alberta

2005 ◽  
Vol 42 (5) ◽  
pp. 1461-1473 ◽  
Author(s):  
B D Smerdon ◽  
C A Mendoza ◽  
A M McCann
Keyword(s):  
Groundwater Flow ◽  
Flow Model ◽  
Hydraulic Head ◽  
Water Levels ◽  
Groundwater Flow Model ◽  
Reservoir Level ◽  
Hydraulic Connection ◽  
Southern Alberta ◽  
Approximate Location ◽  
Pumping Rates

Quantitative investigations, including two aquifer tests and development of a three-dimensional (3D) groundwater flow model, were required to determine the hydraulic connection between an irrigation reservoir and a buried valley aquifer in southern Alberta. Evidence of seepage was detected in the buried valley aquifer 10 km east of the Pine Coulee reservoir at the onset of filling in 1999, when the reservoir level exceeded an elevation of 1035 m above sea level (a.s.l.). Concern for an increase in the local water table and the creation of artesian conditions in the aquifer prompted this study to determine the approximate location of a seepage window that appeared to be connecting the reservoir and aquifer. Observations of hydraulic head in the aquifer during the pumping tests revealed a barrier boundary when the reservoir level was at an elevation of 1035 m a.s.l. and a recharge boundary condition when the elevation exceeded 1039 m a.s.l. These data were used to calibrate a 3D groundwater flow model, which was needed to determine the hydraulic properties and approximate location of the leakage zone. The quantitative investigation showed that seepage likely occurred through the sideslopes of the flooded coulee, rather than through the low-permeability coulee floor sediments or the embankment dam. Further simulations illustrated the expected seepage rates at various reservoir supply levels and the pumping rates required for relief wells installed in the buried valley aquifer to maintain historic aquifer hydraulic head. A brief postanalysis indicated that the forecasted pumping rates were only 15% lower than have been required to maintain preconstruction water levels in the buried valley aquifer.Key words: dams, seepage analysis, groundwater modelling, buried valley aquifer, pumping test.

scholarly journals An objective and efficient method for estimating probabilistic coastal inundation hazards

2019 ◽  
Vol 99 (2) ◽  
pp. 1105-1130 ◽  
Author(s):  
Kun Yang ◽  
Vladimir Paramygin ◽  
Y. Peter Sheng
Keyword(s):  
High Resolution ◽  
Storm Surge ◽  
Computational Cost ◽  
Joint Probability ◽  
Water Levels ◽  
Computational Time ◽  
Probability Method ◽  
Highly Efficient ◽  
Inundation Maps ◽  
Flood Elevation

Abstract The joint probability method (JPM) is the traditional way to determine the base flood elevation due to storm surge, and it usually requires simulation of storm surge response from tens of thousands of synthetic storms. The simulated storm surge is combined with probabilistic storm rates to create flood maps with various return periods. However, the map production requires enormous computational cost if state-of-the-art hydrodynamic models with high-resolution numerical grids are used; hence, optimal sampling (JPM-OS) with a small number of (~ 100–200) optimal (representative) storms is preferred. This paper presents a significantly improved JPM-OS, where a small number of optimal storms are objectively selected, and simulated storm surge responses of tens of thousands of storms are accurately interpolated from those for the optimal storms using a highly efficient kriging surrogate model. This study focuses on Southwest Florida and considers ~ 150 optimal storms that are selected based on simulations using either the low fidelity (with low resolution and simple physics) SLOSH model or the high fidelity (with high resolution and comprehensive physics) CH3D model. Surge responses to the optimal storms are simulated using both SLOSH and CH3D, and the flood elevations are calculated using JPM-OS with highly efficient kriging interpolations. For verification, the probabilistic inundation maps are compared to those obtained by the traditional JPM and variations of JPM-OS that employ different interpolation schemes, and computed probabilistic water levels are compared to those calculated by historical storm methods. The inundation maps obtained with the JPM-OS differ less than 10% from those obtained with JPM for 20,625 storms, with only 4% of the computational time.

scholarly journals rpe v5: An emulator for reduced floating-point precision in large numerical simulations

2016 ◽  
Author(s):  
Andrew Dawson ◽  
Peter Düben
Keyword(s):  
Numerical Simulations ◽  
Climate Models ◽  
Numerical Models ◽  
Computational Cost ◽  
Floating Point ◽  
Climate Modelling ◽  
Acceptable Quality ◽  
Weather And Climate ◽  
Code Changes ◽  
Numerical Precision

Abstract. This paper describes the rpe library which has the capability to emulate the use of arbitrary reduced floating-point precision within large numerical models written in Fortran. The rpe software allows model developers to test how reduced floating-point precision affects the result of their simulations without having to make extensive code changes or port the model onto specialised hardware. The software can be used to identify parts of a program that are problematic for numerical precision and to guide changes to the program to allow a stronger reduction in precision. The development of rpe was motivated by the strong demand for more computing power. If numerical precision can be reduced for an application under consideration while still achieving results of acceptable quality, computational cost can be reduced, since a reduction in numerical precision may allow an increase in performance or a reduction in power consumption. For simulations with weather and climate models, savings due to a reduction in precision could be reinvested to allow model simulations at higher spatial resolution or complexity, or to increase the number of ensemble members to improve predictions. rpe was developed with particular focus on the community of weather and climate modelling, but the software could be used with numerical simulations from other domains.

Sign in / Sign up

Export Citation Format

Share Document