Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reference Model 2”

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 847 records · Page 47

In-situ Condition Monitoring of Components in Small Modular Reactors Using Process and Electrical Signature Analysis. Final report, volume 1. Development of experimental flow control loop, data analysis and plant monitoring

The research and development under this project was focused on the following three major objectives: Objective 1: Identification of critical in-vessel SMR components for remote monitoring and development of their low-order dynamic models, along with a simulation model of an integral pressurized water reactor (iPWR). Objective 2: Development of an experimental flow control loop with motor-driven valves and pumps, incorporating data acquisition and on-line monitoring interface. Objective 3: Development of stationary and transient signal processing methods for electrical signatures, machinery vibration, and for characterizing process variables for equipment monitoring. This objective includes the development of a data analysis toolbox. The following is a summary of the technical accomplishments under this project: - A detailed literature review of various SMR types and electrical signature analysis of motor-driven systems was completed. A bibliography of literature is provided at the end of this report. Assistance was provided by ORNL in identifying some key references. - A review of literature on pump-motor modeling and digital signal processing methods was performed. - An existing flow control loop was upgraded with new instrumentation, data acquisition hardware and software. The upgrading of the experimental loop included the installation of a new submersible pump driven by a three-phase induction motor. All the sensors were calibrated before full-scale experimental runs were performed. - MATLAB-Simulink model of a three-phase induction motor and pump system was completed. The model was used to simulate normal operation and fault conditions in the motor-pump system, and to identify changes in the electrical signatures. - A simulation model of an integral PWR (iPWR) was updated and the MATLAB-Simulink model was validated for known transients. The pump-motor model was interfaced with the iPWR model for testing the impact of primary flow perturbations (upsets) on plant parameters and the pump electrical signatures. Additionally, the reactor simulation is being used to generate normal operation data and data with instrumentation faults and process anomalies. A frequency controller was interfaced with the motor power supply in order to vary the electrical supply frequency. The experimental flow control loop was used to generate operational data under varying motor performance characteristics. Coolant leakage events were simulated by varying the bypass loop flow rate. The accuracy of motor power calculation was improved by incorporating the power factor, computed from motor current and voltage in each phase of the induction motor.- A variety of experimental runs were made for steady-state and transient pump operating conditions. Process, vibration, and electrical signatures were measured using a submersible pump with variable supply frequency. High correlation was seen between motor current and pump discharge pressure signal; similar high correlation was exhibited between pump motor power and flow rate. Wide-band analysis indicated high coherence (in the frequency domain) between motor current and vibration signals. - Wide-band operational data from a PWR were acquired from AMS Corporation and used to develop time-series models, and to estimate signal spectrum and sensor time constant. All the data were from different pressure transmitters in the system, including primary and secondary loops. These signals were pre-processed using the wavelet transform for filtering both low-frequency and high-frequency bands. This technique of signal pre-processing provides minimum distortion of the data, and results in a more optimal estimation of time constants of plant sensors using time-series modeling techniques.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Comparative Study of foF2 Measurements with IRI-2007 Model Predictions During Extended Solar Minimum

The unusually deep and extended solar minimum of cycle 2324 made it very difficult to predict the solar indices 1 or 2 years into the future. Most of the predictions were proven wrong by the actual observed indices. IRI gets its solar, magnetic, and ionospheric indices from an indices file that is updated twice a year. In recent years, due to the unusual solar minimum, predictions had to be corrected downward with every new indices update. In this paper we analyse how much the uncertainties in the predictability of solar activity indices affect the IRI outcome and how the IRI values calculated with predicted and observed indices compared to the actual measurements.Monthly median values of F2 layer critical frequency (foF2) derived from the ionosonde measurements at the mid-latitude ionospheric station Juliusruh were compared with the International Reference Ionosphere (IRI-2007) model predictions. The analysis found that IRIprovides reliable results that compare well with actual measurements, when the definite (observed and adjusted) indices of solar activityare used, while IRI values based on earlier predictions of these indices noticeably overestimated the measurements during the solar minimum.One of the principal objectives of this paper is to direct attention of IRI users to update their solar activity indices files regularly.Use of an older index file can lead to serious IRI overestimations of F-region electron density during the recent extended solar minimum.

Ionosphere↗

Panel discussion: Challenges and improvements in accident dose analysis; Regulatory and industry perspective

Full text of publication follows. Design basis accident (DBA) analyses are one of several ways nuclear power plant licensees demonstrate, and the U.S. Nuclear Regulatory Commission (NRC) staff verifies, that the design and location of a nuclear power plant are appropriate and that the protection of health and safety of the public is maintained. In July 2000, the NRC staff issued Regulatory Guide 1.183, 'Alternative Radiological Source Terms for Evaluating Design Basis Accidents at Nuclear Power Reactors' (RG 1.183) [1]. RG 1.183 provides a method acceptable to the NRC staff for complying with the regulatory requirements contained in Title 10 of the Code of Federal Regulation (10 CFR) Section 50.67 (50.67), 'Accident source term' [2]. RG 1.183 also provides assumptions and parameters used to model postulated DBA radiological analyses for light-water reactors (LWR). Since the initial issuance of RG 1.183, the NRC staff and the nuclear power plant licensees have both gained substantial experience with the implementation of 10 CFR 50.67 and RG 1.183. Based on this experience, comments from licensees, the anticipation of licensing advanced LWRs, and new research, the NRC is developing a proposed revision to RG 1.183. This panel discussion will provide a discussion of that proposed revision, its potential impact on existing and future evaluations, and the industry perspective on the revised guidance. References: 1) USNRC, Regulatory Guide 1.183, 'Alternative Radiological Source Terms for Evaluating Design Basis Accidents at Nuclear Power Reactors,' Revision 0, U.S. Nuclear Regulatory Commission, Washington, D.C. (2000), Electronic copies are also available in ADAMS (http://www.nrc.gov/reading-rm/adams.html), under Accession No. ML003716792. 2) U.S. Code of Federal Regulations, Title 10 (Energy), Part 50 (Domestic Licensing of Production and Utilization Facilities), Section 67, (Accident source term), Office of the Federal Register, Washington, DC (1999). (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Anti-obesity potential of Glycyrrhiza uralensis and licochalcone A through induction of adipocyte browning

Highlights: • G.uralensis or licochalcone A reduced body weight and inguinal fat weights in obese mice. • G.uralensis or licochalcone A increased brown fat population in inguinal white adipose tissue of obese mice. • G.uralensis and LicoA reduce obesity and recover metabolic homeostasis by inducing fat browning. The main function of brown adipose tissue is to dissipate surplus caloric intake into heat energy by thermogenesis, increasing energy expenditure. Inducible brown adipocytes can develop within white adipose tissue (WAT) through a process referred to as browning. Browning of white fat represents a promising strategy for treatment of obesity and the related complications. We investigated whether Glycyrrhiza uralensis and its ingredients modulated adipogenesis through adipocyte browning using 3T3-L1 adipocytes and a high-fat diet (HFD)-induced obesity mice model. Amongst extracts and fractions of G. uralensis, methyl dichloride (MeCl{sub 2}) fraction was the most effective to induce expression of uncoupling protein 1 (UCP1), a fat browning marker, in 3T3-L1 adipocytes. Ingredients of G. uralensis such as licochalcone A (LicoA), isoliquiritigenin, and liquiritigenin induced UCP1 expression in 3T3-L1 adipocytes. After inducing obesity in mice by 6-week HFD, MeCl{sub 2} fraction of G. uralensis or LicoA was intraperitoneally administered for additional 19 days. MeCl{sub 2} fraction or LicoA significantly reduced body weight gain and inguinal fat pad weights. Furthermore, MeCl{sub 2} fraction or LicoA improved metabolic disorders induced by HFD as the treatments decreased serum levels of glucose and cholesterol, and blocked insulin resistance. MeCl{sub 2} fraction or LicoA enhanced expression of brown fat markers such as UCP1, PRDM16, and PGC-1α and increased brown fat phenotype population in inguinal WAT of HFD-fed mice. Our results demonstrate that G. uralensis and LicoA are effective to reduce obesity and to recover metabolic homeostasis by inducing the brown fat phenotype.

60 APPLIED LIFE SCIENCES↗

Scalable aqueous-phase fabrication of reduced graphene oxide nanofiltration membranes by an integrated roll-to-roll (R2R) process

Economic fractionation of multicomponent biomass-derived streams is a key challenge in biobased fuels and chemicals production. Kraft black liquor (BL) is generated at ~1 billion tons/yr globally from biomass pulping and contains ~15wt% total solids including lignin, hemicellulose fragments, and inorganics. Membrane-based BL concentration is attractive but challenged by low solute rejections and poor stability in BL, which combines alkaline pH (~13), high dissolved solids content (15+ wt%) and high temperature (70–85°C). Here we discuss the scaled fabrication and characterization of reduced graphene oxide (rGO) nanofiltration membranes by slot die coating on a roll-to-roll (R2R) with integrated vacuum filtration. Here, we obtain high-quality 90×30 cm rGO membranes (referred to as “R2R-rGO membranes” in this work) with ~100 nm average thickness supported on porous poly(ethersulfone) (PES) sheets, and an effective area of 2700cm 2 . We present characterizations of their microstructure and uniformity as mapped over large areas by scanning electron microscopy (SEM), ellipsometry, adhesion test and X-ray photoelectron spectroscopy (XPS). The constructed ellipsometry model for the rGO membrane/PES substrate composite was successfully applied to 48 different locations of the R2R-rGO membranes. The R2R-rGO membranes showed >98% lignin rejection with a steady state flux of ~12 LMH at 50bar in a 15.7wt% total solids (TS) kraft BL stream at 72°C and 2–3 gal/min crossflow rates. Slot die coating on a R2R platform with integrated vacuum filtration can enable rapid fabrication of high-quality GO membranes at scale without organic solvents or volatile organic compounds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Coupled SAM/Griffin Model of a Reference Fluoride-Salt-Cooled High-Temperature Reactor for Multi-Physics Simulations

A multi-physics coupled simulation model of a reference pebble bed fluoride-salt-cooled high-temperature reactor (PB-FHR) has been developed with SAM and Griffin computer codes for transient safety analysis. The reference problem of a prototypical reactor design serves as the foundation for the U.S. NRC (Nuclear Regulatory Commission) to verify the adequacy of computer codes and evaluation models for specific reactor types. In this work, the previously developed SAM model for PB-FHR primary loop has been updated for the coupled simulation. The updated SAM PB-FHR model includes a 2-D axial symmetric core region and external core components in 0-D/1-D. In addition to the primary loop, a detailed model of the RCCS (reactor cavity cooling system) is added. The 2-D and 1-D domains are tightly coupled using the single-solve approach developed in SAM. In the pebble bed region, the SAM multiscale explicit pebble model is applied to calculate the pebble and TRISO fuel kernel temperatures. The Griffin model used in this work is based on a model developed at Idaho National Laboratory in collaboration with the U.S. NRC. The Griffin neutronics model and SAM thermal hydraulics model is coupled with the Comprehensive Reactor Analysis Bundle (CRAB or alternately BlueCRAB) application. Both steady-state and transient scenarios are simulated to demonstrate the model's suitability for multi-physics simulations of PB-FHR transients.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Global and regional kinematics with GPS

The inherent precision of the doubly differenced phase measurement and the low cost of instrumentation made GPS the space geodetic technique of choice for regional surveys as soon as the constellation reached acceptable geometry in the area of interest: 1985 in western North America, the early 1990's in most of the world. Instrument and site-related errors for horizontal positioning are usually less than 3 mm, so that the dominant source of error is uncertainty in the reference frame defined by the satellites orbits and the tracking stations used to determine them. Prior to about 1992, when the tracking network for most experiments was globally sparse, the number of fiducial sites or the level at which they could be tied to an SLR or VLBI reference frame usually, set the accuracy limit. Recently, with a global network of over 30 stations, the limit is set more often by deficiencies in models for non-gravitational forces acting on the satellites. For regional networks in the northern hemisphere, reference frame errors are currently about 3 parts per billion (ppb) in horizontal position, allowing centimeter-level accuracies over intercontinental distances and less than 1 mm for a 100 km baseline. The accuracy of GPS measurements for monitoring height variations is generally 2-3 times worse than for horizontal motions. As for VLBI, the primary source of error is unmodeled fluctuations in atmospheric water vapor, but both reference frame uncertainties and some instrument errors are more serious for vertical than horizontal measurements. Under good conditions, daily repeatabilities at the level of 10 mm rms were achieved. This paper will summarize the current accuracy of GPS measurements and their implication for the use of SLR to study regional kinematics.

King, Robert W.↗

Multispectral Emission of the Sun during the First Whole Sun Month: Magnetohydrodynamic Simulations

We demonstrate that a three-dimensional magnetohydrodynamic (MHD) simulation of the corona can model its global plasma density and temperature structure with sufficient accuracy to reproduce many of the multispectral properties of the corona observed in extreme ultraviolet (EW) and X-ray emission. The key ingredient to this new type of global MHD model is the inclusion of energy transport processes (coronal heating, anisotropic thermal conduction, and radiative losses) in the energy equation. The calculation of these processes has previously been confined to one-dimensional loop models, idealized two-dimensional computations, and three-dimensional active region models. We refer to this as the thermodynamic MHD model, and we apply it to the time period of Carrington rotation 1913 (1996 August 22 to September 18). The form of the coronal heating term strongly affects the plasma density and temperature of the solutions. We perform our calculation for three different empirical heating models: (1) a heating function exponentially decreasing in radius; (2) the model of Schrijver et al.; and (3) a model reproducing the heating properties of the quiet Sun and active regions. We produce synthetic emission images from the density and temperature calculated with these three heating functions and quantitatively compare them with observations from E W Imaging Telescope on the Solar and Heliospheric Observatory and the soft X-ray telescope on Yohkoh. Although none of the heating models provide a perfect match, heating models 2 and 3 provide a reasonable match to the observations.

Lionello, Roberto↗

Factors Contributing to Pilot Valve Fuel Seal Extrusion in Orbiter PRCS Thrusters

Extrusion of the polytetrafluoroethylene (PTFE) pilot seal used in the monomethylhydrazine (fuel) valve of the Orbiter Primary Reaction Control System (PRCS) thrusters has been implicated in numerous on-orbit thruster failures and on-ground valve failures. Two extrusion mechanisms have been proposed, one or both may be occurring. The first mechanism is attributed to thermal expansion mismatch between adjacent PTFE and metal parts used in the fuel valve, and is referred to as thermal extrusion. The second mechanism is attributed to nitrogen tetroxide (oxidizer) leakage from the adjacent oxidizer valve on the same thruster during ground turnaround, and is referred to as oxidizer-induced extrusion. Model calculations of PTFE pilot seal in an exact pilot valve configuration show that extrusion can be caused by differential thermal expansion, without the intervening influence of oxidizer. Experimental data on semitrapped PTFE and TFM (modified PTFE) specimens simulating a fuel pilot valve configuration show that thermal extrusion 1) is incremental and irreversible, 2) increases with the size of the thermal excursion, 3) decreases with successive thermal cycling, and 4) is accompanied by gap formation. Both PTFE and TFM exhibit a higher affinity for oxidizer than fuel. The property changes associated with oxidizer uptake may explain why oxidizer seals do not exhibit extrusion. Impression replicas of fuel pilot seals removed from the Orbiter fleet show two types of extrusion: extrusion of the entire seal (loaded extrusion), or extrusion of non-sealing surface (unloaded extrusion). Both extrusion types may arise from differences in service history, rather than in failure mechanism. The plausibility oxidizer-induced extrusion was evaluated. Preliminary calculations suggest that enough energy, heat, or gas may be liberated under certain operational scenarios to cause catastrophic extrusion. However, given the lack of supporting data, conclusions implicating oxidizer leakage as a factor in extrusion must be made with caution.

Waller, J.M.↗

Assessing the Performance Limits of Internal Coronagraphs Through End-to-End Modeling

As part of the NASA ROSES Technology Demonstrations for Exoplanet Missions (TDEM) program, we conducted a numerical modeling study of three internal coronagraphs (PIAA, vector vortex, hybrid bandlimited) to understand their behaviors in realistically-aberrated systems with wavefront control (deformable mirrors). This investigation consisted of two milestones: (1) develop wavefront propagation codes appropriate for each coronagraph that are accurate to 1% or better (compared to a reference algorithm) but are also time and memory efficient, and (2) use these codes to determine the wavefront control limits of each architecture. We discuss here how the milestones were met and identify some of the behaviors particular to each coronagraph. The codes developed in this study are being made available for community use. We discuss here results for the HBLC and VVC systems, with PIAA having been discussed in a previous proceeding.

Inner working angle (IWA)↗

Decoupling the Effects of Anthropogenic Emission Reductions from the Meteorology and Natural Emissions in TROPOMI NO 2 Retrievals During the 2020 COVID-19 Lockdowns

Satellite measurements during the COVID-19 lockdowns that began in 2020 revealed unprecedented reductions in NO 2 tropospheric vertical column densities (VCD). These reductions have largely been attributed to reduced anthropogenic emissions associated with abrupt decreases in road traffic and other power consuming business activities. Although decreased emissions tended to be the main contributor to the observed NO 2 VCD reduction, meteorological variability also played a role. Whereas the observed VCD changes were predominantly negative in places where public health policies were strictly enforced, meteorology had both positive and negative effects over short time intervals. Here, we present results from a global study of the NO 2 reductions aimed at disentangling the meteorological and natural emission variability from the anthropogenic emissions over the world’s most populated megacities. For this study, NASA’s TROPOMI NO 2 algorithm was used in conjunction with the Global Modeling Initiative (GMI) chemical transport model to separate the contributions due emissions and meteorology. A priori NO 2 profiles were generated from two GMI simulations performed for 2020 at a resolution of 0.25° longitude x 0.25° latitude. The first simulation used updated, COVID-impacted NOX emissions based on Forster et al. (2020), while the second simulation used the 2019 emissions with the 2020 meteorology, referred to here as 2020BAU. The 2020BAU data set allowed for the decoupling of the emission component from the meteorology. When compared to the same period in 2019, NO 2 column amounts during the lockdowns were reduced in 35 out of 36 cities. While reduced emissions contributed most to the observed total change in NO 2 during the lockdowns, the effects of meteorology were significant, ranging between 40% (Chennai) and 15% (Beijing). In China, an increase in NO 2 levels due to meteorology were observed in five out of the seven cities considered in the study. Use of different a-priori NO 2 profiles from the two simulations in our TROPOMI retrievals allowed for a determination of retrieval errors that ranged from -1.7% to -11.0%. We used the quality assurance flag > 0.75 to select the highest quality scenes in the study period. Using the GMI, we estimated the sampling biases to be in the range -11% to 10% of the total change for the cities in our study.

Brad Fisher↗

Modernizing the Legacy Fission Wire Measurement System for the Advanced Test Reactor-Critical Facility

Operational lifetime extensions of existing research reactors have emphasized the need for refurbishment, replacements, and upgrades to supporting equipment and instrumentation. The Advanced Test Reactor (ATR) at Idaho National Laboratory (INL), which entered service in 1967, has recently completed the sixth core internals change-out and has scheduled operations until at least 2040. Reactor maintenance and operational risk management is critically important in the research reactor community, however supporting measurement systems sometimes get overlooked when maintenance is planned. The Fission Wire Measurement System (FWMS) is a custom measurement system designed in the 1960s to measure the beta-particle activity of irradiated uranium-aluminum fission wires. This measurement is conducted to determine the fission rate profile of the Advanced Reactor Test Critical (ATR-C) facility. The ATR-C is an open-pool, low-power test reactor that was purpose driven to resemble ATR and is used to qualify experiment configurations and verify core models prior to full-power experiment irradiations in ATR. A power distribution measurement in ATR-C uses uranium-aluminum wires that are distributed throughout the ATR-C core to validate simulation and modeling results. These measurements require 340 to 1500 wires to be irradiated and measured within a 12-hour window. The activity of the wires is measured in the required time with the FWMS, which was put into service in 1965 at the Radiation Measurements Laboratory (RML). The system consists of 4 measurement channels and one reference channel, each with a 2-pi proportional gas flow detector and the measurement channels each have an automated sample changer. This legacy system is crucial to the continued operations of ATR and has undergone some minor hardware upgrades since 1965, however the system presently relies on custom control boards, custom gas ion chambers, analog amplifiers/discriminators, and a user interface (UI) for the system written in outdated code. Much of the equipment and software is custom with no commercial replacements or support and limited documentation. The existing control software requires an operating system that is no longer supported, creating more vulnerabilities to continued operations. A project is underway with a third-party vendor to design, build, and document a new control and data acquisition system (CDAS) for the FWMS. The new upgrade will replace the control system, computer, UI, sample changer motors, and main power supply while maintaining the interface with existing detector hardware. The upgraded system will be operated in parallel with the current hardware and software to conduct validation testing. This equipment upgrade demonstrates the commitment at ATR to ensuring successful operations and potential future research reactors at INL.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗

Complex Dependence of Calcite Crack Kinetics on Salinity: The Role of DLVO and Hydration Forces

Abstract Subcritical crack growth (SCG) plays an important role in many geological processes such as delayed earth rupture and rock weathering. The complex dependency of SCG on the in‐crack fluid chemistry, however, is still poorly understood. In this study, we utilize the newly developed surface force‐based fracture theory (SFFT) to elucidate the relative contributions of surface forces and solute transport to the crack growth kinetics of calcite in NaCl solutions. Expanding on Barenblatt's cohesive crack model, SFFT introduces an effective stress intensity at the crack tip that encompasses all the relevant intermolecular forces across the crack in addition to the external far‐field stresses. The nonlinear system of equations portraying the crack opening profile, the solute distribution in a propagating crack, and the crack growth velocity are numerically solved via an implicit scheme. After carefully calibrating the model for calcite‐water systems, the SFFT is used to predict the SCG response of calcite at different NaCl concentrations, based on various hypotheses. These predictions are then compared to existing SCG data from the literature. We demonstrate that the experimentally observed variation of SCG rate with NaCl concentration cannot be explained solely by DLVO forces (electrostatic and Van der Waals interactions). This can be remediated by introducing an exponentially decaying hydration force with a nonlinear, nonmonotonic dependence on NaCl concentration. Furthermore, we demonstrate that accounting for both diffusive and advective transport of ions is important in explaining the absence of a stage‐II SCG response for calcite in electrolyte solutions. Plain Language Summary Subcritical crack growth (SCG) refers to the slow propagation of cracks in materials under a stress below the threshold for catastrophic failure. SCG is a key process in many geological events, for example, delayed earth ruptures and rock weathering. New initiatives such as underground CO 2 and H 2 storage in carbonate reservoirs further call for better understanding of SCG in carbonate minerals subjected to varying fluid chemistry. This study examines the SCG of calcite, a key mineral found in carbonate rocks, intergranular cement in sandstones, and filling material in mineral veins and faults, determining their deformation and strength. A mathematical model is developed to describe how the crack opens and propagates, how solutes (like salts) distribute within the crack, and how the crack surfaces interact with each other. We used the model to predict calcite SCG in water at different salt concentrations and compared it with experimental data. Our results revealed that the hydration force is the dominating factor in determining the complex, non‐linear dependency of SCG on salinity. We also found that both the movement of ions by diffusion and by bulk water flow are crucial for explaining the SCG rates, especially when the cracks grow quickly. Key Points Surface Force‐Based Fracture Theory predicts the complex subcritical crack growth patterns of calcite crystals immersed in NaCl solutions Results highlight the dominant role of hydration forces in altering the fracture behavior of calcite compared to VdW and electric double‐layer forces Advective solute transport explains the absence of stages‐II and ‐III subcritical crack growth responses in solid‐liquid systems

DLVO↗

Modified Closures in Monte Carlo Algorithms for Diffusive Binary Stochastic Media Transport Problems

In a stochastic medium, the material properties at a given spatial location are known only statistically. The most common approach to solving particle transport problems involving binary stochastic media (BSM) is to use the atomic mix approximation in which the transport problem is solved using homogenized (volume-averaged) material properties. A common deterministic model developed for solving BSM particle transport problems is the Levermore-Pomraning (LP) model. Zimmerman and Adams proposed a Monte Carlo algorithm (Algorithm A) that solves the LP equations and another Monte Carlo algorithm (Algorithm B) that locally preserves the sampled material realization; we refer to these Monte Carlo algorithms as the LP and LPLRP (LP local realization preserving) algorithms, respectively. One-dimensional (1D) planar geometry benchmark studies have shown that the LPLRP algorithm is often significantly more accurate than the LP algorithm for problems with an incident angular flux as well as for problems with an interior source. The LPLRP algorithm implementation in these benchmark comparisons made explicit use of the one-dimensional nature of the problem. The LP model is derived by assuming an 'upwind' closure in the coupling term relating the two materials in a binary stochastic medium. Su and Pomraning developed a modified form of this closure by considering the small correlation length limit and requiring the modified closure to produce the correct exponential decay for a source-free half-line albedo problem in rod geometry. They concluded that the modified closure is generally not inferior to the LP closure and in some cases is significantly better. Brantley investigated the use of the Su-Pomraning (SP) closure in the Monte Carlo LP algorithm (LP-SP) for the suite of benchmark problems and concluded that 1) the LP-SP algorithm was somewhat more accurate overall than the LP algorithm and somewhat less accurate overall than the LPLRP algorithm for an incident angular flux benchmark suite, and 2) the LP-SP algorithm was generally the least accurate of the algorithms for an interior source benchmark suite. Larsen et al. performed an asymptotic analysis of the transport equation in 1D planar geometry for the situation in which the physical system is 1) a random binary stochastic medium with material macroscopic total cross sections and mean slab width values of O(1) or smaller and 2) optically thick with weak absorption and sources at each spatial point and therefore globally diffusive. The asymptotic analysis demonstrates that, under these assumptions, the transport equation limits to the conventional diffusion equation with atomically-mixed (volume-averaged) material properties. Larsen et al. also demonstrate that, under these same assumptions, the LP equations asymptotically limit to an atomically-mixed diffusion equation with a diffusion coefficient that is too large, leading to an unphysical flattening of the ensemble-averaged scalar flux distribution. Vasques and Yadav recently performed an asymptotic analysis of an adjusted Levermore-Pomraning closure in which the Markovian transition functions are rescaled such that the equations asymptotically limit to the correct atomically-mixed diffusion equation; we refer to this model as LP-VY. In this paper, we demonstrate that the adjusted LP closure proposed by Vasques and Yadav is a special case of the SP closure obtained by assuming that the absorption cross sections in the two materials are equal. We further demonstrate that the Su-Pomraning closure has the correct asymptotic limit for the diffusive physical system and is therefore an appropriate closure for these problems. We also describe how to incorporate these various closures into the Monte Carlo LP and LPLRP particle transport algorithms, and we present numerical results comparing the accuracy of the Monte Carlo LP, LP-VY, LP-SP, LPLRP, and LPLRP-SP algorithms for the set of diffusive benchmark problems. We demonstrated that the SP closure has the correct asymptotic limit for the diffusive physical system under investigation in this paper and that the VY closure is a special case of the SP closure obtained by assuming equal absorption in the two materials. Both of these closures are readily implemented in the Monte Carlo LP and LPLRP algorithms for particle transport in BSM. Through numerical comparisons to a diffusive benchmark suite, we confirmed that 1) the LP and LPLRP algorithms do not limit to the appropriate diffusion solution and 2) the LP-VY, LP-SP, and LPLRP-SP algorithms limit to the appropriate diffusion solution and are of comparable accuracy. The LPLRP algorithm has been shown to be generally more accurate than the LP algorithm for non-diffusive problems. However, the numerical results in this paper indicate that the LPLRP algorithm does not limit to the correct diffusion equation for the asymptotic limit under consideration. The work in this paper demonstrates that the SP algorithm does limit to the correct diffusion equation for the asymptotic limit considered. Future work will investigate the possibility of whether the LPLRP-SP algorithm is accurate for both diffusive and non-diffusive problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Measurements of azimuthal anisotropy and charged-particle multiplicity in d+Au collisions at √ s NN=200, 62.4, 39, and 19.6 GeV

We present measurements of the elliptic flow (v 2 ) as a function of transverse momentum (p T ), pseudorapidity (η), and centrality in d+Au collisions at √ s NN = 200, 62.4, 39, and 19.6 GeV. The beam-energy scan of d+Au collisions provides a testing ground for the onset of ow signatures in small collision systems. We measure a nonzero v 2 signal at all four collision energies, which, at midrapidity and low p T , is consistent with predictions from viscous hydrodynamic models. Comparisons with calculations from parton transport models (based on the ampt Monte Carlo generator) show good agreement with the data at midrapidity to forward (d-going) rapidities and low p T . At backward (Au-going) rapidities and p T > 1:5 GeV/c, the data diverges from ampt calculations of v 2 relative to the initial geometry, indicating the possible dominance of nongeometry related corre- lations, referred to as nonflow. We also present measurements of the charged-particle multiplicity (dN ch /d ) as a function of η in central d+Au collisions at the same energies. We find that in d+Au collisions at √ s NN = 200 GeV the v 2 scales with dN ch /d over all in the PHENIX acceptance. At √ s NN = 62:4, and 39 GeV, v 2 scales with dN ch /d at midrapidity and forward rapidity, but falls o at backward rapidity. Furthermore, this departure from the dN ch /dη scaling may be a further indication of non ow effects dominating at backward rapidity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Performance and Cost Assessment of a Natural Gas-Fueled Direct sCO 2 Power Plant

This final report details the cost and performance baseline developed at NETL for a natural gas-fueled direct sCO 2 plant, which holds promise as an efficient, cost-effective method of power production with carbon capture and storage (CCS). The steady state plant model includes turbine cooling flows based on a semi-empirical cooling model, and thermal integration with the air separation unit and oxygen compression train. Discussion of the performance of the baseline plant design is followed by details on the costing methodology, overall plant capital cost, and resulting cost of electricity (COE). Performance results yield a baseline sCO 2 plant thermal efficiency of 48.2% (HHV) with 98.2% carbon capture, an improvement over the reference NGCC plant with CCS (45.7% and 90.7% carbon capture). Cost analyses yield a 5% reduction in COE versus the NGCC plant, albeit with significant uncertainty. Finally, the effects of incomplete combustion on the overall plant performance are assessed.

03 NATURAL GAS↗

Lunar motion analysis and laser data management

Work completed in lunar motion analysis and laser data management during the period July 1, 1971 - September 30, 1975 was reported. In this context, analysis refers to theoretical or numerical studies involving real or potential applications of such observations to improvement of the physical model, and data management refers to the process by which observed photon events are turned into observations and are made available to potential users. The data analysis work included: (1) bringing to operational status of computer programs for the numerical integration of the lunar orbit motion and for the application of lunar laser time delays for the improvement of the parameters of the physical model, (2) program improvement and program integrity, (3) three-dimensional ephemeris, and (4) miscellaneous independent studies. The data management work included: (1) data identification, (2) observatory interfaces, and (3) data distribution.

Mulholland, J. D.↗

The Effect of Geocenter Motion on Jason-2 Orbits and the Mean Sea Level

We compute a series of Jason-2 GPS and SLR/DORIS-based orbits using ITRF2005 and the std0905 standards (Lemoine et al. 2010). Our GPS and SLR/DORIS orbit data sets span a period of 2 years from cycle 3 (July 2008) to cycle 74 (July 2010). We extract the Jason-2 orbit frame translational parameters per cycle by the means of a Helmert transformation between a set of reference orbits and a set of test orbits. We compare the annual terms of these time-series to the annual terms of two different geocenter motion models where biases and trends have been removed. Subsequently, we include the annual terms of the modeled geocenter motion as a degree-1 loading displacement correction to the GPS and SLR/DORIS tracking network of the POD process. Although the annual geocenter motion correction would reflect a stationary signal in time, under ideal conditions, the whole geocenter motion is a non-stationary process that includes secular trends. Our results suggest that our GSFC Jason-2 GPS-based orbits are closely tied to the center of mass (CM) of the Earth consistent with our current force modeling, whereas GSFC's SLR/DORIS-based orbits are tied to the origin of ITRF2005, which is the center of figure (CF) for sub-secular scales. We quantify the GPS and SLR/DORIS orbit centering and how this impacts the orbit radial error over the globe, which is assimilated into mean sea level (MSL) error, from the omission of the annual term of the geocenter correction. We find that for the SLR/DORIS std0905 orbits, currently used by the oceanographic community, only the negligence of the annual term of the geocenter motion correction results in a 4.67 plus or minus 3.40 mm error in the Z-component of the orbit frame which creates 1.06 plus or minus 2.66 mm of systematic error in the MSL estimates, mainly due to the uneven distribution of the oceans between the North and South hemisphere.

Melachroinos, S. A.↗