Three-dimensional aspects of fluid flows in channels. I. Meniscus and thin film regimes

We study the forced displacement of a fluid-fluid interface in a three-dimensional channel formed by two parallel solid plates. Using a lattice-Boltzmann method, we study situations in which a slip velocity arises from diffusion effects near the contact line. The difference between the slip and channel velocities determines whether the interface advances as a meniscus or a thin film of fluid is left adhered to the plates. We find that this effect is controlled by the capillary and Péclet numbers. We estimate the crossover from a meniscus to a thin film and find good agreement with numerical results. The penetration regime is examined in the steady state. We find that the occupation fraction of the advancing finger relative to the channel thickness is controlled by the capillary number and the viscosity contrast between the fluids. For high viscosity contrast, lattice-Boltzmann results agree with previous results. For zero viscosity contrast, we observe remarkably narrow fingers. The shape of the finger is found to be universal.

Front dynamics in turbulent media

A study of a stable front propagating in a turbulent medium is presented. The front is generated through a reaction-diffusion equation, and the turbulent medium is statistically modeled using a Langevin equation. Numerical simulations indicate the presence of two different dynamical regimes. These regimes appear when the turbulent flow either wrinkles a still rather sharp propagating interfase or broadens it. Specific dependences of the propagating velocities on stirring intensities appropriate to each case are found and fitted when possible according to theoretically predicted laws. Different turbulent spectra are considered.

Deep crustal electromagnetic structure of Bhuj earthquake region (India) and its implications

The existence of fluids and partial melt in the lower crust of the seismically active Kutch rift basin (on the western continental margin of India) owing to underplating has been proposed in previous geological and geophysical studies. This hypothesis is examined using magnetotelluric (MT) data acquired at 23 stations along two profiles across Kutch Mainland Uplift and Wagad Uplift. A detailed upper crustal structure is also presented using twodimensional inversion of MT data in the Bhuj earthquake (2001) area. The prominent boundaries of reflection in the upper crust at 5, 10 and 20 km obtained in previous seismic reflection profiles correlate with conductive structures in our models. The MT study reveals 1-2 km thick Mesozoic sediments under the Deccan trap cover. The Deccan trap thickness in this region varies from a few meters to 1.5 km. The basement is shallow on the northern side compared to the south and is in good agreement with geological models as well as drilling information. The models for these profiles indicate that the thickness of sediments would further increase southwards into the Gulf of Kutch. Significant findings of the present study indicate 1) the hypocentre region of the earthquake is devoid of fluids, 2) absence of melt (that is emplaced during rifting as suggested from the passive seismological studies) in the lower crust and 3) a low resistive zone in the depth range of 5-20 km. The present MT study rules out fluidsand melt (magma) as the causative factors that triggered the Bhuj earthquake. The estimated porosity value of 0.02% will explain 100-500 ohm·m resistivity values observed in the lower crust. Based on the seismic velocities and geochemical studies, presence of garnet is inferred. The lower crust consists of basalts - probably generated by partial melting of metasomatised garnet peridotite at deeper depths in the lithosphere - and their composition might be modified by reaction with the spinel peridotites.

Complexation to macromolecules with a large number of sites

This paper presents an approach based on the saddle-point approximation to study the equilibrium interactions between small molecules and macromolecules with a large number of sites. For this case, the application of the Darwin–Fowler method results in very simple expressions for the stoichiometric equilibrium constants and their corresponding free energies in terms of integrals of the binding curve plus a correction term which depends on the first derivatives of the binding curve in the points corresponding to an integer value of the mean occupation number. These expressions are simplified when the number of sites tends to infinity, providing an interpretation of the binding curve in terms of the stoichiometric stability constants. The formalism presented is applied to some simple complexation models, obtaining good values for the free energies involved. When heterogeneous complexation is assumed, simple expressions are obtained to relate the macroscopic description of the binding, given by the stoichiomeric constants, with the microscopic description in terms of the intrinsic stability constants or the affinity spectrum. © 1999 American Institute of Physics.

The measurement of ice velocity, mass balance and thinning-rate on Johnsons Glacier, Livingston Island, South Shetland Islands, Antarctica

A network of twenty stakes was set up on Johnsons Glacier in order to determine its dynamics. During the austral summers from 1994-95 to 1997-98, we estimated surface velocities, mass balances and ice thickness variations. Horizontal velocity increased dow nstream from 1 m a- 1 near the ice divides to 40 m a- 1 near the ice terminus. The accumulation zone showed low accumulation rates (maximum of 0,6 m a- 1 (ice)), whereas in the lower part of the glacier, ablation rates were 4,3 m a- 1 (ice). Over the 3-year study period, both in the accumulation and ablation zones, we detected a reduction in the ice surface level ranging from 2 to 10 m from the annual ve rt ical velocities and ice-thinning data, the mass balance was obtained and compared with the mass balance field values, resulting in similar estimates. Flux values were calculated using cross-section data and horizontal velocities, and compared with the results obtained by means of mass balance and ice thinning data using the continuity equation. The two methods gave similar results.

Surface passivation of crystalline silicon by Cat-CVD amorphous and nanocrystalline thin silicon films

In this work, we study the electronic surface passivation of crystalline silicon with intrinsic thin silicon films deposited by Catalytic CVD. The contactless method used to determine the effective surface recombination velocity was the quasi-steady-state photoconductance technique. Hydrogenated amorphous and nanocrystalline silicon films were evaluated as passivating layers on n- and p-type float zone silicon wafers. The best results were obtained with amorphous silicon films, which allowed effective surface recombination velocities as low as 60 and 130 cms -1 on p- and n-type silicon, respectively. To our knowledge, these are the best results ever reported with intrinsic amorphous silicon films deposited by Catalytic CVD. The passivating properties of nanocrystalline silicon films strongly depended on the deposition conditions, especially on the filament temperature. Samples grown at lower filament temperatures (1600 °C) allowed effective surface recombination velocities of 450 and 600 cms -1 on n- and p-type silicon.

Fluid flow processes at basin scale

S u b s u r face fluid flow plays a significant role in many geologic processes and is increasingly being studied in the scale of sedimentary basins and geologic time perspective. Many economic resources such as petroleum and mineral deposits are products of basin scale fluid flow operating over large periods of time. Such ancient flow systems can be studied through analysis of diagenetic alterations and fluid inclusions to constrain physical and chemical conditions of fluids and rocks during their paleohy d r og e o l ogic evolution. Basin simulation models are useful to complement the paleohy d r og e o l ogic record preserved in the rocks and to derive conceptual models on hydraulic basin evolution and generation of economic resources. Different types of fluid flow regimes may evo l ve during basin evolution. The most important with respect to flow rates and capacity for transport of solutes and thermal energy is gr avitational fluid flow driven by the topographic configuration of a basin. Such flow systems require the basin to be elevated above sea level. Consolidational fluid flow is the principal fluid migration process in basins below sea level, caused by loading of compressible rocks. Flow rates of such systems are several orders of magnitude below topogr a p hy driven flow. Howeve r, consolidation may create significant fluid ove rpressure. Episodic dewatering of ove rpressured compart m e n t s m ay cause sudden fluid release with elevated flow velocities and may cause a transient local thermal and chemical disequilibrium betwe e n fluid and rock. This paper gives an ove rv i ew on subsurface fluid flow processes at basin scale and presents examples related to the Pe n e d è s basin in the central Catalan continental margin including the offshore Barcelona half-graben and the compressive South-Pyrenean basin.

Diagenetic processes in a partially dolomitized carbonate reservoir: Casablanca oil field, Valencia Trough, offshore Spain.

Mesozoic and Neogene carbonates located in the Valencia Trough (offshore Spain, western Mediterranean Sea)are oil reservoirs. This paper investigates the diagenetic evolution of the Upper Jurassic limestones, currently dolomitized, that constitute the main reservoir of the Casablanca oil field. Core samples from Casablanca-1A well have been studied to determine the diagenetic products and their relation with porosity evolution, and to reconstruct the fluid flow history prior to and during oil emplacement. On the basis of petrological observations and geochemical analyses (major, minor and trace element composition and oxygen, carbon and strontium isotope composition), a major dolomitization event is recognized postdating subaerial exposure, erosion and karstification. The dolomitization event originated two replacive dolomites (RD1 and RD2) and two dolomite cements (saddle dolomite cement, SDC, and milky-white dolomite cement, MDC)which are partially cogenetic. RD1, RD2 and SDC precipitated at increasing temperatures (over 60ºC and below 110ºC), probably from meteoric water mixed with marine water. The last dolomite type milky-white dolomite cement) precipitated with increasing burial conditions and by arrival of hydrothermal fluids during the Miocene. The post-dolomitization sequence comprises precipitation of calcite cement and partial calcitization of all previous dolomites. The oxygen, carbon and strontium isotope compositions suggest that this calcite cementation occurred from meteoric waters mixed with Burdigalian - Langhian marine waters trapped in the sediments and expelled by compaction in the moderate to deep burial realm. Normal faults were the conduits for upward migration of these fluids as well as for later oil expulsion from the Burdigalian - Langhian source rocks. Late corrosion associated with organic acid-enriched fluids took place prior or simultaneously to oil migration during the Pliocene, enhancing porosity and increasing eservoir quality.

Reaction time asymmetries between expansion and contraction

Different asymmetries between expansion and contraction (radial motions) have been reported in the literature. Often these patterns have been regarded as implying different channels for each type of radial direction (outward versus inwards) operating at a higher level of visual motion processing. In two experiments (detection and discrimination tasks) we report reaction time asymmetries between expansion and contraction. Power functions were fitted to the data. While an exponent of 0.5 accounted for the expansion data better, a value of unity yielded the best fit for the contraction data. Instead of interpreting these differences as corresponding to different higher order motion detectors, we regard these findings as reflecting the fact that expansion and contraction tap two distinct psychophysical input channels underlying the processing of fast and slow velocities respectively.

Reliable computation of robust response tori on the verge of breakdown

We prove the existence and local uniqueness of invariant tori on the verge of breakdown for two systems: the quasi-periodically driven logistic map and the quasi-periodically forced standard map. These systems exemplify two scenarios: the Heagy-Hammel route for the creation of strange non- chaotic attractors and the nonsmooth bifurcation of saddle invariant tori. Our proofs are computer- assisted and are based on a tailored version of the Newton-Kantorovich theorem. The proofs cannot be performed using classical perturbation theory because the two scenarios are very far from the perturbative regime, and fundamental hypotheses such as reducibility or hyperbolicity either do not hold or are very close to failing. Our proofs are based on a reliable computation of the invariant tori and a careful study of their dynamical properties, leading to the rigorous validation of the numerical results with our novel computational techniques.

Effect of dripline flushing on subsurface drip irrigation systems

The velocity of dripline flushing in subsurface drip irrigation (SDI) systems affects system design, cost, management, performance, and longevity. A 30‐day field study was conducted at Kansas State University to analyze the effect of four targeted flushing velocities (0.23, 0.30, 0.46, and 0.61 m/s) for a fixed 15 min duration of flushing and three flushing frequencies (no flushing or flushing every 15 or 30 days) on SDI emitter discharge and sediments within the dripline and removed in the flushing water. At the end of the field experiment (371 h), the amount of solids carried away by the flushing water and retained in every lateral were determined as well as laboratory determination of emitter discharge for every single emitter within each dripline. Greater dripline flushing velocities, which also resulted in greater flushing volumes, tended to result in greater amounts of solids in the flushing water, but the differences were not always statistically significant. Neither the frequency of flushing nor the interaction of flushing frequency and velocity significantly affected the amount of solids in the flushing water. There was a greater concentration of solids in the beginning one‐third of the 90 m laterals, particularly for treatments with no flushing or with slower dripline flushing velocities. As flushing velocity and concurrently flushing volume increased, there was a tendency for greater solids removal and/or more equal distribution within the dripline. At the end of the field study, the average emitter discharge as measured in the laboratory for a total of 3970 emitters was 0.64 L/h. which was significantly less (approximately 2.5%) than the discharge for new and unused emitters. Only six emitters were nearly or fully clogged, with discharges between 0% and 5% of new and unused emitters. Flushing velocity and flushing frequency did not have consistent significant effects on emitter discharge, and those numerical differences that did exist were small (<3%). Emitter discharge was approximately 3% less for the distal ends of the driplines (last 20% of the dripline). Although not a specific factor in the study, the results of solids removals during flushing and solids retention within the different dripline sections suggest that duration of flushing may be a more cost‐effective management option than increasing the dripline flushing velocity through SDI system design. Finally, although microirrigation system components have been improved over the years, the need for flushing to remove solids and reduce clogging potential has not been eliminated

Massive protostars as gamma-ray sources.

Massive protostars have associated bipolar outflows with velocities of hundreds of km s-1. Such outflows can produce strong shocks when they interact with the ambient medium leading to regions of nonthermal radio emission. Aims. We aim at exploring under which conditions relativistic particles are accelerated at the terminal shocks of the protostellar jets and whether they can produce significant gamma-ray emission. Methods. We estimate the conditions necessary for particle acceleration up to very high energies and gamma-ray production in the nonthermal hot spots of jets associated with massive protostars embedded in dense molecular clouds. Results. We show that relativistic bremsstrahlung and proton-proton collisions can make molecular clouds with massive young stellar objects detectable by the Fermi satellite at MeV-GeV energies and by Cherenkov telescope arrays in the GeV-TeV range. Conclusions. Gamma-ray astronomy can be used to probe the physical conditions in star-forming regions and particle acceleration processes in the complex environment of massive molecular clouds.

3D kinematics of the near-IR HH 223 outflow in L723

In this work, we derive the full 3D kinematics of the near-infrared outflow HH 223, located in the dark cloud Lynds 723 (L723), where a well-defined quadrupolar CO outflow is found. HH 223 appears projected on to the two lobes of the eastwest CO outflow. The radio continuum source VLA 2, towards the centre of the CO outflow, harbours a multiple system of low-mass young stellar objects. One of the components has been proposed to be the exciting source of the eastwest CO outflow. From the analysis of the kinematics, we get further evidence on the relationship between the near-infrared and CO outflows and on the location of their exciting source. The proper motions were derived using multi-epoch, narrow-band H2 (2.122 μm line) images. Radial velocities were derived from the 2.122 μm line of the spectra. Because of the extended (∼5 arcmin), S-shaped morphology of the target, the spectra were obtained with the multi-object-spectroscopy (MOS) observing mode using the instrument Long-Slit Intermediate Resolution Infrared Spectrograph (LIRIS) at the 4.2 m William Herschel Telescope. To our knowledge, this work is the first time that MOS observing mode has been successfully used in the near-infrared range for an extended target.

Present-day deformation of the Pyrenees revealed by GPS surveying and earthquake focal mechanisms until 2011

The Pyrenean mountain range is a slowly deforming belt with continuous and moderate seismic activity. To quantify its deformation field, we present the velocity field estimated from a GPS survey of the Pyrenees spanning 18 yr. The PotSis and ResPyr networks, including a total of 85 GPS sites, were installed and first measured in 1992 and 1995 1997, respectively, and remeasured in 2008 and 2010. We obtain a deformation field with velocities less than 1 mm yr−1 across the range. The estimated velocities for individual stations do not differ significantly from zero with 95 per cent confidence. Even so, we estimate a maximum extensional horizontal strain rate of 2.0 ± 1.7 nanostrain per year in a N S direction in the western part of the range. We do not interpret the vertical displacements due to their large uncertainties. In order to compare the horizontal strain rates with the seismic activity, we analyse a set of 194 focal mechanisms using three methods: (i) the 'r' factor relating their P and T axes, (ii) the stress tensors obtained by fault slip inversion and (iii) the strain-rate tensors. Stress and strain-rate tensors are estimated for: (i) the whole data set, (ii) the eastern and western parts of the range separately, and (iii) eight zones, which are defined based on the seismicity and the tectonic patterns of the Pyrenees. Each of these analyses reveals a lateral variation of the deformation style from compression and extension in the east to extension and strike-slip in the west of the range. Although the horizontal components of the strain-rate tensors estimated from the seismic data are slightly smaller in magnitude than those computed from the GPS velocity field, they are consistent within the 2σ uncertainties. Furthermore, the orientations of their principal axes agree with the mapped active faults.

Full instability behavior of N-dimensional dynamical systems with a one-directional nonlinear vector field

We show how certain N-dimensional dynamical systems are able to exploit the full instability capabilities of their fixed points to do Hopf bifurcations and how such a behavior produces complex time evolutions based on the nonlinear combination of the oscillation modes that emerged from these bifurcations. For really different oscillation frequencies, the evolutions describe robust wave form structures, usually periodic, in which selfsimilarity with respect to both the time scale and system dimension is clearly appreciated. For closer frequencies, the evolution signals usually appear irregular but are still based on the repetition of complex wave form structures. The study is developed by considering vector fields with a scalar-valued nonlinear function of a single variable that is a linear combination of the N dynamical variables. In this case, the linear stability analysis can be used to design N-dimensional systems in which the fixed points of a saddle-node pair experience up to N21 Hopf bifurcations with preselected oscillation frequencies. The secondary processes occurring in the phase region where the variety of limit cycles appear may be rather complex and difficult to characterize, but they produce the nonlinear mixing of oscillation modes with relatively generic features