Modelica - A Language for Physical System Modeling, Visualization and Interaction

Modelica is an object-oriented language for modeling of large, complex and heterogeneous physical systems. It is suited for multi-domain modeling, for example for modeling of mechatronics including cars, aircrafts and industrial robots which typically consist of mechanical, electrical and hydraulic subsystems as well as control systems. General equations are used for modeling of the physical phenomena, No particular variable needs to be solved for manually. A Modelica tool will have enough information to do that automatically. The language has been designed to allow tools to generate efficient code automatically. The modeling effort is thus reduced considerably since model components can be reused and tedious and error-prone manual manipulations are not needed. The principles of object-oriented modeling and the details of the Modelica language as well as several examples are presented

Discovery of stars surrounded by iron dust in the LMC

We consider a small sample of oxygen-rich, asymptotic giant branch stars in the Large Magellanic Cloud, observed by the Spitzer Space Telescope, exhibiting a peculiar spectral energy distribution, which can be hardly explained by the common assumption that dust around AGB stars is primarily composed of silicate grains. We suggest that this uncommon class of objects are the progeny of a metal-poor generation of stars, with metallicity $Z \sim 1-2\times 10^{-3}$, formed $\sim 100$ Myr ago. The main dust component in the circumstellar envelope is solid iron. In these stars the poor formation of silicates is set by the strong nucleosynthesis experienced at the base of the envelope, which provokes a scarcity of magnesium atoms and water molecules, required to the silicate formation. The importance of the present results to interpret the data from the incoming James Webb Space Telescope is also discussed.Comment: Accepted for publication in ApJ Letter on 9 January 201

Supernovae data and perturbative deviation from homogeneity

We show that a spherically symmetric perturbation of a dust dominated $\Omega=1$ FRW universe in the Newtonian gauge can lead to an apparent acceleration of standard candles and provide a fit to the magnitude-redshift relation inferred from the supernovae data, while the perturbation in the gravitational potential remains small at all scales. We also demonstrate that the supernovae data does not necessarily imply the presence of some additional non-perturbative contribution by showing that any Lemaitre-Tolman-Bondi model fitting the supernovae data (with appropriate initial conditions) will be equivalent to a perturbed FRW spacetime along the past light cone.Comment: 8 pages, 3 figures; v2: 1 figure added, references added/updated, minor modifications and clarifications, matches published versio

Back-reaction and effective acceleration in generic LTB dust models

We provide a thorough examination of the conditions for the existence of back-reaction and an "effective" acceleration (in the context of Buchert's averaging formalism) in regular generic spherically symmetric Lemaitre-Tolman-Bondi (LTB) dust models. By considering arbitrary spherical comoving domains, we verify rigorously the fulfillment of these conditions expressed in terms of suitable scalar variables that are evaluated at the boundary of every domain. Effective deceleration necessarily occurs in all domains in: (a) the asymptotic radial range of models converging to a FLRW background, (b) the asymptotic time range of non-vacuum hyperbolic models, (c) LTB self-similar solutions and (d) near a simultaneous big bang. Accelerating domains are proven to exist in the following scenarios: (i) central vacuum regions, (ii) central (non-vacuum) density voids, (iii) the intermediate radial range of models converging to a FLRW background, (iv) the asymptotic radial range of models converging to a Minkowski vacuum and (v) domains near and/or intersecting a non-simultaneous big bang. All these scenarios occur in hyperbolic models with negative averaged and local spatial curvature, though scenarios (iv) and (v) are also possible in low density regions of a class of elliptic models in which local spatial curvature is negative but its average is positive. Rough numerical estimates between -0.003 and -0.5 were found for the effective deceleration parameter. While the existence of accelerating domains cannot be ruled out in models converging to an Einstein de Sitter background and in domains undergoing gravitational collapse, the conditions for this are very restrictive. The results obtained may provide important theoretical clues on the effects of back-reaction and averaging in more general non-spherical models.Comment: Final version accepted for publication in Classical and Quantum Gravity. 47 pages in IOP LaTeX macros, 12 pdf figure

Edge Logarithmic Corrections probed by Impurity NMR

Semi-infinite quantum spin chains display spin autocorrelations near the boundary with power-law exponents that are given by boundary conformal field theories. We show that NMR measurements on spinless impurities that break a quantum spin chain lead to a spin-lattice relaxation rate 1/T_1^edge that has a temperature dependence which is a direct probe of the anomalous boundary exponents. For the antiferromagnetic S=1/2 spin chain, we show that 1/T_1^edge behaves as T (log T)^2 instead of (log T)^1/2 for a bulk measurement. We show that, in the case of a one-dimensional conductor described by a Luttinger liquid, a similar measurement leads to a relaxation rate 1/T_1^{edge} behaving as T, independent of the anomalous exponent K_rho.Comment: 4 pages, 1 encapsulated figure, corrected typo

Usefulness of high resolution coastal models for operational oil spill forecast: the "Full City" accident

Oil spill modeling is considered to be an important part of a decision support system (DeSS) for oil spill combatment and is useful for remedial action in case of accidents, as well as for designing the environmental monitoring system that is frequently set up after major accidents. Many accidents take place in coastal areas, implying that low resolution basin scale ocean models are of limited use for predicting the trajectories of an oil spill. In this study, we target the oil spill in connection with the "Full City" accident on the Norwegian south coast and compare operational simulations from three different oil spill models for the area. The result of the analysis is that all models do a satisfactory job. The "standard" operational model for the area is shown to have severe flaws, but by applying ocean forcing data of higher resolution (1.5 km resolution), the model system shows results that compare well with observations. The study also shows that an ensemble of results from the three different models is useful when predicting/analyzing oil spill in coastal areas

Intense Mass Loss from C-rich AGB Stars at low Metallicity?

We argue that the energy injection of pulsations may be of greater importance to the mass-loss rate of AGB stars than metallicity, and that the mass-loss trend with metallicity is not as simple as sometimes assumed. Using our detailed radiation hydrodynamical models that include dust formation, we illustrate the effects of pulsation energy on wind properties. We find that the mass-loss rate scales with the kinetic energy input by pulsations as long as a dust-saturated wind does not occur, and all other stellar parameters are kept constant. This includes the absolute abundance of condensible carbon (not bound in CO), which is more relevant than keeping the C/O-ratio constant when comparing stars of different metallicity. The pressure and temperature gradients in the atmospheres of stars, become steeper and flatter, respectively, when the metallicity is reduced, while the radius where the atmosphere becomes opaque is typically associated with a higher gas pressure. This effect can be compensated for by adjusting the velocity amplitude of the variable inner boundary (piston), which is used to simulate the effects of pulsation, to obtain models with comparable kinetic-energy input. Hence, it is more relevant to compare models with similar energy-injections than of similar velocity amplitude. Since there is no evidence for weaker pulsations in low-metallicity AGB stars, we conclude that it is unlikely that low-metallicity C-stars have a lower mass-loss rate, than their more metal-rich counterparts with similar stellar parameters, as long as they have a comparable amount of condensible carbon.Comment: 4 pages, 3 figures. Accepted for publication in A&A. Updated after language editing. Additional typos fixe

Low-temperature conductivity of quasi-one-dimensional conductors: Luttinger liquid stabilized by impurities

A new non-Fermi-liquid state of quasi-one-dimensional conductors is suggested in which electronic system exists in a form of collection of bounded Luttinger liquids stabilized by impurities. This state is shown to be stable towards interchain electron hopping at low temperatures. Electronic spectrum of the system contains zero modes and collective excitations of the bounded Luttinger liquids in the segments between impurities. Zero modes give rise to randomly distributed localized electronic levels, and long-range interaction generates the Coulomb gap in the density of states at the Fermi energy. Mechanism of conductivity at low temperatures is phonon-assisted hopping via zero-mode states. At higher voltages the excitations of Luttinger liquid are involved in electron transport, and conductivity obeys power-law dependence on voltage. The results provide a qualitative explanation for recent experimental data for NbSe3 and TaS3 crystals.Comment: 12 pages, 1 figur