Numerical evaluation of the damping‐solvent extraction method in the frequency domain

The damping‐solvent extraction method for the analysis of unbounded visco‐elastic media is evaluated numerically in the frequency domain in order to investigate the influence of the computational parameters—domain size, amount of artificial damping, and mesh density—on the accuracy of results. An analytical estimate of this influence is presented, and specific questions regarding the influence of the parameters on the results are answered using the analytical estimate and numerical results for two classical problems: the rigid strip and rigid disc footings on a visco‐elastic half‐space with constant hysteretic material damping. As the domain size is increased, the results become more accurate only at lower frequencies, but are essentially unaffected at higher frequencies. Choosing the domain size to ensure that the static stiffness is computed accurately leads to an unnecessarily large domain for analysis at higher frequencies. The results improve by increasing artificial damping but at a slower rate as the total (material plus artificial) damping ratio ζ_t gets closer to 0.866. However, the results do not deteriorate significantly for the larger amounts of artificial damping, suggesting that ζ_t≈0.6 is appropriate; a larger value is not likely to influence the accuracy of results. Presented results do not support the earlier suggestion that similar accuracy can be achieved by a large bounded domain with small damping or by a small domain with larger damping. Copyright © 2002 John Wiley & Sons, Ltd.     

A response‐based decoupling criterion for multiply‐supported secondary systems

It is often infeasible to carry out coupled analyses of multiply‐supported secondary systems for earthquake excitations. ‘Approximate’ decoupled analyses are then resorted to, unless the response errors due to those are significantly high. This study proposes a decoupling criterion to identify such cases where these errors are likely to be larger than an acceptable level. The proposed criterion is based on the errors in the primary system response due to decoupling and has been obtained by assuming (i) the input excitation to be an ideal white noise process, (ii) cross‐modal correlation to be negligible, and (iii) the combined system to be classically damped. It uses the modal properties of the undamped combined system, and therefore, a perturbation approach has been formulated to determine the combined system properties in case of light to moderately heavy secondary systems. A numerical study has been carried out to illustrate the accuracy achieved with the proposed perturbation formulation. The proposed decoupling criterion has been validated with the help of two example primary‐secondary systems and four example excitation processes. Copyright © 2002 John Wiley & Sons, Ltd.     

Approximate analysis methods for asymmetric plan base‐isolated buildings

An approximate method for linear analysis of asymmetric‐plan, multistorey buildings is specialized for a single‐storey, base‐isolated structure. To find the mode shapes of the torsionally coupled system, the Rayleigh–Ritz procedure is applied using the torsionally uncoupled modes as Ritz vectors. This approach reduces to analysis of two single‐storey systems, each with vibration properties and eccentricities (labelled ‘effective eccentricities’) similar to corresponding properties of the isolation system or the fixed‐base structure. With certain assumptions, the vibration properties of the coupled system can be expressed explicitly in terms of these single‐storey system properties. Three different methods are developed: the first is a direct application of the Rayleigh–Ritz procedure; the second and third use simplifications for the effective eccentricities, assuming a relatively stiff superstructure. The accuracy of these proposed methods and the rigid structure method in determining responses are assessed for a range of system parameters including eccentricity and structure flexibility. For a subset of systems with equal isolation and structural eccentricities, two of the methods are exact and the third is sufficiently accurate; all three are preferred to the rigid structure method. For systems with zero isolation eccentricity, however, all approximate methods considered are inconsistent and should be applied with caution, only to systems with small structural eccentricities or stiff structures. Copyright © 2001 John Wiley & Sons, Ltd.     

Modelling evaporation using an artificial neural network algorithm

This paper investigates the prediction of Class A pan evaporation using the artificial neural network (ANN) technique. The ANN back propagation algorithm has been evaluated for its applicability for predicting evaporation from minimum climatic data. Four combinations of input data were considered and the resulting values of evaporation were analysed and compared with those of existing models. The results from this study suggest that the neural computing technique could be employed successfully in modelling the evaporation process from the available climatic data set. However, an analysis of the residuals from the ANN models developed revealed that the models showed significant error in predictions during the validation, implying loss of generalization properties of ANN models unless trained carefully. The study indicated that evaporation values could be reasonably estimated using temperature data only through the ANN technique. This would be of much use in instances where data availability is limited. Copyright © 2002 John Wiley & Sons, Ltd.     

A double-head photon-counting photometer for fast timing observations

The structure of the electronics of the two-head, high-speed photometer operating on the 152-cm G.D. Cassini telescope of Bologna University has been fully redesigned for high-time-resolution observations of fast variable sources, in particular, of the optical counterparts of X-ray binary systems. The fast photometer now permits sampling times from 0.1 ms to 100 s. The handling of the system clock and I/O configuration have been modified to provide rigorous synchronization of the sampling time with an external high-stability clock and the continuous acquisition of an uninterrupted data stream. The continuous data-acquisition process, which is driven by an external clock, is limited only by the capacity of the hard disk and is supported by a double-buffer SRAM memory designed to overcome any discontinuity in the asynchronous communication with the computer. Real-time monitoring of the data enables continuous evaluation of the weather conditions, and the directly displayed light curve can be used to determine the optimal management of the observations.     

Chlorophyll modulation of mixed layer thermodynamics in a mixed-layer isopycnal General Circulation Model — An example from Arabian Sea and equatorial Pacific

Western tropical Indian Ocean, Arabian Sea, and the equatorial Pacific are known as regions of intense bio-chemical-physical interactions: the Arabian Sea has the largest phytoplankton bloom with seasonal signal, while the equatorial Pacific bloom is perennial with quasi-permanent upwelling. Here, we studied three dimensional ocean thermodynamics comparing recent ocean observation with ocean general circulation model (OPYC) experiment combined with remotely sensed chlorophyll pigment concentrations from the Coastal Zone Color Scanner (CZCS). Using solar radiation parameterization representing observations that a higher abundance of chlorophyll increases absorption of solar irradiance and heating rate in the upper ocean, we showed that the mixed layer thickness decreases more than they would be under clear water conditions. These changes in the model mixed layer were consistent with Joint Global Ocean Flux Study (JGOFS) observations during the 1994-1995 Arabian Sea experiment and epi-fluorescence microscopy (EFM) on samples collected during Equatorial Pacific Ocean Climate Study (EPOCS) in November, 1988. In the Arabian Sea, as the chlorophyll concentrations peak in October (3 mg/m³) after the summer plankton bloom induced by coastal upwelling, the chlorophyll induced biological heating enhanced the sea surface temperature (SST) by as much as 0.6‡C and sub-layer temperature decreases and sub-layer thickness increases. In the equatorial Pacific, modest concentrations of chlorophyll less than 0.3 mg/m³ is enough to introduce a meridional differential heating, which results in reducing the equatorial mixed layer thickness to more than 20 m. The anomalous meridional tilting of the mixed layer bottom enhances off equatorial westward geostrophic currents. Consequently, the equatorial undercurrent transports more water from west to east. We proposed that these numerical model experiments with use of satellite andin situ ocean observations are consistent under three dimensional ocean circulation theory combined with solar radiation transfer process.     

The Mw = 6.0, 7 September 1999 Athens Earthquake

On 7 September 1999 at 11:56 GMT a destructive earthquake (M_w = 6.0) occurred close to Athens (Greece). The rupture process is examined using data from the Cornet local permanent network, as well as teleseismic recordings. Data recorded by a temporary seismological network were analyzed to study the aftershock sequence. The mainshock was relocated at 38.105°N, 23.565°E, about 20 km northwest of Athens. Four foreshocks were also relocated close to the mainshock. The modeling of teleseismic P and SH waves provides a well-constrained focal mechanism of the mainshock (strike = 105°, dip = 55° and rake = -80°) at a depth of 8 km and a seismic moment M₀ = 1.010²⁵ dyn·cm. The obtained fault plane solution represents normal faulting indicating an almost north-south extension. More than 3500 aftershocks were located, 1813 of which present RMS < 0.1 s and ERH, ERZ < 1.0 km. Two main clusters were distinguished, while the depth distribution is concentrated between 2 and 11 km. Over 1000 fault plane solutions of aftershocks were constrained, the majority of which also correspond to N–S extension. No surface breaks were observed but the fault plane solution of the mainshock is in agreement with the tectonics of the area and with the focal mechanisms obtained by aftershocks. The hypocenter of the mainshock is located on the deep western edge of the fault plane. The relocated epicenter coincides with the fringe that represents the highest deformation observed on the differential interferometric image. The calculated source duration is 5 sec, while the estimated dimensions of the fault are 15 km length and 10 km width. The source process is characterized by unilateral eastward rupture propagation, towards the city of Athens. An evident stop phase observed in the recordings of the Cornet local stations is interpreted as a barrier caused by the Aegaleo Mountain.     

An automated strategy for calculation of phase diagram sections and retrieval of rock properties as a function of physical conditions

We formulate an algorithm for the calculation of stable phase relations of a system with constrained bulk composition as a function of its environmental variables. The basis of this algorithm is the approximate representation of the free energy composition surfaces of solution phases by inscribed polyhedra. This representation leads to discretization of high variance phase fields into a continuous mesh of smaller polygonal fields within which the composition and physical properties of the phases are uniquely determined. The resulting phase diagram sections are useful for understanding the phase relations of complex metamorphic systems and for applications in which it is necessary to establish the variations in rock properties such as density, seismic velocities and volatile‐content through a metamorphic cycle. The algorithm has been implemented within a computer program that is general with respect to both the choice of variables and the number of components and phases possible in a system, and is independent of the structure of the equations of state used to describe the phases of the system.     

Up to what temperature is petroleum stable? New insights from a 5200 free radical reactions model

The discovery of crude oils and condensates at ever higher temperatures casts doubt on the validity of the usual geochemical modelling approach, that uses empirical reactions and rate constants. The solution used to account for such a high thermal stability is presently to adjust the rate parameters, but the physical meaning and scientific value of such a strategy can be questioned. We have developed a mechanistic model consisting of 5200 lumped free radical reactions to describe the thermal evolution of a mixture of 52 organic species meant to represent light petroleum. Rate constants used are those available in the literature or estimated using well established thermochemistry-reactivity correlations. Chemical structures included in the model are linear, branched and cyclic hydrocarbons, hydro- and alkyl-aromatics, PAHs, and three heteroatomic compounds. Reactions include cracking and alkylation chains and inhibiting and accelerating reactions from the various reactants. This model has been applied to several mixtures with various proportions of reaction inhibitors and accelerators, and to a composition representing a light mature oil. From the results obtained, we conclude that mature oils will be stable up to 240–260 °C, depending on their composition, and that the thermal cracking of oil to gas is not possible under reasonable basin conditions. The kinetics of petroleum cracking are thus much slower than generally recognized.