Finite element modeling of fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins

In order to use the finite element method for solving fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins effectively and efficiently, we have presented, in this paper, the new concept and numerical algorithms to deal with the fundamental issues associated with the fluid-rock interaction problems. These fundamental issues are often overlooked by some purely numerical modelers. (1) Since the fluid-rock interaction problem involves heterogeneous chemical reactions between reactive aqueous chemical species in the pore-fluid and solid minerals in the rock masses, it is necessary to develop the new concept of the generalized concentration of a solid mineral, so that two types of reactive mass transport equations, namely, the conventional mass transport equation for the aqueous chemical species in the pore-fluid and the degenerated mass transport equation for the solid minerals in the rock mass, can be solved simultaneously in computation. (2) Since the reaction area between the pore-fluid and mineral surfaces is basically a function of the generalized concentration of the solid mineral, there is a definite need to appropriately consider the dependence of the dissolution rate of a dissolving mineral on its generalized concentration in the numerical analysis. (3) Considering the direct consequence of the porosity evolution with time in the transient analysis of fluid-rock interaction problems; we have proposed the term splitting algorithm and the concept of the equivalent source/sink terms in mass transport equations so that the problem of variable mesh Peclet number and Courant number has been successfully converted into the problem of constant mesh Peclet and Courant numbers. The numerical results from an application example have demonstrated the usefulness of the proposed concepts and the robustness of the proposed numerical algorithms in dealing with fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins. (C) 2001 Elsevier Science B.V. All rights reserved.

Minimum-order stable recursive filter design via the genetic algorithm approach

In this paper, the minimum-order stable recursive filter design problem is proposed and investigated. This problem is playing an important role in pipeline implementation sin signal processing. Here, the existence of a high-order stable recursive filter is proved theoretically, in which the upper bound for the highest order of stable filters is given. Then the minimum-order stable linear predictor is obtained via solving an optimization problem. In this paper, the popular genetic algorithm approach is adopted since it is a heuristic probabilistic optimization technique and has been widely used in engineering designs. Finally, an illustrative example is sued to show the effectiveness of the proposed algorithm.

Scope for further analytical solutions for constant flux infiltration into a semi-infinite soil profile or redistribution in a finite soil profile

[1] We attempt to generate new solutions for the moisture content form of the one-dimensional Richards' [1931] equation using the Lisle [1992] equivalence mapping. This mapping is used as no more general set of transformations exists for mapping the one-dimensional Richards' equation into itself. Starting from a given solution, the mapping has the potential to generate an infinite number of new solutions for a series of nonlinear diffusivity and hydraulic conductivity functions. We first seek new analytical solutions satisfying Richards' equation subject to a constant flux surface boundary condition for a semi-infinite dry soil, starting with the Burgers model. The first iteration produces an existing solution, while subsequent iterations are shown to endlessly reproduce this same solution. Next, we briefly consider the problem of redistribution in a finite-length soil. In this case, Lisle's equivalence mapping is generalized to account for arbitrary initial conditions. As was the case for infiltration, however, it is found that new analytical solutions are not generated using the equivalence mapping, although existing solutions are recovered.

An equivalent algorithm for simulating thermal effects of magma intrusion problems in porous rocks

An equivalent algorithm is proposed to simulate thermal effects of the magma intrusion in geological systems, which are composed of porous rocks. Based on the physical and mathematical equivalence, the original magma solidification problem with a moving boundary between the rock and intruded magma is transformed into a new problem without the moving boundary but with a physically equivalent heat source. From the analysis of an ideal solidification model, the physically equivalent heat source has been determined in this paper. The major advantage in using the proposed equivalent algorithm is that the fixed finite element mesh with a variable integration time step can be employed to simulate the thermal effect of the intruded magma solidification using the conventional finite element method. The related numerical results have demonstrated the correctness and usefulness of the proposed equivalent algorithm for simulating the thermal effect of the intruded magma solidification in geological systems. (C) 2003 Elsevier B.V. All rights reserved.

A spectral clustering algorithm for manufacturing cell formation

A graph clustering algorithm constructs groups of closely related parts and machines separately. After they are matched for the least intercell moves, a refining process runs on the initial cell formation to decrease the number of intercell moves. A simple modification of this main approach can deal with some practical constraints, such as the popular constraint of bounding the maximum number of machines in a cell. Our approach makes a big improvement in the computational time. More importantly, improvement is seen in the number of intercell moves when the computational results were compared with best known solutions from the literature. (C) 2009 Elsevier Ltd. All rights reserved.

How does a colony of Apoica flavissima (Hymenoptera: Vespidae, Epiponini) maintain a constant temperature?

The thermal characteristics of a colony of Apoica flavissima, an epiponine wasp, were examined. The nest, with a diameter of slightly less than 30 cm, was built on a twig of an orange tree. The temperature of the roof surface fluctuated greatly, ranging between 19.1 and 41.5 degrees C. However, the temperature in the central cell was kept constant at around 27 degrees C throughout a day. Although heavy rain pelted the nest roof in the morning, the central cell maintained temperatures higher than 25 degrees C. On the contrary, after all immature and adult wasps were removed the temperature in the nest fluctuated considerably. The presence of immature individuals and adult wasps densely covering the under surface of the comb seemed to function as an effective insulator. The smaller temperature fluctuation in the central cell than on the roof surface, when the nest was in the empty state, suggests that the thick spongy tissue of the roof made from curled plant leaf hairs serves as an insulator to prevent the conduction of solar heat into the cells and the outward flow of heat generated in cells, especially at night.

Evaluation of elastic modulus and hardness of thin films by nanoindentation

Simple equations are proposed for determining elastic modulus and hardness properties of thin films on substrates from nanoindentation experiments. An empirical formulation relates the modulus E and hardness H of the film/substrate bilayer to corresponding material properties of the constituent materials via a power-law relation. Geometrical dependence of E and H is wholly contained in the power-law exponents, expressed here as sigmoidal functions of indenter penetration relative to film thickness. The formulation may be inverted to enable deconvolution of film properties from data on the film/substrate bilayers. Berkovich nanoindentation data for dense oxide and nitride films on silicon substrates are used to validate the equations and to demonstrate the film property deconvolution. Additional data for less dense nitride films are used to illustrate the extent to which film properties may depend on the method of fabrication.

Axial Young's modulus prediction of single-walled carbon nanotube arrays with diameters from nanometer to meter scales

Based on a self-similar array model, we systematically investigated the axial Young's modulus (Y-axis) of single-walled carbon nanotube (SWNT) arrays with diameters from nanometer to meter scales by an analytical approach. The results show that the Y-axis of SWNT arrays decreases dramatically with the increases of their hierarchy number (s) and is not sensitive to the specific size and constitution when s is the same, and the specific Young's modulus Y-axis(s) is independent of the packing configuration of SWNTs. Our calculations also show that the Y-axis of SWNT arrays with diameters of several micrometers is close to that of commercial high performance carbon fibers (CFs), but the Y-axis(s) of SWNT arrays is much better than that of high performance CFs. (C) 2005 American Institute of Physics.

A new algorithm and refined bounds for extended gcd computation

Extended gcd computation is interesting itself. It also plays a fundamental role in other calculations. We present a new algorithm for solving the extended gcd problem. This algorithm has a particularly simple description and is practical. It also provides refined bounds on the size of the multipliers obtained.

The Qu-Prolog unification algorithm: Formalisation and correctness

Qu-Prolog is an extension of Prolog which performs meta-level computations over object languages, such as predicate calculi and lambda-calculi, which have object-level variables, and quantifier or binding symbols creating local scopes for those variables. As in Prolog, the instantiable (meta-level) variables of Qu-Prolog range over object-level terms, and in addition other Qu-Prolog syntax denotes the various components of the object-level syntax, including object-level variables. Further, the meta-level operation of substitution into object-level terms is directly represented by appropriate Qu-Prolog syntax. Again as in Prolog, the driving mechanism in Qu-Prolog computation is a form of unification, but this is substantially more complex than for Prolog because of Qu-Prolog's greater generality, and especially because substitution operations are evaluated during unification. In this paper, the Qu-Prolog unification algorithm is specified, formalised and proved correct. Further, the analysis of the algorithm is carried out in a frame-work which straightforwardly allows the 'completeness' of the algorithm to be proved: though fully explicit answers to unification problems are not always provided, no information is lost in the unification process.

Determination of the force constant of a single-beam gradient trap by measurement of backscattered light

A single-beam gradient trap could potentially be used to hold a stylus for scanning force microscopy. With a view to development of this technique, we modeled the optical trap as a harmonic oscillator and therefore characterized it by its force constant. We measured force constants and resonant frequencies for 1-4-mu m-diameter polystyrene spheres in a single-beam gradient trap using measurements of back-scattered light. Force constants were determined with both Gaussian and doughnut laser modes, with powers of 3 and 1 mW, respectively. Typical values for spring constants were measured to be between 10(-6) and 4 x 10(-6) N/m. The resonant frequencies of trapped particles were measured to be between 1 and 10 kHz, and the rms amplitudes of oscillations were estimated to be around 40 nm. Our results confirm that the use of the doughnut mode for single-beam trapping is more efficient in the axial direction. (C) 1996 Optical Society of America.

Ramp and constant power trials produce equivalent critical power estimates

The standard critical power test protocol on the cycle prescribes a series of trials to exhaustion, each at a different but constant power setting. Recently the protocol has been modified and applied to a series of trials to exhaustion each at a different ramp incremental rate. This study was undertaken to compare critical power and anaerobic work capacity estimates in the same group of subjects when derived from the two protocols. Ten male subjects of mixed athletic ability cycled to exhaustion on eight occasions in randomized order over a 3-wk period. Four trials were performed at differing constant power settings and four trials on differing ramp incremental rates. Both critical power and anaerobic work capacity were estimated for each subject by curve fitting of the ramp model and of three versions of the constant power model. After adjusting for inter-subject variability, no significant differences were detected between critical power estimates or between anaerobic work capacity estimates from any model formulation or from the two protocols. It is concluded that both the ramp and constant power protocols produce equivalent estimates for critical power and anaerobic work capacity.

The subcycled Newmark algorithm

The popular Newmark algorithm, used for implicit direct integration of structural dynamics, is extended by means of a nodal partition to permit use of different timesteps in different regions of a structural model. The algorithm developed has as a special case an explicit-explicit subcycling algorithm previously reported by Belytschko, Yen and Mullen. That algorithm has been shown, in the absence of damping or other energy dissipation, to exhibit instability over narrow timestep ranges that become narrower as the number of degrees of freedom increases, making them unlikely to be encountered in practice. The present algorithm avoids such instabilities in the case of a one to two timestep ratio (two subcycles), achieving unconditional stability in an exponential sense for a linear problem. However, with three or more subcycles, the trapezoidal rule exhibits stability that becomes conditional, falling towards that of the central difference method as the number of subcycles increases. Instabilities over narrow timestep ranges, that become narrower as the model size increases, also appear with three or more subcycles. However by moving the partition between timesteps one row of elements into the region suitable for integration with the larger timestep these the unstable timestep ranges become extremely narrow, even in simple systems with a few degrees of freedom. As well, accuracy is improved. Use of a version of the Newmark algorithm that dissipates high frequencies minimises or eliminates these narrow bands of instability. Viscous damping is also shown to remove these instabilities, at the expense of having more effect on the low frequency response.

Simulated-annealing-based genetic algorithm for modeling the optical constants of solids

We propose a simulated-annealing-based genetic algorithm for solving model parameter estimation problems. The algorithm incorporates advantages of both genetic algorithms and simulated annealing. Tests on computer-generated synthetic data that closely resemble optical constants of a metal were performed to compare the efficiency of plain genetic algorithms against the simulated-annealing-based genetic algorithms. These tests assess the ability of the algorithms to and the global minimum and the accuracy of values obtained for model parameters. Finally, the algorithm with the best performance is used to fit the model dielectric function to data for platinum and aluminum. (C) 1997 Optical Society of America.