Non-linear boundary element formulation with tangent operator to analyse crack propagation in quasi-brittle materials

This work deals with analysis of cracked structures using BEM. Two formulations to analyse the crack growth process in quasi-brittle materials are discussed. They are based on the dual formulation of BEM where two different integral equations are employed along the opposite sides of the crack surface. The first presented formulation uses the concept of constant operator, in which the corrections of the nonlinear process are made only by applying appropriate tractions along the crack surfaces. The second presented BEM formulation to analyse crack growth problems is an implicit technique based on the use of a consistent tangent operator. This formulation is accurate, stable and always requires much less iterations to reach the equilibrium within a given load increment in comparison with the classical approach. Comparison examples of classical problem of crack growth are shown to illustrate the performance of the two formulations. (C) 2009 Elsevier Ltd. All rights reserved.

Consolidation of elastic-plastic saturated porous media by the boundary element method

This paper presents a domain boundary element formulation for inelastic saturated porous media with rate-independent behavior for the solid skeleton. The formulation is then applied to elastic-plastic behavior for the solid. Biot`s consolidation theory, extended to include irreversible phenomena is considered and the direct boundary element technique is used for the numerical solution after time discretization by the implicit Euler backward algorithm. The associated nonlinear algebraic problem is solved by the Newton-Raphson procedure whereby the loading/unloading conditions are fully taken into account and the consistent tangent operator defined. Only domain nodes (nodes defined inside the domain) are used to represent all domain values and the corresponding integrals are computed by using an accurate sub-elementation scheme. The developments are illustrated through the Drucker-Prager elastic-plastic model for the solid skeleton and various examples are analyzed with the proposed algorithms. (c) 2008 Elsevier B.V. All rights reserved.

A STUDY OF PENALTY FORMULATIONS USED IN THE NUMERICAL APPROXIMATION OF A RADIALLY SYMMETRIC ELASTICITY PROBLEM

We consider a class of two-dimensional problems in classical linear elasticity for which material overlapping occurs in the absence of singularities. Of course, material overlapping is not physically realistic, and one possible way to prevent it uses a constrained minimization theory. In this theory, a minimization problem consists of minimizing the total potential energy of a linear elastic body subject to the constraint that the deformation field must be locally invertible. Here, we use an interior and an exterior penalty formulation of the minimization problem together with both a standard finite element method and classical nonlinear programming techniques to compute the minimizers. We compare both formulations by solving a plane problem numerically in the context of the constrained minimization theory. The problem has a closed-form solution, which is used to validate the numerical results. This solution is regular everywhere, including the boundary. In particular, we show numerical results which indicate that, for a fixed finite element mesh, the sequences of numerical solutions obtained with both the interior and the exterior penalty formulations converge to the same limit function as the penalization is enforced. This limit function yields an approximate deformation field to the plane problem that is locally invertible at all points in the domain. As the mesh is refined, this field converges to the exact solution of the plane problem.

CONCENTRATION FIELDS NEAR AIR-WATER INTERFACES DURING INTERFACIAL MASS TRANSPORT: OXYGEN TRANSPORT AND RANDOM SQUARE WAVE ANALYSIS

Mass transfer across a gas-liquid interface was studied theoretically and experimentally, using transfer of oxygen into water as the gas-liquid system. The experimental results support the conclusions of a theoretical description of the concentration field that uses random square waves approximations. The effect of diffusion over the concentration records was quantified. It is shown that the peak of the normalized rills concentration fluctuation profiles must be lower than 0.5, and that the position of the peak of the rms value is an adequate measure of the thickness of the diffusive layer. The position of the peak is the boundary between the regions more subject to molecular diffusion or to turbulent transport of dissolved mass.

Verification of a mixed high-order accurate DNS code for laminar turbulent transition by the method of manufactured solutions

This paper presents results on a verification test of a Direct Numerical Simulation code of mixed high-order of accuracy using the method of manufactured solutions (MMS). This test is based on the formulation of an analytical solution for the Navier-Stokes equations modified by the addition of a source term. The present numerical code was aimed at simulating the temporal evolution of instability waves in a plane Poiseuille flow. The governing equations were solved in a vorticity-velocity formulation for a two-dimensional incompressible flow. The code employed two different numerical schemes. One used mixed high-order compact and non-compact finite-differences from fourth-order to sixth-order of accuracy. The other scheme used spectral methods instead of finite-difference methods for the streamwise direction, which was periodic. In the present test, particular attention was paid to the boundary conditions of the physical problem of interest. Indeed, the verification procedure using MMS can be more demanding than the often used comparison with Linear Stability Theory. That is particularly because in the latter test no attention is paid to the nonlinear terms. For the present verification test, it was possible to manufacture an analytical solution that reproduced some aspects of an instability wave in a nonlinear stage. Although the results of the verification by MMS for this mixed-order numerical scheme had to be interpreted with care, the test was very useful as it gave confidence that the code was free of programming errors. Copyright (C) 2009 John Wiley & Sons, Ltd.

An EFG method for the nonlinear analysis of plates undergoing arbitrarily large deformations

The applicability of a meshfree approximation method, namely the EFG method, on fully geometrically exact analysis of plates is investigated. Based on a unified nonlinear theory of plates, which allows for arbitrarily large rotations and displacements, a Galerkin approximation via MLS functions is settled. A hybrid method of analysis is proposed, where the solution is obtained by the independent approximation of the generalized internal displacement fields and the generalized boundary tractions. A consistent linearization procedure is performed, resulting in a semi-definite generalized tangent stiffness matrix which, for hyperelastic materials and conservative loadings, is always symmetric (even for configurations far from the generalized equilibrium trajectory). Besides the total Lagrangian formulation, an updated version is also presented, which enables the treatment of rotations beyond the parameterization limit. An extension of the arc-length method that includes the generalized domain displacement fields, the generalized boundary tractions and the load parameter in the constraint equation of the hyper-ellipsis is proposed to solve the resulting nonlinear problem. Extending the hybrid-displacement formulation, a multi-region decomposition is proposed to handle complex geometries. A criterium for the classification of the equilibrium`s stability, based on the Bordered-Hessian matrix analysis, is suggested. Several numerical examples are presented, illustrating the effectiveness of the method. Differently from the standard finite element methods (FEM), the resulting solutions are (arbitrary) smooth generalized displacement and stress fields. (c) 2007 Elsevier Ltd. All rights reserved.

Discrete approximation to the global spectrum of the tangent operator for flow past a circular cylinder

This paper deals with the calculation of the discrete approximation to the full spectrum for the tangent operator for the stability problem of the symmetric flow past a circular cylinder. It is also concerned with the localization of the Hopf bifurcation in laminar flow past a cylinder, when the stationary solution loses stability and often becomes periodic in time. The main problem is to determine the critical Reynolds number for which a pair of eigenvalues crosses the imaginary axis. We thus present a divergence-free method, based on a decoupling of the vector of velocities in the saddle-point system from the vector of pressures, allowing the computation of eigenvalues, from which we can deduce the fundamental frequency of the time-periodic solution. The calculation showed that stability is lost through a symmetry-breaking Hopf bifurcation and that the critical Reynolds number is in agreement with the value presented in reported computations. (c) 2007 IMACS. Published by Elsevier B.V. All rights reserved.

A stable MPC with zone control

Several MPC applications implement a control strategy in which some of the system outputs are controlled within specified ranges or zones, rather than at fixed set points [J.M. Maciejowski, Predictive Control with Constraints, Prentice Hall, New Jersey, 2002]. This means that these outputs will be treated as controlled variables only when the predicted future values lie outside the boundary of their corresponding zones. The zone control is usually implemented by selecting an appropriate weighting matrix for the output error in the control cost function. When an output prediction is inside its zone, the corresponding weight is zeroed, so that the controller ignores this output. When the output prediction lies outside the zone, the error weight is made equal to a specified value and the distance between the output prediction and the boundary of the zone is minimized. The main problem of this approach, as long as stability of the closed loop is concerned, is that each time an output is switched from the status of non-controlled to the status of controlled, or vice versa, a different linear controller is activated. Thus, throughout the continuous operation of the process, the control system keeps switching from one controller to another. Even if a stabilizing control law is developed for each of the control configurations, switching among stable controllers not necessarily produces a stable closed loop system. Here, a stable M PC is developed for the zone control of open-loop stable systems. Focusing on the practical application of the proposed controller, it is assumed that in the control structure of the process system there is an upper optimization layer that defines optimal targets to the system inputs. The performance of the proposed strategy is illustrated by simulation of a subsystem of an industrial FCC system. (C) 2008 Elsevier Ltd. All rights reserved.

On the Filtering Problem for Continuous-Time Markov Jump Linear Systems with no Observation of the Markov Chain

We consider in this paper the optimal stationary dynamic linear filtering problem for continuous-time linear systems subject to Markovian jumps in the parameters (LSMJP) and additive noise (Wiener process). It is assumed that only an output of the system is available and therefore the values of the jump parameter are not accessible. It is a well known fact that in this setting the optimal nonlinear filter is infinite dimensional, which makes the linear filtering a natural numerically, treatable choice. The goal is to design a dynamic linear filter such that the closed loop system is mean square stable and minimizes the stationary expected value of the mean square estimation error. It is shown that an explicit analytical solution to this optimal filtering problem is obtained from the stationary solution associated to a certain Riccati equation. It is also shown that the problem can be formulated using a linear matrix inequalities (LMI) approach, which can be extended to consider convex polytopic uncertainties on the parameters of the possible modes of operation of the system and on the transition rate matrix of the Markov process. As far as the authors are aware of this is the first time that this stationary filtering problem (exact and robust versions) for LSMJP with no knowledge of the Markov jump parameters is considered in the literature. Finally, we illustrate the results with an example.

Study of the interaction between hydroxymethylnitrofurazone and 2-hydroxypropyl-beta-cyclodextrin

Chagas disease is a serious health problem in Latin America. Hidroxymethylnitrofurazone (NFOH) is a nitrofurazone prodrug more active than nitrofurazone against Trypanosoma cruzi. However, NFOH presents low aqueous solubility, high photodecomposition and high toxicity. The present work is focused on the characterization of an inclusion complex of NFOH in 2-hydroxypropyl-beta-cyclodextrin (HP-beta-CD). The complexation with HP-beta-CD was investigated using reversed-phase liquid chromatography, solubility isotherms and nuclear magnetic resonance. The retention behavior was analyzed on a reversed-phase C-18 column, using acetonitrile-water (20/80, v/v) as the mobile phase, in which HP-beta-CD was incorporated as a mobile phase additive. The decrease in the retention times with increasing concentrations of HP-beta-CD enables the determination of the apparent stability constant of the complex (K = 6.2 +/- 0.3 M-1) by HPLC. The solubility isotherm was studied and the value for the apparent stability constant (K = 7.9 +/- 0.2 M-1) was calculated. The application of continuous variation method indicated the presence of a complex with 1:1 NFOH:HP-beta-CD stoichiometry. The photostability study showed that the formation of an inclusion complex had a destabilizing effect on the photodecomposition of NFOH when compared to that of the ""free"" molecule in solution. The mobility investigation (by NMR longitudinal relaxation time) gives information about the complexation of NFOH with HP-beta-CD. In preliminary toxicity studies, cell viability tests revealed that inclusion complexes were able to decrease the toxic effect (p < 0.01) caused by NFOH. (c) 2008 Elsevier B.V. All rights reserved.

QUADRANT-SEARCHING: A NEW TECHNIQUE TO REDUCE THE SEARCH SPACE OF THE INVERSE PROBLEM OF EIT USING BEM TO THE DIRECT PROBLEM

The image reconstruction using the EIT (Electrical Impedance Tomography) technique is a nonlinear and ill-posed inverse problem which demands a powerful direct or iterative method. A typical approach for solving the problem is to minimize an error functional using an iterative method. In this case, an initial solution close enough to the global minimum is mandatory to ensure the convergence to the correct minimum in an appropriate time interval. The aim of this paper is to present a new, simple and low cost technique (quadrant-searching) to reduce the search space and consequently to obtain an initial solution of the inverse problem of EIT. This technique calculates the error functional for four different contrast distributions placing a large prospective inclusion in the four quadrants of the domain. Comparing the four values of the error functional it is possible to get conclusions about the internal electric contrast. For this purpose, initially we performed tests to assess the accuracy of the BEM (Boundary Element Method) when applied to the direct problem of the EIT and to verify the behavior of error functional surface in the search space. Finally, numerical tests have been performed to verify the new technique.

Prognostic value of left atrial volume in patients who underwent dobutamine stress echocardiography for known or suspected coronary artery disease

Background Left atrial volume indexed (LAVI) has been reported as a predictor of cardiovascular events. We sought to determine the prognostic value of LAVI for predicting the outcome of patients who underwent dobutamine stress echocardiography (DSE) for known or suspected coronary artery disease (CAD). Methods From January 2000 to July 2005, we studied 981 patients who underwent DSE and off-line measurements of LAVI. The value of DSE over clinical and LAVI data was examined using a stepwise log-rank test. Results During a median follow-up of 24 months, 56 (6%) events occurred. By univariate analysis, predictors of events were male sex, diabetes mellitus, previous myocardial infarction, left ventricular ejection fraction (LVEF), left atrial diameter indexed, LAVI, and abnormal DSE. By multivariate analysis, independent predictors were LVEF (relative risk [RR] = 0.98, 95% CI 0.95-1.00), LAVI (RR = 1.04, 95% CI 1.02-1.05), and abnormal DSE (RR = 2.70, 95% CI 1.28-5.69). In an incremental multivariate model, LAVI was additional to clinical data for predicting events (chi(2) 36.8, P < .001). The addition of DSE to clinical and LAVI yielded incremental information (chi(2) 55.3, P < .001). The 3-year event-free survival in patients with normal DSE and LAVI <= 33 mL/m(2) was 96%; with abnormal DSE and LAVI <= 33 mL/m(2), 91%; with normal DSE and LAVI >34 mL/m(2), 83%; and with abnormal DSE and LAVI >34 mL/m(2) 51%. Conclusion Left atrial volume indexed provides independent prognostic information in patients who underwent DSE for known or suspected CAD. Among patients with normal DSE, those with larger LAVI had worse outcome, and among patients with abnormal DSE, LAVI was still predictive. (Am Heart J 2008; 156:1110-6.)

Prognostic value of NDRG1 and SPARC protein expression in breast cancer patients

An increasing number of studies have shown altered expression of secreted protein acidic and rich in cysteine (SPARC) and N-myc down-regulated gene (NDRG1) in several malignancies, including breast carcinoma; however, the role of these potential biomarkers in tumor development and progression is controversial. In this study, NDRG1 and SPARC protein expression was evaluated by immunohistochemistry on tissue microarrays containing breast tumor specimens from patients with 10 years of follow-up. NDRG1 and SPARC protein expression was determined in 596 patients along with other prognostic markers, such as ER, PR, and HER2. The status of NDRG1 and SPARC protein expression was correlated with prognostic variables and patient clinical outcome. Immunostaining revealed that 272 of the 596 cases (45.6%) were positive for NDRG1 and 431 (72.3%) were positive for SPARC. Statistically significant differences were found between the presence of SPARC and NDRG1 protein expression and standard clinicopathological variables. Kaplan-Meier analysis showed that NDRG1 positivity was directly associated with shorter disease-free survival (DFS, P < 0.001) and overall survival (OS, P < 0.001). In contrast, patients expressing low levels of SPARC protein had worse DFS (P = 0.001) and OS (P = 0.001) compared to those expressing high levels. Combined analysis of the two markers indicated that DFS (P < 0.001) and OS rates (P < 0.001) were lowest for patients with NDRG1-positive and SPARC-negative tumors. Furthermore, NDRG1 over-expression and SPARC down-regulation correlated with poor prognosis in patients with luminal A or triple-negative subtype breast cancer. On multivariate analysis using a Cox proportional hazards model, NDRG1 and SPARC protein expression were independent prognostic factors for both DFS and OS of breast cancer patients. These data indicate that NDRG1 over-expression and SPARC down-regulation could play important roles in breast cancer progression and serve as useful biomarkers to better define breast cancer prognosis.

(18)F-FDG PET After 2 Cycles of ABVD Predicts Event-Free Survival in Early and Advanced Hodgkin Lymphoma

Our objective was to assess the prognostic value of (18)F-FDG PET after 2 cycles of chemotherapy using doxorubicin, bleomycin, vinblastine, and dacarbazine (ABVD) in Hodgkin lymphoma (HL) patients overall and in subgroups of patients with early and advanced stages and with low and high risks according to the International Prognostic Score (IPS). Methods: One hundred fifteen patients with newly diagnosed HL were prospectively included in the study. All underwent standard ABVD therapy followed by consolidation radiotherapy in cases of bulky disease. After 2 cycles of ABVD, the patients were evaluated with PET (PET2). Prognostic analysis compared the 3-y event-free survival (EFS) rate to the PET2 results, clinical data, and IPS. Results: Of the 104 evaluated patients, 93 achieved complete remission after first-line therapy. During a median follow-up of 36 mo, relapse or disease progression was seen in 22 patients. Treatment failure was seen in 16 of the 30 PET2-positive patients and in only 6 of the 74 PET2-negative patients. PET2 was the only significant prognostic factor. The 3-y EFS was 53.4% for PET2-positive patients and 90.5% for PET2-negative ones (P < 0.001). When patients were categorized according to low or high IPS risk and according to early or advanced stage of disease, PET2 was also significantly associated with treatment outcome. Conclusion: PET2 is an accurate and independent predictor of EFS in HL. A negative interim (18)F-FDG PET result is highly predictive of treatment success in overall HL patients, as well as in subgroups with early or advanced-stage disease and with low or high IPS risk.

Extent of Pelvic Lymphadenectomy in Women With Squamous Cell Carcinoma of the Uterine Cervix: Is There Any Prognostic Value?

Background and Objectives: Some authors states that the removal of lymph node would only contribute towards assessing the lymph node status and regional disease control, without any benefit for the patients` survival. The aim of this paper was to assess the influence of the number of surgically dissected pelvic lymph nodes (PLN) on disease-free Survival. Methods: Retrospective cohort study on 42 women presenting squamous cell carcinoma (SCC) of the uterine cervix, with metastases in PLN treated by radical surgery. The Cox model was used to identify risk factors for recurrence. The model variables were adjusted for treatment-related factors (year of treatment, surgical margins and postoperative radiotherapy). The cutoff value for classifying the lymphadenectomy as comprehensive (15 PLN or more) or non-comprehensive (<15 PLN) was determined from analysis of the ROC curve. Results: Fourteen recurrences (32.6%) were recorded: three pelvic, eight distant, two both pelvic and distant, and one at an unknown location. The following risk factors for recurrence were identified: invasion of the deep third of the cervix and number of dissected lymph nodes <15. Conclusions: Deep invasion and non-comprehensive pelvic lymphadenectomy are possible risk factors for recurrence of SCC of the uterine cervix with metastases in PLN. J. Surg. Oncol. 2009;100:252-257. (C) 2009 Wiley-Liss, Inc.