The Influence of Friction Stir Welding Process Idealization on Residual Stress and Distortion Predictions for Future Airframe Assembly Simulations

The ability to accurately predict residual stresses and resultant distortions is a key product from process assembly simulations. Assembly processes necessarily consider large structural components potentially making simulations computationally expensive. The objective herein is to develop greater understanding of the influence of friction stir welding process idealization on the prediction of residual stress and distortion and thus determine the minimum required modeling fidelity for future airframe assembly simulations. The combined computational and experimental results highlight the importance of accurately representing the welding forging force and process speed. In addition, the results emphasize that increased CPU simulation times are associated with representing the tool torque, while there is potentially only local increase in prediction fidelity.

Tree-code simulations of proton acceleration from laser-irradiated wire targets

Recent experiments using Terawatt lasers to accelerate protons deposited on thin wire targets are modeled with a new type of gridless plasma simulation code. In contrast to conventional mesh-based methods, this technique offers a unique capability in emulating the complex geometry and open-ended boundary conditions characteristic of contemporary experimental conditions. Comparisons of ion acceleration are made between the tree code and standard particle-in-cell simulations for a typical collisionless

The molecular dynamics of Trypanosoma brucei UDP-galactose 4'-epimerase:a drug target for African sleeping sickness

During the past century, several epidemics of human African trypanosomiasis, a deadly disease caused by the protist Trypanosoma brucei, have afflicted sub-Saharan Africa. Over 10 000 new victims are reported each year, with hundreds of thousands more at risk. As current drug treatments are either highly toxic or ineffective, novel trypanocides are urgently needed. The T. brucei galactose synthesis pathway is one potential therapeutic target. Although galactose is essential for T. brucei survival, the parasite lacks the transporters required to intake galactose from the environment. UDP-galactose 4'-epimerase (TbGalE) is responsible for the epimerization of UDP-glucose to UDP-galactose and is therefore of great interest to medicinal chemists. Using molecular dynamics simulations, we investigate the atomistic motions of TbGalE in both the apo and holo states. The sampled conformations and protein dynamics depend not only on the presence of a UDP-sugar ligand, but also on the chirality of the UDP-sugar C4 atom. This dependence provides important insights into TbGalE function and may help guide future computer-aided drug discovery efforts targeting this protein.

Conformational Dynamics of the Human Propeller Telomeric DNA Quadruplex on a Microsecond Time Scale

The human telomeric DNA sequence with four repeats can fold into a parallel-stranded propeller-type topology. NMR structures solved under molecular crowding experiments correlate with the crystal structures found with crystal-packing interactions that are effectively equivalent to molecular crowding. This topology has been used for rationalization of ligand design and occurs experimentally in a number of complexes with a diversity of ligands, at least in the crystalline state. While G-quartet stems have been well characterised, the interactions of the TTA loop with the G-quartets are much less defined. To better understand the conformational variability and structural dynamics of the propeller-type topology, we performed molecular dynamics simulations in explicit solvent up to 1.5 µs. The analysis provides a detailed atomistic account of the dynamic nature of the TTA loops highlighting their interactions with the G-quartets including formation of an A:A base pair, triad, pentad and hexad. The results present a threshold in quadruplex simulations, with regards to understanding the flexible nature of the sugar-phosphate backbone in formation of unusual architecture within the topology. Furthermore, this study stresses the importance of simulation time in sampling conformational space for this topology.

Prospective control of manual interceptive actions:comparative simulations of extant and new model constructs

Two prospective controllers of hand movements in catching-both based on required velocity control-were simulated. Under certain conditions, this required velocity control led to overshoots of the future interception point. These overshoots were absent in pertinent experiments. To remedy this shortcoming, the required velocity model was reformulated in terms of a neural network, the Vector Integration To Endpoint model, to create a Required Velocity Integration To Endpoint model. Addition of a parallel relative velocity channel, resulting in the Relative and Required Velocity Integration To Endpoint model, provided a better account for the experimentally observed kinematics than the existing, purely behavioral models. Simulations of reaching to intercept decelerating and accelerating objects in the presence of background motion were performed to make distinct predictions for future experiments.

Unsteady Simulations of a Film Cooling Flow from an Inclined Cylindrical Jet

Film cooling is extensively used to provide protection against the severe thermal environment in gas turbine engines. Most of the computational studies on film cooling flow have been done using steady Reynolds-averaged Navier–Stokes calculation procedures. However, the flowfield associated with a jet in a crossflow is highly unsteady and complex with different types of vortical structures. In this paper, a computational investigation about the unsteady phenomena of a jet in a crossflow is performed using detached eddy simulation. Detailed computation of a single row of 35 deg round holes on a flat plate has been obtained for a 1.0 blowing ratio and a 2.0 density ratio. First, time-step size, grid resolution, and computational domain tests for an unsteady simulation have been conducted. Comparison between the results of unsteady Reynolds-averaged Navier–Stokes calculation, detached eddy simulation, and large eddy simulation is also performed. Comparison of the time-averaged detached eddy simulation prediction with the measured film-cooling effectiveness shows that the detached eddy simulation prediction is reasonable. From present detached eddy simulations, the influential coherent vortical structures of a film cooling flow can be seen. The unsteady physics of jet in a crossflow interactions and a jet liftoff in film cooling flows have been explained.

Simulations of magnetic microstructure in thin film elements used for programmable motion of magnetic particles

The results of two-dimensional micromagnetic modeling of magnetization patterns in Permalloy ellipses under the influence of rotating constant-amplitude magnetic fields are discussed. Ellipses of two different lateral sizes have been studied, 0.5m x 1.5m and 1m x 3m. The amplitude of the rotating magnetic field was varied between simulations with the condition that it must be large enough to saturate or nearly saturate the ellipse with the field applied along the long axis of the ellipse. For the smaller ellipse size it is found that the magnetization pattern forms an S state and the direction of the net magnetization lags behind the direction of the applied field. At a critical angle of the rotating magnetic field the direction of the magnetization switches by a large angle to a new S state. Both the critical angle and the angle interval of the switch depend on field amplitude. For this new state, it is instead the applied field direction that lags behind the magnetization direction. The transient magnetization patterns correspond to multi-domain patterns including two vortices, but this state never exists for the equilibrated magnetization patterns. The behavior of the larger ellipse in rotating field is different. With the field applied along the long-axis of the ellipse, the magnetization of the ellipse is nearly saturated with a vortex close to each apex of the ellipse. As the field is rotated, this magnetization pattern remains and the net-magnetization direction lags behind the direction of the field until for a certain angle of the applied field an equilibrium multi-domain state is created. Comparisons are made with corresponding experimental results obtained by performing in-field magnetic force microscopy on Permalloy ellipses.

Combined studies of DFT atomistic modelling and in situ FTIR spectroscopy on surface oxidants and CO oxidation at Ru electrodes

We report the combined studies of density functional theory (DFT) calculations and electrochemical in situ FTIR spectroscopy on surface oxidants and mechanisms of CO oxidation at the Ru(0001) electrodes. It is shown that CO can co-adsorb with both O and OH species at lower potential region where a low coverage of the (2 x 2)-O/OH adlayer formed; the oxidation of CO adsorbates takes place at higher potentials where a high coverage of the (1 x 1)-O/OH adlayer formed. Surface O species are not the active oxidants under all coverages studied, due to the high reaction barriers between CO and O (>1 eV). However, surface OH species with higher coverage are identified as the active oxidants, and CO oxidation takes place via a two-steps' mechanism of CO + 3OH -> COOH + 2OH -> CO2 + H2O + OH, in which three nearby OH species are involved in the CO2 formation: CO reacts with OH, forming COOH; COOH then transfers the H to a nearby OH to form H2O and CO2, at the same time, another H in the H2O transfers to a nearby OH to form a weak adsorbed H2O and a new OH. The reaction barrier of these processes is reduced significantly to around 0.50 eV. These new results not only provide an insight into surface active oxidants on Ru, which is directly relevant to fuel cell catalysis, but also reveals the extra complexity of catalytic reactions taking place at solid/liquid electrochemical interface in comparison to the relatively simpler ones at solid/gas phase. 

Two-dimensional Monte Carlo simulations of H i line formation in massive young stellar object disc winds

Massive young stellar objects (YSOs) are powerful infrared Hi line emitters. It has been suggested that these lines form in an outflow from a disc surrounding the YSO. Here, new two-dimensional Monte Carlo radiative transfer calculations are described which test this hypothesis. Infrared spectra are synthesized for a YSO disc wind model based on earlier hydrodynamical calculations. The model spectra are in qualitative agreement with the observed spectra from massive YSOs, and therefore provide support for a disc wind explanation for the Hi lines. However, there are some significant differences: the models tend to overpredict the Bra/Br? ratio of equivalent widths and produce line profiles which are slightly too broad and, in contrast to typical observations, are double-peaked. The interpretation of these differences within the context of the disc wind picture and suggestions for their resolution via modifications to the assumed disc and outflow structure are discussed. © 2005 RAS.

2D simulations of the double-detonation model for thermonuclear transients from low-mass carbon-oxygen white dwarfs

Thermonuclear explosions may arise in binary star systems in which a carbon-oxygen (CO) white dwarf (WD) accretes helium-rich material from a companion star. If the accretion rate allows a sufficiently large mass of helium to accumulate prior to ignition of nuclear burning, the helium surface layer may detonate, giving rise to an astrophysical transient. Detonation of the accreted helium layer generates shock waves that propagate into the underlying CO WD. This might directly ignite a detonation of the CO WD at its surface (an edge-lit secondary detonation) or compress the core of the WD sufficiently to trigger a CO detonation near the centre. If either of these ignition mechanisms works, the two detonations (helium and CO) can then release sufficient energy to completely unbind the WD. These 'double-detonation' scenarios for thermonuclear explosion of WDs have previously been investigated as a potential channel for the production of Type Ia supernovae from WDs of ~ 1 M . Here we extend our 2D studies of the double-detonation model to significantly less massive CO WDs, the explosion of which could produce fainter, more rapidly evolving transients. We investigate the feasibility of triggering a secondary core detonation by shock convergence in low-mass CO WDs and the observable consequences of such a detonation. Our results suggest that core detonation is probable, even for the lowest CO core masses that are likely to be realized in nature. To quantify the observable signatures of core detonation, we compute spectra and light curves for models in which either an edge-lit or compression-triggered CO detonation is assumed to occur. We compare these to synthetic observables for models in which no CO detonation was allowed to occur. If significant shock compression of the CO WD occurs prior to detonation, explosion of the CO WD can produce a sufficiently large mass of radioactive iron-group nuclei to significantly affect the light curves. In particular, this can lead to relatively slow post-maximum decline. If the secondary detonation is edge-lit, however, the CO WD explosion primarily yields intermediate-mass elements that affect the observables more subtly. In this case, near-infrared observations and detailed spectroscopic analysis would be needed to determine whether a core detonation occurred. We comment on the implications of our results for understanding peculiar astrophysical transients including SN 2002bj, SN 2010X and SN 2005E. © 2012 The Authors Monthly Notices of the Royal Astronomical Society © 2012 RAS.