Towards realistic simulations of lava dome growth using the level set method

The level set method has been implemented in a computational volcanology context. New techniques are presented to solve the advection equation and the reinitialisation equation. These techniques are based upon an algorithm developed in the finite difference context, but are modified to take advantage of the robustness of the finite element method. The resulting algorithm is tested on a well documented Rayleigh–Taylor instability benchmark [19], and on an axisymmetric problem where the analytical solution is known. Finally, the algorithm is applied to a basic study of lava dome growth.

Using the level set method to model endogenous lava dome growth

Modeling volcanic phenomena is complicated by free-surfaces often supporting large rheological gradients. Analytical solutions and analogue models provide explanations for fundamental characteristics of lava flows. But more sophisticated models are needed, incorporating improved physics and rheology to capture realistic events. To advance our understanding of the flow dynamics of highly viscous lava in Peléean lava dome formation, axi-symmetrical Finite Element Method (FEM) models of generic endogenous dome growth have been developed. We use a novel technique, the level-set method, which tracks a moving interface, leaving the mesh unaltered. The model equations are formulated in an Eulerian framework. In this paper we test the quality of this technique in our numerical scheme by considering existing analytical and experimental models of lava dome growth which assume a constant Newtonian viscosity. We then compare our model against analytical solutions for real lava domes extruded on Soufrière, St. Vincent, W.I. in 1979 and Mount St. Helens, USA in October 1980 using an effective viscosity. The level-set method is found to be computationally light and robust enough to model the free-surface of a growing lava dome. Also, by modeling the extruded lava with a constant pressure head this naturally results in a drop in extrusion rate with increasing dome height, which can explain lava dome growth observables more appropriately than when using a fixed extrusion rate. From the modeling point of view, the level-set method will ultimately provide an opportunity to capture more of the physics while benefiting from the numerical robustness of regular grids.

Direct measurement of pulse distortion near the zero-dispersion wavelength in an optical fiber by frequency-resolved optical gating

The technique of frequency-resolved optical gating is used to characterize the intensity and the phase of picosecond pulses after propagation through 700 m of fiber at close to the zero-dispersion wavelength. Using the frequency-resolved optical gating technique, we directly measure the severe temporal distortion resulting from the interplay between self-phase modulation and higher-order dispersion in this regime. The measured intensity and phase of the pulses after propagation are found to be in good agreement with the predictions of numerical simulations with the nonlinear Schrodinger equation. (C) 1997 Optical Society of America.

Ductile fabrics in the zone of active oblique convergence near the Alpine Fault, New Zealand: identifying the neotectonic overprint

The mid-crustal Alpine Schist in central Southern Alps, New Zealand has been exhumed during the past similar to3 m.y. on the hanging wall of the oblique-slip Alpine Fault. These rocks underwent ductile deformation during their passage through the similar to 150-km-wide Pacific-Australia plate boundary zone. Likely to be Cretaceous in age, peak metamorphism predates the largely Pliocene and younger oblique convergence that continues to uplift the Southern Alps today. Late Cenozoic ductile deformation constructively reinforced a pre-existing fabric that was well oriented to accommodate a dextral-transpressive overprint. Quartz microstructures below a recently exhumed brittle-ductile transition zone reflect a late Cenozoic increment of ductile strain that was distributed across deeper levels of the Pacific Plate. Deformation was transpressive, including a dextral-normal shear component that bends and rotates a delaminated panel of Pacific Plate crust onto the oblique footwall ramp of the Alpine Fault. Progressive ductile shear in mylonites at the base of the Pacific Plate overprints earlier fabrics in a dextral-reverse sense, a deformation that accompanies translation of the schists up the Alpine Fault. Ductile shear along that structure affects not only the 12-km-thick section of Alpine mylonites, but is distributed across several kilometres of overlying nonmylonitic rocks. (C) 2001 Elsevier Science Ltd. All rights reserved.

NMR imaging of the diffusion of water at 310 K into semi-IPNs of PEM and Poly(HEMA-co-THFMA)

Magnetic resonance imaging has been used to monitor the diffusion of water at 310 K into a series of semi-IPNs of poly(ethyl methacrylate), PEM, and copolymers of 2-hydroxyethyl methacrylate, HEMA, and tetrahydrofurfuryl methacrylate, THFMA. The diffusion was found to be well described by a Fickian kinetic model in the early stages of the water sorption process, and the diffusion coefficients were found to be slightly smaller than those for the copolymers of HEMA and THFMA, P(HEMA-co-THFMA), containing the same mole fraction of HEMA in the matrix. A second stage sorption process was identified in the later stage of water sorption by the PEM/PTHFMA semi-IPN and for the systems containing a P(HEMA-co-THFMA) component with a mole fraction HEMA of 0.6 or less. This was characterized by the presence of Water near the surface of the cylinders with a longer NMR T-2 relaxation time, which would be characteristic of mobile water, such as water present in large pores or surface fissures. The presence of the drug chlorhexidine in the polymer matrixes at a concentration of 5.625 wt % was found not to modify the properties significantly, but the diffusion coefficients for the water sorption were systematically smaller when the drug was present.

Numerical evaluation of the fields induced by body motion in or near high-field MRI scanners

In modern magnetic resonance imaging, both patients and health care workers are exposed to strong. non-uniform static magnetic fields inside and outside of the scanner. In which body movement may be able to induce electric currents in tissues which could be potentially harmful. This paper presents theoretical investigations into the spatial distribution of induced E-fields in a tissue-equivalent human model when moving at various positions around the magnet. The numerical calculations are based on an efficient. quasi-static, finite-difference scheme. Three-dimensional field profiles from an actively shielded 4 T magnet system are used and the body model projected through the field profile with normalized velocity. The simulation shows that it is possible to induce E-fields/currents near the level of physiological significance under some circumstances and provides insight into the spatial characteristics of the induced fields. The methodology presented herein can be extrapolated to very high field strengths for the evaluation of the effects of motion at a variety of field strengths and velocities. (C) 2004 Elsevier Ltd. All rights reserved.

Intake rates and the functional response in shorebirds (Charadriiformes) eating macro-invertebrates

As field determinations take much effort, it would be useful to be able to predict easily the coefficients describing the functional response of free-living predators, the function relating food intake rate to the abundance of food organisms in the environment. As a means easily to parameterise an individual-based model of shorebird Charadriiformes populations, we attempted this for shorebirds eating macro-invertebrates. Intake rate is measured as the ash-free dry mass (AFDM) per second of active foraging; i.e. excluding time spent on digestive pauses and other activities, such as preening. The present and previous studies show that the general shape of the functional response in shorebirds eating approximately the same size of prey across the full range of prey density is a decelerating rise to a plateau, thus approximating the Holling type 11 ('disc equation') formulation. But field studies confirmed that the asymptote was not set by handling time, as assumed by the disc equation, because only about half the foraging time was spent in successfully or unsuccessfully attacking and handling prey, the rest being devoted to searching. A review of 30 functional responses showed that intake rate in free-living shorebirds varied independently of prey density over a wide range, with the asymptote being reached at very low prey densities (< 150/m(-2)). Accordingly, most of the many studies of shorebird intake rate have probably been conducted at or near the asymptote of the functional response, suggesting that equations that predict intake rate should also predict the asymptote. A multivariate analysis of 468 'spot' estimates of intake rates from 26 shorebirds identified ten variables, representing prey and shorebird characteristics, that accounted for 81 % of the variance in logarithm-transformed intake rate. But four-variables accounted for almost as much (77.3 %), these being bird size, prey size, whether the bird was an oystercatcher Haematopus ostralegus eating mussels Mytilus edulis, or breeding. The four variable equation under-predicted, on average, the observed 30 estimates of the asymptote by 11.6%, but this discrepancy was reduced to 0.2% when two suspect estimates from one early study in the 1960s were removed. The equation therefore predicted the observed asymptote very successfully in 93 % of cases. We conclude that the asymptote can be reliably predicted from just four easily measured variables. Indeed, if the birds are not breeding and are not oystercatchers eating mussels, reliable predictions can be obtained using just two variables, bird and prey sizes. A multivariate analysis of 23 estimates of the half-asymptote constant suggested they were smaller when prey were small but greater when the birds were large, especially in oystercatchers. The resulting equation could be used to predict the half-asymptote constant, but its predictive power has yet to be tested. As well as predicting the asymptote of the functional response, the equations will enable research workers engaged in many areas of shorebird ecology and behaviour to estimate intake rate without the need for conventional time-consuming field studies, including species for which it has not yet proved possible to measure intake rate in the field.

End-to-end modelling of the HEG shock tunnel

We demonstrate an end-to-end computational model of the HEG shock tunnel as a way to extract more precise test flow conditions and as a way of getting predictions of new operating conditions. For a selection of established operating conditions, the L1d program was used to simulate the one-dimensional gas-dynamic processes within the whole of the facility. The program reproduces the compression tube performance reliably and, with the inclusion of a loss factor near the upstream-end of the compression tube, it provides a good estimate of the equilibrium pressure in the shock-reflection region over the set of six standard operating conditions for HEG.

Temperature dependence of the magnetic susceptibility for triangular-lattice antiferromagnets with spatially anisotropic exchange constants

We present the temperature dependence of the uniform susceptibility of spin-half quantum antiferromagnets on spatially anisotropic triangular lattices, using high-temperature series expansions. We consider a model with two exchange constants J1 and J2 on a lattice that interpolates between the limits of a square lattice (J1=0), a triangular lattice (J2=J1), and decoupled linear chains (J2=0). In all cases, the susceptibility, which has a Curie-Weiss behavior at high temperatures, rolls over and begins to decrease below a peak temperature Tp. Scaling the exchange constants to get the same peak temperature shows that the susceptibilities for the square lattice and linear chain limits have similar magnitudes near the peak. Maximum deviation arises near the triangular-lattice limit, where frustration leads to much smaller susceptibility and with a flatter temperature dependence. We compare our results to the inorganic materials Cs2CuCl4 and Cs2CuBr4 and to a number of organic molecular crystals. We find that the former (Cs2CuCl4 and Cs2CuBr4) are weakly frustrated and their exchange parameters determined through the temperature dependence of the susceptibility are in agreement with neutron-scattering measurements. In contrast, the organic materials considered are strongly frustrated with exchange parameters near the isotropic triangular-lattice limit.

On the nonlinear Neumann problem with critical and supercritical nonlinearities

We investigate the solvability of the Neumann problem (1.1) involving a critical Sobolev exponent. In the first part of this work it is assumed that the coeffcients Q and h are at least continuous. Moreover Q is positive on overline Omega and lambda > 0 is a parameter. We examine the common effect of the mean curvature and the shape of the graphs of the coeffcients Q and h on the existence of low energy solutions. In the second part of this work we consider the same problem with Q replaced by - Q. In this case the problem can be supercritical and the existence results depend on integrability conditions on Q and h.

Experimental Analysis of Froude Number Effect on Air Entrainment in the Hydraulic Jump

A hydraulic jump is characterized by strong energy dissipation and mixing, large-scale turbulence, air entrainment, waves and spray. Despite recent pertinent studies, the interaction between air bubbles diffusion and momentum transfer is not completely understood. The objective of this paper is to present experimental results from new measurements performed in rectangular horizontal flume with partially-developed inflow conditions. The vertical distributions of void fraction and air bubbles count rate were recorded for inflow Froude number Fr1 in the range from 5.2 to 14.3. Rapid detrainment process was observed near the jump toe, whereas the structure of the air diffusion layer was clearly observed over longer distances. These new data were compared with previous data generally collected at lower Froude numbers. The comparison demonstrated that, at a fixed distance from the jump toe, the maximum void fraction Cmax increases with the increasing Fr1. The vertical locations of the maximum void fraction and bubble count rate were consistent with previous studies. Finally, an empirical correlation between the upper boundary of the air diffusion layer and the distance from the impingement point was provided.