A Smooth Finite Element Method Based on Reproducing Kernel DMS-Splines

The element-based piecewise smooth functional approximation in the conventional finite element method (FEM) results in discontinuous first and higher order derivatives across element boundaries Despite the significant advantages of the FEM in modelling complicated geometries, a motivation in developing mesh-free methods has been the ease with which higher order globally smooth shape functions can be derived via the reproduction of polynomials There is thus a case for combining these advantages in a so-called hybrid scheme or a `smooth FEM' that, whilst retaining the popular mesh-based discretization, obtains shape functions with uniform C-p (p >= 1) continuity One such recent attempt, a NURBS based parametric bridging method (Shaw et al 2008b), uses polynomial reproducing, tensor-product non-uniform rational B-splines (NURBS) over a typical FE mesh and relies upon a (possibly piecewise) bijective geometric map between the physical domain and a rectangular (cuboidal) parametric domain The present work aims at a significant extension and improvement of this concept by replacing NURBS with DMS-splines (say, of degree n > 0) that are defined over triangles and provide Cn-1 continuity across the triangle edges This relieves the need for a geometric map that could precipitate ill-conditioning of the discretized equations Delaunay triangulation is used to discretize the physical domain and shape functions are constructed via the polynomial reproduction condition, which quite remarkably relieves the solution of its sensitive dependence on the selected knotsets Derivatives of shape functions are also constructed based on the principle of reproduction of derivatives of polynomials (Shaw and Roy 2008a) Within the present scheme, the triangles also serve as background integration cells in weak formulations thereby overcoming non-conformability issues Numerical examples involving the evaluation of derivatives of targeted functions up to the fourth order and applications of the method to a few boundary value problems of general interest in solid mechanics over (non-simply connected) bounded domains in 2D are presented towards the end of the paper

Dilute solution properties of methyl methacrylate—acrylonitrile copolymer (MA2)

This article deals with studies of the dilute solution properties of methyl methacrylate-acrylonitrile (MMA-AN) copolymer of 0.415 mole fraction (mf) of acrylonitrile composition. Mark—Houwink parameters for this copolymer have been evaluated in acetonitrile (MeCN), 2-butanone (MEK), dimethyl formamide (DMF), and γ-butyrolactone (γ-BL). The Mark-Houwink exponent a in all four solvents at all temperatures is larger than the corresponding values of the parent homopolymers. The solvent power is in the order of DMF < γ–BL < MEK < MeCN; [η] decreases with an increase in temperature, which is behavior characteristic of polymers in good solvent. The unperturbed dimensions (K0) values, obtained by the Stockmayer–Fixman method, are lower than those for the parent homopolymers and depend on solvent as well as temperature. The solute—solvent interaction parameter X1 values are close to 0.5; X1 is independent of temperature. The excess interaction parameter XABvalues are negative. The results for this copolymer system in regard to low second virial coefficient A2, large X1, and high a values suggest that the large extension of these copolymer chains is due to the unusual short-range interactions.

Comparative study of different formulations and solution methods for two phase flow in saturated porous media

Six models (Simulators) are formulated and developed with all possible combinations of pressure and saturation of the phases as primary variables. A comparative study between six simulators with two numerical methods, conventional simultaneous and modified sequential methods are carried out. The results of the numerical models are compared with the laboratory experimental results to study the accuracy of the model especially in heterogeneous porous media. From the study it is observed that the simulator using pressure and saturation of the wetting fluid (PW, SW formulation) is the best among the models tested. Many simulators with nonwetting phase as one of the primary variables did not converge when used along with simultaneous method. Based on simulator 1 (PW, SW formulation), a comparison of different solution methods such as simultaneous method, modified sequential and adaptive solution modified sequential method are carried out on 4 test problems including heterogeneous and randomly heterogeneous problems. It is found that the modified sequential and adaptive solution modified sequential methods could save the memory by half and as also the CPU time required by these methods is very less when compared with that using simultaneous method. It is also found that the simulator with PNW and PW as the primary variable which had problem of convergence using the simultaneous method, converged using both the modified sequential method and also using adaptive solution modified sequential method. The present study indicates that pressure and saturation formulation along with adaptive solution modified sequential method is the best among the different simulators and methods tested.

Fractal boundaries of basin of attraction of Newton-Raphson method in helicopter trim

A key problem in helicopter aeroelastic analysis is the enormous computational time required for a numerical solution of the nonlinear system of algebraic equations required for trim, particularly when free wake models are used. Trim requires calculation of the main rotor and tail rotor controls and the vehicle attitude which leads to the six steady forces and moments about the helicopter center of gravity to be zero. An appropriate initial estimate of the trim state is needed for successful helicopter trim. This study aims to determine the control inputs that can have considerable effect on the convergence of trim solution in the aeroelastic analysis of helicopter rotors by investigating the basin of attraction of the nonlinear equations (set of initial guess points from which the nonlinear equations converge). It is illustrated that the three main rotor pitch controls of collective pitch, longitudinal cyclic pitch and lateral cyclic pitch have a significant contribution to the convergence of the trim solution. Trajectories of the Newton iterates are shown and some ideas for accelerating the convergence of a trim solution in the aeroelastic analysis of helicopters are proposed. It is found that the basins of attraction can have fractal boundaries. (C) 2010 Elsevier Ltd. All rights reserved.

Direct Growth of Gold Nanorods on Gold and Indium Tin Oxide Surfaces: Spectral, Electrochemical, and Electrocatalytic Studies

The surfactant-assisted seed-mediated growth method was used for the formation of gold nanorods (GNRs) directly on gold (Au) and indium tin oxide (ITO) surfaces. Citrate-stabilized similar to 2.6 nm spherical gold nanoparticles (AuNPs) were first self-assembled on ITO or Au surfaces modified with (3-mercaptopropyl)-trimethoxysilane (MPTS) sol-gel film and then immersed in a cationic surfactant growth solution to form GNRs. The growth of GNRs on the MPTS sol gel film modified ITO surface was monitored by UV-visible spectroscopy. The ITO surface with the attached spherical AuNPs shows a surface plasmon resonance band at 550 nm. The intensity of this absorption band increases while increasing the immersion time of the AuNP-modified ITO surface into the growth solution, and after 5 h, an additional shoulder band around 680 nm was observed. The intensity of this shoulder band increases, and it was shifted to longer wavelength as the immersion time of the AuNP-modified ITO surface into the growth solution increases. After 20 h, a predominant wave at 720 nm was observed along with a band at 550 nm. Further immersion of the modified ITO surface into the growth solution did not change the absorption characteristics. The bands observed at 550 and 720 nm were characteristics of GNRs, corresponding to transverse and longitudinal waves, respectively. The AFM images showed the presence of GNRs on the surface of the MPTS sol gel modified ITO surface with a typical length of similar to 100-120 nm and a width of similar to 20-22 nm in addition to a few spherical AuNPs, indicating that seeded spherical AuNPs were not completely involved in the GNRs' formation. Finally, the electrocatalytic activity of the surface-grown GNRs on the MPTS sol gel film modified Au electrode toward the oxidation of ascorbic acid (AA) was studied. Unlike a polycrystalline Au electrode, the surface-grown GNR-modified electrode shows two well-defined voltammetric peaks for AA at 0.01 and 0.35 V in alkaline, neutral, and acidic pHs. The cause for the observed two oxidation peaks for AA was due to the presence of both nanorods and spherical nanoparticles on the electrode surface. The presence of spherical AuNPs on the MPTS sol gel film oxidized AA at more positive potential, whereas the GNRs oxidized AA at less positive potential. The observed 340 mV less positive potential shift in the oxidation of AA suggested that GNRs are better electrocatalysts for the oxidation of AA than the spherical AuNPs.

Nonnegative integral solution of linear equations

A method to obtain a nonnegative integral solution of a system of linear equations, if such a solution exists is given. The method writes linear equations as an integer programming problem and then solves the problem using a combination of artificial basis technique and a method of integer forms.

An improved iterative finite element solution for pin joints

Accurate, reliable and economical methods of determining stress distributions are important for fastener joints. In the past the contact stress problems in these mechanically fastened joints using interference or push or clearance fit pins were solved using both inverse and iterative techniques. Inverse techniques were found to be most efficient, but at times inadequate in the presence of asymmetries. Iterative techniques based on the finite element method of analysis have wider applications, but they have the major drawbacks of being expensive and time-consuming. In this paper an improved finite element technique for iteration is presented to overcome these drawbacks. The improved iterative technique employs a frontal solver for elimination of variables not requiring iteration, by creation of a dummy element. This automatically results in a large reduction in computer time and in the size of the problem to be handled during iteration. Numerical results are compared with those available in the literature. The method is used to study an eccentrically located pin in a quasi-isotropic laminated plate under uniform tension.

A New Improved Procedure for the Synthesis of 4-Formylcarane - Precursor to Sandatrile, a Novel Perfumery Chemical

Stereospecific synthesis of 4-formylcarane (2) has been achieved through hydroboration-carbonylation of DELTA-3-carene. Both the reactions are optimised using sodium borohydride. The method is utilised for the synthesis of sandatrile (3), a novel perfumery chemical.

Dislocation characterization in crystals of potash alum grown by seeded solution growth under conditions of low supersaturation

An analysis has been carried out of the genesis and character of growth dislocations present in all growth sectors of single crystals of potash alum. The crystals, grown from seeded solutions by the temperature lowering method under conditions of low supersaturation, presented the well-developed forms: {111} dominant, {100} and {110}. Growth dislocations formed predominately during refacetting of the edges and corners of the seed, rounded during preparation and insertion into the supersaturated solution. From here they become refracted into the {111} sectors which proved to be the most defective. Smaller numbers of dislocations form at the {111}, {100} and {110} seed interfaces and propagate in these sectors. In crystals of inferior quality, a number of inclusions were found predominantly in the fast growing {100} sectors which become the source of additional dislocations. Dislocations present in the original seed did not propagate across the interface into the developing crystal. Dislocations of all characters were observed. The principal Burgers vectors were found to be left angle bracket100right-pointing angle bracket, left angle bracket110right-pointing angle bracket and left angle bracket111right-pointing angle bracket.

Dynamics of activationless reactions in solution

Recent picosecond and subpicosecond laser spectroscopy experiments have revealed several chemically and biologically important reactions in solution in which the reaction potential surface does not present a barrier to the motion along the reaction coordinate.The dynamics of these reactions display diverse and interesting behavior. They include the dependence of relaxation rate on the solvent viscosity, the solvent polarity, the temperature, and the wavelength of the exciting light. In this article we review the recent developments in the theoretical description of activationless processes in solution and compare them with the available experimental results

A parametric method 3 (PM3) study of trans‐stilbene

We report a comparative modified neglect of diatomic overlap (MNDO), Austin method one (AM1), and parametric method 3 (PM3) study of trans‐stilbene (tS) in its ground, excited (singlet and triplet), and ionic (positive and negative polarons and bipolarons) states. We have also calculated the barrier for ring rotation about the backbone single bond. Our results show that PM3 geometries are superior to MNDO and AM1, at least for tS. PM3 predicts, in contrast with MNDO, AM1 and even ab initio 3‐21G, a coplanar structure for tS, in accordance with recent experimental data. Singlet and triplet energies obtained from heats of formation are in surprisingly good agreement with experimental data.

Lipid-amphotericin B complex structure in solution: A possible first step in the aggregation process in cell membranes

The interactions between the polyene antibiotic amphotericin B with dipalmitoylphosphatidylcholine were investigated in vesicles (using circular dichroism) and in chloroform solution (using circular dichroism and IH, I3C, and 31P nuclear magnetic resonance). The results show that amphotericin B readily aggregates in vesicles and that the extent of aggregation depends on the 1ipid:drug concentration ratio. Introduction of sterol molecules into the membrane hastens the process of aggregation of amphotericin B. In chloroform solutions amphotericin B strongly interacts with phospholipid molecules to form a stoichiometric complex. The results suggest that there are interactions between the conjugated heptene stretch of amphotericin B and the methylene groups of lipid acyl chains, while the sugar moiety interacts with the phosphate head group by the formation of a hydrogen bond. A model is proposed for the lipid-amphotericin B complex, in which amphotericin B interacts equally well with the two lipid acyl chains, forming a 1:l complex.

Microwave Irradiation Assisted Method for the Rapid Synthesis of fine Particles of alpha Al2O3 And alpha (Al1-XCrx)(2)O-3 And Their Coatings on Si(100)

Chromium substituted beta diketonate complexes of aluminium have been synthesized and employed as precursors for a novel soft chemistry process wherein microwave irradiation of a solution of the complex yields within minutes well crystallized needles of alpha (Al1 XCrx)(2)O-3 measuring 20 30 nm in diameter and 50 nm long By varying the microwave irradiation parameters and using a surfactant such as polyvinyl pyrrolidone the crystallite size and shape can be controlled and their agglomeration prevented These microstructural parameters as well as the polymorph of the Cr substituted Al2O3 formed may also be controlled by employing a different complex Samples of alpha (Al1 XCrx)(2)O-3 have been characterized by XRD FTIR and TEM The technique results in material of homogeneous metal composition, as shown by EDAX and can be adjusted as desired The technique has been extended to obtain coatings of alpha (Al1 XCrx)(2)O-3 on Si(100)

Synthesis, Reactivity, Structural Aspects, and Solution Dynamics of Cyclopalladated Compounds of N,N `,N `'-Tris(2-anisyl)guanidine

N,N',N `'-Tris(2-anisyl)guanidine, (ArNH)(2)C=NAr (Ar = 2-(MeO)C6H4), was cyclopallaclated with Pd(OC(O)R)(2) (R = Me, CF3) in toluene at 70 degrees C to afford palladacycles Pd{kappa(2)(C,N)-C6H3-(OMe)-3(NHC(NHAr)(=NAr))-2}(mu-OC(O)R)](2)(R = Me (1a) and CF3 (1b)) in 87% and 95% yield, respectively. Palladacycle 1a was subjected to a metathetical reaction with LiBr in aqueous ethanol at 78 degrees C to afford palladacycle Pd{kappa(2)(C,N)-C6H3(OMe)-3(NHC(NHAr)(=NAr))-2}(mu-Br)](2) (2) in 90% yield. Palladacycle 2 was subjected to a bridge-splitting reaction with Lewis bases in CH2Cl2 to afford the monomeric palladacycles Pd{kappa(2)(C,N)-C6H3(OMe)-3(NHC(NHAr)(=NAr))-2}Br(L)] (L = 2,6-Me2C5H3N (3a), 2,4-Me2C5H3N (3b), 3,5-Me2C5H3N (3c), XyNC (Xy = 2,6-Me2C6H3; 4a), (BuNC)-Bu-t (4b), and PPh3 (5)) in 87-95% yield. Palladacycle 2 upon reaction with 2 equiv of XyNC in CH2Cl2 afforded an unanticipated palladacycle, Pd{kappa(2)(C,N)-C(=NXy)(C6H3(OMe)-4)-2(N=C-(NH Ar)(2))-3} Br(CNXy)] (6) in 93% yield, and the driving force for the formation of 6 was ascribed to a ring contraction followed by amine-imine tautomerization. Palladacycles 1 a,b revealed a dimeric transoid in-in conformation with ``open book'' framework in the solid state. In solution, 1 a exhibited a fluxional behavior ascribed to the six-membered ``(C,N)Pd'' ring inversion and partly dissociates to the pincer type and kappa(2)-O,O'-OAc monomeric palladacycles by an anchimerically assisted acetate cleavage process as studied by variable-temperature H-1 NMR data. Palladacycles 3a,b revealed a unique trans configuration around the palladium with lutidine being placed trans to the Pd-C bond, whereas cis stereochemistry was observed between the Pd-C bond and the Lewis base in 4a (as determined by X-ray diffraction data) and 5 (as determined by P-31 and C-13 NMR data). The aforementioned stereochemical difference was explained by invoking relative hardness/softness of the donor atoms around the palladium center. In solution, palladacycles 3a-c exist as a mixture of two interconverting boat conformers via a planar intermediate without any bond breaking due to the six-membered ``(C,N)Pd'' ring inversion, whereas palladacycles 4a,b and 5 exist as a single isomer, as deduced from detailed H-1 NMR studies.

Geometry and quadratic nonlinearity of charge transfer complexes in solution using depolarized hyper-Rayleigh scattering

We report large quadratic nonlinearity in a series of 1:1 molecular complexes between methyl substituted benzene donors and quinone acceptors in solution. The first hyperpolarizability, beta(HRS), which is very small for the individual components, becomes large by intermolecular charge transfer (CT) interaction between the donor and the acceptor in the complex. In addition, we have investigated the geometry of these CT complexes in solution using polarization resolved hyper-Rayleigh scattering (HRS). Using linearly (electric field vector along X direction) and circularly polarized incident light, respectively, we have measured two macroscopic depolarization ratios D = I-2 omega,I-X,I-X/I-2 omega,I-Z,I-X and D' = I-2 omega,I-X,I-C/I-2 omega,I-Z,I-C in the laboratory fixed XYZ frame by detecting the second harmonic scattered light in a polarization resolved fashion. The experimentally obtained first hyperpolarizability, beta(HRS), and the value of macroscopic depolarization ratios, D and D', are then matched with the theoretically deduced values from single and double configuration interaction calculations performed using the Zerner's intermediate neglect of differential overlap self-consistent reaction field technique. In solution, since several geometries are possible, we have carried out calculations by rotating the acceptor moiety around three different axes keeping the donor molecule fixed at an optimized geometry. These rotations give us the theoretical beta(HRS), D and D' values as a function of the geometry of the complex. The calculated beta(HRS), D, and D' values that closely match with the experimental values, give the dominant equilibrium geometry in solution. All the CT complexes between methyl benzenes and chloranil or 1,2-dichloro-4,5-dicyano-p-benzoquinone investigated here are found to have a slipped parallel stacking of the donors and the acceptors. Furthermore, the geometries are staggered and in some pairs, a twist angle as high as 30 degrees is observed. Thus, we have demonstrated in this paper that the polarization resolved HRS technique along with theoretical calculations can unravel the geometry of CT complexes in solution. (C) 2011 American Institute of Physics. doi:10.1063/1.3514922]