Numerical solutions for flow in porous media

A numerical approach is proposed to model the flow in porous media using homogenization theory. The proposed concept involves the analyses of micro‐true flow at pore‐level and macro‐seepage flow at macro‐level. Macro‐seepage and microscopic characteristic flow equations are first derived from the Navier–Stokes equation at low Reynolds number through a two‐scale homogenization method. This homogenization method adopts an asymptotic expansion of velocity and pressure through the micro‐structures of porous media. A slightly compressible condition is introduced to express the characteristic flow through only characteristic velocity. This characteristic flow is then numerically solved using a penalty FEM scheme. Reduced integration technique is introduced for the volumetric term to avoid mesh locking. Finally, the numerical model is examined using two sets of permeability test data on clay and one set of permeability test data on sand. The numerical predictions agree well with the experimental data if constraint water film is considered for clay and two‐dimensional cross‐connection effect is included for sand. Copyright © 2003 John Wiley & Sons, Ltd.     

Metallothioneins and trace metals in the dogwhelk Nucella lapillus (L.) collected from Icelandic coasts

Different sizes of the dogwhelk Nucella lapillus were collected from eight locations along the southwest and north coasts of Iceland. Concentrations of total metallothioneins (MTs) and heavy metals (Cd, Cu, Cr, Mn, Ni, and Zn) were analysed using the silver saturation method and inductively coupled plasma-atomic emission spectrometry. The level of tributyltin (TBT) contamination was also assessed using imposex indices, the vas deferens stage index (VDSI) and relative penis size index (RPSI). Gufunes N. lapillus presented the highest values of VDSI (4.0) and RPSI (11.1), followed by Grenivík individuals (VDSI = 3.0; RPSI = 0.9), while the Strandakirkja population showed the lowest VDSI (0.3) and zero RPSI. At a standardised size (0.25 g dry soft-body weight), Grenivík N. lapillus exhibited significantly higher concentrations of all metals whereas overall metal concentrations were significantly lower in individuals from Strandakirkja and Garðskagi compared to other study sites. Partial correlation analyses with size correction indicated that MT concentrations were better correlated with Cd and Cu concentrations than with other metals. At the standard size, the pattern of MT concentrations in N. lapillus from different sites was, however, very different from those of metal profiles. Such discrepancies between the patterns of MT and metals in N. lapillus might be explained by the fact that MT induction could be influenced by various factors such as temperature, dietary metal intake, growth rate and co-existence of other MT-inducing chemicals.     

Toxicity of antifouling biocides to the intertidal harpacticoid copepod Tigriopus japonicus (Crustacea, Copepoda): effects of temperature and salinity

Intertidal harpacticoid copepods are commonly used in eco-toxicity tests worldwide. They predominately live in mid-high shore rock pools and often experience a wide range of temperature and salinity fluctuation. Most eco-toxicity tests are conducted at fixed temperature and salinity and thus the influence of these environmental factors on chemical toxicity is largely unknown. This study investigated the combined effect of temperature and salinity on the acute toxicity of the copepod Tigriopus japonicus against two common biocides, copper (Cu) and tributyltin (TBT) using a 2 × 3 × 4 factorial design (i.e. two temperatures: 25 and 35 °C; three salinities: 15.0‰, 34.5‰ and 45.0‰; three levels of the biocide plus a control). Copper sulphate and tributyltin chloride were used as the test chemicals while distilled water and acetone were utilised as solvents for Cu and TBT respectively. 96h-LC50s of Cu and TBT were 1024 and 0.149 μg l⁻¹ respectively (at 25 °C; 34.5‰) and, based on these results, nominal biocide concentrations of LC0 (i.e. control), LC30, LC50 and LC70 were employed. Analysis of Covariance using ‘concentration’ as the covariate and both ‘temperature’ and ‘salinity’ as fixed factors, showed a significant interaction between temperature and salinity effects for Cu, mortality increasing with temperature but decreasing with elevated salinity. A similar result was revealed for TBT. Both temperature and salinity are, therefore, important factors affecting the results of acute eco-toxicity tests using these marine copepods. We recommend that such eco-toxicity tests should be conducted at a range of environmentally realistic temperature/salinity regimes, as this will enhance the sensitivity of the test and improve the safety margin in line with the precautionary principle.     

On the relation between the amplitude of light and radial velocity variations of β Canis Majoris stars

The relation between the amplitude of light m and the amplitude of the radial velocity 2K for Canis Majoris stars is investigated. A linear relationship between m and 2K is found. However, the two stars BW Vul and Sco, which have the largest radial velocity variations, do not seem to share this relationship.     

Amplitude ratios of several groups of variable stars

The observed positions of classical cepheids, RR Lyrae stars, Scuti stars and dwarf cepheids in the logg-logT _e plane form a continuous sequence, thereby defining the location of maximum instability. The amplitude ratio (the ratio of radial velocity amplitude to light amplitude) is small for variables at the upper end of the instability strip and increases almost linearly towards the lower end of the strip. The theory of radial pulsation predicts the trend of this correlation.     

ENSO anomalies over the Western United States: present and future patterns in regional climate simulations

Surface temperature, precipitation, specific humidity and wind anomalies associated with the warm and cold phases of ENSO simulated by WRF and HadRM are examined for the present and future decades. WRF is driven by ECHAM5 and CCSM3, respectively, and HadRM is driven by HadCM3. For the current decades, all simulations show some capability in resolving the observed warm-dry and cool-wet teleconnection patterns over the PNW and the Southwest U.S. for warm and cold ENSO. Differences in the regional simulations originate primarily from the respective driving fields. For the future decades, the warm-dry and cool-wet teleconnection patterns in association with ENSO are still represented in ECHAM5-WRF and HadRM. However, there are indications of changes in the ENSO teleconnection patterns for CCSM3-WRF in the future, with wet anomalies dominating in the PNW and the Southwest U.S. for both warm and cold ENSO, in contrast to the canonical patterns of precipitation anomalies. Interaction of anomalous wind flow with local terrain plays a critical role in the generation of anomalous precipitation over the western U.S. Anomalous dry conditions are always associated with anomalous airflow that runs parallel to local mountains and wet conditions with airflow that runs perpendicular to local mountains. Future changes in temperature and precipitation associated with the ENSO events in the regional simulations indicate varying responses depending on the variables examined as well as depending on the phase of ENSO.     

Regional climate model projections for the State of Washington

Global climate models do not have sufficient spatial resolution to represent the atmospheric and land surface processes that determine the unique regional climate of the State of Washington. Regional climate models explicitly simulate the interactions between the large-scale weather patterns simulated by a global model and the local terrain. We have performed two 100-year regional climate simulations using the Weather Research and Forecasting (WRF) model developed at the National Center for Atmospheric Research (NCAR). One simulation is forced by the NCAR Community Climate System Model version 3 (CCSM3) and the second is forced by a simulation of the Max Plank Institute, Hamburg, global model (ECHAM5). The mesoscale simulations produce regional changes in snow cover, cloudiness, and circulation patterns associated with interactions between the large-scale climate change and the regional topography and land-water contrasts. These changes substantially alter the temperature and precipitation trends over the region relative to the global model result or statistical downscaling. To illustrate this effect, we analyze the changes from the current climate (1970–1999) to the mid twenty-first century (2030–2059). Changes in seasonal-mean temperature, precipitation, and snowpack are presented. Several climatological indices of extreme daily weather are also presented: precipitation intensity, fraction of precipitation occurring in extreme daily events, heat wave frequency, growing season length, and frequency of warm nights. Despite somewhat different changes in seasonal precipitation and temperature from the two regional simulations, consistent results for changes in snowpack and extreme precipitation are found in both simulations.     

Linear and nonlinear vibration of non-uniform beams on two-parameter foundations using p-elements

A hierarchical finite element is presented for the geometrically nonlinear free and forced vibration of a non-uniform Timoshenko beam resting on a two-parameter foundation. Legendre orthogonal polynomials are used as enriching shape functions to avoid the shear-locking problem. With the enriching degrees of freedom, the accuracy of the computed results and the computational efficiency are greatly improved. The arc-length iterative method is used to solve the nonlinear eigenvalue equation. The computed results of linear and nonlinear vibration analyses show that the convergence of the proposed element is very fast with respect to the number of Legendre orthogonal polynomials used. Since the elastic foundation and the axial load applied at both ends of the beam affect the ratios of linear frequencies associated with the internal resonance, they influence the nonlinear vibration characteristics of the beam. The axial tensile stress of the beam in nonlinear vibration is investigated in this paper, and attention should be paid to the geometrically nonlinear vibration resulting in considerably large axial tensile stress in the beam.     

A Statistical Upscaling Workflow for Warm Solvent Injection Processes for Heterogeneous Heavy Oil Reservoirs

Natural Resources Research - Hybrid solvent- and thermal-based methods offer an important potential for reducing the environmental impact and enhancing recoveries in heavy oil production processes....