A Staggered Explicit-Implicit Finite Element Formulation for Electroactive Polymers

Electroactive polymers such as dielectric elastomers (DEs) have attracted significant attention in recent years. Computational techniques to solve the coupled electromechanical system of equations for this class of materials have universally centered around fully coupled monolithic formulations, which while generating good accuracy requires significant computational expense. However, this has significantly hindered the ability to solve large scale, fully three-dimensional problems involving complex deformations and electromechanical instabilities of DEs. In this work, we provide theoretical basis for the effectiveness and accuracy of staggered explicit-implicit finite element formulations for this class of electromechanically coupled materials, and elicit the simplicity of the resulting staggered formulation. We demonstrate the stability and accuracy of the staggered approach by solving complex electromechanically coupled problems involving electroactive polymers, where we focus on problems involving electromechanical instabilities such as creasing, wrinkling, and bursting drops. In all examples, essentially identical results to the fully monolithic solution are obtained, showing the accuracy of the staggered approach at a significantly reduced computational cost

Metal-insulator (fermion-boson)-crossover origin of pseudogap phase of cuprates I: anomalous heat conductivity, insulator resistivity boundary, nonlinear entropy

Among all experimental observations of cuprate physics, the metal-insulator-crossover (MIC), seen in the pseudogap (PG) region of the temperature-doping phase diagram of copper-oxides under a strong magnetic field, when the superconductivity is suppressed, is most likely the most intriguing one. Since it was expected that the PG-normal state for these materials, as for conventional superconductors, is conducting. This MIC, revealed in such phenomena as heat conductivity downturn, anomalous Lorentz ratio, insulator resistivity boundary, nonlinear entropy, resistivity temperature upturn, insulating ground state, nematicity- and stripe-phases and Fermi pockets, unambiguously indicates on the insulating normal state, from which the high-temperature superconductivity (HTS) appears. In the present work (article I), we discuss the MIC phenomena mentioned in the title of article. The second work (article II) will be devoted to discussion of other listed above MIC phenomena and also to interpretation of the recent observations in the hidden magnetic order and scanning tunneling microscopy (STM) experiments spin and charge fluctuations as the intra PG and HTS pair ones. We find that all these MIC (called in the literature as non-Fermi liquid) phenomena can be obtained within the Coulomb single boson and single fermion two liquid model, which we recently developed, and the MIC is a crossover of single fermions into those of single bosons. We show that this MIC originates from bosons of Coulomb two liquid model and fermions, whose origin is these bosons. At an increase of doping up to critical value or temperature up to PG boundary temperature, the boson system undegoes bosonic insulator - bosonic metal - fermionic metal transitions.Comment: 13 pages, 3 figure

Implementation of a partitioned algorithm for simulation of large CSI problems

The implementation of a partitioned numerical algorithm for determining the dynamic response of coupled structure/controller/estimator finite-dimensional systems is reviewed. The partitioned approach leads to a set of coupled first and second-order linear differential equations which are numerically integrated with extrapolation and implicit step methods. The present software implementation, ACSIS, utilizes parallel processing techniques at various levels to optimize performance on a shared-memory concurrent/vector processing system. A general procedure for the design of controller and filter gains is also implemented, which utilizes the vibration characteristics of the structure to be solved. Also presented are: example problems; a user's guide to the software; the procedures and algorithm scripts; a stability analysis for the algorithm; and the source code for the parallel implementation

Heat transfer and pressure drop in blade cooling channels with turbulence promoters

Repeated rib roughness elements have been used in advanced turbine cooling designs to enhance the internal heat transfer. Often the ribs are perpendicular to the main flow direction so that they have an angle-of-attack of 90 deg. The objective of the project was to investigate the effect of rib angle-of-attack on the pressure drop and the average heat transfer coefficients in a square duct with two opposite rib-roughned walls for Reynolds number varied from 8000 to 80,000. The rib height-to-equivalent diameter ratio (e/D) was kept at a constant value of 0.063, the rib pitch-to-height ratio (P/e) was varied from 10 to 20, and the rib angle-of-attack (alpha) was varied from 90 deg to 60 deg to 45 deg to 30 deg respectively. Two types of entrance conditions were examined, namely, long duct and sudden contraction. The heat transfer coefficient distribution on the smooth side wall and the rough side wall at the entrance and the fully developed regions were measured. Thermal performance comparison indicated that the pumping power requirement for the rib with an oblique angle to the flow (alpha = 45 deg to 30 deg) was about 20 to 50 percent lower than the rib with a 90 deg angle to the flow for a given heat transfer duty

Computational methods and software systems for dynamics and control of large space structures

Two key areas of crucial importance to the computer-based simulation of large space structures are discussed. The first area involves multibody dynamics (MBD) of flexible space structures, with applications directed to deployment, construction, and maneuvering. The second area deals with advanced software systems, with emphasis on parallel processing. The latest research thrust in the second area involves massively parallel computers

A dynamical mean-field theory study of Nagaoka ferromagnetism

We revisit Nagaoka ferromagnetism in the U=oo Hubbard model within the dynamical mean-field theory (DMFT) using the recently developed continuous time quantum Monte Carlo method as the impurity solver. The stability of Nagaoka ferromagnetism is studied as a function of the temperature, the doping level, and the next-nearest-neighbor lattice hopping t'. We found that the nature of the phase transition as well as the stability of the ferromagnetic state is very sensitive to the t' hopping. Negative t'=-0.1t stabilizes ferromagnetism up to higher doping levels. The paramagnetic state is reached through a first order phase transition. Alternatively, a second order phase transition is observed at t'=0. Very near half-filling, the coherence temperature T_{coh} of the paramagnetic metal becomes very low and ferromagnetism evolves out of an incoherent metal rather than conventional Fermi liquid. We use the DMFT results to benchmark slave-boson method which might be useful in more complicated geometries.Comment: 10 pages, 11 figure

The interaction of small metal particles with refractory oxide supports

Islands and continuous layers of Pd were grown in UHV on Mo and MoO subtrates. As-deposited Pd islands and layers exhibited bulk Pd adsorption properties for CO when the Pd had been deposited at RT and at thicknesses exceeding 3 ML. However, CO adsorption was drastically reduced upon annealing. This deactivation was interpreted in terms of substrate/support interaction involving the diffusion of substrate species toward the Pd surface, using AES, TPD, and work function measurement techniques. A study of the growth and annealing behavior of Pd on Mo(110) was made for thicknesses up to 12 monolayers and substrate temperatures up to 1300K, using AES, XPS, LEED, and work function measurements. At low tempertures Pd formed a monolayer without alloying. In thick layers (12 ML) annealed about 700 K, Mo diffusion into the Pd layer and alloying were noted. Such layers remained continuous up to 1100 K. Thinner Pd layers were less stable and started coalescing upon annealing to as little as 550 K. Significant changes in Pd Auger peak shape, as well as shifts of Pd core levels, were observed during layer growth and annealing