Angular momentum transport modeling: achievements of a gyrokinetic quasi-linear approach

QuaLiKiz, a model based on a local gyrokinetic eigenvalue solver is expanded to include momentum flux modeling in addition to heat and particle fluxes. Essential for accurate momentum flux predictions, the parallel asymmetrization of the eigenfunctions is successfully recovered by an analytical fluid model. This is tested against self-consistent gyrokinetic calculations and allows for a correct prediction of the ExB shear impact on the saturated potential amplitude by means of a mixing length rule. Hence, the effect of the ExB shear is recovered on all the transport channels including the induced residual stress. Including these additions, QuaLiKiz remains ~10 000 faster than non-linear gyrokinetic codes allowing for comparisons with experiments without resorting to high performance computing. The example is given of momentum pinch calculations in NBI modulation experiments

Impact-induced acceleration by obstacles

We explore a surprising phenomenon in which an obstruction accelerates, rather than decelerates, a moving flexible object. It has been claimed that the right kind of discrete chain falling onto a table falls \emph{faster} than a free-falling body. We confirm and quantify this effect, reveal its complicated dependence on angle of incidence, and identify multiple operative mechanisms. Prior theories for direct impact onto flat surfaces, which involve a single constitutive parameter, match our data well if we account for a characteristic delay length that must impinge before the onset of excess acceleration. Our measurements provide a robust determination of this parameter. This supports the possibility of modeling such discrete structures as continuous bodies with a complicated constitutive law of impact that includes angle of incidence as an input.Comment: small changes and corrections, added reference

Automating embedded analysis capabilities and managing software complexity in multiphysics simulation part I: template-based generic programming

An approach for incorporating embedded simulation and analysis capabilities in complex simulation codes through template-based generic programming is presented. This approach relies on templating and operator overloading within the C++ language to transform a given calculation into one that can compute a variety of additional quantities that are necessary for many state-of-the-art simulation and analysis algorithms. An approach for incorporating these ideas into complex simulation codes through general graph-based assembly is also presented. These ideas have been implemented within a set of packages in the Trilinos framework and are demonstrated on a simple problem from chemical engineering

Trust and transactions: Transaction cost analysis with a differential behavioral assumption

Transaction Costs;costs and cost price

MADNESS: A Multiresolution, Adaptive Numerical Environment for Scientific Simulation

MADNESS (multiresolution adaptive numerical environment for scientific simulation) is a high-level software environment for solving integral and differential equations in many dimensions that uses adaptive and fast harmonic analysis methods with guaranteed precision based on multiresolution analysis and separated representations. Underpinning the numerical capabilities is a powerful petascale parallel programming environment that aims to increase both programmer productivity and code scalability. This paper describes the features and capabilities of MADNESS and briefly discusses some current applications in chemistry and several areas of physics

A metabolite-sensitive, thermodynamically-constrained model of\ud cardiac cross-bridge cycling: Implications for force development during ischemia

We present a metabolically regulated model of cardiac active force generation with which we investigate the effects of ischemia on maximum forceproduction. Our model, based on the Rice et al. (2008) model of cross-bridge kinetics, reproduces many of the observed effects of MgATP, MgADP, Pi and H+ on force development while still retaining the force/length/Ca2+ properties of the original model. We introduce three new parameters to account for the competitive binding of H+ to the Ca2+ binding site on troponin C and the binding of MgADP within the cross-bridge cycle. These parameters along with the Pi and H+ regulatory steps within the cross-bridge cycle were constrained using data from the literature and validated using a range of metabolic and sinusoidal length perturbation protocols. The placement of the MgADP binding step between two strongly-bound and force-generating states leads to the emergence of an unexpected effect on the force-MgADP curve, where the trend of the relationship (positive or negative) depends on the concentrations of the other metabolites and [H+]. The model is used to investigate the sensitivity of maximum force production to changes in metabolite concentrations during the development of ischemia