Wei Liu, Yingxue Song, Yanping Chen, Gexian Fan, Pengshan Wang, Kai Li
{"title":"Numerical simulation on staggered grids of three-dimensional brinkman-forchheimer flow and heat transfer in porous media","authors":"Wei Liu, Yingxue Song, Yanping Chen, Gexian Fan, Pengshan Wang, Kai Li","doi":"10.1007/s10596-023-10266-7","DOIUrl":"https://doi.org/10.1007/s10596-023-10266-7","url":null,"abstract":"<p>In this paper, three-dimensional numerical algorithm is constructed to simulate the behavior of the Brinkman-Forchheimer flow and thermal fields. Numerical results of velocity, pressure and temperature are obtained by applying the efficient modified two-grid marker and cell (MAC) algorithm on staggered grids with the second-order backward difference formula (BDF2) time approximation. The modified-upwind idea is introduced to convective heat transfer equations for improving accuracy without any numerical oscillation. The second-order convergence rate can be achieved for pressure, velocity and temperature of considered three-dimensional model. Some numerical experiments are presented to illustrate the efficiency of algorithm. The numerical example with analytical solution is used to validate the effectiveness and accuracy of the algorithm by comparing with the results of traditional MAC algorithm. A time-dependent test is proposed to show a detailed sensitivity analysis to indicate the influence of parameters including the <span>(varepsilon )</span>, Forchheimer number, Brinkman number and thermal diffusivity on the physical properties of Brinkman-Forchheimer flow and heat transfer in porous media.</p>","PeriodicalId":10662,"journal":{"name":"Computational Geosciences","volume":"4 1","pages":""},"PeriodicalIF":2.5,"publicationDate":"2024-01-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139376192","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"XFVM modelling of fracture aperture induced by shear and tensile opening.","authors":"Giulia Conti, Stephan Matthäi, Patrick Jenny","doi":"10.1007/s10596-023-10214-5","DOIUrl":"10.1007/s10596-023-10214-5","url":null,"abstract":"<p><p>In reservoir simulation, it is important to understand the mechanical behaviour of fractured rocks and the effect of shear and tensile displacements of fractures on their aperture. Tensile opening directly enhances the fracture aperture, whereas shear of a preexisting rough-walled fracture creates aperture changes dependent on the local stress state. Since fracture dilatation increases reservoir permeability, both processes must be included in a realistic and consistent manner into the mechanical reservoir simulation model. Here, we use the extended finite volume method (XFVM) to conduct flow and geomechanics simulations. In XFVM, fractures are embedded in a poroelastic matrix and are modelled with discontinuous basis functions. On each fracture segment the tractions and compressive forces are calculated, and one extra degree of freedom is added for both the shear and tensile displacement. In this particular XFVM implementation we assume that linear elasticity and steady state fluid pressure adequately constrain the effective stress. In this paper, shear dilation is not calculated a posteriori, but it enters the equations such that aperture changes directly affect the stress state. This is accomplished by adding shear dilation to the displacement gradients and therefore ascertains a consistent representation in the stress-strain relations and force balances. We illustrate and discuss the influence of this extra term in two simple test cases and in a realistic layer-restricted two-dimensional fracture network subjected to plausible in situ stress and pore pressure conditions.</p>","PeriodicalId":10662,"journal":{"name":"Computational Geosciences","volume":"1 1","pages":"227-239"},"PeriodicalIF":2.5,"publicationDate":"2024-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC11080645/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"52164191","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Power law slip boundary condition for Navier-Stokes equations: Discontinuous Galerkin schemes","authors":"J. K. Djoko, V. S. Konlack, T. Sayah","doi":"10.1007/s10596-023-10265-8","DOIUrl":"https://doi.org/10.1007/s10596-023-10265-8","url":null,"abstract":"<p>This study deals with the numerical analysis of several discontinuous Galerkin (DG) methods for the resolution of the Navier-Stokes equations with power law slip boundary condition. The physical context corresponding to this problem is the description of a flow when a position and the direction slip boundary condition is taken into consideration. The main goal in this work is to examine the solvability, convergence of several DG methods, and to discuss their practical resolution by means of fixed point iterative algorithm. Theoretical findings are backed up by solid computational results.</p>","PeriodicalId":10662,"journal":{"name":"Computational Geosciences","volume":"19 1","pages":""},"PeriodicalIF":2.5,"publicationDate":"2023-12-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139070338","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"On spatially correlated observations in importance sampling methods for subsidence estimation","authors":"Samantha S. R. Kim, Femke C. Vossepoel","doi":"10.1007/s10596-023-10264-9","DOIUrl":"https://doi.org/10.1007/s10596-023-10264-9","url":null,"abstract":"<p>The particle filter is a data assimilation method based on importance sampling for state and parameter estimation. We apply a particle filter in two different quasi-static experiments with models of subsidence caused by a compacting reservoir. The first model considers uncorrelated model state variables and observations, with observed subsidence resulting from a single source of strain. In the second model, subsidence is a summation of subsidence contributions from multiple sources which causes spatial dependencies and correlations in the observed subsidence field. Assimilating these correlated subsidence fields may trigger weight collapse. With synthetic tests, we show in a model of subsidence with 50 independent state variables and spatially correlated subsidence a minimum of <span>(varvec{10^{13}})</span> particles are required to have information in the posterior distribution identical to that in a model with 50 independent and spatially uncorrelated observations. Spatial correlations cause an information loss which can be quantified with mutual information. We illustrate how a stronger spatial correlation results in lower information content in the posterior and we empirically derive the required ensemble size for the importance sampling to remain effective. We furthermore illustrate how this loss of information is reflected in the log likelihood, and how this depends on the number of model state variables. Based on these empirical results, we propose criteria to evaluate the required ensemble size in data assimilation of spatially correlated observation fields.</p>","PeriodicalId":10662,"journal":{"name":"Computational Geosciences","volume":"38 1","pages":""},"PeriodicalIF":2.5,"publicationDate":"2023-12-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138520456","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Parameterizing the fluid forces on limpet shells in unidirectional flow","authors":"Carley Walker, Julian Simeonov, Ian Adams","doi":"10.1007/s10596-023-10263-w","DOIUrl":"https://doi.org/10.1007/s10596-023-10263-w","url":null,"abstract":"<p>Current parameterizations of the hydrodynamic forces on irregular particles consider some shape dependencies, but lack an explicit dependence on the orientation with respect to the flow. In this paper, we propose a new parameterization of the drag and lift forces acting on whole Limpet shells at arbitrary orientations with respect to the direction of flow through the linear regression of fluid forces against the velocity components in an object frame of reference. The fluid forces were estimated using boundary layer-resolving Reynolds-averaged Navier-Stokes (RANS) simulations. We verified the accuracy of the shear stress transport (SST) <span>(k-omega )</span> turbulence model on flat plates with varying angles of attack, and we achieved coefficients of determination versus existing data of approximately 0.95 for both the drag and lift coefficients. From the linear regression of our simulated force data, we developed a model as a function of 3-dimensional orientations to predict the hydrodynamic forces acting on a Limpet shell with coefficients of determination of 0.80 for normal forces and 0.51 for longitudinal forces.</p>","PeriodicalId":10662,"journal":{"name":"Computational Geosciences","volume":"28 7","pages":""},"PeriodicalIF":2.5,"publicationDate":"2023-11-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138520457","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Gillien Latour, Pierre Horgue, François Renard, Romain Guibert, Gérald Debenest
{"title":"Comparative calibration of 1D+2D and 3D hydrogeological watershed models","authors":"Gillien Latour, Pierre Horgue, François Renard, Romain Guibert, Gérald Debenest","doi":"10.1007/s10596-023-10261-y","DOIUrl":"https://doi.org/10.1007/s10596-023-10261-y","url":null,"abstract":"<p>In this work, we study the calibration of the parameters of a hydrogeological watershed model by comparing a 1D+2D approach that combines unsaturated 1D columns and a saturated 2D model, with a full 3D approach. In a first step, a heterogeneous permeability field is estimated by an inversion procedure for each model (2D saturated and 3D unsaturated). The fields obtained are similar but the calculation time is obviously much higher in the case of the 3D model: the 2D model seems therefore sufficient and more efficient to evaluate permeability fields using piezometric measurements in the case of vertically homogeneous aquifers. The second step focuses on the calibration of the hydraulic parameters by adjusting the hydraulic heights either derived from a 1D+2D reference model at several fictitious points distributed over the entire domain, or measured in a dozen real piezometers. Both approaches provide a good fit to the piezometric measurements, but the parameter values differ significantly: the van Genuchten alpha coefficient is unrealistic in the 1D+2D approach, reflecting a poorer consideration of the modeling unsaturated zone, while the porosity value is higher in the 3D approach, which can probably be remedied by developing a more suitable cost function.</p>","PeriodicalId":10662,"journal":{"name":"Computational Geosciences","volume":"11 1","pages":""},"PeriodicalIF":2.5,"publicationDate":"2023-11-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138542803","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Dede Tarwidi, Sri Redjeki Pudjaprasetya, Didit Adytia
{"title":"The effects of dispersion and non-linearity on the simulation of landslide-generated waves using the reduced two-layer non-hydrostatic model","authors":"Dede Tarwidi, Sri Redjeki Pudjaprasetya, Didit Adytia","doi":"10.1007/s10596-023-10262-x","DOIUrl":"https://doi.org/10.1007/s10596-023-10262-x","url":null,"abstract":"<p>This paper revisits the previously developed NH-2LR (reduced two-layer non-hydrostatic) model. The governing equations and numerical schemes are written in terms of normalized variables, with two dimensionless parameters representing dispersion and non-linearity. By utilizing analytical solutions and laboratory experiments, this study aims to validate the numerical NH-2LR model and investigate the effects of dispersion and non-linearity on the resulting waves. The first validation employs the analytical solution of the linear and fully dispersive model of a landslide moving with constant velocity on a flat bottom. The second validation involves a landslide hump sliding over a constant beach slope. A closer look at the run-up height reveals that this case is non-dispersive. Furthermore, we found that the dispersion effect was evident from the beginning of the wave formation process. Finally, we compare our numerical results to experiments on submarine landslides on sloping beaches. We found that dispersion is essential in the early generation and propagation of waves in off-shore regions. Moreover, non-linearity significantly influences the maximum run-up of landslide-generated waves.</p>","PeriodicalId":10662,"journal":{"name":"Computational Geosciences","volume":"130 11","pages":""},"PeriodicalIF":2.5,"publicationDate":"2023-11-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138520443","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Multi-asset closed-loop reservoir management using deep reinforcement learning","authors":"Yusuf Nasir, Louis J. Durlofsky","doi":"10.1007/s10596-023-10255-w","DOIUrl":"https://doi.org/10.1007/s10596-023-10255-w","url":null,"abstract":"","PeriodicalId":10662,"journal":{"name":"Computational Geosciences","volume":"10 8","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-11-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135821136","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Numerical simulation of multiscale fault systems with rate- and state-dependent friction","authors":"Carsten Gräser, Ralf Kornhuber, Joscha Podlesny","doi":"10.1007/s10596-023-10231-4","DOIUrl":"https://doi.org/10.1007/s10596-023-10231-4","url":null,"abstract":"Abstract We consider the deformation of a geological structure with non-intersecting faults that can be represented by a layered system of viscoelastic bodies satisfying rate- and state-depending friction conditions along the common interfaces. We derive a mathematical model that contains classical Dieterich- and Ruina-type friction as special cases and accounts for possibly large tangential displacements. Semi-discretization in time by a Newmark scheme leads to a coupled system of nonsmooth, convex minimization problems for rate and state to be solved in each time step. Additional spatial discretization by a mortar method and piecewise constant finite elements allows for the decoupling of rate and state by a fixed point iteration and efficient algebraic solution of the rate problem by truncated nonsmooth Newton methods. Numerical experiments with a spring slider and a layered multiscale system illustrate the behavior of our model as well as the efficiency and reliability of the numerical solver.","PeriodicalId":10662,"journal":{"name":"Computational Geosciences","volume":"10 5","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-10-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135765679","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A 3D organized point cloud clustering algorithm for seismic fault data based on region growth","authors":"Lihong Zhao, Minghao Cai, Renwei Ding, Yujie Zhang, Shuo Zhao, Jinwei Zhang, Jing Yang","doi":"10.1007/s10596-023-10259-6","DOIUrl":"https://doi.org/10.1007/s10596-023-10259-6","url":null,"abstract":"","PeriodicalId":10662,"journal":{"name":"Computational Geosciences","volume":"21 6","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-10-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"134909119","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}