scholarly journals CHALLENGES OF BOUNDARY CONDITIONS IN OPEN OCEAN MODELS

Author(s):  
Judith Herold ◽  
Eric Lemont ◽  
Stuart Bettington ◽  
Edward Couriel

The development and calibration of coastal numerical models to support engineering design and environmental impact studies is a challenging process and one that requires professional judgement and continual assessment of all aspects of the model makeup. Fundamental to the integrity of the model are appropriate boundary conditions and quality observational data for calibration. Open ocean boundary conditions are typically the most complex and important aspect of a model build. They represent the influence of dynamics occurring beyond the model extent, bridging large-scale dynamics to the small-scale processes in the model. This study discusses the challenges of open ocean boundaries and how we utilised data to achieve an effective model.Recorded Presentation from the vICCE (YouTube Link): https://youtu.be/evL7f_17wZg

Energies ◽  
2018 ◽  
Vol 11 (9) ◽  
pp. 2228 ◽  
Author(s):  
Kaushik Biswas

A combined finite element analysis (FEA) and experimental validation approach to estimating effective edge conductivities of vacuum insulation panels (VIPs) embedded in foam-VIP composites is presented. The edge conductivities were estimated by comparing the simulation results with measurements of small-scale (0.61 × 0.61 m) foam-VIP composites and using an error minimization method. The two composites contained multiple VIPs that were butt-jointed with each other in one composite and separated by foam insulation in the other. Edge conductivities were estimated by considering the neighboring materials, i.e., whether the VIPs were adjacent to other VIPs or foam insulation. Models incorporating the edge conductivities were then used to simulate additional small- and large-scale (2.44 × 1.22 m) composites for validation and evaluation of the overall thermal transmission properties. The simulations used either the same boundary conditions as the experiments or used the experimental parameters to define the appropriate boundary conditions.


2000 ◽  
Vol 663 ◽  
Author(s):  
J. Samper ◽  
R. Juncosa ◽  
V. Navarro ◽  
J. Delgado ◽  
L. Montenegro ◽  
...  

ABSTRACTFEBEX (Full-scale Engineered Barrier EXperiment) is a demonstration and research project dealing with the bentonite engineered barrier designed for sealing and containment of waste in a high level radioactive waste repository (HLWR). It includes two main experiments: an situ full-scale test performed at Grimsel (GTS) and a mock-up test operating since February 1997 at CIEMAT facilities in Madrid (Spain) [1,2,3]. One of the objectives of FEBEX is the development and testing of conceptual and numerical models for the thermal, hydrodynamic, and geochemical (THG) processes expected to take place in engineered clay barriers. A significant improvement in coupled THG modeling of the clay barrier has been achieved both in terms of a better understanding of THG processes and more sophisticated THG computer codes. The ability of these models to reproduce the observed THG patterns in a wide range of THG conditions enhances the confidence in their prediction capabilities. Numerical THG models of heating and hydration experiments performed on small-scale lab cells provide excellent results for temperatures, water inflow and final water content in the cells [3]. Calculated concentrations at the end of the experiments reproduce most of the patterns of measured data. In general, the fit of concentrations of dissolved species is better than that of exchanged cations. These models were later used to simulate the evolution of the large-scale experiments (in situ and mock-up). Some thermo-hydrodynamic hypotheses and bentonite parameters were slightly revised during TH calibration of the mock-up test. The results of the reference model reproduce simultaneously the observed water inflows and bentonite temperatures and relative humidities. Although the model is highly sensitive to one-at-a-time variations in model parameters, the possibility of parameter combinations leading to similar fits cannot be precluded. The TH model of the “in situ” test is based on the same bentonite TH parameters and assumptions as for the “mock-up” test. Granite parameters were slightly modified during the calibration process in order to reproduce the observed thermal and hydrodynamic evolution. The reference model captures properly relative humidities and temperatures in the bentonite [3]. It also reproduces the observed spatial distribution of water pressures and temperatures in the granite. Once calibrated the TH aspects of the model, predictions of the THG evolution of both tests were performed. Data from the dismantling of the in situ test, which is planned for the summer of 2001, will provide a unique opportunity to test and validate current THG models of the EBS.


Energies ◽  
2020 ◽  
Vol 13 (24) ◽  
pp. 6667
Author(s):  
Wenjuan Zhang ◽  
Waleed Diab ◽  
Hadi Hajibeygi ◽  
Mohammed Al Kobaisi

Modeling flow and transport in fractured porous media has been a topic of intensive research for a number of energy- and environment-related industries. The presence of multiscale fractures makes it an extremely challenging task to resolve accurately and efficiently the flow dynamics at both the local and global scales. To tackle this challenge, we developed a computational workflow that adopts a two-level hierarchical strategy based on fracture length partitioning. This was achieved by specifying a partition length to split the discrete fracture network (DFN) into small-scale fractures and large-scale fractures. Flow-based numerical upscaling was then employed to homogenize the small-scale fractures and the porous matrix into an equivalent/effective single medium, whereas the large-scale fractures were modeled explicitly. As the effective medium properties can be fully tensorial, the developed hierarchical framework constructed the discrete systems for the explicit fracture–matrix sub-domains using the nonlinear two-point flux approximation (NTPFA) scheme. This led to a significant reduction of grid orientation effects, thus developing a robust, applicable, and field-relevant framework. To assess the efficacy of the proposed hierarchical workflow, several numerical simulations were carried out to systematically analyze the effects of the homogenized explicit cutoff length scale, as well as the fracture length and orientation distributions. The effect of different boundary conditions, namely, the constant pressure drop boundary condition and the linear pressure boundary condition, for the numerical upscaling on the accuracy of the workflow was investigated. The results show that when the partition length is much larger than the characteristic length of the grid block, and when the DFN has a predominant orientation that is often the case in practical simulations, the workflow employing linear pressure boundary conditions for numerical upscaling give closer results to the full-model reference solutions. Our findings shed new light on the development of meaningful computational frameworks for highly fractured, heterogeneous geological media where fractures are present at multiple scales.


2019 ◽  
Vol 865 ◽  
pp. 681-719
Author(s):  
Catherine A. Vreugdenhil ◽  
Bishakhdatta Gayen ◽  
Ross W. Griffiths

Direct numerical simulations are used to investigate the nature of fully resolved small-scale convection and its role in large-scale circulation in a rotating $f$-plane rectangular basin with imposed surface temperature difference. The large-scale circulation has a horizontal geostrophic component and a deep vertical overturning. This paper focuses on convective circulation with no wind stress, and buoyancy forcing sufficiently strong to ensure turbulent convection within the thermal boundary layer (horizontal Rayleigh numbers $Ra\approx 10^{12}{-}10^{13}$). The dynamics are found to depend on the value of a convective Rossby number, $Ro_{\unicode[STIX]{x0394}T}$, which represents the strength of buoyancy forcing relative to Coriolis forces. Vertical convection shifts from a mean endwall plume under weak rotation ($Ro_{\unicode[STIX]{x0394}T}>10^{-1}$) to ‘open ocean’ chimney convection plus mean vertical plumes at the side boundaries under strong rotation ($Ro_{\unicode[STIX]{x0394}T}<10^{-1}$). The overall heat throughput, horizontal gyre transport and zonally integrated overturning transport are then consistent with scaling predictions for flow constrained by thermal wind balance in the thermal boundary layer coupled to vertical advection–diffusion balance in the boundary layer. For small Rossby numbers relevant to circulation in an ocean basin, vertical heat transport from the surface layer into the deep interior occurs mostly in ‘open ocean’ chimney convection while most vertical mass transport is against the side boundaries. Both heat throughput and the mean circulation (in geostrophic gyres, boundary currents and overturning) are reduced by geostrophic constraints.


2020 ◽  
Vol 104 (2) ◽  
pp. 1581-1596
Author(s):  
Thomas Heinze

Abstract Dynamics of snow avalanches or landslides can be described by rapid granular flow. Experimental investigations of granular flow at laboratory scale are often required to analyze flow behaviour and to develop adequate mathematical and numerical models. Most investigations use image-based analysis, and additional sensors such as pressure gauges are not always possible. Testing various scenarios and parameter variations such as different obstacle shapes and positions as well as basal topography and friction usually requires either the construction of a new laboratory setups for each test or a cumbersome reconstruction. In this work, a highly flexible and modular laboratory setup is presented based on LEGO bricks. The flexibility of the model is demonstrated, and possible extensions for future laboratory tests are outlined. The setup is able to reproduce published laboratory experiments addressing current scientific research topics, such as overflow of a rigid reflector, flow on a bumpy surface and against a rigid wall using standard image-based analysis. This makes the setup applicable for quick scenario testing, e.g. for hypothesis testing or for low-cost testing prior to large-scale experiments, and it can contribute to the validation of external results and to benchmarks of numerical models. Small-scale laboratory setups are also very useful for demonstration purposes such as education and public outreach, both crucial in the context of natural hazards. The presented setup enables variation of parameters such as of slope length, channel width, height and shape, inclination, bed friction, obstacle position and shape, as well as density, composition, amount and grain size of flowing mass. Observable quantities are flow type, flow height, flow path and flow velocity, as well as runout distance, size and shape of the deposited material. Additional sensors allow further quantitative assessments, such as local pressure values.


Author(s):  
M. T. Rahmati ◽  
G. Alfano ◽  
H. Bahai

Flexible risers which are used for transporting oil and gas between the seabed and surface in ultra-deep waters have a very complex internal structure. Therefore, accurate modeling of their behaviour is a great challenge for the oil and gas industry. Constitutive laws based on beam models which allow the large-scale dynamics of pipes to be related to the behaviour of its internal components can be used for multi-scale analysis of flexible risers. An integral part of these models is the small-scale FE model chosen and the detailed implementation of the boundary conditions. The small scale FE analyses are typically carried out on models of up to a few meters length. The computational requirements of these methods limit their applications for only multi-scale structural analysis based on a sequential approach. For nested multi-scale approaches (i.e. the so called FE2 method) and for multi-scale multi-physic analyses, e.g. fluid structure interaction modeling of flexible risers, more efficient methods are required. The emphasis of the present work is on a highly efficient small-scale modelling method for flexible risers. By applying periodic boundary conditions, only a small fraction of a flexible pipe is used for detailed analysis. The computational model is firstly described. Then, the capability of the method in capturing the detailed nonlinear effects and the great advantage in terms of significant CPU time saving entailed by this method are demonstrated. For proof of concept the approach is applied on a simplified 3-layer pipe made of inner and outer polymer layers and an intermediate armour layer made of 40 steel tendons.


2019 ◽  
Author(s):  
Ning Li ◽  
Wolfgang Kinzelbach ◽  
Haitao Li ◽  
Wenpeng Li ◽  
Fei Chen

Abstract. To assess the efficiency of groundwater management of an administrative unit, we propose to decompose the groundwater head changes within an administrative unit into inside and outside contributions by using numerical models. Guantao County of Hebei Province, China, serves as an example to demonstrate the decomposition technique. The groundwater flow model of Guantao was constructed using observed heads as prescribed head boundary conditions. The model was coupled with Hydrus 1D, to calculate the groundwater recharge distribution in time reflecting the delay and damping effects of the soil column on seepage at the surface. The model was calibrated by adjusting parameters such as hydraulic conductivities, recharge infiltration ratios and specific yields. The calibrated parameters are then used in a large model with a boundary at large distance from Guantao administrative boundary to determine the groundwater head changes due to inside drivers. The differences of the two models on the Guantao boundary serve as the specified head values on the boundary for a small scale model, which is used to calculate the groundwater head imposed by outside drivers. To eliminate inconsistencies caused by the different types of boundary conditions of large and small models, the groundwater head changes due to inside drivers must be updated. The results indicate that the groundwater head changes in the centre and south of Guantao County are influenced equally by both inside and outside contributions, while in the north outside contributions have the stronger impact. The sensitivity analysis shows that the groundwater head changes and their decomposition are much more sensitive to infiltration ratios than to the aquifer parameters. The parameters within Guantao have a certain influence on the net groundwater head changes while the parameters outside of Guantao have only an influence on the decomposition.


Author(s):  
Matthias Rempel

Sunspots are central to our understanding of solar (and stellar) magnetism in many respects. On the large scale, they link the magnetic field observable in the photosphere to the dynamo processes operating in the solar interior. Properly interpreting the constraints that sunspots impose on the dynamo process requires a detailed understanding of the processes involved in their formation, dynamical evolution and decay. On the small scale, they give an insight into how convective energy transport interacts with the magnetic field over a wide range of field strengths and inclination angles, leading to sunspot fine structure observed in the form of umbral dots and penumbral filaments. Over the past decade, substantial progress has been made on both observational and theoretical sides. Advanced ground- and space-based observations have resolved, for the first time, the details of umbral dots and penumbral filaments and discovered similarities in their substructures. Numerical models have advanced to the degree that simulations of entire sunspots with sufficient resolution to resolve sunspot fine structure are feasible. A combination of improved helioseismic inversion techniques with seismic forward modelling provides new views on the subsurface structure of sunspots. In this review, we summarize recent progress, with particular focus on numerical modelling.


2012 ◽  
Vol 8 (S294) ◽  
pp. 225-236
Author(s):  
M. Hanasz ◽  
D. Woltanski ◽  
K. Kowalik

AbstractWe review recent developments of amplification models of galactic and intergalactic magnetic field. The most popular scenarios involve variety of physical mechanisms, including turbulence generation on a wide range of physical scales, effects of supernovae, buoyancy as well as the magnetorotational instability. Other models rely on galaxy interaction, which generate galactic and intergalactic magnetic fields during galaxy mergers. We present also global galactic-scale numerical models of the Cosmic Ray (CR) driven dynamo, which was originally proposed by Parker (1992). We conduct a series of direct CR+MHD numerical simulations of the dynamics of the interstellar medium (ISM), composed of gas, magnetic fields and CR components. We take into account CRs accelerated in randomly distributed supernova (SN) remnants, and assume that SNe deposit small-scale, randomly oriented, dipolar magnetic fields into the ISM. The amplification timescale of the large-scale magnetic field resulting from the CR-driven dynamo is comparable to the galactic rotation period. The process efficiently converts small-scale magnetic fields of SN-remnants into galactic-scale magnetic fields. The resulting magnetic field structure resembles the X-shaped magnetic fields observed in edge-on galaxies.


Sign in / Sign up

Export Citation Format

Share Document