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,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
DOE OSTI2010

Critical issues for modeling dust transport in tokamaks

Dust produced in tokamaks is an important issue for fusion. Dust particles can introduce health and safety issues when in the same time can have an impact on reactor performance. Apart from the associated problems there are also potential benefits that make the better understanding of their behavior important. In this work the dust transport code Dust in TOKamakS will be used to explore the effect that variations in the plasma background and the physical model, describing the plasma-dust interaction, have on their predicted trajectories.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Evolutionary optimization of PAW data-sets for accurate high pressure simulations

We examine the challenge of performing accurate electronic structure calculations at high pressures by comparing the results of all-electron full potential linearized augmented-plane-wave calculations, as implemented in the WIEN2k code, with those of the projector augmented wave (PAW) method, as implemented in Quantum ESPRESSO or Abinit code. In particular, we focus on developing an automated and consistent way of generating transferable PAW data-sets that can closely produce the all electron equation of state defined from zero to arbitrary high pressures. The technique we propose is an evolutionary search procedure that exploits the ATOMPAW code to generate atomic data-sets and the Quantum ESPRESSO software suite for total energy calculations. We demonstrate different aspects of its workability by optimizing PAW basis functions of some elements relatively abundant in planetary interiors. In addition, we introduce a new measure of atomic data-set goodness by considering their performance uniformity over an extended pressure range. - Highlights: • An attempt made to achieve the AE-FLAPW accuracy, at the favorable computational efficiency of the PAW approach for an outspread pressure range. • Provided an automatic and consistent route for developing efficient and transferable PAW data-sets. • Explicated an improved measure for the goodness of data-sets by including uniformity of the data-set performance over an extended pressure range.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2011

THE DiskMass SURVEY. III. STELLAR KINEMATICS VIA CROSS-CORRELATION

We describe a new cross-correlation (CC) approach used by our survey to derive stellar kinematics from galaxy-continuum spectroscopy. This approach adopts the formal error analysis derived by Statler, but properly handles spectral masks. Thus, we address the primary concerns regarding application of the CC method to censored data, while maintaining its primary advantage by consolidating kinematic and template-mismatch information toward different regions of the CC function. We identify a systematic error in the nominal CC method of approximately 10% in velocity dispersion incurred by a mistreatment of detector-censored data, which is eliminated by our new method. We derive our approach from first principles, and we use Monte Carlo simulations to demonstrate its efficacy. An identical set of Monte Carlo simulations performed using the well-established penalized-pixel-fitting code of Cappellari and Emsellem compares favorably with the results from our newly implemented software. Finally, we provide a practical demonstration of this software by extracting stellar kinematics from SparsePak spectra of UGC 6918.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

H2OI95: A Stand-Alone Fortran Code for Evaluating the IAPWS-95 Equation-of-State Model for Water (Rev. 1)

H2OI95 is a stand-alone Fortran code for evaluating the IAPWS-95 equation-of-state model for the thermodynamic properties of water. It also evaluates the corresponding thermochemical properties of water consistent with the CODATA recommendations. The IAPWS-95 model is based on a model equation for the dimensionless Helmholtz energy for which the primary variables are the inverse reduced temperature (τ = T cr /T ) and reduced density (δ = ρ/ρ cr ). Here T is the absolute temperature (K), ρ is density (kg/m³), and the subscript “ cr ” refers to the critical point of water (647.096 K and 22.064 MPa pressure in this model, for which ρ cr is 322 kg/m³). The code solves four basic types of problems, distinguished by the specified inputs: 1. Temperature (K) and density (ρ kg/m³) or reduced density (δ) 2. Temperature (K) and pressure (MPa). 3. Temperature (K) on the saturation (liquid-vapor equilibrium) curve 4. Pressure (MPa) on the saturation curve. Each type of problem is run using a corresponding input (text) file. All but the first type of problem require iteration. For example, to solve for desired temperature and pressure, the reduced density must be adjusted to give the desired pressure. Iteration is accomplished using the Newton-Raphson method, though the secant method is also used in solving the fourth type of problem. For the last three types of problems, H2OI95 has been used to conduct numerical studies of convergence and the problem of multiple numerical solutions, some of which are not physically valid. Obtaining valid results depends mainly on appropriate choice of starting values for the density. The default values used in H2OI95 appear to consistently lead to generally desired results. With modification (not addressed here), H2OI95 can be used to support SUPCRT92 and similar codes that compute chemical thermodynamic properties of species and reactions over a wide range of temperature and pressure (273.16-1273K and 0-1000 MPa). Version 1.1 of this software calculates an extended set of water properties, including some derivable solely from IAPWS-95, such as the isobaric thermal expansivity and the isothermal compressibility, and still other properties based on supplementary models sanctioned by IAPWS, including the dynamic viscosity, thermal conductivity, surface tension, and static dielectric constant, and various other properties (notably Debye-Hückel parameters) derivable from the dielectric constant and its partial derivatives with respect to temperature and pressure. A coding error in Version 1.0 was identified while testing the implementation of the viscosity model. This error is described in Appendix C. It mainly affects calculation of “higher order” properties such as compressibility and heat capacity near the critical point.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2010

Simulation of image detectors in radiology for determination of scatter-to-primary ratios using Monte Carlo radiation transport code MCNP/MCNPX

Purpose: The purpose of this study was to compare and validate three methods to simulate radiographic image detectors with the Monte Carlo software MCNP/MCNPX in a time efficient way. Methods: The first detector model was the standard semideterministic radiography tally, which has been used in previous image simulation studies. Next to the radiography tally two alternative stochastic detector models were developed: A perfect energy integrating detector and a detector based on the energy absorbed in the detector material. Validation of three image detector models was performed by comparing calculated scatter-to-primary ratios (SPRs) with the published and experimentally acquired SPR values. Results: For mammographic applications, SPRs computed with the radiography tally were up to 44% larger than the published results, while the SPRs computed with the perfect energy integrating detectors and the blur-free absorbed energy detector model were, on the average, 0.3% (ranging from -3% to 3%) and 0.4% (ranging from -5% to 5%) lower, respectively. For general radiography applications, the radiography tally overestimated the measured SPR by as much as 46%. The SPRs calculated with the perfect energy integrating detectors were, on the average, 4.7% (ranging from -5.3% to -4%) lower than the measured SPRs, whereas for the blur-free absorbed energy detector model, the calculated SPRs were, on the average, 1.3% (ranging from -0.1% to 2.4%) larger than the measured SPRs. Conclusions: For mammographic applications, both the perfect energy integrating detector model and the blur-free energy absorbing detector model can be used to simulate image detectors, whereas for conventional x-ray imaging using higher energies, the blur-free energy absorbing detector model is the most appropriate image detector model. The radiography tally overestimates the scattered part and should therefore not be used to simulate radiographic image detectors.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2015

SU-E-T-02: 90Y Microspheres Dosimetry Calculation with Voxel-S-Value Method: A Simple Use in the Clinic

Purpose: To present a simple and feasible method of voxel-S-value (VSV) dosimetry calculation for daily clinical use in radioembolization (RE) with {sup 90}Y microspheres. Dose distributions are obtained and visualized over CT images. Methods: Spatial dose distributions and dose in liver and tumor are calculated for RE patients treated with Sirtex Medical miscrospheres at our center. Data obtained from the previous simulation of treatment were the basis for calculations: Tc-99m maggregated albumin SPECT-CT study in a gammacamera (Infinia, General Electric Healthcare.). Attenuation correction and ordered-subsets expectation maximization (OSEM) algorithm were applied.For VSV calculations, both SPECT and CT were exported from the gammacamera workstation and registered with the radiotherapy treatment planning system (Eclipse, Varian Medical systems). Convolution of activity matrix and local dose deposition kernel (S values) was implemented with an in-house developed software based on Python code. The kernel was downloaded from www.medphys.it. Final dose distribution was evaluated with the free software Dicompyler. Results: Liver mean dose is consistent with Partition method calculations (accepted as a good standard). Tumor dose has not been evaluated due to the high dependence on its contouring. Small lesion size, hot spots in health tissue and blurred limits can affect a lot the dose distribution in tumors. Extra work includes: export and import of images and other dicom files, create and calculate a dummy plan of external radiotherapy, convolution calculation and evaluation of the dose distribution with dicompyler. Total time spent is less than 2 hours. Conclusion: VSV calculations do not require any extra appointment or any uncomfortable process for patient. The total process is short enough to carry it out the same day of simulation and to contribute to prescription decisions prior to treatment. Three-dimensional dose knowledge provides much more information than other methods of dose calculation usually applied in the clinic.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

photon-plasma: A modern high-order particle-in-cell code

We present the photon-plasma code, a modern high order charge conserving particle-in-cell code for simulating relativistic plasmas. The code is using a high order implicit field solver and a novel high order charge conserving interpolation scheme for particle-to-cell interpolation and charge deposition. It includes powerful diagnostics tools with on-the-fly particle tracking, synthetic spectra integration, 2D volume slicing, and a new method to correctly account for radiative cooling in the simulations. A robust technique for imposing (time-dependent) particle and field fluxes on the boundaries is also presented. Using a hybrid OpenMP and MPI approach, the code scales efficiently from 8 to more than 250.000 cores with almost linear weak scaling on a range of architectures. The code is tested with the classical benchmarks particle heating, cold beam instability, and two-stream instability. We also present particle-in-cell simulations of the Kelvin-Helmholtz instability, and new results on radiative collisionless shocks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

Benchmark specifications for EBR-II shutdown heat removal tests

Argonne National Laboratory (ANL) is hosting an IAEA-coordinated research project on benchmark analyses of sodium-cooled fast reactor passive safety tests performed at the Experimental Breeder Reactor-II (EBR-II). The benchmark project involves analysis of a protected and an unprotected loss of flow tests conducted during an extensive testing program within the framework of the U.S. Integral Fast Reactor program to demonstrate the inherently safety features of EBR-II as a pool-type, sodium-cooled fast reactor prototype. The project is intended to improve the participants' design and safety analysis capabilities for sodium-cooled fast reactors through validation and qualification of safety analysis codes and methods. This paper provides a description of the EBR-II tests included in the program, and outlines the benchmark specifications being prepared to support the IAEA-coordinated research project. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2017

Developmental study of advanced HIS design method for digital I and C+HMIT of PWR plant - 277

A new conceptual frame of how to design and validate a digital HIS (human interface system) on an innovative numerical simulation basis is proposed for the support of plant operators' supervisory control of various types of automated complex NPPs (nuclear power plants). The proposed conceptual framework utilizes the object-oriented software for plant DiD (defense-in depth) risk monitor with the combination of nuclear reactor accident simulation by an advanced nuclear safety analysis code RELAP5/MOD4. The conceptual frame proposed in this paper will be applied for an example practice for the SBLOCA (small break loss of coolant accident) case of passive safety PWR (pressurized water reactor) AP1000. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2015

Preliminary critical experiments validation of a lattice physic code LATC - 14415

Core and System Integrated Engine for design and analysis (COSINE) is an integrated nuclear engineering code package, which is being developed by State Nuclear Power Software Development Center (SNPSDC) in China. A lattice physics code named LATC is included in COSINE.LATC is a 2D, steady-state code for modeling PWR fuel. A WIMS-D 69 group micro library based on ENDF/B-VI is used and there are two transportation methods in LATC, such as two steps method (collision probability method for cell homogenization, discrete coordinate method for assembly homogenization) and one step method (MOC). The UO{sub 2} fuel assembly, which may contain guide tubes, burnable poisons such as Integral Fuel Burnable Absorbers (IFBA) / Wet Annular Burnable Absorbers (WABA)/ Pyrexes or control rods etc, can be calculated in LATC. LATC provides the required cross-section data for the core physics simulator analysis code in COSINE. In order to evaluate the quality of the LATC, a strategy of validation is researched and applied to the LATC. Firstly, a validation requirement analysis is performed, and the corresponding critical experiments are collected, and then a critical experimental matrix considering the variety of fuel enrichments, fuel pitch and water to uranium ratio is generated. Finally, a comparison of results between LATC using collision probability and discrete coordinate method and critical experiments are made to make sure that the code works well. According to the validation strategy discussed above, the preliminary validation is carried out and critical experiments are calculated according to the critical experimental matrix and results from LATC code are compared with reference results. The result shows that the calculated results agree well with critical experiments, and do not indicate any significant trends with respect to the various conditions considered. (authors)

42 ENGINEERING↗
DOE OSTI2010

Target diagnostic control system implementation for the National Ignition Facility (invited)

The extreme physics of targets shocked by NIF's 192-beam laser is observed by a diverse suite of diagnostics. Many diagnostics are being developed by collaborators at other sites, but ad hoc controls could lead to unreliable and costly operations. A diagnostic control system (DCS) framework for both hardware and software facilitates development and eases integration. Each complex diagnostic typically uses an ensemble of electronic instruments attached to sensors, digitizers, cameras, and other devices. In the DCS architecture each instrument is interfaced to a low-cost WINDOWS XP processor and JAVA application. Each instrument is aggregated with others as needed in the supervisory system to form an integrated diagnostic. The JAVA framework provides data management, control services, and operator graphical user interface generation. DCS instruments are reusable by replication with reconfiguration for specific diagnostics in extensible markup language. Advantages include minimal application code, easy testing, and high reliability. Collaborators save costs by assembling diagnostics with existing DCS instruments. This talk discusses target diagnostic instrumentation used on NIF and presents the DCS architecture and framework.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2010

Nuclear Mass Visualization and Analysis at nuclearmasses.org

Nuclear masses play an important role in a wide variety of studies, motivating a tremendous surge in new precision mass measurements and global mass models. To handle this information, and to address inadequacies of existing text-based disseminations of nuclear masses, we have built a suite of codes--online at nuclearmasses.org--to manage, visualize, access, manipulate, share, compare, and analyze nuclear mass datasets. The services at this site are described, along with plans for future development.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2013

A variational multi-symplectic particle-in-cell algorithm with smoothing functions for the Vlasov-Maxwell system

Smoothing functions are commonly used to reduce numerical noise arising from coarse sampling of particles in particle-in-cell (PIC) plasma simulations. When applying smoothing functions to symplectic algorithms, the conservation of symplectic structure should be guaranteed to preserve good conservation properties. In this paper, we show how to construct a variational multi-symplectic PIC algorithm with smoothing functions for the Vlasov-Maxwell system. The conservation of the multi-symplectic structure and the reduction of numerical noise make this algorithm specifically suitable for simulating long-term dynamics of plasmas, such as those in the steady-state operation or long-pulse discharge of a super-conducting tokamak. The algorithm has been implemented in a 6D large scale PIC code. Numerical examples are given to demonstrate the good conservation properties of the multi-symplectic algorithm and the reduction of the noise due to the application of smoothing function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2011

Investigations of laser plasmas dynamics by means of real and virtual Langmuir probes

In this paper we propose a novel technique for LPP-Laser Produced Plasmas investigation, combining high time resolved measurements using compact Langmuir Probes with the output of a theoretical model called HYBLAS developed on purpose, which is able to simulate the charged particles collected by a so-called virtual probe. It will be shown that with an appropriate experimental set-up and with the use of a Matlab software able to accurately analyze the experimental I-V curves, laser plasmas can be investigated properly even if the probe is placed very close to the target surface. This permits not only to study the plume expansion with a high temporal resolution, but also to estimate correctly the self-generated coulomb electric field inside the plume and to detect the inner structure of the the first upcoming expanding plasma. HYBLAS is able to predict and describe the plume expansion at relatively low power densities and is a powerful method to compare directly the experimental current signals with the numerical results if the initial conditions are settled properly. A direct comparison of the theoretical data with the experimental ones realized on different metal targets shows that our method is able to predict properly the overall plasma expansion in the nanosecond laser pulse duration regime. The virtual probe method was moreover tested by comparing the numerical results with another numerical code called MULTI, which simulate the expansion by combining the hydrodynamics equations to a multigroup method in order to include the radiation transport. (authors)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2012

Presentation of FP7 matter project: general overview

The 2010-2012 implementation plan of the European Sustainable Nuclear Industrial Initiative (ESNII), prepared in the frame of the Sustainable Nuclear Energy Technology Platform (SNETP), establishes a very tight time schedule for the start of construction of the European Gen IV prototypes; namely the construction of the LFR ETPP (European Technology Pilot Plant) MYRRHA will start in 2014 and that of the SFR Prototype ASTRID will start in 2017. The GEN IV reactors pose new challenges to the designers and scientists in terms of higher operating temperature, higher fuel burn-up, and in some cases more corrosive environment with respect to the present technologies and which impacts the materials performance. In this frame, the MATTER (Materials Testing and Rules) Project starts well targeted R and D activities to perform careful materials studies in GEN IV operational conditions and to find out criteria for the correct use of these materials in relevant reactor applications. Aim of the MATTER Project (that involved 27 partners and will end in 2015) is to complement the materials researches, in the frame of the European Energy Research Alliance (EERA) guidelines, with the implementation of pre-normative rules. The MATTER Project is divided in 3 technical Domains (called DM): DM1 - Development of test and evaluation guidelines for structural materials: to develop/establish best practice guidelines for testing and evaluation procedures, which are aimed to screen and characterize nuclear materials for innovative nuclear systems. DM2 - Pre-normative R and D for Codes and Standards: Pre-normative activities are performed, comprehensive of experiments, to revise and update the design rules (with an EU level consensus) in order to answer to some short term needs of the two projects ASTRID and MYRRHA with respect to the design and the construction of structural components. DM3 - Joint Program Scheme, implementation and Priorities: to optimise the effectiveness and efficiency of the EERA Joint Program on nuclear materials for innovative reactors and to support specific research activities related to fundamental understanding of ODS steels fabrication. ODS steels are considered candidate materials, in the medium-long term, for high fuel burn-up cladding application. After a brief presentation of DM1 and DM3, this paper mainly focuses on description of Pre-normative R and D activities for Codes and Standards (DM2). (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2011

Measurement and simulation of the cross sections for the production of {sup 148}Gd in thin {sup nat}W and {sup 181}Ta targets irradiated with 0.4- to 2.6-GeV protons

The cross sections for the production of {sup 148}Gd in {sup nat}W and {sup 181}Ta targets irradiated by 0.4-, 0.6-, 0.8-, 1.2-, 1.6-, and 2.6-GeV protons at the ITEP accelerator complex have been measured by direct {alpha} spectrometry without chemical separation. The experimental data have been compared with the data obtained at other laboratories and with the theoretical simulations of the yields on the basis of the BERTINI, ISABEL, CEM03.02, INCL4.2, INCL4.5, CASCADE07, and PHITS codes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Cooperation on Improved Isotopic Identification and Analysis Software for Portable, Electrically Cooled High-Resolution Gamma Spectrometry Systems Final Report

Under a 2006 agreement between the Department of Energy (DOE) of the United States of America and the Institut de Radioprotection et de Sûreté Nucléaire (IRSN) of France, the National Nuclear Security Administration (NNSA) within DOE and IRSN initiated a collaboration to improve isotopic identification and analysis of nuclear material [i.e., plutonium (Pu) and uranium (U)]. The specific aim of the collaborative project was to develop new versions of two types of isotopic identification and analysis software: (1) the fixed-energy response-function analysis for multiple energies (FRAM) codes and (2) multi-group analysis (MGA) codes. The project is entitled Action Sheet 4 – Cooperation on Improved Isotopic Identification and Analysis Software for Portable, Electrically Cooled, High-Resolution Gamma Spectrometry Systems (Action Sheet 4). FRAM and MGA/U235HI are software codes used to analyze isotopic ratios of U and Pu. FRAM is an application that uses parameter sets for the analysis of U or Pu. MGA and U235HI are two separate applications that analyze Pu or U, respectively. They have traditionally been used by safeguards practitioners to analyze gamma spectra acquired with high-resolution gamma spectrometry (HRGS) systems that are cooled by liquid nitrogen. However, it was discovered that these analysis programs were not as accurate when used on spectra acquired with a newer generation of more portable, electrically cooled HRGS (ECHRGS) systems. In response to this need, DOE/NNSA and IRSN collaborated to update the FRAM and U235HI codes to improve their performance with newer ECHRGS systems. Lawrence Livermore National Laboratory (LLNL) and Los Alamos National Laboratory (LANL) performed this work for DOE/NNSA.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗