Search NASA⌕ Search

SEARCH · Search NASA

Results for “analysis 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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

HYSPLIT/MACCS Atmospheric Dispersion Model Technical Documentation and Benchmark Analysis

The nuclear accident consequence analysis code MACCS has traditionally modeled dispersion during downwind transport using a Gaussian plume segment model. MACCS is designed to estimate consequence measures such as air concentrations and ground depositions, radiological doses, and health and economic impacts on a statistical basis over the course of a year to produce annualaveraged output measures. The objective of this work is to supplement the Gaussian atmospheric transport and diffusion (ATD) model currently in MACCS with a new option using the HYSPLIT model. HYSPLIT/MACCS coupling has been implemented, with HYSPLIT as an alternative ATD option. The subsequent calculations in MACCS use the HYSPLIT-generated air concentration, and ground deposition values to calculate the same range of output quantities (dose, health effects, risks, etc.) that can be generated when using the MACCS Gaussian ATD model. Based on the results from the verification test cases, the implementation of the HYSPLIT/MACCS coupling is confirmed. This report contains technical details of the HYSPLIT/MACCS coupling and presents a benchmark analysis using the HYSPLIT/MACCS coupling system. The benchmark analysis, which involves running specific scenarios and sensitivity studies designed to examine how the results generated by the traditional MACCS Gaussian plume segment model compare to the new, higher fidelity HYSPLIT/MACCS modeling option, demonstrates the modeling results that can be obtained by using this new option. The comparisons provided herein can also help decision-makers evaluate the potential benefit of using results based on higher fidelity modeling with the additional computational burden needed to perform the calculations. Three sensitivity studies to investigate the potential impact of alternative modeling options, regarding 1) input meteorological data set, 2) method to estimate stability class, and 3) plume dispersion model for larger distances, on consequence results were also performed. The results of these analyses are provided and discussed in this report.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI1996

Nuclear Safety [Vol. 37, No. 2, April-June 1996]

Nuclear Safety is a journal that covers significant issues in the field of nuclear safety. Its primary scope is safety in the design, construction, operation, and decommissioning of nuclear power reactors worldwide and the research and analysis activities that promote this goal, but it also encompasses the safety aspects of the entire nuclear fuel cycle, including fuel fabrication, spent-fuel processing and handling, and nuclear waste disposal, the handling of fissionable materials and radioisotopes, and the environmental effects of all these activities. Table of Contents for this issue follows. GENERAL SAFETY CONSIDERATIONS: 97 The Nuclear Community and the Public: Cognitive and Cultural Influences on Thinking About Nuclear Risk, M. A. Meyer; 109 Twenty-Third Water Reactor Safety Information Meeting, D. A. Copinger; ACCIDENT ANALYSIS: 126 Analysis of a PWR LBLOCA Without SCRAM, Trevor N. Tyler, Rafael Macian-Juan and John H. Mahaffy; DESIGN FEATURES: 139 Vulnerability of Multiple-Barrier Systems, N. C. Lind; ENVIRONMENTAL EFFECTS: 149 A Study of Wet Catalytic Oxidation of Radioactive Spent Ion Exchange Resin by Hydrogen Peroxide, Xingchao Jian, Tianbao Wu, and Guichun Yun; 157 A Comparison Study and Resolution of Differences Between Emergency Response and Safety Analysis Codes Used at the Savannah River Site, A. A. Simpkins; OPERATING EXPERIENCES: 164 Reactor Shutdown Experience, Compiled by J. W. Cletcher; RECENT DEVELOPMENTS: 167 Reports, Standards, and Safety Guides, D. S. Queener; 172 Proposed Rule Changes as of Dec. 31,1995; ANNOUNCEMENTS: 178 American Nuclear Society 1997 Annual Meeting; 178 American Nuclear Society Nuclear Criticality and Safety Division Topical Meeting; 176 The Authors.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI1996

Nuclear Safety [Vol. 37, No. 2, April-June 1996]

Nuclear Safety is a journal that covers significant issues in the field of nuclear safety. Its primary scope is safety in the design, construction, operation, and decommissioning of nuclear power reactors worldwide and the research and analysis activities that promote this goal, but it also encompasses the safety aspects of the entire nuclear fuel cycle, including fuel fabrication, spent-fuel processing and handling, and nuclear waste disposal, the handling of fissionable materials and radioisotopes, and the environmental effects of all these activities. Table of Contents for this issue follows. GENERAL SAFETY CONSIDERATIONS: 97 The Nuclear Community and the Public: Cognitive and Cultural Influences on Thinking About Nuclear Risk, M. A. Meyer; 109 Twenty-Third Water Reactor Safety Information Meeting, D. A. Copinger; ACCIDENT ANALYSIS: 126 Analysis of a PWR LBLOCA Without SCRAM, Trevor N. Tyler, Rafael Macian-Juan and John H. Mahaffy; DESIGN FEATURES: 139 Vulnerability of Multiple-Barrier Systems, N. C. Lind; ENVIRONMENTAL EFFECTS: 149 A Study of Wet Catalytic Oxidation of Radioactive Spent Ion Exchange Resin by Hydrogen Peroxide, Xingchao Jian, Tianbao Wu, and Guichun Yun; 157 A Comparison Study and Resolution of Differences Between Emergency Response and Safety Analysis Codes Used at the Savannah River Site, A. A. Simpkins; OPERATING EXPERIENCES: 164 Reactor Shutdown Experience, Compiled by J. W. Cletcher; RECENT DEVELOPMENTS: 167 Reports, Standards, and Safety Guides, D. S. Queener; 172 Proposed Rule Changes as of Dec. 31,1995; ANNOUNCEMENTS: 178 American Nuclear Society 1997 Annual Meeting; 178 American Nuclear Society Nuclear Criticality and Safety Division Topical Meeting; 176 The Authors.

05 NUCLEAR FUELS↗
DOE OSTI2015

Initial Probabilistic Evaluation of Reactor Pressure Vessel Fracture with Grizzly and Raven

The Grizzly code is being developed with the goal of creating a general tool that can be applied to study a variety of degradation mechanisms in nuclear power plant components. The first application of Grizzly has been to study fracture in embrittled reactor pressure vessels (RPVs). Grizzly can be used to model the thermal/mechanical response of an RPV under transient conditions that would be observed in a pressurized thermal shock (PTS) scenario. The global response of the vessel provides boundary conditions for local models of the material in the vicinity of a flaw. Fracture domain integrals are computed to obtain stress intensity factors, which can in turn be used to assess whether a fracture would initiate at a pre-existing flaw. These capabilities have been demonstrated previously. A typical RPV is likely to contain a large population of pre-existing flaws introduced during the manufacturing process. This flaw population is characterized stastistically through probability density functions of the flaw distributions. The use of probabilistic techniques is necessary to assess the likelihood of crack initiation during a transient event. This report documents initial work to perform probabilistic analysis of RPV fracture during a PTS event using a combination of the RAVEN risk analysis code and Grizzly. This work is limited in scope, considering only a single flaw with deterministic geometry, but with uncertainty introduced in the parameters that influence fracture toughness. These results are benchmarked against equivalent models run in the FAVOR code. When fully developed, the RAVEN/Grizzly methodology for modeling probabilistic fracture in RPVs will provide a general capability that can be used to consider a wider variety of vessel and flaw conditions that are difficult to consider with current tools. In addition, this will provide access to advanced probabilistic techniques provided by RAVEN, including adaptive sampling and parallelism, which can dramatically decrease run times.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

Development and Implementation of Mechanistic Terry Turbine Models in RELAP-7 to Simulate RCIC Normal Operation Conditions

As part of the efforts to understand the unexpected “self-regulating” mode of the RCIC (Reactor Core Isolation Cooling) systems in Fukushima accidents and extend BWR RCIC and PWR AFW (Auxiliary Feed Water) operational range and flexibility, mechanistic models for the Terry turbine, based on Sandia’s original work [1], have been developed and implemented in the RELAP-7 code to simulate the RCIC system. In 2016, our effort has been focused on normal working conditions of the RCIC system. More complex off-design conditions will be pursued in later years when more data are available. In the Sandia model, the turbine stator inlet velocity is provided according to a reduced-order model which was obtained from a large number of CFD (computational fluid dynamics) simulations. In this work, we propose an alternative method, using an under-expanded jet model to obtain the velocity and thermodynamic conditions for the turbine stator inlet. The models include both an adiabatic expansion process inside the nozzle and a free expansion process outside of the nozzle to ambient pressure. The combined models are able to predict the steam mass flow rate and supersonic velocity to the Terry turbine bucket entrance, which are the necessary input information for the Terry turbine rotor model. The analytical models for the nozzle were validated with experimental data and benchmarked with CFD simulations. The analytical models generally agree well with the experimental data and CFD simulations. The analytical models are suitable for implementation into a reactor system analysis code or severe accident code as part of mechanistic and dynamical models to understand the RCIC behaviors. The newly developed nozzle models and modified turbine rotor model according to the Sandia’s original work have been implemented into RELAP-7, along with the original Sandia Terry turbine model. A new pump model has also been developed and implemented to couple with the Terry turbine model. An input model was developed to test the Terry turbine RCIC system, which generates reasonable results. Both the INL RCIC model and the Sandia RCIC model produce results matching major rated parameters such as the rotational speed, pump torque, and the turbine shaft work for the normal operation condition. The Sandia model is more sensitive to the turbine outlet pressure than the INL model. The next step will be further refining the Terry turbine models by including two-phase flow cases so that off-design conditions can be simulated. The pump model could also be enhanced with the use of the homologous curves.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2012

Approach for Selection of Rayleigh Damping Parameters Used for Time History Analysis

Nonlinearities, whether geometric or material, need to be addressed in seismic analysis. One good analysis method that can address these nonlinearities is direct time integration with Rayleigh damping. Modal damping is the damping typically specified in seismic analysis Codes and Standards. Modal damping is constant for all frequencies where Rayleigh damping varies with frequency. An approach is proposed here for selection of Rayleigh damping coefficients to be used in seismic analyses that is consistent with given Modal damping. The approach uses the difference between the modal damping response and the Rayleigh damping response along with effective mass properties of the model being evaluated to match overall system response levels. This paper provides a simple example problem to demonstrate the approach. It also provides results for a finite element model representing an existing piping system. Displacement, acceleration, and stress results are compared from model runs using modal damping and model runs using Rayleigh damping with coefficients selected using the proposed method.

42 ENGINEERING↗
DOE OSTI2015

Information on the Advanced Plant Experiment (APEX) Test Facility

The purpose of this report provides information related to the design of the Oregon State University Advanced Plant Experiment (APEX) test facility. Information provided in this report have been pulled from the following information sources: Reference 1: R. Nourgaliev and et.al, "Summary Report on NGSAC (Next-Generation Safety Analysis Code) Development and Testing," Idaho National Laboratory, 2011. Note that this is report has not been released as an external report. Reference 2: O. Stevens, Characterization of the Advanced Plant Experiment (APEX) Passive Residual Heat Removal System Heat Exchanger, Master Thesis, June 1996. Reference 3: J. Reyes, Jr., Q. Wu, and J. King, Jr., Scaling Assessment for the Design of the OSU APEX-1000 Test Facility, OSU-APEX-03001 (Rev. 0), May 2003. Reference 4: J. Reyes et al, Final Report of the NRC AP600 Research Conducted at Oregon State University, NUREG/CR-6641, July 1999. Reference 5: K. Welter et al, APEX-1000 Confirmatory Testing to Support AP1000 Design Certification (non-proprietary), NUREG-1826, August 2005.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Simulation of the Fast Reactor Fuel Assembly Duct-Bowing Reactivity Effect Using Monte Carlo Neutron Transport and Finite Element Analysis

This paper discusses a new method of simulating the fuel assembly duct-bowing reactivity coefficient for EBR-II run 138B. Quantification of the fuel assembly duct-bowing reactivity effect in liquid metal–cooled fast reactors has been a persistent problem since they were first designed and operated. Simulation of the duct-bowing reactivity effect is difficult because the level of detail required to simulate the effect has exceeded most modeling capabilities. The new method outlined in this paper utilizes the finite element analysis code ANSYS to analyze the thermal and structural components. Here, the displacement of the fuel assembly duct due to thermal expansion and mechanical interaction was calculated by ANSYS using recorded EBR-II run 138B temperature and power boundary value data. The displacement values were incorporated into to a Monte Carlo model of EBR-II run 138B and keff was calculated. Multiple Monte Carlo calculations were performed with duct displacement values corresponding to different reactor temperatures. Using the calculated keff values associated with the different duct displacement results allowed calculation of the duct-bowing reactivity coefficient. The duct-bowing reactivity coefficient was calculated to be –14.5 × 10 –4 $/°C/ ± 4.4%.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2016

SINGLE PHASE ANALYTICAL MODELS FOR TERRY TURBINE NOZZLE

All BWR RCIC (Reactor Core Isolation Cooling) systems and PWR AFW (Auxiliary Feed Water) systems use Terry turbine, which is composed of the wheel with turbine buckets and several groups of fixed nozzles and reversing chambers inside the turbine casing. The inlet steam is accelerated through the turbine nozzle and impacts on the wheel buckets, generating work to drive the RCIC pump. As part of the efforts to understand the unexpected “self-regulating” mode of the RCIC systems in Fukushima accidents and extend BWR RCIC and PWR AFW operational range and flexibility, mechanistic models for the Terry turbine, based on Sandia National Laboratories’ original work, has been developed and implemented in the RELAP-7 code to simulate the RCIC system. RELAP-7 is a new reactor system code currently under development with the funding support from U.S. Department of Energy. The RELAP-7 code is a fully implicit code and the preconditioned Jacobian-free Newton-Krylov (JFNK) method is used to solve the discretized nonlinear system. This paper presents a set of analytical models for simulating the flow through the Terry turbine nozzles when inlet fluid is pure steam. The implementation of the models into RELAP-7 will be briefly discussed. In the Sandia model, the turbine bucket inlet velocity is provided according to a reduced-order model, which was obtained from a large number of CFD simulations. In this work, we propose an alternative method, using an under-expanded jet model to obtain the velocity and thermodynamic conditions for the turbine bucket inlet. The models include both adiabatic expansion process inside the nozzle and free expansion process out of the nozzle to reach the ambient pressure. The combined models are able to predict the steam mass flow rate and supersonic velocity to the Terry turbine bucket entrance, which are the necessary input conditions for the Terry Turbine rotor model. The nozzle analytical models were validated with experimental data and benchmarked with CFD simulations. The analytical models generally agree well with the experimental data and CFD simulations. The analytical models are suitable for implementation into a reactor system analysis code or severe accident code as part of mechanistic and dynamical models to understand the RCIC behaviors. The cases with two-phase flow at the turbine inlet will be pursued in future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2015

Design and performance testing of a novel far field Gamma system to assay radioactive waste in 400 and 200 litre drums - 15339

A Far Field Gamma ray Measurement system has been developed for measuring and assaying both 400 and 200 litre waste drums with a wide range of activity at nuclear power stations. These include radioactive waste drums with surface dose-rates of less than 2 mSv/h (200 mrem/h) and, at the higher range of activity, 400 litre drums with surface dose-rates of up to 50 mSv/h (5 rem/h). The assay system incorporates an electro-mechanically cooled high purity Germanium profile detector and advanced digital spectroscopic analysis electronics. A novel optimised conical-trapezoidal collimator has been incorporated to reduce the effect of the radioactive background in addition to the conventional detector lead shield with tin and copper graded lining to reduce the effect of lead X-rays. Two features have been incorporated to reduce detector dead time and maintain detector resolution when measuring high activity drums. The first is the inclusion of two automated tungsten filters of different thickness as an integral part of the detector collimator. The filters are used in conjunction with an automated rail system on which the detector platform is mounted. The detector can be positioned at different distances from the drum surface in order to reduce detector dead time. The deployment of the filters and the adjustment of the drum - detector position are automatically controlled based on a user defined dose-rate limit table. The data is provided by two Geiger-Muller dose-rate detectors, which are used to measure the drum surface dose-rate. This data is subsequently used to adjust the detector position and in cases of higher dose-rate deploy one or other of the tungsten filters. The drum rotation platform incorporates a load cell to determine drum weight and hence the drum density. This information is used to calculate an attenuation correction based on the assumption of uniform drum density and uniform distribution of activity within the drum. Where the waste drums have a more complicated and known regular internal structure, such as a small radioactive region surrounded by an annulus of shielding material, the analysis algorithms incorporated into the far field geometry and attenuation correction spectroscopic data analysis code are able to make a more accurate determination of drum activity than would be possible assuming a completely uniform drum. The detector energy and efficiency calibration is achieved using a point source with multiple gamma ray energy peaks. The response of the high purity Germanium detector is also modelled using the Monte Carlo Neutron Photon code and the model is benchmarked using the measured point source calibration data. The performance of the instrument has been determined for a range of uniform density matrix drums containing a set of Eu-152 line or rod sources located in re-entrant tubes positioned on an equal volume basis. When rotated, the test drums with volume distributed line sources simulate waste drums of different density each with a uniformly distributed radioactive source. The measurements have been validated using Monte Carlo simulations. (authors)

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

Implementation and validation of a fully implicit accumulator model in RELAP-7

This paper presents the implementation and validation of an accumulator model in RELAP-7 under the framework of preconditioned Jacobian free Newton Krylov (JFNK) method, based on the similar model used in RELAP5. RELAP-7 is a new nuclear reactor system safety analysis code being developed at the Idaho National Laboratory (INL). RELAP-7 is a fully implicit system code. The JFNK and preconditioning methods used in RELAP-7 is briefly discussed. The slightly modified accumulator model is summarized for completeness. The implemented model was validated with LOFT L3-1 test and benchmarked with RELAP5 results. RELAP-7 and RELAP5 had almost identical results for the accumulator gas pressure and water level, although there were some minor difference in other parameters such as accumulator gas temperature and tank wall temperature. One advantage of the JFNK method is its easiness to maintain and modify models due to fully separation of numerical methods from physical models. It would be straightforward to extend the current RELAP-7 accumulator model to simulate the advanced accumulator design. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

Thermo-Mechanical Analysis of SBLOCA with Loss of One Safety Injection System Using SCDAP/RELAP5

Design Extension Condition (DEC) for nuclear power plants has been being deliberated to bring it within a safety regulation boundary, following the 2011 Fukushima Daiichi accident. Moreover, a new regulation on the LOCA accident analysis designated as 10CFR50.45(c) has been proposed based on a concern of the non-conservatism of current criterion neglecting fuel burn-up effect on the cladding embrittlement. Therefore, a more comprehensive and realistic simulation of fuel rod behavior is required to cope with the stringent environment under LOCA accidents. A SCDAP/RELAP5 system analysis code has been used mainly on the prediction of severe accident phenomena. It has a set of mechanistic models for the fuel deformation, which plays an important role particularly in the early phase progression of the fuel melting. However, limited work has been done focusing on the models only to the fuel failure estimation through the DBA or candidates of DEC. Therefore, some analyses were performed here with a burnup dependent core model to understand the capability of SCDAP/RELAP5 to predict the thermo-mechanical behavior of the fuel cladding such as ballooning, burst and external/internal oxidation. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2012

PWR ENDF/B-VII cross-section libraries for ORIGEN-ARP

New pressurized water reactor (PWR) cross-section libraries were generated for use with the ORIGEN-ARP depletion sequence in the SCALE nuclear analysis code system. These libraries are based on ENDF/B-VII nuclear data and were generated using the two-dimensional depletion sequence, TRITON/NEWT, in SCALE 6.1. The libraries contain multiple burnup-dependent cross-sections for seven PWR fuel designs, with enrichments ranging from 1.5 to 6 wt% {sup 235}U. The burnup range has been extended from the 72 GWd/MTU used in previous versions of the libraries to 90 GWd/MTU. Validation of the libraries using radiochemical assay measurements and decay heat measurements for PWR spent fuel showed good agreement between calculated and experimental data. Verification against detailed TRITON simulations for the considered assembly designs showed that depletion calculations performed in ORIGEN-ARP with the pre-generated libraries provide similar results as obtained with direct TRITON depletion, while greatly reducing the computation time. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2012

AFR-100 safety analyses

The Advanced Fast Reactor-100 (AFR-100) is Argonne National Laboratory's 250 MWth metal-fueled modular sodium-cooled pool-type fast reactor concept. [1] A series of accident sequences that focused on the AFR-100's ability to provide protection against reactor damage during low probability accident sequences resulting from multiple equipment failures were examined. Protected and Unprotected Loss of Flow (PLOF and ULOF) and Unprotected Transient Over-Power (UTOP) accidents were simulated using the SAS4A/SASSYS-1 safety analysis code. The large heat capacity of the sodium in the pool-type reactor allows the AFR-100 to absorb large amounts of energy during a PLOF with relatively small temperature increases throughout the system. During a ULOF with a 25-second flow halving time, coolant and cladding temperatures peak around 720 deg. C within the first minute before reactivity feedback effects decrease power to match the flow. Core radial expansion and fuel Doppler provide the necessary feedback during the UTOP to bring the system back to critical before system temperatures exceed allowable limits. Simulation results indicate that adequate ULOF safety margins exist for the AFR-100 design with flow halving times of twenty-five seconds. Significant safety margins are maintained for PLOF accidents as well as UTOP accidents if a rod stop is used. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2012

Modeling the Structural Response from a Propagating High Explosive Using Smooth Particle Hydrodynamics

This report primarily concerns the use of two massively parallel finite element codes originally written and maintained at Lawrence Livermore National Laboratory. ALE3D is an explicit hydrodynamics code commonly employed to simulate wave propagation from high energy scenarios and the resulting interaction with nearby structures. This coupled response ensures that a structure is accurately applied with a blast loading varying both in space and time. Figure 1 illustrates the radial outward propagation of a pressure wave due to a center detonated spherical explosive originating from the lower left. The radial symmetry seen in this scenario is lost when instead a cylindrocal charge is detonated. Figure 2 indicates that a stronger, faster traveling pressure wave occurs in the direction of the normal axis to the cylinder. The ALE3D name is derived because of the use of arbitrary-Lagrange-Eulerian elements in which the mesh is allowed to advect; a process through which the mesh is modified to alleviate tanlging and general mesh distortion often cuased by high energy scenarios. The counterpart to an advecting element is a Lagrange element, whose mesh moves with the material. Ideally all structural components are kept Lagrange as long as possible to preserve accuracy of material variables and minimize advection related errors. Advection leads to mixed zoning, so using structural Lagrange elements also improves the visualization when post processing the results. A simplified representation of the advection process is shown in Figure 3. First the mesh is distorted due to material motion during the Lagrange step. The mesh is then shifted to an idealized and less distorted state to prevent irregular zones caused by the Lagrange motion. Lastly, the state variables are remapped to the elements of the newly constructed mesh. Note that Figure 3 represents a purely Eulerian mesh relaxation because the mesh is relocated back to the pre-Lagrange position. This is the case when the material flows through a still mesh. This is not typically done in an ALE3D analysis, especially if Lagrange elements exist. Deforming Lagrange elements would certainly tangle with a Eulerian mesh eventually. The best method in this case is to have an advecting mesh positioned as some relaxed version of the pre and post Lagrange step; this gives the best opportunity of modeling a high energy event with a combination of Lagrange and ALE elements. Dyne3D is another explicit dynamic analysis code, ParaDyn being the parallel version. ParaDyn is used for predicting the transient response of three dimensional structures using Lagrangian solid mechanics. Large deformation and mesh tangling is often resolved through the use of an element deletion scheme. This is useful to accommodate component failure, but if it is done purely as a means to preserve a useful mesh it can lead to problems because it does not maintain continuity of the material bulk response. Whatever medium exists between structural components is typically not modeled in ParaDyn. Instead, a structure either has a known loading profile applied or given initial conditions. The many included contact algorithms can calculate the loading response of materials if and when they collide. A recent implementation of an SPH module in which failed or deleted material nodes are converted to independent particles is currently being utilized for a variety of spall related problems and high velocity impact scenarios. Figure 4 shows an example of a projectile, given an initial velocity, and how it fails the first plate which generates SPH particles which then interact with and damage the second plate.

36 MATERIALS SCIENCE↗
DOE OSTI2012

On the prediction of horizontal bubbly flows using the interfacial area transport equation

To solve the two-fluid model utilized in current nuclear reactor system analysis codes, the interfacial area concentration (a i) is estimated through flow regime dependent correlations that rely on static regime transition criteria. This approach does not capture the continuous evolution of the interfacial structures, and thus, it can pose numerical issues near the transition boundaries. The interfacial area transport equation (IATE) can help address these shortcomings by providing a dynamic prediction of a a{sub i} through mechanistic source and sink terms that account for bubble coalescence and breakup. Most of the previous work for this approach has focused on vertical two-phase flow. However, relatively few studies have been performed for horizontal two-phase flows, where buoyancy strongly affects the phase distribution. To develop a one-dimensional, area-averaged form of the IATE for adiabatic, horizontal bubbly flows the following considerations are necessary: (1) pressure drop estimation, (2) bubble velocity/void fraction estimation, (3) determination of bubble interaction mechanisms, and (4) treatment of the asymmetric phase distribution. In the current work, treatment of the asymmetric phase distribution is presented. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2012

Semi-analytical prediction of hydraulic resistance and heat transfer for pipe and channel flows of water at supercritical pressure

Within the Generation IV International Forum the Supercritical Water Reactor is investigated. For its core design and safety analysis the efficient prediction of flow and heat transfer parameters such as the wall-shear stress and the heat-transfer coefficient for pipe and channel flows is needed. For circular pipe flows a numerical model based on the one-dimensional conservation equations of mass, momentum end energy in the radial direction is presented, referred to as a 'semi-analytical' method. An accurate, high-order numerical method is employed to evaluate previously derived analytical solutions of the governing equations. Flow turbulence is modeled using the algebraic approach of Prandtl/van-Karman, including a model for the buffer layer. The influence of wall roughness is taken into account by a new modified numerical damping function of the turbulence model. The thermo-hydraulic properties of water are implemented according to the international standard of 1997. This method has the potential to be used within a sub-channel analysis code and as wall-functions for CFD codes to predict the wall shear stress and the wall temperature. The present study presents a validation of the method with comparison of model results with experiments and multi-dimensional computational (CFD) studies in a wide range of flow parameters. The focus is laid on forced convection flows related to reactor design and near-design conditions. It is found, that the method can accurately predict the wall temperature even under deterioration conditions as they occur in the selected experiments (Yamagata el al. 1972 at 24.5 MPa, Ornatski et al. 1971 at 25.5 and Swenson et al. 1963 at 22.75 MPa). Comparison of the friction coefficient under high heat flux conditions including significant viscosity and density reductions near the wall with various correlations for the hydraulic resistance will be presented; the best agreement is achieve with the correlation of Pioro et al. 2004. It is concluded, that the semi-analytical prediction method is accurate and efficient for forced convection flows even at high heat fluxes. The wall-roughness height is identified as an important model parameter. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗