Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “N CODES”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

Pursuing Dark Energy with Large Galaxy Redshift Surveys: Baryon Acoustic Oscillations and Beyond (Final Technical Report)

We have developed and applied methods of large-scale structure from galaxy redshift surveys to the study of dark energy. This program supported my group's development of the baryon acoustic method and our participation in the Dark Energy Spectroscopic Instrument (DESI) survey. In addition to my leadership role in the DESI collaboration, my group at Harvard focused on several topics addressing key needs of the DESI science program. 1) We continued our design, validation, and documentation of DESI target selection of the luminous red galaxy and emission-line galaxy samples. 2) We applied simulations from our novel and extremely fast N-body code to make DESI catalogs. With these, we will search for possible systematic shifts in the acoustic scale due to clustering bias with an unprecedented volume of simulations. 3) We extended our method for the fast construction of covariance matrices for correlation functions. 4) We developed new large-scale structure methods, including density-field reconstruction and three-point function analysis.

79 ASTRONOMY AND ASTROPHYSICS↗

Improvements in MCNP6.3 (Rev. 1) [Slides]

This presentation outlines improvements in the Monte Carlo N-Particle Code (MCNP) version 6.3. MCNP is an internationally renowned Monte Carlo particle transport code developed by Los Alamos National Laboratory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Computational study of tungsten and depleted uranium photoneutron targets for a 20 MeV electron linear accelerator

Neutron production can be realized with a high energy electron linear accelerator by using Bremsstrahlung and photoneutron converters. In this study, Monte Carlo N-Particle Code (MCNP) was used to evaluate potential photonuclear target designs for a high energy electron linear accelerator for applications such as neutron radiography and neutron resonance spectroscopy. A computational model was developed to inform a target design that would yield a high number of neutrons. It consists of a 20 MeV electron beam incident on a Bremsstrahlung target and a photonuclear target to generate neutrons. This computational model showed that a thickness of 0.75 inches for both tungsten and depleted uranium yields the most neutrons from photoneutron reactions. Saturation in the total number of generated neutrons was observed at over 0.75-inch thickness for both evaluated materials. Depleted uranium yielded approximately twice the number of neutrons overall compared to tungsten. The highest neutron surface flux for Depleted Uranium was 1.06 × 10-4 neutrons/cm2/source electron, and for Tungsten it was 5.12 × 10-5 neutrons/cm2/source electron. The optimal target design for this study’s application would consist of a 0.75 inch-thick block of depleted uranium with the length, width, and/or diameter varying dependent on application.

43 PARTICLE ACCELERATORS↗

Fitness For Service Assessment of a Corroded Heat Exchanger

Within the Fermi National Accelerator complex, there exist various water systems that support accelerator operations. One of these systems is extremely vital to the operation of the machine; that is the cooling system. The cooling system consists of nine relatively large heat exchangers that take untreated pond water and use it to cool the process fluid that further cools machine components. Over the 30 years these heat exchangers have been in operation, they have undergone significant material loss on the channels. This material loss, due to various forms of corrosion such as galvanic and microbiologically influenced corrosion (MIC) and possibly others, has deteriorated more than 80% of the nominal wall thickness of some of the exchangers and placed them in a questionable state. ASME FFS-1 (API 579) has been applied to address the condition of the heat exchangers due to their noncompliance with the governing code, BPVC Sec. VIII Div. 1. The assessments encompassed ASME FFS-1 parts 4: General Metal Loss and 9: Crack Like Flaw using level 1, 2, and 3 analysis techniques based on inspection data obtained by API 510 inspections. Level 1 and 2 assessments were deemed unfit for the corroded regions due to their location relative to a major structural discontinuity (channel to tube-sheet joint), so a level 3 analysis was conducted according to ASME Sec. VIII Div. 2 (design by analysis) rules for pressure vessels. Supplemental information included pond water tests to determine an accurate future corrosion allowance due to lacking inspection history. A leak before break (LBB) route was chosen to evaluate the possibility of leaking prior to the onset of failure. The analysis of one heat exchanger shows that the possibility the channel will develop a pinhole leak over 2.5 more years of operation should not be overlooked, but burst was unlikely from operation. The use of fracture mechanics show, that if a through-wall crack were to develop, it would not propagate further than the channel geometry and cause a leak not greater than 35 GPM. Using ASME Section XI Code Case N-705-1, allowing us to operate with a leak until the next outage given certain operating conditions and developing a leak mitigation procedure, this heat exchanger is deemed fit-for-service.

Humenik, Alex [Fermilab]↗

Analysis of electrophoresis performance

The SAMPLE computer code models electrophoresis separation in a wide range of conditions. Results are included for steady three dimensional continuous flow electrophoresis (CFE), time dependent gel and acetate film experiments in one or two dimensions and isoelectric focusing in one dimension. The code evolves N two dimensional radical concentration distributions in time, or distance down a CFE chamber. For each time or distance increment, there are six stages, successively obtaining the pH distribution, the corresponding degrees of ionization for each radical, the conductivity, the electric field and current distribution, and the flux components in each direction for each separate radical. The final stage is to update the radical concentrations. The model formulation for ion motion in an electric field ignores activity effects, and is valid only for low concentrations; for larger concentrations the conductivity is, therefore, also invalid.

Roberts, G. O.↗

The distribution of dark matter, galaxies, and the intergalactic medium in a cold dark matter dominated universe

The evolution and distribution of galaxies and the intergalactic medium (IGM) have been studied, along with collisionless dark matter in a Universe dominated by cold dark matter. The Einstein-deSitter universe with omega sub 0 = 1 and h = 0.5 was considered (here h = H sub 0 bar 100/kms/Mpc and H sub 0 is the present value of the Hubble constant). It is assumed that initially dark matter composes 90 pct and baryonic matter composes 10 pct of total mass, and that the primordial baryonic matter is comprised of H and He, with the abundance of He equal to 10 pct of H by number. Galaxies are allowed to form out of the IGM, if the total density and baryonic density satisfy an overdensity criterion. Subsequently, the newly formed galaxies release 10 to the 60th ergs of energy into the IGM over a period of 10 to the 8th years. Calculations have been performed with 32 to the 3rd dark matter particles and 32 to the 3rd cells in a cube with comoving side length L = 9.6/h Mpc. Dark matter particles and galaxies have been followed with an N-body code, while the IGM has been followed with a fluid code.

Ryu, Dongsu↗

Long-lived spiral waves in N-body simulations

Results of N-body simulations of disc galaxies using a two-dimensional Cartesian N-body code are presented. Both trailing arm spirals (TAS) and leading arm spirals (LAS) were used with varieties of pitch angles and pattern speeds. LAS perturbations transferred their energy to TAS via swing amplification; TAS perturbations led to TAS arms. In both cases the spiral arms persisted for more than 5 rotation periods, but the maximum amplitude 2-armed spirals were generated by LAS perturbations. The persistence of the trailing arm spiral waves is thought to be caused by the kinematic spiral arm mechanism described by Kalnajs (1973).

Comins, Neil F.↗

A simulation survey of galaxy interactions

Many carefully selected samples of interacting galaxies have been observed extensively in attempts to clarify whether interaction produces activity in galaxies. Because the sample members represent a wide range of encounter parameters and times, one can then study whether there are correlations between observable encounter features and, for example, Seyfert activity. On the other hand, in theoretical studies, simulations typically deal with either time-consuming detailed modelling of single galaxy pairs or tracing a few model encounters over time. The authors extend the observational survey approach by combining it with a simulation survey. The authors are conducting a survey of model encounters, covering the most important encounter parameters over a wide range. Some parameters, such as companion structure and initial velocity, are demonstratably less important and can be ignored in a first pass. The parameter range must be richly enough sampled so that the authors can evaluate the uniqueness of the observable morphology and velocity structure of the resulting simulated pairs to diagnose unobservable companion orbit parameters. They are using a self-gravitating polar n-body code run on the Cray X-MP at the Alabama Supercomputer Network. For each simulation, the authors have stellar and gas distributions predicted over, typically, a billion years, along with information on gas motions within the disk and any material captured by the companion or lost to the system. Features of disturbed spiral galaxies are sensitive enough to time and encounter parameters so that a match of the simulation survey results to observations can be applied as starting points to infer unobservable orbital or system parameters in actual sample members. This should enable them to examine whether interesting observed properties (Seyfert activity, nuclear star-formation rate) are functions of unobservable dynamical properties which characterize each encounter. Any correlations (or lack of some expected ones) will provide strong clues as to how or whether these phenomena are related to interactions. Aside from its use with such observed samples, this survey should greatly speed determination of initial orbital parameters for more detailed subsequent simulations of individual systems.

Byrd, Gene G.↗

Simulations of gas clouds in interacting galaxies

A companion can induce a variety of morphological changes in a galaxy. The author uses N-body simulations to study the effects of different kinds of perturbations on the dynamics of a disk galaxy. The model is two-dimensional, with a disk consisting of about 60,000 particles. Most of the particles (80%) represent the old stellar population with a high velocity dispersion, while the rest (20%) represent gas clouds with a low velocity dispersion. Initially, the velocity dispersion corresponds to Q = 1 for the star particles, and Q = O for the gas particles, where Q is Toomre's (1964) stability parameter. The gas clouds can collide inelastically. The disk is stabilized by a rigid halo potential, and by the random motions of the old star particles. To simulate the effect of an encounter on the disk, a companion galaxy, modelled as a point mass, can move in a co-planar orbit around the disk. A complete description of the N-body code is found in Thomasson (1989). The spiral structures caused by a companion in first a direct and then a retrograde (with respect to the rotation of the disk) parabolic orbit are presented. The associated velocity fields suggest a way to observationally distinguish between leading and trailing spiral arms. The stability of the gas component in a disk in which tidally triggered infall of gas to the center occurs is studied. Finally, the author shows how a ring of gas can form in a disk as a result of a co-planar encounter with another galaxy.

Thomasson, Magnus↗

Heat transfer and pressure measurements and comparison with prediction for the SSME two-stage turbine

Time averaged Stanton number and surface pressure distributions are reported for the first stage vane row, the first stage blade row, and the second stage vane row of the Rocketdyne Space Shuttle Main Engine (SSME) two-stage fuel-side turbine. Unsteady pressure envelope measurements for the first blade are also reported. These measurements were made at 10 percent, 50 percent, and 90 percent span on both the pressure and suction surfaces of the first stage components. Additional Stanton number measurements were made on the first stage blade platform, blade tip, and shroud, and at 50 percent span on the second vane. A shock tube was used as a short duration source of heated and pressurized air to which the turbine was subjected. Platinum thin film heat flux gages were used to obtain the heat flux measurements, while miniature silicon diaphragm flush-mounted pressure transducers were used to obtain the pressure measurements. The first stage vane Stanton number distributions are compared with predictions obtained using a version of STAN5 and quasi-3D Navier-Stokes solution. This same quasi-3D N-S code was also used to obtain predictions for the first blade and the second vane.

Dunn, M. G.↗

The spatial correlation properties of dark galaxy halos in a CDM universe

We use the Hierarchical Particle Mesh (HPM) N-body code written by J. V. Villumsen (Villumsen, 1989) to investigate the two-point spatial correlation function, xi(r), of dark galaxy halos as a function of halo mass and local environment (i.e. high, low, or average mass density). We assume a standard cold dark matter (CDM) universe (omega = 1, delta = 0, H sub 0 = 50,km/sec/Mpc). Because of the large dynamic ranges in mass and length that can be obtained with the HPM code, it is well-suited to an investigation of this sort.

Brainerd, Tereasa G.↗

Effect of tidal fields on star clusters

We follow the dynamical evolution of a star cluster in a galactic tidal field using a restricted N-body code. We find large asymmetric distortions in the outer profile of the cluster in the first 10 or so crossing times as material is lost. Prograde stars escape preferentially and establish a potentially observable retrograde rotation in the halo. We present the rate of particle loss and compare with the prescription proposed by Lee and Ostriker (1987).

Chernoff, David↗

Effects of gas on the global stability of galactic disks - Radial flows

We study numerically the effect of gas on the global stability of a two-component self-gravitating galactic disk embedded in a live halo. The stars are evolved by using a 3D collisionless N-body code, and the gas is represented by an ensemble of finite size inelastic particles. The gravitational interaction of stars and gas is calculated using a TREE method. We find that the evolution of the gaseous distribution in the globally unstable disks can be described by two different regimes. When the gas mass fraction is less than about 10 percent, the gas is channeled toward the galactic center by a growing stellar bar. For higher gas fractions, the gas becomes highly inhomogeneous, and the bar instability in the disk is heavily damped. The gas falls toward the inner kpc due to dynamical friction. Domains of both regimes depend on the efficiency of dissipation in the gas. We also discuss the relevance of the Jeans instability and give an empirical criterion for the global bar instability in a two-component self-gravitating disk.

Shlosman, Isaac↗

Experimental measurements and analytical analysis related to gas turbine heat transfer. Part 1: Time-averaged heat-flux and surface-pressure measurements on the vanes and blades of the SSME fuel-side turbine and comparison with prediction. Part 2: Phase-resolved surface-pressure and heat-flux measurements on the first blade of the SSME fuel-side turbine

Time averaged Stanton number and surface-pressure distributions are reported for the first-stage vane row, the first stage blade row, and the second stage vane row of the Rocketdyne Space Shuttle Main Engine two-stage fuel-side turbine. Unsteady pressure envelope measurements for the first blade are also reported. These measurements were made at 10 percent, 50 percent, and 90 percent span on both the pressure and suction surfaces of the first stage components. Additional Stanton number measurements were made on the first stage blade platform blade tip, and shroud, and at 50 percent span on the second vane. A shock tube was used as a short duration source of heated and pressurized air to which the turbine was subjected. Platinum thin-film heat flux gages were used to obtain the heat flux measurements, while miniature silicon-diaphragm flush-mounted pressure transducers were used to obtain the pressure measurements. The first stage vane Stanton number distributions are compared with predictions obtained using a version of STAN5 and a quasi-3D Navier-Stokes solution. This same quasi-3D N-S code was also used to obtain predictions for the first blade and the second vane.

Source record↗

The jet model of chondrule formation

We estimate the size range of particles that are ejected from a protostellar accretion disk by a protostellar jet. An n-body code is used to determine the subsequent motion of the ejected particles, where the particles are subject to two forces: the gravitational attraction from the protostar; and the gas drag from the halo gas of the accretion disk.

Liffman, K.↗

Self-Consistent Simulation of the Brownian Stage of Dust Growth

It is a widely accepted view that in proto-planetary accretion disks the collision and following sticking of dust particles embedded in the gas eventually leads to the formation of planetesimals (coagulation). For the smallest dust grains, Brownian motion is assumed to be the dominant source of their relative velocities leading to collisions between these dust grains. As the dust grains grow they eventually couple to the turbulent motion of the gas which then drives the coagulation much more efficiently. Many numerical coagulation simulations have been carried out to calculate the fractal dimension of the aggregates, which determines the duration of the ineffective Brownian stage of growth. Predominantly on-lattice and off-lattice methods were used. However, both methods require simplification of the astrophysical conditions. The aggregates found by those methods had a fractal dimension of approximately 2 which is equivalent to a constant, mass-independent friction time. If this value were valid for the conditions in an accretion disk, this would mean that the coagulation process would finally 'freeze out' and the growth of a planetesimal would be impossible within the lifetime of an accretion disk. In order to investigate whether this fractal dimension is model independent, we simulate self-consistently the Brownian stage of the coagulation by an N-particle code. This method has the advantage that no further assumptions about homogeneity of the dust have to be made. In our model, the dust grains are considered as aggregates built up of spheres. The equation of motion of the dust grains is based on the probability density for the diffusive transport within the gas atmosphere. Because of the very low number density of the dust grains, only 2-body-collisions have to be considered. As the Brownian stage of growth is very inefficient, the system is to be simulated over long periods of time. In order to find close particle pairs of the system which are most likely to undergo a collision, we use a particle-in-cell (PIC) method for the early stages of the simulation where the system is still very homogeneous and a tree method later when the particles are more clustered.

Kempf, S.↗

Improved models for the evolution of the Coma cluster of galaxies

The analysis by Fitchett & Webster of the observations of the Coma cluster of galaxies has demonstrated that the center of the Coma Cluster consists of two subclusters. Therefore, it is important to construct realistic dynamical models of a galaxy cluster with two mass centers. Our previous N-body models for the Coma Cluster consisted of point masses or particles with simple interaction properties. In the current paper, we employ a more sophisticated N-body code, which includes dynamical friction, mass exchange, and mergers between galaxies.

Dark↗

Design and Validation of a Gammaray Scanning System for Measuring Irradiated Nuclear Fuel

A prototype gamma-ray scanning system was designed to perform post-irradiation measurements of nuclear fuel at Idaho National Laboratory (INL). The system is composed of a coaxial high-purity germanium (HPGe) detector, a collimator, and mechanical positioning stages that translate the fuel sample across the front of the collimator. This paper documents the design process, modeling and the laboratory characterization performed to validate the system. The modeling of the system was done using the Monte Carlo N-Particle code with vendor-supplied design specifications, dimensional measurements, and x-ray radiographs of the HPGe detector. The characterization included measurements of calibrated radiation sources in different geometries, and calculations for absolute efficiency and dead time. Benchmark measurements were performed by scanning an irradiated fuel rodlet containing ten pellets of 4.9% enriched UO2 in zirconium alloy. These results were compared to simulated spectra to further help characterize the detector model. One goal of the system is to determine the number of fissions per gram of UO2 in the fuel; comparison of the results from this prototype system's assay of the irradiated fuel rodlet are in agreement with previously estimated results.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗