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,585 records · Page 88

Calculation of the Phenix end-of-life test in natural circulation with the CATHARE code

The Inst. of Radiological Protection and Nuclear Safety (IRSN) acts as technical support to French public authorities. As such, IRSN is in charge of safety assessment of operating and under construction reactors, as well as future projects. In this framework, one current objective of IRSN is to evaluate the ability and accuracy of numerical tools to foresee consequences of accidents. One of the advantages pointed up for fast reactors cooled by heavy liquid metal is the possibility of decay heat removal based on natural convection. The promotion of this passive cooling mode in future safety demonstrations will involve the use of adapted and validated numerical codes. After the final shutdown of the Phenix sodium cooled fast reactor in 2009, a set of tests covering different areas was conducted for code validation, including a natural circulation test in the primary circuit. Experimental data were issued by CEA to organize a benchmark exercise in the frame of an IAEA Coordinated Research Project (CRP), with the objective to assess the system-codes capability in simulating the thermal-hydraulics behavior of sodium cooled fast reactors in such accidental conditions. IRSN participated to this benchmark with the CATHARE code. This code, co-developed by CEA, EDF, AREVA and IRSN and widely used for PWR safety studies, was recently extended for sodium applications. This paper presents the CATHARE modeling of the Phenix primary circuit and the results obtained. A relatively good agreement was found with available measurements considering the experimental uncertainties. This work stressed the local aspects of phenomena occurring during the natural convection establishment and the limits of a 0D/1D approach. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Sol-gel deposited transducers characteristics: A weight percent study - 15718

Monitoring functional machinery components that operate at elevated temperatures is costly and time consuming. This work explores leave in place composite transducers for high temperature monitoring of machinery components. A theoretical model was developed utilizing software provided by NASA. This model was then compared to the experimental results. The objective is to determine the characteristics of a set of spray-on composite transducers to identify the ideal transducer for any situation where a certain specification is made: - Operating temperature; - Signal amplitude/Signal to noise ratio; - Specific material property. The transducer is fabricated through a sol-gel deposition process that utilizes spraying the transducer onto the substrate. A sol-gel solution containing bismuth titanate is mixed with a combination of bismuth titanate and lead zirconate titanate(PZT) powders to create a composite transducer. Typical sintering methods have proved to be destructive and difficult to handle for the sol-gel deposition method; however, induction sintering has proved effective. -d{sub 33} decreases with increasing weight percent of bismuth titanate -the maximum usable temperature increases with increasing weight percent bismuth titanate -signal amplitude is also higher with increasing weight percent of bismuth titanate The ideal composition for the best signal amplitude and the highest operating temperature has been shown to be pure bismuth titanate; however, if the properties of the transducer are of utmost importance including conductivity and density then PZT should be added to increase these properties. The theoretical NASA Micromechanical Analysis Code by Generalized Method of Cell model has been validated through the experimental results of these experiments. Now that the composition of the transducer has been optimized to a specific application the future work is comprised of: - Utilizing the information to evaluate wall thickness - Integrating the optimized transducer with a wireless system for remote monitoring of mechanical component - Establishing a field deposition method for the spray-on transducer.

36 MATERIALS SCIENCE↗

Refactoring and Its Benefits

Refactoring is a disciplined technique for restructuring an existing body of code, altering its internal structure without changing its external behavior. It improves internal code structure without altering its external functionality by transforming functions and rethinking algorithms. It is an iterative process. Refactoring include reducing scope, replacing complex instructions with simpler or built-in instructions, and combining multiple statements into one statement. By transforming the code with refactoring techniques it will be faster to change, execute, and download. It is an excellent best practice to adopt for programmers wanting to improve their productivity. Refactoring is similar to things like performance optimizations, which are also behavior- preserving transformations. It also helps us find bugs when we are trying to fix a bug in difficult-to-understand code. By cleaning things up, we make it easier to expose the bug. Refactoring improves the quality of application design and implementation. In general, three cases concerning refactoring. Iterative refactoring, Refactoring when is necessary, Not refactor.Mr. Martin Fowler identifies four key reasons to refractor. Refactoring improves the design of software, makes software easier to understand, helps us find bugs and also helps in executing the program faster. There is an additional benefit of refactoring. It changes the way a developer thinks about the implementation when not refactoring. There are the three types of refactorings. 1) Code refactoring: It often referred to simply as refactoring. This is the refactoring of programming source code. 2) Database refactoring: It is a simple change to a database schema that improves its design while retaining both its behavioral and informational semantics. 3) User interface (UI) refactoring: It is a simple change to the UI which retains its semantics. Finally, we conclude the benefits of Refactoring are: Improves the design of software, Makes software easier to understand, Software gets cleaned up and Helps us to find bugs and Helps us to program faster.

99 GENERAL AND MISCELLANEOUS↗

Studies on the coupling transformer to improve the performance of microwave ion source

A 2.45 GHz microwave ion source has been developed and installed at the Variable Energy Cyclotron Centre to produce high intensity proton beam. It is operational and has already produced more than 12 mA of proton beam with just 350 W of microwave power. In order to optimize the coupling of microwave power to the plasma, a maximally flat matching transformer has been used. In this paper, we first describe an analytical method to design the matching transformer and then present the results of rigorous simulation performed using ANSYS HFSS code to understand the effect of different parameters on the transformed impedance and reflection and transmission coefficients. Based on the simulation results, we have chosen two different coupling transformers which are double ridged waveguides with ridge widths of 24 mm and 48 mm. We have fabricated these transformers and performed experiments to study the influence of these transformers on the coupling of microwave to plasma and extracted beam current from the ion source.

43 PARTICLE ACCELERATORS↗

Theoretical and numerical predictions of hypervelocity impact-generated plasma

The hypervelocity impact generated plasmas (HVIGP) in thermodynamic non-equilibrium state were theoretically analyzed, and a physical model was presented to explore the relationship between plasma ionization degree and internal energy of the system by a group of equations including a chemical reaction equilibrium equation, a chemical reaction rate equation, and an energy conservation equation. A series of AUTODYN 3D (a widely used software in dynamic numerical simulations and developed by Century Dynamic Inc.) numerical simulations of the impacts of hypervelocity Al projectile on its targets at different incident angles were performed. The internal energy and the material density obtained from the numerical simulations were then used to calculate the ionization degree and the electron temperature. Based on a self-developed 2D smooth particle hydrodynamic (SPH) code and the theoretical model, the plasmas generated by 6 hypervelocity impacts were directly simulated and their total charges were calculated. The numerical results are in good agreements with the experimental results as well as the empirical formulas, demonstrating that the theoretical model is justified by the AUTODYN 3D and self-developed 2D SPH simulations and applicable to predict HVIGPs. The study is of significance for astrophysical and cosmonautic researches and safety.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

High Speed Tribology of the Head-Disk Interface by Computer Simulation Final Report CRADA No. TC-297-92-C

The purpose of this project is to develop and validate a molecular dynamics model of nanoscale tribological processes occurring inside the head-disk interface of magnetic recording devices. Within the three years of the project the partners will: write a new molecular dynamics computer simulation code for growth, sliding and impact inside the head-disk interface; identify accurate models for chemical bonding within coatings used to protect magnetic surfaces, especially amorphous carbon thin films; simulate the growth of amorphous carbon thin films and analyse the resulting chemical bonding; simulate friction and wear processes relevant to the head-disk interface; and validate the model through comparison with AFM and nanoindenter studies of friction and wear of amorphous carbon thin films.

36 MATERIALS SCIENCE↗

Experimental validation of CASMO-4E and CASMO-5M for radial fission rate distributions in a westinghouse SVEA-96 Optima2 BWR fuel assembly

Measured and calculated radial total fission rate distributions are compared for the three axial sections of a Westinghouse SVEA-96 Optima2 BWR fuel assembly, comprising 96, 92 and 84 fuel rods, respectively. The measurements were performed on a full-size fuel assembly in the PROTEUS zero-power experimental facility. The measured fission rates are compared to the results of the CASMO-4E and CASMO-5M fuel assembly codes. Detailed measured geometrical data were used in the models, and effects of the surrounding zones of the reactor were taken into account by correction factors derived from MCNPX calculations. The results of the calculations agree well with those of the experiments, with root-mean-square deviations between 1.2% and 1.5% and maximum deviations of 3-4%. The quality of the predictions by CASMO-4E and CASMO-5M is comparable. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

CTF Parallel Performance Improvements

CTF is a modernized version of COBRA-TF (Coolant Boiling in Rod Arrays-Two Fluid), which is a thermal hydraulic (T/H) sub-channel code based on a two-phase, two-fluid model that formulates the conservation equations of mass, energy, and momentum for three fields of vapor, continuous liquid, and entrained liquid droplets. Under the Consortium for Advanced Simulation of Light Water Reactors (CASL) program sponsored by the U.S. Department of Energy, CTF is being further developed and improved for modeling and simulation of light water reactors. The code improvements for predicting T/H responses of pressurized water reactor (PWR) cores are focused on its applications to specific challenge problems deemed significant to the industry, such as departure from nucleate boiling (DNB) and crud-induced power shift. The code is also being used for T/H feedback in normal operating condition depletions in PWR and boiling water reactor designs in the reactor core simulator, Virtual Environment for Reactor Applications Core Simulator (VERA-CS), being developed by CASL. CTF modifications under the CASL program include software optimization, new closure models, and parallelization for modeling full reactor core T/H responses. The purpose of this work was to improve the parallel performance of CTF. CTF was originally parallelized via a domain-decomposition strategy using the Message Passing Interface (MPI), in which the basis of a solution domain was a physical assembly in the model. In other words, a user who wanted to perform a parallel solution with CTF had to provide the code one processor per assembly in the model (e.g., 56 cores for a 56 assembly quarter-core model, or 193 cores for a 193 assembly full-core model). This design was rigid, demanding a very specific number of processors for a model; furthermore, when CTF was run in VERA-CS-which typically uses thousands of processors- many processors sat idle during the CTF solve because CTF was unable to use any more processors than there are physical assemblies. This was not necessarily a problem if the CTF execution speed was very fast, but the average time spent in a CTF solve for prototypical PWR operating conditions was significant (i.e., over 20 % of the simulation wall time was being spent in CTF for a standard core depletion simulation). Larger, more computationally intensive simulations performed in support of the CASL DNB challenge problem demonstrated a need for improvements in this original parallelization approach. The specific goal of this work was to allow CTF to use more processors than there are physical assemblies in the model, which would not only speed up CTF and improve the overall efficiency of VERA-CS but also give users more-flexible options for executing their simulations. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Implementation of a Pressure Correction Term in the CTF Subchannel Code

CTF is a modernized and improved version of the legacy subchannel code COBRA-TF, being jointly developed and maintained by Oak Ridge National Laboratory (ORNL) and North Carolina State University. The code was adopted by ORNL for use in the Consortium for Advanced Simulation of Light Water Reactors (CASL) in 2012 to help address CASL challenge problems. Since that time, activities related to CTF have included implementing software quality assurance measures, implementing new features and models, performing validation and verification testing, establishing and supporting a CTF user group, and developing the code for use in coupled applications. CTF is applicable to single- and two-phase flows in light water reactor geometries at normal and accident operating conditions. The two-phase flow model consists of the six-equation two-phase flow model augmented by equations for the droplet field, and it also includes appropriate source terms. This work focuses on the numerical stability of the six-equation two-phase flow model and is part of a recent effort to improve CTF capabilities for simulating flows in boiling water reactors (BWRs). The six-equation two-phase flow model is not strongly hyperbolic, meaning that some of its eigenvalues are either not real or not distinct. In a Cauchy sense, and for an initial-value problem, strong hyperbolicity was found to be equivalent to well-posedness, which ensures the uniqueness of the solution in the space-time used. This concept is also of importance when designing high-order numerical methods and also when implementing boundary conditions that rely on real eigenvalues such as Riemann solvers. A tremendous amount of research is available in the literature on how to hyperbolize the six-equation two-phase flow model by adding ad hoc terms to the system of equations. Such terms are often referred to as virtual mass or interfacial pressure correction terms. From a mathematical perspective, this approach adds terms to the phasic momentum equation containing first-order partial derivatives of the unknowns to modify the Jacobian matrix, thus the imaginary part of the eigenvalues. For this work, an interfacial pressure correction (IPC) term on the same model as what was used for the CATHARE system code, is implemented in CTF. This allows for recovery of a two-pressure formulation. In the following, the six-equation two-phase flow model is recalled, along with the IPC term and its influence on the eigenvalues. Then, numerical results and convergence studies are presented for a water-faucet problem at a pressure of 10 bar. Finally, conclusions and future work are addressed.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A new code for the design and analysis of the heliostat field layout for power tower system

A new code for the design and analysis of the heliostat field layout for power tower system is developed. In the new code, a new method for the heliostat field layout is proposed based on the edge ray principle of nonimaging optics. The heliostat field boundary is constrained by the tower height, the receiver tilt angle and size and the heliostat efficiency factor which is the product of the annual cosine efficiency and the annual atmospheric transmission efficiency. With the new method, the heliostat can be placed with a higher efficiency and a faster response speed of the design and optimization can be obtained. A new module for the analysis of the aspherical heliostat is created in the new code. A new toroidal heliostat field is designed and analyzed by using the new code. Compared with the spherical heliostat, the solar image radius of the field is reduced by about 30% by using the toroidal heliostat if the mirror shape and the tracking are ideal. In addition, to maximize the utilization of land, suitable crops can be considered to be planted under heliostats. To evaluate the feasibility of the crop growth, a method for calculating the annual distribution of sunshine duration on the land surface is developed as well. (author)

14 SOLAR ENERGY↗

SU-E-T-157: CARMEN: A MatLab-Based Research Platform for Monte Carlo Treatment Planning (MCTP) and Customized System for Planning Evaluation

Purpose: Although there exist several radiotherapy research platforms, such as: CERR, the most widely used and referenced; SlicerRT, which allows treatment plan comparison from various sources; and MMCTP, a full MCTP system; it is still needed a full MCTP toolset that provides users complete control of calculation grids, interpolation methods and filters in order to “fairly” compare results from different TPSs, supporting verification with experimental measurements. Methods: This work presents CARMEN, a MatLab-based platform including multicore and GPGPU accelerated functions for loading RT data; designing treatment plans; and evaluating dose matrices and experimental data.CARMEN supports anatomic and functional imaging in DICOM format, as well as RTSTRUCT, RTPLAN and RTDOSE. Besides, it contains numerous tools to accomplish the MCTP process, managing egs4phant and phase space files.CARMEN planning mode assist in designing IMRT, VMAT and MERT treatments via both inverse and direct optimization. The evaluation mode contains a comprehensive toolset (e.g. 2D/3D gamma evaluation, difference matrices, profiles, DVH, etc.) to compare datasets from commercial TPS, MC simulations (i.e. 3ddose) and radiochromic film in a user-controlled manner. Results: CARMEN has been validated against commercial RTPs and well-established evaluation tools, showing coherent behavior of its multiple algorithms. Furthermore, CARMEN platform has been used to generate competitive complex treatment that has been published in comparative studies. Conclusion: A new research oriented MCTP platform with a customized validation toolset has been presented. Despite of being coded with a high-level programming language, CARMEN is agile due to the use of parallel algorithms. The wide-spread use of MatLab provides straightforward access to CARMEN’s algorithms to most researchers. Similarly, our platform can benefit from the MatLab community scientific developments as filters, registration algorithms etc. Finally, CARMEN arises the importance of grid and filtering control in treatment plan comparison.

60 APPLIED LIFE SCIENCES↗

Electronic, Vibrational, and Elastic Properties of Zn{sub 2}SeTe Crystals with Antihalcopirite Structure

Using the capabilities of the CRYSTAL14 code, the parameters of an equilibrium crystal lattice are obtained and the electronic and vibrational structures of a new chalcopyrite-like Zn{sub 2}SeTe crystal are studied. The energy band structure, long-wave vibrational frequencies at the point Г, and elastic moduli are calculated. The dependence of the contributions from the vibrations of individual atoms to the vibrational modes of a crystal caused by filling the cationic or anionic sublattices in the chalcopyrite structure is established.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Laser scanner data processing and 3D modeling using a free and open source software

The laser scanning is a technology that allows in a short time to run the relief geometric objects with a high level of detail and completeness, based on the signal emitted by the laser and the corresponding return signal. When the incident laser radiation hits the object to detect, then the radiation is reflected. The purpose is to build a three-dimensional digital model that allows to reconstruct the reality of the object and to conduct studies regarding the design, restoration and/or conservation. When the laser scanner is equipped with a digital camera, the result of the measurement process is a set of points in XYZ coordinates showing a high density and accuracy with radiometric and RGB tones. In this case, the set of measured points is called “point cloud” and allows the reconstruction of the Digital Surface Model. Even the post-processing is usually performed by closed source software, which is characterized by Copyright restricting the free use, free and open source software can increase the performance by far. Indeed, this latter can be freely used providing the possibility to display and even custom the source code. The experience started at the Faculty of Engineering in Catania is aimed at finding a valuable free and open source tool, MeshLab (Italian Software for data processing), to be compared with a reference closed source software for data processing, i.e. RapidForm. In this work, we compare the results obtained with MeshLab and Rapidform through the planning of the survey and the acquisition of the point cloud of a morphologically complex statue.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Analysis of direct-drive capsule compression experiments on the Iskra-5 laser facility

We have analyzed and numerically simulated our experiments on the compression of DT-gas-filled glass capsules under irradiation by a small number of beams on the Iskra-5 facility (12 beams) at the second harmonic of an iodine laser ({lambda} = 0.66 {mu}m) for a laser pulse energy of 2 kJ and duration of 0.5 ns in the case of asymmetric irradiation and compression. Our simulations include the construction of a target illumination map and a histogram of the target surface illumination distribution; 1D capsule compression simulations based on the DIANA code corresponding to various target surface regions; and 2D compression simulations based on the NUTCY code corresponding to the illumination conditions. We have succeeded in reproducing the shape of the compressed region at the time of maximum compression and the reduction in neutron yield (compared to the 1D simulations) to the experimentally observed values. For the Iskra-5 conditions, we have considered targets that can provide a more symmetric compression and a higher neutron yield.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

LACEwING: A New Moving Group Analysis Code

We present a new nearby young moving group (NYMG) kinematic membership analysis code, LocAting Constituent mEmbers In Nearby Groups (LACEwING), a new Catalog of Suspected Nearby Young Stars, a new list of bona fide members of moving groups, and a kinematic traceback code. LACEwING is a convergence-style algorithm with carefully vetted membership statistics based on a large numerical simulation of the Solar Neighborhood. Given spatial and kinematic information on stars, LACEwING calculates membership probabilities in 13 NYMGs and three open clusters within 100 pc. In addition to describing the inputs, methods, and products of the code, we provide comparisons of LACEwING to other popular kinematic moving group membership identification codes. As a proof of concept, we use LACEwING to reconsider the membership of 930 stellar systems in the Solar Neighborhood (within 100 pc) that have reported measurable lithium equivalent widths. We quantify the evidence in support of a population of young stars not attached to any NYMGs, which is a possible sign of new as-yet-undiscovered groups or of a field population of young stars.

79 ASTRONOMY AND ASTROPHYSICS↗

Numerical simulation of reactivity measurements in WER-1000 reactor

Reactivity is one of the most used and important concepts in physics and nuclear reactor calculations carried out for the safety analysis. Currently, design calculation of Russian WER reactors are carried out with modern coupled time-dependent neutronic and heat hydraulic codes. They allow to perform numerical simulation of reactivity measurements. However, point kinetic model used for simulation of large reactivity insertion leads to some issues. The paper discusses the numerical simulation of reactivity measurement, shows that the reactivity obtained from the steady state solution does not always correspond to the measured value. Comparison of scram system reactivity worth calculated and measured during physical start-up of unit 3, Kalinin NPP is presented. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Verification of the shift Monte Carlo code with the C5G7 reactor benchmark

Shift is a new hybrid Monte Carlo/deterministic radiation transport code being developed at Oak Ridge National Laboratory. At its current stage of development, Shift includes a parallel Monte Carlo capability for simulating eigenvalue and fixed-source multigroup transport problems. This paper focuses on recent efforts to verify Shift's Monte Carlo component using the two-dimensional and three-dimensional C5G7 NEA benchmark problems. Comparisons were made between the benchmark eigenvalues and those output by the Shift code. In addition, mesh-based scalar flux tally results generated by Shift were compared to those obtained using MCNP5 on an identical model and tally grid. The Shift-generated eigenvalues were within three standard deviations of the benchmark and MCNP5-1.60 values in all cases. The flux tallies generated by Shift were found to be in very good agreement with those from MCNP. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Silo & HDF5 I/O Scaling Improvements on BG/P Systems

Silo and HDF5 are I/O libraries used by many codes important to the LLNL's Weapons and Complex Integration (WCI) mission. In the past year, modest adjustments and tuning of Silo, HDF5 and the I/O configuration of the BG/P platform, Dawn, were undertaken. A key goal of this work was to improve I/O performance without requiring any changes in the application codes themselves. In particular, the application codes have been allowed to continue to use a simplified yet highly flexible I/O paradigm known as 'Poor Man's Parallel I/O', where scalability is achieved through concurrent, serial I/O to multiple files. The results demonstrate substantial performance gains (better than 50x in many cases) at large scale (greater than 64,000 MPI tasks). They describe key enhancements made to Silo, HDF5 and the I/O configuration of our BG/P platform and present very favorable results from scalability studies over a wide range of operating scenarios.

97 MATHEMATICS AND COMPUTING↗