Numerical modelling of chemical effects of magma solidification problems in porous rocks

The solidification of intruded magma in porous rocks can result in the following two consequences: (1) the heat release due to the solidification of the interface between the rock and intruded magma and (2) the mass release of the volatile fluids in the region where the intruded magma is solidified into the rock. Traditionally, the intruded magma solidification problem is treated as a moving interface (i.e. the solidification interface between the rock and intruded magma) problem to consider these consequences in conventional numerical methods. This paper presents an alternative new approach to simulate thermal and chemical consequences/effects of magma intrusion in geological systems, which are composed of porous rocks. In the proposed new approach and algorithm, 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 the proposed mass source and physically equivalent heat source. The major advantage in using the proposed equivalent algorithm is that a fixed mesh of finite elements with a variable integration time-step can be employed to simulate the consequences and effects of the intruded magma solidification using the conventional finite element method. The correctness and usefulness of the proposed equivalent algorithm have been demonstrated by a benchmark magma solidification problem. Copyright (c) 2005 John Wiley & Sons, Ltd.

Numerical simulation of double-diffusion driven convective flow and rock alteration in three-dimensional fluid-saturated geological fault zones

Numerical methods are used to simulate the double-diffusion driven convective pore-fluid flow and rock alteration in three-dimensional fluid-saturated geological fault zones. The double diffusion is caused by a combination of both the positive upward temperature gradient and the positive downward salinity concentration gradient within a three-dimensional fluid-saturated geological fault zone, which is assumed to be more permeable than its surrounding rocks. In order to ensure the physical meaningfulness of the obtained numerical solutions, the numerical method used in this study is validated by a benchmark problem, for which the analytical solution to the critical Rayleigh number of the system is available. The theoretical value of the critical Rayleigh number of a three-dimensional fluid-saturated geological fault zone system can be used to judge whether or not the double-diffusion driven convective pore-fluid flow can take place within the system. After the possibility of triggering the double-diffusion driven convective pore-fluid flow is theoretically validated for the numerical model of a three-dimensional fluid-saturated geological fault zone system, the corresponding numerical solutions for the convective flow and temperature are directly coupled with a geochemical system. Through the numerical simulation of the coupled system between the convective fluid flow, heat transfer, mass transport and chemical reactions, we have investigated the effect of the double-diffusion driven convective pore-fluid flow on the rock alteration, which is the direct consequence of mineral redistribution due to its dissolution, transportation and precipitation, within the three-dimensional fluid-saturated geological fault zone system. (c) 2005 Elsevier B.V. All rights reserved.

Numerical simulation of dendrite growth during solidification

A comprehensive probabilistic model for simulating dendrite morphology and investigating dendritic growth kinetics during solidification has been developed, based on a modified Cellular Automaton (mCA) for microscopic modeling of nucleation, growth of crystals and solute diffusion. The mCA model numerically calculated solute redistribution both in the solid and liquid phases, the curvature of dendrite tips and the growth anisotropy. This modeling takes account of thermal, curvature and solute diffusion effects. Therefore, it can simulate microstructure formation both on the scale of the dendrite tip length. This model was then applied for simulating dendritic solidification of an Al-7%Si alloy. Both directional and equiaxed dendritic growth has been performed to investigate the growth anisotropy and cooling rate on dendrite morphology. Furthermore, the competitive growth and selection of dendritic crystals have also investigated.

Labour flexibility and gender in the service sector: A study of employment practices in law and accounting firms

In seeking to increase the flexibility of their use of employee time, employers can pursue strategies based on the employment of casual and part-time workers (numerical flexibility) or strategies based on ad hoc variation of the working hours of permanent employees (working time flexibility). Patterns of flexibility strategies and their implications are examined in the context of a highly feminised sector of work-clerical and administrative employment in law and accounting firms. We consider whether, as is often assumed, working time flexibility strategies are generally better for employees because they avoid the substitution of core, high quality jobs with the peripheral, relatively insecure employment often associated with casualisation. Analysing data drawn from a survey of law and accounting firms, we argue that there are three distinct flexibility strategies adopted by employers, and that the choice of strategy is influenced by the size of the firm and the extent of feminisation. The quality of employment conditions associated with each strategy is investigated through an analysis of the determinants of training provision for clerical and administrative workers. Rather than an expected simple linear relationship between increasing casualisation and decreasing training provision, we find that firm size and feminisation are implicated. Larger firms that tend to employ at least some men and use a combination of working time and numerical flexibility strategies tend to provide more training than the small, more fully feminised firms that tend to opt for either casualisation or working time flexibility strategies. This suggests that, from an employee perspective, working time flexibility may not be as benevolent as is often thought.

A numerical study of effect of grain boundaries on elastic and plastic properties in nanocomposite materials

Nanocomposite materials have received considerable attention in recent years due to their novel properties. Grain boundaries are considered to play an important role in nanostructured materials. This work focuses on the finite element analysis of the effect of grain boundaries on the overall mechanical properties of aluminium/alumina composites. A grain boundary is incorporated into the commonly used unit cell model to investigate its effect on material properties. By combining the unit cell model with an indentation model, coupled with experimental indentation measurements, the ''effective'' plastic property of the grain boundary is estimated. In addition, the strengthening mechanism is also discussed based on the Estrin-Mecking model.

High-breakdown linear discriminant analysis

The classification rules of linear discriminant analysis are defined by the true mean vectors and the common covariance matrix of the populations from which the data come. Because these true parameters are generally unknown, they are commonly estimated by the sample mean vector and covariance matrix of the data in a training sample randomly drawn from each population. However, these sample statistics are notoriously susceptible to contamination by outliers, a problem compounded by the fact that the outliers may be invisible to conventional diagnostics. High-breakdown estimation is a procedure designed to remove this cause for concern by producing estimates that are immune to serious distortion by a minority of outliers, regardless of their severity. In this article we motivate and develop a high-breakdown criterion for linear discriminant analysis and give an algorithm for its implementation. The procedure is intended to supplement rather than replace the usual sample-moment methodology of discriminant analysis either by providing indications that the dataset is not seriously affected by outliers (supporting the usual analysis) or by identifying apparently aberrant points and giving resistant estimators that are not affected by them.

Numerical modelling of tide-induced beach water table fluctuations

Field studies have shown that the elevation of the beach groundwater table varies with the tide and such variations affect significantly beach erosion or accretion. In this paper, we present a BEM (Boundary Element Method) model for simulating the tidal fluctuation of the beach groundwater table. The model solves the two-dimensional flow equation subject to free and moving boundary conditions, including the seepage dynamics at the beach face. The simulated seepage faces were found to agree with the predictions of a simple model (Turner, 1993). The advantage of the present model is, however, that it can be used with little modification to simulate more complicated cases, e.g., surface recharge from rainfall and drainage in the aquifer may be included (the latter is related to beach dewatering technique). The model also simulated well the field data of Nielsen (1990). In particular, the model replicated three distinct features of local water table fluctuations: steep rising phase versus flat falling phase, amplitude attenuation and phase lagging.

Numerical and experimental study of a homogenizer impinging jet

High-pressure homogenization is a key unit operation used to disrupt cells containing intracellular bioproducts. Modeling and optimization of this unit are restrained by a lack of information on the flow conditions within a homogenizer value. A numerical investigation of the impinging radial jet within a homogenizer value is presented. Results for a laminar and turbulent (k-epsilon turbulent model) jet are obtained using the PHOENICS finite-volume code. Experimental measurement of the stagnation region width and correlation of the cell disruption efficiency with jet stagnation pressure both indicate that the impinging jet in the homogenizer system examined is likely to be laminar under normal operating conditions. Correlation of disruption data with laminar stagnation pressure provides a better description of experimental variability than existing correlations using total pressure drop or the grouping 1/Y(2)h(2).

Anisotropic correlated electron model associated with the Temperley-Lieb algebra

We present an anisotropic correlated electron model on a periodic lattice, constructed from an R-matrix associated with the Temperley-Lieb algebra. By modification of the coupling of the first and last sites we obtain a model with quantum algebra invariance.

Linear IgA dermatosis induced by pregnancy

Linear IgA disease is a rare autoimmune subepidermal bullous disorder characterized by linear IgA deposits at the epidermal basement membrane zone. According to the literature, in patients who have linear IgA disease and become pregnant, the disease tends to improve. We report a case of linear IgA disease induced by pregnancy, successfully treated with dapsone and prednisone with no adverse effects observed in the patient and her newborns.

Leptospirosis pulmonary haemorrhage syndrome is associated with linear deposition of immunoglobulin and complement on the alveolar surface

Leptospirosis is a zoonotic infection associated with severe diseases such as leptospirosis pulmonary haemorrhage syndrome (LPHS). The cause of pulmonary haemorrhage is unclear. Understanding which mechanisms and processes are involved in LPHS will be important in treatment regimens under development for this life-threatening syndrome. In the present study, we evaluated 30 lung specimens from LPHS patients and seven controls using histology and immunohistochemistry (detection of IgM, IgG, IgA and C3) in order to describe the pathological features associated with this syndrome. Immunoglobulin deposits were detected on the alveolar surface in 18/30 LPHS patients. Three staining patterns were observed for the immunoglobulins and C3 in the lung tissues of LPHS patients: AS, delicate linear staining adjacent to the alveolar surface, which was indicative of a membrane covering the luminal surface of type I and II pneumocyte cells; S, heterogeneous staining which was sporadically distributed along the alveolar septum; and IA, weak, focal intra-alveolar granular staining. Human LPHS is associated with individual and unique histological patterns that differ from those of other causes of pulmonary haemorrhage. In the present study, it was found that the linear deposition of immunoglobulins (IgA, IgG and IgM) and complement on the alveolar surface may play a role in the pathogenesis of pulmonary haemorrhage in human leptospirosis.

Comments on the Drinfeld realization of the quantum affine superalgebra U-q[gl(m backslash n)(1)] and its Hopf algebra structure

By generalizing the Reshetikhin and Semenov-Tian-Shansky construction to supersymmetric cases, we obtain the Drinfeld current realization for the quantum affine superalgebra U-q[gl(m\n)((1))]. We find a simple coproduct for the quantum current generators and establish the Hopf algebra structure of this super current algebra.

Direct measurement of instantaneous source speed for a HDR brachytherapy unit using an optical fiber based detector

Purpose: Several attempts to determine the transit time of a high dose rate (HDR) brachytherapy unit have been reported in the literature with controversial results. The determination of the source speed is necessary to accurately calculate the transient dose in brachytherapy treatments. In these studies, only the average speed of the source was measured as a parameter for transit dose calculation, which does not account for the realistic movement of the source, and is therefore inaccurate for numerical simulations. The purpose of this work is to report the implementation and technical design of an optical fiber based detector to directly measure the instantaneous speed profile of a (192)Ir source in a Nucletron HDR brachytherapy unit. Methods: To accomplish this task, we have developed a setup that uses the Cerenkov light induced in optical fibers as a detection signal for the radiation source moving inside the HDR catheter. As the (192)Ir source travels between two optical fibers with known distance, the threshold of the induced signals are used to extract the transit time and thus the velocity. The high resolution of the detector enables the measurement of the transit time at short separation distance of the fibers, providing the instantaneous speed. Results: Accurate and high resolution speed profiles of the 192Ir radiation source traveling from the safe to the end of the catheter and between dwell positions are presented. The maximum and minimum velocities of the source were found to be 52.0 +/- 1.0 and 17.3 +/- 1:2 cm/s. The authors demonstrate that the radiation source follows a uniformly accelerated linear motion with acceleration of vertical bar a vertical bar = 113 cm/s(2). In addition, the authors compare the average speed measured using the optical fiber detector to those obtained in the literature, showing deviation up to 265%. Conclusions: To the best of the authors` knowledge, the authors directly measured for the first time the instantaneous speed profile of a radiation source in a HDR brachytherapy unit traveling from the unit safe to the end of the catheter and between interdwell distances. The method is feasible and accurate to implement on quality assurance tests and provides a unique database for efficient computational simulations of the transient dose. (C) 2010 American Association of Physicists in Medicine. [DOI: 10.1118/1.3483780]

Obesity due to Melanocortin 4 Receptor (MC4R) Deficiency Is Associated with Increased Linear Growth and Final Height, Fasting Hyperinsulinemia, and Incompletely Suppressed Growth Hormone Secretion

Context: Melanocortin receptor 4 (MC4R) deficiency is characterized by increased linear growth greater than expected for the degree of obesity. Objective: The objective of the investigation was to study the somatotroph axis in obese MC4R-deficient patients and equally obese controls. Patients and Methods: We obtained anthropometric measurements and insulin concentrations in 153 MC4R-deficient subjects and 1392 controls matched for age and severity of obesity. We measured fasting IGF-I, IGF-II, IGF binding protein (IGFBP)-1, IGFBP-3, and acid-labile subunit levels in a subset of 33 MC4R-deficient patients and 36 control subjects. We examined pulsatile GH secretion in six adult MC4R-deficient subjects and six obese controls. Results: Height so score was significantly greater in MC4R-deficient children under 5 yr of age compared with controls (mean +/- SEM: 2.3 +/- 0.06 vs. 1.8 +/- 0.04, P < 0.001), an effect that persisted throughout childhood. Final height (cm) was greater in MC4R-deficient men (mean +/- SEM 173 +/- 2.5 vs. 168 +/- 2.1, P < 0.001) and women (mean 165 +/- 2.1 vs. 158 +/- 1.9, P < 0.001). Fasting IGF-I, IGF-II, acid-labile subunit, and IGFBP-3 concentrations were similar in the two groups. GH levels were markedly suppressed in obese controls, but pulsatile GH secretion was retained in MC4R deficiency. The mean maximal GH secretion rate per burst (P < 0.05) and mass per burst (P < 0.05) were increased in MC4R deficiency, consistent with increased pulsatile and total GH secretion. Fasting insulin levels were markedly elevated in MC4R-deficient children. Conclusions: In MC4R deficiency, increased linear growth in childhood leads to increased adult final height, greater than predicted by obesity alone. GH pulsatility is maintained in MC4R deficiency, a finding consistent with animal studies, suggesting a role for MC4R in controlling hypothalamic somatostatinergic tone. Fasting insulin levels are significantly higher in children carrying MC4R mutations. Both of these factors may contribute to the accelerated growth phenotype characteristic of MC4R deficiency. (J Clin Endocrinol Metab 96: E181-E188, 2011)