We examine how cosmological expansion arises in a universe containing a lattice of spinning black holes. We study averaged expansion properties as a function of fundamental properties of the black holes, including the bare mass of the black holes and black hole spin. We then explore how closely the expansion properties correspond to properties of a corresponding matter-dominated FLRW universe. As residual radiation present in the initial data decays, we find good agreement with a matter-dominated FLRW solution, and the effective density in the volume is well-described by the horizon mass of the black hole.

Standard cosmological models rely on an approximate treatment of gravity, utilizing solutions of the linearized Einstein equations as well as physical approximations. In an era of precision cosmology, we should ask: are these approximate predictions sufficiently accurate for comparison to observations, and can we draw meaningful conclusions about properties of our Universe from them? In this work we examine the accuracy of linearized gravity in the presence of collisionless matter and a cosmological constant utilizing fully general relativistic simulations. We observe the gauge-dependence of corrections to linear theory, and note the amplitude of these corrections. For perturbations whose amplitudes are in line with expectations from the standard $\Lambda$CDM model, we find that the full, general relativistic metric is well-described by linear theory in Newtonian and harmonic gauges, while the metric in comoving-synchronous gauge is not. For the largest observed structures in our Universe, our results suggest that corrections to linear gravitational theory can reach or surpass the percent-level.

Using two-dimensional simulations, we numerically explore the dependences of Kelvin-Helmholtz instability upon various physical parameters, including viscosity, width of sheared layer, flow speed, and magnetic field strength. In most cases, a multi-vortex phase exists between the initial growth phase and final single-vortex phase. The parametric study shows that the evolutionary properties, such as phase duration and vortex dynamics, are generally sensitive to these parameters except in certain regimes. An interesting result is that for supersonic flows, the phase durations and saturation of velocity growth approach constant values asymptotically as the sonic Mach number increases. We confirm that the linear coupling between magnetic field and Kelvin-Helmholtz modes is negligible if the magnetic field is weak enough. The morphological behaviour suggests that the multi-vortex coalescence might be driven by the underlying wave-wave interaction. Based on these results, we make a preliminary discussion about several events observed in the solar corona. The numerical models need to be further improved to make a practical diagnostic of the coronal plasma properties.

An attempt at understanding downward overshooting in the convective envelopes of post-main-sequence stars has been made on the basis of three-dimensional large-eddy simulations, using artificially modified OPAL opacity and taking into account radiation and ionization in the equation of state. Two types of star, an intermediate-mass star and a massive star, were considered. To avoid a long thermal relaxation time of the intermediate-mass star, we increased the stellar energy flux artificially while trying to maintain a structure close to the one given by a 1D stellar model. A parametric study of the flux factor was performed. For the massive star, no such process was necessary. Numerical results were analysed when the system reached the statistical steady state. It was shown that the penetration distance in pressure scaleheights is of the order of unity. The scaling relations between penetration distance, input flux and vertical velocity fluctuations studied by Singh et al. were checked. The anisotropy of the turbulent convection and the diffusion models of the third-order moments representing the non-local transport were also investigated. These models are dramatically affected by the velocity fields and no universal constant parameters seem to exist. The limitations of the numerical results were also discussed.

Structures seen in idealized numerical experiments on compressible magnetoconvection in an imposed strong vertical magnetic field show important differences from those detected in observations or realistic numerical simulations of sunspot umbrae. To elucidate the origin of these discrepancies, we present a series of idealized 3D compressible magnetoconvection experiments that differ from previous such experiments in several details, bringing them closer to realistic solar conditions. An initially vertical magnetic field B 0 is imposed on a time snapshot of fully developed solar-like turbulent convection in a layer bounded by a stable layer from above. Upon relaxation to a statistically steady state, the structure of the flow field and magnetic field is examined. Instead of the vigorous granular convection (GRC) well known to take place in magnetized or weakly magnetized convection, for high values of B0 heat is transported by small-scale convection (SSC) in the form of narrow, persistent convective columns consisting of slender upflows accompanied by adjacent downflow patches, which are reminiscent of the 'convectons' identified in earlier semianalytic models. For moderate field strengths, flux separation (FXS) is observed: isolated field-free inclusions of GRC are embedded in a strongly magnetized plasma with SSC. Between the SSC and FXS regimes, a transitional regime (F/S) is identified where convectons dynamically evolve into multiply segmented granular inclusions and back. Our results agree in some aspects more closely with observed umbral structures than earlier idealized models, because they do reproduce the strong localized, patchy downflows immediately adjacent to the narrow convective columns. Based on recent observations of umbral dots, we suggest that in some cases the conditions in sunspot umbrae correspond to the newly identified F/S transitional regime.

We systematically study the renormalizable three-term polynomial inflation in the supersymmetric and non-supersymmetric models. The supersymmetric inflaton potentials can be realized in supergravity theory, and only have two independent parameters. We show that the general renormalizable supergravity model is equivalent to one kind of our supersymmetric models. We find that the spectral index and tensor-to-scalar ratio can be consistent with the Planck and BICEP2 results, but the running of spectral index is always out of the $2\sigma$ range. If we do not consider the BICEP2 experiment, these inflationary models can be highly consistent with the Planck observations and saturate its upper bound on the tensor-to-scalar ratio ($r \le 0.11$). Thus, our models can be tested at the future Planck and QUBIC experiments.

The multidimensional gas-kinetic scheme for the Navier-Stokes equations under gravitational fields [J. Comput. Phys. 226 (2007) 2003-2027] is extended to resistive magnetic flows. The non-magnetic part of the magnetohydrodynamics equations is calculated by a BGK solver modified due to magnetic field. The magnetic part is treated by the flux splitting method based gas-kinetic theory [J. Comput. Phys. 153 (1999) 334-352 ], using a particle distribution function constructed in the BGK solver. To include Lorentz force effects into gas evolution stage is very important to improve the accuracy of the scheme. For some multidimensional problems, the deviations tangential to the cell interface from equilibrium distribution are essential to keep the scheme robust and accurate. Besides implementation of a TVD time discretization scheme, enhancing the dynamic dissipation a little bit is a simply and efficient way to stabilize the calculation. One-dimensional and two-dimensional shock waves tests are calculated to validate this new scheme. A three-dimensional turbulent magneto-convection simulation is used to show the applicability of current scheme to complicated astrophysical flows.

This paper reports an application of gas-kinetic BGK scheme to the computation of turbulent compressible convection in the stellar interior. After incorporating the Sub-grid Scale (SGS) turbulence model into the BGK scheme, we tested the effects of numerical parameters on the quantitative relationships among the thermodynamic variables, their fluctuations and correlations in a very deep, initially gravity-stratified stellar atmosphere. Comparison indicates that the thermal properties and dynamic properties are dominated by different aspects of numerical models separately. An adjustable Deardorff constant in the SGS model $c_\mu=0.25$ and an amplitude of artificial viscosity in the gas-kinetic BGK scheme $C_2=0$ are appropriate for current study. We also calculated the density-weighted auto- and cross-correlation functions in Xiong's (\cite{xiong77}) turbulent stellar convection theories based on which the gradient type of models of the non-local transport and the anisotropy of the turbulence are preliminarily studied. No universal relations or constant parameters were found for these models.