Modeling-Computer Simulations At Northern Basin & Range Region...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Northern Basin & Range Region (Pritchett, 2004) Exploration Activity...
Modeling-Computer Simulations At Central Nevada Seismic Zone...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Central Nevada Seismic Zone Region (Pritchett, 2004) Exploration...
Modeling-Computer Simulations At Geysers Area (Goff & Decker...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Geysers Area (Goff & Decker, 1983) Exploration Activity Details...
Modeling-Computer Simulations At Dixie Valley Geothermal Area...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Dixie Valley Geothermal Area (Wisian & Blackwell, 2004) Exploration...
Modeling-Computer Simulations At Raft River Geothermal Area ...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Raft River Geothermal Area (1980) Exploration Activity Details...
Modeling-Computer Simulations (Lewicki & Oldenburg, 2004) | Open...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations (Lewicki & Oldenburg, 2004) Exploration Activity Details Location...
Modeling-Computer Simulations At Desert Peak Area (Wisian & Blackwell...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Desert Peak Area (Wisian & Blackwell, 2004) Exploration Activity...
Modeling-Computer Simulations (Combs, Et Al., 1999) | Open Energy...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations (Combs, Et Al., 1999) Exploration Activity Details Location Unspecified...
Modeling-Computer Simulations At Yellowstone Region (Laney, 2005...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Yellowstone Region (Laney, 2005) Exploration Activity Details Location...
Modeling-Computer Simulations At Raft River Geothermal Area ...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Raft River Geothermal Area (1979) Exploration Activity Details...
Modeling-Computer Simulations At Raft River Geothermal Area ...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Raft River Geothermal Area (1977) Exploration Activity Details...
Modeling-Computer Simulations (Ozkocak, 1985) | Open Energy Informatio...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations (Ozkocak, 1985) Exploration Activity Details Location Unspecified...
Modeling-Computer Simulations At White Mountains Area (Goff ...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At White Mountains Area (Goff & Decker, 1983) Exploration Activity...
Modeling-Computer Simulations At Stillwater Area (Wisian & Blackwell...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Stillwater Area (Wisian & Blackwell, 2004) Exploration Activity...
Modeling-Computer Simulations At Valles Caldera - Redondo Geothermal...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Valles Caldera - Redondo Geothermal Area (Wilt & Haar, 1986)...
Modeling-Computer Simulations At Dixie Valley Geothermal Area...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Dixie Valley Geothermal Area (Kennedy & Soest, 2006) Exploration...
Modeling-Computer Simulations (Ranalli & Rybach, 2005) | Open...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations (Ranalli & Rybach, 2005) Exploration Activity Details Location...
Modeling-Computer Simulations At Raft River Geothermal Area ...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Raft River Geothermal Area (1983) Exploration Activity Details...
Modeling-Computer Simulations At Dixie Valley Geothermal Area...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Dixie Valley Geothermal Area (Wannamaker, Et Al., 2006) Exploration...
Modeling-Computer Simulations At Obsidian Cliff Area (Hulen,...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Obsidian Cliff Area (Hulen, Et Al., 2003) Exploration Activity Details...
Modeling-Computer Simulations At Walker-Lane Transitional Zone...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Walker-Lane Transitional Zone Region (Laney, 2005) Exploration...
Modeling-Computer Simulations At Valles Caldera - Redondo Geothermal...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Valles Caldera - Redondo Geothermal Area (Roberts, Et Al., 1995)...
Modeling-Computer Simulations At Long Valley Caldera Geothermal...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Long Valley Caldera Geothermal Area (Pribnow, Et Al., 2003)...
Modeling-Computer Simulations At Hawthorne Area (Lazaro, Et Al...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Hawthorne Area (Lazaro, Et Al., 2010) Exploration Activity Details...
Modeling-Computer Simulations At Walker-Lane Transitional Zone...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Walker-Lane Transitional Zone Region (Pritchett, 2004) Exploration...
Modeling-Computer Simulations At Fenton Hill HDR Geothermal Area...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Fenton Hill HDR Geothermal Area (Brown & DuTeaux, 1997) Exploration...
Modeling-Computer Simulations At Coso Geothermal Area (1980)...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Coso Geothermal Area (1980) Exploration Activity Details Location Coso...
Modeling-Computer Simulations At Long Valley Caldera Geothermal...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Long Valley Caldera Geothermal Area (Newman, Et Al., 2006) Exploration...
Modeling-Computer Simulations At The Needles Area (Bell & Ramelli...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At The Needles Area (Bell & Ramelli, 2009) Exploration Activity Details...
Modeling-Computer Simulations At Fenton Hill HDR Geothermal Area...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Fenton Hill HDR Geothermal Area (Goff & Decker, 1983) Exploration...
Modeling-Computer Simulations At Long Valley Caldera Geothermal...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Long Valley Caldera Geothermal Area (Farrar, Et Al., 2003) Exploration...
Modeling-Computer Simulations At Central Nevada Seismic Zone...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Central Nevada Seismic Zone Region (Biasi, Et Al., 2009) Exploration...
Modeling-Computer Simulations At Valles Caldera - Sulphur Springs...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Valles Caldera - Sulphur Springs Geothermal Area (Roberts, Et Al.,...
Modeling-Computer Simulations At Nw Basin & Range Region (Pritchett...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Nw Basin & Range Region (Pritchett, 2004) Exploration Activity Details...
Modeling-Computer Simulations At Long Valley Caldera Geothermal...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Long Valley Caldera Geothermal Area (Tempel, Et Al., 2011) Exploration...
Modeling-Computer Simulations At Nw Basin & Range Region (Biasi...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Nw Basin & Range Region (Biasi, Et Al., 2009) Exploration Activity...
Modeling-Computer Simulations At Coso Geothermal Area (2000)...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Coso Geothermal Area (2000) Exploration Activity Details Location Coso...
Modeling-Computer Simulations At Northern Basin & Range Region...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Northern Basin & Range Region (Biasi, Et Al., 2009) Exploration...
Modeling-Computer Simulations At Valles Caldera - Sulphur Springs...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Valles Caldera - Sulphur Springs Geothermal Area (Wilt & Haar, 1986)...
Modeling-Computer Simulations At Akutan Fumaroles Area (Kolker...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Akutan Fumaroles Area (Kolker, Et Al., 2010) Exploration Activity...
Modeling-Computer Simulations At Walker-Lane Transitional Zone...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Walker-Lane Transitional Zone Region (Biasi, Et Al., 2009) Exploration...
Modeling-Computer Simulations At Coso Geothermal Area (1999)...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Coso Geothermal Area (1999) Exploration Activity Details Location Coso...
Modeling-Computer Simulations At Fish Lake Valley Area (Deymonaz...
Fish Lake Valley Area (Deymonaz, Et Al., 2008) Jump to: navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Fish Lake Valley...
Modeling-Computer Simulations At Nevada Test And Training Range...
Nevada Test And Training Range Area (Sabin, Et Al., 2004) Jump to: navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Nevada...
Modeling-Computer Simulations | Open Energy Information
the risk of inaccurate predictions.1 Potential Pitfalls Uncertainties in initial reservoir conditions and other model inputs can cause inaccuracies in simulations, which...
Modeling-Computer Simulations At U.S. West Region (Sabin, Et...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At U.S. West Region (Sabin, Et Al., 2004) Exploration Activity Details...
Modeling-Computer Simulations At Cove Fort Area (Toksoz, Et Al...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Cove Fort Area (Toksoz, Et Al, 2010) Exploration Activity Details...
Modeling-Computer Simulations At U.S. West Region (Williams ...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At U.S. West Region (Williams & Deangelo, 2008) Exploration Activity...
Modeling-Computer Simulations At Chocolate Mountains Area (Alm...
Simulations Activity Date Usefulness useful DOE-funding Unknown Notes "Shallow temperature gradient drilling began at the CMAGR in January of 2010. 13 temperature...
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Direct Numerical Simulation - Sandia Energy Energy Search Icon Sandia Home Locations Contact Us Employee Locator Energy & Climate Secure & Sustainable Energy Future Stationary Power Energy Conversion Efficiency Solar Energy Wind Energy Water Power Supercritical CO2 Geothermal Natural Gas Safety, Security & Resilience of the Energy Infrastructure Energy Storage Nuclear Power & Engineering Grid Modernization Battery Testing Nuclear Energy Defense Waste Management Programs Advanced
Confidence in Numerical Simulations
Hemez, Francois M.
2015-02-23
This PowerPoint presentation offers a high-level discussion of uncertainty, confidence and credibility in scientific Modeling and Simulation (M&S). It begins by briefly evoking M&S trends in computational physics and engineering. The first thrust of the discussion is to emphasize that the role of M&S in decision-making is either to support reasoning by similarity or to “forecast,” that is, make predictions about the future or extrapolate to settings or environments that cannot be tested experimentally. The second thrust is to explain that M&S-aided decision-making is an exercise in uncertainty management. The three broad classes of uncertainty in computational physics and engineering are variability and randomness, numerical uncertainty and model-form uncertainty. The last part of the discussion addresses how scientists “think.” This thought process parallels the scientific method where by a hypothesis is formulated, often accompanied by simplifying assumptions, then, physical experiments and numerical simulations are performed to confirm or reject the hypothesis. “Confidence” derives, not just from the levels of training and experience of analysts, but also from the rigor with which these assessments are performed, documented and peer-reviewed.
Sandia Energy - Numerical Simulations of Hydrokinetics in the...
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Numerical Simulations of Hydrokinetics in the Roza Canal, Yakima Washington Home Renewable Energy Energy Water Power Computational Modeling & Simulation Numerical Simulations of...
Direct Numerical Simulation of Compressible, Turbulent Flow ...
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
The computational mesh for this direct numerical simulation was over 33 billion cells, and was run on up to 102,400 cores under a DoD HPCMP Frontier Project. Nicholas Bisek and ...
Numerical simulations of strong incompressible magnetohydrodynamic turbulence
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.
Development of Numerical Simulation Capabilities for In Situ...
Office of Scientific and Technical Information (OSTI)
Development of Numerical Simulation Capabilities for In Situ Heating of Oil Shale Citation Details In-Document Search Title: Development of Numerical Simulation Capabilities for In ...
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 ...
Numerical recipes for mold filling simulation
Kothe, D.; Juric, D.; Lam, K.; Lally, B.
1998-07-01
Has the ability to simulate the filling of a mold progressed to a point where an appropriate numerical recipe achieves the desired results? If results are defined to be topological robustness, computational efficiency, quantitative accuracy, and predictability, all within a computational domain that faithfully represents complex three-dimensional foundry molds, then the answer unfortunately remains no. Significant interfacial flow algorithm developments have occurred over the last decade, however, that could bring this answer closer to maybe. These developments have been both evolutionary and revolutionary, will continue to transpire for the near future. Might they become useful numerical recipes for mold filling simulations? Quite possibly. Recent progress in algorithms for interface kinematics and dynamics, linear solution methods, computer science issues such as parallelization and object-oriented programming, high resolution Navier-Stokes (NS) solution methods, and unstructured mesh techniques, must all be pursued as possible paths toward higher fidelity mold filling simulations. A detailed exposition of these algorithmic developments is beyond the scope of this paper, hence the authors choose to focus here exclusively on algorithms for interface kinematics. These interface tracking algorithms are designed to model the movement of interfaces relative to a reference frame such as a fixed mesh. Current interface tracking algorithm choices are numerous, so is any one best suited for mold filling simulation? Although a clear winner is not (yet) apparent, pros and cons are given in the following brief, critical review. Highlighted are those outstanding interface tracking algorithm issues the authors feel can hamper the reliable modeling of today`s foundry mold filling processes.
Direct numerical simulation of turbulent reacting flows
Chen, J.H.
1993-12-01
The development of turbulent combustion models that reflect some of the most important characteristics of turbulent reacting flows requires knowledge about the behavior of key quantities in well defined combustion regimes. In turbulent flames, the coupling between the turbulence and the chemistry is so strong in certain regimes that is is very difficult to isolate the role played by one individual phenomenon. Direct numerical simulation (DNS) is an extremely useful tool to study in detail the turbulence-chemistry interactions in certain well defined regimes. Globally, non-premixed flames are controlled by two limiting cases: the fast chemistry limit, where the turbulent fluctuations. In between these two limits, finite-rate chemical effects are important and the turbulence interacts strongly with the chemical processes. This regime is important because industrial burners operate in regimes in which, locally the flame undergoes extinction, or is at least in some nonequilibrium condition. Furthermore, these nonequilibrium conditions strongly influence the production of pollutants. To quantify the finite-rate chemistry effect, direct numerical simulations are performed to study the interaction between an initially laminar non-premixed flame and a three-dimensional field of homogeneous isotropic decaying turbulence. Emphasis is placed on the dynamics of extinction and on transient effects on the fine scale mixing process. Differential molecular diffusion among species is also examined with this approach, both for nonreacting and reacting situations. To address the problem of large-scale mixing and to examine the effects of mean shear, efforts are underway to perform large eddy simulations of round three-dimensional jets.
Numerical simulation of Rayleigh-Taylor instabilities involving solids
Chang, Chong H.
2015-11-20
This report is a description of research performed by LANL regarding numeric simulations of Rayleigh-Taylor instability.
Numerical simulations of capillary barrier field tests
Morris, C.E.; Stormont, J.C.
1997-12-31
Numerical simulations of two capillary barrier systems tested in the field were conducted to determine if an unsaturated flow model could accurately represent the observed results. The field data was collected from two 7-m long, 1.2-m thick capillary barriers built on a 10% grade that were being tested to investigate their ability to laterally divert water downslope. One system had a homogeneous fine layer, while the fine soil of the second barrier was layered to increase its ability to laterally divert infiltrating moisture. The barriers were subjected first to constant infiltration while minimizing evaporative losses and then were exposed to ambient conditions. The continuous infiltration period of the field tests for the two barrier systems was modelled to determine the ability of an existing code to accurately represent capillary barrier behavior embodied in these two designs. Differences between the field test and the model data were found, but in general the simulations appeared to adequately reproduce the response of the test systems. Accounting for moisture retention hysteresis in the layered system will potentially lead to more accurate modelling results and is likely to be important when developing reasonable predictions of capillary barrier behavior.
Development of Numerical Simulation Capabilities for In Situ...
Office of Scientific and Technical Information (OSTI)
Numerical Simulation Capabilities for In Situ Heating of Oil Shale Hoda, Nazish ExxonMobil Upstream Research Company, Houston, TX, USA; Fang, Chen ExxonMobil Upstream Research...
Development of Numerical Simulation Capabilities for In Situ...
Office of Scientific and Technical Information (OSTI)
for In Situ Heating of Oil Shale Citation Details In-Document Search Title: Development of Numerical Simulation Capabilities for In Situ Heating of Oil Shale Authors: Hoda, ...
Info-Gap Analysis of Truncation Errors in Numerical Simulations...
Office of Scientific and Technical Information (OSTI)
Title: Info-Gap Analysis of Truncation Errors in Numerical Simulations. Authors: Kamm, James R. ; Witkowski, Walter R. ; Rider, William J. ; Trucano, Timothy Guy ; Ben-Haim, Yakov. ...
Numerical Simulation of Ni Grain Growth in a Thermal Gradient
Office of Scientific and Technical Information (OSTI)
665C Numerical Simulation of Ni Grain Growth in a Thermal Gradient Sandia National Laboratories John A. Mitchell and Veena Tikare Sandia National Laboratories, Albuquerque New ...
Direct Numerical Simulations and Robust Predictions of Cloud...
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
cloud. Credit: Computational Science and Engineering Laboratory, ETH Zurich, Switzerland Direct Numerical Simulations and Robust Predictions of Cloud Cavitation Collapse PI Name:...
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 ... Visit OSTI to utilize additional information resources in energy science and technology. A ...
A review of recent advances of numerical simulations of microscale...
Office of Scientific and Technical Information (OSTI)
Due to the reformer's small size, numerical simulations are critical to understand heat and mass transfer phenomena occurring in the systems. This paper reviews the development of ...
Numerical simulation of the environmental impact of hydraulic...
Office of Scientific and Technical Information (OSTI)
Numerical simulation of the environmental impact of hydraulic fracturing of tightshale gas reservoirs on near-surface groundwater: Background, base cases, shallow reservoirs,...
Numerical simulations shed new light on early universe
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Numerical simulations shed new light on early universe Numerical simulations shed new light on early universe The code simulates conditions during the first few minutes of cosmological evolution to model the role of neutrinos, nuclei and other particles in shaping the early universe. April 21, 2016 Los Alamos scientists developed the BURST computer code to predict-to unprecedented precision-the amounts of light nuclei synthesized in the Big Bang. Los Alamos scientists developed the BURST
Numerical simulation experiments on the long-term evolution of...
Office of Scientific and Technical Information (OSTI)
the long-term evolution of a CO2 plume under a sloping caprock Citation Details In-Document Search Title: Numerical simulation experiments on the long-term evolution of a CO2 plume ...
Numerical and laboratory simulations of auroral acceleration
Gunell, H.; De Keyser, J.; Mann, I.
2013-10-15
The existence of parallel electric fields is an essential ingredient of auroral physics, leading to the acceleration of particles that give rise to the auroral displays. An auroral flux tube is modelled using electrostatic Vlasov simulations, and the results are compared to simulations of a proposed laboratory device that is meant for studies of the plasma physical processes that occur on auroral field lines. The hot magnetospheric plasma is represented by a gas discharge plasma source in the laboratory device, and the cold plasma mimicking the ionospheric plasma is generated by a Q-machine source. In both systems, double layers form with plasma density gradients concentrated on their high potential sides. The systems differ regarding the properties of ion acoustic waves that are heavily damped in the magnetosphere, where the ion population is hot, but weakly damped in the laboratory, where the discharge ions are cold. Ion waves are excited by the ion beam that is created by acceleration in the double layer in both systems. The efficiency of this beam-plasma interaction depends on the acceleration voltage. For voltages where the interaction is less efficient, the laboratory experiment is more space-like.
Numerical simulation model for vertical flow in geothermal wells
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.
Direct Numerical Simulations and Robust Predictions of Cloud Cavitation
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Collapse | Argonne Leadership Computing Facility cavitating vapor bubbles above a solid wall This image shows cavitating vapor bubbles above a solid wall. The yellow is a visualization of the pressure peak in the center of the bubble cloud. Credit: Computational Science and Engineering Laboratory, ETH Zurich, Switzerland Direct Numerical Simulations and Robust Predictions of Cloud Cavitation Collapse PI Name: Petros Koumoutsakos PI Email: petros@ethz.ch Institution: ETH Zürich Allocation
Direct Numerical Simulations and Robust Predictions of Cloud Cavitation
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Collapse | Argonne Leadership Computing Facility Initiation of cloud cavitation collapse for 50,000 bubbles Initiation of cloud cavitation collapse for 50,000 bubbles. Jonas Sukys, ETH Zurich Direct Numerical Simulations and Robust Predictions of Cloud Cavitation Collapse PI Name: Petros Koumoutsakos PI Email: petros@ethz.ch Institution: ETH Zurich Allocation Program: INCITE Allocation Hours at ALCF: 72 Million Year: 2016 Research Domain: Engineering Cloud cavitation collapse-the evolution
Direct Numerical Simulations of High Reynolds Number Turbulent Channel Flow
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
| Argonne Leadership Computing Facility Visualization of the spanwise vorticity in a turbulent channel. S. Hoyas and O. Flores while they were at Universidad Politecnica de Madrid Direct Numerical Simulations of High Reynolds Number Turbulent Channel Flow PI Name: Robert Moser PI Email: rmoser@ices.utexas.edu Institution: University of Texas Allocation Program: INCITE Allocation Hours at ALCF: 175 Million Year: 2013 Research Domain: Engineering Approximately 28% of U.S. energy resources are
Numerical simulations of the decay of primordial magnetic turbulence
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.
Numerical Simulations of Small Non-spherical Particles in Turbulence |
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Argonne Leadership Computing Facility Numerical Simulations of Small Non-spherical Particles in Turbulence Event Sponsor: Mathematics and Computer Science Division LANS Seminar Start Date: Aug 31 2016 - 3:00pm Building/Room: Building 240/Room 1404-1405 Location: Argonne National Laboratory Speaker(s): Nimish Pujara Speaker(s) Title: UC Berkeley Motivated by the ubiquity of natural particles in turbulent flows in the natural environment as well as in many industrial processes, we investigate
Transient productivity index for numerical well test simulations
Blanc, G.; Ding, D.Y.; Ene, A.
1997-08-01
The most difficult aspect of numerical simulation of well tests is the treatment of the Bottom Hole Flowing (BHF) Pressure. In full field simulations, this pressure is derived from the Well-block Pressure (WBP) using a numerical productivity index which accounts for the grid size and permeability, and for the well completion. This productivity index is calculated assuming a pseudo-steady state flow regime in the vicinity of the well and is therefore constant during the well production period. Such a pseudo-steady state assumption is no longer valid for the early time of a well test simulation as long as the pressure perturbation has not reached several grid-blocks around the well. This paper offers two different solutions to this problem: (1) The first one is based on the derivation of a Numerical Transient Productivity Index (NTPI) to be applied to Cartesian grids; (2) The second one is based on the use of a Corrected Transmissibility and Accumulation Term (CTAT) in the flow equation. The representation of the pressure behavior given by both solutions is far more accurate than the conventional one as shown by several validation examples which are presented in the following pages.
Numerical simulation of steam injection processes with solvent
Zerpa, L.; Mendez, Z.
1995-12-31
In Venezuela during recent years, gas oil has been evaluated as an additive to increase steam injection process efficiency. The results of laboratory and field tests have shown a significant improvement in the production behavior. Despite these experiences, it is necessary to complement the information with results obtained from numerical simulation studies in order to know injection parameter effects, such as gas oil concentration, schemes and rates of injection, temperatures, etc., and also some mechanisms involved in the process. In this work, the results achieved in the numerical simulation of displacement tests with steam and gas oil are presented. A fully implicit 2-D thermal, three-phase compositional simulator was used to obtain all the data presented in this paper The numerical simulation results show a similar oil production performance to those obtained in the displacement tests with injection of gas oil and steam simultaneously. These results indicate rising of the production rate when the solvent concentration increases. They also reveal that the solvent co-injection scheme improves the productivity in relation to the gas oil pre-injection at low temperature. However, when gas oil is pre-injected at higher temperature, the oil production performance is similar to the co-injection scheme performance. This can attribute to the favorable temperature effect on the diffusion mechanisms. On the other hand, an increase of the gas oil injection rate causes a productivity reduction. In addition, the gas oil capacity to remove more viscous fractions than the original crude was verified. It was determined that the gas oil light fraction volatilization contributes to the process improvement. In general, these results confirm the benefit of using solvent and contribute to the understanding of process mechanisms.
Numerical simulation of multi-layered textile composite reinforcement forming
Wang, P.; Hamila, N.; Boisse, P.
2011-05-04
One important perspective in aeronautics is to produce large, thick or/and complex structural composite parts. The forming stage presents an important role during the whole manufacturing process, especially for LCM processes (Liquid Composites Moulding) or CFRTP (Continuous Fibre Reinforcements and Thermoplastic resin). Numerical simulations corresponding to multi-layered composite forming allow the prediction for a successful process to produce the thick parts, and importantly, the positions of the fibres after forming to be known. This paper details a set of simulation examples carried out by using a semi-discrete shell finite element made up of unit woven cells. The internal virtual work is applied on all woven cells of the element taking into account tensions, in-plane shear and bending effects. As one key problem, the contact behaviours of tool/ply and ply/ply are described in the numerical model. The simulation results not only improve our understanding of the multi-layered composite forming process but also point out the importance of the fibre orientation and inter-ply friction during formability.
Diffusive mesh relaxation in ALE finite element numerical simulations
Dube, E.I.
1996-06-01
The theory for a diffusive mesh relaxation algorithm is developed for use in three-dimensional Arbitary Lagrange/Eulerian (ALE) finite element simulation techniques. This mesh relaxer is derived by a variational principle for an unstructured 3D grid using finite elements, and incorporates hourglass controls in the numerical implementation. The diffusive coefficients are based on the geometric properties of the existing mesh, and are chosen so as to allow for a smooth grid that retains the general shape of the original mesh. The diffusive mesh relaxation algorithm is then applied to an ALE code system, and results from several test cases are discussed.
Numerical simulation of carbon arc discharge for nanoparticle synthesis
Kundrapu, M.; Keidar, M.
2012-07-15
Arc discharge with catalyst-filled carbon anode in helium background was used for the synthesis of carbon nanoparticles. In this paper, we present the results of numerical simulation of carbon arc discharges with arc current varying from 10 A to 100 A in a background gas pressure of 68 kPa. Anode sublimation rate and current voltage characteristics are compared with experiments. Distribution of temperature and species density, which is important for the estimation of the growth of nanoparticles, is obtained. The probable location of nanoparticle growth region is identified based on the temperature range for the formation of catalyst clusters.
Numerical simulation of linear fiction welding (LFW) processes
Fratini, L.; La Spisa, D. [University of Palermo-Dept. of Industrial engineering (Italy)
2011-05-04
Solid state welding processes are becoming increasingly important due to a large number of advantages related to joining ''unweldable'' materials and in particular light weight alloys. Linear friction welding (LFW) has been used successfully to bond non-axisymmetric components of a range of materials including titanium alloys, steels, aluminum alloys, nickel, copper, and also dissimilar material combinations. The technique is useful in the research of quality of the joints and in reducing costs of components and parts of the aeronautic and automotive industries.LFW involves parts to be welded through the relative reciprocating motion of two components under an axial force. In such process the heat source is given by the frictional forces work decaying into heat determining a local softening of the material and proper bonding conditions due to both the temperature increase and the local pressure of the two edges to be welded. This paper is a comparative test between the numerical model in two dimensions, i.e. in plane strain conditions, and in three dimensions of a LFW process of AISI1045 steel specimens. It must be observed that the 3D model assures a faithful simulation of the actual threedimensional material flow, even if the two-dimensional simulation computational times are very short, a few hours instead of several ones as the 3D model. The obtained results were compared with experimental values found out in the scientific literature.
Advanced Numerical Methods and Software Approaches for Semiconductor Device Simulation
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Carey, Graham F.; Pardhanani, A. L.; Bova, S. W.
2000-01-01
In this article we concisely present several modern strategies that are applicable to driftdominated carrier transport in higher-order deterministic models such as the driftdiffusion, hydrodynamic, and quantum hydrodynamic systems. The approaches include extensions of “upwind” and artificial dissipation schemes, generalization of the traditional Scharfetter – Gummel approach, Petrov – Galerkin and streamline-upwind Petrov Galerkin (SUPG), “entropy” variables, transformations, least-squares mixed methods and other stabilized Galerkin schemes such as Galerkin least squares and discontinuous Galerkin schemes. The treatment is representative rather than an exhaustive review and several schemes are mentioned only briefly with appropriate reference to the literature. Some of themore » methods have been applied to the semiconductor device problem while others are still in the early stages of development for this class of applications. We have included numerical examples from our recent research tests with some of the methods. A second aspect of the work deals with algorithms that employ unstructured grids in conjunction with adaptive refinement strategies. The full benefits of such approaches have not yet been developed in this application area and we emphasize the need for further work on analysis, data structures and software to support adaptivity. Finally, we briefly consider some aspects of software frameworks. These include dial-an-operator approaches such as that used in the industrial simulator PROPHET, and object-oriented software support such as those in the SANDIA National Laboratory framework SIERRA.« less
DRIVER TO SUPPORT USE OF NUMERICAL SIMULATION TOOLS
Energy Science and Technology Software Center (OSTI)
2001-02-13
UNIPACK is a computer interface that simplifies and enhances the use of numerical simulation tools to design a primary geometry and/or a forming die for a powder compact and/or to design the pressing process used to shape a powder by compaction. More particularly, it is an interface that utilizes predefined generic geometric configurations to simplify the use of finite element method modeling software to simply and more efficiently design: (1) the shape and size amore » powder compact; (2) a forming die to shape a powder compact; and/or (3) the pressing process used to form a powder compact. UNIPACK is a user interface for a predictive model for powder compaction that incorporates unprecedented flexibility to design powder press tooling and powder pressing processes. UNIPACK works with the Sandia National Laboratories (SNL) Engineering Analysis Cide Access System (SEACAS) to generate a finite element (FE) mesh and automatically perform a FE analysis of powder compaction. UNIPACK was developed to allow a non-expert with minimal training to quickly and easily design/construct a variable dimension component or die in real time on a desktop or laptop personal computer.« less
Numerical simulation of the flow in wire-wrapped pin bundles: Effect of
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
pin-wire contact modeling | Argonne Leadership Computing Facility simulation of the flow in wire-wrapped pin bundles: Effect of pin-wire contact modeling Authors: Merzari, E., Smith, J.G., Tentner, A., Pointer, W.D., Fischer, P. The rapid advancement of numerical techniques and the availability of increasingly powerful supercomputers recently enabled scientists to use large eddy simulation (LES) to simulate numerically the flow in a full subassembly composed of wire-wrapped pins. Because of
Numerical Simulation of Groundwater Withdrawal at the Nevada Test Site
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.
Numerical simulations of altocumulus with a cloud resolving model
Liu, S.; Krueger, S.K.
1996-04-01
Altocumulus and altostratus clouds together cover approximately 22% of the earth`s surface. They play an important role in the earth`s energy budget through their effect on solar and infrared radiation. However, there has been little altocumulus cloud investigation by either modelers or observational programs. Starr and Cox (SC) (1985a,b) simulated an altostratus case as part of the same study in which they modeled a thin layer of cirrus. Although this calculation was originally described as representing altostratus, it probably better represents altocumulus stratiformis. In this paper, we simulate altocumulus cloud with a cloud resolving model (CRM). We simply describe the CRM first. We calculate the same middle-level cloud case as SC to compare our results with theirs. We will look at the role of cloud-scale processes in response to large-scale forcing. We will also discuss radiative effects by simulating diurnal and nocturnal cases. Finally, we discuss the utility of a 1D model by comparing 1D simulations and 2D simulations.
Numerical simulation of the hydrodynamical combustion to strange quark matter
Niebergal, Brian; Ouyed, Rachid; Jaikumar, Prashanth
2010-12-15
We present results from a numerical solution to the burning of neutron matter inside a cold neutron star into stable u,d,s quark matter. Our method solves hydrodynamical flow equations in one dimension with neutrino emission from weak equilibrating reactions, and strange quark diffusion across the burning front. We also include entropy change from heat released in forming the stable quark phase. Our numerical results suggest burning front laminar speeds of 0.002-0.04 times the speed of light, much faster than previous estimates derived using only a reactive-diffusive description. Analytic solutions to hydrodynamical jump conditions with a temperature-dependent equation of state agree very well with our numerical findings for fluid velocities. The most important effect of neutrino cooling is that the conversion front stalls at lower density (below {approx_equal}2 times saturation density). In a two-dimensional setting, such rapid speeds and neutrino cooling may allow for a flame wrinkle instability to develop, possibly leading to detonation.
Numerical thermalization in particle-in-cell simulations with Monte-Carlo collisions
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.
1980-12-01
The second quarterly technical progress report is presented for a program entitled, Application of Numerical Simulation Methodology to Automotive Combustion. The goal of the program is to develop, validate, demonstrate and apply a numerical simulation methodology for in-cylinder reactive flows in internal combustion engines. Previous work on this contract involved the initial development and validation of a finite difference based simulation model for time dependent axisymmetric flows which includes: a generalized coordinate system for arbitrary mesh design and treatment of complex and time dependent boundaries; multiple and interacting chemical species; coupled swirl flow velocity component; and two-equation turbulence closure. In its various stages of development, the model has been used to simulate numerous engine-related problems for validation and demonstration purposes. The technical effort during the current reporting period has concentrated on: reactive flow model development, test and data comparison studies; swirl flow simulations; and in-cylinder compression cycle flow simulations. Results of these studies are discussed.
Numerical simulation of optical feedback on a quantum dot lasers
Al-Khursan, Amin H.; Ghalib, Basim Abdullattif; Al-Obaidi, Sabri J.
2012-02-15
We use multi-population rate equations model to study feedback oscillations in the quantum dot laser. This model takes into account all peculiar characteristics in the quantum dots such as inhomogeneous broadening of the gain spectrum, the presence of the excited states on the quantum dot and the non-confined states due to the presence of wetting layer and the barrier. The contribution of quantum dot groups, which cannot follow by other models, is simulated. The results obtained from this model show the feedback oscillations, the periodic oscillations which evolves to chaos at higher injection current of higher feedback levels. The frequency fluctuation is attributed mainly to wetting layer with a considerable contribution from excited states. The simulation shows that is must be not using simple rate equation models to express quantum dots working at excited state transition.
Three dimensional numerical simulations of the UPS-292-SC engine
O'Rourke, P.J.; Amsden, A.A.
1987-01-01
We present and analyze three-dimensional calculations of the spray, mixing and combustion in the UPS-292 stratified charge engine for three different operating conditions, corresponding to overall air-fuel ratios between 22.4 and 61.0. The numerical calculations are performed with KIVA, a multidimensional arbitrary-mesh, finite-difference hydrodynamics program for internal combustion engine applications. The calculations use a mesh of 10,000 computational cells, which conform to the shape of the piston bowl and cylinder and move to follow piston motion. Each operating condition is calculated from intake valve closure at 118/sup 0/ BTDC to 90/sup 0/ ATDC and requires approximately three hours of CRAY-XMP computer time.
Numerical Simulations of Leakage from Underground LPG Storage Caverns
Yamamoto, Hajime; Pruess, Karsten
2004-09-01
To secure a stable supply of petroleum gas, underground storage caverns for liquified petroleum gas (LPG) are commonly used in many countries worldwide. Storing LPG in underground caverns requires that the surrounding rock mass remain saturated with groundwater and that the water pressure be higher than the liquid pressure inside the cavern. In previous studies, gas containment criteria for underground gas storage based on hydraulic gradient and pressure have been discussed, but these studies do not consider the physicochemical characteristics and behavior of LPG such as vaporization and dissolution in groundwater. Therefore, while these studies are very useful for designing storage caverns, they do not provide better understanding of the either the environmental effects of gas contamination or the behavior of vaporized LPG. In this study, we have performed three-phase fluid flow simulations of gas leakage from underground LPG storage caverns, using the multiphase multicomponent nonisothermal simulator TMVOC (Pruess and Battistelli, 2002), which is capable of solving the three-phase nonisothermal flow of water, gas, and a multicomponent mixture of volatile organic chemicals (VOCs) in multidimensional heterogeneous porous media. A two-dimensional cross-sectional model resembling an actual underground LPG facility in Japan was developed, and gas leakage phenomena were simulated for three different permeability models: (1) a homogeneous model, (2) a single-fault model, and (3) a heterogeneous model. In addition, the behavior of stored LPG was studied for the special case of a water curtain suddenly losing its function because of operational problems, or because of long-term effects such as clogging of boreholes. The results of the study indicate the following: (1) The water curtain system is a very powerful means for preventing gas leakage from underground storage facilities. By operating with appropriate pressure and layout, gas containment can be ensured. (2
Direct Numerical Simulation of Autoiginition of a Hydrogen Jet in a
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Preheated Cross Flow | Argonne Leadership Computing Facility Numerical Simulation of Autoiginition of a Hydrogen Jet in a Preheated Cross Flow Authors: Abdilghanie, A., Frouzakis, C.E., Fischer, P Autoignition of a nitrogen-diluted hydrogen mixture issuing from a round nozzle into a cross-flowing turbulent stream of preheated air flowing in a channel at a friction Reynolds number Re = 180 is inves- tigated via 3-D direct numerical simulations (DNS) at two crossflow stream temperatures (930
Numerical simulation of plasma heating of a composite powder particle
Demetriou, M.D.; Lavine, A.S.; Ghoniem, N.M.
1999-07-01
The use of fine composite powder particles (composed of a ceramic core and a metallic coating) in plasma spraying processes is desirable in developing thin film coatings that possess high abrasion as well as high fracture resistance. Quantitative knowledge of the thermal behavior of a composite particle in a plasma beam is essential in optimizing the process variables to achieve uniform melting of the coating material. In this work, a numerical model is developed to analyze the in-flight thermal behavior of a spherically symmetric WC-Co composite particle travelling in an argon arc-jet DC plasma under strongly unsteady plasma conditions. The model gives quantitative as well as qualitative information about the thermal response of the heated particle. The important features that are addressed are the temperature response of the particle; the history of the location of the melting and vaporization fronts; and the physical state of the particle at the end of its flight. For the conditions investigated, it was determined that the internal conduction resistance is negligible as compared to the net external resistance. However, the presence of the ceramic base was found to affect the transient heating process since its content in the particle composition determines the time constant of the process. Another interesting observation is that proper selection of the particle injection speed and injection location can be effective means for optimizing the heating process and achieving uniform melting of the coating material.
On Improving Analytical Models of Cosmic Reionization for Matching Numerical Simulations
Kaurov, Alexander A.
2016-01-01
The methods for studying the epoch of cosmic reionization vary from full radiative transfer simulations to purely analytical models. While numerical approaches are computationally expensive and are not suitable for generating many mock catalogs, analytical methods are based on assumptions and approximations. We explore the interconnection between both methods. First, we ask how the analytical framework of excursion set formalism can be used for statistical analysis of numerical simulations and visual representation of the morphology of ionization fronts. Second, we explore the methods of training the analytical model on a given numerical simulation. We present a new code which emerged from this study. Its main application is to match the analytical model with a numerical simulation. Then, it allows one to generate mock reionization catalogs with volumes exceeding the original simulation quickly and computationally inexpensively, meanwhile reproducing large scale statistical properties. These mock catalogs are particularly useful for CMB polarization and 21cm experiments, where large volumes are required to simulate the observed signal.
Numerical simulation of tectonic plates motion and seismic process in Central Asia
Peryshkin, A. Yu.; Makarov, P. V. Eremin, M. O.
2014-11-14
An evolutionary approach proposed in [1, 2] combining the achievements of traditional macroscopic theory of solid mechanics and basic ideas of nonlinear dynamics is applied in a numerical simulation of present-day tectonic plates motion and seismic process in Central Asia. Relative values of strength parameters of rigid blocks with respect to the soft zones were characterized by the ? parameter that was varied in the numerical experiments within ? = 1.11.8 for different groups of the zonal-block divisibility. In general, the numerical simulations of tectonic block motion and accompanying seismic process in the model geomedium indicate that the numerical solutions of the solid mechanics equations characterize its deformation as a typical behavior of a nonlinear dynamic system under conditions of self-organized criticality.
Numerical simulations of a vertical tail of a commercial aircraft with
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
active flow control | Argonne Leadership Computing Facility simulations of a vertical tail of a commercial aircraft with active flow control Authors: Rasquin, M., Martin, J., Jansen, K. A series of numerical simulations of a realistic vertical tail of a commercial aircraft, with a tapered swept stabilizer and a rudder, is considered in this work with application of flow control. Flow control is known to have the capacity to augment the streamwise momentum near the rudder suction peak where
Numerical Simulation of Horizontal Continuous Casting Process of C194 Copper Alloy
Huang Guojie; Xie Shuisheng; Cheng Lei; Cheng Zhenkang [State Key Laboratory for Fabrication and Processing of Nonferrous Metals, Beijing General Research Institute for Non-ferrous Metals, China, 100088 (China)
2007-05-17
Horizontal Continuous Casting (H.C.C) is an important method to cast C194 copper ingot. In this paper, numerical simulation is adopted to investigate the casting process in order to optimize the H.C.C technical parameters, such as the casting temperature, casting speed and cooling intensity. According to the numerical results, the reasonable parameters are that the casting temperature is between 1383K{approx}1463K, the casting speed is between 7.2m/h{approx}10.8m/h and the speed of cooling water is between 3.6m/s{approx}4.6m/s. The results of numerical simulation provide the significant reference to the subsequent experiments.
Dynamical properties of fractal networks: Scaling, numerical simulations, and physical realizations
Nakayama, T.; Yakubo, K. ); Orbach, R.L. )
1994-04-01
This article describes the advances that have been made over the past ten years on the problem of fracton excitations in fractal structures. The relevant systems to this subject are so numerous that focus is limited to a specific structure, the percolating network. Recent progress has followed three directions: scaling, numerical simulations, and experiment. In a happy coincidence, large-scale computations, especially those involving array processors, have become possible in recent years. Experimental techniques such as light- and neutron-scattering experiments have also been developed. Together, they form the basis for a review article useful as a guide to understanding these developments and for charting future research directions. In addition, new numerical simulation results for the dynamical properties of diluted antiferromagnets are presented and interpreted in terms of scaling arguments. The authors hope this article will bring the major advances and future issues facing this field into clearer focus, and will stimulate further research on the dynamical properties of random systems.
2D numerical simulation of the MEP energy-transport model with a finite difference scheme
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.
Rossi, Tuomas P. Sakko, Arto; Puska, Martti J.; Lehtola, Susi; Nieminen, Risto M.
2015-03-07
We present an approach for generating local numerical basis sets of improving accuracy for first-principles nanoplasmonics simulations within time-dependent density functional theory. The method is demonstrated for copper, silver, and gold nanoparticles that are of experimental interest but computationally demanding due to the semi-core d-electrons that affect their plasmonic response. The basis sets are constructed by augmenting numerical atomic orbital basis sets by truncated Gaussian-type orbitals generated by the completeness-optimization scheme, which is applied to the photoabsorption spectra of homoatomic metal atom dimers. We obtain basis sets of improving accuracy up to the complete basis set limit and demonstrate that the performance of the basis sets transfers to simulations of larger nanoparticles and nanoalloys as well as to calculations with various exchange-correlation functionals. This work promotes the use of the local basis set approach of controllable accuracy in first-principles nanoplasmonics simulations and beyond.
Analysis of the flamelet concept in the numerical simulation of laminar partially premixed flames
Consul, R.; Oliva, A.; Perez-Segarra, C.D.; Carbonell, D.; de Goey, L.P.H.
2008-04-15
The aim of this work is to analyze the application of flamelet models based on the mixture fraction variable and its dissipation rate to the numerical simulation of partially premixed flames. Although the main application of these models is the computation of turbulent flames, this work focuses on the performance of flamelet concept in laminar flame simulations removing, in this way, turbulence closure interactions. A well-known coflow methane/air laminar flame is selected. Five levels of premixing are taken into account from an equivalence ratio {phi}={infinity} (nonpremixed) to {phi}=2.464. Results obtained using the flamelet approaches are compared to data obtained from the detailed solution of the complete transport equations using primitive variables. Numerical simulations of a counterflow flame are also presented to support the discussion of the results. Special emphasis is given to the analysis of the scalar dissipation rate modeling. (author)
Wimmer, Thomas Srimathveeravalli, Govindarajan; Gutta, Narendra; Ezell, Paula C.; Monette, Sebastien; Maybody, Majid; Erinjery, Joseph P.; Durack, Jeremy C.; Coleman, Jonathan A.; Solomon, Stephen B.
2015-02-15
PurposeNumerical simulations are used for treatment planning in clinical applications of irreversible electroporation (IRE) to determine ablation size and shape. To assess the reliability of simulations for treatment planning, we compared simulation results with empiric outcomes of renal IRE using computed tomography (CT) and histology in an animal model.MethodsThe ablation size and shape for six different IRE parameter sets (70–90 pulses, 2,000–2,700 V, 70–100 µs) for monopolar and bipolar electrodes was simulated using a numerical model. Employing these treatment parameters, 35 CT-guided IRE ablations were created in both kidneys of six pigs and followed up with CT immediately and after 24 h. Histopathology was analyzed from postablation day 1.ResultsAblation zones on CT measured 81 ± 18 % (day 0, p ≤ 0.05) and 115 ± 18 % (day 1, p ≤ 0.09) of the simulated size for monopolar electrodes, and 190 ± 33 % (day 0, p ≤ 0.001) and 234 ± 12 % (day 1, p ≤ 0.0001) for bipolar electrodes. Histopathology indicated smaller ablation zones than simulated (71 ± 41 %, p ≤ 0.047) and measured on CT (47 ± 16 %, p ≤ 0.005) with complete ablation of kidney parenchyma within the central zone and incomplete ablation in the periphery.ConclusionBoth numerical simulations for planning renal IRE and CT measurements may overestimate the size of ablation compared to histology, and ablation effects may be incomplete in the periphery.
A new dipolar potential for numerical simulations of polar fluids on the 4D hypersphere
Caillol, Jean-Michel; Trulsson, Martin
2014-09-28
We present a new method for Monte Carlo or Molecular Dynamics numerical simulations of three-dimensional polar fluids. The simulation cell is defined to be the surface of the northern hemisphere of a four-dimensional (hyper)sphere. The point dipoles are constrained to remain tangent to the sphere and their interactions are derived from the basic laws of electrostatics in this geometry. The dipole-dipole potential has two singularities which correspond to the following boundary conditions: when a dipole leaves the northern hemisphere at some point of the equator, it reappears at the antipodal point bearing the same dipole moment. We derive all the formal expressions needed to obtain the thermodynamic and structural properties of a polar liquid at thermal equilibrium in actual numerical simulation. We notably establish the expression of the static dielectric constant of the fluid as well as the behavior of the pair correlation at large distances. We report and discuss the results of extensive numerical Monte Carlo simulations for two reference states of a fluid of dipolar hard spheres and compare these results with previous methods with a special emphasis on finite size effects.
Rutqvist, Jonny; Rutqvist, J.; Moridis, G.J.
2008-06-01
In this paper, we describe the development and application of a numerical simulator that analyzes the geomechanical performance of hydrate-bearing sediments, which may become an important future energy supply. The simulator is developed by coupling a robust numerical simulator of coupled fluid flow, hydrate thermodynamics, and phase behavior in geologic media (TOUGH+HYDRATE) with an established geomechanical code (FLAC3D). We demonstrate the current simulator capabilities and applicability for two examples of geomechanical responses of hydrate bearing sediments during production-induced hydrate dissociation. In these applications, the coupled geomechanical behavior within hydrate-bearing seducements are considered through a Mohr-Coulomb constitutive model, corrected for changes in pore-filling hydrate and ice content, based on laboratory data. The results demonstrate how depressurization-based gas production from oceanic hydrate deposits may lead to severe geomechanical problems unless care is taken in designing the production scheme. We conclude that the coupled simulator can be used to design production strategies for optimizing production, while avoiding damaging geomechanical problems.
Direct numerical simulations of fluid flow, heat transfer and phase changes
Juric, D.; Tryggvason, G.; Han, J.
1997-04-01
Direct numerical simulations of fluid flow, heat transfer, and phase changes are presented. The simulations are made possible by a recently developed finite difference/front tracking method based on the one-field formulation of the governing equations where a single set of conservation equations is written for all the phases involved. The conservation equations are solved on a fixed rectangular grid, but the phase boundaries are kept sharp by tracking them explicitly by a moving grid of lower dimension. The method is discussed and applications to boiling heat transfer and the solidification of drops colliding with a wall are shown.
Numerical simulation of alumina spraying in argon-helium plasma jet
Chang, C.H.
1992-01-01
A new numerical model is described for simulating thermal plasmas containing entrained particles, with emphasis on plasma spraying applications. The plasma is represented as a continuum multicomponent chemically reacting ideal gas, while the particles are tracked as discrete Lagrangian entities coupled to the plasma. Computational results are presented from a transient simulation of alumina spraying in a turbulent argon-helium plasma jet in air environment, including torch geometry, substrate, and multiple species with chemical reactions. Particle-plasma interactions including turbulent dispersion have been modeled in a fully self-consistent manner. Interactions between the plasma and the torch and substrate walls are modeled using wall functions. (15 refs.)
Numerical simulation of alumina spraying in argon-helium plasma jet
Chang, C.H.
1992-08-01
A new numerical model is described for simulating thermal plasmas containing entrained particles, with emphasis on plasma spraying applications. The plasma is represented as a continuum multicomponent chemically reacting ideal gas, while the particles are tracked as discrete Lagrangian entities coupled to the plasma. Computational results are presented from a transient simulation of alumina spraying in a turbulent argon-helium plasma jet in air environment, including torch geometry, substrate, and multiple species with chemical reactions. Particle-plasma interactions including turbulent dispersion have been modeled in a fully self-consistent manner. Interactions between the plasma and the torch and substrate walls are modeled using wall functions. (15 refs.)
Holladay, Jamelyn D.; Wang, Yong
2015-05-01
Microscale (<5W) reformers for hydrogen production have been investigated for over a decade. These devices are intended to provide hydrogen for small fuel cells. Due to the reformer’s small size, numerical simulations are critical to understand heat and mass transfer phenomena occurring in the systems. This paper reviews the development of the numerical codes and details the reaction equations used. The majority of the devices utilized methanol as the fuel due to methanol’s low reforming temperature and high conversion, although, there are several methane fueled systems. As computational power has decreased in cost and increased in availability, the codes increased in complexity and accuracy. Initial models focused on the reformer, while more recently, the simulations began including other unit operations such as vaporizers, inlet manifolds, and combustors. These codes are critical for developing the next generation systems. The systems reviewed included, plate reactors, microchannel reactors, annulus reactors, wash-coated, packed bed systems.
Rider, William; Kamm, J. R.; Tomkins, C. D.; Zoldi, C. A.; Prestridge, K. P.; Marr-Lyon, M.; Rightley, P. M.; Benjamin, R. F.
2002-01-01
We consider the detailed structures of mixing flows for Richtmyer-Meshkov experiments of Prestridge et al. [PRE 00] and Tomkins et al. [TOM 01] and examine the most recent measurements from the experimental apparatus. Numerical simulations of these experiments are performed with three different versions of high resolution finite volume Godunov methods. We compare experimental data with simulations for configurations of one and two diffuse cylinders of SF{sub 6} in air using integral measures as well as fractal analysis and continuous wavelet transforms. The details of the initial conditions have a significant effect on the computed results, especially in the case of the double cylinder. Additionally, these comparisons reveal sensitive dependence of the computed solution on the numerical method.
Mikellides, Ioannis G.; Goebel, Dan M.; Snyder, John Steven; Katz, Ira; Herman, Daniel A.
2010-12-01
Numerical simulations of neutralizer hollow cathodes at various operating conditions and orifice sizes are presented. The simulations were performed using a two-dimensional axisymmetric model that solves numerically an extensive system of conservation laws for the partially ionized gas in these devices. The results for the plasma are compared directly with Langmuir probe measurements. The computed keeper voltages are also compared with the observed values. Whenever model inputs and/or specific physics of the cathode discharge were uncertain or unknown additional sensitivity calculations have been performed to quantify the uncertainties. The model has also been employed to provide insight into recent ground test observations of the neutralizer cathode in NASA's evolutionary xenon thruster. It is found that a likely cause of the observed keeper voltage drop in a long duration test of the engine is cathode orifice erosion.
Numerical simulations of stripping effects in high-intensity hydrogen ion linacs
Carneiro, J.-P.; Mustapha, B.; Ostroumov, P.N.; /Argonne
2008-12-01
Numerical simulations of H{sup -} stripping losses from blackbody radiation, electromagnetic fields, and residual gas have been implemented into the beam dynamics code TRACK. Estimates of the stripping losses along two high-intensity H{sup -} linacs are presented: the Spallation Neutron Source linac currently being operated at Oak Ridge National Laboratory and an 8 GeV superconducting linac currently being designed at Fermi National Accelerator Laboratory.
Bae, Y. Y.; Hong, S. D.; Kim, Y. W.
2012-07-01
A number of computational works have been performed so far for the simulation of heat transfer in a supercritical fluid. The simulations, however, faced a lot of difficulties when heat transfer deteriorates due either to buoyancy or by acceleration. When the bulk temperature approaches the pseudo-critical temperature the fluid experiences a severe axial density gradient on top of a severe radial one. Earlier numerical calculations showed, without exception, unrealistic over-predictions, as soon as the bulk temperature exceeded the pseudo-critical temperature. The over-predictions might have been resulted from an inapplicability of widely-used turbulence models. One of the major causes for the difficulties may probably be an assumption of a constant turbulent Prandtl number. Recent research, both numerical and experimental, indicates that the turbulent Prandtl number is never a constant when the gradient of physical properties is significant. This paper describes the applicability of a variable turbulent Prandtl number to the numerical simulation of heat transfer in supercritical fluids flowing in narrow vertical tubes. (authors)
Numerical Simulations of Subscale Wind Turbine Rotor Inboard Airfoils at Low Reynolds Number
Blaylock, Myra L.; Maniaci, David Charles; Resor, Brian R.
2015-04-01
New blade designs are planned to support future research campaigns at the SWiFT facility in Lubbock, Texas. The sub-scale blades will reproduce specific aerodynamic characteristics of utility-scale rotors. Reynolds numbers for megawatt-, utility-scale rotors are generally above 2-8 million. The thickness of inboard airfoils for these large rotors are typically as high as 35-40%. The thickness and the proximity to three-dimensional flow of these airfoils present design and analysis challenges, even at the full scale. However, more than a decade of experience with the airfoils in numerical simulation, in the wind tunnel, and in the field has generated confidence in their performance. Reynolds number regimes for the sub-scale rotor are significantly lower for the inboard blade, ranging from 0.7 to 1 million. Performance of the thick airfoils in this regime is uncertain because of the lack of wind tunnel data and the inherent challenge associated with numerical simulations. This report documents efforts to determine the most capable analysis tools to support these simulations in an effort to improve understanding of the aerodynamic properties of thick airfoils in this Reynolds number regime. Numerical results from various codes of four airfoils are verified against previously published wind tunnel results where data at those Reynolds numbers are available. Results are then computed for other Reynolds numbers of interest.
An Implicit Algorithm for the Numerical Simulation of Shape-Memory Alloys
Becker, R; Stolken, J; Jannetti, C; Bassani, J
2003-10-16
Shape-memory alloys (SMA) have the potential to be used in a variety of interesting applications due to their unique properties of pseudoelasticity and the shape-memory effect. However, in order to design SMA devices efficiently, a physics-based constitutive model is required to accurately simulate the behavior of shape-memory alloys. The scope of this work is to extend the numerical capabilities of the SMA constitutive model developed by Jannetti et. al. (2003), to handle large-scale polycrystalline simulations. The constitutive model is implemented within the finite-element software ABAQUS/Standard using a user defined material subroutine, or UMAT. To improve the efficiency of the numerical simulations, so that polycrystalline specimens of shape-memory alloys can be modeled, a fully implicit algorithm has been implemented to integrate the constitutive equations. Using an implicit integration scheme increases the efficiency of the UMAT over the previously implemented explicit integration method by a factor of more than 100 for single crystal simulations.
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.
Some effects of data base variations on numerical simulations of uranium migration
Carnahan, C.L.
1987-12-01
Numerical simulations of migration of chemicals in the geosphere depend on knowledge of identities of chemical species and on values of chemical equilibrium constants supplied to the simulators. In this work, some effects of variability in assumed speciation and in equilibrium constants were examined, using migration of uranium as an example. Various simulations were done of uranium migration in systems with varying oxidation potential, pH, and mator component content. A simulation including formation of aqueous species UO/sub 2//sup 2 +/, UO/sub 2/CO/sub 3//sup 0/, UO/sub 2/(CO/sub 3/)/sub 2//sup 2 -/, UO/sub 2/(CO/sub 3/)/sub 3//sup 4 -/, (UO/sub 2/)/sub 2/CO/sub 3/(OH)/sub 3//sup -/, UO/sub 2//sup +/, U(OH)/sub 4//sup 0/, and U(OH)/sub 5//sup -/ is compared to simulation excluding formation of UO/sub 2//sup +/ and U(OH)/sub 5//sup -/. These simulations relied on older data bases, and they are compared to a further simulation using recently published data on formation of U(OH)/sub 4//sup 0/, (UO/sub 2/)/sub 2/CO/sub 3/(OH)/sub 3//sup -/, UO/sub 2/(CO/sub 3/)/sub 5//sup 5 -/, and U(CO/sub 3/)/sub 5//sup 6 -/. Significant differences in dissolved uranium concentrations are noted among the simulations. Differences are noted also in precipitation of two solids, USiO/sub 4/(c) (coffinite) and CaUO/sub 4/(c) (calcium uranate), although the solubility products of the solids were not varied in the simulations. 18 refs., 9 figs., 2 tabs.
Two-dimensional numerical simulation of boron diffusion for pyramidally textured silicon
Ma, Fa-Jun Duttagupta, Shubham; Shetty, Kishan Devappa; Meng, Lei; Hoex, Bram; Peters, Ian Marius; Samudra, Ganesh S.
2014-11-14
Multidimensional numerical simulation of boron diffusion is of great relevance for the improvement of industrial n-type crystalline silicon wafer solar cells. However, surface passivation of boron diffused area is typically studied in one dimension on planar lifetime samples. This approach neglects the effects of the solar cell pyramidal texture on the boron doping process and resulting doping profile. In this work, we present a theoretical study using a two-dimensional surface morphology for pyramidally textured samples. The boron diffusivity and segregation coefficient between oxide and silicon in simulation are determined by reproducing measured one-dimensional boron depth profiles prepared using different boron diffusion recipes on planar samples. The established parameters are subsequently used to simulate the boron diffusion process on textured samples. The simulated junction depth is found to agree quantitatively well with electron beam induced current measurements. Finally, chemical passivation on planar and textured samples is compared in device simulation. Particularly, a two-dimensional approach is adopted for textured samples to evaluate chemical passivation. The intrinsic emitter saturation current density, which is only related to Auger and radiative recombination, is also simulated for both planar and textured samples. The differences between planar and textured samples are discussed.
Godfrey, Brendan B.; Vay, Jean-Luc
2013-09-01
Rapidly growing numerical instabilities routinely occur in multidimensional particle-in-cell computer simulations of plasma-based particle accelerators, astrophysical phenomena, and relativistic charged particle beams. Reducing instability growth to acceptable levels has necessitated higher resolution grids, high-order field solvers, current filtering, etc. except for certain ratios of the time step to the axial cell size, for which numerical growth rates and saturation levels are reduced substantially. This paper derives and solves the cold beam dispersion relation for numerical instabilities in multidimensional, relativistic, electromagnetic particle-in-cell programs employing either the standard or the ColeKarkkainnen finite difference field solver on a staggered mesh and the common Esirkepov current-gathering algorithm. Good overall agreement is achieved with previously reported results of the WARP code. In particular, the existence of select time steps for which instabilities are minimized is explained. Additionally, an alternative field interpolation algorithm is proposed for which instabilities are almost completely eliminated for a particular time step in ultra-relativistic simulations.
Figueroa, Aldo [Facultad de Ciencias, Universidad Autnoma del Estado de Morelos, Cuernavaca, Morelos 62209 (Mexico)] [Facultad de Ciencias, Universidad Autnoma del Estado de Morelos, Cuernavaca, Morelos 62209 (Mexico); Meunier, Patrice; Villermaux, Emmanuel [Aix-Marseille Univ., CNRS, Centrale Marseille, IRPHE, Marseille F-13384 (France)] [Aix-Marseille Univ., CNRS, Centrale Marseille, IRPHE, Marseille F-13384 (France); Cuevas, Sergio; Ramos, Eduardo [Instituto de Energas Renovables, Universidad Nacional Autnoma de Mxico, A.P. 34, Temixco, Morelos 62580 (Mexico)] [Instituto de Energas Renovables, Universidad Nacional Autnoma de Mxico, A.P. 34, Temixco, Morelos 62580 (Mexico)
2014-01-15
We present a combination of experiment, theory, and modelling on laminar mixing at large Pclet number. The flow is produced by oscillating electromagnetic forces in a thin electrolytic fluid layer, leading to oscillating dipoles, quadrupoles, octopoles, and disordered flows. The numerical simulations are based on the Diffusive Strip Method (DSM) which was recently introduced (P. Meunier and E. Villermaux, The diffusive strip method for scalar mixing in two-dimensions, J. Fluid Mech. 662, 134172 (2010)) to solve the advection-diffusion problem by combining Lagrangian techniques and theoretical modelling of the diffusion. Numerical simulations obtained with the DSM are in reasonable agreement with quantitative dye visualization experiments of the scalar fields. A theoretical model based on log-normal Probability Density Functions (PDFs) of stretching factors, characteristic of homogeneous turbulence in the Batchelor regime, allows to predict the PDFs of scalar in agreement with numerical and experimental results. This model also indicates that the PDFs of scalar are asymptotically close to log-normal at late stages, except for the large concentration levels which correspond to low stretching factors.
On a framework for generating PoD curves assisted by numerical simulations
Subair, S. Mohamed Agrawal, Shweta Balasubramaniam, Krishnan Rajagopal, Prabhu; Kumar, Anish; Rao, Purnachandra B.; Tamanna, Jayakumar
2015-03-31
The Probability of Detection (PoD) curve method has emerged as an important tool for the assessment of the performance of NDE techniques, a topic of particular interest to the nuclear industry where inspection qualification is very important. The conventional experimental means of generating PoD curves though, can be expensive, requiring large data sets (covering defects and test conditions), and equipment and operator time. Several methods of achieving faster estimates for PoD curves using physics-based modelling have been developed to address this problem. Numerical modelling techniques are also attractive, especially given the ever-increasing computational power available to scientists today. Here we develop procedures for obtaining PoD curves, assisted by numerical simulation and based on Bayesian statistics. Numerical simulations are performed using Finite Element analysis for factors that are assumed to be independent, random and normally distributed. PoD curves so generated are compared with experiments on austenitic stainless steel (SS) plates with artificially created notches. We examine issues affecting the PoD curve generation process including codes, standards, distribution of defect parameters and the choice of the noise threshold. We also study the assumption of normal distribution for signal response parameters and consider strategies for dealing with data that may be more complex or sparse to justify this. These topics are addressed and illustrated through the example case of generation of PoD curves for pulse-echo ultrasonic inspection of vertical surface-breaking cracks in SS plates.
Van Eerten, Hendrik J.; MacFadyen, Andrew I.
2012-06-01
We discuss jet dynamics for narrow and wide gamma-ray burst (GRB) afterglow jets and the observational implications of numerical simulations of relativistic jets in two dimensions. We confirm earlier numerical results that sideways expansion of relativistic jets during the bulk of the afterglow emission phase is logarithmic in time and find that this also applies to narrow jets with half opening angle of 0.05 rad. As a result, afterglow jets remain highly nonspherical until after they have become nonrelativistic. Although sideways expansion steepens the afterglow light curve after the jet break, the jet edges becoming visible dominates the jet break, which means that the jet break is sensitive to the observer angle even for narrow jets. Failure to take the observer angle into account can lead to an overestimation of the jet energy by up to a factor of four. This weakens the challenge posed to the magneter energy limit by extreme events such as GRB090926A. Late-time radio calorimetry based on a spherical nonrelativistic outflow model remains relevant when the observer is approximately on-axis and where differences of a few in flux level between the model and the simulation are acceptable. However, this does not imply sphericity of the outflow and therefore does not translate to high observer angles relevant to orphan afterglows. For more accurate calorimetry and in order to model significant late-time features such as the rise of the counterjet, detailed jet simulations remain indispensable.
Terascale direct numerical simulations of turbulent combustion using S3D
Chen, Jackie; Klasky, Scott A; Hawkes, Evatt R; Sankaran, Ramanan; Choudhary, Alok; Yoo, Chun S; Liao, Wei-keng; Podhorszki, Norbert
2009-01-01
Computational science is paramount to the understanding of underlying processes in internal combustion engines of the future that will utilize non-petroleum-based alternative fuels, including carbon-neutral biofuels, and burn in new combustion regimes that will attain high efficiency while minimizing emissions of particulates and nitrogen oxides. Next-generation engines will likely operate at higher pressures, with greater amounts of dilution and utilize alternative fuels that exhibit a wide range of chemical and physical properties. Therefore, there is a significant role for high-fidelity simulations, direct numerical simulations (DNS), specifically designed to capture key turbulence-chemistry interactions in these relatively uncharted combustion regimes, and in particular, that can discriminate the effects of differences in fuel properties. In DNS, all of the relevant turbulence and flame scales are resolved numerically using high-order accurate numerical algorithms. As a consequence terascale DNS are computationally intensive, require massive amounts of computing power and generate tens of terabytes of data. Recent results from terascale DNS of turbulent flames are presented here, illustrating its role in elucidating flame stabilization mechanisms in a lifted turbulent hydrogen/air jet flame in a hot air coflow, and the flame structure of a fuel-lean turbulent premixed jet flame. Computing at this scale requires close collaborations between computer and combustion scientists to provide optimized scaleable algorithms and software for terascale simulations, efficient collective parallel I/O, tools for volume visualization of multiscale, multivariate data and automating the combustion workflow. The enabling computer science, applied to combustion science, is also required in many other terascale physics and engineering simulations. In particular, performance monitoring is used to identify the performance of key kernels in the DNS code, S3D and especially memory
Numerical simulations of impulsively generated Alfvn waves in solar magnetic arcades
Chmielewski, P.; Murawski, K.; Musielak, Z. E.; Srivastava, A. K.
2014-09-20
We perform numerical simulations of impulsively generated Alfvn waves in an isolated solar arcade, which is gravitationally stratified and magnetically confined. We study numerically the propagation of Alfvn waves along the magnetic structure that extends from the lower chromosphere, where the waves are generated, to the solar corona, and analyze the influence of the arcade size and the width of the initial pulses on the wave propagation and reflection. Our model of the solar atmosphere is constructed by adopting the temperature distribution based on the semi-empirical VAL-C model and specifying the curved magnetic field lines that constitute the asymmetric magnetic arcade. The propagation and reflection of Alfvn waves in this arcade is described by 2.5-dimensional magnetohydrodynamic equations that are numerically solved by the FLASH code. Our numerical simulations reveal that the Alfvn wave amplitude decreases as a result of a partial reflection of Alfvn waves in the solar transition region, and that the waves that are not reflected leak through the transition region and reach the solar corona. We also find the decrement of the attenuation time of Alfvn waves for wider initial pulses. Moreover, our results show that the propagation of Alfvn waves in the arcade is affected by the spatial dependence of the Alfvn speed, which leads to phase mixing that is stronger for more curved and larger magnetic arcades. We discuss the processes that affect the Alfvn wave propagation in an asymmetric solar arcade and conclude that besides phase mixing in the magnetic field configuration, the plasma properties of the arcade, the size of the initial pulse, and the structure of the solar transition region all play a vital role in the Alfvn wave propagation.
THE REBOUND CONDITION OF DUST AGGREGATES REVEALED BY NUMERICAL SIMULATION OF THEIR COLLISIONS
Wada, Koji; Tanaka, Hidekazu; Yamamoto, Tetsuo; Suyama, Toru; Kimura, Hiroshi
2011-08-10
Collisional growth of dust aggregates is a plausible root of planetesimals forming in protoplanetary disks. However, a rebound of colliding dust aggregates prevents dust from growing into planetesimals. In fact, rebounding aggregates are observed in laboratory experiments but not in previous numerical simulations. Therefore, the condition of rebound between dust aggregates should be clarified to better understand the processes of dust growth and planetesimal formation. We have carried out numerical simulations of aggregate collisions for various types of aggregates and succeeded in reproducing a rebound of colliding aggregates under specific conditions. Our finding is that in the rebound process, the key factor of the aggregate structure is the coordination number, namely, the number of particles in contact with a particle. A rebound is governed by the energy dissipation along with restructuring of the aggregates and a large coordination number inhibits the restructuring at collisions. Results of our numerical simulation for various aggregates indicate that they stick to each other when the mean coordination number is less than 6, regardless of their materials and structures, as long as their collision velocity is less than the critical velocity for fragmentation. This criterion of the coordination number would correspond to a filling factor of {approx}0.3, which is somewhat larger than that reported in laboratory experiments. In protoplanetary disks, dust aggregates are expected to have low bulk densities (<0.1 g cm{sup -3}) during their growth, which would prevent dust aggregates from rebounding. This result supports the formation of planetesimals with direct dust growth in protoplanetary disks.
Naziar, J.; Couch, R.; Davis, M.
1996-01-01
Traditionally, aeropropulsion structural performance and aerodynamic performance have been designed separately and later mated together via flight testing. In today`s atmosphere of declining resources, it is imperative that more productive ways of designing and verifying aeropropulsion performance and structural interaction be made available to the aerospace industry. One method of obtaining a more productive design and evaluation capability is through the use of numerical simulations. Currently, Lawrence Livermore National Laboratory has developed a generalized fluid/structural interaction code known as ALE3D. This code is capable of characterizing fluid and structural interaction for components such as the combustor, fan/stators, inlet and/or nozzles. This code solves the 3D Euler equations and has been applied to several aeropropulsion applications such as a supersonic inlet and a combustor rupture simulation. To characterize aerodynamic-structural interaction for rotating components such as the compressor, appropriate turbomachinery simulations would need to be implemented within the ALE3D structure. The Arnold Engineering Development Center is currently developing a three-dimensional compression system code known as TEACC (Turbine Engine Analysis Compressor Code). TEACC also solves the 3D Euler equations and is intended to simulate dynamic behavior such as inlet distortion, surge or rotating stall. The technology being developed within the TEACC effort provides the necessary turbomachinery simulation for implementation into ALE3D. This paper describes a methodology to combine three-dimensional aerodynamic turbomachinery technology into the existing aerodynamic-structural interaction simulation, ALE3D to obtain the desired aerodynamic and structural integrated simulation for an aeropropulsion system.
Numerical simulation of fracture rocks and wave propagation by means of fractal theory
Valle G., R. del
1994-12-31
A numerical approach was developed for the dynamic simulation of fracture rocks and wave propagation. Based on some ideas of percolation theory and fractal growth, a network of particles and strings represent the rock model. To simulate an inhomogeneous medium, the particles and springs have random distributed elastic parameters and are implemented in the dynamic Navier equation. Some of the springs snap with criteria based on the confined stress applied, therefore creating a fractured rock consistent with the physical environment. The basic purpose of this research was to provide a method to construct a fractured rock with confined stress conditions as well as the wave propagation imposed in the model. Such models provide a better understanding of the behavior of wave propagation in fractured media. The synthetic seismic data obtained henceforth, can be used as a tool to develop methods for characterizing fractured rocks by means of geophysical inference.
Direct numerical simulation of turbulent flow in a rotating square duct
Dai, Yi-Jun; Huang, Wei-Xi Xu, Chun-Xiao; Cui, Gui-Xiang
2015-06-15
A fully developed turbulent flow in a rotating straight square duct is simulated by direct numerical simulations at Re{sub ?} = 300 and 0 ? Ro{sub ?} ? 40. The rotating axis is parallel to two opposite walls of the duct and normal to the main flow. Variations of the turbulence statistics with the rotation rate are presented, and a comparison with the rotating turbulent channel flow is discussed. Rich secondary flow patterns in the cross section are observed by varying the rotation rate. The appearance of a pair of additional vortices above the pressure wall is carefully examined, and the underlying mechanism is explained according to the budget analysis of the mean momentum equations.
Numerical simulation to study the transient self focusing of laser beam in plasma
Sharma, R. P.; Hussain, Saba Gaur, Nidhi
2015-02-15
In this paper, we present the numerical simulation for the coupled system of equations governing the dynamics of laser and Ion Acoustic Wave (IAW) in a collisionless plasma, when the coupling between the waves is through ponderomotive non-linearity. The nonlinear evolution of the laser beam is studied when the pump laser is perturbed by a periodic perturbation. By changing the perturbation wave number, we have studied its effect on the nonlinear evolution pattern of laser beam. In order to have a physical insight into the nonlinear dynamics of laser beam evolution in time and space, we have studied the laser and IAW spectra containing spatial harmonics. The magnitude of these harmonics changes with time and leads to time dependent localization of laser beam in spatial domain. The nonlinear dynamics of this localization is investigated in detail by using simulation and a semi-analytical model.
Romero-Talam?s, C A; Hooper, E B; Hill, D N; Cohen, B I; McLean, H S; Wood, R D; Moller, J M
2006-03-15
Data from a recently installed insertable magnetic probe array in the Sustained Spheromak Physics Experiment (SSPX) [E. B. Hooper et al., Nucl. Fusion 39, 863 (1999)] is compared against NIMROD [C. R. Sovinec et al., J. Comp. Phys. 195, 355 (2004)], a full 3D resistive magnetohydrodynamic code that is used to simulate SSPX plasmas. The experiment probe consists of a linear array of chip inductors arranged in clusters that are spaced every 2 cm, and spans the entire machine radius at the flux conserver midplane. Both the experiment and the numerical simulations show the appearance, shortly after breakdown, of a column with a hollow current profile that precedes magnetic reconnection, a process essential to the formation of closed magnetic flux surfaces. However, there are differences between the experiment and the simulation in how the column evolves after it is formed. These differences are studied to help identify the mechanisms that eventually lead to closed-flux surfaces (azimuthally averaged) and flux amplification, which occur in both the experiment and the simulation.
Oostrom, Martinus; Wietsma, Thomas W.; Strickland, Christopher E.; Freedman, Vicky L.; Truex, Michael J.
2012-02-01
Soil desiccation, in conjunction with surface infiltration control, is considered at the Hanford Site as a potential technology to limit the flux of technetium and other contaminants in the vadose zone to the groundwater. An intermediate-scale experiment was conducted to test the response of a series of instruments to desiccation and subsequent rewetting of porous media. The instruments include thermistors, thermocouple psychrometers, dual-probe heat pulse sensors, heat dissipation units, and humidity probes. The experiment was simulated with the multifluid flow simulator STOMP, using independently obtained hydraulic and thermal porous medium properties. All instrument types used for this experiment were able to indicate when the desiccation front passed a certain location. In most cases the changes were sharp, indicating rapid changes in moisture content, water potential, or humidity. However, a response to the changing conditions was recorded only when the drying front was very close to a sensor. Of the tested instruments, only the heat dissipation unit and humidity probes were able to detect rewetting. The numerical simulation results reasonably match the experimental data, indicating that the simulator captures the pertinent gas flow and transport processes related to desiccation and rewetting and may be useful in the design and analysis of field tests.
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Austin, Ryan A.; Barton, Nathan R.; Reaugh, John E.; Fried, Laurence E.
2015-05-14
A numerical model is developed to study the shock wave ignition of HMX crystal. The model accounts for the coupling between crystal thermal/mechanical responses and chemical reactions that are driven by the temperature field. This allows for the direct numerical simulation of decomposition reactions in the hot spots formed by shock/impact loading. The model is used to simulate intragranular pore collapse under shock wave loading. In a reference case: (i) shear-enabled micro-jetting is responsible for a modest extent of reaction in the pore collapse region, and (ii) shear banding is found to be an important mode of localization. The shearmore » bands, which are filled with molten HMX, grow out of the pore collapse region and serve as potential ignition sites. The model predictions of shear banding and reactivity are found to be quite sensitive to the respective flow strengths of the solid and liquid phases. In this regard, it is shown that reasonable assumptions of liquid-HMX viscosity can lead to chemical reactions within the shear bands on a nanosecond time scale.« less
Austin, Ryan A.; Barton, Nathan R.; Reaugh, John E.; Fried, Laurence E.
2015-05-14
A numerical model is developed to study the shock wave ignition of HMX crystal. The model accounts for the coupling between crystal thermal/mechanical responses and chemical reactions that are driven by the temperature field. This allows for the direct numerical simulation of decomposition reactions in the hot spots formed by shock/impact loading. The model is used to simulate intragranular pore collapse under shock wave loading. In a reference case: (i) shear-enabled micro-jetting is responsible for a modest extent of reaction in the pore collapse region, and (ii) shear banding is found to be an important mode of localization. The shear bands, which are filled with molten HMX, grow out of the pore collapse region and serve as potential ignition sites. The model predictions of shear banding and reactivity are found to be quite sensitive to the respective flow strengths of the solid and liquid phases. In this regard, it is shown that reasonable assumptions of liquid-HMX viscosity can lead to chemical reactions within the shear bands on a nanosecond time scale.
NUMERICAL SIMULATIONS OF THE EFFECTS OF CHANGING FUEL FOR TURBINES FIRED BY NATURAL GAS AND SYNGAS
Sabau, Adrian S; Wright, Ian G
2007-01-01
Gas turbines in integrated gasification combined cycle (IGCC) power plants burn a fuel gas (syngas) in which the proportions of hydrocarbons, H2, CO, water vapor, and minor impurity levels may vary significantly from those in natural gas, depending on the input feed to the gasifier and the gasification process. A data structure and computational methodology is presented for the numerical simulation of a turbine thermodynamic cycle for various fuel types, air/fuel ratios, and coolant flow rates. The approach used allowed efficient handling of turbine components and different variable constraints due to fuel changes. Examples are presented for a turbine with four stages and cooled blades. The blades were considered to be cooled in an open circuit, with air provided from appropriate compressor stages. Results are presented for the temperatures of the hot gas, alloy surface (coating-superalloy interface), and coolant, as well as for cooling flow rates. Based on the results of the numerical simulations, values were calculated for the fuel flow rates, airflow ratios, and coolant flow rates required to maintain the superalloy in the first stage blade at the desired temperature when the fuel was changed from natural gas (NG) to syngas (SG). One NG case was conducted to assess the effect of coolant pressure matching between the compressor extraction points and corresponding turbine injection points. It was found that pressure matching is a feature that must be considered for high combustion temperatures. The first series of SG simulations was conducted using the same inlet mass flow and pressure ratios as those for the NG case. The results showed that higher coolant flow rates and a larger number of cooled turbine rows were needed for the SG case. Thus, for this first case, the turbine size would be different for SG than for NG. In order to maintain the original turbine configuration (i.e., geometry, diameters, blade heights, angles, and cooling circuit characteristics) for
SCATTERING OF THE f-MODE BY SMALL MAGNETIC FLUX ELEMENTS FROM OBSERVATIONS AND NUMERICAL SIMULATIONS
Felipe, T.; Braun, D.; Crouch, A.; Birch, A.
2012-10-01
The scattering of f-modes by magnetic tubes is analyzed using three-dimensional numerical simulations. An f-mode wave packet is propagated through a solar atmosphere embedded with three different flux tube models that differ in radius and total magnetic flux. A quiet-Sun simulation without a tube present is also performed as a reference. Waves are excited inside the flux tube and propagate along the field lines, and jacket modes are generated in the surroundings of the flux tube, carrying 40% as much energy as the tube modes. The resulting scattered wave is mainly an f-mode composed of a mixture of m = 0 and m = {+-}1 modes. The amplitude of the scattered wave approximately scales with the magnetic flux. A small amount of power is scattered into the p{sub 1}-mode. We have evaluated the absorption and phase shift from a Fourier-Hankel decomposition of the photospheric vertical velocities. They are compared with the results obtained from the ensemble average of 3400 small magnetic elements observed in high-resolution MDI Doppler datacubes. The comparison shows that the observed dependence of the phase shift with wavenumber can be matched reasonably well with the simulated flux tube model. The observed variation of the phase shifts with the azimuthal order m appears to depend on details of the ensemble averaging, including possible motions of the magnetic elements and asymmetrically shaped elements.
Numerical simulation of gas flow through unsaturated fractured rock at Yucca Mountain, Nevada
Cooper, C.A.
1990-01-01
Numerical analysis is used to identify the physical phenomena associated with barometrically driven gas (air and water vapor) flow through unsaturated fractured rock at Yucca Mountain, Nevada. Results from simple finite difference simulations indicate that for a fractured rock scenario, the maximum velocity of air out of an uncased 10 cm borehole is 0.002 m s{sub {minus}1}. An equivalent porous medium (EPM) model was incorporated into a multiphase, multicomponent simulator to test more complex conceptual models. Results indicate that for a typical June day, a diurnal pressure wave propagates about 160 m into the surrounding Tiva Canyon hydrogeologic unit. Dry air that enters the formation evaporates water around the borehole which reduces capillary pressure. Multiphase countercurrent flow develops in the vicinity of the hole; the gas phase flows into the formation while the liquid phase flows toward the borehole. The effect occurs within 0.5 m of the borehole. The amount of water vapor leaving the formation during 1 day is 900 cm{sup 3}. This is less than 0.1% of the total recharge into the formation, suggesting that the barometric effect may be insignificant in drying the unsaturated zone. However, gas phase velocities out of the borehole (3 m s{sup {minus}1}), indicating that observed flow rates from wells along the east flank of Yucca Mountain were able to be simulated with a barometric model.
Dong, S.
2015-02-15
We present a family of physical formulations, and a numerical algorithm, based on a class of general order parameters for simulating the motion of a mixture of N (N⩾2) immiscible incompressible fluids with given densities, dynamic viscosities, and pairwise surface tensions. The N-phase formulations stem from a phase field model we developed in a recent work based on the conservations of mass/momentum, and the second law of thermodynamics. The introduction of general order parameters leads to an extremely strongly-coupled system of (N−1) phase field equations. On the other hand, the general form enables one to compute the N-phase mixing energy density coefficients in an explicit fashion in terms of the pairwise surface tensions. We show that the increased complexity in the form of the phase field equations associated with general order parameters in actuality does not cause essential computational difficulties. Our numerical algorithm reformulates the (N−1) strongly-coupled phase field equations for general order parameters into 2(N−1) Helmholtz-type equations that are completely de-coupled from one another. This leads to a computational complexity comparable to that for the simplified phase field equations associated with certain special choice of the order parameters. We demonstrate the capabilities of the method developed herein using several test problems involving multiple fluid phases and large contrasts in densities and viscosities among the multitude of fluids. In particular, by comparing simulation results with the Langmuir–de Gennes theory of floating liquid lenses we show that the method using general order parameters produces physically accurate results for multiple fluid phases.
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Bansal, Gaurav; Mascarenhas, Ajith; Chen, Jacqueline H.
2014-10-01
In our paper, two- and three-dimensional direct numerical simulations (DNS) of autoignition phenomena in stratified dimethyl-ether (DME)/air turbulent mixtures are performed. A reduced DME oxidation mechanism, which was obtained using rigorous mathematical reduction and stiffness removal procedure from a detailed DME mechanism with 55 species, is used in the present DNS. The reduced DME mechanism consists of 30 chemical species. This study investigates the fundamental aspects of turbulence-mixing-autoignition interaction occurring in homogeneous charge compression ignition (HCCI) engine environments. A homogeneous isotropic turbulence spectrum is used to initialize the velocity field in the domain. Moreover, the computational configuration corresponds to amore » constant volume combustion vessel with inert mass source terms added to the governing equations to mimic the pressure rise due to piston motion, as present in practical engines. DME autoignition is found to be a complex three-staged process; each stage corresponds to a distinct chemical kinetic pathway. The distinct role of turbulence and reaction in generating scalar gradients and hence promoting molecular transport processes are investigated. Then, by applying numerical diagnostic techniques, the different heat release modes present in the igniting mixture are identified. In particular, the contribution of homogeneous autoignition, spontaneous ignition front propagation, and premixed deflagration towards the total heat release are quantified.« less
Bansal, Gaurav; Mascarenhas, Ajith; Chen, Jacqueline H.
2014-10-01
In our paper, two- and three-dimensional direct numerical simulations (DNS) of autoignition phenomena in stratified dimethyl-ether (DME)/air turbulent mixtures are performed. A reduced DME oxidation mechanism, which was obtained using rigorous mathematical reduction and stiffness removal procedure from a detailed DME mechanism with 55 species, is used in the present DNS. The reduced DME mechanism consists of 30 chemical species. This study investigates the fundamental aspects of turbulence-mixing-autoignition interaction occurring in homogeneous charge compression ignition (HCCI) engine environments. A homogeneous isotropic turbulence spectrum is used to initialize the velocity field in the domain. Moreover, the computational configuration corresponds to a constant volume combustion vessel with inert mass source terms added to the governing equations to mimic the pressure rise due to piston motion, as present in practical engines. DME autoignition is found to be a complex three-staged process; each stage corresponds to a distinct chemical kinetic pathway. The distinct role of turbulence and reaction in generating scalar gradients and hence promoting molecular transport processes are investigated. Then, by applying numerical diagnostic techniques, the different heat release modes present in the igniting mixture are identified. In particular, the contribution of homogeneous autoignition, spontaneous ignition front propagation, and premixed deflagration towards the total heat release are quantified.
Benchmark of numerical tools simulating beam propagation and secondary particles in ITER NBI
Sartori, E. Veltri, P.; Serianni, G.; Dlougach, E.; Hemsworth, R.; Singh, M.
2015-04-08
Injection of high energy beams of neutral particles is a method for plasma heating in fusion devices. The ITER injector, and its prototype MITICA (Megavolt ITER Injector and Concept Advancement), are large extrapolations from existing devices: therefore numerical modeling is needed to set thermo-mechanical requirements for all beam-facing components. As the power and charge deposition originates from several sources (primary beam, co-accelerated electrons, and secondary production by beam-gas, beam-surface, and electron-surface interaction), the beam propagation along the beam line is simulated by comprehensive 3D models. This paper presents a comparative study between two codes: BTR has been used for several years in the design of the ITER HNB/DNB components; SAMANTHA code was independently developed and includes additional phenomena, such as secondary particles generated by collision of beam particles with the background gas. The code comparison is valuable in the perspective of the upcoming experimental operations, in order to prepare a reliable numerical support to the interpretation of experimental measurements in the beam test facilities. The power density map calculated on the Electrostatic Residual Ion Dump (ERID) is the chosen benchmark, as it depends on the electric and magnetic fields as well as on the evolution of the beam species via interaction with the gas. Finally the paper shows additional results provided by SAMANTHA, like the secondary electrons produced by volume processes accelerated by the ERID fringe-field towards the Cryopumps.
Numerical construction and flow simulation in networks of fractures using fractals
Yortsos, Y.C.; Acuna, J.A.
1991-11-01
Present models for the representation of naturally fractured systems rely on the double-porosity Warren-Root model or on random arrays of fractures. However, field observation in outcrops has demonstrated the existence of multiple length scales in many naturally fractured media. The existing models fail to capture this important fractal property. In this paper, we use concepts from the theory of fragmentation and from fractal geometry for the numerical construction of networks of fractures that have fractal characteristics. The method is based mainly on the work of Barnsley (1) and allows for great flexibility in the development of patterns. Numerical techniques are developed for the simulation of unsteady single phase flow in such networks. It is found that the pressure transient response of finite fractals behaves according to the analytical predictions of Chang and Yortsos (6), provided that there exists a power law in the mass-radius relationship around the test well location. Otherwise, the finite size effects become significant and interfere severely with the identification of the underlying fractal structure. 21 refs., 13 figs.
THREE-DIMENSIONAL NUMERICAL SIMULATIONS OF FAST-TO-ALFVEN CONVERSION IN SUNSPOTS
Felipe, T.
2012-10-20
The conversion of fast waves to the Alfven mode in a realistic sunspot atmosphere is studied through three-dimensional numerical simulations. An upward propagating fast acoustic wave is excited in the high-{beta} region of the model. The new wave modes generated at the conversion layer are analyzed from the projections of the velocity and magnetic field in their characteristic directions, and the computation of their wave energy and fluxes. The analysis reveals that the maximum efficiency of the conversion to the slow mode is obtained for inclinations of 25 Degree-Sign and low azimuths, while the Alfven wave conversions peak at high inclinations and azimuths between 50 Degree-Sign and 120 Degree-Sign . Downward propagating Alfven waves appear at the regions of the sunspot where the orientation of the magnetic field is in the direction opposite to the wave propagation, since at these locations the Alfven wave couples better with the downgoing fast magnetic wave which is reflected due to the gradients of the Alfven speed. The simulations show that the Alfven energy at the chromosphere is comparable to the acoustic energy of the slow mode, being even higher at high inclined magnetic fields.
Terascale Direct Numerical Simulations of Turbulent Combustion: Capabilities and Limits (PReSS Talk)
Yoo, Chun Sang
2009-03-26
The rapid growth in computational capabilities has provided great opportunities for direct numerical simulations (DNS) of turbulent combustion, a type of simulations without any turbulence model. With the help of terascale high performance supercomputing (HPC) resources, we are now able to provide fundamental insight into turbulence-chemistry interaction in simple laboratory-scale turbulent flames with detailed chemistry using three-dimensional (3D) DNS. However, the actual domain size of 3D-DNS is still limited within {approx} O(10 cm{sup 3}) due to its tremendously high grid resolution required to resolve the smallest turbulent length scale as well as flame structures. Moreover, 3D-DNS will require more computing powers to investigate next-generation engines, of which operating conditions will be characterized by higher pressures, lower temperatures, and higher levels of dilution. In this talk, I will discuss the capabilities and limits of DNS of turbulent combustion and present some results of ignition/extinction characteristics of a highly diluted hydrogen flame counter-flowing against heated air. The results of our recent 3D-DNS of a spatially-developing turbulent lifted hydrogen jet flame in heated coflow will also be presented. The 3D-DNS was performed at a jet Reynolds number of 11,000 with {approx} 1 billion grid points, which required 3.5 million CPU hours on Cray XT3/XT4 at Oak Ridge National Laboratories.
GPU accelerated flow solver for direct numerical simulation of turbulent flows
Salvadore, Francesco; Botti, Michela
2013-02-15
Graphical processing units (GPUs), characterized by significant computing performance, are nowadays very appealing for the solution of computationally demanding tasks in a wide variety of scientific applications. However, to run on GPUs, existing codes need to be ported and optimized, a procedure which is not yet standardized and may require non trivial efforts, even to high-performance computing specialists. In the present paper we accurately describe the porting to CUDA (Compute Unified Device Architecture) of a finite-difference compressible Navier–Stokes solver, suitable for direct numerical simulation (DNS) of turbulent flows. Porting and validation processes are illustrated in detail, with emphasis on computational strategies and techniques that can be applied to overcome typical bottlenecks arising from the porting of common computational fluid dynamics solvers. We demonstrate that a careful optimization work is crucial to get the highest performance from GPU accelerators. The results show that the overall speedup of one NVIDIA Tesla S2070 GPU is approximately 22 compared with one AMD Opteron 2352 Barcelona chip and 11 compared with one Intel Xeon X5650 Westmere core. The potential of GPU devices in the simulation of unsteady three-dimensional turbulent flows is proved by performing a DNS of a spatially evolving compressible mixing layer.
Numerical simulation of laminar plasma dynamos in a cylindrical von Karman flow
Khalzov, I. V.; Brown, B. P.; Schnack, D. D.; Forest, C. B. [University of Wisconsin, 1150 University Avenue, Madison, Wisconsin 53706 (United States); Ebrahimi, F. [University of New Hampshire, 8 College Road, Durham, New Hampshire 03824 (United States)
2011-03-15
The results of a numerical study of the magnetic dynamo effect in cylindrical von Karman plasma flow are presented with parameters relevant to the Madison Plasma Couette Experiment. This experiment is designed to investigate a broad class of phenomena in flowing plasmas. In a plasma, the magnetic Prandtl number Pm can be of order unity (i.e., the fluid Reynolds number Re is comparable to the magnetic Reynolds number Rm). This is in contrast to liquid metal experiments, where Pm is small (so, Re>>Rm) and the flows are always turbulent. We explore dynamo action through simulations using the extended magnetohydrodynamic NIMROD code for an isothermal and compressible plasma model. We also study two-fluid effects in simulations by including the Hall term in Ohm's law. We find that the counter-rotating von Karman flow results in sustained dynamo action and the self-generation of magnetic field when the magnetic Reynolds number exceeds a critical value. For the plasma parameters of the experiment, this field saturates at an amplitude corresponding to a new stable equilibrium (a laminar dynamo). We show that compressibility in the plasma results in an increase of the critical magnetic Reynolds number, while inclusion of the Hall term in Ohm's law changes the amplitude of the saturated dynamo field but not the critical value for the onset of dynamo action.
Numerical simulations of the two-dimensional multimode Richtmyer-Meshkov instability
Thornber, B.; Zhou, Y.
2015-03-15
The two-dimensional Richtmyer-Meshkov instability occurs as shock waves pass through a perturbed material interface, triggering transition to an inhomogeneous turbulence variable density flow. This paper presents a series of large-eddy-simulations of the two dimensional turbulent RM instability and compares the results to the fully three dimensional simulations. There are two aims for this paper, the first is to explore what numerical resolution is required for a statistically converged solution for a two dimensional inhomogeneous flow field. The second aim is to elucidate the key differences in flow physics between the two dimensional and three dimensional Richtmyer-Meshkov instabilities, particularly their asymptotic self-similar regime. Convergence is achieved using 64 independent realisations and grid resolutions up to 4096{sup 2} in the plane. It is shown that for narrowband cases the growth rate θ = 0.48 which is substantially higher than the three-dimensional equivalent. Mix measures are consistently lower compared to three-dimensional, and the kinetic energy distribution is homogeneous at late time. The broadband case has a similar initial growth rate as the three-dimensional case, with a marginally lower θ = 0.63. Mix is similar in magnitude, but is reducing at late time. The spectra in both cases exhibit the dual-cascade expected from two-dimensional turbulence.
Numerical simulation of viscoelastic layer rearrangement in polymer melts using OpenFOAM®
Köpplmayr, Thomas Mayrhofer, Elias
2015-05-22
In addition to their shear-thinning behavior, polymer melts are characterized by first and second normal stress differences, which cause secondary motions. Polymer coextrusion processes involve viscoelastic two-phase flows that influence layer formation. Using polymer melts with different pigmentation makes visible the layers deformed by second normal stress differences. We used a new solver for the OpenFOAM CFD toolbox which handles viscoelastic two-phase flows. A derivative of the volume-of-fluid (VoF) methodology was employed to describe the interface. Different types of polymer melt, such as polyethylene (PE), polypropylene (PP) and polyethylene terephthalate (PET) were investigated. In a coextrusion process, the less viscous phase usually tends to encapsulate the more viscous one. However, the different viscoelastic properties of the melts also influence interface deformation. The materials were characterized by small-amplitude oscillatory-shear rheometry, and a multimode Giesekus model was used to fit shear viscosity, storage and loss modulus. Our simulations also took interfacial tension into account. Experimental observations and corresponding numerical simulations were found to be in good accordance.
Zhao Xinghai; Mathews, Grant J.
2011-01-15
General relativistic corrections to the expansion rate of the Universe arise when the Einstein equations are averaged over a spatial volume in a locally inhomogeneous cosmology. It has been suggested that they may contribute to the observed cosmic acceleration. In this paper, we propose a new scheme that utilizes numerical simulations to make a realistic estimate of the magnitude of these corrections for general inhomogeneities in (3+1) spacetime. We then quantitatively calculate the volume averaged expansion rate using N-body large-scale structure simulations and compare it with the expansion rate in a standard FRW cosmology. We find that in the weak gravitational field limit, the converged corrections are slightly larger than the previous claimed 10{sup -5} level, but not large enough nor even of the correct sign to drive the current cosmic acceleration. Nevertheless, the question of whether the cumulative effect can significantly change the expansion history of the Universe needs to be further investigated with strong-field relativity.
A review of direct numerical simulations of astrophysical detonations and their implications
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Parete-Koon, Suzanne T.; Smith, Christopher R.; Papatheodore, Thomas L.; Bronson Messer, O. E.
2013-04-11
Multi-dimensional direct numerical simulations (DNS) of astrophysical detonations in degenerate matter have revealed that the nuclear burning is typically characterized by cellular structure caused by transverse instabilities in the detonation front. Type Ia supernova modelers often use one- dimensional DNS of detonations as inputs or constraints for their whole star simulations. While these one-dimensional studies are useful tools, the true nature of the detonation is multi-dimensional. The multi-dimensional structure of the burning influences the speed, stability, and the composition of the detonation and its burning products, and therefore, could have an impact on the spectra of Type Ia supernovae. Considerablemore » effort has been expended modeling Type Ia supernovae at densities above 1x107 g∙cm-3 where the complexities of turbulent burning dominate the flame propagation. However, most full star models turn the nuclear burning schemes off when the density falls below 1x107 g∙cm-3 and distributed burning begins. The deflagration to detonation transition (DDT) is believed to occur at just these densities and consequently they are the densities important for studying the properties of the subsequent detonation. In conclusion, this work reviews the status of DNS studies of detonations and their possible implications for Type Ia supernova models. It will cover the development of Detonation theory from the first simple Chapman-Jouguet (CJ) detonation models to the current models based on the time-dependent, compressible, reactive flow Euler equations of fluid dynamics.« less
Numerical simulation of jet mixing concepts in Tank 241-SY-101
Trent, D.S.; Michener, T.E.
1993-03-01
The episodic gas release events (GRES) that have characterized the behavior of Tank 241-SY-101 for the past several years are thought to result from gases generated by the waste material in it that become trapped in the layer of settled solids at the bottom of the tank. Several concepts for mitigating the GREs have been proposed. One concept involves mobilizing the solid particles with mixing jets. The rationale behind this idea is to prevent formation of a consolidated layer of settled solids at the bottom of the tank, thus inhibiting the accumulation of gas bubbles in this layer. Numerical simulations were conducted using the TEMPEST computer code to assess the viability and effectiveness of the proposed jet discharge concepts and operating parameters. Before these parametric studies were commenced, a series of turbulent jet studies were conducted that established the adequacy of the TEMPEST code for this application. Configurations studied for Tank 241-SY-101 include centrally located downward discharging jets, draft tubes, and horizontal jets that are either stationary or rotating. Parameter studies included varying the jet discharge velocity, jet diameter, discharge elevation, and material properties. A total of 18 simulations were conducted and are reported in this document. The effect of gas bubbles on the mixing dynamics was not included within the scope of this study.
Numerical simulation of the compressor coil of the plasma dynamic accelerator
Thomas, P.
1997-01-01
The plasma dynamic accelerator accelerates a plasma to very high velocities in a coaxial accelerator and then compresses it in a compressor coil, achieving high densities. The axial component of the current distribution, extending from the tip of the coaxial accelerator`s center electrode to the coil turns, causes compressing forces, the radial component yields accelerating forces. The rapid change of the coil current induces azimuthal eddy currents in the plasma that interact with the coil`s magnetic field, again yielding Lorentz forces. Aerodynamic compression may also be an important effect. A new two-dimensional magnetohydrodynamics code is used to investigate which of these effects are really important for the compression. The code allows one to simulate all effects mentioned separately and in combination. In a first step only aerodynamic compression is considered. Then each electromagnetic effect is imposed on the system. Finally, a complete simulation of the compressor coil is performed. The analysis of the results provides new insights in the way the coil operates. This paper presents important aspects of the mathematical model and of the numerical implementation and reports results.
TOUGH2: A general-purpose numerical simulator for multiphase nonisothermal flows
Pruess, K.
1991-06-01
Numerical simulators for multiphase fluid and heat flows in permeable media have been under development at Lawrence Berkeley Laboratory for more than 10 yr. Real geofluids contain noncondensible gases and dissolved solids in addition to water, and the desire to model such `compositional` systems led to the development of a flexible multicomponent, multiphase simulation architecture known as MULKOM. The design of MULKOM was based on the recognition that the mass-and energy-balance equations for multiphase fluid and heat flows in multicomponent systems have the same mathematical form, regardless of the number and nature of fluid components and phases present. Application of MULKOM to different fluid mixtures, such as water and air, or water, oil, and gas, is possible by means of appropriate `equation-of-state` (EOS) modules, which provide all thermophysical and transport parameters of the fluid mixture and the permeable medium as a function of a suitable set of primary thermodynamic variables. Investigations of thermal and hydrologic effects from emplacement of heat-generating nuclear wastes into partially water-saturated formations prompted the development and release of a specialized version of MULKOM for nonisothermal flow of water and air, named TOUGH. TOUGH is an acronym for `transport of unsaturated groundwater and heat` and is also an allusion to the tuff formations at Yucca Mountain, Nevada. The TOUGH2 code is intended to supersede TOUGH. It offers all the capabilities of TOUGH and includes a considerably more general subset of MULKOM modules with added capabilities. The paper briefly describes the simulation methodology and user features.
Yang, Xiaofan; Scheibe, Timothy D.; Richmond, Marshall C.; Perkins, William A.; Vogt, Sarah J.; Codd, Sarah L.; Seymour, Joseph D.; Mckinley, Matthew I.
2013-04-01
A significant body of current research is aimed at developing methods for numerical simulation of flow and transport in porous media that explicitly resolve complex pore and solid geometries, and at utilizing such models to study the relationships between fundamental pore-scale processes and macroscopic manifestations at larger (i.e., Darcy) scales. A number of different numerical methods for pore-scale simulation have been developed, and have been extensively tested and validated for simplified geometries. However, validation of pore-scale simulations of fluid velocity for complex, three-dimensional (3D) pore geometries that are representative of natural porous media is challenging due to our limited ability to measure pore-scale velocity in such systems. Recent advances in magnetic resonance imaging (MRI) offer the opportunity to measure not only the pore geometry, but also local fluid velocities under steady-state flow conditions in 3D and with high spatial resolution. In this paper, we present a 3D velocity field measured at sub-pore resolution (tens of micrometers) over a centimeter-scale 3D domain using MRI methods. We have utilized the measured pore geometry to perform 3D simulations of Navier-Stokes flow over the same domain using direct numerical simulation techniques. We present a comparison of the numerical simulation results with the measured velocity field. It is shown that the numerical results match the observed velocity patterns well overall except for a variance and small systematic scaling which can be attributed to the known experimental error in the MRI measurements. The comparisons presented here provide strong validation of the pore-scale simulation methods and new insights for interpretation of uncertainty in MRI measurements of pore-scale velocity. This study also provides a potential benchmark for future comparison of other pore-scale simulation methods.
Numerical simulations of quiet sun magnetism: On the contribution from a small-scale dynamo
Rempel, M.
2014-07-10
We present a series of radiative MHD simulations addressing the origin and distribution of the mixed polarity magnetic field in the solar photosphere. To this end, we consider numerical simulations that cover the uppermost 2-6 Mm of the solar convection zone and we explore scales ranging from 2 km to 25 Mm. We study how the strength and distribution of the magnetic field in the photosphere and subsurface layers depend on resolution, domain size, and boundary conditions. We find that 50% of the magnetic energy at the τ = 1 level comes from fields with the less than 500 G strength and that 50% of the energy resides on scales smaller than about 100 km. While the probability distribution functions are essentially independent of resolution, properly describing the spectral energy distribution requires grid spacings of 8 km or smaller. The formation of flux concentrations in the photosphere exceeding 1 kG requires a mean vertical field strength greater than 30-40 G at τ = 1. The filling factor of kG flux concentrations increases with overall domain size as the magnetic field becomes organized by larger, longer-lived flow structures. A solution with a mean vertical field strength of around 85 G at τ = 1 requires a subsurface rms field strength increasing with depth at the same rate as the equipartition field strength. We consider this an upper limit for the quiet Sun field strength, which implies that most of the convection zone is magnetized close to the equipartition. We discuss these findings in view of recent high-resolution spectropolarimetric observations of quiet Sun magnetism.
Numerical simulation of solar heat absorption within indoor space by means of composite grid method
Omori, Toshiaki; Murakami, Shuzo; Kato, Shinsuke
1997-12-31
This paper describes the method for numerical simulation of solar radiation entering indoor spaces through fenestration. The proposed method can systematically deal with the interception of sunlight by buildings in the outdoor space and obstacles in the indoor space by tracing a large number of particles directed toward the sun. Configuration factors from the fenestration to the sky are also three-dimensionally treated by accounting for outdoor geometries. Distribution of the solar heat absorption in the indoor space is calculated by assuming radiation equilibrium. After the solar heat absorption analysis is carried out, heat transfer analysis in the space is conducted taking account of longwave radiation, convective heat transfer, thermal conduction, and cooling/heating by air conditioning. Then, the indoor thermal environment is evaluated using the resulting temperature distribution of air and indoor surfaces. To evaluate the applicability of these procedures, the thermal environment in a model hall with large glass windows and an overhang is predicted. The analyzed hall is assumed to be located near a tall building.
Pan Yi; Buonanno, Alessandra; Buchman, Luisa T.; Chu, Tony; Scheel, Mark A.; Kidder, Lawrence E.; Pfeiffer, Harald P.
2010-04-15
We present the first attempt at calibrating the effective-one-body (EOB) model to accurate numerical relativity simulations of spinning, nonprecessing black-hole binaries. Aligning the EOB and numerical waveforms at low frequency over a time interval of 1000M, we first estimate the phase and amplitude errors in the numerical waveforms and then minimize the difference between numerical and EOB waveforms by calibrating a handful of EOB-adjustable parameters. In the equal-mass, spin aligned case, we find that phase and fractional amplitude differences between the numerical and EOB (2,2) mode can be reduced to 0.01 radian and 1%, respectively, over the entire inspiral waveforms. In the equal-mass, spin antialigned case, these differences can be reduced to 0.13 radian and 1% during inspiral and plunge, and to 0.4 radian and 10% during merger and ringdown. The waveform agreement is within numerical errors in the spin aligned case while slightly over numerical errors in the spin antialigned case. Using Enhanced LIGO and Advanced LIGO noise curves, we find that the overlap between the EOB and the numerical (2,2) mode, maximized over the initial phase and time of arrival, is larger than 0.999 for binaries with total mass 30M{sub {center_dot}-}200M{sub {center_dot}}. In addition to the leading (2,2) mode, we compare four subleading modes. We find good amplitude and frequency agreements between the EOB and numerical modes for both spin configurations considered, except for the (3,2) mode in the spin antialigned case. We believe that the larger difference in the (3,2) mode is due to the lack of knowledge of post-Newtonian spin effects in the higher modes.
Direct Numerical Simulation of Interfacial Flows: Implicit Sharp-Interface Method (I-SIM)
Robert Nourgaliev; Theo Theofanous; HyeongKae Park; Vincent Mousseau; Dana Knoll
2008-01-01
In recent work (Nourgaliev, Liou, Theofanous, JCP in press) we demonstrated that numerical simulations of interfacial flows in the presence of strong shear must be cast in dynamically sharp terms (sharp interface treatment or SIM), and that moreover they must meet stringent resolution requirements (i.e., resolving the critical layer). The present work is an outgrowth of that work aiming to overcome consequent limitations on the temporal treatment, which become still more severe in the presence of phase change. The key is to avoid operator splitting between interface motion, fluid convection, viscous/heat diffusion and reactions; instead treating all these non-linear operators fully-coupled within a Newton iteration scheme. To this end, the SIMs cut-cell meshing is combined with the high-orderaccurate implicit Runge-Kutta and the recovery Discontinuous Galerkin methods along with a Jacobian-free, Krylov subspace iteration algorithm and its physics-based preconditioning. In particular, the interfacial geometry (i.e., markers positions and volumes of cut cells) is a part of the Newton-Krylov solution vector, so that the interface dynamics and fluid motions are fully-(non-linearly)-coupled. We show that our method is: (a) robust (L-stable) and efficient, allowing to step over stability time steps at will while maintaining high-(up to the 5th)-order temporal accuracy; (b) fully conservative, even near multimaterial contacts, without any adverse consequences (pressure/velocity oscillations); and (c) highorder-accurate in spatial discretization (demonstrated here up to the 12th-order for smoothin-the-bulk-fluid flows), capturing interfacial jumps sharply, within one cell. Performance is illustrated with a variety of test problems, including low-Mach-number manufactured solutions, shock dynamics/tracking with slow dynamic time scales, and multi-fluid, highspeed shock-tube problems. We briefly discuss preconditioning, and we introduce two physics-based preconditioners
Daeva, S.G.; Setukha, A.V.
2015-03-10
A numerical method for solving a problem of diffraction of acoustic waves by system of solid and thin objects based on the reduction the problem to a boundary integral equation in which the integral is understood in the sense of finite Hadamard value is proposed. To solve this equation we applied piecewise constant approximations and collocation methods numerical scheme. The difference between the constructed scheme and earlier known is in obtaining approximate analytical expressions to appearing system of linear equations coefficients by separating the main part of the kernel integral operator. The proposed numerical scheme is tested on the solution of the model problem of diffraction of an acoustic wave by inelastic sphere.
Goldberg, L.F.
1992-04-01
Aspects of the information propagation modeling behavior of integral machine computer simulation programs are investigated in terms of a transmission line. In particular, the effects of pressure-linking and temporal integration algorithms on the amplitude ratio and phase angle predictions are compared against experimental and closed-form analytic data. It is concluded that the discretized, first order conservation balances may not be adequate for modeling information propagation effects at characteristic numbers less than about 24. An entropy transport equation suitable for generalized use in Stirling machine simulation is developed. The equation is evaluated by including it in a simulation of an incompressible oscillating flow apparatus designed to demonstrate the effect of flow oscillations on the enhancement of thermal diffusion. Numerical false diffusion is found to be a major factor inhibiting validation of the simulation predictions with experimental and closed-form analytic data. A generalized false diffusion correction algorithm is developed which allows the numerical results to match their analytic counterparts. Under these conditions, the simulation yields entropy predictions which satisfy Clausius' inequality.
Martinez-Tossas, Luis A.; Churchfield, Matthew J.; Meneveau, Charles
2015-06-18
In this work we report on results from a detailed comparative numerical study from two Large Eddy Simulation (LES) codes using the Actuator Line Model (ALM). The study focuses on prediction of wind turbine wakes and their breakdown when subject to uniform inflow. Previous studies have shown relative insensitivity to subgrid modeling in the context of a finite-volume code. The present study uses the low dissipation pseudo-spectral LES code from Johns Hopkins University (LESGO) and the second-order, finite-volume OpenFOAMcode (SOWFA) from the National Renewable Energy Laboratory. When subject to uniform inflow, the loads on the blades are found to be unaffected by subgrid models or numerics, as expected. The turbulence in the wake and the location of transition to a turbulent state are affected by the subgrid-scale model and the numerics.
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Martinez-Tossas, Luis A.; Churchfield, Matthew J.; Meneveau, Charles
2015-06-18
In this work we report on results from a detailed comparative numerical study from two Large Eddy Simulation (LES) codes using the Actuator Line Model (ALM). The study focuses on prediction of wind turbine wakes and their breakdown when subject to uniform inflow. Previous studies have shown relative insensitivity to subgrid modeling in the context of a finite-volume code. The present study uses the low dissipation pseudo-spectral LES code from Johns Hopkins University (LESGO) and the second-order, finite-volume OpenFOAMcode (SOWFA) from the National Renewable Energy Laboratory. When subject to uniform inflow, the loads on the blades are found to bemore » unaffected by subgrid models or numerics, as expected. The turbulence in the wake and the location of transition to a turbulent state are affected by the subgrid-scale model and the numerics.« less
Andrea Prosperetti
2006-03-24
The report briefly describes the activities carried out in the course of the project. A first line of research was the development of systematic closure relations for averaged equations for disperse multiphase flow. A second line was the development of efficient numerical methods for the simulation of Navier-Stokes flows with many suspended particles. The report also lists the 21 journal articles in which this work is more fully decsribed.
Ohshima, Hiroyuki; Uwaba, Tomoyuki; Hashimoto, Akihiko; Imai, Yasutomo; Ito, Masahiro
2015-12-31
A numerical simulation system, which consists of a deformation analysis program and three kinds of thermal-hydraulics analysis programs, is being developed in Japan Atomic Energy Agency in order to offer methodologies to clarify thermal-hydraulic phenomena in fuel assemblies of sodium-cooled fast reactors under various operating conditions. This paper gives the outline of the system and its applications to fuel assembly analyses as a validation study.
Not Available
2012-02-01
New code will help accelerate design improvements by providing a high-fidelity simulation tool to study power performance, structural loading, and the interactions between devices in arrays.
Modeling-Computer Simulations At Nw Basin & Range Region (Laney...
previous seismic experiments and earthquake-monitoring projects, and data donated from mining, geothermal, and petroleum companies. We also collected (May 2002 and August 2004) two...
Modeling-Computer Simulations At Nw Basin & Range Region (Blackwell...
generic Basin and Range systems based on Dixie Valley data that help to understand the nature of large scale constraints on the location and characteristics of the geothermal...
Modeling-Computer Simulations At San Juan Volcanic Field Area...
San Juan region, to further investigate both the thermal history of the region and the nature of the influence of the San Juan volcanic field thermal source on the thermal history...
Modeling-Computer Simulations At Long Valley Caldera Geothermal...
(2003) The Mechanics of Unrest at Long Valley Caldera, California. 2. Constraining the Nature of the Source Using Geodetic and Micro-Gravity Data John O. Langbein (2003)...
Modeling-Computer Simulations At Central Nevada Seismic Zone...
generic Basin and Range systems based on Dixie Valley data that help to understand the nature of large scale constraints on the location and characteristics of the geothermal...
Modeling-Computer Simulations At Northern Basin & Range Region...
generic Basin and Range systems based on Dixie Valley data that help to understand the nature of large scale constraints on the location and characteristics of the geothermal...
Modeling-Computer Simulations (Laney, 2005) | Open Energy Information
in the near surface: Available technologies for monitoring CO2 in the near-surface environment include (1) the infrared gas analyzer (IRGA) for measurement of concentrations at...
Modeling-Computer Simulations At Dixie Valley Geothermal Area...
vein structure associated with ore deposits. References David D. Blackwell, Richard P. Smith, Al Waibel, Maria C. Richards, Patrick Stepp (2009) Why Basin and Range Systems are...
Modeling-Computer Simulations (Gritto & Majer) | Open Energy...
are shown in Figure 1. The parameters of the fault were modeled after Coates and Schoenberg (1995), where the orientation of the fault relative to the finite-difference grid...
Modeling-Computer Simulations (Walker, Et Al., 2005) | Open Energy...
occurrence model for geothermal systems based on fundamental geologic data. References J. D. Walker, A. E. Sabin, J. R. Unruh, J. Combs, F. C. Monastero (2005) Development Of...
Modeling-Computer Simulations At Kilauea East Rift Geothermal...
importance of water convection for distributing heat in the East Rift Zone. References Albert J. Rudman, David Epp (1983) Conduction Models Of The Temperature Distribution In The...
Numerical simulation of transient, incongruent vaporization induced by high power laser
Tsai, C.H.
1981-01-01
A mathematical model and numerical calculations were developed to solve the heat and mass transfer problems specifically for uranum oxide subject to laser irradiation. It can easily be modified for other heat sources or/and other materials. In the uranium-oxygen system, oxygen is the preferentially vaporizing component, and as a result of the finite mobility of oxygen in the solid, an oxygen deficiency is set up near the surface. Because of the bivariant behavior of uranium oxide, the heat transfer problem and the oxygen diffusion problem are coupled and a numerical method of simultaneously solving the two boundary value problems is studied. The temperature dependence of the thermal properties and oxygen diffusivity, as well as the highly ablative effect on the surface, leads to considerable non-linearities in both the governing differential equations and the boundary conditions. Based on the earlier work done in this laboratory by Olstad and Olander on Iron and on Zirconium hydride, the generality of the problem is expanded and the efficiency of the numerical scheme is improved. The finite difference method, along with some advanced numerical techniques, is found to be an efficient way to solve this problem.
Numerical simulation of the shock-tip leakage vortex interaction in a HPC front stage
Hoeger, M.; Fritsch, G.; Bauer, D.
1999-07-01
For a single-stage transonic compressor rig at the TU Darmstadt, three-dimensional viscous simulations are compared to L2F measurements and data from the EGV leading edge instrumentation to demonstrate the predictive capability of the Navier-Stokes code TRACE{_}S. In a second step the separated regions at the blade tip are investigated in detail to gain insight into the mechanisms of tip leakage vortex-shock interaction at operating points close to stall, peak efficiency, and choke. At the casing the simulations reveal a region with axially reversed flow, leading to a rotationally asymmetric displacement of the outermost stream surface and a localized additional pitch-average blockage of approximately 2 percent. Loss mechanisms and streamline patterns deduced from the simulation are also discussed. Although the flow is essentially three-dimensional, a simple model for local blockage from tip leakage is demonstrated to significantly improve two-dimensional simulations on S1-surfaces.
Numerical Simulation of Combustion and Rotor-Stator Interaction in a Turbine Combustor
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Isvoranu, Dragos D.; Cizmas, Paul G. A.
2003-01-01
This article presents the development of a numerical algorithm for the computation of flow and combustion in a turbine combustor. The flow and combustion are modeled by the Reynolds-averaged Navier-Stokes equations coupled with the species-conservation equations. The chemistry model used herein is a two-step, global, finite-rate combustion model for methane and combustion gases. The governing equations are written in the strong conservation form and solved using a fully implicit, finite-difference approximation. The gas dynamics and chemistry equations are fully decoupled. A correction technique has been developed to enforce the conservation of mass fractions. The numerical algorithm developed herein has beenmore » used to investigate the flow and combustion in a one-stage turbine combustor.« less
Lavergne, F.; Sab, K.; Sanahuja, J.; Bornert, M.; Toulemonde, C.
2015-05-15
Prestress losses due to creep of concrete is a matter of interest for long-term operations of nuclear power plants containment buildings. Experimental studies by Granger (1995) have shown that concretes with similar formulations have different creep behaviors. The aim of this paper is to numerically investigate the effect of size distribution and shape of elastic inclusions on the long-term creep of concrete. Several microstructures with prescribed size distribution and spherical or polyhedral shape of inclusions are generated. By using the 3D numerical homogenization procedure for viscoelastic microstructures proposed by Šmilauer and Bažant (2010), it is shown that the size distribution and shape of inclusions have no measurable influence on the overall creep behavior. Moreover, a mean-field estimate provides close predictions. An Interfacial Transition Zone was introduced according to the model of Nadeau (2003). It is shown that this feature of concrete's microstructure can explain differences between creep behaviors.
Numerical simulation: Toward the design of high-efficiency planar perovskite solar cells
Liu, Feng; Zhu, Jun E-mail: sydai@ipp.ac.cn; Wei, Junfeng; Li, Yi; Lv, Mei; Yang, Shangfeng; Zhang, Bing; Yao, Jianxi; Dai, Songyuan E-mail: sydai@ipp.ac.cn
2014-06-23
Organo-metal halide perovskite solar cells based on planar architecture have been reported to achieve remarkably high power conversion efficiency (PCE, >16%), rendering them highly competitive to the conventional silicon based solar cells. A thorough understanding of the role of each component in solar cells and their effects as a whole is still required for further improvement in PCE. In this work, the planar heterojunction-based perovskite solar cells were simulated with the program AMPS (analysis of microelectronic and photonic structures)-1D. Simulation results revealed a great dependence of PCE on the thickness and defect density of the perovskite layer. Meanwhile, parameters including the work function of the back contact as well as the hole mobility and acceptor density in hole transport materials were identified to significantly influence the performance of the device. Strikingly, an efficiency over 20% was obtained under the moderate simulation conditions.
Pore-Scale and Multiscale Numerical Simulation of Flow and Transport in a Laboratory-Scale Column
Scheibe, Timothy D.; Perkins, William A.; Richmond, Marshall C.; McKinley, Matthey I.; Romero Gomez, Pedro DJ; Oostrom, Martinus; Wietsma, Thomas W.; Serkowski, John A.; Zachara, John M.
2015-02-01
Pore-scale models are useful for studying relationships between fundamental processes and phenomena at larger (i.e., Darcy) scales. However, the size of domains that can be simulated with explicit pore-scale resolution is limited by computational and observational constraints. Direct numerical simulation of pore-scale flow and transport is typically performed on millimeter-scale volumes at which X-ray computed tomography (XCT), often used to characterize pore geometry, can achieve micrometer resolution. In contrast, the scale at which a continuum approximation of a porous medium is valid is usually larger, on the order of centimeters to decimeters. Furthermore, laboratory experiments that measure continuum properties are typically performed on decimeter-scale columns. At this scale, XCT resolution is coarse (tens to hundreds of micrometers) and prohibits characterization of small pores and grains. We performed simulations of pore-scale processes over a decimeter-scale volume of natural porous media with a wide range of grain sizes, and compared to results of column experiments using the same sample. Simulations were conducted using high-performance codes executed on a supercomputer. Two approaches to XCT image segmentation were evaluated, a binary (pores and solids) segmentation and a ternary segmentation that resolved a third category (porous solids with pores smaller than the imaged resolution). We used a mixed Stokes-Darcy simulation method to simulate the combination of Stokes flow in large open pores and Darcy-like flow in porous solid regions. Simulations based on the ternary segmentation provided results that were consistent with experimental observations, demonstrating our ability to successfully model pore-scale flow over a column-scale domain.
NUMERICAL SIMULATIONS OF CHROMOSPHERIC ANEMONE JETS ASSOCIATED WITH MOVING MAGNETIC FEATURES
Yang, Liping; He, Jiansen; Tu, Chuanyi; Zhang, Lei; Peter, Hardi; Feng, Xueshang; Zhang, Shaohua
2013-11-01
Observations with the space-based solar observatory Hinode show that small-scale magnetic structures in the photosphere are found to be associated with a particular class of jets of plasma in the chromosphere called anemone jets. The goal of our study is to conduct a numerical experiment of such chromospheric anemone jets related to the moving magnetic features (MMFs). We construct a 2.5 dimensional numerical MHD model to describe the process of magnetic reconnection between the MMFs and the pre-existing ambient magnetic field, which is driven by the horizontal motion of the magnetic structure in the photosphere. We include thermal conduction parallel to the magnetic field and optically thin radiative losses in the corona to account for a self-consistent description of the evaporation process during the heating of the plasma due to the reconnection process. The motion of the MMFs leads to the expected jet and our numerical results can reproduce many observed characteristics of chromospheric anemone jets, topologically and quantitatively. As a result of the tearing instability, plasmoids are generated in the reconnection process that are consistent with the observed bright moving blobs in the anemone jets. An increase in the thermal pressure at the base of the jet is also driven by the reconnection, which induces a train of slow-mode shocks propagating upward. These shocks are a secondary effect, and only modulate the outflow of the anemone jet. The jet itself is driven by the energy input due to the reconnection of the MMFs and the ambient magnetic field.
Martin, V.; Egido, J.L.; Khoo, T.L.; Lauritsen, T.
1993-11-11
The electromagnetic decay of the nuclei {sup 152-154-156}Dy is analyzed using microscopic Hartree-Fock calculations at finite temperature. The theoretical collective transition probabilities are implemented in numerical simulations to produce theoretical espectra. Thermal shape fluctuations are also taken into account. The inclusion of these correlation is crucial in order to understand the main features of the collective E2 spectra of these isotopes at different energies. The theoretical calculations suggest a shape change as responsible for the unusual features of the spectrum of the nucleus {sup 154}Dy at high energy.
Karakulov, Valerii V.; Smolin, Igor Yu. E-mail: skrp@ftf.tsu.ru; Skripnyak, Vladimir A. E-mail: skrp@ftf.tsu.ru
2014-11-14
Mechanical behavior of stochastic metal-ceramic composites with the aluminum matrix under high-rate deformation at shock-wave loading is numerically simulated with consideration for structural evolution. Effective values of mechanical parameters of metal-ceramic composites Al
Oliva, A.; Costa, M.; Perez Segarra, C.D. )
1991-01-01
A numerical model has been developed for determination of thermal behavior of solar collector. The model takes into account the multidimensional and transient aspects that characterize the phenomenon of heat transfer in a collector. The modelization carried out allows the analysis of the influence of such aspects as: flow nonuniformity distribution, areas of shadow, and variations in dimension and properties of the different elements. These aspects can be analyzed equally for steady and nonsteady outdoor conditions. Illustrative situations of the influence on the collector performance of the different aspects previously mentioned are shown.
Soker, N.; Sarazin, C.L.; O'Dea, C.P.
1988-04-01
Three-dimensional numerical hydrodynamic simulations are used to study the bending of radio jets. The simulations are compared with observations of jets in narrow-angle-tail radio sources. Two mechanisms for the observed bending are considered: direct bending of quasi-continuous jets by ram pressure from intergalactic gas and bending by pressure gradients in the interstellar gas of the host galaxy, the pressure gradients themselves being the result of ram pressure by intergalactic gas. It is shown that the pressure gradients are much less effective in bending jets, implying that the jets have roughly 30 times lower momentum fluxes if they are bent by this mechanism. Ram-pressure bending produces jets with kidney-shaped cross sections; when observed from the side, these jets appear to have diffuse extensions on the downstream side. On the other hand, pressure-gradient bending causes the jets to be densest near their upstream side. 31 references.
Mazzeo, Brian A.; Patil, Anjali N.; Klis, Jeffrey M.; Hurd, Randy C.; Truscott, Tadd T.; Guthrie, W. Spencer
2014-02-18
Delaminations in bridge decks typically result from corrosion of the top mat of reinforcing steel, which leads to a localized separation of the concrete cover from the underlying concrete. Because delaminations cannot be detected using visual inspection, rapid, large-area interrogation methods are desired to characterize bridge decks without disruption to traffic, without the subjectivity inherent in existing methods, and with increased inspector safety. To this end, disposable impactors such as water droplets or ice chips can be dropped using automatic dispensers onto concrete surfaces to excite mechanical vibrations while acoustic responses can be recorded using air-coupled microphones. In this work, numerical simulations are used to characterize the flexural response of a model concrete bridge deck subject to both steel and ice impactors, and the results are compared with similar experiments performed in the laboratory on a partially delaminated concrete bridge deck slab. The simulations offer greater understanding of the kinetics of impacts and the responses of materials.
Final Report: A Model Management System for Numerical Simulations of Subsurface Processes
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.
Egorov, I.
2014-06-15
This paper describes the development of a computation model of a pulsed voltage generator for a repetitive electron accelerator. The model is based on a principle circuit of the generator, supplemented with the parasitics elements of the construction. Verification of the principle model was achieved by comparison of simulation with experimental results, where reasonable agreement was demonstrated for a wide range of generator load resistance.
Sergeeva, Ekaterina A; Katichev, A R; Kirillin, M Yu
2011-01-24
Using the radiative transfer theory and Monte Carlo simulations, we analyse the effect of scattering in a medium and of the size of the detector pinhole on the formation of the fluorescent signal in standard two-photon fluorescence microscopy (TPFM) systems. The theoretical analysis is based on a small-angle diffusion approximation of the radiative transfer equation, adapted to calculate the propagation of focused infrared radiation in media similar to the biological tissues in their optical properties. The accuracy of the model is evaluated by comparing the calculated excitation intensity in a highly scattering medium with the results of Monte Carlo simulations. To simulate a tightly focused Gaussian beam by the Monte Carlo method, the so called 'ray-optics' approach that correctly takes into account the finite size and shape of the beam waist is applied. It is shown that in the combined confocal and two-photon scanning microscopy systems not equipped with an external 'nondescanned' detector, the scattering significantly affects both the nonlinear excitation efficiency in the medium and the fluorescence collection efficiency of the system. In such systems, the rate of the useful TPFM signal in-depth decay is 1.5 - 2 times higher than in systems equipped with a 'nondescanned' detector. (application of lasers and laser-optical methods in life sciences)
Pahn, T.; Jonkman, J.; Rolges, R.; Robertson, A.
2012-11-01
Physically measuring the dynamic responses of wind turbine support structures enables the calculation of the applied loads using an inverse procedure. In this process, inverse means deriving the inputs/forces from the outputs/responses. This paper presents results of a numerical verification of such an inverse load calculation. For this verification, the comprehensive simulation code FAST is used. FAST accounts for the coupled dynamics of wind inflow, aerodynamics, elasticity and turbine controls. Simulations are run using a 5-MW onshore wind turbine model with a tubular tower. Both the applied loads due to the instantaneous wind field and the resulting system responses are known from the simulations. Using the system responses as inputs to the inverse calculation, the applied loads are calculated, which in this case are the rotor thrust forces. These forces are compared to the rotor thrust forces known from the FAST simulations. The results of these comparisons are presented to assess the accuracy of the inverse calculation. To study the influences of turbine controls, load cases in normal operation between cut-in and rated wind speed, near rated wind speed and between rated and cut-out wind speed are chosen. The presented study shows that the inverse load calculation is capable of computing very good estimates of the rotor thrust. The accuracy of the inverse calculation does not depend on the control activity of the wind turbine.
Perez, Jean Carlos; Chandran, Benjamin D. G.
2013-10-20
We present direct numerical simulations of inhomogeneous reduced magnetohydrodynamic (RMHD) turbulence between the Sun and the Alfvn critical point. These are the first such simulations that take into account the solar-wind outflow velocity and the radial inhomogeneity of the background solar wind without approximating the nonlinear terms in the governing equations. RMHD turbulence is driven by outward-propagating Alfvn waves (z {sup +} fluctuations) launched from the Sun, which undergo partial non-WKB reflection to produce sunward-propagating Alfvn waves (z {sup } fluctuations). We present 10 simulations with different values of the correlation time ?{sub c{sub sun}{sup +}} and perpendicular correlation length L{sub ?} of outward-propagating Alfvn waves at the coronal base. We find that between 15% and 33% of the z {sup +} energy launched into the corona dissipates between the coronal base and Alfvn critical point. Between 33% and 40% of this input energy goes into work on the solar-wind outflow, and between 22% and 36% escapes as z {sup +} fluctuations through the simulation boundary at r = r{sub A}. The z {sup } power spectra scale like k{sub perpendicular}{sup -?{sup }}, where k is the wavenumber in the plane perpendicular to B{sub 0}. In our simulation with the smallest value of ?{sub c{sub sun}{sup +}} (?2 minutes) and largest value of L{sub ?} (2 10{sup 4} km), we find that ?{sup +} decreases approximately linearly with increasing ln (r), reaching a value of 1.3 at r = 11.1 R{sub ?}. Our simulations with larger values of ?{sub c{sub sun}{sup +}} exhibit alignment between the contours of constant ?{sup +}, ?{sup }, ?{sub 0}{sup +}, and ?{sub 0}{sup -}, where ?{sup } are the Elssser potentials and ?{sub 0}{sup } are the outer-scale parallel Elssser vorticities.
Numerical simulation of a thermoacoustic refrigerator. 2: Stratified flow around the stack
Worlikar, A.S.; Knio, O.M.; Klein, R.
1998-08-10
The unsteady, two-dimensional, thermally stratified flow in the neighborhood of an idealized thermoacoustic stack is analyzed using a low-Mach-number model that extends the adiabatic flow scheme developed in part 1 (Journal of Computational Physics 127, 424 (1996)). The extension consists of incorporation of numerical solvers for the energy equations in the fluid and the stack plates, and construction and implementation of fast Poisson solver for the velocity potential based on a domain decomposition/boundary Green`s function technique. The unsteady computations are used to predict the steady-state, acoustically generated temperature gradient across a two-dimensional couple and to analyze its dependence on the amplitude of the prevailing resonant wave. Computed results are compared to theoretical predictions and experimental data.
Numerical simulation study on fluid dynamics of plasma window using argon
Huang, S.; Zhu, K.; Shi, B. L.; Lu, Y. R.; Hershcovitch, A.; Yang, L.; Zhang, X. Y.; Wei, G. D.
2013-07-15
In this paper, a numerical 2D FLUENT-based magneto-hydrodynamic model has been developed to investigate the arc and flow field of plasma window, which is used as a windowless vacuum sealing device. The gas inlet, arc creation-developing and plasma expansion segments are all incorporated together in the integral model. An axis-symmetry cathode structure (hollow cathode) is used in the model. Current distribution of the arc is presented and discussed. The temperature, velocity, and pressure field are presented to show the physical mechanisms for the high pressure gap within the plasma window. Flow acceleration and viscosity effect are concluded as the main reasons for the pressure drop. The result for the pressure distribution in the cylindrical tube section has a good agreement with the analytical model. The validation for the sealing ability of plasma window is verified.
Real time control and numerical simulation of pipeline subjected to landslide
Cuscuna, S.; Giusti, G.; Gramola, C.
1984-06-01
This paper describes SNAM research activity in the study of behaviour and real-time control of pipelines in landslide areas. The subject can be delt considering three different aspects: 1. Geotechnical characterization of unstable soils. The mechanical parameters of soil and the landslide types are defined; 2. Structural analysis of pipe-soil system. By means of a finite element program it's possible to study the pipe-soil interaction; in this numerical code the soil parameters attend by the non-linear elastic behaviour of pipe restraints. The results of this analysis are the location of the expected most stressed sections of pipe and the global behaviour of pipe inside the soil. 3. Instrumental control. The adoption of a suitable appliance of vibrating wire strain gauges allows the strain control of pipe in time. The aim is to make possible timely interventions in order to guarantee the installation safety.
Three dimensional numerical simulations of the UPS-292 stratified charge engine
O'Rourke, P.J.; Amsden, A.A.
1987-01-01
The authors present and analyze three-dimensional calculations of the spray, mixing and combustion in the UPS-292 stratified charge engine for three different operating conditions, corresponding to overall air-fuel ratios between 22.4 and 61.0. The numerical calculations are performed with KIVA, a multidimensional arbitrary-mesh, finite-difference hydrodynamics program for internal combustion engine applications. The calculations use a mesh of 10,000 computational cells. Each operating condition is calculated from intake valve closure at 118/sup 0/ BTDC to 90/sup 0/ ATDC and requires approximately three hours of CRAY-XMP computer time. Combustion occurs primarily in the wake of the spark plug, and to include the effects of the spark plug on the flow field, we use a novel internal obstacle treatment. The methodology, in which internal obstacles are represented by computational particles, promises to be applicable to the calculation of the flows around intake and exhaust valves.
Hakan Ozaltun & Herman Shen
2011-11-01
This article presents assessment of the mechanical behavior of U-10wt% Mo (U10Mo) alloy based monolithic fuel plates subject to irradiation. Monolithic, plate-type fuel is a new fuel form being developed for research and test reactors to achieve higher uranium densities within the reactor core to allow the use of low-enriched uranium fuel in high-performance reactors. Identification of the stress/strain characteristics is important for understanding the in-reactor performance of these plate-type fuels. For this work, three distinct cases were considered: (1) fabrication induced residual stresses (2) thermal cycling of fabricated plates; and finally (3) transient mechanical behavior under actual operating conditions. Because the temperatures approach the melting temperature of the cladding during the fabrication and thermal cycling, high temperature material properties were incorporated to improve the accuracy. Once residual stress fields due to fabrication process were identified, solution was used as initial state for the subsequent simulations. For thermal cycling simulation, elasto-plastic material model with thermal creep was constructed and residual stresses caused by the fabrication process were included. For in-service simulation, coupled fluid-thermal-structural interaction was considered. First, temperature field on the plates was calculated and this field was used to compute the thermal stresses. For time dependent mechanical behavior, thermal creep of cladding, volumetric swelling and fission induced creep of the fuel foil were considered. The analysis showed that the stresses evolve very rapidly in the reactor. While swelling of the foil increases the stress of the foil, irradiation induced creep causes stress relaxation.
Numerical simulation and design of a fluxset sensor by finite element method
Preis, K.; Bardi, I.; Biro, O.; Richter, K.R.; Pavo, J.; Gasparics, A.; Ticar, I.
1998-09-01
A 3D model of a fluxset sensor serving to measure magnetic fields arising in Eddy Current Nondestructive Testing applications is analyzed by the finite element method. The voltage induced in the pick-up coil is obtained by computing the flux of the core of the sensor for several values of the exciting current at various external fields. It is shown that the time shift of the ensuing voltage impulse depends linearly on the external field in a wide range. The behavior of the sensor is furthermore simulated in a real nondestructive testing arrangement consisting of an exciting coil located above a conducting plate with a crack.
Numerical Simulation of Earth Pressure on Head Chamber of Shield Machine with FEM
Li Shouju; Kang Chengang [State Key Laboratory of structural analysis for industrial equipment, Dalian University of Technology, Dalian 116023 (China); Sun, Wei [School of Mechanical Engineering, Dalian University of Technology, Dalian 116023 (China); Shangguan Zichang [School of Civil and Hydraulic Engineering, Dalian University of Technology, Dalian 116023 (China); Institute of Civil Engineering, Dalian Fishery University, Dalian 116023 (China)
2010-05-21
Model parameters of conditioned soils in head chamber of shield machine are determined based on tree-axial compression tests in laboratory. The loads acting on tunneling face are estimated according to static earth pressure principle. Based on Duncan-Chang nonlinear elastic constitutive model, the earth pressures on head chamber of shield machine are simulated in different aperture ratio cases for rotating cutterhead of shield machine. Relationship between pressure transportation factor and aperture ratio of shield machine is proposed by using aggression analysis.
Numerical simulation of cathode plasma dynamics in magnetically insulated vacuum transmission lines
Thoma, C.; Genoni, T. C.; Welch, D. R.; Rose, D. V.; Clark, R. E.; Miller, C. L.; Stygar, W. A.; Kiefer, M. L.
2015-03-15
A novel algorithm for the simulation of cathode plasmas in particle-in-cell codes is described and applied to investigate cathode plasma evolution in magnetically insulated transmission lines (MITLs). The MITL electron sheath is modeled by a fully kinetic electron species. Electron and ion macroparticles, both modeled as fluid species, form a dense plasma which is initially localized at the cathode surface. Energetic plasma electron particles can be converted to kinetic electrons to resupply the electron flux at the plasma edge (the “effective” cathode). Using this model, we compare results for the time evolution of the cathode plasma and MITL electron flow with a simplified (isothermal) diffusion model. Simulations in 1D show a slow diffusive expansion of the plasma from the cathode surface. But in multiple dimensions, the plasma can expand much more rapidly due to anomalous diffusion caused by an instability due to the strong coupling of a transverse magnetic mode in the electron sheath with the expanding resistive plasma layer.
Bocchi, M.; Ummels, B.; Chittenden, J. P.; Lebedev, S. V.; Frank, A.; Blackman, E. G.
2013-04-10
The physics of accretion disks is of fundamental importance for understanding of a wide variety of astrophysical sources that includes protostars, X-ray binaries, and active galactic nuclei. The interplay between hydrodynamic flows and magnetic fields and the potential for turbulence-producing instabilities is a topic of active research that would benefit from the support of dedicated experimental studies. Such efforts are in their infancy, but in an effort to push the enterprise forward we propose an experimental configuration which employs a modified cylindrical wire array Z-pinch to produce a rotating plasma flow relevant to accretion disks. We present three-dimensional resistive magnetohydrodynamic simulations which show how this approach can be implemented. In the simulations, a rotating plasma cylinder or ring is formed, with typical rotation velocity {approx}30 km s{sup -1}, Mach number {approx}4, and Reynolds number in excess of 10{sup 7}. The plasma is also differentially rotating. Implementation of different external magnetic field configurations is discussed. It is found that a modest uniform vertical field of 1 T can affect the dynamics of the system and could be used to study magnetic field entrainment and amplification through differential rotation. A dipolar field potentially relevant to the study of accretion columns is also considered.
Numerical simulations of optically thick accretion onto a black hole. II. Rotating flow
Fragile, P. Chris; Olejar, Ally; Anninos, Peter
2014-11-20
In this paper, we report on recent upgrades to our general relativistic radiation magnetohydrodynamics code, Cosmos++, including the development of a new primitive inversion scheme and a hybrid implicit-explicit solver with a more general M {sub 1} closure relation for the radiation equations. The new hybrid solver helps stabilize the treatment of the radiation source terms, while the new closure allows for a much broader range of optical depths to be considered. These changes allow us to expand by orders of magnitude the range of temperatures, opacities, and mass accretion rates, and move a step closer toward our goal of performing global simulations of radiation-pressure-dominated black hole accretion disks. In this work, we test and validate the new method against an array of problems. We also demonstrate its ability to handle super-Eddington, quasi-spherical accretion. Even with just a single proof-of-principle simulation, we already see tantalizing hints of the interesting phenomenology associated with the coupling of radiation and gas in super-Eddington accretion flows.
Marxen, Olaf, E-mail: olaf.marxen@vki.ac.be [Center for Turbulence Research, Building 500, Stanford University, Stanford, CA 94305-3035 (United States) [Center for Turbulence Research, Building 500, Stanford University, Stanford, CA 94305-3035 (United States); Aeronautics and Aerospace Department, von Karman Institute for Fluid Dynamics, Chausse de Waterloo, 72, 1640 Rhode-St-Gense (Belgium); Magin, Thierry E. [Aeronautics and Aerospace Department, von Karman Institute for Fluid Dynamics, Chausse de Waterloo, 72, 1640 Rhode-St-Gense (Belgium)] [Aeronautics and Aerospace Department, von Karman Institute for Fluid Dynamics, Chausse de Waterloo, 72, 1640 Rhode-St-Gense (Belgium); Shaqfeh, Eric S.G.; Iaccarino, Gianluca [Center for Turbulence Research, Building 500, Stanford University, Stanford, CA 94305-3035 (United States)] [Center for Turbulence Research, Building 500, Stanford University, Stanford, CA 94305-3035 (United States)
2013-12-15
A new numerical method is presented here that allows to consider chemically reacting gases during the direct numerical simulation of a hypersonic fluid flow. The method comprises the direct coupling of a solver for the fluid mechanical model and a library providing the physio-chemical model. The numerical method for the fluid mechanical model integrates the compressible NavierStokes equations using an explicit time advancement scheme and high-order finite differences. This NavierStokes code can be applied to the investigation of laminar-turbulent transition and boundary-layer instability. The numerical method for the physio-chemical model provides thermodynamic and transport properties for different gases as well as chemical production rates, while here we exclusively consider a five species air mixture. The new method is verified for a number of test cases at Mach 10, including the one-dimensional high-temperature flow downstream of a normal shock, a hypersonic chemical reacting boundary layer in local thermodynamic equilibrium and a hypersonic reacting boundary layer with finite-rate chemistry. We are able to confirm that the diffusion flux plays an important role for a high-temperature boundary layer in local thermodynamic equilibrium. Moreover, we demonstrate that the flow for a case previously considered as a benchmark for the investigation of non-equilibrium chemistry can be regarded as frozen. Finally, the new method is applied to investigate the effect of finite-rate chemistry on boundary layer instability by considering the downstream evolution of a small-amplitude wave and comparing results with those obtained for a frozen gas as well as a gas in local thermodynamic equilibrium.
PROBABILISTIC SIMULATION OF SUBSURFACE FLUID FLOW: A STUDY USING A NUMERICAL SCHEME
Buscheck, Timothy Eric
1980-03-01
There has been an increasing interest in probabilistic modeling of hydrogeologic systems. The classical approach to groundwater modeling has been deterministic in nature, where individual layers and formations are assumed to be uniformly homogeneous. Even in the case of complex heterogeneous systems, the heterogeneities describe the differences in parameter values between various layers, but not within any individual layer. In a deterministic model a single-number is assigned to each hydrogeologic parameter, given a particular scale of interest. However, physically there is no such entity as a truly uniform and homogeneous unit. Single-number representations or deterministic predictions are subject to uncertainties. The approach used in this work models such uncertainties with probabilistic parameters. The resulting statistical distributions of output variables are analyzed. A numerical algorithm, based on axiomatic principles of probability theory, performs arithmetic operations between probability distributions. Two subroutines are developed from the algorithm and incorporated into the computer program TERZAGI, which solves groundwater flow problems in saturated, multi-dimensional systems. The probabilistic computer program is given the name, PROGRES. The algorithm has been applied to study the following problems: one-dimensional flow through homogeneous media, steady-state and transient flow conditions, one-dimensional flow through heterogeneous media, steady-state and transient flow conditions, and two-dimensional steady-stte flow through heterogeneous media. The results are compared with those available in the literature.
Crandall, Dustin; Bromhal, Grant; Karpyn, Zuleima T.
2010-07-01
Understanding how fracture wall-roughness affects fluid flow is important when modeling many subsurface transport problems. Computed tomography scanning provides a unique view of rock fractures, allowing the measurement of fracture wall-roughness, without destroying the initial rock sample. For this computational fluid dynamics study, we used several different methods to obtain three-dimensional meshes of a computed tomography scanned fracture in Berea sandstone. These volumetric meshes had different wall-roughnesses, which we characterized using the Joint Roughness Coefficient and the fractal dimension of the fracture profiles. We then related these macroscopic roughness parameters to the effective flow through the fractures, as determined from Navier-Stokes numerical models. Thus, we used our fracture meshes to develop relationships between the observed roughness properties of the fracture geometries and flow parameters that are of importance for modeling flow through fractures in field scale models. Fractures with high Joint Roughness Coefficients and fractal dimensions were shown to exhibit tortuous flow paths, be poorly characterized by the mean geometric aperture, and have a fracture transmissivity 35 times smaller than the smoother modeled fracture flows.
Two-dimensional numerical simulation of a Stirling engine heat exchanger
Ibrahim, M.B.; Tew, R.C.; Dudenhoefer, J.E.
1994-09-01
This paper describes the first phase of an effort to develop multidimensional models of Stirling engine components; the ultimate goal is to model an entire engine working space. More specifically, this paper describes parallel plate and tubular heat exchanger models with emphasis on the central part of the channel (i.e., ignoring hydrodynamic and thermal end effects). The model assumes: Laminar, incompressible flow with constant thermophysical properties. In addition, a constant axial temperature gradient is imposed. The governing equations, describing the model, have been solved Crack-Nicloson finite-difference scheme. Model predictions have been compared with analytical solutions for oscillating/reversing flow and heat transfer in order to check numerical accuracy. The simplifying assumptions will later be relaxed to permit modeling of incompressible, laminar/turbulent flow that occurs in Stirling heat exchanger. Excellent agreement has been obtained for the model predictions with analytical solutions available for both flow in circular tubes and between parallel plates. Also the heat transfer computational results are in good agreement with the heat transfer analytical results for parallel plates.
New Particle-in-Cell Code for Numerical Simulation of Coherent Synchrotron Radiation
Balsa Terzic, Rui Li
2010-05-01
We present a first look at the new code for self-consistent, 2D simulations of beam dynamics affected by the coherent synchrotron radiation. The code is of the particle-in-cell variety: the beam bunch is sampled by point-charge particles, which are deposited on the grid; the corresponding forces on the grid are then computed using retarded potentials according to causality, and interpolated so as to advance the particles in time. The retarded potentials are evaluated by integrating over the 2D path history of the bunch, with the charge and current density at the retarded time obtained from interpolation of the particle distributions recorded at discrete timesteps. The code is benchmarked against analytical results obtained for a rigid-line bunch. We also outline the features and applications which are currently being developed.
Stevens, B.; Cotton, W.R.; Feingold, G.
1996-04-01
Over the past decade or so the evolution and equilibria of persistent decks of stratocumulus climatologically clinging to the edge of summertime subtropical highs have been an issue of increased scientific inquiry. The particular interest in the microphysical structure of these clouds stems from a variety of hypotheses which suggest that anthropogenic influences or biogenic feedbacks may alter the structure of these clouds in a climatically significant manner. Most of these hypotheses are quite tentative, based as they are on simple formulations of boundary layer structures and interactions between drops and aerosols. This work is concerned with an assessment of the microphysical structure of marine stratocumulus as simulated by an LES-EM model.
Numerical simulation of fiber and wire array Z-pinches with Trac-II
Reisman, D
1998-09-01
Trac-II is a two dimensional axisymmetric resistive MHD code. It simulates all three spatial components (r, z, φ) of the magnetic field and fluid velocity vectors, and the plasma is treated as a single fluid with two temperatures (T_{e},T_{i}). In addition, it can optionally include a self-consistent external circuit. Recent modifications to the code include the addition of the 3-T radiation model, a 4-phase (solid-liquid-vapor-plasma) equation of state model (QEOS), a 4-phase electrical/thermal conductivity model, and an implicit solution of poloidal B_{z},B_{r}) magnetic field diffusion. These changes permit a detailed study of fiber and wire array Z-pinches. Specifically, Trac-II is used to study the wire array Z-pinch at the PBFA-Z pulse power generator at Sandia National Laboratory. First, in 1-D we examine the behavior of a single wire in the Z-pinch. Then, using these results as initial radial conditions in 2-D, we investigate the dynamics of wire array configurations in the r-z and r-θ plane. In the r-z plane we examine the growth of the m=0 or "sausage" instability in single wires within the array. In the r-θ plane we examine the merging behavior between neighboring wires. Special emphasis is placed on trying to explain how instability growth affects the performance of the Z-pinch. Lastly, we introduce Trac-III, a 3-D MHD code, and illustrate the m=1 or "kink" instability. We also discuss how Trac-III can be modified to simulate the wire array Z-pinch.
Numerical simulation of a wave-guide mixing layer on a Cray C-90
Greenough, J.A.; Crutchfield, W.Y.; Rendleman, C.A.
1995-05-19
The development of a three-dimensional spatially evolving compressible mixing layer is investigated numerically using a parallel implementation of Adaptive Mesh Refinement (AMR) on a Cray C-90. The parallel implementation allowed the flow to be highly resolved while significantly reducing the wall-clock runtime. A sustained computation rate of 5.3 Gigaflops including I/O was obtained for a typical production run on a 16 processor machine. A novel mixing layer configuration is investigated where a pressure mismatch is maintained between the two inlet streams. In addition, the sonic character of the two streams is sufficiently different so that the pressure relief wave is trapped in the high speed stream. The trapped wave forces the mixing layer to form a characteristic cellular pattern. The cellular structure introduces curvature into the mixing layer that excites centrifugal instabilities characterized by large-scale counter-rotating vortical pairs embedded within the mixing layer. These are the dominant feature of the flow. Visualizations of these structures in cross-section show the pumping action which lifts dense fluid up into light gas. This effect has a strong impact on mixing enhancement as monitored by a conserved scalar formulation. Once the large-scale structures axe well established in the flow and undergo intensification from favorable velocity gradients, the time-averaged integrated product shows almost a four-fold increase. A spectral analysis of the flow-field over the cellular structures, as part of a full space-time analysis, shows these structures to be zero-frequency modes that develop from low level essentially broad-banded noise. This characterization of the vortical structures and their appearance is consistent with a recent linear stability analysis, of a mixing layer over a curved wall that predicts the most unstable modes to be zero frequency streamwise vortices.
Watts, C.A.
1993-09-01
In this dissertation the possibility that chaos and simple determinism are governing the dynamics of reversed field pinch (RFP) plasmas is investigated. To properly assess this possibility, data from both numerical simulations and experiment are analyzed. A large repertoire of nonlinear analysis techniques is used to identify low dimensional chaos in the data. These tools include phase portraits and Poincare sections, correlation dimension, the spectrum of Lyapunov exponents and short term predictability. In addition, nonlinear noise reduction techniques are applied to the experimental data in an attempt to extract any underlying deterministic dynamics. Two model systems are used to simulate the plasma dynamics. These are the DEBS code, which models global RFP dynamics, and the dissipative trapped electron mode (DTEM) model, which models drift wave turbulence. Data from both simulations show strong indications of low dimensional chaos and simple determinism. Experimental date were obtained from the Madison Symmetric Torus RFP and consist of a wide array of both global and local diagnostic signals. None of the signals shows any indication of low dimensional chaos or low simple determinism. Moreover, most of the analysis tools indicate the experimental system is very high dimensional with properties similar to noise. Nonlinear noise reduction is unsuccessful at extracting an underlying deterministic system.
Numerical simulation of a thermoacoustic refrigerator. I. Unsteady adiabatic flow around the stack
Worlikar, A.S.; Knio, O.M.
1996-09-01
A low Mach-number compressible flow model for the simulation of acoustically driven flow in a thermoacoustic stack is constructed. The model is based on the assumption that the acoustic wavelength is much larger than the characteristic hydrodynamic lengthscale. Thus, a simplified description of the flow is obtained which still retains the essential features of acoustically induced velocity oscillations near solid boundaries. A vorticity-based formulation of the governing equation is derived which relies on the Helmholtz decomposition of the velocity vector into irrotational and divergence-free components. Irrotational motion is used to represent the action of acoustic waves. Meanwhile the divergence-free velocity component is used to capture the nonlinear vortical perturbations due to no-slip boundaries. A simplified version of the model is applied to analyze unsteady flow in the neighborhood of an idealized thermo-acoustic stack which consists of a periodic array of thin plates placed in an acoustic standing wave. Computed results are used to analyze, for different stack configurations, the nonlinear response of the flow to different acoustic driving amplitudes and frequencies. In particular, it is shown that the flow is dominated by the motion of vortices which result from the shedding of boundary layers from the edges of the stack. The dependence of energy losses on stack configuration and operating conditions is also examined. 28 refs., 23 figs., 2 tabs.
A STELLAR WIND ORIGIN FOR THE G2 CLOUD: THREE-DIMENSIONAL NUMERICAL SIMULATIONS
De Colle, Fabio; Raga, A. C.; Contreras-Torres, Flavio F.; Toledo-Roy, Juan C.
2014-07-10
We present three-dimensional, adaptive mesh refinement simulations of G2, a cloud of gas moving in a highly eccentric orbit toward the galactic center. We assume that G2 originates from a stellar wind interacting with the environment of the Sgr A* black hole. The stellar wind forms a cometary bubble which becomes increasingly elongated as the star approaches periastron. A few months after periastron passage, streams of material begin to accrete on the central black hole with accretion rates M-dot ?10{sup ?8} M {sub ?}yr{sup 1}. Predicted Br? emission maps and position-velocity diagrams show an elongated emission resembling recent observations of G2. A large increase in luminosity is predicted by the emission coming from the shocked wind region during periastron passage. The observations, showing a constant Br? luminosity, remain puzzling, and are explained here assuming that the emission is dominated by the free-wind region. The observed Br? luminosity (?8 10{sup 30}ergs{sup 1}) is reproduced by a model with a v{sub w} = 50kms{sup 1} wind velocity and a 10{sup 7} M {sub ?}yr{sup 1} mass-loss rate if the emission comes from the shocked wind. A faster and less dense wind reproduces the Br? luminosity if the emission comes from the inner, free-wind region. The extended cometary wind bubble, largely destroyed by the tidal interaction with the black hole, reforms a few years after periastron passage. As a result, the Br? emission is more compact after periastronpassage.
Masada, Youhei; Sano, Takayoshi E-mail: sano@ile.osaka-u.ac.jp
2014-10-10
The mechanism of large-scale dynamos in rigidly rotating stratified convection is explored by direct numerical simulations (DNS) in Cartesian geometry. A mean-field dynamo model is also constructed using turbulent velocity profiles consistently extracted from the corresponding DNS results. By quantitative comparison between the DNS and our mean-field model, it is demonstrated that the oscillatory α{sup 2} dynamo wave, excited and sustained in the convection zone, is responsible for large-scale magnetic activities such as cyclic polarity reversal and spatiotemporal migration. The results provide strong evidence that a nonuniformity of the α-effect, which is a natural outcome of rotating stratified convection, can be an important prerequisite for large-scale stellar dynamos, even without the Ω-effect.
Neeraj Gupta
2008-03-31
A series of numerical simulations of carbon dioxide (CO{sub 2}) injection were conducted as part of a program to assess the potential for geologic sequestration in deep geologic reservoirs (the Rose Run and Copper Ridge formations), at the American Electric Power (AEP) Mountaineer Power Plant outside of New Haven, West Virginia. The simulations were executed using the H{sub 2}O-CO{sub 2}-NaCl operational mode of the Subsurface Transport Over Multiple Phases (STOMP) simulator (White and Oostrom, 2006). The objective of the Rose Run formation modeling was to predict CO{sub 2} injection rates using data from the core analysis conducted on the samples. A systematic screening procedure was applied to the Ohio River Valley CO{sub 2} storage site utilizing the Features, Elements, and Processes (FEP) database for geological storage of CO{sub 2} (Savage et al., 2004). The objective of the screening was to identify potential risk categories for the long-term geological storage of CO{sub 2} at the Mountaineer Power Plant in New Haven, West Virginia. Over 130 FEPs in seven main classes were assessed for the project based on site characterization information gathered in a geological background study, testing in a deep well drilled on the site, and general site conditions. In evaluating the database, it was apparent that many of the items were not applicable to the Mountaineer site based its geologic framework and environmental setting. Nine FEPs were identified for further consideration for the site. These FEPs generally fell into categories related to variations in subsurface geology, well completion materials, and the behavior of CO{sub 2} in the subsurface. Results from the screening were used to provide guidance on injection system design, developing a monitoring program, performing reservoir simulations, and other risk assessment efforts. Initial work indicates that the significant FEPs may be accounted for by focusing the storage program on these potential issues. The
Wang, X.H.; Zhao, D.Q.; Jiang, L.Q.; Yang, W.B.
2009-07-01
Numerical analysis was used to study the deposition and burning characteristics of combining co-combustion with slagging combustion technologies in this paper. The pyrolysis and burning kinetic models of different fuels were implanted into the WBSF-PCC2 (wall burning and slag flow in pulverized co-combustion) computation code, and then the slagging and co-combustion characteristics (especially the wall burning mechanism of different solid fuels and their effects on the whole burning behavior in the cylindrical combustor at different mixing ratios under the condition of keeping the heat input same) were simulated numerically. The results showed that adding wood powder at 25% mass fraction can increase the temperature at the initial stage of combustion, which is helpful to utilize the front space of the combustor. Adding wood powder at a 25% mass fraction can increase the reaction rate at the initial combustion stage; also, the coal ignitability is improved, and the burnout efficiency is enhanced by about 5% of suspension and deposition particles, which is helpful for coal particles to burn entirely and for combustion devices to minimize their dimensions or sizes. The results also showed that adding wood powder at a proper ratio is helpful to keep the combustion stability, not only because of the enhancement for the burning characteristics, but also because the running slag layer structure can be changed more continuously, which is very important for avoiding the abnormal slag accumulation in the slagging combustor. The theoretic analysis in this paper proves that unification of co-combustion and slagging combustion technologies is feasible, though more comprehensive and rigorous research is needed.
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Infrastructure Energy Storage Nuclear Power & Engineering Grid Modernization Battery Testing ... Heavy Duty Fuels DISI Combustion HCCISCCI Fundamentals Spray Combustion Modeling ...
Kwon, Kyung; Fan, Liang-Shih; Zhou, Qiang; Yang, Hui
2014-09-30
A new and efficient direct numerical method with second-order convergence accuracy was developed for fully resolved simulations of incompressible viscous flows laden with rigid particles. The method combines the state-of-the-art immersed boundary method (IBM), the multi-direct forcing method, and the lattice Boltzmann method (LBM). First, the multi-direct forcing method is adopted in the improved IBM to better approximate the no-slip/no-penetration (ns/np) condition on the surface of particles. Second, a slight retraction of the Lagrangian grid from the surface towards the interior of particles with a fraction of the Eulerian grid spacing helps increase the convergence accuracy of the method. An over-relaxation technique in the procedure of multi-direct forcing method and the classical fourth order Runge-Kutta scheme in the coupled fluid-particle interaction were applied. The use of the classical fourth order Runge-Kutta scheme helps the overall IB-LBM achieve the second order accuracy and provides more accurate predictions of the translational and rotational motion of particles. The preexistent code with the first-order convergence rate is updated so that the updated new code can resolve the translational and rotational motion of particles with the second-order convergence rate. The updated code has been validated with several benchmark applications. The efficiency of IBM and thus the efficiency of IB-LBM were improved by reducing the number of the Lagragian markers on particles by using a new formula for the number of Lagrangian markers on particle surfaces. The immersed boundary-lattice Boltzmann method (IBLBM) has been shown to predict correctly the angular velocity of a particle. Prior to examining drag force exerted on a cluster of particles, the updated IB-LBM code along with the new formula for the number of Lagrangian markers has been further validated by solving several theoretical problems. Moreover, the unsteadiness of the drag force is examined when a
Danish, Mohammad Suman, Sawan Srinivasan, Balaji
2014-12-15
The pressure Hessian tensor plays a key role in shaping the behavior of the velocity gradient tensor, and in turn, that of many incumbent non-linear processes in a turbulent flow field. In compressible flows, the role of pressure Hessian is even more important because it represents the level of fluid-thermodynamic coupling existing in the flow field. In this work, we first perform a direct numerical simulation-based study to clearly identify, isolate, and understand various important inviscid mechanisms that govern the evolution of the pressure Hessian tensor in compressible turbulence. The ensuing understanding is then employed to introduce major improvements to the existing Lagrangian model of the pressure Hessian tensor (the enhanced Homogenized Euler equation or EHEE) in terms of (i) non-symmetric, non-isentropic effects and (ii) improved representation of the anisotropic portion of the pressure Hessian tensor. Finally, we evaluate the new model extensively by comparing the new model results against known turbulence behavior over a range of Reynolds and Mach numbers. Indeed, the new model shows much improved performance as compared to the EHEE model.
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Schilling, Oleg; Mueschke, Nicholas J.
2010-10-18
Data from a 1152X760X1280 direct numerical simulation (DNS) of a transitional Rayleigh-Taylor mixing layer modeled after a small Atwood number water channel experiment is used to comprehensively investigate the structure of mean and turbulent transport and mixing. The simulation had physical parameters and initial conditions approximating those in the experiment. The budgets of the mean vertical momentum, heavy-fluid mass fraction, turbulent kinetic energy, turbulent kinetic energy dissipation rate, heavy-fluid mass fraction variance, and heavy-fluid mass fraction variance dissipation rate equations are constructed using Reynolds averaging applied to the DNS data. The relative importance of mean and turbulent production, turbulent dissipationmoreand destruction, and turbulent transport are investigated as a function of Reynolds number and across the mixing layer to provide insight into the flow dynamics not presently available from experiments. The analysis of the budgets supports the assumption for small Atwood number, Rayleigh/Taylor driven flows that the principal transport mechanisms are buoyancy production, turbulent production, turbulent dissipation, and turbulent diffusion (shear and mean field production are negligible). As the Reynolds number increases, the turbulent production in the turbulent kinetic energy dissipation rate equation becomes the dominant production term, while the buoyancy production plateaus. Distinctions between momentum and scalar transport are also noted, where the turbulent kinetic energy and its dissipation rate both grow in time and are peaked near the center plane of the mixing layer, while the heavy-fluid mass fraction variance and its dissipation rate initially grow and then begin to decrease as mixing progresses and reduces density fluctuations. All terms in the transport equations generally grow or decay, with no qualitative change in their profile, except for the pressure flux contribution to the total turbulent kinetic
Diwekar, Urmila; Shastri, Yogendra (Vishwamitra Research Institute Clarendon Hills, IL); Subrmanyan, Karthik; Zitney, S.E.
2007-11-04
APECS (Advanced Process Engineering Co-Simulator) is an integrated software suite that combines the power of process simulation with high-fidelity, computational fluid dynamics (CFD) for improved design, analysis, and optimization of process engineering systems. The APECS system uses commercial process simulation (e.g., Aspen Plus) and CFD (e.g., FLUENT) software integrated with the process-industry standard CAPE-OPEN (CO) interfaces. This breakthrough capability allows engineers to better understand and optimize the fluid mechanics that drive overall power plant performance and efficiency. The focus of this paper is the CAPE-OPEN complaint stochastic modeling and reduced order model computational capability around the APECS system. The usefulness of capabilities is illustrated with coal fired, gasification based, FutureGen power plant simulation. These capabilities are used to generate efficient reduced order models and optimizing model complexities.
Rosenberg, Duane L; Pouquet, Dr. Annick; Mininni, Dr. Pablo D.; Marino, Dr. Raffaele
2015-01-01
We report results on rotating stratified turbulence in the absence of forcing, with large-scale isotropic initial conditions, using direct numerical simulations computed on grids of up to $4096^3$ points. The Reynolds and Froude numbers are respectively equal to $Re=5.4\\times 10^4$ and $Fr=0.0242$. The ratio of the Brunt-V\\"ais\\"al\\"a to the inertial wave frequency, $N/f$, is taken to be equal to 5, a choice appropriate to model the dynamics of the southern abyssal ocean at mid latitudes. This gives a global buoyancy Reynolds number $R_B=ReFr^2=32$, a value sufficient for some isotropy to be recovered in the small scales beyond the Ozmidov scale, but still moderate enough that the intermediate scales where waves are prevalent are well resolved. We concentrate on the large-scale dynamics and confirm that the Froude number based on a typical vertical length scale is of order unity, with strong gradients in the vertical. Two characteristic scales emerge from this computation, and are identified from sharp variations in the spectral distribution of either total energy or helicity. A spectral break is also observed at a scale at which the partition of energy between the kinetic and potential modes changes abruptly, and beyond which a Kolmogorov-like spectrum recovers. Large slanted layers are ubiquitous in the flow in the velocity and temperature fields, and a large-scale enhancement of energy is also observed, directly attributable to the effect of rotation.
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Schilling, Oleg; Mueschke, Nicholas J.
2010-10-18
Data from a 1152X760X1280 direct numerical simulation (DNS) of a transitional Rayleigh-Taylor mixing layer modeled after a small Atwood number water channel experiment is used to comprehensively investigate the structure of mean and turbulent transport and mixing. The simulation had physical parameters and initial conditions approximating those in the experiment. The budgets of the mean vertical momentum, heavy-fluid mass fraction, turbulent kinetic energy, turbulent kinetic energy dissipation rate, heavy-fluid mass fraction variance, and heavy-fluid mass fraction variance dissipation rate equations are constructed using Reynolds averaging applied to the DNS data. The relative importance of mean and turbulent production, turbulent dissipationmore » and destruction, and turbulent transport are investigated as a function of Reynolds number and across the mixing layer to provide insight into the flow dynamics not presently available from experiments. The analysis of the budgets supports the assumption for small Atwood number, Rayleigh/Taylor driven flows that the principal transport mechanisms are buoyancy production, turbulent production, turbulent dissipation, and turbulent diffusion (shear and mean field production are negligible). As the Reynolds number increases, the turbulent production in the turbulent kinetic energy dissipation rate equation becomes the dominant production term, while the buoyancy production plateaus. Distinctions between momentum and scalar transport are also noted, where the turbulent kinetic energy and its dissipation rate both grow in time and are peaked near the center plane of the mixing layer, while the heavy-fluid mass fraction variance and its dissipation rate initially grow and then begin to decrease as mixing progresses and reduces density fluctuations. All terms in the transport equations generally grow or decay, with no qualitative change in their profile, except for the pressure flux contribution to the total turbulent kinetic
Chen, Jacqueline H.; Hawkes, Evatt R.
2004-08-01
Direct numerical simulation (DNS) with complex chemistry was used to study statistics of displacement and consumption speeds in turbulent lean premixed methane-air flames. The main focus of the study is an evaluation of the extent to which a turbulent flame in the thin reaction zones regime can be described by an ensemble of strained laminar flames. Conditional averages with respect to strain for displacement and consumption speeds are presented over a wide range of strain typically encountered in a turbulent flame, compared with previous studies that either made local pointwise comparisons or conditioned the data on small strain and curvature. The conditional averages for positive strains are compared with calculated data from two different canonical strained laminar configurations to determine which is the optimal representation of a laminar flame structure embedded in a turbulent flame: the reactant-to-product (R-to-P) configuration or the symmetric twin flame configuration. Displacement speed statistics are compared for the progress-variable isosurface of maximum reaction rate and an isosurface toward the fresh gases, which are relevant for both modeling and interpretation of experiment results. Displacement speeds in the inner reaction layer are found to agree very well with the laminar R-to-P calculations over a wide range of strain for higher Damkhler number conditions, well beyond the regime in which agreement was expected. For lower Damkhler numbers, a reduced response to strain is observed, consistent with previous studies and theoretical expectations. Compared with the inner layer, broader and shifted probability density functions (PDFs) of displacement speed were observed in the fresh gases, and the agreement with the R-to-P calculations deteriorated. Consumption speeds show a poorer agreement with strained laminar calculations, which is attributed to multidimensional effects and a more attenuated unsteady response to strain fluctuations; however, they
Scheibe, Timothy D.; Richmond, Marshall C.
2002-01-30
This paper describes a numerical model of juvenile salmonid migration in the Columbia and Snake Rivers. The model, called the Fish Individual-based Numerical Simulator or FINS, employs a discrete, particle-based approach to simulate the migration and history of exposure to dissolved gases of individual fish. FINS is linked to a two-dimensional (vertically-averaged) hydrodynamic simulator that quantifies local water velocity, temperature, and dissolved gas levels as a function of river flow rates and dam operations. Simulated gas exposure histories can be input to biological mortality models to predict the effects of various river configurations on fish injury and mortality due to dissolved gas supersaturation. Therefore, FINS serves as a critical linkage between hydrodynamic models of the river system and models of biological impacts. FINS was parameterized and validated based on observations of individual fish movements collected using radiotelemetry methods during 1997 and 1998. A quasi-inverse approach was used to decouple fish swimming movements from advection with the local water velocity, allowing inference of time series of non-advective displacements of individual fish from the radiotelemetry data. Statistical analyses of these displacements are presented, and confirm that strong temporal correlation of fish swimming behavior persists in some cases over several hours. A correlated random-walk model was employed to simulate the observed migration behavior, and parameters of the model were estimated that lead to close correspondence between predictions and observations.
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.
Not Available
1993-12-01
The long-range goal of the Numerical Tokamak Project (NTP) is the reliable prediction of tokamak performance using physics-based numerical tools describing tokamak physics. The NTP is accomplishing the development of the most advanced particle and extended fluid model`s on massively parallel processing (MPP) environments as part of a multi-institutional, multi-disciplinary numerical study of tokamak core fluctuations. The NTP is a continuing focus of the Office of Fusion Energy`s theory and computation program. Near-term HPCC work concentrates on developing a predictive numerical description of the core plasma transport in tokamaks driven by low-frequency collective fluctuations. This work addresses one of the greatest intellectual challenges to our understanding of the physics of tokamak performance and needs the most advanced computational resources to progress. We are conducting detailed comparisons of kinetic and fluid numerical models of tokamak turbulence. These comparisons are stimulating the improvement of each and the development of hybrid models which embody aspects of both. The combination of emerging massively parallel processing hardware and algorithmic improvements will result in an estimated 10**2--10**6 performance increase. Development of information processing and visualization tools is accelerating our comparison of computational models to one another, to experimental data, and to analytical theory, providing a bootstrap effect in our understanding of the target physics. The measure of success is the degree to which the experimentally observed scaling of fluctuation-driven transport may be predicted numerically. The NTP is advancing the HPCC Initiative through its state-of-the-art computational work. We are pushing the capability of high performance computing through our efforts which are strongly leveraged by OFE support.
Skibinski, Jakub; Wejrzanowski, Tomasz; Caban, Piotr; Kurzydlowski, Krzysztof J.
2014-10-06
In the present study numerical simulations of epitaxial growth of gallium nitride in Metal Organic Vapor Phase Epitaxy reactor AIX-200/4RF-S is addressed. Epitaxial growth means crystal growth that progresses while inheriting the laminar structure and the orientation of substrate crystals. One of the technological problems is to obtain homogeneous growth rate over the main deposit area. Since there are many agents influencing reaction on crystal area such as temperature, pressure, gas flow or reactor geometry, it is difficult to design optimal process. According to the fact that it's impossible to determine experimentally the exact distribution of heat and mass transfer inside the reactor during crystal growth, modeling is the only solution to understand the process precisely. Numerical simulations allow to understand the epitaxial process by calculation of heat and mass transfer distribution during growth of gallium nitride. Including chemical reactions in numerical model allows to calculate the growth rate of the substrate and estimate the optimal process conditions for obtaining the most homogeneous product.
Martin, V.; Egido, J.L.; Khoo, T.L.; Lauritsen, T.
1995-06-01
The electromagnetic decay of the nuclei {sup 152,154,156}Dy is analyzed using microscopic Hartree-Fock calculations at finite temperature. The theoretical collective transition probabilities are implemented in numerical simulations to produce theoretical spectra. Thermal shape fluctuations are also taken into account. The inclusion of these correlation is crucial in order to understand the main features of the collective {ital E}2 spectra of these isotopes at different energies. The theoretical calculations suggest a shape change as responsible for the unusual features of the spectrum of the nucleus {sup 154}Dy at high energy.
Balakrishna, Jayashree; Bondarescu, Ruxandra; Daues, Gregory; Bondarescu, Mihai
2008-01-15
Excited state soliton stars are studied numerically for the first time. The stability of spherically symmetric S-branch excited state oscillatons under radial perturbations is investigated using a 1D code. We find that these stars are inherently unstable either migrating to the ground state or collapsing to black holes. Higher excited state configurations are observed to cascade through intermediate excited states during their migration to the ground state. This is similar to excited state boson stars [J. Balakrishna, E. Seidel, and W.-M. Suen, Phys. Rev. D 58, 104004 (1998).]. Ground state oscillatons are then studied in full 3D numerical relativity. Finding the appropriate gauge condition for the dynamic oscillatons is much more challenging than in the case of boson stars. Different slicing conditions are explored, and a customized gauge condition that approximates polar slicing in spherical symmetry is implemented. Comparisons with 1D results and convergence tests are performed. The behavior of these stars under small axisymmetric perturbations is studied and gravitational waveforms are extracted. We find that the gravitational waves damp out on a short time scale, enabling us to obtain the complete waveform. This work is a starting point for the evolution of real scalar field systems with arbitrary symmetries.
Rosa, M.P.; Podowski, M.Z.
1995-09-01
This paper is concerned with the analysis of dynamics and stability of boiling channels and systems. The specific objectives are two-fold. One of them is to present the results of a study aimed at analyzing the effects of various modeling concepts and numerical approaches on the transient response and stability of parallel boiling channels. The other objective is to investigate the effect of closed-loop feedback on stability of a boiling water reactor (BWR). Various modeling and computational issues for parallel boiling channels are discussed, such as: the impact of the numerical discretization scheme for the node containing the moving boiling boundary on the convergence and accuracy of computations, and the effects of subcooled boiling and other two-phase flow phenomena on the predictions of marginal stability conditions. Furthermore, the effects are analyzed of local loss coefficients around the recirculation loop of a boiling water reactor on stability of the reactor system. An apparent paradox is explained concerning the impact of changing single-phase losses on loop stability. The calculations have been performed using the DYNOBOSS computer code. The results of DYNOBOSS validation against other computer codes and experimental data are shown.
Brandao, C. S. S.; De Araujo, J. C. N., E-mail: claudiosoriano.uesc@gmail.com, E-mail: jcarlos.dearaujo@inpe.br [Divisao de Astrofisica, Instituto Nacional de Pesquisas Espaciais, S. J. Campos, SP 12227-010 (Brazil)
2012-05-01
A way to probe alternative theories of gravitation is to study if they could account for the structures of the universe. We therefore modified the well-known Gadget-2 code to probe alternative theories of gravitation through galactic dynamics. As an application, we simulate the evolution of spiral galaxies to probe alternative theories of gravitation whose weak field limits have a Yukawa-like gravitational potential. These simulations show that galactic dynamics can be used to constrain the parameters associated with alternative theories of gravitation. It is worth stressing that the recipe given in this study can be applied to any other alternative theory of gravitation in which the superposition principle is valid.
Mauger, Sarah; Colin de Verdière, Guillaume; Bergé, Luc; Skupin, Stefan; Friedrich Schiller University, Institute of Condensed Matter Theory and Optics, 07743 Jena
2013-02-15
A computer cluster equipped with Graphics Processing Units (GPUs) is used for simulating nonlinear optical wave packets undergoing Kerr self-focusing and stimulated Brillouin scattering in fused silica. We first recall the model equations in full (3+1) dimensions. These consist of two coupled nonlinear Schrödinger equations for counterpropagating optical beams closed with a source equation for light-induced acoustic waves seeded by thermal noise. Compared with simulations on a conventional cluster of Central Processing Units (CPUs), GPU-based computations allow us to use a significant (16 times) larger number of mesh points within similar computation times. Reciprocally, simulations employing the same number of mesh points are between 3 and 20 times faster on GPUs than on the same number of classical CPUs. Performance speedups close to 45 are reported for isolated functions evaluating, e.g., the optical nonlinearities. Since the field intensities may reach the ionization threshold of silica, the action of a defocusing electron plasma is also addressed.
Gupta, A.; Moridis, G.J.; Kneafsey, T.J.; Sloan, Jr., E.D.
2009-08-15
The numerical simulator TOUGH+HYDRATE (T+H) was used to predict the transient pure methane hydrate (no sediment) dissociation data. X-ray computed tomography (CT) was used to visualize the methane hydrate formation and dissociation processes. A methane hydrate sample was formed from granular ice in a cylindrical vessel, and slow depressurization combined with thermal stimulation was applied to dissociate the hydrate sample. CT images showed that the water produced from the hydrate dissociation accumulated at the bottom of the vessel and increased the hydrate dissociation rate there. CT images were obtained during hydrate dissociation to confirm the radial dissociation of the hydrate sample. This radial dissociation process has implications for dissociation of hydrates in pipelines, suggesting lower dissociation times than for longitudinal dissociation. These observations were also confirmed by the numerical simulator predictions, which were in good agreement with the measured thermal data during hydrate dissociation. System pressure and sample temperature measured at the sample center followed the CH{sub 4} hydrate L{sub w}+H+V equilibrium line during hydrate dissociation. The predicted cumulative methane gas production was within 5% of the measured data. Thus, this study validated our simulation approach and assumptions, which include stationary pure methane hydrate-skeleton, equilibrium hydrate-dissociation and heat- and mass-transfer in predicting hydrate dissociation in the absence of sediments. It should be noted that the application of T+H for the pure methane hydrate system (no sediment) is outside the general applicability limits of T+H.
Oudini, N.; Taccogna, F.; Aanesland, A.
2014-06-15
Laser photo-detachment is used as a method to measure or determine the negative ion density and temperature in electronegative plasmas. In essence, the method consists of producing an electropositive channel (negative ion free region) via pulsed laser photo-detachment within an electronegative plasma bulk. Electrostatic probes placed in this channel measure the change in the electron density. A second pulse might be used to track the negative ion recovery. From this, the negative ion density and temperature can be determined. We study the formation and relaxation of the electropositive channel via a two-dimensional Particle-In-Cell/Mote Carlo collision model. The simulation is mainly carried out in a Hydrogen plasma with an electronegativity of α = 1, with a parametric study for α up to 20. The temporal and spatial evolution of the plasma potential and the electron densities shows the formation of a double layer (DL) confining the photo-detached electrons within the electropositive channel. This DL evolves into two fronts that move in the opposite directions inside and outside of the laser spot region. As a consequence, within the laser spot region, the background and photo-detached electron energy distribution function relaxes/thermalizes via collisionless effects such as Fermi acceleration and Landau damping. Moreover, the simulations show that collisional effects and the DL electric field strength might play a non-negligible role in the negative ion recovery within the laser spot region, leading to a two-temperature negative ion distribution. The latter result might have important effects in the determination of the negative ion density and temperature from laser photo detachment diagnostic.
Richard W. Johnson
2005-09-01
The final design of the very high temperature reactor (VHTR) of the fourth generation of nuclear power plants (Gen IV) has not yet been established. The VHTR may be either a prismatic (block) or pebble bed type. It may be either gas-cooled or cooled with an as yet unspecified molten salt. However, a conceptual design of a gas-cooled VHTR, based on the General Atomics GT-MHR, does exist and is called the prismatic VHTR reference design, MacDonald et al [2003], General Atomics [1996]. The present validation studies are based on the prismatic VHTR reference design. In the prismatic VHTR reference design, the flow in the lower plenum will be introduced by dozens of turbulent jets issuing into a large crossflow that must negotiate dozens of cylindrical support columns as it flows toward the exit duct of the reactor vessel. The jets will not all be at the same temperature due to the radial variation of power density expected in the core. However, it is important that the coolant be well mixed when it enters the power conversion unit to ensure proper operation and long life of the power conversion machinery. Hence, it is deemed important to be able to accurately model the flow and mixing of the variable temperature coolant in the lower plenum and exit duct. Accurate flow modeling involves determining modeling strategies including the fineness of the grid needed, iterative convergence tolerance, numerical discretization method used, whether the flow is steady or unsteady, and the turbulence model and wall treatment employed. It also involves validation of the computer code and turbulence model against a series of separate and combined flow phenomena and selection of the data used for the validation. The present report describes progress made to date for the task entitled ‘CFD software validation of jets in crossflow’ which was designed to investigate the issues pertaining to the validation process.
Gokaltun, Seckin; Munroe, Norman; Subramaniam, Shankar
2014-12-31
This study presents a new drag model, based on the cohesive inter-particle forces, implemented in the MFIX code. This new drag model combines an existing standard model in MFIX with a particle-based drag model based on a switching principle. Switches between the models in the computational domain occur where strong particle-to-particle cohesion potential is detected. Three versions of the new model were obtained by using one standard drag model in each version. Later, performance of each version was compared against available experimental data for a fluidized bed, published in the literature and used extensively by other researchers for validation purposes. In our analysis of the results, we first observed that standard models used in this research were incapable of producing closely matching results. Then, we showed for a simple case that a threshold is needed to be set on the solid volume fraction. This modification was applied to avoid non-physical results for the clustering predictions, when governing equation of the solid granular temperate was solved. Later, we used our hybrid technique and observed the capability of our approach in improving the numerical results significantly; however, improvement of the results depended on the threshold of the cohesive index, which was used in the switching procedure. Our results showed that small values of the threshold for the cohesive index could result in significant reduction of the computational error for all the versions of the proposed drag model. In addition, we redesigned an existing circulating fluidized bed (CFB) test facility in order to create validation cases for clustering regime of Geldart A type particles.
Michael David Petersen
2001-05-01
Using computer simulations, the performance of several CdTe based photovoltaic structures has been studied. The advantages and disadvantages of band gap grading, through the use of (Zn,Cd)Te, have also been investigated in these structures. Grading at the front interface between a CdS window layer and a CdTe absorber layer, can arise due to interdiffusion between the materials during growth or due to the intentional variation of the material composition. This grading has been shown to improve certain performance metrics, such as the open-circuit voltage, while degrading others, such as the fill factor, depending on the amount and distance of the grading. The presence of a Schottky barrier as the back contact has also been shown to degrade the photovoltaic performance of the device, resulting in a characteristic IV curve. However, with the appropriate band gap grading at the back interface, it has been shown that the performance can be enhanced through more efficient carrier collection. These results were then correlated with experimental observations of the performance degradation in devices subjected to light and heat stress.
Michael David Petersen
2001-06-27
Using computer simulations, the performance of several CdTe based photovoltaic structures has been studied. The advantages and disadvantages of band gap grading, through the use of (Zn, Cd)Te, have also been investigated in these structures. Grading at the front interface between a CdS window layer and a CdTe absorber layer, can arise due to interdiffusion between the materials during growth or due to the intentional variation of the material composition. This grading has been shown to improve certain performance metrics, such as the open-circuit voltage, while degrading others, such as the fill factor, depending on the amount and distance of the grading. The presence of a Schottky barrier as the back contact has also been shown to degrade the photovoltaic performance of the device, resulting in a characteristic IV curve. However, with the appropriate band gap grading at the back interface, it has been shown that the performance can be enhanced through more efficient carrier collection. These results were then correlated with experimental observations of the performance degradation in devices subjected to light and heat stress.
Pruess, K.; Nordbotten, J.
2010-12-28
We have used the TOUGH2-MP/ECO2N code to perform numerical simulation studies of the long-term behavior of CO{sub 2} stored in an aquifer with a sloping caprock. This problem is of great practical interest, and is very challenging due to the importance of multi-scale processes. We find that the mechanism of plume advance is different from what is seen in a forced immiscible displacement, such as gas injection into a water-saturated medium. Instead of pushing the water forward, the plume advances because the vertical pressure gradients within the plume are smaller than hydrostatic, causing the groundwater column to collapse ahead of the plume tip. Increased resistance to vertical flow of aqueous phase in anisotropic media leads to reduced speed of updip plume advancement. Vertical equilibrium models that ignore effects of vertical flow will overpredict the speed of plume advancement. The CO{sub 2} plume becomes thinner as it advances, yet the speed of advancement remains constant over the entire simulation period of up to 400 years, with migration distances of more than 80 km. Our simulations include dissolution of CO{sub 2} into the aqueous phase and associated density increase, and molecular diffusion. However, no convection develops in the aqueous phase because it is suppressed by the relatively coarse (sub-) horizontal gridding required in a regional-scale model. A first crude sub-grid-scale model was developed to represent convective enhancement of CO{sub 2} dissolution. This process is found to greatly reduce the thickness of the CO{sub 2} plume, but, for the parameters used in our simulations, does not affect the speed of plume advancement.
Kaindl, W.; Lades, M.; Wachutka, G.; Kaminski, N.; Niemann, E.
1999-03-01
Silicon carbide (SiC) receives strong attention for high-power, high-temperature, and high-frequency device applications due to its promising material properties. In order to measure the ionization time constants of dopants in 4H/6H-SiC within a wide range of temperature, nitrogen (N), aluminum (Al), and boron (B) have been characterized using thermal admittance spectroscopy (AS) and deep level transient spectroscopy (DLTS). The temperature extrapolation of the results obtained by AS shows excellent agreement with those obtained by DLTS, yielding the base for an evaluation of incomplete ionization effects in SiC devices within usual operation ranges. The measured data has been analyzed using numerical drift-diffusion simulations based on the method of finite-elements. A numerical investigation of the different freeze-out characteristics of free carriers in p{sup +}n, n{sup +}p, and Schottky diodes shows that unlike in the case of B, the ionization time constant of Al can be exclusively measured in the highly doped region of a p{sup +}n diode.
Kanevce, Ana; Kuciauskas, Darius; Levi, Dean H.; Johnston, Steven W.; Allende Motz, Alyssa M.
2015-07-28
We use two-dimensional numerical simulations to analyze high spatial resolution time-resolved spectroscopy data. This analysis is applied to two-photon excitation time-resolved photoluminescence (2PE-TRPL) but is broadly applicable to all microscopic time-resolved techniques. By solving time-dependent drift-diffusion equations, we gain insight into carrier dynamics and transport characteristics. Accurate understanding of measurement results establishes the limits and potential of the measurement and enhances its value as a characterization method. Diffusion of carriers outside of the collection volume can have a significant impact on the measured decay but can also provide an estimate of carrier mobility as well as lifetime. In addition to material parameters, the experimental conditions, such as spot size and injection level, can impact the measurement results. Although small spot size provides better resolution, it also increases the impact of diffusion on the decay; if the spot size is much smaller than the diffusion length, it impacts the entire decay. By reproducing experimental 2PE-TRPL decays, the simulations determine the bulk carrier lifetime from the data. The analysis is applied to single-crystal and heteroepitaxial CdTe, material important for solar cells, but it is also applicable to other semiconductors where carrier diffusion from the excitation volume could affect experimental measurements.
Modeling-Computer Simulations At U.S. West Region (Laney, 2005...
Hydrothermal Systems in the Basin and Range and Cordilleran United States, Moore, Nash, Nemcok, Lutz, Norton, Kaspereit, Berard, van de Putte, Johnson and Deymonaz....
Fairchild, B.T.
1987-01-01
These proceedings contain papers on simulators with artificial intelligence, and the human decision making process; visuals for simulators: human factors, training, and psycho-physical impacts; the role of institutional structure on simulation projects; maintenance trainers for economic value and safety; biomedical simulators for understanding nature, for medical benefits, and the physiological effects of simulators; the mathematical models and numerical techniques that drive today's simulators; and the demography of simulators, with census papers identifying the population of real-time simulator training devices; nuclear reactors.
Sun, Rui; Ismail, Tamer M.; Ren, Xiaohan; Abd El-Salam, M.
2015-05-15
Highlights: • The effects of moisture content on the burning process of MSW are investigated. • A two-dimensional mathematical model was built to simulate the combustion process. • Temperature distributions, process rates, gas species were measured and simulated. • The The conversion ratio of C/CO and N/NO in MSW are inverse to moisture content. - Abstract: In order to reveal the features of the combustion process in the porous bed of a waste incinerator, a two-dimensional unsteady state model and experimental study were employed to investigate the combustion process in a fixed bed of municipal solid waste (MSW) on the combustion process in a fixed bed reactor. Conservation equations of the waste bed were implemented to describe the incineration process. The gas phase turbulence was modeled using the k–ε turbulent model and the particle phase was modeled using the kinetic theory of granular flow. The rate of moisture evaporation, devolatilization rate, and char burnout was calculated according to the waste property characters. The simulation results were then compared with experimental data for different moisture content of MSW, which shows that the incineration process of waste in the fixed bed is reasonably simulated. The simulation results of solid temperature, gas species and process rate in the bed are accordant with experimental data. Due to the high moisture content of fuel, moisture evaporation consumes a vast amount of heat, and the evaporation takes up most of the combustion time (about 2/3 of the whole combustion process). The whole bed combustion process reduces greatly as MSW moisture content increases. The experimental and simulation results provide direction for design and optimization of the fixed bed of MSW.
Garcia, C. Amanda; Halford, Keith J.; Laczniak, Randell J.
2010-02-12
Hydraulic conductivities of volcanic and carbonate lithologic units at the Nevada Test Site were estimated from flow logs and aquifer-test data. Borehole flow and drawdown were integrated and interpreted using a radial, axisymmetric flow model, AnalyzeHOLE. This integrated approach is used because complex well completions and heterogeneous aquifers and confining units produce vertical flow in the annular space and aquifers adjacent to the wellbore. AnalyzeHOLE simulates vertical flow, in addition to horizontal flow, which accounts for converging flow toward screen ends and diverging flow toward transmissive intervals. Simulated aquifers and confining units uniformly are subdivided by depth into intervals in which the hydraulic conductivity is estimated with the Parameter ESTimation (PEST) software. Between 50 and 150 hydraulic-conductivity parameters were estimated by minimizing weighted differences between simulated and measured flow and drawdown. Transmissivity estimates from single-well or multiple-well aquifer tests were used to constrain estimates of hydraulic conductivity. The distribution of hydraulic conductivity within each lithology had a minimum variance because estimates were constrained with Tikhonov regularization. AnalyzeHOLE simulated hydraulic-conductivity estimates for lithologic units across screened and cased intervals are as much as 100 times less than those estimated using proportional flow-log analyses applied across screened intervals only. Smaller estimates of hydraulic conductivity for individual lithologic units are simulated because sections of the unit behind cased intervals of the wellbore are not assumed to be impermeable, and therefore, can contribute flow to the wellbore. Simulated hydraulic-conductivity estimates vary by more than three orders of magnitude across a lithologic unit, indicating a high degree of heterogeneity in volcanic and carbonate-rock units. The higher water transmitting potential of carbonate-rock units relative
LaVenue, A.M.; Haug, A.; Kelley, V.A.
1988-03-01
This hydrogeologic modeling study has been performed as part of the regional hydrologic characterization of the Waste Isolation Pilot Plant (WIPP) Site in southeastern New Mexico. The study resulted in an estimation of the transmissivity distrubution, hydraulic potentials, flow field, and fluid densities in the Culebra Dolomite Member of the Permian Rustler Formation at the WIPP site. The three-dimensional finite-difference code SWIFT-II was employed for the numerical modeling, using variable-fluid-density and a single-porosity formulation. The modeled area includes and extends beyond the WIPP controlled zone (Zone 3). The work performed consisted of modeling the hydrogeology of the Culebra using two approaches: (1) steady-state modeling to develop the best estimate of the undisturbed head distribution, i.e., of the situation before sinking if the WIPP shafts, which began in 1981; and (2) superimposed transient modeling of local hydrologic responses to excavation of the three WIPP shafts at the center of the WIPP site, as well as to various well tests. Boundary conditions (prescribed constant fluid pressures and densities) were estimated using hydraulic-head and fluid-density data obtained from about 40 wells at and near the WIPP site. The transient modeling used the calculated steady-state freshwater heads as initial conditions. 107 refs., 112 figs., 22 tabs.
Rundle, John B.
2004-12-31
Physical Review Letters (Tiampo et al., in press), demonstrate that the Southern California system is ergodic in the same way that is seen in the models. These results will be discussed in more detail below. However, the point that needs to be emphasized is that it was the combination of model investigation via theory and simulation coupled with assimilation and classification of real data and applying the methods of statistical mechanics to real fault systems that led to both a successful forecasting algorithm and a deeper understanding of the nature of earthquake fault systems. This paper describes in some detail the results obtained in the previous funding period. We present these in three groups. (A) Investigation of statistical physics models and applications. (B) Earthquake fault systems and Greens functions for complex sources and (C) Space time patterns, data analysis and forecasting.