The effects of numerical resolution on hydrodynamical surface convection simulations and spectral line formation

The computationally demanding nature of radiative-hydrodynamical simulations of stellar surface convection warrants an investigation of the sensitivity of the convective structure and spectral synthesis to the numerical resolution and dimension of the simulations, which is presented here. With too coarse a resolution the predicted spectral lines tend to be too narrow, reflecting insufficient Doppler broadening from the convective motions, while at the currently highest affordable resolution the line shapes have converged essentially perfectly to the observed profiles. Similar conclusions are drawn from the line asymmetries and shifts. In terms of abundances, weak FeI and FeII lines show a very small dependence (~0.02 dex) while for intermediate strong lines with significant non-thermal broadening the sensitivity increases (~0.10 dex). Problems arise when using 2D convection simulations to describe an inherent 3D phenomenon, which translates to inaccurate atmospheric velocity fields and temperature and pressure structures. In 2D the theoretical line profiles tend to be too shallow and broad compared with the 3D calculations and observations, in particular for intermediate strong lines. In terms of abundances, the 2D results are systematically about 0.1 dex lower than for the 3D case for FeI lines. Furthermore, the predicted line asymmetries and shifts are much inferior in 2D. Given these shortcomings and computing time considerations it is better to use 3D simulations of even modest resolution than high-resolution 2D simulations.Comment: Accepted for A&

Realistic Magnetohydrodynamical Simulation of Solar Local Supergranulation

Three-dimensional numerical simulations of solar surface magnetoconvection using realistic model physics are conducted. The thermal structure of convective motions into the upper radiative layers of the photosphere, the main scales of convective cells and the penetration depths of convection are investigated. We take part of the solar photosphere with size of 60x60 Mm in horizontal direction and by depth 20 Mm from level of the visible solar surface. We use a realistic initial model of the Sun and apply equation of state and opacities of stellar matter. The equations of fully compressible radiation magnetohydrodynamics with dynamical viscosity and gravity are solved. We apply: 1) conservative TVD difference scheme for the magnetohydrodynamics, 2) the diffusion approximation for the radiative transfer, 3) dynamical viscosity from subgrid scale modeling. In simulation we take uniform two-dimesional grid in gorizontal plane and nonuniform grid in vertical direction with number of cells 600x600x204. We use 512 processors with distributed memory multiprocessors on supercomputer MVS-100k in the Joint Computational Centre of the Russian Academy of Sciences.Comment: 6 pages, 5 figures, submitted to the proceedings of the GONG 2008 / SOHO XXI conferenc

Relationships Between the Performance of Time/Frequency Standards and Navigation/Communication Systems

The relationship between system performance and clock or oscillator performance is discussed. Tradeoffs discussed include: short term stability versus bandwidth requirements; frequency accuracy versus signal acquisition time; flicker of frequency and drift versus resynchronization time; frequency precision versus communications traffic volume; spectral purity versus bit error rate, and frequency standard stability versus frequency selection and adjustability. The benefits and tradeoffs of using precise frequency and time signals are various levels of precision and accuracy are emphasized

Solar Oscillations and Convection: II. Excitation of Radial Oscillations

Solar p-mode oscillations are excited by the work of stochastic, non-adiabatic, pressure fluctuations on the compressive modes. We evaluate the expression for the radial mode excitation rate derived by Nordlund and Stein (Paper I) using numerical simulations of near surface solar convection. We first apply this expression to the three radial modes of the simulation and obtain good agreement between the predicted excitation rate and the actual mode damping rates as determined from their energies and the widths of their resolved spectral profiles. We then apply this expression for the mode excitation rate to the solar modes and obtain excellent agreement with the low l damping rates determined from GOLF data. Excitation occurs close to the surface, mainly in the intergranular lanes and near the boundaries of granules (where turbulence and radiative cooling are large). The non-adiabatic pressure fluctuations near the surface are produced by small instantaneous local imbalances between the divergence of the radiative and convective fluxes near the solar surface. Below the surface, the non-adiabatic pressure fluctuations are produced primarily by turbulent pressure fluctuations (Reynolds stresses). The frequency dependence of the mode excitation is due to effects of the mode structure and the pressure fluctuation spectrum. Excitation is small at low frequencies due to mode properties -- the mode compression decreases and the mode mass increases at low frequency. Excitation is small at high frequencies due to the pressure fluctuation spectrum -- pressure fluctuations become small at high frequencies because they are due to convection which is a long time scale phenomena compared to the dominant p-mode periods.Comment: Accepted for publication in ApJ (scheduled for Dec 10, 2000 issue). 17 pages, 27 figures, some with reduced resolution -- high resolution versions available at http://www.astro.ku.dk/~aake/astro-ph/0008048

Coupling of Transport and Chemical Processes in Catalytic Combustion

Catalytic combustors have demonstrated the ability to operate efficiently over a much wider range of fuel air ratios than are imposed by the flammability limits of conventional combustors. Extensive commercial use however needs the following: (1) the design of a catalyst with low ignition temperature and high temperature stability, (2) reducing fatigue due to thermal stresses during transient operation, and (3) the development of mathematical models that can be used as design optimization tools to isolate promising operating ranges for the numerous operating parameters. The current program of research involves the development of a two dimensional transient catalytic combustion model and the development of a new catalyst with low temperature light-off and high temperature stablity characteristics

On the Fourier transform of the characteristic functions of domains with C1C^1 -smooth boundary

We consider domains $D\subseteq\mathbb R^n$ with $C^1$ -smooth boundary and study the following question: when the Fourier transform $\hat{1_D}$ of the characteristic function $1_D$ belongs to $L^p(\mathbb R^n)$?Comment: added two references; added footnotes on pages 6 and 1

Computing automorphic forms on Shimura curves over fields with arbitrary class number

We extend methods of Greenberg and the author to compute in the cohomology of a Shimura curve defined over a totally real field with arbitrary class number. Via the Jacquet-Langlands correspondence, we thereby compute systems of Hecke eigenvalues associated to Hilbert modular forms of arbitrary level over a totally real field of odd degree. We conclude with two examples which illustrate the effectiveness of our algorithms.Comment: 15 pages; final submission to ANTS I

Red Blood Cells and Turbulence

Measurements were made of the turbulence intensity of blood of various hematocrits (volume percentage of red cells in blood) flowing through an orifice. The maximum relative turbulence intensity was found to occur in the hematocrit range of 20% - 30%

Magnetoelastic coupling in RETiO3 (RE = La, Nd, Sm, Gd, Y)

A detailed analysis of the crystal structure in RETiO3 with RE = La, Nd, Sm, Gd, and Y reveals an intrinsic coupling between orbital degrees of freedom and the lattice which cannot be fully attributed to the structural deformation arising from bond-length mismatch. The TiO6 octahedra in this series are all irregular with the shape of the distortion depending on the RE ionic radius. These octahedron distortions vary more strongly with temperature than the tilt and rotation angles. Around the Ti magnetic ordering all compounds exhibit strong anomalies in the thermal-expansion coefficients, these anomalies exhibit opposite signs for the antiferromagnetic and ferromagnetic compounds. Furthermore the strongest effects are observed in the materials close to the magnetic cross-over from antiferromagnetic to ferromagnetic order

High-resolution models of solar granulation: the 2D case

Using grid refinement, we have simulated solar granulation in 2D. The refined region measures 1.97*2.58 Mm (vertical*horizontal). Grid spacing there is 1.82*2.84 km. The downflows exhibit strong Kelvin-Helmholtz instabilities. Below the photosphere, acoustic pulses are generated. They proceed laterally (in some cases distances of at least the size of our refined domain) and may be enhanced when transversing downflows) as well as upwards where, in the photosphere they contribute significantly to 'turbulence' (velocity gradients, etc.) The acoustic pulses are ubiquitous in that at any time several of them are seen in our high-resolution domain. Their possible contributions to p-mode excitation or heating of the chromosphere needs to be investigated