Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computerized simulation”

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 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Electrochemical Corrosion Studies for Modeling Metallic Waste Form Release Rates

The isotope 99 Tc is an important fission product generated from nuclear power production. Because of its long half-life (t 1/2 = 2.13 ∙ 10 5 years) and beta-radiotoxicity (β⁻ = 292 keV), it is a major concern in the long-term management of spent nuclear fuel. In the spent nuclear fuel, Tc is present as an alloy with Mo, Ru, Rh, and Pd called the epsilon-phase, the relative amount of which increases with fuel burn-up. In some separation schemes for spent nuclear fuel, Tc would be separated from the spent fuel and disposed of in a durable waste form. Technetium waste forms under consideration include metallic alloys, oxide ceramics and borosilicate glass. In the development of a metallic waste form, after separation from the spent fuel, Tc would be converted to the metal, incorporated into an alloy and the resulting waste form stored in a repository. Metallic alloys under consideration include Tc–Zr alloys, Tc–stainless steel alloys and Tc–Inconel alloys (Inconel is an alloy of Ni, Cr and iron which is resistant to corrosion). To predict the long-term behavior of the metallic Tc waste form, understanding the corrosion properties of Tc metal and Tc alloys in various chemical environments is needed, but efforts to model the behavior of Tc metallic alloys are limited. One parameter that should also be considered in predicting the long-term behavior of the Tc waste form is the ingrowth of stable Ru that occurs from the radioactive decay of 99 Tc ( 99 Tc → 99 Ru + β⁻). After a geological period of time, significant amounts of Ru will be present in the Tc and may affect its corrosion properties. Studying the effect of Ru on the corrosion behavior of Tc is also of importance. In this context, we studied the electrochemical behavior of Tc metal, Tc-Ni alloys (to model Tc-Inconel alloy) and Tc-Ru alloys in acidic media. The study of Tc-U alloys has also been performed in order to better understand the nature of Tc in metallic spent fuel. Computational modeling and simulations were performed to shed light on experimental results and explain structural and kinetics trends.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2016

ARES Modeling of High-foot Implosions (NNSA Milestone #5466)

ARES “capsule only” simulations demonstrated results of applying an ASC code to a suite of high-foot ICF implosion experiments. While a capability to apply an asymmetric FDS drive to the capsule-only model using add-on Python routines exists, it was not exercised here. The ARES simulation results resemble the results from HYDRA simulations documented in A. Kritcher, et al., Phys. Plasmas, 23, 052709 (2016); namely, 1D simulation and data are in reasonable agreement for the lowest velocity experiments, but diverge from each other at higher velocities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Nonparametric estimators for incomplete surveys

Nonparametric estimators, such as the 1/V{sub max} estimator and the C{sup −} estimator, have been applied extensively to estimate luminosity functions (LFs) of astronomical sources from complete, truncated survey data sets. Application of such estimators to incomplete data sets typically requires further truncation of data, separation into subsets of constant completeness, and/or correction for incompleteness-induced bias. In this paper, we derive generalizations of the above estimators designed for use with incomplete, truncated data sets. We compare these generalized nonparametric estimators, investigate some of their simple statistical properties, and validate them using Monte Carlo simulation methods. We apply a nonparametric estimator to data obtained from the extended Baryon Oscillation Spectroscopic Survey to estimate the QSO LF for redshifts 0.68<z<4.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

The Effect of Corner Modes in the Initial Conditions of Cosmological Simulations

In view of future high-precision large-scale structure surveys, it is important to quantify the percent and subpercent level effects in cosmological N-body simulations from which theoretical predictions are drawn. One such effect involves deciding whether to zero all modes above the one-dimensional Nyquist frequency, the so-called “corner” modes, in the initial conditions. We investigate this effect by comparing power spectra, density distribution functions, halo mass functions, and halo profiles in simulations with and without these modes. For a simulation with a mass resolution of m{sub p}∼10{sup 11} h{sup −1} M{sub ⊙}, we find that at z>6, the difference in the matter power spectrum is large at wavenumbers above ∼80% of k{sub Ny}, reducing to below 2% at all scales by z∼3. Including corner modes results in a better match between low- and high-resolution simulations at wavenumbers around the Nyquist frequency of the low-resolution simulation, but the effect of the corner modes is smaller than the effect of particle discreteness. The differences in mass functions are 3% for the smallest halos at z = 6 for the m{sub p}∼10{sup 11} h{sup −1} M{sub ⊙} simulation, but we find no significant difference in the stacked profiles of well-resolved halos at z⩽6. Thus removing power at |k|>k{sub Ny} in the initial conditions of cosmological simulations has a small effect on small scales and high redshifts, typically below a few percent.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Implications of Galaxy Buildup for Putative IMF Variations in Massive Galaxies

Recent observational evidence for initial mass function (IMF) variations in massive quiescent galaxies at z = 0 challenges the long-established paradigm of a universal IMF. While a few theoretical models relate the IMF to birth cloud conditions, the physical driver underlying these putative IMF variations is still largely unclear. Here we use post-processing analysis of the Illustris cosmological hydrodynamical simulation to investigate possible physical origins of IMF variability with galactic properties. We do so by tagging stellar particles in the simulation (each representing a stellar population of ≈10{sup 6} M{sub ⊙}) with individual IMFs that depend on various physical conditions, such as velocity dispersion, metallicity, or star formation rate, at the time and place in which the stars are formed. We then follow the assembly of these populations throughout cosmic time and reconstruct the overall IMF of each z = 0 galaxy from the many distinct IMFs it is composed of. Our main result is that applying the observed relations between IMF and galactic properties to the conditions at the star formation sites does not result in strong enough IMF variations between z = 0 galaxies. Steeper physical IMF relations are required for reproducing the observed IMF trends, and some stellar populations must form with more extreme IMFs than those observed. The origin of this result is the hierarchical nature of massive galaxy assembly, and it has implications for the reliability of the strong observed trends, for the ability of cosmological simulations to capture certain physical conditions in galaxies, and for theories of star formation aiming to explain the physical origin of a variable IMF.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

The Populations of Carina. I. Decoding the Color–Magnitude Diagram

We investigate the color–magnitude diagram (CMD) of the Carina dwarf spheroidal galaxy using data of Stetson et al. and synthetic CMDs based on isochrones of Dotter et al., in terms of the parameters [Fe/H], age, and [ α /Fe], for the cases when (i) [ α /Fe] is held constant and (ii) [ α /Fe] is varied. The data are well described by four basic epochs of star formation, having [Fe/H] = −1.85, −1.5, −1.2, and ∼−1.15 and ages ∼13, 7, ∼3.5, and ∼1.5 Gyr, respectively (for [ α /Fe] = 0.1, constant [ α /Fe], and [ α /Fe] = 0.2, 0.1, −0.2, −0.2, variable [ α /Fe]), with small spreads in [Fe/H] and age of order 0.1 dex and 1–3 Gyr. Within an elliptical radius 13.′1, the mass fractions of the populations, at their times of formation, were (in decreasing age order) 0.34, 0.39, 0.23, and 0.04. This formalism reproduces five observed CMD features (two distinct subgiant branches of old and intermediate-age populations, two younger, main-sequence components, and the small color dispersion on the red giant branch (RGB). The parameters of the youngest population are less certain than those of the others, and given it is less centrally concentrated, it may not be directly related to them. High-resolution spectroscopically analyzed RGB samples appear statistically incomplete compared with those selected using radial velocity, which contain bluer stars comprising ∼5%–10% of the samples. We conjecture these objects may, at least in part, be members of the youngest population. We use the CMD simulations to obtain insight into the population structure of Carina's upper RGB.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Development of an Integrated Hydrological Model for Simulation of Surface Runoff and Stream Flow in Tims Branch Watershed, Savannah River Site, South Carolina - 16203

This research is part of continued efforts to develop an integrated transport model to predict the distribution of tin within the overland flow in Tims Branch Watershed (TBW) at Savannah River Site (SRS), South Carolina. Tin was introduced into TBW during the application of an innovative mercury remediation technology implemented by the U.S. Department of Energy's Office of Environmental Management, which involved the injection of stannous (tin) chloride into mercury contaminated groundwater. Understanding the fate and transport of tin and its compounds is of primary importance due to their potential impact on the environment. In this study, we developed an overland hydrology model (MIKE SHE) capable of simulating surface hydrology throughout the TBW. The modeling application used historical precipitation, groundwater levels, geological data, and river discharges that were retrieved from government databases and input to the model. The model was developed to simulate spatio-temporal distribution of flow discharges, flow duration, and water levels in the TBW. The preliminary simulation results indicated that the model was capable of predicting flow depth and velocity within the study area during extreme climate events. The developed hydrology model will later on be coupled with transport model to simulate tin distribution within Tims Branch stream during seasonal storm events. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2019

An Integrated Hydrological Model for Long-Term Monitoring in Tims Branch Watershed, SC - 19435

The Tims Branch watershed in Savannah River Site (SRS), South Carolina is a small-scale watershed that has been impacted by contaminated runoff and discharge wastewater from nuclear research facilities located within the watershed over the past century. Innovative remediation technologies, such as tin injection and pump and treat, have been implemented to remove heavy metals and volatile organic compounds (VOCs) from the watershed. Today, the Tims Branch watershed serves as a test-bed for monitoring the performance of remediation activities. This study describes the development of the two-dimensional integrated surface water/groundwater model for Tims Branch watershed to assess the impact of cleanup activities on hydrological processes in this watershed. This model will serve as a tool to simulate the effect of climatic events on overland flow distribution and further on will be used to assess contaminant transport. This paper illustrates the application of the developed model for simulation of overland flow depth and velocity in Tims Branch watershed and streamflow in A-014 outfall tributary. The watershed includes the second order Tims Branch stream, several outfalls, and a network of tributaries. Simulation of non-calibrated overland flow was conducted for the years 2016 using available stream flow data within the watershed. Recent flow data from a newly deployed stage gauge was also used to assess the model results. The results show the model ability to simulate overland flow and stream flow and can be used as a tool to represent the hydrological conditions in the Tims Branch watershed. The results of this model simulation will be further used to develop contaminant transport model. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2016

Evaluation of fission product releases obtained from KueFA safety tests in light of the setup's accuracy - 18692

The Cold Finger Apparatus (KuehlFinger-Apparatur - KueFA) at the European Commission's JRC-Karlsruhe has been designed to assess the effects of Depressurization Loss of Forced Circulation (D-LOFC) accident scenarios on spherical irradiated High Temperature Reactor (HTR) fuel. At approximately 1600 deg. C, the reference maximum temperature for D-LOFC accidents, it has been observed that high-quality German HTR fuel pebbles release only a limited quantity of Cs but a significant amount of Ag. This paper revisits fission product release data obtained from the recently KueFA-tested fuel pebbles HFR-K5/3 and HFR-EU1/3 and presents associated investigations to improve the accuracy of the experimental setup. In light of these findings, implications for previously conducted KueFA heating tests are discussed, considering both experimental data as well as simulated release curves calculated with the Source Term Analysis Code system (STACY). (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2018

Rossi and Feynman-alpha formulas including prompt neutron decay and detector die away time constants - 25396

A simple Rossi-alpha formula and the corresponding Feynman-alpha formula have been derived that take into account two time constants; one characterizing the prompt neutron decay in the fissile source domain and the other the detector die-away time. The derivation of the formula is an extension of the classic heuristic method. A test problem is defined and simulated using the Monte Carlo MCNPX code in order to illustrate the applicability of the formulas. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2018

Synergistic Effect of Nitrogen and Refractory Material on TiN Formation and Equiaxed Grain Structure of Ferritic Stainless Steel

The effect of nitrogen content on the formation of an equiaxed solidification structure of Fe-16Cr steel was investigated. Moreover, two different kinds of refractory materials, i.e., alumina and magnesia, were employed to control the type of oxide inclusion. The characteristics of TiN(-oxide) inclusions were quantitatively analyzed in both molten steel and solidified samples. When the melting was carried out in the alumina refractory, the grain size continuously decreased with increasing nitrogen content. However, a minimum grain size was observed at a specific nitrogen content (approx. 150 ppm) when the steel was melted in the magnesia refractory. Most of the single TiN particles had a cuboidal shape and fine irregularly shaped particles were located along the grain boundary due to the microsegregation of Ti at the grain boundary during solidification. The type of TiN-oxide hybrid inclusion was strongly affected by the refractory material where Al{sub 2}O{sub 3}-TiN and MgAl{sub 2}O{sub 4}-TiN hybrid-type inclusions were obtained in the alumina and magnesia refractory experiments, respectively. The formation of oxide inclusions was well predicted by thermochemical computations and it was commonly found that oxide particles were initially formed, followed by the nucleation and growth of TiN. When the nitrogen content increased, the number density of TiN linearly increased in the alumina refractory experiments. However, the number of TiN exhibits a maximum at about [N] = 150 ppm, at which a minimum grain size was obtained in the magnesia refractory experiments. Therefore, the larger the number density of TiN, the smaller the primary grain size after solidification. The number density of TiN in the steel melted in the magnesia refractory was greater than that in the steel melted in the alumina refractory at given Ti and N contents, which was due to the lower planar lattice disregistry of MgAl{sub 2}O{sub 4}-TiN interface rather than that of Al{sub 2}O{sub 3}-TiN interface. When ∆T{sub TiN} (= difference between the TiN precipitation temperature and the liquidus of the steel) was 20 K to 40 K, the number density of effective TiN was maximized and thus, the grain size was minimized after solidification. Finally, although most of the TiN particles were smaller than 1 μm in the molten steel samples irrespective of the nitrogen content, TiN particles larger than 10 μm were observed in the solidified samples when the nitrogen content was greater than 150 ppm. The growth of TiN particles during melting and solidification was well predicted by the combinatorial simulation of the ‘Ostwald ripening model’ based on the Lifshitz–Slyozov–Wagner theory in conjunction with the ‘Diffusion controlled model’ using Ohnaka’s microsegregation equation.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Removal of Bifilms from Al Melts by Stirring in Unbaffled and Baffled Crucibles

The free surface of aluminum melt, during stirring under different rotation speeds in unbaffled and baffled crucibles, was simulated numerically using a RANS turbulence model coupled with a particle-tracking method. The effect of mechanical stirring of Al melts containing 0.3, 0.7, and 4.5 wt pct Mg in unbaffled and baffled crucibles on the rate of elimination of bifilms from the melt was also investigated using a reduced pressure test. The results showed that the stirring increased the rate of removal of bifilms from the melt due to an increase in the rate of formation of cracks on the bifilms. Using baffles in the crucible accelerated the rate of removal of bifilms by changing the pattern of movement of bifilms in the melt. The critical stirring speeds, above which the oxide layer on the Al melt could submerge into the bulk liquid, were determined by the simulation. It was confirmed by the experimental results to be about 158 and 420 RPM for the unbaffled and baffled crucibles, respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2013

A proposal for the realization of universal quantum gates via superconducting qubits inside a cavity

A family of quantum logic gates is proposed via superconducting (SC) qubits coupled to a SC-cavity. The Hamiltonian for SC-charge qubits inside a single mode cavity is considered. Three- and two-qubit operations are generated by applying a classical magnetic field with the flux. Therefore, a number of quantum logic gates are realized. Numerical simulations and calculation of the fidelity are used to prove the success of these operations for these gates. -- Highlights: •A family of quantum logic gates is proposed via SC-qubits coupled to a cavity. •Three- and two-qubit operations are generated via a classical field with the flux. •Numerical simulations and calculation of the fidelity are used to prove the success of these operations for these gates.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

A mesh partitioning algorithm for preserving spatial locality in arbitrary geometries

Highlights: •An algorithm for partitioning computational meshes is proposed. •The Morton order space-filling curve is modified to achieve improved locality. •A spatial locality metric is defined to compare results with existing approaches. •Results indicate improved performance of the algorithm in complex geometries. -- Abstract: A space-filling curve (SFC) is a proximity preserving linear mapping of any multi-dimensional space and is widely used as a clustering tool. Equi-sized partitioning of an SFC ignores the loss in clustering quality that occurs due to inaccuracies in the mapping. Often, this results in poor locality within partitions, especially for the conceptually simple, Morton order curves. We present a heuristic that improves partition locality in arbitrary geometries by slicing a Morton order curve at points where spatial locality is sacrificed. In addition, we develop algorithms that evenly distribute points to the extent possible while maintaining spatial locality. A metric is defined to estimate relative inter-partition contact as an indicator of communication in parallel computing architectures. Domain partitioning tests have been conducted on geometries relevant to turbulent reactive flow simulations. The results obtained highlight the performance of our method as an unsupervised and computationally inexpensive domain partitioning tool.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Nonadiabatic electron response in the Hasegawa-Wakatani equations

Tokamak edge turbulence is strongly influenced by parallel electron physics, which relaxes density and potential fluctuations towards electron adiabatic response. Beginning with the paradigmatic Hasegawa-Wakatani equations (HWEs) for resistive tokamak edge turbulence, a unique decomposition of the electric potential (φ) into adiabatic (a) and nonadiabatic (b) portions is derived, based on the requirement that a neither drive nor respond to the parallel current j{sub ∥}. The form of the decomposition clarifies that, at perpendicular scales large relative to the sound radius, the electron adiabatic response controls the nonzonal φ, not the fluctuating density n. Simple energy balance arguments allow one to rigorously bound the ratio of rms nonzonal nonadiabatic fluctuations (b(tilde sign)) relative to adiabatic ones (ã). The role of the vorticity nonlinearity in transferring energy between adiabatic and nonadiabatic fluctuations aids intuitive understanding of self-sustained turbulence in the HWEs. When the normalized parallel resistivity is weak, b(tilde sign) becomes effectively slaved, allowing the reduction to an approximate one-field model that remains valid for strong turbulence. In addition to guiding physical intuition, the one-field reduction should greatly ease further analytical manipulations. Direct numerical simulation of the 2D HWEs confirms the convergence of the asymptotic formula for b(tilde sign)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

A comprehensive comparison for simulations of cosmic-ray muons underground

The two leading simulation frameworks used for the simulation of cosmic-ray muons underground are FLUKA and Geant4. There have been in the past various questions raised as to the equivalence of these codes regarding cosmogenically produced neutrons and radioactivity in an underground environment. Many experiments choose one of these frameworks, and because they typically have different geometries or locations, the issues relating to code comparison are compounded. We report on an effort to compare the results of each of these codes in simulations which have simple geometry that is consistent between the two codes. It is seen that in terms of integrated neutron flux and neturon capture statistics the codes agree well in a broad sense. There are, however, differences that will be subject of further study. Comparisons of the simulations to available data are considered and the difficulties of such comparisons are pointed out.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2014

Spiral-like multi-beam emission via transformation electromagnetics

Transformation electromagnetics offers an unconventional approach for the design of novel radiating devices. Here, we propose an electromagnetic metamaterial able to split an isotropic radiation into multiple directive beams. By applying transformations that modify distance and angles, we show how the multiple directive beams can be steered at will. We describe transformation of the metric space and the calculation of the material parameters. Different transformations are proposed for a possible physical realization through the use of engineered artificial metamaterials. Full wave simulations are performed to validate the proposed approach. The idea paves the way to interesting applications in various domains in microwave and optical regimes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Self-generation of dissipative solitons in magnonic quasicrystal active ring resonator

Self-generation of dissipative solitons in the magnonic quasicrystal (MQC) active ring resonator is studied theoretically and experimentally. The developed magnonic crystal has quasiperiodic Fibonacci type structure. Frequency selectivity of the MQC together with the parametric three-wave decay of magnetostatic surface spin wave (MSSW) leads to the dissipative soliton self-generation. The transfer matrix method is used to describe MQC transmission responses. Besides, the model of MQC active ring resonator is suggested. The model includes three coupled differential equations describing the parametric decay of MSSW and two differential equations of linear oscillators describing the frequency selectivity of MQC. Numerical simulation results of dissipative soliton self-generation are in a fair agreement with experimental data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗