Construction of three-phase data to model multiphase flow in porous media: Comparing an optimization approach to the finite element approach

Multiphase flow modelling is a major issue in the assessment of groundwater pollution. Three-phase flows are commonly governed by mathematical models that associate a pressure equation with two saturation equations. These equations involve a number of secondary variables that reflect the fluid behaviour in a porous medium. To improve the computational efficiency of multiphase flow simulators, several simplified reformulations of three-phase flow equations have been proposed. However, they require the construction of new secondary variables adapted to the reformulated flow equations. In this article, two different approaches are compared to quantify these variables. A numerical example is given for a typical fine sand.     

Evidence for enhanced hydration on the northern flank of Olympus Mons, Mars

In this paper we report about a small region on the northern scarp of Olympus Mons showing an increase of the 3 μm hydration band in the OMEGA spectra, together with low superficial temperatures. Although water ice clouds can occurs on the flank of big martian volcanoes, radiative transfer modeling indicates that atmospheric water ice alone cannot justify the shape of the observed band. A fit of the 1.9–3 μm absorption features is obtained by hypothesizing that the study region consists of a mixture of dust and water ice covered by an optically thin (τ=0.08 at 3 μm) layer of dust. Thermal modeling also suggests that water ice in this region may be stable during most of the martian year due to the saturation of the atmosphere. If water ice is responsible for the observed spectral behavior, it might consist of a number of ice or snow patches possibly deposited in small depressions.     

Zur geomagnetischen Tiefensondierung

Zusammenfassung Eine Kurvenvergleichsmethode zur Interpretation räumlicher Induktions-pfeilprofile wird unter der Voraussetzung entwickelt, dass die Verteilung der von den geomagnetischen Variationen im Erdinnern induzierten Ströme linear ist. Sie gründet sich auf die lineare Abhängigkeit der Länge des Induktionspfeils von der vertikalen Komponente des Biot-Savartschen Gesetzes. Die Frequenzabhängigkeit des Induktionspfeiles ermöglicht die Angabe eines für die Leitfähigkeitsstruktur charakteristischen Widerstandes und einer Induktivität. Wird zur Bestimmung von Induktionspfeilen das gesamte Horizontalfeld benutzt, so geht die Linearität der Beziehung zwischen Pfeillänge und Ortsfunktion verloren. Für eine Reihe von Strommodellen werden derartige Ortsfunktionenp _z berechnet und diskutiert. Sie dienen als Beispiele von Musterkurven für die entwickelte Interpretationsmethode.

---

Summary One curve-matching method is developed in order to interprete the geomagnetic induction arrow profiles when the distribution of telluric currents in the underground is assumed to be linear. It bases on the linear dependence of the lengths of the induction arrows from the vertical component of the Biot-Savarts Law, which is a function of the radius vector. This is valid when the inducing field has only a horizontal component and the induced horizontal field is zero. In case that the model approximates the current distribution one gets a resistivity and an inductivity from the dependence of frequency of the induction arrows. Some model curves concerning several current distributions are drawn and discussed.

---

---

Mitteilung aus dem Geomagnetischen Institut der DAW, Potsdam, Nr. 212.     

Die Darstellung des geomagnetischen Potentials zur Epoche 1945,0 durch eine Entwicklung nach Kugelfunktionen bis zur 15. Ordnung

Zusammenfassung Es werden die Ergebnisse einer Berechnung des geomagnetischen Potentials in Form einer Reihenentwicklung nach Kugelfunktionen bis zur 15. Ordnung auf Grund der magnetischen Weltkarten für die Epoche 1945,0 vonE. H. Vestine und andern vorgelegt. Bei der Ableitung des Potentialausdrucks nach der Methode der kleinsten Quadrate wurden die Orthogonalitätseigenschaften der Kugelflächenfunktionen vollständig ausgenutzt.

---

Summary The results of the analysis of the main geomagnetic field for the epoch 1945.0 in a series of spherical harmonics to the 15th degree are presented. The analysis is based on the world magnetic charts derived byVestine et al. The coefficients of the potential are calculated by least-squares approximation taking advantage of the orthogonality of the spherical harmonics over discrete ranges.

     

Strain determination from the measurement of pebble shapes

Cross-sections of pebbles or other near-ellipsoidal bodies in rocks can be measured in three planes, conveniently chosen orthogonal to each other, but otherwise oriented arbitrarily. The three sets of measurements, referred to a single Cartesian coordinate system, can in principle be used to calculate the average pebble shape. Some of these measurements, however, are necessarily redundant and, generally, not compatible without adjustments for the distribution of an error. The calculation of the average ellipsoid is developed, as is that of the strain which transforms an original sphere into such an ellipsoid. This calculation involves the distribution by a weighted least-squares method of the error due to redundancy and leads to the assignment of an estimated error to each of the second-rank tensor components of the strain. The assumption is made that, after being referred to the common coordinate system, the measured values for each tensor component have a normal distribution. This assumption may not always be correct but could be tested before this method is used.