National Library of Energy BETA

Sample records for lcoe numerical modeling

  1. numerical modeling | OpenEI Community

    Open Energy Info (EERE)

    Submitted by Ocop(5) Member 15 July, 2014 - 07:07 MHK LCOE Reporting Guidance Draft Cost Current DOE LCOE numerical modeling Performance Tidal Wave To normalize competing...

  2. LCOEs and Renewables

    Gasoline and Diesel Fuel Update (EIA)

    Victor Niemeyer Program Manager, Energy and Environmental Policy Analysis and Company Strategy Program EIA LCOE/LACE Workshop July 25, 2013 LCOEs and Renewables 2 © 2013 Electric Power Research Institute, Inc. All rights reserved. EPRI Generation Options Report Provides Excellent Example of LCOE Use By Robin Bedillion of EPRI's Strategic Energy Analysis Group Reference: EPRI Report 1026656 (free from EPRI.com, search for "1026656") 3 © 2013 Electric Power Research Institute, Inc. All

  3. LCOE | OpenEI Community

    Open Energy Info (EERE)

    LCOE Home Kch's picture Submitted by Kch(24) Member 15 July, 2014 - 07:07 MHK Cost Breakdown Structure Draft CBS current energy GMREC LCOE levelized cost of energy marine energy...

  4. Economic Competitiveness of U.S. Utility-Scale Photovoltaics Systems in 2015: Regional Cost Modeling of Installed Cost ($/W) and LCOE ($/kWh)

    SciTech Connect (OSTI)

    Fu, Ran; James, Ted L.; Chung, Donald; Gagne, Douglas; Lopez, Anthony; Dobos, Aron

    2015-06-14

    Utility-scale photovoltaics (PV) system growth is largely driven by the economic metrics of total installed costs and levelized cost of electricity (LCOE), which differ by region. This study details regional cost factors, including environment (wind speed and snow loads), labor costs, material costs, sales taxes, and permitting costs using a new system-level bottom-up cost modeling approach. We use this model to identify regional all-in PV installed costs for fixed-tilt and one-axis tracker systems in the United States with consideration of union and non-union labor costs in 2015. LCOEs using those regional installed costs are then modeled and spatially presented. Finally, we assess the cost reduction opportunities of increasing module conversion efficiencies on PV system costs in order to indicate the possible economic impacts of module technology advancements and help future research and development (R&D) effects in the context of U.S. SunShot targets.

  5. LCOE and LACE: Seeking Simplicity

    Gasoline and Diesel Fuel Update (EIA)

    LCOE and LACE Seeking Simplicity Michael Leifman GE Power & Water Marketing Presented to the EIA LCOE/LACE Workshop July 25, 2013 © 2013 General Electric Company. All Rights Reserved 2 GE Power & Water July 25, 2013 A love-hate relationship * We love Levelized Cost of Electricity (LCOE) because it is simple * Everything is boiled down to one number * Provides an easy point of comparison across technologies * Serves as a shorthand for many dimensions BUT * LCOE masks many site-specific

  6. LCOE Uncertainty Analysis for Hydropower using Monte Carlo Simulations

    SciTech Connect (OSTI)

    Chalise, Dol Raj; O'Connor, Patrick W; DeNeale, Scott T; Uria Martinez, Rocio; Kao, Shih-Chieh

    2015-01-01

    Levelized Cost of Energy (LCOE) is an important metric to evaluate the cost and performance of electricity production generation alternatives, and combined with other measures, can be used to assess the economics of future hydropower development. Multiple assumptions on input parameters are required to calculate the LCOE, which each contain some level of uncertainty, in turn affecting the accuracy of LCOE results. This paper explores these uncertainties, their sources, and ultimately the level of variability they introduce at the screening level of project evaluation for non-powered dams (NPDs) across the U.S. Owing to site-specific differences in site design, the LCOE for hydropower varies significantly from project to project unlike technologies with more standardized configurations such as wind and gas. Therefore, to assess the impact of LCOE input uncertainty on the economics of U.S. hydropower resources, these uncertainties must be modeled across the population of potential opportunities. To demonstrate the impact of uncertainty, resource data from a recent nationwide non-powered dam (NPD) resource assessment (Hadjerioua et al., 2012) and screening-level predictive cost equations (O Connor et al., 2015) are used to quantify and evaluate uncertainties in project capital and operations & maintenance costs, and generation potential at broad scale. LCOE dependence on financial assumptions is also evaluated on a sensitivity basis to explore ownership/investment implications on project economics for the U.S. hydropower fleet. The results indicate that the LCOE for U.S. NPDs varies substantially. The LCOE estimates for the potential NPD projects of capacity greater than 1 MW range from 40 to 182 $/MWh, with average of 106 $/MWh. 4,000 MW could be developed through projects with individual LCOE values below 100 $/MWh. The results also indicate that typically 90 % of LCOE uncertainty can be attributed to uncertainties in capital costs and energy production; however

  7. New Report Shows Downward Trend in LCOE for Wind | Department...

    Office of Environmental Management (EM)

    New Report Shows Downward Trend in LCOE for Wind New Report Shows Downward Trend in LCOE for Wind May 18, 2015 - 2:48pm Addthis A new report recently published by the U.S. ...

  8. Concentrating Solar Power: Advanced Projects Offering Low LCOE...

    Broader source: Energy.gov (indexed) [DOE]

    Advanced Projects Offering Low LCOE Opportunities (CSP: APOLLO) Building upon the successful outcomes of the 2012 SunShot Concentrating Solar Power (CSP) Research & Development ...

  9. 51-Mile Hydroelectric Power Project Demonstration of new methodologies to reduce the LCOE for small, hydropower development

    Broader source: Energy.gov [DOE]

    51-Mile Hydroelectric Power Project Demonstration of new methodologies to reduce the LCOE for small, hydropower development

  10. Development of Reference Models and Design Tools (LCOE Models...

    Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site

    FY 09 Lab Call: Research & Assessment for MHK Development 2014 Water Power Program Peer Review Compiled Presentations: Marine and Hydrokinetic Technologies Effects on the Physical ...

  11. Levelized cost of energy (LCOE) metric to characterize solar absorber coatings for the CSP industry

    SciTech Connect (OSTI)

    Boubault, Antoine; Ho, Clifford K.; Hall, Aaron; Lambert, Timothy N.; Ambrosini, Andrea

    2015-07-08

    The contribution of each component of a power generation plant to the levelized cost of energy (LCOE) can be estimated and used to increase the power output while reducing system operation and maintenance costs. The LCOE is used in order to quantify solar receiver coating influence on the LCOE of solar power towers. Two new parameters are introduced: the absolute levelized cost of coating (LCOC) and the LCOC efficiency. Depending on the material properties, aging, costs, and temperature, the absolute LCOC enables quantifying the cost-effectiveness of absorber coatings, as well as finding optimal operating conditions. The absolute LCOC is investigated for different hypothetic coatings and is demonstrated on Pyromark 2500 paint. Results show that absorber coatings yield lower LCOE values in most cases, even at significant costs. Optimal reapplication intervals range from one to five years. At receiver temperatures greater than 700 C, non-selective coatings are not always worthwhile while durable selective coatings consistently reduce the LCOEup to 12% of the value obtained for an uncoated receiver. Moreover the absolute LCOC is a powerful tool to characterize and compare different coatings, not only considering their initial efficiencies but also including their durability.

  12. LBNL/NREL Analysis Predicts Record Low LCOE for Wind Energy in 2012-2013 |

    Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site

    Department of Energy LBNL/NREL Analysis Predicts Record Low LCOE for Wind Energy in 2012-2013 LBNL/NREL Analysis Predicts Record Low LCOE for Wind Energy in 2012-2013 February 24, 2012 - 11:27am Addthis This is an excerpt from the First Quarter 2012 edition of the Wind Program R&D Newsletter. A recent analysis conducted by the Lawrence Berkley National Laboratory (LBNL) and the National Renewable Energy Laboratory (NREL) suggests that lower capital costs and continued increases in wind

  13. Levelized cost of energy (LCOE) metric to characterize solar absorber coatings for the CSP industry

    SciTech Connect (OSTI)

    Boubault, Antoine; Ho, Clifford K.; Hall, Aaron; Lambert, Timothy N.; Ambrosini, Andrea

    2015-07-08

    The contribution of each component of a power generation plant to the levelized cost of energy (LCOE) can be estimated and used to increase the power output while reducing system operation and maintenance costs. The LCOE is used in order to quantify solar receiver coating influence on the LCOE of solar power towers. Two new parameters are introduced: the absolute levelized cost of coating (LCOC) and the LCOC efficiency. Depending on the material properties, aging, costs, and temperature, the absolute LCOC enables quantifying the cost-effectiveness of absorber coatings, as well as finding optimal operating conditions. The absolute LCOC is investigated for different hypothetic coatings and is demonstrated on Pyromark 2500 paint. Results show that absorber coatings yield lower LCOE values in most cases, even at significant costs. Optimal reapplication intervals range from one to five years. At receiver temperatures greater than 700 °C, non-selective coatings are not always worthwhile while durable selective coatings consistently reduce the LCOE—up to 12% of the value obtained for an uncoated receiver. Moreover the absolute LCOC is a powerful tool to characterize and compare different coatings, not only considering their initial efficiencies but also including their durability.

  14. Levelized cost of energy (LCOE) metric to characterize solar absorber coatings for the CSP industry

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Boubault, Antoine; Ho, Clifford K.; Hall, Aaron; Lambert, Timothy N.; Ambrosini, Andrea

    2015-07-08

    The contribution of each component of a power generation plant to the levelized cost of energy (LCOE) can be estimated and used to increase the power output while reducing system operation and maintenance costs. The LCOE is used in order to quantify solar receiver coating influence on the LCOE of solar power towers. Two new parameters are introduced: the absolute levelized cost of coating (LCOC) and the LCOC efficiency. Depending on the material properties, aging, costs, and temperature, the absolute LCOC enables quantifying the cost-effectiveness of absorber coatings, as well as finding optimal operating conditions. The absolute LCOC is investigatedmore » for different hypothetic coatings and is demonstrated on Pyromark 2500 paint. Results show that absorber coatings yield lower LCOE values in most cases, even at significant costs. Optimal reapplication intervals range from one to five years. At receiver temperatures greater than 700 °C, non-selective coatings are not always worthwhile while durable selective coatings consistently reduce the LCOE—up to 12% of the value obtained for an uncoated receiver. Moreover the absolute LCOC is a powerful tool to characterize and compare different coatings, not only considering their initial efficiencies but also including their durability.« less

  15. MHK LCOE Reporting Guidance Draft | OpenEI Community

    Open Energy Info (EERE)

    details on the Reference Models can be found here. The presentation was created in Microsoft Powerpoint, therefore it is recommended that Microsoft Powerpoint or Powerpoint...

  16. Funding Opportunity Announcement: Concentrating Solar Power: Advanced Projects Offering Low LCOE Opportunities

    Broader source: Energy.gov [DOE]

    The SunShot Initiative's Concentrating Solar Power: Advanced Projects Offering Low LCOE Opportunities (CSP: APOLLO) funding opportunity announcement (FOA) seeks transformative projects targeting all components of a concentrating solar power (CSP) plant. Projects should seek to meet the targets set out in the SunShot Vision Study , enabling CSP to become fully cost-competitive with traditional forms of electric power generation. Projects can address challenges in any technical system of the plant, including solar collectors, receivers and heat transfer fluids, thermal energy storage, power cycles, as well as operations and maintenance.

  17. Numerical modeling of water injection into vapor-dominatedgeothermal...

    Office of Scientific and Technical Information (OSTI)

    Technical Report: Numerical modeling of water injection into vapor-dominatedgeothermal reservoirs Citation Details In-Document Search Title: Numerical modeling of water injection ...

  18. Numerical Modeling Studies of The Dissolution-Diffusion-Convection...

    Office of Scientific and Technical Information (OSTI)

    Numerical Modeling Studies of The Dissolution-Diffusion-Convection ProcessDuring CO2 Storage in Saline Aquifers Citation Details In-Document Search Title: Numerical Modeling ...

  19. Numerical Modeling At Dixie Valley Geothermal Area (McKenna ...

    Open Energy Info (EERE)

    Numerical Modeling At Dixie Valley Geothermal Area (McKenna & Blackwell, 2003) Jump to: navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Numerical Modeling...

  20. Wave | OpenEI Community

    Open Energy Info (EERE)

    Submitted by Ocop(5) Member 15 July, 2014 - 07:07 MHK LCOE Reporting Guidance Draft Cost Current DOE LCOE numerical modeling Performance Tidal Wave To normalize competing...

  1. Performance | OpenEI Community

    Open Energy Info (EERE)

    Submitted by Ocop(5) Member 15 July, 2014 - 07:07 MHK LCOE Reporting Guidance Draft Cost Current DOE LCOE numerical modeling Performance Tidal Wave To normalize competing...

  2. Current | OpenEI Community

    Open Energy Info (EERE)

    Submitted by Ocop(5) Member 15 July, 2014 - 07:07 MHK LCOE Reporting Guidance Draft Cost Current DOE LCOE numerical modeling Performance Tidal Wave To normalize competing...

  3. Tidal | OpenEI Community

    Open Energy Info (EERE)

    Submitted by Ocop(5) Member 15 July, 2014 - 07:07 MHK LCOE Reporting Guidance Draft Cost Current DOE LCOE numerical modeling Performance Tidal Wave To normalize competing...

  4. Advanced Numerical Model for Irradiated Concrete

    SciTech Connect (OSTI)

    Giorla, Alain B.

    2015-03-01

    In this report, we establish a numerical model for concrete exposed to irradiation to address these three critical points. The model accounts for creep in the cement paste and its coupling with damage, temperature and relative humidity. The shift in failure mode with the loading rate is also properly represented. The numerical model for creep has been validated and calibrated against different experiments in the literature [Wittmann, 1970, Le Roy, 1995]. Results from a simplified model are shown to showcase the ability of numerical homogenization to simulate irradiation effects in concrete. In future works, the complete model will be applied to the analysis of the irradiation experiments of Elleuch et al. [1972] and Kelly et al. [1969]. This requires a careful examination of the experimental environmental conditions as in both cases certain critical information are missing, including the relative humidity history. A sensitivity analysis will be conducted to provide lower and upper bounds of the concrete expansion under irradiation, and check if the scatter in the simulated results matches the one found in experiments. The numerical and experimental results will be compared in terms of expansion and loss of mechanical stiffness and strength. Both effects should be captured accordingly by the model to validate it. Once the model has been validated on these two experiments, it can be applied to simulate concrete from nuclear power plants. To do so, the materials used in these concrete must be as well characterized as possible. The main parameters required are the mechanical properties of each constituent in the concrete (aggregates, cement paste), namely the elastic modulus, the creep properties, the tensile and compressive strength, the thermal expansion coefficient, and the drying shrinkage. These can be either measured experimentally, estimated from the initial composition in the case of cement paste, or back-calculated from mechanical tests on concrete. If some

  5. NUMERICAL MODELING OF CATHODE CONTACT MATERIAL DENSIFICATION

    SciTech Connect (OSTI)

    Koeppel, Brian J.; Liu, Wenning N.; Stephens, Elizabeth V.; Khaleel, Mohammad A.

    2011-11-01

    Numerical modeling was used to simulate the constrained sintering process of the cathode contact layer during assembly of solid oxide fuel cells (SOFCs). A finite element model based on the continuum theory for sintering of porous bodies was developed and used to investigate candidate low-temperature cathode contact materials. Constitutive parameters for various contact materials under investigation were estimated from dilatometry screening tests, and the influence of processing time, processing temperature, initial grain size, and applied compressive stress on the free sintering response was predicted for selected candidate materials. The densification behavior and generated stresses within a 5-cell planar SOFC stack during sintering, high temperature operation, and room temperature shutdown were predicted. Insufficient constrained densification was observed in the stack at the proposed heat treatment, but beneficial effects of reduced grain size, compressive stack preload, and reduced thermal expansion coefficient on the contact layer densification and stresses were observed.

  6. An Updated Numerical Model Of The Larderello-Travale Geothermal...

    Open Energy Info (EERE)

    Numerical Model Of The Larderello-Travale Geothermal System, Italy Jump to: navigation, search OpenEI Reference LibraryAdd to library Journal Article: An Updated Numerical Model Of...

  7. Numerical Modeling Of Basin And Range Geothermal Systems | Open...

    Open Energy Info (EERE)

    for extensional geothermal systems that include structure, heat input, and permeability distribution have been established using numerical models. Extensional geothermal...

  8. Numerical Modelling of Geothermal Systems a Short Introduction...

    Open Energy Info (EERE)

    Modelling of Geothermal Systems a Short Introduction Jump to: navigation, search OpenEI Reference LibraryAdd to library General: Numerical Modelling of Geothermal Systems a Short...

  9. Numerical Modeling of PCCI Combustion | Department of Energy

    Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site

    PCCI Combustion Numerical Modeling of PCCI Combustion 2004 Diesel Engine Emissions Reduction (DEER) Conference Presentation: Lawrence Livermore National Laboratory/University of Michigan 2004_deer_flowers.pdf (252.97 KB) More Documents & Publications Modeling of HCCI and PCCI Combustion Processes Bridging the Gap between Fundamental Physics and Chemistry and Applied Models for HCCI Engines Numerical Modeling of HCCI Combustion

  10. Numerical Modeling At Raft River Geothermal Area (1983) | Open...

    Open Energy Info (EERE)

    Raft River Geothermal Area (1983) Jump to: navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Numerical Modeling At Raft River Geothermal Area (1983)...

  11. Final Report on Experimental and Numerical Modeling Activities...

    Office of Scientific and Technical Information (OSTI)

    Technical Report: Final Report on Experimental and Numerical Modeling Activities for the Newark Basin Citation Details In-Document Search Title: Final Report on Experimental and ...

  12. Numerical Modeling At Coso Geothermal Area (1995) | Open Energy...

    Open Energy Info (EERE)

    transform is employed to characterize guided-wave's velocity-frequency dispersion, and numerical methods are used to simulate the guided-wave propagation. The modeling...

  13. Numerical Modeling At Lightning Dock Geothermal Area (O'Brien...

    Open Energy Info (EERE)

    Basin Additional References Retrieved from "http:en.openei.orgwindex.php?titleNumericalModelingAtLightningDockGeothermalArea(O%27Brien,EtAl.,1984)&oldid762871...

  14. Mathematical and Numerical Analyses of Peridynamics for Multiscale Materials Modeling

    SciTech Connect (OSTI)

    Gunzburger, Max

    2015-02-17

    We have treated the modeling, analysis, numerical analysis, and algorithmic development for nonlocal models of diffusion and mechanics. Variational formulations were developed and finite element methods were developed based on those formulations for both steady state and time dependent problems. Obstacle problems and optimization problems for the nonlocal models were also treated and connections made with fractional derivative models.

  15. Use of ARM observations and numerical models to determine radiative...

    Office of Scientific and Technical Information (OSTI)

    We investigated whether the West African anvil clouds connected with squall line MCSs passing over the Niamey ARM site could be simulated in a numerical model by comparing the ...

  16. Rapid installation of numerical models in multiple parent codes

    SciTech Connect (OSTI)

    Brannon, R.M.; Wong, M.K.

    1996-10-01

    A set of``model interface guidelines``, called MIG, is offered as a means to more rapidly install numerical models (such as stress-strain laws) into any parent code (hydrocode, finite element code, etc.) without having to modify the model subroutines. The model developer (who creates the model package in compliance with the guidelines) specifies the model`s input and storage requirements in a standardized way. For portability, database management (such as saving user inputs and field variables) is handled by the parent code. To date, NUG has proved viable in beta installations of several diverse models in vectorized and parallel codes written in different computer languages. A NUG-compliant model can be installed in different codes without modifying the model`s subroutines. By maintaining one model for many codes, MIG facilitates code-to-code comparisons and reduces duplication of effort potentially reducing the cost of installing and sharing models.

  17. Numerical Modeling of HCCI Combustion | Department of Energy

    Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site

    HCCI Combustion Numerical Modeling of HCCI Combustion Presentation given at DEER 2006, August 20-24, 2006, Detroit, Michigan. Sponsored by the U.S. DOE's EERE FreedomCar and Fuel Partnership and 21st Century Truck Programs. 2006_deer_aceves.pdf (840 KB) More Documents & Publications High Fidelity Modeling of Premixed Charge Compression Ignition Engines New Methodologies for Analysis of Premixed Charge Compression Ignition Engines Modeling of HCCI and PCCI Combustion Processes

  18. Numerical Modeling At Neal Hot Springs Geothermal Area (U.S....

    Open Energy Info (EERE)

    Area Exploration Technique Numerical Modeling Activity Date 2011 - 2011 Usefulness useful DOE-funding Unknown Exploration Basis A numerical reservoir model was created to...

  19. AEETES---A solar reflux receiver thermal performance numerical model

    SciTech Connect (OSTI)

    Hogan, R.E. Jr.

    1991-01-01

    Reflux solar receivers for dish-Stirling electric power generation systems are currently being investigated by several companies and laboratories. In support of these efforts, the AEETES thermal performance numerical model has been developed to predict thermal performance of pool-boiler and heat-pipe reflux receivers. The formulation of the AEETES numerical model, which is applicable to axisymmetric geometries with asymmetric incident fluxes, is presented in detail. Thermal efficiency predictions agree to within 4.1% with test data from on-sun tests of a pool-boiler reflux receiver. Predicted absorber and sidewall temperatures agree with thermocouple data to within 3.3.% and 7.3%, respectively. The importance of accounting for the asymmetric incident fluxes is demonstrated in comparisons with predictions using azimuthally averaged variables. The predicted receiver heat losses are characterized in terms of convective, solar and infrared radiative, and conductive heat transfer mechanisms. 27 refs., 9 figs., 4 tabs.

  20. ASSIMILATION OF DOPPLER RADAR DATA INTO NUMERICAL WEATHER MODELS

    SciTech Connect (OSTI)

    Chiswell, S.; Buckley, R.

    2009-01-15

    During the year 2008, the United States National Weather Service (NWS) completed an eight fold increase in sampling capability for weather radars to 250 m resolution. This increase is expected to improve warning lead times by detecting small scale features sooner with increased reliability; however, current NWS operational model domains utilize grid spacing an order of magnitude larger than the radar data resolution, and therefore the added resolution of radar data is not fully exploited. The assimilation of radar reflectivity and velocity data into high resolution numerical weather model forecasts where grid spacing is comparable to the radar data resolution was investigated under a Laboratory Directed Research and Development (LDRD) 'quick hit' grant to determine the impact of improved data resolution on model predictions with specific initial proof of concept application to daily Savannah River Site operations and emergency response. Development of software to process NWS radar reflectivity and radial velocity data was undertaken for assimilation of observations into numerical models. Data values within the radar data volume undergo automated quality control (QC) analysis routines developed in support of this project to eliminate empty/missing data points, decrease anomalous propagation values, and determine error thresholds by utilizing the calculated variances among data values. The Weather Research and Forecasting model (WRF) three dimensional variational data assimilation package (WRF-3DVAR) was used to incorporate the QC'ed radar data into input and boundary conditions. The lack of observational data in the vicinity of SRS available to NWS operational models signifies an important data void where radar observations can provide significant input. These observations greatly enhance the knowledge of storm structures and the environmental conditions which influence their development. As the increase in computational power and availability has made higher

  1. Mathematical and Numerical Analyses of Peridynamics for Multiscale Materials Modeling

    SciTech Connect (OSTI)

    Du, Qiang

    2014-11-12

    The rational design of materials, the development of accurate and efficient material simulation algorithms, and the determination of the response of materials to environments and loads occurring in practice all require an understanding of mechanics at disparate spatial and temporal scales. The project addresses mathematical and numerical analyses for material problems for which relevant scales range from those usually treated by molecular dynamics all the way up to those most often treated by classical elasticity. The prevalent approach towards developing a multiscale material model couples two or more well known models, e.g., molecular dynamics and classical elasticity, each of which is useful at a different scale, creating a multiscale multi-model. However, the challenges behind such a coupling are formidable and largely arise because the atomistic and continuum models employ nonlocal and local models of force, respectively. The project focuses on a multiscale analysis of the peridynamics materials model. Peridynamics can be used as a transition between molecular dynamics and classical elasticity so that the difficulties encountered when directly coupling those two models are mitigated. In addition, in some situations, peridynamics can be used all by itself as a material model that accurately and efficiently captures the behavior of materials over a wide range of spatial and temporal scales. Peridynamics is well suited to these purposes because it employs a nonlocal model of force, analogous to that of molecular dynamics; furthermore, at sufficiently large length scales and assuming smooth deformation, peridynamics can be approximated by classical elasticity. The project will extend the emerging mathematical and numerical analysis of peridynamics. One goal is to develop a peridynamics-enabled multiscale multi-model that potentially provides a new and more extensive mathematical basis for coupling classical elasticity and molecular dynamics, thus enabling next

  2. Numerical simulation model for vertical flow in geothermal wells

    SciTech Connect (OSTI)

    Tachimori, M.

    1982-01-01

    A numerical simulation model for vertical flow in geothermal wells is presented. The model consists of equations for the conservation of mass, momentum, and energy, for thermodynamic state of water, for friction losses, for slip velocity relations, and of the criteria for various flow regimes. A new set of correlations and criteria is presented for two-phase flow to improve the accuracy of predictions; bubbly flow - Griffith and Wallis correlation, slug flow - Nicklin et al. one, annular-mist flow - Inoue and Aoki and modified by the author. The simulation method was verified by data from actual wells.

  3. Numerical Model for Conduction-Cooled Current Lead Heat Loads

    SciTech Connect (OSTI)

    White, M.J.; Wang, X.L.; Brueck, H.D.; /DESY

    2011-06-10

    Current leads are utilized to deliver electrical power from a room temperature junction mounted on the vacuum vessel to a superconducting magnet located within the vacuum space of a cryostat. There are many types of current leads used at laboratories throughout the world; however, conduction-cooled current leads are often chosen for their simplicity and reliability. Conduction-cooled leads have the advantage of using common materials, have no superconducting/normal state transition, and have no boil-off vapor to collect. This paper presents a numerical model for conduction-cooled current lead heat loads. This model takes into account varying material and fluid thermal properties, varying thicknesses along the length of the lead, heat transfer in the circumferential and longitudinal directions, electrical power dissipation, and the effect of thermal intercepts. The model is validated by comparing the numerical model results to ideal cases where analytical equations are valid. In addition, the XFEL (X-Ray Free Electron Laser) prototype current leads are modeled and compared to the experimental results from testing at DESY's XFEL Magnet Test Stand (XMTS) and Cryomodule Test Bench (CMTB).

  4. Comparison of numerical models of a pyrotechnic device

    SciTech Connect (OSTI)

    Pierce, K.G.

    1986-01-01

    The predictions of two numerical models of a hot-wire initiated pyrotechnic device are compared to each other and to experimental results. Both models employ finite difference approximations to the heat diffusion equation in cylindrical coordinates. The temperature dependence of the thermal properties of the pyrotechnic materials and of the bridgewire are modeled. An Arrhenius' model is used to describe the exothermic reaction in the powder. One model employs a single radial coordinate and predicts the radial temperature distribution in the bridgewire and surrounding powder mass. In addition to the radial coordinate, the other model also employs a longitudinal coordinate to predict the temperature distribution parallel to the axis of the bridgewire. The predictions of the two-dimensional model concerning the energy requirements for ignition and the energy losses from the ends of the bridgewire are presented. A comparison of the predictions of the two models and the development of thermal gradients are employed to define the regime where the assumption, in the one-dimensional model, of negligible heat transfer axial to the bridgewire does not lead to significant error. The general problems associated with predicting ignition from a diffusion model are also discussed.

  5. On Numerical Considerations for Modeling Reactive Astrophysical Shocks

    SciTech Connect (OSTI)

    Papatheodore, Thomas L; Messer, Bronson

    2014-01-01

    Simulating detonations in astrophysical environments is often complicated by numerical approximations to shock structure. A common prescription to ensure correct detonation speeds (and associated quantities) is to prohibit burning inside the numerically broadened shock (Fryxell et al. 1989). We have performed a series of simulations to verify the efficacy of this approximation and to understand how resolution and dimensionality might affect its use. Our results show that, in one dimension, prohibiting burning in the shock is important wherever the carbon burning length is not resolved, in keeping with the results of Fryxell et al. (1989). In two dimensions, we find that the prohibition of shock burning effectively inhibits the development of cellular structure for all but the most highly-resolved cases. We discuss the possible impacts this outcome may have on sub-grid models and detonation propagation in Type Ia supernovae.

  6. Numerical modeling of solar magnetostatic structures bounded by current sheets

    SciTech Connect (OSTI)

    Pizzo, V.J. )

    1990-12-01

    A numerical method for efficiently determining the magnetostatic equilibrium configuration of erupted solar flux concentrations, such as sunspots and flux tubes, is presented. The magnetic structures are taken to be approximately vertically oriented and axisymmetric in the surface layers and are assumed to be isolated from the surrounding photosphere by a vanishingly thin current sheet. Since the location of the current sheet is initially unknown, the final structure is generated iteratively as a free-surface problem, with the magnetic configuration for each iterate being obtained from the horizontal force balance equation, subject to the appropriate boundary conditions. Multigrid methods are used at each stage to solve the equilibrium equation, which is mapped algebraically into a body-fitted coordinate system via transfinite interpolation techniques. Several model flux tubes and sunspots are computed to illustrate the procedure, and the accuracy of the numerical method is assessed against exact analytic solutions. 32 refs.

  7. Cost | OpenEI Community

    Open Energy Info (EERE)

    Cost Home Ocop's picture Submitted by Ocop(5) Member 15 July, 2014 - 07:07 MHK LCOE Reporting Guidance Draft Cost Current DOE LCOE numerical modeling Performance Tidal Wave To...

  8. Progress report on LBL's numerical modeling studies on Cerro Prieto

    SciTech Connect (OSTI)

    Halfman-Dooley, S.E.; Lippman, M.J.; Bodvarsson, G.S.

    1989-04-01

    An exploitation model of the Cerro Prieto geothermal system is needed to assess the energy capacity of the field, estimate its productive lifetime and develop an optimal reservoir management plan. The model must consider the natural state (i.e., pre-exploitation) conditions of the system and be able to predict changes in the reservoir thermodynamic conditions (and fluid chemistry) in response to fluid production (and injection). This paper discusses the results of a three-dimensional numerical simulation of the natural state conditions of the Cerro Prieto field and compares computed and observed pressure and temperature/enthalpy changes for the 1973--1987 production period. 16 refs., 24 figs., 2 tabs.

  9. LCOE and LACE

    Gasoline and Diesel Fuel Update (EIA)

    728 386 519 519 420 341 1981-2014 Natural Gas Nonassociated, Wet After Lease Separation 215 279 468 391 332 273 1981-2014 Natural Gas Associated-Dissolved, Wet After Lease Separation 513 107 51 128 88 68 1981-2014 Dry Natural Gas 701 371 502 502 402 327 1981-2014 Natural Gas Liquids (Million Barrels) 1981 Lease Separation

    15 279 468 391 332 273 1981-2014 Adjustments -8 115 53 158 -57 -34 1981-2014 Revision Increases 45 46 34 65 59 4 1981-2014 Revision Decreases 40 64 12 209 19 41

  10. Numerical modeling of the SNS H{sup ?} ion source

    SciTech Connect (OSTI)

    Veitzer, Seth A.; Beckwith, Kristian R. C.; Kundrapu, Madhusudhan; Stoltz, Peter H.

    2015-04-08

    Ion source rf antennas that produce H- ions can fail when plasma heating causes ablation of the insulating coating due to small structural defects such as cracks. Reducing antenna failures that reduce the operating capabilities of the Spallation Neutron Source (SNS) accelerator is one of the top priorities of the SNS H- Source Program at ORNL. Numerical modeling of ion sources can provide techniques for optimizing design in order to reduce antenna failures. There are a number of difficulties in developing accurate models of rf inductive plasmas. First, a large range of spatial and temporal scales must be resolved in order to accurately capture the physics of plasma motion, including the Debye length, rf frequencies on the order of tens of MHz, simulation time scales of many hundreds of rf periods, large device sizes on tens of cm, and ion motions that are thousands of times slower than electrons. This results in large simulation domains with many computational cells for solving plasma and electromagnetic equations, short time steps, and long-duration simulations. In order to reduce the computational requirements, one can develop implicit models for both fields and particle motions (e.g. divergence-preserving ADI methods), various electrostatic models, or magnetohydrodynamic models. We have performed simulations using all three of these methods and have found that fluid models have the greatest potential for giving accurate solutions while still being fast enough to perform long timescale simulations in a reasonable amount of time. We have implemented a number of fluid models with electromagnetics using the simulation tool USim and applied them to modeling the SNS H- ion source. We found that a reduced, single-fluid MHD model with an imposed magnetic field due to the rf antenna current and the confining multi-cusp field generated increased bulk plasma velocities of > 200 m/s in the region of the antenna where ablation is often observed in the SNS source. We report

  11. Numerical Modeling of Impact Initiation of High Explosives

    SciTech Connect (OSTI)

    Wu, C J; Piggott, T; Yoh, J; Reaugh, J

    2006-05-31

    We performed continuum mechanics simulations to examine the behavior of energetic materials in Ballistic Chamber Impact (BIC) experiments, using an Arbitrary Lagrangian-Eulerian code (ALE3D). Our simulations revealed that interface friction plays an important role in inducing the formation of shear bands, which result in 'hot spots' for ignition. The temperature localization during BIC impact was found to be significant in materials with high yield strength. In those materials, there are multiple locations inside shear bands can achieve temperatures exceeding the threshold temperature for reaction. In addition, we investigated the relevant parameters influencing the pressure profile of a BIC test by numerical analysis from a simple phenomenological model. To our surprise, we found that the peaks of BIC pressure profiles not only can be a result of multi-center chemical reactions, but can also arise from factors associated apparatus configuration.

  12. An integrated analytical framework for quantifying the LCOE of waste-to-energy facilities for a range of greenhouse gas emissions policy and technical factors

    SciTech Connect (OSTI)

    Townsend, Aaron K.; Webber, Michael E.

    2012-07-15

    This study presents a novel integrated method for considering the economics of waste-to-energy (WTE) facilities with priced greenhouse gas (GHG) emissions based upon technical and economic characteristics of the WTE facility, MSW stream, landfill alternative, and GHG emissions policy. The study demonstrates use of the formulation for six different policy scenarios and explores sensitivity of the results to ranges of certain technical parameters as found in existing literature. The study shows that details of the GHG emissions regulations have large impact on the levelized cost of energy (LCOE) of WTE and that GHG regulations can either increase or decrease the LCOE of WTE depending on policy choices regarding biogenic fractions from combusted waste and emissions from landfills. Important policy considerations are the fraction of the carbon emissions that are priced (i.e. all emissions versus only non-biogenic emissions), whether emissions credits are allowed due to reducing fugitive landfill gas emissions, whether biogenic carbon sequestration in landfills is credited against landfill emissions, and the effectiveness of the landfill gas recovery system where waste would otherwise have been buried. The default landfill gas recovery system effectiveness assumed by much of the industry yields GHG offsets that are very close to the direct non-biogenic GHG emissions from a WTE facility, meaning that small changes in the recovery effectiveness cause relatively larger changes in the emissions factor of the WTE facility. Finally, the economics of WTE are dependent on the MSW stream composition, with paper and wood being advantageous, metal and glass being disadvantageous, and plastics, food, and yard waste being either advantageous or disadvantageous depending upon the avoided tipping fee and the GHG emissions price.

  13. Numerical modeling of self-limiting and self-enhancing caprock...

    Office of Scientific and Technical Information (OSTI)

    modeling of self-limiting and self-enhancing caprock alteration induced by CO2 storage in a depleted gas reservoir Citation Details In-Document Search Title: Numerical modeling of ...

  14. Numerical Modeling At Coso Geothermal Area (1997) | Open Energy...

    Open Energy Info (EERE)

    velocity structure have been estimated. It is suggested that the identification and modeling of guided waves is an effective tool to locate fracture-induced, low-velocity...

  15. Numerical Modeling At Coso Geothermal Area (2010) | Open Energy...

    Open Energy Info (EERE)

    model was developed using Poly3D to simulate the distribution and magnitude of stress concentration in the vicinity of the borehole floor, and determine the conditions...

  16. Numerical Modeling At Coso Geothermal Area (2006) | Open Energy...

    Open Energy Info (EERE)

    transport and exchange analysis Notes Finite element models of single-phase, variable-density fluid flow, conductive- convective heat transfer, fluid-rock isotope exchange, and...

  17. Numerical modeling of aerial bursts and ablation melting of Libyan...

    Office of Scientific and Technical Information (OSTI)

    Sponsoring Org: USDOE Country of Publication: United States Language: English Subject: 58 GEOSCIENCES; ABLATION; DESERTS; GLASS; MELTING; MATHEMATICAL MODELS; LIBYAN ARAB ...

  18. Numerically Solvable Model for Resonant Collisions of Electronswith Diatomic Molecules

    SciTech Connect (OSTI)

    Houfek, Karel; Rescigno, T.N.; McCurdy, C.W.

    2006-01-27

    We describe a simple model for electron-molecule collisions that has one nuclear and one electronic degree of freedom and that can be solved to arbitrarily high precision, without making the Born-Oppenheimer approximation, by employing a combination of the exterior complex scaling method and a finite-element implementation of the discrete variable representation. We compare exact cross sections for vibrational excitation and dissociative attachment with results obtained using the local complex potential approximation as commonly applied in the ''boomerang'' model, and suggest how this two-dimensional model can be used to test the underpinnings of contemporary nonlocal approximations to resonant collisions.

  19. Numerical Modeling of the Nucleation Conditions of Petal-Centerline...

    Open Energy Info (EERE)

    model using Poly3D has been developed to investigate the conditions in which the stress concentration below the floor of a borehole can cause tensile stress necessary to...

  20. Numerical Modeling At Dixie Valley Geothermal Area (Benoit, 1999...

    Open Energy Info (EERE)

    large geothermal flow test was performed where there was 6 geothermal wells flowing at once and 8 idle wells being monitored. The conceptual model developed from this flow test...

  1. Advances in the numerical modeling of field-reversed configurations

    SciTech Connect (OSTI)

    Belova, Elena V.; Davidson, Ronald C.; Ji, Hantao; Yamada, Masaaki

    2006-05-15

    The field-reversed configuration (FRC) is a compact torus with little or no toroidal magnetic field. A theoretical understanding of the observed FRC equilibrium and stability properties presents significant challenges due to the high plasma beta, plasma flows, large ion gyroradius, and the stochasticity of the particle orbits. Advanced numerical simulations are generally required to describe and understand the detailed behavior of FRC plasmas. Results of such simulations are presented in this paper. It is shown that 3D nonlinear hybrid simulations using the HYM code [E. V. Belova et al., Phys. Plasmas 7, 4996 (2000)] reproduce all major experimentally observed stability properties of elongated (theta-pinch-formed) FRCs. Namely, the scaling of the growth rate of the n=1 tilt mode with the S*/E parameter (S* is the FRC kinetic parameter, E is elongation, and n is toroidal mode number), the nonlinear saturation of the tilt mode, ion toroidal spin-up, and the growth of the n=2 rotational mode have been demonstrated and studied in detail. The HYM code has also been used to study stability properties of FRCs formed by the counterhelicity spheromak merging method. A new stability regime has been found for FRCs with elongation E{approx}1, which requires a close-fitting conducting shell and energetic beam ion stabilization.

  2. Numerical modeling of an all vanadium redox flow battery.

    SciTech Connect (OSTI)

    Clausen, Jonathan R.; Brunini, Victor E.; Moffat, Harry K.; Martinez, Mario J.

    2014-01-01

    We develop a capability to simulate reduction-oxidation (redox) flow batteries in the Sierra Multi-Mechanics code base. Specifically, we focus on all-vanadium redox flow batteries; however, the capability is general in implementation and could be adopted to other chemistries. The electrochemical and porous flow models follow those developed in the recent publication by [28]. We review the model implemented in this work and its assumptions, and we show several verification cases including a binary electrolyte, and a battery half-cell. Then, we compare our model implementation with the experimental results shown in [28], with good agreement seen. Next, a sensitivity study is conducted for the major model parameters, which is beneficial in targeting specific features of the redox flow cell for improvement. Lastly, we simulate a three-dimensional version of the flow cell to determine the impact of plenum channels on the performance of the cell. Such channels are frequently seen in experimental designs where the current collector plates are borrowed from fuel cell designs. These designs use a serpentine channel etched into a solid collector plate.

  3. Numerical simulations of altocumulus with a cloud resolving model

    SciTech Connect (OSTI)

    Liu, S.; Krueger, S.K.

    1996-04-01

    Altocumulus and altostratus clouds together cover approximately 22% of the earth`s surface. They play an important role in the earth`s energy budget through their effect on solar and infrared radiation. However, there has been little altocumulus cloud investigation by either modelers or observational programs. Starr and Cox (SC) (1985a,b) simulated an altostratus case as part of the same study in which they modeled a thin layer of cirrus. Although this calculation was originally described as representing altostratus, it probably better represents altocumulus stratiformis. In this paper, we simulate altocumulus cloud with a cloud resolving model (CRM). We simply describe the CRM first. We calculate the same middle-level cloud case as SC to compare our results with theirs. We will look at the role of cloud-scale processes in response to large-scale forcing. We will also discuss radiative effects by simulating diurnal and nocturnal cases. Finally, we discuss the utility of a 1D model by comparing 1D simulations and 2D simulations.

  4. Joint physical and numerical modeling of water distribution networks.

    SciTech Connect (OSTI)

    Zimmerman, Adam; O'Hern, Timothy John; Orear, Leslie Jr.; Kajder, Karen C.; Webb, Stephen Walter; Cappelle, Malynda A.; Khalsa, Siri Sahib; Wright, Jerome L.; Sun, Amy Cha-Tien; Chwirka, J. Benjamin; Hartenberger, Joel David; McKenna, Sean Andrew; van Bloemen Waanders, Bart Gustaaf; McGrath, Lucas K.; Ho, Clifford Kuofei

    2009-01-01

    This report summarizes the experimental and modeling effort undertaken to understand solute mixing in a water distribution network conducted during the last year of a 3-year project. The experimental effort involves measurement of extent of mixing within different configurations of pipe networks, measurement of dynamic mixing in a single mixing tank, and measurement of dynamic solute mixing in a combined network-tank configuration. High resolution analysis of turbulence mixing is carried out via high speed photography as well as 3D finite-volume based Large Eddy Simulation turbulence models. Macroscopic mixing rules based on flow momentum balance are also explored, and in some cases, implemented in EPANET. A new version EPANET code was developed to yield better mixing predictions. The impact of a storage tank on pipe mixing in a combined pipe-tank network during diurnal fill-and-drain cycles is assessed. Preliminary comparison between dynamic pilot data and EPANET-BAM is also reported.

  5. Numerical approaches to combustion modeling. Progress in Astronautics and Aeronautics. Vol. 135

    SciTech Connect (OSTI)

    Oran, E.S.; Boris, J.P. )

    1991-01-01

    Various papers on numerical approaches to combustion modeling are presented. The topics addressed include; ab initio quantum chemistry for combustion; rate coefficient calculations for combustion modeling; numerical modeling of combustion of complex hydrocarbons; combustion kinetics and sensitivity analysis computations; reduction of chemical reaction models; length scales in laminar and turbulent flames; numerical modeling of laminar diffusion flames; laminar flames in premixed gases; spectral simulations of turbulent reacting flows; vortex simulation of reacting shear flow; combustion modeling using PDF methods. Also considered are: supersonic reacting internal flow fields; studies of detonation initiation, propagation, and quenching; numerical modeling of heterogeneous detonations, deflagration-to-detonation transition to reactive granular materials; toward a microscopic theory of detonations in energetic crystals; overview of spray modeling; liquid drop behavior in dense and dilute clusters; spray combustion in idealized configurations: parallel drop streams; comparisons of deterministic and stochastic computations of drop collisions in dense sprays; ignition and flame spread across solid fuels; numerical study of pulse combustor dynamics; mathematical modeling of enclosure fires; nuclear systems.

  6. Foundation heat exchangers for residential ground source heat pump systems Numerical modeling and experimental validation

    SciTech Connect (OSTI)

    Xing, Lu; Cullin, James; Spitler, Jeffery; Im, Piljae; Fisher, Daniel

    2011-01-01

    A new type of ground heat exchanger that utilizes the excavation often made for basements or foundations has been proposed as an alternative to conventional ground heat exchangers. This article describes a numerical model that can be used to size these foundation heat exchanger (FHX) systems. The numerical model is a two-dimensional finite-volume model that considers a wide variety of factors, such as soil freezing and evapotranspiration. The FHX numerical model is validated with one year of experimental data collected at an experimental house located near Oak Ridge, Tennessee. The model shows good agreement with the experimental data-heat pump entering fluid temperatures typically within 1 C (1.8 F) - with minor discrepancies due to approximations, such as constant moisture content throughout the year, uniform evapotranspiration over the seasons, and lack of ground shading in the model.

  7. ARRAY OPTIMIZATION FOR TIDAL ENERGY EXTRACTION IN A TIDAL CHANNEL A NUMERICAL MODELING ANALYSIS

    SciTech Connect (OSTI)

    Yang, Zhaoqing; Wang, Taiping; Copping, Andrea

    2014-04-18

    This paper presents an application of a hydrodynamic model to simulate tidal energy extraction in a tidal dominated estuary in the Pacific Northwest coast. A series of numerical experiments were carried out to simulate tidal energy extraction with different turbine array configurations, including location, spacing and array size. Preliminary model results suggest that array optimization for tidal energy extraction in a real-world site is a very complex process that requires consideration of multiple factors. Numerical models can be used effectively to assist turbine siting and array arrangement in a tidal turbine farm for tidal energy extraction.

  8. Numerical Investigation of Flapwise-Torsional Vibration Model of a Smart Section Blade with Microtab

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Li, Nailu; Balas, Mark J.; Yang, Hua; Jiang, Wei; Magar, Kaman T.

    2015-01-01

    This study presents a method to develop an aeroelastic model of a smart section blade equipped with microtab. The model is suitable for potential passive vibration control study of the blade section in classic flutter. Equations of the model are described by the nondimensional flapwise and torsional vibration modes coupled with the aerodynamic model based on the Theodorsen theory and aerodynamic effects of the microtab based on the wind tunnel experimental data. The aeroelastic model is validated using numerical data available in the literature and then utilized to analyze the microtab control capability on flutter instability case and divergence instabilitymore » case. The effectiveness of the microtab is investigated with the scenarios of different output controllers and actuation deployments for both instability cases. The numerical results show that the microtab can effectively suppress both vibration modes with the appropriate choice of the output feedback controller.« less

  9. Proceedings of the Numerical Modeling for Underground Nuclear Test Monitoring Symposium

    SciTech Connect (OSTI)

    Taylor, S.R.; Kamm, J.R.

    1993-11-01

    The purpose of the meeting was to discuss the state-of-the-art in numerical simulations of nuclear explosion phenomenology with applications to test ban monitoring. We focused on the uniqueness of model fits to data, the measurement and characterization of material response models, advanced modeling techniques, and applications of modeling to monitoring problems. The second goal of the symposium was to establish a dialogue between seismologists and explosion-source code calculators. The meeting was divided into five main sessions: explosion source phenomenology, material response modeling, numerical simulations, the seismic source, and phenomenology from near source to far field. We feel the symposium reached many of its goals. Individual papers submitted at the conference are indexed separately on the data base.

  10. OpenEI Community - LCOE

    Open Energy Info (EERE)

    Cost Breakdown Structure Draft http:en.openei.orgcommunitydocumentmhk-cost-breakdown-structure-draft

    The generalized Cost Breakdown Structure (CBS) for marine and...

  11. NUMERICAL VERIFICATION OF THE RELAP-7 CORE CHANNEL SINGLE-PHASE MODEL

    SciTech Connect (OSTI)

    Haihua Zhao; Ling Zou; Hongbin Zhang; Richard Martineau

    2014-06-01

    The RELAP-7 code is the next generation of nuclear reactor system safety analysis code being developed at the Idaho National Laboratory (INL). All the physics in RELAP-7 are fully coupled and the errors resulted from the traditional operator-splitting approach are eliminated. By using 2nd order methods in both time and space and eliminating operator-splitting errors, the numerical error of RELAP-7 can be minimized. Numerical verification is the process to verify the orders of numerical methods. It is an important part of modern verification and validation process. The core channel component in RELAP-7 is designed to simulate coolant flow as well as the conjugated heat transfer between coolant flow and the fuel rod. A special treatment at fuel centerline to avoid numerical singularity for the cylindrical heat conduction in the continuous finite element mesh is discussed. One steady state test case and one fast power up transient test case are utilized for the verification of the core channel model with single-phase flow. Analytical solution for the fuel pin temperature and figures of merit such as peak clad temperature and peak fuel temperature are used to define numerical errors. These cases prove that the mass and energy are well conserved and 2nd order convergence rates for both time and space are achieved in the core channel model.

  12. On Improving Analytical Models of Cosmic Reionization for Matching Numerical Simulations

    SciTech Connect (OSTI)

    Kaurov, Alexander A.

    2016-01-01

    The methods for studying the epoch of cosmic reionization vary from full radiative transfer simulations to purely analytical models. While numerical approaches are computationally expensive and are not suitable for generating many mock catalogs, analytical methods are based on assumptions and approximations. We explore the interconnection between both methods. First, we ask how the analytical framework of excursion set formalism can be used for statistical analysis of numerical simulations and visual representation of the morphology of ionization fronts. Second, we explore the methods of training the analytical model on a given numerical simulation. We present a new code which emerged from this study. Its main application is to match the analytical model with a numerical simulation. Then, it allows one to generate mock reionization catalogs with volumes exceeding the original simulation quickly and computationally inexpensively, meanwhile reproducing large scale statistical properties. These mock catalogs are particularly useful for CMB polarization and 21cm experiments, where large volumes are required to simulate the observed signal.

  13. CHEMICAL TRANSPORT IN A FISSURED BOCK: VERIFICATION OF A NUMERICAL MODEL

    SciTech Connect (OSTI)

    Rasmuson, A.; Narasimhan, T.N.; Neretnieks, I.

    1982-04-01

    Numerical models for simulating chemical transport in fissured rocks constitute powerful tools for evaluating the acceptability of geological nuclear waste repositories. Due to the very long-term, high toxicity of some nuclear waste products, the models are required to predict, in certain cases, the spatial and temporal distribution of chemical concentration less than 0.001% of the concentration released from the repository. Whether numerical models can provide such accuracies is a major question addressed in the present work. To this end, we have verified a numerical model, TRUMP, which solves the advective diffusion equation in general three dimensions with or without decay and source terms. The method is based on an integrated finite-difference approach. The model was verified against known analytic solution of the one-dimensional advection-diffusion problem as well as the problem of advection-diffusion in a system of parallel fractures separated by spherical particles. The studies show that as long as the magnitude of advectance is equal to or less than that of conductance for the closed surface bounding any volume element in the region (that is, numerical Peclet number <2), the numerical method can indeed match the analytic solution within errors of ±10{sup -3} % or less. The realistic input parameters used in the sample calculations suggest that such a range of Peclet numbers is indeed likely to characterize deep groundwater systems in granitic and ancient argillaceous systems. Thus TRUMP in its present form does provide a viable tool for use in nuclear waste evaluation studies. A sensitivity analysis based on the analytic solution suggests that the errors in prediction introduced due to uncertainties in input parameters is likely to be larger than the computational inaccuracies introduced by the numerical model. Currently, a disadvantage in the TRUMP model is that the iterative method of solving the set of simultaneous equations is rather slow when time

  14. 2D numerical simulation of the MEP energy-transport model with a finite difference scheme

    SciTech Connect (OSTI)

    Romano, V. . E-mail: romano@dmi.unict.it

    2007-02-10

    A finite difference scheme of Scharfetter-Gummel type is used to simulate a consistent energy-transport model for electron transport in semiconductors devices, free of any fitting parameters, formulated on the basis of the maximum entropy principle. Simulations of silicon n{sup +}-n-n{sup +} diodes, 2D-MESFET and 2D-MOSFET and comparisons with the results obtained by a direct simulation of the Boltzmann transport equation and with other energy-transport models, known in the literature, show the validity of the model and the robustness of the numerical scheme.

  15. Recent numerical and algorithmic advances within the volume tracking framework for modeling interfacial flows

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    François, Marianne M.

    2015-05-28

    A review of recent advances made in numerical methods and algorithms within the volume tracking framework is presented. The volume tracking method, also known as the volume-of-fluid method has become an established numerical approach to model and simulate interfacial flows. Its advantage is its strict mass conservation. However, because the interface is not explicitly tracked but captured via the material volume fraction on a fixed mesh, accurate estimation of the interface position, its geometric properties and modeling of interfacial physics in the volume tracking framework remain difficult. Several improvements have been made over the last decade to address these challenges.more » In this study, the multimaterial interface reconstruction method via power diagram, curvature estimation via heights and mean values and the balanced-force algorithm for surface tension are highlighted.« less

  16. A Numerical Model of the Temperature Field of the Cast and Solidified Ceramic Material

    SciTech Connect (OSTI)

    Kavicka, Frantisek; Sekanina, Bohumil; Stransky, Karel; Stetina, Josef [Brno University of Technology, Brno, Technicka 2 (Czech Republic); Dobrovska, Jana [Technical University of Ostrava, Ostrava, Tr. 17.listopadu 17 (Czech Republic)

    2010-06-15

    Corundo-baddeleyit material (CBM)--EUCOR--is a heat- and wear-resistant material even at extreme temperatures. This article introduces a numerical model of solidification and cooling of this material in a non-metallic mould. The model is capable of determining the total solidification time of the casting and also the place of the casting which solidifies last. Furthermore, it is possible to calculate the temperature gradient in any point and time, and also determine the local solidification time and the solidification interval of any point. The local solidification time is one of the input parameters for the cooperating model of chemical heterogeneity. This second model and its application on samples of EUCOR prove that the applied method of measurement of chemical heterogeneity provides detailed quantitative information on the material structure and makes it possible to analyse the solidification process. The analysis of this process entails statistical processing of the results of the measurements of the heterogeneity of the components of EUCOR and performs correlation of individual components during solidification. The crystallisation process seems to be very complicated, where the macro- and microscopic segregations differ significantly. The verification of both numerical models was conducted on a real cast 350x200x400 mm block.

  17. The shallow water equations as a hybrid flow model for the numerical and experimental analysis of hydro power stations

    SciTech Connect (OSTI)

    Ostermann, Lars; Seidel, Christian

    2015-03-10

    The numerical analysis of hydro power stations is an important method of the hydraulic design and is used for the development and optimisation of hydro power stations in addition to the experiments with the physical submodel of a full model in the hydraulic laboratory. For the numerical analysis, 2D and 3D models are appropriate and commonly used.The 2D models refer mainly to the shallow water equations (SWE), since for this flow model a large experience on a wide field of applications for the flow analysis of numerous problems in hydraulic engineering already exists. Often, the flow model is verified by in situ measurements. In order to consider 3D flow phenomena close to singularities like weirs, hydro power stations etc. the development of a hybrid fluid model is advantageous to improve the quality and significance of the global model. Here, an extended hybrid flow model based on the principle of the SWE is presented. The hybrid flow model directly links the numerical model with the experimental data, which may originate from physical full models, physical submodels and in-situ measurements. Hence a wide field of application of the hybrid model emerges including the improvement of numerical models and the strong coupling of numerical and experimental analysis.

  18. Numerical Modeling of the Lake Mary Road Bridge for Foundation Reuse Assessment

    SciTech Connect (OSTI)

    Sitek, M. A.; Bojanowski, C.; Lottes, S. A.

    2015-04-01

    This project uses numerical techniques to assess the structural integrity and capacity of the bridge foundations and, as a result, reduces the risk associated with reusing the same foundation for a new superstructure. Nondestructive test methods of different types were used in combination with the numerical modeling and analysis. The onsite tests included visual inspection, tomography, ground penetrating radar, drilling boreholes and coreholes, and the laboratory tests on recovered samples. The results were utilized to identify the current geometry of the structure with foundation, including the hidden geometry of the abutments and piers, and soil and foundation material properties. This data was used to build the numerical models and run computational analyses on a high performance computer cluster to assess the structural integrity of the bridge and foundations including the suitability of the foundation for reuse with a new superstructure and traffic that will increase the load on the foundations. Computational analysis is more cost-effective and gives an advantage of getting more detailed knowledge about the structural response. It also enables to go beyond non-destructive testing and find the failure conditions without destroying the structure under consideration.

  19. Numerical models for static denting and dynamic puncture of gas transmission linepipe and their validation

    SciTech Connect (OSTI)

    Zarea, M.F.; Toumbas, D.N.; Philibert, C.E.; Deo, I.

    1996-12-31

    Gas transmission pipe resistance to external damage is a subject of great attention at Gaz de France and in Europe. Existing results cover part of the necessary criteria for the residual life of damaged pipelines, but more knowledge is needed on defect creation. The authors propose to complement existing experimental work which is limited to the explored range of parameters by validated numerical models. The first, simple static denting model aims at optimizing the conditions for calculating the residual stress distribution needed to assess the fatigue life of dents and dents and gouges. The second, more complex dynamic puncture model calculates both the puncture force and the puncture energy for a given pipe, excavator and tooth geometry. These models can contribute to enhance the external damage prevention policies of transmission pipeline operators.

  20. Monitoring and Numerical Modeling of Shallow CO{sub 2} Injection, Greene County, Missouri

    SciTech Connect (OSTI)

    Rovey, Charles; Gouzie, Douglas; Biagioni, Richard

    2013-09-30

    The project titled Monitoring and Numerical Modeling of Shallow CO{sub 2} Injection, Greene County, Missouri provided training for three graduate students in areas related to carbon capture and storage. Numerical modeling of CO{sub 2} injection into the St. Francois aquifer at the Southwest Power Plant Site in Greene County, Missouri indicates that up to 4.1 x 10{sup 5} metric tons of CO{sub 2} per year could be injected for 30 years without exceeding a 3 MPa differential injection pressure. The injected CO{sub 2} would remain sequestered below the top of the overlying caprock (St. Francois confining unit) for more than 1000 years. Geochemical modeling indicates that portions of the injected CO{sub 2} will react rapidly with trace minerals in the aquifer to form various solid carbonate mineral phases. These minerals would store significant portions of injected CO{sub 2} over geologic time scales. Finally, a GIS data base on the pore-fluid chemistry of the overlying aquifer system in Missouri, the Ozark aquifer, was compiled from many sources. This data base could become useful in monitoring for leakage from future CO{sub 2} sequestration sites.

  1. Numeric-modeling sensitivity analysis of the performance of wind turbine arrays

    SciTech Connect (OSTI)

    Lissaman, P.B.S.; Gyatt, G.W.; Zalay, A.D.

    1982-06-01

    An evaluation of the numerical model created by Lissaman for predicting the performance of wind turbine arrays has been made. Model predictions of the wake parameters have been compared with both full-scale and wind tunnel measurements. Only limited, full-scale data were available, while wind tunnel studies showed difficulties in representing real meteorological conditions. Nevertheless, several modifications and additions have been made to the model using both theoretical and empirical techniques and the new model shows good correlation with experiment. The larger wake growth rate and shorter near wake length predicted by the new model lead to reduced interference effects on downstream turbines and hence greater array efficiencies. The array model has also been re-examined and now incorporates the ability to show the effects of real meteorological conditions such as variations in wind speed and unsteady winds. The resulting computer code has been run to show the sensitivity of array performance to meteorological, machine, and array parameters. Ambient turbulence and windwise spacing are shown to dominate, while hub height ratio is seen to be relatively unimportant. Finally, a detailed analysis of the Goodnoe Hills wind farm in Washington has been made to show how power output can be expected to vary with ambient turbulence, wind speed, and wind direction.

  2. Solar Advisor Model; Session: Modeling and Analysis (Presentation)

    SciTech Connect (OSTI)

    Blair, N.

    2008-04-01

    This project supports the Solar America Initiative by: (1) providing a consistent framework for analyzing and comparing power system costs and performance across the range of solar technologies and markets, PV, solar heat systems, CSP, residential, commercial and utility markets; (2) developing and validating performance models to enable accurate calculation of levelized cost of energy (LCOE); (3) providing a consistent modeling platform for all TPP's; and (4) supporting implementation and usage of cost models.

  3. Numerical modeling of the action of an explosion on an iron slab

    SciTech Connect (OSTI)

    Sugak, S.G.; Fortov, V.E.; Kanel', G.I.; Ni, A.L.; Stel'makh, V.G.

    1983-09-01

    This article examines the explosion of a condensed high explosive (HE) by a flat Armco iron slab. The fundamental physical processes accompanying the nonstationary two-dimensional compression pulse propagation in a metal are analyzed, a mathematical model of the process is formulated, and its numerical characteristics are determined. A cylindrical charge of trinitrotoluene of 20 mm diameter and height, initiated at the axial point on the free surface, was mounted on a 10-mm-thick and 120-mm-diameter armco iron disk. A finite-difference method using a quadrangular computational mesh is employed. It is assumed that the process of detonation transformation occurs instantaneously and is characterized only by singularities in the equation of state. Numerical modeling of the loading of an iron slab by an explosion showed that the simple kinetic relationship proposed for the fracture computation will assure a realistic description of the fracture process. It is determined that the progress of a reversible polymorphic transformation in the specimen material substantially influences the nature of its fracture by explosion. The magnitude of the strain on the strength properties of the material must be taken into account in order to achieve a detailed description of the fracture.

  4. Numerical modeling of carrier gas flow in atomic layer deposition vacuum reactor: A comparative study of lattice Boltzmann models

    SciTech Connect (OSTI)

    Pan, Dongqing; Chien Jen, Tien [Department of Mechanical Engineering, University of Wisconsin-Milwaukee, Milwaukee, Wisconsin 53201 (United States); Li, Tao [School of Mechanical Engineering, Dalian University of Technology, Dalian 116024 (China); Yuan, Chris, E-mail: cyuan@uwm.edu [Department of Mechanical Engineering, University of Wisconsin-Milwaukee, 3200 North Cramer Street, Milwaukee, Wisconsin 53211 (United States)

    2014-01-15

    This paper characterizes the carrier gas flow in the atomic layer deposition (ALD) vacuum reactor by introducing Lattice Boltzmann Method (LBM) to the ALD simulation through a comparative study of two LBM models. Numerical models of gas flow are constructed and implemented in two-dimensional geometry based on lattice BhatnagarGrossKrook (LBGK)-D2Q9 model and two-relaxation-time (TRT) model. Both incompressible and compressible scenarios are simulated and the two models are compared in the aspects of flow features, stability, and efficiency. Our simulation outcome reveals that, for our specific ALD vacuum reactor, TRT model generates better steady laminar flow features all over the domain with better stability and reliability than LBGK-D2Q9 model especially when considering the compressible effects of the gas flow. The LBM-TRT is verified indirectly by comparing the numerical result with conventional continuum-based computational fluid dynamics solvers, and it shows very good agreement with these conventional methods. The velocity field of carrier gas flow through ALD vacuum reactor was characterized by LBM-TRT model finally. The flow in ALD is in a laminar steady state with velocity concentrated at the corners and around the wafer. The effects of flow fields on precursor distributions, surface absorptions, and surface reactions are discussed in detail. Steady and evenly distributed velocity field contribute to higher precursor concentration near the wafer and relatively lower particle velocities help to achieve better surface adsorption and deposition. The ALD reactor geometry needs to be considered carefully if a steady and laminar flow field around the wafer and better surface deposition are desired.

  5. Full-Scale Numerical Modeling of Turbulent Processes in the Earth's Ionosphere

    SciTech Connect (OSTI)

    Eliasson, B.; Stenflo, L.; Shukla, P. K.

    2008-10-15

    We present a full-scale simulation study of ionospheric turbulence by means of a generalized Zakharov model based on the separation of variables into high-frequency and slow time scales. The model includes realistic length scales of the ionospheric profile and of the electromagnetic and electrostatic fields, and uses ionospheric plasma parameters relevant for high-latitude radio facilities such as Eiscat and HAARP. A nested grid numerical method has been developed to resolve the different length-scales, while avoiding severe restrictions on the time step. The simulation demonstrates the parametric decay of the ordinary mode into Langmuir and ion-acoustic waves, followed by a Langmuir wave collapse and short-scale caviton formation, as observed in ionospheric heating experiments.

  6. Precarious Rock Methodology for Seismic Hazard: Physical Testing, Numerical Modeling and Coherence Studies

    SciTech Connect (OSTI)

    Anooshehpoor, Rasool; Purvance, Matthew D.; Brune, James N.; Preston, Leiph A.; Anderson, John G.; Smith, Kenneth D.

    2006-09-29

    This report covers the following projects: Shake table tests of precarious rock methodology, field tests of precarious rocks at Yucca Mountain and comparison of the results with PSHA predictions, study of the coherence of the wave field in the ESF, and a limited survey of precarious rocks south of the proposed repository footprint. A series of shake table experiments have been carried out at the University of Nevada, Reno Large Scale Structures Laboratory. The bulk of the experiments involved scaling acceleration time histories (uniaxial forcing) from 0.1g to the point where the objects on the shake table overturned a specified number of times. The results of these experiments have been compared with numerical overturning predictions. Numerical predictions for toppling of large objects with simple contact conditions (e.g., I-beams with sharp basal edges) agree well with shake-table results. The numerical model slightly underpredicts the overturning of small rectangular blocks. It overpredicts the overturning PGA for asymmetric granite boulders with complex basal contact conditions. In general the results confirm the approximate predictions of previous studies. Field testing of several rocks at Yucca Mountain has approximately confirmed the preliminary results from previous studies, suggesting that he PSHA predictions are too high, possibly because the uncertainty in the mean of the attenuation relations. Study of the coherence of wavefields in the ESF has provided results which will be very important in design of the canisters distribution, in particular a preliminary estimate of the wavelengths at which the wavefields become incoherent. No evidence was found for extreme focusing by lens-like inhomogeneities. A limited survey for precarious rocks confirmed that they extend south of the repository, and one of these has been field tested.

  7. Sensitivity of injection costs to input petrophysical parameters in numerical geologic carbon sequestration models

    SciTech Connect (OSTI)

    Cheng, C. L.; Gragg, M. J.; Perfect, E.; White, Mark D.; Lemiszki, P. J.; McKay, L. D.

    2013-08-24

    Numerical simulations are widely used in feasibility studies for geologic carbon sequestration. Accurate estimates of petrophysical parameters are needed as inputs for these simulations. However, relatively few experimental values are available for CO2-brine systems. Hence, a sensitivity analysis was performed using the STOMP numerical code for supercritical CO2 injected into a model confined deep saline aquifer. The intrinsic permeability, porosity, pore compressibility, and capillary pressure-saturation/relative permeability parameters (residual liquid saturation, residual gas saturation, and van Genuchten alpha and m values) were varied independently. Their influence on CO2 injection rates and costs were determined and the parameters were ranked based on normalized coefficients of variation. The simulations resulted in differences of up to tens of millions of dollars over the life of the project (i.e., the time taken to inject 10.8 million metric tons of CO2). The two most influential parameters were the intrinsic permeability and the van Genuchten m value. Two other parameters, the residual gas saturation and the residual liquid saturation, ranked above the porosity. These results highlight the need for accurate estimates of capillary pressure-saturation/relative permeability parameters for geologic carbon sequestration simulations in addition to measurements of porosity and intrinsic permeability.

  8. A Numerical Algorithm for the Solution of a Phase-Field Model of Polycrystalline Materials

    SciTech Connect (OSTI)

    Dorr, M R; Fattebert, J; Wickett, M E; Belak, J F; Turchi, P A

    2008-12-04

    We describe an algorithm for the numerical solution of a phase-field model (PFM) of microstructure evolution in polycrystalline materials. The PFM system of equations includes a local order parameter, a quaternion representation of local orientation and a species composition parameter. The algorithm is based on the implicit integration of a semidiscretization of the PFM system using a backward difference formula (BDF) temporal discretization combined with a Newton-Krylov algorithm to solve the nonlinear system at each time step. The BDF algorithm is combined with a coordinate projection method to maintain quaternion unit length, which is related to an important solution invariant. A key element of the Newton-Krylov algorithm is the selection of a preconditioner to accelerate the convergence of the Generalized Minimum Residual algorithm used to solve the Jacobian linear system in each Newton step. Results are presented for the application of the algorithm to 2D and 3D examples.

  9. Large Eddy Simulation of Wind Turbine Wakes. Detailed Comparisons of Two Codes Focusing on Effects of Numerics and Subgrid Modeling

    SciTech Connect (OSTI)

    Martinez-Tossas, Luis A.; Churchfield, Matthew J.; Meneveau, Charles

    2015-06-18

    In this work we report on results from a detailed comparative numerical study from two Large Eddy Simulation (LES) codes using the Actuator Line Model (ALM). The study focuses on prediction of wind turbine wakes and their breakdown when subject to uniform inflow. Previous studies have shown relative insensitivity to subgrid modeling in the context of a finite-volume code. The present study uses the low dissipation pseudo-spectral LES code from Johns Hopkins University (LESGO) and the second-order, finite-volume OpenFOAMcode (SOWFA) from the National Renewable Energy Laboratory. When subject to uniform inflow, the loads on the blades are found to be unaffected by subgrid models or numerics, as expected. The turbulence in the wake and the location of transition to a turbulent state are affected by the subgrid-scale model and the numerics.

  10. Large Eddy Simulation of Wind Turbine Wakes. Detailed Comparisons of Two Codes Focusing on Effects of Numerics and Subgrid Modeling

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Martinez-Tossas, Luis A.; Churchfield, Matthew J.; Meneveau, Charles

    2015-06-18

    In this work we report on results from a detailed comparative numerical study from two Large Eddy Simulation (LES) codes using the Actuator Line Model (ALM). The study focuses on prediction of wind turbine wakes and their breakdown when subject to uniform inflow. Previous studies have shown relative insensitivity to subgrid modeling in the context of a finite-volume code. The present study uses the low dissipation pseudo-spectral LES code from Johns Hopkins University (LESGO) and the second-order, finite-volume OpenFOAMcode (SOWFA) from the National Renewable Energy Laboratory. When subject to uniform inflow, the loads on the blades are found to bemore » unaffected by subgrid models or numerics, as expected. The turbulence in the wake and the location of transition to a turbulent state are affected by the subgrid-scale model and the numerics.« less

  11. NUMERICAL MODELING OF CONTAMINANT TRANSPORT IN FRACTURED POROUS MEDIA USING MIXED FINITE ELEMENT AND FINITE VOLUME METHODS

    SciTech Connect (OSTI)

    Taylor, G.; Dong, C.; Sun, S.

    2010-03-18

    A mathematical model for contaminant species passing through fractured porous media is presented. In the numerical model, we combine two locally conservative methods, i.e. mixed finite element (MFE) and the finite volume methods. Adaptive triangle mesh is used for effective treatment of the fractures. A hybrid MFE method is employed to provide an accurate approximation of velocities field for both the fractures and matrix which are crucial to the convection part of the transport equation. The finite volume method and the standard MFE method are used to approximate the convection and dispersion terms respectively. The model is used to investigate the interaction of adsorption with transport and to extract information on effective adsorption distribution coefficients. Numerical examples in different fractured media illustrate the robustness and efficiency of the proposed numerical model.

  12. Numerical Modeling of Geomechanical Processes Related to CO{sub 2} Injection within Generic Reservoirs

    SciTech Connect (OSTI)

    Eckert, Andreas

    2013-05-31

    In this project generic anticline structures have been used for numerical modeling analyses to study the influence of geometrical parameters, fluid flow boundary conditions, in situ stress regime and inter-bedding friction coefficient on geomechanical risks such as fracture reactivation and fracture generation. The resulting stress states for these structures are also used to determine safe drilling directions and a methodology for wellbore trajection optimization is developed that is applicable for non-Andersonian stress states. The results of the fluid flow simulation show that the type of fluid flow boundary condition is of utmost importance and has significant impact on all injection related parameters. It is recommended that further research is conducted to establish a method to quantify the fluid flow boundary conditions for injection applications. The results of the geomechanical simulation show that in situ stress regime is a crucial, if not the most important, factor determining geomechanical risks. For extension and strike slip stress regimes anticline structures should be favored over horizontally layered basin as they feature higher ΔP{sub c} magnitudes. If sedimentary basins are tectonically relaxed and their state of stress is characterized by the uni-axial strain model the basin is in exact frictional equilibrium and fluids should not be injected. The results also show that low inter bedding friction coefficients effectively decouple layers resulting in lower ΔP{sub c} magnitudes, especially for the compressional stress regime.

  13. NUMERICAL MODELING OF THE COAGULATION AND POROSITY EVOLUTION OF DUST AGGREGATES

    SciTech Connect (OSTI)

    Okuzumi, Satoshi; Sakagami, Masa-aki [Graduate School of Human and Environmental Studies, Kyoto University, Yoshida-nihonmatsu-cho, Sakyo-ku, Kyoto 606-8501 (Japan); Tanaka, Hidekazu, E-mail: satoshi.okuzumi@ax2.ecs.kyoto-u.ac.j [Institute of Low Temperature Science, Hokkaido University, Sapporo 060-0819 (Japan)

    2009-12-20

    Porosity evolution of dust aggregates is crucial in understanding dust evolution in protoplanetary disks. In this study, we present useful tools to study the coagulation and porosity evolution of dust aggregates. First, we present a new numerical method for simulating dust coagulation and porosity evolution as an extension of the conventional Smoluchowski equation. This method follows the evolution of the mean porosity for each aggregate mass simultaneously with the evolution of the mass distribution function. This method reproduces the results of previous Monte Carlo simulations with much less computational expense. Second, we propose a new collision model for porous dust aggregates on the basis of our N-body experiments on aggregate collisions. As the first step, we focus on 'hit-and-stick' collisions, which involve neither compression nor fragmentation of aggregates. We first obtain empirical data on porosity changes between the classical limits of ballistic cluster-cluster and particle-cluster aggregation. Using the data, we construct a recipe for the porosity change due to general hit-and-stick collisions as well as formulae for the aerodynamical and collisional cross sections. Our collision model is thus more realistic than a previous model of Ormel et al. based on the classical aggregation limits only. Simple coagulation simulations using the extended Smoluchowski method show that our collision model explains the fractal dimensions of porous aggregates observed in a full N-body simulation and a laboratory experiment. By contrast, similar simulations using the collision model of Ormel et al. result in much less porous aggregates, meaning that this model underestimates the porosity increase upon unequal-sized collisions. Besides, we discover that aggregates at the high-mass end of the distribution can have a considerably small aerodynamical cross section per unit mass compared with aggregates of lower masses. This occurs when aggregates drift under uniform

  14. Numerical method to test a theoretical model of the quantum interferen...

    Office of Scientific and Technical Information (OSTI)

    A numerical method is provided to fit the experimental conductivity to the complicated conductivity expression for the quantum interference effect of Anderson localization. This ...

  15. Numerical modeling of multi-GeV laser wakefield electron acceleration inside a dielectric capillary tube

    SciTech Connect (OSTI)

    Paradkar, B. S.; Cros, B.; Maynard, G.; Mora, P.

    2013-08-15

    Numerical modeling of laser wakefield electron acceleration inside a gas filled dielectric capillary tube is presented. Guiding of a short pulse laser inside a dielectric capillary tube over a long distance (∼1 m) and acceleration of an externally injected electron bunch to ultra-relativistic energies (∼5-10 GeV) are demonstrated in the quasi-linear regime of laser wakefield acceleration. Two dimensional axisymmetric simulations were performed with the code WAKE-EP (Extended Performances), which allows computationally efficient simulations of such long scale plasma. The code is an upgrade of the quasi-static particle code, WAKE [P. Mora and T. M. Antonsen, Jr., Phys. Plasmas 4, 217 (1997)], to simulate the acceleration of an externally injected electron bunch (including beam loading effect) and propagation of the laser beam inside a dielectric capillary. The influence of the transverse electric field of the plasma wake on the radial loss of the accelerated electrons to the dielectric wall is investigated. The stable acceleration of electrons to multi-GeV energy with a non-resonant laser pulse with a large spot-size is demonstrated.

  16. GETEM -Geothermal Electricity Technology Evaluation Model

    Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site

    1 GETEM -Geothermal Electricity Technology Evaluation Model Background: GETEM was originally developed for the Department of Energy's Geothermal Technologies Program to provide both a method for quantifying the power generation cost from geothermal energy, and a means of assessing how technology advances might impact those generation costs. Generation cost is determined as the Levelized-Cost-of-Electricity (LCOE). The model is intended to provide representative estimates of cost and performance

  17. Numerical research of the optimal control problem in the semi-Markov inventory model

    SciTech Connect (OSTI)

    Gorshenin, Andrey K.

    2015-03-10

    This paper is devoted to the numerical simulation of stochastic system for inventory management products using controlled semi-Markov process. The results of a special software for the systems research and finding the optimal control are presented.

  18. HELIOSPHERIC PROPAGATION OF CORONAL MASS EJECTIONS: COMPARISON OF NUMERICAL WSA-ENLIL+CONE MODEL AND ANALYTICAL DRAG-BASED MODEL

    SciTech Connect (OSTI)

    Vrnak, B.; ic, T.; Dumbovi?, M.; Temmer, M.; Mstl, C.; Veronig, A. M.; Taktakishvili, A.; Mays, M. L.; Odstr?il, D. E-mail: tzic@geof.hr E-mail: manuela.temmer@uni-graz.at E-mail: astrid.veronig@uni-graz.at E-mail: m.leila.mays@nasa.gov

    2014-08-01

    Real-time forecasting of the arrival of coronal mass ejections (CMEs) at Earth, based on remote solar observations, is one of the central issues of space-weather research. In this paper, we compare arrival-time predictions calculated applying the numerical ''WSA-ENLIL+Cone model'' and the analytical ''drag-based model'' (DBM). Both models use coronagraphic observations of CMEs as input data, thus providing an early space-weather forecast two to four days before the arrival of the disturbance at the Earth, depending on the CME speed. It is shown that both methods give very similar results if the drag parameter ? = 0.1 is used in DBM in combination with a background solar-wind speed of w = 400 km s{sup 1}. For this combination, the mean value of the difference between arrival times calculated by ENLIL and DBM is ?-bar =0.099.0 hr with an average of the absolute-value differences of |?|-bar =7.1 hr. Comparing the observed arrivals (O) with the calculated ones (C) for ENLIL gives O C = 0.3 16.9 hr and, analogously, O C = +1.1 19.1 hr for DBM. Applying ? = 0.2 with w = 450 km s{sup 1} in DBM, one finds O C = 1.7 18.3 hr, with an average of the absolute-value differences of 14.8 hr, which is similar to that for ENLIL, 14.1 hr. Finally, we demonstrate that the prediction accuracy significantly degrades with increasing solar activity.

  19. A Numerical Model Without Truncation Error for a Steady-State Analysis of a Once-Through Steam Generator

    SciTech Connect (OSTI)

    Sim, Yoon Sub; Kim, Eui Kwang; Eoh, Jae Hyuk [Korea Atomic Energy Research Institute (Korea, Republic of)

    2005-06-15

    To overcome the drawbacks of conventional schemes for a numerical analysis of a steam generator (SG), an efficient numerical model has been developed to analyze the steady state of a once-through-type SG where the feedwater is heated to superheated steam. In the developed model, the temperature and enthalpy are defined at the boundary of a calculation cell, and the exact solutions for the temperature distribution in a calculation cell are utilized. This feature of the developed model frees calculation from the undesirable effects of numerical diffusion, and only a small number of nodes are required. Also, the developed model removes the ambiguity from the parameter values at the inlet and exit of a calculation.The BoSupSG-SS computer code was developed by using the analysis model, and it performed well with only three calculation nodes to analyze a superheated SG. The developed model can be effectively used for the cases where a fast one-dimensional calculation is required such as an SG or system design analysis.

  20. Numerical modeling of the groundwater contaminant transport for the Lake Karachai Area: The methodological approach and the basic two- dimensional regional model

    SciTech Connect (OSTI)

    Petrov, A.V.; Samsonova, L.M.; Vasil`kova, N.A.; Zinin, A.I.; Zinina, G.A. |

    1994-06-01

    Methodological aspects of the numerical modeling of the groundwater contaminant transport for the Lake Karachay area are discussed. Main features of conditions of the task are the high grade of non-uniformity of the aquifer in the fractured rock massif and the high density of the waste solutions, and also the high volume of the input data: both on the part of parameters of the aquifer (number of pump tests) and on the part of observations of functions of processes (long-time observations by the monitoring well grid). The modeling process for constructing the two dimensional regional model is described, and this model is presented as the basic model for subsequent full three-dimensional modeling in sub-areas of interest. Original powerful mathematical apparatus and computer codes for finite-difference numerical modeling are used.

  1. Coupling a Mesoscale Numerical Weather Prediction Model with Large-Eddy Simulation for Realistic Wind Plant Aerodynamics Simulations (Poster)

    SciTech Connect (OSTI)

    Draxl, C.; Churchfield, M.; Mirocha, J.; Lee, S.; Lundquist, J.; Michalakes, J.; Moriarty, P.; Purkayastha, A.; Sprague, M.; Vanderwende, B.

    2014-06-01

    Wind plant aerodynamics are influenced by a combination of microscale and mesoscale phenomena. Incorporating mesoscale atmospheric forcing (e.g., diurnal cycles and frontal passages) into wind plant simulations can lead to a more accurate representation of microscale flows, aerodynamics, and wind turbine/plant performance. Our goal is to couple a numerical weather prediction model that can represent mesoscale flow [specifically the Weather Research and Forecasting model] with a microscale LES model (OpenFOAM) that can predict microscale turbulence and wake losses.

  2. Experimental and numerical investigations of beryllium strength models using the Rayleigh-Taylor instability

    SciTech Connect (OSTI)

    Henry de Frahan, M. T.; Belof, J. L.; Cavallo, R. M.; Raevsky, V. A.; Ignatova, O. N.; Lebedev, A.; Ancheta, D. S.; El-dasher, B. S.; Florando, J. N.; Gallegos, G. F.; Johnsen, E.; LeBlanc, M. M.

    2015-06-14

    A recent collaboration between LLNL and VNIIEF has produced a set of high explosive driven Rayleigh-Taylor strength data for beryllium. Design simulations using legacy strength models from Steinberg-Lund and Preston-Tonks-Wallace (PTW) suggested an optimal design that would delineate between not just different strength models, but different parameters sets of the PTW model. Application of the models to the post-shot results, however, shows close to classical growth. We characterize the material properties of the beryllium tested in the experiments. We also discuss recent efforts to simulate the data using the legacy strength models as well as the more recent RING relaxation model developed at VNIIEF. Finally, we present shock and ramp-loading recovery experiments conducted as part of the collaboration.

  3. Numerical modeling of pulsed laser-material interaction and of laser plume dynamics

    SciTech Connect (OSTI)

    Zhao, Qiang; Shi, Yina

    2015-03-10

    We have developed two-dimensional Arbitrary Lagrangian Eulerian (ALE) code which is used to study the physical processes, the plasma absorption, the crater profile, and the temperature distribution on metallic target and below the surface. The ALE method overcomes problems with Lagrangian moving mesh distortion by mesh smoothing and conservative quantities remapping from Lagrangian mesh to smoothed one. A new second order accurate diffusion solver has been implemented for the thermal conduction and radiation transport on distorted mesh. The results of numerical simulation of pulsed laser ablation are presented. The influences of different processes, such as time evolution of the surface temperature, interspecies interactions (elastic collisions, recombination-dissociation reaction), interaction with an ambient gas are examined. The study presents particular interest for the analysis of experimental results obtained during pulsed laser ablation.

  4. Numerical modeling of the effects of wave energy converter characteristics on nearshore wave conditions

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Chang, G.; Ruehl, K.; Jones, C. A.; Roberts, J.; Chartrand, C.

    2015-12-24

    Modeled nearshore wave propagation was investigated downstream of simulated wave energy converters (WECs) to evaluate overall near- and far-field effects of WEC arrays. Model sensitivity to WEC characteristics and WEC array deployment scenarios was evaluated using a modified version of an industry standard wave model, Simulating WAves Nearshore (SWAN), which allows the incorporation of device-specific WEC characteristics to specify obstacle transmission. The sensitivity study illustrated that WEC device type and subsequently its size directly resulted in wave height variations in the lee of the WEC array. Wave heights decreased up to 30% between modeled scenarios with and without WECs formore » large arrays (100 devices) of relatively sizable devices (26 m in diameter) with peak power generation near to the modeled incident wave height. Other WEC types resulted in less than 15% differences in modeled wave height with and without WECs, with lesser influence for WECs less than 10 m in diameter. Wave directions and periods were largely insensitive to changes in parameters. Furthermore, additional model parameterization and analysis are required to fully explore the model sensitivity of peak wave period and mean wave direction to the varying of the parameters.« less

  5. Numerical modeling of the effects of wave energy converter characteristics on nearshore wave conditions

    SciTech Connect (OSTI)

    Chang, G.; Ruehl, K.; Jones, C. A.; Roberts, J.; Chartrand, C.

    2015-12-24

    Modeled nearshore wave propagation was investigated downstream of simulated wave energy converters (WECs) to evaluate overall near- and far-field effects of WEC arrays. Model sensitivity to WEC characteristics and WEC array deployment scenarios was evaluated using a modified version of an industry standard wave model, Simulating WAves Nearshore (SWAN), which allows the incorporation of device-specific WEC characteristics to specify obstacle transmission. The sensitivity study illustrated that WEC device type and subsequently its size directly resulted in wave height variations in the lee of the WEC array. Wave heights decreased up to 30% between modeled scenarios with and without WECs for large arrays (100 devices) of relatively sizable devices (26 m in diameter) with peak power generation near to the modeled incident wave height. Other WEC types resulted in less than 15% differences in modeled wave height with and without WECs, with lesser influence for WECs less than 10 m in diameter. Wave directions and periods were largely insensitive to changes in parameters. Furthermore, additional model parameterization and analysis are required to fully explore the model sensitivity of peak wave period and mean wave direction to the varying of the parameters.

  6. A numerical investigation of transonic axial compressor rotor flow using a low-Reynolds-number {kappa}-{epsilon} turbulence model

    SciTech Connect (OSTI)

    Arima, T.; Sonoda, T.; Shirotori, M.; Tamura, A.; Kikuchi, K.

    1999-01-01

    The authors have developed a computer simulation code for three-dimensional viscous flow in turbomachinery based on the time-averaged compressible Navier-Stokes equations and a low-Reynolds-number {kappa}-{epsilon} turbulence model. It is described in detail in this paper. The code is used to compute the flow fields for two types of rotor (a transonic fan NASA Rotor 67 and a transonic axial compressor NASA rotor 37), and numerical results are compared to experimental data based on aerodynamic probe and laser anemometer measurements. In the case of Rotor 67, calculated and experimental results are compared under the design speed to validate the code. The calculated results show good agreement with the experimental data, such as the rotor performance map and the spanwise distribution of total pressure, total temperature, and flow angle downstream of the rotor. In the case of Rotor 37, detailed comparisons between the numerical results and the experimental data are made under the design speed condition to assess the overall quality of the numerical solution. Furthermore, comparisons under the part-speed condition are used to investigate a flow field without passage shock. The results are well predicted qualitatively. However, considerable quantitative discrepancies remain in predicting the flow near the tip. In order to assess the predictive capabilities of the developed code, computed flow structures are presented with the experimental data for each rotor and the cause of the discrepancies is discussed.

  7. A numerical test of KPZ scaling; Potts models coupled to two-dimensional quantum gravity

    SciTech Connect (OSTI)

    Baillie, C.F. , 91 - Nozay ); Johnston, D.A. . Dept. of Mathematics)

    1992-06-07

    In this paper, the authors perform Monte-Carlo simulations using the Wolff cluster algorithm of the q = 2 (Ising), 3, 4 and q = 10 Potts models on dynamical phi-cubed graphs of spherical topology with up to 5000 nodes. The authors find that the measured critical exponents are in reasonable agreement with those from the exact solution of the Ising model and with those calculated from KPZ scaling for q = 3, 4 where no exact solution is available. Using Binder's cumulant, the authors find that the q = 10 Potts model displays a first order phase transition on a dynamical graph, as it does on a fixed lattice. The authors also examine the internal geometry of the graphs generated in the simulation, finding a linear relationship between ring length probabilities and the central charge of the Potts model.

  8. The value of numerical modeling in understanding the complete load/deformation behavior of cohesive-frictional materials

    SciTech Connect (OSTI)

    Potyondy, D.O.; Fairhurst, C.E.

    1999-07-01

    The post-peak load/deformation behavior of cohesive-frictional materials is an integral part of the overall response of a specimen to compressive loading. A more comprehensive understanding of the pre- and post-peak behavior is necessary. Recent developments in numerical modeling that allow study of the overall response of a synthetic material containing discrete heterogeneities and discontinuities both at the micro (particle) scale and at the larger scale of jointed rock masses can greatly aid the interpretation and application of laboratory test results on these materials.

  9. Final Report: A Model Management System for Numerical Simulations of Subsurface Processes

    SciTech Connect (OSTI)

    Zachmann, David

    2013-10-07

    The DOE and several other Federal agencies have committed significant resources to support the development of a large number of mathematical models for studying subsurface science problems such as groundwater flow, fate of contaminants and carbon sequestration, to mention only a few. This project provides new tools to help decision makers and stakeholders in subsurface science related problems to select an appropriate set of simulation models for a given field application.

  10. Numerical modeling of immiscible two-phase flow in micro-models using a commercial CFD code

    SciTech Connect (OSTI)

    Crandall, Dustin; Ahmadia, Goodarz; Smith, Duane H.

    2009-01-01

    Off-the-shelf CFD software is being used to analyze everything from flow over airplanes to lab-on-a-chip designs. So, how accurately can two-phase immiscible flow be modeled flowing through some small-scale models of porous media? We evaluate the capability of the CFD code FLUENT{trademark} to model immiscible flow in micro-scale, bench-top stereolithography models. By comparing the flow results to experimental models we show that accurate 3D modeling is possible.

  11. Synthesis of Numerical Methods for Modeling Wave Energy Converter-Point Absorbers: Preprint

    SciTech Connect (OSTI)

    Li, Y.; Yu, Y. H.

    2012-05-01

    During the past few decades, wave energy has received significant attention among all ocean energy formats. Industry has proposed hundreds of prototypes such as an oscillating water column, a point absorber, an overtopping system, and a bottom-hinged system. In particular, many researchers have focused on modeling the floating-point absorber as the technology to extract wave energy. Several modeling methods have been used such as the analytical method, the boundary-integral equation method, the Navier-Stokes equations method, and the empirical method. However, no standardized method has been decided. To assist the development of wave energy conversion technologies, this report reviews the methods for modeling the floating-point absorber.

  12. Numerical model for the vacuum pyrolysis of scrap tires in batch reactors

    SciTech Connect (OSTI)

    Yang, J.; Tanguy, P.A.; Roy, C.

    1995-06-01

    A quantitative model for scrap tire pyrolysis in a batch scale reactor developed comprises the following basic phenomena: conduction inside tire particles; conduction, convection, and radiation between the feedstock particles or between the fluids and the particles; tire pyrolysis reaction; exothermicity and endothermicity caused by tire decomposition and volatilization; and the variation of the composition and the thermal properties of tire particles. This model was used to predict the transient temperature and density distributions in the bed of particles, the volatile product evolution rate, the mass change, the energy consumption during the pyrolysis process, and the pressure history in a tire pyrolysis reactor with a load of 1 kg. The model predictions agree well with independent experimental data.

  13. Deep geological isolation of nuclear waste: numerical modeling of repository scale hydrology

    SciTech Connect (OSTI)

    Dettinger, M.D.

    1980-04-01

    The Scope of Work undertaken covers three main tasks, described as follows: (Task 1) CDM provided consulting services to the University on modeling aspects of the study having to do with transport processes involving the local groundwater system near the repository and the flow of fluids and vapors through the various porous media making up the repository system. (Task 2) CDM reviewed literature related to repository design, concentrating on effects of the repository geometry, location and other design factors on the flow of fluids within the repository boundaries, drainage from the repository structure, and the eventual transport of radionucldies away from the repository site. (Task 3) CDM, in a joint effort with LLL personnel, identified generic boundary and initial conditions, identified processes to be modeled, and recommended a modeling approach with suggestions for appropriate simplifications and approximations to the problem and identifiying important parameters necessary to model the processes. This report consists of two chapters and an appendix. The first chapter (Chapter III of the LLL report) presents a detailed description and discussion of the modeling approach developed in this project, its merits and weaknesses, and a brief review of the difficulties anticipated in implementing the approach. The second chapter (Chapter IV of the LLL report) presents a summary of a survey of researchers in the field of repository performance analysis and a discussion of that survey in light of the proposed modeling approach. The appendix is a review of the important physical processes involved in the potential hydrologic transport of radionuclides through, around and away from deep geologic nuclear waste repositories.

  14. Numerical modeling and experiments of creep crack growth under cyclic loading

    SciTech Connect (OSTI)

    Brust, F.W.

    1995-12-31

    This paper presents a summary of some recent studies of creep crack growth under history dependent load conditions. The effect of a proper constitutive law is illustrated. Moreover, the asymptotic fields are reconsidered under cyclic creep conditions. In addition, several experiments are modeled and the behavior of integral parameters is discussed.

  15. Numerical Study of Coal Gasification Using Eulerian-Eulerian Multiphase Model

    SciTech Connect (OSTI)

    Shi, S.; Guenther, C.; Orsino, S.

    2007-09-01

    Gasification converts the carbon-containing material into a synthesis gas (syngas) which can be used as a fuel to generate electricity or used as a basic chemical building block for a large number of uses in the petrochemical and refining industries. Based on the mode of conveyance of the fuel and the gasifying medium, gasification can be classified into fixed or moving bed, fluidized bed, and entrained flow reactors. Entrained flow gasifiers normally feature dilute flow with small particle size and can be successfully modeled with the Discrete Phase Method (DPM). For the other types, the Eulerian-Eulerian (E-E) or the so called two-fluid multiphase model is a more appropriate approach. The E-E model treats the solid phase as a distinct interpenetrating granular fluid and it is the most general-purposed multi-fluid model. This approach provides transient, three-dimensional, detailed information inside the reactor which would otherwise be unobtainable through experiments due to the large scale, high pressure and/or temperature. In this paper, a transient, three-dimensional model of the Power Systems Development Facility (PSDF) transport gasifier will be presented to illustrate how Computational Fluid Dynamics (CFD) can be used for large-scale complicated geometry with detailed physics and chemistry. In the model, eleven species are included in the gas phase while four pseudo-species are assumed in the solid phase. A total of sixteen reactions, both homogeneous (involving only gas phase species) and heterogeneous (involving species in both gas and solid phases), are used to model the coal gasification chemistry. Computational results have been validated against PSDF experimental data from lignite to bituminous coals under both air and oxygen blown conditions. The PSDF gasifier geometry was meshed with about 70,000, hexahedra-dominated cells. A total of six cases with different coal, feed gas, and/or operation conditions have been performed. The predicted and

  16. Modeling and numerical simulation of oscillatory two-phase flows, with application to boiling water nuclear reactors

    SciTech Connect (OSTI)

    Rosa, M.P.; Podowski, M.Z.

    1995-09-01

    This paper is concerned with the analysis of dynamics and stability of boiling channels and systems. The specific objectives are two-fold. One of them is to present the results of a study aimed at analyzing the effects of various modeling concepts and numerical approaches on the transient response and stability of parallel boiling channels. The other objective is to investigate the effect of closed-loop feedback on stability of a boiling water reactor (BWR). Various modeling and computational issues for parallel boiling channels are discussed, such as: the impact of the numerical discretization scheme for the node containing the moving boiling boundary on the convergence and accuracy of computations, and the effects of subcooled boiling and other two-phase flow phenomena on the predictions of marginal stability conditions. Furthermore, the effects are analyzed of local loss coefficients around the recirculation loop of a boiling water reactor on stability of the reactor system. An apparent paradox is explained concerning the impact of changing single-phase losses on loop stability. The calculations have been performed using the DYNOBOSS computer code. The results of DYNOBOSS validation against other computer codes and experimental data are shown.

  17. One-dimensional numerical fluid dynamics model of the spreading of liquefied gaseous fuel (LGF) on water

    SciTech Connect (OSTI)

    Stein, W.; Ermak, D.L.

    1980-11-04

    A computer model has been developed to simulate the spreading of an evaporating liquefied gaseous fuel that has been spilled on the surface of a denser liquid. This would correspond, for example, to the spreading of liquefied natural gas spilled onto water. The model is based on the one-dimensional, time-dependent equations of conservation of mass and momentum, with the assumption that the pool of liquid fuel spreads in a radially symmetric manner. It includes the effects of vaporization, shear at the fuel-liquid interface, and buoyancy due to the density difference between the fuel and the liquid onto which it is spilled. Both instantaneous and continuous spills of finite volume are treated. The height and spreading velocity of the pool of spilled fuel are calculated as functions of time and radius by numerically solving the conservation equations with a finite difference method. Output of the calculations is presented in both tabular and graphical form.

  18. One-dimensional numerical fluid dynamics model of the spreading of liquefied gaseous fuel (LGF) on water

    SciTech Connect (OSTI)

    Stein, W.; Ermak, D.L.

    1981-01-01

    A computer model has been developed to simulate the spreading of an evaporating liquefied gaseous fuel that has been spilled on the surface of a denser liquid. This would correspond, for example, to the spreading of liquefied natural gas spilled onto water. The model is based on the one-dimensional, time-dependent equations of conservation of mass and momentum, with the assumption that the pool of liquid fuel spreads in a radially symmetric manner. It includes the effects of vaporization, shear at the fuel-liquid interface, and buoyancy due to the density difference between the fuel and the liquid onto which it is spilled. Both instantaneous and continuous spills of finite volume are treated. The height and spreading velocity of the pool of spilled fuel are calculated as functions of time and radius by numerically solving the conservation equations with a finite difference method.Output of the calculations is presented in both tabular and graphical form.

  19. MEAN-FIELD MODELING OF AN α{sup 2} DYNAMO COUPLED WITH DIRECT NUMERICAL SIMULATIONS OF RIGIDLY ROTATING CONVECTION

    SciTech Connect (OSTI)

    Masada, Youhei; Sano, Takayoshi E-mail: sano@ile.osaka-u.ac.jp

    2014-10-10

    The mechanism of large-scale dynamos in rigidly rotating stratified convection is explored by direct numerical simulations (DNS) in Cartesian geometry. A mean-field dynamo model is also constructed using turbulent velocity profiles consistently extracted from the corresponding DNS results. By quantitative comparison between the DNS and our mean-field model, it is demonstrated that the oscillatory α{sup 2} dynamo wave, excited and sustained in the convection zone, is responsible for large-scale magnetic activities such as cyclic polarity reversal and spatiotemporal migration. The results provide strong evidence that a nonuniformity of the α-effect, which is a natural outcome of rotating stratified convection, can be an important prerequisite for large-scale stellar dynamos, even without the Ω-effect.

  20. A direct numerical simulation-based investigation and modeling of pressure Hessian effects on compressible velocity gradient dynamics

    SciTech Connect (OSTI)

    Danish, Mohammad Suman, Sawan Srinivasan, Balaji

    2014-12-15

    The pressure Hessian tensor plays a key role in shaping the behavior of the velocity gradient tensor, and in turn, that of many incumbent non-linear processes in a turbulent flow field. In compressible flows, the role of pressure Hessian is even more important because it represents the level of fluid-thermodynamic coupling existing in the flow field. In this work, we first perform a direct numerical simulation-based study to clearly identify, isolate, and understand various important inviscid mechanisms that govern the evolution of the pressure Hessian tensor in compressible turbulence. The ensuing understanding is then employed to introduce major improvements to the existing Lagrangian model of the pressure Hessian tensor (the enhanced Homogenized Euler equation or EHEE) in terms of (i) non-symmetric, non-isentropic effects and (ii) improved representation of the anisotropic portion of the pressure Hessian tensor. Finally, we evaluate the new model extensively by comparing the new model results against known turbulence behavior over a range of Reynolds and Mach numbers. Indeed, the new model shows much improved performance as compared to the EHEE model.

  1. Investigation of the effect of chemistry models on the numerical predictions of the supersonic combustion of hydrogen

    SciTech Connect (OSTI)

    Kumaran, K.; Babu, V.

    2009-04-15

    In this numerical study, the influence of chemistry models on the predictions of supersonic combustion in a model combustor is investigated. To this end, 3D, compressible, turbulent, reacting flow calculations with a detailed chemistry model (with 37 reactions and 9 species) and the Spalart-Allmaras turbulence model have been carried out. These results are compared with earlier results obtained using single step chemistry. Hydrogen is used as the fuel and three fuel injection schemes, namely, strut, staged (i.e., strut and wall) and wall injection, are considered to evaluate the impact of the chemistry models on the flow field predictions. Predictions of the mass fractions of major species, minor species, dimensionless stagnation temperature, dimensionless static pressure rise and thrust percentage along the combustor length are presented and discussed. Overall performance metrics such as mixing efficiency and combustion efficiency are used to draw inferences on the nature (whether mixing- or kinetic-controlled) and the completeness of the combustion process. The predicted values of the dimensionless wall static pressure are compared with experimental data reported in the literature. The calculations show that multi step chemistry predicts higher and more wide spread heat release than what is predicted by single step chemistry. In addition, it is also shown that multi step chemistry predicts intricate details of the combustion process such as the ignition distance and induction distance. (author)

  2. Numerical modeling of laser tunneling ionization in explicit particle-in-cell codes

    SciTech Connect (OSTI)

    Chen, M.; Cormier-Michel, E.; Geddes, C.G.R.; Bruhwiler, D.L.; Yu, L.L.; Esarey, E.; Schroeder, C.B.; Leemans, W.P.

    2013-03-01

    Methods for the calculation of laser tunneling ionization in explicit particle-in-cell codes used for modeling laserplasma interactions are compared and validated against theoretical predictions. Improved accuracy is obtained by using the direct current form for the ionization rate. Multi level ionization in a single time step and energy conservation have been considered during the ionization process. The effects of grid resolution and number of macro-particles per cell are examined. Implementation of the ionization algorithm in two different particle-in-cell codes is compared for the case of ionization-based electron injection in a laserplasma accelerator.

  3. Application of numerical modeling in a clean-coal demonstration project

    SciTech Connect (OSTI)

    Latham, C.E.; Laursen, T.A.; Bellanca, C.; Duong, H.

    1992-11-01

    Currently, utility boilers equipped with cell burners comprise 13% of pre-NSPS coal-fired generating capacity. The cell burner rapidly mixes the pulverized coal and combustion air resulting in rapid combustion and high NO{sub x} generation. A US Department of Energy (DOE) Clean-Coal Technology Demonstration project is underway at Dayton Power & Light`s J. M. Stuart Station to demonstrate the Low-NO{sub x} Cell{trademark} burner (LNCB{trademark}) on a 605-MWe utility boiler originally equipped with cell burners. The LNCB{trademark} is designed to reduce NO{sub x} emissions by delaying the mixing of the coal and the combustion air without boiler pressure part modifications. Preliminary post-retrofit testing results showed unexpectedly high carbon monoxide (CO) and hydrogen sulfide (H{sub 2}S) concentrations below the lowest burner row. The substoichiometric operation of the lowest burner row caused the relatively high concentrations in the lower furnace. Babcock & Wilcox`s flow, combustion, and heat transfer models were used to predict the CO concentrations in the lower furnace. The predictions were compared to field measurements for three different operating conditions. Based on this validation, the models were used to evaluate several methods for mitigating the CO concentrations. The results of this analysis are presented and discussed. The most attractive alternative was selected and will be implemented during the spring of 1992. The effectiveness of the new alternative will be available when the installation is complete and testing resumes.

  4. Numerical comparison between a gyrofluid and gyrokinetic model investigating collisionless magnetic reconnection

    SciTech Connect (OSTI)

    Zacharias, O.; Kleiber, R.; Borchardt, M.; Comisso, L.; Grasso, D.; Hatzky, R.

    2014-06-15

    The first detailed comparison between gyrokinetic and gyrofluid simulations of collisionless magnetic reconnection has been carried out. Both the linear and nonlinear evolution of the collisionless tearing mode have been analyzed. In the linear regime, we have found a good agreement between the two approaches over the whole spectrum of linearly unstable wave numbers, both in the drift kinetic limit and for finite ion temperature. Nonlinearly, focusing on the small-Δ′ regime, with Δ′ indicating the standard tearing stability parameter, we have compared relevant observables such as the evolution and saturation of the island width, as well as the island oscillation frequency in the saturated phase. The results are basically the same, with small discrepancies only in the value of the saturated island width for moderately high values of Δ′. Therefore, in the regimes investigated here, the gyrofluid approach can describe the collisionless reconnection process as well as the more complete gyrokinetic model.

  5. Mixture Preparation and Nitric Oxide Formation in a GDI Engine studied by Combined Laser Diagnostics and Numerical Modeling

    SciTech Connect (OSTI)

    Volker Sick; Dennis N. Assanis

    2002-11-27

    Through the combination of advanced imaging laser diagnostics with multi-dimensional computer models, a new understanding of the performance of direct-injection gasoline engines is pursuit. The work focuses on the fuel injection process, the breakup of the liquid into a fine spray and the mixing of the fuel with the in-cylinder gases. Non-intrusive laser diagnostics will be used to measure the spatial distribution of droplets and vaporized fuel with very high temporal resolution. These data along with temperature measurements will be used to validate a new spray breakup model for gasoline direct-injection. Experimental data on near wall fuel distributions will be used for comparison with a model that predicts the spray-wall interaction and the dynamics of the liquid film on the surface. Quantitative measurements of local nitric oxide concentrations inside the combustion chamber will provide a critical test for a numerical simulation of the nitric oxide formation process. This model is based on a modified flamelet approach and will be used to study the effects of exhaust gas recirculation.

  6. An improved neutral diffusion model and numerical solution of the two dimensional edge plasma fluid equations. Final report

    SciTech Connect (OSTI)

    Prinja, A.K.

    1998-09-01

    relatively smooth as a consequence of the less localized recycling, leading to an improved convergence rate of the numerical algorithm. Peak plasma density is lower and the temperature correspondingly higher than those predicted by the standard diffusion model. It is believed that the FFCD model is more accurate. With both the TP continuation and multigrid methods, the author has demonstrated the robustness of these two methods. A mutually beneficial hybridization between the TP method and multigrid methods is clearly an alternative for edge plasma simulation. While the fundamental transport model considered in this work has ignored important physics such as drifts and currents, he has nevertheless demonstrated the versatility and robustness of the numerical scheme to handle such new physics. The application of gaseous-radiative divertor model in this work is just a beginning and up to this point numerically, the future is exciting.

  7. Tall Tower Wind Energy Monitoring and Numerical Model Validation in Northern Nevada

    SciTech Connect (OSTI)

    Koracin, D.; Kaplan, M.; Smith, C.; McCurdy, G.; Wolf, A.; McCord, T.; King, K.; Belu, R.; Horvath, K.

    2015-10-01

    The main objectives of this project were to conduct a tall-tower and sodar field campaign in complex terrain, investigate wind properties relevant to wind energy assessment, and evaluate high-resolution models with fixed and adaptive grid structures. Two 60-m towers at Virginia Peak ridges near Washoe Valley, Nevada, were instrumented with cup and vane anemometers as well as sonic anemometers, and an acoustic sounder (hereafter sodar) was installed near one of the towers. The towers were located 2,700 m apart with a vertical distance of 140 m elevation between their bases. Each tower had a downhill exposure of rolling complex terrain, with the nearby valley floor 3,200 m to the west and 800 m below the summit. Cup anemometers were installed at both towers at 20, 40, and 60 m, wind vanes at 20 and 60 m, and sonic anemometers at 20 and 60 m. The sodar measurements were nominally provided every 10 m in vertical distance from 40 to 200 m with the quality of the data generally decreasing with height. Surface air temperature, atmospheric pressure, and radiation measurements were conducted at 1.5 m AGL at both of the towers. Although the plan was to conduct a 1-year period of data collection, we extended the period (October 5, 2012 through February 24, 2014) to cover for possible data loss from instrument or communication problems. We also present a preliminary analysis of the towers and sodar data, including a detailed inventory of available and missing data as well as outliers. The analysis additionally includes calculation of the Weibull parameters, turbulence intensity, and initial computation of wind power density at various heights.

  8. ON THE IMPACT OF SUPER RESOLUTION WSR-88D DOPPLER RADAR DATA ASSIMILATION ON HIGH RESOLUTION NUMERICAL MODEL FORECASTS

    SciTech Connect (OSTI)

    Chiswell, S

    2009-01-11

    Assimilation of radar velocity and precipitation fields into high-resolution model simulations can improve precipitation forecasts with decreased 'spin-up' time and improve short-term simulation of boundary layer winds (Benjamin, 2004 & 2007; Xiao, 2008) which is critical to improving plume transport forecasts. Accurate description of wind and turbulence fields is essential to useful atmospheric transport and dispersion results, and any improvement in the accuracy of these fields will make consequence assessment more valuable during both routine operation as well as potential emergency situations. During 2008, the United States National Weather Service (NWS) radars implemented a significant upgrade which increased the real-time level II data resolution to 8 times their previous 'legacy' resolution, from 1 km range gate and 1.0 degree azimuthal resolution to 'super resolution' 250 m range gate and 0.5 degree azimuthal resolution (Fig 1). These radar observations provide reflectivity, velocity and returned power spectra measurements at a range of up to 300 km (460 km for reflectivity) at a frequency of 4-5 minutes and yield up to 13.5 million point observations per level in super-resolution mode. The migration of National Weather Service (NWS) WSR-88D radars to super resolution is expected to improve warning lead times by detecting small scale features sooner with increased reliability; however, current operational mesoscale model domains utilize grid spacing several times larger than the legacy data resolution, and therefore the added resolution of radar data is not fully exploited. The assimilation of super resolution reflectivity and velocity data into high resolution numerical weather model forecasts where grid spacing is comparable to the radar data resolution is investigated here to determine the impact of the improved data resolution on model predictions.

  9. Canyon Disposal Initiative - Numerical Modeling of Contaminant Transport from Grouted Residual Waste in the 221-U Facility (U Plant)

    SciTech Connect (OSTI)

    Rockhold, Mark L.; White, Mark D.; Freeman, Eugene J.

    2004-10-12

    This letter report documents initial numerical analyses conducted by PNNL to provide support for a feasibility study on decommissioning of the canyon buildings at Hanford. The 221-U facility is the first of the major canyon buildings to be decommissioned. The specific objective of this modeling effort was to provide estimates of potential rates of migration of residual contaminants out of the 221-U facility during the first 40 years after decommissioning. If minimal contaminant migration is predicted to occur from the facility during this time period, then the structure may be deemed to provide a level of groundwater protection that is essentially equivalent to the liner and leachate collection systems that are required at conventional landfills. The STOMP code was used to simulate transport of selected radionuclides out of a canyon building, representative of the 221-U facility after decommissioning, for a period of 40 years. Simulation results indicate that none of the selected radionuclides that were modeled migrated beyond the concrete structure of the facility during the 40-year period of interest. Jacques (2001) identified other potential contaminants in the 221-U facility that were not modeled, however, including kerosene, phenol, and various metals. Modeling of these contaminants was beyond the scope of this preliminary effort due to increased complexity. Simulation results indicate that contaminant release from the canyon buildings will be diffusion controlled at early times. Advection is expected to become much more important at later times, after contaminants have diffused out of the facility and into the surrounding soil environment. After contaminants have diffused out of the facility, surface infiltration covers will become very important for mitigating further transport of contaminants in the underlying vadose zone and groundwater.

  10. The deposition and burning characteristics during slagging co-firing coal and wood: modeling and numerical simulation

    SciTech Connect (OSTI)

    Wang, X.H.; Zhao, D.Q.; Jiang, L.Q.; Yang, W.B.

    2009-07-01

    Numerical analysis was used to study the deposition and burning characteristics of combining co-combustion with slagging combustion technologies in this paper. The pyrolysis and burning kinetic models of different fuels were implanted into the WBSF-PCC2 (wall burning and slag flow in pulverized co-combustion) computation code, and then the slagging and co-combustion characteristics (especially the wall burning mechanism of different solid fuels and their effects on the whole burning behavior in the cylindrical combustor at different mixing ratios under the condition of keeping the heat input same) were simulated numerically. The results showed that adding wood powder at 25% mass fraction can increase the temperature at the initial stage of combustion, which is helpful to utilize the front space of the combustor. Adding wood powder at a 25% mass fraction can increase the reaction rate at the initial combustion stage; also, the coal ignitability is improved, and the burnout efficiency is enhanced by about 5% of suspension and deposition particles, which is helpful for coal particles to burn entirely and for combustion devices to minimize their dimensions or sizes. The results also showed that adding wood powder at a proper ratio is helpful to keep the combustion stability, not only because of the enhancement for the burning characteristics, but also because the running slag layer structure can be changed more continuously, which is very important for avoiding the abnormal slag accumulation in the slagging combustor. The theoretic analysis in this paper proves that unification of co-combustion and slagging combustion technologies is feasible, though more comprehensive and rigorous research is needed.

  11. A numerical study of short residence time FCC riser flows with a new flow/kinetics modeling technique.

    SciTech Connect (OSTI)

    Chang, S. L.

    1998-08-25

    Fluid Catalytic Cracking (FCC) technology is the most important process used by the refinery industry to convert crude oil to valuable lighter products such as gasoline. New and modified processes are constantly developed by refinery companies to improve their global competitiveness and meet more stringent environmental regulations. Short residence time FCC riser reactor is one of the advanced processes that the refining industry is actively pursuing because it can improve the yield selectivity and efficiency of an FCC unit. However, as the residence time becomes shorter, the impact of the mixing between catalyst and feed oil at the feed injection region on the product yield becomes more significant. Currently, most FCC computer models used by the refineries perform sophisticated kinetic calculations on simplified flow field and can not be used to evaluate the impact of fluid mixing on the performance of an FCC unit. Argonne National Laboratory (AFL) is developing a computational fluid dynamic (CFD) code ICRKFLO for FCC riser flow modeling. The code, employing hybrid hydrodynamic-chemical kinetic coupling techniques, is used to investigate the effect of operating and design conditions on the product yields of FCC riser reactors. Numerical calculations were made using the code to examine the impacts of the operating and design conditions on the product yields. The controlling parameters under investigation include the residence time, reaction temperature, and catalyst/oil ratio. This paper describes the CFD code, presents computation results, and discusses the effects of operating conditions on the performance of short residence time FCC riser reactors.

  12. Reference Model 6 (RM6): Oscillating Wave Energy Converter.

    SciTech Connect (OSTI)

    Bull, Diana L; Smith, Chris; Jenne, Dale Scott; Jacob, Paul; Copping, Andrea; Willits, Steve; Fontaine, Arnold; Brefort, Dorian; Gordon, Margaret Ellen; Copeland, Robert; Jepsen, Richard A.

    2014-10-01

    This report is an addendum to SAND2013-9040: Methodology for Design and Economic Analysis of Marine Energy Conversion (MEC) Technologies. This report describes an Oscillating Water Column Wave Energy Converter reference model design in a complementary manner to Reference Models 1-4 contained in the above report. In this report, a conceptual design for an Oscillating Water Column Wave Energy Converter (WEC) device appropriate for the modeled reference resource site was identified, and a detailed backward bent duct buoy (BBDB) device design was developed using a combination of numerical modeling tools and scaled physical models. Our team used the methodology in SAND2013-9040 for the economic analysis that included costs for designing, manufacturing, deploying, and operating commercial-scale MEC arrays, up to 100 devices. The methodology was applied to identify key cost drivers and to estimate levelized cost of energy (LCOE) for this RM6 Oscillating Water Column device in dollars per kilowatt-hour (%24/kWh). Although many costs were difficult to estimate at this time due to the lack of operational experience, the main contribution of this work was to disseminate a detailed set of methodologies and models that allow for an initial cost analysis of this emerging technology. This project is sponsored by the U.S. Department of Energy's (DOE) Wind and Water Power Technologies Program Office (WWPTO), within the Office of Energy Efficiency & Renewable Energy (EERE). Sandia National Laboratories, the lead in this effort, collaborated with partners from National Laboratories, industry, and universities to design and test this reference model.

  13. Solutions of the two-dimensional Hubbard model: Benchmarks and results from a wide range of numerical algorithms

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    LeBlanc, J. P. F.; Antipov, Andrey E.; Becca, Federico; Bulik, Ireneusz W.; Chan, Garnet Kin-Lic; Chung, Chia -Min; Deng, Youjin; Ferrero, Michel; Henderson, Thomas M.; Jiménez-Hoyos, Carlos A.; et al

    2015-12-14

    Numerical results for ground-state and excited-state properties (energies, double occupancies, and Matsubara-axis self-energies) of the single-orbital Hubbard model on a two-dimensional square lattice are presented, in order to provide an assessment of our ability to compute accurate results in the thermodynamic limit. Many methods are employed, including auxiliary-field quantum Monte Carlo, bare and bold-line diagrammatic Monte Carlo, method of dual fermions, density matrix embedding theory, density matrix renormalization group, dynamical cluster approximation, diffusion Monte Carlo within a fixed-node approximation, unrestricted coupled cluster theory, and multireference projected Hartree-Fock methods. Comparison of results obtained by different methods allows for the identification ofmore » uncertainties and systematic errors. The importance of extrapolation to converged thermodynamic-limit values is emphasized. Furthermore, cases where agreement between different methods is obtained establish benchmark results that may be useful in the validation of new approaches and the improvement of existing methods.« less

  14. Solutions of the two-dimensional Hubbard model: Benchmarks and results from a wide range of numerical algorithms

    SciTech Connect (OSTI)

    LeBlanc, J. P. F.; Antipov, Andrey E.; Becca, Federico; Bulik, Ireneusz W.; Chan, Garnet Kin-Lic; Chung, Chia -Min; Deng, Youjin; Ferrero, Michel; Henderson, Thomas M.; Jiménez-Hoyos, Carlos A.; Kozik, E.; Liu, Xuan -Wen; Millis, Andrew J.; Prokof’ev, N. V.; Qin, Mingpu; Scuseria, Gustavo E.; Shi, Hao; Svistunov, B. V.; Tocchio, Luca F.; Tupitsyn, I. S.; White, Steven R.; Zhang, Shiwei; Zheng, Bo -Xiao; Zhu, Zhenyue; Gull, Emanuel

    2015-12-14

    Numerical results for ground-state and excited-state properties (energies, double occupancies, and Matsubara-axis self-energies) of the single-orbital Hubbard model on a two-dimensional square lattice are presented, in order to provide an assessment of our ability to compute accurate results in the thermodynamic limit. Many methods are employed, including auxiliary-field quantum Monte Carlo, bare and bold-line diagrammatic Monte Carlo, method of dual fermions, density matrix embedding theory, density matrix renormalization group, dynamical cluster approximation, diffusion Monte Carlo within a fixed-node approximation, unrestricted coupled cluster theory, and multireference projected Hartree-Fock methods. Comparison of results obtained by different methods allows for the identification of uncertainties and systematic errors. The importance of extrapolation to converged thermodynamic-limit values is emphasized. Furthermore, cases where agreement between different methods is obtained establish benchmark results that may be useful in the validation of new approaches and the improvement of existing methods.

  15. Thermal effects of groundwater flow through subarctic fens: A case study based on field observations and numerical modeling

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Sjöberg, Ylva; Coon, Ethan; K. Sannel, A. Britta; Pannetier, Romain; Harp, Dylan; Frampton, Andrew; Painter, Scott L.; Lyon, Steve W.

    2016-02-05

    Modeling and observation of ground temperature dynamics are the main tools for understanding current permafrost thermal regimes and projecting future thaw. Until recently, most studies on permafrost have focused on vertical ground heat fluxes. Groundwater can transport heat in both lateral and vertical directions but its influence on ground temperatures at local scales in permafrost environments is not well understood. In this paper, we combine field observations from a subarctic fen in the sporadic permafrost zone with numerical simulations of coupled water and thermal fluxes. At the Tavvavuoma study site in northern Sweden, ground temperature profiles and groundwater levels weremore » observed in boreholes. These observations were used to set up one- and two-dimensional simulations down to 2 m depth across a gradient of permafrost conditions within and surrounding the fen. Two-dimensional scenarios representing the fen under various hydraulic gradients were developed to quantify the influence of groundwater flow on ground temperature. Our observations suggest that lateral groundwater flow significantly affects ground temperatures. This is corroborated by modeling results that show seasonal ground ice melts 1 month earlier when a lateral groundwater flux is present. Further, although the thermal regime may be dominated by vertically conducted heat fluxes during most of the year, isolated high groundwater flow rate events such as the spring freshet are potentially important for ground temperatures. Finally, as sporadic permafrost environments often contain substantial portions of unfrozen ground with active groundwater flow paths, knowledge of this heat transport mechanism is important for understanding permafrost dynamics in these environments.« less

  16. A mixture-energy-consistent six-equation two-phase numerical model for fluids with interfaces, cavitation and evaporation waves

    SciTech Connect (OSTI)

    Pelanti, Marica; Shyue, Keh-Ming

    2014-02-15

    We model liquidgas flows with cavitation by a variant of the six-equation single-velocity two-phase model with stiff mechanical relaxation of SaurelPetitpasBerry (Saurel et al., 2009) [9]. In our approach we employ phasic total energy equations instead of the phasic internal energy equations of the classical six-equation system. This alternative formulation allows us to easily design a simple numerical method that ensures consistency with mixture total energy conservation at the discrete level and agreement of the relaxed pressure at equilibrium with the correct mixture equation of state. Temperature and Gibbs free energy exchange terms are included in the equations as relaxation terms to model heat and mass transfer and hence liquidvapor transition. The algorithm uses a high-resolution wave propagation method for the numerical approximation of the homogeneous hyperbolic portion of the model. In two dimensions a fully-discretized scheme based on a hybrid HLLC/Roe Riemann solver is employed. Thermo-chemical terms are handled numerically via a stiff relaxation solver that forces thermodynamic equilibrium at liquidvapor interfaces under metastable conditions. We present numerical results of sample tests in one and two space dimensions that show the ability of the proposed model to describe cavitation mechanisms and evaporation wave dynamics.

  17. Use of ARM observations and numerical models to determine radiative and latent heating profiles of mesoscale convective systems for general circulation models

    SciTech Connect (OSTI)

    Houze, Jr., Robert A.

    2013-11-13

    We examined cloud radar data in monsoon climates, using cloud radars at Darwin in the Australian monsoon, on a ship in the Bay of Bengal in the South Asian monsoon, and at Niamey in the West African monsoon. We followed on with a more in-depth study of the continental MCSs over West Africa. We investigated whether the West African anvil clouds connected with squall line MCSs passing over the Niamey ARM site could be simulated in a numerical model by comparing the observed anvil clouds to anvil structures generated by the Weather Research and Forecasting (WRF) mesoscale model at high resolution using six different ice-phase microphysical schemes. We carried out further simulations with a cloud-resolving model forced by sounding network budgets over the Niamey region and over the northern Australian region. We have devoted some of the effort of this project to examining how well satellite data can determine the global breadth of the anvil cloud measurements obtained at the ARM ground sites. We next considered whether satellite data could be objectively analyzed to so that their large global measurement sets can be systematically related to the ARM measurements. Further differences were detailed between the land and ocean MCS anvil clouds by examining the interior structure of the anvils with the satellite-detected the CloudSat Cloud Profiling Radar (CPR). The satellite survey of anvil clouds in the Indo-Pacific region was continued to determine the role of MCSs in producing the cloud pattern associated with the MJO.

  18. Fish Individual-based Numerical Simulator (FINS): A particle-based model of juvenile salmonid movement and dissolved gas exposure history in the Columbia River Basin

    SciTech Connect (OSTI)

    Scheibe, Timothy D.; Richmond, Marshall C.

    2002-01-30

    This paper describes a numerical model of juvenile salmonid migration in the Columbia and Snake Rivers. The model, called the Fish Individual-based Numerical Simulator or FINS, employs a discrete, particle-based approach to simulate the migration and history of exposure to dissolved gases of individual fish. FINS is linked to a two-dimensional (vertically-averaged) hydrodynamic simulator that quantifies local water velocity, temperature, and dissolved gas levels as a function of river flow rates and dam operations. Simulated gas exposure histories can be input to biological mortality models to predict the effects of various river configurations on fish injury and mortality due to dissolved gas supersaturation. Therefore, FINS serves as a critical linkage between hydrodynamic models of the river system and models of biological impacts. FINS was parameterized and validated based on observations of individual fish movements collected using radiotelemetry methods during 1997 and 1998. A quasi-inverse approach was used to decouple fish swimming movements from advection with the local water velocity, allowing inference of time series of non-advective displacements of individual fish from the radiotelemetry data. Statistical analyses of these displacements are presented, and confirm that strong temporal correlation of fish swimming behavior persists in some cases over several hours. A correlated random-walk model was employed to simulate the observed migration behavior, and parameters of the model were estimated that lead to close correspondence between predictions and observations.

  19. Numerical implementation, verification and validation of two-phase flow four-equation drift flux model with Jacobian-free Newton–Krylov method

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Zou, Ling; Zhao, Haihua; Zhang, Hongbin

    2016-08-24

    This study presents a numerical investigation on using the Jacobian-free Newton–Krylov (JFNK) method to solve the two-phase flow four-equation drift flux model with realistic constitutive correlations (‘closure models’). The drift flux model is based on Isshi and his collaborators’ work. Additional constitutive correlations for vertical channel flow, such as two-phase flow pressure drop, flow regime map, wall boiling and interfacial heat transfer models, were taken from the RELAP5-3D Code Manual and included to complete the model. The staggered grid finite volume method and fully implicit backward Euler method was used for the spatial discretization and time integration schemes, respectively. Themore » Jacobian-free Newton–Krylov method shows no difficulty in solving the two-phase flow drift flux model with a discrete flow regime map. In addition to the Jacobian-free approach, the preconditioning matrix is obtained by using the default finite differencing method provided in the PETSc package, and consequently the labor-intensive implementation of complex analytical Jacobian matrix is avoided. Extensive and successful numerical verification and validation have been performed to prove the correct implementation of the models and methods. Code-to-code comparison with RELAP5-3D has further demonstrated the successful implementation of the drift flux model.« less

  20. Study on the well-posedness, convergence and the stability of the semi-implicit upwind numerical solver for the multi-fluid model

    SciTech Connect (OSTI)

    Lee, S. Y.; Park, C. E.; Hibiki, T.; Ishii, M.; Ransom, V. H.

    2012-07-01

    The well-posedness, convergence and the stability of the two-fluid code has been studied for a long time. Most of the investigations concern the semi-implicit upwind solution scheme for the six equation two-fluid model such as used in RELAP5 3 or TRACE 21. Since the system code, SPACE 2, adopts one more field, a droplet field, it consists of nine equations (3 mass, 3 momentum and 3 energy balance equations) and thus more involved investigations are necessary to confirm the stability and convergence. For this objective, the old issue of the well-posedness, convergence and the stability is revisited and some general guidelines to develop a well-posed numerical multi-fluid model are derived as follows; (1) Hyperbolicity of the corresponding system of partial differential equations is not a necessary condition for the development of a numerical model for multi-phase flow, but whether or not it is hyperbolic can provide guidance relative to initial conditions, boundary conditions, and expected high frequency behavior of the model. (2) A necessary condition for a well-posed numerical model is stability in the von Neumann sense, i.e. growth factor less than 1.0 for the shortest wave-length, 2{Delta}x. (3) The smallest node size used for convergence studies should be of the order of the characteristic dimension of the average description, i.e. smaller nodes can be used so long as they do not result in unphysical growth of wave-lengths less than the characteristic dimension. The usual mathematical definition of convergence i.e. the behavior of the solution as the node size approaches zero, is not appropriate for the discrete averaged numerical model, since there is diminished physical meaning to behavior at wavelengths less than the characteristic dimension of the average description. Under these guidelines, dispersion analysis and von Neumann stability analysis are performed for the three field multi-fluid, semi-implicit, upwind numerical model to show that the necessary

  1. Positive field-cooled dc susceptibility in granular superconductors interpreted through numerical simulations on a simple Josephson-junction-array model

    SciTech Connect (OSTI)

    Auletta, C.; Raiconi, G.; De Luca, R.; Pace, S.

    1995-05-01

    We have performed numerical simulations of a field-cooled dc susceptibility experiment carried out for granular superconductors by modeling these systems with a simple Josephson-junction array proposed by the authors. By this analysis the temperature dependence of the positive field-cooled susceptibility at very low values of the applied magnetic field, observed by Braunisch {ital et} {ital al}. [Phys. Rev. Lett. 68, 1908 (1992)] for some ceramic superonductors, has been reproduced and interpreted.

  2. Addressing numerical challenges in introducing a reactive transport code into a land surface model: a biogeochemical modeling proof-of-concept with CLM–PFLOTRAN 1.0

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Tang, Guoping; Yuan, Fengming; Bisht, Gautam; Hammond, Glenn E.; Lichtner, Peter C.; Kumar, Jitendra; Mills, Richard T.; Xu, Xiaofeng; Andre, Ben; Hoffman, Forrest M.; et al

    2016-03-04

    We explore coupling to a configurable subsurface reactive transport code as a flexible and extensible approach to biogeochemistry in land surface models. A reaction network with the Community Land Model carbon–nitrogen (CLM-CN) decomposition, nitrification, denitrification, and plant uptake is used as an example. We implement the reactions in the open-source PFLOTRAN (massively parallel subsurface flow and reactive transport) code and couple it with the CLM. To make the rate formulae designed for use in explicit time stepping in CLMs compatible with the implicit time stepping used in PFLOTRAN, the Monod substrate rate-limiting function with a residual concentration is used to represent the limitation ofmore » nitrogen availability on plant uptake and immobilization. We demonstrate that CLM–PFLOTRAN predictions (without invoking PFLOTRAN transport) are consistent with CLM4.5 for Arctic, temperate, and tropical sites.Switching from explicit to implicit method increases rigor but introduces numerical challenges. Care needs to be taken to use scaling, clipping, or log transformation to avoid negative concentrations during the Newton iterations. With a tight relative update tolerance (STOL) to avoid false convergence, an accurate solution can be achieved with about 50 % more computing time than CLM in point mode site simulations using either the scaling or clipping methods. The log transformation method takes 60–100 % more computing time than CLM. The computing time increases slightly for clipping and scaling; it increases substantially for log transformation for half saturation decrease from 10−3 to 10−9 mol m−3, which normally results in decreasing nitrogen concentrations. The frequent occurrence of very low concentrations (e.g. below nanomolar) can increase the computing time for clipping or scaling by about 20 %, double for log transformation. Overall, the log transformation method is accurate and robust, and the clipping and scaling

  3. Addressing numerical challenges in introducing a reactive transport code into a land surface model: A biogeochemical modeling proof-of-concept with CLM PFLOTRAN 1.0

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Tang, Guoping; Yuan, Fengming; Bisht, Gautam; Hammond, Glenn E.; Lichtner, Peter C.; Kumar, Jitendra; Mills, Richard T.; Xu, Xiaofeng; Andre, Ben; Hoffman, Forrest M.; et al

    2016-03-04

    Here, we explore coupling to a configurable subsurface reactive transport code as a flexible and extensible approach to biogeochemistry in land surface models. A reaction network with the Community Land Model carbon nitrogen (CLM-CN) decomposition, nitrification, denitrification, and plant uptake is used as an example. We implement the reactions in the open-source PFLOTRAN (massively parallel subsurface flow and reactive transport) code and couple it with the CLM. To make the rate formulae designed for use in explicit time stepping in CLMs compatible with the implicit time stepping used in PFLOTRAN, the Monod substrate rate-limiting function with a residual concentration ismore » used to represent the limitation of nitrogen availability on plant uptake and immobilization. We demonstrate that CLM PFLOTRAN predictions (without invoking PFLOTRAN transport) are consistent with CLM4.5 for Arctic, temperate, and tropical sites. Switching from explicit to implicit method increases rigor but introduces numerical challenges. Care needs to be taken to use scaling, clipping, or log transformation to avoid negative concentrations during the Newton iterations. With a tight relative update tolerance (STOL) to avoid false convergence, an accurate solution can be achieved with about 50 % more computing time than CLM in point mode site simulations using either the scaling or clipping methods. The log transformation method takes 60–100 % more computing time than CLM. The computing time increases slightly for clipping and scaling; it increases substantially for log transformation for half saturation decrease from 10–3 to 10–9 mol m–3, which normally results in decreasing nitrogen concentrations. The frequent occurrence of very low concentrations (e.g. below nanomolar) can increase the computing time for clipping or scaling by about 20 %, double for log transformation. Overall, the log transformation method is accurate and robust, and the clipping and scaling methods are

  4. Use of ARM observations and numerical models to determine radiative and latent heating profiles of mesoscale convective systems for general circulation models

    SciTech Connect (OSTI)

    Tao, Wei-Kuo; Houze, Robert, A., Jr.; Zeng, Xiping

    2013-03-14

    observations to evaluate model simulations In cooperation with Profs. Bob Houze at University of Washington and Steven Rutledge at Colorado State University, numerical model results were evaluated with observations from W- and C-band radars and CloudSat/TRMM satellites. These studies exhibited some shortcomings of current numerical models, such as too little of thin anvil clouds, directing the future improvement of cloud microphysics parameterization in CRMs. Two papers of Powell et al (2012) and Zeng et al. (2013), summarizing these studies, were published in the Journal of the Atmospheric Sciences. 4. Analyzed the water budgets of MCSs Using ARM data from TWP-ICE, ARM-SGP and other field campaigns, the Goddard CRM simulations were carried out to analyze the water budgets of clouds from TWP-ICE and AMMA. The simulations generated a set of datasets on clouds and radiation, which are available http://cloud.gsfc.nasa.gov/. The cloud datasets were available for modelers and other researchers aiming to improve the representation of cloud processes in multi-scale modeling frameworks, GCMs and climate models. Special datasets, such as 3D cloud distributions every six minutes for TWP-ICE, were requested and generated for ARM/ASR investigators. Data server records show that 86,206 datasets were downloaded by 120 users between April of 2010 and January of 2012. 5. MMF simulations The Goddard MMF (multi-scale modeling framework) has been improved by coupling with the Goddard Land Information System (LIS) and the Goddard Earth Observing System Model, Version 5 (GOES5). It has also been optimized on NASA HEC supercomputers and can be run over 4000 CPUs. The improved MMF with high horizontal resolution (1 x 1 degree) is currently being applied to cases covering 2005 and 2006. The results show that the spatial distribution pattern of precipitation rate is well simulated by the MMF through comparisons with satellite retrievals from the CMOPRH and GPCP data sets. In addition, the MMF results

  5. A three-dimensional numerical model of predevelopment conditions in the Death Valley regional ground-water flow system, Nevada and California

    SciTech Connect (OSTI)

    D'Agnese, F.A.; O'Brien, G.M.; Faunt, C.C.; Belcher, W.R.; San Juan, Carma

    2002-11-22

    In the early 1990's, two numerical models of the Death Valley regional ground-water flow system were developed by the U.S. Department of Energy. In general, the two models were based on the same basic hydrogeologic data set. In 1998, the U.S. Department of Energy requested that the U.S. Geological Survey develop and maintain a ground-water flow model of the Death Valley region in support of U.S. Department of Energy programs at the Nevada Test Site. The purpose of developing this ''second-generation'' regional model was to enhance the knowledge and understanding of the ground-water flow system as new information and tools are developed. The U.S. Geological Survey also was encouraged by the U.S. Department of Energy to cooperate to the fullest extent with other Federal, State, and local entities in the region to take advantage of the benefits of their knowledge and expertise. The short-term objective of the Death Valley regional ground-water flow system project was to develop a steady-stat e representation of the predevelopment conditions of the ground-water flow system utilizing the two geologic interpretations used to develop the previous numerical models. The long-term objective of this project was to construct and calibrate a transient model that simulates the ground-water conditions of the study area over the historical record that utilizes a newly interpreted hydrogeologic conceptual model. This report describes the result of the predevelopment steady-state model construction and calibration.

  6. Numerical modeling of regional ground-water flow in the deep-basin brine aquifer of the Palo Duro Basin, Texas Panhandle

    SciTech Connect (OSTI)

    Wirojanagud, P.; Kreitler, C.W.; Smith, D.A.

    1986-01-01

    Bedded Permian-age evaporite sequences in the Palo Duro Basin are being considered for a permanent nuclear waste repository by the U.S. Department of Energy. The purpose of this modeling study is to provide an understanding of regional ground-water flow in the formations beneath the Permian evaporite section. From this understanding, more detailed, smaller scale studies can be designed. This study is also intended to provide a better understanding of the boundary conditions and permeabilities of the aquifer and aquitard system as well as provide estimates of ground-water travel times across the basin. Numerical simulations were made of the Wolfcamp aquifer modeled as a single layer and of the entire Deep-Basin Brine aquifer system, including the Wolfcamp aquifer, modeled as a single layer.

  7. Numerical modeling of a 2K J-T heat exchanger used in Fermilab Vertical Test Stand VTS-1

    SciTech Connect (OSTI)

    Gupta, Prabhat Kumar; Rabehl, Roger

    2014-07-01

    Fermilab Vertical Test Stand-1 (VTS-1) is in operation since 2007 for testing the superconducting RF cavities at 2 K. This test stand has single layer coiled finned tubes heat exchanger before J-T valve. A finite difference based thermal model has been developed in Engineering Equation Solver (EES) to study its thermal performance during filling and refilling to maintain the constant liquid level of test stand. The model is also useful to predict its performance under other various operating conditions and will be useful to design the similar kind of heat exchanger for future needs. Present paper discusses the different operational modes of this heat exchanger and its thermal characteristics under these operational modes. Results of this model have also been compared with the experimental data gathered from the VTS-1 heat exchanger and they are in good agreement with the present model.

  8. DRSPALL: Impact of the Modification of the Numerical Spallings Model on Waste Isolation Pilot Plant Performance Assessment.

    SciTech Connect (OSTI)

    Kicker, Dwayne Curtis; Herrick, Courtney G.; Zeitler, Todd; Malama, Bwalya; Rudeen, David Keith; Gilkey, Amy P.

    2016-01-01

    The numerical code DRSPALL (from direct release spallings) is written to calculate the volume of Waste Isolation Pilot Plant (WIPP) solid waste subject to material failure and transport to the surface as a result of a hypothetical future inadvertent drilling intrusion. An error in the implementation of the DRSPALL finite difference equations was discovered as documented in Software Problem Report (SPR) 13-001. The modifications to DRSPALL to correct the finite difference equations are detailed, and verification and validation testing has been completed for the modified DRSPALL code. The complementary cumulative distribution function (CCDF) of spallings releases obtained using the modified DRSPALL is higher compared to that found in previous WIPP performance assessment (PA) calculations. Compared to previous PAs, there was an increase in the number of vectors that result in a nonzero spallings volume, which generally translates to an increase in spallings releases. The overall mean CCDFs for total releases using the modified DRSPALL are virtually unchanged, thus the modification to DRSPALL did not impact WIPP PA calculation results.

  9. Evaluation of cloud prediction and determination of critical relative humidity for a mesoscale numerical weather prediction model

    SciTech Connect (OSTI)

    Seaman, N.L.; Guo, Z.; Ackerman, T.P.

    1996-04-01

    Predictions of cloud occurrence and vertical location from the Pennsylvannia State University/National Center for Atmospheric Research nonhydrostatic mesoscale model (MM5) were evaluated statistically using cloud observations obtained at Coffeyville, Kansas, as part of the Second International satellite Cloud Climatology Project Regional Experiment campaign. Seventeen cases were selected for simulation during a November-December 1991 field study. MM5 was used to produce two sets of 36-km simulations, one with and one without four-dimensional data assimilation (FDDA), and a set of 12-km simulations without FDDA, but nested within the 36-km FDDA runs.

  10. Numerical Modeling of Complex Targets for High-Energy- Density Experiments with Ion Beams and other Drivers

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Koniges, Alice; Liu, Wangyi; Lidia, Steven; Schenkel, Thomas; Barnard, John; Friedman, Alex; Eder, David; Fisher, Aaron; Masters, Nathan

    2016-03-01

    We explore the simulation challenges and requirements for experiments planned on facilities such as the NDCX-II ion accelerator at LBNL, currently undergoing commissioning. Hydrodynamic modeling of NDCX-II experiments include certain lower temperature effects, e.g., surface tension and target fragmentation, that are not generally present in extreme high-energy laser facility experiments, where targets are completely vaporized in an extremely short period of time. Target designs proposed for NDCX-II range from metal foils of order one micron thick (thin targets) to metallic foam targets several tens of microns thick (thick targets). These high-energy-density experiments allow for the study of fracture as wellmore » as the process of bubble and droplet formation. We incorporate these physics effects into a code called ALE-AMR that uses a combination of Arbitrary Lagrangian Eulerian hydrodynamics and Adaptive Mesh Refinement. Inclusion of certain effects becomes tricky as we must deal with non-orthogonal meshes of various levels of refinement in three dimensions. A surface tension model used for droplet dynamics is implemented in ALE-AMR using curvature calculated from volume fractions. Thick foam target experiments provide information on how ion beam induced shock waves couple into kinetic energy of fluid flow. Although NDCX-II is not fully commissioned, experiments are being conducted that explore material defect production and dynamics.« less

  11. A summary of methods for approximating salt creep and disposal room closure in numerical models of multiphase flow

    SciTech Connect (OSTI)

    Freeze, G.A.; Larson, K.W.; Davies, P.B.

    1995-10-01

    Eight alternative methods for approximating salt creep and disposal room closure in a multiphase flow model of the Waste Isolation Pilot Plant (WIPP) were implemented and evaluated: Three fixed-room geometries three porosity functions and two fluid-phase-salt methods. The pressure-time-porosity line interpolation method is the method used in current WIPP Performance Assessment calculations. The room closure approximation methods were calibrated against a series of room closure simulations performed using a creep closure code, SANCHO. The fixed-room geometries did not incorporate a direct coupling between room void volume and room pressure. The two porosity function methods that utilized moles of gas as an independent parameter for closure coupling. The capillary backstress method was unable to accurately simulate conditions of re-closure of the room. Two methods were found to be accurate enough to approximate the effects of room closure; the boundary backstress method and pressure-time-porosity line interpolation. The boundary backstress method is a more reliable indicator of system behavior due to a theoretical basis for modeling salt deformation as a viscous process. It is a complex method and a detailed calibration process is required. The pressure lines method is thought to be less reliable because the results were skewed towards SANCHO results in simulations where the sequence of gas generation was significantly different from the SANCHO gas-generation rate histories used for closure calibration. This limitation in the pressure lines method is most pronounced at higher gas-generation rates and is relatively insignificant at lower gas-generation rates. Due to its relative simplicity, the pressure lines method is easier to implement in multiphase flow codes and simulations have a shorter execution time.

  12. Numerical modeling of heat-transfer and the influence of process parameters on tailoring the grain morphology of IN718 in electron beam additive manufacturing

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Raghavan, Narendran; Dehoff, Ryan; Pannala, Sreekanth; Simunovic, Srdjan; Kirka, Michael; Turner, John; Carlson, Neil; Babu, Sudarsanam S.

    2016-04-26

    The fabrication of 3-D parts from CAD models by additive manufacturing (AM) is a disruptive technology that is transforming the metal manufacturing industry. The correlation between solidification microstructure and mechanical properties has been well understood in the casting and welding processes over the years. This paper focuses on extending these principles to additive manufacturing to understand the transient phenomena of repeated melting and solidification during electron beam powder melting process to achieve site-specific microstructure control within a fabricated component. In this paper, we have developed a novel melt scan strategy for electron beam melting of nickel-base superalloy (Inconel 718) andmore » also analyzed 3-D heat transfer conditions using a parallel numerical solidification code (Truchas) developed at Los Alamos National Laboratory. The spatial and temporal variations of temperature gradient (G) and growth velocity (R) at the liquid-solid interface of the melt pool were calculated as a function of electron beam parameters. By manipulating the relative number of voxels that lie in the columnar or equiaxed region, the crystallographic texture of the components can be controlled to an extent. The analysis of the parameters provided optimum processing conditions that will result in columnar to equiaxed transition (CET) during the solidification. Furthermore, the results from the numerical simulations were validated by experimental processing and characterization thereby proving the potential of additive manufacturing process to achieve site-specific crystallographic texture control within a fabricated component.« less

  13. Self-Consistant Numerical Modeling of E-Cloud Driven Instability of a Bunch Train in the CERN SPS

    SciTech Connect (OSTI)

    Vay, J-L.; Furman, M.A.; Secondo, R.; Venturini, M.; Fox, J.D.; Rivetta, C.H,

    2010-09-01

    The simulation package WARP-POSINST was recently upgraded for handling multiple bunches and modeling concurrently the electron cloud buildup and its effect on the beam, allowing for direct self-consistent simulation of bunch trains generating, and interacting with, electron clouds. We have used the WARP-POSINST package on massively parallel supercomputers to study the growth rate and frequency patterns in space-time of the electron cloud driven transverse instability for a proton bunch train in the CERN SPS accelerator. Results suggest that a positive feedback mechanism exists between the electron buildup and the e-cloud driven transverse instability, leading to a net increase in predicted electron density. Comparisons to selected experimental data are also given. Electron clouds have been shown to trigger fast growing instabilities on proton beams circulating in the SPS and other accelerators. So far, simulations of electron cloud buildup and their effects on beam dynamics have been performed separately. This is a consequence of the large computational cost of the combined calculation due to large space and time scale disparities between the two processes. We have presented the latest improvements of the simulation package WARP-POSINST for the simulation of self-consistent ecloud effects, including mesh refinement, and generation of electrons from gas ionization and impact at the pipe walls. We also presented simulations of two consecutive bunches interacting with electrons clouds in the SPS, which included generation of secondary electrons. The distribution of electrons in front of the first beam was initialized from a dump taken from a preceding buildup calculation using the POSINST code. In this paper, we present an extension of this work where one full batch of 72 bunches is simulated in the SPS, including the entire buildup calculation and the self-consistent interaction between the bunches and the electrons. Comparisons to experimental data are also given.

  14. Direct Numerical Modeling of E-Cloud Driven Instability of a Bunch Train in the CERN SPS

    SciTech Connect (OSTI)

    Vay, J-L.; Furman, M.A.; Venturini, M.

    2011-03-01

    The simulation package WARP-POSINST was recently upgraded for handling multiple bunches and modeling concurrently the electron cloud buildup and its effect on the beam, allowing for direct self-consistent simulation of bunch trains generating, and interacting with, electron clouds. We have used the WARP-POSINST package on massively parallel supercomputers to study the buildup and interaction of electron clouds with a proton bunch train in the CERN SPS accelerator. Results suggest that a positive feedback mechanism exists between the electron buildup and the e-cloud driven transverse instability, leading to a net increase in predicted electron density. Electron clouds have been shown to trigger fast growing instabilities on proton beams circulating in the SPS and other accelerators. So far, simulations of electron cloud buildup and their effects on beam dynamics have been performed separately. This is a consequence of the large computational cost of the combined calculation due to large space and time scale disparities between the two processes. We have presented the latest improvements of the simulation package WARP-POSINST for the simulation of self-consistent ecloud effects, including mesh refinement, and generation of electrons from gas ionization and impact at the pipe walls. We also presented simulations of two consecutive bunches interacting with electrons clouds in the SPS, which included generation of secondary electrons. The distribution of electrons in front of the first beam was initialized from a dump taken from a preceding buildup calculation using the POSINST code. In this paper, we present an extension of this work where one full batch of 72 bunches is simulated in the SPS, including the entire buildup calculation and the self-consistent interaction between the bunches and the electrons.

  15. Numerical Modeling of 90Sr and 137Cs Transport from a Spill in the B-Cell of the 324 Building, Hanford Site 300 Area

    SciTech Connect (OSTI)

    Rockhold, Mark L.; Bacon, Diana H.; Freedman, Vicky L.; Lindberg, Michael J.; Clayton, Ray E.

    2012-03-19

    To characterize the extent of contamination under the 324 Building, a pit was excavated on the north side of the building in 2010 by Washington Closure Hanford LLC (WCH). Horizontal closed-end steel access pipes were installed under the foundation of the building from this pit and were used for measuring temperatures and exposure rates under the B-Cell. The deployed sensors measured elevated temperatures of up to 61 C (142 F) and exposure rates of up to 8,900 R/hr. WCH suspended deactivation of the facility because it recognized that building safety systems and additional characterization data might be needed for remediation of the contaminated material. The characterization work included additional field sampling, laboratory measurements, and numerical flow and transport modeling. Laboratory measurements of sediment physical, hydraulic, and geochemical properties were performed by Pacific Northwest National Laboratory (PNNL) and others. Geochemical modeling and subsurface flow and transport modeling also were performed by PNNL to evaluate the possible extent of contamination in the unsaturated sand and gravel sediments underlying the building. Historical records suggest that the concentrated 137Cs- and 90Sr-bearing liquid wastes that were spilled in B-Cell were likely from a glass-waste repository testing program associated with the Federal Republic of Germany (FRG). Incomplete estimates of the aqueous chemical composition (no anion data provided) of the FRG waste solutions were entered into a geochemical speciation model and were charge balanced with nitrate to estimate waste composition. Additional geochemical modeling was performed to evaluate reactions of the waste stream with the concrete foundation of the building prior to the stream entering the subsurface.

  16. Concentrating Solar Power: Advanced Projects Offering Low LCOE

    Office of Environmental Management (EM)

    Energy Compliance Order issued to Los Alamos National Laboratory Compliance Order issued to Los Alamos National Laboratory Pursuant to the authority of the Secretary of Energy under section 234B of the Atomic Energy Act of 1954, as amended, and 10 C.F.R. 5 824.4(b) of the Department's Procedural Rules for the Assessment of Civil Penalties for Classzjied Information Security Violations, I am today issuing the enclosed Compliance Order to Los Alamos National Security, LLC (LANS). Compliance

  17. Numerical Modeling | Open Energy Information

    Open Energy Info (EERE)

    4.0 4.1 Jerome Sacks,William Welch,Toby Mitchell,Henry Wynn. 1989. Design and Analysis of Computer Experiments. Statistical Science. . Page Area Activity Start Date Activity End...

  18. Numerical simulations for low energy nuclear reactions including...

    Office of Scientific and Technical Information (OSTI)

    Numerical simulations for low energy nuclear reactions including direct channels to validate statistical models Citation Details In-Document Search Title: Numerical simulations for ...

  19. Sandia Energy - Numerical Simulations of Hydrokinetics in the...

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    Numerical Simulations of Hydrokinetics in the Roza Canal, Yakima Washington Home Renewable Energy Energy Water Power Computational Modeling & Simulation Numerical Simulations of...

  20. NUMERICAL MODELING OF THE DISRUPTION OF COMET D/1993 F2 SHOEMAKER-LEVY 9 REPRESENTING THE PROGENITOR BY A GRAVITATIONALLY BOUND ASSEMBLAGE OF RANDOMLY SHAPED POLYHEDRA

    SciTech Connect (OSTI)

    Movshovitz, Naor; Asphaug, Erik; Korycansky, Donald

    2012-11-10

    We advance the modeling of rubble-pile solid bodies by re-examining the tidal breakup of comet Shoemaker-Levy 9, an event that occurred during a 1.33 R encounter with Jupiter in 1992 July. Tidal disruption of the comet nucleus led to a chain of sub-nuclei {approx}100-1000 m diameter; these went on to collide with the planet two years later. They were intensively studied prior to and during the collisions, making SL9 the best natural benchmark for physical models of small-body disruption. For the first time in the study of this event, we use numerical codes treating rubble piles as collections of polyhedra. This introduces forces of dilatation and friction, and inelastic response. As in our previous studies we conclude that the progenitor must have been a rubble pile, and we obtain approximately the same pre-breakup diameter ({approx}1.5 km) in our best fits to the data. We find that the inclusion of realistic fragment shapes leads to grain locking and dilatancy, so that even in the absence of friction or other dissipation we find that disruption is overall more difficult than in our spheres-based simulations. We constrain the comet's bulk density at {rho}{sub bulk} {approx} 300-400 kg m{sup -3}, half that of our spheres-based predictions and consistent with recent estimates derived from spacecraft observations.

  1. Confidence in Numerical Simulations

    SciTech Connect (OSTI)

    Hemez, Francois M.

    2015-02-23

    This PowerPoint presentation offers a high-level discussion of uncertainty, confidence and credibility in scientific Modeling and Simulation (M&S). It begins by briefly evoking M&S trends in computational physics and engineering. The first thrust of the discussion is to emphasize that the role of M&S in decision-making is either to support reasoning by similarity or to “forecast,” that is, make predictions about the future or extrapolate to settings or environments that cannot be tested experimentally. The second thrust is to explain that M&S-aided decision-making is an exercise in uncertainty management. The three broad classes of uncertainty in computational physics and engineering are variability and randomness, numerical uncertainty and model-form uncertainty. The last part of the discussion addresses how scientists “think.” This thought process parallels the scientific method where by a hypothesis is formulated, often accompanied by simplifying assumptions, then, physical experiments and numerical simulations are performed to confirm or reject the hypothesis. “Confidence” derives, not just from the levels of training and experience of analysts, but also from the rigor with which these assessments are performed, documented and peer-reviewed.

  2. Cost and Performance Assumptions for Modeling Electricity Generation Technologies

    SciTech Connect (OSTI)

    Tidball, Rick; Bluestein, Joel; Rodriguez, Nick; Knoke, Stu

    2010-11-01

    The goal of this project was to compare and contrast utility scale power plant characteristics used in data sets that support energy market models. Characteristics include both technology cost and technology performance projections to the year 2050. Cost parameters include installed capital costs and operation and maintenance (O&M) costs. Performance parameters include plant size, heat rate, capacity factor or availability factor, and plant lifetime. Conventional, renewable, and emerging electricity generating technologies were considered. Six data sets, each associated with a different model, were selected. Two of the data sets represent modeled results, not direct model inputs. These two data sets include cost and performance improvements that result from increased deployment as well as resulting capacity factors estimated from particular model runs; other data sets represent model input data. For the technologies contained in each data set, the levelized cost of energy (LCOE) was also evaluated, according to published cost, performance, and fuel assumptions.

  3. Reference Inflow Characterization for River Resource Reference Model (RM2)

    SciTech Connect (OSTI)

    Neary, Vincent S

    2011-12-01

    Sandia National Laboratory (SNL) is leading an effort to develop reference models for marine and hydrokinetic technologies and wave and current energy resources. This effort will allow the refinement of technology design tools, accurate estimates of a baseline levelized cost of energy (LCoE), and the identification of the main cost drivers that need to be addressed to achieve a competitive LCoE. As part of this effort, Oak Ridge National Laboratory was charged with examining and reporting reference river inflow characteristics for reference model 2 (RM2). Published turbulent flow data from large rivers, a water supply canal and laboratory flumes, are reviewed to determine the range of velocities, turbulence intensities and turbulent stresses acting on hydrokinetic technologies, and also to evaluate the validity of classical models that describe the depth variation of the time-mean velocity and turbulent normal Reynolds stresses. The classical models are found to generally perform well in describing river inflow characteristics. A potential challenge in river inflow characterization, however, is the high variability of depth and flow over the design life of a hydrokinetic device. This variation can have significant effects on the inflow mean velocity and turbulence intensity experienced by stationary and bottom mounted hydrokinetic energy conversion devices, which requires further investigation, but are expected to have minimal effects on surface mounted devices like the vertical axis turbine device designed for RM2. A simple methodology for obtaining an approximate inflow characterization for surface deployed devices is developed using the relation umax=(7/6)V where V is the bulk velocity and umax is assumed to be the near-surface velocity. The application of this expression is recommended for deriving the local inflow velocity acting on the energy extraction planes of the RM2 vertical axis rotors, where V=Q/A can be calculated given a USGS gage flow time

  4. Numerical computation of Pop plot

    SciTech Connect (OSTI)

    Menikoff, Ralph

    2015-03-23

    The Pop plot — distance-of-run to detonation versus initial shock pressure — is a key characterization of shock initiation in a heterogeneous explosive. Reactive burn models for high explosives (HE) must reproduce the experimental Pop plot to have any chance of accurately predicting shock initiation phenomena. This report describes a methodology for automating the computation of a Pop plot for a specific explosive with a given HE model. Illustrative examples of the computation are shown for PBX 9502 with three burn models (SURF, WSD and Forest Fire) utilizing the xRage code, which is the Eulerian ASC hydrocode at LANL. Comparison of the numerical and experimental Pop plot can be the basis for a validation test or as an aid in calibrating the burn rate of an HE model. Issues with calibration are discussed.

  5. Modeling pure methane hydrate dissociation using a numerical simulator from a novel combination of X-ray computed tomography and macroscopic data

    SciTech Connect (OSTI)

    Gupta, A.; Moridis, G.J.; Kneafsey, T.J.; Sloan, Jr., E.D.

    2009-08-15

    The numerical simulator TOUGH+HYDRATE (T+H) was used to predict the transient pure methane hydrate (no sediment) dissociation data. X-ray computed tomography (CT) was used to visualize the methane hydrate formation and dissociation processes. A methane hydrate sample was formed from granular ice in a cylindrical vessel, and slow depressurization combined with thermal stimulation was applied to dissociate the hydrate sample. CT images showed that the water produced from the hydrate dissociation accumulated at the bottom of the vessel and increased the hydrate dissociation rate there. CT images were obtained during hydrate dissociation to confirm the radial dissociation of the hydrate sample. This radial dissociation process has implications for dissociation of hydrates in pipelines, suggesting lower dissociation times than for longitudinal dissociation. These observations were also confirmed by the numerical simulator predictions, which were in good agreement with the measured thermal data during hydrate dissociation. System pressure and sample temperature measured at the sample center followed the CH{sub 4} hydrate L{sub w}+H+V equilibrium line during hydrate dissociation. The predicted cumulative methane gas production was within 5% of the measured data. Thus, this study validated our simulation approach and assumptions, which include stationary pure methane hydrate-skeleton, equilibrium hydrate-dissociation and heat- and mass-transfer in predicting hydrate dissociation in the absence of sediments. It should be noted that the application of T+H for the pure methane hydrate system (no sediment) is outside the general applicability limits of T+H.

  6. 2011 Cost of Wind Energy Review

    SciTech Connect (OSTI)

    Tegen, S.; Lantz, E.; Hand, M.; Maples, B.; Smith, A.; Schwabe, P.

    2013-03-01

    This report describes the levelized cost of energy (LCOE) for a typical land-based wind turbine installed in the United States in 2011, as well as the modeled LCOE for a fixed-bottom offshore wind turbine installed in the United States in 2011. Each of the four major components of the LCOE equation are explained in detail, such as installed capital cost, annual energy production, annual operating expenses, and financing, and including sensitivity ranges that show how each component can affect LCOE. These LCOE calculations are used for planning and other purposes by the U.S. Department of Energy's Wind Program.

  7. Numerical Simulation of Inter-basin Groundwater Flow into Northern Yucca Flat, Nevada National Security Site, Using the Death Valley Regional Flow System Model

    SciTech Connect (OSTI)

    Pohlmann Karl,Ye Ming

    2012-03-01

    Models of groundwater flow for the Yucca Flat area of the Nevada National Security Site (NNSS) are under development by the U.S. Department of Energy (DOE) for corrective action investigations of the Yucca Flat-Climax Mine Corrective Action Unit (CAU). One important aspect of these models is the quantity of inter-basin groundwater flow from regional systems to the north. This component of flow, together with its uncertainty, must be properly accounted for in the CAU flow models to provide a defensible regional framework for calculations of radionuclide transport that will support determinations of the Yucca Flat-Climax Mine contaminant boundary. Because characterizing flow boundary conditions in northern Yucca Flat requires evaluation to a higher level of detail than the scale of the Yucca Flat-Climax Mine CAU model can efficiently provide, a study more focused on this aspect of the model was required.

  8. Disruptive Innovation in Numerical Hydrodynamics

    SciTech Connect (OSTI)

    Waltz, Jacob I.

    2012-09-06

    We propose the research and development of a high-fidelity hydrodynamic algorithm for tetrahedral meshes that will lead to a disruptive innovation in the numerical modeling of Laboratory problems. Our proposed innovation has the potential to reduce turnaround time by orders of magnitude relative to Advanced Simulation and Computing (ASC) codes; reduce simulation setup costs by millions of dollars per year; and effectively leverage Graphics Processing Unit (GPU) and future Exascale computing hardware. If successful, this work will lead to a dramatic leap forward in the Laboratory's quest for a predictive simulation capability.

  9. Numerical Prediction of Experimentally Observed Behavior of a Scale Model of an Offshore Wind Turbine Supported by a Tension-Leg Platform: Preprint

    SciTech Connect (OSTI)

    Prowell, I.; Robertson, A.; Jonkman, J.; Stewart, G. M.; Goupee, A. J.

    2013-01-01

    Realizing the critical importance the role physical experimental tests play in understanding the dynamics of floating offshore wind turbines, the DeepCwind consortium conducted a one-fiftieth-scale model test program where several floating wind platforms were subjected to a variety of wind and wave loading condition at the Maritime Research Institute Netherlands wave basin. This paper describes the observed behavior of a tension-leg platform, one of three platforms tested, and the systematic effort to predict the measured response with the FAST simulation tool using a model primarily based on consensus geometric and mass properties of the test specimen.

  10. Addressing numerical challenges in introducing a reactive transport code into a land surface model: A biogeochemical modeling proof-of-concept with CLM-PFLOTRAN 1.0: Modeling Archive

    DOE Data Explorer [Office of Scientific and Technical Information (OSTI)]

    R. T. Mills; G. Bisht; G. E. Hammond; J. Kumar; P. C. Lichtner; F. M. Hoffman; X. Xu; F. Yuan; S. L. Painter; G. Tang; and P.E. Thornton; B. Andre

    2016-04-19

    This Modeling Archive is in support of an NGEE Arctic discussion paper under review and available at doi:10.5194/gmd-9-927-2016. The purpose is to document the simulations to allow verification, reproducibility, and follow-up studies. This dataset contains shell scripts to create the CLM-PFLOTRAN cases, specific input files for PFLOTRAN and CLM, outputs, and python scripts to make the figures using the outputs in the publication. Through these results, we demonstrate that CLM-PFLOTRAN can approximately reproduce CLM results in selected cases for the Arctic, temperate and tropic sites. In addition, the new framework facilitates mechanistic representations of soil biogeochemistry processes in the land surface model.

  11. Direct Numerical Simulation

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    Direct Numerical Simulation - Sandia Energy Energy Search Icon Sandia Home Locations Contact Us Employee Locator Energy & Climate Secure & Sustainable Energy Future Stationary Power Energy Conversion Efficiency Solar Energy Wind Energy Water Power Supercritical CO2 Geothermal Natural Gas Safety, Security & Resilience of the Energy Infrastructure Energy Storage Nuclear Power & Engineering Grid Modernization Battery Testing Nuclear Energy Defense Waste Management Programs Advanced

  12. The effect of variable atmospheric forcing on oceanic subduction of a passive tracer in a numerical model: Implications for global warming

    SciTech Connect (OSTI)

    Horsfall, F.; Bleck, R.; Hanson, H.P.

    1997-11-01

    This study addresses the issue of the ocean`s response to the changing climate. The objectives is to determine the effect of variable atmospheric forcing on the ocean on decadal time scales, specifically on the subduction of a passive tracer. In the context of the model used in this study, this tracer is {open_quotes}tagged{close_quotes} water that is subducted into the thermocline and into the deep ocean. The model used in this study is the Miami Isopycnic Coordinate Ocean Model which has a realistic Atlantic domain from 20{degrees}S to 60{degrees}N. There are twelve model layers, the first (top) layer being the thermodynamically active mixed layer and the lower eleven layers all having constant potential density ({sigma}{sub {theta}}). The atmospheric forcing changes vary latitudinally, allowing for a maximum increase in wind at midlatitudes and a maximum increase in temperature at the poles. In these experiments, it was found that wind speed and temperature effects dominate in bringing about changes in mixed-layer depth and in tracer penetration at high latitudes, with wind speed effects having the greater weight. It is apparent from the results that the weakening of the North Atlantic thermohaline circulation is dependent on the atmospheric changes in air temperature and in the wind field. 11 refs., 2 figs.

  13. Renewable Energy Cost Modeling. A Toolkit for Establishing Cost-Based Incentives in the United States

    SciTech Connect (OSTI)

    Gifford, Jason S.; Grace, Robert C.; Rickerson, Wilson H.

    2011-05-01

    This report serves as a resource for policymakers who wish to learn more about levelized cost of energy (LCOE) calculations, including cost-based incentives. The report identifies key renewable energy cost modeling options, highlights the policy implications of choosing one approach over the other, and presents recommendations on the optimal characteristics of a model to calculate rates for cost-based incentives, FITs, or similar policies. These recommendations shaped the design of NREL's Cost of Renewable Energy Spreadsheet Tool (CREST), which is used by state policymakers, regulators, utilities, developers, and other stakeholders to assist with analyses of policy and renewable energy incentive payment structures. Authored by Jason S. Gifford and Robert C. Grace of Sustainable Energy Advantage LLC and Wilson H. Rickerson of Meister Consultants Group, Inc.

  14. Mechanical diode: Comparing numerical and experimental characterizations

    SciTech Connect (OSTI)

    Sagartz, M.J.; Segalman, D.; Simmermacher, T.

    1998-02-01

    In this introductory work, joint compliance is studied in both a numerical and experimental setting. A simple bolted interface is used as the test article and compliance is measured for the joint in both compression and in tension. This simple interface is shown to exhibit a strong non-linearity near the transition from compression to tension (or vice-versa). Modeling issues pertaining to numerically solving for the compliance are addressed. It is shown that the model predictions, in spite of convergence being very sensitive to numerical artifacts of the interface model, are in good agreement with experimentally measured strains and joint compliances. The joint behavior is a mechanical analogy to a diode, i.e., in compression, the joint is very stiff, acting almost as a rigid link, while in tension the joint is relatively soft, acting as a spring.

  15. Direct Numerical Simulation

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    Infrastructure Energy Storage Nuclear Power & Engineering Grid Modernization Battery Testing ... Heavy Duty Fuels DISI Combustion HCCISCCI Fundamentals Spray Combustion Modeling ...

  16. Rythmos Numerical Integration Package

    Energy Science and Technology Software Center (OSTI)

    2006-09-01

    Rythmos numerically integrates transient differential equations. The differential equations can be explicit or implicit ordinary differential equations ofr formulated as fully implicit differential-algebraic equations. Methods include backward Euler, forward Euler, explicit Runge-Kutta, and implicit BDF at this time. Native support for operator split methods and strict modularity are strong design goals. Forward sensitivity computations will be included in the first release with adjoint sensitivities coming in the near future. Rythmos heavily relies on Thyra formore » linear algebra and nonlinear solver interfaces to AztecOO, Amesos, IFPack, and NOX in Tilinos. Rythmos is specially suited for stiff differential equations and thos applictions where operator split methods have a big advantage, e.g. Computational fluid dynamics, convection-diffusion equations, etc.« less

  17. Numerical simulations of strong incompressible magnetohydrodynamic turbulence

    SciTech Connect (OSTI)

    Mason, J.; Cattaneo, F.; Perez, J. C.; Boldyrev, S.

    2012-05-15

    Magnetised plasma turbulence pervades the universe and is likely to play an important role in a variety of astrophysical settings. Magnetohydrodynamics (MHD) provides the simplest theoretical framework in which phenomenological models for the turbulent dynamics can be built. Numerical simulations of MHD turbulence are widely used to guide and test the theoretical predictions; however, simulating MHD turbulence and accurately measuring its scaling properties is far from straightforward. Computational power limits the calculations to moderate Reynolds numbers and often simplifying assumptions are made in order that a wider range of scales can be accessed. After describing the theoretical predictions and the numerical approaches that are often employed in studying strong incompressible MHD turbulence, we present the findings of a series of high-resolution direct numerical simulations. We discuss the effects that insufficiencies in the computational approach can have on the solution and its physical interpretation.

  18. An integrated experimental and numerical study: Developing a reaction

    Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site

    transport model that couples chemical reactions of mineral dissolution/precipitation with spatial and temporal flow variations in CO2/brine/rock systems | Department of Energy An integrated experimental and numerical study: Developing a reaction transport model that couples chemical reactions of mineral dissolution/precipitation with spatial and temporal flow variations in CO2/brine/rock systems An integrated experimental and numerical study: Developing a reaction transport model that

  19. Simple intrinsic defects in GaAs : numerical supplement.

    SciTech Connect (OSTI)

    Schultz, Peter Andrew

    2012-04-01

    This Report presents numerical tables summarizing properties of intrinsic defects in gallium arsenide, GaAs, as computed by density functional theory. This Report serves as a numerical supplement to the results published in: P.A. Schultz and O.A. von Lilienfeld, 'Simple intrinsic defects in GaAs', Modelling Simul. Mater. Sci Eng., Vol. 17, 084007 (2009), and intended for use as reference tables for a defect physics package in device models. The numerical results for density functional theory calculations of properties of simple intrinsic defects in gallium arsenide are presented.

  20. Numerical

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    ... saturation mechanism depends sensitively on the viscosity parameterized by Pm. ... A more thorough treatment of the dynamic variation of dominant eigenmodes can be found for ...

  1. RELAP-7 Numerical Stabilization: Entropy Viscosity Method

    SciTech Connect (OSTI)

    R. A. Berry; M. O. Delchini; J. Ragusa

    2014-06-01

    The RELAP-7 code is the next generation nuclear reactor system safety analysis code being developed at the Idaho National Laboratory (INL). The code is based on the INL's modern scientific software development framework, MOOSE (Multi-Physics Object Oriented Simulation Environment). The overall design goal of RELAP-7 is to take advantage of the previous thirty years of advancements in computer architecture, software design, numerical integration methods, and physical models. The end result will be a reactor systems analysis capability that retains and improves upon RELAP5's capability and extends the analysis capability for all reactor system simulation scenarios. RELAP-7 utilizes a single phase and a novel seven-equation two-phase flow models as described in the RELAP-7 Theory Manual (INL/EXT-14-31366). The basic equation systems are hyperbolic, which generally require some type of stabilization (or artificial viscosity) to capture nonlinear discontinuities and to suppress advection-caused oscillations. This report documents one of the available options for this stabilization in RELAP-7 -- a new and novel approach known as the entropy viscosity method. Because the code is an ongoing development effort in which the physical sub models, numerics, and coding are evolving, so too must the specific details of the entropy viscosity stabilization method. Here the fundamentals of the method in their current state are presented.

  2. LBNL/NREL Analysis Predicts Record Low LCOE for Wind Energy in...

    Broader source: Energy.gov (indexed) [DOE]

    turbine productivity will drive down the levelized cost of wind energy for U.S. wind projects constructed in 2012 - 2013. Estimates focused on changes in capital costs and turbine ...

  3. Stephanie von Numers | Department of Energy

    Broader source: Energy.gov (indexed) [DOE]

    Stephanie von Numers - Communications and Web Coordinator, Education & Workforce Development Stephanie von Numers joined the U.S. Department of Energy's Office of Energy Efficiency ...

  4. Recent advances in two-phase flow numerics

    SciTech Connect (OSTI)

    Mahaffy, J.H.; Macian, R.

    1997-07-01

    The authors review three topics in the broad field of numerical methods that may be of interest to individuals modeling two-phase flow in nuclear power plants. The first topic is iterative solution of linear equations created during the solution of finite volume equations. The second is numerical tracking of macroscopic liquid interfaces. The final area surveyed is the use of higher spatial difference techniques.

  5. Numerical simulations shed new light on early universe

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    Numerical simulations shed new light on early universe Numerical simulations shed new light on early universe The code simulates conditions during the first few minutes of cosmological evolution to model the role of neutrinos, nuclei and other particles in shaping the early universe. April 21, 2016 Los Alamos scientists developed the BURST computer code to predict-to unprecedented precision-the amounts of light nuclei synthesized in the Big Bang. Los Alamos scientists developed the BURST

  6. Numerical recipes for mold filling simulation

    SciTech Connect (OSTI)

    Kothe, D.; Juric, D.; Lam, K.; Lally, B.

    1998-07-01

    Has the ability to simulate the filling of a mold progressed to a point where an appropriate numerical recipe achieves the desired results? If results are defined to be topological robustness, computational efficiency, quantitative accuracy, and predictability, all within a computational domain that faithfully represents complex three-dimensional foundry molds, then the answer unfortunately remains no. Significant interfacial flow algorithm developments have occurred over the last decade, however, that could bring this answer closer to maybe. These developments have been both evolutionary and revolutionary, will continue to transpire for the near future. Might they become useful numerical recipes for mold filling simulations? Quite possibly. Recent progress in algorithms for interface kinematics and dynamics, linear solution methods, computer science issues such as parallelization and object-oriented programming, high resolution Navier-Stokes (NS) solution methods, and unstructured mesh techniques, must all be pursued as possible paths toward higher fidelity mold filling simulations. A detailed exposition of these algorithmic developments is beyond the scope of this paper, hence the authors choose to focus here exclusively on algorithms for interface kinematics. These interface tracking algorithms are designed to model the movement of interfaces relative to a reference frame such as a fixed mesh. Current interface tracking algorithm choices are numerous, so is any one best suited for mold filling simulation? Although a clear winner is not (yet) apparent, pros and cons are given in the following brief, critical review. Highlighted are those outstanding interface tracking algorithm issues the authors feel can hamper the reliable modeling of today`s foundry mold filling processes.

  7. Numerical uncertainty in computational engineering and physics

    SciTech Connect (OSTI)

    Hemez, Francois M

    2009-01-01

    Obtaining a solution that approximates ordinary or partial differential equations on a computational mesh or grid does not necessarily mean that the solution is accurate or even 'correct'. Unfortunately assessing the quality of discrete solutions by questioning the role played by spatial and temporal discretizations generally comes as a distant third to test-analysis comparison and model calibration. This publication is contributed to raise awareness of the fact that discrete solutions introduce numerical uncertainty. This uncertainty may, in some cases, overwhelm in complexity and magnitude other sources of uncertainty that include experimental variability, parametric uncertainty and modeling assumptions. The concepts of consistency, convergence and truncation error are overviewed to explain the articulation between the exact solution of continuous equations, the solution of modified equations and discrete solutions computed by a code. The current state-of-the-practice of code and solution verification activities is discussed. An example in the discipline of hydro-dynamics illustrates the significant effect that meshing can have on the quality of code predictions. A simple method is proposed to derive bounds of solution uncertainty in cases where the exact solution of the continuous equations, or its modified equations, is unknown. It is argued that numerical uncertainty originating from mesh discretization should always be quantified and accounted for in the overall uncertainty 'budget' that supports decision-making for applications in computational physics and engineering.

  8. Oceanic Trace Gases Numeric Data Packages from the Carbon Dioxide Information Analysis Center (CDIAC)

    DOE Data Explorer [Office of Scientific and Technical Information (OSTI)]

    CDIAC products include numeric data packages, publications, trend data, atlases, models, etc. and can be searched for by subject area, keywords, authors, product numbers, time periods, collection sites, spatial references, etc. Most data sets or packages, many with numerous data files, are free to download from CDIAC's ftp area. The collection provides access to the Oceanographic Numeric Data Packages (NDPs).

  9. Numerical Simulation of Groundwater Withdrawal at the Nevada Test Site

    SciTech Connect (OSTI)

    Carroll, Rosemary; Giroux, Brian; Pohll, Greg; Hershey, Ronald; Russell, Charles; Howcroft, William

    2004-01-28

    Alternative uses of the Nevada Test Site (NTS) may require large amounts of water to construct and/or operate. The only abundant source of water at the NTS is groundwater. This report describes preliminary modeling to quantify the amount of groundwater available for development from three hydrographic areas at the NTS. Modeling was conducted with a three-dimensional transient numerical groundwater flow model.

  10. Numerical simulations of capillary barrier field tests

    SciTech Connect (OSTI)

    Morris, C.E.; Stormont, J.C.

    1997-12-31

    Numerical simulations of two capillary barrier systems tested in the field were conducted to determine if an unsaturated flow model could accurately represent the observed results. The field data was collected from two 7-m long, 1.2-m thick capillary barriers built on a 10% grade that were being tested to investigate their ability to laterally divert water downslope. One system had a homogeneous fine layer, while the fine soil of the second barrier was layered to increase its ability to laterally divert infiltrating moisture. The barriers were subjected first to constant infiltration while minimizing evaporative losses and then were exposed to ambient conditions. The continuous infiltration period of the field tests for the two barrier systems was modelled to determine the ability of an existing code to accurately represent capillary barrier behavior embodied in these two designs. Differences between the field test and the model data were found, but in general the simulations appeared to adequately reproduce the response of the test systems. Accounting for moisture retention hysteresis in the layered system will potentially lead to more accurate modelling results and is likely to be important when developing reasonable predictions of capillary barrier behavior.

  11. Numerical Studies of Impurities in Fusion Plasmas

    DOE R&D Accomplishments [OSTI]

    Hulse, R. A.

    1982-09-01

    The coupled partial differential equations used to describe the behavior of impurity ions in magnetically confined controlled fusion plasmas require numerical solution for cases of practical interest. Computer codes developed for impurity modeling at the Princeton Plasma Physics Laboratory are used as examples of the types of codes employed for this purpose. These codes solve for the impurity ionization state densities and associated radiation rates using atomic physics appropriate for these low-density, high-temperature plasmas. The simpler codes solve local equations in zero spatial dimensions while more complex cases require codes which explicitly include transport of the impurity ions simultaneously with the atomic processes of ionization and recombination. Typical applications are discussed and computational results are presented for selected cases of interest.

  12. Interagency mechanical operations group numerical systems group

    SciTech Connect (OSTI)

    1997-09-01

    This report consists of the minutes of the May 20-21, 1971 meeting of the Interagency Mechanical Operations Group (IMOG) Numerical Systems Group. This group looks at issues related to numerical control in the machining industry. Items discussed related to the use of CAD and CAM, EIA standards, data links, and numerical control.

  13. Simple intrinsic defects in InAs : numerical predictions.

    SciTech Connect (OSTI)

    Schultz, Peter Andrew

    2013-03-01

    This Report presents numerical tables summarizing properties of intrinsic defects in indium arsenide, InAs, as computed by density functional theory using semi-local density functionals, intended for use as reference tables for a defect physics package in device models.

  14. Direct numerical simulation of turbulent reacting flows

    SciTech Connect (OSTI)

    Chen, J.H.

    1993-12-01

    The development of turbulent combustion models that reflect some of the most important characteristics of turbulent reacting flows requires knowledge about the behavior of key quantities in well defined combustion regimes. In turbulent flames, the coupling between the turbulence and the chemistry is so strong in certain regimes that is is very difficult to isolate the role played by one individual phenomenon. Direct numerical simulation (DNS) is an extremely useful tool to study in detail the turbulence-chemistry interactions in certain well defined regimes. Globally, non-premixed flames are controlled by two limiting cases: the fast chemistry limit, where the turbulent fluctuations. In between these two limits, finite-rate chemical effects are important and the turbulence interacts strongly with the chemical processes. This regime is important because industrial burners operate in regimes in which, locally the flame undergoes extinction, or is at least in some nonequilibrium condition. Furthermore, these nonequilibrium conditions strongly influence the production of pollutants. To quantify the finite-rate chemistry effect, direct numerical simulations are performed to study the interaction between an initially laminar non-premixed flame and a three-dimensional field of homogeneous isotropic decaying turbulence. Emphasis is placed on the dynamics of extinction and on transient effects on the fine scale mixing process. Differential molecular diffusion among species is also examined with this approach, both for nonreacting and reacting situations. To address the problem of large-scale mixing and to examine the effects of mean shear, efforts are underway to perform large eddy simulations of round three-dimensional jets.

  15. Numerical evaluation of effective unsaturated hydraulic properties...

    Office of Scientific and Technical Information (OSTI)

    unsaturated hydraulic properties for fractured rocks Citation Details In-Document Search Title: Numerical evaluation of effective unsaturated hydraulic properties for ...

  16. Category:Modeling Techniques | Open Energy Information

    Open Energy Info (EERE)

    Analytical Modeling 1 pages C Conceptual Model 1 pages M Modeling-Computer Simulations 1 pages N Numerical Modeling 1 pages P Portfolio Risk...

  17. Numerical simulations of the decay of primordial magnetic turbulence

    SciTech Connect (OSTI)

    Kahniashvili, Tina [McWilliams Center for Cosmology and Department of Physics, Carnegie Mellon University, 5000 Forbes Ave, Pittsburgh, Pennsylvania 15213 (United States); Department of Physics, Laurentian University, Ramsey Lake Road, Sudbury, ON P3E 2C (Canada); Abastumani Astrophysical Observatory, Ilia State University, 2A Kazbegi Ave, Tbilisi, GE-0160 (Georgia); Brandenburg, Axel [Nordita, AlbaNova University Center, Roslagstullsbacken 23, 10691 Stockholm (Sweden); Department of Astronomy, Stockholm University, SE 10691 Stockholm (Sweden); Tevzadze, Alexander G. [Abastumani Astrophysical Observatory, Ilia State University, 2A Kazbegi Ave, Tbilisi, GE-0160 (Georgia); Faculty of Exact and Natural Sciences, Tbilisi State University, 1 Chavchavadze Avenue Tbilisi, GE-0128 (Georgia); Ratra, Bharat [Department of Physics, Kansas State University, 116 Cardwell Hall, Manhattan, Kansas 66506 (United States)

    2010-06-15

    We perform direct numerical simulations of forced and freely decaying 3D magnetohydrodynamic turbulence in order to model magnetic field evolution during cosmological phase transitions in the early Universe. Our approach assumes the existence of a magnetic field generated either by a process during inflation or shortly thereafter, or by bubble collisions during a phase transition. We show that the final configuration of the magnetic field depends on the initial conditions, while the velocity field is nearly independent of initial conditions.

  18. Improvements to the RELAP5-3D Nearly-Implicit Numerical Scheme

    SciTech Connect (OSTI)

    Richard A. Riemke; Walter L. Weaver; RIchard R. Schultz

    2005-05-01

    The RELAP5-3D computer program has been improved with regard to its nearly-implicit numerical scheme for twophase flow and single-phase flow. Changes were made to the nearly-implicit numerical scheme finite difference momentum equations as follows: (1) added the velocity flip-flop mass/energy error mitigation logic, (2) added the modified Henry-Fauske choking model, (3) used the new time void fraction in the horizontal stratification force terms and gravity head, and (4) used an implicit form of the artificial viscosity. The code modifications allow the nearly-implicit numerical scheme to be more implicit and lead to enhanced numerical stability.

  19. Direct Numerical Simulation of Compressible, Turbulent Flow ...

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    The computational mesh for this direct numerical simulation was over 33 billion cells, and was run on up to 102,400 cores under a DoD HPCMP Frontier Project. Nicholas Bisek and ...

  20. Numerical routines for predicting ignition in pyrotechnic devices

    SciTech Connect (OSTI)

    Pierce, K.G.

    1986-06-01

    Two numerical models of the thermal processes leading to ignition in a pyrotechnic device have been developed. These models are based on finite difference approximations to the heat diffusion equation, with temperature-dependent thermal properties, in a single spatial coordinate. The derivation of the finite difference equations is discussed and the methods employed at boundaries and interfaces are given. The sources of the thermal-properties data are identified and how these data are used is explained. The program structure is explained and example runs of the programs are given.

  1. Numerical simulation of the flow in wire-wrapped pin bundles: Effect of

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    pin-wire contact modeling | Argonne Leadership Computing Facility simulation of the flow in wire-wrapped pin bundles: Effect of pin-wire contact modeling Authors: Merzari, E., Smith, J.G., Tentner, A., Pointer, W.D., Fischer, P. The rapid advancement of numerical techniques and the availability of increasingly powerful supercomputers recently enabled scientists to use large eddy simulation (LES) to simulate numerically the flow in a full subassembly composed of wire-wrapped pins. Because of

  2. Numerical Modeling At Coso Geothermal Area (2000) | Open Energy...

    Open Energy Info (EERE)

    porosity, which are then used to constrain and delineate possible zones of intense heat, fracture accumulation and fluid saturation. Poisson's ratio at Coso ranges from 0.15 to...

  3. Numerical Modeling At Dixie Valley Geothermal Area (Iovenitti...

    Open Energy Info (EERE)

    H. Ibser, Jennifer Lewicki, B. Mack. Kennedy, Michael Swyer (2013) Egs Exploration Methodology Project Using the Dixie Valley Geothermal System, Nevada, Status Update Christoph...

  4. Final Report on Experimental and Numerical Modeling Activities...

    Office of Scientific and Technical Information (OSTI)

    Authors: Mukhopadhyay, Sumit ; Spycher, Nicolas ; Pester, Nick ; Saldi, Giuseppe ; Beyer, John ; Houseworth, Jim ; Knauss, Kevin Publication Date: 2014-09-04 OSTI Identifier: ...

  5. A NUMERICAL MODEL OF STANDARD TO BLOWOUT JETS

    SciTech Connect (OSTI)

    Archontis, V.; Hood, A. W.

    2013-06-01

    We report on three-dimensional (3D) MHD simulations of the formation of jets produced during the emergence and eruption of solar magnetic fields. The interaction between an emerging and an ambient magnetic field in the solar atmosphere leads to (external) reconnection and the formation of ''standard'' jets with an inverse Y-shaped configuration. Eventually, low-atmosphere (internal) reconnection of sheared fieldlines in the emerging flux region produces an erupting magnetic flux rope and a reconnection jet underneath it. The erupting plasma blows out the ambient field and, moreover, it unwinds as it is ejected into the outer solar atmosphere. The fast emission of the cool material that erupts together with the hot outflows due to external/internal reconnection form a wider ''blowout'' jet. We show the transition from ''standard'' to ''blowout'' jets and report on their 3D structure. The physical plasma properties of the jets are consistent with observational studies.

  6. Validation of Numerical Two-Fluid and Kinetic Plasma Models

    SciTech Connect (OSTI)

    Daniel Barnes

    2011-03-25

    This was a four year grant commencing October 1, 2003 and finishing September 30, 2007. The funding was primarily used to support the work of the Principal Investigator, who collaborated with Profs. Scott Parker and John Cary at U. Colorado, and with two students, N. Xiang and J. Cheng also of U. Colorado. The technical accomplishments of this grant can be found in the publications listed in the final Section here. The main accomplishments of the grant work were: (1) Development and implementation of time-implicit two-fluid simulation methods in collaboration with the NIMROD team; and (2) Development and testing of a new time-implicit delta-f, energy-conserving method The basic two-fluid method, with many improvements is used in present NIMROD calculations. The energy-conserving delta-f method is under continuing development under contract between Coronado Consulting, a New Mexico sole proprietorship and the Oak Ridge National Laboratory.

  7. Numerical Modeling At Coso Geothermal Area (2007) | Open Energy...

    Open Energy Info (EERE)

    the Coso Range. Notes A finite element analysis is used to establish the 3D state of stress within the tectonic setting of the Coso Range. The mean and differential stress...

  8. Numerical Modeling At Coso Geothermal Area (1999) | Open Energy...

    Open Energy Info (EERE)

    microseismic travel time data Lees, J.M.; Wu, H. (1 August 1999) P wave anisotropy, stress, and crack distribution at Coso geothermal field, California Additional References...

  9. Analytical-Numerical Modeling Of Komatiite Lava Emplacement And...

    Open Energy Info (EERE)

    submarine tuff. Flow distances must have been long (tens to hundreds of kilometers) and flow volumes must have been very high (hundreds to thousands of km3). Lava...

  10. Numerical Modeling of Transient Basin and Range Extensional Geothermal...

    Open Energy Info (EERE)

    R. McKenna and D. D. Blackwell Conference PROCEEDINGS, TOUGH Symposium 2003 Lawrence Berkeley National Laboratory; Berkeley, California; 2003 Published PROCEEDINGS, TOUGH...

  11. A Hydro-Thermo-Mechanical Numerical Model For Hdr Geothermal...

    Open Energy Info (EERE)

    on the fluid flow distribution in an HDR geothermal reservoir. Authors T. W. Hicks, R. J. Pine, J. Willis-Richards, S. Xu, A. J. Jupe and N. E. V. Rodrigues Published Journal...

  12. Hydrodynamic and numerical modeling of a spherical homogeneous.pdf

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

  13. Numerical Modeling of WECS at Ecole Centrale de Nantes

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    LHEEA Lab. : Hydrodynamics, Energetics and Atmospheric Environment Staff : 100, Director : ... windwave coupling, freak wave statistics, wave propagation over complex ...

  14. Numerical method for shear bands in ductile metal with inclusions

    SciTech Connect (OSTI)

    Plohr, Jee Yeon N [Los Alamos National Laboratory; Plohr, Bradley J [Los Alamos National Laboratory

    2010-01-01

    A numerical method for mesoscale simulation of high strain-rate loading of ductile metal containing inclusions is described. Because of small-scale inhomogeneities, such a composite material is prone to localized shear deformation (adiabatic shear bands). The modeling framework is the Generalized Method of Cells of Paley and Aboudi [Mech. Materials, vol. 14, pp. /27-139, 1992], which ensures that the micromechanical response of the material is reflected in the behavior of the composite at the mesoscale. To calculate the effective plastic strain rate when shear bands are present, the analytic and numerical analysis of shear bands by Glimm, Plohr, and Sharp [Mech. Materials, vol. 24, pp. 31-41, 1996] is adapted and extended.

  15. Numerical thermalization in particle-in-cell simulations with Monte-Carlo collisions

    SciTech Connect (OSTI)

    Lai, P. Y.; Lin, T. Y.; Lin-Liu, Y. R.; Chen, S. H.

    2014-12-15

    Numerical thermalization in collisional one-dimensional (1D) electrostatic (ES) particle-in-cell (PIC) simulations was investigated. Two collision models, the pitch-angle scattering of electrons by the stationary ion background and large-angle collisions between the electrons and the neutral background, were included in the PIC simulation using Monte-Carlo methods. The numerical results show that the thermalization times in both models were considerably reduced by the additional Monte-Carlo collisions as demonstrated by comparisons with Turner's previous simulation results based on a head-on collision model [M. M. Turner, Phys. Plasmas 13, 033506 (2006)]. However, the breakdown of Dawson's scaling law in the collisional 1D ES PIC simulation is more complicated than that was observed by Turner, and the revised scaling law of the numerical thermalization time with numerical parameters are derived on the basis of the simulation results obtained in this study.

  16. Modeling

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    Modeling & Analysis, News, News & Events, Photovoltaic, Renewable Energy, Research & Capabilities, Solar, Solar Newsletter, SunShot, Systems Analysis Sandia Develops Stochastic ...

  17. Modeling

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    Monte Carlo modeling it was found that for noisy signals with a significant background component, accuracy is improved by fitting the total emission data which includes the...

  18. Modeling

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    Science and Actuarial Practice" Read More Permalink New Project Is the ACME of Computer Science to Address Climate Change Analysis, Climate, Global Climate & Energy, Modeling, ...

  19. Modeling

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    Solar Sandia Labs Releases New Version of PVLib Toolbox Sandia has released version 1.3 of PVLib, its widely used Matlab toolbox for modeling photovoltaic (PV) power ...

  20. Modeling

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    ... Sandia Will Host PV Bankability Workshop at Solar Power International (SPI) 2013 Computational Modeling & Simulation, Distribution Grid Integration, Energy, Facilities, Grid ...

  1. Modeling

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    Though adequate for modeling mean transport, this approach does not address ... Microphysics such as diffusive transport and chemical kinetics are represented by ...

  2. Modeling of geothermal systems

    SciTech Connect (OSTI)

    Bodvarsson, G.S.; Pruess, K.; Lippmann, M.J.

    1985-03-01

    During the last decade the use of numerical modeling for geothermal resource evaluation has grown significantly, and new modeling approaches have been developed. In this paper we present a summary of the present status in numerical modeling of geothermal systems, emphasizing recent developments. Different modeling approaches are described and their applicability discussed. The various modeling tasks, including natural-state, exploitation, injection, multi-component and subsidence modeling, are illustrated with geothermal field examples. 99 refs., 14 figs.

  3. Microsoft PowerPoint - 1-A-3-Beaudry-LOSIQUE-Brief to Blade Workshop...

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    ... Energy Capture Efficiency Tower 2020 LCOE Target FY10 Budget for COE reductions ... Modeled cost of energy to reflect impact of DOE R&D; lower than actual market prices Rotor ...

  4. PowerPoint Presentation

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    ... This suggests that any of the three model options are suitable for estimates of a PV system's output or cost. Annual Output (kWhyr) LCOE (kWh) Sandia Module 360,000 11.8 CEC ...

  5. Numerical prediction of window condensation potential

    SciTech Connect (OSTI)

    McGowan, A.G.

    1995-08-01

    Although a substantial amount of effort has been expended to develop numerical methods for determining windows U-factors (EE 1983; Goss and Curcija 1994; Standaert 1985; CSA 1993a; NFRC 1991), there has been little work to data on using numerical methods to predict condensation potential. It is perhaps of direct interest to most ASHRAE members to determine heat loss and solar gains through windows as a precursor to sizing heating and cooling equipment, but condensation has long been recognized as an extremely important issue for consumers (and, consequently, for window manufacturers). Moreover, building scientists recognize the link between condensation and increased energy consumption (due to latent loads), reduced occupant comfort and indoor air quality (from the presence of bacteria and mold), and structural damage (where accumulated condensation is absorbed by the building material, thus reducing their structural stability). The National Fenestration Rating Council (NFRC) is developing a rating method for condensation potential in fenestration products as part of its mandate from the US Department of Energy (DOE). A rating method would benefit from the use of simulation as a supplement to physical condensation resistance testing to reduce the cost and time required for implementation and increase the flexibility of the rating method. This paper outlines one of the necessary components in the application of numerical methods for evaluating condensation in fenestration products. The theoretical approach and its practical application are discussed, as well as some comparisons between numerical prediction and physical test results for a sample of products.

  6. Numerical likelihood analysis of cosmic ray anisotropies

    SciTech Connect (OSTI)

    Carlos Hojvat et al.

    2003-07-02

    A numerical likelihood approach to the determination of cosmic ray anisotropies is presented which offers many advantages over other approaches. It allows a wide range of statistically meaningful hypotheses to be compared even when full sky coverage is unavailable, can be readily extended in order to include measurement errors, and makes maximum unbiased use of all available information.

  7. Numerical and experimental analysis of a retrievable offshore loading facility

    SciTech Connect (OSTI)

    Sterndorff, M.J.; O`Brien, P.

    1995-12-31

    ROLF (Retrievable Offshore Loading Facility) has been proposed as an alternative offshore oil export tanker loading system for the North Sea. The system consists of a flexible riser ascending from the seabed in a lazy wave configuration to the bow of a dynamically positioned tanker. In order to supplant and support the numerical analyses performed to design the system, an extensive model test program was carried out in a 3D offshore basin at scale 1:50. A model riser with properties equivalent to the properties of the oil filled prototype riser installed in seawater was tested in several combinations of waves and current. During the tests the forces at the bow of the tanker and at the pipeline end manifold were measured together with the motions of the tanker and the riser. The riser motions were measured by means of a video based 3D motion monitoring system. Of special importance was accurate determination of the minimum bending radius for the riser. This was derived based on the measured riser motions. The results of the model tests were compared to numerical analyses by an MCS proprietary riser analysis program.

  8. Error Estimation for Fault Tolerance in Numerical Integration...

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    Error Estimation for Fault Tolerance in Numerical Integration Solvers Event Sponsor: ... In numerical integration solvers, approximation error can be estimated at a low cost. We ...

  9. Numerical Verification of Bounce Harmonic Resonances in Neoclassical

    Office of Scientific and Technical Information (OSTI)

    for Tokamaks Kimin Kim, Jong-Kyu Park and Allen H. Boozer 70 PLASMA PHYSICS AND FUSION TECHNOLOGY Tokamaks, Numerical Verification Tokamaks, Numerical Verification This...

  10. Development of Numerical Simulation Capabilities for In Situ...

    Office of Scientific and Technical Information (OSTI)

    Development of Numerical Simulation Capabilities for In Situ Heating of Oil Shale Citation Details In-Document Search Title: Development of Numerical Simulation Capabilities for In ...

  11. Stress-dependent permeability of fractured rock masses: A numerical...

    Office of Scientific and Technical Information (OSTI)

    permeability of fractured rock masses: A numerical study Citation Details In-Document Search Title: Stress-dependent permeability of fractured rock masses: A numerical study We ...

  12. Toward portable programming of numerical linear algebra on manycore...

    Office of Scientific and Technical Information (OSTI)

    Toward portable programming of numerical linear algebra on manycore nodes. Citation Details In-Document Search Title: Toward portable programming of numerical linear algebra on ...

  13. Modeling

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    with application in modeling NDCX-II experiments Wangyi Liu 1 , John Barnard 2 , Alex Friedman 2 , Nathan Masters 2 , Aaron Fisher 2 , Alice Koniges 2 , David Eder 2 1 LBNL, USA, 2...

  14. Modeling

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    NASA Earth at Night Video EC, Energy, Energy Efficiency, Global, Modeling, News & Events, Solid-State Lighting, Videos NASA Earth at Night Video Have you ever wondered what the ...

  15. On the Numerical Dispersion of Electromagnetic Particle-In-Cell Code : Finite Grid Instability

    SciTech Connect (OSTI)

    Meyers, Michael David; Huang, Chengkun; Zeng, Yong; Yi, Sunghwan; Albright, Brian James

    2014-07-15

    The Particle-In-Cell (PIC) method is widely used in relativistic particle beam and laser plasma modeling. However, the PIC method exhibits numerical instabilities that can render unphysical simulation results or even destroy the simulation. For electromagnetic relativistic beam and plasma modeling, the most relevant numerical instabilities are the finite grid instability and the numerical Cherenkov instability. We review the numerical dispersion relation of the electromagnetic PIC algorithm to analyze the origin of these instabilities. We rigorously derive the faithful 3D numerical dispersion of the PIC algorithm, and then specialize to the Yee FDTD scheme. In particular, we account for the manner in which the PIC algorithm updates and samples the fields and distribution function. Temporal and spatial phase factors from solving Maxwell's equations on the Yee grid with the leapfrog scheme are also explicitly accounted for. Numerical solutions to the electrostatic-like modes in the 1D dispersion relation for a cold drifting plasma are obtained for parameters of interest. In the succeeding analysis, we investigate how the finite grid instability arises from the interaction of the numerical 1D modes admitted in the system and their aliases. The most significant interaction is due critically to the correct representation of the operators in the dispersion relation. We obtain a simple analytic expression for the peak growth rate due to this interaction.

  16. Numerical simulation of scattering of acoustic waves by inelastic bodies using hypersingular boundary integral equation

    SciTech Connect (OSTI)

    Daeva, S.G.; Setukha, A.V.

    2015-03-10

    A numerical method for solving a problem of diffraction of acoustic waves by system of solid and thin objects based on the reduction the problem to a boundary integral equation in which the integral is understood in the sense of finite Hadamard value is proposed. To solve this equation we applied piecewise constant approximations and collocation methods numerical scheme. The difference between the constructed scheme and earlier known is in obtaining approximate analytical expressions to appearing system of linear equations coefficients by separating the main part of the kernel integral operator. The proposed numerical scheme is tested on the solution of the model problem of diffraction of an acoustic wave by inelastic sphere.

  17. Hydrologic Modeling Capabilities

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    Understanding complex hydrologic systems requires the ability to develop, utilize, and interpret both numerical and analytical models. The Defense Waste Management Programs has both experience and technical knowledge to use and develop Earth systems models. Hydrological Modeling Models are simplified representations of reality, which we accept do not capture every detail of reality. Mathematical and numerical models can be used to rigorously test geologic and hydrologic assumptions, determine

  18. Numerical calculation of two-phase turbulent jets

    SciTech Connect (OSTI)

    Saif, A.A.

    1995-05-01

    Two-phase turbulent round jets were numerically simulated using a multidimensional two-phase CFD code based on the two-fluid model. The turbulence phenomena were treated with the standard k-{epsilon} model. It was modified to take into account the additional dissipation of turbulent kinetic energy by the dispersed phase. Within the context of the two-fluid model it is more appropriate and physically justified to treat the diffusion by an interfacial force in the momentum equation. In this work, the diffusion force and the additional dissipation effect by the dispersed phase were modeled starting from the classical turbulent energy spectrum analysis. A cut-off frequency was proposed to decrease the dissipation effect by the dispersed phase when large size particles are introduced in the flow. The cut-off frequency combined with the bubble-induced turbulence effect allows for an increase in turbulence for large particles. Additional care was taken in choosing the right kind of experimental data from the literature so that a good separate effect test was possible for their models. The models predicted the experimental data very closely and they were general enough to predict extreme limit cases: water-bubble and air-droplet jets.

  19. Advanced numerical methods in mesh generation and mesh adaptation

    SciTech Connect (OSTI)

    Lipnikov, Konstantine; Danilov, A; Vassilevski, Y; Agonzal, A

    2010-01-01

    Numerical solution of partial differential equations requires appropriate meshes, efficient solvers and robust and reliable error estimates. Generation of high-quality meshes for complex engineering models is a non-trivial task. This task is made more difficult when the mesh has to be adapted to a problem solution. This article is focused on a synergistic approach to the mesh generation and mesh adaptation, where best properties of various mesh generation methods are combined to build efficiently simplicial meshes. First, the advancing front technique (AFT) is combined with the incremental Delaunay triangulation (DT) to build an initial mesh. Second, the metric-based mesh adaptation (MBA) method is employed to improve quality of the generated mesh and/or to adapt it to a problem solution. We demonstrate with numerical experiments that combination of all three methods is required for robust meshing of complex engineering models. The key to successful mesh generation is the high-quality of the triangles in the initial front. We use a black-box technique to improve surface meshes exported from an unattainable CAD system. The initial surface mesh is refined into a shape-regular triangulation which approximates the boundary with the same accuracy as the CAD mesh. The DT method adds robustness to the AFT. The resulting mesh is topologically correct but may contain a few slivers. The MBA uses seven local operations to modify the mesh topology. It improves significantly the mesh quality. The MBA method is also used to adapt the mesh to a problem solution to minimize computational resources required for solving the problem. The MBA has a solid theoretical background. In the first two experiments, we consider the convection-diffusion and elasticity problems. We demonstrate the optimal reduction rate of the discretization error on a sequence of adaptive strongly anisotropic meshes. The key element of the MBA method is construction of a tensor metric from hierarchical edge

  20. Technical progress report for application of numerical simulation methodology to automotive combustion

    SciTech Connect (OSTI)

    1980-12-01

    The second quarterly technical progress report is presented for a program entitled, Application of Numerical Simulation Methodology to Automotive Combustion. The goal of the program is to develop, validate, demonstrate and apply a numerical simulation methodology for in-cylinder reactive flows in internal combustion engines. Previous work on this contract involved the initial development and validation of a finite difference based simulation model for time dependent axisymmetric flows which includes: a generalized coordinate system for arbitrary mesh design and treatment of complex and time dependent boundaries; multiple and interacting chemical species; coupled swirl flow velocity component; and two-equation turbulence closure. In its various stages of development, the model has been used to simulate numerous engine-related problems for validation and demonstration purposes. The technical effort during the current reporting period has concentrated on: reactive flow model development, test and data comparison studies; swirl flow simulations; and in-cylinder compression cycle flow simulations. Results of these studies are discussed.

  1. Effect of virtual mass on the characteristics and the numerical stability in two-phase flows

    SciTech Connect (OSTI)

    No, H.C.; Kazimi, M.S.

    1981-04-01

    It is known that the typical six equation two-fluid model of the two-phase flow possesses complex characteristics, exhibits unbounded instabilities in the short-wavelength limit and constitutes an ill-posed initial value problem. Among the suggestions to overcome these difficulties, one model for the virtual mass force terms were studied here, because the virtual mass represents real physical effects to accomplish the dissipation for numerical stability. It was found that the virtual mass has a profound effect upon the mathematical characteristic and numerical stability. Here a quantitative bound on the coefficient of the virtual mass terms was suggested for mathematical hyperbolicity and numerical stability. It was concluded that the finite difference scheme with the virtual mass model is restricted only by the convective stability conditions with the above suggested value.

  2. Modeling

    SciTech Connect (OSTI)

    Loth, E.; Tryggvason, G.; Tsuji, Y.; Elghobashi, S. E.; Crowe, Clayton T.; Berlemont, A.; Reeks, M.; Simonin, O.; Frank, Th; Onishi, Yasuo; Van Wachem, B.

    2005-09-01

    Slurry flows occur in many circumstances, including chemical manufacturing processes, pipeline transfer of coal, sand, and minerals; mud flows; and disposal of dredged materials. In this section we discuss slurry flow applications related to radioactive waste management. The Hanford tank waste solids and interstitial liquids will be mixed to form a slurry so it can be pumped out for retrieval and treatment. The waste is very complex chemically and physically. The ARIEL code is used to model the chemical interactions and fluid dynamics of the waste.

  3. Numerical simulation of tectonic plates motion and seismic process in Central Asia

    SciTech Connect (OSTI)

    Peryshkin, A. Yu.; Makarov, P. V. Eremin, M. O.

    2014-11-14

    An evolutionary approach proposed in [1, 2] combining the achievements of traditional macroscopic theory of solid mechanics and basic ideas of nonlinear dynamics is applied in a numerical simulation of present-day tectonic plates motion and seismic process in Central Asia. Relative values of strength parameters of rigid blocks with respect to the soft zones were characterized by the ? parameter that was varied in the numerical experiments within ? = 1.11.8 for different groups of the zonal-block divisibility. In general, the numerical simulations of tectonic block motion and accompanying seismic process in the model geomedium indicate that the numerical solutions of the solid mechanics equations characterize its deformation as a typical behavior of a nonlinear dynamic system under conditions of self-organized criticality.

  4. Advanced numerical methods for three dimensional two-phase flow calculations

    SciTech Connect (OSTI)

    Toumi, I.; Caruge, D.

    1997-07-01

    This paper is devoted to new numerical methods developed for both one and three dimensional two-phase flow calculations. These methods are finite volume numerical methods and are based on the use of Approximate Riemann Solvers concepts to define convective fluxes versus mean cell quantities. The first part of the paper presents the numerical method for a one dimensional hyperbolic two-fluid model including differential terms as added mass and interface pressure. This numerical solution scheme makes use of the Riemann problem solution to define backward and forward differencing to approximate spatial derivatives. The construction of this approximate Riemann solver uses an extension of Roe`s method that has been successfully used to solve gas dynamic equations. As far as the two-fluid model is hyperbolic, this numerical method seems very efficient for the numerical solution of two-phase flow problems. The scheme was applied both to shock tube problems and to standard tests for two-fluid computer codes. The second part describes the numerical method in the three dimensional case. The authors discuss also some improvements performed to obtain a fully implicit solution method that provides fast running steady state calculations. Such a scheme is not implemented in a thermal-hydraulic computer code devoted to 3-D steady-state and transient computations. Some results obtained for Pressurised Water Reactors concerning upper plenum calculations and a steady state flow in the core with rod bow effect evaluation are presented. In practice these new numerical methods have proved to be stable on non staggered grids and capable of generating accurate non oscillating solutions for two-phase flow calculations.

  5. A hybrid numerical fluid dynamics code for resistive magnetohydrodynamics

    SciTech Connect (OSTI)

    2006-04-01

    Spasmos is a computational fluid dynamics code that uses two numerical methods to solve the equations of resistive magnetohydrodynamic (MHD) flows in compressible, inviscid, conducting media[1]. The code is implemented as a set of libraries for the Python programming language[2]. It represents conducting and non-conducting gases and materials with uncomplicated (analytic) equations of state. It supports calculations in 1D, 2D, and 3D geometry, though only the 1D configuation has received significant testing to date. Because it uses the Python interpreter as a front end, users can easily write test programs to model systems with a variety of different numerical and physical parameters. Currently, the code includes 1D test programs for hydrodynamics (linear acoustic waves, the Sod weak shock[3], the Noh strong shock[4], the Sedov explosion[5], magnetic diffusion (decay of a magnetic pulse[6], a driven oscillatory "wine-cellar" problem[7], magnetic equilibrium), and magnetohydrodynamics (an advected magnetic pulse[8], linear MHD waves, a magnetized shock tube[9]). Spasmos current runs only in a serial configuration. In the future, it will use MPI for parallel computation.

  6. Numerical simulation of multi-layered textile composite reinforcement forming

    SciTech Connect (OSTI)

    Wang, P.; Hamila, N.; Boisse, P.

    2011-05-04

    One important perspective in aeronautics is to produce large, thick or/and complex structural composite parts. The forming stage presents an important role during the whole manufacturing process, especially for LCM processes (Liquid Composites Moulding) or CFRTP (Continuous Fibre Reinforcements and Thermoplastic resin). Numerical simulations corresponding to multi-layered composite forming allow the prediction for a successful process to produce the thick parts, and importantly, the positions of the fibres after forming to be known. This paper details a set of simulation examples carried out by using a semi-discrete shell finite element made up of unit woven cells. The internal virtual work is applied on all woven cells of the element taking into account tensions, in-plane shear and bending effects. As one key problem, the contact behaviours of tool/ply and ply/ply are described in the numerical model. The simulation results not only improve our understanding of the multi-layered composite forming process but also point out the importance of the fibre orientation and inter-ply friction during formability.

  7. A hybrid numerical fluid dynamics code for resistive magnetohydrodynamics

    Energy Science and Technology Software Center (OSTI)

    2006-04-01

    Spasmos is a computational fluid dynamics code that uses two numerical methods to solve the equations of resistive magnetohydrodynamic (MHD) flows in compressible, inviscid, conducting media[1]. The code is implemented as a set of libraries for the Python programming language[2]. It represents conducting and non-conducting gases and materials with uncomplicated (analytic) equations of state. It supports calculations in 1D, 2D, and 3D geometry, though only the 1D configuation has received significant testing to date. Becausemore » it uses the Python interpreter as a front end, users can easily write test programs to model systems with a variety of different numerical and physical parameters. Currently, the code includes 1D test programs for hydrodynamics (linear acoustic waves, the Sod weak shock[3], the Noh strong shock[4], the Sedov explosion[5], magnetic diffusion (decay of a magnetic pulse[6], a driven oscillatory "wine-cellar" problem[7], magnetic equilibrium), and magnetohydrodynamics (an advected magnetic pulse[8], linear MHD waves, a magnetized shock tube[9]). Spasmos current runs only in a serial configuration. In the future, it will use MPI for parallel computation.« less

  8. Subtask 2.2 - Creating A Numerical Technique for Microseismic Data Inversion

    SciTech Connect (OSTI)

    Anastasia Dobroskok; Yevhen Holubnyak; James Sorensen

    2009-05-01

    Geomechanical and geophysical monitoring are the techniques which can complement each other and provide enhancement in the solutions of many problems of geotechnical engineering. One of the most promising geophysical techniques is passive seismic monitoring. The essence of the technique is recording the acoustic signals produced in the subsurface, either naturally or in response to human activity. The acoustic signals are produced by mechanical displacements on the contacts of structural elements (e.g., faults, boundaries of rock blocks, natural and induced fractures). The process can be modeled by modern numerical techniques developed in geomechanics. The report discusses a study that was aimed at the unification of the passive seismic monitoring and numerical modeling for the monitoring of the hydraulic fracture propagation. The approach adopted in the study consisted of numerical modeling of the seismicity accompanying hydraulic fracture propagation and defining seismic attributes and patterns characterizing the process and fracture parameters. Numerical experiments indicated that the spatial distribution of seismic events is correlated to geometrical parameters of hydrofracture. Namely, the highest density of the events is observed along fracture contour, and projection of the events to the fracture plane makes this effect most pronounced. The numerical experiments also showed that dividing the totality of the events into groups corresponding to the steps of fracture propagation allows for reconstructing the geometry of the resulting fracture more accurately than has been done in the majority of commercial applications.

  9. Numerical description of cavitation on axisymmetric bodies

    SciTech Connect (OSTI)

    Hickox, C.E.; Hailey, C.E.; Wolfe, W.P.; Watts, H.A.; Gross, R.J.; Ingber, M.S.

    1988-01-01

    This paper reports on ongoing studies which are directed toward the development of predictive techniques for the modeling of steady cavitation on axisymmetric bodies. The primary goal of the modeling effort is the prediction of cavity shape and pressure distribution from which forces and moments can be calculated. Here we present an overview of the modeling techniques developed and compare predictions with experimental data obtained from water tunnel tests for both limited and supercavitation. 14 refs., 4 figs.

  10. Open-source extreme conditions modeling tools

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    ... HomeEnergy, Modeling & Analysis, News, Renewable Energy, Water PowerOpen-source extreme ... numerical tools for use in modeling extreme conditions of wave energy converters (WECs). ...