Numerical simulation of comet nuclei I. Water-ice comets

A one-dimensional simulation of pure water-ice cometary nuclei is presented, and the effect of the nucleus as a heat reservoir is considered. The phase transition from amorphous to crystalline ice is studied for two cases: (1) where the released latent heat goes entirely into heating adjacent layers and (2) where the released latent heat goes entirely into sublimation. For a Halley-like orbit it was found that for case 1 the phase boundary penetrates about 15 m on the first orbit and does not advance until sublimation brings the surface to some 10 m from the phase boundary. For case 2 the phase boundary penetrates about 1 m below the surface and remains at this depth as the surface sublimates. For an orbit like that of Schwassmann-Wachmann 1 the phase boundary penetrates about 50 m initially for case 1 and about 1 m for case 2. There is no further transformation until the entire comet is heated slowly to near the transition temperature, after which the entire nucleus is converted to crystalline ice. For an Encke-type orbit case 1 gives a nearly continuous transition of the entire nucleus to crystalline ice, while for case 2 the initial penetration is about 8 m and remains at this depth relative to the surface as sublimation decreases the cometary radius. Thus the entire comet is converted to crystalline ice just before it is completely dissipated.     

Fast iterative solution of large undrained soil‐structure interaction problems

In view of rapid developments in iterative solvers, it is timely to re‐examine the merits of using mixed formulation for incompressible problems. This paper presents extensive numerical studies to compare the accuracy of undrained solutions resulting from the standard displacement formulation with a penalty term and the two‐field mixed formulation. The standard displacement and two‐field mixed formulations are solved using both direct and iterative approaches to assess if it is cost‐effective to achieve more accurate solutions. Numerical studies of a simple footing problem show that the mixed formulation is able to solve the incompressible problem ‘exactly’, does not create pressure and stress instabilities, and obviate the need for an ad hoc penalty number. In addition, for large‐scale problems where it is not possible to perform direct solutions entirely within available random access memory, it turns out that the larger system of equations from mixed formulation also can be solved much more efficiently than the smaller system of equations arising from standard formulation by using the symmetric quasi‐minimal residual (SQMR) method with the generalized Jacobi (GJ) preconditioner. Iterative solution by SQMR with GJ preconditioning also is more elegant, faster, and more accurate than the popular Uzawa method. Copyright © 2003 John Wiley & Sons, Ltd.     

Numerical simulations of comet nuclei: II. The effect of the dust mantle

The insulating effect of an evolving dust mantle is examined. The role of this mantle in determining the surface temperature of the ice core is studied as a function of the mass fraction of the dust in the ice-dust mixture and the thermal conductivity of the nucleus. Using the so-called “looselattice” model of D.A. Mendis and G.D. Brin (1977, Moon17, 359–372) (which was also extended to include cracks and pores in the mantle), it was found that both high dust to ice ratios and high core conductivities inhibit mantle blowoff. Indeed, it is often possible to build an essentially permanent dust mantle around an ice nucleus, so that the nucleus will take on an asteroidal appearance.     

Transient biouptake flux and accumulation by microorganisms: The case of two types of sites with Langmuir adsorption

The uptake of a chemical species by an aquatic microorganism is modelled considering two kinds of sites where Langmuirian adsorption is followed by first order internalisation kinetics. Simpler models, such as only one internalisation route (while most of the adsorption takes place on non-internalising sites) or a linear isotherm for adsorption on one or both sites, become limiting cases of this double-Langmuirian model. The model considers the sites located on the spherical (or semi-spherical) surface of the organism, and takes diffusion from the medium into explicit account. The numerical solution for the internalisation flux shows a maximum. We provide an estimate for the time needed to reach a certain proximity to steady state. The transient solution confirms that the analytical expressions for the steady-state flux are usually valid and that the accumulated amounts reflect the impact of the short-time uptake. The Instantaneous Steady-State Approximation (ISSA), where an intercept of the linear regression of accumulated amount as a function of time is interpreted as an adsorbed amount, can be critically assessed with the transient numerical code for two cases: (i) when the total burden of metal on the cell is the input data and (ii) when an extraction procedure provides further information on the adsorbed and internalised amount.     

Stochastic Modelling of the Impact of Flood Protection Measures Along the River Waal in the Netherlands

River flooding is a problem of international interest. In the past few years many countries suffered from severe floods. A large part of the Netherlands is below sea level and river levels. The Dutch flood defences along the river Rhine are designed for water levels with a probability of exceedance of 1/1250 per year. These water levels are computed with a hydrodynamic model using a deterministic bed level and a deterministic design discharge. Traditionally, the safety against flooding in the Netherlands is obtained by building and reinforcing dikes. Recently, a new policy was proposed to cope with increasing design discharges in the Rhine and Meuse rivers. This policy is known as the Room for the River (RfR) policy, in which a reduction of flood levels is achieved by measures creating space for the river, such as dike replacement, side channels and floodplain lowering. As compared with dike reinforcement, these measures may have a stronger impact on flow and sediment transport fields, probably leading to stronger morphological effects. As a result of the latter the flood conveyance capacity may decrease over time. An a priori judgement of safety against flooding on the basis of an increased conveyance capacity of the river can be quite misleading. Therefore, the determination of design water levels using a fixed-bed hydrodynamic model may not be justified and the use of a mobile-bed approach may be more appropriate. This problem is addressed in this paper, using a case study of the river Waal (one of the Rhine branches in the Netherlands). The morphological response of the river Waal to a flood protection measure (floodplain lowering in combination with summer levee removal) is analysed. The effect of this measure is subject to various sources of uncertainty. Monte Carlo simulations are applied to calculate the impact of uncertainties in the river discharge on the bed levels. The impact of the “uncertain” morphological response on design flood level predictions is analysed for three phenomena, viz. the impact of the spatial morphological variation over years, the impact of the seasonal morphological variation and the impact of the morphological variability around bifurcation points. The impact of seasonal morphological variations turns out to be negligible, but the other two phenomena appear to have each an appreciable impact (order of magnitude 0.05–0.1 m) on the computed design water levels. We have to note however, that other sources of uncertainty (e.g. uncertainty in hydraulic roughness predictor), which may be of influence, are not taken into consideration. In fact, the present investigation is limited to the sensitivity of the design water levels to uncertainties in the predicted bed level.     

The Dynamics of Small Postglacial Lake Basins and the Nature of Their Archaeological Record: A Case Study of the Middle Palaeolithic Site Neumark‐Nord 2, Germany

Due to favorable conditions of preservation, sedimentary basins provide rich records of human behavior and its environmental context. The conditions for the preservation of archaeological material however vary between basin types (large, river‐fed or small, closed basins), while conditions also differ within a particular basin environment. The goal of this paper is to understand how the dynamics of a small postglacial basin such as Neumark‐Nord 2, a context that dominates the Eemian archaeological record, affected the archaeology situated at its basin margin. The approach used here is to correlate the archaeological record with reconstructions of patterns of deposition and the water conditions within the basin, using lithology, micromorphology, pollen, and macroremains from a transect running from the basin center to the margin. The results show that (1) find levels were exposed to overland flow‐induced winnowing, which vertically concentrated finds but did not cause significant transport, (2) find levels correspond to phases of increased water presence in the basin, and (3) lateral shifts in hominin activity areas may reflect adjustments in the water level. The research shows the importance of large‐scale archaeological excavations and a multidisciplinary sampling strategy that covers both the basin center and the margins, when studying postglacial basin localities like Neumark‐Nord 2.