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.

Power Secant Method applied to natural frequency extraction of Timoshenko beam structures

This work deals with an improved plane frame formulation whose exact dynamic stiffness matrix (DSM) presents, uniquely, null determinant for the natural frequencies. In comparison with the classical DSM, the formulation herein presented has some major advantages: local mode shapes are preserved in the formulation so that, for any positive frequency, the DSM will never be ill-conditioned; in the absence of poles, it is possible to employ the secant method in order to have a more computationally efficient eigenvalue extraction procedure. Applying the procedure to the more general case of Timoshenko beams, we introduce a new technique, named ""power deflation"", that makes the secant method suitable for the transcendental nonlinear eigenvalue problems based on the improved DSM. In order to avoid overflow occurrences that can hinder the secant method iterations, limiting frequencies are formulated, with scaling also applied to the eigenvalue problem. Comparisons with results available in the literature demonstrate the strength of the proposed method. Computational efficiency is compared with solutions obtained both by FEM and by the Wittrick-Williams algorithm.

An oracle approach for interaction neighborhood estimation in random fields

We consider the problem of interaction neighborhood estimation from the partial observation of a finite number of realizations of a random field. We introduce a model selection rule to choose estimators of conditional probabilities among natural candidates. Our main result is an oracle inequality satisfied by the resulting estimator. We use then this selection rule in a two-step procedure to evaluate the interacting neighborhoods. The selection rule selects a small prior set of possible interacting points and a cutting step remove from this prior set the irrelevant points. We also prove that the Ising models satisfy the assumptions of the main theorems, without restrictions on the temperature, on the structure of the interacting graph or on the range of the interactions. It provides therefore a large class of applications for our results. We give a computationally efficient procedure in these models. We finally show the practical efficiency of our approach in a simulation study.

Efficiency issues of evolutionary k-means

One of the top ten most influential data mining algorithms, k-means, is known for being simple and scalable. However, it is sensitive to initialization of prototypes and requires that the number of clusters be specified in advance. This paper shows that evolutionary techniques conceived to guide the application of k-means can be more computationally efficient than systematic (i.e., repetitive) approaches that try to get around the above-mentioned drawbacks by repeatedly running the algorithm from different configurations for the number of clusters and initial positions of prototypes. To do so, a modified version of a (k-means based) fast evolutionary algorithm for clustering is employed. Theoretical complexity analyses for the systematic and evolutionary algorithms under interest are provided. Computational experiments and statistical analyses of the results are presented for artificial and text mining data sets. (C) 2010 Elsevier B.V. All rights reserved.

Bayesian network classifiers: Beyond classification accuracy

This work proposes and discusses an approach for inducing Bayesian classifiers aimed at balancing the tradeoff between the precise probability estimates produced by time consuming unrestricted Bayesian networks and the computational efficiency of Naive Bayes (NB) classifiers. The proposed approach is based on the fundamental principles of the Heuristic Search Bayesian network learning. The Markov Blanket concept, as well as a proposed ""approximate Markov Blanket"" are used to reduce the number of nodes that form the Bayesian network to be induced from data. Consequently, the usually high computational cost of the heuristic search learning algorithms can be lessened, while Bayesian network structures better than NB can be achieved. The resulting algorithms, called DMBC (Dynamic Markov Blanket Classifier) and A-DMBC (Approximate DMBC), are empirically assessed in twelve domains that illustrate scenarios of particular interest. The obtained results are compared with NB and Tree Augmented Network (TAN) classifiers, and confinn that both proposed algorithms can provide good classification accuracies and better probability estimates than NB and TAN, while being more computationally efficient than the widely used K2 Algorithm.

Efficient synthesis of benzothiazine and acrylamide compounds

This article describes the synthesis of the new (2Z)-2-(4-methoxybenzylidene)-6-nitro-4H -benzo[1,4]thiazin-3-one, (2Z)-2-(4-methoxybenzylidene)-4-methyl-6-nitro-4H-benzo[1,4]thiazin-3-one, (2Z)-6-amino-2-(4-methoxybenzylidene)-4H -benzo[1,4]thiazin-3-one, (2Z)-6-butylamino-2-(4-methoxybenzylidene)-4-methyl-4H-benzo[1,4]-thiazin-3-one and (2E)-N-alkyl-N-(2-hydroxy-5-nitrophenyl)-3-phenylacrylamides and the spectroscopic data. The arylidenebenzothiazine compounds were prepared using the Knoevenagel condensation with substituted benzaldehydes in the presence of sodium methoxide in DMF. The presence of a nitro substituent in the 4-position, water and a slightly acid reaction medium in this condensation caused the rupture of the benzothiazine ring and subsequent formation of the phenylacrylamide compounds. A crystallographic data was presented for (2E)-3-(4-bromophenyl)-N-dodecyl-N -(2-hydroxy-5-nitrophenyl) acrylamide.

Africanized honey bees are efficient at detecting, uncapping and removing dead brood

The hygienic behavior of honey bees is based on a two-step process, including uncapping and removing diseased, dead, damaged, or parasitized brood inside the cell. We evaluated during periods of 1 h the time that hygienic and non-hygienic colonies of Africanized honey bees spend to detect, uncap and remove pin-killed brood using comb inserts with transparent walls placed in observation hives. We observed that hygienic colonies are significantly faster in detecting, uncapping and removing dead brood in the cells (P < 0.001).

Determination of the Most Efficient Temperature for Radiofrequency Ablation of Renal Cells: A Prospective Study in Dogs

Background and Purpose: Radiofrequency (RF) ablation of renal tumors is a major technique for tumor cell destruction while preserving healthy renal parenchyma. There is no consensus in the literature regarding the optimal temperature, impedance, and time for RF application for effective cell destruction. This study investigated two variables while keeping time unchanged: Temperature for RF cell destruction and tissue impedance in dog kidneys. Materials and Methods: Sixteen dogs had renal punctures through videolaparoscopy for RF interstitial tissue ablation. A RF generator was applied for 10 minutes to the dog's kidney at different target temperatures: 80 degrees C, 90 degrees C, and 100 degrees C. On postoperative day14, the animals were sacrificed and nephrectomized. All lesions were macroscopically and microscopically examined. The bioelectrical impedance was evaluated at three different temperatures. Results: Renal injuries were wider and deeper at 90 degrees C (P < 0.001), and they were similar at 80 degrees C and 100 degrees C. The bioelectrical impedance was lower at 90 degrees C than at the temperatures of 80 degrees C and 100 degrees C (P < 0.001). Viable cells in the RF ablation tissue area were not found in the microscopic examination. Conclusion: The most effective cell destruction in terms of width and depth was achieved at 90 degrees C, which was also the optimal temperature for tissue impedance. RF ablation of renal cells eliminated all viable cells.

Organic molecule assembled between carbon nanotubes: A highly efficient switch device

We report a highly efficient switch built from an organic molecule assembled between single-wall carbon nanotube electrodes. We theoretically show that changes in the distance between the electrodes alter the molecular conformation within the gap, affecting in a dramatic way the electronic and charge transport properties, with an on/off ratio larger than 300. This opens up the perspective of combining molecular electronics with carbon nanotubes, bringing great possibilities for the design of nanodevices.

Lossless filter for multiple repeats with bounded edit distance

Background: Identifying local similarity between two or more sequences, or identifying repeats occurring at least twice in a sequence, is an essential part in the analysis of biological sequences and of their phylogenetic relationship. Finding such fragments while allowing for a certain number of insertions, deletions, and substitutions, is however known to be a computationally expensive task, and consequently exact methods can usually not be applied in practice. Results: The filter TUIUIU that we introduce in this paper provides a possible solution to this problem. It can be used as a preprocessing step to any multiple alignment or repeats inference method, eliminating a possibly large fraction of the input that is guaranteed not to contain any approximate repeat. It consists in the verification of several strong necessary conditions that can be checked in a fast way. We implemented three versions of the filter. The first is simply a straightforward extension to the case of multiple sequences of an application of conditions already existing in the literature. The second uses a stronger condition which, as our results show, enable to filter sensibly more with negligible (if any) additional time. The third version uses an additional condition and pushes the sensibility of the filter even further with a non negligible additional time in many circumstances; our experiments show that it is particularly useful with large error rates. The latter version was applied as a preprocessing of a multiple alignment tool, obtaining an overall time (filter plus alignment) on average 63 and at best 530 times smaller than before (direct alignment), with in most cases a better quality alignment. Conclusion: To the best of our knowledge, TUIUIU is the first filter designed for multiple repeats and for dealing with error rates greater than 10% of the repeats length.

omega-Transaminases as efficient biocatalysts to obtain novel chiral selenium-amine ligands for Pd-catalysis

omega-Transaminases have been evaluated as biocatalysts in the reductive amination of organoselenium acetophenones to the corresponding amines, and in the kinetic resolution of racemic organoselenium amines. Kinetic resolution proved to be more efficient than the asymmetric reductive amination. By using these methodologies we were able to obtain both amine enantiomers in high enantiomeric excess (up to 99%). Derivatives of the obtained optically pure o-selenium 1-phenylethyl amine were evaluated as ligands in the palladium-catalyzed asymmetric alkylation, giving the alkylated product in up to 99% ee.

Gold is not a Faradaic-Efficient Borohydride Oxidation Electrocatalyst: An Online Electrochemical Mass Spectrometry Study

Direct borohydride fuel cells are promising high energy density portable generators. However, their development remains limited by the complexity of the anodic reaction: The borohydride oxidation reaction (BOR) kinetics is slow and occurs at high overvoltages, while it may compete with the heterogeneous hydrolysis of BH(4)(-). Nevertheless, one usually admits that gold is rather inactive toward the heterogeneous hydrolysis of BH(4)(-) and presents some activity regarding the BOR, therefore yielding to the complete eight-electron BOR. In the present paper, by coupling online mass spectrometry to electrochemistry, we in situ monitored the H(2) yield during BOR experiments on sputtered gold electrodes. Our results show non-negligible H(2) generation on Au on the whole BOR potential range (0-0.8 V vs reversible hydrogen electrode), thus revealing that gold cannot be considered as a faradaic-efficient BOR electrocatalyst. We further propose a relevant reaction pathway for the BOR on gold that accounts for these findings.

Mechanistic studies of the 'blue' Cu enzyme, bilirubin oxidase, as a highly efficient electrocatalyst for the oxygen reduction reaction

The 'blue copper' enzyme bilirubin oxidase from Myrothecium verrucaria shows significantly enhanced adsorption on a pyrolytic graphite 'edge' (PGE) electrode that has been covalently modified with naphthyl-2-carboxylate functionalities by diazonium coupling. Modified electrodes coated with bilirubin oxidase show electrocatalytic voltammograms for the direct, four-electron reduction of O(2) by bilirubin oxidase with up to four times the current density of an unmodified PGE electrode. Electrocatalytic voltammograms measured with a rapidly rotating electrode (to remove effects of O(2) diffusion limitation) have a complex shape (an almost linear dependence of current on potential below pH 6) that is similar regardless of how PGE is chemically modified. Importantly, the same waveform is observed if bilirubin oxidase is adsorbed on Au(111) or Pt(111) single-crystal electrodes (at which activity is short-lived). The electrocatalytic behavior of bilirubin oxidase, including its enhanced response on chemically-modified PGE, therefore reflects inherent properties that do not depend on the electrode material. The variation of voltammetric waveshapes and potential-dependent (O(2)) Michaelis constants with pH and analysis in terms of the dispersion model are consistent with a change in rate-determining step over the pH range 5-8: at pH 5, the high activity is limited by the rate of interfacial redox cycling of the Type 1 copper whereas at pH 8 activity is much lower and a sigmoidal shape is approached, showing that interfacial electron transfer is no longer a limiting factor. The electrocatalytic activity of bilirubin oxidase on Pt(111) appears as a prominent pre-wave to electrocatalysis by Pt surface atoms, thus substantiating in a single, direct experiment that the minimum overpotential required for O(2) reduction by the enzyme is substantially smaller than required at Pt. At pH 8, the onset of O(2) reduction lies within 0.14 V of the four-electron O(2)/2H(2)O potential.

Efficient method for obtaining Lep(ob)/Lep(ob)-derived animal models using adipose tissue transplantations

Background: Leptin-deficient mice (Lep(ob)/Lep(ob), also known as ob/ob) are of great importance for studies of obesity, diabetes and other correlated pathologies. Thus, generation of animals carrying the Lep(ob) gene mutation as well as additional genomic modifications has been used to associate genes with metabolic diseases. However, the infertility of Lep(ob)/Lep(ob) mice impairs this kind of breeding experiment. Objective: To propose a new method for production of Lep(ob)/Lep(ob) animals and Lep(ob)/Lep(ob)-derived animal models by restoring the fertility of Lep(ob)/Lep(ob) mice in a stable way through white adipose tissue transplantations. Methods: For this purpose, 1 g of peri-gonadal adipose tissue from lean donors was used in subcutaneous transplantations of Lep(ob)/Lep(ob) animals and a crossing strategy was established to generate Lep(ob)/Lep(ob)-derived mice. Results: The presented method reduced by four times the number of animals used to generate double transgenic models (from about 20 to 5 animals per double mutant produced) and minimized the number of genotyping steps (from 3 to 1 genotyping step, reducing the number of Lep gene genotyping assays from 83 to 6). Conclusion: The application of the adipose transplantation technique drastically improves both the production of Lep(ob)/Lep(ob) animals and the generation of Lep(ob)/Lep(ob)-derived animal models. International Journal of Obesity (2009) 33, 938-944; doi: 10.1038/ijo.2009.95; published online 16 June 2009

An efficient heuristic for total flowtime minimisation in no-wait flowshops

In this paper, we address the problem of scheduling jobs in a no-wait flowshop with the objective of minimising the total completion time. This problem is well-known for being nondeterministic polynomial-time hard, and therefore, most contributions to the topic focus on developing algorithms able to obtain good approximate solutions for the problem in a short CPU time. More specifically, there are various constructive heuristics available for the problem [such as the ones by Rajendran and Chaudhuri (Nav Res Logist 37: 695-705, 1990); Bertolissi (J Mater Process Technol 107: 459-465, 2000), Aldowaisan and Allahverdi (Omega 32: 345-352, 2004) and the Chins heuristic by Fink and Voa (Eur J Operat Res 151: 400-414, 2003)], as well as a successful local search procedure (Pilot-1-Chins). We propose a new constructive heuristic based on an analogy with the two-machine problem in order to select the candidate to be appended in the partial schedule. The myopic behaviour of the heuristic is tempered by exploring the neighbourhood of the so-obtained partial schedules. The computational results indicate that the proposed heuristic outperforms existing ones in terms of quality of the solution obtained and equals the performance of the time-consuming Pilot-1-Chins.