Tidal evolution of hierarchical and inclined systems

We investigate the dynamical evolution of hierarchical three-body systems under the effect of tides, when the ratio of the orbital semi-major axes is small and the mutual inclination is relatively large (greater than 20°). Using the quadrupolar non-restricted approximation for the gravitational interactions and the viscous linear model for tides, we derive the averaged equations of motion in a vectorial formalism which is suitable to model the long-term evolution of a large variety of exoplanetary systems in very eccentric and inclined orbits. In particular, it can be used to derive constraints for stellar spin-orbit misalignment, capture in Cassini states, tidal-Kozai migration, or damping of the mutual inclination. Because our model is valid for the non-restricted problem, it can be used to study systems of identical mass or for the outer restricted problem, such as the evolution of a planet around a binary of stars. Here, we apply our model to various situations in the HD 11964, HD 80606, and HD 98800 systems.     

Data Assimilation (4D-VAR) for Shallow-Water Flow: The Case of the Chicoutimi River

This is a web presentation of the work presented at the 10th Annual Conference of the CFD Society of Canada, “CFD 2002”, at the University of Windsor on June 9-11, 2002. This discussion paper presents the four-dimensional variational data assimilation (4D-VAR) technique as a tool to forecast floods. This discussion will be limited to hydrological forecast. We assume that the weather, here a large rainstorm, had already been forecasted by the meteorological services. In the 4D-VAR technique, we need to minimize, in the sense of Lagrange, a cost function which measures the difference between the forecast and the observations. The physical equations acts as a set of constraints. Here, the model is the shallow-water equations modified to include sediment transport. The minimum was found by using the steepest descent algorithm. This is made possible because the gradient of the cost function can be calculated analytically by using the adjoint equations of the model. To illustrate the 4D-VAR technique, the bypass of a simple theoretical dam as well as the more complex overflowing of the Chicoutimi River at the Chute-Garneau dam (during the 1996 flood) are investigated.     

Phytoliths as indicators of pedogenesis and paleoenvironmental changes in the Brazilian cerrado

This paper presents the preliminary results of phytolith analyses of a peat located in the cerrado of the Uberaba municipality, State of Minas Gerais, Brazil. The phytolith soil sequence is discussed by comparison with phytolith assemblages extracted from dominant plants (Cyperaceae and Poaceae). Increasing mean age of phytolith assemblages with depth is assumed. Poaceae Cuneiform bulliform cell, Parallelepiped bulliform cell and Elongate smooth long cell types dominate up to 80 cm, rapidly drop to 42% at 80 cm and regularly decrease from 42% to 2% upwards. Cyperaceae Rondel concave type shows the inverse trend, being dominant in the upper part of the profile. This pattern can be assigned to increasing selective dissolution of the Cyperaceae phytolith type with depth, or/and to a decrease of water stress suffered by the grasses, leading to a decrease of bulliform cell silicification. Soil processes and paleo-environmental changes hypotheses are discussed.     

Numerical simulation of percussive drilling

This paper presents a novel dynamical model to analyze the long‐term response of a percussive drilling system. This departs from existing approaches that usually consider a single activation and bit/rock interaction cycle for the analysis of the process performance. The proposed model integrates the axial dynamics of an elastic piston and an elastic drill bit, a motion‐dependent pressure law to drive the piston, and a generalized bit/rock interaction law representative of the dynamic indentation taking place at the bit/rock interface. It applies to down‐the‐hole percussive drilling as well as top‐hole, with minor modifications. The model does not account for the angular motion or the hole cleaning, however. The model is first formulated mathematically; then, a finite‐dimensional approximation is proposed for computations. Numerical analyses of the model response, for a low‐size down‐the‐hole percussive system, follow. The period‐1 stationary response for the reference configuration is studied in detail, and parametric analyses assessing the influence on the rate of penetration of the bit/rock interaction parameters, the feed force, and the percussive activation parameters are conducted. These analyses reveal that the multiscale nature of the process is well captured by the model and recover expected trends for the influence of the parameters. They also suggest that a significant increase of the penetration rate can be achieved by increasing the percussive frequency. Copyright © 2015 John Wiley & Sons, Ltd.     

Granitic coastal geomorphology: applying integrated terrestrial and bathymetric LiDAR with multibeam sonar to examine coastal landscape evolution

Coasts composed of resistant lithologies such as granite are generally highly resistant to erosion. They tend to evolve over multiple sea level cycles with highstands acting to remove subaerially weathered material. This often results in a landscape dominated by plunging cliffs with shore platforms rarely occurring. The long‐term evolution of these landforms means that throughout the Quaternary these coasts have been variably exposed to different sea level elevations which means erosion may have been concentrated at different elevations from today. Investigations of the submarine landscape of granitic coasts have however been hindered by an inability to accurately image the nearshore morphology. Only with the advent of multibeam sonar and aerial laser surveying can topographic data now be seamlessly collected from above and below sea level. This study tests the utility of these techniques and finds that very accurate measurements can be made of the nearshore thereby allowing researchers to study the submarine profile with the same accuracy as the subaerial profile. From a combination of terrestrial and marine LiDAR data with multibeam sonar data, it is found that the morphology of granite domes is virtually unaffected by erosion at sea level. It appears that evolution of these landscapes on the coast is a very slow process with modern sea level acting only to remove subaerially weathered debris. The size and orientation of the joints determines the erosional potential of the granite. Where joints are densely spaced (<2 m apart) or the bedrock is highly weathered can semi‐horizontal surfaces form. Copyright © 2014 John Wiley & Sons, Ltd.     

Effects of substratum type on fish assemblages in shallow areas of a tropical estuary

While there is already a comprehensive understanding of the effects of environmental variables, such as dissolved oxygen, temperature and salinity, on the structure, biomass and metabolism of aquatic biota in estuarine habitats, the effect of sedimentation, a harmful underlying factor, remains unclear. The aim of this study was to assess the differences in fish assemblages along the freshwater to salt water gradient of a large tropical estuary, and to evaluate the effects on them of habitat disturbance associated with shallow water sedimentation in the intertidal areas. Fish and environmental variables were recorded in the upper, middle and lower salinity zones of the estuary at ebb tide in both the dry and wet seasons. Sediment samples associated with different levels of habitat disturbance were analysed using granulometry, and their organic content and dissolved oxygen levels were quantified. Water temperature, salinity, pH and dissolved oxygen levels were also measured. Habitat disturbance was found to be correlated with the density, biomass and richness of fish assemblages. A total of 77 species were recorded, forming two distinct fish assemblages, with the Eleotridae family dominating in the upper zone, and Gerreidae, Gobiidae and Tetraodontidae the most common in the middle and lower estuary. Changes in the structure of fish assemblages, including reductions in density, biomass and richness, were associated with disturbance to natural features, where muddy sediment was replaced by sandy sediment and the quantity of organic matter fell. Atherinella brasiliensis was the species which showed a preference for the disturbed areas in the middle and lower zones, while Dormitator maculatus showed a preference for them in the upper estuary. They may be taken as indicators of habitat disturbance due to shallow sedimentation.     

Modelling in-situ upgrading of heavy oil using operator splitting method

The in-situ upgrading (ISU) of bitumen and oil shale is a very challenging process to model numerically because of the large number of components that need to be modelled using a system of equations that are both highly non-linear and strongly coupled. Operator splitting methods are one way of potentially improving computational performance. Each numerical operator in a process is modelled separately, allowing the best solution method to be used for the given numerical operator. A significant drawback to the approach is that decoupling the governing equations introduces an additional source of numerical error, known as the splitting error. The best splitting method for modelling a given process minimises the splitting error whilst improving computational performance compared to a fully implicit approach. Although operator splitting has been widely used for the modelling of reactive-transport problems, it has not yet been applied to the modelling of ISU. One reason is that it is not clear which operator splitting technique to use. Numerous such techniques are described in the literature and each leads to a different splitting error. While this error has been extensively analysed for linear operators for a wide range of methods, the results cannot be extended to general non-linear systems. It is therefore not clear which of these techniques is most appropriate for the modelling of ISU. In this paper, we investigate the application of various operator splitting techniques to the modelling of the ISU of bitumen and oil shale. The techniques were tested on a simplified model of the physical system in which a solid or heavy liquid component is decomposed by pyrolysis into lighter liquid and gas components. The operator splitting techniques examined include the sequential split operator (SSO), the Strang-Marchuk split operator (SMSO) and the iterative split operator (ISO). They were evaluated on various test cases by considering the evolution of the discretization error as a function of the time-step size compared with the results obtained from a fully implicit simulation. We observed that the error was least for a splitting scheme where the thermal conduction was performed first, followed by the chemical reaction step and finally the heat and mass convection operator (SSO-CKA). This method was then applied to a more realistic model of the ISU of bitumen with multiple components, and we were able to obtain a speed-up of between 3 and 5.     

Towards a hierarchical parametrization to address prior uncertainty in ensemble-based data assimilation

Ensemble-based methods are becoming popular assisted history matching techniques with a growing number of field applications. These methods use an ensemble of model realizations, typically constructed by means of geostatistics, to represent the prior uncertainty. The performance of the history matching is very dependent on the quality of the initial ensemble. However, there is a significant level of uncertainty in the parameters used to define the geostatistical model. From a Bayesian viewpoint, the uncertainty in the geostatistical modeling can be represented by a hyper-prior in a hierarchical formulation. This paper presents the first steps towards a general parametrization to address the problem of uncertainty in the prior modeling. The proposed parametrization is inspired in Gaussian mixtures, where the uncertainty in the prior mean and prior covariance is accounted by defining weights for combining multiple Gaussian ensembles, which are estimated during the data assimilation. The parametrization was successfully tested in a simple reservoir problem where the orientation of the major anisotropic direction of the permeability field was unknown.