Smooth finite strain plasticity with non-local pressure support

The aim of this work is to introduce an alternative framework to solve problems of finite strain elastoplasticity including anisotropy and kinematic hardening coupled with any isotropic hyperelastic law. After deriving the constitutive equations and inequalities without any of the customary simplifications, we arrive at a new general elasto-plastic system. We integrate the elasto-plastic algebraico-differential system and replace the loading–unloading condition by a Chen–Mangasarian smooth function to obtain a non-linear system solved by a trust region method. Despite being non-standard, this approach is advantageous, since quadratic convergence is always obtained by the non-linear solver and very large steps can be used with negligible effect in the results. Discretized equilibrium is, in contrast with traditional approaches, smooth and well behaved. In addition, since no return mapping algorithm is used, there is no need to use a predictor. The work follows our previous studies of element technology and highly non-linear visco-elasticity. From a general framework, with exact linearization, systematic particularization is made to prototype constitutive models shown as examples. Our element with non-local pressure support is used. Examples illustrating the generality of the method are presented with excellent results

VoroCrust: Voronoi Meshing Without Clipping

Polyhedral meshes are increasingly becoming an attractive option with particular advantages over traditional meshes for certain applications. What has been missing is a robust polyhedral meshing algorithm that can handle broad classes of domains exhibiting arbitrarily curved boundaries and sharp features. In addition, the power of primal-dual mesh pairs, exemplified by Voronoi-Delaunay meshes, has been recognized as an important ingredient in numerous formulations. The VoroCrust algorithm is the first provably-correct algorithm for conforming polyhedral Voronoi meshing for non-convex and non-manifold domains with guarantees on the quality of both surface and volume elements. A robust refinement process estimates a suitable sizing field that enables the careful placement of Voronoi seeds across the surface circumventing the need for clipping and avoiding its many drawbacks. The algorithm has the flexibility of filling the interior by either structured or random samples, while preserving all sharp features in the output mesh. We demonstrate the capabilities of the algorithm on a variety of models and compare against state-of-the-art polyhedral meshing methods based on clipped Voronoi cells establishing the clear advantage of VoroCrust output.Comment: 18 pages (including appendix), 18 figures. Version without compressed images available on https://www.dropbox.com/s/qc6sot1gaujundy/VoroCrust.pdf. Supplemental materials available on https://www.dropbox.com/s/6p72h1e2ivw6kj3/VoroCrust_supplemental_materials.pd

Randomized controlled trial of a home-based action observation intervention to improve walking in Parkinson disease

Published in final edited form as: Arch Phys Med Rehabil. 2016 May ; 97(5): 665–673. doi:10.1016/j.apmr.2015.12.029.OBJECTIVE: To examine the feasibility and efficacy of a home-based gait observation intervention for improving walking in Parkinson disease (PD). DESIGN: Participants were randomly assigned to an intervention or control condition. A baseline walking assessment, a training period at home, and a posttraining assessment were conducted. SETTING: The laboratory and participants' home and community environments. PARTICIPANTS: Nondemented individuals with PD (N=23) experiencing walking difficulty. INTERVENTION: In the gait observation (intervention) condition, participants viewed videos of healthy and parkinsonian gait. In the landscape observation (control) condition, participants viewed videos of moving water. These tasks were completed daily for 8 days. MAIN OUTCOME MEASURES: Spatiotemporal walking variables were assessed using accelerometers in the laboratory (baseline and posttraining assessments) and continuously at home during the training period. Variables included daily activity, walking speed, stride length, stride frequency, leg swing time, and gait asymmetry. Questionnaires including the 39-item Parkinson Disease Questionnaire (PDQ-39) were administered to determine self-reported change in walking, as well as feasibility. RESULTS: At posttraining assessment, only the gait observation group reported significantly improved mobility (PDQ-39). No improvements were seen in accelerometer-derived walking data. Participants found the at-home training tasks and accelerometer feasible to use. CONCLUSIONS: Participants found procedures feasible and reported improved mobility, suggesting that observational training holds promise in the rehabilitation of walking in PD. Observational training alone, however, may not be sufficient to enhance walking in PD. A more challenging and adaptive task, and the use of explicit perceptual learning and practice of actions, may be required to effect change

Finitary Topos for Locally Finite, Causal and Quantal Vacuum Einstein Gravity

Previous work on applications of Abstract Differential Geometry (ADG) to discrete Lorentzian quantum gravity is brought to its categorical climax by organizing the curved finitary spacetime sheaves of quantum causal sets involved therein, on which a finitary (:locally finite), singularity-free, background manifold independent and geometrically prequantized version of the gravitational vacuum Einstein field equations were seen to hold, into a topos structure. This topos is seen to be a finitary instance of both an elementary and a Grothendieck topos, generalizing in a differential geometric setting, as befits ADG, Sorkin's finitary substitutes of continuous spacetime topologies. The paper closes with a thorough discussion of four future routes we could take in order to further develop our topos-theoretic perspective on ADG-gravity along certain categorical trends in current quantum gravity research.Comment: 49 pages, latest updated version (errata corrected, references polished) Submitted to the International Journal of Theoretical Physic

Computational Eulerian Hydrodynamics and Galilean Invariance

Eulerian hydrodynamical simulations are a powerful and popular tool for modeling fluids in astrophysical systems. In this work, we critically examine recent claims that these methods violate Galilean invariance of the Euler equations. We demonstrate that Eulerian hydrodynamics methods do converge to a Galilean-invariant solution, provided a well-defined convergent solution exists. Specifically, we show that numerical diffusion, resulting from diffusion-like terms in the discretized hydrodynamical equations solved by Eulerian methods, accounts for the effects previously identified as evidence for the Galilean non-invariance of these methods. These velocity-dependent diffusive terms lead to different results for different bulk velocities when the spatial resolution of the simulation is kept fixed, but their effect becomes negligible as the resolution of the simulation is increased to obtain a converged solution. In particular, we find that Kelvin-Helmholtz instabilities develop properly in realistic Eulerian calculations regardless of the bulk velocity provided the problem is simulated with sufficient resolution (a factor of 2-4 increase compared to the case without bulk flows for realistic velocities). Our results reiterate that high-resolution Eulerian methods can perform well and obtain a convergent solution, even in the presence of highly supersonic bulk flows.Comment: Version accepted by MNRAS Oct 2, 2009. Figures degraded. For high-resolution color figures and movies of the numerical simulations, please visit http://www.astro.caltech.edu/~brant/Site/Computational_Eulerian_Hydrodynamics_and_Galilean_Invariance.htm