Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computer 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 1,675 records · Page 93

Magnetic reconnection under anisotropic magnetohydrodynamic approximation

We study the formation of slow-mode shocks in collisionless magnetic reconnection by using one- and two-dimensional collisionless MHD codes based on the double adiabatic approximation and the Landau closure model. We bridge the gap between the Petschek-type MHD reconnection model accompanied by a pair of slow shocks and the observational evidence of the rare occasion of in-situ slow shock observations. Our results showed that once magnetic reconnection takes place, a firehose-sense (p{sub ∥}>p{sub ⊥}) pressure anisotropy arises in the downstream region, and the generated slow shocks are quite weak comparing with those in an isotropic MHD. In spite of the weakness of the shocks, however, the resultant reconnection rate is 10%–30% higher than that in an isotropic case. This result implies that the slow shock does not necessarily play an important role in the energy conversion in the reconnection system and is consistent with the satellite observation in the Earth's magnetosphere.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

AN ANGLE-DEPENDENT SYNCHROTRON SELF-COMPTON MODEL FOR RELATIVISTIC JET SOURCES

We report on the development of a numerical code to calculate the angle-dependent synchrotron + synchrotron self-Compton radiation from relativistic jet sources with partially ordered magnetic fields and anisotropic particle distributions. Using a multi-zone radiation transfer approach, we can simulate magnetic-field configurations ranging from perfectly ordered (unidirectional) to randomly oriented (tangled). We demonstrate that synchrotron self-Compton model fits to the spectral energy distributions (SEDs) of extragalactic jet sources may be possible with a wide range of magnetic-field values, depending on their orientation with respect to the jet axis and the observer. This is illustrated with the example of a spectral fit to the SED of Mrk 421 from multiwavelength observations in 2006, where acceptable fits are possible with magnetic-field values varying within a range of an order of magnitude for different degrees of B-field alignment and orientation.

79 ASTRONOMY AND ASTROPHYSICS↗

CONSTRAINTS ON THE SHAPE OF THE MILKY WAY DARK MATTER HALO FROM JEANS EQUATIONS APPLIED TO SLOAN DIGITAL SKY SURVEY DATA

We search for evidence of dark matter in the Milky Way by utilizing the stellar number density distribution and kinematics measured by the Sloan Digital Sky Survey (SDSS) to heliocentric distances exceeding {approx}10 kpc. We employ the cylindrically symmetric form of Jeans equations and focus on the morphology of the resulting acceleration maps, rather than the normalization of the total mass as done in previous, mostly local, studies. Jeans equations are first applied to a mock catalog based on a cosmologically derived N-body+SPH simulation, and the known acceleration (gradient of gravitational potential) is successfully recovered. The same simulation is also used to quantify the impact of dark matter on the total acceleration. We use Galfast, a code designed to quantitatively reproduce SDSS measurements and selection effects, to generate a synthetic stellar catalog. We apply Jeans equations to this catalog and produce two-dimensional maps of stellar acceleration. These maps reveal that in a Newtonian framework, the implied gravitational potential cannot be explained by visible matter alone. The acceleration experienced by stars at galactocentric distances of {approx}20 kpc is three times larger than what can be explained by purely visible matter. The application of an analytic method for estimating the dark matter halo axis ratio to SDSS data implies an oblate halo with q{sub DM} = 0.47 {+-} 0.14 within the same distance range. These techniques can be used to map the dark matter halo to much larger distances from the Galactic center using upcoming deep optical surveys, such as LSST.

79 ASTRONOMY AND ASTROPHYSICS↗

Calculating uncertainty on K-effective with MONK10 - 14583

Criticality safety assessments require a demonstration that a particular configuration of fissile material has an adequate sub-critical margin (k-effective sufficiently below unity) to ensure that the risk of criticality under normal operation and accident conditions is acceptable. The required sub-critical margin depends upon the uncertainty in the estimated value of k-effective. The uncertainty in the calculated value of k-effective arises from a number of sources, including: manufacturing tolerances on input data to the code (affecting geometry, compositions and densities), uncertainty in the nuclear data used by the code, stochastic uncertainty resulting from Monte Carlo simulation and modelling approximations/errors, including the inevitable bugs in the software. The ANSWERS Software Service, in collaboration with industrial partners, is developing a number of techniques to better understand and quantify uncertainty on predicted values of k-effective, using MONK. The SPRUCE utility code has been developed to allow uncertainty to be estimated using sampling methods. This can include the sampling of input parameters (including dimensions, compositions and densities) from statistical distributions. It can also include sampling different nuclear data libraries. A set of nuclear data libraries has been generated for this purpose by sampling from statistical distributions that represent the uncertainties in the published nuclear data evaluated files; a set of libraries has been produced for Latin Hypercube Sampling. By varying the input data and nuclear data, separate and combined uncertainties due to manufacturing tolerances and nuclear data can be derived. By performing least squares fitting on the results it is also possible to estimate the contribution of each of the uncertain inputs and a sensitivity method in MONK can break down the nuclear data uncertainty. (authors)

42 ENGINEERING↗

EFFICIENT THREE-DIMENSIONAL NLTE DUST RADIATIVE TRANSFER WITH SKIRT

We present an updated version of SKIRT, a three-dimensional (3D) Monte Carlo radiative transfer code developed to simulate dusty galaxies. The main novel characteristics of the SKIRT code are the use of a stellar foam to generate random positions, an efficient combination of eternal forced scattering and continuous absorption, and a new library approach that links the radiative transfer code to the DustEM dust emission library. This approach enables a fast, accurate, and self-consistent calculation of the dust emission of arbitrary mixtures of transiently heated dust grains and polycyclic aromatic hydrocarbons, even for full 3D models containing millions of dust cells. We have demonstrated the accuracy of the SKIRT code through a set of simulations based on the edge-on spiral galaxy UGC 4754. The models we ran were gradually refined from a smooth, two-dimensional, local thermal equilibrium (LTE) model to a fully 3D model that includes non-LTE (NLTE) dust emission and a clumpy structure of the dusty interstellar medium. We find that clumpy models absorb UV and optical radiation less efficiently than smooth models with the same amount of dust, and that the dust in clumpy models is on average both cooler and less luminous. Our simulations demonstrate that, given the appropriate use of optimization techniques, it is possible to efficiently and accurately run Monte Carlo radiative transfer simulations of arbitrary 3D structures of several million dust cells, including a full calculation of the NLTE emission by arbitrary dust mixtures.

79 ASTRONOMY AND ASTROPHYSICS↗

NORMAL TYPE Ia SUPERNOVAE FROM VIOLENT MERGERS OF WHITE DWARF BINARIES

One of the most important questions regarding the progenitor systems of Type Ia supernovae (SNe Ia) is whether mergers of two white dwarfs can lead to explosions that reproduce observations of normal events. Here we present a fully three-dimensional simulation of a violent merger of two carbon-oxygen white dwarfs with masses of 0.9 M{sub Sun} and 1.1 M{sub Sun} combining very high resolution and exact initial conditions. A well-tested combination of codes is used to study the system. We start with the dynamical inspiral phase and follow the subsequent thermonuclear explosion under the plausible assumption that a detonation forms in the process of merging. We then perform detailed nucleosynthesis calculations and radiative transfer simulations to predict synthetic observables from the homologously expanding supernova ejecta. We find that synthetic color light curves of our merger, which produces about 0.62 M{sub Sun} of {sup 56}Ni, show good agreement with those observed for normal SNe Ia in all wave bands from U to K. Line velocities in synthetic spectra around maximum light also agree well with observations. We conclude that violent mergers of massive white dwarfs can closely resemble normal SNe Ia. Therefore, depending on the number of such massive systems available these mergers may contribute at least a small fraction to the observed population of normal SNe Ia.

79 ASTRONOMY AND ASTROPHYSICS↗

Evaluation of fuel cycle scenarios on MOX fuel recycling in PWRs and SFRs

Prospects on advanced fuel cycle scenario are considered for achieving a progressive integration of Sodium Fast Reactor (SFR) technology within the current French Pressurized Water Reactor (PWR) nuclear fleet, in a view to benefit from fissile material multi-recycling capability. A step by step process is envisioned, and emphasis is put on its potential implementation through the nuclear mass inventory calculations with the COSAC code. The overall time scale is not optimized. The first step, already implemented in several countries, the plutonium coming from the reprocessing of used Light Water Reactor (LWR) fuels is recycled into a small number of LWRs. The second step is the progressive introduction of the first SFRs, in parallel with the continuation of step 1. This second step lets to prepare the optimized multi recycling of MOX fuel which is considered in step 3. Step 3 is characterized by the introduction of a greater number of SFR and MOX management between EPR reactors and SFRs. In the final step 4, all the fleet is formed with SFRs. This study assesses the viability of each step of the overall scenario. The switch from one step to the other one could result from different constrains related to issues such as resources, waste, experience feedback, public acceptance, country policy, etc.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

3D Visualisation of Un-georeferenced, Coordinate-poor Spatial Data in Complex D and D and WM Environments - 18294

This paper describes the development of an algorithm-based utility tool using Visual Basic that examines un-georeferenced or GPS-deprived spatial datasets (3D datasets lacking one or more of the key x, y and z components). Taking into account any known geometrical parameters which are input via a user-friendly front-end, the new code automatically evaluates attributes, including any spatial coordinate components, detects invalid survey data points and fills in bad or missing data points utilising data-specific interpolation (e.g. cubic or linear spline) of valid data points, such that local extrema are not created. Previously troublesome 2D data can be easily converted to realistic 3D contour display, integrated with other time-series data to create 4D movies, and incorporated into other existing 3D datasets. More importantly those processed and 'more complete' data can be easily stored in 3D or 4D spatial database management systems for future interrogation by using clear, realistic and intuitive 3D visualisation. Test cases are presented to demonstrate the value and versatility of the newly developed utility tool which include the 3D visualisation of: surveyed wall thicknesses of a helical coil, radiation mapping of a plant and aerial survey mapping of radiation in a region near Fukushima. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

ASPCAP: THE APOGEE STELLAR PARAMETER AND CHEMICAL ABUNDANCES PIPELINE

The Apache Point Observatory Galactic Evolution Experiment (APOGEE) has built the largest moderately high-resolution ( R ≈ 22,500) spectroscopic map of the stars across the Milky Way, and including dust-obscured areas. The APOGEE Stellar Parameter and Chemical Abundances Pipeline (ASPCAP) is the software developed for the automated analysis of these spectra. ASPCAP determines atmospheric parameters and chemical abundances from observed spectra by comparing observed spectra to libraries of theoretical spectra, using χ {sup 2} minimization in a multidimensional parameter space. The package consists of a fortran90 code that does the actual minimization and a wrapper IDL code for book-keeping and data handling. This paper explains in detail the ASPCAP components and functionality, and presents results from a number of tests designed to check its performance. ASPCAP provides stellar effective temperatures, surface gravities, and metallicities precise to 2%, 0.1 dex, and 0.05 dex, respectively, for most APOGEE stars, which are predominantly giants. It also provides abundances for up to 15 chemical elements with various levels of precision, typically under 0.1 dex. The final data release (DR12) of the Sloan Digital Sky Survey III contains an APOGEE database of more than 150,000 stars. ASPCAP development continues in the SDSS-IV APOGEE-2 survey.

79 ASTRONOMY AND ASTROPHYSICS↗

Design studies and commissioning plans for plasma acceleration research station experimental program

Plasma acceleration research station is an electron beam driven plasma wakefield acceleration test stand proposed for CLARA facility in Daresbury Laboratory. In this paper, the interaction between the electron beam and the plasma is numerically characterised via 2D numerical studies by using VSIM code. The wakefields induced by a single bunch travelling through the plasma were found to vary from 200 MV/m to 3 GV/m for a range of bunch length, bunch radius, and plasma densities. Energy gain for the particles populating the bunch tail through the wakefields driven by the head of the bunch was demonstrated. After determining the achievable field for various beams and plasma configurations, a reference setting was determined for further studies. Considering this reference setting, the beam quality studies were performed for a two-bunch acceleration case. The maximum energy gain as well as the energy spread mitigation by benefiting from the beam loading was investigated by positioning the witness and driver bunches with respect to each other. Emittance growth mechanisms were studied considering the beam-plasma and beam-wakefield interactions. Eventually, regarding the findings, the initial commissioning plans and the aims for the later stages were summarised.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Transition to CASMO-5M and SIMULATE-3K for stability analyses of the Swiss BWRs

For stability analyses of the Swiss operating Boiling-Water-Reactors (BWRs), the methodology employed and validated so far at the Paul Scherrer Inst. (PSI) was based on the RAMONA-3 code with a hybrid upstream static lattice/core analysis approach using CASMO-4 and PRESTO-2. More recently, steps were undertaken towards a new methodology based on the SIMULATE-3K (S3K) code for the dynamical analyses combined with the CMSYS system relying on the CASMO/SIMULATE-3 suite of codes and which was established at PSI to serve as framework for the development and validation of reference core models of all the Swiss reactors and operated cycles. This paper presents a first validation of the new methodology on the basis of a benchmark recently organised by a Swiss utility and including the participation of several international organisations with various codes/methods. Now in parallel, a transition from CASMO-4E (C4E) to CASMO-5M (C5M) as basis for the CMSYS core models was also recently initiated at PSI. Consequently, it was considered adequate to address the impact of this transition both for the steady-state core analyses as well as for the stability calculations and to achieve thereby, an integral approach for the validation of the new S3K methodology. Therefore, a comparative assessment of C4 versus C5M is also presented in this paper with particular emphasis on the void coefficients and their impact on the downstream stability analysis results. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Analysis of BWR OPRM plant data and detection algorithms with DSSPP

All U.S. BWRs are required to have licensed stability solutions that satisfy General Design Criteria (GDC) 10 and 12 of 10 CFR 50 Appendix A. Implemented solutions are either detect and suppress or preventive in nature. Detection and suppression of power oscillations is accomplished by specialized hardware and software such as the Oscillation Power Range Monitor (OPRM) utilized in Option III and Detect and Suppress Solution - Confirmation Density (DSS-CD) stability Long-Term Solutions (LTSs). The detection algorithms are designed to recognize a Thermal-Hydraulic Instability (THI) event and initiate control rod insertion before the power oscillations increase much higher above the noise level that may threaten the fuel integrity. Option III is the most widely used long-term stability solution in the US and has more than 200 reactor years of operational history. DSS-CD represents an evolutionary step from the stability LTS Option III and its licensed domain envelopes the Maximum Extended Load Line Limit Analysis Plus (MELLLA +) domain. In order to enhance the capability to investigate the sensitivity of key parameters of stability detection algorithms, GEH has developed a new engineering analysis code, namely DSSPP (Detect and Suppress Solution Post Processor), which is introduced in this paper. The DSSPP analysis tool represents a major advancement in the method for diagnosing the design of stability detection algorithms that enables designers to perform parametric studies of the key parameters relevant for THI events and to fine tune these system parameters such that a potential spurious scram might be avoided. Demonstrations of DSSPPs application are also presented in this paper utilizing actual plant THI data. A BWR/6 plant had a plant transient that included unplanned recirculation pump transfer from fast to slow speed resulting in about 100% to {approx}40% rated power decrease and about 99% to {approx}30% rated core flow decrease. As the feedwater temperature is reduced to equilibrium conditions, the power increased from about {approx}40% to about {approx}60% with little change inflow. A THI event developed and subsequently, an OPRM initiated scram occurred with Option III. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A GENERAL CIRCULATION MODEL FOR GASEOUS EXOPLANETS WITH DOUBLE-GRAY RADIATIVE TRANSFER

We present a new version of our code for modeling the atmospheric circulation on gaseous exoplanets, now employing a 'double-gray' radiative transfer scheme, which self-consistently solves for fluxes and heating throughout the atmosphere, including the emerging (observable) infrared flux. We separate the radiation into infrared and optical components, each with its own absorption coefficient, and solve standard two-stream radiative transfer equations. We use a constant optical absorption coefficient, while the infrared coefficient can scale as a power law with pressure; however, for simplicity, the results shown in this paper use a constant infrared coefficient. Here we describe our new code in detail and demonstrate its utility by presenting a generic hot Jupiter model. We discuss issues related to modeling the deepest pressures of the atmosphere and describe our use of the diffusion approximation for radiative fluxes at high optical depths. In addition, we present new models using a simple form for magnetic drag on the atmosphere. We calculate emitted thermal phase curves and find that our drag-free model has the brightest region of the atmosphere offset by {approx}12 Degree-Sign from the substellar point and a minimum flux that is 17% of the maximum, while the model with the strongest magnetic drag has an offset of only {approx}2 Degree-Sign and a ratio of 13%. Finally, we calculate rates of numerical loss of kinetic energy at {approx}15% for every model except for our strong-drag model, where there is no measurable loss; we speculate that this is due to the much decreased wind speeds in that model.

79 ASTRONOMY AND ASTROPHYSICS↗

Mechanistic source terms for the XE-100 reactor XSTERM Code Development Approach - 18587

X-Energy is developing the Xe-100 pebble bed modular High Temperature Gas Reactor (HTGR) as an advanced nuclear power source. To achieve its safety objectives, the Xe-100 reactor relies on materials inherent characteristics and passive safety features. Key features include the inherent high temperature characteristics of TRISO-coated fuel particles, graphite moderator, and helium coolant, along with passive heat removal capability of a low-power-density core. This paper will describe the approach adopted by X-Energy to develop and validate mechanistic source terms for the Xe-100. The Mechanistic Source Term (MST) approach for the Xe-100 reactor will build upon that of the prior DOE Next Generation Nuclear Program (NGNP) (Refs. 1, 2, 3 and 4). MST development is closely related to Fuel Qualification (FQ), Licensing Basis Event (LBE) selection and Probabilistic Risk Assessment (PRA). This paper provides a brief overview of the MST road map process. Software plays a crucial role in establishing and quantifying MSTs required for licensing the Xe-100. An insight into X-Energy s comprehensive, integrated software development program for calculating time-dependent radionuclide releases from the Xe-100 over a wide range of events is provided. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Simulation of radiation energy release in air showers

A simulation study of the energy released by extensive air showers in the form of MHz radiation is performed using the CoREAS simulation code. We develop an efficient method to extract this radiation energy from air-shower simulations. We determine the longitudinal profile of the radiation energy release and compare it to the longitudinal profile of the energy deposit by the electromagnetic component of the air shower. We find that the radiation energy corrected for the geometric dependence of the geomagnetic emission scales quadratically with the energy in the electromagnetic component of the air shower with a second-order dependence on the atmospheric density at the position of the maximum shower development X {sub max}. In a measurement where X {sub max} is not accessible, this second order dependence can be approximated using the zenith angle of the incoming direction of the air shower with only a minor loss in accuracy. Our method results in an intrinsic uncertainty of 4% in the determination of the energy in the electromagnetic air-shower component, which is well below current experimental uncertainties.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Formation of a strong electric field resulting in the excitation of microplasma discharges at the edge of a dielectric film on a metal in a plasma flow

Results are presented from experimental and analytical studies of the processes resulting in the excitation of microplasma discharges (MPDs) on a metal surface partially covered with a thin dielectric film under the action of an external plasma flow in vacuum. It is shown experimentally that MPDs are excited at the interface between the open metal surface and the region covered by the dielectric film. The probability of MPD excitation is investigated as a function of the thickness of the dielectric film deposited on the metal. It is found that, for a film thickness of 1 μm, the probability of MPD excitation is close to unity. As the film thickness decreases below ~10 nm or increases above ~10 μm, the probability of MPD excitation is reduced by more than two orders of magnitude. A two-dimensional kinetic numerical code is developed that allows one to model the processes of Debye sheath formation and generation of a strong electric field near the edge of a finite-thickness dielectric film on a metal surface in a plasma flow for different configurations of the film edge. It is shown that the maximum value of the tangential component of the electric field is reached at the film edge and amounts to E{sub max} ≈ |φ{sub 0}|/2d (where φ{sub 0} < 0 is the electric potential applied to the metal and d is the film thickness), which for typical conditions of experiments on the excitation of MPDs on metal surfaces (φ{sub 0} ≈–400 V, d ≈ 1 μm) yields E{sub max} ≈ 2 MV/cm. The results of kinetic simulations confirm the qualitative idea about the mechanism of the formation of a strong electric field resulting in the excitation of MPDs at the edge of a dielectric film on a metal surface in a plasma flow and agree with experimental data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Improvement of three-field based safety analysis code, SPACE, through verification and validation

To verify and validate the state-of-the-art safety analysis code, SPACE, features of three field governing equations were arranged and the comparative analysis with both analytic solution and experimental benchmarks were performed. With the world wide streamline on the safety analysis, SPACE uses two-fluid, three-field governing equations which aim to implement multi-dimensional and multi-scaled analysis. Through the Verification and Validation (V and V), developing code based on the three-field governing equation have not only been confirmed but also improved. V and V activities on SPACE have been conducted as to interfacial drag model in vertical turbulent flow. The verification of the interfacial drag model was based on cell size sensitivity test for confirming feasible calculation range of cell size. Furthermore, through the validation with analytic solution, the modification of the interfacial drag model in vertical turbulent annulus flow was carried out and the nominated model was incorporated. Modified interfacial drag model, Asali's correlation, yields improved result than previous interfacial drag model. The validation with experimental benchmarks had also been performed on the nearly horizontal counter current limitation. To confirm the possibility of the further application of the CCFL model in case with horizontal pipe, validation with nearly horizontal pipe experiment was carried out. V and V activities of SPACE with analytic solution and experimental benchmarks allow us to confirm/identify the tailorability for safety analysis. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Software development infrastructure for the HYBRID modeling and simulation project

One of the goals of the HYBRID modeling and simulation project is to assess the economic viability of hybrid systems in a market that contains renewable energy sources like wind. The idea is that it is possible for the nuclear plant to sell non-electric energy cushions, which absorb (at least partially) the volatility introduced by the renewable energy sources. This system is currently modeled in the Modelica programming language. To assess the economics of the system, an optimization procedure is trying to find the minimal cost of electricity production. The RAVEN code is used as a driver for the whole problem. It is assumed that at this stage, the HYBRID modeling and simulation framework can be classified as non-safety “research and development” software. The associated quality level is Quality Level 3 software. This imposes low requirements on quality control, testing and documentation. The quality level could change as the application development continues.Despite the low quality requirement level, a workflow for the HYBRID developers has been defined that include a coding standard and some documentation and testing requirements. The repository performs automated unit testing of contributed models. The automated testing is achieved via an open-source python script called BuildingsP from Lawrence Berkeley National Lab. BuildingsPy runs Modelica simulation tests using Dymola in an automated manner and generates and runs unit tests from Modelica scripts written by developers. In order to assure effective communication between the different national laboratories a biweekly videoconference has been set-up, where developers can report their progress and issues. In addition, periodic face-face meetings are organized intended to discuss high-level strategy decisions with management. A second means of communication is the developer email list. This is a list to which everybody can send emails that will be received by the collective of the developers and managers involved in the project. Thirdly, to exchange documents quickly, a SharePoint directory has been set-up. SharePoint allows teams and organizations to intelligently share, and collaborate on content from anywhere.

97 MATHEMATICS AND COMPUTING↗