National Library of Energy BETA

Sample records for lcoe numerical modeling

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

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

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

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

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

    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. Department...

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

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

    LBNLNREL Analysis Predicts Record Low LCOE for Wind Energy in 2012-2013 LBNLNREL Analysis Predicts Record Low LCOE for Wind Energy in 2012-2013 February 24, 2012 - 11:27am...

  4. MHK LCOE Reporting Guidance Draft | OpenEI Community

    Open Energy Info (EERE)

    AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on QA:QA J-E-1 SECTION J APPENDIX E LIST OFAMERICA'SHeavyAgency (IRENA) JumpLiterature Review HomeM-7MHK LCOE Reporting

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

    Office of Energy Efficiency and Renewable Energy (EERE)

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

  6. Supplementary Information Potential for Electricity Generation from Renewable Resources and Levelized Cost of Electricity (LCOE)

    E-Print Network [OSTI]

    Suo, Zhigang

    Supplementary Information Potential for Electricity Generation from Renewable Resources and Levelized Cost of Electricity (LCOE) Electrical energy can be generated from renewable resources the potential to meet the worldwide demand of electricity and they contribute to the total generation

  7. Numerical wind speed simulation model

    SciTech Connect (OSTI)

    Ramsdell, J.V.; Athey, G.F.; Ballinger, M.Y.

    1981-09-01

    A relatively simple stochastic model for simulating wind speed time series that can be used as an alternative to time series from representative locations is described in this report. The model incorporates systematic seasonal variation of the mean wind, its standard deviation, and the correlation speeds. It also incorporates systematic diurnal variation of the mean speed and standard deviation. To demonstrate the model capabilities, simulations were made using model parameters derived from data collected at the Hanford Meteorology Station, and results of analysis of simulated and actual data were compared.

  8. Numerical Model Construction with Closed Observables

    E-Print Network [OSTI]

    Felix Dietrich; Gerta Köster; Hans-Joachim Bungartz

    2015-10-18

    Performing analysis, optimization and control using simulations of many-particle systems is computationally demanding when no macroscopic model for the dynamics of the variables of interest is available. In case observations on the macroscopic scale can only be produced via legacy simulator code or live experiments, finding a model for these macroscopic variables is challenging. In this paper, we employ time-lagged embedding theory to construct macroscopic numerical models from output data of a black box, such as a simulator or live experiments. Since the state space variables of the constructed, coarse model are dynamically closed and observable by an observation function, we call these variables closed observables. The approach is an online-offline procedure, as model construction from observation data is performed offline and the new model can then be used in an online phase, independent of the original. We illustrate the theoretical findings with numerical models constructed from time series of a two-dimensional ordinary differential equation system, and from the density evolution of a transport-diffusion system. Applicability is demonstrated in a real-world example, where passengers leave a train and the macroscopic model for the density flow onto the platform is constructed with our approach. If only the macroscopic variables are of interest, simulation runtimes with the numerical model are three orders of magnitude lower compared to simulations with the original fine scale model. We conclude with a brief discussion of possibilities of numerical model construction in systematic upscaling, network optimization and uncertainty quantification.

  9. High performance computing and numerical modelling

    E-Print Network [OSTI]

    ,

    2014-01-01

    Numerical methods play an ever more important role in astrophysics. This is especially true in theoretical works, but of course, even in purely observational projects, data analysis without massive use of computational methods has become unthinkable. The key utility of computer simulations comes from their ability to solve complex systems of equations that are either intractable with analytic techniques or only amenable to highly approximative treatments. Simulations are best viewed as a powerful complement to analytic reasoning, and as the method of choice to model systems that feature enormous physical complexity such as star formation in evolving galaxies, the topic of this 43rd Saas Fee Advanced Course. The organizers asked me to lecture about high performance computing and numerical modelling in this winter school, and to specifically cover the basics of numerically treating gravity and hydrodynamics in the context of galaxy evolution. This is still a vast field, and I necessarily had to select a subset ...

  10. Lattice Boltzmann Model for Numerical Relativity

    E-Print Network [OSTI]

    Ilseven, E

    2015-01-01

    In the Bona-Masso formulation, Einstein equations are written as a set of flux conservative first order hyperbolic equations that resemble fluid dynamics equations. Based on this formulation, we construct a lattice Boltzmann model for Numerical Relativity. Our model is validated with well-established tests, showing good agreement with analytical solutions. Furthermore, we show that by increasing the relaxation time, we gain stability at the cost of losing accuracy, and by decreasing the lattice spacings while keeping a constant numerical diffusivity, the accuracy and stability of our simulations improves. Finally, in order to show the potential of our approach a linear scaling law for parallelisation with respect to number of CPU cores is demonstrated. Our model represents the first step in using lattice kinetic theory to solve gravitational problems.

  11. Numerical Modeling of Eastern Connecticut's Visual Resources1

    E-Print Network [OSTI]

    Numerical Modeling of Eastern Connecticut's Visual Resources1 Daniel L. Civco 2/ l Submitted Conservation, The University of Connecticut, Storrs, Connecticut. / Abstract: A numerical model capable of accurately predicting the preference for landscape photographs of selected points in eastern Connecticut

  12. 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.

  13. INCORPORATION OF GROUNDWATER FLOW INTO NUMERICAL MODELS AND DESIGN MODELS

    E-Print Network [OSTI]

    -coupled, ground-source heat pumps, groundwater, heat pump, heat exchanger, heat transfer, numerical models-loop ground-coupled heat exchangers. Green and Perry (1961) demonstrated that the value of effective thermal on the design and performance of vertical closed-loop ground heat exchangers. Based on the investigation results

  14. Contributions Motivation Literature Models Synergy Measure Numerical Examples Conclusions Supply Chain Network Models for Humanitarian

    E-Print Network [OSTI]

    Nagurney, Anna

    Contributions Motivation Literature Models Synergy Measure Numerical Examples Conclusions Supply Motivation Literature Models Synergy Measure Numerical Examples Conclusions 1 Contributions 2 Motivation 3 University of Massachusetts Amherst #12;Contributions Motivation Literature Models Synergy Measure Numerical

  15. NUMERICAL APPROXIMATION OF A MACROSCOPIC MODEL OF PEDESTRIAN FLOWS

    E-Print Network [OSTI]

    Chalons, Christophe

    NUMERICAL APPROXIMATION OF A MACROSCOPIC MODEL OF PEDESTRIAN FLOWS CHRISTOPHE CHALONS Abstract for the description of the flow of pedestrians. Solutions of the associated Riemann problem are known to be possibly in obtaining an efficient numerical scheme. Numerical evidences are proposed. Key words. macroscopic pedestrian

  16. Validation of Eddy-renewal model by numerical simulation

    E-Print Network [OSTI]

    Garbe, Christoph S.

    Validation of Eddy-renewal model by numerical simulation Li-Ping Hung1,2 , Christoph S. Garbe1, E-mail: wttsai@ncu.edu.tw Abstract. The eddy-renewal model proposes that the dominant vortical flows as the diffusivity effect in the thin diffusive sublayer beneath surface. Key Words: eddy-renewal model, numerical

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

    Open Energy Info (EERE)

    A Hydro-Thermo-Mechanical Numerical Model For Hdr Geothermal Reservoir Evaluation Jump to: navigation, search OpenEI Reference LibraryAdd to library Journal Article: A...

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

    Open Energy Info (EERE)

    Numerical Modeling of the Nucleation Conditions of Petal-Centerline Fractures below a Borehole Floor, A Sensitivity Study and Application to the Coso Geothermal Field Jump to:...

  19. 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 LaboratoryUniversity of...

  20. Starting laminar plumes: Comparison of laboratory and numerical modeling

    E-Print Network [OSTI]

    van Keken, Peter

    Starting laminar plumes: Comparison of laboratory and numerical modeling Judith Vatteville Institut, France (davaille@fast.u-psud.fr) [1] A detailed comparison of starting laminar plumes in viscous fluids. Davaille (2009), Starting laminar plumes: Comparison of laboratory and numerical modeling, Geochem. Geophys

  1. Evolution of the Bohemian Massif: Insights from numerical modeling

    E-Print Network [OSTI]

    Cerveny, Vlastislav

    Evolution of the Bohemian Massif: Insights from numerical modeling Petra Maierová Supervisor: Doc of Geophysics Faculty of Mathematics and Physics Charles University in Prague #12;February 4, 2013Evolution Conclusions Outline #12;February 4, 2013Evolution of the Bohemian Massif: Insights from numerical modeling 3

  2. NUMERICAL MODELING OF LOW FREQUENCY HYDRO-ACOUSTIC WAVES

    E-Print Network [OSTI]

    Kirby, James T.

    NUMERICAL MODELING OF LOW FREQUENCY HYDRO-ACOUSTIC WAVES GENERATED BY SUBMARINE TSUNAMIGENIC#al to increase the reliability of the system · Can we use precursors of tsunami? Hydro numerical models applicable on an oceanic scale #12;Index · Introduc#on on hydro

  3. Numerical models of phosphate esters in the Chattahoochee River

    E-Print Network [OSTI]

    Haffey, Samuel Fraad, 1973-

    2004-01-01

    A numerical model was constructed to assess the magnitude of organophosphoric acid triester sinks in the Chattahoochee River and to identify concentration patterns downstream of point source discharges. The model was built ...

  4. A numerical model simulation of longshore transport for Galveston Island 

    E-Print Network [OSTI]

    Gilbreath, Stephen Alexander

    1995-01-01

    The shoreline changes, deposition patterns, and longshore transport rates were calculated for the coast of Galveston Island using a numerical model simulation. The model only simulated changes due to waves creating longshore currents. East Beach...

  5. Category:Numerical Modeling | Open Energy Information

    Open Energy Info (EERE)

    AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on QA:QA J-E-1 SECTION JEnvironmentalBowerbank,CammackFLIR Jump to: navigation,GroundNear Infrared SurveysNumerical

  6. Data Assimilation for Idealised Mathematical Models of Numerical Weather Prediction

    E-Print Network [OSTI]

    Wirosoetisno, Djoko

    Data Assimilation for Idealised Mathematical Models of Numerical Weather Prediction Supervisors). Background: Numerical Weather Prediction (NWP) has seen significant gains in accuracy in recent years due in weather dynamics, e.g., the asymptotic balance seen in high and low pressure systems. Aims of the project

  7. Fully kinetic numerical modeling of a plasma thruster

    E-Print Network [OSTI]

    Szabo, James Joseph, 1969-

    2001-01-01

    A Hall effect plasma thruster with conductive acceleration channel walls was numerically modeled using 2D3V Particle-in-Cell (PIC) and Monte-Carlo Collision (MCC) methodolo- gies. Electron, ion, and neutral dynamics were ...

  8. Stress and diffusion induced interface motion: Modelling and numerical simulations

    E-Print Network [OSTI]

    Styles, Vanessa

    Stress and diffusion induced interface motion: Modelling and numerical simulations Harald Garcke of Mathematics, University of Sussex, Brighton, BN1 9QH, U.K. Abstract We propose a phase field model for stress stress effects. In this paper we will demonstrate that the model can also be used to describe other

  9. Numerical modelling of tsunami mitigation by mangroves Putu Harry Gunawan

    E-Print Network [OSTI]

    Mancini, Simona

    Numerical modelling of tsunami mitigation by mangroves Putu Harry Gunawan LAMA (Laboratoire d'Analyse et de Mathmatiques Appliques) UPEM putu-harry.gunawan@univ.paris-est.fr Abstract Figure 1: Mangrove-Tsunami Model. The role of mangroves (coastal forests) in the mitigation of tsunami impacts is a debated topic

  10. Diffusion and Dispersion Characterization of a Numerical Tsunami Model

    E-Print Network [OSTI]

    Tolkova, Elena

    Diffusion and Dispersion Characterization of a Numerical Tsunami Model David Burwell, Elena Tolkova and Arun Chawla February 6, 2007 1 Introduction Method Of Splitting Tsunami (MOST) is a depth averaged long wave tsunami inundation model that was originally developed by Titov and Synolakis (1995) for 1D

  11. EXPERIMENTAL CHARACTERIZATION AND NUMERICAL MODELING OF THERMAL AND ELECTROCHEMISTRY EFFECTS IN

    E-Print Network [OSTI]

    Kassegne, Samuel Kinde

    EXPERIMENTAL CHARACTERIZATION AND NUMERICAL MODELING OF THERMAL AND ELECTROCHEMISTRY EFFECTS IN 3D and Numerical Modeling of Thermal and Electrochemistry Effects in 3D Bionanoelectronics Platform by Neha

  12. 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...

  13. Numerical

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

    AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Homesum_a_epg0_fpd_mmcf_m.xls" ,"Available from WebQuantityBonneville Power Administration wouldMass map shinesSolarNewsusceptometer under pressureNavyNumerical simulations of current

  14. Small Buoys for Energy Harvesting : Experimental and Numerical Modeling Studies

    E-Print Network [OSTI]

    Grilli, Stéphan T.

    Small Buoys for Energy Harvesting : Experimental and Numerical Modeling Studies St´ephan T. Grilli for wave energy harvesting (free-floating or slackly moored), to produce about 1 KW per unit at full scale-contained (water tight) resonating multiple-spar buoy (or Starspar), in which a longer central spar houses the LEG

  15. Thermodynamics cycle analysis and numerical modeling of thermoelastic cooling systems

    E-Print Network [OSTI]

    Rubloff, Gary W.

    Thermodynamics cycle analysis and numerical modeling of thermoelastic cooling systems Suxin Qian level. However, a thermoelastic cooling system integrated with heat transfer fluid loops have not been;2012) (a.k.a. elastocaloric cooling). These solid-state cooling systems offer us alternatives to eliminate

  16. Small Buoys for Energy Harvesting : Experimental and Numerical Modeling Studies

    E-Print Network [OSTI]

    Grilli, Stéphan T.

    Small Buoys for Energy Harvesting : Experimental and Numerical Modeling Studies St´ephan T. Grilli, Inc., Cranston, RI, USA ABSTRACT We present the development and application of small buoy systems climates, of two new types of buoy systems equipped with an embedded Linear Electric Generator (LEG; made

  17. Numerical Modeling of Periodic Composite Media for Electromagnetic Shielding Application

    E-Print Network [OSTI]

    Koledintseva, Marina Y.

    on a conventional mixing theory, have served as the fundamentals for these techniques. In these formulationsNumerical Modeling of Periodic Composite Media for Electromagnetic Shielding Application Dagang Wu-difference time-domain (FDTD) method. The results are compared with conventional mixing theories and 3D Fourier

  18. Direct Numerical Simulations and Modeling of Jets in Crossflow

    E-Print Network [OSTI]

    Mahesh, Krishnan

    Direct Numerical Simulations and Modeling of Jets in Crossflow A THESIS SUBMITTED TO THE FACULTY. i #12;To my parents and my grandparents, and to Ramnath ii #12;Abstract Jets in crossflow are used to study the different aspects of round jets in a crossflow. The first problem studies

  19. Validation of Erosion Modeling: Physical and Numerical Mehrad Kamalzare1

    E-Print Network [OSTI]

    Franklin, W. Randolph

    -3590 ABSTRACT The overall intent of this research is to develop numerical models of erosion of levees, dams is necessary for emergency plans for levee or dam breaches. Griffis, 2007 addressed the overall design and hydraulic shear stress. #12;2 Xu and Zhang (2009) found that in addition to soil type, the degree

  20. Integrating Numerical Computation into the Modeling Instruction Curriculum

    E-Print Network [OSTI]

    Caballero, Marcos D; Aiken, John M; Douglas, Scott S; Scanlon, Erin M; Thoms, Brian; Schatz, Michael F

    2012-01-01

    We describe a way to introduce physics high school students with no background in programming to computational problem-solving experiences. Our approach builds on the great strides made by the Modeling Instruction reform curriculum. This approach emphasizes the practices of "Developing and using models" and "Computational thinking" highlighted by the NRC K-12 science standards framework. We taught 9th-grade students in a Modeling-Instruction-based physics course to construct computational models using the VPython programming environment. Numerical computation within the Modeling Instruction curriculum provides coherence among the curriculum's different force and motion models, links the various representations which the curriculum employs, and extends the curriculum to include real-world problems that are inaccessible to a purely analytic approach.

  1. On numerical considerations for modeling reactive astrophysical shocks

    SciTech Connect (OSTI)

    Papatheodore, Thomas L.; Messer, O. E. Bronson, E-mail: tpapathe@utk.edu, E-mail: bronson@ornl.gov [Department of Physics and Astronomy, University of Tennessee, Knoxville, TN 37996 (United States)

    2014-02-10

    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. 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. 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 models of Type Ia supernovae, including potential impacts on observables.

  2. A high-resolution, cloud-assimilating numerical weather prediction model for solar irradiance forecasting

    E-Print Network [OSTI]

    Mathiesen, Patrick; Collier, Craig; Kleissl, Jan

    2013-01-01

    of numerical weather prediction solar irradiance forecasts numerical weather prediction model for solar irradiance weather prediction for intra?day solar  forecasting in the 

  3. 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 resolution real-time model simulations possible, the need to obtain observations to both initialize numerical models and verify their output has become increasingly important. The assimilation of high resolution radar observations therefore provides a vital component in the development and utility of numerical model forecasts for both weather forecasting and contaminant transport, including future opportunities to improve wet deposition computations explicitly.

  4. 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 generation atomistic-to-continuum multiscale simulations. In addition, a rigorous studyof nite element discretizations of peridynamics will be considered. Using the fact that peridynamics is spatially derivative free, we will also characterize the space of admissible peridynamic solutions and carry out systematic analyses of the models, in particular rigorously showing how peridynamics encompasses fracture and other failure phenomena. Additional aspects of the project include the mathematical and numerical analysis of peridynamics applied to stochastic peridynamics models. In summary, the project will make feasible mathematically consistent multiscale models for the analysis and design of advanced materials.

  5. 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.

  6. 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).

  7. 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.

  8. 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.

  9. Numerical model of Electron Cyclotron Resonance Ion Source

    E-Print Network [OSTI]

    Mironov, V; Bondarchenko, A; Efremov, A; Loginov, V

    2015-01-01

    Important features of Electron Cyclotron Resonance Ion Source (ECRIS) operation are accurately reproduced with a numerical code. The code uses the particle-in-cell technique to model a dynamics of ions in ECRIS plasma. It is shown that gas dynamical ion confinement mechanism is sufficient to provide the ion production rates in ECRIS close to the experimentally observed values. Extracted ion currents are calculated and compared to the experiment for few sources. Changes in the extracted ion currents are obtained with varying the gas flow into the source chamber and the microwave power. Empirical scaling laws for ECRIS design are studied and the underlying physical effects are discussed.

  10. Numerical Modeling of Submarine Mass-Movement Generated Waves Using RANS Model

    E-Print Network [OSTI]

    Yim, Solomon C.

    Numerical Modeling of Submarine Mass-Movement Generated Waves Using RANS Model D. YUK and S. YIM model for predicting waves generated by nearshore submarine mass-movements is described. The model to coastal scientists and engineers for many years. Coastal wave generation due to submarine mass movement

  11. 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.

  12. Fracture Density Estimation Using Spectral Analysis of Reservoir Reflections: A Numerical Modeling Approach

    E-Print Network [OSTI]

    Pearce, Fred

    2003-01-01

    We use a 3-D finite difference numerical model to generate synthetic seismograms from a simple fractured reservoir

  13. Numeric modeling of fire suppression by organophosphorous inhibitors

    E-Print Network [OSTI]

    Makhviladze, G M; Zykov, A P

    2008-01-01

    Numerical calculations of the effect of organophosphorous inhibitor (CF3CH2O)3P and its mixtures with carbon dioxide on propane flames are carried out using the three dimensional Reynolds-averaged Navier-Stokes (RANS) equations in the low Mach number approximation. The k-e model of turbulence, the EDC combustion model and the weighted-sum-of-gray-gases model of radiation are used. The Westbrook global-kinetic scheme with fractional order of reaction was used for the calculation of chemical reaction rate of propane combustion. The empirical expression for the correction factor for the chemical reaction rate was used to model the effect of organophosphorous inhibitor no the reaction. Two series of test calculations for different values of the correction factor are carried out. Dependences of the minimum extinguishing concentration of the inhibitor per carbon dioxide volume concentration in the extinguishing mixtures were obtained. The results of test calculations are shown to agree reasonably with the experimen...

  14. A numerical model of aerosol scavenging: Part 1, Microphysics parameterization

    SciTech Connect (OSTI)

    Molenkamp, C.R.; Bradley, M.M.

    1991-09-01

    We have developed a three-dimensional numerical model (OCTET) to simulate the dynamics and microphysics of clouds and the transport, diffusion and precipitation scavenging of aerosol particles. In this paper we describe the cloud microphysics and scavenging parameterizations. The representation of cloud microphysics is a bulk- water parameterization which includes water vapor and five types of hydrometeors (cloud droplets, rain drops, ice crystals, snow, and graupel). A parallel parameterization represents the scavenging interactions between pollutant particles and hydrometeors including collection of particles because of condensation nucleation, Brownian and phoretic attachment, and inertial capture, resuspension because of evaporation and sublimation; and transfer interactions where particles collected by one type of hydrometeor are transferred to another type of freezing, melting, accretion, riming and autoconversion.

  15. Numerical Methods for the Bogoliubov-Tolmachev-Shirkov model in superconductivity theory

    E-Print Network [OSTI]

    Zhihao Ge; Ruihua Li

    2015-03-08

    In the work, the numerical methods are designed for the Bogoliubov-Tolmachev-Shirkov model in superconductivity theory. The numerical methods are novel and effective to determine the critical transition temperature and approximate to the energy gap function of the above model. Finally, a numerical example confirming the theoretical results is presented.

  16. NUMERICAL SOLUTION OF RESERVOIR FLOW MODELS BASED ON LARGE TIME STEP OPERATOR SPLITTING ALGORITHMS

    E-Print Network [OSTI]

    NUMERICAL SOLUTION OF RESERVOIR FLOW MODELS BASED ON LARGE TIME STEP OPERATOR SPLITTING ALGORITHMS. Special focus is posed on the numerical solution algorithms for the saturation equation, which. The general background for the reservoir ow model is reviewed, and the main features of the numerical

  17. NUMERICAL SOLUTION OF RESERVOIR FLOW MODELS BASED ON LARGE TIME STEP OPERATOR SPLITTING ALGORITHMS

    E-Print Network [OSTI]

    NUMERICAL SOLUTION OF RESERVOIR FLOW MODELS BASED ON LARGE TIME STEP OPERATOR SPLITTING ALGORITHMS focus is posed on the numerical solution algorithms for the saturation equation, which is a convection­ eral background for the reservoir flow model is reviewed, and the main features of the numerical

  18. Passive millimeter-wave retrieval of global precipitation utilizing satellites and a numerical weather prediction model

    E-Print Network [OSTI]

    Surussavadee, Chinnawat

    2007-01-01

    This thesis develops and validates the MM5/TBSCAT/F([lambda]) model, composed of a mesoscale numerical weather prediction (NWP) model (MM5), a two-stream radiative transfer model (TBSCAT), and electromagnetic models for ...

  19. Numerical study of energy diffusion in King models

    E-Print Network [OSTI]

    Tom Theuns

    1995-11-07

    The energy diffusion coefficients D_n(E) (n=1,2) for a system of equal mass particles moving self-consistently in an N-body realisation of a King model are computed from the probability per unit time, P(E, Delta E), that a star with initial energy E will undergo an energy change Delta E. In turn, P is computed from the number of times during the simulation that a particle in a state of given energy undergoes a transition to another state. These particle states are defined directly from the time evolution of E by identifying them with the event occuring between two local maxima in the E(t) curve. If one assumes next that energy changes are uncorrelated between different states, one can use diffusion theory to compute D_n(E). The simulations employ N=512, 2048,... , 32768 particles and are performed using an implementation of Aarseth's direct integrator N-body1 on a massively parallel computer. The more than seven million transitions measured in the largest N simulation provide excellent statistics. The numerically determined D(E)'s are compared against their theoretical counterparts which are computed from phase-space averaged rates of energy change due to independent binary encounters. The overall agreement between them is impressive over most of the energy range, notwithstanding the very different type of approximations involved, giving considerable support to the valid usage of these theoretical expressions to simulate dynamical evolution in Fokker-Planck type calculations.

  20. A regional numerical ocean model of the circulation in the Bay of Biscay

    E-Print Network [OSTI]

    Drijfhout, Sybren

    A regional numerical ocean model of the circulation in the Bay of Biscay Y. Friocourt,1,2,3 B Peninsula and in the Bay of Biscay is investigated by means of a regional ocean model. In particular numerical ocean model of the circulation in the Bay of Biscay, J. Geophys. Res., 112, C09008, doi:10

  1. Experimental and Numerical Investigations of Nonlinearity in soils using Advanced Laboratory-Scaled models

    E-Print Network [OSTI]

    ENINALS Experimental and Numerical Investigations of Nonlinearity in soils using Advanced of Nonlinearity in Soils Using Advanced Laboratory-scaled Models Concluding Workshop Ispra, May 28-30 2013 IFSTTAR-geological constraints - centrifuge modelling - numerical modelling MILESTONES - Experiencing techniques for realising

  2. NUMERICAL MODELLING OF IMPACT NOISE Cristina Daz Cereceda | Laboratori de Clcul Numric | Universitat Politcnica de Catalunya

    E-Print Network [OSTI]

    Huerta, Antonio

    NUMERICAL MODELLING OF IMPACT NOISE Cristina Díaz Cereceda | Laboratori de Càlcul Numèric Modelling of Impact Noise. In Congreso de Métodos Numéricos en Ingeniería 2009. [2] J. Hetherington, J. Poblet-Puig, A. Rodríguez-Ferran. Numerical Modelling of Impact Noise in Lightweight Floors: a Modal

  3. Numerical Modeling of CIGS Solar Cells: Definition of the Baseline and

    E-Print Network [OSTI]

    Sites, James R.

    Thesis Numerical Modeling of CIGS Solar Cells: Definition of the Baseline and Explanation our supervision by Markus Gloeckler entitled "Numerical Modeling of CIGS Solar Cells: Definition. A three-layer structure, simulating a Cu(InGa)Se2 (CIGS) heterojunction solar cell, was set up using

  4. NUMERICAL MODELING OF CIGS AND CdTe SOLAR CELLS: SETTING THE BASELINE

    E-Print Network [OSTI]

    Sites, James R.

    NUMERICAL MODELING OF CIGS AND CdTe SOLAR CELLS: SETTING THE BASELINE M. Gloeckler, A important complications that are often found in experimental CIGS and CdTe solar cells. 1. INTRODUCTION Numerical modeling of polycrystalline thin-film solar cells is an important strategy to test the viability

  5. Detection of microseismic compressional (P) body waves aided by numerical modeling of oceanic noise sources

    E-Print Network [OSTI]

    Stutzmann, Eléonore

    Detection of microseismic compressional (P) body waves aided by numerical modeling of oceanic noise), Detection of microseismic compressional (P) body waves aided by numerical modeling of oceanic noise sources in seismic noise, body waves present appealing properties but are still challenging to extract. Here we first

  6. CONVERGENCE OF NUMERICAL APPROXIMATIONS TO A PHASE FIELD BENDING ELASTICITY MODEL OF

    E-Print Network [OSTI]

    CONVERGENCE OF NUMERICAL APPROXIMATIONS TO A PHASE FIELD BENDING ELASTICITY MODEL OF MEMBRANE, Number 1, Pages 1­18 CONVERGENCE OF NUMERICAL APPROXIMATIONS TO A PHASE FIELD BENDING ELASTICITY MODEL OF MEMBRANE DEFORMATIONS QIANG DU AND XIAOQIANG WANG This paper is dedicated to Prof.Max Gunzburger

  7. Penetration of solar radiation in the upper ocean: A numerical model for oceanic and coastal waters

    E-Print Network [OSTI]

    Lee, Zhongping

    Penetration of solar radiation in the upper ocean: A numerical model for oceanic and coastal waters (2005), Penetration of solar radiation in the upper ocean: A numerical model for oceanic and coastal; Siegel et al., 1995] have demonstrated that the penetration of EVIS in the upper layer of the ocean plays

  8. Numerical modeling of fluid flow and time-lapse seismics to monitor CO2

    E-Print Network [OSTI]

    Santos, Juan

    Numerical modeling of fluid flow and time-lapse seismics to monitor CO2 Sequestration in aquifers J, ITALY). IMAL, 30/5/2014 Numerical modeling of fluid flow and time-lapse seismics to monitor CO2 Sequestration in aquifers ­ p. #12;Introduction. I Storage of CO2 in geological formations is a procedure

  9. Numerical Model of a Tensioner System and Flex Joint 

    E-Print Network [OSTI]

    Huang, Han

    2013-07-27

    D will be integrated into a numerical code, known as COUPLE, for the simulation of the dynamic interaction among the hull of a floating structure, such as SPAR or TLP, its mooring system and riser system under the impact of wind, current and waves...

  10. Efficient numerical algorithms for surface formulations of mathematical models for biomolecule analysis and design

    E-Print Network [OSTI]

    Bardhan, Jaydeep Porter, 1978-

    2006-01-01

    This thesis presents a set of numerical techniques that extend and improve computational modeling approaches for biomolecule analysis and design. The presented research focuses on surface formulations of modeling problems ...

  11. Mathematical Models and Numerical Solutions of Liquid-Solid and Solid-Liquid Phase Change

    E-Print Network [OSTI]

    Surana, Karan S.; Joy, Aaron; Quiros, Luis; Reddy, JN

    2015-04-01

    This paper presents numerical simulations of liquid-solid and solid-liquid phase change processes using mathematical models in Lagrangian and Eulerian descriptions. The mathematical models are derived by assuming a smooth ...

  12. Thermoelectric Generators for Automotive Waste Heat Recovery Systems Part I: Numerical Modeling

    E-Print Network [OSTI]

    Xu, Xianfan

    Thermoelectric Generators for Automotive Waste Heat Recovery Systems Part I: Numerical Modeling (TEG) designed for automotive waste heat recovery systems. This model is capable of computing telluride TEMs. Key words: Thermoelectric generators, waste heat recovery, automotive exhaust, skutterudites

  13. Fractional Calculus in Hydrologic Modeling: A Numerical Perspective

    SciTech Connect (OSTI)

    David A. Benson; Mark M. Meerschaert; Jordan Revielle

    2012-01-01

    Fractional derivatives can be viewed either as a handy extension of classical calculus or, more deeply, as mathematical operators defined by natural phenomena. This follows the view that the diffusion equation is defined as the governing equation of a Brownian motion. In this paper, we emphasize that fractional derivatives come from the governing equations of stable Levy motion, and that fractional integration is the corresponding inverse operator. Fractional integration, and its multi-dimensional extensions derived in this way, are intimately tied to fractional Brownian (and Levy) motions and noises. By following these general principles, we discuss the Eulerian and Lagrangian numerical solutions to fractional partial differential equations, and Eulerian methods for stochastic integrals. These numerical approximations illuminate the essential nature of the fractional calculus.

  14. LCOE | OpenEI Community

    Open Energy Info (EERE)

    AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on QA:QA J-E-1 SECTION J APPENDIX E LISTStar2-0057-EAInvervar HydroElectric Cooperative(Redirected from RenewableLCOLCOE

  15. OpenEI Community - LCOE

    Open Energy Info (EERE)

    AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on QA:QA J-E-1 SECTION J APPENDIXsourceII Jump to:InformationInformationOorjaen TheGeothermalCost Breakdown

  16. Turbulent Supersonic Channel Flow: Direct Numerical Simulation and Modeling

    E-Print Network [OSTI]

    Heinz, Stefan

    modeling: the turbulence frequency production mechanism, wall damping effects on turbulence model frequency production mechanisms and wall damping effects may be explained very well on the basis, Chik w = wall viscosity = kinematic viscosity, = T = turbulent kinematic viscosity, Ck d = pressure

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

    Open Energy Info (EERE)

    is needed to preserve the geochemical signature of the reservoir and . Finally, a new stress model is planned to be used for Dixie Valley, the model will utilize a boundary...

  18. Numerically Efficient Water Quality Modeling and Security Applications 

    E-Print Network [OSTI]

    Mann, Angelica

    2013-02-04

    . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13 1.4 Summary of contributions . . . . . . . . . . . . . . . . . . . . . . . . 17 2 THE MERLION WATER QUALITY MODEL . . . . . . . . . . . . . . . . 20 2.1 Model development . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20 2.2 Origin... tracking algorithm . . . . . . . . . . . . . . . . . . . . . . . . 24 2.3 Merlion water quality model . . . . . . . . . . . . . . . . . . . . . . . 31 2.4 Summary . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34 3 EFFICIENT MULTI...

  19. Numerical Simulation of Pulse-Tube Refrigerators: 1D model I.A. Lyulina1

    E-Print Network [OSTI]

    Eindhoven, Technische Universiteit

    of a piston, an aftercooler (AC), a regenerator, a cold heat exchanger (CHX), a tube, a hot heat exchanger numerical model has been introduced to study steady oscillatory heat and mass transfer in the tube section, numerical simulation, high resolution scheme 1 Introduction The pulse tube is a relatively new type

  20. Numerical modeling of extreme rogue waves generated by directional energy focusing

    E-Print Network [OSTI]

    Grilli, Stéphan T.

    Numerical modeling of extreme rogue waves generated by directional energy focusing Christophe that contributes to the generation of extreme waves, also known as rogue waves, in the ocean. To simulate and analyze this phenomenon, we generate extreme waves in a 3D numerical wave tank (NWT), by specifying

  1. Estimate of recharge from radiocarbon dating of groundwater and numerical flow and transport modeling

    E-Print Network [OSTI]

    Zhu, Chen

    of radioactive waste repositories. A case study of a regional aquifer in northeastern Arizona shows the rechargeEstimate of recharge from radiocarbon dating of groundwater and numerical flow and transport ages using a linked numerical 14 C transport and flow model while hydraulic conductivity values

  2. NUMERICAL MODELING OF COASTAL TSUNAMI DISSIPATION AND Stphan T. Grilli1

    E-Print Network [OSTI]

    Grilli, Stéphan T.

    1 NUMERICAL MODELING OF COASTAL TSUNAMI DISSIPATION AND IMPACT Stéphan T. Grilli1 , Jeffrey C Tehranirad2 Recent observations of the coastal impact of large tsunamis (e.g., Indian Ocean 2004; Tohoku 2011) and related numerical and theoretical works have made it increasingly clear that tsunami waves arrive

  3. NUMERICAL MODELING OF COASTAL TSUNAMI IMPACT DISSIPATION Stphan T. Grilli1

    E-Print Network [OSTI]

    Kirby, James T.

    1 NUMERICAL MODELING OF COASTAL TSUNAMI IMPACT DISSIPATION AND IMPACT Stéphan T. Grilli1 , Jeffrey Tehranirad2 Recent observations of the coastal impact of large tsunamis (e.g., Indian Ocean 2004; Tohoku 2011) and related numerical and theoretical works have made it increasingly clear that tsunami waves arrive

  4. Numerical model for steel catenary riser on seafloor support 

    E-Print Network [OSTI]

    You, Jung Hwan

    2007-04-25

    of development of a simplified seafloor support model. This model simulates the seafloor-pipe interaction as a flexible pipe supported on a bed of springs. Constants for the soil springs were derived from finite element studies performed in a separate, parallel...

  5. Numerically Estimating Internal Models of Dynamic Virtual Objects

    E-Print Network [OSTI]

    Sekuler, Robert

    human subjects to manipulate a computer-animated virtual object. This virtual object (vO) was a high, human cognition, human information processing, ideal performer, internal model, virtual object, virtual, specifically how humans acquire an internal model of a dynamic virtual object. Our methodology minimizes

  6. Decision-Making to Reduce Manufacturing Greenhouse Gas Emissions

    E-Print Network [OSTI]

    Reich-Weiser, Corinne

    2010-01-01

    Concentrator design to minimize LCOE,” Proceedings of thetool at SolFocus [122]. The LCOE cost model provides theinstallation variables, the LCOE metric is a step beyond the

  7. An Analysis of the Costs, Benefits, and Implications of Different Approaches to Capturing the Value of Renewable Energy Tax Incentives

    E-Print Network [OSTI]

    Bolinger, Mark

    2014-01-01

    levelized cost of energy (“LCOE”). Tax Equity Yield (after-power closer to achieving LCOE goals (and at no additionallevelized cost of energy (“LCOE”). 3. Model Descriptions and

  8. Seismic scattering attributes to estimate reservoir fracture density : a numerical modeling study

    E-Print Network [OSTI]

    Pearce, Frederick D. (Frederick Douglas), 1978-

    2003-01-01

    We use a 3-D finite difference numerical model to generate synthetic seismograms from a simple fractured reservoir containing evenly-spaced, discrete, vertical fracture zones. The fracture zones are represented using a ...

  9. NUMERICAL MODELING OF 3D ORGANIC SOLAR CELLS Presented to the

    E-Print Network [OSTI]

    Kassegne, Samuel Kinde

    .2 Physical Process in Polymer: Fullerene Build Heterojunction Solar CellsNUMERICAL MODELING OF 3D ORGANIC SOLAR CELLS _______________ A Thesis Presented to the Faculty Solar Cells by Anurag Kaushik Master of Science in Electrical Engineering San Diego State University

  10. Numerical modeling of elastic wave scattering by near-surface heterogeneities

    E-Print Network [OSTI]

    Al Muhaidib, Abdulaziz

    2013-01-01

    A perturbation method for elastic waves and numerical forward modeling are used to calculate the effects of seismic wave scattering from arbitrary shape shallow subsurface heterogeneities. Wave propagation is simulated ...

  11. Numerical Modelling of a Pulse Combustion Burner: Limiting Conditions of Stable

    E-Print Network [OSTI]

    Vuik, Kees

    Numerical Modelling of a Pulse Combustion Burner: Limiting Conditions of Stable Operation P.A. van in the burner system. Self-sustained pulse combustion and high-intensity sound waves result if the system

  12. Mass and charge flow in nanopores: numerical simulation via mesoscale models

    E-Print Network [OSTI]

    Cecconi, Fabio

    Mass and charge flow in nanopores: numerical simulation via mesoscale models Mauro Chinappi1 at nanoscale is here addressed via a recent developed mesoscale approach. In particular the flow

  13. 2D-Modelling of pellet injection in the poloidal plane: results of numerical tests

    E-Print Network [OSTI]

    Paris-Sud XI, Université de

    2D-Modelling of pellet injection in the poloidal plane: results of numerical tests P. Lalousis developed for computing the expansion of pellet-produced clouds in the poloidal plane. The expansion

  14. Observational and Numerical Modeling Studies of Turbulence on the Texas-Louisiana Continental Shelf 

    E-Print Network [OSTI]

    Zhang, Zheng

    2013-05-24

    ), dissipation rate of temperature variance (?), eddy diffusivity of temperature (?'t), and eddy diffusivity of density (?'p). Numerical models were also applied to test their capability of simulating these turbulence quantities. At site D, TKE, E, and ? were...

  15. A comparative analysis of numerical simulation and analytical modeling of horizontal well cyclic steam injection 

    E-Print Network [OSTI]

    Ravago Bastardo, Delmira Cristina

    2005-08-29

    The main objective of this research is to compare the performance of cyclic steam injection using horizontal wells based on the analytical model developed by Gunadi against that based on numerical simulation. For comparison, a common reservoir...

  16. P. Douillet S. Ouillon E. Cordier A numerical model for ne suspended sediment transport

    E-Print Network [OSTI]

    transport á Numerical model á Deposition á Erosion á New Caledonia Introduction In New Caledonia, where open-cast for cohesive particles (e.g., Nicholson and O'Connor 1986; Teisson 1991; Brenon and Le Hir 1999), around sand

  17. Numerical modeling and experimental validation of uniform microchamber filling in centrifugal microfluidics

    E-Print Network [OSTI]

    Banerjee, Debjyoti

    Numerical modeling and experimental validation of uniform microchamber filling in centrifugal analysis of microchamber filling in centrifugal microfluidics is presented. In the development of micro microchambers, such as those needed for nucleic acid amplification or detection. With centrifugal devices

  18. Seismic Scattering Attributes to Estimate Reservoir Fracture Density: A Numerical Modeling Study

    E-Print Network [OSTI]

    Pearce, Frederick Douglas

    We use a 3-D finite difference numerical model to generate synthetic seismograms from a simple fractured reservoir containing evenly-spaced, discrete, vertical fracture zones. The fracture zones are represented using a ...

  19. Numerical modeling of hydraulic fracture initiation and development

    E-Print Network [OSTI]

    2007-05-25

    Calculation scheme for modeling a hydraulic fracturing process: horizontal section of a ...... Jr., “Overview of current hydraulic fracturing design and treatment technology. .... A. A. Dobroskok, A. Ghassemi, and A. M. Linkov, “Extended structural ...

  20. NUMERICAL MODELING OF TURBULENT FLOW IN A COMBUSTION TUNNEL

    E-Print Network [OSTI]

    Ghoniem, A.F.

    2013-01-01

    1VJcDona·ld, H. (1979) Combustion r 1 iodeJ·ing in Two and1979) Practical Turbulent-Combustion Interaction Models forInternation on Combustors. Combustion The 17th Symposium

  1. Ductile fracture modeling : theory, experimental investigation and numerical verification

    E-Print Network [OSTI]

    Xue, Liang, 1973-

    2007-01-01

    The fracture initiation in ductile materials is governed by the damaging process along the plastic loading path. A new damage plasticity model for ductile fracture is proposed. Experimental results show that fracture ...

  2. A robotic crawler exploiting directional frictional interactions: experiments, numerics, and derivation of a reduced model

    E-Print Network [OSTI]

    Giovanni Noselli; Antonio DeSimone

    2014-08-26

    We present experimental and numerical results for a model crawler which is able to extract net positional changes from reciprocal shape changes, i.e. 'breathing-like' deformations, thanks to directional, frictional interactions with a textured solid substrate, mediated by flexible inclined feet. We also present a simple reduced model that captures the essential features of the kinematics and energetics of the gait, and compare its predictions with the results from experiments and from numerical simulations.

  3. Development and validation of a vertically two-dimensional mesoscale numerical model 

    E-Print Network [OSTI]

    Walters, Michael Kent

    1985-01-01

    DEVELOPMENT AND VALIDATION OF A VERTICALLY TWO-DIMENSIONAL MESOSCALE NUMERICAL MODEL A Thesis by MICHAEL KENT WALTERS Submitted to the Graduate College of Texas AsM University in partial fulfillment of the requirement for the degree of MASTER... OF SCIENCE August 1985 Major Subject: Meteorology DEVELOPMENT AND VALIDATION OF A VERTICALLY TWO-DIMENSIONAL MESOSCALE NUMERICAL MODEL A Thesis by MICHAEL KENT WALTERS Approved as to style and content by: Dusan Djuric (Chair of Committee) WP...

  4. Moist processes and the quasi-hydrostatic approximation in a mesoscale numerical model 

    E-Print Network [OSTI]

    Kennedy, Charles Joseph

    1987-01-01

    MOIST PROCESSES AND THE QUASI-HYDROSTATIC APPROXIMATION IN A MESOSCALE NUMERICAL MODEL A Thesis by CHARLES JOSEPH KENNEDY Submitted to the Graduate College of Texas A&M University in partial fulfillment of the requirement for the degree... of MASTER OF SCIENCE December 1987 Major Subject: Meteorology MOIST PROCESSES AND THE QUASI-HYDROSTATIC APPROXIMATION IN A MESOSCALE NUMERICAL MODEL A Thesis by CHARLES JOSEPH KENNEDY Approved as to style and content by: Dusan Djuric (Chair...

  5. Numerical modeling of the wind flow over a transverse dune

    E-Print Network [OSTI]

    Ascânio D. Araújo; Eric J. R. Parteli; Thorsten Poeschel; José S. Andrade Jr.; Hans J. Herrmann

    2013-09-30

    Transverse dunes, which form under unidirectional winds and have fixed profile in the direction perpendicular to the wind, occur on all celestial objects of our solar system where dunes have been detected. Here we perform a numerical study of the average turbulent wind flow over a transverse dune by means of computational fluid dynamics simulations. We find that the length of the zone of recirculating flow at the dune lee --- the {\\em{separation bubble}} --- displays a surprisingly strong dependence on the wind shear velocity, $u_{\\ast}$: it is nearly independent of $u_{\\ast}$ for shear velocities within the range between $0.2\\,$m$$s and $0.8\\,$m$$s but increases linearly with $u_{\\ast}$ for larger shear velocities. Our calculations show that transport in the direction opposite to dune migration within the separation bubble can be sustained if $u_{\\ast}$ is larger than approximately $0.39\\,$m$$s, whereas a larger value of $u_{\\ast}$ (about $0.49\\,$m$$s) is required to initiate this reverse transport.

  6. Numerical modeling of the wind flow over a transverse dune

    E-Print Network [OSTI]

    Araújo, Ascânio D; Poeschel, Thorsten; Andrade, José S; Herrmann, Hans J

    2013-01-01

    Transverse dunes, which form under unidirectional winds and have fixed profile in the direction perpendicular to the wind, occur on all celestial objects of our solar system where dunes have been detected. Here we perform a numerical study of the average turbulent wind flow over a transverse dune by means of computational fluid dynamics simulations. We find that the length of the zone of recirculating flow at the dune lee --- the {\\em{separation bubble}} --- displays a surprisingly strong dependence on the wind shear velocity, $u_{\\ast}$: it is nearly independent of $u_{\\ast}$ for shear velocities within the range between $0.2\\,$m$$s and $0.8\\,$m$$s but increases linearly with $u_{\\ast}$ for larger shear velocities. Our calculations show that transport in the direction opposite to dune migration within the separation bubble can be sustained if $u_{\\ast}$ is larger than approximately $0.39\\,$m$$s, whereas a larger value of $u_{\\ast}$ (about $0.49\\,$m$$s) is required to initiate this reverse transport.

  7. Numerical models of pressure pulse generation by imploding metal liners

    SciTech Connect (OSTI)

    Humphries, S. Jr. [Acceleration Associates, Albuquerque, NM (United States)] [Acceleration Associates, Albuquerque, NM (United States); [Univ. of New Mexico, Albuquerque, NM (United States). Dept. of Electrical Engineering; Ekdahl, C.A. [Los Alamos National Lab., NM (United States)] [Los Alamos National Lab., NM (United States)

    1996-12-01

    The authors describe numerical calculations of pressure pulse generation using imploding liners. Liners are metal cylinders that are magnetically compressed by an intense axial current flow from a high-power pulse generator. The simulations cover the acceleration of the liner, collision with an internal diagnostic target, followed by compression and shock wave heating of the target. With the projected current waveform of the Atlas capacitor bank (in development at Los Alamos National Laboratory), initial results suggest that it may be possible to achieve pressures exceeding 3,000 Gpa (30 Mbar) in a 4 mm diameter sample over an interval of 100--200 ns. The simulations were carried out with Crunch, a new one-dimensional hydrodynamics package for advanced personal computers. The program uses finite-element techniques to solve the coupled problems of hydrodynamics and magnetic diffusion. Crunch fully supports loading and interpolating Sesame equation-of-state tables. The program exhibits excellent stability, even for collisions between material shells and shock convergence on axis. In contrast to previous work, the present studies follow the full process through solid target collision and compression. The work supports the High-Energy Density Physics Program of Los Alamos National Laboratory (LANL), a component of the US Department of Energy Stockpile Stewardship Program. The purpose of this program is maintenance of the nuclear weapons stockpile through improved computational ability and above-ground experiments. Imploding liners driven by conventional capacitor banks constitute a portion of the program to study matter at high pressure.

  8. 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.

  9. IMPROVED NUMERICAL METHODS FOR MODELING RIVER-AQUIFER INTERACTION.

    SciTech Connect (OSTI)

    Tidwell, Vincent C.; Sue Tillery; Phillip King

    2008-09-01

    A new option for Local Time-Stepping (LTS) was developed to use in conjunction with the multiple-refined-area grid capability of the U.S. Geological Survey's (USGS) groundwater modeling program, MODFLOW-LGR (MF-LGR). The LTS option allows each local, refined-area grid to simulate multiple stress periods within each stress period of a coarser, regional grid. This option is an alternative to the current method of MF-LGR whereby the refined grids are required to have the same stress period and time-step structure as the coarse grid. The MF-LGR method for simulating multiple-refined grids essentially defines each grid as a complete model, then for each coarse grid time-step, iteratively runs each model until the head and flux changes at the interfacing boundaries of the models are less than some specified tolerances. Use of the LTS option is illustrated in two hypothetical test cases consisting of a dual well pumping system and a hydraulically connected stream-aquifer system, and one field application. Each of the hypothetical test cases was simulated with multiple scenarios including an LTS scenario, which combined a monthly stress period for a coarse grid model with a daily stress period for a refined grid model. The other scenarios simulated various combinations of grid spacing and temporal refinement using standard MODFLOW model constructs. The field application simulated an irrigated corridor along the Lower Rio Grande River in New Mexico, with refinement of a small agricultural area in the irrigated corridor.The results from the LTS scenarios for the hypothetical test cases closely replicated the results from the true scenarios in the refined areas of interest. The head errors of the LTS scenarios were much smaller than from the other scenarios in relation to the true solution, and the run times for the LTS models were three to six times faster than the true models for the dual well and stream-aquifer test cases, respectively. The results of the field application show that better estimates of daily stream leakage can be made with the LTS simulation, thereby improving the efficiency of daily operations for an agricultural irrigation system. ACKNOWLEDGEMENTSThe authors appreciatively acknowledge support for Sue Tillery provided by Sandia National Laboratories' through a Campus Executive Laboratory Directed Research and Development (LDRD) research project.Funding for this study was provided by Directed Research and Development (LDRD) research project.

  10. Numerical Solution of a Free-Boundary Problem for Percussive Deep Drilling Modeling by BEM

    E-Print Network [OSTI]

    Mikhailov, Sergey

    Numerical Solution of a Free-Boundary Problem for Percussive Deep Drilling Modeling by BEM S to a stationary-periodic quasi-static model of rock percus- sive deep drilling is presented. The rock is modeled. An auxiliary problem of stationary inden- tation of a rigid drill bit is considered first, where it is assumed

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

    SciTech Connect (OSTI)

    Xing, Lu [Oklahoma State University; Cullin, James [Oklahoma State University; Spitler, Jeffery [Oklahoma State University; Im, Piljae [ORNL; Fisher, Daniel [Oklahoma State University

    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.

  12. Acoustics 2000 1 The Two Dimensional Numerical Modeling

    E-Print Network [OSTI]

    -difference methods have often been used. This report terms the wave equations suited to waves in fluids, acoustic, Sven Treitel, and Alford, 1976) but the acoustic wave equations have also been used for geophysical. Key features of the model at present are: (i) The use of acoustic wave equation (ii) Time domain

  13. 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.

  14. Numerical modelling and analysis of a room temperature magnetic

    E-Print Network [OSTI]

    and analysis of a room temperature magnetic refrigeration system Department: Fuel Cells and Solid State are separated by channels of a heat transfer fluid. The time-dependent model solves the momentum and continuity equations of the flow of the heat transfer fluid and the coupled energy equations of the heat transfer

  15. A Two Pressure Numerical Model of Two Fluid Mixing \\Lambda

    E-Print Network [OSTI]

    New York at Stoney Brook, State University of

    fluid. The model is designed for the study of acceleration driven mixing layers in a chunk mix regime to the description of natural phenomena, as in meteorological temperature in­ version driven updrafts quantities directly, to close the system of equations. Closure is a physics level identity, independent

  16. NUMERICAL MODELING OF FLUID FLOW AND TIME-LAPSE ...

    E-Print Network [OSTI]

    gabriela

    CO2 injection operation at the Sleipner gas field in the North Sea, operated by Statoil ... The simultaneous flow of brine and CO2 is modeled with the Black-Oil formulation for ..... As water saturation is reduced, and the larger pores drained first, ...

  17. Numerical Modeling of Seafloor Interation with Steel Catenary Riser 

    E-Print Network [OSTI]

    You, Jung Hwan

    2012-10-19

    , and moment. Primary outputs from this model include the deflected shape of the riser pipe and bending moments along riser length. The code also provides the location of maximum trench depth and the position where the maximum bending moment occurs and any...

  18. THEACCURACY OF NUMERICAL MODELS FOR CONTINUUM STANLY STEINBERG

    E-Print Network [OSTI]

    Steinberg, Stanly

    Water Rain Evaporation Figure 1.1. A Lake These notes are directed towards students who have some modest of water in the lake. Water can enter and leave the lake by various means: rain; evaporation; rivers and how much error is there is solving the mathematical model. #12; 2 Lake Ground Ground Water Ground

  19. 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.

  20. 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.

  1. Mathematical formulation and numerical modeling of wax deposition in pipelines from enthalpyporosity approach

    E-Print Network [OSTI]

    Firoozabadi, Abbas

    Mathematical formulation and numerical modeling of wax deposition in pipelines from enthalpy in the layer composed of the two-phase wax and oil (that is the gel layer), and (3) use Fick's law to describe to temperature gradient. In this work, a rigorous math- ematical model for the prediction of wax deposition

  2. Intercomparison of Single-Column Numerical Models for the Prediction of Radiation Fog

    E-Print Network [OSTI]

    Intercomparison of Single-Column Numerical Models for the Prediction of Radiation Fog THIERRY-term forecasting of fog is a difficult issue that can have a large societal impact. Radiation fog appears layers of the atmosphere. Current NWP models poorly forecast the life cycle of fog, and improved NWP

  3. Modeling and Design of RF MEMS Structures Using Computationally Efficient Numerical Techniques

    E-Print Network [OSTI]

    Tentzeris, Manos

    Modeling and Design of RF MEMS Structures Using Computationally Efficient Numerical Techniques N. A Abstract The modeling of MEMS structures using MRTD is presented. Many complex RF structures have been communication systems efficiently and accurately. Specifically, micromachined structures such as MEMS

  4. NUMERICAL COUPLING OF ELECTRIC CIRCUIT EQUATIONS AND ENERGY-TRANSPORT MODELS FOR SEMICONDUCTORS

    E-Print Network [OSTI]

    Jüngel, Ansgar

    NUMERICAL COUPLING OF ELECTRIC CIRCUIT EQUATIONS AND ENERGY-TRANSPORT MODELS FOR SEMICONDUCTORS effects is proposed. The charged particle flow in the semiconductor devices is described by the energy-transport equations for the electrons and the drift-diffusion equations for the holes. The electric circuit is modeled

  5. Effects of winds and Caribbean eddies on the frequency of Loop Current eddy shedding: A numerical model study

    E-Print Network [OSTI]

    Effects of winds and Caribbean eddies on the frequency of Loop Current eddy shedding: A numerical of Mexico, Caribbean Current, Loop Current, eddy shedding, winds and eddies, numerical ocean model Citation of Loop Current eddy shedding: A numerical model study, J. Geophys. Res., 108(C10), 3324, doi:10

  6. Characterization of Texas lignite and numerical modeling of its in-situ gasification 

    E-Print Network [OSTI]

    Wang, Yih-Jy

    1983-01-01

    Modeling Site selection for in-situ gasification projects normally involves application of site screen1ng criteria. Some of these cr1teria were discussed by Russell et al. (1983). Numerical simulation may play an important role in s1te selection...CHARACTERIZATION OF TEXAS LIGNITE AND NUMERICAL MODELING OF ITS IN-SITU GASIFICATION A Thesis by YIH-JY WANG Submitted to the Graduate College of Texas A&M University in partial fulfillment of the requirements for the degree of MASTER...

  7. 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

  8. The NINJA-2 project: Detecting and characterizing gravitational waveforms modelled using numerical binary black hole simulations

    E-Print Network [OSTI]

    :,; Abbott, B P; Abbott, R; Abbott, T; Abernathy, M R; Accadia, T; Acernese, F; Ackley, K; Adams, C; Adams, T; Addesso, P; Adhikari, R X; Affeldt, C; Agathos, M; Aggarwal, N; Aguiar, O D; Ain, A; Ajith, P; Alemic, A; Allen, B; Allocca, A; Amariutei, D; Andersen, M; Anderson, R; Anderson, S B; Anderson, W G; Arai, K; Araya, M C; Arceneaux, C; Areeda, J; Aston, S M; Astone, P; Aufmuth, P; Aulbert, C; Austin, L; Aylott, B E; Babak, S; Baker, P T; Ballardin, G; Ballmer, S W; Barayoga, J C; Barbet, M; Barish, B C; Barker, D; Barone, F; Barr, B; Barsotti, L; Barsuglia, M; Barton, M A; Bartos, I; Bassiri, R; Basti, A; Batch, J C; Bauchrowitz, J; Bauer, Th S; Behnke, B; Bejger, M; Beker, M G; Belczynski, C; Bell, A S; Bell, C; Bergmann, G; Bersanetti, D; Bertolini, A; Betzwieser, J; Beyersdorf, P T; Bilenko, I A; Billingsley, G; Birch, J; Biscans, S; Bitossi, M; Bizouard, M A; Black, E; Blackburn, J K; Blackburn, L; Blair, D; Bloemen, S; Blom, M; Bock, O; Bodiya, T P; Boer, M; Bogaert, G; Bogan, C; Bond, C; Bondu, F; Bonelli, L; Bonnand, R; Bork, R; Born, M; Boschi, V; Bose, Sukanta; Bosi, L; Bradaschia, C; Brady, P R; Braginsky, V B; Branchesi, M; Brau, J E; Briant, T; Bridges, D O; Brillet, A; Brinkmann, M; Brisson, V; Brooks, A F; Brown, D A; Brown, D D; Brückner, F; Buchman, S; Bulik, T; Bulten, H J; Buonanno, A; Burman, R; Buskulic, D; Buy, C; Cadonati, L; Cagnoli, G; Bustillo, J Calderón; Calloni, E; Camp, J B; Campsie, P; Cannon, K C; Canuel, B; Cao, J; Capano, C D; Carbognani, F; Carbone, L; Caride, S; Castiglia, A; Caudill, S; Cavaglià, M; Cavalier, F; Cavalieri, R; Celerier, C; Cella, G; Cepeda, C; Cesarini, E; Chakraborty, R; Chalermsongsak, T; Chamberlin, S J; Chao, S; Charlton, P; Chassande-Mottin, E; Chen, X; Chen, Y; Chincarini, A; Chiummo, A; Cho, H S; Chow, J; Christensen, N; Chu, Q; Chua, S S Y; Chung, S; Ciani, G; Clara, F; Clark, J A; Cleva, F; Coccia, E; Cohadon, P -F; Colla, A; Collette, C; Colombini, M; Cominsky, L; Constancio, M; Conte, A; Cook, D; Corbitt, T R; Cordier, M; Cornish, N; Corpuz, A; Corsi, A; Costa, C A; Coughlin, M W; Coughlin, S; Coulon, J -P; Countryman, S; Couvares, P; Coward, D M; Cowart, M; Coyne, D C; Coyne, R; Craig, K; Creighton, J D E; Crowder, S G; Cumming, A; Cunningham, L; Cuoco, E; Dahl, K; Canton, T Dal; Damjanic, M; Danilishin, S L; D'Antonio, S; Danzmann, K; Dattilo, V; Daveloza, H; Davier, M; Davies, G S; Daw, E J; Day, R; Dayanga, T; Debreczeni, G; Degallaix, J; Deléglise, S; Del Pozzo, W; Denker, T; Dent, T; Dereli, H; Dergachev, V; De Rosa, R; DeRosa, R T; DeSalvo, R; Dhurandhar, S; Díaz, M; Di Fiore, L; Di Lieto, A; Di Palma, I; Di Virgilio, A; Donath, A; Donovan, F; Dooley, K L; Doravari, S; Dossa, S; Douglas, R; Downes, T P; Drago, M; Drever, R W P; Driggers, J C; Du, Z; Dwyer, S; Eberle, T; Edo, T; Edwards, M; Effler, A; Eggenstein, H; Ehrens, P; Eichholz, J; Eikenberry, S S; Endr?czi, G; Essick, R; Etzel, T; Evans, M; Evans, T; Factourovich, M; Fafone, V; Fairhurst, S; Fang, Q; Farinon, S; Farr, B; Farr, W M; Favata, M; Fehrmann, H; Fejer, M M; Feldbaum, D; Feroz, F; Ferrante, I; Ferrini, F; Fidecaro, F; Finn, L S; Fiori, I; Fisher, R P; Flaminio, R; Fournier, J -D; Franco, S; Frasca, S; Frasconi, F; Frede, M; Frei, Z; Freise, A; Frey, R; Fricke, T T; Fritschel, P; Frolov, V V; Fulda, P; Fyffe, M; Gair, J; Gammaitoni, L; Gaonkar, S; Garufi, F; Gehrels, N; Gemme, G; Genin, E; Gennai, A; Ghosh, S; Giaime, J A; Giardina, K D; Giazotto, A; Gill, C; Gleason, J; Goetz, E; Goetz, R; Gondan, L; González, G; Gordon, N; Gorodetsky, M L; Gossan, S; Goßler, S; Gouaty, R; Gräf, C; Graff, P B; Granata, M; Grant, A; Gras, S; Gray, C; Greenhalgh, R J S; Gretarsson, A M; Groot, P; Grote, H; Grover, K; Grunewald, S; Guidi, G M; Guido, C; Gushwa, K; Gustafson, E K; Gustafson, R; Hammer, D; Hammond, G; Hanke, M; Hanks, J; Hanna, C; Hanson, J; Harms, J; Harry, G M; Harry, I W; Harstad, E D; Hart, M; Hartman, M T; Haster, C -J; Haughian, K; Heidmann, A; Heintze, M; Heitmann, H; Hello, P; Hemming, G; Hendry, M; Heng, I S; Heptonstall, A W; Heurs, M; Hewitson, M; Hild, S; Hoak, D; Hodge, K A; Holt, K; Hooper, S; Hopkins, P; Hosken, D J; Hough, J; Howell, E J; Hu, Y; Hughey, B; Husa, S; Huttner, S H; Huynh, M; Huynh-Dinh, T; Ingram, D R; Inta, R; Isogai, T; Ivanov, A; Iyer, B R; Izumi, K; Jacobson, M; James, E; Jang, H; Jaranowski, P; Ji, Y; Jiménez-Forteza, F; Johnson, W W; Jones, D I; Jones, R; Jonker, R J G; Ju, L; K, Haris; Kalmus, P; Kalogera, V; Kandhasamy, S; Kang, G; Kanner, J B; Karlen, J; Kasprzack, M; Katsavounidis, E; Katzman, W; Kaufer, H; Kawabe, K; Kawazoe, F; Kéfélian, F; Keiser, G M; Keitel, D; Kelley, D B; Kells, W; Khalaidovski, A; Khalili, F Y; Khazanov, E A; Kim, C; Kim, K; Kim, N; Kim, N G; Kim, Y -M; King, E J; King, P J; Kinzel, D L; Kissel, J S; Klimenko, S; Kline, J; Koehlenbeck, S; Kokeyama, K; Kondrashov, V; Koranda, S; Korth, W Z; Kowalska, I

    2014-01-01

    The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational-wave astrophysics communities. The purpose of NINJA is to study the ability to detect gravitational waves emitted from merging binary black holes and recover their parameters with next-generation gravitational-wave observatories. We report here on the results of the second NINJA project, NINJA-2, which employs 60 complete binary black hole hybrid waveforms consisting of a numerical portion modelling the late inspiral, merger, and ringdown stitched to a post-Newtonian portion modelling the early inspiral. In a "blind injection challenge" similar to that conducted in recent LIGO and Virgo science runs, we added 7 hybrid waveforms to two months of data recolored to predictions of Advanced LIGO and Advanced Virgo sensitivity curves during their first observing runs. The resulting data was analyzed by gravitational-wave detection algorithms and 6 of the waveforms were recovered w...

  9. Numerical model for vibration damping resulting from the first-order phase transformations

    E-Print Network [OSTI]

    Melnik, Roderick

    : Martensite transformation; Thermo-mechanical coupling; Vibration damping; Ginzburg­Landau theory 1Numerical model for vibration damping resulting from the first-order phase transformations L Landau­Ginzburg theory that couples nonlinear mechanical and thermal fields. The free energy function

  10. Nonlinear inverse problem for a model of ion-exchange filter: numerical recovery of parameters

    E-Print Network [OSTI]

    methods of ion exchange is based on passing of source water or partially treated water through a filter-exchange filters results in highly mineralized, acidic and alkaline waste water [9]. With the continuous1 Nonlinear inverse problem for a model of ion-exchange filter: numerical recovery of parameters

  11. Numerical modeling of fluid flow and time-lapse seismograms applied to

    E-Print Network [OSTI]

    Santos, Juan

    ; and CO2 and CO2 are the CO2 mole fraction and the CO2 mass fraction in the brine phase. This conversionNumerical modeling of fluid flow and time-lapse seismograms applied to CO2 storage and monitoring G and time-lapse seismograms applied to CO2 storage and monitoring ­ p. #12;Introduction · Fossil

  12. Numerical analysis of electric field formulations of the eddy current model

    E-Print Network [OSTI]

    Rodríguez, Rodolfo

    Numerical analysis of electric field formulations of the eddy current model Alfredo Berm´udez1 methods for the numeri- cal solution of the eddy current problem in a bounded conducting domain crossed): 78M10, 65N30 Key words Low-frequency harmonic Maxwell equations, eddy currents, finite elements

  13. Mathematical and numerical analysis of a transient non-linear axisymmetric eddy current model

    E-Print Network [OSTI]

    Rodríguez, Rodolfo

    Mathematical and numerical analysis of a transient non-linear axisymmetric eddy current model the theoretically predicted behavior of the method, are reported. Keywords transient eddy current · axisymmetric is the accurate computation of power losses in the ferromagnetic components of the core due to hysteresis and eddy-current

  14. NUMERICAL SOLUTION OF A TRANSIENT NON-LINEAR AXISYMMETRIC EDDY CURRENT MODEL WITH NON-LOCAL

    E-Print Network [OSTI]

    Rodríguez, Rodolfo

    NUMERICAL SOLUTION OF A TRANSIENT NON-LINEAR AXISYMMETRIC EDDY CURRENT MODEL WITH NON@ing-mat.udec.cl This paper deals with an axisymmetric transient eddy current problem in conductive nonlinear magnetic media of the proposed scheme. Keywords: transient eddy current problem; electromagnetic losses; nonlinear magnetic

  15. A simple numerical model of the apparent loss of eddy current conductivity due to surface roughness

    E-Print Network [OSTI]

    Nagy, Peter B.

    A simple numerical model of the apparent loss of eddy current conductivity due to surface roughness of eddy current conductivity has been suggested as a possible means to allow the nondestructive evaluation, the path of the eddy current must follow a more tortuous route in the material, which produces a reduction

  16. 792 / JOURNAL OF ENVIRONMENTAL ENGINEERING / SEPTEMBER 1999 NUMERICAL MODEL OF SEDIMENTATION/THICKENING

    E-Print Network [OSTI]

    Wells, Scott A.

    792 / JOURNAL OF ENVIRONMENTAL ENGINEERING / SEPTEMBER 1999 NUMERICAL MODEL OF SEDIMENTATION with the ASCE Manager of Journals. The manuscript for this paper was submitted for review and possible publication on July 20, 1998. This paper is part of the Journal of Environmental Engineering, Vol. 125, No. 9

  17. Numerical and analytical modelling of entropy noise in a supersonic nozzle with a shock

    E-Print Network [OSTI]

    Nicoud, Franck

    Numerical and analytical modelling of entropy noise in a supersonic nozzle with a shock M. Leyko a these sources, the noise coming from the turbulent flame within the combustor is already identified in a supersonic nozzle with a shock, Journal of Sound and Vibration (2011), doi:10.1016/j.jsv.2011.01.025 #12;Two

  18. DROP-SCALE NUMERICAL MODELING OF CHEMICAL PARTITIOING DURING CLOUD HYDROMETEOR FREEZING

    E-Print Network [OSTI]

    Stuart, Amy L.

    DROP-SCALE NUMERICAL MODELING OF CHEMICAL PARTITIOING DURING CLOUD HYDROMETEOR FREEZING A.L. Stuart freezing provide greatly varying estimates of the retention efficiency of volatile solutes (e.g., Lamb understanding of the dependence of partitioning on chemical properties and freezing conditions (Stuart

  19. International Series of Numerical Mathematics, Vol. 160, 521540 Freezing of Living Cells: Mathematical Models

    E-Print Network [OSTI]

    Turova, Varvara

    International Series of Numerical Mathematics, Vol. 160, 521­540 Freezing of Living Cells, stresses arising due to non-simultaneous freezing of water in- side and outside of cells are modeled and outside of living cells during freezing is derived by applying an appropriate averaging technique

  20. COMPUTATIONAL CHALLENGES IN THE NUMERICAL TREATMENT OF LARGE AIR POLLUTION MODELS

    E-Print Network [OSTI]

    Ostromsky, Tzvetan

    COMPUTATIONAL CHALLENGES IN THE NUMERICAL TREATMENT OF LARGE AIR POLLUTION MODELS I. DIMOV , K. GEORGIEVy, TZ. OSTROMSKY , R. J. VAN DER PASz, AND Z. ZLATEVx Abstract. The air pollution, and especially the reduction of the air pollution to some acceptable levels, is an important environmental problem, which

  1. Ensemble Kalman Filter Data Assimilation in a 1D Numerical Model Used for Fog Forecasting

    E-Print Network [OSTI]

    Ensemble Kalman Filter Data Assimilation in a 1D Numerical Model Used for Fog Forecasting SAMUEL RE, a need exists for accurate and updated fog and low-cloud forecasts. Couche Brouillard Eau Liquide (COBEL for the very short-term forecast of fog and low clouds. This forecast system assimilates local observations

  2. A faster numerical scheme for a coupled system modelling soil erosion and sediment transport

    E-Print Network [OSTI]

    d'Orléans, Université

    A faster numerical scheme for a coupled system modelling soil erosion and sediment transport M flow and the bed sediment, are classically described by a well-established system coupling the shallow states and the positivity of both water depth and sediment concentration. Recently, finite volume schemes

  3. Mathematical, physical and numerical principles essential for models of turbulent mixing

    SciTech Connect (OSTI)

    Sharp, David Howland [Los Alamos National Laboratory; Lim, Hyunkyung [STONY BROOK UNIV; Yu, Yan [STONY BROOK UNIV; Glimm, James G [STONY BROOK UNIV

    2009-01-01

    We propose mathematical, physical and numerical principles which are important for the modeling of turbulent mixing, especially the classical and well studied Rayleigh-Taylor and Richtmyer-Meshkov instabilities which involve acceleration driven mixing of a fluid discontinuity layer, by a steady accerleration or an impulsive force.

  4. Modeling and Numerical Simulation of Bioheat Transfer and Biomechanics in Soft Tissue

    E-Print Network [OSTI]

    Zhang, Jun

    Modeling and Numerical Simulation of Bioheat Transfer and Biomechanics in Soft Tissue #3; Wensheng techniques are eÆcient. Key words: Bioheat transfer, biomechanics, discretization, iterative solver. 1 do in engineering area by solving constitutive equations. One of the major diÆculties in biomechanics

  5. NUMERICAL MODELLING OF THE DEEP IMPACT MISSION EXPERIMENT K. Wnnemann1 , and H. J. Melosh1

    E-Print Network [OSTI]

    Collins, Gareth

    NUMERICAL MODELLING OF THE DEEP IMPACT MISSION EXPERIMENT K. Wünnemann1 , G. S. Collins2 , and H. J@lpl.arizona.edu, 2 Department of Earth Science and Engineering, Imperial College London, London SW7 2AZ, UK, g.collins@imperial.ac.uk . Introduction: NASA's Deep Impact Mission (launched January 2005) will provide, for the first time ever

  6. Numerical Simulation of the December 26, 2004 Indian Ocean Tsunami using a Boussinesq model

    E-Print Network [OSTI]

    Kirby, James T.

    tsunami education or tsunami warning system in the re- gion exaccerbated the number of fatalities, even. Scientists had been warning of the growing exposure of coastal residents to tsunami hazards for yearsNumerical Simulation of the December 26, 2004 Indian Ocean Tsunami using a Boussinesq model Philip

  7. September 25, 2006 Numerical modeling of the effect of carbon dioxide

    E-Print Network [OSTI]

    Paris-Sud XI, Université de

    September 25, 2006 Numerical modeling of the effect of carbon dioxide sequestration on the rate souterrain de dioxyde de carbone sur la déformation des calcaires par dissolution sous contrainte: résultats@obs.ujf- grenoble.fr, marielle.collombet@ujf-grenoble.fr, yleguen@lgit.obs.ujf-grenoble.fr. #12;Abstract When carbon

  8. LABORATORY-NUMERICAL MODEL COMPARISONS OF CANYON FLOWS: A PARAMETER STUDY.

    E-Print Network [OSTI]

    LABORATORY-NUMERICAL MODEL COMPARISONS OF CANYON FLOWS: A PARAMETER STUDY. DON L. BOYER Cite Scientifique ­ 59655 - Villeneuve d'Ascq cedex (France) ABSTRACT An integrated set of laboratory surrounding a coastal canyon and, further, to explore the degree to which laboratory experiments can provide

  9. CONCEPTUAL AND NUMERICAL MODELS OF OXYGEN DIFFUSION, SULPHIDE OXIDATION AND ACID MINE DRAINAGE WITHIN

    E-Print Network [OSTI]

    Aubertin, Michel

    CONCEPTUAL AND NUMERICAL MODELS OF OXYGEN DIFFUSION, SULPHIDE OXIDATION AND ACID MINE DRAINAGE The generation and transport of acid mine drainage (AMD) through discretely fractured porous media is simulated. RÉSUMÉ La génération et le transport des produits issus du drainage minier acide (DMA) dans un milieu

  10. January 2, 2008 Numerical modeling of the effect of carbon dioxide

    E-Print Network [OSTI]

    Boyer, Edmond

    January 2, 2008 Numerical modeling of the effect of carbon dioxide sequestration on the rate souterrain de dioxyde de carbone sur la déformation des calcaires par dissolution sous contrainte: résultats;Abstract When carbon dioxide (CO2) is injected into an aquifer or a depleted geological reservoir, its

  11. Numerical modelling of the blowing phase in the production of glass containers

    E-Print Network [OSTI]

    Eindhoven, Technische Universiteit

    Numerical modelling of the blowing phase in the production of glass containers W. Dijkstraa , R to simulate the blowing of glass containers that are not rotationally symmetric. The contact between glass illustrates the results. Keywords: Boundary element method, blowing phase, glass, Stokes equations 1

  12. MODELING AND ADAPTIVE NUMERICAL TECHNIQUES FOR OXI-DATION OF CERAMIC COMPOSITES

    E-Print Network [OSTI]

    Adjerid, Slimane

    . INTRODUCTION Oxidation shortens the life of ceramic matrix composites by, e.g., chang- ing the elasticMODELING AND ADAPTIVE NUMERICAL TECHNIQUES FOR OXI- DATION OF CERAMIC COMPOSITES S. Adjerid, M. Ai reaction 1-3 . Composite materials are protected by coatings; however, cracks that form as a result

  13. NUMERICAL MODELLING OF AUTOGENOUS HEALING AND RECOVERY OF MECHANICAL PROPERTIES IN ULTRA-HIGH

    E-Print Network [OSTI]

    Boyer, Edmond

    in the cementitious matrix can react with carbon dioxide dissolved in the water filling the crack. Autogenous healingNUMERICAL MODELLING OF AUTOGENOUS HEALING AND RECOVERY OF MECHANICAL PROPERTIES IN ULTRA into the crack and leads to a partial recovery of mechanical properties (Young's modulus, tensile strength

  14. Numerical Passage from Radiative Heat Transfer to Nonlinear Diffusion Models \\Lambda

    E-Print Network [OSTI]

    Schmeiser, Christian

    Numerical Passage from Radiative Heat Transfer to Nonlinear Diffusion Models \\Lambda A. Klar y C. Schmeiser z Abstract Radiative heat transfer equations including heat conduction are consid­ ered situations are presented. Keywords. radiative heat transfer, asymptotic analysis, nonlinear diffusion limit

  15. LABORATORY OBSERVATIONS AND NUMERICAL MODELING OF THE EFFECTS OF AN ARRAY OF WAVE ENERGY CONVERTERS

    E-Print Network [OSTI]

    Haller, Merrick

    1 LABORATORY OBSERVATIONS AND NUMERICAL MODELING OF THE EFFECTS OF AN ARRAY OF WAVE ENERGY of wave energy converters (WECs) on water waves through the analysis of extensive laboratory experiments absorption is a reasonable predictor of the effect of WECs on the far field. Keywords: wave- energy; spectral

  16. Numerical Modeling of Hydro-acoustic Waves In Weakly Compressible Fluid Ali Abdolali1,2

    E-Print Network [OSTI]

    Kirby, James T.

    Numerical Modeling of Hydro-acoustic Waves In Weakly Compressible Fluid Ali Abdolali1,2 , James T of Civil Engineering, University of Roma Tre Low-frequency hydro-acoustic waves are precursors of tsunamis. Detection of hydro-acoustic waves generated due to the water column compression triggered by sudden seabed

  17. Numerical simulation of breaking waves by a multi-scale turbulence model

    E-Print Network [OSTI]

    Zhao, Qun

    and diffusion are of the same order at the trough level. Above the trough level, turbulent convection dominates-dimensional MAC type finite difference method. The third-order upwind scheme proposed by Kawamura and KawaharaNumerical simulation of breaking waves by a multi-scale turbulence model Qun Zhaoa,*, Steve

  18. On the application of robust numerical methods to a complete-flow wave-current model

    E-Print Network [OSTI]

    Madden, Niall

    -current interaction in the presence of weakly turbulent flow leading to an Orr-Sommerfeld type problem, and a system of two singularly perturbed reaction-diffusion equations from a k- turbulence model. The numerical of modifying the scheme to stabilize it, a mesh tailored to the specific problem is used. In this study we

  19. On the application of robust numerical methods to a complete-flow wave-current model

    E-Print Network [OSTI]

    Madden, Niall

    -current interaction in the presence of weakly turbulent flow leading to an Orr-Sommerfeld type problem and a system of two singularly perturbed reaction-diffusion equations from a k- turbulence model. The numerical layers. Instead of modifying the scheme to stabilize it, a mesh tailored to the specific problem is used

  20. On the seismic response of deep-seated rock slope instabilities --Insights from numerical modeling

    E-Print Network [OSTI]

    On the seismic response of deep-seated rock slope instabilities -- Insights from numerical modeling 2015 Accepted 3 April 2015 Available online 15 April 2015 Keywords: Seismic landslide hazard Spectral component of seismic hazard in mountainous regions. While many seismic slope stability analysis methods

  1. Numerical modeling of heat transfer and fluid flow in rotor-stator cavities with throughflow

    E-Print Network [OSTI]

    Boyer, Edmond

    Numerical modeling of heat transfer and fluid flow in rotor-stator cavities with throughflow S in a rotor-stator cavity subjected to a superimposed throughflow with heat transfer. Nu- merical predictions field from the heat transfer process. The turbulent flux is approximated by a gradient hypothesis

  2. Numerical modeling of heat transfer and fluid flow in rotor-stator cavities with throughflow

    E-Print Network [OSTI]

    Boyer, Edmond

    Numerical modeling of heat transfer and fluid flow in rotor-stator cavities with throughflow S the dynamical effects from the heat transfer process. The fluid flow in an enclosed disk system with axial with heat transfer along the stator, which corresponds to the experiment of Djaoui et al. [2]. Our results

  3. Precipitation sensitivity to autoconversion rate in a Numerical Weather Prediction model

    E-Print Network [OSTI]

    Marsham, John

    1 Precipitation sensitivity to autoconversion rate in a Numerical Weather Prediction model Céline;2 Summary Aerosols are known to significantly affect cloud and precipitation patterns and intensity. The impact of changing cloud droplet number concentration (CDNC), on cloud and precipitation evolution can

  4. 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.

  5. Numerical investigation of high-pressure combustion in rocket engines using Flamelet/Progress-variable models

    E-Print Network [OSTI]

    Coclite, A; De Palma, P; Pascazio, G

    2015-01-01

    The present paper deals with the numerical study of high pressure LOx/H2 or LOx/hydrocarbon combustion for propulsion systems. The present research effort is driven by the continued interest in achieving low cost, reliable access to space and more recently, by the renewed interest in hypersonic transportation systems capable of reducing time-to-destination. Moreover, combustion at high pressure has been assumed as a key issue to achieve better propulsive performance and lower environmental impact, as long as the replacement of hydrogen with a hydrocarbon, to reduce the costs related to ground operations and increase flexibility. The current work provides a model for the numerical simulation of high- pressure turbulent combustion employing detailed chemistry description, embedded in a RANS equations solver with a Low Reynolds number k-omega turbulence model. The model used to study such a combustion phenomenon is an extension of the standard flamelet-progress-variable (FPV) turbulent combustion model combined ...

  6. Modelling and Numerical Simulation of Gas Migration in a Nuclear Waste Repository

    E-Print Network [OSTI]

    Bourgeat, Alain; Smai, Farid

    2010-01-01

    We present a compositional compressible two-phase, liquid and gas, flow model for numerical simulations of hydrogen migration in deep geological radioactive waste repository. This model includes capillary effects and the gas diffusivity. The choice of the main variables in this model, Total or Dissolved Hydrogen Mass Concentration and Liquid Pressure, leads to a unique and consistent formulation of the gas phase appearance and disappearance. After introducing this model, we show computational evidences of its adequacy to simulate gas phase appearance and disappearance in different situations typical of underground radioactive waste repository.

  7. A numerical model for the gamma-ray emission of the microquasar LS 5039

    E-Print Network [OSTI]

    V. Bosch-Ramon; J. M. Paredes

    2004-07-01

    The possible association between the microquasar LS 5039 and the EGRET source 3EG J1824-1514 suggests that microquasars could also be sources of high energy gamma-rays. In this paper, we explore, with a detailed numerical model, if this system can produce the emission detected by EGRET (>100 MeV) through inverse Compton (IC) scattering. Our numerical approach considers a population of relativistic electrons entrained in a cylindrical inhomogeneous jet, interacting with both the radiation and the magnetic fields, taking into account the Thomson and Klein-Nishina regimes of interaction. The computed spectrum reproduces the observed spectral characteristics at very high energy.

  8. 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.

  9. Numerical upscaling for the eddy-current model with stochastic magnetic materials

    SciTech Connect (OSTI)

    Eberhard, Jens P. [Computer Simulation Technology, Bad Nauheimer Strasse, 19, D-64289 Darmstadt (Germany)], E-mail: jens.eberhard@cst.com; Popovic, Dan [Simulation in Technology, University of Heidelberg, Im Neuenheimer Feld 368, D-69120 Heidelberg (Germany)], E-mail: dan.popovic@stud.uni-heidelberg.de; Wittum, Gabriel [Simulation in Technology, University of Heidelberg, Im Neuenheimer Feld 368, D-69120 Heidelberg (Germany)], E-mail: wittum@uni-hd.de

    2008-04-01

    This paper deals with the upscaling of the time-harmonic Maxwell equations for heterogeneous media. We analyze the eddy-current approximation of Maxwell's equations to describe the electric field for heterogeneous, isotropic magnetic materials. The magnetic permeability of the materials is assumed to have random heterogeneities described by a Gaussian random field. We apply the so-called Coarse Graining method to develop a numerical upscaling of the eddy-current model. The upscaling uses filtering and averaging procedures in Fourier space which results in a formulation of the eddy-current model on coarser resolution scales where the influence of sub-scale fluctuations is modeled by effective scale- and space-dependent reluctivity tensors. The effective reluctivity tensors can be obtained by solving local partial differential equations which contain a Laplacian as well as a curl-curl operator. We present a computational method how the equation of the combined operators can be discretized and solved numerically using an extended variational formulation compared to standard discretizations. We compare the results of the numerical upscaling of the eddy-current model with theoretical results of Eberhard [J.P. Eberhard, Upscaling for the time-harmonic Maxwell equations with heterogeneous magnetic materials, Physical Review E 72 (3), (2005)] and obtain a very good agreement.

  10. Numerical study of a slip-link model for polymer melts and nanocomposites

    E-Print Network [OSTI]

    Diego Delbiondo; Elian Masnada; Samy Merabia; Marc Couty; Jean-Louis Barrat

    2013-06-10

    We present a numerical study of the slip link model introduced by Likhtman for describing the dy- namics of dense polymer melts. After reviewing the technical aspects associated with the implemen- tation of the model, we extend previous work in several directions. The dependence of the relaxation modulus with the slip link density and the slip link stiffness is reported. Then the nonlinear rheolog- ical properties of the model, for a particular set of parameters, are explored. Finally, we introduce excluded volume interactions in a mean field such as manner in order to describe inhomogeneous systems, and we apply this description to a simple nanocomposite model. With this extension, the slip link model appears as a simple and generic model of a polymer melt, that can be used as an alternative to molecular dynamics for coarse grained simulations of complex polymeric systems.

  11. 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

  12. Numerical-Model Investigation of the Hydrothermal Regime of a Straight-Through Shallow Cooling Pond

    SciTech Connect (OSTI)

    Sokolov, A. S. [JSC 'VNIIG im. B. E. Vedeneeva' (Russian Federation)] [JSC 'VNIIG im. B. E. Vedeneeva' (Russian Federation)

    2013-11-15

    A mathematic model based on solution of hydrodynamics and heat-transfer equations by the finite-element method is constructed to predict the hydrothermal regime of a straight-through shallow cooling pond, which provides cooling circulating water to a repository of spent nuclear fuel. Numerical experiments made it possible to evaluate the influence exerted by wind conditions and flow rate of water in the river on the temperature of the circulating water.

  13. Numerical modeling of hydraulic fracture problem in permeable medium using cohesive zone model

    E-Print Network [OSTI]

    Paris-Sud XI, Université de

    or radioactive waste [30], geothermal Corresponding author Email address: benoit.carrier@enpc.fr (Benoit Carrier processes. During the last sixty years, numerous papers [3, 7, 14, 21, 35, 34, 41, 25, 27, 22, 38, 36, 1. In the recent years, a scaling and asymptotic framework was built to determine the influence of the physical

  14. 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.

  15. gamma-ray emission from microquasars: a numerical model for LSI+61 303

    E-Print Network [OSTI]

    Valenti Bosch-Ramon; Josep M. Paredes

    2004-07-01

    We explore the possible association between the microquasar LSI+61 303 and the EGRET source 2CG 135+01/3EG J0241$+$6103 by studying, with a detailed numerical model, whether this system can produce the emission and the variability detected by EGRET ($>$100 MeV) through inverse Compton (IC) scattering. Our numerical approach considers a population of relativistic electrons entrained in a cylindrical inhomogeneous jet, interacting with both the radiation and the magnetic fields, taking into account the Thomson and Klein-Nishina regimes of interaction. Our results reproduce the observed spectral characteristics and variability at $\\gamma$-rays, thus strengthening the identification of LSI+61 303 as a high-energy $\\gamma$-ray source.

  16. Numerical modelling of sandstone uniaxial compression test using a mix-mode cohesive fracture model

    E-Print Network [OSTI]

    Gui, Yilin; Kodikara, Jayantha

    2015-01-01

    A mix-mode cohesive fracture model considering tension, compression and shear material behaviour is presented, which has wide applications to geotechnical problems. The model considers both elastic and inelastic displacements. Inelastic displacement comprises fracture and plastic displacements. The norm of inelastic displacement is used to control the fracture behaviour. Meantime, a failure function describing the fracture strength is proposed. Using the internal programming FISH, the cohesive fracture model is programmed into a hybrid distinct element algorithm as encoded in Universal Distinct Element Code (UDEC). The model is verified through uniaxial tension and direct shear tests. The developed model is then applied to model the behaviour of a uniaxial compression test on Gosford sandstone. The modelling results indicate that the proposed cohesive fracture model is capable of simulating combined failure behaviour applicable to rock.

  17. The NINJA-2 project: Detecting and characterizing gravitational waveforms modelled using numerical binary black hole simulations

    E-Print Network [OSTI]

    The LIGO Scientific Collaboration; the Virgo Collaboration; the NINJA-2 Collaboration; :; J. Aasi; B. P. Abbott; R. Abbott; T. Abbott; M. R. Abernathy; T. Accadia; F. Acernese; K. Ackley; C. Adams; T. Adams; P. Addesso; R. X. Adhikari; C. Affeldt; M. Agathos; N. Aggarwal; O. D. Aguiar; A. Ain; P. Ajith; A. Alemic; B. Allen; A. Allocca; D. Amariutei; M. Andersen; R. Anderson; S. B. Anderson; W. G. Anderson; K. Arai; M. C. Araya; C. Arceneaux; J. Areeda; S. M. Aston; P. Astone; P. Aufmuth; C. Aulbert; L. Austin; B. E. Aylott; S. Babak; P. T. Baker; G. Ballardin; S. W. Ballmer; J. C. Barayoga; M. Barbet; B. C. Barish; D. Barker; F. Barone; B. Barr; L. Barsotti; M. Barsuglia; M. A. Barton; I. Bartos; R. Bassiri; A. Basti; J. C. Batch; J. Bauchrowitz; Th. S. Bauer; B. Behnke; M. Bejger; M. G. Beker; C. Belczynski; A. S. Bell; C. Bell; G. Bergmann; D. Bersanetti; A. Bertolini; J. Betzwieser; P. T. Beyersdorf; I. A. Bilenko; G. Billingsley; J. Birch; S. Biscans; M. Bitossi; M. A. Bizouard; E. Black; J. K. Blackburn; L. Blackburn; D. Blair; S. Bloemen; M. Blom; O. Bock; T. P. Bodiya; M. Boer; G. Bogaert; C. Bogan; C. Bond; F. Bondu; L. Bonelli; R. Bonnand; R. Bork; M. Born; V. Boschi; Sukanta Bose; L. Bosi; C. Bradaschia; P. R. Brady; V. B. Braginsky; M. Branchesi; J. E. Brau; T. Briant; D. O. Bridges; A. Brillet; M. Brinkmann; V. Brisson; A. F. Brooks; D. A. Brown; D. D. Brown; F. Brückner; S. Buchman; T. Bulik; H. J. Bulten; A. Buonanno; R. Burman; D. Buskulic; C. Buy; L. Cadonati; G. Cagnoli; J. Calderón Bustillo; E. Calloni; J. B. Camp; P. Campsie; K. C. Cannon; B. Canuel; J. Cao; C. D. Capano; F. Carbognani; L. Carbone; S. Caride; A. Castiglia; S. Caudill; M. Cavaglià; F. Cavalier; R. Cavalieri; C. Celerier; G. Cella; C. Cepeda; E. Cesarini; R. Chakraborty; T. Chalermsongsak; S. J. Chamberlin; S. Chao; P. Charlton; E. Chassande-Mottin; X. Chen; Y. Chen; A. Chincarini; A. Chiummo; H. S. Cho; J. Chow; N. Christensen; Q. Chu; S. S. Y. Chua; S. Chung; G. Ciani; F. Clara; J. A. Clark; F. Cleva; E. Coccia; P. -F. Cohadon; A. Colla; C. Collette; M. Colombini; L. Cominsky; M. Constancio Jr.; A. Conte; D. Cook; T. R. Corbitt; M. Cordier; N. Cornish; A. Corpuz; A. Corsi; C. A. Costa; M. W. Coughlin; S. Coughlin; J. -P. Coulon; S. Countryman; P. Couvares; D. M. Coward; M. Cowart; D. C. Coyne; R. Coyne; K. Craig; J. D. E. Creighton; S. G. Crowder; A. Cumming; L. Cunningham; E. Cuoco; K. Dahl; T. Dal Canton; M. Damjanic; S. L. Danilishin; S. D'Antonio; K. Danzmann; V. Dattilo; H. Daveloza; M. Davier; G. S. Davies; E. J. Daw; R. Day; T. Dayanga; G. Debreczeni; J. Degallaix; S. Deléglise; W. Del Pozzo; T. Denker; T. Dent; H. Dereli; V. Dergachev; R. De Rosa; R. T. DeRosa; R. DeSalvo; S. Dhurandhar; M. Díaz; L. Di Fiore; A. Di Lieto; I. Di Palma; A. Di Virgilio; A. Donath; F. Donovan; K. L. Dooley; S. Doravari; S. Dossa; R. Douglas; T. P. Downes; M. Drago; R. W. P. Drever; J. C. Driggers; Z. Du; S. Dwyer; T. Eberle; T. Edo; M. Edwards; A. Effler; H. Eggenstein; P. Ehrens; J. Eichholz; S. S. Eikenberry; G. Endr?czi; R. Essick; T. Etzel; M. Evans; T. Evans; M. Factourovich; V. Fafone; S. Fairhurst; Q. Fang; S. Farinon; B. Farr; W. M. Farr; M. Favata; H. Fehrmann; M. M. Fejer; D. Feldbaum; F. Feroz; I. Ferrante; F. Ferrini; F. Fidecaro; L. S. Finn; I. Fiori; R. P. Fisher; R. Flaminio; J. -D. Fournier; S. Franco; S. Frasca; F. Frasconi; M. Frede; Z. Frei; A. Freise; R. Frey; T. T. Fricke; P. Fritschel; V. V. Frolov; P. Fulda; M. Fyffe; J. Gair; L. Gammaitoni; S. Gaonkar; F. Garufi; N. Gehrels; G. Gemme; E. Genin; A. Gennai; S. Ghosh; J. A. Giaime; K. D. Giardina; A. Giazotto; C. Gill; J. Gleason; E. Goetz; R. Goetz; L. Gondan; G. González; N. Gordon; M. L. Gorodetsky; S. Gossan; S. Goßler; R. Gouaty; C. Gräf; P. B. Graff; M. Granata; A. Grant; S. Gras; C. Gray; R. J. S. Greenhalgh; A. M. Gretarsson; P. Groot; H. Grote; K. Grover; S. Grunewald; G. M. Guidi; C. Guido; K. Gushwa; E. K. Gustafson; R. Gustafson; D. Hammer; G. Hammond; M. Hanke; J. Hanks; C. Hanna; J. Hanson; J. Harms; G. M. Harry; I. W. Harry; E. D. Harstad; M. Hart; M. T. Hartman; C. -J. Haster; K. Haughian; A. Heidmann; M. Heintze; H. Heitmann; P. Hello; G. Hemming; M. Hendry; I. S. Heng; A. W. Heptonstall; M. Heurs; M. Hewitson; S. Hild; D. Hoak; K. A. Hodge; K. Holt; S. Hooper; P. Hopkins; D. J. Hosken; J. Hough; E. J. Howell; Y. Hu; B. Hughey; S. Husa; S. H. Huttner; M. Huynh; T. Huynh-Dinh; D. R. Ingram; R. Inta; T. Isogai; A. Ivanov; B. R. Iyer; K. Izumi; M. Jacobson; E. James; H. Jang; P. Jaranowski; Y. Ji; F. Jiménez-Forteza; W. W. Johnson; D. I. Jones; R. Jones; R. J. G. Jonker; L. Ju; Haris K; P. Kalmus; V. Kalogera; S. Kandhasamy; G. Kang; J. B. Kanner; J. Karlen; M. Kasprzack; E. Katsavounidis; W. Katzman; H. Kaufer; K. Kawabe; F. Kawazoe; F. Kéfélian; G. M. Keiser; D. Keitel; D. B. Kelley; W. Kells; A. Khalaidovski

    2014-01-05

    The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational-wave astrophysics communities. The purpose of NINJA is to study the ability to detect gravitational waves emitted from merging binary black holes and recover their parameters with next-generation gravitational-wave observatories. We report here on the results of the second NINJA project, NINJA-2, which employs 60 complete binary black hole hybrid waveforms consisting of a numerical portion modelling the late inspiral, merger, and ringdown stitched to a post-Newtonian portion modelling the early inspiral. In a "blind injection challenge" similar to that conducted in recent LIGO and Virgo science runs, we added 7 hybrid waveforms to two months of data recolored to predictions of Advanced LIGO and Advanced Virgo sensitivity curves during their first observing runs. The resulting data was analyzed by gravitational-wave detection algorithms and 6 of the waveforms were recovered with false alarm rates smaller than 1 in a thousand years. Parameter estimation algorithms were run on each of these waveforms to explore the ability to constrain the masses, component angular momenta and sky position of these waveforms. We also perform a large-scale monte-carlo study to assess the ability to recover each of the 60 hybrid waveforms with early Advanced LIGO and Advanced Virgo sensitivity curves. Our results predict that early Advanced LIGO and Advanced Virgo will have a volume-weighted average sensitive distance of 300Mpc (1Gpc) for $10M_{\\odot}+10M_{\\odot}$ ($50M_{\\odot}+50M_{\\odot}$) binary black hole coalescences. We demonstrate that neglecting the component angular momenta in the waveform models used in matched-filtering will result in a reduction in sensitivity for systems with large component angular momenta. [Abstract abridged for ArXiv, full version in PDF

  18. 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.

  19. A phase screen model for simulating numerically the propagation of a laser beam in rain

    SciTech Connect (OSTI)

    Lukin, I P; Rychkov, D S; Falits, A V [Institute of Atmospheric Optics, Siberian Branch, Russian Academy of Sciences, Tomsk (Russian Federation); Lai, Kin S; Liu, Min R [DSO National Laboratories 20 (Singapore)

    2009-09-30

    The method based on the generalisation of the phase screen method for a continuous random medium is proposed for simulating numerically the propagation of laser radiation in a turbulent atmosphere with precipitation. In the phase screen model for a discrete component of a heterogeneous 'air-rain droplet' medium, the amplitude screen describing the scattering of an optical field by discrete particles of the medium is replaced by an equivalent phase screen with a spectrum of the correlation function of the effective dielectric constant fluctuations that is similar to the spectrum of a discrete scattering component - water droplets in air. The 'turbulent' phase screen is constructed on the basis of the Kolmogorov model, while the 'rain' screen model utiises the exponential distribution of the number of rain drops with respect to their radii as a function of the rain intensity. Theresults of the numerical simulation are compared with the known theoretical estimates for a large-scale discrete scattering medium. (propagation of laser radiation in matter)

  20. Universitt StuttgartInstitut fr Wasserbau, Lehrstuhl fr Hydromechanik und Hydrosystemmodellierung Workshop on Numerical Models for Carbon Dioxide Storage in Geological Formations

    E-Print Network [OSTI]

    Cirpka, Olaf Arie

    Hydrosystemmodellierung Workshop on Numerical Models for Carbon Dioxide Storage in Geological Formations 1/16 Modelling April 2008 Workshop on Numerical Models for Carbon Dioxide Storage in Geological Formations #12 on Numerical Models for Carbon Dioxide Storage in Geological Formations 2/16 CO2 leakage mitigation using

  1. Stochastic models and numerical algorithms for a class of regulatory gene networks

    E-Print Network [OSTI]

    Thomas Fournier; Jean-Pierre Gabriel; Christian Mazza; Jerome Pasquier; Jose Galbete; Nicolas Mermod

    2008-10-01

    Regulatory gene networks contain generic modules like those involving feedback loops, which are essential for the regulation of many biological functions. We consider a class of self-regulated genes which are the building blocks of many regulatory gene networks, and study the steady state distributions of the associated Gillespie algorithm by providing efficient numerical algorithms. We also study a regulatory gene network of interest in synthetic biology and in gene therapy, using mean-field models with time delays. Convergence of the related time-nonhomogeneous Markov chain is established for a class of linear catalytic networks with feedback loops

  2. Numerical modeling of the high-temperature geothermal system of Amatitlan, Guatemala

    SciTech Connect (OSTI)

    Pham, M.; Menzies, A.J.; Sanyal, S.K. [GeothermEx, Inc., Richmond, CA (United States)] [and others

    1996-12-31

    Using the conceptual model of the Amatitlin geothermal system as a guide, a 14 by 10 km numerical model was set up with five layers and a total of 1,220 grid blocks. The number of grid blocks varies from layer to layer, with finer resolution in the two production layers. Boundary blocks for conductive heat transfer, heat and mass recharge, shallow discharge and deep regional fluid flow were added to the model as appropriate. With an inflow of hot fluid (336{degrees}C) at a rate of 450 tones per hour, an inflow of cool fluid (50{degrees}C) at a rate of 620 tons per hour and permeabilities varying from 1.5 to 50 md, the initial state temperatures on all five levels were well matched by the numerical model. The initial state results indicate that highest permeability is present to the southeast of the productive wells, suggesting that considerable additional production potential exists in this area. A two-phase region exists in the upper levels of the reservoir and in the outflow plume; the latter is consistent with the presence of fumaroles along the shoreline of Lago de Amatitlan. To further calibrate the numerical model, enthalpy and pressure data collected during tests of wells AMF-1 and -2 were matched, primarily by varying reservoir storage capacity. Good matches to the observed data were obtained for both wells, although sharp transients could not be matched. Modeling proceeded to the prediction stage under two scenarios: production of 12 and 25 MW (gross) for 30 years. 100% injection of separated water and condensate was assumed to take place in the vicinity of well AMF-3 in both cases. The results show an increase in enthalpy caused by expansion of the two-phase zone during the first year of production, followed by an enthalpy decline as reservoir liquid begins to migrate toward the production area. The two existing production wells (AMF-1 and -2) are predicted to be capable of maintaining output in the 12 MW case throughout the 30-year project life.

  3. Numerical modeling of mixed sediment resuspension, transport, and deposition during the March 1998 episodic events in southern Lake

    E-Print Network [OSTI]

    Numerical modeling of mixed sediment resuspension, transport, and deposition during the March 1998 sediment resuspension of mixed (cohesive plus noncohesive) sediment is developed and applied to quantitatively simulate the March 1998 resuspension events in southern Lake Michigan. Some characteristics

  4. Tsunami inundation modeling in constructed environments: A physical and numerical comparison of free-surface elevation, velocity, and momentum flux

    E-Print Network [OSTI]

    Lynett, Patrick

    Tsunami inundation modeling in constructed environments: A physical and numerical comparison April 2013 Available online 17 May 2013 Keywords: Tsunami Inundation Macro-roughness Benchmark COULWAVE Friction factor A laboratory benchmark test for tsunami inundation through an urban waterfront including

  5. 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 Bhatnagar–Gross–Krook (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.

  6. Langevin power curve analysis for numerical WEC models with new insights on high frequency power performance

    E-Print Network [OSTI]

    Mücke, Tanja A; Milan, Patrick; Peinke, Joachim

    2015-01-01

    Based on the Langevin equation it has been proposed to obtain power curves for wind turbines from high frequency data of wind speed measurements u(t) and power output P (t). The two parts of the Langevin approach, power curve and drift field, give a comprehensive description of the conversion dynamic over the whole operating range of the wind turbine. The method deals with high frequent data instead of 10 min means. It is therefore possible to gain a reliable power curve already from a small amount of data per wind speed. Furthermore, the method is able to visualize multiple fixed points, which is e.g. characteristic for the transition from partial to full load or in case the conversion process deviates from the standard procedures. In order to gain a deeper knowledge it is essential that the method works not only for measured data but also for numerical wind turbine models and synthetic wind fields. Here, we characterize the dynamics of a detailed numerical wind turbine model and calculate the Langevin power...

  7. Numerical modeling of roll structures in mesoscale vortexes over the Black Sea

    E-Print Network [OSTI]

    Iarova, D A

    2014-01-01

    This paper is a case study of horizontal atmospheric rolls that formed over the Black Sea on 16 August 2007. The rolls were discovered in WRF modeling results for a mesoscale cyclone that originated over the sea on 15 August 2007. The roll formation mechanisms, such as Rayleigh-Benard convective instability, dynamic instability, advection and stretching of vertical velocity field inhomogeneities, are considered. It is shown that indeed convective instability played an important role in the roll formation but dynamic instability did not occur. In order to distinguish other possible mechanisms of the roll formation numerical experiments were performed. In these experiments sea surface temperature in the initial conditions was decreased in order to prevent convective instability. Even though convective instability was suppressed roll-like structures still appeared in the modeling results, although their height and circulation velocity were smaller than in the control run. It was found that these structures were ...

  8. Final Report for "Accurate Numerical Models of the Secondary Electron Yield from Grazing-incidence Collisions".

    SciTech Connect (OSTI)

    Seth A Veitzer

    2008-10-21

    Effects of stray electrons are a main factor limiting performance of many accelerators. Because heavy-ion fusion (HIF) accelerators will operate in regimes of higher current and with walls much closer to the beam than accelerators operating today, stray electrons might have a large, detrimental effect on the performance of an HIF accelerator. A primary source of stray electrons is electrons generated when halo ions strike the beam pipe walls. There is some research on these types of secondary electrons for the HIF community to draw upon, but this work is missing one crucial ingredient: the effect of grazing incidence. The overall goal of this project was to develop the numerical tools necessary to accurately model the effect of grazing incidence on the behavior of halo ions in a HIF accelerator, and further, to provide accurate models of heavy ion stopping powers with applications to ICF, WDM, and HEDP experiments.

  9. Two-dimensional numerical model of underground oil-shale retorting

    SciTech Connect (OSTI)

    Travis, B.J.; Hommert, P.J.; Tyner, C.E.

    1983-01-01

    A two-dimensional numerical model of underground oil shale retorting, which fully couples retorting chemistry with fluid and heat flow, has been developed. The model solves the time-dependent, two-dimensional mass, momentum, and energy balance equations for a nine-component fluid (O/sub 2/, N/sub 2/, H/sub 2/, CO/sub 2/, CO, CH/sub 4/, CH/sub x/, H/sub 2/O, and oil). Water and oil can flow in the liquid and/or vapor phases. Retort chemistry includes kerogen pyrolysis, carbonate decomposition, char reactions, and combustion. Also, detailed modeling of heat flow and chemistry inside shale particles allows large rubble sizes as well as small sizes to be considered. The model is compared to one-dimensional experimental data obtained from Lawrence Livermore National Laboratory. The model can be used to examine the effect of two-dimensional variations in shale grade, rubble size, permeability, porosity, geometry, inflow gas composition, etc. on retorting efficiency and process optimization. A sample calculation is presented.

  10. 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.

  11. The effect of simplifying the building description on the numerical modeling of its thermal performance

    SciTech Connect (OSTI)

    Stetiu, C.

    1993-07-01

    A thermal building simulation program is a numerical model that calculates the response of the building envelopes to weather and human activity, simulates dynamic heating and cooling loads, and heating and cooling distribution systems, and models building equipment operation. The scope of the research is to supply the users of such programs with information about the dangers and benefits of simplifying the input to their models. The Introduction describes the advantages of modeling the heat transfer mechanisms in a building. The programs that perform this type of modeling have, however, limitations. The user is therefore often put in the situation of simplifying the floor plans of the building under study, but not being able to check the effects that this approximation introduces in the results of the simulation. Chapter 1 is a description of methods. It also introduces the floor plans for the office building under study and the ``reasonable`` floor plans simplifications. Chapter 2 presents DOE-2, the thermal building simulation program used in the sensitivity study. The evaluation of the accuracy of the DOE-2 program itself is also presented. Chapter 3 contains the sensitivity study. The complicated nature of the process of interpreting the temperature profile inside a space leads to the necessity of defining different building modes. The study compares the results from the model of the detailed building description with the results from the models of the same building having simplified floor plans. The conclusion is reached that a study of the effects of simplifying the floor plans of a building is important mainly for defining the cases in which this approximation is acceptable. Different results are obtained for different air conditioning/load regimes of the building. 9 refs., 24 figs.

  12. Numerical Modeling of Thermal EOR: Comprehensive Coupling of an AMR-Based Model

    E-Print Network [OSTI]

    Paris-Sud XI, Université de

    Flow and Geomechanics N. Guy*, G. Enchéry and G. Renard IFP Energies nouvelles, 1-4 avenue de Bois of Thermal EOR: Comprehensive Coupling of an AMR-Based Model of Thermal Fluid Flow and Geomechanics when both thermal fluid flow and geomechanics are coupled in order to take into account variations

  13. 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.

  14. Tidal Downsizing model. I. Numerical methods: saving giant planets from tidal disruptions

    E-Print Network [OSTI]

    Nayakshin, Sergei

    2014-01-01

    Tidal Downsizing (TD) is a recently developed planet formation theory that supplements the classical Gravitational disc Instability (GI) model with planet migration inward and tidal disruptions of GI fragments in the inner regions of the disc. Numerical methods for a detailed population synthesis of TD planets are presented here. As an example application, the conditions under which GI fragments collapse faster than they migrate into the inner $a\\sim$ few AU disc are considered. It is found that most gas fragments are tidally or thermally disrupted unless (a) their opacity is $\\sim 3$ orders of magnitude less than the interstellar dust opacity at metallicities typical of the observed giant planets, or (b) the opacity is high but the fragments accrete large dust grains (pebbles) from the disc. Case (a) models produce very low mass solid cores ($M_{\\rm core} < 0.1$ Earth masses) and follow a negative correlation of giant planet frequency with host star metallicity. In contrast, case (b) models produce massiv...

  15. Numerical modeling of mixed sediment resuspension, transport, and deposition during the March 1998 episodic events in southern Lake

    E-Print Network [OSTI]

    Stroud, Jonathan

    Numerical modeling of mixed sediment resuspension, transport, and deposition during the March 1998 2006; published 17 February 2007. [1] A two-dimensional sediment transport model capable of simulating sediment resuspension of mixed (cohesive plus noncohesive) sediment is developed and applied

  16. Introduction Flattening the Earth Continuation procedure Flat Earth Numerical simulations Continuation from a flat to a round Earth model

    E-Print Network [OSTI]

    Boyer, Edmond

    Introduction Flattening the Earth Continuation procedure Flat Earth Numerical simulations Continuation from a flat to a round Earth model in the coplanar orbit transfer problem M. Cerf1, T. Haberkorn, SADCO 2011, March 2nd M. Cerf, T. Haberkorn, E. Tr´elat Continuation from a flat to a round Earth model

  17. 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.

  18. Asymptotics and numerical efficiency of the Allen-Cahn model for phase interfaces with low energy in solids

    E-Print Network [OSTI]

    Hans-Dieter Alber

    2015-05-20

    The accurate simulation of phase interfaces in solids requires small model error and small numerical error. If a phase field model is used and the interface carries low interface energy, then the model error is only small if the interface width in the model is chosen small. Yet, for effective numerical computation the interface width should be large. Choosing the parameters, which determine the width, is therefore an optimality problem. We study this problem for the Allen-Cahn equation coupled to the elasticity equations by constructing an asymptotic solution of second order, which yields an expansion for the kinetic relation of the model. This expansion determines the choice of the parameters, however only if the difference between the expansion and the exact kinetic relation is uniformly small with respect to a second parameter controlling the interface energy. To show this uniformity we determine the asymptotics with respect to this second parameter by scaling of the model equations. Our investigations are formal.

  19. 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 acceleration (e.g., gravity) and their collision is induced by the difference in their terminal velocities. We point out an important implication of this discovery for dust growth in protoplanetary disks.

  20. Numerical modeling of the initial state and matching of well test data from the Zunil geothermal field, Guatemala

    SciTech Connect (OSTI)

    Menzies, A.J.; Granados, E.E.; Sanyal, .K.; Merida-I., L.; Caicedo-A, A.

    1991-01-01

    A significant amount of geoscientific and reservoir engineering data have been collected from the Zunil geothermal field since 1973. The data have been used to define a conceptual model for the field which has formed the basis for the construction of a three dimensional numerical simulation model. The numerical model has successfully matched both the initial state of the reservoir, as indicated by subsurface temperature and pressure distributions within the presently drilled area, and available well test data. The well test data include short and long term discharge tests and a comprehensive pressure interference test. Calibration of the model will continue during 1991 when the results from drilling and testing of three additional deep wells are available. The model will then be used to study various long term production scenarios for the proposed 15 MW power development.

  1. The VOLNA code for the numerical modelling of tsunami waves: generation, propagation and

    E-Print Network [OSTI]

    of the finite volume scheme implemented in the code. We explain the numerical treatment of the wet is decided on the base of inundation maps which are produced with this type of numerical tools. Finally we and the perspectives for future research presented. Key words: tsunami waves, shallow water equations, tsunami

  2. Numerical modeling of alongshore sediment transport and shoreline change along the Galveston coast 

    E-Print Network [OSTI]

    Sitanggang, Khairil Irfan

    2005-02-17

    coast for the period of 2001-2011. The first assessment shows that GENESIS has a numerical error which tends to be large for low energy wave (small breaking wave height) and large breaking wave angle. This numerical inaccuracy cannot be neglected...

  3. Numerical simulation of the 2011 Tohoku tsunami: Comparison with field observations and sensitivity to model parameters

    E-Print Network [OSTI]

    Kirby, James T.

    Numerical simulation of the 2011 Tohoku tsunami: Comparison with field observations and sensitivity history, created a major tsunami that caused numerous deaths and enormous destruction on the nearby Hon- shu coast. Various tsunami sources were developed for this event, based on inverting seismic or GPS

  4. A controlled distributed parameter model for a fluid-flexible structure system: numerical simulations and experiment validations

    E-Print Network [OSTI]

    Baudouin, Lucie

    A controlled distributed parameter model for a fluid-flexible structure system: numerical consider the problem of active reduction of vibrations in a fluid-flexible structure system and the sloshing of the fuel inside the wing's tank. The control is performed using piezoelectric patches

  5. 2250 IEEE TRANSACTIONS ON PLASMA SCIENCE, VOL. 37, NO. 11, NOVEMBER 2009 A Numerical Model of Plasma-Actuator Effects

    E-Print Network [OSTI]

    Huang, Xun

    therein. Most recent ap- plications include high-speed flow control using localized arc filament plasma2250 IEEE TRANSACTIONS ON PLASMA SCIENCE, VOL. 37, NO. 11, NOVEMBER 2009 A Numerical Model of Plasma-Actuator Effects in Flow-Induced Noise Control Ed Peers, Xun Huang, and Xinfu Luo Abstract

  6. IMA Journal of Numerical Analysis (1988) 8, 415-433 Finite Element Methods for a Model for Full Waveform Acoustic

    E-Print Network [OSTI]

    Harrison, Mark

    1988-01-01

    IMA Journal of Numerical Analysis (1988) 8, 415-433 Finite Element Methods for a Model for Full are given and then a discrete-time, explicit finite element procedure is defined and analysed, with finite on the existence, uniqueness and finite element approximation of the solution of Biot's equations were given in [15

  7. LETTER Earth Planets Space, 63, 847851, 2011 Three-dimensional numerical modeling of tsunami-related internal gravity

    E-Print Network [OSTI]

    Occhipinti, Giovanni "Ninto"

    LETTER Earth Planets Space, 63, 847­851, 2011 Three-dimensional numerical modeling of tsunami, 2011; Accepted June 30, 2011; Online published September 27, 2011) The tremendous tsunami following, to reproduce the tsunami signature observed in the airglow by the imager located in Hawaii and clearly showing

  8. Segmentation and Tracking of Mesoscale Eddies in Numeric Ocean Models Vishal Sood, Bin John, Ramprasad Balasubramanian and Amit Tandon*

    E-Print Network [OSTI]

    Tandon, Amit

    Segmentation and Tracking of Mesoscale Eddies in Numeric Ocean Models Vishal Sood, Bin John suggested that the mesoscale eddies and mesoscale features play a strong role in carrying heat poleward oceanographers an invaluable tool to assess mesoscale eddies and the Lagrangian characteristics of this mesoscale

  9. Numerical study of the mass spectrum in the 2D O(3) sigma model with a theta term

    E-Print Network [OSTI]

    B. Alles; A. Papa

    2007-11-12

    It has been conjectured that the mass spectrum of the O(3) non-linear sigma model with a theta term in 2 dimensions may possess an excited state, which decays when theta is lowered from pi below a critical value. Since the direct numerical investigation of the model is prevented by a sign problem, we try to infer some information on the mass spectrum at real theta by studying the model at imaginary theta via analytic continuation. A modified Swendsen-Wang cluster algorithm has been introduced to simulate the model with the theta term.

  10. 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.

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

    SciTech Connect (OSTI)

    Vršnak, B.; Žic, T.; Dumbovi?, M.; Temmer, M.; Möstl, 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.09±9.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.

  12. A high-resolution, cloud-assimilating numerical weather prediction model for solar irradiance forecasting

    E-Print Network [OSTI]

    Mathiesen, Patrick; Collier, Craig; Kleissl, Jan

    2013-01-01

    cycle:  The RUC.  Monthly Weather  Review.   132, 495?518.  th  Conference on  Numerical Weather Prediction.   American closure schemes.   Monthly Weather Review.   122, 927?945.  

  13. Soil Liquefaction–Induced Uplift of Underground Structures: Physical and Numerical Modeling

    E-Print Network [OSTI]

    Chian, Siau Chen; Tokimatsu, Kohji; M.ASCE; Madabhushi, Santana Phani Gopal

    2014-07-11

    in the liquefiable soil deposit. These consistencies with field observations clearly demonstrate and pave the prospects of applying numerical and/or experimental analyses for geotechnical problems associated with the floatation of underground structures...

  14. Numerical modeling of species transport in turbulent flow and experimental study on aerosol sampling 

    E-Print Network [OSTI]

    Vijayaraghavan, Vishnu Karthik

    2007-04-25

    Numerical simulations were performed to study the turbulent mixing of a scalar species in straight tube, single and double elbow flow configurations. Different Reynolds Averaged Navier Stokes (RANS) and Large Eddy Simulation ...

  15. Grid-Based Surrogate Reservoir Modeling (SRM) for Fast Track Analysis of Numerical Reservoir Simulation Models at the Grid block Level

    E-Print Network [OSTI]

    Mohaghegh, Shahab

    capability of being able to replicate the pressure and saturation distribution throughout the reservoirSPE 153844 Grid-Based Surrogate Reservoir Modeling (SRM) for Fast Track Analysis of Numerical Reservoir Simulation Models at the Grid block Level Shahab D. Mohaghegh, West Virginia University

  16. Numerical method toward fast and accurate calculation of dilute quantum gas using Uehling-Uhlenbeck model equation

    E-Print Network [OSTI]

    Yano, Ryosuke

    2015-01-01

    We consider a fast and accurate numerical method to calculate the dilute quantum gas. To attain our aim, the Uehing-Uhlenbeck (U-U) model equation is calculated using the direct simulation Monte Carlo (DSMC) method. The proposed algorithm in the DSMC method enables us to obtain the accurate thermalization using a small number of sample particles and calculate the dilute quantum gas dynamics in practical time. Finally, we apply the proposed numerical method to the calculation of the viscosity coefficient on the basis of Green-Kubo expression.

  17. Qualification of the ITER CS Quench Detection System using Numerical Modeling

    SciTech Connect (OSTI)

    Martovetsky, Nicolai N [ORNL; Radovinsky, Alexey L [ORNL

    2013-01-01

    Abstract The ITER Central Solenoid (CS) magnet needs to be protected against overheating of the conductor in the event of the occurrence of a normal zone (NZ). Due to a large amount of stored energy and slow NZ propagation, the NZ needs to be detected and the switchyard needs to open the breakers within 2 s after detection of the NZ. The CS will be discharged on a dump resistor with a time constant of 7.5 s. During operation of the CS and its interaction with the poloidal field (PF) coils and plasma current, the CS experiences large inductive voltages from multiple sources, including nonlinear signals from eddy currents in the vacuum vessel and plasma current variation, that make the task of detecting the resistive signal even more difficult. This inductive voltage needs to be cancelled by quench detection (QD) hardware (e.g., bridges, converters, filters, processors) and appropriate processing of the QD signals to reliably detect NZ initiation and propagation. Two redundant schemes are proposed as the baseline for the CS QD System: 1) A scheme with Regular Voltage Taps (RVT) from triads of Double Pancakes (DP) supplemented by Central Difference Averaging (CDA) and by digital suppression of the inductive voltage from all active coils (the CS and PF coils). Voltage taps are taken from helium outlets at the CS outer diameter. 2)A scheme with Cowound Voltage Taps (CVT) taken from cowound wires routed from the helium inlet at the CS inner diameter. Summary of results of the numerical modeling of the performance of both baseline CS QD systems is presented in this paper. Index Terms Quench detection, Central Solenoid, ITER

  18. 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.

  19. Observation and numerical modeling of chromospheric evaporation during the impulsive phase of a solar flare

    E-Print Network [OSTI]

    Imada, Shinsuke; Watanabe, Tetsuya

    2015-01-01

    We have studied the chromospheric evaporation flow during the impulsive phase of the flare by using the Hinode/EIS observation and 1D hydrodynamic numerical simulation coupled to the time-dependent ionization. The observation clearly shows that the strong redshift can be observed at the base of the flaring loop only during the impulsive phase. We performed two different numerical simulations to reproduce the strong downflows in FeXII and FeXV during the impulsive phase. By changing the thermal conduction coefficient, we carried out the numerical calculation of chromospheric evaporation in the thermal conduction dominant regime (conductivity coefficient kappa0 = classical value) and the enthalpy flux dominant regime (kappa0 = 0.1 x classical value). The chromospheric evaporation calculation in the enthalpy flux dominant regime could reproduce the strong redshift at the base of the flare during the impulsive phase. This result might indicate that the thermal conduction can be strongly suppressed in some cases o...

  20. A three-dimensional numerical model of a micro laminar flow fuel cell with a bridge-shaped microchannel cross-section

    E-Print Network [OSTI]

    Kenis, Paul J. A.

    A three-dimensional numerical model of a micro laminar flow fuel cell with a bridge: Membraneless fuel cell Laminar flow fuel cell Numerical model Convection-diffusion equations Electrode kinetics equations COMSOL a b s t r a c t The operation of a laminar flow fuel cell (LFFC) involves complex interplay

  1. Remarks on numerical simulation of the PKN model of hydrofracturing in proper variables. Various leak-off regimes

    E-Print Network [OSTI]

    Kusmierczyk, P; Wrobel, M

    2012-01-01

    The problem of hydraulic fracture for the PKN model is considered within the framework of approach presented recently by Linkov (2011). The modified formulation is further enhanced by employing an improved regularized boundary condition near the crack tip. This increases solution accuracy especially for singular leak-off regimes. A new dependent variable having clear physical sense is introduced. A comprehensive analysis of numerical algorithms based on various dependent variables is provided.

  2. Numerical study on transient heat transfer under soil with plastic mulch in agriculture applications using a nonlinear finite element model

    E-Print Network [OSTI]

    De Castro, Carlos Armando

    2011-01-01

    In this paper is developed a simple mathematical model of transient heat transfer under soil with plastic mulch in order to determine with numerical studies the influence of different plastic mulches on the soil temperature and the evolutions of temperatures at different depths with time. The governing differential equations are solved by a Galerkin Finite Element Model, taking into account the nonlinearities due to radiative heat exchange between the soil surface, the plastic mulch and the atmosphere. The model was validated experimentally giving good approximation of the model to the measured data. Simulations were run with the validated model in order to determine the optimal combination of mulch optical properties to maximize the soil temperature with a Taguchi's analysis, proving that the material most used nowadays in Colombia is not the optimal and giving quantitative results of the properties the optimal mulch must possess.

  3. THE STATE OF THE ART OF NUMERICAL MODELING OF THERMOHYDROLOGIC FLOW IN FRACTURED ROCK MASSES

    E-Print Network [OSTI]

    Wang, J.S.Y.

    2013-01-01

    disposal of radioactive waste: The Sandia waste isolation flow and transport (SWIFT) model: Sandia Laboratories Report

  4. Mathematical Models and Numerical Simulations of phase change in Lagrangian and Eulerian descriptions

    E-Print Network [OSTI]

    Quiros Fonseca, Luis Alonso

    2012-12-31

    This thesis presents development of mathematical models for liquid-solid and solid-liquid phase change phenomena in Lagrangian and Eulerian descriptions. The mathematical models are derived by assuming a smooth interface ...

  5. THE STATE OF THE ART OF NUMERICAL MODELING OF THERMOHYDROLOGIC FLOW IN FRACTURED ROCK MASSES

    E-Print Network [OSTI]

    Wang, J.S.Y.

    2013-01-01

    improving production by hydraulic fracturing 8 the focus otfor fractures. (d) Hydraulic Fracturing: The model has been

  6. Atmospheric Test Models and Numerical Experiments for the Simulation of the Global Distributions of Weather Data Transponders III. Horizontal Distributions

    SciTech Connect (OSTI)

    Molenkamp, C.R.; Grossman, A.

    1999-12-20

    A network of small balloon-borne transponders which gather very high resolution wind and temperature data for use by modern numerical weather predication models has been proposed to improve the reliability of long-range weather forecasts. The global distribution of an array of such transponders is simulated using LLNL's atmospheric parcel transport model (GRANTOUR) with winds supplied by two different general circulation models. An initial study used winds from CCM3 with a horizontal resolution of about 3 degrees in latitude and longitude, and a second study used winds from NOGAPS with a 0.75 degree horizontal resolution. Results from both simulations show that reasonable global coverage can be attained by releasing balloons from an appropriate set of launch sites.

  7. 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.

  8. Atmospheric test models and numerical experiments for the simulation of the global distribution of weather data transponders

    SciTech Connect (OSTI)

    Grossman, A; Molenkamp, C R

    1999-08-25

    A proposal has been made to establish a high density global network of atmospheric micro transponders to record time, temperature, and wind data with time resolution of {le} 1 minute, temperature accuracy of {+-} 1 K, spatial resolution no poorer than {approx}3km horizontally and {approx}0.1km vertically, and 2-D speed accuracy of {le} 1m/s. This data will be used in conjunction with advanced numerical weather prediction models to provide increases in the reliability of long range weather forecasts. Major advances in data collection technology will be required to provide the proposed high-resolution data collection network. Systems studies must be undertaken to determine insertion requirements, spacing, and evolution of the transponder ensemble, which will be used to collect the data. Numerical models which provide realistic global weather pattern simulations must be utilized in order to perform these studies. A global circulation model with a 3{sup o} horizontal resolution has been used for initial simulations of the generation and evolution of transponder distributions. These studies indicate that reasonable global coverage of transponders can be achieved by a launch scenario consisting of the sequential launch of transponders at specified heights from a globally distributed set of launch sites.

  9. A high-resolution, cloud-assimilating numerical weather prediction model for solar irradiance forecasting

    E-Print Network [OSTI]

    Mathiesen, Patrick; Collier, Craig; Kleissl, Jan

    2013-01-01

    of the WRF model solar irradiance forecasts in Andalusia (Beyer, H. , 2009.    Irradiance forecasting for the power dependent probabilistic irradiance  forecasts for coastal 

  10. THE STATE OF THE ART OF NUMERICAL MODELING OF THERMOHYDROLOGIC FLOW IN FRACTURED ROCK MASSES

    E-Print Network [OSTI]

    Wang, J.S.Y.

    2013-01-01

    thermal calculations for the WIPP site in southeastern NewWaste Isolation Pilot Plant (WIPP) in bedded salt. The codepersonal communication 1980). WIPP The modeling for WIPP is

  11. A Numerical Model For The Dynamics Of Pyroclastic Flows At Galeras...

    Open Energy Info (EERE)

    a topographic barrier located more than 16 km from the source. Two initial solid volumetric fractions are modeled. For both cases, some of the structures located more distant...

  12. 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.

  13. A unified numerical model of collisional depolarization and broadening rates due to hydrogen atom collisions

    E-Print Network [OSTI]

    Derouich, M; Barklem, P S

    2015-01-01

    Interpretation of solar polarization spectra accounting for partial or complete frequency redistribution requires data on various collisional processes. Data for depolarization and polarization transfer are needed but often missing, while data for collisional broadening are usually more readily available. Recent work by Sahal-Br\\'echot and Bommier concluded that despite underlying similarities in the physics of collisional broadening and depolarization processes, relationships between them are not possible to derive purely analytically. We aim to derive accurate numerical relationships between the collisional broadening rates and the collisional depolarization and polarization transfer rates due to hydrogen atom collisions. Such relationships would enable accurate and efficient estimation of collisional data for solar applications. Using earlier results for broadening and depolarization processes based on general (i.e. not specific to a given atom), semi-classical calculations employing interaction potentials...

  14. Interstellar hydrogen fluxes measured by IBEX-Lo in 2009: numerical modeling and comparison with the data

    E-Print Network [OSTI]

    Katushkina, O A; Alexashov, D B; Schwadron, N A; McComas, D J

    2015-01-01

    In this paper, we perform numerical modeling of the interstellar hydrogen fluxes measured by IBEX-Lo during orbit 23 (spring 2009) using a state-of-the-art kinetic model of the interstellar neutral hydrogen distribution in the heliosphere. This model takes into account the temporal and heliolatitudinal variations of the solar parameters as well as non-Maxwellian kinetic properties of the hydrogen distribution due to charge exchange in the heliospheric interface. We found that there is a qualitative difference between the IBEX-Lo data and the modeling results obtained with the three-dimensional, time-dependent model. Namely, the model predicts a larger count rate in energy bin~2 (20-41 eV) than in energy bin~1 (11-21 eV), while the data shows the opposite case. We perform study of the model parameter effects on the IBEX-Lo fluxes and the ratio of fluxes in two energy channels. We shown that the most important parameter, which has a major influence on the ratio of the fluxes in the two energy bins, is the solar...

  15. / free-slip / free-slip Numerical models of slab deformation in the upper mantle

    E-Print Network [OSTI]

    Cerveny, Vlastislav

    of the plate. Model: Model domain: Results: Effect of decoupling mechanism 410 km 660 km oceánska doska nadlozná doska Results: Effect of boundary condition Effect of yield stress Plate velocity Kinematic, crust.s 1020 Pa.s 1021 Pa.s Free-slip, crust = 1020 Pas, y = 108 Pa Kinematic, crust = 1020 Pas, y = 109 Pa

  16. NUMERICAL MODELING FOR THE FORMATION MECHANISM OF 3D TOPOGRAPHY ON MICROBIAL MAT SURFACES 

    E-Print Network [OSTI]

    Patel, Harsh Jay

    2013-09-27

    a node in the fluid to solid (aggregate or the particle) ? b number of discrete velocities 14 In order to stochastically solve the stochastic solute diffusion equation (Eq. 6), initial and boundary conditions are specified before solving...) Model........................................... 1.2.2 The Diffusion Limited Aggregation-Cellular Automata (DLA-CA) Model................................................................................................ 1.2.3 The Reaction-Diffusion...

  17. Simulation of a Polar Low Case in the North Atlantic with different regional numerical models

    E-Print Network [OSTI]

    Zahn, Matthias

    by the DWD (German Weather Service) by means of their forecast model HRM (High Resolution Model) and another University Press, Cambridge. (a) CLM (b) REMO (c) HRM, DWD (d) BWK Figure 1: 1(a)- 1(c)10m wind velocity pressure from CLM and REMO simulations and HRM analysis, DWD, respectively, at 15/10/93, 6:00, 1(d) surface

  18. Controls on advance of tidewater glaciers: results from numerical modeling applied to Columbia Glacier

    E-Print Network [OSTI]

    Nick, F. M.; van der Veen, Cornelis J.; Oerlemans, J.

    2007-07-11

    at the glacier terminus and presented the flotation model. In the flotation model the terminus position is defined as the point where the ice thickness exceeds the flotation thickness by an amount H0. If the glacier thins, the terminus will retreat to a point... where this condition is again satisfied. Vieli et al. [2001] modified the flotation criterion and defined the thickness in excess of flotation H0 as a fraction of the flotation thickness. A recent modeling study by Nick and Oerlemans [2006] compared both...

  19. A simplified numerical model of coronal energy dissipation based on reduced MHD

    E-Print Network [OSTI]

    E. Buchlin; V. Aletti; S. Galtier; M. Velli; G. Einaudi; J. -C. Vial

    2003-03-27

    A 3D model intermediate between cellular automata (CA) models and the reduced magnetohydrodynamic (RMHD) equations is presented to simulate solar impulsive events generated along a coronal magnetic loop. The model consists of a set of planes distributed along a magnetic loop between which the information propagates through Alfven waves. Statistical properties in terms of power-laws for energies and durations of dissipative events are obtained, and their agreement with X-ray and UV flares observations is discussed. The existence of observational biases is also discussed.

  20. Numerical Modeling of the Transient Thermal Interference of Vertical U-Tube Haet Exchangers 

    E-Print Network [OSTI]

    Muraya, Norman K.

    1994-01-01

    Non-linear finite element models were developed to simulate transient heat and mass transfer in the soil surrounding the ground heat exchangers of ground-coupled heat pumps (GCHPs) operating in the cooling mode. Parametric ...

  1. Numerical and analytical modeling of heat transfer between fluid and fractured rocks

    E-Print Network [OSTI]

    Li, Wei, S.M. Massachusetts Institute of Technology

    2014-01-01

    Modeling of heat transfer between fluid and fractured rocks is of particular importance for energy extraction analysis in EGS, and therefore represents a critical component of EGS design and performance evaluation. In ...

  2. Numerical Modeling of Fractured Shale-Gas and Tight-Gas Reservoirs Using Unstructured Grids 

    E-Print Network [OSTI]

    Olorode, Olufemi Morounfopefoluwa

    2012-02-14

    Various models featuring horizontal wells with multiple induced fractures have been proposed to characterize flow behavior over time in tight gas and shale gas systems. Currently, there is little consensus regarding the effects of non...

  3. An Efficient Numerical Scheme for Simulating Unidirectional Irregular Waves Based on a Hybrid Wave Model 

    E-Print Network [OSTI]

    Jia, Dongxing 1984-

    2012-11-15

    The Unidirectional Hybrid Wave Model (UHWM) predicts irregular wave kinematics and pressure accurately in comparison with its linear counterpart and modification, especially near the free surface. Hence, in using the Morrison equation it has been...

  4. 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.

  5. 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.

  6. 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.

  7. 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.

  8. Modeling and numerical techniques for high-speed digital simulation of nuclear power plants

    SciTech Connect (OSTI)

    Wulff, W.; Cheng, H.S.; Mallen, A.N.

    1987-01-01

    Conventional computing methods are contrasted with newly developed high-speed and low-cost computing techniques for simulating normal and accidental transients in nuclear power plants. Six principles are formulated for cost-effective high-fidelity simulation with emphasis on modeling of transient two-phase flow coolant dynamics in nuclear reactors. Available computing architectures are characterized. It is shown that the combination of the newly developed modeling and computing principles with the use of existing special-purpose peripheral processors is capable of achieving low-cost and high-speed simulation with high-fidelity and outstanding user convenience, suitable for detailed reactor plant response analyses.

  9. 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)

  10. Model Validation and Spatial Interpolation by Combining Observations with Outputs from Numerical

    E-Print Network [OSTI]

    Washington at Seattle, University of

    ""r,c,rn The authors are for hel]JfuI #12;Abstract Constructing maps of pollution levels is vital for air quality concentrations. Key tlJords: air pollution, Ba~yesian inference, change of support, likelihood approaches, Matern Resolutions 2.5 Modeling a Nonstationary Covariance . 3 Estimation 3.1 Algorithm 4 Application: Air Pollution

  11. NUMERICAL MODEL OF PLASMA CONFINEMENT J. C. Sprott and E. J. Strait

    E-Print Network [OSTI]

    Sprott, Julien Clinton

    , a fast, simple, time-dependent, zero-dimensional computer model for plasmas confined in cylindrical chamber wall temperature (eV) B = average magnetic field strength (kG) B o = peak magnetic field strength) For the levitated toroidal oetupole, we take the values: T w = .025 eV a = 50 em L = 800 em Q = " 104 o S = 103

  12. 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.

  13. A numerical procedure to model and monitor CO2 sequestration in ...

    E-Print Network [OSTI]

    santos

    Sep 7, 2012 ... The model of the formation is based on the porosity and clay .... The porosity and dry bulk and shear modulus of the formation were obtained using a pore ... The viscosity, density and bulk modulus of CO2 needed for the flow ...

  14. Numerical modeling of spray cooling-assisted dermatologic laser surgery for treatment of port wine stains

    E-Print Network [OSTI]

    Aguilar, Guillermo

    . The objectives of this paper are: to improve the thermal modeling of skin undergoing CSC-assisted DLS for PWS that incorporate all the complexity of the problem are a valuable and fundamental research tool. Keywords: Spray Engineering, University of California, Riverside, CA 92521 bBeckman Laser Institute Medical Clinic, University

  15. A nanoscale numerical model of calcium silicate hydrate P.C. Fonseca a,

    E-Print Network [OSTI]

    Andrade, Jose

    to pre- dict the bulk properties of cement and concrete, such as shrinkage, creep, permeability and the properties of structural concrete are not fully under- stood. Models are becoming increasingly important, and cracking. C­S­H is responsible for much of the cohesive proper- ties in concrete but the chemical origin

  16. An XFEM Model for Carbon Sequestration Journal: International Journal for Numerical Methods in Engineering

    E-Print Network [OSTI]

    Gracie, Robert

    the temperatures and pressures are such that the CO2 exists in a supercritical state. The supercritical CO2 is less method (XFEM) model is presented to analyze the injection and sequestration of carbon dioxide (CO2 method (SUPG-FEM-FDM) to approximate the distribution of CO2 in the aquifer. Near well enrichment

  17. Malcolm Roberts: Research Statement My research interests are mathematical modelling, numerical analysis, and

    E-Print Network [OSTI]

    Roberts, Malcolm

    analysis, and high-performance computing. I have experience working in the areas of compu- tational fluid mechanics, subgrid models for turbulence, and flows in complex geometries. 1 Introduction Consider, which is the product of the characteristic length and velocity of the flow divided by the kinematic

  18. Numerical study of three-dimensional PIC for the surface plasmon excitation based on Drude model

    E-Print Network [OSTI]

    Liu, La-Qun; Wang, Hui-Hui; Liu, Da-Gang

    2015-01-01

    This paper explores the time-domain equations of noble metals, in which Drude model is adopted to describe the dielectric constant, to implement three-dimensional particle-in-cell (PIC) simulations for the surface plasmon excitation with the finite-difference time-domain method (FDTD). A three-dimensional model for an electron bunch movement near the metal film is constructed, and particle-in-cell (PIC) simulations are carried out with various metal films of different thicknesses. The frequency of surface plasmon obtained from PIC simulation is agreed with that from theory. Furthermore, the surface plasmon wave properties of excitation and propagation with the metal film is summarized by PIC results.

  19. Non-smooth Dynamics Using Differential-algebraic Equations Perspective: Modeling and Numerical Solutions 

    E-Print Network [OSTI]

    Gotika, Priyanka

    2012-02-14

    : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : 73 ix LIST OF FIGURES FIGURE Page 1 Constitutive relation for Bingham-type uid. Force due to dash- pot cannot be written as a function of kinematic variables (that is, velocity of the mass). On the other hand, as shown in the gure, velocity... of constitutive relation between velocity and force due to dash-pot for a dash-pot governed by visco-elastic model. : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : 9 5 A pictorial representation of constitutive relation between velocity...

  20. Modeling electrospinning process and a numerical scheme using Lattice Boltzmann method to simulate viscoelastic fluid flows 

    E-Print Network [OSTI]

    Karra, Satish

    2009-05-15

    [35] developed a stability analysis for the ax- isymmetric Rayleigh instability as well as the non-axisymmetric whipping instability in Newtonian fluids. According to them, the Rayleigh instability due to electrical forces is equivalent to the surface... for viscoelasticity between the beads. . . . . . . . 19 7 Typical result for the discrete particle model showing the bending loop in the jet. Top view shows that the envelope of the jet trajectory is a cone. The number of beads for this simulation N = 100 and non...

  1. Numerical Model for the Deformation of Nucleated Cells by Optical Stretchers

    E-Print Network [OSTI]

    Ihab Sraj; Joshua Francois; David W. M. Marr; Charles D. Eggleton

    2015-01-15

    In this paper, we seek to model the deformation of nucleated cells by single diode-laser bar optical stretchers. We employ a recently developed computational model, the Dynamic Ray-Tracing method, to determine the stress distribution induced by the applied optical forces on a capsule encapsulating a nucleus of different optical properties. These forces are shape dependent and can deform real non-rigid objects; thus resulting in a dynamically changing optical stress distribution with cell and nucleus deformation. Chinese hamster ovary cell is a common biological cell that is of interest to the biomedical community because of their use in recombinant protein therapeutics and is an example of a nucleated cell. To this end, we model chinese hamster ovary cells as two three-dimensional elastic capsules of variable inner capsule size immersed in a fluid where the hydrodynamic forces are calculated using the Immersed Boundary Method. Our results show that the presence of a nucleus has a major effect on the force distribution on the cell surface and the net deformation. Scattering and gradient forces are reported for different nucleus sizes and the effect of nucleus size on the cell deformation is discussed.

  2. Numerical Model for the Deformation of Nucleated Cells by Optical Stretchers

    E-Print Network [OSTI]

    Sraj, Ihab; Marr, David W M; Eggleton, Charles D

    2015-01-01

    In this paper, we seek to model the deformation of nucleated cells by single diode-laser bar optical stretchers. We employ a recently developed computational model, the Dynamic Ray-Tracing method, to determine the stress distribution induced by the applied optical forces on a capsule encapsulating a nucleus of different optical properties. These forces are shape dependent and can deform real non-rigid objects; thus resulting in a dynamically changing optical stress distribution with cell and nucleus deformation. Chinese hamster ovary cell is a common biological cell that is of interest to the biomedical community because of their use in recombinant protein therapeutics and is an example of a nucleated cell. To this end, we model chinese hamster ovary cells as two three-dimensional elastic capsules of variable inner capsule size immersed in a fluid where the hydrodynamic forces are calculated using the Immersed Boundary Method. Our results show that the presence of a nucleus has a major effect on the force dis...

  3. Modelling the high-energy emission from gamma-ray binaries using numerical relativistic hydrodynamics

    E-Print Network [OSTI]

    Dubus, Guillaume; Fromang, Sébastien

    2015-01-01

    Detailed modeling of the high-energy emission from gamma-ray binaries has been propounded as a path to pulsar wind physics. Fulfilling this ambition requires a coherent model of the flow and its emission in the region where the pulsar wind interacts with the stellar wind of its companion. We developed a code that follows the evolution and emission of electrons in the shocked pulsar wind based on inputs from a relativistic hydrodynamical simulation. The code is used to model the well-documented spectral energy distribution and orbital modulations from LS 5039. The pulsar wind is fully confined by a bow shock and a back shock. The particles are distributed into a narrow Maxwellian, emitting mostly GeV photons, and a power law radiating very efficiently over a broad energy range from X-rays to TeV gamma rays. Most of the emission arises from the apex of the bow shock. Doppler boosting shapes the X-ray and VHE lightcurves, constraining the system inclination to $i\\approx 35^{\\rm o}$. There is a tension between th...

  4. 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.

  5. 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

    In this work, it has been shown that, for the given sets of parameters (transport coefficients), the Tangent-Predictor (TP) continuation method, which was used in the coarsest grid, works remarkably well. The problems in finding an initial guess that resides well within Newton`s method radius of convergence are alleviated by correcting the initial guess by the predictor step of the TP method. The TP method works well also in neutral gas puffing and impurity simulations. The neutral gas puffing simulation is performed by systematically increasing the fraction of puffing rate according to the TP method until it reaches a desired condition. Similarly, the impurity simulation characterized by using the fraction of impurity density as the continuation parameter, is carried out in line with the TP method. Both methods show, as expected, a better performance than the classical embedding (CE) method. The convergence criteria {epsilon} is set to be 10{sup {minus}9} based on the fact that lower value of {epsilon} does not alter the solution significantly. Correspondingly, the number of Newton`s iterations in the corrector step of the TP method decrease substantially, an extra point in terms of code speed. The success of the TP method enlarges the possibility of including other sets of parameters (operations and physics). With the availability of the converged coarsest grid solution, the next forward step to the multigrid cycle becomes possible. The multigrid method shows that the memory storage problems that plagued the application of Newton`s method on fine grids, are of no concern. An important result that needs to be noted here is the performance of the FFCD model. The FFCD model is relatively simple and is based on the overall results the model has shown to predict different divertor plasma parameters. The FFCD model treats exactly the implementation of the deep penetration of energetic neutrals emerging from the divertor plate. The resulting ionization profiles are 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.

  6. Mathematical and numerical models to achieve high speed with special-purpose parallel processors

    SciTech Connect (OSTI)

    Cheng, H.S.; Wulff, W.; Mallen, A.N.

    1986-07-01

    One simulation facility that has been developed is the BNL Plant Analyzer, currently set up for BWR plant simulations at up to seven times faster than real-time process speeds. The principal hardware components of the BNL Plant Analyzer are two units of special-purpose parallel processors, the AD10 of Applied Dynamics International and a PDP-11/34 host computer. The AD10 is specifically designed for time-critical system simulations, utilizing the modern parallel processing technology with pipeline architecture. The simulator employs advanced modeling techniques and efficient integration techniques in conjunction with the parallel processors to achieve high speed performance.

  7. Numerical modeling of boundary layer flow under shoaling and breaking waves 

    E-Print Network [OSTI]

    Pattipawaej, Olga Catherina

    1998-01-01

    velocities and shear stresses in the bottom boundary layer under breaking waves. The experiment wss conducted in the Precision Wave Tank located in the Ocean Engineering Laboratory at the University of Delaware. The flume was 33 m long, 0. 6 m wide, and 1... modeling of surf zone hydrodynamics. " Rep. No. CACR-95-97, Center for Applied Coastal Research, University Delaware, Newark, Delaware. Cox, D. T. , Kobayashi, N. , and Okayasu, A. (1996). "Bottom shear stress in the surf zone. " J. Geophys. Res. , 101(C...

  8. Numerical modeling of two-phase behavior in the PEFC gas diffusion layer

    SciTech Connect (OSTI)

    Mukherjee, Partha Pa223876 [Los Alamos National Laboratory; Kang, Qinjun [Los Alamos National Laboratory; Mukundan, Rangachary [Los Alamos National Laboratory; Borup, Rod L [Los Alamos National Laboratory

    2009-01-01

    A critical performance limitation in the polymer electrolye fuel cell (PEFC) is attributed to the mass transport loss originating from suboptimal liquid water transport and flooding phenomena. Liquid water can block the porous pathways in the fibrous gas diffusion layer (GDL) and the catalyst layer (CL), thus hindering oxygen transport from the flow field to the electrochemically actives sites in the catalyst layer. In this paper, the study of the two phase behavior and the durability implications due to the wetting characteristics in the carbon paper GDL are presented using a pore-scale modeling framework.

  9. 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 laser–plasma 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 laser–plasma accelerator.

  10. Numerical study of a binary Yukawa model in regimes characteristic of globular proteins in solutions

    SciTech Connect (OSTI)

    Giacometti, Achille; Gazzillo, Domenico; Pastore, Giorgio; Das, Tushar Kanti [Istituto Nazionale per la Fisica della Materia and Dipartimento di Chimica Fisica, Universita di Venezia, S. Marta DD 2137, I-30123 Venice (Italy); Dipartimento di Fisica Teorica, Universita di Trieste and INFM-DEMOCRITOS, National Simulation Center, Strada Costiera 11, Miramare, I-34100 Trieste (Italy); ICTP, Diploma Course, Strada Costiera 11, Miramare P.O Box 586, I-34100 Trieste (Italy)

    2005-03-01

    The main goal of this paper is to assess the limits of validity, in the regime of low concentration and strong Coulomb coupling (high molecular charges), of a simple perturbative approximation to the radial distribution functions (RDF's), based upon a low-density expansion of the potential of mean force and proposed to describe protein-protein interactions in a recent small-angle-scattering (SAS) experimental study. A highly simplified Yukawa (screened Coulomb) model of monomers and dimers of a charged globular protein ({beta}-lactoglobulin) in solution is considered. We test the accuracy of the RDF approximation, as a necessary complementary part of the previous experimental investigation, by comparison with the fluid structure predicted by approximate integral equations and exact Monte Carlo (MC) simulations. In the MC calculations, an Ewald construction for Yukawa potentials has been used to take into account the long-range part of the interactions in the weakly screened cases. Our results confirm that the perturbative first-order approximation is valid for this system even at strong Coulomb coupling, provided that the screening is not too weak (i.e., for Debye length smaller than monomer radius). A comparison of the MC results with integral equation calculations shows that both the hypernetted-chain (HNC) and Percus-Yevick closures have a satisfactory behavior under these regimes, with the HNC being superior throughout. The relevance of our findings for interpreting SAS results is also discussed.

  11. Fast Track Analysis of Shale Numerical Models A. Kalantari-Dahaghi ,SPE, S. Esmaili, SPE, West Virginia University, S.D. Mohaghegh, SPE, Intelligent Solution

    E-Print Network [OSTI]

    Mohaghegh, Shahab

    SPE 162699 Fast Track Analysis of Shale Numerical Models A. Kalantari-Dahaghi ,SPE, S. Esmaili, SPE of SPE copyright. Abstract Latest advances in shale gas reservoir simulation and modeling have made it possible to optimize and enhance the production from organic rich shale gas reservoirs. Reservoir simulator

  12. Motivation Literature Trans. Network Models Relative Total Cost Indices Numerical Examples Summary A Relative Total Cost Index for the Evaluation of

    E-Print Network [OSTI]

    Nagurney, Anna

    A Relative Total Cost Index for the Evaluation of Transportation Network Robustness in the Presence Index #12;Motivation Literature Trans. Network Models Relative Total Cost Indices Numerical Examples, Qiang Qiang Relative Total Cost Index #12;Motivation Literature Trans. Network Models Relative Total

  13. 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.

  14. 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.

  15. IEEE TRANSACTIONS ON COMPUTER-AIDED DESIGN OF INTEGRATED CIRCUITS AND SYSTEMS, VOL. 17, NO. 4, APRIL 1998 305 Numerically Stable Green Function for Modeling and

    E-Print Network [OSTI]

    California at Berkeley, University of

    IEEE TRANSACTIONS ON COMPUTER-AIDED DESIGN OF INTEGRATED CIRCUITS AND SYSTEMS, VOL. 17, NO. 4, APRIL 1998 305 Numerically Stable Green Function for Modeling and Analysis of Substrate Coupling. Meyer, Fellow, IEEE Abstract-- The Green function over a multilayer substrate is derived by solving

  16. Numerical analysis of a locking-free mixed finite element method for a bending moment formulation of Reissner-Mindlin plate model

    E-Print Network [OSTI]

    Rodríguez, Rodolfo

    ´on, Chile. This paper deals with the approximation of the bending of a clamped plate modeled by ReissnerNumerical analysis of a locking-free mixed finite element method for a bending moment formulation the thickness t is small. Here, we propose a new mixed formulation in terms of the bending moments, shear stress

  17. An analytical and numerical model to determine stresses in a Rock Melt Drill produced glass liner for potential use on Mars 

    E-Print Network [OSTI]

    McConnell, Joshua B

    2000-01-01

    A numerical and analytical model was constructed to determine the resultant stress state imposed on a Rock Melt Drill produced liner. The purpose of this study was to determine if the liner produced would possess the strength required to prevent a...

  18. Numerical Transform Inversion to Analyze Teletraffic Models Gagan L. Choudhury, a David M. Lucantoni a and Ward Whitt b

    E-Print Network [OSTI]

    Whitt, Ward

    . Lucantoni a and Ward Whitt b a AT&T Bell Laboratories, Holmdel, NJ 07733­3030, USA b AT&T Bell Laboratories, Murray Hill, NJ 07974­0636, USA We describe recently developed algorithms for numerically inverting

  19. Numerical failure analysis of a stretch-bending test on dual-phase steel sheets using a phenomenological fracture model

    E-Print Network [OSTI]

    Luo, Meng

    Advanced High Strength Steels (AHSS) are increasingly used in automotive industry due to their superior strength and substantial weight advantage. However, their compromised ductility gives rise to numerous manufacturing ...

  20. 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.

  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. 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

    This three-year project, in cooperation with Professor Bob Houze at University of Washington, has been successfully finished as planned. Both ARM (the Atmospheric Radiation Measurement Program) data and cloud-resolving model (CRM) simulations were used to identify the water budgets of clouds observed in two international field campaigns. The research results achieved shed light on several key processes of clouds in climate change (or general circulation models), which are summarized below. 1. Revealed the effect of mineral dust on mesoscale convective systems (MCSs) Two international field campaigns near a desert and a tropical coast provided unique data to drive and evaluate CRM simulations, which are TWP-ICE (the Tropical Warm Pool International Cloud Experiment) and AMMA (the African Monsoon Multidisciplinary Analysis). Studies of the two campaign data were contrasted, revealing that much mineral dust can bring about large MCSs via ice nucleation and clouds. This result was reported as a PI presentation in the 3rd ASR Science Team meeting held in Arlington, Virginia in March 2012. A paper on the studies was published in the Journal of the Atmospheric Sciences (Zeng et al. 2013). 2. Identified the effect of convective downdrafts on ice crystal concentration Using the large-scale forcing data from TWP-ICE, ARM-SGP (the Southern Great Plains) and other field campaigns, Goddard CRM simulations were carried out in comparison with radar and satellite observations. The comparison between model and observations revealed that convective downdrafts could increase ice crystal concentration by up to three or four orders, which is a key to quantitatively represent the indirect effects of ice nuclei, a kind of aerosol, on clouds and radiation in the Tropics. This result was published in the Journal of the Atmospheric Sciences (Zeng et al. 2011) and summarized in the DOE/ASR Research Highlights Summaries (see http://www.arm.gov/science/highlights/RMjY5/view). 3. Used radar 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 were compared with three reanalyses (MERRA, ERA-Interim and CFSR). Although the MMF tends

  3. Multiscale Constitutive Modeling and Numerical Simulations of the Thermomechanical Response of Polycrystalline NiTi Shape Memory Alloy

    E-Print Network [OSTI]

    Sengupta, Arkaprabha

    2010-01-01

    Phase transformation in Nitinol polycrystals 3.1 Martensiticvector pairs (m ? , b ? ) for Nitinol . . . . . . . . . . .4 Constitutive model of Nitinol 4.1 Constitutive models for

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

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

    foam infiltrated with a phase change material (PCM). The project will extend the graphite foamPCM LHTES system to make it compatible with supercritical CO2 power cycle...

  5. Concentrating Solar Power: Advanced Projects Offering Low LCOE

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

    AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Homesum_a_epg0_fpd_mmcf_m.xls" ,"Available from WebQuantity of Natural GasAdjustmentsShirleyEnergy A plug-in electric vehicle (PEV)Day-June 22,FresnoSky EnergyOpportunities (CSP:

  6. 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.

  7. Florian SEITZ: Atmospheric and oceanic impacts to Earth rotations numerical studies with a dynamic Earth system model

    E-Print Network [OSTI]

    Schuh, Harald

    with a dynamic Earth system model (completed in October 2004) Variations of Earth rotation are caused Earth system model DyMEG has been developed. It is based on the balance of angular momentum

  8. A Study to Verify the Material Surface Concept of Water Table by Examining Analytical and Numerical Models

    E-Print Network [OSTI]

    Dadi, Sireesh Kumar

    2011-10-21

    : the Neuman model, which assumes instantaneous drainage from the unsaturated zone; the Moench model, which considered gradual drainage from the unsaturated zone using a series of exponential terms in the water table boundary condition; and the Mathias...

  9. Atmospheric Test Models and Numerical Experiments for the Simulation of the Global Distribution of Weather Data Transponders II. Vertical Transponder Motion Considerations

    SciTech Connect (OSTI)

    Grossman, A.; Errico, R.M.

    1999-11-29

    The vertical motion of constant density atmospheric balloons has been considered via an equation of motion for the vertical displacement of a balloon, due to vertical air motion, which can be numerically solved for balloon positions. Initial calculations are made for a constant density atmosphere. Various vertical wind models with relatively large amplitudes are applied to the model to determine how tightly the balloons are coupled to the reference level and the time scale for the balloons to change to the wind driven reference altitude. A surface launch of a balloon to a 6 km reference altitude is modeled using a detailed atmospheric pressure-density-temperature profile in the equation of motion. The results show the balloons to be relatively tightly coupled ({approx} 50-100 m) to the reference altitude.

  10. Numerical Simulation of PulseTube Refrigerators: 1D model I.A. Lyulina 1 , R.M.M. Mattheij 1 , A.S. Tijsseling 1 , A.T.A.M. de Waele 2

    E-Print Network [OSTI]

    Eindhoven, Technische Universiteit

    of a piston, an aftercooler (AC), a regenerator, a cold heat exchanger (CHX), a tube, a hot heat exchanger ( numerical model has been introduced to study steady oscillatory heat and mass transfer in the tube section, numerical simulation, high resolution scheme 1 Introduction The pulse tube is a relatively new type

  11. Numerical analysis of a finite element method for the axisymmetric eddy current model of an induction furnace

    E-Print Network [OSTI]

    Rodríguez, Rodolfo

    of an induction furnace ALFREDO BERM ´UDEZ Departamento de Matem´atica Aplicada, Universidade de Santiago de problem arising from the modeling of an induction furnace. By taking advantage of the cylindrical symmetry a formulation of the eddy current problem arising from the modeling of an induction furnace, which is based

  12. Numerical modeling of time-lapse seismic data from fractured reservoirs including fluid flow and geochemical processes 

    E-Print Network [OSTI]

    Shekhar, Ravi

    2009-05-15

    and amplitude variation with offset (AVO) results for our example model predicts that CO2 is easier to detect than brine in the fractured reservoirs. The effects of geochemical processes on seismics are simulated by time-lapse modeling for t = 1000 years. My...

  13. Comparison of numerical, analytical, and empirical models to estimate wetting patterns for surface and subsurface drip irrigation

    E-Print Network [OSTI]

    Kandelous, Maziar M.; Šim?nek, Ji?í

    2010-01-01

    and clay (33.5%), and water contents for pressures of -33times, antecedent water contents, and emitter locations, andthe saturated soil water content as one of the model

  14. 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.

  15. Development of analytical and numerical models predicting the deposition rate of electrically charged particles in turbulent channel flows 

    E-Print Network [OSTI]

    Ko, Hanseo

    1994-01-01

    An analytical model is established to predict an electrostatically charged particle deposition as a function of particle size in fully-developed turbulent pipe flow. The convectivediffusion flux equation is solved for the particle concentration as a...

  16. Numerical Modeling Studies of The Dissolution-Diffusion-Convection Process During CO2 Storage in Saline Aquifers

    E-Print Network [OSTI]

    Pruess, Karsten

    2008-01-01

    Properties of the Supercritical CO2–Water Pure Interface, J.and pressure conditions, supercritical CO2 (scCO2) has lowerthe overlying free supercritical CO2 phase. This is modeled

  17. Experimental measurements and numerical modeling of fast-ion distributions in the Alcator C-Mod Tokamak

    E-Print Network [OSTI]

    Bader, Aaron Craig

    2012-01-01

    In this thesis we discuss measurements and modeling of minority heated fast-ion distributions in the Ion Cyclotron Range of Frequencies (ICRF) on the Alcator C-Mod tokamak. Analysis of fast-ions >100Te is important for ...

  18. PoPe (Projection on Proper elements) for code control: verification, numerical convergence and reduced models. Application to plasma turbulence simulations

    E-Print Network [OSTI]

    Cartier-Michaud, T; Sarazin, Y; Abiteboul, J; Bufferand, H; Dif-Pradalier, G; Garbet, X; Grandgirard, V; Latu, G; Norscini, C; Passeron, C; Tamain, P

    2015-01-01

    The Projection on Proper elements (PoPe) is a novel method of code control dedicated to 1) checking the correct implementation of models, 2) determining the convergence of numerical methods and 3) characterizing the residual errors of any given solution at very low cost. The basic idea is to establish a bijection between a simulation and a set of equations that generate it. Recovering equations is direct and relies on a statistical measure of the weight of the various operators. This method can be used in any dimensions and any regime, including chaotic ones. This method also provides a procedure to design reduced models and quantify the ratio costs to benefits. PoPe is applied to a kinetic and a fluid code of plasma turbulence.

  19. Numerical modeling of mixed sediment resuspension, transport, and deposition during the March 1998 episodic events in southern Lake Michigan

    SciTech Connect (OSTI)

    Lee, Cheegwan; Schwab, David J.; Beletsky, Dmitry; Stroud, Jonathan; Lesht, B. M.

    2007-02-17

    A two-dimensional sediment transport model capable of simulating sediment resuspension of mixed (cohesive+noncohesive) sediment is developed and applied to quantitatively simulate the March 1998 resuspension events in southern Lake Michigan. Some characteristics of the model are the capability to incorporate several floc size classes, a physically-based settling velocity formula, bed armoring, and sediment availability limitation. Important resuspension parameters were estimated from field and laboratory measurement data. The model reproduced the resuspension plume (observed by the SeaWIFS satellite and field instruments) and recently measured sedimentation rate distribution (using radiotracer techniques) fairly well. Model results were verified with field measurements of suspended sediment concentration and settling flux (by ADCPs and sediment traps). Both wave conditions and sediment bed properties (critical shear stress, fine sediment fraction, and limited sediment availability or source) are the critical factors that determine the concentration distribution and width of the resuspension plume. The modeled sedimentation pattern shows preferential accumulation of sediment on the eastern side of the lake, which agrees with the observed sedimentation pattern despite a predominance of particle sources from the western shoreline. The main physical mechanisms determining the sedimentation pattern are 1) the two counter-rotating circulation gyres producing offshore mass transport along the southeastern coast during northerly wind and 2) the settling velocity of sediment flocs which controls the deposition location.

  20. Numerical modeling of mixed sediment resuspension, transport, and deposition during the March 1998 episodic events in southern Lake Michigan.

    SciTech Connect (OSTI)

    Lee, C.; Schwab, D. J.; Beletsky, D.; Stroud, J.; Lesht, B.; PNNL; NOAA; Univ. of Michigan; Univ. of Pennsylvania

    2007-02-17

    A two-dimensional sediment transport model capable of simulating sediment resuspension of mixed (cohesive plus noncohesive) sediment is developed and applied to quantitatively simulate the March 1998 resuspension events in southern Lake Michigan. Some characteristics of the model are the capability to incorporate several floc size classes, a physically based settling velocity formula, bed armoring, and sediment availability limitation. Important resuspension parameters were estimated from field and laboratory measurement data. The model reproduced the resuspension plume (observed by the SeaWIFS satellite and field instruments) and recently measured sedimentation rate distribution (using radiotracer techniques) fairly well. Model results were verified with field measurements of suspended sediment concentration and settling flux (by ADCPs and sediment traps). Both wave conditions and sediment bed properties (critical shear stress, fine sediment fraction, and limited sediment availability or source) are the critical factors that determine the concentration distribution and width of the resuspension plume. The modeled sedimentation pattern shows preferential accumulation of sediment on the eastern side of the lake, which agrees with the observed sedimentation pattern despite a predominance of particle sources from the western shoreline. The main physical mechanisms determining the sedimentation pattern are (1) the two counter-rotating circulation gyres producing offshore mass transport along the southeastern coast during northerly wind and (2) the settling velocity of sediment flocs which controls the deposition location.

  1. Modelling and numerical approximation of a 2.5D set of equations for mesoscale atmospheric processes

    E-Print Network [OSTI]

    Kalise, Dante

    2011-01-01

    The set of 3D inviscid primitive equations for the atmosphere is dimensionally reduced by a Discontinuous Galerkin discretization in one horizontal direction. The resulting model is a 2D system of balance laws where with a source term depending on the layering procedure and the choice of coupling fluxes, which is established in terms of upwind considerations. The "2.5D" system is discretized via a WENO-TVD scheme based in a flux limiter centered approach. We study four tests cases related to atmospheric phenomena to analyze the physical validity of the model.

  2. Three-dimensional numerical models with varied material properties and erosion rates: Implications for the mechanics and kinematics

    E-Print Network [OSTI]

    Mueller, Karl

    for the mechanics and kinematics of compressive wedges Phaedra Upton,1,2 Karl Mueller,3 and Yue-Gau Chen4 Received develop three-dimensional mechanical models of a compressive wedge and investigate how the form and kinematics of the outboard wedge are affected by variation in initial topography, material properties

  3. Preprint of the paper "A Numerical Refraction-Diffraction Model in the Ray-front Coordinate System"

    E-Print Network [OSTI]

    Colominas, Ignasi

    theory was adapted to water waves3 , allowing engineers and scientists to assess the diffractive effects" (CD-ROM), Parte IV: "Fluid Mechanics", Sección 6: "Shallow Water and Water Waves". S.R. Idelsohn, E, Ray Methods. Abstract. A model for the propagation of water waves is presented, which retrieves

  4. Numerical studies of a one-dimensional three-spin spin-glass model with long-range interactions 

    E-Print Network [OSTI]

    Larson, Derek; Katzgraber, Helmut G.; Moore, M. A.; Young, A. P.

    2010-01-01

    are presented in Sec. V and our conclusions are summarized in Sec. VI. II. THEORETICAL BACKGROUND The field theory associated with p-spin models is a cubic field theory1?3 with the following Ginzburg-Landau-Wilson Hamiltonian: HGLW =? ddr? t2 ??? q??2 ?r...

  5. A TWO PRESSURE NUMERICAL MODEL OF TWO FLUID MIXING JAMES GLIMM, HYEONSEONG JIN, MARC LAFOREST, FOLKERT TANGERMAN,

    E-Print Network [OSTI]

    Zhang, Yongmin

    . The model is designed for the study of acceleration driven mixing layers in a chunk mix regime dominated, as in meteorological temperature inversion driven updrafts and the geological evolution of salt dome formations approxi- mately as a function of the averaged flow quantities directly, to close the system of equations

  6. NUMERICAL ANALYSIS: This refers to the analysis of mathematical problems by numerical means, es-

    E-Print Network [OSTI]

    Atkinson, Kendall

    NUMERICAL ANALYSIS: This refers to the analysis of mathematical problems by numerical means, es- pecially mathematical problems arising from models based on calculus. Effective numerical analysis requires such as rootfinding and numerical integration; but we will also look at the structure of computers and the impli

  7. Numerical Modeling Studies of The Dissolution-Diffusion-Convection ProcessDuring CO2 Storage in Saline Aquifers

    SciTech Connect (OSTI)

    Pruess, Karsten; Zhang, Keni

    2008-11-17

    For purposes of geologic storage, CO2 would be injected into saline formations at supercritical temperature and pressure conditions, and would form a separate phase that is immiscible with the aqueous phase (brine). At typical subsurface temperature and pressure conditions, supercritical CO2 (scCO2) has lower density than the aqueous phase and would experience an upward buoyancy force. Accordingly, the CO2 is expected to accumulate beneath the caprock at the top of the permeable interval, and could escape from the storage formation wherever (sub-)vertical pathways are available, such as fractures or faults through the caprock, or improperly abandoned wells. Over time, an increasing fraction of CO2 may dissolve in the aqueous phase, and eventually some of the aqueous CO2 may react with rock minerals to form poorly soluble carbonates. Dissolution into the aqueous phase and eventual sequestration as carbonates are highly desirable processes as they would increase permanence and security of storage. Dissolution of CO2 will establish phase equilibrium locally between the overlying CO2 plume and the aqueous phase beneath. If the aqueous phase were immobile, CO2 dissolution would be limited by the rate at which molecular diffusion can remove dissolved CO2 from the interface between CO2-rich and aqueous phases. This is a slow process. However, dissolution of CO2 is accompanied by a small increase in the density of the aqueous phase, creating a negative buoyancy force that can give rise to downward convection of CO2-rich brine, which in turn can greatly accelerate CO2 dissolution. This study explores the process of dissolution-diffusion-convection (DDC), using high-resolution numerical simulation. We find that geometric features of convection patterns are very sensitive to small changes in problem specifications, reflecting self-enhancing feedbacks and the chaotic nature of the process. Total CO2 dissolution rates on the other hand are found to be quite robust against modest changes in problem parameters, and are essentially constant as long as no dissolved CO2 reaches the lower boundary of the system.

  8. Quasi-periodic oscillations in accreting magnetic white dwarfs II. The asset of numerical modelling for interpreting observations

    E-Print Network [OSTI]

    Busschaert, C; Michaut, C; Bonnet-Bidaud, J -M; Mouchet, M

    2015-01-01

    Magnetic cataclysmic variables are close binary systems containing a strongly magnetized white dwarf that accretes matter coming from an M-dwarf companion. High-energy radiation coming from those objects is emitted from the accretion column close to the white dwarf photosphere at the impact region. Its properties depend on the characteristics of the white dwarf and an accurate accretion column model allows the properties of the binary system to be inferred, such as the white dwarf mass, its magnetic field, and the accretion rate. We study the temporal and spectral behaviour of the accretion region and use the tools we developed to accurately connect the simulation results to the X-ray and optical astronomical observations. The radiation hydrodynamics code Hades was adapted to simulate this specific accretion phenomena. Classical approaches were used to model the radiative losses of the two main radiative processes: bremsstrahlung and cyclotron. The oscillation frequencies and amplitudes in the X-ray and optic...

  9. Numerical examination of plasmoid-induced reconnection model for solar flares: the relation between plasmoid velocity and reconnection rate

    E-Print Network [OSTI]

    Keisuke Nishida; Masaki Shimizu; Daikou Shiota; Hiroyuki Takasaki; Tetsuya Magara; Kazunari Shibata

    2008-09-04

    The plasmoid-induced-reconnection model explaining solar flares based on bursty reconnection produced by an ejecting plasmoid suggests a possible relation between the ejection velocity of a plasmoid and the rate of magnetic reconnection. In this study, we focus on the quantitative description of this relation. We performed magnetohydrodynamic (MHD) simulations of solar flares by changing the values of resistivity and the plasmoid velocity. The plasmoid velocity has been changed by applying an additional force to the plasmoid to see how the plasmoid velocity affects the reconnection rate. An important result is that the reconnection rate has a positive correlation with the plasmoid velocity, which is consistent with the plasmoid-induced-reconnection model for solar flares. We also discuss an observational result supporting this positive correlation.

  10. Numerical Modeling of AGN Jets: Formation of Magnetically Dominated Lobes and Stability Properties of Current-carrying Jets

    E-Print Network [OSTI]

    Masanori Nakamura; Hui Li; Steven Diehl; Shengtai Li

    2007-07-10

    We argue the behavior of Poynting flux-dominated outflows from AGN in the galactic cluster systems by performing three-dimensional MHD simulations within the framework of the "magnetic tower" model. Of particular interests are the structure of MHD waves, the cylindrical radial force balance, the (de)collimation, and the stability properties of magnetic tower jets. Transition between the jet/lobe and the formation of wiggling jet by growing current-driven instability are discussed.

  11. 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...

  12. Numerical modelling of VLF radio wave propagation through earth-ionosphere waveguide and its application to sudden ionospheric disturbances

    E-Print Network [OSTI]

    Pal, Sujay

    2015-01-01

    In this thesis, we theoretically predict the normal characteristics of Very Low Frequency (3~30 kHz) radio wave propagation through Earth-ionosphere waveguide corresponding to normal behavior of the D-region ionosphere. We took the VLF narrow band data from the receivers of Indian Centre for Space Physics (ICSP) to validate our models. Detection of sudden ionospheric disturbances (SIDs) are common to all the measurements. We apply our theoretical models to infer the D-region characteristics and to reproduce the observed VLF signal behavior corresponding to such SIDs. We develop a code based on ray theory to simulate the diurnal behavior of VLF signals over short propagation paths (2000~3000 km). The diurnal variation from this code are comparable to the variation obtained from a more general Long Wave Propagation Capability (LWPC) code which is based on mode theory approach. We simulate the observational results obtained during the Total Solar Eclipse of July 22, 2009 in India. We also report and simulate a h...

  13. Numerical study on convection diffusion for gasification agent in underground coal gasification. Part I: establishment of mathematical models and solving method

    SciTech Connect (OSTI)

    Yang, L.H.; Ding, Y.M. [China University of Mining & Technology, Xuzhou (China). College of Resources and Geoscience

    2009-07-01

    The aim of this article is to discuss the distribution law of the gasification agent concentration in a deep-going way during underground coal gasification and the new method of solving the problem for the convection diffusion of the gas. In this paper, the basic features of convection diffusion for the gas produced in underground coal gasification are studied. On the basis of the model experiment, through the analysis of the distribution and patterns of variation for the fluid concentration field in the process of the combustion and gasification of the coal seams within the gasifier, the 3-D non-linear unstable mathematical models on the convection diffusion for oxygen are established. In order to curb such pseudo-physical effects as numerical oscillation and surfeit which frequently occurred in the solution of the complex mathematical models, the novel finite unit algorithm, the upstream weighted multi-cell balance method is advanced in this article, and its main derivation process is introduced.

  14. 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.

  15. Advanced Numerical Weather Prediction Techniques for Solar Irradiance Forecasting : : Statistical, Data-Assimilation, and Ensemble Forecasting

    E-Print Network [OSTI]

    Mathiesen, Patrick James

    2013-01-01

    of numerical weather prediction solar irradiance forecastsof numerical weather prediction solar irradiance forecastsnumerical weather prediction model for solar irradiance

  16. FiltrationinPorousMediaandIndustrialApplications(Cetraro,Italy1998),9-77,LectureNotesinMath.,1734,Springer,Berlin,2000 NUMERICAL SOLUTION OF RESERVOIR FLOW MODELS BASED ON

    E-Print Network [OSTI]

    Karlsen, Kenneth Hvistendahl

    ,Springer,Berlin,2000 NUMERICAL SOLUTION OF RESERVOIR FLOW MODELS BASED ON LARGE TIME STEP OPERATOR for the saturation equation, which is a convection dominated, degenerate convection-diffusion equation. Both theory and applications are discussed. The general background for the reservoir flow model is reviewed, and the main

  17. Cost and Performance Assumptions for Modeling Electricity Generation Technologies

    SciTech Connect (OSTI)

    Tidball, R.; Bluestein, J.; Rodriguez, N.; Knoke, S.

    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.

  18. Gas permeation carbon capture --- Process modeling and optimization

    SciTech Connect (OSTI)

    Morinelly, Juan; Miller, David

    2011-01-01

    A multi-staged gas permeation carbon capture process model was developed in Aspen Custom Modeler{reg_sign} (ACM) and optimized in the context of the retrofit of a 550 MW subcritical pulverized coal (PC) power plant. The gas permeation stages in the process are described by a custom multi-component, hollowfiber membrane model. Gas transport across the asymmetric membrane was modeled according to the solution-diffusion model for the selective skin layer and the assumption of negligible flux resistance by the porous support. Counter-current, one-dimensional plug flow was assumed with permeate pressure drop in the fiber lumen side due to capillary constrained flow. A modular optimization framework was used to minimize the levelized cost of electricity (LCOE) by optimizing a set of key process variables. The framework allows the external control of multiple simulation modules from different software packages from a common interface.

  19. Numerical Modeling | Open Energy Information

    Open Energy Info (EERE)

    AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on QA:QA J-E-1 SECTION J APPENDIXsourceII Jump to:Information 3rd|Northfork ElectricName Nueces

  20. Numerical Simulation of Free Standing Hybrid Risers 

    E-Print Network [OSTI]

    Hou, Tiancong

    2014-08-13

    and offloading (FPSO). The numerical scheme is mainly based on an existing in-house numerical code, known as COUPLE. In using COUPLE, URA is modeled as a rigid body. The URA connects a vertical steel riser and a flexible jumper, both of which are modeled by beam...

  1. High performance computing and numerical Volker Springel

    E-Print Network [OSTI]

    Masci, Frank

    High performance computing and numerical modelling Volker Springel Lecture Notes 43rd Saas-mail: volker.springel@h-its.org 1 arXiv:1412.5187v1[astro-ph.GA]16Dec2014 #12;#12;Contents High performance computing and numerical modelling . . . . . . . . . . . . . . . 1 Volker Springel 1 Preamble

  2. Numerical wind tunnels

    E-Print Network [OSTI]

    Souza, Paulo Victor Santos

    2015-01-01

    Flow of viscous fluids are not usually discussed in detail in general and basic courses of physics. This is due in part to the fact that the Navier-Stokes equation has analytical solution only for a few restricted cases, while more sophisticated problems can only be solved by numerical methods. In this text, we present a computer simulation of wind tunnel, i.e., we present a set of programs to solve the Navier-Stokes equation for an arbitrary object inserted in a wind tunnel. The tunnel enables us to visualize the formation of vortices behind object, the so-called von K\\'arm\\'an vortices, and calculate the drag force on the object. We believe that this numerical wind tunnel can support the teacher and allow a more elaborate discussion of viscous flow. The potential of the tunnel is exemplified by the study of the drag on a simplified model of wing whose angle of attack can be controlled. A link to download the programs that make up the tunnel appears at the end.

  3. www.eprg.group.cam.ac.uk EPRGWORKINGPAPERNON-TECHNICALSUMMARY

    E-Print Network [OSTI]

    de Gispert, Adrià

    as uncertainties are resolved reduces significantly its expected Levelised Cost Of generating Electricity (LCOE deterministic economic model to identify its LCOE without considering uncertainty. In the second step #12;www.eprg.group.cam.ac.uk EPRGWORKINGPAPERNON-TECHNICALSUMMARY we identified sources of uncertainty in the design that would affect its LCOE

  4. Numerical modeling of gas migration into and through faulted sand reservoirs in Pabst Field (Main Pass East Block 259), northern Gulf of Mexico 

    E-Print Network [OSTI]

    Li, Yuqian

    2006-08-16

    allow gas communication among the sands. Meanwhile, three fault families break up the three sands into numerous compartments. A primary fault and large synthetic and antithetic faults act as gas migration pathways: the synthetic and antithetic faults...

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

    E-Print Network [OSTI]

    Pelanti, Marica

    flows is relevant in numerous areas of engineering, from naval and submarine systems design to aerospace and nuclear power plants technologies. Cavitating fluids are multiphase mixtures that often involve complex

  6. A QUARTER-CENTURY OF OBSERVATIONS OF COMET 10P/TEMPEL 2 AT LOWELL OBSERVATORY: CONTINUED SPIN-DOWN, COMA MORPHOLOGY, PRODUCTION RATES, AND NUMERICAL MODELING

    SciTech Connect (OSTI)

    Knight, Matthew M.; Schleicher, David G.; Schwieterman, Edward W.; Christensen, Samantha R. [Lowell Observatory, 1400 W. Mars Hill Road, Flagstaff, AZ 86001 (United States); Farnham, Tony L., E-mail: knight@lowell.edu [Department of Astronomy, University of Maryland, College Park, MD 20742-2421 (United States)

    2012-11-01

    We report on photometry and imaging of Comet 10P/Tempel 2 obtained at Lowell Observatory from 1983 through 2011. We measured a nucleus rotation period of 8.950 {+-} 0.002 hr from 16 nights of imaging acquired between 2010 September and 2011 January. This rotation period is longer than the period we previously measured in 1999, which was itself longer than the period measured in 1988, and demonstrates that Tempel 2 is continuing to spin down, presumably due to torques caused by asymmetric outgassing. A nearly linear jet was observed which varied little during a rotation cycle in both R and CN images acquired during the 1999 and 2010 apparitions. We measured the projected direction of this jet throughout the two apparitions and, under the assumption that the source region of the jet was near the comet's pole, determined a rotational pole direction of R.A./decl. = 151 Degree-Sign /+59 Degree-Sign from CN measurements and R.A./decl. = 173 Degree-Sign /+57 Degree-Sign from dust measurements (we estimate a circular uncertainty of 3 Degree-Sign for CN and 4 Degree-Sign for dust). Different combinations of effects likely bias both gas and dust solutions and we elected to average these solutions for a final pole direction of R.A./decl. = 162 Degree-Sign {+-} 11 Degree-Sign /+58 Degree-Sign {+-} 1 Degree-Sign . Photoelectric photometry was acquired on 3 nights in 1983, 2 nights in 1988, 19 nights in 1999/2000, and 10 nights in 2010/2011. The activity exhibited a steep 'turn-on' {approx}3 months prior to perihelion (the exact timing of which varies) and a relatively smooth decline after perihelion. The activity during the 1999 and 2010 apparitions was similar; limited data in 1983 and 1988 (along with IUE data from the literature) were systematically higher and the difference cannot be explained entirely by the smaller perihelion distance. We measured a 'typical' composition, in agreement with previous investigators. Monte Carlo numerical modeling with our pole solution best replicated the observed coma morphology for a source region located near a comet latitude of +80 Degree-Sign and having a radius of {approx}10 Degree-Sign . Our model reproduced the seasonal changes in activity, suggesting that the majority of Tempel 2's activity originates from a small active region located near the pole. We also find that a cosine-squared solar angle function gives the best fit as compared to a standard cosine function.

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

    SciTech Connect (OSTI)

    Neary, Vincent S [ORNL

    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-series and stage vs. cross-section area rating relationship.

  8. 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.

  9. Numerical Simulation of Laminar Reacting Flows with Complex Chemistry

    E-Print Network [OSTI]

    Bell, John B.

    Numerical Simulation of Laminar Reacting Flows with Complex Chemistry M S Day and J B Bell Lawrence: 47.40.Fw, 82.40.Py Submitted to: Combust. Theory Modelling #12;Numerical Simulation of Laminar

  10. Lateral boundary errors in regional numerical weather

    E-Print Network [OSTI]

    ?umer, Slobodan

    Lateral boundary errors in regional numerical weather prediction models Author: Ana Car Advisor, they describe evolution of atmospher - weather forecast. Every NWP model solves the same system of equations (1: assoc. prof. dr. Nedjeljka Zagar January 5, 2015 Abstract Regional models are used in many national

  11. Numerical Solution of 3D Poisson-Nernst-Planck Equations Coupled with Classical Density Functional Theory for Modeling Ion and Electron Transport in a Confined Environment

    SciTech Connect (OSTI)

    Meng, Da; Zheng, Bin; Lin, Guang; Sushko, Maria L.

    2014-08-29

    We have developed efficient numerical algorithms for the solution of 3D steady-state Poisson-Nernst-Planck equations (PNP) with excess chemical potentials described by the classical density functional theory (cDFT). The coupled PNP equations are discretized by finite difference scheme and solved iteratively by Gummel method with relaxation. The Nernst-Planck equations are transformed into Laplace equations through the Slotboom transformation. Algebraic multigrid method is then applied to efficiently solve the Poisson equation and the transformed Nernst-Planck equations. A novel strategy for calculating excess chemical potentials through fast Fourier transforms is proposed which reduces computational complexity from O(N2) to O(NlogN) where N is the number of grid points. Integrals involving Dirac delta function are evaluated directly by coordinate transformation which yields more accurate result compared to applying numerical quadrature to an approximated delta function. Numerical results for ion and electron transport in solid electrolyte for Li ion batteries are shown to be in good agreement with the experimental data and the results from previous studies.

  12. Freund Publishing House Ltd. International Journal of Nonlinear Sciences and Numerical Simulation 4, 239-250, 2003 Traffic Flow CA Model in Which Only the Cars Following

    E-Print Network [OSTI]

    Simulation 4, 239-250, 2003 Traffic Flow CA Model in Which Only the Cars Following the Trail of the Ahead Car from the NS model in that only the cars following the trail of the ahead car can be delayed. In other words, a car with spacing ahead longer than the car velocity limit M can not be delayed in the new model

  13. Multigrid methods for improving the variational data assimilation in numerical weather prediction

    E-Print Network [OSTI]

    Kwak, Do Young

    Multigrid methods for improving the variational data assimilation in numerical weather prediction conditions are needed to solve numerical weather prediction models: initial condition and boundary condition: numerical weather prediction, variational data assimilation, minimization procedure, multigrid methods, cell

  14. 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.

  15. Proceedings of the Workshop on Numerical Modeling of Thermohydrological Flow in Fractured Rock Masses, Feb. 19-20, 1980, Berkeley, CA

    E-Print Network [OSTI]

    Witherspoon, P.A.

    2010-01-01

    n The variable is the hydraulic head Energy T r a n s p o rhydraulic fractur­ ing process requires the modeling of the energy

  16. 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.

  17. Introduction to exact numerical computation

    E-Print Network [OSTI]

    Escardó, Martín

    Introduction to exact numerical computation Notes 2 Floating-point computation 4 3 Exact numerical (called the operational semantics) for computing a syntactical representative of the mathematical entity

  18. Threedimensional numerical simulation for various geometries

    E-Print Network [OSTI]

    Herbin, Raphaèle

    modelling and numerical simulation of natural gas­fed solid oxide cells (Solid Oxide Fuel Cell, SOFC) at a stationary regime. The principle of a Solid Oxide Fuel Cell (SOFC) is based on the conversion of the chemical is taken into account in the present model. The SOFC systems seem to be of great interest for use

  19. 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.

  20. 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.

  1. 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.

  2. Numerical modeling of the Snowmass Creek paleoglacier, Colorado, and climate in the Rocky Mountains during the Bull Lake glaciation (MIS 6)

    SciTech Connect (OSTI)

    Eric M. Leonard; Mitchell A. Plummer; Paul E. Carrara

    2014-04-01

    Well-preserved moraines from the penultimate, or Bull Lake, glaciation of Snowmass Creek Valley in the Elk Range of Colorado present an opportunity to examine the character of the high-altitude climate in the Rocky Mountains during Marine Oxygen Isotope Stage 6. This study employs a 2-D coupled mass/energy balance and flow model to assess the magnitudes of temperature and precipitation change that could have sustained the glacier in mass-balance equilibrium at its maximum extent during the Bull Lake glaciation. Variable substrate effects on glacier flow and ice thickness make the modeling somewhat more complex than in geologically simpler settings. Model results indicate that a temperature depression of about 6.7°C compared with the present (1971–2000 AD) would have been necessary to sustain the Snowmass Creek glacier in mass-balance equilibrium during the Bull Lake glaciation, assuming no change in precipitation amount or seasonality. A 50% increase or decrease from modern precipitation would have been coupled with 5.2°C and 9.1°C Bull Lake temperature depressions respectively. Uncertainty in these modeled temperature depressions is about 1°C.

  3. 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.

  4. Advanced Numerical Methods for Complex Environmental Models: Needs and Availability, 2013, 301-326 301 Istvn Farag, gnes Havasi and Zahari Zlatev (Eds.)

    E-Print Network [OSTI]

    Dimov, Ivan

    investigations related to air pollution in Bulgaria ([6, 7]), Denmark ([8-11]), England ([12]), Europe ([2, 4 been used in some inter-comparisons of European large-scale air pollution models ([20, 21]). Recently Publishers CHAPTER 4 Part C: Sensitivity of European Pollution Levels to Changes of Human-Made Emissions

  5. HGSYSTEM/UF{sub 6} model enhancements for plume rise and dispersion around buildings, lift-off of buoyant plumes, and robustness of numerical solver

    SciTech Connect (OSTI)

    Hanna, S.R.; Chang, J.C.

    1997-01-01

    The HGSYSTEM/UF{sub 6} model was developed for use in preparing Safety Analysis Reports (SARs) by estimating the consequences of possible accidental releases of UF{sub 6} to the atmosphere at the gaseous diffusion plants (GDPs) located in Portsmouth, Ohio, and Paducah, Kentucky. Although the latter report carries a 1996 date, the work that is described was completed in late 1994. When that report was written, the primary release scenarios of interest were thought to be gas pipeline and liquid tank ruptures over open terrain away from the influence of buildings. However, upon further analysis of possible release scenarios, the developers of the SARs decided it was necessary to also consider accidental releases within buildings. Consequently, during the fall and winter of 1995-96, modules were added to HGSYSTEM/UF{sub 6} to account for flow and dispersion around buildings. The original HGSYSTEM/UF{sub 6} model also contained a preliminary method for accounting for the possible lift-off of ground-based buoyant plumes. An improved model and a new set of wind tunnel data for buoyant plumes trapped in building recirculation cavities have become available that appear to be useful for revising the lift-off algorithm and modifying it for use in recirculation cavities. This improved lift-off model has been incorporated in the updated modules for dispersion around buildings.

  6. Freund Publishing House Ltd. International Journal of Nonlinear Sciences and Numerical Simulation 5(1) 5-8, 2004 Thermo-electro-hydrodynamic model for electrospinning process

    E-Print Network [OSTI]

    Pan, Ning

    5(1) 5-8, 2004 Thermo-electro-hydrodynamic model for electrospinning process Yu-Qin Wana , Qian Guob-fine polymer fibers, obtained by electrospinning, have a wide range of potential applications such as fluid behind electrospinning. When electrical force is dominant over the other forces, Bratu equation

  7. NUMERICAL SIMULATIONS OF TRANSVERSE COMPRESSION AND DENSIFICATION IN WOOD

    E-Print Network [OSTI]

    Nairn, John A.

    NUMERICAL SIMULATIONS OF TRANSVERSE COMPRESSION AND DENSIFICATION IN WOOD John A. Nairn1 Professor- terials is a useful tool for stress analysis and for failure modeling. Although FEA of wood as an anisotropic continuum is used, numerical modeling of realistic wood structures, including details of wood

  8. Experimental and numerical analysis of a deepwater mini-TLP 

    E-Print Network [OSTI]

    Guichard, Aurelien

    2001-01-01

    As the quest for oil and gas resources drives the industry to ever deeper waters, model testing still represents an essential step after numerical modeling when designing offshore platforms in these hostile environments. ...

  9. Numerical Weather Forecasting at the Savannah River Site

    SciTech Connect (OSTI)

    Buckley, R.L. [Westinghouse Savannah River Company, AIKEN, SC (United States)

    1998-08-01

    This paper discusses the use of an advanced three-dimensional prognostic numerical model to provide space and time-dependent meteorological data for use in the WIND System dispersion models.

  10. Observational and Numerical of the Adriatic Bora

    E-Print Network [OSTI]

    Gohm, Alexander

    and radiosoundings) 2. Mesoscale numerical model RAMS: · 5 or 6 nested grids · Horizontal mesh size as low as 800 m RAMS domain 5 #12;The event ­ 28 March 2002: · Deep bora: deep NE flow, no critical level field: · The bora flow weakens and separates during daytime. RAMS vertical transects across Vratnik Pass

  11. Numerical simulation of an electroweak oscillon

    SciTech Connect (OSTI)

    Graham, N. [Department of Physics, Middlebury College, Middlebury, Vermont 05753 (United States)

    2007-10-15

    Numerical simulations of the bosonic sector of the SU(2)xU(1) electroweak standard model in 3+1 dimensions have demonstrated the existence of an oscillon--an extremely long-lived, localized, oscillatory solution to the equations of motion--when the Higgs mass is equal to twice the W{sup {+-}} boson mass. It contains total energy roughly 30 TeV localized in a region of radius 0.05 fm. A detailed description of these numerical results is presented.

  12. 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.

  13. Status of NINJA: the Numerical INJection Analysis project

    E-Print Network [OSTI]

    Benjamin Aylott; John G. Baker; William D. Boggs; Michael Boyle; Patrick R. Brady; Duncan A. Brown; Bernd Brügmann; Luisa T. Buchman; Alessandra Buonanno; Laura Cadonati; Jordan Camp; Manuela Campanelli; Joan Centrella; Shourov Chatterjis; Nelson Christensen; Tony Chu; Peter Diener; Nils Dorband; Zachariah B. Etienne; Joshua Faber; Stephen Fairhurst; Benjamin Farr; Sebastian Fischetti; Gianluca Guidi; Lisa M. Goggin; Mark Hannam; Frank Herrmann; Ian Hinder; Sascha Husa; Vicky Kalogera; Drew Keppel; Lawrence E. Kidder; Bernard J. Kelly; Badri Krishnan; Pablo Laguna; Carlos O. Lousto; Ilya Mandel; Pedro Marronetti; Richard Matzner; Sean T. McWilliams; Keith D. Matthews; R. Adam Mercer; Satyanarayan R. P. Mohapatra; Abdul H. Mroué; Hiroyuki Nakano; Evan Ochsner; Yi Pan; Larne Pekowsky; Harald P. Pfeiffer; Denis Pollney; Frans Pretorius; Vivien Raymond; Christian Reisswig; Luciano Rezzolla; Oliver Rinne; Craig Robinson; Christian Röver; Lucía Santamaría; Bangalore Sathyaprakash; Mark A. Scheel; Erik Schnetter; Jennifer Seiler; Stuart L. Shapiro; Deirdre Shoemaker; Ulrich Sperhake; Alexander Stroeer; Riccardo Sturani; Wolfgang Tichy; Yuk Tung Liu; Marc van der Sluys; James R. van Meter; Ruslan Vaulin; Alberto Vecchio; John Veitch; Andrea Viceré; John T. Whelan; Yosef Zlochower

    2009-05-26

    The 2008 NRDA conference introduced the Numerical INJection Analysis project (NINJA), a new collaborative effort between the numerical relativity community and the data analysis community. NINJA focuses on modeling and searching for gravitational wave signatures from the coalescence of binary system of compact objects. We review the scope of this collaboration and the components of the first NINJA project, where numerical relativity groups shared waveforms and data analysis teams applied various techniques to detect them when embedded in colored Gaussian noise.

  14. Supplementary Information Potential for Electricity Generation from Renewable Resources and Levelized Cost of Electricity (LCOE)

    E-Print Network [OSTI]

    Suo, Zhigang

    the potential to meet the worldwide demand of electricity and they contribute to the total generation of providing enough energy to meet the world demand of electricity, the current amount of electricitySupplementary Information Potential for Electricity Generation from Renewable Resources

  15. New Report Shows Downward Trend in LCOE for Wind | Department of Energy

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

    AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Homesum_a_epg0_fpd_mmcf_m.xls" ,"Available from WebQuantity of Natural GasAdjustmentsShirleyEnergyTher i nAand DOEDepartmentNew Jersey is homeAdvanced Stages of Development ||report

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

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

    AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on Delicious RankADVANCED MANUFACTURING OFFICE INDUSTRIALU.S. DepartmentJeanKeyLANL Sustainable Design GuideON THE

  17. Astrophysical jets: Observations, numerical simulations, and laboratory experiments

    E-Print Network [OSTI]

    Bellan, Paul M.

    Astrophysical jets: Observations, numerical simulations, and laboratory experiments P. M. Bellan,1; published online 22 April 2009 This paper provides summaries of ten talks on astrophysical jets given of observation, numerical modeling, and laboratory experiment. One essential feature of jets, namely

  18. MULTISCALE NUMERICAL STUDY OF TURBULENT FLOW AND BUBBLE ENTRAINMENT

    E-Print Network [OSTI]

    Kirby, James T.

    MULTISCALE NUMERICAL STUDY OF TURBULENT FLOW AND BUBBLE ENTRAINMENT IN THE SURF ZONE BY GANGFENG MA . . . . . . . . . . . . . . . . . . . . . 3 1.1.2 Numerical Investigations . . . . . . . . . . . . . . . . . . . . . 5 1.2 Bubble Entrainment . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15 2.4 Bubble Entrainment Model . . . . . . . . . . . . . . . . . . . . . . . 18 2.5 Bubble Breakup

  19. Gnie mcanique Numerical Modeling of Water Splitting

    E-Print Network [OSTI]

    Psaltis, Demetri

    reaction takes place.. Laboratory of Renewable Energy Science and Engineering Parametric Study Annual & Tom Bosserez from KU Leuven , Mikaël Dumortier from EPFL Motivation & Objectives It has become accepted in wide ranges of societies that alternative sources of energy need to be developped. Using

  20. Numerical and experimental modelling of microwave applicators

    E-Print Network [OSTI]

    Dibben, David

    1995-10-24

    This thesis presents a time domain finite element method for the solution of microwave heating problems. This is the first time that this particular technique has been applied to microwave heating. It is found that the standard frequency domain...

  1. numerical modeling | OpenEI Community

    Open Energy Info (EERE)

    AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on QA:QA J-E-1 SECTION J APPENDIX ECoop IncIowa (UtilityMichigan)data book Homefuelleasing Homemapsmockupsnotes

  2. OpenEI Community - numerical modeling

    Open Energy Info (EERE)

    AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on QA:QA J-E-1 SECTION J APPENDIXsourceII JumpQuarterly Smart Grid Data available for download onst,/0Thoughts

  3. Numerical study of high heat ux pool boiling heat transfer Ying He a,*, Masahiro Shoji b

    E-Print Network [OSTI]

    Maruyama, Shigeo

    Numerical study of high heat ¯ux pool boiling heat transfer Ying He a,*, Masahiro Shoji b , Shigeo simulation model of boiling heat transfer is proposed based on a numerical macrolayer model [S. Maruyama, M. Shoji, S. Shimizu, A numerical simulation of transition boiling heat transfer, in: Proceedings

  4. Numerical Simulation of Laminar Reacting Flows with Complex Chemistry

    E-Print Network [OSTI]

    Bell, John B.

    Numerical Simulation of Laminar Reacting Flows with Complex Chemistry M S Day and J B Bell Lawrence Simulation of Laminar Reacting Flows 2 1. Introduction Detailed modelling of time-dependent reacting ows

  5. An investigation of the numerical treatment of condensation

    E-Print Network [OSTI]

    Sasson, Joseph

    The simulation of complete condensation continues to challenge the numerical methods currently used for multi-phase flow modeling; especially at low pressures, the change of phase process from a two-phase mixture to liquid ...

  6. NUMERICAL APPROXIMATIONS OF ALLEN-CAHN AND CAHN ...

    E-Print Network [OSTI]

    2010-03-18

    schemes that satisfy a corresponding discrete energy law, or in other words, energy ... Our first objective is to design stabilized semi-implicit schemes that satisfy an energy ...... Numerical analysis of a continuum model of phase transition.

  7. Evaluation of numerical weather prediction for intra-day solar forecasting in the continental United States

    E-Print Network [OSTI]

    Mathiesen, Patrick; Kleissl, Jan

    2011-01-01

    transport and  numerical weather modeling.   J.  Applied cross correlations.    Weather and Forecasting, 8:4, 401?of radiation for numerical weather prediction and climate 

  8. A numerical study of pseudoscalar inflation

    E-Print Network [OSTI]

    Cheng, Shu-Lin; Ng, Kin-Wang

    2015-01-01

    A numerical study of a pseudoscalar inflation having an axion-photon-like coupling is performed by solving numerically the coupled differential equations of motion for inflaton and photon mode functions from the onset of inflation to the end of reheating. The backreaction due to particle production is also included self-consistently. We find that this particular inflation model realizes the idea of a warm inflation in which a steady thermal bath is established by the particle production. In most cases this thermal bath exceeds the amount of radiation released in the reheating process. In the strong coupling regime, the transition from the inflationary to the radiation-dominated phase does not involve either a preheating or reheating process. In addition, energy density peaks produced near the end of inflation may lead to the formation of primordial black holes.

  9. A numerical study of pseudoscalar inflation

    E-Print Network [OSTI]

    Shu-Lin Cheng; Wolung Lee; Kin-Wang Ng

    2015-08-02

    A numerical study of a pseudoscalar inflation having an axion-photon-like coupling is performed by solving numerically the coupled differential equations of motion for inflaton and photon mode functions from the onset of inflation to the end of reheating. The backreaction due to particle production is also included self-consistently. We find that this particular inflation model realizes the idea of a warm inflation in which a steady thermal bath is established by the particle production. In most cases this thermal bath exceeds the amount of radiation released in the reheating process. In the strong coupling regime, the transition from the inflationary to the radiation-dominated phase does not involve either a preheating or reheating process. In addition, energy density peaks produced near the end of inflation may lead to the formation of primordial black holes.

  10. An Interactive Course in Numerical Methods for the Earth Susan E. Allen

    E-Print Network [OSTI]

    Stockie, John

    are typically in physics, mathematics, chemistry, biology or geography, and exposure to numerical computing science are becoming more and more dependent on numerical simulation, modelling, and computationally courses in numerical methods using MatlabTM as the programming language, or a graduate-level mathematics

  11. 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.

  12. Numerical simulation of water flow around a rigid fishing net

    E-Print Network [OSTI]

    Roger Lewandowski; Géraldine Pichot

    2006-12-20

    This paper is devoted to the simulation of the flow around and inside a rigid axisymmetric net. We describe first how experimental data have been obtained. We show in detail the modelization. The model is based on a Reynolds Averaged Navier-Stokes turbulence model penalized by a term based on the Brinkman law. At the out-boundary of the computational box, we have used a "ghost" boundary condition. We show that the corresponding variational problem has a solution. Then the numerical scheme is given and the paper finishes with numerical simulations compared with the experimental data.

  13. Merging of coronal and heliospheric numerical two-dimensional D. Odstrcil,1,5

    E-Print Network [OSTI]

    California at Berkeley, University of

    spatial and temporal scales between the Sun and Earth. Specialized physically based numerical models have been developed to address particular aspects of the entire system. However, an integrated modeling physically based numerical models have been developed to address particular aspects of the entire system

  14. Numerical simulation of tsunami waves generated by deformable submarine landslides

    E-Print Network [OSTI]

    Kirby, James T.

    wave model Tsunami wave Numerical modeling a b s t r a c t This paper presents a new submarine energy is mostly concentrated on a narrow band of the dominant slide direction for the waves generated-up along the coast. For example, submarine mass failure is considered as one of the major sources

  15. Numerical Approximation of Vortex Density Evolution in a Superconductor.

    E-Print Network [OSTI]

    Styles, Vanessa

    Numerical Approximation of Vortex Density Evolution in a Superconductor. C.M. Elliott & V. Styles Abstract A #12;nite volume/element approximation of a mean #12;eld model of superconducting vortices in one approximations of a two-dimensional version of the mean #12;eld model of superconducting vortices considered

  16. Testing gravitational-wave searches with numerical relativity waveforms: Results from the first Numerical INJection Analysis (NINJA) project

    E-Print Network [OSTI]

    Benjamin Aylott; John G. Baker; William D. Boggs; Michael Boyle; Patrick R. Brady; Duncan A. Brown; Bernd Brügmann; Luisa T. Buchman; Alessandra Buonanno; Laura Cadonati; Jordan Camp; Manuela Campanelli; Joan Centrella; Shourov Chatterji; Nelson Christensen; Tony Chu; Peter Diener; Nils Dorband; Zachariah B. Etienne; Joshua Faber; Stephen Fairhurst; Benjamin Farr; Sebastian Fischetti; Gianluca Guidi; Lisa M. Goggin; Mark Hannam; Frank Herrmann; Ian Hinder; Sascha Husa; Vicky Kalogera; Drew Keppel; Lawrence E. Kidder; Bernard J. Kelly; Badri Krishnan; Pablo Laguna; Carlos O. Lousto; Ilya Mandel; Pedro Marronetti; Richard Matzner; Sean T. McWilliams; Keith D. Matthews; R. Adam Mercer; Satyanarayan R. P. Mohapatra; Abdul H. Mroué; Hiroyuki Nakano; Evan Ochsner; Yi Pan; Larne Pekowsky; Harald P. Pfeiffer; Denis Pollney; Frans Pretorius; Vivien Raymond; Christian Reisswig; Luciano Rezzolla; Oliver Rinne; Craig Robinson; Christian Röver; Lucía Santamaría; Bangalore Sathyaprakash; Mark A. Scheel; Erik Schnetter; Jennifer Seiler; Stuart L. Shapiro; Deirdre Shoemaker; Ulrich Sperhake; Alexander Stroeer; Riccardo Sturani; Wolfgang Tichy; Yuk Tung Liu; Marc van der Sluys; James R. van Meter; Ruslan Vaulin; Alberto Vecchio; John Veitch; Andrea Viceré; John T. Whelan; Yosef Zlochower

    2009-07-09

    The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational-wave data analysis communities. The purpose of NINJA is to study the sensitivity of existing gravitational-wave search algorithms using numerically generated waveforms and to foster closer collaboration between the numerical relativity and data analysis communities. We describe the results of the first NINJA analysis which focused on gravitational waveforms from binary black hole coalescence. Ten numerical relativity groups contributed numerical data which were used to generate a set of gravitational-wave signals. These signals were injected into a simulated data set, designed to mimic the response of the Initial LIGO and Virgo gravitational-wave detectors. Nine groups analysed this data using search and parameter-estimation pipelines. Matched filter algorithms, un-modelled-burst searches and Bayesian parameter-estimation and model-selection algorithms were applied to the data. We report the efficiency of these search methods in detecting the numerical waveforms and measuring their parameters. We describe preliminary comparisons between the different search methods and suggest improvements for future NINJA analyses.

  17. Planning numerical approximations Richard Power

    E-Print Network [OSTI]

    Williams, Sandra

    Planning numerical approximations Richard Power Sandra Williams 21st September 2009 #12;Table proportions (e.g., more than a quarter, 25.9 per cent) Proportions are a convenient well-defined subproblem Common in factual discourse (e.g., newspaper articles) Important for generating from data (but neglected

  18. Graduate Survey of Numerical Methods Background material

    E-Print Network [OSTI]

    Corless, Robert M.

    Numerical Methods, Numerical Analysis, Scientific Computing, Com- putational Mathematics, Computational. Numerical Meth- ods are what one uses to solve a problem from continuous mathematics (vide Nick Trefethen), and in what sense. Scientific Computing is the use of numerical methods to solve problems of scientific

  19. Numerical semigroups Easy bounds on ng

    E-Print Network [OSTI]

    Elizalde, Sergi

    Numerical semigroups Easy bounds on ng Improved bounds on ng The number of numerical semigroups #12;Numerical semigroups Easy bounds on ng Improved bounds on ng Motivation Definitions The tree on ng Improved bounds on ng Motivation Definitions The tree T of numerical semigroups Succession rules

  20. 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.

  1. Numerical Simulations of Gamma-Ray Burst Explosions

    E-Print Network [OSTI]

    Lazzati, Davide; López-Cámara, Diego

    2015-01-01

    Gamma-ray bursts are a complex, non-linear system that evolves very rapidly through stages of vastly different conditions. They evolve from scales of few hundred kilometers where they are very dense and hot to cold and tenuous on scales of parsecs. As such, our understanding of such a phenomenon can truly increase by combining theoretical and numerical studies adopting different numerical techniques to face different problems and deal with diverse conditions. In this review, we will describe the tremendous advancement in our comprehension of the bursts phenomenology through numerical modeling. Though we will discuss studies mainly based on jet dynamics across the progenitor star and the interstellar medium, we will also touch upon other problems such as the jet launching, its acceleration, and the radiation mechanisms. Finally, we will describe how combining numerical results with observations from Swift and other instruments resulted in true understanding of the bursts phenomenon and the challenges still lyi...

  2. The effective viscosity of rocksalt: implementation of steady-state creep laws in numerical

    E-Print Network [OSTI]

    van Keken, Peter

    of salt tectonics have received renewed attention in the last few years. This has been stimulated in numerical models of salt diapirism P.E. van Keken, C.J. Spiers, A.P. van den Berg and E.J. Muyzert and fluid-enhanced grain-boundary diffusion creep, has been used in numerical models of salt diapirism

  3. Numerical Methods for a Nonlinear BVP Arising in Physical Oceanography

    E-Print Network [OSTI]

    Riccardo Fazio; Alessandra Jannelli

    2013-10-08

    In this paper we report and compare the numerical results for an ocean circulation model obtained by the classical truncated boundary formulation, the free boundary approach and a quasi-uniform grid treatment of the problem. We apply a shooting method to the truncated boundary formulation and finite difference methods to both the free boundary approach and the quasi-uniform grid treatment. Using the shooting method, supplemented by the Newton's iterations, we show that the ocean circulation model cannot be considered as a simple test case. In fact, for this method we are forced to use as initial iterate a value close to the correct missing initial condition in order to be able to get a convergent numerical solution. The reported numerical results allow us to point out how the finite difference method with a quasi-uniform grid is the less demanding approach and that the free boundary approach provides a more reliable formulation than the classical truncated boundary formulation.

  4. NUMERICAL SIMULATION OF NATURAL GAS-SWIRL BURNER

    SciTech Connect (OSTI)

    Ala Qubbaj

    2005-03-01

    A numerical simulation of a turbulent natural gas jet diffusion flame at a Reynolds number of 9000 in a swirling air stream is presented. The numerical computations were carried out using the commercially available software package CFDRC. The instantaneous chemistry model was used as the reaction model. The thermal, composition, flow (velocity), as well as stream function fields for both the baseline and air-swirling flames were numerically simulated in the near-burner region, where most of the mixing and reactions occur. The results were useful to interpret the effects of swirl in enhancing the mixing rates in the combustion zone as well as in stabilizing the flame. The results showed the generation of two recirculating regimes induced by the swirling air stream, which account for such effects. The present investigation will be used as a benchmark study of swirl flow combustion analysis as a step in developing an enhanced swirl-cascade burner technology.

  5. Numerical Simulations of Bouncing Jets

    E-Print Network [OSTI]

    Bonito, Andrea; Lee, Sanghyun

    2015-01-01

    Bouncing jets are fascinating phenomenons occurring under certain conditions when a jet impinges on a free surface. This effect is observed when the fluid is Newtonian and the jet falls in a bath undergoing a solid motion. It occurs also for non-Newtonian fluids when the jets falls in a vessel at rest containing the same fluid. We investigate numerically the impact of the experimental setting and the rheological properties of the fluid on the onset of the bouncing phenomenon. Our investigations show that the occurrence of a thin lubricating layer of air separating the jet and the rest of the liquid is a key factor for the bouncing of the jet to happen. The numerical technique that is used consists of a projection method for the Navier-Stokes system coupled with a level set formulation for the representation of the interface. The space approximation is done with adaptive finite elements. Adaptive refinement is shown to be very important to capture the thin layer of air that is responsible for the bouncing.

  6. Class Generation for Numerical Wind Atlases

    E-Print Network [OSTI]

    Class Generation for Numerical Wind Atlases Risø National Laboratory Wind Energy Department Constructing a Numerical Wind Atlas 5 2.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . 13 2.4 The existing procedure at Risø . . . . . . . . . . . . . . . . . . . 13 3 Representing a Wind

  7. On the numerical solution of some problems of environmental pollution

    E-Print Network [OSTI]

    Ehrhardt, Matthias

    On the numerical solution of some problems of environmental pollution Quang A Dang1 , Matthias world and their importance will even increase in the future. High pollution of air, water and soil may to the mathematical formula- tion of air pollution models. Let G be a cylindrical domain in the three

  8. 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.

  9. RisR1252(EN) The Numerical Wind Atlas

    E-Print Network [OSTI]

    Risø­R­1252(EN) The Numerical Wind Atlas -- the KAMM/WAsP Method Helmut P. Frank, Ole Rathmann The method of combining the Karlsruhe Atmospheric Mesoscale Model, KAMM, with the Wind Atlas Analysis and Application Program, WAsP, to make local predictions of the wind resource is presented. It combines

  10. Capturing Incomplete Information in Resource Allocation Problems through Numerical

    E-Print Network [OSTI]

    Powell, Warren B.

    Capturing Incomplete Information in Resource Allocation Problems through Numerical Patterns Arun information in resource allocation models. We use a popular goodness-of-fit measure known as the Cramer-Von Mises metric based on the empirical distribution function (EDF) as the foundation of our approach. We

  11. Hydraulic Fractures: multiscale phenomena, asymptotic and numerical solutions

    E-Print Network [OSTI]

    Peirce, Anthony

    Hydraulic Fractures: multiscale phenomena, asymptotic and numerical solutions SANUM Conference (UMN) Eduard Siebrits (SLB) #12;2 Outline · Examples of hydraulic fractures · Governing equations well stimulation Fracturing Fluid Proppant #12;5 Quarries #12;6 Magma flow Tarkastad #12;7 Model EQ 1

  12. Hydraulic Fractures: multiscale phenomena, asymptotic and numerical solutions

    E-Print Network [OSTI]

    Peirce, Anthony

    Hydraulic Fractures: multiscale phenomena, asymptotic and numerical solutions CSIRO CSS TCP Detournay (UMN) Eduard Siebrits (SLB) #12;2 Outline · Examples of hydraulic fractures · Governing equations well stimulation Fracturing Fluid Proppant #12;5 Quarries #12;6 Magma flow Tarkastad #12;7 Model EQ 1

  13. Landolt-Bornstein Numerical Data and Functional Relationships

    E-Print Network [OSTI]

    Landolt-B¨ornstein Numerical Data and Functional Relationships NEW SERIES 4.1.1.4: Solar, the temperature structure of theoretical models is obtained by solving an energy equation. In the solar atmosphere Potsdam, Germany msteffen@aip.de October 29, 2009 #12;#12;Ref. p. 44] 4.1.1.4 Solar photosphere

  14. NUMERICAL SOLUTION OF THE POLYMER SYSTEM BY FRONT TRACKING

    E-Print Network [OSTI]

    , models two­phase, three­component polymer flooding in an oil reservoir. Two examples are presented, where. Key words and phrases. Front tracking, nonstrictly hyperbolic systems, polymer flooding, dimensionalNUMERICAL SOLUTION OF THE POLYMER SYSTEM BY FRONT TRACKING V. HAUGSE, K. H. KARLSEN, K.--A. LIE

  15. Numerical Convergence Properties of Option Pricing PDEs with Uncertain Volatility

    E-Print Network [OSTI]

    Forsyth, Peter A.

    is then simply the cost of this hedge. The principal source of risk is the price of the underlying assetNumerical Convergence Properties of Option Pricing PDEs with Uncertain Volatility D. M. Pooley , P, 2001 Abstract The pricing equations derived from uncertain volatility models in finance are often cast

  16. An implicit numerical algorithm general relativistic hydrodynamics

    E-Print Network [OSTI]

    A. Hujeirat

    2008-01-09

    An implicit numerical algorithm general relativistic hydrodynamics This article has been replaced by arXiv:0801.1017

  17. Petroleum Engineering 301 Petroleum Engineering Numerical Methods

    E-Print Network [OSTI]

    10 Petroleum Engineering 301 Petroleum Engineering Numerical Methods Credit 3: (2-3) Required for Juniors Catalog Description: Use of numerical methods in a variety of petroleum engineering problems methods. Prerequisites(s): MATH 308 Textbook Required: Numerical Methods for Engineers, 6th Edition

  18. From Numerical Analysis to Computational Science

    E-Print Network [OSTI]

    Li, Tiejun

    . Numerical computing has, of course, been part of mathematics for a very long time. Al- gorithms by the namesFrom Numerical Analysis to Computational Science Bj¨orn Engquist · Gene Golub 1. Introduction The modern development of numerical computing is driven by the rapid in- crease in computer performance

  19. Simulating Reionization in Numerical Cosmology

    E-Print Network [OSTI]

    Aaron Sokasian; Tom Abel; Lars E. Hernquist

    2001-05-10

    The incorporation of radiative transfer effects into cosmological hydrodynamical simulations is essential for understanding how the intergalactic medium (IGM) makes the transition from a neutral medium to one that is almost fully ionized. Here, we present an approximate numerical method designed to study in a statistical sense how a cosmological density field is ionized by a set of discrete point sources. A diffuse background radiation field is also computed self-consistently in our procedure. The method requires relatively few time steps and can be employed with simulations having high resolution. We describe the details of the algorithm and provide a description of how the method can be applied to the output from a pre-existing cosmological simulation to study the systematic reionization of a particular ionic species. As a first application, we compute the reionization of He II by quasars in the redshift range 3 to 6.

  20. Area E Numerics and Scientific Computing Foundation in Numerical Mathematics F4E1

    E-Print Network [OSTI]

    Sturm, Karl-Theodor

    Area E ­ Numerics and Scientific Computing Foundation in Numerical Mathematics F4E1: · Scientific Selected Topics in Scientific Computing · V5E3 Advanced Topics in Numerical Methods in Science matrices) · Computational Finance (e.g. option pricing, fast numerical methods) · Visualization Methods (e

  1. NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. 2005; 12:683

    E-Print Network [OSTI]

    NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. 2005; 12:683 Published Numerical Linear Algebra and its Applications The fourth workshop of the ERCIM Working Group on `Matrix Computations and Statistics' and the First International workshop on `Numerical Linear Algebra and its

  2. Turbulent Combustion in Type Ia Supernova Models

    E-Print Network [OSTI]

    F. K. Roepke; W. Hillebrandt

    2006-09-15

    We review the astrophysical modeling of type Ia supernova explosions and describe numerical methods to implement numerical simulations of these events. Some results of such simulations are discussed.

  3. Comparison of integral equation and physical scale modelling of the electromagnetic response of models

    E-Print Network [OSTI]

    Farquharson, Colin G.

    Comparison of integral equation and physical scale modelling of the electromagnetic response history of EM numerical modelling in geophysics. · Another integral equation modelling program;Introduction: a brief history · Two main approaches to numerical modelling: integral equation; finite

  4. Direct numerical simulation of turbulent reacting flows

    SciTech Connect (OSTI)

    Chen, J.H. [Sandia National Laboratories, Livermore, CA (United States)

    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.

  5. Non-Gaussian numerical errors versus mass hierarchy

    E-Print Network [OSTI]

    Y. Meurice; M. B. Oktay

    2000-05-12

    We probe the numerical errors made in renormalization group calculations by varying slightly the rescaling factor of the fields and rescaling back in order to get the same (if there were no round-off errors) zero momentum 2-point function (magnetic susceptibility). The actual calculations were performed with Dyson's hierarchical model and a simplified version of it. We compare the distributions of numerical values obtained from a large sample of rescaling factors with the (Gaussian by design) distribution of a random number generator and find significant departures from the Gaussian behavior. In addition, the average value differ (robustly) from the exact answer by a quantity which is of the same order as the standard deviation. We provide a simple model in which the errors made at shorter distance have a larger weight than those made at larger distance. This model explains in part the non-Gaussian features and why the central-limit theorem does not apply.

  6. NUMERICAL STUDY OF A TURBULENT HYDRAULIC JUMP Qun Zhao 1 Shubhra K. Misra1

    E-Print Network [OSTI]

    Zhao, Qun

    . Hydraulic jumps are commonly used as energy dissipators and they have been studied intensively by hydraulicNUMERICAL STUDY OF A TURBULENT HYDRAULIC JUMP Qun Zhao 1 Shubhra K. Misra1 Ib A. Svendsen 1 (Member of a turbulent hydraulic jump. The numerical model is based on RIPPLE (Kothe et al., 1994) with two turbulence

  7. Numerical study on optimal Stirling engine regenerator matrix designs taking into account the effects of matrix

    E-Print Network [OSTI]

    Numerical study on optimal Stirling engine regenerator matrix designs taking into account matrix design that improves the efficiency of a Stirling engine has been developed in a numerical study of the existing SM5 Stirling engine. A new, detailed, one-dimensional Stirling engine model that delivers results

  8. Physical mechanisms involved in grooved flat heat pipes: experimental and numerical analyses

    E-Print Network [OSTI]

    Paris-Sud XI, Université de

    1 Physical mechanisms involved in grooved flat heat pipes: experimental and numerical analyses.lefevre@insa-lyon.fr) Abstract An experimental database, obtained with flat plate heat pipes (FPHP) with longitudinal grooves to results of numerical models. Physical mechanisms involved in grooved heat pipes are discussed, including

  9. NUMERICAL SIMULATIONS OF LONG TERM UNSATURATED FLOW AND ACID MINE DRAINAGE AT WASTE ROCK PILES

    E-Print Network [OSTI]

    Aubertin, Michel

    NUMERICAL SIMULATIONS OF LONG TERM UNSATURATED FLOW AND ACID MINE DRAINAGE AT WASTE ROCK PILES Omar present a numerical modeling study of unsaturated water flow and acid mine drainage in idealized (but of oxygen diffusion and acid mine drainage through the waste rock piles showed that oxygen is generally

  10. Impact of Turbulence Closures and Numerical Errors for the Optimization of Flow Control Devices

    E-Print Network [OSTI]

    Paris-Sud XI, Université de

    Impact of Turbulence Closures and Numerical Errors for the Optimization of Flow Control Devices J the use of a Kriging-based global optimization method to determine optimal control parameters conduct an optimization process and measure the impact of numerical and modeling errors on the optimal

  11. Observations and Numerical Simulations of Urban Heat Island and Sea Breeze Circulations over New York City

    E-Print Network [OSTI]

    Raman, Sethu

    (HIC) associated with an urban area can significantly alter lower tropospheric winds and lowObservations and Numerical Simulations of Urban Heat Island and Sea Breeze Circulations over New disaster on September 11, 2001. An ARPS model numerical simulation was conducted to explore the complex

  12. 1, 303350, 2007 A numerical

    E-Print Network [OSTI]

    Paris-Sud XI, Université de

    -induced snow transport for a high5 alpine area in the Berchtesgaden National Park (Germany), and we discuss sublimation from the surface15 of the trees, and a model for the simulation of wind-induced snow transport heat of fusion (melting or freezing) and the latent heat of vaporization (evaporation or condensation

  13. Numerical Weather Forecasting at the Savannah River Site

    SciTech Connect (OSTI)

    Buckley, R.L. [Westinghouse Savannah River Company, AIKEN, SC (United States)

    1998-11-01

    Weather forecasts at the Savannah River Site (SRS) are important for applications to emergency response. The fate of accidentally-released radiological materials and toxic chemicals can be determined by providing wind and turbulence input to atmospheric transport models. This operation has been routinely performed at SRS using the WIND System, a system of computer models and monitors which collect data from towers situated throughout the SRS. However, the information provided to these models is spatially homogeneous (in one or two dimensions) with an elementary forecasting capability. This paper discusses the use of an advanced three-dimensional prognostic numerical model to provide space and time-dependent meteorological data for use in the WIND System dispersion models. The extensive meteorological data collection at SRS serves as a ground truth for further model development as well as for use in other applications.

  14. Numerical Study of a Turbulent Hydraulic Jump

    E-Print Network [OSTI]

    Zhao, Qun

    Numerical Study of a Turbulent Hydraulic Jump Qun Zhao, Shubhra Misra, Ib. A. Svendsen and James T of a Turbulent Hydraulic Jump ­ p.1/14 #12;Objective Our ultimate goal is to study the breaking waves. Numerical Study of a Turbulent Hydraulic Jump ­ p.2/14 #12;A moving bore Qiantang Bore China (Courtesy of Dr J

  15. 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.

  16. Direct numerical simulation of electroconvective instability and hysteretic current-voltage response of a permselective membrane

    E-Print Network [OSTI]

    Pham, Van Sang

    We present a systematic, multiscale, fully detailed numerical modeling for dynamics of fluid flow and ion transport covering Ohmic, limiting, and overlimiting current regimes in conductance of ion-selective membrane. By ...

  17. Numerical Investigation of the Effect of Chirality of Carbon Nanotube on the Interfacial Thermal Resistance 

    E-Print Network [OSTI]

    Hu, Yuzhu

    2014-06-05

    . In this study, the interfacial thermal resistance between a carbon nanotube (CNT) and carbonate molten salt eutectics were calculated by using numerical models that were then implemented in Molecular Dynamics (MD) simulations. The estimates for “R_(k)” obtained...

  18. Three-dimensional numerical manifold method simulations for blocky rock analysis

    E-Print Network [OSTI]

    Shentu, Longfei

    2011-01-01

    After decades of development, people realize that there are wider and more various applications of numerical modeling and analysis. However, current feasible software tools cannot satisfy engineering and commercial goals. ...

  19. t -software package for numerical simulations of radioactive contaminant transport in groundwater

    E-Print Network [OSTI]

    Frolkovic, Peter

    in groundwater Peter Frolkovic Michael Lampe Gabriel Wittum September 16, 2005 Abstract The software package r3t or as numerical solution of some groundwater flow model, e.g., the density driven flow problem. The matrix Di = Di

  20. On the numerical simulation of the instationary free fall of a solid in a fluid.

    E-Print Network [OSTI]

    from meteorology, sedimentology and aerospace engineering to biology. We present a new numerical method], sedimentology [1], aerospace engineering [21] and biology (e.g. models for animal flights [2]). Recently

  1. A numerical study of compressible flow in labyrinth seals with rub grooves 

    E-Print Network [OSTI]

    Adams, Richard Gordon

    1997-01-01

    A labyrinth sea] uses a series of knife blades and cavities to create a high frictional flow path that turns a pressure potential into kinetic energy, which is then dissipated as heat. An existing numerical model is extended ...

  2. NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. 2011; 18:961980

    E-Print Network [OSTI]

    De Sterck, Hans

    NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. 2011; 18:961­980 Published aggregation method of [1], and the Markov chain algebraic multigrid (MCAMG) method of [3] using the OTF

  3. CHAPTER 1 -NUMERICAL INTEGRATION There are many techniques for performing numerical integrations, although no

    E-Print Network [OSTI]

    Boal, David

    . (1992) for further reading. There are several approaches to the integration of differential equationsCHAPTER 1 - NUMERICAL INTEGRATION There are many techniques for performing numerical integrations or numerical calculation of derivatives. Another approach, which we show for completeness in Sec. 1.1, is based

  4. Barchan dune asymmetry: Numerical investigation

    E-Print Network [OSTI]

    Parteli, Eric J R; Bourke, Mary C; Tsoar, Haim; Poeschel, Thorsten; Herrmann, Hans J

    2013-01-01

    Barchan dunes --- crescent-shaped dunes that form in areas of unidirectional winds and low sand availability --- commonly display an asymmetric shape, with one limb extended downwind. Several factors have been identified as potential causes for barchan dune asymmetry on Earth and Mars: asymmetric bimodal wind regime, topography, influx asymmetry and dune collision. However, the dynamics and potential range of barchan morphologies emerging under each specific scenario that leads to dune asymmetry are far from being understood. In the present work, we use dune modeling in order to investigate the formation and evolution of asymmetric barchans. We find that a bimodal wind regime causes limb extension when the divergence angle between primary and secondary winds is larger than $90^{\\circ}$, whereas the extended limb evolves into a seif dune if the ratio between secondary and primary transport rates is larger than 25%. Calculations of dune formation on an inclined surface under constant wind direction also lead to...

  5. NUMERICAL SIMULATIONS OF CHROMOSPHERIC MICROFLARES

    SciTech Connect (OSTI)

    Jiang, R. L.; Fang, C.; Chen, P. F., E-mail: fangc@nju.edu.c [Department of Astronomy, Nanjing University, Nanjing 210093 (China)

    2010-02-20

    With gravity, ionization, and radiation being considered, we perform 2.5 dimensional (2.5D) compressible resistive magnetohydrodynamic (MHD) simulations of chromospheric magnetic reconnection using the CIP-MOCCT scheme. The temperature distribution of the quiet-Sun atmospheric model VALC and the helium abundance (10%) are adopted. Our 2.5D MHD simulation reproduces qualitatively the temperature enhancement observed in chromospheric microflares. The temperature enhancement DELTAT is demonstrated to be sensitive to the background magnetic field, whereas the total evolution time DELTAt is sensitive to the magnitude of the anomalous resistivity. Moreover, we found a scaling law, which is described as DELTAT/DELTAt {approx} n{sub H} {sup -1.5} B {sup 2.1}eta{sub 0} {sup 0.88}. Our results also indicate that the velocity of the upward jet is much greater than that of the downward jet, and the X-point may move up or down.

  6. Connecting Numerical Relativity and Data Analysis of Gravitational Wave Detectors

    E-Print Network [OSTI]

    Deirdre Shoemaker; Karan Jani; Lionel London; Larne Pekowsky

    2015-03-09

    Gravitational waves deliver information in exquisite detail about astrophysical phenomena, among them the collision of two black holes, a system completely invisible to the eyes of electromagnetic telescopes. Models that predict gravitational wave signals from likely sources are crucial for the success of this endeavor. Modeling binary black hole sources of gravitational radiation requires solving the Eintein equations of General Relativity using powerful computer hardware and sophisticated numerical algorithms. This proceeding presents where we are in understanding ground-based gravitational waves resulting from the merger of black holes and the implications of these sources for the advent of gravitational-wave astronomy.

  7. 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.

  8. On Splitting-Based Numerical Methods for Convection-Diffusion Equations

    E-Print Network [OSTI]

    Chertock, Alina

    the performance of the splitting-based methods on a number of numerical exam- ples including the polymer system arising in modeling of the flooding processes in enhanced oil recovery, systems modeling the propagation mechanics, astrophysics, meteorology, multiphase flow in oil reservoirs, polymer flow, financial modeling

  9. COMBINED MODELING OF ACCELERATION, TRANSPORT, AND HYDRODYNAMIC...

    Office of Scientific and Technical Information (OSTI)

    COMBINED MODELING OF ACCELERATION, TRANSPORT, AND HYDRODYNAMIC RESPONSE IN SOLAR FLARES. I. THE NUMERICAL MODEL Citation Details In-Document Search Title: COMBINED MODELING OF...

  10. SAND2013-7207

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

    to calculate the LCOE of various devices. 20 This model has lead Oceantec to select a new platform (an offshore oscillating water column in a point absorber configuration) to...

  11. Initial Conditions for Numerical Relativity -- Introduction to numerical methods for solving elliptic PDEs

    E-Print Network [OSTI]

    Okawa, Hirotada

    2013-01-01

    Numerical relativity became a powerful tool to investigate the dynamics of binary problems with black holes or neutron stars as well as the very structure of General Relativity. Although public numerical relativity codes are available to evolve such systems, a proper understanding of the methods involved is quite important. Here we focus on the numerical solution of elliptic partial differential equations. Such equations arise when preparing initial data for numerical relativity, but also for monitoring the evolution of black holes. Because such elliptic equations play an important role in many branches of physics, we give an overview of the topic, and show how to numerically solve them with simple examples and sample codes written in C++ and Fortran90 for beginners in numerical relativity or other fields requiring numerical expertise.

  12. Initial Conditions for Numerical Relativity -- Introduction to numerical methods for solving elliptic PDEs

    E-Print Network [OSTI]

    Hirotada Okawa

    2013-08-15

    Numerical relativity became a powerful tool to investigate the dynamics of binary problems with black holes or neutron stars as well as the very structure of General Relativity. Although public numerical relativity codes are available to evolve such systems, a proper understanding of the methods involved is quite important. Here we focus on the numerical solution of elliptic partial differential equations. Such equations arise when preparing initial data for numerical relativity, but also for monitoring the evolution of black holes. Because such elliptic equations play an important role in many branches of physics, we give an overview of the topic, and show how to numerically solve them with simple examples and sample codes written in C++ and Fortran90 for beginners in numerical relativity or other fields requiring numerical expertise.

  13. Numerical weather forecasting at the Savannah River Site

    SciTech Connect (OSTI)

    Buckley, R.L. [Westinghouse Savannah River Site, Aiken, SC (United States)

    1998-12-31

    Weather forecasts at the Savannah River Site (SRS) are important for applications to emergency response. The fate of accidentally released radiological materials and toxic chemicals can be determined by providing wind and turbulence input to atmospheric transport models. This operation has been routinely performed at SRS using the WIND system, a system of computer models and monitors that collects data from towers situated throughout the SRS. However, the information provided to these models is spatially homogeneous (in one or two dimensions) with an elementary forecasting capability. This paper discusses the use of an advanced three-dimensional prognostic numerical model to provide space- and time-dependent meteorological data for use in the WIND system dispersion models. The extensive meteorological data collection at SRS serves as a ground truth for further model development as well as for use in other applications. A prognostic mesoscale model, the regional atmospheric modeling system (RAMS), is used to provide these forecasts. Use of RAMS allows for incorporation of mesoscale features such as the sea breeze, which has been shown to affect local weather conditions. This paper discusses the mesoscale model and its configuration for the operational simulation, as well as an application using a dispersion model at the SRS.

  14. 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.

  15. The Numerical Simulation Of A Transitional Flow In The VKI-GENOA Turbine Cascade

    E-Print Network [OSTI]

    Yershov, Sergiy; Yakovlev, Viktor; Gryzun, Maria

    2015-01-01

    This study presents a numerical simulation of a 3D viscous flow in the VKI-Genoa cascade that takes into account the laminar-turbulent transition. The numerical simulation is performed using the Reynolds-averaged Navier-Stokes equations and the two-equation k-omega SST turbulence model. The algebraic Production Term Modification model is used for modeling the laminar-turbulent transition. Computations of both fully turbulent and transitional flows are carried out. The contours of the Mach number, the turbulence kinetic energy, the entropy function, as well as limiting streamlines are presented. The analysis of the numerical results demonstrates the influence of the laminar-turbulent transition on the secondary flow pattern. The comparison between the present computational results and the existing experimental and numerical data shows that the proposed approach reflects sufficiently the physics of the laminar-turbulent transition in turbine cascades.

  16. Numerical determination of the material properties of porous dust cakes

    E-Print Network [OSTI]

    D. Paszun; C. Dominik

    2008-02-13

    The formation of planetesimals requires the growth of dust particles through collisions. Micron-sized particles must grow by many orders of magnitude in mass. In order to understand and model the processes during this growth, the mechanical properties, and the interaction cross sections of aggregates with surrounding gas must be well understood. Recent advances in experimental (laboratory) studies now provide the background for pushing numerical aggregate models onto a new level. We present the calibration of a previously tested model of aggregate dynamics. We use plastic deformation of surface asperities as the physical model to bring critical velocities for sticking into accordance with experimental results. The modified code is then used to compute compression strength and the velocity of sound in the aggregate at different densities. We compare these predictions with experimental results and conclude that the new code is capable of studying the properties of small aggregates.

  17. Spherically symmetric cosmological spacetimes with dust and radiation — numerical implementation

    SciTech Connect (OSTI)

    Lim, Woei Chet [Department of Mathematics, University of Waikato, Private Bag 3105, Hamilton 3240 (New Zealand); Regis, Marco [Dipartimento di Fisica, Università di Torino and INFN, Torino (Italy); Clarkson, Chris, E-mail: wclim@waikato.ac.nz, E-mail: regis@to.infn.it, E-mail: chris.clarkson@gmail.com [Astrophysics, Cosmology and Gravity Centre, and Department of Mathematics and Applied Mathematics, University of Cape Town, Rondebosch 7701, Cape Town (South Africa)

    2013-10-01

    We present new numerical cosmological solutions of the Einstein Field Equations. The spacetime is spherically symmetric with a source of dust and radiation approximated as a perfect fluid. The dust and radiation are necessarily non-comoving due to the inhomogeneity of the spacetime. Such a model can be used to investigate non-linear general relativistic effects present during decoupling or big-bang nucleosynthesis, as well as for investigating void models of dark energy with isocurvature degrees of freedom. We describe the full evolution of the spacetime as well as the redshift and luminosity distance for a central observer. After demonstrating accuracy of the code, we consider a few example models, and demonstrate the sensitivity of the late time model to the degree of inhomogeneity of the initial radiation contrast.

  18. Toward standard testbeds for numerical relativity

    E-Print Network [OSTI]

    Miguel Alcubierre; Gabrielle Allen; Carles Bona; David Fiske; Tom Goodale; F. Siddharta Guzman; Ian Hawke; Scott H. Hawley; Sascha Husa; Michael Koppitz; Christiane Lechner; Denis Pollney; David Rideout; Marcelo Salgado; Erik Schnetter; Edward Seidel; Hisa-aki Shinkai; Bela Szilagyi; Deirdre Shoemaker; Ryoji Takahashi; Jeffrey Winicour

    2003-05-06

    In recent years, many different numerical evolution schemes for Einstein's equations have been proposed to address stability and accuracy problems that have plagued the numerical relativity community for decades. Some of these approaches have been tested on different spacetimes, and conclusions have been drawn based on these tests. However, differences in results originate from many sources, including not only formulations of the equations, but also gauges, boundary conditions, numerical methods, and so on. We propose to build up a suite of standardized testbeds for comparing approaches to the numerical evolution of Einstein's equations that are designed to both probe their strengths and weaknesses and to separate out different effects, and their causes, seen in the results. We discuss general design principles of suitable testbeds, and we present an initial round of simple tests with periodic boundary conditions. This is a pivotal first step toward building a suite of testbeds to serve the numerical relativists and researchers from related fields who wish to assess the capabilities of numerical relativity codes. We present some examples of how these tests can be quite effective in revealing various limitations of different approaches, and illustrating their differences. The tests are presently limited to vacuum spacetimes, can be run on modest computational resources, and can be used with many different approaches used in the relativity community.

  19. 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.

  20. Numerical Estimation of Frictional Torques with Rate and State Friction

    E-Print Network [OSTI]

    Arun K. Singh; T. N. Singh

    2015-01-20

    In this paper, numerical estimation of frictional torques is carried out of a rotary elastic disc on a hard and rough surface under different rotating conditions. A one dimensional spring- mass rotary system is numerically solved under the quasistatic condition with the rate and state dependent friction model. It is established that torque of frictional strength as well as torque of steady dynamic stress increases with radius and found to be maximum at the periphery of the disc. Torque corresponding to frictional strength estimated using the analytical solution matches closely with the simulation only in the case of high stiffness of the connecting spring. In steady relaxation simulation, a steadily rotating disc is suddenly stopped and relaxational angular velocity and corresponding frictional torque decreases with both steady angular velocity and stiffness of the connecting spring in the velocity strengthening regime. In velocity weakening regime, in contrast, torque of relaxation stress deceases but relaxation velocity increases. The reason for the contradiction is explained.

  1. A destructuration theory and its application to SANICLAY model

    E-Print Network [OSTI]

    Taiebat, Mahdi; Dafalias, Yannis; Peek, Ralf

    2010-01-01

    AND ANALYTICAL METHODS IN GEOMECHANICS Int. J. Numer. Anal.a¨ a Numerical Models in Geomechanics, Proceedings of NUMOGAnalytical Methods in Geomechanics 2006; 30(12):1231–1257.

  2. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids (2012)

    E-Print Network [OSTI]

    Jameson, Antony

    2012-01-01

    dynamics, transitional or massively detached flows, turbulent mixing and combustion, or aerodynamic noise. CASTONGUAY AND A. JAMESON numerical methods available within the framework of industrial CFD applications

  3. 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.

  4. Numerical Simulation on Laser Fusion in China

    SciTech Connect (OSTI)

    Zhu Shaoping; Pei Wenbing; Xu Yan; Gu Peijun; Lan Ke; Ye Wenhua; Wu Junfeng; Li Jinghong; Gao Yaoming; Zheng Chunyang; Li Shuanggui; Mo Zeyao; Yan Jun [Institute of Applied Physics and Computational Mathematics, Beijing 100088 (China); Zhang Weiyan [National High-Tech Inertial Confinement Fusion Committee of China, Beijing 100088 (China)

    2009-05-02

    Numerical simulation is a powerful tool to get insight into the physics of laser fusion. Much effort has been devoted to develop the numerical simulation code series named LARED in China. The code series LARED are composed of six parts and enable us to have the simulation capability for the key processes in laser fusion. In recent years, a number of numerical simulations using LARED have been carried out and the simulation is checked by experiments done at the laser facility SG-II and SG-III prototype. In the present talk, some details of LARED code series will be introduced, and some simulation results, especially recent work on the opacities, will be shown.

  5. Numerical Modeling of Hydraulic Fracturing in Oil Sands

    E-Print Network [OSTI]

    2008-11-16

    as follows. First, a neat uid, such as water (called. 1. ... Department of Civil and Environmental Engineering,. University ... thick deposits in deep and semi-deep underground layers. ... hydrofracturing, in which hot water/steam is injected into the ...

  6. A Numeric Predictive Failure Model for Percussive Excavation

    E-Print Network [OSTI]

    Green, Alex Nicholas

    2011-01-01

    Journal of the Soil mechanics and Foundations Division, ASCEConference on Soil Mechanics and Foundation Engineering,Journal of the Soil Mechanics and Foundations Division,

  7. A UNIFIED NUMERICAL MODEL FOR SATURATED-UNSATURATED GROUNDWATER FLOW

    E-Print Network [OSTI]

    Narasimhan, T.N.

    2011-01-01

    between the disciplines of soil mechanics, soil physics andbetween the fields mechanics, soil physics and hydrogeology.fields of soil mechanics, hydrogeology,and soil physics. In

  8. NUMERICAL MODEL FOR LAND SUBSIDENCE IN SHALLOW GROUNDWATER SYSTEMS

    E-Print Network [OSTI]

    Narasimhan, T.N.

    2010-01-01

    1948, Fundamentals of Soil Mechanics, John Wiley and Sons,from the field of soil mechanics. A doubly-draining claySolutions for Soil and Rock Mechanics, John Wiley and Sons,

  9. Numerical issues in plasticity models for granular flows

    E-Print Network [OSTI]

    has a long history, largely due to its connection to soil mechan- ics in civil and military 27695-8205, USA Abstract Friction plays a fundamental role in the Mechanics of granular materials. Two material mechanics also plays a fun- damental role in countless geophysical phenomena. While most

  10. Rate-dependent interface models: formulation and numerical applications

    E-Print Network [OSTI]

    Corigliano, Alberto

    degradation in polymer matrix composites, the most meaningful example of application being the simulation may, in the ®rst approximation, be conducted using conventional elements to represent the composite, the phenomenon of delamination in polymer matrix composites is concerned, the viscous properties of the polymer

  11. Numerical models of Rayleigh-Taylor instabilities superimposed upon convection

    E-Print Network [OSTI]

    Schmeling, Harro

    one R.T.-overturn to multiple convective with salt tectonics (R.T.), plutonism (R.T.), plate overturns. Secondly, how is the temperature field tectonics and mantle flows (convection) (see e.g. affected

  12. Numerical modeling of methane venting from lake sediments

    E-Print Network [OSTI]

    Scandella, Benjamin P. (Benjamin Paul)

    2010-01-01

    The dynamics of methane transport in lake sediments control the release of methane into the water column above, and the portion that reaches the atmosphere may contribute significantly to the greenhouse effect. The observed ...

  13. NUMERICAL MODEL FOR LAND SUBSIDENCE IN SHALLOW GROUNDWATER SYSTEMS

    E-Print Network [OSTI]

    Narasimhan, T.N.

    2010-01-01

    and R. L. Klausing, 1969, Land subsidence due to groundwater7612-10874 Fig. S. Land subsidence at Pixley, California:Symposium on Land Subsidence, Anaheim, CA, December 10-

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

    Open Energy Info (EERE)

    form View source History View New Pages Recent Changes All Special Pages Semantic SearchQuerying Get Involved Help Apps Datasets Community Login | Sign Up Search Page Edit with...

  15. 1-Dimensional Numerical Model of Thermal Conduction and Vapor Diffusion

    E-Print Network [OSTI]

    Schörghofer, Norbert

    developed by Samar Khatiwala, 2001 extended to variable thermal properties and irregular grid by Norbert Sch for c. Upper boundary condition: a) Radiation Q + k T z z=0 = T4 z=0 Q is the incoming solar flux of Water Vapor with Phase Transitions developed by Norbert Sch¨orghofer, 2003­2004 3 phases: vapor, free

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

    Open Energy Info (EERE)

    reservoir pressure much lower than hydrostatic. This is a consequence of its natural evolution from an initial liquid-dominated to the current steam-dominated system. Beneath a...

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

    Open Energy Info (EERE)

    erosion becomes more effective for substrates that are increasingly unconsolidated and water rich. For thermo-mechanical erosion to be responsible for the formation of the...

  18. Numerical Modeling of Hydraulic Fracturing in Oil Sands

    E-Print Network [OSTI]

    2008-11-16

    Hydraulic fracturing is a widely used and e cient technique for enhancing oil extraction from heavy oil sands ..... phenomenon are the main issues involved in hydraulic fracturing. ..... energy ux due to conduction and convection: Lei = @T. @xi.

  19. Mathematical Theory and Numerical Methods for Biomolecular Modeling /

    E-Print Network [OSTI]

    White, Michael Robert

    2013-01-01

    J. A. de la Torre and P. Espa˜ nol. Coarse-graining BrownianPhys. , 69(4), 1978. [27] P. Espa˜ nol. Einstein–Helfand80:061113, 2009. [28] P. Espa˜ nol and I. Z´ u n ˜ iga. On

  20. Numerical and analytical modeling of sanding onset prediction 

    E-Print Network [OSTI]

    Yi, Xianjie

    2004-09-30

    To provide technical support for sand control decision-making, it is necessary to predict the production condition at which sand production occurs. Sanding onset prediction involves simulating the stress state on the surface of an oil/gas producing...

  1. NUMERICAL MODELING OF TURBULENT FLOW IN A COMBUSTION TUNNEL

    E-Print Network [OSTI]

    Ghoniem, A.F.

    2013-01-01

    of turbulent combustion in a 1 ean propane~a1 of high sp r mturbulent combustion stabilized behind a step in a propane~combustion behind a step at R while Su corresponding to a propane-

  2. Discussions and Closures Discussion of "Numerical Modeling of

    E-Print Network [OSTI]

    Horvath, John S.

    similar to other geosynthetic product categories such as geomembrane and geotextile. As with other types

  3. Numerical modelling of current transfer in nonlinear anisotropic conductive media

    E-Print Network [OSTI]

    Baranowski, Robert Paul

    on the nature of current transport. The main motivation for this work was the desire for a better understanding of the conceptually difficult behaviour of current transport in superconducting bodies and examines current transfer quantitatively for a number...

  4. Numerical modeling of magnetohydrodynamic activity in the Swarthmore Spheromak Experiment

    E-Print Network [OSTI]

    Brown, Michael R.

    of Physics. DOI: 10.1063/1.1362294 I. OVERVIEW There are several unresolved mysteries related to solar and astrophysical magnetism. First is the origin of solar and astrophysical magnetic fields: the dynamo problem Alfve´nic outflow jets from a solar flare reconnection site have been observed using a Doppler technique

  5. A UNIFIED NUMERICAL MODEL FOR SATURATED-UNSATURATED GROUNDWATER FLOW

    E-Print Network [OSTI]

    Narasimhan, T.N.

    2011-01-01

    A.W. , 11 The measurement of pore pressure in the triaxialtest 11 9 in Pore Pressure and Suction in Soils,9 1973 A.W. , 11 The pore pressure coefficients, A and B

  6. NUMERICAL MODEL FOR LAND SUBSIDENCE IN SHALLOW GROUNDWATER SYSTEMS

    E-Print Network [OSTI]

    Narasimhan, T.N.

    2010-01-01

    loading: of ultimate pore-pressures generated. contours LOOdcoefficient relating pore pressure change to X' effectiveof 3.05 m of excess pore-pressure within each-bed. The clay

  7. Modeling and numerical approximation of two-phase ...

    E-Print Network [OSTI]

    2011-02-09

    Feb 9, 2011 ... separated by an interface, these methods represent the interfacial tension as a body-force or bulk-stress spreading over a narrow region ...

  8. Geology 775 NUMERICAL MODELING FOR EARTH SCIENCE APPLICATIONS

    E-Print Network [OSTI]

    and engineers. Brian D. Hahn, 2002. [On reserve] Matlab, Second Edition: A Practical Introduction to Programming with Applications. Amos Gilat, 2010. [The 2008 edition is pretty reasonable on Amazon] Course Goals and Learning to address them 3. Write computer programs that use finite difference and finite element methods to solve

  9. Geology 575 NUMERICAL MODELING FOR EARTH SCIENCE APPLICATIONS

    E-Print Network [OSTI]

    and engineers. Brian D. Hahn, 2002. [On reserve] Matlab, Second Edition: A Practical Introduction to Programming with Applications. Amos Gilat, 2010. [The 2008 edition is pretty reasonable on Amazon] Course Goals and Learning to address them 3. Write computer programs that use finite difference and finite element methods to solve

  10. A UNIFIED NUMERICAL MODEL FOR SATURATED-UNSATURATED GROUNDWATER FLOW

    E-Print Network [OSTI]

    Narasimhan, T.N.

    2011-01-01

    Saturated-Unsaturated Groundwater Flow Ph.D. Dissertation in~ " Fundamental principles of groundwater flow uv e in Flowunsaturated flow in a groundwater basi.n 11 9 Hater

  11. Numerical modeling of fluid flow and time-lapse seismograms ...

    E-Print Network [OSTI]

    gabriela

    is separated from natural gas produced and is currently being injected into the Utsira Sand, a saline aquifer. • Injection started in 1996 and is planned to ...

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

    Office of Scientific and Technical Information (OSTI)

    for presentation at the Geological Society of America, October 22-25, 2006, Philadelphia, PA. Research Org: Sandia National Laboratories Sponsoring Org: USDOE Country of...

  13. 1 NUMERICAL MODELS OF SALT MARSH 2 EVOLUTION: ECOLOGICAL, GEOMORPHIC,

    E-Print Network [OSTI]

    Mudd, Simon Marius

    ]. They typically form in sheltered 45environments where fine sediments can accumulate, such 46as in estuaries et al., 2004a; Kirwan and Murray, 1 Department of Earth Sciences and Marine Program, Boston

  14. Mathematical Modelling and Numerical Simulation of Oil Pollution Problems

    E-Print Network [OSTI]

    due to convective and storm energy, and by major currents that include the Loop Current and its associ- ated rings. Shelf processes are also strongly influenced by seasonal river discharge, winds, and storms from active and passive radars, as well as from visible and infrared sensors were used to de- termine

  15. Numerical modeling of hydraulic fracture initiation and development

    E-Print Network [OSTI]

    2007-05-25

    Studying initiation and propagation of hydraulic fractures is carried out based on the hypersingular ... and for geothermal energy ...... means of enhancing oil and gas recovery, and it serves its purpose if propagates for tens or hundreds meters.

  16. Performance monitoring and numerical modelling of a deep circular excavation

    E-Print Network [OSTI]

    Schwamb, Tina

    2014-07-01

    : King post wall: Vertical posts are driven into the ground prior to excavation. Horizontal timber sheets or pre-cast concrete elements are then placed between the post with increasing excavation depth and the posts are anchored into the ground behind... sinking: Rings of pre-cast concrete segments are erected and driven into the ground with a steel cutting edge installed at the lowest ring (Figure 2.1c). A lubricant (usually bentonite) is filled into the annulus behind the wall, a calculated jacking force...

  17. Numerical models of the geodynamo and observational constraints

    E-Print Network [OSTI]

    Dormy, Emmanuel

    , and paleomagnetic observations have been interpreted as actual geomagnetic features and used as constraints. Keywords: Geodynamo; Earth core; geomagnetic field; magnetic records. Index terms: Dynamo theory; time much more recently. Papers coauthored by dynamo specialists and geophysicists observing the field over

  18. Realistic Numerical Modeling of Solar Magnetoconvection and Oscillations

    E-Print Network [OSTI]

    Kitiashvili, Irina; Kosovichev, Alexander; Wray, Alan; Mansour, Nagi

    2009-01-01

    We have developed 3D, compressible, non-linear radiative MHD simulations to study the influence of the magnetic field of various strength and geometry on the turbulent convective cells and on the excitation mechanisms of the acoustic oscillations. The results reveal substantial changes of the granulation structure with increased magnetic field, and a frequency-dependent reduction in the oscillation power. These simulation results reproduce the enhanced high-frequency acoustic emission observed at the boundaries of active region ("acoustic halo" phenomenon). In the presence of inclined magnetic field the solar convection develops filamentary structure with flows concentrated along magnetic filaments, and also exhibits behavior of running magnetoconvective waves, resembling recent observations of the sunspot penumbra dynamics from Hinode/SOT.

  19. University of Reading A Nonhydrostatic Numerical Model in -coordinates

    E-Print Network [OSTI]

    Xue, Ming

    scheme, to spatial resolution and to diffusion are also examined. Deep orographic convection is also variables, where the flux-corrected transport scheme (FCT) is employed. FCT eliminates the problem. Subgrid scale turbulence is parameterizd using the deformation and Richardson number dependent formulation

  20. University of Reading A Nonhydrostatic Numerical Model in scoordinates

    E-Print Network [OSTI]

    Xue, Ming

    scheme, to spatial resolution and to diffusion are also examined. Deep orographic convection is also of thermodynamical variables, where the flux­corrected transport scheme (FCT) is employed. FCT eliminates the problem. Subgrid scale turbulence is parameterizd using the deformation and Richardson number dependent formulation

  1. Numerical modelling of flows involving submerged bodies and free surfaces 

    E-Print Network [OSTI]

    Topper, Mathew Bernard Robert

    2011-06-28

    Kinetic energy extraction devices for ocean and river flows are often located in the vicinity of the fluid free surface. This differs from wind turbines where the atmosphere may be considered to extend to infinity for ...

  2. 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.

  3. Numerical Modeling of HCCI Combustion | Department of Energy

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

    Detroit, Michigan. Sponsored by the U.S. DOE's EERE FreedomCar and Fuel Partnership and 21st Century Truck Programs. 2006deeraceves.pdf More Documents & Publications High...

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

    Open Energy Info (EERE)

    fromthe valley floor) separated by a thick sequence (about4 km) of clastic sediments derived from the adjacentranges, and a relatively permeable, high angle faultthat functions...

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

    Open Energy Info (EERE)

    AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on QA:QAsource History ViewMayo, Maryland:NPI VenturesNewSt. Louis, Minnesota:Nulato, Alaska: Energy Resources5) Jump to:

  6. Numerical Modeling At Dixie Valley Geothermal Area (Benoit, 1999) | Open

    Open Energy Info (EERE)

    AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on QA:QAsource History ViewMayo, Maryland:NPI VenturesNewSt. Louis, Minnesota:Nulato, Alaska: Energy Resources5) Jump

  7. Numerical Modeling At Dixie Valley Geothermal Area (Iovenitti, Et Al.,

    Open Energy Info (EERE)

    AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on QA:QAsource History ViewMayo, Maryland:NPI VenturesNewSt. Louis, Minnesota:Nulato, Alaska: Energy Resources5) Jump2013)

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

    Open Energy Info (EERE)

    AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on QA:QAsource History ViewMayo, Maryland:NPI VenturesNewSt. Louis, Minnesota:Nulato, Alaska: Energy

  9. Numerical Modeling of Transient Basin and Range Extensional Geothermal

    Open Energy Info (EERE)

    AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on QA:QAsource History ViewMayo, Maryland:NPI VenturesNewSt. Louis, Minnesota:Nulato, Alaska: EnergySystems | Open Energy

  10. Numerical Modelling of Geothermal Systems a Short Introduction | Open

    Open Energy Info (EERE)

    AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on QA:QAsource History ViewMayo, Maryland:NPI VenturesNewSt. Louis, Minnesota:Nulato, Alaska: EnergySystems | Open

  11. Analytical-Numerical Modeling Of Komatiite Lava Emplacement And Thermal

    Open Energy Info (EERE)

    AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on QA:QA J-E-1 SECTION JEnvironmental Jump to:EAandAmminex A S Jump to: navigation,Inof Ground Source Heat| Open1986)

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

    Open Energy Info (EERE)

    AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on QA:QA J-E-1 SECTION J APPENDIXsourceII Jump to:Information 3rd|Northfork ElectricName Nueces Electric7) Jump to:

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

    Open Energy Info (EERE)

    AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on QA:QA J-E-1 SECTION J APPENDIXsourceII Jump to:Information 3rd|Northfork ElectricName Nueces Electric7) Jump to:9)

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

    Open Energy Info (EERE)

    AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on QA:QA J-E-1 SECTION J APPENDIXsourceII Jump to:Information 3rd|Northfork ElectricName Nueces Electric7) Jump

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

    Open Energy Info (EERE)

    AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on QA:QA J-E-1 SECTION J APPENDIXsourceII Jump to:Information 3rd|Northfork ElectricName Nueces Electric7) Jump6) Jump

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

    Open Energy Info (EERE)

    AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on QA:QA J-E-1 SECTION J APPENDIXsourceII Jump to:Information 3rd|Northfork ElectricName Nueces Electric7) Jump6)

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

    Open Energy Info (EERE)

    AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on QA:QA J-E-1 SECTION J APPENDIXsourceII Jump to:Information 3rd|Northfork ElectricName Nueces Electric7) Jump6)Coso

  18. Numerical Modeling Of Basin And Range Geothermal Systems | Open Energy

    Open Energy Info (EERE)

    AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on QA:QA J-E-1 SECTION J APPENDIXsourceII Jump to:Information 3rd|Northfork ElectricName Nueces Electric7)

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

    Open Energy Info (EERE)

    AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on QA:QA J-E-1 SECTION J APPENDIXsourceII Jump to:Information 3rd|Northfork ElectricName Nueces Electric7)Fractures

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

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

    AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Homesum_a_epg0_fpd_mmcf_m.xls" ,"Available from WebQuantityBonneville Power Administration would likeUniverse (JournalvivoHighHussein Khalil Hussein KhalilStatistical Self-Similarity