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,423 records · Page 79
All sourcesNASA NTRSDOE OSTI
DOE OSTI2012

CONVECTIVE BABCOCK-LEIGHTON DYNAMO MODELS

We present the first global, three-dimensional simulations of solar/stellar convection that take into account the influence of magnetic flux emergence by means of the Babcock-Leighton (BL) mechanism. We have shown that the inclusion of a BL poloidal source term in a convection simulation can promote cyclic activity in an otherwise steady dynamo. Some cycle properties are reminiscent of solar observations, such as the equatorward propagation of toroidal flux near the base of the convection zone. However, the cycle period in this young sun (rotating three times faster than the solar rate) is very short ({approx}6 months) and it is unclear whether much longer cycles may be achieved within this modeling framework, given the high efficiency of field generation and transport by the convection. Even so, the incorporation of mean-field parameterizations in three-dimensional convection simulations to account for elusive processes such as flux emergence may well prove useful in the future modeling of solar and stellar activity cycles.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

DISTRIBUTION OF MAGNETIC DISCONTINUITIES IN THE SOLAR WIND AND IN MAGNETOHYDRODYNAMIC TURBULENCE

The statistical properties of magnetic discontinuities in the solar wind are investigated by measuring fluctuations in the magnetic field direction, given by the rotation {Delta}{theta} that the magnetic field vector undergoes during time interval {Delta}t. We show that the probability density function, P({Delta}{theta}), can be described by a simple model in which the magnetic field vector purely rotates with a relative increment {Delta}B/B that is lognormally distributed. We find that the probability density function of increments, P({Delta}B/B), has a remarkable scaling property: the normalized variable x = ({Delta}B/B) {center_dot} ({Delta}t/{Delta}t{sub 0}){sup -{alpha}} has a universal lognormal distribution for all time intervals {Delta}t. We then compare measurements from the solar wind with those from direct numerical simulations of magnetohydrodynamic (MHD) turbulence. We find good agreement for P({Delta}{theta}) obtained in the two cases when the magnetic guide field to fluctuations ratio B{sub 0}/b{sub rms} is chosen accordingly. However, the scale invariance of P({Delta}B/B) is broken in the MHD simulations with relatively limited inertial interval, which causes P({Delta}{theta}) to scale with measurement interval differently than in the solar wind.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

GLOBAL NON-SPHERICAL OSCILLATIONS IN THREE-DIMENSIONAL 4π SIMULATIONS OF THE H-INGESTION FLASH

We performed three-dimensional simulations of proton-rich material entrainment into {sup 12}C-rich He-shell flash convection and the subsequent H-ingestion flash that took place in the post-asymptotic giant branch star Sakurai's object. Observations of the transient nature and anomalous abundance features are available to validate our method and assumptions, with the aim of applying them to very low-metallicity stars in the future. We include nuclear energy feedback from H burning and cover the full 4π geometry of the shell. Runs on 768{sup 3} and 1536{sup 3} grids agree well with each other and have been followed for 1500 minutes and 1200 minutes. After an 850 minute long quiescent entrainment phase, the simulations enter into a global non-spherical oscillation that is launched and sustained by individual ignition events of H-rich fluid pockets. Fast circumferential flows collide at the antipode and cause the formation and localized ignition of the next H-overabundant pocket. The cycle repeats for more than a dozen times while its amplitude decreases. During the global oscillation, the entrainment rate increases temporarily by a factor of ≈100. Entrained entropy quenches convective motions in the upper layer until the burning of entrained H establishes a separate convection zone. The lower-resolution run hints at the possibility that another global oscillation, perhaps even more violent, will follow. The location of the H-burning convection zone agrees with a one-dimensional model in which the mixing efficiency is calibrated to reproduce the light curve. The simulations have been performed at the NSF Blue Waters supercomputer at NCSA.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

Two-dimensional s-polarized solitary waves in relativistic plasmas. I. The fluid plasma model

The properties of two-dimensional linearly s-polarized solitary waves are investigated by fluid-Maxwell equations and particle-in-cell (PIC) simulations. These self-trapped electromagnetic waves appear during laser-plasma interactions, and they have a dominant electric field component E{sub z}, normal to the plane of the wave, that oscillates at a frequency below the electron plasma frequency {omega}{sub pe}. A set of equations that describe the waves are derived from the plasma fluid model in the case of cold or warm plasma and then solved numerically. The main features, including the maximum value of the vector potential amplitude, the total energy, the width, and the cavitation radius are presented as a function of the frequency. The amplitude of the vector potential increases monotonically as the frequency of the wave decreases, whereas the width reaches a minimum value at a frequency of the order of 0.82 {omega}{sub pe}. The results are compared with a set of PIC simulations where the solitary waves are excited by a high-intensity laser pulse.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Experimental evaluation and numerical modelling of timber-framed walls

Timber-framed wall buildings are seen all over Europe, especially in seismic regions, given its adequacy to resist earthquakes. The pombaline buildings, developed after the big 1755 earthquake that destroyed Lisbon, constitute one of the best examples of ancient seismic-resistant structures based on timber-framed masonry walls. But in addition to suffering natural degradation and frequent inadequate interventions, these buildings were not prepared to resist the severe seismic actions envisaged by today’s structural codes. In this context, they increasingly have to undergo rehabilitation and reinforcing works. The success of these projects strongly depends on the accuracy of the structural modelling, which in turn depends on knowledge of the behaviour of the existing structure. Aiming at providing reliable data to be used in structural design, a series of tests was carried out on the base element that constitutes the timber-framed masonry walls: the Saint Andrew’s cross. The experimental tests were numerically simulated in order to determine the parameters that lead to the most accurate results when using linear-elastic models.

42 ENGINEERING↗
DOE OSTI2012

The potential, limitations, and challenges of divide and conquer quantum electronic structure calculations on energetic materials

High explosives are an important class of energetic materials used in many weapons applications. Even with modern computers, the simulation of the dynamic chemical reactions and energy release is exceedingly challenging. While the scale of the detonation process may be macroscopic, the dynamic bond breaking responsible for the explosive release of energy is fundamentally quantum mechanical. Thus, any method that does not adequately describe bonding is destined to lack predictive capability on some level. Performing quantum mechanics calculations on systems with more than dozens of atoms is a gargantuan task, and severe approximation schemes must be employed in practical calculations. We have developed and tested a divide and conquer (DnC) scheme to obtain total energies, forces, and harmonic frequencies within semi-empirical quantum mechanics. The method is intended as an approximate but faster solution to the full problem and is possible due to the sparsity of the density matrix in many applications. The resulting total energy calculation scales linearly as the number of subsystems, and the method provides a path-forward to quantum mechanical simulations of millions of atoms.

36 MATERIALS SCIENCE↗
DOE OSTI2015

X-Ray Detector Simulations - Oral Presentation

The free-electron laser at LCLS produces X-Rays that are used in several facilities. This light source is so bright and quick that we are capable of producing movies of objects like proteins. But making these movies would not be possible without a device that can detect the X-Rays and produce images. We need X-Ray cameras. The challenges LCLS faces include the X-Rays’ high repetition rate of 120 Hz, short pulses that can reach 200 femto-seconds, and extreme peak brightness. We need detectors that are compatible with this light source, but before they can be used in the facilities, they must first be characterized. My project was to do just that, by making a computer simulation program. My presentation discusses the individual detectors I simulated, the details of my program, and how my project will help determine which detector is most useful for a specific experiment.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
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 OSTI2012

TRACE and TRAC-BF1 benchmark against Leibstadt plant data during the event inadvertent opening of relief valves

In framework of introducing TRACE code to transient analyses system codes for Leibstadt Power Plant (KKL), a conversion process of existing TRAC-BF1 model to TRACE has been started within KKL. In the first step, TRACE thermal-hydraulic model for KKL has been developed based on existing TRAC-BF1 model. In order to assess the code models a simulation of plant transient event is required. In this matter simulations of inadvertent opening of 8 relief valves event have been performed. The event occurs at KKL during normal operation, and it started when 8 relief valves open resulting in depressurization of the Reactor Pressure Vessel (RPV). The reactor was shutdown safely by SCRAM at low level. The high pressure core spray (HPCS) and the reactor core isolation cooling (RCIC) have been started manually in order to compensate the level drop. The remaining water in the feedwater (FW) lines flashes due to saturation conditions originated from RPV depressurization and refills the reactor downcomer. The plant boundary conditions have been used in the simulations and the FW flow rate has been adjusted for better prediction. The simulations reproduce the plant data with good agreement. It can be concluded that the TRAC-BF1 existing model has been used successfully to develop the TRACE model and the results of the calculations have shown good agreement with plant recorded data. Beside the modeling assessment, the TRACE and TRAC-BF1 capabilities to reproduce plant physical behavior during the transient have shown satisfactory results. The first step of developing KKL model for TRACE has been successfully achieved and this model is further developed in order to simulate more complex plant behavior such as Turbine Trip. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
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↗