We present new theoretical predictions for the galaxy three-point correlation function (3PCF) using high-resolution dissipationless cosmologicalsimulations of a flat {Lambda}CDM Universe which resolve galaxy-size halos and subhalos. We create realistic mock galaxy catalogs by assigning luminosities and colors to dark matter halos and subhalos, and we measure the reduced 3PCF as a function of luminosity and color in both real and redshift space. As galaxy luminosity and color are varied, we find small differences in the amplitude and shape dependence of the reduced 3PCF, at a level qualitatively consistent with recent measurements from the SDSS and 2dFGRS. We confirm that discrepancies between previous 3PCF measurements can be explained in part by differences in binning choices. We explore the degree to which a simple local bias model can fit the simulated 3PCF. The agreement between the model ...
We compare the existing observational data on type Ia supernovae with the evolutions of the Universe predicted by a one-parameter family of tachyon models which we have introduced recently [Phys. Rev. D 69, 123512 (2004)]. Among the set of the trajectories of the model which are compatible with the data there is a consistent subset for which the Universe ends up in a new type of soft cosmological singularity dubbed big brake. This opens up yet another scenario for the future history of the Universe besides the one predicted by the standard #LAMBDA#CDM model.
This paper reports on Monte Carlo path integral simulations of a Bianchi Type IX cosmology which yield a wave function which reflects the known chaotic dynamics of the classical system and predicts a significant probability for a correlation between large universe volume and large anisotropy.
Anti-proton and positron Galactic cosmic ray (GCR) spectra are among the key targets for indirect detection of dark matter (DM). The boost factors, corresponding to an enhancement of the signal|linked to the clumpiness properties of the dark matter distribution|, have been taken as high as thousands in the past. The dramatic impact of these boost factors for indirect detection of antiparticles, for instance with the PAMELA satellite or the coming AMS-02 experiment, asks for their detailed calculation. We take into account the state-of-the-art results of high resolution N-body dark matter simulations to calculate the most likely energy dependent boost factors|linked to the GCR propagation properties|, for anti-protons and positrons. The results from extreme, but still possible, configurations of the clumpy dark matter component is also discussed. Starting from the mass and space distributions of sub-halos, the anti-proton and positron propagators are used to ...
(B204) The meeting will bring together observers and theorists in a highly interactive format, to further connect the local and cosmological star formation communities. Forward looking talks, aimed at the other communities, will survey terminology, achievements, problems and aspirations. Discussion will focus on the definition of the key questions, how the different communities can help each other, and preparations for the incorporation of realistic star formation into cosmologicalsimulations.
A Monte Carlo simulation of the vacuum Bianchi type-IX (mixmaster) cosmology yields a significant correlation between large universe volume and high anisotropy. An analog of the model's chaotic classical behavior is seen in the break up of the universe wave function at large volume into fingers in the corners of the minisuperspace anisotropy potential.
An estimate for the number of ionizing photons per baryon as a function of redshift is computed based on the plausible extrapolation of the observed galaxy UV luminosity function and the latest results on the properties of the escape fraction of ionizing radiation. It is found that, if the escape fraction for low mass galaxies (Mtot<10^{11}Msun) is assumed to be negligibly small, as indicated by numerical simulations, then there are not enough ionizing photons to reionize the universe by z=6 for the cosmology favored by the WMAP 3rd year results, while the WMAP 1st year cosmology is marginally consistent with the reionization requirement. The escape fraction as a function of galaxy mass would have to be constant to within a factor of two for the whole mass range of galaxies for reionization to be possible within the WMAP 3rd year cosmology.
A Monte Carlo simulation of the vacuum Bianchi type-IX (mixmaster) cosmology yields a significant correlation between large universe volume and high anisotropy. An analog of the model's chaotic classical behavior is seen in the break up of the universe wave function at large volume into fingers in the corners of the minisuperspace anisotropy potential.
Cosmological N-body and hydrodynamic simulations start with a realization of a random density fluctuation field representing a cosmological model at an early epoch. The density field is often replaced by a set of particles whose positions and velocities are set to conform to the desired density field. Each particle represents a cloud of huge number of real particles. Positions and velocities of particles are subsequently integrated by various numerical codes. We have simulated a set of collisionless collapses of Gaussian density peaks by using the PM and P(3) M codes. We find that in cosmologicalsimulations the physics at scales below the mean particle separation(MPS) is dominated by inaccuracies in describing the initial density field, and cannot be studied even by the high force-resolution codes. Since density fluctuations are ill-defined at scales smaller ...
We provide a quantitative assessment of the probability distribution function of the concentration parameter of galaxy clusters. We do so by using the probability distribution function of halo formation times, calculated by means of the excursion set formalism, and a formation redshift-concentration scaling derived from results of N-body simulations. Our results suggest that the observed high concentrations of several clusters are quite unlikely in the standard Lambda CDM cosmological model, but that due to various inherent uncertainties, the statistical range of the predicted distribution may be significantly wider than commonly acknowledged. In addition, the probability distribution function of the Einstein radius of A1689 is evaluated, confirming that the observed value of ~45" +/- 5" is very improbable in the currently favoured cosmological model. If, however, a variance of ~20% in the theoretically predicted value of ...
The Fourier transform of cosmological density perturbations can be represented in terms of amplitudes and phases for each Fourier mode. We investigate the phase evolution of these modes using a mixture of analytical and numerical techniques. Using a toy model of one-dimensional perturbations evolving under the Zel'dovich approximation as an initial motivation, we develop a statistic that quantifies the information content of the distribution of phases. Using numerical simulations beginning with more realistic Gaussian random-phase initial conditions, we show that the information content of the phases grows from zero in the initial conditions, first slowly and then rapidly when structures become non-linear. This growth of phase information can be expressed in terms of an effective entropy: Gaussian initial conditions are a maximum entropy realisation of the initial power spectrum, gravitational evolution decreases the phase entropy. We show that ...
Suppose an astronomer is equipped with a device capable of detecting emissions -- whether they be electromagnetic, gravitational, or neutrino -- from transient sources distributed throughout the cosmos. Because of source rate density evolution and variation of cosmological volume elements, the sources first detected when the machine is switched on are likely to be ones in the high-redshift universe; as observation time increases, rarer, more local, events will be found. We characterize the observer's evolving record of events in terms of a "probability event horizon", converging on the observer from great distances at enormous speed, and illustrate it by simulating neutron star birth events distributed throughout the cosmos. As an initial application of the concept, we determine the approach of this horizon for gamma-ray bursts (GRBs) by fitting to redshift data. The event rates required to fit the model are consistent with the proposed link ...
A new N-body and hydrodynamical code, called RAMSES, is presented. It has been designed to study structure formation in the universe with high spatial resolution. The code is based on Adaptive Mesh Refinement (AMR) technique, with a tree based data structure allowing recursive grid refinements on a cell-by-cell basis. The N-body solver is very similar to the one developed for the ART code (Kravtsov et al. 97), with minor differences in the exact implementation. The hydrodynamical solver is based on a second-order Godunov method, a modern shock-capturing scheme known to compute accurately the thermal history of the fluid component. The accuracy of the code is carefully estimated using various test cases, from pure gas dynamical tests to cosmological ones. The specific refinement strategy used in cosmologicalsimulations is described, and potential spurious effects associated to shock waves propagation in the resulting AMR ...
We present a showcase for deriving bounds on the neutrino masses from laboratory experiments and cosmological observations. We compare the frequentist and Bayesian bounds on the effective electron neutrino mass m_beta which the KATRIN neutrino mass experiment is expected to obtain, using both an analytical likelihood function and Monte Carlo simulations of KATRIN. Assuming a uniform prior in m_beta, we find that a null result yields an upper bound of about 0.17 eV at 90% confidence in the Bayesian analysis, to be compared with the frequentist KATRIN reference value of 0.20 eV. This is a significant difference when judged relative to the systematic and statistical uncertainties of the experiment. On the other hand, an input m_beta=0.35 eV, which is the KATRIN 5sigma detection threshold, would be detected at virtually the same level. Finally, we combine the simulated KATRIN results with cosmological data ...
Lyman-limit absorption systems can play many important roles during and after cosmological reionization. Unfortunately, due to the prohibitively large dynamic range required, it is impossible to self-consistently include these systems in cosmologicalsimulations. Using fast and versatile semi-numeric simulations, we systematically explore the spatial distribution of absorption systems during and following reionization. We self-calibrate the resulting number of absorbers to the mean free path (mfp) of the ionizing ultraviolet background (UVB), and present results at a given mfp and neutral hydrogen fraction. We use a simple optical depth criterion to identify the locations of absorbers. Our approach is fairly robust to uncertainties such as missing subgrid structure. Unlike at lower redshifts where the UVB is relatively uniform, at higher redshifts the fluctuations in the UVB and the HII morphology of ...
Dark matter particles form halos that contribute the major part of the mass of galaxy clusters. The formation of these cosmological structures have been investigated both observationally and in numerical simulations, which have confirmed the existence of a universal mass profile. However, the dynamic behaviour of dark matter in halos is not as well understood. We have used observations of 16 equilibrated galaxy clusters to show that the random velocities of dark matter particles are larger on average along the radial direction than along the tangential, and that the magnitude of this velocity anisotropy is radially varying. Our measurement implies that the collective behaviour of dark matter particles is fundamentally different from that of normal particles and the radial variation of the anisotropy velocity agrees with the predictions of numerical simulation.
The central engine of Gamma Ray Bursts is hidden from direct probing with photons mainly due to the high densities involved. Inferences on their properties are thus made from their cosmological setting, energetics, low-energy counterparts and variability. If GRBs are powered by hypercritical accretion onto compact objects, on small spatial scales the flow will exhibit fluctuations, which could in principle be reflected in the power output of the central engine and ultimately in the high energy prompt emission. Here we address this issue by characterizing the variability in neutrino cooled accretion flows through local shearing box simulations with magnetic fields, and then convolving them on a global scale with large scale dynamical simulations of accretion disks. The resulting signature is characteristic, and sensitive to the details of the cooling mechanism, providing in principle a discriminant for GRB central engine ...
We create realistic, full-sky, half-arcminute resolution simulations of the microwave sky matched to the most recent astrophysical observations. The primary purpose of these simulations is to test the data reduction pipeline for the Atacama Cosmology Telescope (ACT) experiment; however, we have widened the frequency coverage beyond the ACT bands and utilized the easily accessible HEALPix map format to make these simulations applicable to other current and near future microwave background experiments. Some of the novel features of these simulations are that the radio and infrared galaxy populations are correlated with the galaxy cluster and group populations, the primordial microwave background is lensed by the dark matter structure in the simulation via a ray-tracing code, the contribution to the thermal and kinetic Sunyaev-Zel'dovich (SZ) signals from ...
Cosmological hydrodynamical simulations as well as observations indicate that spiral galaxies are comprised of five different components: dark matter halo, stellar disc, stellar bulge, gaseous disc and gaseous halo. While the first four components have been extensively considered in numerical simulations of binary galaxy mergers, the effect of a hot gaseous halo has usually been neglected even though it can contain up to 80% of the total gas within the galaxy virial radius. We present a series of hydrodynamic simulations of major mergers of disc galaxies, that for the first time include a diffuse, rotating, hot gaseous halo. Through cooling and accretion, the hot halo can dissipate and refuel the cold gas disc before and after a merger. This cold gas can subsequently form stars, thus impacting the morphology and kinematics of the remnant. Simulations of isolated systems with total ...
We investigate the formation and evolution of dwarf galaxies in a high resolution, hydrodynamical cosmologicalsimulation of a Milky Way sized halo and its environment. Our simulation includes gas cooling, star formation, supernova feedback, metal enrichment and UV heating. In total, 90 satellites and more than 400 isolated dwarf galaxies are formed in the simulation, allowing a systematic study of the internal and environmental processes that determine their evolution. We find that 95% of satellite galaxies are gas-free at z=0, and identify three mechanisms for gas loss: supernova feedback, tidal stripping, and photo-evaporation due to re-ionization. Gas-rich satellite galaxies are only found with total masses above ~ 5x10^9 solar masses. In contrast, for isolated dwarf galaxies, a total mass of ~ 10^9 solar masses constitutes a sharp transition; less massive galaxies are predominantly gas-free at z=0, ...
We examine the late-time behavior of initially expanding homogeneous cosmological models satisfying Einstein's equation with a positive cosmological constant ..lambda... It is shown that such models of all Bianchi types except IX exponentially evolve toward the de Sitter solution, with time scale (3/..lambda..)/sup 1/2/. The behavior of Bianchi type-IX universes is similar, provided that ..lambda.. is sufficiently large compared with spatial-curvature terms. Thus, a positive cosmological constant provides an effective means of isotropizing homogeneous universes.
We examine the late-time behavior of initially expanding homogeneous cosmological models satisfying Einstein's equation with a positive cosmological constant #LAMBDA#. It is shown that such models of all Bianchi types except IX exponentially evolve toward the de Sitter solution, with time scale (3/#LAMBDA#)/sup 1/2/. The behavior of Bianchi type-IX universes is similar, provided that #LAMBDA# is sufficiently large compared with spatial-curvature terms. Thus, a positive cosmological constant provides an effective means of isotropizing homogeneous universes.
Abstract We present the first numerical,-N-body, hydrodynamical, chemical simulations of cosmic structure formation in the framework of non-Gaussian models. We study the impact of primordial non-Gaussianities on early chemistry (e-, H, H+, H-, He, He+, He++, H2, H+2, D, D+, HD, HeH+), molecular and atomic gas cooling, star formation, metal (C, O, Si, Fe, Mg, S) enrichment, Population-III (popIII) and Population-II-I (popII) transition and on the evolution of -visible- objects. We find that non-Gaussianities can have some consequences on baryonic structure formation at very early epochs, but the subsequent evolution at later times washes out any difference among the various models. When assuming reasonable values for primordial non-Gaussian perturbations, it turns out that they are responsi...
One hope to solve the cosmological constant problem is to identify a symmetry principle, based on which the cosmological constant can be reduced either to zero, or to a tiny value. Here, we note that requiring that the vacuum state is Lorentz invariant significantly reduces the theoretical value of the vacuum energy density. Hence, this also reduces the discrepancy between the observed value of the cosmological constant and its theoretical expectation, down from 123 orders of magnitude to 56 orders of magnitude. We find that, at one loop level, massless particles do not yield any contribution to the cosmological constant. Another important consequence of Lorentz symmetry is stabilization of the gravitational hierarchy: the cosmological constant (divided by Newton's constant) does not run as the quartic power of the renormalization group scale, but instead only logarithmically.
We propose an X-ray mission called Xenia to search for decaying superweakly interacting Dark Matter particles (super-WIMP) with a mass in the keV range. The mission and its observation plan are capable of providing a major break through in our understanding of the nature of Dark Matter (DM). It will confirm, or reject, predictions of a number of particle physics models by increasing the sensitivity of the search for decaying DM by about two orders of magnitude through a wide-field imaging X-ray spectrometer in combination with a dedicated observation program. The proposed mission will provide unique limits on the mixing angle and mass of neutral leptons, right handed partners of neutrinos, which are important Dark Matter candidates. The existence of these particles is strongly motivated by observed neutrino flavor oscillations and the problem of baryon asymmetry of the Universe. In super-WIMP models, the details of the formation of the cosmic web are different from those of ...
We study the role of rolling tachyons in the cosmological model with dilatonic gravity. In the string frame, flat space solutions of both initial-stage and late-time are obtained in closed form. In the Einstein frame, we show that every expanding solution is decelerating.
We study the impact of non-Abelian T-duality transformations on a string based cosmological model. The implementation of the pre-big-bang scenario is investigated. We found a region of the dual phase where such a picture is possible.
We show the existence of a rather general class of closed cosmological models of Bianchi type IX that do not exhibit recollapse but expand for all times. This is despite the fact that these models satisfy the strong energy condition by a wide margin.
Solutions are obtained for the Einstein field equations governing 40 rotating and time-dependent Bianchi type-IX cosmologies with thermal and nonthermal perfectly fluid sources. Some general properties of the solutions are discussed, and the exact computational details are given in a companion paper. All the cosmologies have spacelike, timelike or null-like homogeneous hypersurfaces depending on whether the constant alpha is less than or greater than one. The pressure of the cosmologies is assumed to be zero, and the energy density is less than zero. 9 references.
We discuss the asymptotic dynamical evolution of spatially homogeneous brane-world cosmological models close to the initial singularity. We find that generically the cosmological singularity is isotropic in Bianchi type IX brane-world models and consequently these models do not exhibit Mixmaster or chaotic-like behaviour close to the initial singularity. We argue that this is typical of more general cosmological models in the brane-world scenario. In particular, we show that an isotropic singularity is a past-attractor in all orthogonal Bianchi models and is a local past-attractor in a class of inhomogeneous brane-world models. (letter to the editor)
We discuss various superstring effective actions and, in particular, their common sector which leads to the so-called pre-big-bang cosmology (cosmology in a weak coupling limit of heterotic superstring theory). Using the conformal relationship between these two theories we present Kasner asymptotic solutions of Bianchi type IX geometries within these theories and make predictions about possible emergence of chaos. Finally, we present a possible method of generating Horava-Witten cosmological solutions out of the well-known general relativistic or pre-big-bang solutions.
We discuss various superstring effective actions and, in particular, their common sector which leads to the so-called pre-big-bang cosmology (cosmology in a weak coupling limit of heterotic superstring theory. Using the conformal relationship between these two theories we present Kasner asymptotic solutions of Bianchi type IX geometries within these theories and make predictions about possible emergence of chaos. Finally, we present a possible method of generating Horava-Witten cosmological solutions out of the well-known general relativistic or pre-big-bang solutions.
A Bianchi type-IX cosmological model has been found as a solution of Einstien's vacuum field equations with a cosmological constant. The solution represents a rotating generalziation fo the de Sitter universe. This universe shows a transition to exponential expansion and the vorticity begins to decay exponentially at the grand-unified-theory time. The point of time for this transition is independent of the magnitude of the vorticity. During the Guth inflationary era the vorticity decays by a factor of the order 10 US.
A Bianchi type-IX cosmological model has been found as a solution of Einstien's vacuum field equations with a cosmological constant. The solution represents a rotating generalziation fo the de Sitter universe. This universe shows a transition to exponential expansion and the vorticity begins to decay exponentially at the grand-unified-theory time. The point of time for this transition is independent of the magnitude of the vorticity. During the Guth inflationary era the vorticity decays by a factor of the order 10"-"1"4"2.
We describe exact cosmological solutions with rotation and expansion in the low-energy effective string theory. These models are spatially homogeneous (closed Bianchi type IX) and they belong to the family of shear-free metrics which are causal (no closed timelike curves are allowed), admit no parallax effects and do not disturb the isotropy of the background radiation. The dilaton and the axion fields are nontrivial, in general, and we consider both cases with and without the central charge (effective cosmological constant)
We describe exact cosmological solutions with rotation and expansion in the low-energy effective string theory. These models are spatially homogeneous (closed Bianchi type IX) and they belong to the family of shear-free metrics which are causal (no closed timelike curves are allowed), admit no parallax effects and do not disturb the isotropy of the background radiation. The dilaton and the axion fields are nontrivial, in general, and we consider both cases with and without the central charge (effective cosmological constant).
The propriety of the cosmic no-hair conjecture to the Bianchi-type-IX spacetime is discussed from a quantum cosmological point of view. It is shown that most, but not all, classical universes which are created quantum cosmologically are inflationary. The probability of inflation among such universes is also discussed.
The propriety of the cosmic no-hair conjecture to the Bianchi-type-IX spacetime is discussed from a quantum cosmological point of view. It is shown that most, but not all, classical universes which are created quantum cosmologically are inflationary. The probability of inflation among such universes is also discussed.
A higher-dimensional homogeneous spacetime is investigated satisfying the vacuum Einstein equations. It is assumed that the algebra of Killing vectors L admits a non-trivial Levi decomposition L=N+so(3), i.e. that the subalgebras N and so(3) do not commute. It is found that the model behaves in a non-chaotic way and cosmological dimensional reduction inevitably occurs. This model completes all the possible types within the class of higher-dimensional extensions of Bianchi type-IX cosmology.
The problem on change of deuterium abundance in the process of galactic evolution (star evolution, supernova explosions, nucleosynthesis in supermassive objects) is considered. It is shown that the observable deuterium quantity in the interstellar medium must correspond to its cosmological abundance. This conclusion is independent of the rate of accretion of intergalactic gas by Galaxy. The effect of hypothetical pregalactic active objects on cosmological deuterium is small. It is poind out that observations of interstellar deuterium in absorbtion at lambda=91.6 cm are significant.
Nordstroem's first theory of gravitation, which is a Lorentz covariant scalar theory, is discussed, and the cosmological solution to the field equation is derived. It is shown that there are two physically equivalent representations of the theory that differ in the system of units with which measurements are made. The cosmological solution is found to be an infinite, Euclidean universe with a big bang and a big crunch.
A higher-dimensional homogeneous spacetime is investigated satisfying the vacuum Einstein equations. It is assumed that the algebra of Killing vectors L admits a non-trivial Levi decomposition L=N+so(3), i.e. that the subalgebras N and so(3) do not commute. It is found that the model behaves in a non-chaotic way and cosmological dimensional reduction inevitably occurs. This model completes all the possible types within the class of higher-dimensional extensions of Bianchi type-IX cosmology. (orig.).
Recent measurements of cosmic ray leptons by PAMELA, ATIC, HESS and Fermi revealed interesting excesses. Many authors suggested particle Dark Matter (DM) annihilations could be at the origin of these effects. In this paper, we critically assess this interpretation by reviewing some results questioning the naturalness and robustness of such an interpretation. Natural values for the DM particle parameters lead to a poor leptons production so that models often require signal enhancement effects that we constrain here. Considering DM annihilations are likely to produce antiprotons as well, we use the PAMELA antiproton to proton ratio measurements to constrain a possible exotic contribution. We also consider the possibility of an enhancement due to a nearby clump of DM. This scenario appears unlikely when compared to the state-of-the-art cosmological N-body simulations. We conclude that the bulk of the observed signals most likely has no link with ...
We consider the simplest scenario when black strings (cigars) penetrate the cosmological brane. As a result, the brane has a Swiss-cheese structure, with Schwarzschild black holes immersed in a Friedmann-Lema\\^{\\i}tre-Robertson-Walker brane. There is no dark radiation in the model, the cosmological regions of the brane are characterized by a cosmological constant $\\Lambda$ and flat spatial sections. Regardless of the value of $\\Lambda$, these brane-world universes forever expand and forever decelerate. The totality of source terms in the modified Einstein equation sum up to a dust, establishing a formal equivalence with the general relativistic Einstein-Straus model. However in this brane-world scenario with black strings the evolution of the cosmological fluid strongly depends on $\\Lambda$. For $\\Lambda$ less or equal to zero it has positive energy density $\\rho$ and negative pressure $p$ and at ...
Recent work on [ital N]=2 supersymmetric Bianchi type IX cosmologies coupled to a scalar field is extended to a general treatment of homogeneous quantum cosmologies with explicitly solvable momentum constraints, i.e., Bianchi types I, II, VII, VIII in addition to the Bianchi type IX, and special cases, namely, the Freidmann universes, the Kantowski-Sachs space, and Taub-NUT space. In addition to the earlier explicit solution of the Wheeler-DeWitt equation for Bianchi type IX, describing a virtual wormhole fluctuation, an additional explicit solution is given and identified with the no-boundary state.''
Recent work on N=2 supersymmetric Bianchi type IX cosmologies coupled to a scalar field is extended to a general treatment of homogeneous quantum cosmologies with explicitly solvable momentum constraints, i.e., Bianchi types I, II, VII, VIII in addition to the Bianchi type IX, and special cases, namely, the Freidmann universes, the Kantowski-Sachs space, and Taub-NUT space. In addition to the earlier explicit solution of the Wheeler-DeWitt equation for Bianchi type IX, describing a virtual wormhole fluctuation, an additional explicit solution is given and identified with the ''no-boundary state.''.
We perform a detailed phase-space analysis of various phantom cosmological models, where the dark energy sector interacts with the dark matter one. We examine whether there exist late-time scaling attractors, corresponding to an accelerating universe and possessing dark energy and dark matter densities of the same order. We find that all the examined models, although accepting stable late-time accelerated solutions, cannot alleviate the coincidence problem, unless one imposes a form of fine-tuning in the model parameters. It seems that interacting phantom cosmology cannot fulfill the basic requirement that led to its construction.
In the inflationary scenario of loop quantum cosmology (LQC) in the presence of inverse-volume corrections, we give analytic formulas for the power spectra of scalar and tensor perturbations convenient to confront with observations. Since inverse-volume corrections can provide strong contributions to the running spectral indices, inclusion of terms higher than the second-order runnings in the power spectra is crucially important. Using the recent data of cosmic microwave background (CMB) and other cosmological experiments, we place bounds on the quantum corrections for a quadratic inflaton potential.
This is a design study into the capabilities of the Australian Square Kilometre Array Pathfinder (ASKAP) in performing a full-sky low redshift neutral hydrogen survey, termed WALLABY, and the potential cosmological constraints one can attain from measurement of the matter power spectrum. We find that the full sky survey will likely attain 600,000 redshifts which, when combined with expected Planck CMB data, will constrain the Dark Energy equation of state to 20%, for the first time making cosmological constraints from radio observations competitive with the best existing optical surveys.
A previously studied theory of gravitation in flat space-time is applied to homogeneous and isotropic cosmological models. There exist two different classes of models without singularities: (i) ever-expanding models, (ii) oscillating models. The first class contains models with hot big bang. For these models there exist at the beginning of the universe-in contrast to Einstein's theory-very high but finite densities of matter and radiation with a big bang of very short duration. After short time these models pass into the homogeneous and isotropic models of Einstein's theory with spatial curvature equal to zero and cosmological constant ALPHA >= O.
If quantum fields exist in extra compact dimensions, they will give rise to a quantum vacuum or Casimir energy. That vacuum energy will manifest itself as a cosmological constant. The fact that supernova and cosmic microwave background data indicate that the cosmological constant is of the same order as the critical mass density to close the universe supplies a lower bound on the size of the extra dimensions. Recent laboratory constraints on deviations from Newton's law place an upper limit. The allowed region is so small as to suggest that either extra compact dimensions do not exist, or their number is about to be tightly constrained by experimental data.
If quantum fields exist in extra compact dimensions, they will give rise to a quantum vacuum or Casimir energy. That vacuum energy will manifest itself as a cosmological constant. The fact that supernova and cosmic microwave background data indicate that the cosmological constant is of the same order as the critical mass density to close the universe supplies a lower bound on the size of the extra dimensions. Recent laboratory constraints on deviations from Newton's law place an upper limit. The allowed region is so small as to suggest that either extra compact dimensions do not exist, or their properties are about to be tightly constrained by experimental data.
The discovery that the cosmic expansion is accelerating has been followed by an intense theoretical and experimental response in physics and astronomy. The discovery implies that our most basic notions about how gravity work are violated on cosmological distance scales. One simple fix is the introduction of a cosmological constant into the field equations for general relativity. However, the extremely small value of the cosmological constant, relative to theoretical expectations, has led theorists to explore a wide variety of alternative explanations that involve the introduction of an exotic negative-pressure fluid or a modification of general relativity. Here we briefly review the evidence for cosmic acceleration. We then survey some of the theoretical attempts to account for it, including the cosmological constant, quintessence and its variants, mass-varying neutrinos, and modifications of general ...
Versions of parameterized pseudo-Newtonian gravity theories specially designed for cosmology have been introduced in recent cosmology literature. The modifications demand a zero-pressure fluid in the context of versions of modified Poisson-like equation with two different gravitational potentials. We consider such modifications in the context of relativistic gravity theories where the action is a general algebraic function of the scalar curvature, the scalar field, and the kinetic term of the field. In general it is not possible to isolate the zero-pressure fluid component simultaneously demanding a modification in the Poisson-like equation. Only in the small-scale limit we can realize some special forms of the attempted modifications. We address some loopholes in the possibility of showing non-Einstein gravity nature based on pseudo-Newtonian modifications in the cosmological context. We point out that future observations ...
The next decade promises an observational revolution which will change cosmology forever. The precise measurement of the angular anisotropy of the cosmic microwave background should specify to a few percent all of the parameters of the cosmological model which effect astrophysics. The growth of structure will then be determined (but not yet observed) until gravitational collapse becomes highly non-linear and stars, galaxies and active galactic nuclei (AGN) form. These processes are hard to model with basic physics because they are complex and allow a rich variety of expression. Instead observations will determine when the first stars and quasars formed, and how and when galaxies assembled. If we can reconcile the numerous contradictions which characterize the subject today, cosmology will become a mature subject, founded on the agreement between detailed, inclusive and realistic models, which make precise predictions, and ...
We investigate the quantum cosmology of spatially homogeneous models with compact spatial sections admitting a u(2) isometry algebra. The metric ansatz in these models is that of Bianchi type IX with two scale factors set to be equal. We apply the Hartle-Hawking no-boundary path integral prescription and find the semi-classical contributions to the wave function. Exact formulae are obtainable for certain contributions and otherwise the limits of large and small anisotropy (for the pure vacuum case) and large spatial volume or small anisotropy (for the case with a positive cosmological constant) are considered. For the pure vacuum case we find no semiclassical components which would correspond to Lorentzian universes. For the case with a cosmological constant the Hartle-Hawking boundary conditions formally constrain one of the parameters in the Lorentzian solutions to be purely imaginary. Possible interpretations of this ...
We investigate the quantum cosmology of spatially homogeneous models with compact spatial sections admitting a u(2) isometry algebra. The metric ansatz in these models is that of Bianchi type IX with two scale factors set to be equal. We apply the Hartle-Hawking no-boundary path integral prescription and find the semi-classical contributions to the wave function. Exact formulae are obtainable for certain contributions and otherwise the limits of large and small anisotropy (for the pure vacuum case) and large spatial volume or small anisotropy (for the case with a positive cosmological constant) are considered. For the pure vacuum case we find no semiclassical components which would correspond to Lorentzian universes. For the case with a cosmological constant the Hartle-Hawking boundary conditions formally constrain one of the parameters in the Lorentzian solutions to be purely imaginary. Possible interpretations of this ...
Warm-intermediate inflationary universe models in the context of braneworld cosmologies, are studied. This study is done in the weak and strong dissipative regimes. We find that, the scalar potentials and dissipation coefficients in terms of the scalar field, evolves as type-power-law and powers of logarithms, respectively. General conditions required for these models to be realizable are derived and discussed. We also study the scalar and tensor perturbations for each regime. We use recent astronomical observations to constraint the parameters appearing in the braneworld models.
The Thomsen and Frandsen (1983) cosmological test employing the observed correlation between elliptical galaxy surface brightness and scale size in order to avoid problems due to dynamical evolution is presently extended to make use of arbitrary galaxy samples. The explicit dependence on stellar evolution is also demonstrated. On the assumption that this evolution is calculable, an equation for the deceleration parameter entirely in terms of observables is derived. The test is applied to two available samples. 25 references.
The loop quantum cosmology 'improved dynamics' of the Bianchi type IX model are studied. The action of the Hamiltonian constraint operator is obtained via techniques developed for the Bianchi type I and type II models, no new input is required. It is shown that the big bang and big crunch singularities are resolved by quantum gravity effects. We also present effective equations which provide quantum geometry corrections to the classical equations of motion.
The possibility of having inflation in a renormalizable cosmological model is investigated. The cosmic no-hair conjecture is proved to hold for all Bianchi types except Bianchi type IX. By the use of a conformal transformation on the metric we show that these models are equivalent to the ones described by the Einstein-Hilbert action for gravity minimally coupled to a set of scalar fields with inflationary potentials. Henceforth, we prove that inflationary solutions behave as attractors in solution space, making it a natural event in the evolution of such models.
The possibility of having inflation in a renormalizable cosmological model is investigated. The cosmic no-hair conjecture is proved to hold for all Bianchi types except Bianchi type IX. By the use of a conformal transformation on the metric we show that these models are equivalent to the ones described by the Einstein-Hilbert action for gravity minimally coupled to a set of scalar fields with inflationary potentials. Henceforth, we prove that inflationary solutions behave as attractors in solution space, making it a natural event in the evolution of such models.
The t<0 branch of pre-big bang cosmological scenarios is subject to a gravitational wave instability. The unstable behaviour of tensor perturbations is derived in a very simple way in Hwang's covariant and gauge-invariant formalism developed for extended theories of gravity. A simple interpretation of this instability as the effect of an ''antifriction'' is given, and it is argued that a universe must eventually enter the expanding phase. (orig.)
A consistent combination of quantum geometry effects rules out a large class of models of loop quantum cosmology and their critical densities as they have been used in the recent literature. In particular, the critical density at which an isotropic universe filled with a free, massless scalar field would bounce must be well below the Planck density. In the presence of anisotropy, no model of the Schwarzschild black hole interior analyzed so far is consistent.
The dimensionally reduced effective action of the bosonic sector of the heterotic string in critical dimensions is employed to derive a Wheeler-DeWitt equation for the Bianchi type-IX cosmology. An exact solution is found that becomes strongly peaked around the isotropic limit as the volume of the three-geometry increases. In principle the global O(6,6) symmetry of the effective action can be employed to generate new solutions from the one presented here.
The dimensionally reduced effective action of the bosonic sector of the heterotic string in critical dimensions is employed to derive a Wheeler-DeWitt equation for the Bianchi type-IX cosmology. An exact solution is found that becomes strongly peaked around the isotropic limit as the volume of the three-geometry increases. In principle the global O(6,6) symmetry of the effective action can be employed to generate new solutions from the one presented here.
Ekpyrotic and cyclic cosmologies provide theories of the very early and of the very late universe. In these models, the big bang is described as a collision of branes - and thus the big bang is not the beginning of time. Before the big bang, there is an ekpyrotic phase with equation of state w=P/({rho}) >>1 (where P is the average pressure and {rho} the average energy density) during which the universe slowly contracts. This phase resolves the standard cosmological puzzles and generates a nearly scale-invariant spectrum of cosmological perturbations containing a significant non-Gaussian component. At the same time it produces small-amplitude gravitational waves with a blue spectrum. The dark energy dominating the present-day cosmological evolution is reinterpreted as a small attractive force between our brane and a parallel one. This force eventually induces a new ekpyrotic phase and a new ...
The red sequence is an important feature of galaxy clusters and plays a crucial role in optical cluster detection. Measurement of the slope and scatter of the red sequence are affected both by selection of red sequence galaxies and measurement errors. In this paper, we describe a new error corrected Gaussian Mixture Model for red sequence galaxy identification. Using this technique, we can remove the effects of measurement error and extract unbiased information about the intrinsic properties of the red sequence. We use this method to select red sequence galaxies in each of the 13,823 clusters in the maxBCG catalog, and measure the red sequence ridgeline location and scatter of each. These measurements provide precise constraints on the variation of the average red galaxy populations in the observed frame with redshift. We find that the scatter of the red sequence ridgeline increases mildly with redshift, and that the slope decreases with redshift. We also observe that the slope does ...
The observation that Type Ia supernovae are fainter than expected given their red shifts has led to the conclusion that the expansion of the universe is accelerating. The widely accepted hypothesis is that this acceleration is caused by a cosmological constant or, more generally, some dark energy field that pervades the universe. This hypothesis presents a challenge to physics so severe that one is motivated to explore alternative explanations. In this paper, we explore whether the data from Type Ia supernovae can be explained with an idea that is almost as old as that of the cosmological constant, namely, that the strength of gravity varies on a cosmic timescale. This topic is an ideal one for investigation by an undergraduate physics major because the entire chain of reasoning from models to data analysis is well within the mathematical and conceptual sophistication of a motivated undergraduate.
We study the dynamics of states perturbatively expanded about a harmonic system of loop quantum cosmology, exhibiting a bounce. In particular, the evolution equations for the first and second order moments of the system are analyzed. These moments back-react on the trajectories of the expectation values of the state and hence alter the energy density at the bounce. This analysis is performed for isotropic loop quantum cosmology coupled to a scalar field with a small but non-zero constant potential, hence in a regime in which the kinetic energy of matter dominates. Analytic restrictions on the existence of dynamical coherent states and the meaning of semi-classicality within these systems are discussed. A numerical investigation of the trajectories of states that remain semi-classical across the bounce demonstrates that, at least for such states, the bounce persists and that its properties are similar to the standard case, in which the moments ...
Cosmological constant a.k.a. dark energy problem is considered to be one major challenge in modern cosmology. Here we present a model where large scale structure formation causes spatially-flat FRW universe to fragment into numerous `FRW islands' surrounded by vacuum. We show that this mechanism can explain the origin of dark energy as well as the late time cosmic acceleration. This explanation of dark energy does not require any exotic matter source nor an extremely fine-tuned cosmological constant. This explanation is given within classical general relativity and relies on the fact that our universe has been undergoing structure formation since its recent past.
The full theory and the semiclassical description of loop quantum cosmology (LQC) have been studied in the Friedmann-Robertson-Walker and Bianchi I models. As an extension to include both anisotropy and intrinsic curvature, this paper investigates the cosmological model of Kantowski-Sachs spacetime with a free massless scalar field at the level of phenomenological dynamics with the LQC discreteness corrections. The LQC corrections are implemented in two different improved quantization schemes. In both schemes, the big bang and big crunch singularities of the classical solution are resolved and replaced by the big bounces when the area or volume scale factor approaches the critical values in the Planck regime measured by the reference of the scalar field momentum. Symmetries of scaling are also noted and suggest that the fundamental spatial scale (area gap) may give rise to a temporal scale. The bouncing scenarios are in an analogous fashion of ...
General relativity postulates the Minkowski space-time to be the standard flat geometry against which we compare all curved space-times and the gravitational ground state where particles, quantum fields and their vacuum states are primarily conceived. On the other hand, experimental evidences show that there exists a non-zero cosmological constant, which implies in a deSitter space-time, not compatible with the assumed Minkowski structure. Such inconsistency is shown to be a consequence of the lack of a application independent curvature standard in Riemann's geometry, leading eventually to the cosmological constant problem in general relativity. We show how the curvature standard in Riemann's geometry can be fixed by Nash's theorem on locally embedded Riemannian geometries, which imply in the existence of extra dimensions. The resulting gravitational theory is more general than general relativity, similar to brane-world gravity, but where the ...
The Lorentz and coordinate covariant calculus of spinors in Riemannian spacetime, which is the mathematical model for the description of the quantum mechanics of elementary particles with spin interacting with the classical gravitation field, is explored. The Dirac equation describing the interaction of neutrinos with the gravitational fields of the Robertson-Walker cosmological world models is separated, and the spectrum of eigenfunctions and eigenvalues for particular choices of the set of quantum numbers is given explicitly for the k = 0 and k = +1 models, although only the radial equations determining the final quantum number are given for the k = -1 model. The mathematical theory of the motion of a perfect fluid whose elements interact via long-range neutrino-exchange forces, as well as gravitationally, is developed. The formalism for calculating, by calculating the Bogoliubov transformation of the Fock space operators that instantaneously diagonalize their ...
After mathematicians and physicists had learned that the structure of physical space was not necessarily Euclidean, it became conceivable that the global topological structure of space was non-trivial. In the context of the late 19th century debates on physical space this speculation gave rise to the problem of classifying spaces of constant curvature from a topological point of view. William Kingdon Clifford, Felix Klein and Wilhelm Killing, the latter of whom devoted a substantial amount of work to the topic in the early 1890s, clearly perceived this problem as relevant for both mathematics and natural philosophy (i.e., physics or cosmology). To some extent, a cosmological interest may even be found among those authors who restated the space form problem in more modern terms in the early 20th century, such as Heinz Hopf.
Obtaining the set of cosmological parameters consistent with observational data is an important exercise in current cosmological research. It involves finding the global maximum of the likelihood function in the multi-dimensional parameter space. Currently sampling based methods, which are in general stochastic in nature, like Markov-Chain Monte Carlo(MCMC), are being commonly used for parameter estimation. The beauty of stochastic methods is that the computational cost grows, at the most, linearly in place of exponentially (as in grid based approaches) with the dimensionality of the search space. MCMC methods sample the full joint probability distribution (posterior) from which one and two dimensional probability distributions, best fit (average) values of parameters and then error bars can be computed. In the present work we demonstrate the application of another stochastic method, named Particle Swarm Optimization (PSO), that is widely used ...
Gravitational lensing distorts the cosmic microwave background (CMB) temperature and polarization fields and encodes valuable information on distances and growth rates at intermediate redshifts into the lensed power spectra. The non-Gaussian bandpower covariance induced by the lenses is negligible to l=2000 for all but the B polarization field where it increases the net variance by up to a factor of 10 and favors an observing strategy with 3 times more area than if it were Gaussian. To quantify the cosmological information, we introduce two lensing observables, characterizing nearly all of the information, which simplify the study of non-Gaussian impact, parameter degeneracies, dark energy models, and complementarity with other cosmological probes. Information on the intermediate redshift parameters rapidly becomes limited by constraints on the cold dark matter density and initial amplitude of fluctuations as observations improve. Extraction of ...
An exact cosmological solution of Einstein's equations which has time-dependent rotation is presented. The t-constant sections are of Bianchi type II. The source of this geometry is a fluid which has not been thermalized. (Author).
The current status of cosmological observations is presented; and the light curves and radiation spectra from supernova SN1987A are used in comparison between expected and observed universal nucleosynthesis and star evolution data./aip/.
The canonical quantum theory of gravity-quantum geometrodynamics (QG)-is applied to the homogeneous Bianchi type IX cosmological model. As a result, a framework for the quantum theory of homogeneous cosmologies is developed. We show that the theory is internally consistent and prove that it possesses the correct classical limit (the theory of general relativity). To emphasize the special role that the constraints play in this new theory, we compare it to the traditional ADM square-root and Wheeler-DeWitt quantization schemes. We show that, unlike traditional approaches, QG leads to a well-defined Schroedinger equation for the wavefunction of the universe that is inherently coupled to the expectation value of the constraint equations. This coupling to the constraints is responsible for the appearance of a coherent spacetime picture. Thus, the physical meaning of the constraints of the theory is quite different from Dirac's ...
The canonical quantum theory of gravity-quantum geometrodynamics (QG)-is applied to the homogeneous Bianchi type IX cosmological model. As a result, a framework for the quantum theory of homogeneous cosmologies is developed. We show that the theory is internally consistent and prove that it possesses the correct classical limit (the theory of general relativity). To emphasize the special role that the constraints play in this new theory, we compare it to the traditional ADM square-root and Wheeler-DeWitt quantization schemes. We show that, unlike traditional approaches, QG leads to a well-defined Schroedinger equation for the wavefunction of the universe that is inherently coupled to the expectation value of the constraint equations. This coupling to the constraints is responsible for the appearance of a coherent spacetime picture. Thus, the physical meaning of the constraints of the theory is quite different from Dirac's interpretation. In ...
It is proven that there do not exist any diagonal Bianchi type-IX universes which expand for an infinite time, provided only that the matter satisfies the dominant energy condition and has non-negative average principal pressures.
It is proven that there do not exist any diagonal Bianchi type-IX universes which expand for an infinite time, provided only that the matter satisfies the dominant energy condition and has non-negative average principal pressures.
indian academy of sciences - journal of astrophysics and astronomy on astronomy (eds. v. k. kapahi, n. k. dadhich, g. swarup and special issue in honour of s. chandrasekhar big bang and alternative cosmologies; ...
The energy-momentum tensor of a massless spinor field is constructed and studied based on the previously proposed interpretation of quantum effects of such a field in the anisotropic metric of Bianchi type IX. The characteristic properties of the energy-momentum tensor in the mixed universe model are discussed.
We investigate the quantum cosmology of spatially homogeneous models with compact spatial sections admitting a u(2) isometry algebra. The metric ansatz in these models is that of Bianchi type IX with two scale factors set to be equal. We apply the Hartle-Hawking no-boundary path integral prescription and find the semi-classical contributions to the wave function. Exact formulae are obtainable for certain contributions and otherwise the limits of large and small anisotropy (for the pure vacuum case) and large spatial volume or small anisotropy (for the case with a positive cosmological constant) are considered. For the pure vacuum case we find no rapidly oscillating semiclassical components in the wave function, and hence do not recover lorentzian space-time as a prediction of the no-boundary proposal. For the case with a cosmological constant the wave function does contain rapidly oscillating components and thus predicts ...
We investigate the quantum cosmology of spatially homogeneous models with compact spatial sections admitting a u(2) isometry algebra. The metric ansatz in these models is that of Bianchi type IX with two scale factors set to be equal. We apply the Hartle-Hawking no-boundary path integral prescription and find the semi-classical contributions to the wave function. Exact formulae are obtainable for certain contributions and otherwise the limits of large and small anisotropy (for the pure vacuum case) and large spatial volume or small anisotropy (for the case with a positive cosmological constant) are considered. For the pure vacuum case we find no rapidly oscillating semiclassical components in the wave function, and hence do not recover lorentzian space-time as a prediction of the no-boundary proposal. For the case with a cosmological constant the wave function does contain rapidly oscillating components and thus predicts ...
Observational work on quasars, galaxies, and stars is summarized. Theoretical studies covering stars and stellar evolution, galaxies, clusters and cosmology, high energy astrophysics the solar system and the Sun are described. (ESA)
We investigate several varying-mass dark-matter particle models in the framework of phantom cosmology. We examine whether there exist late-time cosmological solutions, corresponding to an accelerating universe and possessing dark energy and dark matter densities of the same order. Imposing exponential or power-law potentials and exponential or power-law mass dependence, we conclude that the coincidence problem cannot be solved or even alleviated. Thus, if dark energy is attributed to the phantom paradigm, varying-mass dark matter models cannot fulfill the basic requirement that led to their construction.
By employing the higher dimensional version of the Wu-Yang Ansatz we obtain magnetically charged new black hole solutions in the Einstein-Yang-Mills-Lovelock (EYML) theory with second ($\\alpha2$) and third ($\\alpha3$)order parameters. These parameters, where $\\alpha2$ is also known as the Gauss-Bonnet parameter, modify the horizons (and the resulting thermodynamical properties) of the black holes. It is shown also that asymptotically, these parameters contribute to an effective cosmological constant -without cosmological constant- so that the solution behaves de-Sitter (Anti de-Sitter) like.
The geodesic equation for space-like objects moving along a circular trajectory in the expanding universe is considered. Our analysis leads to the conclusion that ancient cosmological tachyons may exist in the present-day world and may play an important role in (i) the internal structure of hadrons conceived as nonlocal objects called strings, (ii) the T-symmetry violation observed in the weak K-decays, (iii) the multidimensional unified field theories of Kaluza-Klein type, and in (iv) the classical models of charged particles which combine ordinary electromagnetism with a self-interacting version of Newtonian gravity. 18 refs.
The geodesic equation for space-like objects moving along a circular trajectory in the expanding universe is considered. Our analysis leads to the conclusion that ancient cosmological tachyons may exist in the present-day world and may play an important role in (i) the internal structure of hadrons conceived as nonlocal objects called strings, (ii) the T-symmetry violation observed in the weak K-decays, (iii) the multidimensional unified field theories of Kaluza-Klein type, and in (iv) the classical models of charged particles which combine ordinary electromagnetism with a self-interacting version of Newtonian gravity. 18 refs.
A (n +1)-dimensional Einstein-Gauss-Bonnet (EGB) model is considered. For diagonal cosmological type metrics the equations of motion are reduced to a set of Lagrange equations. The effective Lagrangian contains two minisuperspace metrics on R^n. The first one is the well-known 2-metric of pseudo-Euclidean signature and the second one is the Finslerian 4-metric that is proportional to n-dimensional Berwald-Moor 4-metric. When a synchronous-like time gauge is considered the equations of motion are reduced to an autonomous system of first order differential equations. For the case of the pure Gauss-Bonnet model two exact solutions with power-law and exponential dependence of scale factors (w.r.t. synchronous-like variable) are obtained. In cosmological case the first solution was obtained earlier by A. Toporensky and P. Tretyakov for n = 5 (and verified for n = 6,7). A generalization of the effective Lagrangian to the Lowelock case is conjectured. ...
A sum-over-histories generalized quantum theory is developed for homogeneous minisuperspace type A Bianchi cosmological models, focusing on the particular example of the classically recollapsing Bianchi type-IX universe. The decoherence functional for such universes is exhibited. We show how the probabilities of decoherent sets of alternative, coarse-grained histories of these model universes can be calculated. We consider in particular the probabilities for classical evolution defined by a suitable coarse graining. For a restricted class of initial conditions and coarse grainings we exhibit the approximate decoherence of alternative histories in which the universe behaves classically and those in which it does not. For these situations we show that the probability is near unity for the universe to recontract classically if it expands classically. We also determine the relative probabilities of quasiclassical trajectories for initial states of WKB form, recovering ...
Some properties of the universe are fixed by physics derived from mathematical symmetries, others may have been selected from an ensemble of possibilities. Some successes and failures of anthropic reasoning in this context are reviewed in the light of recent developments in astrobiology, cosmology and unification physics. Specific issues raised include our spacetime location (including the reason for the present age of the universe), the timescale of biological evolution, the tuning of global cosmological parameters, the origin of the Large Numbers of astrophysics, and the parameters of the Standard Model. Out of the twenty parameters of the Standard Model,the basic behavior and structures of the world (nucleons, nuclei,atoms, molecules, planets, stars, galaxies) depend mainly on five of them: $m_e,m_u,m_d,\\alpha,\\alpha_G$, three of which are independent in the context of Grand Unified Theories (that is, not related by any known symmetry). ...
We follow the evolution of the curvaton on superhorizon scales and check that the spectral tilt of the curvaton perturbations is unchanged as the curvaton becomes non-relativistic. Both inflation and pre-big bang cosmology can be treated since the curvaton mechanism within the two scenarios works the same way. We also discuss the amplitude of the density perturbations, which leads to some interesting constrains on the pre-big bang scenario. It is shown that within a SL(3,R) non-linear sigma model one of the three axions has the right coupling to the dilaton and moduli to yield a flat spectrum with a high string scale, if a quadratic non-perturbative potential is generated and an intermediate string phase lasts long enough.
Peculiar velocities induce apparent line of sight displacements of galaxies in redshift space, distorting the pattern of clustering in the radial versus transverse directions. On large scales, the amplitude of the distortion yields a measure of the dimensionless linear growth rate \\ff of fluctuations, which is related to the cosmological density \\Omega and the linear bias factor b in linearly biassed standard cosmology by \\ff \\approx \\Omega^{0.6} /b. To make the maximum statistical use of the data in a wide angle redshift survey, and for the greatest accuracy, the spherical character of the distortion needs to be treated properly, rather than in the simpler plane parallel approximation. In the linear regime, the redshift space correlation function is described by a spherical distortion operator acting on the true correlation function. It is pointed out here that there exists an operator, which is essentially the logarithmic derivative with ...
The existence of a cosmic neutrino background - the analogue of the cosmic microwave background - is a fundamental prediction of standard big bang cosmology. Up to now, the observational evidence for its existence is rather indirect and rests entirely on cosmological observations of, e.g., the light elemental abundances, the anisotropies in the cosmic microwave background, and the large scale distribution of matter. Here, we review more direct, weak interaction based detection techniques for the cosmic neutrino background in the present epoch and in our local neighbourhood. We show that, with current technology, all proposals are still off by some orders of magnitude in sensitivity to lead to a guaranteed detection of the relic neutrinos. The most promising laboratory search, based on neutrino capture on beta decaying nuclei, may be done in future experiments designed to measure the neutrino mass through decay kinematics. (orig.)
In this work we investigate the multivariate statistical description of the matter distribution in the nonlinear regime. We introduce the multivariate Edgeworth expansion of the lognormal distribution to model the cosmological matter field. Such a technique could be useful to generate and reconstruct three-dimensional nonlinear cosmological density fields with the information of higher order correlation functions. We explicitly calculate the expansion up to third order in perturbation theory making use of the multivariate Hermite polynomials up to sixth order. The probability distribution function for the matter field includes at this level the two-point, the three-point and the four-point correlation functions. We use the hierarchical model to formulate the higher order correlation functions based on combinations of the two-point correlation function. This permits us to find compact expressions for the skewness and kurtosis terms of the ...
The Kosambi-Cartan-Chern (KCC) theory represents a powerful mathematical method for the analysis of dynamical systems. In this approach one describes the evolution of a dynamical system in geometric terms, by considering it as a geodesic in a Finsler space. By associating a non-linear connection and a Berwald type connection to the dynamical system, five geometrical invariants are obtained, with the second invariant giving the Jacobi stability of the system. The Jacobi (in)stability is a natural generalization of the (in)stability of the geodesic flow on a differentiable manifold endowed with a metric (Riemannian or Finslerian) to the non-metric setting. In the present paper we review the basic mathematical formalism of the KCC theory, and present some specific applications of this method in general relativity, cosmology and astrophysics. In particular we investigate the Jacobi stability of the general relativistic static fluid sphere with a linear barotropic ...
Within the framework of the Bianchi type-IX homogeneous space, we set up a system of coupled equations for the cosmic scale factors, scalar field, and radiative energy density. At the tree level, the equations are written in a self-consistent, Hartree-Fock form. For phi/sup 4/ theory, the system of nine first-order differential equations is solved numerically for a varying ratio of the energy of anisotropy to the vacuum energy. As the vacuum energy increases, there appears to be less reheating, since the energy of anisotropy is more efficiently converted into isotropic expansion. If the energy of anisotropy is large enough, the inflationary phase is prevented. In this case, a series of cosmological phase transitions will take place each time the square of the effective mass changes its sign.
Within the framework of the Bianchi type-IX homogeneous space, we set up a system of coupled equations for the cosmic scale factors, scalar field, and radiative energy density. At the tree level, the equations are written in a self-consistent, Hartree-Fock form. For phi"4 theory, the system of nine first-order differential equations is solved numerically for a varying ratio of the energy of anisotropy to the vacuum energy. As the vacuum energy increases, there appears to be less reheating, since the energy of anisotropy is more efficiently converted into isotropic expansion. If the energy of anisotropy is large enough, the inflationary phase is prevented. In this case, a series of cosmological phase transitions will take place each time the square of the effective mass changes its sign.
A new mathematical framework is formulated to derive the effective equations of motion for the constrained quantum system which possesses an internal clock. In the realm close to classical behavior, the quantum evolution is approximated by a finite system of coupled but ordinary differential equations adhered to the weakly imposed Hamiltonian constraint. For the simplified version of loop quantum cosmology in the Bianchi I model with a free massless scalar filed, the resulting effective equations of motion affirm the bouncing scenario predicted by the previous studies: The big bang singularity is resolved and replaced by the big bounces, which take place up to three times, once in each diagonal direction, whenever the directional density approaches the critical value in the regime of Planckian density. It is also revealed that back-reaction arises from the quantum corrections and modifies the precise value of the directional density at the bouncing epoch. ...
Model selection aims to determine which theoretical models are most plausible given some data, without necessarily asking about the preferred values of the model parameters. A common model selection question is to ask when new data require introduction of an additional parameter, describing a newly-discovered physical effect. We review several model selection statistics, and then focus on use of the Bayesian evidence, which implements the usual Bayesian analysis framework at the level of models rather than parameters. We describe our CosmoNest code, which is the first computationally-efficient implementation of Bayesian model selection in a cosmological context. We apply it to recent WMAP satellite data, examining the need for a perturbation spectral index differing from the scale-invariant (Harrison-Zel'dovich) case.
In Randall-Sundrum-type brane-world cosmologies, the dynamical equations on the three-brane differ from the general relativity equations by terms that carry the effects of imbedding and of the free gravitational field in the five-dimensional bulk. Instead of starting from an ansatz for the metric, we derive the covariant nonlinear dynamical equations for the gravitational and matter fields on the brane. The local energy-momentum corrections are significant only at very high energies, and in this regime we show that fluid world-lines have a non-gravitational acceleration off the brane. The imprint on the brane of the nonlocal gravitational field in the bulk is more subtle, and we provide a careful decomposition of this effect. The nonlocal energy density determines the tidal acceleration in the off-brane direction, and can oppose singularity formation via the generalized Raychaudhuri equation. Unlike the nonlocal energy density and flux, the nonlocal anisotropic ...
The relative merits and difficlties of the primary and secondry origin hypotheses for the observed cosmic-ray antiprotons, including the new low-energy measurement of Buffington, et al are discussed. We conclude that the cosmic-ray antiproton data may be evidence for antimatter galaxies and baryon symmetric cosmology. The present bar P data are consistent with a primary extragalactic component having /p /equiv 1+/- 3.2/0.7x10 to the -4 independent of energy. We propose that the primary extragalactic cosmic ray antiprotons are most likely from active galaxies and that expected disintegration of bar alpha/alpha ban alpha/alpha. We further predict a value for ban alpha/alpha /equiv 10 to the -5, within range of future cosmic ray detectors.
We investigate the extent to which the pure magnification effect of gravitational lensing can be extracted from galaxy clustering statistics, by a nulling method which aims to eliminate terms arising from the intrinsic clustering of galaxies. The aim is to leave statistics which are free from the uncertainties of galaxy bias. We find that nulling can be done effectively, leaving data which are relatively insensitive to uncertainties in galaxy bias and its evolution, leading to cosmological parameter estimation which is effectively unbiased. This advantage comes at the expense of increased statistical errors, which are in some cases large, but it offers a robust alternative analysis method to cosmic shear for cosmological imaging surveys designed for weak lensing studies, or to full modelling of the clustering signal including magnification effects.
We investigate the asymptotic properties of the Bianchi type IX cosmological model in the brane-world scenario. The matter content is assumed to be a combination of a perfect fluid and a minimimally coupled scalar field that is restricted to the brane. A detailed qualitative analysis of the Bianchi type IX brane-world containing a scalar field having an exponential potential is undertaken. It is found that the brane-Robertson-Walker solution is a local source for the expanding Bianchi type IX models, and if k"2<2 the ever-expanding Bianchi IX models asymptote to the power-law inflationary solution. The only other local sink is the contracting brane-Robertson Walker solution. An analysis of the Bianchi type IX models with a scalar field with a general potential is discussed, and it is shown that in the case of expanding models, for physical scalar field potentials close to the initial singularity, the scalar field is effectively massless, and the solution is ...
We use the bulge Sersic index n and bulge-to-total mass ratio (B/T) to explore the fundamental question of how bulges form. We perform two-dimensional bulge-disk-bar decomposition on H-band images of 143 bright, high-mass (M _* #>=# 1.0 x 10"1"0 M _s_u_n) low-to-moderately inclined (i < 70"0) spirals. Our results are as follows. (1) Our H-band bar fraction (#approx#58%) is consistent with that from ellipse fits. (2) 70% of the stellar mass is in disks, 10% in bars, and 20% in bulges. (3) A large fraction (#approx#69%) of bright spirals have B/T#<=# 0.2, and #approx#76% have low n #<=# 2 bulges. These bulges exist in barred and unbarred galaxies across a wide range of Hubble types. (4) About 65% (68%) of bright spirals with n #<=# 2 (B/T #<=# 0.2) bulges host bars, suggesting a possible link between bars and bulges. (5) We compare the results with predictions from a set of #LAMBDA#CDM models. In the models, a high-mass spiral can have a bulge with ...
The propagating of tachyons in an expanding universe is discussed. It is shown that a primordial tachyon in the big-bang universe cannot survive unless it had very large energy initially. In an indefinitely expanding universe the tachyon trajectory turns back in time. This time barrier is found to exist even in the quantum mechanical discussion of tachyons. This property is used to set limits on the mass of a tachyon. The possible astronomical checks on the hypothesis that neutrinos or photons may be tachyonic are also discussed. (author).
The propagation of tachyons in an expanding universe is discussed. It is shown that a primordial tachyon in the big bang universe cannot survive unless it had very large energy initially. In an indefinitely expanding universe the tachyon trajectory turns back in time. This time barrier is found to exist even in the quantum mechanical discussion of tachyons. This property is used to set limits on the mass of a tachyon. The possible astronomical checks on the hypothesis that neutrinos or photons may be tachyonic are also discussed.
We demonstrate the existence of shear-free cosmological models with rotation and expansion which support inflationary scenarios. The corresponding metrics belong to the family of spatially homogeneous models with the geometry of the closed universe (Bianchi type IX). We show that the global vorticity does not prevent inflation and can even accelerate it.
In this thesis we study the relic density n{sub {chi}} of non--relativistic long--lived or stable particles {chi} in various non--standard cosmological scenarios. First, we discuss the relic density in the non--standard cosmological scenario in which the temperature is too low for the particles {chi} to achieve full chemical equilibrium. We also investigated the case where {chi} particles are non--thermally produced from the decay of heavier particles in addition to the usual thermal production. In low temperature scenario, we calculate the relic abundance starting from arbitrary initial temperatures T{sub 0} of the radiation--dominated epoch and derive approximate solutions for the temperature dependence of the relic density which can accurately reproduces numerical results when full thermal equilibrium is not achieved. If full equilibrium is reached, our ansatz no longer reproduces the correct temperature dependence of the {chi} number ...
The baryogenesis process in the early hot universe is investigated by means of relativistic kinetic theory. An exact solution to the kinetic equations for supermassive bosons serves to refine previous results: the optimum baryon-production domain is now complemented by bosons of low mass, thus removing the cosmological lower bound that had limited the mass of superheavy bosons. 14 references.
It is proven that there do not exist any Bianchi type-IX universes which expand for an infinite time, provided only that the matter satisfies the dominant energy condition and has non-negative average pressure (i.e., a non-negative trace of the spatial projection of the stress-energy tensor). This generalizes a proof recently given by us for the case of diagonal Bianchi type-IX solutions.
It is proven that there do not exist any Bianchi type-IX universes which expand for an infinite time, provided only that the matter satisfies the dominant energy condition and has non-negative average pressure (i.e., a non-negative trace of the spatial projection of the stress-energy tensor). This generalizes a proof recently given by us for the case of diagonal Bianchi type-IX solutions.
We study the dynamics of the scalar tensor cosmological model in the presence of tachyon field. In an alternative approach, in two exponential and power law form of the scalar field functions in the model, field equations are solved by simultaneously best fitting the model parameters with the most recent observational data. This approach gives us an observationally verified interpretation of the dynamics of the universe. We then discuss the best fitted of equation of state parameter, the statefinder parameters and the reconstructed scalar field in the model.
It is shown that observations of the lithium isotope ratio in high surface temperature Population II stars may be critical to cosmological nucleosynthesis models. In particular, decaying particle scenarios as derived in some supersymmetric models may stand or fall with such observations. 18 references.
We analyze the behavior of some Bianchi type-IX, mixmaster cosmological models possessing extra dimensions. We find that, unlike the three-dimensional case, they do not exhibit chaotic behavior on approach to their initial singularity. A finite sequence of stochastic mixmaster oscillations is, in general, replaced by monotonic contraction of the three-space scale factors on approach to a singularity whenever additional spatial dimensions exist with the metric form we consider.
We analyze the behavior of some Bianchi type-IX, mixmaster cosmological models possessing extra dimensions. We find that, unlike the three-dimensional case, they do not exhibit chaotic behavior on approach to their initial singularity. A finite sequence of stochastic mixmaster oscillations is, in general, replaced by monotonic contraction of the three-space scale factors on approach to a singularity whenever additional spatial dimensions exist with the metric form we consider.
We apply Hawking's proposal for the initial state of the Universe to a rotating Bianchi type-IX minisuperspace model in order to examine the problem of rotation in a closed cosmology. Using the Feynman path-integral techniques instead of the Wheeler-DeWitt equation, we are able to say that for a small angular velocity the closed Friedmann-Robertson-Walker model is indeed the most probable initial state for the Universe.
We apply Hawking's proposal for the initial state of the Universe to a rotating Bianchi type-IX minisuperspace model in order to examine the problem of rotation in a closed cosmology. Using the Feynman path-integral techniques instead of the Wheeler-DeWitt equation, we are able to say that for a small angular velocity the closed Friedmann-Robertson-Walker model is indeed the most probable initial state for the Universe.
Spatially homogeneous space-times of Bianchi type IX are considered. A general scheme for the derivation of exact solutions of Einstein's equations corresponding to perfect fluid plus pure radiation fields is outlined. Some simple rotating Bianchi type IX cosmological models are presented. The details of these solutions are also discussed. 9 refs. (author).
We consider a generalization of the DGP model, by adding a second brane with localized curvature, and allowing for a bulk cosmological constant and brane tensions. We study radion and graviton fluctuations in detail, enabling us to check for ghosts and tachyons. By tuning our parameters accordingly, we find bigravity models that are free from ghosts and tachyons. These models will lead to large distance modifications of gravity that could be observable in the near future.
The determination of relative abundances of helium and deuterium in Jupiter and Saturn provides estimates of the primordial abundances of these elements. The values so inferred can be compared with theoretical predictions of the standard model of the big-bang theory as a test of the model and suggest either that deuterium is more efficiently consumed in subsequent nucleosynthesis than hitherto assumed, or that the standard big-bang model should be revised in order for there to have been less production of primordial helium.
We argue that the occurrence of late-time acceleration can conveniently be described by first-order general relativity covariantly coupled to fermions. Dark energy arises as a gravitationally driven BCS condensate of fermions which forms in the early universe. At late times, the gap and chemical potential evolve to have an equation of state with effective negative pressure, thus naturally leading to acceleration.
The relative merits and difficulties of the primary and secondary origin hypotheses for the observed cosmic ray antiprotons, including the low energy measurement of Buffington, were examined. It is concluded that the cosmic ray antiproton data may be strong evidence for antimatter galaxies and baryon symmetric cosmology. The present antiproton data are consistent with a primary extragalactic component having antiproton/proton approximately equal to .0032 + or - 0.7.
In this short note we report our finding that within the context of alternative version of the Brans-Dicke theory (for #omega# #>=# -3/2, where #omega# is the Brans-Dicke parameter) the anisotropic Bianchi type cosmological models evolve towards the de Sitter isotropic universe. In short it is shown that during inflation there is no difference between the Brans-Dicke theory and General Relativity. Our result can thus be viewed as a generalization of the Wald's theorem for General Relativity. (author). 5 refs.
We study the high energy behaviour of Yang-Mills theory under the inclusion of gravity. In the weak-gravity limit, the running gauge coupling receives no contribution from the gravitational sector, if all symmetries are preserved. This holds true with and without cosmological constant. We also show that asymptotic freedom persists in general field-theory-based gravity scenarios including gravitational shielding as well as asymptotically safe gravity.
We present a set of exact solutions for quantum Bianchi type-IX anisotropic cosmological models (including the Taub model) of the form {Psi}={ital We}{sup {minus}{ital S}}. These solutions are spread over all values of anisotropy near the singularity, but at larger values of the radius of the universe they are strongly peaked around the {ital k}=+1 Friedmann-Robertson-Walker model.
We consider the pre-Big-Bang scenario with a massive axion field which starts to dominate energy density when oscillating in an instanton-induced potential and subsequently reheats the universe as it decays into photons, thus creating adiabatic CMB perturbations. We find that the fluctuations in the axion field can give rise to a nearly flat spectrum of adiabatic perturbations with a spectral tilt {delta}n in the range -0.1< or approx. {delta}n < or approx. 1.
We consider the problem of rotation in a homogeneous Bianchi type IX cosmological model. Studying the Wheeler-De Witt equation corresponding to this minisuperspace, and adopting a particular choice of the factor ordering, we are able to find a particular solution which is strongly peaked about isotropy. This result confirms all the previous investigations in this field, and suggests the conclusion that the machian nature of our universe is the natural, i.e. most likely, outcome of the Planck epoch. (orig.).
We consider the problem of rotation in a homogeneous Bianchi type IX cosmological model. Studying the Wheeler-De Witt equation corresponding to this minisuperspace, and adopting a particular choice of the factor ordering, we are able to find a particular solution which is strongly peaked about isotropy. This result confirms all the previous investigations in this field, and suggests the conclusion that the machian nature of our universe is the natural, i.e. most likely, outcome of the Planck epoch. (orig.).
This series of animations depicts the results of a computer simulation of the effects of deforestation on the surface temperature and precipitation in the Amazon region. Six animations are shown - a control simulation, a deforestation simulation, and the difference between the two for both surface temperature and rainfall. The period simulated is the second half of 1987.
... DARPA UUV and the UUV support simulator. ... Modular connections for robot arm manipulators ... Intelligent mobile robots will be performing many ...
In the next decade Type Ia supernovae (SNe Ia) will be used to test theories predicting changes in the Dark Energy equation of state with time. Ultimately this requires a dedicated space mission like JDEM. SNe Ia are mature cosmological probes --- their limitations are well characterized, and a path to improvement is clear. Dominant systematic errors include photometric calibration, selection effects, reddening, and population-dependent differences. Building on past lessons, well-controlled new surveys are poised to make strides in these areas: the Palomar Transient Factory, Skymapper, La Silla QUEST, Pan-STARRS, the Dark Energy Survey, LSST, and JDEM. They will obviate historical calibrations and selection biases, and allow comparisons via large subsamples. Some systematics follow from our ignorance of SN Ia progenitors, which there is hope of determining with SN Ia rate studies from 0cosmology, SNe Ia regulate galactic ...
We consider spatially homogeneous (but generally non-isotropic) cosmologies in the recently proposed Horava-Lifshitz gravity and compare them to those of general relativity using Hamiltonian methods. In all cases, the problem is described by an effective point particle moving in a potential well with exponentially steep walls. Focusing on the closed-space cosmological model (Bianchi type IX), the mixmaster dynamics is now completely dominated by the quadratic Cotton tensor potential term for a very small volume of the universe. Unlike general relativity, where the evolution toward the initial singularity always exhibits chaotic behavior with alternating Kasner epochs, the anisotropic universe in Horava-Lifshitz gravity (with parameter lambda > 1/3) is described by a particle moving in a frozen potential well with fixed (but arbitrary) energy E. Alternating Kasner epochs still provide a good description of the early universe for very ...
We consider spatially homogeneous (but generally non-isotropic) cosmologies in the recently proposed Horava-Lifshitz gravity and compare them to those of general relativity using Hamiltonian methods. In all cases, the problem is described by an effective point particle moving in a potential well with exponentially steep walls. Focusing on the closed-space cosmological model (Bianchi type IX), the mixmaster dynamics is now completely dominated by the quadratic Cotton tensor potential term for a very small volume of the universe. Unlike general relativity, where the evolution toward the initial singularity always exhibits chaotic behavior with alternating Kasner epochs, the anisotropic universe in Horava-Lifshitz gravity (with parameter #lambda# > 1/3) is described by a particle moving in a frozen potential well with fixed (but arbitrary) energy E. Alternating Kasner epochs still provide a good description of the early universe for very large ...
Weak lensing measurements are entering a precision era to statistically map the distribution of matter in the universe. The most common measurement has been of the variance of the projected surface density of matter, which corresponds to the induced correlation in alignments of background galaxies. This measurement of the fluctuations is insensitive to the total mass content, like using waves on the ocean to measure its depths. But when the depth is shallow as happens near a beach, waves become skewed. Similarly, a measurement of skewness in the projected matter distribution directly measures the total matter content of the universe. While skewness has already been convincingly detected, its constraint on cosmology is still weak. We address optimal analyses for the CFHT Legacy Survey in the presence of noise. We show that a compensated Gaussian filter with a width of 2.5 arc minutes optimizes the cosmological constraint, yielding $\\Delta ...
To systematically analyze the dynamical implications of the matter content in cosmology, we generalize earlier dynamical systems approaches so that perfect fluids with a general barotropic equation of state can be treated. We focus on locally rotationally symmetric Bianchi type IX and Kantowski-Sachs orthogonal perfect fluid models, since such models exhibit a particularly rich dynamical structure and also illustrate typical features of more general cases. For these models, we recast Einstein's field equations into a regular system on a compact state space, which is the basis for our analysis. We prove that models expand from a singularity and recollapse to a singularity when the perfect fluid satisfies the strong energy condition. When the matter source admits Einstein's static model, we present a comprehensive dynamical description, which includes the qualitative asymptotic behavior, of models in the neighborhood of the Einstein model; the results refute earlier ...
We re-examine the process of loop quantization for flat isotropic models in cosmology. In particular, we contrast different inequivalent `loop quantizations' of these simple models through their respective successes and limitations and assess whether they can lead to any viable physical description. We propose three simple requirements which any such admissible quantum model should satisfy: i) independence from any auxiliary structure, such as a fiducial interval/cell introduced to define the phase space when integrating over non-compact manifolds; ii) existence of a well defined classical limit and iii) provide a sensible "Planck scale" where quantum gravitational effects become manifest. We show that even when it may seem that one can have several possible loop quantizations, these physical requirements considerably narrow down the consistent choices. Apart for the so called improved dynamics of LQC, none of the other available inequivalent loop quantizations ...
The quantum behavior of the vacuum Bianchi type-IX universe with the cosmological constant is investigated in terms of the Ashtekar variables. An exact solution to the quantum Hamiltonian constraint in the holomorphic representation is given. This solution reduces to the Hartle-Hawking wave function in the spatially isotropic sector and extends in the triad representation to the classically forbidden region where the determinant of the spatial metric becomes negative. The analysis of the quantum Robertson-Walker universe indicates that if the superspace is extended to such a classically forbidden region, the holomorphic representation picks up some restricted class of solutions in general. This observation leads to a new ansatz on the boundary condition of the Universe. In particular, the behavior of the Lorentzian and Euclidean WKB orbits corresponding to the solution suggests a new picture on the semiclassical behavior of the quantum Universe: that the Universe ...
The quantum behavior of the vacuum Bianchi type-IX universe with the cosmological constant is investigated in terms of the Ashtekar variables. An exact solution to the quantum Hamiltonian constraint in the holomorphic representation is given. This solution reduces to the Hartle-Hawking wave function in the spatially isotropic sector and extends in the triad representation to the classically forbidden region where the determinant of the spatial metric becomes negative. The analysis of the quantum Robertson-Walker universe indicates that if the superspace is extended to such a classically forbidden region, the holomorphic representation picks up some restricted class of solutions in general. This observation leads to a new ansatz on the boundary condition of the Universe. In particular, the behavior of the Lorentzian and Euclidean WKB orbits corresponding to the solution suggests a new picture on the semiclassical behavior of the quantum Universe: that the Universe ...
Gravitational lensing by the dual cusp catastrophes of the cold dark matter (CDM) caustic rings at cosmological distances may provide the tantalizing opportunity to detect CDM indirectly, and discriminate between axions and weakly interacting massive particles (WIMPs). Caustics are places where the CDM particles are naturally focussed. Our focus is upon the caustic rings which are closed tubes whose cross-section is an elliptic umbilic catastrophe with three dual cusps. A caustic ring has a specific density profile, a specific geometry and, therefore, precisely calculable gravitational lensing signatures. The magnification monotonically increases as the line of sight approaches to the cusps where it diverges in the limit of zero velocity dispersion. In this limit, we find 37% magnification at a sample point near the outer cusps of the CDM caustic rings at cosmological distances. In the presence of finite velocity dispersion, the lower and upper ...
A diagrammatic approach to calculate n-point correlators of the primordial curvature perturbation \\zeta was developed a few years ago following the spirit of the Feynman rules in Quantum Field Theory. The methodology is very useful and time-saving, as it is for the case of the Feynman rules in the particle physics context, but, unfortunately, is not very well known by the cosmology community. In the present work, we extend such an approach in order to include not only scalar field perturbations as the generators of \\zeta, but also vector field perturbations. The purpose is twofold: first, we would like the diagrammatic approach (which we would call the Feynman-like rules) to become widespread among the cosmology community; second, we intend to give an easy tool to formulate any correlator of \\zeta for those cases that involve vector field perturbations and that, therefore, may generate prolonged stages of anisotropic expansion and/or ...
Some models for quantum gravity (QG) violate Lorentz invariance and predict an energy dependence of the speed of light, leading to a dispersion of high-energy gamma-ray signals that travel over cosmological distances. Limits on the dispersion from short-duration substructures observed in gamma-rays emitted by gamma-ray bursts (GRBs) at cosmological distances have provided interesting bounds on Lorentz invariance violation (LIV). Recent observations of unprecedentedly fast flares in the very-high energy gamma-ray emission of the active galactic nuclei (AGNs) Mkn 501 in 2005 and PKS 2155-304 in 2006 resulted in the most constraining limits on LIV from light-travel observations, approaching the Planck mass scale, at which QG effects are assumed to become important. I review the current status of LIV searches using GRBs and AGN flare events, and discuss limitations of light-travel time analyses and prospects for future instruments in the gamma-ray ...
The detailed formulation for loop quantum cosmology (LQC) in the Bianchi I model with a scalar massless field has been constructed. In this paper, its effective dynamics is studied in two improved strategies for implementing the LQC discreteness corrections. Both schemes show that the big bang is replaced by the big bounces, which take place up to three times, once in each diagonal direction, when the area or volume scale factor approaches the critical values in the Planck regime measured by the reference of the scalar field momentum. These two strategies give different evolutions: In one scheme, the effective dynamics is independent of the choice of the finite sized cell prescribed to make Hamiltonian finite; in the other, the effective dynamics reacts to the macroscopic scales introduced by the boundary conditions. Both schemes reveal interesting symmetries of scaling, which are reminiscent of the relational interpretation of quantum mechanics and also suggest ...
We propose to use alternative cosmic tracers to measure the dark energy equation of state and the matter content of the Universe [w(z) & Omega_m]. Our proposed method consists of two components: (a) tracing the Hubble relation using HII galaxies which can be detected up to very large redshifts, z~4, as an alternative to supernovae type Ia, and (b) measuring the clustering pattern of X-ray selected AGN at a median redshift of z~1. Each component of the method can in itself provide interesting constraints on the cosmological parameters, especially under our anticipation that we will reduce the corresponding random and systematic errors significantly. However, by joining their likelihood functions we will be able to put stringent cosmological constraints and break the known degeneracies between the dark energy equation of state (whether it is constant or variable) and the matter content of the universe and provide a powerful and alternative ...
We find instanton/cosmological solutions with biaxial Bianchi-IX symmetry, involving nontrivial spatial dependence of the CP{sup 1}-CP{sup 2}-sigma-models coupled to gravity. Such manifolds arise in N = 1, d 4 supergravity with supermatter actions and hence the solutions can be embedded in supergravity. There is a natural way in which the standard coordinates of these manifolds can be mapped into the four-dimensional physical space. Due to its special symmetry, we start with CP{sup 2} with its corresponding scalar ansatz; this further requires the spacetime to be SU(2) x U(1)-invariant. The problem then reduces to a set of ordinary differential equations whose analytical properties and solutions are discussed. Among the solutions there is a surprising, special family of exact solutions which owe their existence to the nontrivial topology of CP{sup 2} and are in 1-1 correspondence with matter-free Bianchi-IX metrics. These solutions can also be found by coupling ...
We find instanton/cosmological solutions with biaxial Bianchi-IX symmetry, involving nontrivial spatial dependence of the CP"1-CP"2-sigma-models coupled to gravity. Such manifolds arise in N = 1, d 4 supergravity with supermatter actions and hence the solutions can be embedded in supergravity. There is a natural way in which the standard coordinates of these manifolds can be mapped into the four-dimensional physical space. Due to its special symmetry, we start with CP"2 with its corresponding scalar ansatz; this further requires the spacetime to be SU(2) x U(1)-invariant. The problem then reduces to a set of ordinary differential equations whose analytical properties and solutions are discussed. Among the solutions there is a surprising, special family of exact solutions which owe their existence to the nontrivial topology of CP"2 and are in 1-1 correspondence with matter-free Bianchi-IX metrics. These solutions can also be found by coupling CP"1 to gravity. The ...
We propose to use alternative cosmic tracers to measure the dark energy equation of state and the matter content of the Universe [w(z) & \\Omega_m]. Our proposed method consists of two components: (a) tracing the Hubble relation using HII-like starburst galaxies, as an alternative to SNIa, which can be detected up to very large redshifts, z~4, and (b) measuring the clustering pattern of X-ray selected AGN at a median redshift of ~1. Each component of the method can in itself provide interesting constraints on the cosmological parameters, especially under our anticipation that we will reduce the corresponding random and systematic errors significantly. However, by joining their likelihood functions we will be able to put stringent cosmological constraints and break the known degeneracies between the dark energy equation of state (whether it is constant or variable) and the matter content of the universe and provide a powerful and alternative ...
A decreasing gravitational constant, G, coupled with angular momentum conservation is expected to increrase a planetary semimajor axis, a, as \\dot a/a=-\\dot G/G. Analysis of lunar laser ranging data strongly limits such temporal variations and constrains a local (~1 AU) scale expansion of the solar system as \\dot a/a=-\\dot G/G =-(4\\pm9)\\times10^{-13} yr^{-1}, including that due to cosmological effects.
A basic problem in gravitational physics is the resolution of spacetime singularities where general relativity breaks down. The simplest such singularities are conical singularities arising from orbifold identifications of flat space, and the most challenging are spacelike singularities inside black holes (and in cosmology). Topology changing processes also require evolution through classically singular spacetimes. I briefly review how a phase of closed string tachyon condensate replaces, and helps to resolve, basic singularities of each of these types. Finally I discuss some interesting features of singularities arising in the small volume limit of compact negatively curved spaces and the emerging zoology of spacelike singularities.
We study an analytical solution to the Einstein's equations in 2 + 1-dimensions. The space-time is dynamical and has a line symmetry. The matter content is a minimally coupled, massless, scalar field. Depending on the value of certain parameters, this solution represents three distinct space-times. The first one is at space-time. Then, we have a big bang model with a negative curvature scalar and a real scalar field. The last case is a big bang model with event horizons where the curvature scalar vanishes and the scalar field changes from real to purely imaginary. (author)
In this paper, the authors use free field techniques in D = 2 string theory t calculate the perturbation of the special state algebras when the cosmological constant is turned on. In particular, the authors find that the 'ground cone' preserved by the ring structure is promoted to a three-dimensional hyperboloid as conjectured by Witten. On the other hand, the perturbed (1,1) current algebra of moduli deformations is computed completely, and no simple geometrical interpretation is found. The authors also quote some facts concerning the Liouville matrix a model dictionary in this class of theories.
This paper reviews and elaborates on the issue of the dilaton transformation under the usual {tau} {yields} {alpha}{prime}/{tau} target space duality and its non-static generalization (or {sigma}-model duality). It is found that the transformation law {tau} {yields} {alpha}{prime}/{tau}, {phi} {yields} {phi} In ({tau}/{alpha}{prime}) which guarantees duality at the one-loop {sigma}-model level should be modified at two (and higher) loop order. The non-static duality is illustrated on the example of cosmological solutions in D {ge} 2 with time-dependent radii of space torus.
It is shown that all contracting, spatially homogeneous, orthogonal Bianchi cosmologies that are sourced by an ultra-stiff fluid with an arbitrary and, in general, varying equation of state asymptote to the spatially flat and isotropic universe in the neighbourhood of the big crunch singularity. This result is employed to investigate the asymptotic dynamics of a collapsing Bianchi type IX universe sourced by a scalar field rolling down a steep, negative exponential potential. A toroidally compactified version of M*-theory that leads to such a potential is discussed and it is shown that the isotropic attractor solution for a collapsing Bianchi type IX universe is supersymmetric when interpreted in an 11-dimensional context.
The diagonal and symmetric Bianchi type IX models are coupled to a homogeneous spinor field. An action for the combined fields is constructed, where the orthonormal basis used is given explicitly in terms of the metric. This allows one to vary the action with respect to the metric and the spinor fields only. Next, a Hamiltonian formulation is given, and a qualitative solution for the problem is presented. We also show that the k = +1 FRW (Friedmann--Robertson--Walker) model is not compatible with a homogeneous spinor field, while the more complicated models are.
The diagonal and symmetric Bianchi type IX models are coupled to a homogeneous spinor field. An action for the combined fields is constructed, where the orthonormal basis used is given explicitly in terms of the metric. This allows one to vary the action with respect to the metric and the spinor fields only. Next, a Hamiltonian formulation is given, and a qualitative solution for the problem is presented. We also show that the k = +1 FRW (Friedmann--Robertson--Walker) model is not compatible with a homogeneous spinor field, while the more complicated models are.
The authors draw the story of the concept of infinity in sciences from the ancient Greek civilization to modern times, from the Greek idea of atom to quantum physics and cosmology. It is shown how what was considered as limits in ancient times, have been progressively seen as unknown worlds, and are now being investigated with the use of big machines like particle accelerators. This book is composed of 3 main chapters: 1) the quest for the infinitely big, 2) the quest for the infinitely small, and 3) the big-bang theory that appears as the meeting point of the 2 infinites. (A.C.)
We investigate the dynamics of an 11-dimensional homogeneous cosmological model. We assume that the t = const hypersurfaces are products of a 3-dimensional Bianchi type-IX space and a 7-dimensional torus. Most results of our investigation hold when the 7-dimensional torus is replaced by an m-dimensional torus T/sup m/. We show that for a large class of vacuum solutions the physical space expands while the microspace contracts providing a natural mechanism of dimensional reduction. Matter satisfying a simple barotropic equation of state always breaks the process of dynamical dimensional reduction. With special attention we study the behavior of our model close to the initial singularity. In contrast with the 4-dimensional Bianchi type-IX cosmological model the Kasner solution always describes an approach to the initial singularity. We study the transition from the Kasner regime to the oscillatory regime. We show that matter does not ...
We investigate the dynamics of an 11-dimensional homogeneous cosmological model. We assume that the t = const hypersurfaces are products of a 3-dimensional Bianchi type-IX space and a 7-dimensional torus. Most results of our investigation hold when the 7-dimensional torus is replaced by an m-dimensional torus T/sup m/. We show that for a large class of vacuum solutions the physical space expands while the microspace contracts providing a natural mechanism of dimensional reduction. Matter satisfying a simple barotropic equation of state always breaks the process of dynamical dimensional reduction. With special attention we study the behavior of our model close to the initial singularity. In contrast with the 4-dimensional Bianchi type-IX cosmological model the Kasner solution always describes an approach to the initial singularity. We study the transition from the Kasner regime to the oscillatory regime. We show that matter does not ...
Our Universe is ruled by quantum mechanics and its extension quantum field theory. However, the explanations for a number of cosmological phenomena such as inflation, dark energy, symmetry breakings, and phase transitions need the presence of classical scalar fields. Although the process of condensation of scalar fields in the lab is fairly well understood, the extension of results to a cosmological context is not trivial. Here we investigate the formation of a condensate--a classical scalar field--after reheating of the Universe. We assume a light quantum scalar field produced by the decay of a heavy particle, which for simplicity is assumed to be another scalar. We show that during the radiation domination epoch under certain conditions, the decay of the heavy particle alone is sufficient for the production of a condensate. This process is very similar to preheating--the exponential particle production at the end of inflation. During the ...
The authors review two types of simulators for the analysis of cosmic ray neutron-induced soft errors (SE's). One of them is the neutron-induced soft error simulator (NISES). A recently proposed nuclear reaction theory forms the foundation for the nuclear reaction database used in NISES. The other simulator, the simplified simulator MBGR, is based on a modified version of the burst generation rate (BGR) model. Both simulators accurately simulate neutron-induced SE rates (SER's). MBGR actually provides an easier and quicker estimation of neutron-induced SER's than NISES. On the other hand, NISES covers more applications; it simulates neutron-induced charge collection, multiple-bit SE, and [alpha]-induced SE analysis.
The authors review two types of simulators for the analysis of cosmic ray neutron-induced soft errors (SE's). One of them is the neutron-induced soft error simulator (NISES). A recently proposed nuclear reaction theory forms the foundation for the nuclear reaction database used in NISES. The other simulator, the simplified simulator MBGR, is based on a modified version of the burst generation rate (BGR) model. Both simulators accurately simulate neutron-induced SE rates (SER's). MBGR actually provides an easier and quicker estimation of neutron-induced SER's than NISES. On the other hand, NISES covers more applications; it simulates neutron-induced charge collection, multiple-bit SE, and #alpha#-induced SE analysis
A special purpose hardware machine, which simulates up to one half-million gates and 2m byte RAM ICS at a 5 millisecond clock speed is described. This is accomplished with a hardware logic (HAL) simulator. This performance is achieved with 32 distributed special parallel processors, which utilize block oriented simulation technique. The technique promises a good cost hardware logic simulator. 7 references.
A computer program 'CPA simulator' has been designed to simulate the response parameters of a typical charge preamplifier (CPA) used in nuclear front-end electronics. The response parameters are studied considering three common connection methods between the detector and the CPA. The simulator predicts and illustrates the role and influence of the different components of the CPA and detector circuits.
A stochastic wind simulation for VAWTs (VSTOC) has been developed which yields turbulent wind-velocity fluctuations for rotationally sampled points. This allows three-component wind-velocity fluctuations to be simulated at specified nodal points on the wi...
The impacts of climate change on crop productivity are often assessed using simulations from a numerical climate model as an input to a crop simulation model. The precision of these predictions reflects...Full Text Available
Computer simulations of agricultural finance and management may be helpful to students preparing for agricultural careers. Simulations would provide problem definition and problem solving in realistic interactions with complex and ambiguous situations similar to those encountered in farming. The paper describes a major financial problem for farmers, explores the implications of that problem for students, and outlines a computer-based simulation model which would give students experience in dealing with the problem. Data sources for simulations, potentials of interactive exercises among students, and advantages of using real data with existing ambiguities are described. The board game "AGRICULTURE" is discussed as a prototype for simulation. (LFL)
This paper presents an approach to comparing computer run time of building simulation programs. The computing run time of a simulation program depends on several key factors, including the calculation algorithm and modeling capabilities of the program, the run period, the simulation time step, the complexity of the energy models, the run control settings, and the software and hardware configurations of the computer that is used to make the simulation runs. To demonstrate the approach, simulation runs are performed for several representative DOE-2.1E and EnergyPlus energy models. The computer run time of these energy models are then compared and analyzed.
Sep 1, 2011 ... 2005. Computer simulation modeling of recreation use: current status, case studies, and future directions. Gen. Tech. Rep. RMRS-GTR-143. ...
An application that provides a flexible and easy to use interface to the GAMMA spectral simulation package is described that is targeted at investigations using in vivo MR spectroscopic methods....Full Text Available
We give relations for the embedding of spatially-flat Friedmann-Robertson-Walker cosmological models of Einstein's theory in flat manifolds of the type used in Kaluza-Klein theory. We present embedding diagrams that depict different 4D universes as hypersurfaces in a higher-dimensional flat manifold. The morphology of the hypersurfaces is found to depend on the equation of state of the matter. The hypersurfaces possess a line-like curvature singularity infinitesimally close to the t 0{sup +} 3-surface, where t is the time expired since the big bang. The family of timelike comoving geodesics on any given hypersurface is found to have a caustic on the singular line, which we conclude is the 5D position of the point-like big bang.
In a first paper (Forni & Aghanim 1999), we developed several statistical discriminators to test the non-gaussian nature of a signal. These tests are based on the study of the coefficients in a wavelet decomposition basis. In this paper, we apply them in a cosmological context, to the study of the nature of the Cosmic Microwave Background (CMB) anisotropies. The latter represent the superposition of primary anisotropy imprints of the initial density perturbations and secondary ones due to photon interactions after recombination. In an inflationary scenario (standard Cold Dark Matter) with gaussian distributed fluctuations, we study the statistical signature of the secondary effects. More specifically, we investigate the dominant effects arising from the Compton scattering of CMB photons in ionised regions of the Universe: the Sunyaev-Zel'dovich effect of galaxy clusters and the effects of a spatially inhomogeneous re-ionisation of the Universe. Our study ...
The primary aims of the BESS-Polar program are precise measurements of the low-energy antiproton spectrum and search for cosmologically significant antimatter, which would provide new clues to understand the early Universe. The second flight (BESS-Polar II) over Antarctica was successfully carried out in December 2007 - January 2008. We performed 24.5 days scientific observation just at the solar minimum. The payload worked well during the flight and 4.7 billion cosmic-ray events were collected, which corresponds to 10-20 times statistics of the BESS data taken in the previous solar minimum period (1995 and 1997). Based on the BESS-Polar II data, we will present recent preliminary results of cosmic-ray antiproton measurements and sensitive search for antimatter. (authors)
We study the effects of quantum production of open strings on the relativistic scattering of D-branes. We find strong corrections to the brane trajectory from copious production of highly-excited open strings, whose typical oscillator level is proportional to the square of the rapidity. In the corrected trajectory, the branes rapidly coincide and remain trapped in a configuration with enhanced symmetry. This is a purely stringy effect which makes relativistic brane collisions exceptionally inelastic. We trace this effect to velocity-dependent corrections to the open string mass, which render open strings between relativistic D-branes surprisingly light. We observe that pair-creation of open strings could play an important role in cosmological scenarios in which branes approach each other at very high speeds. (author)
Relativistic strongly magnetized winds outflowing from fast-rotating compact objects like millisecond pulsars with surface magnetic fields of $\\sim 10^{15}-10^{16}$ G are plausible sources of cosmological $\\gamma$-ray bursts. In such winds, there are at least three regions where extremely powerful X-ray and $\\gamma$-ray emission may be generated. The first radiating region is the wind photosphere that is at a distance of $\\sim 10^9$ cm from the compact object. The second radiating region is at a distance of $\\sim 10^{13}-10^{14}$ cm. In this region, the striped component of the wind field is transformed into large-amplitude electromagnetic waves. The third radiating region is at a distance of $\\sim 10^{16}-10^{17}$ cm, where deceleration of the wind due to its interaction with an ambient medium becomes important. Radiation from all these regions is considered.
It is suggested that the model of a neutrino ball described in terms of the standard model extended by adding right-handed neutrinos and the Majorana scalar field can be presented in order to explain a body of weak interacting neutrinos. Neutrino interaction with the scalar Majorana field violates the lepton number and produces the mass splitting of neutrino due to the sea-saw mechanism. In this model a neutrino ball is an object which appears as a result of a first order cosmological phase transition. It can be regarded as a ball filled with Dirac neutrinos and can be treated as a remnant of the phase transition with unbroken global lepton symmetry. In this paper we study the macroscopic parameters of such a configuration. In the result the mass-radius curve M(R) for this object is obtained. (orig.).
This essay aims to summarize the main physical features arising from a new supersymmetric theory of gravitation. Based on preliminary discussions about classical field theory, cosmology, algebra and group theory, and taking formal results and theoretical considerations in comparison with several contributions from great authors, present work deals with gravity inside the limits of a meta-field theory, that is, a non-quantized but consistent representation of supergravity, the supersymmetry between gravitons and gravitinos. The introduction of meta-fields furnishes an independent framework for the study of gravity despite of constraints of quantization, treating the supersymmetric partners as deterministic actors of gravitation and not simply probabilistic entities. I explain my belief that gravitational field, by its own nature, is not quantizable in the same foot as the other fields, what does not means that we can not understand gravity by similar formal veins. ...
We review and extend in several directions recent results on the asymptotic safety approach to quantum gravity. The central issue in this approach is the search of a Fixed Point having suitable properties, and the tool that is used is a type of Wilsonian renormalization group equation. We begin by discussing various cutoff schemes, i.e. ways of implementing the Wilsonian cutoff procedure. We compare the beta functions of the gravitational couplings obtained with different schemes, studying first the contribution of matter fields and then the so-called Einstein-Hilbert truncation, where only the cosmological constant and Newton's constant are retained. In this context we make connection with old results, in particular we reproduce the results of the epsilon expansion and the perturbative one loop divergences. We then apply the Renormalization Group to higher derivative gravity. In the case of a general action quadratic in curvature we recover, within certain ...
ABSTRACT New observations from the Hubble Ultra Deep Field suggest that the star formation rate at Formula Not Shown drops off faster than previously thought. Using a newly determined star formation rate for the normal mode of Population II/I (PopII/I) stars, including this new constraint, we compute the Thomson scattering optical depth and find a result that is marginally consistent with Wilkinson Microwave Anisotropy Probe 5 results. We also reconsider the role of Population III (PopIII) stars in light of cosmological and stellar evolution constraints. While this input may be needed for reionization, we show that it is essential in order to account for cosmic chemical evolution in the early universe. We investigate the consequences of PopIII stars on the local metallicity distribution fu...
In this letter we discuss the connection between so-called homoclinic chaos and the violation of energy conditions in locally rotationally symmetric Bianchi type IX models, where the matter is assumed to be nontilted dust and a positive cosmological constant. We show that homoclinic chaos in these models is an artifact of unphysical assumptions: it requires that there exist solutions with positive matter energy density #rho#>0 that evolve through the singularity and beyond as solutions with negative matter energy density #rho#<0. Homoclinic chaos is absent when it is assumed that the dust particles always retain their positive mass. In addition, we discuss more general models: for solutions that are not locally rotationally symmetric we demonstrate that the construction of extensions through the singularity, which is required for homoclinic chaos, is not possible in general.
The behaviour of the wave function of the universe under the barrier for the anisotropic cosmological Bianchi type-IX model taking account of the scalar field is explored. In view of the known difficulties with the interpretation of multidimensional ones is offered. For this purpose in the frameworks of the semiclassical approach the system of characteristics equations relative to one variable is written out. This system describes a bundle of the characteristics along which the multidimensional problem is reduced to a one-dimensional one that allows to utilize the standard interpretation of the wave function as well as the usual Schroedinger equation. The obtained results for the Bianchi type-IX model are reduced to the following statement: the universe tunnels through the barrier from an isotropic state with small anisotropy that is necessary for providing a ling-lived inflation to derive the universe.
Disorder on the string theory landscape may significantly affect dynamics of eternal inflation leading to the possibility for some vacua on the landscape to become dynamically preferable over others. We systematically study effects of a generic disorder on the landscape starting by identifying a sector with built-in disorder -- a set of de Sitter vacua corresponding to compactifications of the Type IIB string theory on Calabi-Yau manifolds with a number of warped Klebanov-Strassler throats attached randomly to the bulk part of the Calabi-Yau. Further, we derive continuum limit of the vacuum dynamics equations on the landscape. Using methods of dynamical renormalization group we determine the late time behavior of the probability distribution for an observer to measure a given value of the cosmological constant. We find the diffusion of the probability distribution to significantly slow down in sectors of the landscape where the number of nearest neighboring vacua ...
The detection of a stochastic background of gravitational waves could significantly impact our understanding of the physical processes that shaped the early Universe. The challenge lies in separating the cosmological signal from other stochastic processes such as instrument noise and astrophysical foregrounds. One approach is to build two or more detectors and cross correlate their output, thereby enhancing the common gravitational wave signal relative to the uncorrelated instrument noise. When only one detector is available, as will likely be the case with the Laser Interferometer Space Antenna (LISA), alternative analysis techniques must be developed. Here we show that models of the noise and signal transfer functions can be used to tease apart the gravitational and instrument noise contributions. We discuss the role of gravitational wave insensitive "null channels" formed from particular combinations of the time delay interferometry, and derive a new combination ...
The relativistic time dilation is reviewed in a cosmological context. We show that a clock or twin paradox does not arise if cosmic time is properly taken into account. The receding galaxy background provides a unique frame of reference, and the proper times of geodesic as well as accelerated observers can be linked to the universal cosmic time parameter. This suggests to compare the proper time differentials of the respective observers by determining their state of motion in the galaxy grid. In this way, each observer can figure out whether his proper time is dilated or contracted relative to any other. In particular one can come to unambiguous conclusions on the aging of uniformly moving observers, without reference to asymmetries in measurement procedures or accelerations they may have undergone.
We consider perturbations of a particular type (homogeneous of Bianchi type IX) from a closed Friedmann-Robertson-Walker model. By considering a perfect fluid, we obtain closed-form solutions for the perturbations in three interesting cases: stiff matter, radiation, and dust. As a cosmological application, we obtain an upper limit on the relative shear, (sigma/theta)_0 < b/sub 0exp/, if the microwave background radiation was last scattered at a mean red-shift element of [7,10"3] and the density parameter is #OMEGA# element of (1,10), where b/sub 0exp/ is the relative root-mean-square deviation of the temperature over the whole sky.
The Arnowitt-Deser-Misner canonical formulation of general relativity is extended to the covariant brane-world theory in arbitrary dimensions. The exclusive probing of the extra dimensions makes a substantial difference, allowing for the construction of a non-constrained canonical theory. The quantum states of the brane-world geometry are defined by the Tomonaga-Schwinger equation, whose integrability conditions are determined by the classical perturbations of submanifolds contained in the Nash's differentiable embedding theorem. In principle, quantum brane-world theory can be tested by current experiments in astrophysics and by near future laboratory experiments at Tev energy. The implications to the black-hole information loss problem, to the accelerating cosmology, and to a quantum mathematical theory of four-sub manifolds are briefly commented.
A fully consistent linear perturbation theory for cosmology is derived in the presence of quantum corrections as they are suggested by properties of inverse volume operators in loop quantum gravity. The underlying constraints present a consistent deformation of the classical system, which shows that the discreteness in loop quantum gravity can be implemented in effective equations without spoiling space-time covariance. Nevertheless, non-trivial quantum corrections do arise in the constraint algebra. Since correction terms must appear in tightly controlled forms to avoid anomalies, detailed insights for the correct implementation of constraint operators can be gained. The procedures of this article thus provide a clear link between fundamental quantum gravity and phenomenology.
A fundamental Doppler-like but asymmetric wave effect that shifts received signals in frequency in proportion to their respective source distances, was recently described as means for a whole new generation of communication technology using angle and distance, potentially replacing TDM, FDM or CDMA, for multiplexing. It is equivalent to wave packet compression by scaling of time at the receiver, converting path-dependent phase into distance-dependent shifts, and can multiply the capacity of physical channels. The effect was hitherto unsuspected in physics, appears to be responsible for both the cosmological acceleration and the Pioneer 10/11 anomaly, and is exhibited in audio data. This paper discusses how it may be exploited for instant, passive ranging of signal sources, for verification, rescue and navigation; incoherent aperture synthesis for smaller, yet more accurate radars; universal immunity to jamming or interference; and precision frequency scaling of ...
A variety of forecasting approaches have been used in the energy area. Forecast simulations provide a common basis for comparing the relative precision of different forecasting methods. This paper presents a forecast simulation made using smoothing procedures for forecasting time series. The simulation design, statistical analysis, and conclusions are provided.
The paper sets up a simulation of a geothermal heat pump with a non-azeotropic mixture. The model is modified and verified with experimental data. The results of the simulation show that the systematic model can predict the performance within {+-}12% of the experimental data. As a result, some improvements can be provided on the basis of the simulation platform. (Author)
University research group with research areas: * Land based and submersible autonomous robots, (UUVs: AUVs and ROVs); * Controllers, electronics, sensor design and fusion, motion control; * Guidance and navigation of underwater vehicles; * AI, neural networks, fuzzy logic, subsumption control, behaviour based control; * Optical fibre and ultrasonic sensors for proximal object detection; * Robot arm control, visual servoing; * Imaging sonar applications; * Simulator development: UUV simulator; imaging sonar simulator; Aircraft/flight simulator.
In connection with the ERP project 'Dynamic modelling of staged gasification processes' a gasification simulator has been constructed. The simulator consists of: a mathematical model of the gasification process developed at Technical University of Denmark, a user interface programme, IGSS, and a communication interface between the two programmes. (BA)
2D simulation modeling of a river basin is based on notion of a River Basin as an opened non equilibrium nature system. All components of the system (elevation, water, soil properties and others) are linked to single multi - functional model. The objective of flood simulation is one the model functions. The objectives of the paper are as follows: a mean scale basin floods computing and the model validation; a large scale basin simulation and demands for its validation; results and demands for flood hazard assessment. (orig.)
A three-dimensional numerical simulation of mammatus-like clouds is presented. A portion of a cirrus outflow anvil cloud is simulated including cloud ice and snow microphysical representations. The simulated mammatus clouds appear in a cellular pattern and are compared with the few available previously published physical observations of mammatus. Copyright Copyright 2006 Royal Meteorological Society
The magnicon is a highly attractive candidate to be the RF source for a future multi-Tev linear collider. Physical models and computer codes have been developed which can provide start-to-end self-consistent simulations of a magnicon, including precise simulations of the high-convergence electron gun, RF-system, magnetic system, and beam collector. The 3-D beam dynamics simulations include realistic fields, finite beam size and transverse space charge effects. The codes allow one to provide steady-state simulations of the entire tube, so as to evaluate transient process of magnicon excitation, parasitic mode self-excitation, stability analysis, and tolerance analysis. The results of the simulations are found to be in good agreement with magnicon experiments. A brief description of the physical models and simulation codes employed will be given.
In this paper we present results of using parallel supercomputers to simulate beam dynamics in next-generation high intensity ion linacs. Our approach uses a three-dimensional space charge calculation with six types of boundary conditions. The simulations use a hybrid approach involving transfer maps to treat externally applied fields (including rf cavities) and parallel particle-in-cell techniques to treat the space-charge fields. The large-scale simulation results presented here represent a three order of magnitude improvement in simulation capability, in terms of problem size and speed of execution, compared with typical two-dimensional serial simulations. Specific examples will be presented, including simulation of the spallation neutron source (SNS) linac and the Low Energy Demonstrator Accelerator (LEDA) beam halo experiment.
;Contents: Comparison of Simulation Modules of TRANSYT and INTEGRATION Models; Evaluation of SCATSIM-RTA Adaptive Traffic Network Simulation Model; Comparison NETSIM, NETFLO I, and NETFLO II Traffic Simulation Models for Fixed-Time Signal Control; Traffic Flow Simulation Through Parallel Processing; Cluster Analysis as Tool in Traffic Engineering; Traffic Platoon Dispersion Modeling on Arterial Streets; Hybrid Model for Estimating Permitted Left-Turn Saturations Flow Rate; and Passing Sight Distance and Overtaking Dilemma on Two-Lane Roads.
The full scope simulator in Daya Bay Nuclear Power Station is imported from THOMSON-CSF company in France. It is a large facility, using digital simulation technology, that simulate the main control room of Daya Bay Nuclear Power Station with the high fidelity. It is used for training operators; system analyzing including operation procedure verification. The input output interface system of the simulator is the bridge of the information exchange between instrumentation and the computer. From the point of view of data transfer, this paper discusses the operation principle of the I/O interface system.
This paper describes a simulation model for predicting end-use energy consumption in residential sectors of a city or region. In this model, the annual energy consumption of a dwelling is simulated from the occupants' schedule of living activities, weather data and energy efficiencies of appliances and dwellings. By summing up the simulation results for various household categories, total energy consumption for the residential sector in a region can be estimated. In this paper, energy consumption for Osaka City is simulated. The result is compared with statistical data. The effects of energy efficiency standards and urban heat island phenomena are examined. (author)
The author describes the purpose, classification and main process of acceptance test of full scope simulator of Daya Bay Nuclear Power Plant, including the correction of non-conformance items which are discovered during the performance of acceptance tests. The results of the acceptance tests show that the model accepted by the full scope simulator of Daya Bay Nuclear Power Plant is fully able to cope with the simulation of normal transients and incidental transients and the performance of the simulator indeed compiled with the technical specifications which are defined n the relevant contracts.
We evaluate the neutrino fluxes to be expected from neutralino lightest supersymmetric particle (LSP) annihilations inside the Sun, within the minimal supersymmetric extension of the standard model with supersymmetry-breaking scalar and gaugino masses constrained to be universal at the grand unified theory scale [the constrained minimal supersymmetric standard model (CMSSM)]. We find that there are large regions of typical CMSSM (m_1_/_2,m_0) planes where the LSP density inside the Sun is not in equilibrium, so that the annihilation rate may be far below the capture rate. We show that neutrino fluxes are dependent on the solar model at the 20% level, and adopt the AGSS09 model of Serenelli et al. for our detailed studies. We find that there are large regions of the CMSSM (m_1_/_2,m_0) planes where the capture rate is not dominated by spin-dependent LSP-proton scattering, e.g., at large m_1_/_2 along the CMSSM coannihilation strip. We calculate neutrino fluxes above various threshold ...
The stability properties of a class of spacetimes with quasiregular singularities is discussed. Quasiregular singularities are the end points of incomplete, inextendible geodesics at which the Riemann tensor and its derivatives remain at least bounded in all parallel-propagated orthonormal (PPON) frames; observers approaching such a singularity would find that their world lines come to an end in a finite proper time. The Taub-NUT (Newman-Unti-Tamburino)-type cosmologies investigated are R/sup 1/ x T/sup 3/ and R/sup 3/ x S/sup 1/ flat Kasner spacetimes, the two-parameter family of spatially homogeneous but anisotropic Bianchi type-IX Taub-NUT spacetimes, and an infinite-dimensional family of Einstein-Rosen-Gowdy spacetimes studied by Moncrief. The behavior of matter near the quasiregular singularity in each of these spacetimes is explored through an examination of the behavior of the stress-energy tensors and scalars for conformally coupled and minimally coupled, ...
The stability properties of a class of spacetimes with quasiregular singularities is discussed. Quasiregular singularities are the end points of incomplete, inextendible geodesics at which the Riemann tensor and its derivatives remain at least bounded in all parallel-propagated orthonormal (PPON) frames; observers approaching such a singularity would find that their world lines come to an end in a finite proper time. The Taub-NUT (Newman-Unti-Tamburino)-type cosmologies investigated are R"1 x T"3 and R"3 x S"1 flat Kasner spacetimes, the two-parameter family of spatially homogeneous but anisotropic Bianchi type-IX Taub-NUT spacetimes, and an infinite-dimensional family of Einstein-Rosen-Gowdy spacetimes studied by Moncrief. The behavior of matter near the quasiregular singularity in each of these spacetimes is explored through an examination of the behavior of the stress-energy tensors and scalars for conformally coupled and minimally coupled, massive and massless ...
In this paper, numerical simulation of a concentric tube heat exchanger is presented to determine the convective heat transfer coefficient and friction factor in a smooth tube. Increasing the convective heat transfer coefficient can increase heat transfer rate in a concentric tube heat exchanger from a given tubular surface area. This can be achieved by using heat transfer augmentation devices. This work constitutes the initial phase of the numerical simulation of heat transfer from tubes employing augmentation devices, such as twisted tapes, wire-coil inserts, for heat transfer enhancement. A computational fluid dynamics (CFD) simulation tool was developed with CFX software and the results obtained from the simulations are validated with the empirical correlations for a smooth tube heat exchanger. The difficulties associated with the simulation of a heat exchanger augmented with ...
Details of the interaction of photons with tissue phantoms are elucidated using Monte Carlo simulations. In particular, photon sampling volumes and photon pathlengths are determined for a variety of scattering and absorption parameters. The Monte Carlo simulations are specifically designed to model light delivery and collection geometries relevant to clinical applications of optical biopsy techniques. The Monte Carlo simulations assume that light is delivered and collected by two, nearly-adjacent optical fibers and take into account the numerical aperture of the fibers as well as reflectance and refraction at interfaces between different media. To determine the validity of the Monte Carlo simulations for modeling the interactions between the photons and the tissue phantom in these geometries, the simulations were compared to measurements of aqueous suspensions of polystyrene ...
Field-grown snap bean plants were treated with simulated acidic rain applied either alone or in combination with ambient rain and the effects on growth and yield were determined. In plots where ambient rain was excluded, a retractable canopy was activated to shield the crop. Four levels of acidity at pH values of 5.0, 4.2, 3.4 and 2.6 were applied in four replicate treatments and the experiment was conducted in two successive years (1981 and 1982). In plots that received only simulated rain, yield was not adversely affected by acidic rain; in 1981, a positive linear relationship was present between acidity of simulated rain and yield, but in 1981, no effect was found. In contrast, in plots that received both simulated and ambient rain, a negative linear relationship between acidity in simulated rain and yield was observed in both years.
We have carried out numerical simulations of strongly gravitating systems based on the Einstein equations coupled to the relativistic hydrodynamic equations using adaptive mesh refinement (AMR) techniques. AMR coalescences of neutron stars can now be simulated with sufficient resolution covering the neutron stars while having the computational domain extend to the local wave zone. We show an AMR simulation carried out with a workstation having an accuracy equivalent to that of a 1025^3 regular uni-grid simulation, which is, to the best of our knowledge, larger than all previous simulations of similar NS systems on supercomputers. We believe the capability opens new possibilities in general relativistic simulations.
Although computer simulation has played a central role in the study of nucleation and growth since the earliest molecular dynamics simulations almost 50 years ago, confusion surrounding the effect of finite size on such simulations have limited their applicability. Modeling solidification in molten tantalum on the BlueGene/L computer, we report here on the first atomistic simulation of solidification that verifies independence from finite size effects during the entire nucleation and growth process, up to the onset of coarsening. We show that finite size scaling theory explains the observed maximal grain sizes for systems up to about 8,000,000 atoms. For larger simulations, a cross-over from finite size scaling to more physical size-independent behavior is observed.
The traditional process for developing large-scale simulations is cumbersome, time consuming, costly, and in some cases, inadequate. The topics of software components and component-based software engineering are being explored by software professionals in academic and industrial settings. A component is a well-delineated, relatively independent, and replaceable part of a software system that performs a specific function. Many researchers have addressed the potential to derive a component-based approach to simulations in general, and a few have focused on military simulations in particular. In a component-based approach, functional or logical blocks of the simulation entities are represented as coherent collections of components satisfying explicitly defined interface requirements. A simulation is a top-level aggregate comprised of a collection of components that interact with each ...
This paper describes TENT, a system designed for integrating the various tools and parts of a CFD simulation system (CFD - Computational Fluid Dynamics) and the simulation code itself in one interactive simulation environment. Initially the system was intended to be a TEstbed for Numerical Turbine simulation but it turned out that the approach was general enough to be suitable for more complex simulation tasks in the field of aircraft simulation. The paper will give a brief description of a CFD simulation system and its component interactions, which are the basis for describing the Component Architecture based on CORBA. This includes the control architecture using an event model and data flow considerations. Because most CFD codes use parallelization techniques based on Message Passing, a so called MIMD Object Abstraction is described in ...
Many simulation studies in biomedicine are based on a similar sequence of processing steps, starting from images and running through geometric model generation, assignment of tissue properties, numerical...Full Text Available
For the recruitment, training and position qualification of the simulator instructors and feedback of training effect, the management approaches are formulated in 'The System for Simulator Training and Licensing Examination of Daya Bay Nuclear Power Station Operators'. The concrete requirements on the professional knowledge, work experience and foreign language ability of a simulator instructor are put forward. The process of instructor training is designed. The training items include the trainer training, pedagogy training, time management training, operation activities training during outage of unit, 'shadow' training and on-the-jot training on simulator courses. Job rotation is realized between simulator instructor and licensing personnel on site. New simulator instructor must pass the qualification identification. After a duration of 2 years, ...
Solid molecular orthohydrogen exhibits orientational order at low temperatures. The orthohydrogen molecules, which are quadrupoles, order in the Pa3 structure. We have simulated this ordering, and explored the behaviour under dilution by spherical parahydrogen molecules.
We examined the effect of spatial iconicity (a perceptual simulation of canonical locations of objects) and word-order frequency on language processing and episodic memory of orientation. Participants...Full Text Available
Combined QM(PM3)/MM molecular dynamics simulations together with QM(DFT)/MM optimizations for key configurations have been performed to elucidate the enzymatic catalysis mechanism on the detoxification...Full Text Available
The percentages of total airflows over the nasal respiratory and olfactory epithelium of female rabbits were calculated from computational fluid dynamics (CFD) simulations of steady-state inhalation....Full Text Available
The Subseabed Disposal Program has recently completed a 30-day, 0.286-scale laboratory simulation of an In Situ Heat Transfer Experiment (ISHTE), scheduled for eventual deployment in the deep ocean 600 nautical miles north of Hawaii. 10 figures.
The Subseabed Disposal Program has recently completed a 30-day, 0.286-scale laboratory simulation of an In Situ Heat Transfer Experiment (ISHTE), scheduled for eventual deployment in the deep ocean 600 nautical miles north of Hawaii. 10 figures.
This paper overviews the general features of computer program 'HEC-5, Simulation of Flood Control and Conservation Systems', with emphasis on the capabilities of the most recent release of HEC-5, Version 7.2, dated March 1991. HEC-5 can simulate the essential features and operation goals and constraints of simple or complex systems with simulation intervals ranging from minutes to one month. Single event flood analysis and period of record conservation analysis may be accomplished with the model. Flood control analysis includes balanced system operation for downstream damage centers with consideration of forecasted local flows and hydrologic routing. In addition, induced surcharge operation based on spillway gate regulation schedules can be simulated. Hydropower analysis may include run-of-river, peaking, and pumped storage plants as well as system power operation. Water supply ...
The effect of feedback and materials on perceptual learning was examined in normal hearing listeners exposed to cochlear implant simulations. Generalization was most robust when feedback paired...Full Text Available
Purpose.The hypothesis that drivers with homonymous hemianopia (HH) would take a lane position that increased the safety margin on their blind side was tested with a driving simulator.Methods.Twelve...Full Text Available
A simulation is developed that qualitatively describes the small-zone-gel-filtration behaviour of a reversibly associating protein. The results reflect the dependence of the apparent molecular weight...Full Text Available
In recent years the modelling of interannual climate variability has been studied, the atmospheric energy and water cycles, and climate simulations with the ECHAM3 model. In addition, the climate simulations of several models have been compared with special emphasis in the area of northern Europe
... program incorporates the proven concepts of advanced distributed learning, modeling and simulation, and electronic digital libraries (Figure 1) to ...
A previous paper introduced the use of wave digital filters as a basic building block for power system simulation, particularly suitable for real-time applications. This paper stresses the simulation of non-linear and switching elements, emphasizing the advantages of the wave filters implementation. The digital structure is maintained even when non-linear components change their characteristics or power electronic devices switch their states. As a very important by-product, the suppression of numerical oscillations related to the trapezoidal rule is achieved in a rather simple way, with no effects on simulation results.
The interaction between a pair of counter-rotating vortices and a lean premixed CH{sub 4}/O{sub 2}/N{sub 2} flame ({Phi} = + 0.55) has been studied by direct numerical simulations using detailed and reduced chemical reaction schemes. Results from the complex chemistry simulation are discussed with respect to earlier experiments and differences in the simulations using detailed and reduces chemistry are investigated. Transient evolutions of the flame surface and the total heat release rate are compared and modifications in the evolution of the local flame structure are displayed. (authors) 22 refs.
This paper describes a signal simulator which is designed to calibrate eddy current probes. By using simultaneously an excitation coil and an electrically conducting plane, the simulator can emulate signals which correspond to spot welds with different flaws. The presence of a conducting plane allows a wide variation of the amplitude and phase of the excitation coil voltage during the calibration of a probe. A mathematical model for the simulator is derived and used to study the variation of the impedance change upon the parameters of the probe. Numerical computation shows that the impedance change depends in an important way on the frequency of the excitation coil current.
... used their Phoenics code to simulate the effect of additives, such as talcum powder which was impregnated in combustible cartridge cases, on the ...
... operation, experimental data reduction, and application of PES results to studies related to estuarine and coastal sediment transport problems are ...
... that defines hydrometeor categories for the duration of the simulation ... and to observe the movement of any particular system ... particle settling occurred ...
12], the Predictive Capability Maturity Model (PCMM) [ref. 13], the Uncertainty Structure scale [ref. 11], and the Simulation Readiness Level (SRL) scale ...
This paper reports that, to obtain better simulation results for a Canada deuterium uranium (CANDU) reactor operation, a new simulation method is developed that uses actual detector readings as a correction factor. Detector readings from a CANDU reactor are used to correct the calculated flux distribution during core calculation iterations. A suitable function is found to describe the relationship between the detector flux and the fluxes of mesh points around the detector. The new simulation method is tested by performing numerical calculations for the Wolsung reactor (a CANDU-600). The results show that the new method predicts the core state more accurately with fewer iterations.
the progressive simplification of composite sandwich panels. The equations ... The thermal and structural behavior of composite sandwich panels with a ...
In this paper we report on multi-physics, multi-billion macroparticle simulation of beam transport in a free electron laser (FEL) linac for future light source applications. The simulation includes a self-consistent calculation of 3D space-charge effects, short-range geometry wakefields, longitudinal coherent synchrotron radiation (CSR) wakefields, and detailed modeling of RF acceleration and focusing. We discuss the need for and the challenges associated with such large-scale simulation. Applications to the study of the microbunching instability in an FEL linac are also presented.
Most of the investigations have a theoretical sustenance based on molecular simulation. The area of application of molecular simulation is very wide, in the Materials Technology Department assigned to the Applied Sciences Management have been treated problems about metallic nano structures, glasses, interfaces, and molecules, to sustain and to explain some of the experimental results. Energy calculations are carried out to determine minimum energy structures, for later on to carry out calculations of some of their properties; as well as the images simulation of Electron microscopy and X-ray diffraction. (Author)
This document is a reference guide to the Xyce Parallel Electronic Simulator, and is a companion document to the Xyce Users Guide. The focus of this document is (to the extent possible) exhaustively list device parameters, solver options, parser options, and other usage details of Xyce. This document is not intended to be a tutorial. Users who are new to circuit simulation are better served by the Xyce Users Guide.
This document is a reference guide to the Xyce Parallel Electronic Simulator, and is a companion document to the Xyce Users' Guide. The focus of this document is (to the extent possible) exhaustively list device parameters, solver options, parser options, and other usage details of Xyce. This document is not intended to be a tutorial. Users who are new to circuit simulation are better served by the Xyce Users' Guide.
The utilization of electrically heated rods for the simulation of nuclear fuel rods represents an universally adopted method by the nuclear industry to study thermalhydraulic problems. The present work represents the development of a method to obtain the time variation of the electric linear power necessary to simulate a given nuclear power transient in order to yield the same temperature and heat flux conditions in the surface of the electrical heater that would be obtained by the nuclear fuel rod. (Author).
Simulation of the skim-off method for radon measurement by activated charcoal has been carried out by the finite elements method. The variation of radon concentration with peak appearance has been simulated. The thickness of the layer of charcoal removed after exposure and the time of peak appearance have been varied. For the charcoal and canister geometry used, the best results were found when the thickness of the layer removed was 1.72 mm.
2D MEDICI simulator is used to investigate hardening solutions to single-event burnout (SEB). SEB parametric dependencies such as carrier lifetime reduction, base enlargement, and emitter doping decrease have been verified and a p"+ plug modification approach for SEB hardening of power MOSFETs is validated with simulations on actual device structures.
Motivation is given for a technique to more thoroughly test semiconductor equipment control systems. A description is given of a simulator-based control system testing technique. Potential benefits that could be realized by using this technique in the semiconductor industry as well as benefits documented by using this technique in other industries are described. Specific requirements for using the technique in the semiconductor industry are outlined. A summary of a survey of nine commercial simulation systems is given. Finally, the outcome of the survey is compared with the requirements for using the technique.
A phantom simulation imaging quality control device is described that effectively simulates one centimeter lesions, using steel ball bearings as gamma ray attenuators. The bearings are mounted in a synthetic resinous sheet in an orthogonal pattern. The phantom can provide uniformity, resolution, linearity, distortion and field size checks, all with a single exposure.
The upgraded CDF II detector has collected first data during the initial operation of the Tevatron accelerator in Run II. The simulation of the CDF electromagnetic and hadronic central and upgraded plug (forward) calorimeter is based on the Gflash calorimeter parameterization package used within the GEANT based detector simulation of the Run II CDF detector. We present the results of tuning the central and plug calorimeter response to test beam data.
Accurate modeling of system response and scatter distribution is crucial for image reconstruction in emission tomography. Monte Carlo simulations are very well suited to calculate these quantities. However, Monte Carlo simulations are also slow and many simulated counts are needed to provide a sufficiently exact estimate of the detection probabilities. In order to overcome these problems, we propose to split the simulation into two parts, the detection system and the object to be imaged (the patient). A so-called 'virtual boundary' that separates these two parts is introduced. Within the patient, particles are simulated conventionally. Whenever a photon reaches the virtual boundary, its detection probability is calculated analytically by evaluating a multi-dimensional B-spline that depends on the photon position, direction and energy. The unknown B-spline knot values that define ...
The Finnish simulation code APROS and especially its 5-equation model is applied to modelling the horizontal steam generator. Different nodalizations are used in the secondary side of different models. Simulation results of the stationary state run are compared with results of RELAP5/MOD2 calculations and with an experimental plant data. (2 refs., 3 figs., 4 tabs.).
The Finnish simulation code APROS and especially its 5-equation model is applied to modelling the horizontal steam generator. Different nodalizations are used in the secondary side of different models. Simulation results of the stationary state run are compared with results of RELAP5/MOD2 calculations and with an experimental plant data. (2 refs., 3 figs., 4 tabs.).
Both the accuracy and real-time are quite strictly required for the RCS thermal-hydraulic model of the simulation of nuclear power station. So, a large amount of reasonable simplifications should be carried on. The author introduces the RCS thermal-hydraulic model of the full scope simulator of Daya Bay Nuclear Power Station, and emphasizes on the fundamental hypotheses and computation method of its kernel part--DEFI programme.
It is shown that a simple algorithm which exactly segregates between adiabatic and non-adiabatic electrons in particle-in-cell simulations of drift modes yields excellent conservation properties (e.g. particle number, energy) compared to the conventional df scheme. The removal of the free streaming term in the evolution of the marker weight is shown to be responsible for the improved linear and nonlinear properties of the simulated plasma.
This document describes the High Performance Electrical Modeling and Simulation (HPEMS) Global Verification Test Suite (VERTS). The VERTS is a regression test suite used for verification of the electrical circuit simulation codes currently being developed by the HPEMS code development team. This document contains descriptions of the Tier I test cases.
This document describes the High Performance Electrical Modeling and Simulation (HPEMS) Global Verification Test Suite (VERTS). The VERTS is a regression test suite used for verification of the electrical circuit simulation codes currently being developed by the HPEMS code development team. This document contains descriptions of the Tier I test cases.
The SciDAC2 accelerator project at SLAC aims to simulate an entire three-cryomodule radio frequency (RF) unit of the International Linear Collider (ILC) main Linac. Petascale computing resources supported by advances in Applied Mathematics (AM) and Computer Science (CS) and INCITE Program are essential to enable such very large-scale electromagnetic accelerator simulations required by the ILC Global Design Effort. This poster presents the recent advances and achievements in the areas of CS/AM through collaborations.
The mainly nonstationary operation of a solar-heated receiver can be simulated with sufficient accuracy only if data about the dynamic behavior are available. For this reason, the dynamic behavior of a solar cavity receiver with parabolic dish collector is investigated. The development of a mathematical simulation considering heat transfer and storage processes is presented and the procedure for a numerical solution is illustrated. The performance of the calculation method is finally demonstrated by simulating the passage of a cloud.
The mainly instationary operation of a solar heated receiver can be simulated with sufficient accuracy only if data about the dynamic behavior are available. For this reason, the dynamic behavior of a solar cavity receiver with parabolic dish collector is investigated. The development of a mathematical simulation considering heat transfer and storage processes is presented and the procedure for a numerical solution is illustrated. The performance of the calculation method is finally demonstrated by simulating the passage of a cloud.
The characteristics of martensitic transformations in solids are largely determined by the elastic strain that develops as martensite particles grow and interact. To study the development of microstructure, a finite-element computer simulation model was constructed to mimic the transformation process. The transformation is athermal and simulated at each incremental step by transforming the cell which maximizes the decrease in the free energy. To determine the free energy change, the elastic energy developed during martensite growth is calculated from the theory of linear elasticity for elastically homogeneous media, and updated as the transformation proceeds.
This demonstration project uses modern simulation techniques to illustrate the important technologies and design variables that an auto-designer would consider in production a high efficiency, low emissions vehicle. Simulation and modeling techniques use the idea of capturing the relationships between real components of the systems with mathematical equations. These equations are then solved on a computer to simulate the behavior or performance of the system under various conditions. In the current demonstration project, we focus on many variations of a hydrogen-powered vehicle.
Abstract This paper tries to investigate the problem of memory through one of its most intriguing patterns - chiasmus - reflected in old poetry, sacred texts, philosophy and theology, visual arts, as well as biology. It aims to search for some provisory explanation of why man was able once to acquire such excellence in memorizing internally thousands of lines of poetry, whereas now memory is expelled outside the human body and mind in a mere digital file. Contrasting the so-called "wish-dream for immortality" of the contemporary post-human body, this paper takes another path and looks into some old cosmologies and visions in which chiasmus constantly emerged as an enduring cultural paradigm with ontological relevance. Spanning from Plato to Christian theology, up to the contemporary Neo-Platonism of archetypal psychology, this paper hopes to put forth, if not a theory, at least a vision about man and cosmos, cosmos and man, a chiastic epiphany in which the body and ...
Observations of redshift-space distortions in spectroscopic galaxy surveys offer an attractive method for measuring the build-up of cosmological structure, which depends both on the expansion rate of the Universe and our theory of gravity. Galaxies occupy dark matter halos, whose redshift space clustering has a complex dependence on bias that cannot be inferred from the behavior of matter. We identify two distinct corrections on quasilinear scales (~ 30-80 Mpc/h): the non-linear mapping between real and redshift space positions, and the non-linear suppression of power in the velocity divergence field. We model the first non-perturbatively using the scale-dependent Gaussian streaming model, which we show is accurate at the 10 (s>25) Mpc/h for the monopole (quadrupole) halo correlation functions. We use perturbation theory to predict the real space pairwise halo velocity statistics. Our fully analytic model is accurate at the 2 per cent level only on scales s > ...
Magnetic fields correlated on several kiloparsec scales are seen in spiral galaxies. Their origin could be due to the winding up of a primordial cosmological field or due to amplification of a small seed field by a turbulent galactic dynamo. Both options have difficulties: There is no known battery mechanism for producing the required primordial field. Equally the turbulent dynamo may self destruct before being able to produce the large scale field, due to excess generation of small scale power. The current status of these difficulties is discussed. The resolution could depend on the nature of the saturated field produced by the small scale dynamo. We argue that the small scale fields do not fill most of the volume of the fluid and instead concentrate into intermittent ropes, with their peak value of order equipartition fields, and radii much smaller than their lengths. In this case these fields neither drain significant energy from the turbulence nor convert eddy ...
Towards the end of the 19th century, Kelvin pronounced as the "clouds of physics" 1) the failure of the Michelson-Morely experiment to detect an ether wind, 2) the violation of the classical mechanical equipartition theorem in statistical thermodynamics. And he believed that the removal of these clouds would bring physics to an end. But as we know, the removal of these clouds led to the two great breakthoughts of modern physics: 1) The theory of relativity, and 2) to quantum mechanics. Towards the end of the 20th century more clouds of physics became apparent. They are 1) the riddle of quantum gravity, 2) the superluminal quantum correlations, 3) the small cosmological constant. Furthermore, there is the riddle of dark energy making up 70% of the physical universe, the non-baryonic cold dark matter making up 26% and the very small initial entropy of the universe. An attempt is made to explain the importance of these clouds for the future of physics. Conjectures for ...
The origin of life and the origin of the universe represent two of the most important problems of science. Both are resolved by hydro-gravitational dynamics (HGD) cosmology (Gibson 1996, Schild 1996, Gibson 2009ab), which predicts frozen primordial hydrogen-helium gas planets in clumps as the dark matter of galaxies. Merging Earth-mass planets formed stars, moons and comets to incubate and cosmically seed the first life. Cometary panspermia (Hoyle and Wickramasinghe 1981, 1982; Wickramasinghe et al. 2009) occurs naturally by HGD mechanisms. Comets and moons are fragments from mergers of stardust covered frozen gas planets in their step-wise growth to star mass. Supernovae from stellar over-accretion of planets produce stardust (C, N, O, P etc.) chemical fertilizer. Planets collect this infected radioactive dust gravitationally, to provide liquid water domains in contact with life nutrients seeded with life prototypes. The first mutating, evolving, life from HGD ...
The matter spectrum of the MSSM, including three right-handed neutrino supermultiplets and one pair of Higgs-Higgs conjugate superfields, can be obtained by compactifying the E_{8} x E_{8} heterotic string and M-theory on Calabi-Yau manifolds with specific SU(4) vector bundles. These theories have the standard model gauge group augmented by an additional gauged U(1)_{B-L}. Their minimal content requires that the B-L gauge symmetry be spontaneously broken by a vacuum expectation value of at least one right-handed sneutrino. In previous papers, we presented the results of a quasi-analytic renormalization group analysis showing that B-L gauge symmetry is indeed radiatively broken with an appropriate B-L/electroweak hierarchy. In this paper, we extend these results by 1) enlarging the initial parameter space and 2) explicitly calculating all renormalization group equations numerically, without approximation. The regions of the initial parameter space leading to realistic vacua are ...
Seeing convolved two-dimensional de Vaucouleurs profiles have been fitted to digital images of the brightest elliptical galaxies in two nearby groups of galaxies and three clusters of galaxies with redshifts in the range 0.06 to 0.15. The digital images were created by scanning electronographic plates obtained by the use of several versions of the RGO electronographic camera. The linear relation between effective surface brightness #mu#/sub e/ and the logarithm of the effective radius r/sub e/, first discovered by Kormendy to be valid for nearby elliptical galaxies, is found to be equally valid for the largest and brightest galaxies in rich clusters of galaxies. The apparent universal nature of that relation indicates that the same type of processes might have been involved during the formation of all ellipticals from cD galaxies in clusters to fainter ellipticals in small groups. It is suggested that the (I/sub e/,r/sub e/) relation is not changed by the expected dynamical evolution ...
We study the sensitivity to neutrino masses of a Galactic supernova neutrino signal as could be measured with the detectors presently in operation and with future large volume water \\v{C}erencov and scintillator detectors. The analysis uses the full statistics of neutrino events. The method proposed uses the principles of Bayesian inference reasoning and has shown a remarkably independence of astrophysical assumptions. We show that, after accounting for the uncertainties in the detailed astrophysical description of the neutrino signal and taking into account the effects of neutrino oscillations in the supernova mantle, detectors presently in operation can have enough sensitivity to reveal a neutrino mass (or to set upper limits) at the level of 1 eV. This is sensibly better than present results from tritium $\\beta$-decay experiments, competitive with the most conservative limits from neutrinoless double $\\beta$-decay and less precise but remarkably less dependent from prior ...
Within the context of space-time (D-particle) foam in string/brane-theory it is demonstrated that it is possible to generate non-extensive statistics. The D-particle foam model involves point-like brane defects (D-particles), which provide the topologically non-trivial foamy structures of space-time. The D-particles can capture and emit stringy matter and this leads to a recoil of D-particles. It is indicated how one effect of such a recoil of D-particles is a back reaction on the space-time metric of Finsler type which is stochastic. We show that such a type of stochastic space-time foam can lead to cosmological effects similar to those induced by modifications of particle distributions within the framework of Tsallis entropies. The restrictions placed on the free parameters of the Finsler type metric are obtained from solving the Boltzmann equation in this background for relic abundances of a Lightest Supersymmetric Particle (LSP) dark matter candidate. It is ...
Measurement of redshifted 21-cm emission from neutral hydrogen promises to be the most effective method for studying the reionisation history of hydrogen and, indirectly, the first galaxies. These studies will be limited not by raw sensitivity to the signal, but rather, by bright foreground radiation from Galactic and extragalactic radio sources and the Galactic continuum. In addition, leakage due to gain errors and non-ideal feeds conspire to further contaminate low-frequency radio obsevations. This leakage leads to a portion of the complex linear polarisation signal finding its way into Stokes I, and inhibits the detection of the non-polarised cosmological signal from the epoch of reionisation. In this work, we show that rotation measure synthesis can be used to recover the signature of cosmic hydrogen reionisation in the presence of contamination by polarised foregrounds. To achieve this, we apply the rotation measure synthesis technique to the Stokes I ...
Primordial black holes are unique probes of cosmology, general relativity, quantum gravity and non standard particle physics. They can be considered as the ultimate particle accelerator in their last (explosive) moments since they are supposed to reach, very briefly, the Planck temperature. Upper limits on the primordial black hole number density of mass $M_{\\star} = 5 10^{14}$ g, the Hawking mass (born in the big-bang terminating their life presently), is determined comparing their predicted cumulative $\\gamma$-ray emission, galaxy-wise, to the one observed by the EGRET satellite, once corrected for non thermal $\\gamma$-ray background emission induced by cosmic ray protons and electrons interacting with light and matter in the Milky Way. A model with free gas emissivities is used to map the Galaxy in the 100 MeV photon range, where the peak of the primordial black hole emission is expected. The best gas emissivities and additional model parameters are obtained ...
The extraction of cosmological parameters from microwave background observations relies on specific assumptions about the statistical properties of the data, in particular that the p-point distributions of temperature fluctuations are jointly-normal. Using a battery of statistical tests, we assess the multivariate Gaussian nature of the Wilkinson Microwave Anisotropy Probe (WMAP) 1st year data. The statistics we use fall into three classes which test different aspects of joint-normality: the first set assess the normality of marginal (one-point) distributions using familiar univariate methods; the second involves statistics that directly assess joint-normality; and the third explores the evidence of non-linearity in the relationship between variates. We applied these tests to frequency maps, `foreground-cleaned' assembly maps and all-sky CMB-only maps. The assembly maps are of particular interest as when combined with the kp2 mask, we recreate the region used in ...
NGC 1407 is the central elliptical in a nearby evolved galaxy group apparently destined to become a cluster core. We use the kinematics of globular clusters to probe the dynamics and mass profile of the group's center, out to 60 kpc (~10 R_eff) -- the most extended data set to date around an early-type galaxy. This sample consists of 172 GC velocities, most of them newly obtained using Keck/DEIMOS, with a few additional objects identified as DGTOs or as IGCs. We find weak rotation in the GC system's outer parts, with the metal-poor and metal-rich GCs misaligned. The RMS velocity profile declines rapidly to a radius of ~20 kpc, and then becomes flat or rising to ~60 kpc. There is evidence that the GC orbits have a tangential bias that is strongest for the metal-poor GCs -- possibly contradicting theoretical expectations. We construct cosmologically-motivated galaxy+dark halo dynamical models and infer a mass within 60 kpc of ~3x10^12 M_Sun, extrapolating to a virial ...
The dynamics of a general Bianchi type-IX model with three scale factors is examined. The matter content of the model is assumed to be comoving dust plus a positive cosmological constant. The model presents a critical point of saddle-center-center type in the finite region of phase space. This critical point engenders in the phase space dynamics the topology of stable and unstable four dimensional tubes RxS"3, where R is a saddle direction and S"3 is the manifold of unstable periodic orbits in the center-center sector. A general characteristic of the dynamical flow is an oscillatory mode about orbits of an invariant plane of the dynamics which contains the critical point and a Friedmann-Robertson-Walker (FRW) singularity. We show that a pair of tubes (one stable, one unstable) emerging from the neighborhood of the critical point towards the FRW singularity have homoclinic transversal crossings. The homoclinic intersection manifold has topology RxS"2 and is ...
The perturbations of the L.R.S. class A spatially homogeneous spacetimes are treated using Hamiltonian methods in conjunction with techniques from the theory of Lie group harmonic analysis. These latter techniques lead to a simple way of handling any set of tensor equations on these background spacetimes which has the same symmetry group as the space-time metric. Once this approach is developed, the Hamiltonian formulation is used to recover in a clean way the Bonanos equations for the perturbations of the perfect fluid models of the class of spacetimes under consideration. The conserved quantities associated with the four-dimensional symmetry group are evaluated and their role in the linearized Hamiltonian dynamics is discussed. The time-dependent linear canonical transformation of the linearized vacuum gravitational phase space adapted to Moncrief's gauge invariant decomposition is described in general for these models and evaluated explicitly for a class of lower dimensional ...
Particle physics is driven by five great topics. Neutrino oscillations and masses are now at the fore. The standard model with extensions to supersymmetry and a Higgs to generate mass explains much of the field. The origins of CP violation are not understood. The possibility of extra dimensions has raised tantalizing new questions. A fifth topic lurking in the background is the possibility of something totally different. Many of the questions raised by these topics require powerful new accelerators. It is not an overstatement to say that for some of the issues, the accelerator is almost the experiment. Indeed some of the questions require machines beyond our present capability. As this volume attests, there are parts of the particle physics program that have been significantly advanced without the use of accelerators such as the subject of neutrino oscillations and many aspects of the particle-cosmology interface. At this stage in the development of physics, both ...
Very soon a new generation of reactor and accelerator neutrino oscillation experiments - Double Chooz, Daya Bay, Reno and T2K - will seek for oscillation signals generated by the mixing parameter theta_13. The knowledge of this angle is a fundamental milestone to optimize further experiments aimed at detecting CP violation in the neutrino sector. Leptonic CP violation is a key phenomenon that has profound implications in particle physics and cosmology but it is clearly out of reach for the aforementioned experiments. Since late 90's, a world-wide activity is in progress to design facilities that can access CP violation in neutrino oscillation and perform high precision measurements of the lepton counterpart of the Cabibbo-Kobayashi-Maskawa matrix. In this paper the status of these studies will be summarized, focusing on the options that are best suited to exploit existing European facilities (firstly CERN and the INFN Gran Sasso Laboratories) or technologies where ...
If high-redshift QSOs are ejected from the nuclei of low-redshift galaxies, as some have claimed, a large portion of their redshift must be intrinsic (non-Doppler). If these intrinsic components have preferred values, redshifts will tend to cluster around these preferred values and produce peaks in the redshift distribution. Doppler ejection and Hubble flow components will broaden each peak. Because ejection velocities are randomly directed and Hubble flow components are always positive, in this model all peaks are expected to show an asymmetry, extending further out in the red wing. If peaks are present showing this predicted asymmetry, it can lead directly to an estimate of quasar distances. Using two quasar samples, one with high redshifts and one with low, it is shown here that not only do all peaks in these two redshift distributions occur at previously predicted preferred values, they also all show the predicted extra extension in the red wing. For the low and high redshift ...
Modified gravity theories may provide an alternative to dark energy to explain cosmic acceleration. We argue that the observational program developed to test dark energy needs to be augmented to capture new tests of gravity on astrophysical scales. Several distinct signatures of gravity theories exist outside the linear regime, especially owing to the screening mechanism that operates inside halos like the Milky Way to ensure that gravity tests in the solar system are satisfied. This opens up several decades in length scale and new classes of galaxies at low-redshift that can be exploited by surveys. While theoretical work on models of gravity is in the early stages, we can already identify new regimes which cosmological surveys could target to test gravity. These include: 1. A small scale component that focuses on the interior and vicinity of galaxy and cluster halos. 2. Spectroscopy of low redshift galaxies, especially galaxies smaller than the Milky Way, in ...
We conjecture that all CP violations (both Dirac and Majorana types) arise from a common origin in neutrino seesaw. With this conceptually attractive and simple conjecture, we deduce that mu-tau breaking shares the common origin with all CP violations. We study the common origin of mu-tau and CP breaking in the Dirac mass matrix of seesaw Lagrangian (with right-handed neutrinos being mu-tau blind), which uniquely leads to inverted mass-ordering of light neutrinos. We then predict different correlations between the two small mu-tau breaking observables theta_{13} - 0 and theta_{23} -45, which can be tested against our previous normal mass-ordering scheme by the on-going oscillation experiments. We also analyze the correlations of theta_{13} with Jarlskog invariant and neutrinoless double-beta-decay observable. From the common origin of CP and mu-tau breaking in the neutrino seesaw, we establish a direct link between the low energy CP violations and the cosmological ...
Gamma-Ray Bursts (GRBs) are thought to originate at cosmological distances from the most powerful explosions in the Universe. If GRBs are not beamed then the distribution of their number as a function of gamma-ray flux implies that they occur once per (0.3-40) million years per bright galaxy and that they deposit >10^{53} ergs into their surrounding interstellar medium. The blast wave generated by a GRB explosion would be washed out by interstellar turbulence only after tens of millions of years when it finally slows down to a velocity of 10 km/s. This rather long lifetime implies that there could be up to several tens of active GRB remnants in each galaxy at any given time. For many years, radio observations have revealed the enigmatic presence of expanding neutral-hydrogen (HI) supershells of kpc radius in the Milky Way and in other nearby galaxies. The properties of some supershells cannot be easily explained in terms of conventional sources such as stellar ...
If gamma-ray bursts originate in galaxies at cosmological distances, the host galaxy should be detected if a burst error box is searched deep enough; are the host galaxies present? We present and implement a statistical methodology which evaluates whether the observed galaxy detections in a burst's error box are consistent with the presence of the host galaxy, or whether all the detections can be attributed to unrelated background galaxies. This methodology requires the model-dependent distribution of host galaxy fluxes. While our methodology was derived for galaxies in burst error boxes, it can be applied to other candidate host objects (e.g., active galaxies) and to other types of error boxes. As examples, we apply this methodology to two published studies of burst error boxes. We find that the nine error boxes observed by Larson and McLean (1997) are too large to discriminate between the presence or absence of host galaxies, while the absence of bright galaxies ...
We adapt the Jain-Taylor (2003) shear-ratio geometric lensing method to measure the dark energy equation of state, w = pv/?v and its time derivative from dark matter haloes in cosmologies with arbitrary spatial curvature. The full shear-ratio covariance matrix is calculated for lensed sources, including the intervening large-scale structure and photometric redshift errors as additional sources of noise, and a maximum likelihood method for applying the test is presented. Decomposing the lensing matter distribution into dark matter haloes we calculate the parameter covariance matrix for an arbitrary experiment. Combining with the expected results from the cosmic microwave background (CMB) we design an optimal survey for probing dark energy. This shows that a targeted survey imaging 60 of the largest clusters in a hemisphere with five-band optical photometric redshifts to a median galaxy depth of zm = 0.9 could measure w0 ? w(z = 0) to a marginal 1? error of ?w0 = ...
A distributed watershed hydrologic model soil and water assessment tool (SWAT, 2000 Version) was applied to simulate stream flow and nutrient loadings (including total nitrogen (TN) and total phosphorus (TP)) in Gucheng Lake Basin, which is located in the lower reach of Yangtze River Basin. Model parameters related to stream flow and nutrient loadings were calibrated and validated by the observed value, and the simulation showed that the simulated values were reasonably comparable to the observed data, suggesting the validity of SWAT model. The simulated surface runoff and nutrient yield between two stages, viz. 1981-2000 and 1951-1960, were compared, basically considering land reclamation in the mid-1970s. The result indicated that, the average annual runoff and output of TN and TP in 198...
During the last years, large-scale simulations of realistic physical environments which support the interaction of multiple participants over the Internet have become increasingly available and economically significant, most notably in the computer gaming industry. Such systems, commonly called networked virtual environments (NVEs), are usually based on a client-server architecture where for performance reasons and bandwidth restrictions, the simulation is partially deferred to the clients. This inevitable architectural choice renders the simulation vulnerable to attacks against the semantic integrity of the simulation: malicious clients may attempt to compromise the physical and logical laws governing the simulation, or to alter the causality of events a posteriori. In this paper, we initiate the systematic study of semantic integrity in NVEs from a security point of view. We argue ...
A compact proton accelerator for medical applications is being developed at Lawrence Livermore National Laboratory. The accelerator architecture is based on the dielectric wall accelerator (DWA) concept. One critical area to consider is the switch region. Electric field simulations and thermal calculations of the switch area were performed to help determine the operating limits of rmed SiC switches. Different geometries were considered for the field simulation including the shape of the thin Indium solder meniscus between the electrodes and SiC. Electric field simulations were also utilized to demonstrate how the field stress could be reduced. Both transient and steady steady-state thermal simulations were analyzed to find the average power capability of the switches.
This project is about the simulation design of an engine control unit (ECU) for an Otto cycle engine with electronic fuel injection (EFI). The simulation includes a model for the ECU as well as physical parameters of the engine, which allows closed-loop control and monitoring of various systems.This simulation has been realized using Simulink and Stateflow, which are components of Mathworks' MATLAB software.The program allows control of various parameters of the ECU, as well as the simulation of failures to verify that the designed ECU is fault-tolerant and can control the engine using an open loop control. The main function the ECU provides is fuel metering.Subsequently, this program could be used as a tool to quickly develop and test models of ECU in order to control an engine in laborat...
We review our recent work on an atomistic approach to the development of predictive process simulation tools. First principles methods, molecular dynamics simulations, and experimental results are used to construct a database of defect and dopant energetics in Si. This is used as input for kinetic Monte Carlo simulations. C and B trapping of the Si self- interstitial is shown to help explain the enormous disparity in its measured diffusivity. Excellent agreement is found between experiments and simulations of transient enhanced diffusion following 20-80 keV B implants into Si, and with those of 50 keV Si implants into complex B-doped structures. Our simulations predict novel behavior of the time evolution of the electrically active B fraction during annealing.
Piping in nuclear power plants is vital to the proper operation and safety of these facilities. To assure safety in the unlikely event of a pipe break, it is necessary to evaluate the consequences from the resulting whipping pipe on neighboring components and structures. Numerical simulations allow for rapid evaluation of these consequences. Before simulations can be accepted, however, the methodology and computer codes must be validated against experimental results. This paper uses a probabilistic approach to validate pipe whip simulations against limited experimental results. Probabilistic analysis software was developed and coupled to existing deterministic finite element software. An example of a whipping pipe impacting against a reinforced concrete slab was simulated. The described pr...
This paper presents a trial of probabilistic simulation for performance assessment of high-level waste (HLW) disposal using the Monte Carlo method. Based on the Reference Case conceptual model in H12, a new integrated simulation system which allowed us to rapidly evaluate the effect of data uncertainty was developed. The doses to hypothetical exposure group were compared with the results of H12 that were performed by using a point-wise approach, in terms of maximum total dose. This study showed that H12 results were consistent with results of probabilistic simulation and also showed that transmissivity had a strong influence for the uncertainty of the system performance in all simulation time. (author)
The software package VIBA-lab1, which incorporates PIXE and RBS energy spectra simulation has now been extended to include the simulation of elemental maps from 3D structures. VIBA-lab1 allows the user to define a wide variety of experimental parameters, e.g. energy and species of incident ions, excitation and detection geometry, etc. When the relevant experimental parameters as well as target composition are defined, the program can then simulate the corresponding PIXE and RBS spectra. VIBA-LAB2 has been written with applications in nuclear microscopy in mind. A set of drag-and-drop tools has been incorporated to allow the user to define a three-dimensional sample object of mixed elemental composition. PIXE energy spectra simulations are then carried out on pixel-by-pixel basis and the corresponding intensity distributions or elemental maps can be computed. Several simulated ...
The influence of the velocity profile at the inlet boundary on the simulation of air velocity distribution inside an electrostatic precipitator is presented in this study. Measurements and simulations were performed in a duct and an electrostatic precipitator (ESP). A four-hole cobra probe was used for the measurement of velocity distribution. The flow simulation was performed by using the computational fluid dynamics (CFD) code FLUENT. Numerical calculations for the air flow were carried out by solving the Reynolds-averaged Navier-Stokes equations coupled with the realizable k-{epsilon} turbulence model equations. Simulations were performed with two different velocity profiles at the inlet boundary - one with a uniform (ideal) velocity profile and the other with a non-uniform (real) velocity profile to demonstrate the effect of velocity inlet boundary condition on the flow ...
Maps of wind resources were used to develop a method of evaluating climatic changes and their potential impacts on wind energy resources. Global IPCC climate change scenarios were used to predict climatic conditions for the future, while past wind resource availability was simulated and validated using NCEP and NCAR reanalysis data as well as observed meteorological data from Environment Canada. The simulations were used to compare each 5-year period with a 50-year reference period. Regional scale climate change impacts were evaluated using a statistical dynamic down-scaling method. Advanced meteorological models were used to predict wind flow patterns across specific landscapes. The evolution of past wind resource availability was then simulated. Five-year wind resource simulations for a 50-year period were simulated at 25 km"2 wind speeds at 80 m above the ground. Average wind ...
We examine the observable properties of simulated barred galaxies including radial mass profiles, edge-on structure and kinematics, bar lengths and pattern speed evolution for detailed comparison to real systems. We have run several simulations in which bars are created through inherent instabilities in self-consistent simulations of a realistic disc+halo galaxy model with a disc-dominated, flat rotation curve. These simulations were run at high (N=20M particles) and low (N=500K) resolution to test numerical convergence. We determine the pattern speeds in simulations directly from the phase angle of the bar versus time and the Tremaine-Weinberg method. Fundamental dynamics do not change between the high and low resolution, suggesting that convergence has been reached in this case. We find the higher resolution is needed to simulate structural and kinematic ...
Piping in nuclear power plants is vital to the proper operation and safety of these facilities. To assure safety in the unlikely event of a pipe break, it is necessary to evaluate the consequences from the resulting whipping pipe on neighboring components and structures. Numerical simulations allow for rapid evaluation of these consequences. Before simulations can be accepted, however, the methodology and computer codes must be validated against experimental results. This paper uses a probabilistic approach to validate pipe whip simulations against limited experimental results. Probabilistic analysis software was developed and coupled to existing deterministic finite element software. An example of a whipping pipe impacting against a reinforced concrete slab was simulated. The described probabilistic approach was used to validate the numerical simulations. The conclusions obtained ...
A nuclear power plant engineering (NPP) simulator may include many system functions such as thermal-hydraulic calculation, 3D reactor neutron kinetics model, control and protection system, display and operational human-machine interface, intelligent alarm system, etc. If all those functions are simulated using a single code, the size of the program will cause structural or managerial problems, and the hardware requirement could be tremendous. Thus, any local function errors or changes may affect the whole code, which will make development and maintenance extremely costly. So the nondistributed code is not considered to be flexible and feasible. Currently, distributed and parallel technologies have been applied to develop NPPs engineering simulator for safety analysis, verification of advanced main control room display, operation and intelligent alarm design. The distributed simulation divides a huge ...
The plasma plume generated by a new type of high energy Janitor known as the railplug, is examined. The railplug is a miniaturized railgun that has the potential for improving ignition characteristics of combustible mixtures in engines. The objective of the study is to gain an uderstanding of the characteristics of the plasma created by a transparent railplug, and to validate a multidimensional computer simulation of the plasma and shock fronts. The nature of the plume emitted by the railplug was examined for three levels of electrical energy while firing into air at a pressure of 1 atm. The computer model is to be used to predict trends in railplug performance for various railplug designs, energies, and ambient conditions. The velocity of the plasma movement inside a transparent railplug was measured, as well as the velocity of the plume ejected from the cavity. A shock is produced at the initiation point of the arc and propagates down the cavity, eventually ...
The evolution of FRI jets has been long studied in the framework of the FRI-FRII dichotomy. In this paper, we test the present theoretical and observational models via a relativistic numerical simulation of the jets in the radio galaxy 3C 31. We use the parameters derived from the modelling presented by \\cite{lb02a,lb02b} as input parameters for the simulation of the evolution of the source, thus assuming that they have not varied over the lifetime of the source. We simulate about 10 % of the total lifetime of the jets in 3C 31. Realistic density and pressure gradients for the atmosphere are used. The simulation includes an equation of state for a two-component relativistic gas that allows a separate treatment of leptonic and baryonic matter. We compare our results with the modelling of the observational data of the source. Our results show that the bow shock evolves self-similarly at a quasi-constant ...
Ion implantation and subsequent annealing are essential stages in today's advanced CMOS processing. Although the dopant implanted profile can be accurately predicted by analytical fits calibrated with SIMS profiles, the damage has to be estimated with a binary collision approximation implant simulator. Some models have been proposed, like the '+n', in an attempt to simplify the anneal simulation. We have used the atomistic kinetic Monte Carlo dados to elucidate which are the implant modeling features most relevant in the simulation of transient enhanced diffusion (TED). For the experimental conditions studied we find that the spatial correlation of the I, V Frenkel pairs is not critical in order to yield the correct I supersaturation, that can be simulated just taking into account the net I-V excess distribution. In contrast to, simulate impurity clustering/deactivation when there ...
We present a dynamical spectral model for Large Eddy Simulation of the incompressible magnetohydrodynamic (MHD) equations based on the Eddy Damped Quasi Normal Markovian approximation. This model extends classical spectral Large Eddy Simulations for the Navier-Stokes equations to incorporate general (non Kolmogorovian) spectra as well as eddy noise. We derive the model for MHD and show that introducing a new eddy-damping time for the dynamics of spectral tensors in the absence of equipartition between the velocity and magnetic fields leads to better agreement with direct numerical simulations, an important point for dynamo computations.
Energy deposition process by relativistic fast electrons produced by ultra-intense laser pulses is discussed. The process is calculated with a two dimensional Fokker-Planck simulation code including binary and collective collisions coupled with electromagnetic field. We focused on Velocity Distribution Function (VDF) dependence in the simulation. The results show that the spread angle of the fast electrons distribution affects energy deposition area and deposited energy is concentrated in the vicinity of the propagation axis of the fast electrons. It may be also suggested that self-pinch effect of a fast electron beam causes large deposition energy. (author)
The statistical {gamma}-spectra of highly excited even-even rare earth nuclei are simulated applying appropriate level density and strength function to a given nucleus. Hindrance effects due to K-conservation are taken into account. Simulations are compared to experimental data from the {sup 163}Dy({sup 3}He,{alpha}){sup 162}Dy and {sup 173}Yb({sup 3}He,{alpha}){sup 172}Yb reactions. The influence of the K quantum number at higher energies is discussed. 21 refs., 7 figs., 2 tabs.
We investigate the formation of nanostructures in 2D strained alloys on face centered cubic (111) surfaces by means of equilibrium Monte Carlo simulations. In the framework of an off-lattice model, we consider one monolayer of two bulk-immiscible adsorbates A and B with negative and positive misfit relative to the substrate, respectively. Simulations show that the adsorbates partly self-organize into island or stripe-like patterns. We show how these structures depend on the relative misfits, interaction, and concentration of components. The morphology is quite different for phase separation and intermixing regimes.
Numerical simulation of dynamic fracture of concrete slabs, impacted by steel projectiles, was carried out in this study. The concrete response was described by a continuum damage model. This continuum damage model was originally developed to study rock fragmentation and was modified in the present study with an emphasis on the post-limit structural response. The model was implemented into a transient dynamic explicit finite element code LS-DYNA2D and the code was then used for the numerical simulations. The specific impact configuration of this study follows the experiment series conducted by Hanchak et al. Comparisons between calculated results and measured data were made. Good agreements were found.
Tests for the applicability of a CFD (Computational Fluid Dynamics) code for simulating activity transport in PNA (Pulsed Neutron Activation) fluid measurements have been performed. The CFD code was combined with a Monte Carlo code used for the calculation of the initial activity distribution. The results from the calculations show that it is possible to use CFD for calculation of the activity distribution in PNA. The mainly qualitative results in this work are encouraging and suggest further work. In the continuation of this work a response function for the gamma detector will be calculated so that a PNA time spectrum can be simulated. A more accurate comparison with experimental data can then be performed
This report provides test results on the thermal performance of an Owens-Illinois Sunpak liquid, evacuated tube, solar collector under simulated conditions. The test was conducted using the Marshall Space Flight Center Solar Simulator in accordance with the test requirements specified in ASHRAE 93-77 (Method of Testing to Determine the Thermal Performance of Solar Collectors) and the procedures contained in MTCP-FA-SHAC-400 (Procedure for Operation of the MSFC Solar Simulator Facility). The tests were performed on a module used on the early demonstration projects. A current production module is undergoing tests with results to be in a subsequent report.
This revised edition discusses a large number of renewable energy systems. Future trends are outlined. Economic efficiency calculations are explained and analyzed. Many exemplary calculations and graphical representations illustrate the various technologies and methods of calculation. This standard book addresses university students, scientists and engineers in research and industry. The accompanying DVD is one of the most extensive compilations of demo programs and full versions for simulation of renewable energy systems, all illustrations of the book in electronic format, additional text and a user-friendly navigation aid. (orig.)
The present paper is concerned with the prediction method for the motions and tensions of mooring lines on slack-moored floating oceanic architectural building in coastal zone. The motions and the tensions of mooring lines were estimated by the time domain numerical simulation in the present method. The present method includes the nonlinearity of mooring system, slow varying wave drift force. The authors performed a time domain numerical simulation for motions and the tensions of mooring lines on floating barge in two component waves. Also, the results obtained by numerical simulation were compared with the experimental results. The validity of the present method was confirmed.
The validation of the forecast model for early emergency response to nuclear accidents is evaluated by trace tests in atmosphere in Daya Bay nuclear power site. The simulation experiment of the Daya Bay nuclear power site shows that the particle spreading image and the time-integrated concentration distribution given by plume concentration prediction model can perform the variation of pathway of the pollutant transport, as well as the effects of topography on transport and diffusion of pollutants. The simulation of five trace tests in field shows that 59.1% of ratios between predicted results and observed results are within the range of 10, and 41% of ratios are within the range of 5 approximately. (authors)
Three-dimensional magnetohydrodynamic simulations are executed in a full toroidal geometry to clarify the physical mechanisms of the Internal Reconnection Event (IRE), which is observed in the spherical tokamak experiments. The simulation results reproduce several main properties of IRE. Comparison between the numerical results and experimental observation indicates fairly good agreements regarding nonlinear behavior, such as appearance of localized helical distortion, appearance of characteristic conical shape in the pressure profile during thermal quench, and subsequent appearance of the m=2/n=1 type helical distortion of the torus. (author)
A new scheme to accurately retain kinetic electron effects in particle-in-cell (PIC) simulations for the case of electrostatic drift waves is presented. The splitting scheme, which is based on exact separation between adiabatic and on adiabatic electron responses, is shown to yield more accurate linear growth rates than the standard df scheme. The linear and nonlinear elliptic problems that arise in the splitting scheme are solved using a multi-grid solver. The multi-grid particle-in-cell approach offers an attractive path, both from the physics and numerical points of view, to simulate kinetic electron dynamics in global toroidal plasmas.
The solvation properties of a chlorine ion in small water clusters are investigated using state-of-the-art statistical mechanics. The simulations employ the polarizable water model developed recently by Dang [J. Chem. Phys. [bold 97], 2659 (1992)]. The ion--water interaction potentials are defined such that the successive binding energies for the ionic clusters, and the solvation enthalpy, bulk vertical binding energy, and structural properties of the aqueous solution agree with the best available results obtained from experiments. Simulated vertical electron binding energies of the ionic clusters Cl[sup [minus
An accelerator experiment was performed using a low-energy antiproton beam to measure antiproton detection efficiency of BESS, a balloon-borne spectrometer with a superconducting solenoid. Measured efficiencies showed good agreement with calculated ones derived from the BESS Monte Carlo simulation based on GEANT/GHEISHA. With detailed verification of the BESS simulation, the relative systematic error of detection efficiency derived from the BESS simulation has been determined to be {+-}5%, compared with the previous estimation of {+-}15% which was the dominant uncertainty for measurements of cosmic-ray antiproton flux.
In the framework of the carbon storage two forestry methods are analyzed by the author: the standing capitalization (slowing of the crops) and the spices substitution (slow growing species by fast growing species). The construction of a model taking into account the all stage of the carbon implication (biomass, soil, forest products) offers simulation and quantitative results on these methods. (A.L.B.)
The paper is concerned with the simulation of practical testing situations which are too difficult and/or expensive to replicate in a laboratory environment. Numerical experiments are described which simulate the differential eddy current probe response to the build-up and chemical flushing of magnetite in the crevice gap of a PWR steam generator unit. The simulation results agree well with the only experimental data available to the authors and lead to the conclusion that conventional differential eddy current probes should be capable of characterizing crevice gap conditions with respect to the presence of magnetite. (author).
The hot uniform aluminium plasma was produced by irradiating thin aluminium dotted foil smoothly with the 9th 0.53 ?m laser on Shenguang II laser facility. The emitted spectrum was measured from the front and tangential direction of the target with two crystal spectrometers, and the quantitative spectrum from the front of the target was obtained. The state of laser- produced plasma was simulated with the radiation hydrodynamics code MULTI-1D, and the emitted spectrum was calculated with the spectrum code of Collision-Radiation model under the simulated plasma state. The experimental spectrum accords with the simulated one. (authors)
Thermal expansions of a simulated fuel (SS-1) were measured by using a Dilatometer (DIL402C) from room temperature to 1900 K. The main procedure of an uncertainty evaluation followed the strategy of the UO{sub 2} fuel. Referring to the ISO (International Organization for Standardization) guide, the uncertainties of the thermal expansion were quantified in three parts - the initial length, the length variation, and the system calibration factor. The uncertainty of the thermal expansion for a simulated fuel was also compared with those of UO{sub 2} fuel.
The results of this study show that Monte Carlo simulations of diffusion in homogeneous turbulence can be formulated in terms of the Eulerian space-time velocity autocorrelation function. Numerical results obtained using this approach agree with results obtained by Taylor (1921) using the Lagrangian autocorrelation function. Use was made of the equivalence of the Lagrangian and Eulerian Monte Carlo approaches to derive analytic relations between the Lagrangian integral time scale and the Eulerian integral space and time scales. These analytic results were verified by comparison to Monte Carlo simulations and to other theoretical results. They are in general agreement with many existing theories and semi-empirical relations.
It has been issued that ROP(Regional Overpower Protection) for Wolsong Unit 1 needed to be reanalyzed due to the aging effect. Thermo-hydraulics and core simulation have to be performed for calculation of the fuel bundle power, channel power and detector signal production. PPV/MULTICELL/RFSP code system was used to calculate the power distribution for the ROP analysis. In this study, 232 cases out of 926 scenarios which include postulated accidents such as Startup after Short Shutdown, Shim Cases, Stepback, Insertion and Withdrawal of Reactivity Control Rods were simulated.
The simulation of processes in engineering and the physical sciences has progressed rapidly over the last several years. With rapid developments in supercomputers, parallel processing, numerical algorithms and software, scientists and engineers are now positioned to quantitatively simulate systems requiring many billions of arithmetic operations. The need to understand and assimilate such massive amounts of data has been a driving force in the development of both hardware and software to create visual representations of the underling physical systems. In this paper, and the accompanying videotape, the evolution and development of the visualization process in scientific computing will be reviewed. Specific applications and associated imaging hardware and software technology illustrate both the computational needs and the evolving trends. 6 refs.
An accelerator beam experiment was performed using a low-energy antiproton beam to measure antiproton detection efficiency of the BESS detector. Measured and calculated efficiencies derived from the BESS Monte Carlo simulation based on GRANT/GHEISHA showed good agreement. With detailed verification of the BESS simulation, the relative systematic error of detection efficiency derived from the BESS simulation has been determined to be {+-}5%, compared with the previous estimation of {+-}15% which was the dominant uncertainty for measurements of cosmic-ray antiproton flux. (author)
Electrostatic plasma simulation code for 2D rectangular geometry is presented. Main distinguishing feature of the code is its orientation on the beam-plasma interaction. The code and its graphical interface were developed using MATLAB programming language. Simulation results of inhomogeneous plasma interaction with modulated electron beams of different width are compared. In case of wide beam the front of Langmuir waves generated in point of local plasma resonance is planar and in case of thin beam (or ribbon beam) the front has approximately half-circular form.
Engine flow simulation was provided by four separately mounted air ejectors connected to a high-pressure air supply. The engine nacelle center lines were ...
A stochastic wind simulation for VAWTs (VSTOC) has been developed which yields turbulent wind-velocity fluctuations for rotationally sampled points. This allows three-component wind-velocity fluctuations to be simulated at specified nodal points on the wind-turbine rotor. A first-order convection scheme is used which accounts for the decrease in streamwise velocity as the flow passes through the wind-turbine rotor. The VSTOC simulation is independent of the particular analytical technique used to predict the aerodynamic and performance characteristics of the turbine. The VSTOC subroutine may be used simply as a subroutine in a particular VAWT prediction code or it may be used as a subroutine in an independent processor. The independent processor is used to interact with a version of the VAWT prediction code which is segmented into deterministic and stochastic modules. Using VSTOC in this fashion is very efficient with ...
Description of the current status of the developments of a simulation of the Darlington Nuclear Generating Station shutdown and regulating systems, DARSIM. The DARSIM program simulates the spatial neutron dynamics, the regulations of the reactor power, and shutdown system 1 and shutdown system 2 software. The DARSIM program operates in the interactive simulation (INSIM) program environment. DARSIM was installed on the APOLLO computer at the Atomic Energy Control Board (AECB) and a version for an IBM-PC was also provided for the exclusive use of the AECB. Shutdown system software was updated to incorporate the latest revisions in the functional specifications. Additional developments were provided to assist in the use and interpretation of the DARSIM results.
... UGV), in the air (UAV), or underwater (UUV) -across ... Figure 2 Tactical Mobile Robot MOUT concept: ... 10 robots to clear 20 rooms in 60 minutes ...
A catalytic converter of a ceramic monolith honeycomb substrate, coated with a washcoat of catalyst and attached to a natural gas/diesel dual fuel engine was simulated and studied experimentally. The paper describes the application of one-dimensional finite element model for the transient and steady state operation. Laminar flow was approximated using a dispersed plug flow model, and chemical kinetics were simulated using LHHW (Langmuir/ Hinshelwood/ Hougan/ Watson) type expressions. Simulation results were compared with experimental results for heating and cooling cycles which resulted from speed and load changes on the engine. The comparison showed a maximum difference between the two sets of emission levels of about 10 per cent, showing that the one-dimensional model is acceptable model for this dual fuel engine converter combination. 50 refs., 3 tabs., 13 figs.
The antimicrobial activity of simulated solar disinfection (SODIS) in the presence and absence of riboflavin against various protozoa and helminth organisms was investigated in this study. Assays were conducted in transparent 12 well microtitre plates containing a suspension of test organisms in the presence or absence of 250mM riboflavin. Plates were exposed to simulated sunlight at an optical irradiance of 550Wm^-^2 (watts per square metre) delivered from a SUNTEST(TM) CPS+ solar simulator. Aliquots of the test suspensions were taken at set time points and the viability of the test organisms was determined by either culture, microscopy or flow cytometry where applicable. With Acanthamoeba, Naegleria, Entamoeba and Giardia exposure to SODIS at an optical irradiance of 550Wm^-^2 for up to ...
Inductive electromagnetic launchers, or coilguns, use discrete solenoidal coils to accelerate a coaxial conductive armature. To date, Sandia has been using an internally developed code, SLINGSHOT, as a point-mass lumped circuit element simulation tool for modeling coilgun behavior for design and verification purposes. This code has shortcomings in terms of accurately modeling gun performance under stressful electromagnetic propulsion environments. To correct for these limitations, it was decided to attempt to closely couple two Sandia simulation codes, Xyce and ALEGRA, to develop a more rigorous simulation capability for demanding launch applications. This report summarizes the modifications made to each respective code and the path forward to completing interfacing between them.
Abstract Suspension of solids in stirred reactor is widely used for catalytic reactions, dissolution, etc. Quality of solid suspension is an important parameter required for the reliable design, optimum performance, and scale up of the system. Quality of suspension depends on local characteristics of solid velocity and hold up profiles. The present work was focused on investigating quality of solid suspension using ultrasound velocity profiler (UVP) measurements and CFD simulations. The slip velocity measurements carried out with UVP were used to evaluate different drag correlations used in CFD simulations. Results discussed in this work would be useful for extending the applications of CFD models for simulating large stirred slurry reactors.
This paper describes the simulation of the operating environment of the New York Power Pool (NYPP) within the Dispatcher Training Simulator facility at the New York State Electric and Gas Corporation (NYSEG). The NYSEG control center uses real-time data received from the NYPP control center to control its internal generation. In addition, the NYSEG dispatchers work very closely with the NYPP dispatchers. The exact emulation of the Generation Control functions of the NYPP EMS, and the associated data exchange with the NYSEG ECS, contribute greatly to the high degree of realism achieved in the NYSEG DTS facility. The accuracy of the power system simulation in the DTS has been validated against actual power system disturbances recorded at NYSEG.
In this paper the simulation of hybrid processes containing membrane units will be discussed. For this purpose a user-defined module for simulation and design of membrane processes was implemented into the simulation program AspenPlus. The advantages can be summarized as follows: - any combination of membrane processes with all other units already implemented in AspenPlus is possible, including internal recycle streams, - utilization of the physical property models and data bases of AspenPlus is possible, - cost and sensitivity analysis can be performed. These benefits are demonstrated in detail for a membrane vapor recovery unit for the treatment of tank farm off-gas, for a two-stage reverse osmosis plant for organic/-organic separations and for a combination of distillation and pervaporation for the separation of a dimethylcarbonate/methanol mixture. (orig.)
In this book are found technical and scientific papers on the main works of the Direction of the Fuel Cycle (DCC) in France. The study fields are: the up-side of the nuclear fuel cycle with theoretical studies (plasma simulation) and technological developments and instrumentation (lasers diodes, carbides plasma projection, carbon 13 enrichment); the down-side nuclear fuel cycle with theoretical studies (ion Eu{sup 3+} complexation simulation, decay simulation, uranium and plutonium diffusion study, electrolyser operating simulation), scenario studies ( recycling, wastes management), experimental studies; dismantling and cleaning (soils cleaning, surface-active agent for decontamination, fault tree analysis); analysis with expert systems and mass spectrometry. (A.L.B.)
This guide describes the simulator`s governing equations, constitutive functions and numerical solution algorithms of the STOMP (Subsurface Transport Over Multiple Phases) simulator, a scientific tool for analyzing multiple phase subsurface flow and transport. The STOMP simulator`s fundamental purpose is to produce numerical predictions of thermal and hydrologic flow and transport phenomena in variably saturated subsurface environments, which are contaminated with volatile or nonvolatile organic compounds. Auxiliary applications include numerical predictions of solute transport processes including radioactive chain decay processes. In writing these guides for the STOMP simulator, the authors have assumed that the reader comprehends concepts and theories associated with multiple-phase hydrology, heat transfer, thermodynamics, radioactive chain decay, and nonhysteretic relative permeability, ...
The impact of the azimuthal wave refraction in the middle atmosphere on the distribution of gravity wave amplitudes, propagation azimuths, and other wave parameters is investigated using a numerical ray-tracing model of gravity wave propagation through a representative zonal mean reference model of geostrophic winds and temperature in the middle atmosphere. Simulations are first performed with only a single type of gravity wave to help explain some important effects which occur during the refraction process. Then, a multiray simulation is performed which traces a crude spectrum of waves from different altitudes through the atmosphere for every month of a climatological year. The simulated wave climatologies are compared with observations. 108 refs.
Abstract The electrochemical properties of a perfluorosulfonic acid (PFSA) membrane are estimated using a combination of molecular dynamics simulation and statistical thermodynamic model. We obtain all parameters in an ionic conductivity model from an atomistic simulation and remove all adjusted model parameters. From a microscopic point of view, the hydrated PFSA membrane shows micro-phase segregation which separated into hydrophilic and hydrophobic phases. Our present work originates with this phenomenon and we treat this phase segregation as if it is a continuous phase for each of which the proton (H+) is transported inside the PFSA membrane/solvent (water and alcohols) mixture. The chemical potential for a given system is estimated using a molecular simulation technique to predict the ...
terms of 11) numerical models, 2D SPH simulations, and analytic models. The 1D model was useful for determining the evolution over very long timescales. ...
... one for each node, and an i/o queue to manage ... ecuted as if it came from a terminal attached to a ... In its default mode PC-NETSIM is stochastic, that is ...
RODbaffle heat exchanger is a kind of shell-and-tube heat exchanger with eminent performance. Because of the characteristics of tube bundle support structure, fluid flow on shellside is longitudinal and periodical, and fluid flow is of symmetry in landscape orientation. According to the fluid flow characteristics on shellside and ignoring the impact of inlet, outlet and shell wall, a periodic flow unit duct was taken as the simplified model of the shellside to perform numerical simulation by using CFD code, FLUENT. It was found that the both errors in magnitude of the main fluid velocities and heat transfer coefficients between results of simulation and that of experiment or correlations are less than 10%, and the errors of pressure drop between simulation and correlation are no more than 20%. The detailed characteristics and relation of fluid flow and heat transfer on shellside of the RODbaffle heat exchanger were analyzed ...
... that range from relatively low frequency strong ground motion modeling [Vidale et al., 1985] to high frequency reverse time migration of seismic ...
... being pulled out of water, the movement is restricted ... This heating depended on the particle size ... hydrometeor injection into the upper tropopause. ...
Data on French nuclear electricity generation sites, nuclear power plant operations personnel, operation simulators, nuclear training centers and training statistics are presented.
A non-dissipative drift kinetic simulation scheme, which rigorously satisfies the time-reversibility, is applied to the three-mode coupling problem of the ion temperature gradient (ITG) instability. It is found from the simulation that the three-mode ITG system repeats growth and decay with a period which shows a logarithmic divergence for infinitesimal initial perturbations. Accordingly, time average of the mode amplitude vanishes, as the initial amplitude approaches to zero. An exact solution is analytically given for a class of initial conditions. An excellent agreement is confirmed between the analytical solution and numerical results. The results obtained here provide a useful reference for basic benchmarking of theories and simulation of the ITG modes. (author)
To evaluate soil erodibility under different land uses and to study the applicability of nomograph for estimation of soil erodibility a field experiment was conducted under both natural and simulated rainfall conditions under four land uses viz. barren, cultivated, grassland, and forest in the sub-mountainous tract of Punjab (India). Measured soil erodibility (K) values varied from 033 to 067 under natural rainfall conditions and from 023 to 040 under simulated rainfall conditions. Among different land uses, measured K was in the order of barren > cultivated > grassland > forest soils. The values of the K estimated by nomograph were very low as compared to the observed values. The trends were also in contrast to these observed values of K under simulated and natural rainfall conditions. To...
Due to size and costs of big slewing rings, the common methods of experimental trial and error together with destructive test methods for verification of the inductive hardening process are limited. Thereby it is essential to simulate the hardening process in addition to the tests. With these simulations, it should be possible to get conclusions for the hardening pattern, the residual stress and distortion. With actually existing methods a fully three dimensional simulation is extremely time consuming. For that reason, a new method was developed that combines two dimensional with three dimensional simulation to shorten the time for calculation tremendously. (orig.)
Stochastic models for the simulation of global radiation are discussed. Thermal transients in the ground are analyzed. The performance of buried-pipe storage and a space heating system with long-term storage is described.
The utilization of electrically heated rods for the simulation of nuclear fuel pins represents a generally adopted method by the nuclear industry to study thermalhydraulic problems. Usually its is necessary to determine the time variation of the electric linear power to simulate a given nuclear power transient in order to yield the same temperature and heat flux conditions in the surface of the electrical heater that would be observed in the nuclear fuel pin. The present work analyzes the limits of the usually adopted simulation methods and shows a manner to obtain the required electrical linear power that reduces oscillations and yields accurate results for the thermal conditions of the rod surface wall. (author). 5 refs, 5 figs, 1 tab.
The utilization of electrically heated rods for the simulation of a nuclear fuel pins represents a generally adopted method by the nuclear industry to study thermalhydraulic problems. Usually, it is necessary to determine the time variation of the electric linear power to simulate a given nuclear power transient in order to yield the same temperature and heat flux conditions in the surface of the electrical heater that would be observed in the nuclear fuel pin. The present work analyses the limits of the usually adopted simulation methods and shows a manner to obtain the required electrical linear power that reduces oscillations and yields accurate results for the thermal conditions of the rod surface wall. (author) 5 refs., 5 figs., 1 tab.
... AGARD AR-319, Volume 2. Knight, D., Zhou ... a Turbulent Boundary Layer in a Supersonic Flow. ... of Development of Separated Flows in Compression ...
... turbine blade heat-transfer and study the heat-transfer augmentation in idealized geometry and its ... responsible for heat transfer augmentation. ...
Aug 18, 2004 ... The Arabian Gulf region presents a challenge to meteorologists trying to simulate weather with computer models because sea-surface and land ...
The results obtained from experimental investigations and mathematical simulation of horizontal steam generators are considered. Recommendations for continuing these works are given.
Mar 15, 2011 ... Information Technology. NCCS Discover Supercomputer. Science mission data collection, analysis and simulations pose unprecedented computer ...
We have simulated the transient enhanced diffusion (TED) of boron fluoride (BF_2"+) implanted in crystalline and germanium amorphized silicon. Based on recently published models, the effect of fluorine on boron diffusion in silicon has been introduced and a modelling has been suggested. In order to simulate the boron experimental profiles, we have assumed that fluorine forms clusters involving interstitial boron which reduces the junction depth. Experimental results indicate that fluorine behaviour depends on amorphization energy. Moreover, even no germanium preamorphization is performed, silicon is still amorphized by fluorine species. Hence, BF_2"+ implantation leads to an amorphous/crystalline (a/c) interface near the surface. An improvement of published models is suggested taking into account fluorine effects. The simulations satisfactory reproduce the SIMS experimental profiles for a large scale of experimental ...
Jan 12, 2010 ... Visual simulation of offshore liquefied natural gas (lng) terminals in a decision-making context. Baird, Brian E.; Sheppard, Stephen R. J.; ...
In this paper we present a model and simulator for many clusters of heterogeneous PCs belonging to a local network. These clusters are assumed to be connected to each other through a global network and each cluster is managed via a local scheduler which is shared by many users. We validate our simulator by comparing the experimental and analytical results of a M/M/4 queuing system. These studies indicate that the simulator is consistent. Next, we do the comparison with a real batch system and we obtain an average error of 10.5\\% for the response time and 12\\% for the makespan. We conclude that the simulator is realistic and well describes the behaviour of a large-scale system. Thus we can study the scheduling of our system called \\dirac in a high throughput context. We justify our decentralized, adaptive and opportunistic approach in comparison to a centralized approach in such a context.
An impingement jet system was used to study flow-assisted corrosion (FAC) of 3003 aluminum (Al) alloy in ethylene glycol-water solutions that simulates the automotive coolant by corrosion potential and electrochemical impedance spectroscopy (EIS) measurements as well as computational fluid dynamics (CFD) simulation. The effects of solution pH and fluid impact angle on Al FAC were determined. An increase of solution pH enhances the activity of Al due to dissolution of Al oxide film in alkaline environment. Moreover, Al activity decreases with the increasing fluid impact angle to the specimen. A CFD simulation shows that, with the increase of impact angle, the electrode area under high-velocity flow field decreases and that under low-velocity flow field increases. Consequently, the shear str...
Dynamic clamp is a powerful method that allows the introduction of artificial electrical components into target cells to simulate ionic conductances and synaptic inputs. This method is based...Full Text Available
A complete analysis is performed to simulate biomass production incorporated into a realistic whole farm situation, including or replacing a typical crop mix. Representative farms are constructed to accommodate such simulation. Four management systems are simulated for each firm, with each simulation depicting a different crop mix and/or use of different farming technologies and production methods. The first simulation was a base farm plan in which the operator would maintain the historical crop mix for the area, participate in all price support programs, and not participate in either a conservative reserve or a biomass production program. In the second simulation, the operator would again maintain the historical crop mix, would not participate in a conservation reserve or biomass production program, and would be ineligible to participate in any price support ...
A complete analysis is performed to simulate biomass production incorporated into a realistic whole farm situation, including or replacing a typical crop mix. Representative farms are constructed to accommodate such simulation. Four management systems are simulated for each firm, with each simulation depicting a different crop mix and/or use of different farming technologies and production methods. The first simulation was a base farm plan in which the operator would maintain the historical crop mix for the area, participate in all price support programs, and not participate in either a conservative reserve or a biomass production program. In the second simulation, the operator would again maintain the historical crop mix, would not participate in a conservation reserve or biomass production program, and would be ineligible to participate in any price support ...
The aims of the present work were to investigate possible compatibility between sodium zirconium phosphate (NZP) and Synroc titanate phases, to prepare NZP-based waste forms by hot-pressing rather than sintering, and to investigate the incorporation in NZP of (a) Cs/Sr as simulated heat-generating nuclides; (b) simulated actinides; and (c) simulated Purex waste. The NZP samples were prepared by methods similar to those used for Synroc. The precursor NZP phase was formed from tetrabutyl zirconate Zr(OC{sub 4}H{sub 9}){sub 4}, sodium nitrate, and 85% orthophosphoric acid. Simulated waste nitrate solutions were then mixed with the liquid precursor. After stir drying of the precursor, calcination was carried out at 700{degree}C to remove nitrates and organics.
7.5.6 Conclusions of the Generic Metroplex Queueing Simulation . ...... the systemwide arrival delays incurred at the metroplex terminal-area ...... assignment rules, estimated nominal transit speeds or times by domain, and stochastic ...
We present results of multiple-time-scale simulations of 5, 10 and 15 keV low temperature ion implantation of arsenic on silicon (100), followed by high temperature anneals. The simulations start with a molecular dynamics (MD) calculation of the primary state of damage after 10ps. The results are then coupled to a kinetic Monte Carlo (MC) simulation of bulk defect diffusion and clustering. Dose accumulation is achieved considering that at low temperatures the damage produced in the lattice is stable. After the desired dose is accumulated, the system is annealed at 800 degrees C for several seconds. The results provide information on the evolution for the damage microstructure over macroscopic length and time scales and affords direct comparison to experimental results. We discuss the database of inputs to the MC model and how it affects the diffusion process.
Orientational constraints obtained from solid state NMR experiments on anisotropic samples are used here in molecular dynamics (MD) simulations for determining the structure and dynamics of several different membrane-bound molecules. The new MD technique is based on the inclusion of orientation dependent pseudo-forces in the COSMOS-NMR force field. These forces drive molecular rotations and re-orientations in the simulation, such that the motional time-averages of the tensorial NMR properties approach the experimentally measured parameters. The orientational-constraint-driven MD simulations are universally applicable to all NMR interaction tensors, such as chemical shifts, dipolar couplings and quadrupolar interactions. The strategy does not depend on the initial choice of coordinates, and...
... which may have embedded remote terminals, and (4 ... model the time-variant or stochastic aspects of ... on such system parameters as queue size, time ...
This paper presents a new approach to improve the in silico modeling of ultrasound transducer arrays. While current simulation tools accurately predict the theoretical element...Full Text Available
Engine model: An engine and governor model adapted from a heavy lift helicopter simulation provides a restistic time delay between aerodynamic rotor ...
... 1 1 if on-line printout is desired (data will also be stored on tape) = 0 or blank, reduced data will be stored on a listable tape P3 M Codes B - Corivert ...
joint task force (JTF), division (OIV) and regimental combat team (RCT) with ..... Headquarters, Department of the Army. (HQOA), (2010). Field Manual 2-0 ...
The Stockpile Stewardship Program (SSP) is a single, highly integrated technical program for maintaining the surety and reliability of the U.S. nuclear stockpile. The SSP uses past nuclear test data along with current and future non-nuclear test data, computational modeling and simulation, and experimental facilities to advance understanding of nuclear weapons. It includes stockpile surveillance, experimental research, development and engineering (D&E) programs, and an appropriately scaled production capability to support stockpile requirements. This integrated national program requires the continued use of current facilities and programs along with new experimental facilities and computational enhancements to support these programs. The Advanced Simulation and Computing Program (ASC) is a cornerstone of the SSP, providing simulation capabilities and computational resources to support the annual stockpile assessment ...
The Stockpile Stewardship Program (SSP) is a single, highly integrated technical program for maintaining the safety and reliability of the U.S. nuclear stockpile. The SSP uses past nuclear test data along with current and future nonnuclear test data, computational modeling and simulation, and experimental facilities to advance understanding of nuclear weapons. It includes stockpile surveillance, experimental research, development and engineering programs, and an appropriately scaled production capability to support stockpile requirements. This integrated national program requires the continued use of current facilities and programs along with new experimental facilities and computational enhancements to support these programs. The Advanced Simulation and Computing Program (ASC)1 is a cornerstone of the SSP, providing simulation capabilities and computational resources to support the annual stockpile assessment and ...
To further optimize thermo-hydraulic capability of HCSB-TBM, two different thermo-hydraulic schemes would be proposed and investigated by relevant calculating and simulating of thermo-dynamics in normal and peak heat flux condition. With computational fluid dynamics (CFD) code FLUENT, the actual temperature distribution of key components of HCSB-TBM is simulated. (authors)
The rare earth antimonates RE_3Sb_5O_1_2 constitute an isostructural series, where the rare earth site symmetry is S_4. The spectroscopic properties of pure or doped compounds of this family are studied in order to derive their energy level schemes and to simulate them through the crystal field theory, by using the approximated D_2_d potential. The paramagnetic susceptibilities are calculated with the derived wavefunctions and compared to experiment.
We demonstrate a two-dimensional device simulator for MOSFET structures that incorporates models for defect distributions and show predicted effects on device switching performance for various spatial distributions of defects in amorphous and polycrystalline silicon.
This report describes the test of a dynamic phase compensation unit for a wind turbine with directly connected induction generators. The compensation unit is based on thyristor switched capacitors, where conventional wind turbine compensations use mechanical contactors to switch the capacitors. The influence on power quality analysed, and influence on component lifetime is discussed. Besides, simulation models in Matlab/Simulink are presented, including a flicker meter model. (au)
Symplectic integrators designed for simulating soft matter at constant temperature and constant pressure (or constant surface tension) are presented. In addition to the well-known merits of symplectic integrators, such as long time stability, these methods allow the estimation of thermal heat production/absorption and thus entropy differences through a simple function of the time scaling factor ?. The relative entropy estimation around crystal-liquid phase transition is given for systems with soft core repulsive potential.
Gas and stars in spiral galaxies are modelled with the DUAL code, using hydrodynamic and N-body techniques. The simulations reveal morphological differences mirroring the dual morphologies seen in B and K' band observations of many spiral galaxies. In particular, the gaseous images are more flocculent with lower pitch angles than the stellar images, and the stellar arm-interarm contrast correlates with the degree of morphological decoupling.
The peculiarities of pulsed laser melting and evaporation of the superconducting ceramics are analyzed by means of numerical simulation. The appearance of the overheated metastable states in solid and liquid phases is shown as a result of the phase front dynamics and volume nature of laser energy release. A method of dynamic adaptation for the multifront Stefan problem is proposed.
The model of transient enhanced diffusion of ion-implanted As is formulated and the finite-difference method for numerical solution of the system of equations obtained is developed. The nonuniform distribution of point defects near the interface and more accurate description of arsenic clustering are simultaneously taken into account. Simulation of As diffusion during rapid annealing gives a reasonable agreement with the experimental data. (authors)
Computer simulations and experiments indicate that we can shape the infrared picosecond pulses of the Mark III FEL in amplitude, frequency, and phase. Strongly modulated fundamental and second harmonic pulses have been generated by operating the Mark III FEL in the regime of strong sideband growth. In this paper, we present the results of simulations and experiments for second harmonic generation with fundamental inputs from 2 to 3 {mu}m.
The paper deals with a model in which the motion of a gas in plasma accerators and high-curent discharges in the present of a skin effect is treated as expulsion by a 'magnetic piston' under the action of surface current flow. Specifically examined is the situation where the initial gas pressure is negligible in comparison with the magnetic pressure, and the motion of the gas may be treated as a self-simulating one. A system of hydrodynamic equations of the problem is derived and integrated by Adam's method. Results are plotted and discussed.
A reliable inventory control system was developed at the Fast Flux Test Facility (FFTF) to keep track of the occupancy of 900 refueling facility locations, to compile historical data on the movement of each reactor assembly, and to simulate assembly moves. The simulate capability is valuable because it allows verification of documents before they are issued for use in the plant, and eliminates the possibility of planning illegal or impossible moves. The system is installed on a UNIVAC 1100 computer and is maintained using a data base management system by Sperry Univac called MAPPER.
There is need of pressure surge analyses when valves or pumps are activated or piping systems fail (pipe rupture). Based on actual problems the influences of boundary conditions upon fluid simulation results are discussed. Hints concerning realistic dynamic analyses of piping systems are presented. Some of the simulations results are compared with measurements. (orig.)
This is an inquiry-based module created specifically for use with the PhET simulation Geometric Optics. It includes a lesson plan for teaching about convex lenses, step-by-step student directions, and a set of Power Point "clicker" questions for use as warm-up questions or informal assessment. This resource is part of a larger collection of interactive Java simulations for students of physics, developed by the Physics Education Technology project at the University of Colorado.
Performance predictions have been made for attached-sunspace types of passively solar heated buildings. The predictions are based on hour-by-hour computer simulations using computer models developed in the framework of PASOLE, the Los Alamos Scientific Laboratory (LASL) passive solar energy simulation program. The models have been validated by detailed comparison with actual hourly temperature measurements taken in attached-sunspace test rooms at LASL.
A physical model of free-electron laser (FEL) amplifier with variable-parameter wiggler magnets for one-dimensional numerical simulation is presented and a numerical example is given. The wiggler parameters, efficiency of energy conversion between electron beam and laser field, laser intensity, phase-space distributions and energy spectrum of electrons are computed. The period of synchronous oscillation and saturation value of laser intensity agree with estimated one.
The Monte Carlo simulation method is applied to calculate the nucleon transfer part of the imaginary optical-model potential in colliding heavy ions. In the process of transfer of nucleons from one nucleus to the other, forward and backward neutrons and protons are treated separately by taking into consideration the effect of the Coulomb forces and the potential energy surfaces (driving forces). The Pauli exclusion principle is also incorporated properly in the calculation.
A Siemens Mevatron KV2 accelerator installed at the Royal Adelaide Hospital employs cylindrical solid-walled electron cones for some electron collimation. The cones being used at present result in treatment fields that do not always conform with the International Electrotechnical Commission (IEC) Standards (particularly at high energies). The aim of this project was to simulate the existing cones using Monte Carlo methods in order to evaluate potential cone modifications required to overcome the field irregularities. Simulations were performed using the EGS4 (Electron Gamma Shower version 4, distribution II) Monte Carlo code installed on a DEC Alpha workstation at the University of South Australia. To rigorously simulate the existing electron cones it was necessary to also simulate various components within the treatment head of the linear accelerator. Results of simulations for ...
Molecular simulation aims at simulating particles in interaction, describing a physico-chemical system. When considering Markov Chain Monte Carlo sampling in this context, we often meet the same problem of statistical efficiency as with Molecular Dynamics for the simulation of complex molecules (polymers for example). The search for a correct sampling of the space of possible configurations with respect to the Boltzmann-Gibbs distribution is directly related to the statistical efficiency of such algorithms (i.e. the ability of rapidly providing uncorrelated states covering all the configuration space). We investigated how to improve this efficiency with the help of Artificial Evolution (AE). AE algorithms form a class of stochastic optimization algorithms inspired by Darwinian evolution. Efficiency measures that can be turned into efficiency criteria have been first searched before identifying parameters that could be ...
Molecular dynamics computer simulations of a dense nanocrystalline Nickel sample are performed on the parallel Cray T3 (EPFL-Lausanne). The sample contains 50 grains with an average size of 3 nm. By perturbing the relaxed configuration elastic and plastic properties are studied. (author) 1 fig., 1 ref.
The point-defect-impurity pair diffusion model proposed recently by Mulvaney and Richardson is adopted and modified to simulate the coupled diffusion of phosphorus and self-interstitials in phosphorus-implanted silicon. The assumption of implantation-induced, but empirically determined initial interstitial distributions of Gaussian shape allows a simulation of the net effect of transient enhanced diffusion. As a result an improved modeling of phosphorus diffusion in silicon is achieved for a broad range of ion-implantation and annealing conditions. (author).
We study the dynamical evolution of globular clusters containing primordial binaries, including full single and binary stellar evolution using our Monte Carlo cluster evolution code updated with an adaptation of the single and binary stellar evolution codes SSE and BSE from Hurley et al. We describe the modifications that we have made to the code. We present several test calculations and comparisons with existing studies to illustrate the validity of the code. We show that our code finds very good agreement with direct N-body simulations including primordial binaries and stellar evolution. We find significant differences in the evolution of the global properties of the simulated clusters using stellar evolution compared with simulations without any stellar evolution. In particular, we find that the mass loss from the stellar evolution acts as a significant energy production channel simply by reducing the total gravitational ...
A series of first principles molecular dynamics and Monte Carlo simulations were carried out for liquid water to assess the validity and reproducibility of different sampling approaches. These simulations include Car-Parrinello molecular dynamics simulations using the program CPMD with different values of the fictitious electron mass in the microcanonical and canonical ensembles, Born-Oppenheimer molecular dynamics using the programs CPMD and CP2K in the microcanonical ensemble, and Metropolis Monte Carlo using CP2K in the canonical ensemble. With the exception of one simulation for 128 water molecules, all other simulations were carried out for systems consisting of 64 molecules. It is found that the structural and thermodynamic properties of these simulations are in excellent agreement with each other as long as adiabatic sampling is maintained in the ...
In order to investigate the extrusion characteristics of magnesium alloy (AZ31), a computer simulation was attempted. Tensile properties of as-cast billet with different strain rates were incorporated into the simulation. The results showed a great change in distribution of stress and strain at near die region by the ram travel. It was found that the average grain size at longitudinal edge of the extruded sheet was smaller than that at center region, which might be attributed to recrystallization caused by severe plastic deformation. (orig.)
Two TNT high explosives simulants have been developed. Small scale testing has shown them to be insensitive to: impact, spark, friction, temperature, and shock. The materials have been scaled to 0.5 kg quantities and samples given to the Protective Services Department for field evaluation using explosives detecting canines.
The continuous bioprocessing of simulated steam-stripped sulfite pulp mill evaporator condensates using Candida utilis to produce single cell protein has been studied. At 32 degrees C, processing can be conducted at space velocities of up to 0.43 hr/sup -1/ with over 90% removal of acetic acid present at an initial concentration of 2 g/liter. (Refs. 16).
A new algorithm for the treatment of sliding interfaces between solids with or without friction in an Eulerian wavecode is described. The algorithm has been implemented in the two-dimensional version of the CTH code. The code was used to simulate penetration and perforation of aluminum plates by rigid, conical-nosed tungsten projectiles. Comparison with experimental data is provided.
Interaction of drilling fluids with a geothermal reservoir formation can result in significant permeability impairment and therefore reduced well productivity. This interaction is studied under simulated in situ geothermal conditions of overburden stress, pore fluid pressure, temperature, and pore fluid chemistry. Permeability impairment of an East Mesa KGRA reservoir material is evaluated as a function of stagnation time, drilling fluid, and temperature. Results indicate that all of these parameters contribute significantly to the magnitude and the reversibility of the impairment.
Simulations of young stellar systems suggest that planets embedded in a circumstellar disk can produce many distinctive structures, including rings, gaps and spiral arms. This video compares computer simulations of hypothetical systems to an image of system SAO 206462 taken by the Subaru Telescope and its HiCIAO instrument. n nCredit: NASA Goddard Space Flight Center/NCSA
Power devices such as MOSFETSs and IGBTs, include parasitic structures that can give rise to destructive failures such as breakdown and latch-up. To determine a suitable strategy for device radiation hardening, simulation software like MEDICI-2D can be used to model the effects of technological modifications and device parameters that are difficult to measure experimentally. (authors).
The corrosion rate of low alloy steel SA-508 and carbon steel A-410b in simulated operation and shutdown conditions of pressurized water reactor has been determined Moreover potentiodynamic polarization curves and galvanic effect through coupling of AISI-304 have been carried out under shutdown simulated condition. (Author) 8 refs.
Taking both white noise and colored environmental noise into account, a general stochastic logistic model under regime switching is proposed and studied. Sufficient conditions for extinction, nonpersistence in the mean, weak persistence, stochastic permanence and global attractivity are established. The critical number between weak persistence and extinction is obtained. Moreover, some simulation figures are introduced to illustrate the main results.
This work deals with a time-resolved optical study of the avalanche and streamer formation phases leading to breakdown in atmospheric nitrogen. The authors present the results obtained for nitrogen, from experiments and two-dimensional model simulations. This model is used to obtain a better insight in the relevant mechanisms and processes by a comparison of measurements and simulation data. The trends of externally measured quantities correspond with those predicted by the model.
A simple and continuous model for the on-state current of polysilicon thin-film transistors, suitable for implementation in circuit simulators, is presented. The model includes the potential barrier at the grain boundaries, the channel length modulation and the excess current due to impact ionization. Comparison between measured output characteristics and the model shows excellent agreement over wide range of bias voltages and for devices with different gate lengths.
In this paper, finite-difference time-domain method is used to research the transient electromagnetic pulse (EMP) coupling to the cylinder object which is shielded in another shielding cavity with a hole numerically. The different coupling effectiveness of three transient EMP is simulated. The distribution of the electric field in the shielding cavity and the current on the cavity are researched. The shielding effectiveness is evaluated. (authors)
This document is a reference guide to the Xyce Parallel Electronic Simulator, and is a companion document to the Xyce Users Guide. The focus of this document is (to the extent possible) exhaustively list device parameters, solver options, parser options, and other usage details of Xyce. This document is not intended to be a tutorial. Users who are new to circuit simulation are better served by the Xyce Users Guide. The Xyce Parallel Electronic Simulator has been written to support, in a rigorous manner, the simulation needs of the Sandia National Laboratories electrical designers. It is targeted specifically to run on large-scale parallel computing platforms but also runs well on a variety of architectures including single processor workstations. It also aims to support a variety of devices and models specific to Sandia needs. This document is intended to complement the Xyce Users Guide. It contains ...
The He-shell flash convection in AGB stars is the site for the high-temperature component of the s-process in low- and intermediate mass giants, driven by the Ne22 neutron source. [...] The upper convection boundary plays a critical role during the H-ingestion episode that may lead to neutron-bursts in the most metal-poor AGB stars. We address these problems through global 3-dimensional hydrodynamic simulations including the entire spherical He-shell flash convection zone (as oposed to the 3D box-in-a-star simulations). An important aspect of our current effort is to establish the feasibility of our appoach. We explain why we favour the explicit treatment over the anelastic approximation for this problem. The simulations presented in this paper use a Cartesian grid of 512^3 cells and have been run on four 8-core workstations for four days to simulate ~5000s, which corresponds to almost ten convective ...
We develop a 3-D model to simulate the synthetic aperture radar (SAR) image formation process of an undulated vegetation canopy such as corn grown in fields with large periodic drainage reliefs. We explain how the simulated SAR image of undulated vegetation medium is obtained by the convolution of a 2-D slice of the 3-D simulated SAR system point spread function [(PSF), emulating the SAR beam modeled by a cosine modulated Gaussian], with the 2-D projection of the observed undulated vegetation canopy (modeled with scatterers randomly distributed in 3-D undulated space) followed by the extraction of each look envelope, the summation of looks, and sampling in azimuthal and range directions. Our model is useful to study the parameters involved in the formation and the analysis of SAR images of undulated vegetation medium. Validation of simulations made with actual SAR images shows that undulated corn crop ...
We present results from a large volume simulation of Hydrogen reionization. We combine 3d radiative transfer calculations and an N-body simulation, describing structure formation in the IGM, to detail the growth of HII regions around high redshift galaxies. Our simulation tracks 1024^3 dark matter particles, in a box of side length 65.6 Mpc/h. This large volume allows us to accurately characterize the size distribution of HII regions throughout most of the reionization process. At the same time, our simulation resolves many of the small galaxies likely responsible for reionization. It confirms a picture anticipated by analytic models: HII regions grow collectively around highly-clustered sources, and have a well-defined characteristic size, evolving from a sub-Mpc scale at the beginning of reionization to R>10 Mpc towards the end. We present a detailed statistical description of our results, and ...
This paper describes the graphical simulation of a traffic environment. The environment includes streets leading to an intersection, the intersection, vehicle traffic, and signal lights in the intersection controlled by different methods. The simulation allows for the study of parameters affecting traffic environments and the study of different control strategies for traffic signal lights, including conventional, fuzzy, and adaptive control methods. Realistic traffic environments are simulated including a cross intersection, with one or more lanes of traffic in each direction, with and without turn lanes. Vehicle traffic patterns are a mixture of cars going straight and making right or left turns. The free velocities of vehicles follow a normal distribution with a mean of the ``posted`` speed limit. Actual velocities depend on such factors as the proximity and velocity of surrounding traffic, approaches to intersections, ...
Numerical simulations predict that a large number of baryons reside in intergalactic space at temperatures between 10^5-10^7 K. Highly-ionized metals, such as O VII and O VIII, are good tracers of this ``warm-hot intergalactic medium'', or WHIM. For collisionally-ionized gas, the ionization fraction of each ion peaks at some particular temperature (``peak temperatures''), so different ions can therefore trace the IGM at different temperatures. We performed a hydrodynamic simulation to study the metal distributions in the IGM. We then draw random lines-of-sight across the simulated region and synthesize resonance absorption line spectra in a similar way to simulating the Ly-alpha forest. By studying the distribution functions of H- and He-like O, Si and Fe in a collisionally-ionized IGM and comparing with semi-analytic results based on the Press-Schechter formalism, we find: (1) ions with higher peak ...
The improvement in the functions of the viscous flow calculation method VEGA-SHIP around a ship and the expansion of application range were described as the numerical flow simulation in ship and ocean engineering and at the same time application examples to the ocean engineering by the general-purpose flow simulation code FLOW-3D handling the non-steady flow with a free surface were introduced as the numerical simulation regarding such products as a water gate and a dam. In the VEGA-SHIP, water surface was handled as a fixed wall so that wave could not be calculated. Therefore, an algorithm for calculating wave on the water surface was added to the VEGA-SHIP and a calculation method simultaneously considering the creation of wave around the ship and viscosity was developed. The FLOW-3D was used to calculate the phenomenon where inside liquid moved greatly due to the oscillation of a tank and hit against and damaged the tank ...
The three-dimensional (3-D), multiphase, computational fluid dynamic (CFD) code FLUENT is used to simulated two-phase flow behaviour in a CANDU header manifold under low (natural circulation) flow conditions. This behaviour was previously inferred from experimental data. The CFD simulations reported here are being used to support these inferences and to obtain a better understanding of phase distribution in the header manifold. The simulations seem to show that the vapor-water mixture models in the FLUENT code do not capture properly phase separation in the header and proper phase branching at the header-feeder connections that have been observed in experiments at low flows. The simulations using discrete-phase model in FLUENT, which tracks the pathlines of the individual vapor bubbles in the water continuum phase, show interesting, complicated and, in some cases, unexpected bubble trajectories from the ...
The Stockpile Stewardship Program (SSP) is a single, highly integrated technical program for maintaining the safety and reliability of the U.S. nuclear stockpile. The SSP uses past nuclear test data along with future non-nuclear test data, computational modeling and simulation, and experimental facilities to advance understanding of nuclear weapons. It includes stockpile surveillance, experimental research, development and engineering programs, and an appropriately scaled production capability to support stockpile requirements. This integrated national program will require the continued use of current facilities and programs along with new experimental facilities and computational enhancements to support these programs. The Advanced Simulation and Computing program (ASC) is a cornerstone of the SSP, providing simulation capabilities and computational resources to support the annual stockpile assessment and certification, to ...
(Abridged) We use 3D SPH calculations with higher resolution, as well as with more realistic viscosity and sound-speed prescriptions than previous work to examine the eccentric instability which underlies the superhump phenomenon in semi-detached binaries. We illustrate the importance of the two-armed spiral mode in the generation of superhumps. Differential motions in the fluid disc cause converging flows which lead to strong spiral shocks once each superhump cycle. The dissipation associated with these shocks powers the superhump. We compare 2D and 3D results, and conclude that 3D simulations are necessary to faithfully simulate the disc dynamics. We ran our simulations for unprecedented durations, so that an eccentric equilibrium is established except at high mass ratios where the growth rate of the instability is very low. Our improved simulations give a closer match to the observed relationship ...
In a Pb-Bi cooled direct contact steam generation fast reactor water is injected directly above the core, the produced steam is separated at the top and is send to the turbine. Neither the direct contact phenomenon nor the two-phase flow simulations in CFD have been thoroughly described yet. A first attempt in simulating such two-phase flow in 2D using the CFD code Fluent is presented in this paper. The volume of fluid explicit model was used. Other important simulation parameters were: pressure velocity relation PISO, discretization scheme body force weighted for pressure, second order upwind for momentum and CISCAM for void fraction. Boundary conditions were mass flow inlet (Pb-Bi 0 kg/s and steam 0.07 kg/s) and pressure outlet. The effect of mesh size (0.5 mm and 0.2 mm cells) was investigated as well as the effect of the turbulent model. It was found that using a fine mesh is very important in order to achieve larger ...
Orientational constraints obtained from solid state NMR experiments on anisotropic samples are used here in molecular dynamics (MD) simulations for determining the structure and dynamics of several different membrane-bound molecules. The new MD technique is based on the inclusion of orientation dependent pseudo-forces in the COSMOS-NMR force field. These forces drive molecular rotations and re-orientations in the simulation, such that the motional time-averages of the tensorial NMR properties approach the experimentally measured parameters. The orientational-constraint-driven MD simulations are universally applicable to all NMR interaction tensors, such as chemical shifts, dipolar couplings and quadrupolar interactions. The strategy does not depend on the initial choice of coordinates, and is in principle suitable for any flexible molecule. To test the method on three systems of increasing complexity, we used as constraints ...
Orientational constraints obtained from solid state NMR experiments on anisotropic samples are used here in molecular dynamics (MD) simulations for determining the structure and dynamics of several different membrane-bound molecules. The new MD technique is based on the inclusion of orientation dependent pseudo-forces in the COSMOS-NMR force field. These forces drive molecular rotations and re-orientations in the simulation, such that the motional time-averages of the tensorial NMR properties approach the experimentally measured parameters. The orientational-constraint-driven MD simulations are universally applicable to all NMR interaction tensors, such as chemical shifts, dipolar couplings and quadrupolar interactions. The strategy does not depend on the initial choice of coordinates, and is in principle suitable for any flexible molecule. To test the method on three systems of increasing complexity, we used as constraints ...
Orientational constraints obtained from solid state NMR experiments on anisotropic samples are used here in molecular dynamics (MD) simulations for determining the structure and dynamics of several different membrane-bound molecules. The new MD technique is based on the inclusion of orientation dependent pseudo-forces in the COSMOS-NMR force field. These forces drive molecular rotations and re-orientations in the simulation, such that the motional time-averages of the tensorial NMR properties approach the experimentally measured parameters. The orientational-constraint-driven MD simulations are universally applicable to all NMR interaction tensors, such as chemical shifts, dipolar couplings and quadrupolar interactions. The strategy does not depend on the initial choice of coordinates, and is in principle suitable for any flexible molecule. To test the method on three systems of increasing complexity, we used as constraints ...
This paper presents a new methodology to create realistic 3D microstructures of polycrystals. The virtual microstructures are based on statistical data describing the morphological and crystallographic textures of a sample, obtained from an EBSD analysis. In addition, the methodology can reproduce the observed surface on top of the simulated microstructure. This feature allows finite element calculations on these virtual aggregates to be compared to experimental results of mechanical tests. Such a comparison leads to the identification of the mechanical parameters of constitutive laws, such as critical resolved shear stress and strain hardening, using an optimization algorithm. Two materials were simulated in this study: TiAl and grade 702 zirconium. The first one presents twins inside the microstructure and the second one has an anisotropic texture. Based on 2D simulations, the important parameters necessary to describe a ...
This thesis - Resource Management in Broadband Communication Networks - deals with different ways of optimizing the available resources of data- or telecommunication networks. Especially topics like optimal routing, load balancing and fast recovery of routes in case of link failures are covered. The first part gives a brief description of some of the existing protocols for routing and controlling resources, such as RSVP, OSPF, BGP, PNNI, etc. The remaining part concerns the following two fundamentally different approaches to resource management etc.: - Software agents - Simulated ants In the beginning of the part concerning software agents a description of what exactly constitutes software agents, according to the scientists and according to the organization, FIPA is given. After this, the main results from the IMPACT project are presented. The IMPACT project is an EU-project that aimed at developing a demonstration platform, where software agents handled virtually ...
The authors analyze the consequences of models of structure formation for higher-order (n-point) galaxy correlation functions in the mildly non-linear regime. Several variations of the standard {Omega} = 1 cold dark matter model with scale-invariant primordial perturbations have recently been introduced to obtain more power on large scales, R{sub p} {approximately}20 h{sup {minus}1} Mpc, e.g., low-matter-density (non-zero cosmological constant) models, {open_quote}tilted{close_quote} primordial spectra, and scenarios with a mixture of cold and hot dark matter. They also include models with an effective scale-dependent bias, such as the cooperative galaxy formation scenario of Bower, et al. The authors show that higher-order (n-point) galaxy correlation functions can provide a useful test of such models and can discriminate between models with true large-scale power in the density field and those where the galaxy power arises from scale-dependent bias: a bias with ...
We use deep HST/ACS observations to calculate the star formation history (SFH) of the Cetus dwarf spheroidal (dSph) galaxy. Our photometry reaches below the oldest main-sequence turnoffs, which allows us to estimate the age and duration of the main episode of star formation in Cetus. This is well approximated by a single episode that peaked roughly 12 #+-# 0.5 Gyr ago and lasted no longer than about 1.9 #+-# 0.5 Gyr (FWHM). Our solution also suggests that essentially no stars formed in Cetus during the past 8 Gyr. This makes Cetus' SFH comparable to that of the oldest Milky Way dSphs. Given the current isolation of Cetus in the outer fringes of the Local Group, the dominant old population implies that Cetus is a clear outlier in the morphology-Galactocentric distance relation that holds for the majority of the Milky Way dwarf satellites. Our results also show that Cetus continued forming stars until z#approx =# 1, long after the universe was reionized, and that there is no clear ...
This Report of activity displays the work developed at Centre d`Etude Nucleaires de Bordeaux-Gradignan (CENBG) in the period 1995-1996. A number of 68 short progress notes are presented in the fields of Cosmic particles (3), Nuclear Astrophysics (4), High Spins-Large Deformations (8), Basic Interactions (9), Exotic Nuclei (10), Hybrid Systems (4), Theoretical Physics (28), Cellular Microanalysis by means of Nuclear Probes (3) and Technical Development (3). The main problems attacked at CENBG, established as orientations by scientific board of IN2P3 (CENBG is a UMR of Bordeaux-1 Univ. and CNRS-IN2P3) are the following: the research of Majorana/Dirac nature of neutrino (in the framework of NEMO Experiment) in relation with the problem of origin of mass in the Universe and the cosmological puzzle of dark matter, the investigation of the origin of high energy (20 to 200 GeV) cosmic radiation (in the framework of CELESTE Experiment), the study of extreme issues in ...
The main goals of nuclear astrophysics have been to probe the interiors of stars, stellar explosions, the early moments of cosmic expansion, and the formation and evolution of galaxies and cosmic structure by measurement and application of the relevant nuclear physics. The approach to these goals have generally been from three directions: 1) Careful measurements of the relevant nuclear reactions; 2) Detailed computer models of the relevant astrophysical environments; and 3) Observations of the relevant terrestrial and extra-terrestrial atomic and isotopic abundances. These approaches provide not only insight into the formation and evolution of the elements, but are also pillars upon which a variety of cosmological models as well as models for physics beyond the standard model of particle physics can stand or fall. At present there is a very exciting frontier on all three of these approaches. The development and applications of radioactive-ion-beam and ...
We have developed a kinetic Monte Carlo (kMC) simulator that links atomic migration and binding energies determined primarily from first principles calculations to macroscopic phenomena and laboratory time scales. Input for the kMC simulation is obtained from a combination of ab initio planewave pseudopotential calculations, molecular dynamics simulations, and experimental data. The simulator is validated against an extensive series of experimental studies of the diffusion of B spikes in self-implanted Si. The implant energy, dose, and dose rate, as well as the detailed thermal history of the sample, are included. Good agreement is obtained with the experimental data for temperatures between 750 and 950 C and times from 15 to 255 s. At 1050o C we predict too little diffusion after 105 s compared to experiment: apparently, some mechanism which is not adequately represented by our model becomes important ...
This report covers application of Argonne National Laboratory`s (ANL`s) computer codes to simulation and analysis of components of the magnetohydrodynamic (MHD) power train system at the Component Development and Integration Facility (CDIF). Major components of the system include a 50-MWt coal-fired, two-stage combustor and an MHD channel. The combustor, designed and built by TRW, includes a deswirl section between the first and the second-stage combustor and a converging nozzle following the second-stage combustor, which connects to the MHD channel. ANL used computer codes to simulate and analyze flow characteristics in various components of the MHD system. The first-stage swirl combustor was deemed a mature technology and, therefore, was not included in the computer simulation. Several versions of the ICOMFLO computer code were used for the deswirl section and second-stage combustor. The MGMHD code, upgraded with a slag ...