Manufacturing Algebra. Part II: aggregation, control and simulation

Manufacturing Algebra provides a set of mathematical entities together with some composition rules, that are specially conceived for modelling and controlling a manufacturing system. Here only the modelling capabilities are outlined together with a simple case study. The scope of the paper is to familiarize the reader with the proposed methodology, and to highlight some peculiarities. Formulation is reduced to a minimum. Among the algebra peculiarities, both manufacturing process and the factory layout are defined in their basic elements, and the link between them is given. Specifically the Manufacturing Model (parts, operations) includes time and space coordinates in order it could be employed by factory elements like Production units and Control Units. This calls for the definition of event and event sequence and of the relevant discrete-event elements and operators. A further peculiarity to be clarified in the second part, is the capability of aggregating algebra elements into higher level components, thus favoring hierarchical description and control of manufacturing system

A streamwise-constant model of turbulent pipe flow

A streamwise-constant model is presented to investigate the basic mechanisms responsible for the change in mean flow occuring during pipe flow transition. Using a single forced momentum balance equation, we show that the shape of the velocity profile is robust to changes in the forcing profile and that both linear non-normal and nonlinear effects are required to capture the change in mean flow associated with transition to turbulence. The particularly simple form of the model allows for the study of the momentum transfer directly by inspection of the equations. The distribution of the high- and low-speed streaks over the cross-section of the pipe produced by our model is remarkably similar to one observed in the velocity field near the trailing edge of the puff structures present in pipe flow transition. Under stochastic forcing, the model exhibits a quasi-periodic self-sustaining cycle characterized by the creation and subsequent decay of "streamwise-constant puffs", so-called due to the good agreement between the temporal evolution of their velocity field and the projection of the velocity field associated with three-dimensional puffs in a frame of reference moving at the bulk velocity. We establish that the flow dynamics are relatively insensitive to the regeneration mechanisms invoked to produce near-wall streamwise vortices and that using small, unstructured background disturbances to regenerate the streamwise vortices is sufficient to capture the formation of the high- and low-speed streaks and their segregation leading to the blunting of the velocity profile characteristic of turbulent pipe flow

A Comparison of Measured Crab and Vela Glitch Healing Parameters with Predictions of Neutron Star Models

There are currently two well-accepted models that explain how pulsars exhibit glitches, sudden changes in their regular rotational spin-down. According to the starquake model, the glitch healing parameter, Q, which is measurable in some cases from pulsar timing, should be equal to the ratio of the moment of inertia of the superfluid core of a neutron star (NS) to its total moment of inertia. Measured values of the healing parameter from pulsar glitches can therefore be used in combination with realistic NS structure models as one test of the feasibility of the starquake model as a glitch mechanism. We have constructed NS models using seven representative equations of state of superdense matter to test whether starquakes can account for glitches observed in the Crab and Vela pulsars, for which the most extensive and accurate glitch data are available. We also present a compilation of all measured values of Q for Crab and Vela glitches to date which have been separately published in the literature. We have computed the fractional core moment of inertia for stellar models covering a range of NS masses and find that for stable NSs in the realistic mass range 1.4 +/- 0.2 solar masses, the fraction is greater than 0.55 in all cases. This range is not consistent with the observational restriction Q < 0.2 for Vela if starquakes are the cause of its glitches. This confirms results of previous studies of the Vela pulsar which have suggested that starquakes are not a feasible mechanism for Vela glitches. The much larger values of Q observed for Crab glitches (Q > 0.7) are consistent with the starquake model predictions and support previous conclusions that starquakes can be the cause of Crab glitches.Comment: 8 pages, including 3 figures and 1 table. Accepted for publication in Ap

Linear and Nonlinear Evolution and Diffusion Layer Selection in Electrokinetic Instability

In the present work fournontrivial stages of electrokinetic instability are identified by direct numerical simulation (DNS) of the full Nernst-Planck-Poisson-Stokes (NPPS) system: i) The stage of the influence of the initial conditions (milliseconds); ii) 1D self-similar evolution (milliseconds-seconds); iii) The primary instability of the self-similar solution (seconds); iv) The nonlinear stage with secondary instabilities. The self-similar character of evolution at intermediately large times is confirmed. Rubinstein and Zaltzman instability and noise-driven nonlinear evolution to over-limiting regimes in ion-exchange membranes are numerically simulated and compared with theoretical and experimental predictions. The primary instability which happens during this stage is found to arrest self-similar growth of the diffusion layer and specifies its characteristic length as was first experimentally predicted by Yossifon and Chang (PRL 101, 254501 (2008)). A novel principle for the characteristic wave number selection from the broadbanded initial noise is established.Comment: 13 pages, 8 figure

Holographic Magnetic Star

A warm fermionic AdS star under a homogeneous magnetic field is explored. We obtain the relativistic Landau levels by using Dirac equation and use the Tolman-Oppenheimer-Volkoff (TOV) equation to study the physical profiles of the star. Bulk properties such as sound speed, adiabatic index, and entropy density within the star are calculated analytically and numerically. Bulk temperature increases the mass limit of the AdS star but external magnetic field has the opposite effect. The results are partially interpreted in terms of the pre-thermalization process of the gauge matter at the AdS boundary after the mass injection. The entropy density is found to demonstrate similar temperature dependence as the magnetic black brane in the AdS in certain limits regardless of the different nature of the bulk and Hawking temperatures. Total entropy of the AdS star is also found to be an increasing function of the bulk temperature and a decreasing function of the magnetic field, similar behaviour to the mass limit. Since both total entropy and mass limit are global quantities, they could provide some hints to the value of entropy and energy of the dual gauge matter before and during the thermalization.Comment: 39 pages, 14 figures, 1 table, comments and references added, to appear in JHE

Equation of state and opacities for hydrogen atmospheres of magnetars

The equation of state and radiative opacities of partially ionized, strongly magnetized hydrogen plasmas, presented in a previous paper [ApJ 585, 955 (2003), astro-ph/0212062] for the magnetic field strengths 8.e11 G < B < 3.e13 G, are extended to the field strengths 3.e13 G < B < 1.e15 G, relevant for magnetars. The first- and second-order thermodynamic functions and radiative opacities are calculated and tabulated for 5.e5 < T < 4.e7 K in a wide range of densities. We show that bound-free transitions give an important contribution to the opacities in the considered range of B in the outer neutron-star atmosphere layers. Unlike the case of weaker fields, bound-bound transitions are unimportant.Comment: 7 pages, 6 figures, LaTeX using emulateapj.cls (included). Accepted by Ap

A Spectral Method for Elliptic Equations: The Neumann Problem

Let $\Omega$ be an open, simply connected, and bounded region in $\mathbb{R}^{d}$, $d\geq2$, and assume its boundary $\partial\Omega$ is smooth. Consider solving an elliptic partial differential equation $-\Delta u+\gamma u=f$ over $\Omega$ with a Neumann boundary condition. The problem is converted to an equivalent elliptic problem over the unit ball $B$, and then a spectral Galerkin method is used to create a convergent sequence of multivariate polynomials $u_{n}$ of degree $\leq n$ that is convergent to $u$. The transformation from $\Omega$ to $B$ requires a special analytical calculation for its implementation. With sufficiently smooth problem parameters, the method is shown to be rapidly convergent. For $u\in C^{\infty}(\overline{\Omega})$ and assuming $\partial\Omega$ is a $C^{\infty}$ boundary, the convergence of $\Vert u-u_{n}\Vert_{H^{1}}$ to zero is faster than any power of $1/n$. Numerical examples in $\mathbb{R}^{2}$ and $\mathbb{R}^{3}$ show experimentally an exponential rate of convergence.Comment: 23 pages, 11 figure

Numerical simulations of two dimensional magnetic domain patterns

I show that a model for the interaction of magnetic domains that includes a short range ferromagnetic and a long range dipolar anti-ferromagnetic interaction reproduces very well many characteristic features of two-dimensional magnetic domain patterns. In particular bubble and stripe phases are obtained, along with polygonal and labyrinthine morphologies. In addition, two puzzling phenomena, namely the so called `memory effect' and the `topological melting' observed experimentally are also qualitatively described. Very similar phenomenology is found in the case in which the model is changed to be represented by the Swift-Hohenberg equation driven by an external orienting field.Comment: 8 pages, 8 figures. Version to appear in Phys. Rev.