Coarse-grid higher order finite-difference time-domain algorithm with low dispersion errors

Higher order (2,4) FDTD schemes used for numerical solutions of Maxwell`s equations are focused on diminishing the truncation errors caused by the Taylor series expansion of the spatial derivatives. These schemes use a larger computational stencil, which generally makes use of the two constant coefficients, C-1 and C-2, for the four-point central-difference operators. In this paper we propose a novel way to diminish these truncation errors, in order to obtain more accurate numerical solutions of Maxwell`s equations. For such purpose, we present a method to individually optimize the pair of coefficients, C-1 and C-2, based on any desired grid size resolution and size of time step. Particularly, we are interested in using coarser grid discretizations to be able to simulate electrically large domains. The results of our optimization algorithm show a significant reduction in dispersion error and numerical anisotropy for all modeled grid size resolutions. Numerical simulations of free-space propagation verifies the very promising theoretical results. The model is also shown to perform well in more complex, realistic scenarios.

Three-Dimensional Stress Distribution in the Human Periodontal Ligament in Masticatory, Parafunctional, and Trauma Loads: Finite Element Analysis

Background: The presence of the periodontal ligament (PDL) makes it possible to absorb and distribute loads produced during masticatory function and other tooth contacts into the alveolar process via the alveolar bone proper. However, several factors affect the integrity of periodontal structures causing the destruction of the connective matrix and cells, the loss of fibrous attachment, and the resorption of alveolar bone. Methods: The purpose of this study was to evaluate the stress distribution by finite element analysis in a PDL in three-dimensional models of the upper central incisor under three different load conditions: 100 N occlusal loading at 45 degrees (model 1: masticatory load); 500 N at the incisal edge at 45 degrees (model 2: parafunctional habit); and 800 N at the buccal surface at 90 degrees (model 3: trauma case). The models were built from computed tomography scans. Results: The stress distribution was quite different among the models. The most significant values (harmful) of tensile and compressive stresses were observed in models 2 and 3, with similarly distinct patterns of stress distributions along the PDL. Tensile stresses were observed along the internal and external aspects of the PDL, mostly at the cervical and middle thirds. Conclusions: The stress generation in these models may affect the integrity of periodontal structures. A better understanding of the biomechanical behavior of the PDL under physiologic and traumatic loading conditions might enhance the understanding of the biologic reaction of the PDL in health and disease. J Periodontol 2009;80:1859-1867.

Determination of pore size distributions by regularization and finite element collocation

The problem of extracting pore size distributions from characterization data is solved here with particular reference to adsorption. The technique developed is based on a finite element collocation discretization of the adsorption integral, with fitting of the isotherm data by least squares using regularization. A rapid and simple technique for ensuring non-negativity of the solutions is also developed which modifies the original solution having some negativity. The technique yields stable and converged solutions, and is implemented in a package RIDFEC. The package is demonstrated to be robust, yielding results which are less sensitive to experimental error than conventional methods, with fitting errors matching the known data error. It is shown that the choice of relative or absolute error norm in the least-squares analysis is best based on the kind of error in the data. (C) 1998 Elsevier Science Ltd. All rights reserved.

Finite element modelling of temperature gradient driven rock alteration and mineralization in porous rock masses

We present finite element simulations of temperature gradient driven rock alteration and mineralization in fluid saturated porous rock masses. In particular, we explore the significance of production/annihilation terms in the mass balance equations and the dependence of the spatial patterns of rock alteration upon the ratio of the roll over time of large scale convection cells to the relaxation time of the chemical reactions. Special concepts such as the gradient reaction criterion or rock alteration index (RAI) are discussed in light of the present, more general theory. In order to validate the finite element simulation, we derive an analytical solution for the rock alteration index of a benchmark problem on a two-dimensional rectangular domain. Since the geometry and boundary conditions of the benchmark problem can be easily and exactly modelled, the analytical solution is also useful for validating other numerical methods, such as the finite difference method and the boundary element method, when they are used to dear with this kind of problem. Finally, the potential of the theory is illustrated by means of finite element studies related to coupled flow problems in materially homogeneous and inhomogeneous porous rock masses. (C) 1998 Elsevier Science S.A. All rights reserved.

New techniques in nutritional assessment: body composition methods

New techniques in air-displacement plethysmography seem to have overcome many of the previous problems of poor reproducibility and validity. These have made body-density measurements available to a larger range of individuals, including children, elderly and sick patients who often have difficulties in being submerged underwater in hydrodensitometry systems. The BOD POD air-displacement system (BOD POD body composition system; Life Measurement Instruments, Concord, CA, USA) is more precise than hydrodensitometry, is simple and rapid to operate (approximately 1 min measurements) and the results agree closely with those of hydrodensitometry (e.g. +/-3.4% for estimation of body fat). Body line scanners employing the principles of three-dimensional photography are potentially able to measure the surface area and volume of the body and its segments even more rapidly (approximately 10 s), but the validity of the measurements needs to be established. Advances in i.r. spectroscopy and mathematical modelling for calculating the area under the curve have improved precision for measuring enrichment of (H2O)-H-2 in studies of water dilution (CV 0.1-0.9% within the range of 400-1000 mu l/l) in saliva, plasma and urine. The technique is rapid and compares closely with mass spectrometry (bias 1 (SD 2) %). Advances in bedside bioelectrical-impedance techniques are making possible potential measurements of skinfold thicknesses and limb muscle mass electronically. Preliminary results suggest that the electronic method is more reproducible (intra-and inter-individual reproducibility for measuring skinfold thicknesses) and associated with less bias (+ 12%), than anthropometry (+ 40%). In addition to these selected examples, the 'mobility' or transfer of reference methods between centres has made the distinction between reference and bedside or field techniques less distinct than in the past.

Finite element analysis of flow patterns near geological lenses in hydrodynamic and hydrothermal systems

We use theoretical and numerical methods to investigate the general pore-fluid flow patterns near geological lenses in hydrodynamic and hydrothermal systems respectively. Analytical solutions have been rigorously derived for the pore-fluid velocity, stream function and excess pore-fluid pressure near a circular lens in a hydrodynamic system. These analytical solutions provide not only a better understanding of the physics behind the problem, but also a valuable benchmark solution for validating any numerical method. Since a geological lens is surrounded by a medium of large extent in nature and the finite element method is efficient at modelling only media of finite size, the determination of the size of the computational domain of a finite element model, which is often overlooked by numerical analysts, is very important in order to ensure both the efficiency of the method and the accuracy of the numerical solution obtained. To highlight this issue, we use the derived analytical solutions to deduce a rigorous mathematical formula for designing the computational domain size of a finite element model. The proposed mathematical formula has indicated that, no matter how fine the mesh or how high the order of elements, the desired accuracy of a finite element solution for pore-fluid flow near a geological lens cannot be achieved unless the size of the finite element model is determined appropriately. Once the finite element computational model has been appropriately designed and validated in a hydrodynamic system, it is used to examine general pore-fluid flow patterns near geological lenses in hydrothermal systems. Some interesting conclusions on the behaviour of geological lenses in hydrodynamic and hydrothermal systems have been reached through the analytical and numerical analyses carried out in this paper.

Formal derivation of finite state machines for class testing

Previous work on generating state machines for the purpose of class testing has not been formally based. There has also been work on deriving state machines from formal specifications for testing non-object-oriented software. We build on this work by presenting a method for deriving a state machine for testing purposes from a formal specification of the class under test. We also show how the resulting state machine can be used as the basis for a test suite developed and executed using an existing framework for class testing. To derive the state machine, we identify the states and possible interactions of the operations of the class under test. The Test Template Framework is used to formally derive the states from the Object-Z specification of the class under test. The transitions of the finite state machine are calculated from the derived states and the class's operations. The formally derived finite state machine is transformed to a ClassBench testgraph, which is used as input to the ClassBench framework to test a C++ implementation of the class. The method is illustrated using a simple bounded queue example.

New methods for determining quasi-stationary distributions for Markov chains

We shall be concerned with the problem of determining quasi-stationary distributions for Markovian models directly from their transition rates Q. We shall present simple conditions for a mu-invariant measure m for Q to be mu-invariant for the transition function, so that if m is finite, it can be normalized to produce a quasi-stationary distribution. (C) 2000 Elsevier Science Ltd. All rights reserved.

Finite element modelling of heat transfer through permeable cracks in hydrothermal systems with upward throughflow

We use the finite element method to model the heat transfer phenomenon through permeable cracks in hydrothermal systems with upward throughflow. Since the finite element method is an approximate numerical method, the method must be validated before it is used to soh,e any new, kind of problem. However, the analytical solution, which can be used to validate the finite element method and other numerical methods, is rather limited in the literature, especially, for the problem considered here. Keeping this in mind, we have derived analytical solutions for the temperature distribution along the vertical axis of a crack in a fluid-saturated porous layer. After the finite element method is validated by comparing the numerical solution with the analytical solution for the same benchmark problem, it is used to investigate the pore-fluid flow and heat transfer in layered hydrothermal systems with vertical permeable cracks. The related analytical and numerical results have demonstrated that vertical cracks are effective and efficient members to transfer heat energy from the bottom section to the top section in hydrothermal systems with upward throughflow.

Finite-difference time-domain-based studies of MRI pulsed field gradient-induced eddy currents inside the human body

In modern magnetic resonance imaging (MRI), patients are exposed to strong, rapidly switching magnetic gradient fields that, in extreme cases, may be able to elicit nerve stimulation. This paper presents theoretical investigations into the spatial distribution of induced current inside human tissues caused by pulsed z-gradient fields. A variety of gradient waveforms have been studied. The simulations are based on a new, high-definition, finite-difference time-domain method and a realistic inhomogeneous 10-mm resolution human body model with appropriate tissue parameters. it was found that the eddy current densities are affected not only by the pulse sequences but by many parameters such as the position of the body inside the gradient set, the local biological material properties and the geometry of the body. The discussion contains a comparison of these results with previous results found in the literature. This study and the new methods presented herein will help to further investigate the biological effects caused by the switched gradient fields in a MRI scan. (C) 2002 Wiley Periodicals, Inc.

The effect of tidal volume and peep on CO2 and oxygenation during resuscitation of very premature lambs

Background: Over-ventilation causing low arterial carbon dioxide levels (PaCO2) has been associated with the development of neonatal chronic lung disease and adverse outcomes. This may occur very soon after birth. Aim: To investigate the effect on PaCO2 and oxygenation of very premature lambs resuscitated with different tidal volumes and PEEP. Methods: Anaesthetised lambs delivered at 126 days gestation were randomised to 15 min resuscitation with 3 regimes: (1) Laerdal resuscitation bag (B) with 100% oxygen and no PEEP, (2) fixed tidal volume (VT) of 5 mL/kg, or (3) VT of 10 mL/kg, both delivered with a Babylog 8000 ventilator in volume guarantee mode with 8 cm H2O PEEP and variable FiO2. Frequent blood gases were measured and VT, mean airway pressure (Paw), minute volume (MV), ventilation rate (VR), respiratory system compliance (Crs) and alveolar/arterial oxygen difference (AaDO2) were recorded. Results: Twenty lambs were studied. B (1) was associated with more variable VT and peak inspiratory pressures (PIP) compared to fixed tidal volumes (2 and 3). The lambs ventilated with 10 mL/kg were over-ventilated, those ventilated with 5 mL/kg were slightly under-ventilated. Those ventilated with the Laerdal bag had a mean VT of 7.5 mL/kg and were normocarbic. The different tidal volumes had little effect on oxygenation. PEEP improved oxygenation. The table shows the values at 15 minutes expressed as mean and SEM. TABLE. No caption av... TABLE. No caption av... Image Tools Conclusion: Very premature lambs can be effectively resuscitated from birth using volume guarantee ventilation. Within minutes of birth different tidal volumes had a large effect on PaCO2 and no effect on oxygenation. Studies are needed to determine the appropriate tidal volume for resuscitating very premature infants to maintain acceptable levels of PaCO2. © International Pediatrics Research Foundation, Inc. 2004. All Rights Reserved.

Effects of tidal volume and positive end-expiratory pressure during resuscitation of very premature lambs

Background: Guidelines recommend neonatal resuscitation without controlling tidal volume or positive end-expiratory pressure (PEEP). However, these may improve gas exchange, lung volume and outcome. Aim: To investigate resuscitation of very premature lambs with a Laerdal bag without PEEP versus volume guarantee ventilation with PEEP. Methods: Anaesthetized lambs (n = 20) delivered at 125 d gestation were randomized to three groups receiving 15 min resuscitation: (1) Laerdal bag and no PEEP; (2) ventilation with a tidal volume of 5 ml/kg and 8 cm H2O PEEP; (3) ventilation with 10 ml/kg and 8 cm H2O PEEP. They were then all ventilated for 2 h with tidal volumes of 5 or 10 ml/kg, and 8 cm H2O PEEP. Ventilation parameters and blood gases were recorded. Results: Different tidal volumes affected PaCO2 within minutes, with 10 ml/kg causing severe hypocarbia. PEEP had little effect on PaCO2. Oxygenation improved significantly with PEEP of 8 cm H2O, irrespective of tidal volume. Conclusion: Very premature lambs can be resuscitated effectively using volume-guarantee ventilation and PEEP. Tidal volumes affected PaCO2 within minutes but had little effect on oxygenation. PEEP halved the oxygen requirement compared with no PEEP. Resuscitating premature babies with controlled tidal volumes and PEEP might improve their outcome.

A comparison of methods to identify open-lung PEEP

Purpose: Many methods exist in the literature for identifying PEEP to set in ARDS patients following a lung recruitment maneuver (RM). We compared ten published parameters for setting PEEP following a RM. Methods: Lung injury was induced by bilateral lung lavage in 14 female Dorset sheep, yielding a PaO(2) 100-150 mmHg at F(I)O(2) 1.0 and PEEP 5 cmH(2)O. A quasi-static P-V curve was then performed using the supersyringe method; PEEP was set to 20 cmH(2)O and a RM performed with pressure control ventilation (inspiratory pressure set to 40-50 cmH(2)O), until PaO(2) + PaCO(2) > 400 mmHg. Following the RM, a decremental PEEP trial was performed. The PEEP was decreased in 1 cmH(2)O steps every 5 min until 15 cmH(2)O was reached. Parameters measured during the decremental PEEP trial were compared with parameters obtained from the P-V curve. Results: For setting PEEP, maximum dynamic tidal respiratory compliance, maximum PaO(2), maximum PaO(2) + PaCO(2), and minimum shunt calculated during the decremental PEEP trial, and the lower Pflex and point of maximal compliance increase on the inflation limb of the P-V curve (Pmci,i) were statistically indistinguishable. The PEEP value obtained using the deflation upper Pflex and the point of maximal compliance decrease on the deflation limb were significantly higher, and the true inflection point on the inflation limb and minimum PaCO(2) were significantly lower than the other variables. Conclusion: In this animal model of ARDS, dynamic tidal respiratory compliance, maximum PaO(2), maximum PaO(2) + PaCO(2), minimum shunt, inflation lower Pflex and Pmci,i yield similar values for PEEP following a recruitment maneuver.

Effectiveness of the CAGE questionnaire, gamma-glutamyltransferase and mean corpuscular volume of red blood cells as markers for alcohol-related problems in the workplace

Objective: To evaluate the usefulness of gamma-glutamyltransferase (GGT) and mean corpuscular volume (MCV), as well as that of the CAGE questionnaire, in workplace screening for alcohol abuse/dependence. Methods: A total of 183 male employees were submitted to structured interviews (Structured Clinical Interview for DSM-IV 2.0 and CAGE questionnaire). Blood samples were collected. Diagnostic accuracy and odds ratio were determined for the CAGE, GGT and MCV. Results: The CAGE questionnaire presented the best sensitivity for alcohol dependence (91%; specificity, 87.8%) and for alcohol abuse (87.5%, specificity, 80.9%), which increased when the questionnaire was used in combination with GGT (sensitivity, 100% and 87.5%, respectively; specificity, 68% and 61.5, respectively). CAGE positive results and/or alterations in GGT were less likely to occur among employees not presenting alcohol abuse/ dependence than among those presenting such abuse (OR for CAGE = 13, p < 0.05; OR for CAGE-GGT = 11, p < 0.05) or dependence (OR for CAGE = 76, p < 0.0 1; OR for GGT = 5, p < 0.0 1). Employees not presenting alcohol abuse/dependence were also several times more likely to present negative CAGE or GGT results. Conclusions: The use short, simple questionnaires, combined with that of low-cost biochemical markers, such as GGT, can serve as an initial screening for alcohol-related problems, especially for employees in hazardous occupations. The data provided can serve to corroborate clinical findings. (C) 2008 Elsevier Ltd. All rights reserved.

Urodynamic evaluation of patients with lower urinary tract symptoms and small prostate volume

Introduction: The relevance of prostate size in the pathophysiology of lower urinary tract symptoms (LUTS) is controversial. We evaluated the urodynamic findings in patients with LUTS and small prostate volumes. Materials and Methods: 84 patients aged >= 50 years with LUTS and prostates < 40 ml were evaluated. All had an International Prostate Symptom Score (IPSS) >= 8. Average age was 62.0 +/- 8.1 years. We evaluated the impact of bladder outlet obstruction (BOO) and detrusor overactivity (DO) on the voiding symptoms and urodynamic findings. Results: Mean prostate volume was 29.2 +/- 7.2 ml and mean IPSS was 13.5 +/- 4.6. BOO was the main finding, affecting 42 (50.0%) patients, followed by detrusor underactivity (DU) in 41 (48.8%) and DO in 28 (33.3%) patients. Patients without BOO were significantly older than the obstructed (64.0 +/- 8.8 and 60.1 +/- 6.9 years, respectively; p = 0.026) and had an increased prevalence of DU (76.2 and 21.4%, respectively; p < 0.001). Comparison of patients with and without DO showed reduced bladder capacity and compliance in the DO group (p < 0.001). No other comparisons were significant. Conclusion: Half of the patients with LUTS and small prostates are not obstructed and may have DO or decreased detrusor contractility as the basis for their voiding symptoms. Our results emphasize the value of urodynamics in this population, especially when invasive treatments are being considered. Copyright (c) 2008 S. Karger AG, Basel.