Three-dimensional, fully adaptive simulations of phase-field fluid models

We present an efficient numerical methodology for the 31) computation of incompressible multi-phase flows described by conservative phase-field models We focus here on the case of density matched fluids with different viscosity (Model H) The numerical method employs adaptive mesh refinements (AMR) in concert with an efficient semi-implicit time discretization strategy and a linear, multi-level multigrid to relax high order stability constraints and to capture the flow`s disparate scales at optimal cost. Only five linear solvers are needed per time-step. Moreover, all the adaptive methodology is constructed from scratch to allow a systematic investigation of the key aspects of AMR in a conservative, phase-field setting. We validate the method and demonstrate its capabilities and efficacy with important examples of drop deformation, Kelvin-Helmholtz instability, and flow-induced drop coalescence (C) 2010 Elsevier Inc. All rights reserved

A Robust, Fully Adaptive Hybrid Level-Set/Front-Tracking Method for Two-Phase Flows with an Accurate Surface Tension Computation

We present a variable time step, fully adaptive in space, hybrid method for the accurate simulation of incompressible two-phase flows in the presence of surface tension in two dimensions. The method is based on the hybrid level set/front-tracking approach proposed in [H. D. Ceniceros and A. M. Roma, J. Comput. Phys., 205, 391400, 2005]. Geometric, interfacial quantities are computed from front-tracking via the immersed-boundary setting while the signed distance (level set) function, which is evaluated fast and to machine precision, is used as a fluid indicator. The surface tension force is obtained by employing the mixed Eulerian/Lagrangian representation introduced in [S. Shin, S. I. Abdel-Khalik, V. Daru and D. Juric, J. Comput. Phys., 203, 493-516, 2005] whose success for greatly reducing parasitic currents has been demonstrated. The use of our accurate fluid indicator together with effective Lagrangian marker control enhance this parasitic current reduction by several orders of magnitude. To resolve accurately and efficiently sharp gradients and salient flow features we employ dynamic, adaptive mesh refinements. This spatial adaption is used in concert with a dynamic control of the distribution of the Lagrangian nodes along the fluid interface and a variable time step, linearly implicit time integration scheme. We present numerical examples designed to test the capabilities and performance of the proposed approach as well as three applications: the long-time evolution of a fluid interface undergoing Rayleigh-Taylor instability, an example of bubble ascending dynamics, and a drop impacting on a free interface whose dynamics we compare with both existing numerical and experimental data.

A Study on Unconditionally Stable FDTD Methods for the Modeling of Metamaterials

We assess the performance of three unconditionally stable finite-difference time-domain (FDTD) methods for the modeling of doubly dispersive metamaterials: 1) locally one-dimensional FDTD; 2) locally one-dimensional FDTD with Strang splitting; and (3) alternating direction implicit FDTD. We use both double-negative media and zero-index media as benchmarks.

Distinct conformational properties determined by implicit and explicit representation of protein-solvent interactions. An analytical and computer simulation study

In the protein folding problem, solvent-mediated forces are commonly represented by intra-chain pairwise contact energy. Although this approximation has proven to be useful in several circumstances, it is limited in some other aspects of the problem. Here we show that it is possible to achieve two models to represent the chain-solvent system. one of them with implicit and other with explicit solvent, such that both reproduce the same thermodynamic results. Firstly, lattice models treated by analytical methods, were used to show that the implicit and explicitly representation of solvent effects can be energetically equivalent only if local solvent properties are time and spatially invariant. Following, applying the same reasoning Used for the lattice models, two inter-consistent Monte Carlo off-lattice models for implicit and explicit solvent are constructed, being that now in the latter the solvent properties are allowed to fluctuate. Then, it is shown that the chain configurational evolution as well as the globule equilibrium conformation are significantly distinct for implicit and explicit solvent systems. Actually, strongly contrasting with the implicit solvent version, the explicit solvent model predicts: (i) a malleable globule, in agreement with the estimated large protein-volume fluctuations; (ii) thermal conformational stability, resembling the conformational hear resistance of globular proteins, in which radii of gyration are practically insensitive to thermal effects over a relatively wide range of temperatures; and (iii) smaller radii of gyration at higher temperatures, indicating that the chain conformational entropy in the unfolded state is significantly smaller than that estimated from random coil configurations. Finally, we comment on the meaning of these results with respect to the understanding of the folding process. (C) 2009 Elsevier B.V. All rights reserved.

Approximate controllability of second-order distributed implicit functional systems

In this paper we study the approximate controllability of control systems with states and controls in Hilbert spaces, and described by a second-order semilinear abstract functional differential equation with infinite delay. Initially we establish a characterization for the approximate controllability of a second-order abstract linear system and, in the last section, we compare the approximate controllability of a semilinear abstract functional system with the approximate controllability of the associated linear system. (C) 2008 Elsevier Ltd. All rights reserved.

THE IMPLICIT FUNCTION THEOREM FOR CONTINUOUS FUNCTIONS

In the present paper we obtain a new homological version of the implicit function theorem and some versions of the Darboux theorem. Such results are proved for continuous maps on topological manifolds. As a consequence. some versions of these classic theorems are proved when we consider differenciable (not necessarily C-1) maps.

An implicit technique for solving 3D low Reynolds number moving free surface flows

This paper describes the development of an implicit finite difference method for solving transient three-dimensional incompressible free surface flows. To reduce the CPU time of explicit low-Reynolds number calculations, we have combined a projection method with an implicit technique for treating the pressure on the free surface. The projection method is employed to uncouple the velocity and the pressure fields, allowing each variable to be solved separately. We employ the normal stress condition on the free surface to derive an implicit technique for calculating the pressure at the free surface. Numerical results demonstrate that this modification is essential for the construction of methods that are more stable than those provided by discretizing the free surface explicitly. In addition, we show that the proposed method can be applied to viscoelastic fluids. Numerical results include the simulation of jet buckling and extrudate swell for Reynolds numbers in the range [0.01, 0.5]. (C) 2008 Elsevier Inc. All rights reserved.

Efficient solutions to robust, semi-implicit discretizations of the immersed boundary method

The immersed boundary method is a versatile tool for the investigation of flow-structure interaction. In a large number of applications, the immersed boundaries or structures are very stiff and strong tangential forces on these interfaces induce a well-known, severe time-step restriction for explicit discretizations. This excessive stability constraint can be removed with fully implicit or suitable semi-implicit schemes but at a seemingly prohibitive computational cost. While economical alternatives have been proposed recently for some special cases, there is a practical need for a computationally efficient approach that can be applied more broadly. In this context, we revisit a robust semi-implicit discretization introduced by Peskin in the late 1970s which has received renewed attention recently. This discretization, in which the spreading and interpolation operators are lagged. leads to a linear system of equations for the inter-face configuration at the future time, when the interfacial force is linear. However, this linear system is large and dense and thus it is challenging to streamline its solution. Moreover, while the same linear system or one of similar structure could potentially be used in Newton-type iterations, nonlinear and highly stiff immersed structures pose additional challenges to iterative methods. In this work, we address these problems and propose cost-effective computational strategies for solving Peskin`s lagged-operators type of discretization. We do this by first constructing a sufficiently accurate approximation to the system`s matrix and we obtain a rigorous estimate for this approximation. This matrix is expeditiously computed by using a combination of pre-calculated values and interpolation. The availability of a matrix allows for more efficient matrix-vector products and facilitates the design of effective iterative schemes. We propose efficient iterative approaches to deal with both linear and nonlinear interfacial forces and simple or complex immersed structures with tethered or untethered points. One of these iterative approaches employs a splitting in which we first solve a linear problem for the interfacial force and then we use a nonlinear iteration to find the interface configuration corresponding to this force. We demonstrate that the proposed approach is several orders of magnitude more efficient than the standard explicit method. In addition to considering the standard elliptical drop test case, we show both the robustness and efficacy of the proposed methodology with a 2D model of a heart valve. (C) 2009 Elsevier Inc. All rights reserved.

Comparasion of two minimally invasive methods on the longevity of glass ionomer cement restorations: short-term results of a pilot study

The purpose of this study was to evaluate the clinical performance of glass ionomer cement (GIC) restorations comparing two minimally invasive methods in permanent teeth after 12 months. Fifty pregnant women (second trimester of pregnancy), mean age 22 ± 5.30 years, were treated by two previously trained operators. The treatment approaches tested were: chemomechanical method (CarisolvTM; MediTeam) and atraumatic restorative treatment (ART). A split-mouth study design was used in which the two treatments were randomly placed in 50 matched pairs of permanent teeth. The chemomechanical method (CM) was the test group and the ART was the control group. The treatments were performed in Public Health Centers. The tested restorative material was a high-strength GIC (Ketac Molar; 3M/ESPE). The restorations were placed according to the ART guidelines. Two calibrated independent examiners evaluated the restorations in accordance with ART criteria. The inter-examiner kappa was 0.97. Data were analyzed using 95% confidence interval on the binomial distribution and Fisher's exact test at 5% significance level. In a 12-month follow-up, 86% of the restorations were evaluated. In the test group (CM), 100% (CI=93.3-100%) of the restorations were considered successful. In the control group (ART) 97.6% (CI=87.4-99.9%) of the restorations were considered successful and 2.4% unsuccessful (marginal defect >0.5 mm). There was no statistically significant difference between the 12-mounth success rate for both groups (Fisher's exact test: P=0.49) and between the two operators (Fisher's exact test: P=1.00). Both minimally invasive methods, chemomechanical method and ART, showed a similar clinical performance after 12 months of follow up.

Histopathological evaluation of different methods of experimental induction of periapical periodontitis

This study evaluated histopathologically different methods of experimental induction of periapical periodontitis. The radiographic and microbiological evaluations have been performed in a previous investigation. Fifty-seven root canals from dogs' teeth were assigned to 4 groups. In GI (n=14) and GII (n=14), the root canals were exposed to oral environment for 180 days; in GIII (n=14) and GIV (n=15) the root canals were exposed for 7 days and then the access cavities were restored and remained sealed for 53 days. The root apices of GI and GIII were perforated, whilst those of GII and GIV remained intact. After induction of periapical periodontitis, the dogs were euthanized. Serial sections were obtained and stained with hematoxylin and eosin. Data of the histopathological evaluation were submitted to Kruskal-Wallis and Dunn's tests at 5% significance level. The inflammatory periapical reaction and resorption of mineralized tissues were less intense in GII than in the other groups (p<0.05). There was no histopathological difference among the experimentally induced periapical lesions in the teeth with coronal sealing. On the other hand, when coronal sealing was not performed, greater intensity of induced periapical periodontitis was observed in the teeth with apical perforation.

Microbial distribution in the root canal system after periapical lesion induction using different methods

The aim of this study was to evaluate the microbial distribution in the root canal system after periapical lesion induction in dogs' teeth using different methods. Fifty-two root canals were assigned to 4 groups (n=13). Groups I and II: root canals were exposed to the oral cavity for 180 days; groups III and IV: root canals were exposed for 7 days and then the coronal openings were sealed for 53 days. The root apices of groups I and III were perforated, while those of groups II and IV remained intact. After the experimental periods, the animals were euthanized and the anatomic pieces containing the roots were processed and stained with the Brown & Brenn method to assess the presence and distribution of microorganisms. The incidence of microorganisms at different sites of the roots and periapical lesions was analyzed statistically by the chi-square test at 5% significance level. All groups presented microorganisms in the entire root canal system. A larger number of microorganisms was observed on the root canal walls, apical delta and dentinal tubules (p<0.05), followed by cementum and cemental resorption areas. In spite of the different periods of exposure to the oral environment, the methods used for induction of periapical periodontitis yielded similar distribution of microorganisms in the root canal system.

Comparative study of two commercially pure titanium casting methods

The interest in using titanium to fabricate removable partial denture (RPD) frameworks has increased, but there are few studies evaluating the effects of casting methods on clasp behavior. OBJECTIVE: This study compared the occurrence of porosities and the retentive force of commercially pure titanium (CP Ti) and cobalt-chromium (Co-Cr) removable partial denture circumferential clasps cast by induction/centrifugation and plasma/vacuum-pressure. MATERIAL AND METHODS: 72 frameworks were cast from CP Ti (n=36) and Co-Cr alloy (n=36; control group). For each material, 18 frameworks were casted by electromagnetic induction and injected by centrifugation, whereas the other 18 were casted by plasma and injected by vacuum-pressure. For each casting method, three subgroups (n=6) were formed: 0.25 mm, 0.50 mm, and 0.75 mm undercuts. The specimens were radiographed and subjected to an insertion/removal test simulating 5 years of framework use. Data were analyzed by ANOVA and Tukey's to compare materials and cast methods (α=0.05). RESULTS: Three of 18 specimens of the induction/centrifugation group and 9 of 18 specimens of plasma/vacuum-pressure cast presented porosities, but only 1 and 7 specimens, respectively, were rejected for simulation test. For Co-Cr alloy, no defects were found. Comparing the casting methods, statistically significant differences (p<0.05) were observed only for the Co-Cr alloy with 0.25 mm and 0.50 mm undercuts. Significant differences were found for the 0.25 mm and 0.75 mm undercuts dependent on the material used. For the 0.50 mm undercut, significant differences were found when the materials were induction casted. CONCLUSION: Although both casting methods produced satisfactory CP Ti RPD frameworks, the occurrence of porosities was greater in the plasma/vacuum-pressure than in the induction/centrifugation method, the latter resulting in higher clasp rigidity, generating higher retention force values.

Comparative study of intermaxillary relationships of manual and swallowing methods

This study compared the mandibular displacement from three methods of centric relation record using an anterior jig associated with (A) chin point guidance, (B) swallowing (control group) and (C) bimanual manipulation. Ten patients aged 25-39 years were selected if they met the following inclusion criteria: complete dentition (up to the second molars), Angle class I and absence of signs and symptoms of temporomandibular disorders and diagnostic casts showing stability in the maximum intercuspation (MI) position. Impressions of maxillary and mandibular arches were made with an irreversible hydrocolloid impression material. Master casts of each patient were obtained, mounted on a microscope table in MI as a reference position and 5 records of each method were made per patient. The mandibular casts were then repositioned with records interposed and new measurements were obtained. The difference between the two readings allowed measuring the displacement of the mandible in the anteroposterior and lateral axes. Data were analyzed statistically by ANOVA and Tukey's test at 5% significance level. There was no statistically significant differences (p>0.05) among the three methods for measuring lateral displacement (A=0.38 ± 0.26, B=0.32 ± 0.25 and C=0.32 ± 0.23). For the anteroposterior displacement (A=2.76 ± 1.43, B=2.46 ± 1.48 and C=2.97 ± 1.51), the swallowing method (B) differed significantly from the others (p<0.05), but no significant difference (p>0.05) was found between chin point guidance (A) and bimanual manipulation (C). In conclusion, the swallowing method produced smaller mandibular posterior displacement than the other methods.

Validity and reliability of methods for the detection of secondary caries around amalgam restorations in primary teeth

Secondary caries has been reported as the main reason for restoration replacement. The aim of this in vitro study was to evaluate the performance of different methods - visual inspection, laser fluorescence (DIAGNOdent), radiography and tactile examination - for secondary caries detection in primary molars restored with amalgam. Fifty-four primary molars were photographed and 73 suspect sites adjacent to amalgam restorations were selected. Two examiners evaluated independently these sites using all methods. Agreement between examiners was assessed by the Kappa test. To validate the methods, a caries-detector dye was used after restoration removal. The best cut-off points for the sample were found by a Receiver Operator Characteristic (ROC) analysis, and the area under the ROC curve (Az), and the sensitivity, specificity and accuracy of the methods were calculated for enamel (D2) and dentine (D3) thresholds. These parameters were found for each method and then compared by the McNemar test. The tactile examination and visual inspection presented the highest inter-examiner agreement for the D2 and D3 thresholds, respectively. The visual inspection also showed better performance than the other methods for both thresholds (Az = 0.861 and Az = 0.841, respectively). In conclusion, the visual inspection presented the best performance for detecting enamel and dentin secondary caries in primary teeth restored with amalgam.

Reliability and discriminatory power of methods for dental plaque quantification

OBJECTIVE: This in situ study evaluated the discriminatory power and reliability of methods of dental plaque quantification and the relationship between visual indices (VI) and fluorescence camera (FC) to detect plaque. MATERIAL AND METHODS: Six volunteers used palatal appliances with six bovine enamel blocks presenting different stages of plaque accumulation. The presence of plaque with and without disclosing was assessed using VI. Images were obtained with FC and digital camera in both conditions. The area covered by plaque was assessed. Examinations were done by two independent examiners. Data were analyzed by Kruskal-Wallis and Kappa tests to compare different conditions of samples and to assess the inter-examiner reproducibility. RESULTS: Some methods presented adequate reproducibility. The Turesky index and the assessment of area covered by disclosed plaque in the FC images presented the highest discriminatory powers. CONCLUSION: The Turesky index and images with FC with disclosing present good reliability and discriminatory power in quantifying dental plaque.