Models for the 3-D axisymmetric gravitational potential of the Milky Way Galaxy - A detailed modelling of the Galactic disk

Aims. Galaxy mass models based on simple and analytical functions for the density and potential pairs have been widely proposed in the literature. Disk models constrained by kinematic data alone give information on the global disk structure only very near the Galactic plane. We attempt to circumvent this issue by constructing disk mass models whose three-dimensional structures are constrained by a recent Galactic star counts model in the near-infrared and also by observations of the hydrogen distribution in the disk. Our main aim is to provide models for the gravitational potential of the Galaxy that are fully analytical but also with a more realistic description of the density distribution in the disk component. Methods. From the disk model directly based on the observations (here divided into the thin and thick stellar disks and the HI and H$_2$ disks subcomponents), we produce fitted mass models by combining three Miyamoto-Nagai disk profiles of any "model order" (1, 2, or 3) for each disk subcomponent. The Miyamoto-Nagai disks are combined with models for the bulge and "dark halo" components and the total set of parameters is adjusted by observational kinematic constraints. A model which includes a ring density structure in the disk, beyond the solar Galactic radius, is also investigated. Results. The Galactic mass models return very good matches to the imposed observational constraints. In particular, the model with the ring density structure provides a greater contribution of the disk to the rotational support inside the solar circle. The gravitational potential models and their associated force-fields are described in analytically closed forms, and in addition, they are also compatible with our best knowledge of the stellar and gas distributions in the disk component. The gravitational potential models are suited for investigations of orbits in the Galactic disk.Comment: 22 pages, 13 figures, 11 tables, accepted for publication in A&

An embedded formulation with conforming

Use of strong discontinuities with satisfaction of compatibilit

Arbitrary bi-dimensional finite strain crack propagation

In the past two decades numerous numerical procedures for crack propagation have been developed. Lately, enrichment methods (either local, such as SDA or global, such as XFEM) have been applied with success to simple problems, typically involving some intersections. For arbitrary finite strain propagation, numerous difficulties are encountered: modeling of intersection and coalescence, step size dependence and the presence of distorted finite elements. In order to overcome these difficulties, an approach fully capable of dealing with multiple advancing cracks and self-contact is presented (see Fig.1). This approach makes use of a coupled Arbitrary Lagrangian-Eulerian method (ALE) and local tip remeshing. This is substantially less costly than a full remeshing while retaining its full versatility. Compared to full remeshing, angle measures and crack paths are superior. A consistent continuationbased linear control is used to force the critical tip to be exactly critical, while moving around the candidate set. The critical crack front is identified and propagated when one of the following criteria reaches a material limiting value: (i) the stress intensity factor; or (ii) the element-ahead tip stress. These are the control equations. The ability to solve crack intersection and coalescence problems is shown. Additionally, the independence from crack tip and step size and the absence of blade and dagger-shaped finite elements is observed. Classic benchmarks are computed leading to excellent crack path and load-deflection results, where convergence rate is quadratic