Training and Research on Probabilistic Hydro-Thermo-Mechanical...
Office of Scientific and Technical Information (OSTI)
Training and Research on Probabilistic Hydro-Thermo-Mechanical Modeling of Carbon Dioxide Geological Sequestration in Fractured Porous Rocks Citation Details In-Document Search...
A Hydro-Thermo-Mechanical Numerical Model For Hdr Geothermal...
Published Journal International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, 1996 DOI 10.10160148-9062(96)00002-2 Citation T. W. Hicks,R. J. Pine,J....
A Hydro-Thermo-Mechanical Numerical Model For Hdr Geothermal Reservoir
AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on DeliciousPlasmaP a gHigh4-FD-a < RAPIDâ€Žcommunication facilities | Open33 CFREnergyEnergyResultsAonEvaluation |
Yu, Dunji [ORNL; An, Ke [ORNL; Gao, Carrie Y [ORNL; Heller, William T [ORNL; Chen, Xu [Tianjin UNiversity, China
2013-01-01T23:59:59.000Z
A portable hydro-thermo-mechanical loading cell has been designed to enable in-situ small angle neutron scattering (SANS) studies of proton exchange membranes (PEM) under immersed tensile loadings at different temperatures. The cell consists of three main parts as follows: a letter-paper-size motor-driven mechanical load frame, a SANS friendly reservoir that provides stable immersed and thermal sample conditions, and a data acquisition & control system. The ex-situ tensile tests of Nafion 212 membranes demonstrated a satisfactory thermo-mechanical testing performance of the cell for either dry or immersed conditions at elevated temperatures. The in-situ SANS tensile measurements on the Nafion 212 membranes immersed in D2O at 70oC proved the feasibility and capability of the cell for small angle scattering study on deformation behaviors of PEM and other polymer materials under hydro-thermo-mechanical loading.
Training and Research on Probabilistic Hydro-Thermo-Mechanical Modeling of
Office of Scientific and Technical Information (OSTI)
AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE:1 First Use of Energy for All Purposes (Fuel and Nonfuel), 2002; Level:5 TablesExports to3,1,50022,3,,0,,6,1,Separation 23Tribal EnergyCatalytic CobySpeeding accessProjectDecadal/Centennial Scales (Technical Report)
Training and Research on Probabilistic Hydro-Thermo-Mechanical Modeling of
Office of Scientific and Technical Information (OSTI)
AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE:1 First Use of Energy for All Purposes (Fuel and Nonfuel), 2002; Level:5 TablesExports to3,1,50022,3,,0,,6,1,SeparationConnect1.08]Te[subscriptM-PACEResponses to aConnect Towards an explicitCarbon Dioxide Geological
Numerical wind speed simulation model
Ramsdell, J.V.; Athey, G.F.; Ballinger, M.Y.
1981-09-01T23:59:59.000Z
A relatively simple stochastic model for simulating wind speed time series that can be used as an alternative to time series from representative locations is described in this report. The model incorporates systematic seasonal variation of the mean wind, its standard deviation, and the correlation speeds. It also incorporates systematic diurnal variation of the mean speed and standard deviation. To demonstrate the model capabilities, simulations were made using model parameters derived from data collected at the Hanford Meteorology Station, and results of analysis of simulated and actual data were compared.
INCORPORATION OF GROUNDWATER FLOW INTO NUMERICAL MODELS AND DESIGN MODELS
03/10/99 1 INCORPORATION OF GROUNDWATER FLOW INTO NUMERICAL MODELS AND DESIGN MODELS Jeffrey D-coupled, ground-source heat pumps, groundwater, heat pump, heat exchanger, heat transfer, numerical models transport of heat by moving groundwater may be an important factor in reducing the necessary size of closed
High performance computing and numerical modelling
,
2014-01-01T23:59:59.000Z
Numerical methods play an ever more important role in astrophysics. This is especially true in theoretical works, but of course, even in purely observational projects, data analysis without massive use of computational methods has become unthinkable. The key utility of computer simulations comes from their ability to solve complex systems of equations that are either intractable with analytic techniques or only amenable to highly approximative treatments. Simulations are best viewed as a powerful complement to analytic reasoning, and as the method of choice to model systems that feature enormous physical complexity such as star formation in evolving galaxies, the topic of this 43rd Saas Fee Advanced Course. The organizers asked me to lecture about high performance computing and numerical modelling in this winter school, and to specifically cover the basics of numerically treating gravity and hydrodynamics in the context of galaxy evolution. This is still a vast field, and I necessarily had to select a subset ...
Lattice Boltzmann Model for Numerical Relativity
Ilseven, E
2015-01-01T23:59:59.000Z
In the Bona-Masso formulation, Einstein equations are written as a set of flux conservative first order hyperbolic equations that resemble fluid dynamics equations. Based on this formulation, we construct a lattice Boltzmann model for Numerical Relativity. Our model is validated with well-established tests, showing good agreement with analytical solutions. Furthermore, we show that by increasing the relaxation time, we gain stability at the cost of losing accuracy, and by decreasing the lattice spacings while keeping a constant numerical diffusivity, the accuracy and stability of our simulations improves. Finally, in order to show the potential of our approach a linear scaling law for parallelisation with respect to number of CPU cores is demonstrated. Our model represents the first step in using lattice kinetic theory to solve gravitational problems.
A numerical model of perturbation gas chromatography
DeBarro, Marc Joseph
1985-01-01T23:59:59.000Z
for in the model presented by Glover and Lsu. These mass transfer resistances are accommodated through an overall gss side mass transfer coefficient. Numerically generated peaks are compared to experimental profiles to con- firm the validity to the equations... and the method. Various solvent/polybutadiene systems are explored and individual profiles of the species are examined to under- stand the interaction which occurs between components for multicomponent sys- tems. Both ternary and quaternary systems...
Numerical modeling of vertical cavity semiconductor lasers
Chow, W.W.; Hadley, G.R.
1996-08-01T23:59:59.000Z
A vertical cavity surface emitting laser (VCSEL) is a diode laser whose optical cavity is formed by growing or depositing DBR mirror stacks that sandwich an active gain region. The resulting short cavity supports lasing into a single longitudinal mode normal to the wafer, making these devices ideal for a multitude of applications, ranging from high-speed communication to high-power sources (from 2D arrays). This report describes the development of a numerical VCSEL model, whose goal is to both further their understanding of these complex devices and provide a tool for accurate design and data analysis.
NUMERICAL MODELING OF CATHODE CONTACT MATERIAL DENSIFICATION
Koeppel, Brian J.; Liu, Wenning N.; Stephens, Elizabeth V.; Khaleel, Mohammad A.
2011-11-01T23:59:59.000Z
Numerical modeling was used to simulate the constrained sintering process of the cathode contact layer during assembly of solid oxide fuel cells (SOFCs). A finite element model based on the continuum theory for sintering of porous bodies was developed and used to investigate candidate low-temperature cathode contact materials. Constitutive parameters for various contact materials under investigation were estimated from dilatometry screening tests, and the influence of processing time, processing temperature, initial grain size, and applied compressive stress on the free sintering response was predicted for selected candidate materials. The densification behavior and generated stresses within a 5-cell planar SOFC stack during sintering, high temperature operation, and room temperature shutdown were predicted. Insufficient constrained densification was observed in the stack at the proposed heat treatment, but beneficial effects of reduced grain size, compressive stack preload, and reduced thermal expansion coefficient on the contact layer densification and stresses were observed.
Numerical Modeling At Dixie Valley Geothermal Area (McKenna ...
models were made to analyze the transient behavior of extensional basin and range geothermal systems. Notes Numerical models were made utilizing a range of bulk rock...
A numerical model of aerosol scavenging
Bradley, M.M.; Molenkamp, C.R.
1991-10-01T23:59:59.000Z
Using a three-dimensional numerical cloud/smoke-plume model, we have simulated the burning of a large, mid-latitude city following a nuclear exchange. The model includes 18 dynamic and microphysical equations that predict the fire-driven airflow, cloud processes, and smoke-cloud interactions. In the simulation, the intense heating from the burning city produces a firestorm with updraft velocities exceeding 60 m/s. Within 15 minutes of ignition, the smoke plume penetrates the tropopause. The updraft triggers a cumulonimbus cloud that produces significant quantities of ice, snow, and hail. These solid hydrometeors, as well as cloud droplets and rain, interact with the smoke particles from the fire. At the end of the one-hour simulation, over 20% of the smoke is in slowly falling snowflakes. If the snow reaches the ground before the flakes completely sublimate (or melt and then evaporate), then only approximately 50% of the smoke will survive the scavenging processes and remain in the atmosphere to affect the global climate.
Numerical Modeling of the Nucleation Conditions of Petal-Centerline...
Coso Geothermal Field Jump to: navigation, search OpenEI Reference LibraryAdd to library Conference Proceedings: Numerical Modeling of the Nucleation Conditions of...
Evolution of the Bohemian Massif: Insights from numerical modeling
Cerveny, Vlastislav
Evolution of the Bohemian Massif: Insights from numerical modeling Petra Maierová Supervisor: Doc of Geophysics Faculty of Mathematics and Physics Charles University in Prague #12;February 4, 2013Evolution Conclusions Outline #12;February 4, 2013Evolution of the Bohemian Massif: Insights from numerical modeling 3
Numerical Modeling of Eastern Connecticut's Visual Resources1
Conservation, The University of Connecticut, Storrs, Connecticut. / Abstract: A numerical model capableNumerical Modeling of Eastern Connecticut's Visual Resources1 Daniel L. Civco 2/ l Submitted of accurately predicting the preference for landscape photographs of selected points in eastern Connecticut
Buoyant melting instabilities beneath extending lithosphere: 1. Numerical models
Tackley, Paul J.
Buoyant melting instabilities beneath extending lithosphere: 1. Numerical models John W. Hernlund,1,2 Paul J. Tackley,1,3 and David J. Stevenson4 Received 18 November 2006; revised 18 October 2007 diffusely extending lithosphere is studied using numerical convection models covering a wide range
Starting laminar plumes: Comparison of laboratory and numerical modeling
van Keken, Peter
Starting laminar plumes: Comparison of laboratory and numerical modeling Judith Vatteville Institut, France (davaille@fast.u-psud.fr) [1] A detailed comparison of starting laminar plumes in viscous fluids. Davaille (2009), Starting laminar plumes: Comparison of laboratory and numerical modeling, Geochem. Geophys
A numerical model simulation of longshore transport for Galveston Island
Gilbreath, Stephen Alexander
1995-01-01T23:59:59.000Z
The shoreline changes, deposition patterns, and longshore transport rates were calculated for the coast of Galveston Island using a numerical model simulation. The model only simulated changes due to waves creating longshore currents. East Beach...
NUMERICAL MODELING FOR MULTIPHASE INCOMPRESSIBLE FLOW WITH PHASE CHANGE
Abdou, Mohamed
NUMERICAL MODELING FOR MULTIPHASE INCOMPRESSIBLE FLOW WITH PHASE CHANGE Xiao-Yong Luo, Ming-Jiu Ni for multiphase flows. A con- tinuum surface force (CSF) tension model is used in the present cases. Phase change
Numerical models of phosphate esters in the Chattahoochee River
Haffey, Samuel Fraad, 1973-
2004-01-01T23:59:59.000Z
A numerical model was constructed to assess the magnitude of organophosphoric acid triester sinks in the Chattahoochee River and to identify concentration patterns downstream of point source discharges. The model was built ...
Data Assimilation for Idealised Mathematical Models of Numerical Weather Prediction
Wirosoetisno, Djoko
Data Assimilation for Idealised Mathematical Models of Numerical Weather Prediction Supervisors). Background: Numerical Weather Prediction (NWP) has seen significant gains in accuracy in recent years due in weather dynamics, e.g., the asymptotic balance seen in high and low pressure systems. Aims of the project
Fully kinetic numerical modeling of a plasma thruster
Szabo, James Joseph, 1969-
2001-01-01T23:59:59.000Z
A Hall effect plasma thruster with conductive acceleration channel walls was numerically modeled using 2D3V Particle-in-Cell (PIC) and Monte-Carlo Collision (MCC) methodolo- gies. Electron, ion, and neutral dynamics were ...
Numerical studies of a simple Coulomb blockade model
Shao, Jianfeng
1991-01-01T23:59:59.000Z
NUMERICAL STUDIES OF A SIMPLE COULOMB BLOCKADE MODEL A Thesis by JIANFENG SHAO Submitted to the Office of Graduate Studies of Texas ASM University in partial fulfillment of the requirements for the degree of MASTER OF SCIENCE August 1991... Major Subject: Physics NUMERICAL STUDIES OF A SIMPLE COULOMB BLOCKADE MODEL A Thesis by JIANFENG SHAO Approved as to style and content by: Roland E, Allen (Chair of Committee) /, 1 r oseph H. R s ( Member) Chin B. Su (Member) Richard L...
Numerical Modeling of Brain Dynamics in Traumatic Situations -Impulsive Translations
Burtscher, Martin
.S.A. Abstract We numerically model the brain dy- namics during and after impulsive head translations using brain injuries appear among boxers and shaken babies despite minimal rotations of their heads. Modeling head translations also helps understand the brain dynamics during head rotations about an arbitrary
Numerical modelling and analysis of a room temperature magnetic
Numerical modelling and analysis of a room temperature magnetic refrigeration system Thomas Frank and analysis of a room temperature magnetic refrigeration system Department: Fuel Cells and Solid State-dimensional mathematical model of an Active Magnetic Regenerator (AMR) system which is used for magnetic refrigeration
Numerical modelling of tsunami mitigation by mangroves Putu Harry Gunawan
Mancini, Simona
Numerical modelling of tsunami mitigation by mangroves Putu Harry Gunawan LAMA (Laboratoire d'Analyse et de Mathmatiques Appliques) UPEM putu-harry.gunawan@univ.paris-est.fr Abstract Figure 1: Mangrove-Tsunami Model. The role of mangroves (coastal forests) in the mitigation of tsunami impacts is a debated topic
Diffusion and Dispersion Characterization of a Numerical Tsunami Model
Tolkova, Elena
Diffusion and Dispersion Characterization of a Numerical Tsunami Model David Burwell, Elena Tolkova and Arun Chawla February 6, 2007 1 Introduction Method Of Splitting Tsunami (MOST) is a depth averaged long wave tsunami inundation model that was originally developed by Titov and Synolakis (1995) for 1D
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE:1 First Use of Energy for All Purposes (Fuel and Nonfuel), 2002; Level:5(Million Cubic Feet) Oregon (Including Vehicle Fuel) (Million Cubic Feet) Natural Gas DeliveredStents -(NEI) ScienceNumerical simulations of current
Direct Numerical Simulations and Modeling of Jets in Crossflow
Mahesh, Krishnan
Direct Numerical Simulations and Modeling of Jets in Crossflow A THESIS SUBMITTED TO THE FACULTY. i #12;To my parents and my grandparents, and to Ramnath ii #12;Abstract Jets in crossflow are used to study the different aspects of round jets in a crossflow. The first problem studies
Validation of Erosion Modeling: Physical and Numerical Mehrad Kamalzare1
Franklin, W. Randolph
-3590 ABSTRACT The overall intent of this research is to develop numerical models of erosion of levees, dams is necessary for emergency plans for levee or dam breaches. Griffis, 2007 addressed the overall design and hydraulic shear stress. #12;2 Xu and Zhang (2009) found that in addition to soil type, the degree
Friction versus dilation revisited: Insights from theoretical and numerical models
Einat, Aharonov
Friction versus dilation revisited: Insights from theoretical and numerical models N. Makedonska,1 controlled by the frictional strength of the fault gouge, a granular layer that accumulates between the fault friction coefficient) of such granular layers is the systems resistance to dilation, a byprocess
Numerical Modeling of the 2011 Tohoku Earthquake Tsunami
Frandsen, Jannette B.
Evolution of Ocean WAVEs) to reconstruct the tsunami across the Pacific and its transformation aroundNumerical Modeling of the 2011 Tohoku Earthquake Tsunami Yoshiki Yamazaki Post-doctoral Research 2011 Tohoku earthquake (Mw 9.0) generated a massive tsunami devastated the entire Pacific coast
Small Buoys for Energy Harvesting : Experimental and Numerical Modeling Studies
Grilli, Stéphan T.
Small Buoys for Energy Harvesting : Experimental and Numerical Modeling Studies St´ephan T. Grilli for wave energy harvesting (free-floating or slackly moored), to produce about 1 KW per unit at full scale-contained (water tight) resonating multiple-spar buoy (or Starspar), in which a longer central spar houses the LEG
Integrating Numerical Computation into the Modeling Instruction Curriculum
Caballero, Marcos D; Aiken, John M; Douglas, Scott S; Scanlon, Erin M; Thoms, Brian; Schatz, Michael F
2012-01-01T23:59:59.000Z
We describe a way to introduce physics high school students with no background in programming to computational problem-solving experiences. Our approach builds on the great strides made by the Modeling Instruction reform curriculum. This approach emphasizes the practices of "Developing and using models" and "Computational thinking" highlighted by the NRC K-12 science standards framework. We taught 9th-grade students in a Modeling-Instruction-based physics course to construct computational models using the VPython programming environment. Numerical computation within the Modeling Instruction curriculum provides coherence among the curriculum's different force and motion models, links the various representations which the curriculum employs, and extends the curriculum to include real-world problems that are inaccessible to a purely analytic approach.
Refined numerical models for multidimensional Type Ia supernova simulations
Reinecke, M; Niemeyer, J C
2002-01-01T23:59:59.000Z
Following up on earlier work on this topic (Reinecke et al. 1999, A&A 347, pp. 724 and 739), we present an improved set of numerical models for simulations of white dwarfs exploding as Type Ia supernovae (SNe Ia). Two-dimensional simulations were used to test the reliability and numerical robustness of these algorithms; the results indicate that integral quantities like the total energy release are insensitive to changes of the grid resolution (above a certain threshold), which was not the case for our former code. The models were further enhanced to allow fully three-dimensional simulations of SNe Ia. A direct comparison of a 2D and a 3D calculation with identical initial conditions shows that the explosion is considerably more energetic in three dimensions; this is most likely caused by the assumption of axisymmetry in 2D, which inhibits the growth of flame instabilities in the azimuthal direction and thereby decreases the flame surface.
Refined numerical models for multidimensional Type Ia supernova simulations
M. Reinecke; W. Hillebrandt; J. C. Niemeyer
2001-11-26T23:59:59.000Z
Following up on earlier work on this topic (Reinecke et al. 1999, A&A 347, pp. 724 and 739), we present an improved set of numerical models for simulations of white dwarfs exploding as Type Ia supernovae (SNe Ia). Two-dimensional simulations were used to test the reliability and numerical robustness of these algorithms; the results indicate that integral quantities like the total energy release are insensitive to changes of the grid resolution (above a certain threshold), which was not the case for our former code. The models were further enhanced to allow fully three-dimensional simulations of SNe Ia. A direct comparison of a 2D and a 3D calculation with identical initial conditions shows that the explosion is considerably more energetic in three dimensions; this is most likely caused by the assumption of axisymmetry in 2D, which inhibits the growth of flame instabilities in the azimuthal direction and thereby decreases the flame surface.
On numerical considerations for modeling reactive astrophysical shocks
Papatheodore, Thomas L.; Messer, O. E. Bronson, E-mail: tpapathe@utk.edu, E-mail: bronson@ornl.gov [Department of Physics and Astronomy, University of Tennessee, Knoxville, TN 37996 (United States)
2014-02-10T23:59:59.000Z
Simulating detonations in astrophysical environments is often complicated by numerical approximations to shock structure. A common prescription to ensure correct detonation speeds and associated quantities is to prohibit burning inside the numerically broadened shock. We have performed a series of simulations to verify the efficacy of this approximation and to understand how resolution and dimensionality might affect its use. Our results show that in one dimension, prohibiting burning in the shock is important wherever the carbon burning length is not resolved, in keeping with the results of Fryxell et al. In two dimensions, we find that the prohibition of shock burning effectively inhibits the development of cellular structure for all but the most highly resolved cases. We discuss the possible impacts this outcome may have on sub-grid models and detonation propagation in models of Type Ia supernovae, including potential impacts on observables.
ASSIMILATION OF DOPPLER RADAR DATA INTO NUMERICAL WEATHER MODELS
Chiswell, S.; Buckley, R.
2009-01-15T23:59:59.000Z
During the year 2008, the United States National Weather Service (NWS) completed an eight fold increase in sampling capability for weather radars to 250 m resolution. This increase is expected to improve warning lead times by detecting small scale features sooner with increased reliability; however, current NWS operational model domains utilize grid spacing an order of magnitude larger than the radar data resolution, and therefore the added resolution of radar data is not fully exploited. The assimilation of radar reflectivity and velocity data into high resolution numerical weather model forecasts where grid spacing is comparable to the radar data resolution was investigated under a Laboratory Directed Research and Development (LDRD) 'quick hit' grant to determine the impact of improved data resolution on model predictions with specific initial proof of concept application to daily Savannah River Site operations and emergency response. Development of software to process NWS radar reflectivity and radial velocity data was undertaken for assimilation of observations into numerical models. Data values within the radar data volume undergo automated quality control (QC) analysis routines developed in support of this project to eliminate empty/missing data points, decrease anomalous propagation values, and determine error thresholds by utilizing the calculated variances among data values. The Weather Research and Forecasting model (WRF) three dimensional variational data assimilation package (WRF-3DVAR) was used to incorporate the QC'ed radar data into input and boundary conditions. The lack of observational data in the vicinity of SRS available to NWS operational models signifies an important data void where radar observations can provide significant input. These observations greatly enhance the knowledge of storm structures and the environmental conditions which influence their development. As the increase in computational power and availability has made higher resolution real-time model simulations possible, the need to obtain observations to both initialize numerical models and verify their output has become increasingly important. The assimilation of high resolution radar observations therefore provides a vital component in the development and utility of numerical model forecasts for both weather forecasting and contaminant transport, including future opportunities to improve wet deposition computations explicitly.
Mathematical and Numerical Analyses of Peridynamics for Multiscale Materials Modeling
Du, Qiang
2014-11-12T23:59:59.000Z
The rational design of materials, the development of accurate and efficient material simulation algorithms, and the determination of the response of materials to environments and loads occurring in practice all require an understanding of mechanics at disparate spatial and temporal scales. The project addresses mathematical and numerical analyses for material problems for which relevant scales range from those usually treated by molecular dynamics all the way up to those most often treated by classical elasticity. The prevalent approach towards developing a multiscale material model couples two or more well known models, e.g., molecular dynamics and classical elasticity, each of which is useful at a different scale, creating a multiscale multi-model. However, the challenges behind such a coupling are formidable and largely arise because the atomistic and continuum models employ nonlocal and local models of force, respectively. The project focuses on a multiscale analysis of the peridynamics materials model. Peridynamics can be used as a transition between molecular dynamics and classical elasticity so that the difficulties encountered when directly coupling those two models are mitigated. In addition, in some situations, peridynamics can be used all by itself as a material model that accurately and efficiently captures the behavior of materials over a wide range of spatial and temporal scales. Peridynamics is well suited to these purposes because it employs a nonlocal model of force, analogous to that of molecular dynamics; furthermore, at sufficiently large length scales and assuming smooth deformation, peridynamics can be approximated by classical elasticity. The project will extend the emerging mathematical and numerical analysis of peridynamics. One goal is to develop a peridynamics-enabled multiscale multi-model that potentially provides a new and more extensive mathematical basis for coupling classical elasticity and molecular dynamics, thus enabling next generation atomistic-to-continuum multiscale simulations. In addition, a rigorous study of nite element discretizations of peridynamics will be considered. Using the fact that peridynamics is spatially derivative free, we will also characterize the space of admissible peridynamic solutions and carry out systematic analyses of the models, in particular rigorously showing how peridynamics encompasses fracture and other failure phenomena. Additional aspects of the project include the mathematical and numerical analysis of peridynamics applied to stochastic peridynamics models. In summary, the project will make feasible mathematically consistent multiscale models for the analysis and design of advanced materials.
Numerical Model for Conduction-Cooled Current Lead Heat Loads
White, M.J.; Wang, X.L.; /Fermilab; Brueck, H.D.; /DESY
2011-06-10T23:59:59.000Z
Current leads are utilized to deliver electrical power from a room temperature junction mounted on the vacuum vessel to a superconducting magnet located within the vacuum space of a cryostat. There are many types of current leads used at laboratories throughout the world; however, conduction-cooled current leads are often chosen for their simplicity and reliability. Conduction-cooled leads have the advantage of using common materials, have no superconducting/normal state transition, and have no boil-off vapor to collect. This paper presents a numerical model for conduction-cooled current lead heat loads. This model takes into account varying material and fluid thermal properties, varying thicknesses along the length of the lead, heat transfer in the circumferential and longitudinal directions, electrical power dissipation, and the effect of thermal intercepts. The model is validated by comparing the numerical model results to ideal cases where analytical equations are valid. In addition, the XFEL (X-Ray Free Electron Laser) prototype current leads are modeled and compared to the experimental results from testing at DESY's XFEL Magnet Test Stand (XMTS) and Cryomodule Test Bench (CMTB).
On Numerical Considerations for Modeling Reactive Astrophysical Shocks
Papatheodore, Thomas L [ORNL] [ORNL; Messer, Bronson [ORNL] [ORNL
2014-01-01T23:59:59.000Z
Simulating detonations in astrophysical environments is often complicated by numerical approximations to shock structure. A common prescription to ensure correct detonation speeds (and associated quantities) is to prohibit burning inside the numerically broadened shock (Fryxell et al. 1989). We have performed a series of simulations to verify the efficacy of this approximation and to understand how resolution and dimensionality might affect its use. Our results show that, in one dimension, prohibiting burning in the shock is important wherever the carbon burning length is not resolved, in keeping with the results of Fryxell et al. (1989). In two dimensions, we find that the prohibition of shock burning effectively inhibits the development of cellular structure for all but the most highly-resolved cases. We discuss the possible impacts this outcome may have on sub-grid models and detonation propagation in Type Ia supernovae.
Numerical Modeling At Coso Geothermal Area (1995) | Open Energy Information
AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on DeliciousPlasmaPLawrence County,1980) | OpenAl., 2001) |NewNorthNorwalk,Novomer Jump to:NuPowerNumerical Modeling At
Numerical Modeling At Dixie Valley Geothermal Area (Benoit, 1999) | Open
AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on DeliciousPlasmaPLawrence County,1980) | OpenAl., 2001) |NewNorthNorwalk,Novomer Jump to:NuPowerNumerical Modeling
Numerical Modeling At Dixie Valley Geothermal Area (Iovenitti, Et Al.,
AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on DeliciousPlasmaPLawrence County,1980) | OpenAl., 2001) |NewNorthNorwalk,Novomer Jump to:NuPowerNumerical Modeling2013)
Numerical Modeling At Coso Geothermal Area (2000) | Open Energy Information
AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on DeliciousPlasmaP a gHigh Plains WindInformationNVN-079666New look toNotes from the CallNuclearNumerical Modeling At
Numerical Modeling At Coso Geothermal Area (2006) | Open Energy Information
AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on DeliciousPlasmaP a gHigh Plains WindInformationNVN-079666New look toNotes from the CallNuclearNumerical Modeling
Pearce, Fred
2003-01-01T23:59:59.000Z
We use a 3-D finite difference numerical model to generate synthetic seismograms from a simple fractured reservoir
A numerical model of aerosol scavenging: Part 1, Microphysics parameterization
Molenkamp, C.R.; Bradley, M.M.
1991-09-01T23:59:59.000Z
We have developed a three-dimensional numerical model (OCTET) to simulate the dynamics and microphysics of clouds and the transport, diffusion and precipitation scavenging of aerosol particles. In this paper we describe the cloud microphysics and scavenging parameterizations. The representation of cloud microphysics is a bulk- water parameterization which includes water vapor and five types of hydrometeors (cloud droplets, rain drops, ice crystals, snow, and graupel). A parallel parameterization represents the scavenging interactions between pollutant particles and hydrometeors including collection of particles because of condensation nucleation, Brownian and phoretic attachment, and inertial capture, resuspension because of evaporation and sublimation; and transfer interactions where particles collected by one type of hydrometeor are transferred to another type of freezing, melting, accretion, riming and autoconversion.
Effective viscosity of active suspensions: Three-dimensional numerical modeling
Levan Jibuti; Walter Zimmermann; Salima Rafaï; Philippe Peyla
2014-12-10T23:59:59.000Z
A three-dimensional model is proposed for Chlamydomonas Reinhardtii swimming with a breaststroke-like beating of its two flagella. The model reveals unusual angular orbits of the active swimmer under a linear shear flow. Namely, the swimmer sustains orientation transiently across the flow when flagella plane is perpendicular to the shear plane, and amplify the shear-induced rotation along the flow. Such behavior is a result of the interplay between shear-induced deformation and swimmer's periodic beating motion that exerts internal torques on the torque-free swimmer. This particular behavior has some significant consequences on the rheological properties of the suspension that tends to confirm previous experimental results [Phys. Rev. Lett. 104, 098102 (2010)]. We calculated the intrinsic viscosity of the suspension with such isolated modeled microswimmers (dilute case) in shear flow using numerical simulations based on Rotne-Prager approximation. The results show an increased intrinsic viscosity for active swimmer suspensions in comparison to non-active ones in accordance with previous experimental measurements. A major enhancement of the active swimmer viscosity occurs due to the effectively extended shape of the deformable swimming cells. We also recover the experimentally observed shear thinning behavior.
Numerical Methods for the Bogoliubov-Tolmachev-Shirkov model in superconductivity theory
Zhihao Ge; Ruihua Li
2015-03-08T23:59:59.000Z
In the work, the numerical methods are designed for the Bogoliubov-Tolmachev-Shirkov model in superconductivity theory. The numerical methods are novel and effective to determine the critical transition temperature and approximate to the energy gap function of the above model. Finally, a numerical example confirming the theoretical results is presented.
Surussavadee, Chinnawat
2007-01-01T23:59:59.000Z
This thesis develops and validates the MM5/TBSCAT/F([lambda]) model, composed of a mesoscale numerical weather prediction (NWP) model (MM5), a two-stream radiative transfer model (TBSCAT), and electromagnetic models for ...
AI-Based Simulation: An Alternative to Numerical Simulation and Modeling
Mohaghegh, Shahab
: Numerical Modeling, Simulation, Artificial Intelligence, Data Min- ing, Reservoir Modeling, Reservoir data for brown fields. The run-time of AI-Based reservoir models that provide complete field responses Computational Fluid Dynamics (CFD) to Numer- ical Reservoir Simulation (NRS) most of the computational modeling
Numerical study of energy diffusion in King models
Tom Theuns
1995-11-07T23:59:59.000Z
The energy diffusion coefficients D_n(E) (n=1,2) for a system of equal mass particles moving self-consistently in an N-body realisation of a King model are computed from the probability per unit time, P(E, Delta E), that a star with initial energy E will undergo an energy change Delta E. In turn, P is computed from the number of times during the simulation that a particle in a state of given energy undergoes a transition to another state. These particle states are defined directly from the time evolution of E by identifying them with the event occuring between two local maxima in the E(t) curve. If one assumes next that energy changes are uncorrelated between different states, one can use diffusion theory to compute D_n(E). The simulations employ N=512, 2048,... , 32768 particles and are performed using an implementation of Aarseth's direct integrator N-body1 on a massively parallel computer. The more than seven million transitions measured in the largest N simulation provide excellent statistics. The numerically determined D(E)'s are compared against their theoretical counterparts which are computed from phase-space averaged rates of energy change due to independent binary encounters. The overall agreement between them is impressive over most of the energy range, notwithstanding the very different type of approximations involved, giving considerable support to the valid usage of these theoretical expressions to simulate dynamical evolution in Fokker-Planck type calculations.
Numerical Modeling of Nonlinear Coupling between Lines/Beams with Multiple Floating Bodies
Yang, Chan K.
2010-07-14T23:59:59.000Z
STUDY 1: HYDRO-PNEUMATIC TENSIONER MODEL FOR SPAR GLOBAL PERFORMANCE ANALYSIS ................................................102 4.1 Introduction... ....................................................................................................102 4.2 Principal Dimensions of Spar System............................................................102 4.3 Numerical Model............................................................................................104 4.4 Validation of the Model...
An efficient numerical terrestrial scheme (ENTS) for fast earth system modelling
Williamson, Mark
An efficient numerical terrestrial scheme (ENTS) for fast earth system modelling Mark Williamson Working Paper 83 #12;An efficient numerical terrestrial scheme (ENTS) for fast earth system modelling Mark for long time period simulations and large ensemble studies in Earth system models of intermediate
Numerical modeling of fluid flow and time-lapse seismograms ...
gabriela
1. Inst. del Gas y del Petr´oleo - FI, Univ. de Buenos Aires, ARGENTINA. 2 .... 1 + i??s ). (9) where ?e > ?s are relaxation times and Mr = KG,µm. Numerical ...
Numerical modelling of hyperbolic conservation laws using bicharacteristics
Hanke-Bourgeois, Martin
volume methods -1- #12;Overview I. Hyperbolic Conservation Laws Theory of bicharacteristics and evolution, bicharacteritsics stability, accuracy, error analysis III. Numerical Experiments: Wave equation system, Euler eqs finite volume methods -2- #12;· airflow · hydraulic schock · meteorological flow Source: efluid
Numerical Modeling At Coso Geothermal Area (2010) | Open Energy Information
AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on DeliciousPlasmaP a gHigh Plains WindInformationNVN-079666New look toNotes from the CallNuclearNumericalNumerical
Numerical analysis of a one-dimensional elastodynamic model of dry friction and unilateral contact
Renard, Yves - Pôle de Mathématiques, Institut National des Sciences Appliquées de Lyon
.e. hyperbolic) model with dry friction. Since we consider a Coulomb friction law with a slip velocity dependentNumerical analysis of a one-dimensional elastodynamic model of dry friction and unilateral contact in the numerical analysis of more elaborated dynamic purely elastic problems with dry friction. Ó 2001 Elsevier
Penetration of solar radiation in the upper ocean: A numerical model for oceanic and coastal waters
Lee, Zhongping
Penetration of solar radiation in the upper ocean: A numerical model for oceanic and coastal waters (2005), Penetration of solar radiation in the upper ocean: A numerical model for oceanic and coastal; Siegel et al., 1995] have demonstrated that the penetration of EVIS in the upper layer of the ocean plays
Numerical Modeling of CIGS Solar Cells: Definition of the Baseline and
Sites, James R.
Thesis Numerical Modeling of CIGS Solar Cells: Definition of the Baseline and Explanation our supervision by Markus Gloeckler entitled "Numerical Modeling of CIGS Solar Cells: Definition. A three-layer structure, simulating a Cu(InGa)Se2 (CIGS) heterojunction solar cell, was set up using
NUMERICAL MODELING OF CIGS AND CdTe SOLAR CELLS: SETTING THE BASELINE
Sites, James R.
NUMERICAL MODELING OF CIGS AND CdTe SOLAR CELLS: SETTING THE BASELINE M. Gloeckler, A important complications that are often found in experimental CIGS and CdTe solar cells. 1. INTRODUCTION Numerical modeling of polycrystalline thin-film solar cells is an important strategy to test the viability
Wells, Scott A.
792 / JOURNAL OF ENVIRONMENTAL ENGINEERING / SEPTEMBER 1999 NUMERICAL MODEL OF SEDIMENTATION publication on July 20, 1998. This paper is part of the Journal of Environmental Engineering, Vol. 125, No. 9 would enable engineers to optimize the process design. A numerical model of the physics of gravity
NUMERICAL SOLUTION OF RESERVOIR FLOW MODELS BASED ON LARGE TIME STEP OPERATOR SPLITTING ALGORITHMS
NUMERICAL SOLUTION OF RESERVOIR FLOW MODELS BASED ON LARGE TIME STEP OPERATOR SPLITTING ALGORITHMS and analysis of large time step operator splitting algorithms for the numerical simulation of multiphase flow the main ideas behind these novel operator splitting algorithms for a basic twophase flow model. Special
Numerical Model of a Tensioner System and Flex Joint
Huang, Han
2013-07-27T23:59:59.000Z
Top Tensioned Riser (TTR) and Steel Catenary Riser (SCR) are often used in a floating oil/gas production system deployed in deep water for oil transport. This study focuses on the improvements to the existing numerical code, known as CABLE3D...
Numerical Modeling of Human Effect on Indoor Propagation
Paris-Sud XI, Université de
several types of obstacles such a brick enclosure walls, one metal heater, two metal wardrobes, two metal antenna placed 85 cm above the floor level. The spatial step is /10, or is the wavelength, chosen by the numerical dispersions. The choice of the spatial step is a compromise between the minimization of inaccuracy
Numerical modeling of short pulse laser interaction with Au nanoparticle surrounded by water
Zhigilei, Leonid V.
Numerical modeling of short pulse laser interaction with Au nanoparticle surrounded by water Alexey, University of Virginia, USA Available online 3 February 2007 Abstract Short pulse laser interaction modeling; Nanoparticles; Cell targeting; Laser damage 1. Introduction Short pulse laser irradiation
Thermoelectric Generators for Automotive Waste Heat Recovery Systems Part I: Numerical Modeling
Xu, Xianfan
Thermoelectric Generators for Automotive Waste Heat Recovery Systems Part I: Numerical Modeling (TEG) designed for automotive waste heat recovery systems. This model is capable of computing telluride TEMs. Key words: Thermoelectric generators, waste heat recovery, automotive exhaust, skutterudites
Bardhan, Jaydeep Porter, 1978-
2006-01-01T23:59:59.000Z
This thesis presents a set of numerical techniques that extend and improve computational modeling approaches for biomolecule analysis and design. The presented research focuses on surface formulations of modeling problems ...
Numerical modelling of a radio-frequency micro ion thruster
Tsay, Michael Meng-Tsuan
2006-01-01T23:59:59.000Z
A simple performance model is developed for an inductively-coupled radio-frequency micro ion thruster. Methods of particle and energy balance are utilized for modeling the chamber plasma discharge. A transformer model is ...
Fractional Calculus in Hydrologic Modeling: A Numerical Perspective
David A. Benson; Mark M. Meerschaert; Jordan Revielle
2012-01-01T23:59:59.000Z
Fractional derivatives can be viewed either as a handy extension of classical calculus or, more deeply, as mathematical operators defined by natural phenomena. This follows the view that the diffusion equation is defined as the governing equation of a Brownian motion. In this paper, we emphasize that fractional derivatives come from the governing equations of stable Levy motion, and that fractional integration is the corresponding inverse operator. Fractional integration, and its multi-dimensional extensions derived in this way, are intimately tied to fractional Brownian (and Levy) motions and noises. By following these general principles, we discuss the Eulerian and Lagrangian numerical solutions to fractional partial differential equations, and Eulerian methods for stochastic integrals. These numerical approximations illuminate the essential nature of the fractional calculus.
Analysis and verification of a three-dimensional hydrodynamic numerical model
Abraham, David Daniel
1991-01-01T23:59:59.000Z
of Advisory Committee: Dr. W. P. James A three-dimensional hydrodynamic numerical model (RMA10) was analyized and verified. The verification process involved the comparison of model and measured water surface elevations, velocities, and depth averaged... velocities. The verification (measured) data were obtained from a large (200 cfs) flume. The numerical grid was set up such that the spacial location of the computational points was as close as practical to that of the measured data. The model computed...
Numerical Modeling At Raft River Geothermal Area (1983) | Open Energy
AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on DeliciousPlasmaPLawrence County,1980) | OpenAl., 2001) |NewNorthNorwalk,Novomer Jump to:NuPowerNumerical1984)
Numerical Modeling of Transient Basin and Range Extensional Geothermal
AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on DeliciousPlasmaPLawrence County,1980) | OpenAl., 2001) |NewNorthNorwalk,Novomer Jump to:NuPowerNumerical1984)Systems |
Numerical Modelling of Geothermal Systems a Short Introduction | Open
AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on DeliciousPlasmaPLawrence County,1980) | OpenAl., 2001) |NewNorthNorwalk,Novomer Jump to:NuPowerNumerical1984)Systems
Numerical Modeling At Coso Geothermal Area (2007) | Open Energy Information
AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on DeliciousPlasmaP a gHigh Plains WindInformationNVN-079666New look toNotes from the CallNuclearNumerical
Turbulent Supersonic Channel Flow: Direct Numerical Simulation and Modeling
Heinz, Stefan
modeling: the turbulence frequency production mechanism, wall damping effects on turbulence model frequency production mechanisms and wall damping effects may be explained very well on the basis, Chik w = wall viscosity = kinematic viscosity, = T = turbulent kinematic viscosity, Ck d = pressure
Numerically Efficient Water Quality Modeling and Security Applications
Mann, Angelica
2013-02-04T23:59:59.000Z
. . . . . . . . . . . . . . . . . . . . . . . . . . . . 13 1.4 Summary of contributions . . . . . . . . . . . . . . . . . . . . . . . . 17 2 THE MERLION WATER QUALITY MODEL . . . . . . . . . . . . . . . . 20 2.1 Model development . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20 2.2 Origin... tracking algorithm . . . . . . . . . . . . . . . . . . . . . . . . 24 2.3 Merlion water quality model . . . . . . . . . . . . . . . . . . . . . . . 31 2.4 Summary . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34 3 EFFICIENT MULTI...
Bürger, Raimund
-dimensional model of sedimentation of suspensions of small solid particles dispersed in a viscous fluid. This model accepted spatially one-dimensional sedimentation model [35] gives rise to one scalar, nonlinear hyperbolicINTERNATIONAL JOURNAL OF c 2011 Institute for Scientific NUMERICAL ANALYSIS AND MODELING Computing
Bürger, Raimund
-dimensional model of sedimentation of suspensions of small solid particles dispersed in a viscous fluid. This model accepted spatially one-dimensional sedimentation model [35] gives rise to one scalar, nonlinear hyperbolicINTERNATIONAL JOURNAL OF c 2012 Institute for Scientific NUMERICAL ANALYSIS AND MODELING Computing
NUMERICAL MODELING OF COASTAL TSUNAMI DISSIPATION AND Stphan T. Grilli1
Grilli, Stéphan T.
1 NUMERICAL MODELING OF COASTAL TSUNAMI DISSIPATION AND IMPACT Stéphan T. Grilli1 , Jeffrey C Tehranirad2 Recent observations of the coastal impact of large tsunamis (e.g., Indian Ocean 2004; Tohoku 2011) and related numerical and theoretical works have made it increasingly clear that tsunami waves arrive
NUMERICAL MODELING OF COASTAL TSUNAMI IMPACT DISSIPATION Stphan T. Grilli1
Kirby, James T.
1 NUMERICAL MODELING OF COASTAL TSUNAMI IMPACT DISSIPATION AND IMPACT Stéphan T. Grilli1 , Jeffrey Tehranirad2 Recent observations of the coastal impact of large tsunamis (e.g., Indian Ocean 2004; Tohoku 2011) and related numerical and theoretical works have made it increasingly clear that tsunami waves arrive
Numerical Simulation of Pulse-Tube Refrigerators: 1D model I.A. Lyulina1
Eindhoven, Technische Universiteit
of a piston, an aftercooler (AC), a regenerator, a cold heat exchanger (CHX), a tube, a hot heat exchanger numerical model has been introduced to study steady oscillatory heat and mass transfer in the tube section, numerical simulation, high resolution scheme 1 Introduction The pulse tube is a relatively new type
USING LEARNING MACHINES TO CREATE SOLAR RADIATION MAPS FROM NUMERICAL WEATHER PREDICTION MODELS,
Paris-Sud XI, UniversitÃ© de
USING LEARNING MACHINES TO CREATE SOLAR RADIATION MAPS FROM NUMERICAL WEATHER PREDICTION MODELS to develop a methodology to generate solar radiation maps using information from different sources. First with conclusions and next works in the last section. Keywords: Solar Radiation maps, Numerical Weather Predictions
Numerical models of caldera deformation: Effects of multiphase and multicomponent hydrothermal studies addressing the effects of multiphase flow on crustal mechanics have been attempted. Recent numerical simulations of multiphase (liquid-gas), multicomponent (H2OÂCO2) hydrothermal fluid flow
An Updated Numerical Model Of The Larderello-Travale Geothermal...
been used to conduct a regional modeling study to investigate the production mechanism of superheated steam, the interactions between the geothermal field and the surrounding deep...
Numerical Modeling At Dixie Valley Geothermal Area (Iovenitti...
inversion of gravity data, and thermal conductivity contrasts between basin-filling sediments and range-basement rocks. The three dimensional conductive steady state model...
Numerically Estimating Internal Models of Dynamic Virtual Objects
Sekuler, Robert
human subjects to manipulate a computer-animated virtual object. This virtual object (vO) was a high, human cognition, human information processing, ideal performer, internal model, virtual object, virtual, specifically how humans acquire an internal model of a dynamic virtual object. Our methodology minimizes
Numerical model for steel catenary riser on seafloor support
You, Jung Hwan
2007-04-25T23:59:59.000Z
of development of a simplified seafloor support model. This model simulates the seafloor-pipe interaction as a flexible pipe supported on a bed of springs. Constants for the soil springs were derived from finite element studies performed in a separate, parallel...
Ravago Bastardo, Delmira Cristina
2005-08-29T23:59:59.000Z
The main objective of this research is to compare the performance of cyclic steam injection using horizontal wells based on the analytical model developed by Gunadi against that based on numerical simulation. For comparison, ...
Numerical modeling of elastic wave scattering by near-surface heterogeneities
Al Muhaidib, Abdulaziz
2013-01-01T23:59:59.000Z
A perturbation method for elastic waves and numerical forward modeling are used to calculate the effects of seismic wave scattering from arbitrary shape shallow subsurface heterogeneities. Wave propagation is simulated ...
Numerical modeling of species transport in turbulent flow and experimental study on aerosol sampling
Vijayaraghavan, Vishnu Karthik
2007-04-25T23:59:59.000Z
Numerical simulations were performed to study the turbulent mixing of a scalar species in straight tube, single and double elbow flow configurations. Different Reynolds Averaged Navier Stokes (RANS) and Large Eddy Simulation (LES) models were used...
Mass and charge flow in nanopores: numerical simulation via mesoscale models
Cecconi, Fabio
Mass and charge flow in nanopores: numerical simulation via mesoscale models Mauro Chinappi1 at nanoscale is here addressed via a recent developed mesoscale approach. In particular the flow
2D-Modelling of pellet injection in the poloidal plane: results of numerical tests
Paris-Sud XI, Université de
2D-Modelling of pellet injection in the poloidal plane: results of numerical tests P. Lalousis developed for computing the expansion of pellet-produced clouds in the poloidal plane. The expansion
Seismic Scattering Attributes to Estimate Reservoir Fracture Density: A Numerical Modeling Study
Pearce, Frederick Douglas
We use a 3-D finite difference numerical model to generate synthetic seismograms from a simple fractured reservoir containing evenly-spaced, discrete, vertical fracture zones. The fracture zones are represented using a ...
P. Douillet S. Ouillon E. Cordier A numerical model for ne suspended sediment transport
transport á Numerical model á Deposition á Erosion á New Caledonia Introduction In New Caledonia, where open-cast for cohesive particles (e.g., Nicholson and O'Connor 1986; Teisson 1991; Brenon and Le Hir 1999), around sand
Seismic scattering attributes to estimate reservoir fracture density : a numerical modeling study
Pearce, Frederick D. (Frederick Douglas), 1978-
2003-01-01T23:59:59.000Z
We use a 3-D finite difference numerical model to generate synthetic seismograms from a simple fractured reservoir containing evenly-spaced, discrete, vertical fracture zones. The fracture zones are represented using a ...
Development and validation of a vertically two-dimensional mesoscale numerical model
Walters, Michael Kent
1985-01-01T23:59:59.000Z
DEVELOPMENT AND VALIDATION OF A VERTICALLY TWO-DIMENSIONAL MESOSCALE NUMERICAL MODEL A Thesis by MICHAEL KENT WALTERS Submitted to the Graduate College of Texas AsM University in partial fulfillment of the requirement for the degree of MASTER... OF SCIENCE August 1985 Major Subject: Meteorology DEVELOPMENT AND VALIDATION OF A VERTICALLY TWO-DIMENSIONAL MESOSCALE NUMERICAL MODEL A Thesis by MICHAEL KENT WALTERS Approved as to style and content by: Dusan Djuric (Chair of Committee) WP...
Moist processes and the quasi-hydrostatic approximation in a mesoscale numerical model
Kennedy, Charles Joseph
1987-01-01T23:59:59.000Z
MOIST PROCESSES AND THE QUASI-HYDROSTATIC APPROXIMATION IN A MESOSCALE NUMERICAL MODEL A Thesis by CHARLES JOSEPH KENNEDY Submitted to the Graduate College of Texas A&M University in partial fulfillment of the requirement for the degree... of MASTER OF SCIENCE December 1987 Major Subject: Meteorology MOIST PROCESSES AND THE QUASI-HYDROSTATIC APPROXIMATION IN A MESOSCALE NUMERICAL MODEL A Thesis by CHARLES JOSEPH KENNEDY Approved as to style and content by: Dusan Djuric (Chair...
Giovanni Noselli; Antonio DeSimone
2014-08-26T23:59:59.000Z
We present experimental and numerical results for a model crawler which is able to extract net positional changes from reciprocal shape changes, i.e. 'breathing-like' deformations, thanks to directional, frictional interactions with a textured solid substrate, mediated by flexible inclined feet. We also present a simple reduced model that captures the essential features of the kinematics and energetics of the gait, and compare its predictions with the results from experiments and from numerical simulations.
NUMERICAL MODELING OF TURBULENT FLOW IN A COMBUSTION TUNNEL
Ghoniem, A.F.
2013-01-01T23:59:59.000Z
1VJcDona·ld, H. (1979) Combustion r 1 iodeJ·ing in Two and1979) Practical Turbulent-Combustion Interaction Models forInternation on Combustors. Combustion The 17th Symposium
Numerical Modeling of Seafloor Interation with Steel Catenary Riser
You, Jung Hwan
2012-10-19T23:59:59.000Z
degradation it is possible to simulate the trench formation process and estimate deflections and moments along the riser length. The seabed model is used to perform parametric studies to assess the effects of stiffness, soil strength, amplitude of pipe...
Ductile fracture modeling : theory, experimental investigation and numerical verification
Xue, Liang, 1973-
2007-01-01T23:59:59.000Z
The fracture initiation in ductile materials is governed by the damaging process along the plastic loading path. A new damage plasticity model for ductile fracture is proposed. Experimental results show that fracture ...
Numerical models of pressure pulse generation by imploding metal liners
Humphries, S. Jr. [Acceleration Associates, Albuquerque, NM (United States)] [Acceleration Associates, Albuquerque, NM (United States); [Univ. of New Mexico, Albuquerque, NM (United States). Dept. of Electrical Engineering; Ekdahl, C.A. [Los Alamos National Lab., NM (United States)] [Los Alamos National Lab., NM (United States)
1996-12-01T23:59:59.000Z
The authors describe numerical calculations of pressure pulse generation using imploding liners. Liners are metal cylinders that are magnetically compressed by an intense axial current flow from a high-power pulse generator. The simulations cover the acceleration of the liner, collision with an internal diagnostic target, followed by compression and shock wave heating of the target. With the projected current waveform of the Atlas capacitor bank (in development at Los Alamos National Laboratory), initial results suggest that it may be possible to achieve pressures exceeding 3,000 Gpa (30 Mbar) in a 4 mm diameter sample over an interval of 100--200 ns. The simulations were carried out with Crunch, a new one-dimensional hydrodynamics package for advanced personal computers. The program uses finite-element techniques to solve the coupled problems of hydrodynamics and magnetic diffusion. Crunch fully supports loading and interpolating Sesame equation-of-state tables. The program exhibits excellent stability, even for collisions between material shells and shock convergence on axis. In contrast to previous work, the present studies follow the full process through solid target collision and compression. The work supports the High-Energy Density Physics Program of Los Alamos National Laboratory (LANL), a component of the US Department of Energy Stockpile Stewardship Program. The purpose of this program is maintenance of the nuclear weapons stockpile through improved computational ability and above-ground experiments. Imploding liners driven by conventional capacitor banks constitute a portion of the program to study matter at high pressure.
IMPROVED NUMERICAL METHODS FOR MODELING RIVER-AQUIFER INTERACTION.
Tidwell, Vincent C.; Sue Tillery; Phillip King
2008-09-01T23:59:59.000Z
A new option for Local Time-Stepping (LTS) was developed to use in conjunction with the multiple-refined-area grid capability of the U.S. Geological Survey's (USGS) groundwater modeling program, MODFLOW-LGR (MF-LGR). The LTS option allows each local, refined-area grid to simulate multiple stress periods within each stress period of a coarser, regional grid. This option is an alternative to the current method of MF-LGR whereby the refined grids are required to have the same stress period and time-step structure as the coarse grid. The MF-LGR method for simulating multiple-refined grids essentially defines each grid as a complete model, then for each coarse grid time-step, iteratively runs each model until the head and flux changes at the interfacing boundaries of the models are less than some specified tolerances. Use of the LTS option is illustrated in two hypothetical test cases consisting of a dual well pumping system and a hydraulically connected stream-aquifer system, and one field application. Each of the hypothetical test cases was simulated with multiple scenarios including an LTS scenario, which combined a monthly stress period for a coarse grid model with a daily stress period for a refined grid model. The other scenarios simulated various combinations of grid spacing and temporal refinement using standard MODFLOW model constructs. The field application simulated an irrigated corridor along the Lower Rio Grande River in New Mexico, with refinement of a small agricultural area in the irrigated corridor.The results from the LTS scenarios for the hypothetical test cases closely replicated the results from the true scenarios in the refined areas of interest. The head errors of the LTS scenarios were much smaller than from the other scenarios in relation to the true solution, and the run times for the LTS models were three to six times faster than the true models for the dual well and stream-aquifer test cases, respectively. The results of the field application show that better estimates of daily stream leakage can be made with the LTS simulation, thereby improving the efficiency of daily operations for an agricultural irrigation system. ACKNOWLEDGEMENTSThe authors appreciatively acknowledge support for Sue Tillery provided by Sandia National Laboratories' through a Campus Executive Laboratory Directed Research and Development (LDRD) research project.Funding for this study was provided by Directed Research and Development (LDRD) research project.
Efficient Dynamic Modeling, Numerical Optimal Control and Experimental Results for Various Gaits
Stryk, Oskar von
Efficient Dynamic Modeling, Numerical Optimal Control and Experimental Results for Various Gaits. A fully three- dimensional dynamical model of Sony's four-legged robot is used to state an optimal control model and the algorithm for evaluating the dynamics. The formulation of the optimal control problem
Numerical Modelling of a Pulse Combustion Burner: Limiting Conditions of Stable
Vuik, Kees
Numerical Modelling of a Pulse Combustion Burner: Limiting Conditions of Stable Operation P.A. van a mathematical analysis of a simple model for thermal pulse combustion and determines conditions under which this model can describe stable pulse operation. 1 Introduction Compared to conventional combustion, pulse
A CONSISTENT MODELLING METHODOLOGY FOR SECONDARY1 SETTLING TANKS: A RELIABLE NUMERICAL METHOD2
Bürger, Raimund
relations for hindered settling, compression and dispersion can be used within the model, allowing the user, continuous sedimentation, secondary clarifier, simulation5 model, partial differential equation6 NomenclatureA CONSISTENT MODELLING METHODOLOGY FOR SECONDARY1 SETTLING TANKS: A RELIABLE NUMERICAL METHOD2
Numerical Solution of a Free-Boundary Problem for Percussive Deep Drilling Modeling by BEM
Mikhailov, Sergey
Numerical Solution of a Free-Boundary Problem for Percussive Deep Drilling Modeling by BEM S to a stationary-periodic quasi-static model of rock percus- sive deep drilling is presented. The rock is modeled. An auxiliary problem of stationary inden- tation of a rigid drill bit is considered first, where it is assumed
Numerical Experiments of Some Krylov Subspace Methods for Black Oil Model
Lai, Choi-Hong
Numerical Experiments of Some Krylov Subspace Methods for Black Oil Model Jianwen Cao #3; Choi of linear systems originated from the black oil model in oil reservoir simulation. There exists some Krylov subspace algorithms and pre- conditioning techniques for the black oil model as appeared in the literature
Xing, Lu [Oklahoma State University; Cullin, James [Oklahoma State University; Spitler, Jeffery [Oklahoma State University; Im, Piljae [ORNL; Fisher, Daniel [Oklahoma State University
2011-01-01T23:59:59.000Z
A new type of ground heat exchanger that utilizes the excavation often made for basements or foundations has been proposed as an alternative to conventional ground heat exchangers. This article describes a numerical model that can be used to size these foundation heat exchanger (FHX) systems. The numerical model is a two-dimensional finite-volume model that considers a wide variety of factors, such as soil freezing and evapotranspiration. The FHX numerical model is validated with one year of experimental data collected at an experimental house located near Oak Ridge, Tennessee. The model shows good agreement with the experimental data-heat pump entering fluid temperatures typically within 1 C (1.8 F) - with minor discrepancies due to approximations, such as constant moisture content throughout the year, uniform evapotranspiration over the seasons, and lack of ground shading in the model.
Numerical Integration Numerical Summation
Cohen, Henri
Numerical Integration Numerical Summation Numerical Extrapolation Numerical Recipes for Multiprecision Computations #12;Numerical Integration Numerical Summation Numerical Extrapolation Multiprecision, integration, summation, extrapolation, evaluation of continued fractions, Euler products and sums, complete
Lattice Percolation Approach to Numerical Modeling of Tissue Aging
Privman, Vladimir; Libert, Sergiy
2015-01-01T23:59:59.000Z
We describe a percolation-type approach to modeling of the processes of aging and certain other properties of tissues analyzed as systems consisting of interacting cells. Tissues are considered as structures made of regular healthy, senescent, dead (apoptotic) cells, and studied dynamically, with the ongoing processes including regular cell division to fill vacant sites left by dead cells, healthy cells becoming senescent or dying, and other processes. Statistical-mechanics description can provide patterns of time dependence and snapshots of morphological system properties. An illustrative application of the developed theoretical modeling approach is reported, confirming recent experimental findings that inhibition of senescence can lead to extended lifespan.
Small Buoys for Energy Harvesting : Experimental and Numerical Modeling Studies
Grilli, Stéphan T.
-contained (water tight) resonating multiple-spar buoy (or Starspar), in which a longer central spar houses the LEG and is surrounded by shallower, satellite spars, providing both form stability and a reduced overall average draft peak spectral period. For spar buoys, the former is simply controlled by buoy draft. Scale model
THEACCURACY OF NUMERICAL MODELS FOR CONTINUUM STANLY STEINBERG
Steinberg, Stanly
Water Rain Evaporation Figure 1.1. A Lake These notes are directed towards students who have some modest of water in the lake. Water can enter and leave the lake by various means: rain; evaporation; rivers and how much error is there is solving the mathematical model. #12; 2 Lake Ground Ground Water Ground
Numerical Modeling of Seafloor Interation with Steel Catenary Riser
You, Jung Hwan
2012-10-19T23:59:59.000Z
, and moment. Primary outputs from this model include the deflected shape of the riser pipe and bending moments along riser length. The code also provides the location of maximum trench depth and the position where the maximum bending moment occurs and any...
EFFICIENT NUMERICAL SOLUTION TECHNIQUES IN COMPOSITION MODEL 1
for the study of enhanced oil recovery techniques which lower the surface tension between phase and viscosity the compositional effects and mass transfer, must be used to accurately simulate enhanced oil recovery processes Ewing, 2 and Magne Espedal 4 Abstract. A generalized isothermal multiphase model for enhanced oil
Numerical modeling of an all vanadium redox flow battery.
Clausen, Jonathan R.; Brunini, Victor E.; Moffat, Harry K.; Martinez, Mario J.
2014-01-01T23:59:59.000Z
We develop a capability to simulate reduction-oxidation (redox) flow batteries in the Sierra Multi-Mechanics code base. Specifically, we focus on all-vanadium redox flow batteries; however, the capability is general in implementation and could be adopted to other chemistries. The electrochemical and porous flow models follow those developed in the recent publication by [28]. We review the model implemented in this work and its assumptions, and we show several verification cases including a binary electrolyte, and a battery half-cell. Then, we compare our model implementation with the experimental results shown in [28], with good agreement seen. Next, a sensitivity study is conducted for the major model parameters, which is beneficial in targeting specific features of the redox flow cell for improvement. Lastly, we simulate a three-dimensional version of the flow cell to determine the impact of plenum channels on the performance of the cell. Such channels are frequently seen in experimental designs where the current collector plates are borrowed from fuel cell designs. These designs use a serpentine channel etched into a solid collector plate.
Joint physical and numerical modeling of water distribution networks.
Zimmerman, Adam; O'Hern, Timothy John; Orear, Leslie Jr.; Kajder, Karen C.; Webb, Stephen Walter; Cappelle, Malynda A.; Khalsa, Siri Sahib; Wright, Jerome L.; Sun, Amy Cha-Tien; Chwirka, J. Benjamin; Hartenberger, Joel David; McKenna, Sean Andrew; van Bloemen Waanders, Bart Gustaaf; McGrath, Lucas K.; Ho, Clifford Kuofei
2009-01-01T23:59:59.000Z
This report summarizes the experimental and modeling effort undertaken to understand solute mixing in a water distribution network conducted during the last year of a 3-year project. The experimental effort involves measurement of extent of mixing within different configurations of pipe networks, measurement of dynamic mixing in a single mixing tank, and measurement of dynamic solute mixing in a combined network-tank configuration. High resolution analysis of turbulence mixing is carried out via high speed photography as well as 3D finite-volume based Large Eddy Simulation turbulence models. Macroscopic mixing rules based on flow momentum balance are also explored, and in some cases, implemented in EPANET. A new version EPANET code was developed to yield better mixing predictions. The impact of a storage tank on pipe mixing in a combined pipe-tank network during diurnal fill-and-drain cycles is assessed. Preliminary comparison between dynamic pilot data and EPANET-BAM is also reported.
ARRAY OPTIMIZATION FOR TIDAL ENERGY EXTRACTION IN A TIDAL CHANNEL – A NUMERICAL MODELING ANALYSIS
Yang, Zhaoqing; Wang, Taiping; Copping, Andrea
2014-04-18T23:59:59.000Z
This paper presents an application of a hydrodynamic model to simulate tidal energy extraction in a tidal dominated estuary in the Pacific Northwest coast. A series of numerical experiments were carried out to simulate tidal energy extraction with different turbine array configurations, including location, spacing and array size. Preliminary model results suggest that array optimization for tidal energy extraction in a real-world site is a very complex process that requires consideration of multiple factors. Numerical models can be used effectively to assist turbine siting and array arrangement in a tidal turbine farm for tidal energy extraction.
Modeling and Design of RF MEMS Structures Using Computationally Efficient Numerical Techniques
Tentzeris, Manos
Modeling and Design of RF MEMS Structures Using Computationally Efficient Numerical Techniques N. A Abstract The modeling of MEMS structures using MRTD is presented. Many complex RF structures have been communication systems efficiently and accurately. Specifically, micromachined structures such as MEMS
Numerical modeling of fiber lasers with long and ultra-long ring cavity
Turitsyn, Sergei K.
Numerical modeling of fiber lasers with long and ultra-long ring cavity I.A. Yarutkina,1,2, O. S. Kobtsev, S. Kukarin, and Y. Fedotov, "Ultra-low repetition rate mode-locked fiber laser with high highlight two important aspects related to a mathematical modeling of pulsed fiber lasers with long
A consistent modelling methodology for secondary settling tanks: a reliable numerical method
Bürger, Raimund
accounting for sediment compressibility, and a dispersion term for turbulence. In addition, the solutionA consistent modelling methodology for secondary settling tanks: a reliable numerical method modelling methodology for secondary settling tanks (SSTs) leads to a partial differential equation (PDE
URBAN AEROSOLS SURVEY USING LIDAR AND NUMERICAL MODEL S. GEFFROY1
Boyer, Edmond
URBAN AEROSOLS SURVEY USING LIDAR AND NUMERICAL MODEL S. GEFFROY1 , L. SOULHAC2 , E. FREJAFON3 , R technologique ALATA BP2, F-60550 Verneuil-en-Halatte, France. Keywords: LIDAR, URBAN AEROSOLS, MODEL, IMPACT SURVEY. INTRODUCTION The impact of particulate matters and aerosols on environment and on radiative
Numerical modeling of heat transfer and fluid flow in rotor-stator cavities with throughflow
Boyer, Edmond
Numerical modeling of heat transfer and fluid flow in rotor-stator cavities with throughflow S modeling of the turbulent flow in a rotor-stator cavity subjected to a superimposed throughflow with heat the dynamical effects from the heat transfer process. The fluid flow in an enclosed disk system with axial
NUMERICAL SOLUTION OF RESERVOIR FLOW MODELS BASED ON LARGE TIME STEP OPERATOR SPLITTING ALGORITHMS
processes. A black-oil model is commonly used to describe water injection. This model works well. Special focus is posed on the numerical solution algorithms for the saturation equation, which is a convection dominated, degenerate convection-di#11;usion equation. Both theory and applications are discussed
Gracie, Robert
PeerReview Only An XFEM Model for Carbon Sequestration Journal: International Journal for Numerical method, Carbon Sequestration, Multiphase flow, XFEM, Multifield systems, Petrov-Galerkin httpScience (www.interscience.wiley.com). DOI: 10.1002/nme An XFEM Model for Carbon Sequestration Chris Ladubec
Intercomparison of Single-Column Numerical Models for the Prediction of Radiation Fog
Intercomparison of Single-Column Numerical Models for the Prediction of Radiation Fog THIERRY-term forecasting of fog is a difficult issue that can have a large societal impact. Radiation fog appears layers of the atmosphere. Current NWP models poorly forecast the life cycle of fog, and improved NWP
Nick, F. M.; van der Veen, Cornelis J.; Oerlemans, J.
2007-07-11T23:59:59.000Z
A one-dimensional numerical ice flow model is used to study the advance of a tidewater glacier into deep water. Starting with ice-free conditions, the model simulates glacier growth at higher elevations followed by advance on land to the head...
FOUNDATION, ANALYSIS, AND NUMERICAL INVESTIGATION OF A VARIATIONAL NETWORK-BASED MODEL FOR RUBBER
Boyer, Edmond
FOUNDATION, ANALYSIS, AND NUMERICAL INVESTIGATION OF A VARIATIONAL NETWORK-BASED MODEL FOR RUBBER, many models based on polymer chain statistics have been proposed to describe rubber elasticity. Recently, Alicandro, Cicalese, and the first author rigorously derived a continuum theory of rubber
Effects of winds and Caribbean eddies on the frequency of Loop Current eddy shedding: A numerical of Mexico, Caribbean Current, Loop Current, eddy shedding, winds and eddies, numerical ocean model Citation of Loop Current eddy shedding: A numerical model study, J. Geophys. Res., 108(C10), 3324, doi:10
Incorporation of the capillary hysteresis model HYSTR into the numerical code TOUGH
Niemi, A.; Bodvarsson, G.S.; Pruess, K.
1991-11-01T23:59:59.000Z
As part of the work performed to model flow in the unsaturated zone at Yucca Mountain Nevada, a capillary hysteresis model has been developed. The computer program HYSTR has been developed to compute the hysteretic capillary pressure -- liquid saturation relationship through interpolation of tabulated data. The code can be easily incorporated into any numerical unsaturated flow simulator. A complete description of HYSTR, including a brief summary of the previous hysteresis literature, detailed description of the program, and instructions for its incorporation into a numerical simulator are given in the HYSTR user`s manual (Niemi and Bodvarsson, 1991a). This report describes the incorporation of HYSTR into the numerical code TOUGH (Transport of Unsaturated Groundwater and Heat; Pruess, 1986). The changes made and procedures for the use of TOUGH for hysteresis modeling are documented.
Characterization of Texas lignite and numerical modeling of its in-situ gasification
Wang, Yih-Jy
1983-01-01T23:59:59.000Z
Modeling Site selection for in-situ gasification projects normally involves application of site screen1ng criteria. Some of these cr1teria were discussed by Russell et al. (1983). Numerical simulation may play an important role in s1te selection...CHARACTERIZATION OF TEXAS LIGNITE AND NUMERICAL MODELING OF ITS IN-SITU GASIFICATION A Thesis by YIH-JY WANG Submitted to the Graduate College of Texas A&M University in partial fulfillment of the requirements for the degree of MASTER...
Proceedings of the Numerical Modeling for Underground Nuclear Test Monitoring Symposium
Taylor, S.R.; Kamm, J.R. [eds.
1993-11-01T23:59:59.000Z
The purpose of the meeting was to discuss the state-of-the-art in numerical simulations of nuclear explosion phenomenology with applications to test ban monitoring. We focused on the uniqueness of model fits to data, the measurement and characterization of material response models, advanced modeling techniques, and applications of modeling to monitoring problems. The second goal of the symposium was to establish a dialogue between seismologists and explosion-source code calculators. The meeting was divided into five main sessions: explosion source phenomenology, material response modeling, numerical simulations, the seismic source, and phenomenology from near source to far field. We feel the symposium reached many of its goals. Individual papers submitted at the conference are indexed separately on the data base.
Numerical Investigation of Flapwise-Torsional Vibration Model of a Smart Section Blade with Microtab
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Li, Nailu; Balas, Mark J.; Yang, Hua; Jiang, Wei; Magar, Kaman T.
2015-01-01T23:59:59.000Z
This study presents a method to develop an aeroelastic model of a smart section blade equipped with microtab. The model is suitable for potential passive vibration control study of the blade section in classic flutter. Equations of the model are described by the nondimensional flapwise and torsional vibration modes coupled with the aerodynamic model based on the Theodorsen theory and aerodynamic effects of the microtab based on the wind tunnel experimental data. The aeroelastic model is validated using numerical data available in the literature and then utilized to analyze the microtab control capability on flutter instability case and divergence instabilitymore »case. The effectiveness of the microtab is investigated with the scenarios of different output controllers and actuation deployments for both instability cases. The numerical results show that the microtab can effectively suppress both vibration modes with the appropriate choice of the output feedback controller.« less
Bruneau, Steve
in Pack Ice Roelof C. Dragt Offshore Engineering Faculty of Mechanical, Maritime and Material Engineering of experiments to validate a Graphics Processing Unit based numerical modelling of ship operations in 2D pack ice interaction, 2D Model Experiments, Image Processing. I. INTRODUCTION A ship travelling through pack ice
:,; Abbott, B P; Abbott, R; Abbott, T; Abernathy, M R; Accadia, T; Acernese, F; Ackley, K; Adams, C; Adams, T; Addesso, P; Adhikari, R X; Affeldt, C; Agathos, M; Aggarwal, N; Aguiar, O D; Ain, A; Ajith, P; Alemic, A; Allen, B; Allocca, A; Amariutei, D; Andersen, M; Anderson, R; Anderson, S B; Anderson, W G; Arai, K; Araya, M C; Arceneaux, C; Areeda, J; Aston, S M; Astone, P; Aufmuth, P; Aulbert, C; Austin, L; Aylott, B E; Babak, S; Baker, P T; Ballardin, G; Ballmer, S W; Barayoga, J C; Barbet, M; Barish, B C; Barker, D; Barone, F; Barr, B; Barsotti, L; Barsuglia, M; Barton, M A; Bartos, I; Bassiri, R; Basti, A; Batch, J C; Bauchrowitz, J; Bauer, Th S; Behnke, B; Bejger, M; Beker, M G; Belczynski, C; Bell, A S; Bell, C; Bergmann, G; Bersanetti, D; Bertolini, A; Betzwieser, J; Beyersdorf, P T; Bilenko, I A; Billingsley, G; Birch, J; Biscans, S; Bitossi, M; Bizouard, M A; Black, E; Blackburn, J K; Blackburn, L; Blair, D; Bloemen, S; Blom, M; Bock, O; Bodiya, T P; Boer, M; Bogaert, G; Bogan, C; Bond, C; Bondu, F; Bonelli, L; Bonnand, R; Bork, R; Born, M; Boschi, V; Bose, Sukanta; Bosi, L; Bradaschia, C; Brady, P R; Braginsky, V B; Branchesi, M; Brau, J E; Briant, T; Bridges, D O; Brillet, A; Brinkmann, M; Brisson, V; Brooks, A F; Brown, D A; Brown, D D; Brückner, F; Buchman, S; Bulik, T; Bulten, H J; Buonanno, A; Burman, R; Buskulic, D; Buy, C; Cadonati, L; Cagnoli, G; Bustillo, J Calderón; Calloni, E; Camp, J B; Campsie, P; Cannon, K C; Canuel, B; Cao, J; Capano, C D; Carbognani, F; Carbone, L; Caride, S; Castiglia, A; Caudill, S; Cavaglià, M; Cavalier, F; Cavalieri, R; Celerier, C; Cella, G; Cepeda, C; Cesarini, E; Chakraborty, R; Chalermsongsak, T; Chamberlin, S J; Chao, S; Charlton, P; Chassande-Mottin, E; Chen, X; Chen, Y; Chincarini, A; Chiummo, A; Cho, H S; Chow, J; Christensen, N; Chu, Q; Chua, S S Y; Chung, S; Ciani, G; Clara, F; Clark, J A; Cleva, F; Coccia, E; Cohadon, P -F; Colla, A; Collette, C; Colombini, M; Cominsky, L; Constancio, M; Conte, A; Cook, D; Corbitt, T R; Cordier, M; Cornish, N; Corpuz, A; Corsi, A; Costa, C A; Coughlin, M W; Coughlin, S; Coulon, J -P; Countryman, S; Couvares, P; Coward, D M; Cowart, M; Coyne, D C; Coyne, R; Craig, K; Creighton, J D E; Crowder, S G; Cumming, A; Cunningham, L; Cuoco, E; Dahl, K; Canton, T Dal; Damjanic, M; Danilishin, S L; D'Antonio, S; Danzmann, K; Dattilo, V; Daveloza, H; Davier, M; Davies, G S; Daw, E J; Day, R; Dayanga, T; Debreczeni, G; Degallaix, J; Deléglise, S; Del Pozzo, W; Denker, T; Dent, T; Dereli, H; Dergachev, V; De Rosa, R; DeRosa, R T; DeSalvo, R; Dhurandhar, S; Díaz, M; Di Fiore, L; Di Lieto, A; Di Palma, I; Di Virgilio, A; Donath, A; Donovan, F; Dooley, K L; Doravari, S; Dossa, S; Douglas, R; Downes, T P; Drago, M; Drever, R W P; Driggers, J C; Du, Z; Dwyer, S; Eberle, T; Edo, T; Edwards, M; Effler, A; Eggenstein, H; Ehrens, P; Eichholz, J; Eikenberry, S S; Endr?czi, G; Essick, R; Etzel, T; Evans, M; Evans, T; Factourovich, M; Fafone, V; Fairhurst, S; Fang, Q; Farinon, S; Farr, B; Farr, W M; Favata, M; Fehrmann, H; Fejer, M M; Feldbaum, D; Feroz, F; Ferrante, I; Ferrini, F; Fidecaro, F; Finn, L S; Fiori, I; Fisher, R P; Flaminio, R; Fournier, J -D; Franco, S; Frasca, S; Frasconi, F; Frede, M; Frei, Z; Freise, A; Frey, R; Fricke, T T; Fritschel, P; Frolov, V V; Fulda, P; Fyffe, M; Gair, J; Gammaitoni, L; Gaonkar, S; Garufi, F; Gehrels, N; Gemme, G; Genin, E; Gennai, A; Ghosh, S; Giaime, J A; Giardina, K D; Giazotto, A; Gill, C; Gleason, J; Goetz, E; Goetz, R; Gondan, L; González, G; Gordon, N; Gorodetsky, M L; Gossan, S; Goßler, S; Gouaty, R; Gräf, C; Graff, P B; Granata, M; Grant, A; Gras, S; Gray, C; Greenhalgh, R J S; Gretarsson, A M; Groot, P; Grote, H; Grover, K; Grunewald, S; Guidi, G M; Guido, C; Gushwa, K; Gustafson, E K; Gustafson, R; Hammer, D; Hammond, G; Hanke, M; Hanks, J; Hanna, C; Hanson, J; Harms, J; Harry, G M; Harry, I W; Harstad, E D; Hart, M; Hartman, M T; Haster, C -J; Haughian, K; Heidmann, A; Heintze, M; Heitmann, H; Hello, P; Hemming, G; Hendry, M; Heng, I S; Heptonstall, A W; Heurs, M; Hewitson, M; Hild, S; Hoak, D; Hodge, K A; Holt, K; Hooper, S; Hopkins, P; Hosken, D J; Hough, J; Howell, E J; Hu, Y; Hughey, B; Husa, S; Huttner, S H; Huynh, M; Huynh-Dinh, T; Ingram, D R; Inta, R; Isogai, T; Ivanov, A; Iyer, B R; Izumi, K; Jacobson, M; James, E; Jang, H; Jaranowski, P; Ji, Y; Jiménez-Forteza, F; Johnson, W W; Jones, D I; Jones, R; Jonker, R J G; Ju, L; K, Haris; Kalmus, P; Kalogera, V; Kandhasamy, S; Kang, G; Kanner, J B; Karlen, J; Kasprzack, M; Katsavounidis, E; Katzman, W; Kaufer, H; Kawabe, K; Kawazoe, F; Kéfélian, F; Keiser, G M; Keitel, D; Kelley, D B; Kells, W; Khalaidovski, A; Khalili, F Y; Khazanov, E A; Kim, C; Kim, K; Kim, N; Kim, N G; Kim, Y -M; King, E J; King, P J; Kinzel, D L; Kissel, J S; Klimenko, S; Kline, J; Koehlenbeck, S; Kokeyama, K; Kondrashov, V; Koranda, S; Korth, W Z; Kowalska, I
2014-01-01T23:59:59.000Z
The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational-wave astrophysics communities. The purpose of NINJA is to study the ability to detect gravitational waves emitted from merging binary black holes and recover their parameters with next-generation gravitational-wave observatories. We report here on the results of the second NINJA project, NINJA-2, which employs 60 complete binary black hole hybrid waveforms consisting of a numerical portion modelling the late inspiral, merger, and ringdown stitched to a post-Newtonian portion modelling the early inspiral. In a "blind injection challenge" similar to that conducted in recent LIGO and Virgo science runs, we added 7 hybrid waveforms to two months of data recolored to predictions of Advanced LIGO and Advanced Virgo sensitivity curves during their first observing runs. The resulting data was analyzed by gravitational-wave detection algorithms and 6 of the waveforms were recovered w...
Banerjee, Debjyoti
microfluidics Jonathan Siegrist,*a Mary Amasia,a Navdeep Singh,b Debjyoti Banerjeeb and Marc Madoua Received 1st analysis of microchamber filling in centrifugal microfluidics is presented. In the development of micro on centrifugal microfluidic platforms, numerical modeling using the Volume of Fluids method is performed
Numerical modelling of avalanches based on Saint-Venant equations using a kinetic scheme
Paris-Sud XI, Université de
avalanches are treated here as a dry granular flow with Coulomb-type behavior. The numerical finite volume of an avalanche over simplified topography. Coulomb-type behavior with constant and variable friction angle modelling, Coulomb friction, Saint-Venant equations, finite volume kinetic scheme. 1 Introduction Granular
Numerical methods for the simulation of a corrosion model in a nuclear waste deep repository $
Paris-Sud XI, Université de
Numerical methods for the simulation of a corrosion model in a nuclear waste deep repository $ C of the French nuclear waste management agency ANDRA, investigations are conducted to optimize and finalize by the Nuclear Waste Management Agency ANDRA Corresponding author. Phone: +49 30 20372 560, Fax: +49 30 2044975
Numerical analysis of electric field formulations of the eddy current model
RodrÃguez, Rodolfo
Numerical analysis of electric field formulations of the eddy current model Alfredo BermÂ´udez1 methods for the numeri- cal solution of the eddy current problem in a bounded conducting domain crossed): 78M10, 65N30 Key words Low-frequency harmonic Maxwell equations, eddy currents, finite elements
Mathematical and numerical analysis of a transient non-linear axisymmetric eddy current model
RodrÃguez, Rodolfo
Mathematical and numerical analysis of a transient non-linear axisymmetric eddy current model the theoretically predicted behavior of the method, are reported. Keywords transient eddy current Â· axisymmetric is the accurate computation of power losses in the ferromagnetic components of the core due to hysteresis and eddy-current
NUMERICAL SOLUTION OF A TRANSIENT NON-LINEAR AXISYMMETRIC EDDY CURRENT MODEL WITH NON-LOCAL
RodrÃguez, Rodolfo
NUMERICAL SOLUTION OF A TRANSIENT NON-LINEAR AXISYMMETRIC EDDY CURRENT MODEL WITH NON@ing-mat.udec.cl This paper deals with an axisymmetric transient eddy current problem in conductive nonlinear magnetic media of the proposed scheme. Keywords: transient eddy current problem; electromagnetic losses; nonlinear magnetic
A simple numerical model of the apparent loss of eddy current conductivity due to surface roughness
Nagy, Peter B.
A simple numerical model of the apparent loss of eddy current conductivity due to surface roughness of eddy current conductivity has been suggested as a possible means to allow the nondestructive evaluation, the path of the eddy current must follow a more tortuous route in the material, which produces a reduction
Ensemble Kalman Filter Data Assimilation in a 1D Numerical Model Used for Fog Forecasting
Ribes, Aurélien
Ensemble Kalman Filter Data Assimilation in a 1D Numerical Model Used for Fog Forecasting SAMUEL RE significant. This led to the implementation of an ensemble Kalman filter (EnKF) within COBEL-ISBA. The new by using an ensemble Kalman filter (EnKF; Evensen 1994, 2003). Theoreti- cally, ensemble filters
NUMERICAL MODELLING OF THE DEEP IMPACT MISSION EXPERIMENT K. Wnnemann1 , and H. J. Melosh1
Collins, Gareth
NUMERICAL MODELLING OF THE DEEP IMPACT MISSION EXPERIMENT K. Wünnemann1 , G. S. Collins2 , and H. J@lpl.arizona.edu, 2 Department of Earth Science and Engineering, Imperial College London, London SW7 2AZ, UK, g.collins@imperial.ac.uk . Introduction: NASA's Deep Impact Mission (launched January 2005) will provide, for the first time ever
Boyer, Edmond
A LARGE SCALE CONTINUUM-DISCRETE NUMERICAL MODELLING: APPLICATION TO OVERBURDEN DAMAGE OF A SALT damage on top of an underground solution mining, an in-situ experiment is undertaken above a salt cavity in the Lorraine region (NE of France). The overburden overlying the salt cavity is characterized by a competent
NUMERICAL MODELING OF 3D ORGANIC SOLAR CELLS Presented to the
Kassegne, Samuel Kinde
NUMERICAL MODELING OF 3D ORGANIC SOLAR CELLS _______________ A Thesis Presented to the Faculty and mass equilibrium reactions that occur in the 3D configuration of the organic solar cell. #12;vi TABLE.......................................................................................3 1.3 3-D Organic Solar Cell
Mathematical, physical and numerical principles essential for models of turbulent mixing
Sharp, David Howland [Los Alamos National Laboratory; Lim, Hyunkyung [STONY BROOK UNIV; Yu, Yan [STONY BROOK UNIV; Glimm, James G [STONY BROOK UNIV
2009-01-01T23:59:59.000Z
We propose mathematical, physical and numerical principles which are important for the modeling of turbulent mixing, especially the classical and well studied Rayleigh-Taylor and Richtmyer-Meshkov instabilities which involve acceleration driven mixing of a fluid discontinuity layer, by a steady accerleration or an impulsive force.
A numerical modeling study on desert oasis self-supporting mechanisms
Chu, Peter C.
A numerical modeling study on desert oasis self-supporting mechanisms Peter C. Chua, *, Shihua Lub February 2005 Abstract Oasis self-supporting mechanisms due to oasis breeze circulation (OBC) are proposed from the oasis makes the oasis surface colder than the surrounding desert surface. The sensible heat
ENINALS Experimental and Numerical Investigations of Nonlinearity in soils using Advanced of Nonlinearity in Soils Using Advanced Laboratory-scaled Models Concluding Workshop Ispra, May 28-30 2013 IFSTTAR linearity in heterogeneous soil deposits due to strong motion shaking METODOLOGY - engineering
January 2, 2008 Numerical modeling of the effect of carbon dioxide
Boyer, Edmond
January 2, 2008 Numerical modeling of the effect of carbon dioxide sequestration on the rate souterrain de dioxyde de carbone sur la déformation des calcaires par dissolution sous contrainte: résultats;Abstract When carbon dioxide (CO2) is injected into an aquifer or a depleted geological reservoir, its
The Hamiltonian Particle-Mesh (HPM) method for numerical modeling of atmospheric flows.
Kim, Guebuem
The Hamiltonian Particle-Mesh (HPM) method for numerical modeling of atmospheric flows. Seoleun Shin 15. Feb. 2011 Abstract The Hamiltonian Particle-Mesh (HPM) method is an interesting alternative have developed schemes based on the HPM method for the shallow-water equations on the sphere, nonhydro
Numerical modeling of extreme rogue waves generated by directional energy focusing
Grilli, Stéphan T.
Numerical modeling of extreme rogue waves generated by directional energy focusing Christophe angle of directional energy focusing. We find that an over- turning rogue wave can have different are characterized by their brief occurrence in space and time, resulting from a local focusing of wave energy
Numerical Modelling of Unsaturated Flow in Uniform and Heterogeneous Waste Rock Piles
Aubertin, Michel
Numerical Modelling of Unsaturated Flow in Uniform and Heterogeneous Waste Rock Piles O Fala1 , M Aubertin1,3 , J Molson1 , B Bussière2,3 , G W Wilson4 , R Chapuis1 and V Martin1 ABSTRACT Waste rock piles these piles, many physical, geochemical and biological processes can contribute to the production of AMD
Numerical Passage from Radiative Heat Transfer to Nonlinear Diffusion Models \\Lambda
Schmeiser, Christian
Numerical Passage from Radiative Heat Transfer to Nonlinear Diffusion Models \\Lambda A. Klar y C. Schmeiser z Abstract Radiative heat transfer equations including heat conduction are considÂ ered situations are presented. Keywords. radiative heat transfer, asymptotic analysis, nonlinear diffusion limit
NUMERICAL MODELLING OF AUTOGENOUS HEALING AND RECOVERY OF MECHANICAL PROPERTIES IN ULTRA-HIGH
Boyer, Edmond
in the cementitious matrix can react with carbon dioxide dissolved in the water filling the crack. Autogenous healingNUMERICAL MODELLING OF AUTOGENOUS HEALING AND RECOVERY OF MECHANICAL PROPERTIES IN ULTRA into the crack and leads to a partial recovery of mechanical properties (Young's modulus, tensile strength
Numerical modeling of heat transfer and fluid flow in rotor-stator cavities with throughflow
Boyer, Edmond
Numerical modeling of heat transfer and fluid flow in rotor-stator cavities with throughflow S in a rotor-stator cavity subjected to a superimposed throughflow with heat transfer. Nu- merical predictions field from the heat transfer process. The turbulent flux is approximated by a gradient hypothesis
LABORATORY OBSERVATIONS AND NUMERICAL MODELING OF THE EFFECTS OF AN ARRAY OF WAVE ENERGY CONVERTERS
Haller, Merrick
1 LABORATORY OBSERVATIONS AND NUMERICAL MODELING OF THE EFFECTS OF AN ARRAY OF WAVE ENERGY of wave energy converters (WECs) on water waves through the analysis of extensive laboratory experiments absorption is a reasonable predictor of the effect of WECs on the far field. Keywords: wave- energy; spectral
Numerical Simulation of the December 26, 2004 Indian Ocean Tsunami using a Boussinesq model
Kirby, James T.
tsunami education or tsunami warning system in the re- gion exaccerbated the number of fatalities, even. Scientists had been warning of the growing exposure of coastal residents to tsunami hazards for yearsNumerical Simulation of the December 26, 2004 Indian Ocean Tsunami using a Boussinesq model Philip
September 25, 2006 Numerical modeling of the effect of carbon dioxide
Paris-Sud XI, Université de
September 25, 2006 Numerical modeling of the effect of carbon dioxide sequestration on the rate souterrain de dioxyde de carbone sur la déformation des calcaires par dissolution sous contrainte: résultats@obs.ujf- grenoble.fr, marielle.collombet@ujf-grenoble.fr, yleguen@lgit.obs.ujf-grenoble.fr. #12;Abstract When carbon
Numerical model for vibration damping resulting from the first-order phase transformations
Melnik, Roderick
: Martensite transformation; Thermo-mechanical coupling; Vibration damping; GinzburgLandau theory 1Numerical model for vibration damping resulting from the first-order phase transformations L LandauGinzburg theory that couples nonlinear mechanical and thermal fields. The free energy function
NUMERICAL MODELLING OF THERMAL-ELECTRICAL PHENOMENA IN SPARK PLASMA SINTERING
Boyer, Edmond
NUMERICAL MODELLING OF THERMAL-ELECTRICAL PHENOMENA IN SPARK PLASMA SINTERING P. Mondaleka , L'Etudes Structurales), France c Université Paul Sabatier, Toulouse, France Abstract. Spark Plasma Sintering belongs: Finite element method, Spark plasma sintering, powder compaction. INTRODUCTION Spark Plasma Sintering
DROP-SCALE NUMERICAL MODELING OF CHEMICAL PARTITIOING DURING CLOUD HYDROMETEOR FREEZING
Stuart, Amy L.
DROP-SCALE NUMERICAL MODELING OF CHEMICAL PARTITIOING DURING CLOUD HYDROMETEOR FREEZING A.L. Stuart freezing provide greatly varying estimates of the retention efficiency of volatile solutes (e.g., Lamb understanding of the dependence of partitioning on chemical properties and freezing conditions (Stuart
Turova, Varvara
International Series of Numerical Mathematics, Vol. 160, 521540 Freezing of Living Cells, stresses arising due to non-simultaneous freezing of water in- side and outside of cells are modeled and outside of living cells during freezing is derived by applying an appropriate averaging technique
COMPUTATIONAL CHALLENGES IN THE NUMERICAL TREATMENT OF LARGE AIR POLLUTION MODELS
Ostromsky, Tzvetan
COMPUTATIONAL CHALLENGES IN THE NUMERICAL TREATMENT OF LARGE AIR POLLUTION MODELS I. DIMOV , K. GEORGIEVy, TZ. OSTROMSKY , R. J. VAN DER PASz, AND Z. ZLATEVx Abstract. The air pollution, and especially the reduction of the air pollution to some acceptable levels, is an important environmental problem, which
Nonlinear inverse problem for a model of ion-exchange filter: numerical recovery of parameters
) and pressure (15 MPa) of hot steam [2]. Some units are made of cheap corrosion and heat-resistant steel which1 Nonlinear inverse problem for a model of ion-exchange filter: numerical recovery of parameters]. Power-generating units of TPP operate under severe corrosive conditions: high temperature (515 - 530°C
NUMERICAL VERIFICATION OF THE RELAP-7 CORE CHANNEL SINGLE-PHASE MODEL
Haihua Zhao; Ling Zou; Hongbin Zhang; Richard Martineau
2014-06-01T23:59:59.000Z
The RELAP-7 code is the next generation of nuclear reactor system safety analysis code being developed at the Idaho National Laboratory (INL). All the physics in RELAP-7 are fully coupled and the errors resulted from the traditional operator-splitting approach are eliminated. By using 2nd order methods in both time and space and eliminating operator-splitting errors, the numerical error of RELAP-7 can be minimized. Numerical verification is the process to verify the orders of numerical methods. It is an important part of modern verification and validation process. The core channel component in RELAP-7 is designed to simulate coolant flow as well as the conjugated heat transfer between coolant flow and the fuel rod. A special treatment at fuel centerline to avoid numerical singularity for the cylindrical heat conduction in the continuous finite element mesh is discussed. One steady state test case and one fast power up transient test case are utilized for the verification of the core channel model with single-phase flow. Analytical solution for the fuel pin temperature and figures of merit such as peak clad temperature and peak fuel temperature are used to define numerical errors. These cases prove that the mass and energy are well conserved and 2nd order convergence rates for both time and space are achieved in the core channel model.
Coclite, A; De Palma, P; Pascazio, G
2015-01-01T23:59:59.000Z
The present paper deals with the numerical study of high pressure LOx/H2 or LOx/hydrocarbon combustion for propulsion systems. The present research effort is driven by the continued interest in achieving low cost, reliable access to space and more recently, by the renewed interest in hypersonic transportation systems capable of reducing time-to-destination. Moreover, combustion at high pressure has been assumed as a key issue to achieve better propulsive performance and lower environmental impact, as long as the replacement of hydrogen with a hydrocarbon, to reduce the costs related to ground operations and increase flexibility. The current work provides a model for the numerical simulation of high- pressure turbulent combustion employing detailed chemistry description, embedded in a RANS equations solver with a Low Reynolds number k-omega turbulence model. The model used to study such a combustion phenomenon is an extension of the standard flamelet-progress-variable (FPV) turbulent combustion model combined ...
Modelling and Numerical Simulation of Gas Migration in a Nuclear Waste Repository
Bourgeat, Alain; Smai, Farid
2010-01-01T23:59:59.000Z
We present a compositional compressible two-phase, liquid and gas, flow model for numerical simulations of hydrogen migration in deep geological radioactive waste repository. This model includes capillary effects and the gas diffusivity. The choice of the main variables in this model, Total or Dissolved Hydrogen Mass Concentration and Liquid Pressure, leads to a unique and consistent formulation of the gas phase appearance and disappearance. After introducing this model, we show computational evidences of its adequacy to simulate gas phase appearance and disappearance in different situations typical of underground radioactive waste repository.
Numerical upscaling for the eddy-current model with stochastic magnetic materials
Eberhard, Jens P. [Computer Simulation Technology, Bad Nauheimer Strasse, 19, D-64289 Darmstadt (Germany)], E-mail: jens.eberhard@cst.com; Popovic, Dan [Simulation in Technology, University of Heidelberg, Im Neuenheimer Feld 368, D-69120 Heidelberg (Germany)], E-mail: dan.popovic@stud.uni-heidelberg.de; Wittum, Gabriel [Simulation in Technology, University of Heidelberg, Im Neuenheimer Feld 368, D-69120 Heidelberg (Germany)], E-mail: wittum@uni-hd.de
2008-04-01T23:59:59.000Z
This paper deals with the upscaling of the time-harmonic Maxwell equations for heterogeneous media. We analyze the eddy-current approximation of Maxwell's equations to describe the electric field for heterogeneous, isotropic magnetic materials. The magnetic permeability of the materials is assumed to have random heterogeneities described by a Gaussian random field. We apply the so-called Coarse Graining method to develop a numerical upscaling of the eddy-current model. The upscaling uses filtering and averaging procedures in Fourier space which results in a formulation of the eddy-current model on coarser resolution scales where the influence of sub-scale fluctuations is modeled by effective scale- and space-dependent reluctivity tensors. The effective reluctivity tensors can be obtained by solving local partial differential equations which contain a Laplacian as well as a curl-curl operator. We present a computational method how the equation of the combined operators can be discretized and solved numerically using an extended variational formulation compared to standard discretizations. We compare the results of the numerical upscaling of the eddy-current model with theoretical results of Eberhard [J.P. Eberhard, Upscaling for the time-harmonic Maxwell equations with heterogeneous magnetic materials, Physical Review E 72 (3), (2005)] and obtain a very good agreement.
Numerical study of a slip-link model for polymer melts and nanocomposites
Diego Delbiondo; Elian Masnada; Samy Merabia; Marc Couty; Jean-Louis Barrat
2013-06-10T23:59:59.000Z
We present a numerical study of the slip link model introduced by Likhtman for describing the dy- namics of dense polymer melts. After reviewing the technical aspects associated with the implemen- tation of the model, we extend previous work in several directions. The dependence of the relaxation modulus with the slip link density and the slip link stiffness is reported. Then the nonlinear rheolog- ical properties of the model, for a particular set of parameters, are explored. Finally, we introduce excluded volume interactions in a mean field such as manner in order to describe inhomogeneous systems, and we apply this description to a simple nanocomposite model. With this extension, the slip link model appears as a simple and generic model of a polymer melt, that can be used as an alternative to molecular dynamics for coarse grained simulations of complex polymeric systems.
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
François, Marianne M.
2015-01-01T23:59:59.000Z
A review of recent advances made in numerical methods and algorithms within the volume tracking framework is presented. The volume tracking method, also known as the volume-of-fluid method has become an established numerical approach to model and simulate interfacial flows. Its advantage is its strict mass conservation. However, because the interface is not explicitly tracked but captured via the material volume fraction on a fixed mesh, accurate estimation of the interface position, its geometric properties and modeling of interfacial physics in the volume tracking framework remain difficult. Several improvements have been made over the last decade to address these challenges.more »In this paper, the multimaterial interface reconstruction method via power diagram, curvature estimation via heights and mean values and the balanced-force algorithm for surface tension are highlighted.« less
Numerical-Model Investigation of the Hydrothermal Regime of a Straight-Through Shallow Cooling Pond
Sokolov, A. S. [JSC 'VNIIG im. B. E. Vedeneeva' (Russian Federation)] [JSC 'VNIIG im. B. E. Vedeneeva' (Russian Federation)
2013-11-15T23:59:59.000Z
A mathematic model based on solution of hydrodynamics and heat-transfer equations by the finite-element method is constructed to predict the hydrothermal regime of a straight-through shallow cooling pond, which provides cooling circulating water to a repository of spent nuclear fuel. Numerical experiments made it possible to evaluate the influence exerted by wind conditions and flow rate of water in the river on the temperature of the circulating water.
An investigation of analytical and numerical sucker rod pumping mathematical models
Schafer, Donald Joseph
1987-01-01T23:59:59.000Z
of MASTER OF SCIENCE May 1987 Major Subject: Petroleum Engineering AN INVESTIGATION OF ANALYTICAL AND NUMERICAL SUCKER ROD PUMPING MATHEMATICAL MODELS A Thesis by DONALD JOSEPH SCHAFER Approved as to style and content by: 7d JW. J ni (Chai... to Sucker Rod Pumping Research, Inc. , developed a method for computing downhole forces and displacements using an analog computer simulation. The procedure that results from this work, commonly called the API method, considers the total sucker rod...
Numerical modeling of hydraulic fracture problem in permeable medium using cohesive zone model
Paris-Sud XI, Université de
or radioactive waste [30], geothermal Corresponding author Email address: benoit.carrier@enpc.fr (Benoit Carrier processes. During the last sixty years, numerous papers [3, 7, 14, 21, 35, 34, 41, 25, 27, 22, 38, 36, 1. In the recent years, a scaling and asymptotic framework was built to determine the influence of the physical
Moist processes and the quasi-hydrostatic approximation in a mesoscale numerical model
Kennedy, Charles Joseph
1987-01-01T23:59:59.000Z
of Committee) James P. McGuirk (Member) J'ohn M. Klinck (Member) James R. Sco ns (Head of Department) December 1987 ABSTRACT Moist Processes and the Ouasi-Hydrostatic Approximation in a Mesoscale Numerical Model. (December 1987) Charles Joseph...HV)ds' ? gHp s + gHps a dg 1 gt = (gt), s 1 1 (19) the pressure tendency at the model top equation: g f V ~ (pHV)ds' ? VS Vp Q ( el 1 + 0 Yp CpT ? V (H0) ds' )' ? ) (20) Richardson's equation for vertical motion: s f , , f Id d'D &DDVdd ' ? 0 D, 0...
Thermo--inertial bouncing of a relativistic collapsing sphere: A numerical model
L. Herrera; A. Di Prisco; W. Barreto
2005-12-05T23:59:59.000Z
We present a numerical model of a collapsing radiating sphere, whose boundary surface undergoes bouncing due to a decreasing of its inertial mass density (and, as expected from the equivalence principle, also of the ``gravitational'' force term) produced by the ``inertial'' term of the transport equation. This model exhibits for the first time the consequences of such an effect, and shows that under physically reasonable conditions this decreasing of the gravitational term in the dynamic equation may be large enough as to revert the collapse and produce a bouncing of the boundary surface of the sphere.
Bammann, Douglas J.; Johnson, G. C. (University of California, Berkeley, CA); Marin, Esteban B.; Regueiro, Richard A. (University of Colorado, Boulder, CO)
2006-01-01T23:59:59.000Z
In this report we present the formulation of the physically-based Evolving Microstructural Model of Inelasticity (EMMI) . The specific version of the model treated here describes the plasticity and isotropic damage of metals as being currently applied to model the ductile failure process in structural components of the W80 program . The formulation of the EMMI constitutive equations is framed in the context of the large deformation kinematics of solids and the thermodynamics of internal state variables . This formulation is focused first on developing the plasticity equations in both the relaxed (unloaded) and current configurations. The equations in the current configuration, expressed in non-dimensional form, are used to devise the identification procedure for the plasticity parameters. The model is then extended to include a porosity-based isotropic damage state variable to describe the progressive deterioration of the strength and mechanical properties of metals induced by deformation . The numerical treatment of these coupled plasticity-damage constitutive equations is explained in detail. A number of examples are solved to validate the numerical implementation of the model.
W. Schmidt; J. C. Niemeyer; W. Hillebrandt
2006-01-23T23:59:59.000Z
We present a one-equation subgrid scale model that evolves the turbulence energy corresponding to unresolved velocity fluctuations in large eddy simulations. The model is derived in the context of the Germano consistent decomposition of the hydrodynamical equations. The eddy-viscosity closure for the rate of energy transfer from resolved toward subgrid scales is localised by means of a dynamical procedure for the computation of the closure parameter. Therefore, the subgrid scale model applies to arbitrary flow geometry and evolution. For the treatment of microscopic viscous dissipation a semi-statistical approach is used, and the gradient-diffusion hypothesis is adopted for turbulent transport. A priori tests of the localised eddy-viscosity closure and the gradient-diffusion closure are made by analysing data from direct numerical simulations. As an a posteriori testing case, the large eddy simulation of thermonuclear combustion in forced isotropic turbulence is discussed. We intend the formulation of the subgrid scale model in this paper as a basis for more advanced applications in numerical simulations of complex astrophysical phenomena involving turbulence.
Cirpka, Olaf Arie
on Numerical Models for Carbon Dioxide Storage in Geological Formations 6/16 Transport Bacteria in water Hydrosystemmodellierung Workshop on Numerical Models for Carbon Dioxide Storage in Geological Formations 1/16 Modelling April 2008 Workshop on Numerical Models for Carbon Dioxide Storage in Geological Formations #12
Schmidt, W; Niemeyer, J C
2006-01-01T23:59:59.000Z
We present a one-equation subgrid scale model that evolves the turbulence energy corresponding to unresolved velocity fluctuations in large eddy simulations. The model is derived in the context of the Germano consistent decomposition of the hydrodynamical equations. The eddy-viscosity closure for the rate of energy transfer from resolved toward subgrid scales is localised by means of a dynamical procedure for the computation of the closure parameter. Therefore, the subgrid scale model applies to arbitrary flow geometry and evolution. For the treatment of microscopic viscous dissipation a semi-statistical approach is used, and the gradient-diffusion hypothesis is adopted for turbulent transport. A priori tests of the localised eddy-viscosity closure and the gradient-diffusion closure are made by analysing data from direct numerical simulations. As an a posteriori testing case, the large eddy simulation of thermonuclear combustion in forced isotropic turbulence is discussed. We intend the formulation of the sub...
The LIGO Scientific Collaboration; the Virgo Collaboration; the NINJA-2 Collaboration; :; J. Aasi; B. P. Abbott; R. Abbott; T. Abbott; M. R. Abernathy; T. Accadia; F. Acernese; K. Ackley; C. Adams; T. Adams; P. Addesso; R. X. Adhikari; C. Affeldt; M. Agathos; N. Aggarwal; O. D. Aguiar; A. Ain; P. Ajith; A. Alemic; B. Allen; A. Allocca; D. Amariutei; M. Andersen; R. Anderson; S. B. Anderson; W. G. Anderson; K. Arai; M. C. Araya; C. Arceneaux; J. Areeda; S. M. Aston; P. Astone; P. Aufmuth; C. Aulbert; L. Austin; B. E. Aylott; S. Babak; P. T. Baker; G. Ballardin; S. W. Ballmer; J. C. Barayoga; M. Barbet; B. C. Barish; D. Barker; F. Barone; B. Barr; L. Barsotti; M. Barsuglia; M. A. Barton; I. Bartos; R. Bassiri; A. Basti; J. C. Batch; J. Bauchrowitz; Th. S. Bauer; B. Behnke; M. Bejger; M. G. Beker; C. Belczynski; A. S. Bell; C. Bell; G. Bergmann; D. Bersanetti; A. Bertolini; J. Betzwieser; P. T. Beyersdorf; I. A. Bilenko; G. Billingsley; J. Birch; S. Biscans; M. Bitossi; M. A. Bizouard; E. Black; J. K. Blackburn; L. Blackburn; D. Blair; S. Bloemen; M. Blom; O. Bock; T. P. Bodiya; M. Boer; G. Bogaert; C. Bogan; C. Bond; F. Bondu; L. Bonelli; R. Bonnand; R. Bork; M. Born; V. Boschi; Sukanta Bose; L. Bosi; C. Bradaschia; P. R. Brady; V. B. Braginsky; M. Branchesi; J. E. Brau; T. Briant; D. O. Bridges; A. Brillet; M. Brinkmann; V. Brisson; A. F. Brooks; D. A. Brown; D. D. Brown; F. Brückner; S. Buchman; T. Bulik; H. J. Bulten; A. Buonanno; R. Burman; D. Buskulic; C. Buy; L. Cadonati; G. Cagnoli; J. Calderón Bustillo; E. Calloni; J. B. Camp; P. Campsie; K. C. Cannon; B. Canuel; J. Cao; C. D. Capano; F. Carbognani; L. Carbone; S. Caride; A. Castiglia; S. Caudill; M. Cavaglià; F. Cavalier; R. Cavalieri; C. Celerier; G. Cella; C. Cepeda; E. Cesarini; R. Chakraborty; T. Chalermsongsak; S. J. Chamberlin; S. Chao; P. Charlton; E. Chassande-Mottin; X. Chen; Y. Chen; A. Chincarini; A. Chiummo; H. S. Cho; J. Chow; N. Christensen; Q. Chu; S. S. Y. Chua; S. Chung; G. Ciani; F. Clara; J. A. Clark; F. Cleva; E. Coccia; P. -F. Cohadon; A. Colla; C. Collette; M. Colombini; L. Cominsky; M. Constancio Jr.; A. Conte; D. Cook; T. R. Corbitt; M. Cordier; N. Cornish; A. Corpuz; A. Corsi; C. A. Costa; M. W. Coughlin; S. Coughlin; J. -P. Coulon; S. Countryman; P. Couvares; D. M. Coward; M. Cowart; D. C. Coyne; R. Coyne; K. Craig; J. D. E. Creighton; S. G. Crowder; A. Cumming; L. Cunningham; E. Cuoco; K. Dahl; T. Dal Canton; M. Damjanic; S. L. Danilishin; S. D'Antonio; K. Danzmann; V. Dattilo; H. Daveloza; M. Davier; G. S. Davies; E. J. Daw; R. Day; T. Dayanga; G. Debreczeni; J. Degallaix; S. Deléglise; W. Del Pozzo; T. Denker; T. Dent; H. Dereli; V. Dergachev; R. De Rosa; R. T. DeRosa; R. DeSalvo; S. Dhurandhar; M. Díaz; L. Di Fiore; A. Di Lieto; I. Di Palma; A. Di Virgilio; A. Donath; F. Donovan; K. L. Dooley; S. Doravari; S. Dossa; R. Douglas; T. P. Downes; M. Drago; R. W. P. Drever; J. C. Driggers; Z. Du; S. Dwyer; T. Eberle; T. Edo; M. Edwards; A. Effler; H. Eggenstein; P. Ehrens; J. Eichholz; S. S. Eikenberry; G. Endr?czi; R. Essick; T. Etzel; M. Evans; T. Evans; M. Factourovich; V. Fafone; S. Fairhurst; Q. Fang; S. Farinon; B. Farr; W. M. Farr; M. Favata; H. Fehrmann; M. M. Fejer; D. Feldbaum; F. Feroz; I. Ferrante; F. Ferrini; F. Fidecaro; L. S. Finn; I. Fiori; R. P. Fisher; R. Flaminio; J. -D. Fournier; S. Franco; S. Frasca; F. Frasconi; M. Frede; Z. Frei; A. Freise; R. Frey; T. T. Fricke; P. Fritschel; V. V. Frolov; P. Fulda; M. Fyffe; J. Gair; L. Gammaitoni; S. Gaonkar; F. Garufi; N. Gehrels; G. Gemme; E. Genin; A. Gennai; S. Ghosh; J. A. Giaime; K. D. Giardina; A. Giazotto; C. Gill; J. Gleason; E. Goetz; R. Goetz; L. Gondan; G. González; N. Gordon; M. L. Gorodetsky; S. Gossan; S. Goßler; R. Gouaty; C. Gräf; P. B. Graff; M. Granata; A. Grant; S. Gras; C. Gray; R. J. S. Greenhalgh; A. M. Gretarsson; P. Groot; H. Grote; K. Grover; S. Grunewald; G. M. Guidi; C. Guido; K. Gushwa; E. K. Gustafson; R. Gustafson; D. Hammer; G. Hammond; M. Hanke; J. Hanks; C. Hanna; J. Hanson; J. Harms; G. M. Harry; I. W. Harry; E. D. Harstad; M. Hart; M. T. Hartman; C. -J. Haster; K. Haughian; A. Heidmann; M. Heintze; H. Heitmann; P. Hello; G. Hemming; M. Hendry; I. S. Heng; A. W. Heptonstall; M. Heurs; M. Hewitson; S. Hild; D. Hoak; K. A. Hodge; K. Holt; S. Hooper; P. Hopkins; D. J. Hosken; J. Hough; E. J. Howell; Y. Hu; B. Hughey; S. Husa; S. H. Huttner; M. Huynh; T. Huynh-Dinh; D. R. Ingram; R. Inta; T. Isogai; A. Ivanov; B. R. Iyer; K. Izumi; M. Jacobson; E. James; H. Jang; P. Jaranowski; Y. Ji; F. Jiménez-Forteza; W. W. Johnson; D. I. Jones; R. Jones; R. J. G. Jonker; L. Ju; Haris K; P. Kalmus; V. Kalogera; S. Kandhasamy; G. Kang; J. B. Kanner; J. Karlen; M. Kasprzack; E. Katsavounidis; W. Katzman; H. Kaufer; K. Kawabe; F. Kawazoe; F. Kéfélian; G. M. Keiser; D. Keitel; D. B. Kelley; W. Kells; A. Khalaidovski
2014-01-05T23:59:59.000Z
The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational-wave astrophysics communities. The purpose of NINJA is to study the ability to detect gravitational waves emitted from merging binary black holes and recover their parameters with next-generation gravitational-wave observatories. We report here on the results of the second NINJA project, NINJA-2, which employs 60 complete binary black hole hybrid waveforms consisting of a numerical portion modelling the late inspiral, merger, and ringdown stitched to a post-Newtonian portion modelling the early inspiral. In a "blind injection challenge" similar to that conducted in recent LIGO and Virgo science runs, we added 7 hybrid waveforms to two months of data recolored to predictions of Advanced LIGO and Advanced Virgo sensitivity curves during their first observing runs. The resulting data was analyzed by gravitational-wave detection algorithms and 6 of the waveforms were recovered with false alarm rates smaller than 1 in a thousand years. Parameter estimation algorithms were run on each of these waveforms to explore the ability to constrain the masses, component angular momenta and sky position of these waveforms. We also perform a large-scale monte-carlo study to assess the ability to recover each of the 60 hybrid waveforms with early Advanced LIGO and Advanced Virgo sensitivity curves. Our results predict that early Advanced LIGO and Advanced Virgo will have a volume-weighted average sensitive distance of 300Mpc (1Gpc) for $10M_{\\odot}+10M_{\\odot}$ ($50M_{\\odot}+50M_{\\odot}$) binary black hole coalescences. We demonstrate that neglecting the component angular momenta in the waveform models used in matched-filtering will result in a reduction in sensitivity for systems with large component angular momenta. [Abstract abridged for ArXiv, full version in PDF
Numerical modelling of sandstone uniaxial compression test using a mix-mode cohesive fracture model
Gui, Yilin; Kodikara, Jayantha
2015-01-01T23:59:59.000Z
A mix-mode cohesive fracture model considering tension, compression and shear material behaviour is presented, which has wide applications to geotechnical problems. The model considers both elastic and inelastic displacements. Inelastic displacement comprises fracture and plastic displacements. The norm of inelastic displacement is used to control the fracture behaviour. Meantime, a failure function describing the fracture strength is proposed. Using the internal programming FISH, the cohesive fracture model is programmed into a hybrid distinct element algorithm as encoded in Universal Distinct Element Code (UDEC). The model is verified through uniaxial tension and direct shear tests. The developed model is then applied to model the behaviour of a uniaxial compression test on Gosford sandstone. The modelling results indicate that the proposed cohesive fracture model is capable of simulating combined failure behaviour applicable to rock.
A phase screen model for simulating numerically the propagation of a laser beam in rain
Lukin, I P; Rychkov, D S; Falits, A V [Institute of Atmospheric Optics, Siberian Branch, Russian Academy of Sciences, Tomsk (Russian Federation); Lai, Kin S; Liu, Min R [DSO National Laboratories 20 (Singapore)
2009-09-30T23:59:59.000Z
The method based on the generalisation of the phase screen method for a continuous random medium is proposed for simulating numerically the propagation of laser radiation in a turbulent atmosphere with precipitation. In the phase screen model for a discrete component of a heterogeneous 'air-rain droplet' medium, the amplitude screen describing the scattering of an optical field by discrete particles of the medium is replaced by an equivalent phase screen with a spectrum of the correlation function of the effective dielectric constant fluctuations that is similar to the spectrum of a discrete scattering component - water droplets in air. The 'turbulent' phase screen is constructed on the basis of the Kolmogorov model, while the 'rain' screen model utiises the exponential distribution of the number of rain drops with respect to their radii as a function of the rain intensity. Theresults of the numerical simulation are compared with the known theoretical estimates for a large-scale discrete scattering medium. (propagation of laser radiation in matter)
Monitoring and Numerical Modeling of Shallow CO{sub 2} Injection, Greene County, Missouri
Rovey, Charles; Gouzie, Douglas; Biagioni, Richard
2013-09-30T23:59:59.000Z
The project titled Monitoring and Numerical Modeling of Shallow CO{sub 2} Injection, Greene County, Missouri provided training for three graduate students in areas related to carbon capture and storage. Numerical modeling of CO{sub 2} injection into the St. Francois aquifer at the Southwest Power Plant Site in Greene County, Missouri indicates that up to 4.1 x 10{sup 5} metric tons of CO{sub 2} per year could be injected for 30 years without exceeding a 3 MPa differential injection pressure. The injected CO{sub 2} would remain sequestered below the top of the overlying caprock (St. Francois confining unit) for more than 1000 years. Geochemical modeling indicates that portions of the injected CO{sub 2} will react rapidly with trace minerals in the aquifer to form various solid carbonate mineral phases. These minerals would store significant portions of injected CO{sub 2} over geologic time scales. Finally, a GIS data base on the pore-fluid chemistry of the overlying aquifer system in Missouri, the Ozark aquifer, was compiled from many sources. This data base could become useful in monitoring for leakage from future CO{sub 2} sequestration sites.
Stochastic models and numerical algorithms for a class of regulatory gene networks
Thomas Fournier; Jean-Pierre Gabriel; Christian Mazza; Jerome Pasquier; Jose Galbete; Nicolas Mermod
2008-10-01T23:59:59.000Z
Regulatory gene networks contain generic modules like those involving feedback loops, which are essential for the regulation of many biological functions. We consider a class of self-regulated genes which are the building blocks of many regulatory gene networks, and study the steady state distributions of the associated Gillespie algorithm by providing efficient numerical algorithms. We also study a regulatory gene network of interest in synthetic biology and in gene therapy, using mean-field models with time delays. Convergence of the related time-nonhomogeneous Markov chain is established for a class of linear catalytic networks with feedback loops
Numerical modeling of the high-temperature geothermal system of Amatitlan, Guatemala
Pham, M.; Menzies, A.J.; Sanyal, S.K. [GeothermEx, Inc., Richmond, CA (United States)] [and others
1996-12-31T23:59:59.000Z
Using the conceptual model of the Amatitlin geothermal system as a guide, a 14 by 10 km numerical model was set up with five layers and a total of 1,220 grid blocks. The number of grid blocks varies from layer to layer, with finer resolution in the two production layers. Boundary blocks for conductive heat transfer, heat and mass recharge, shallow discharge and deep regional fluid flow were added to the model as appropriate. With an inflow of hot fluid (336{degrees}C) at a rate of 450 tones per hour, an inflow of cool fluid (50{degrees}C) at a rate of 620 tons per hour and permeabilities varying from 1.5 to 50 md, the initial state temperatures on all five levels were well matched by the numerical model. The initial state results indicate that highest permeability is present to the southeast of the productive wells, suggesting that considerable additional production potential exists in this area. A two-phase region exists in the upper levels of the reservoir and in the outflow plume; the latter is consistent with the presence of fumaroles along the shoreline of Lago de Amatitlan. To further calibrate the numerical model, enthalpy and pressure data collected during tests of wells AMF-1 and -2 were matched, primarily by varying reservoir storage capacity. Good matches to the observed data were obtained for both wells, although sharp transients could not be matched. Modeling proceeded to the prediction stage under two scenarios: production of 12 and 25 MW (gross) for 30 years. 100% injection of separated water and condensate was assumed to take place in the vicinity of well AMF-3 in both cases. The results show an increase in enthalpy caused by expansion of the two-phase zone during the first year of production, followed by an enthalpy decline as reservoir liquid begins to migrate toward the production area. The two existing production wells (AMF-1 and -2) are predicted to be capable of maintaining output in the 12 MW case throughout the 30-year project life.
Numerical modeling of mixed sediment resuspension, transport, and deposition during the March 1998 sediment resuspension of mixed (cohesive plus noncohesive) sediment is developed and applied to quantitatively simulate the March 1998 resuspension events in southern Lake Michigan. Some characteristics
Lynett, Patrick
Tsunami inundation modeling in constructed environments: A physical and numerical comparison April 2013 Available online 17 May 2013 Keywords: Tsunami Inundation Macro-roughness Benchmark COULWAVE Friction factor A laboratory benchmark test for tsunami inundation through an urban waterfront including
Pan, Dongqing; Chien Jen, Tien [Department of Mechanical Engineering, University of Wisconsin-Milwaukee, Milwaukee, Wisconsin 53201 (United States); Li, Tao [School of Mechanical Engineering, Dalian University of Technology, Dalian 116024 (China); Yuan, Chris, E-mail: cyuan@uwm.edu [Department of Mechanical Engineering, University of Wisconsin-Milwaukee, 3200 North Cramer Street, Milwaukee, Wisconsin 53211 (United States)
2014-01-15T23:59:59.000Z
This paper characterizes the carrier gas flow in the atomic layer deposition (ALD) vacuum reactor by introducing Lattice Boltzmann Method (LBM) to the ALD simulation through a comparative study of two LBM models. Numerical models of gas flow are constructed and implemented in two-dimensional geometry based on lattice Bhatnagar–Gross–Krook (LBGK)-D2Q9 model and two-relaxation-time (TRT) model. Both incompressible and compressible scenarios are simulated and the two models are compared in the aspects of flow features, stability, and efficiency. Our simulation outcome reveals that, for our specific ALD vacuum reactor, TRT model generates better steady laminar flow features all over the domain with better stability and reliability than LBGK-D2Q9 model especially when considering the compressible effects of the gas flow. The LBM-TRT is verified indirectly by comparing the numerical result with conventional continuum-based computational fluid dynamics solvers, and it shows very good agreement with these conventional methods. The velocity field of carrier gas flow through ALD vacuum reactor was characterized by LBM-TRT model finally. The flow in ALD is in a laminar steady state with velocity concentrated at the corners and around the wafer. The effects of flow fields on precursor distributions, surface absorptions, and surface reactions are discussed in detail. Steady and evenly distributed velocity field contribute to higher precursor concentration near the wafer and relatively lower particle velocities help to achieve better surface adsorption and deposition. The ALD reactor geometry needs to be considered carefully if a steady and laminar flow field around the wafer and better surface deposition are desired.
Numerical modeling of roll structures in mesoscale vortexes over the Black Sea
Iarova, D A
2014-01-01T23:59:59.000Z
This paper is a case study of horizontal atmospheric rolls that formed over the Black Sea on 16 August 2007. The rolls were discovered in WRF modeling results for a mesoscale cyclone that originated over the sea on 15 August 2007. The roll formation mechanisms, such as Rayleigh-Benard convective instability, dynamic instability, advection and stretching of vertical velocity field inhomogeneities, are considered. It is shown that indeed convective instability played an important role in the roll formation but dynamic instability did not occur. In order to distinguish other possible mechanisms of the roll formation numerical experiments were performed. In these experiments sea surface temperature in the initial conditions was decreased in order to prevent convective instability. Even though convective instability was suppressed roll-like structures still appeared in the modeling results, although their height and circulation velocity were smaller than in the control run. It was found that these structures were ...
Temperature distributions in the laser-heated diamond anvil cell from 3-D numerical modeling
Rainey, E. S. G.; Kavner, A. [Department of Earth and Space Sciences, University of California, Los Angeles, California 90095 (United States); Hernlund, J. W. [Department of Earth and Planetary Science, University of California, Berkeley, California 94720 (United States); Earth-Life Science Institute, Megoro, Tokyo 152-8551 (Japan)
2013-11-28T23:59:59.000Z
We present TempDAC, a 3-D numerical model for calculating the steady-state temperature distribution for continuous wave laser-heated experiments in the diamond anvil cell. TempDAC solves the steady heat conduction equation in three dimensions over the sample chamber, gasket, and diamond anvils and includes material-, temperature-, and direction-dependent thermal conductivity, while allowing for flexible sample geometries, laser beam intensity profile, and laser absorption properties. The model has been validated against an axisymmetric analytic solution for the temperature distribution within a laser-heated sample. Example calculations illustrate the importance of considering heat flow in three dimensions for the laser-heated diamond anvil cell. In particular, we show that a “flat top” input laser beam profile does not lead to a more uniform temperature distribution or flatter temperature gradients than a wide Gaussian laser beam.
Seth A Veitzer
2008-10-21T23:59:59.000Z
Effects of stray electrons are a main factor limiting performance of many accelerators. Because heavy-ion fusion (HIF) accelerators will operate in regimes of higher current and with walls much closer to the beam than accelerators operating today, stray electrons might have a large, detrimental effect on the performance of an HIF accelerator. A primary source of stray electrons is electrons generated when halo ions strike the beam pipe walls. There is some research on these types of secondary electrons for the HIF community to draw upon, but this work is missing one crucial ingredient: the effect of grazing incidence. The overall goal of this project was to develop the numerical tools necessary to accurately model the effect of grazing incidence on the behavior of halo ions in a HIF accelerator, and further, to provide accurate models of heavy ion stopping powers with applications to ICF, WDM, and HEDP experiments.
Full-Scale Numerical Modeling of Turbulent Processes in the Earth's Ionosphere
Eliasson, B. [Department of Physics, Umeaa University, SE-901 87 Umeaa (Sweden); Stenflo, L. [Department of Physics, Umeaa University, SE-901 87 Umeaa (Sweden); Department of Physics, Linkoeping University, SE-581 83 Linkoeping (Sweden); Shukla, P. K. [Institut fuer Theoretische Physik IV, Fakultaet fuer Physik und Astronomie, Ruhr-Universitaet Bochum, D-44780 Bochum (Germany)
2008-10-15T23:59:59.000Z
We present a full-scale simulation study of ionospheric turbulence by means of a generalized Zakharov model based on the separation of variables into high-frequency and slow time scales. The model includes realistic length scales of the ionospheric profile and of the electromagnetic and electrostatic fields, and uses ionospheric plasma parameters relevant for high-latitude radio facilities such as Eiscat and HAARP. A nested grid numerical method has been developed to resolve the different length-scales, while avoiding severe restrictions on the time step. The simulation demonstrates the parametric decay of the ordinary mode into Langmuir and ion-acoustic waves, followed by a Langmuir wave collapse and short-scale caviton formation, as observed in ionospheric heating experiments.
place. 2. The basic model can be varied to suit foragers that optimise either their rate of net energy uptake or their foraging ef®ciency. 3. The model requires speci®cation of the time and energy budgetsAn economic model of the limits to foraging range in central place foragers with numerical
Two-dimensional numerical model of underground oil-shale retorting
Travis, B.J.; Hommert, P.J.; Tyner, C.E.
1983-01-01T23:59:59.000Z
A two-dimensional numerical model of underground oil shale retorting, which fully couples retorting chemistry with fluid and heat flow, has been developed. The model solves the time-dependent, two-dimensional mass, momentum, and energy balance equations for a nine-component fluid (O/sub 2/, N/sub 2/, H/sub 2/, CO/sub 2/, CO, CH/sub 4/, CH/sub x/, H/sub 2/O, and oil). Water and oil can flow in the liquid and/or vapor phases. Retort chemistry includes kerogen pyrolysis, carbonate decomposition, char reactions, and combustion. Also, detailed modeling of heat flow and chemistry inside shale particles allows large rubble sizes as well as small sizes to be considered. The model is compared to one-dimensional experimental data obtained from Lawrence Livermore National Laboratory. The model can be used to examine the effect of two-dimensional variations in shale grade, rubble size, permeability, porosity, geometry, inflow gas composition, etc. on retorting efficiency and process optimization. A sample calculation is presented.
Anooshehpoor, Rasool; Purvance, Matthew D.; Brune, James N.; Preston, Leiph A.; Anderson, John G.; Smith, Kenneth D.
2006-09-29T23:59:59.000Z
This report covers the following projects: Shake table tests of precarious rock methodology, field tests of precarious rocks at Yucca Mountain and comparison of the results with PSHA predictions, study of the coherence of the wave field in the ESF, and a limited survey of precarious rocks south of the proposed repository footprint. A series of shake table experiments have been carried out at the University of Nevada, Reno Large Scale Structures Laboratory. The bulk of the experiments involved scaling acceleration time histories (uniaxial forcing) from 0.1g to the point where the objects on the shake table overturned a specified number of times. The results of these experiments have been compared with numerical overturning predictions. Numerical predictions for toppling of large objects with simple contact conditions (e.g., I-beams with sharp basal edges) agree well with shake-table results. The numerical model slightly underpredicts the overturning of small rectangular blocks. It overpredicts the overturning PGA for asymmetric granite boulders with complex basal contact conditions. In general the results confirm the approximate predictions of previous studies. Field testing of several rocks at Yucca Mountain has approximately confirmed the preliminary results from previous studies, suggesting that he PSHA predictions are too high, possibly because the uncertainty in the mean of the attenuation relations. Study of the coherence of wavefields in the ESF has provided results which will be very important in design of the canisters distribution, in particular a preliminary estimate of the wavelengths at which the wavefields become incoherent. No evidence was found for extreme focusing by lens-like inhomogeneities. A limited survey for precarious rocks confirmed that they extend south of the repository, and one of these has been field tested.
Stetiu, C.
1993-07-01T23:59:59.000Z
A thermal building simulation program is a numerical model that calculates the response of the building envelopes to weather and human activity, simulates dynamic heating and cooling loads, and heating and cooling distribution systems, and models building equipment operation. The scope of the research is to supply the users of such programs with information about the dangers and benefits of simplifying the input to their models. The Introduction describes the advantages of modeling the heat transfer mechanisms in a building. The programs that perform this type of modeling have, however, limitations. The user is therefore often put in the situation of simplifying the floor plans of the building under study, but not being able to check the effects that this approximation introduces in the results of the simulation. Chapter 1 is a description of methods. It also introduces the floor plans for the office building under study and the ``reasonable`` floor plans simplifications. Chapter 2 presents DOE-2, the thermal building simulation program used in the sensitivity study. The evaluation of the accuracy of the DOE-2 program itself is also presented. Chapter 3 contains the sensitivity study. The complicated nature of the process of interpreting the temperature profile inside a space leads to the necessity of defining different building modes. The study compares the results from the model of the detailed building description with the results from the models of the same building having simplified floor plans. The conclusion is reached that a study of the effects of simplifying the floor plans of a building is important mainly for defining the cases in which this approximation is acceptable. Different results are obtained for different air conditioning/load regimes of the building. 9 refs., 24 figs.
Cheng, C. L.; Gragg, M. J.; Perfect, E.; White, Mark D.; Lemiszki, P. J.; McKay, L. D.
2013-08-24T23:59:59.000Z
Numerical simulations are widely used in feasibility studies for geologic carbon sequestration. Accurate estimates of petrophysical parameters are needed as inputs for these simulations. However, relatively few experimental values are available for CO2-brine systems. Hence, a sensitivity analysis was performed using the STOMP numerical code for supercritical CO2 injected into a model confined deep saline aquifer. The intrinsic permeability, porosity, pore compressibility, and capillary pressure-saturation/relative permeability parameters (residual liquid saturation, residual gas saturation, and van Genuchten alpha and m values) were varied independently. Their influence on CO2 injection rates and costs were determined and the parameters were ranked based on normalized coefficients of variation. The simulations resulted in differences of up to tens of millions of dollars over the life of the project (i.e., the time taken to inject 10.8 million metric tons of CO2). The two most influential parameters were the intrinsic permeability and the van Genuchten m value. Two other parameters, the residual gas saturation and the residual liquid saturation, ranked above the porosity. These results highlight the need for accurate estimates of capillary pressure-saturation/relative permeability parameters for geologic carbon sequestration simulations in addition to measurements of porosity and intrinsic permeability.
Numerical Modeling of Thermal EOR: Comprehensive Coupling of an AMR-Based Model
Paris-Sud XI, Université de
Flow and Geomechanics N. Guy*, G. Enchéry and G. Renard IFP Energies nouvelles, 1-4 avenue de Bois of Thermal EOR: Comprehensive Coupling of an AMR-Based Model of Thermal Fluid Flow and Geomechanics when both thermal fluid flow and geomechanics are coupled in order to take into account variations
A numerical procedure to model and monitor CO2 sequestration in
Santos, Juan
Oceanografia e di Geofisica Sperimentale (OGS), Trieste, ITALY IC-MSQUARE 2012, September 3-7, 2012A numerical
Tidal Downsizing model. I. Numerical methods: saving giant planets from tidal disruptions
Nayakshin, Sergei
2014-01-01T23:59:59.000Z
Tidal Downsizing (TD) is a recently developed planet formation theory that supplements the classical Gravitational disc Instability (GI) model with planet migration inward and tidal disruptions of GI fragments in the inner regions of the disc. Numerical methods for a detailed population synthesis of TD planets are presented here. As an example application, the conditions under which GI fragments collapse faster than they migrate into the inner $a\\sim$ few AU disc are considered. It is found that most gas fragments are tidally or thermally disrupted unless (a) their opacity is $\\sim 3$ orders of magnitude less than the interstellar dust opacity at metallicities typical of the observed giant planets, or (b) the opacity is high but the fragments accrete large dust grains (pebbles) from the disc. Case (a) models produce very low mass solid cores ($M_{\\rm core} < 0.1$ Earth masses) and follow a negative correlation of giant planet frequency with host star metallicity. In contrast, case (b) models produce massiv...
Khangaonkar, Tarang; Yang, Zhaoqing; Kim, Tae Yun; Roberts, Mindy
2011-07-20T23:59:59.000Z
Through extensive field data collection and analysis efforts conducted since the 1950s, researchers have established an understanding of the characteristic features of circulation in Puget Sound. The pattern ranges from the classic fjordal behavior in some basins, with shallow brackish outflow and compensating inflow immediately below, to the typical two-layer flow observed in many partially mixed estuaries with saline inflow at depth. An attempt at reproducing this behavior by fitting an analytical formulation to past data is presented, followed by the application of a three-dimensional circulation and transport numerical model. The analytical treatment helped identify key physical processes and parameters, but quickly reconfirmed that response is complex and would require site-specific parameterization to include effects of sills and interconnected basins. The numerical model of Puget Sound, developed using unstructured-grid finite volume method, allowed resolution of the sub-basin geometric features, including presence of major islands, and site-specific strong advective vertical mixing created by bathymetry and multiple sills. The model was calibrated using available recent short-term oceanographic time series data sets from different parts of the Puget Sound basin. The results are compared against (1) recent velocity and salinity data collected in Puget Sound from 2006 and (2) a composite data set from previously analyzed historical records, mostly from the 1970s. The results highlight the ability of the model to reproduce velocity and salinity profile characteristics, their variations among Puget Sound subbasins, and tidally averaged circulation. Sensitivity of residual circulation to variations in freshwater inflow and resulting salinity gradient in fjordal sub-basins of Puget Sound is examined.
Boyer, Edmond
Introduction Flattening the Earth Continuation procedure Flat Earth Numerical simulations Continuation from a flat to a round Earth model in the coplanar orbit transfer problem M. Cerf1, T. Haberkorn, SADCO 2011, March 2nd M. Cerf, T. Haberkorn, E. Tr´elat Continuation from a flat to a round Earth model
Beaumont, Christopher
Numerical Modeling of Salt Tectonics on Passive Continental Margins: Preliminary Assessment Sciences The University of Leeds LS2 9JT Leeds United Kingdom Abstract Salt tectonics in passive model of frictional-plastic sedimentary overburden overlying a linear viscous salt layer. We present
Paris-Sud XI, Université de
Stochastic reduced-order model for an automotive vehicle in presence of numerous local elastic a high modal density in the low-frequency range, such as an automotive vehicle. This type of structure is applied on a complex computational model of an automotive vehicle. 1 INTRODUCTION This work is performed
Paris-Sud XI, Université de
1 Numerical modelling of hybrid arc/laser welding: a Level Set approach for weld bead formation.Bellet@mines-paristech.fr ABSTRACT The joining of high thickness steel sheets by means of hybrid Laser/GMAW welding processes of the workpiece borders. Two finite elements models are presented to illustrate: (i) A hybrid arc/laser welding
An efficient numerical model for incompressible two-phase flow in fractured media Hussein Hoteit a,1
Firoozabadi, Abbas
in fractured hydrocarbon reservoirs [16]. In this model, the matrixfracture mass transfer is describedAn efficient numerical model for incompressible two-phase flow in fractured media Hussein Hoteit a,1 , Abbas Firoozabadi a,b,* a Reservoir Engineering Research Institute, Palo Alto, CA, USA b Yale
A numerical model of aerosol scavenging. Part 2, Simulation of a large city fire
Bradley, M.M.; Molenkamp, C.R.
1991-10-01T23:59:59.000Z
Using a three-dimensional numerical cloud/smoke-plume model, we have simulated the burning of a large, mid-latitude city following a nuclear exchange. The model includes 18 dynamic and microphysical equations that predict the fire-driven airflow, cloud processes, and smoke-cloud interactions. In the simulation, the intense heating from the burning city produces a firestorm with updraft velocities exceeding 60 m/s. Within 15 minutes of ignition, the smoke plume penetrates the tropopause. The updraft triggers a cumulonimbus cloud that produces significant quantities of ice, snow, and hail. These solid hydrometeors, as well as cloud droplets and rain, interact with the smoke particles from the fire. At the end of the one-hour simulation, over 20% of the smoke is in slowly falling snowflakes. If the snow reaches the ground before the flakes completely sublimate (or melt and then evaporate), then only approximately 50% of the smoke will survive the scavenging processes and remain in the atmosphere to affect the global climate.
Eckert, Andreas
2013-05-31T23:59:59.000Z
In this project generic anticline structures have been used for numerical modeling analyses to study the influence of geometrical parameters, fluid flow boundary conditions, in situ stress regime and inter-bedding friction coefficient on geomechanical risks such as fracture reactivation and fracture generation. The resulting stress states for these structures are also used to determine safe drilling directions and a methodology for wellbore trajection optimization is developed that is applicable for non-Andersonian stress states. The results of the fluid flow simulation show that the type of fluid flow boundary condition is of utmost importance and has significant impact on all injection related parameters. It is recommended that further research is conducted to establish a method to quantify the fluid flow boundary conditions for injection applications. The results of the geomechanical simulation show that in situ stress regime is a crucial, if not the most important, factor determining geomechanical risks. For extension and strike slip stress regimes anticline structures should be favored over horizontally layered basin as they feature higher ?P{sub c} magnitudes. If sedimentary basins are tectonically relaxed and their state of stress is characterized by the uni-axial strain model the basin is in exact frictional equilibrium and fluids should not be injected. The results also show that low inter bedding friction coefficients effectively decouple layers resulting in lower ?P{sub c} magnitudes, especially for the compressional stress regime.
Hans-Dieter Alber
2015-05-20T23:59:59.000Z
The accurate simulation of phase interfaces in solids requires small model error and small numerical error. If a phase field model is used and the interface carries low interface energy, then the model error is only small if the interface width in the model is chosen small. Yet, for effective numerical computation the interface width should be large. Choosing the parameters, which determine the width, is therefore an optimality problem. We study this problem for the Allen-Cahn equation coupled to the elasticity equations by constructing an asymptotic solution of second order, which yields an expansion for the kinetic relation of the model. This expansion determines the choice of the parameters, however only if the difference between the expansion and the exact kinetic relation is uniformly small with respect to a second parameter controlling the interface energy. To show this uniformity we determine the asymptotics with respect to this second parameter by scaling of the model equations. Our investigations are formal.
NUMERICAL MODELING OF THE COAGULATION AND POROSITY EVOLUTION OF DUST AGGREGATES
Okuzumi, Satoshi; Sakagami, Masa-aki [Graduate School of Human and Environmental Studies, Kyoto University, Yoshida-nihonmatsu-cho, Sakyo-ku, Kyoto 606-8501 (Japan); Tanaka, Hidekazu, E-mail: satoshi.okuzumi@ax2.ecs.kyoto-u.ac.j [Institute of Low Temperature Science, Hokkaido University, Sapporo 060-0819 (Japan)
2009-12-20T23:59:59.000Z
Porosity evolution of dust aggregates is crucial in understanding dust evolution in protoplanetary disks. In this study, we present useful tools to study the coagulation and porosity evolution of dust aggregates. First, we present a new numerical method for simulating dust coagulation and porosity evolution as an extension of the conventional Smoluchowski equation. This method follows the evolution of the mean porosity for each aggregate mass simultaneously with the evolution of the mass distribution function. This method reproduces the results of previous Monte Carlo simulations with much less computational expense. Second, we propose a new collision model for porous dust aggregates on the basis of our N-body experiments on aggregate collisions. As the first step, we focus on 'hit-and-stick' collisions, which involve neither compression nor fragmentation of aggregates. We first obtain empirical data on porosity changes between the classical limits of ballistic cluster-cluster and particle-cluster aggregation. Using the data, we construct a recipe for the porosity change due to general hit-and-stick collisions as well as formulae for the aerodynamical and collisional cross sections. Our collision model is thus more realistic than a previous model of Ormel et al. based on the classical aggregation limits only. Simple coagulation simulations using the extended Smoluchowski method show that our collision model explains the fractal dimensions of porous aggregates observed in a full N-body simulation and a laboratory experiment. By contrast, similar simulations using the collision model of Ormel et al. result in much less porous aggregates, meaning that this model underestimates the porosity increase upon unequal-sized collisions. Besides, we discover that aggregates at the high-mass end of the distribution can have a considerably small aerodynamical cross section per unit mass compared with aggregates of lower masses. This occurs when aggregates drift under uniform acceleration (e.g., gravity) and their collision is induced by the difference in their terminal velocities. We point out an important implication of this discovery for dust growth in protoplanetary disks.
Menzies, A.J.; Granados, E.E.; Sanyal, .K.; Merida-I., L.; Caicedo-A, A.
1991-01-01T23:59:59.000Z
A significant amount of geoscientific and reservoir engineering data have been collected from the Zunil geothermal field since 1973. The data have been used to define a conceptual model for the field which has formed the basis for the construction of a three dimensional numerical simulation model. The numerical model has successfully matched both the initial state of the reservoir, as indicated by subsurface temperature and pressure distributions within the presently drilled area, and available well test data. The well test data include short and long term discharge tests and a comprehensive pressure interference test. Calibration of the model will continue during 1991 when the results from drilling and testing of three additional deep wells are available. The model will then be used to study various long term production scenarios for the proposed 15 MW power development.
Kirby, James T.
Numerical simulation of the 2011 Tohoku tsunami: Comparison with field observations and sensitivity history, created a major tsunami that caused numerous deaths and enormous destruction on the nearby Hon- shu coast. Various tsunami sources were developed for this event, based on inverting seismic or GPS
The VOLNA code for the numerical modelling of tsunami waves: generation, propagation and
of the finite volume scheme implemented in the code. We explain the numerical treatment of the wet is decided on the base of inundation maps which are produced with this type of numerical tools. Finally we and the perspectives for future research presented. Key words: tsunami waves, shallow water equations, tsunami
Harrison, Mark
1988-01-01T23:59:59.000Z
IMA Journal of Numerical Analysis (1988) 8, 415-433 Finite Element Methods for a Model for Full are given and then a discrete-time, explicit finite element procedure is defined and analysed, with finite on the existence, uniqueness and finite element approximation of the solution of Biot's equations were given in [15
Occhipinti, Giovanni "Ninto"
LETTER Earth Planets Space, 63, 847851, 2011 Three-dimensional numerical modeling of tsunami, 2011; Accepted June 30, 2011; Online published September 27, 2011) The tremendous tsunami following, to reproduce the tsunami signature observed in the airglow by the imager located in Hawaii and clearly showing
RodrÃguez, Rodolfo
Numerical analysis of a finite element method for the axisymmetric eddy current model, 27002, Lugo, Spain The aim of this paper is to analyze a finite element method to solve an eddy current of the method are reported. Keywords: low-frequency harmonic Maxwell equations, eddy current problems, finite
Demouchy, Sylvie
Numerical modelling of erosion processes in the Himalayas of Nepal: effects of spatial variations of the morphology of mountain belts. Here we investigate the modalities of defor- mation in Central Nepal on a c that the pattern of uplift in Nepal is mainly dependent on both erodability and fault geometry, rather than
Dadi, Sireesh Kumar
2011-10-21T23:59:59.000Z
-Butler model, which obtained a new drainage function based on a linearized Richard’s equation but limited the variation of soil moisture and hydraulic conductivity in the unsaturated zone to exponential functions. Numerical analysis was conducted with VS2DT...
Tandon, Amit
Segmentation and Tracking of Mesoscale Eddies in Numeric Ocean Models Vishal Sood, Bin John suggested that the mesoscale eddies and mesoscale features play a strong role in carrying heat poleward oceanographers an invaluable tool to assess mesoscale eddies and the Lagrangian characteristics of this mesoscale
Baudouin, Lucie
A controlled distributed parameter model for a fluid-flexible structure system: numerical consider the problem of active reduction of vibrations in a fluid-flexible structure system and the sloshing of the fuel inside the wing's tank. The control is performed using piezoelectric patches
Sim, Yoon Sub; Kim, Eui Kwang; Eoh, Jae Hyuk [Korea Atomic Energy Research Institute (Korea, Republic of)
2005-06-15T23:59:59.000Z
To overcome the drawbacks of conventional schemes for a numerical analysis of a steam generator (SG), an efficient numerical model has been developed to analyze the steady state of a once-through-type SG where the feedwater is heated to superheated steam. In the developed model, the temperature and enthalpy are defined at the boundary of a calculation cell, and the exact solutions for the temperature distribution in a calculation cell are utilized. This feature of the developed model frees calculation from the undesirable effects of numerical diffusion, and only a small number of nodes are required. Also, the developed model removes the ambiguity from the parameter values at the inlet and exit of a calculation.The BoSupSG-SS computer code was developed by using the analysis model, and it performed well with only three calculation nodes to analyze a superheated SG. The developed model can be effectively used for the cases where a fast one-dimensional calculation is required such as an SG or system design analysis.
Numerical study of the mass spectrum in the 2D O(3) sigma model with a theta term
B. Alles; A. Papa
2007-11-12T23:59:59.000Z
It has been conjectured that the mass spectrum of the O(3) non-linear sigma model with a theta term in 2 dimensions may possess an excited state, which decays when theta is lowered from pi below a critical value. Since the direct numerical investigation of the model is prevented by a sign problem, we try to infer some information on the mass spectrum at real theta by studying the model at imaginary theta via analytic continuation. A modified Swendsen-Wang cluster algorithm has been introduced to simulate the model with the theta term.
Vršnak, B.; Žic, T.; Dumbovi?, M. [Hvar Observatory, Faculty of Geodesy, University of Zagreb, Ka??eva 26, HR-10000 Zagreb (Croatia); Temmer, M.; Möstl, C.; Veronig, A. M. [Kanzelhöhe Observatory—IGAM, Institute of Physics, University of Graz, Universittsplatz 5, A-8010 Graz (Austria); Taktakishvili, A.; Mays, M. L. [NASA Goddard Space Flight Center, Greenbelt, MD 20771 (United States); Odstr?il, D., E-mail: bvrsnak@geof.hr, E-mail: tzic@geof.hr, E-mail: mdumbovic@geof.hr, E-mail: manuela.temmer@uni-graz.at, E-mail: christian.moestl@uni-graz.at, E-mail: astrid.veronig@uni-graz.at, E-mail: aleksandre.taktakishvili-1@nasa.gov, E-mail: m.leila.mays@nasa.gov, E-mail: dusan.odstrcil@nasa.gov [George Mason University, Fairfax, VA 22030 (United States)
2014-08-01T23:59:59.000Z
Real-time forecasting of the arrival of coronal mass ejections (CMEs) at Earth, based on remote solar observations, is one of the central issues of space-weather research. In this paper, we compare arrival-time predictions calculated applying the numerical ''WSA-ENLIL+Cone model'' and the analytical ''drag-based model'' (DBM). Both models use coronagraphic observations of CMEs as input data, thus providing an early space-weather forecast two to four days before the arrival of the disturbance at the Earth, depending on the CME speed. It is shown that both methods give very similar results if the drag parameter ? = 0.1 is used in DBM in combination with a background solar-wind speed of w = 400 km s{sup –1}. For this combination, the mean value of the difference between arrival times calculated by ENLIL and DBM is ?-bar =0.09±9.0 hr with an average of the absolute-value differences of |?|-bar =7.1 hr. Comparing the observed arrivals (O) with the calculated ones (C) for ENLIL gives O – C = –0.3 ± 16.9 hr and, analogously, O – C = +1.1 ± 19.1 hr for DBM. Applying ? = 0.2 with w = 450 km s{sup –1} in DBM, one finds O – C = –1.7 ± 18.3 hr, with an average of the absolute-value differences of 14.8 hr, which is similar to that for ENLIL, 14.1 hr. Finally, we demonstrate that the prediction accuracy significantly degrades with increasing solar activity.
Mathiesen, Patrick; Collier, Craig; Kleissl, Jan
2013-01-01T23:59:59.000Z
cycle: The RUC. Monthly Weather Review. 132, 495?518. th Conference on Numerical Weather Prediction. American closure schemes. Monthly Weather Review. 122, 927?945.
Two-dimensional numerical models of open-top hydrothermal convection at high Rayleigh and Nusselt
Wilcock, William
for Nu up to 5060. Solutions are characterized by an unstable bottom thermal boundary layer where equation. To avoid classical numerical artifacts such as nonphysical oscillatory behavior and artificial
An evaluation of pocket-model, numerical readout breath alcohol testing instruments
Van Tassel, William Edward
2004-11-15T23:59:59.000Z
Eight small-scale breath alcohol measurement devices were tested for accuracy, precision and the ability to not yield false positive and false negative readings. These pocket-sized breath testers (PMBTs), which provided numerical readout of Br...
Qualification of the ITER CS Quench Detection System using Numerical Modeling
Martovetsky, Nicolai N [ORNL; Radovinsky, Alexey L [ORNL
2013-01-01T23:59:59.000Z
Abstract The ITER Central Solenoid (CS) magnet needs to be protected against overheating of the conductor in the event of the occurrence of a normal zone (NZ). Due to a large amount of stored energy and slow NZ propagation, the NZ needs to be detected and the switchyard needs to open the breakers within 2 s after detection of the NZ. The CS will be discharged on a dump resistor with a time constant of 7.5 s. During operation of the CS and its interaction with the poloidal field (PF) coils and plasma current, the CS experiences large inductive voltages from multiple sources, including nonlinear signals from eddy currents in the vacuum vessel and plasma current variation, that make the task of detecting the resistive signal even more difficult. This inductive voltage needs to be cancelled by quench detection (QD) hardware (e.g., bridges, converters, filters, processors) and appropriate processing of the QD signals to reliably detect NZ initiation and propagation. Two redundant schemes are proposed as the baseline for the CS QD System: 1) A scheme with Regular Voltage Taps (RVT) from triads of Double Pancakes (DP) supplemented by Central Difference Averaging (CDA) and by digital suppression of the inductive voltage from all active coils (the CS and PF coils). Voltage taps are taken from helium outlets at the CS outer diameter. 2)A scheme with Cowound Voltage Taps (CVT) taken from cowound wires routed from the helium inlet at the CS inner diameter. Summary of results of the numerical modeling of the performance of both baseline CS QD systems is presented in this paper. Index Terms Quench detection, Central Solenoid, ITER
A new strategy for discrete element numerical models: 2. Sandbox applications
Sandiford, Mike
and visualized by the modeler. Reliable modeling of geomechanic processes enables the structural interpreter
Draxl, C.; Churchfield, M.; Mirocha, J.; Lee, S.; Lundquist, J.; Michalakes, J.; Moriarty, P.; Purkayastha, A.; Sprague, M.; Vanderwende, B.
2014-06-01T23:59:59.000Z
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.
Imada, Shinsuke; Watanabe, Tetsuya
2015-01-01T23:59:59.000Z
We have studied the chromospheric evaporation flow during the impulsive phase of the flare by using the Hinode/EIS observation and 1D hydrodynamic numerical simulation coupled to the time-dependent ionization. The observation clearly shows that the strong redshift can be observed at the base of the flaring loop only during the impulsive phase. We performed two different numerical simulations to reproduce the strong downflows in FeXII and FeXV during the impulsive phase. By changing the thermal conduction coefficient, we carried out the numerical calculation of chromospheric evaporation in the thermal conduction dominant regime (conductivity coefficient kappa0 = classical value) and the enthalpy flux dominant regime (kappa0 = 0.1 x classical value). The chromospheric evaporation calculation in the enthalpy flux dominant regime could reproduce the strong redshift at the base of the flare during the impulsive phase. This result might indicate that the thermal conduction can be strongly suppressed in some cases o...
Kenis, Paul J. A.
A three-dimensional numerical model of a micro laminar flow fuel cell with a bridge: Membraneless fuel cell Laminar flow fuel cell Numerical model Convection-diffusion equations Electrode kinetics equations COMSOL a b s t r a c t The operation of a laminar flow fuel cell (LFFC) involves complex interplay
Bradley, M.M.; Peterson, K.R.; Rodriguez, D.J.
1988-11-17T23:59:59.000Z
During the hours immediately following a nuclear exchange, large fires could inject enormous quantities of smoke into the atmosphere. This smoke, together with dust from surface bursts, would severely restrict visibilities and darken the skies over large areas for days. This, in turn, could impact surface and air operations and systems. These effects could be mitigated by various scavenging mechanisms within the convective clouds that form above the fires. In order to evaluate impacts of post-nuclear-exchange smoke injection, we are developing a three-dimensional numerical smoke plume model (OCTET) to simulate the dynamics and microphysical processes within smoke plumes and convective clouds above large fires. This model is based on the dynamic framework of the Klemp-Wilhelmson (1978) convective storm model and includes parameterizations of scavenging processes. In addition, we are combining results of laboratory research, field experiments, and detailed numerical modeling of cloud microphysical processes in order to better understand smoke scavenging mechanisms. In this brief demonstration of capabilities, we present results from the OCTET model and from a three-dimensional mesoscale model. The smoke plume and fire-induced cloud simulations demonstrate the effects of nucleation scavenging, seasonal variation of atmospheric stability, and various fuel sources. The mesoscale simulations (that use the plume model output as input) demonstrate the mesoscale transport and diffusion of smoke and predict optical depths over the hypothetical target area. No dust effects have been included in these simulations. 2 refs., 8 figs.
Kiryukhin, A.V.; Sugrobov, V.M.
1986-01-21T23:59:59.000Z
The application of the two-dimensional numerical heat-transfer model to the Pauzhetka hydrothermal system allowed us to establish that: (1) a shallow magma body with the anomalous temperature of 700-1000 C and with a volume of 20-30 km{sup 3} may be a heat source for the formation of the Pauzhetka hydrothermal system. (2) The water feeding source of the Pauzhetka hydrothermal system may be meteoric waters which are infiltrated at an average rate of 5-10 kg/s {center_dot} km{sup 2}. The coupling of the numerical heat-transfer model with hydroisotopic data (D,T,{sup 18}O) obtained from the results of testing of exploitation wells, rivers and springs is the basis to understand more clearly the position of recharge areas and the structure of water flows in the hydrothermal system.
A numerical procedure to model and monitor CO2 sequestration in
Santos, Juan
area of research. o We present a methodology integrating numerical simulation of CO2 -brine o The simultaneous flow of brine and CO2 is described by the well-known Black in the brine but the brine is not allowed to vaporize into the CO2 phase. o This formulation uses
LABORATORY-NUMERICAL MODEL COMPARISONS OF CANYON FLOWS: A PARAMETER STUDY.
, but the enhanced viscosities needed to obtain numerical stability give boundary layers that are too wide along length scales, one the fluid depth and another a more narrow boundary-layer-like thickness [O(RoBu-1 is the interaction of an oscillatory, along-slope background current with an isolated canyon incised in an otherwise
Numerical and analytical modelling of entropy noise in a supersonic nozzle with a shock
Boyer, Edmond
and numerical assessments of the indirect noise generated through a nozzle are presented. The configuration. R¨ohle, F. Thiele, B. Noll, The Entropy Wave Generator (EWG): A reference case on entropy noise, Journal of Sound and Vibration 326 (2009) 574-598] where an entropy wave is generated upstream of a nozzle
McConnell, Joshua B
2000-01-01T23:59:59.000Z
to produce a temperature profile in the liner thickness. An analytical stress model, using the results of the derived equations and the numerical thermal model, was constructed to determine the magnitude of the stresses the liner is subjected to after...
Scalar-Scalar Ladder Model in the Unequal-Mass Case. III - Numerical Studies of the P-Wave Case -
Ichio Fukui; Noriaki Setoh
1999-01-21T23:59:59.000Z
The eigenvalue problem for the p-wave bound states formed by two unequal-mass scalar particles through the massive scalar particle exchange is analyzed numerically in the framework of the Bethe-Salpeter ladder model. As in the s-wave case, the eigenvalues of the coupling constant are found to become complex for some mass configurations in some range of the bound state mass. The Bethe-Salpeter amplitudes of the low-lying bound states are also investigated.
De Castro, Carlos Armando
2011-01-01T23:59:59.000Z
In this paper is developed a simple mathematical model of transient heat transfer under soil with plastic mulch in order to determine with numerical studies the influence of different plastic mulches on the soil temperature and the evolutions of temperatures at different depths with time. The governing differential equations are solved by a Galerkin Finite Element Model, taking into account the nonlinearities due to radiative heat exchange between the soil surface, the plastic mulch and the atmosphere. The model was validated experimentally giving good approximation of the model to the measured data. Simulations were run with the validated model in order to determine the optimal combination of mulch optical properties to maximize the soil temperature with a Taguchi's analysis, proving that the material most used nowadays in Colombia is not the optimal and giving quantitative results of the properties the optimal mulch must possess.
Field observations and numerical model experiments for the snowmelt process at a field site
Ohara, N; Kawas, M L
2006-01-01T23:59:59.000Z
a one dimensional heat transfer model is introduced in orderone dimensional heat transfer model is con- venient forHowever, the simple heat transfer model is a useful tool to
THE STATE OF THE ART OF NUMERICAL MODELING OF THERMOHYDROLOGIC FLOW IN FRACTURED ROCK MASSES
Wang, J.S.Y.
2013-01-01T23:59:59.000Z
disposal of radioactive waste: The Sandia waste isolation flow and transport (SWIFT) model: Sandia Laboratories Report
NUMERICAL MODELING FOR THE FORMATION MECHANISM OF 3D TOPOGRAPHY ON MICROBIAL MAT SURFACES
Patel, Harsh Jay
2013-09-27T23:59:59.000Z
determined. A deterministic model includes both a dynamic model, and a static model. Compared to a deterministic solute diffusion model, a stochastic diffusion model has one or more stochastic elements. Concentration diffusion system formulated... solute diffusion equation, which merely appears in the form of a number line, the quasi-2D solute diffusion equation forms a Cartesian grid system. Also, for the explicit quasi-2D solute diffusion equation (Eq. 4), concentration variation is measured...
Numerical Modeling At Dixie Valley Geothermal Area (McKenna & Blackwell,
AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on DeliciousPlasmaPLawrence County,1980) | OpenAl., 2001) |NewNorthNorwalk,Novomer Jump to:NuPowerNumerical
Numerical Modeling At Lightning Dock Geothermal Area (O'Brien, Et Al.,
AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on DeliciousPlasmaPLawrence County,1980) | OpenAl., 2001) |NewNorthNorwalk,Novomer Jump to:NuPowerNumerical1984) | Open
Numerical Modeling At Neal Hot Springs Geothermal Area (U.S. Geothermal
AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on DeliciousPlasmaPLawrence County,1980) | OpenAl., 2001) |NewNorthNorwalk,Novomer Jump to:NuPowerNumerical1984) |
A numerical model for the coupled long-term evolution of salt marshes and tidal flats
Fagherazzi, Sergio
-shore mudflat model that takes into account tidal effects; Waeles et al. [2004] incor- porated in the same
Molenkamp, C.R.; Grossman, A.
1999-12-20T23:59:59.000Z
A network of small balloon-borne transponders which gather very high resolution wind and temperature data for use by modern numerical weather predication models has been proposed to improve the reliability of long-range weather forecasts. The global distribution of an array of such transponders is simulated using LLNL's atmospheric parcel transport model (GRANTOUR) with winds supplied by two different general circulation models. An initial study used winds from CCM3 with a horizontal resolution of about 3 degrees in latitude and longitude, and a second study used winds from NOGAPS with a 0.75 degree horizontal resolution. Results from both simulations show that reasonable global coverage can be attained by releasing balloons from an appropriate set of launch sites.
Grossman, A; Molenkamp, C R
1999-08-25T23:59:59.000Z
A proposal has been made to establish a high density global network of atmospheric micro transponders to record time, temperature, and wind data with time resolution of {le} 1 minute, temperature accuracy of {+-} 1 K, spatial resolution no poorer than {approx}3km horizontally and {approx}0.1km vertically, and 2-D speed accuracy of {le} 1m/s. This data will be used in conjunction with advanced numerical weather prediction models to provide increases in the reliability of long range weather forecasts. Major advances in data collection technology will be required to provide the proposed high-resolution data collection network. Systems studies must be undertaken to determine insertion requirements, spacing, and evolution of the transponder ensemble, which will be used to collect the data. Numerical models which provide realistic global weather pattern simulations must be utilized in order to perform these studies. A global circulation model with a 3{sup o} horizontal resolution has been used for initial simulations of the generation and evolution of transponder distributions. These studies indicate that reasonable global coverage of transponders can be achieved by a launch scenario consisting of the sequential launch of transponders at specified heights from a globally distributed set of launch sites.
Not Available
1980-09-01T23:59:59.000Z
Nineteen papers were presented at the workshop on modeling thermohydrologic flow in fractured masses. This workshop was a result of the interest currently being given to the isolation of nuclear wastes in geologic formations. Included in these proceedings are eighteen of the presentations, one abstract and summaries of the panel discussions. The papers are listed under the following categories: introduction; overviews; fracture modelings; repository studies; geothermal models; and recent developments. Eighteen of the papers have been abstracted and indexed.
/ free-slip / free-slip Numerical models of slab deformation in the upper mantle
Cerveny, Vlastislav
of the plate. Model: Model domain: Results: Effect of decoupling mechanism 410 km 660 km oceánska doska nadlozná doska Results: Effect of boundary condition Effect of yield stress Plate velocity Kinematic, crust.s 1020 Pa.s 1021 Pa.s Free-slip, crust = 1020 Pas, y = 108 Pa Kinematic, crust = 1020 Pas, y = 109 Pa
A simplified numerical model of coronal energy dissipation based on reduced MHD
E. Buchlin; V. Aletti; S. Galtier; M. Velli; G. Einaudi; J. -C. Vial
2003-03-27T23:59:59.000Z
A 3D model intermediate between cellular automata (CA) models and the reduced magnetohydrodynamic (RMHD) equations is presented to simulate solar impulsive events generated along a coronal magnetic loop. The model consists of a set of planes distributed along a magnetic loop between which the information propagates through Alfven waves. Statistical properties in terms of power-laws for energies and durations of dissipative events are obtained, and their agreement with X-ray and UV flares observations is discussed. The existence of observational biases is also discussed.
Numerical Modeling of Fractured Shale-Gas and Tight-Gas Reservoirs Using Unstructured Grids
Olorode, Olufemi Morounfopefoluwa
2012-02-14T23:59:59.000Z
Various models featuring horizontal wells with multiple induced fractures have been proposed to characterize flow behavior over time in tight gas and shale gas systems. Currently, there is little consensus regarding the effects of non...
A Numerical Model For The Dynamics Of Pyroclastic Flows At Galeras...
model is applied to one of the slopes of Galeras Volcano to show: (1) the temperature evolution with the time; (2) dynamic pressure change; and (3) particle concentration along the...
Numerical Modeling of Fractured Shale-Gas and Tight-Gas Reservoirs Using Unstructured Grids
Olorode, Olufemi Morounfopefoluwa
2012-02-14T23:59:59.000Z
Various models featuring horizontal wells with multiple induced fractures have been proposed to characterize flow behavior over time in tight gas and shale gas systems. Currently, there is little consensus regarding the effects of non...
Numerical and analytical modeling of heat transfer between fluid and fractured rocks
Li, Wei, S.M. Massachusetts Institute of Technology
2014-01-01T23:59:59.000Z
Modeling of heat transfer between fluid and fractured rocks is of particular importance for energy extraction analysis in EGS, and therefore represents a critical component of EGS design and performance evaluation. In ...
Numerical Modeling of the Transient Thermal Interference of Vertical U-Tube Haet Exchangers
Muraya, Norman K.
1994-01-01T23:59:59.000Z
Non-linear finite element models were developed to simulate transient heat and mass transfer in the soil surrounding the ground heat exchangers of ground-coupled heat pumps (GCHPs) operating in the cooling mode. Parametric ...
Numerical modeling of alongshore sediment transport and shoreline change along the Galveston coast
Sitanggang, Khairil Irfan
2005-02-17T23:59:59.000Z
An alongshore sediment transport and shoreline change analysis on Galveston Island in the period of 1990-2001 is conducted in this study using the Generalized Model for Simulating Shoreline Change (GENESIS). The study is divided into three main...
NUMERICAL MODELING FOR THE FORMATION MECHANISM OF 3D TOPOGRAPHY ON MICROBIAL MAT SURFACES
Patel, Harsh Jay
2013-09-27T23:59:59.000Z
, an innovative modeling approach was employed that focuses on the interface growth of the microbial mat surfaces using a combined stochastic and deterministic approach. A range of different initial conditions were simulated to evaluate the 3D topography evolution...
Airborne observations and numerical modeling of fetch- limited waves in the Gulf of Tehuantepec
Romero, Leonel
2008-01-01T23:59:59.000Z
velocity scaling in wind wave generation. Boundary-Layerlinear theory of of wind wave generation applied to waveSource terms in a third-generation wind wave model. J. Phys.
Thermodynamic Modeling and Numerical Simulation of Single-Shaft Microturbine Performance
Hao, X.; Zhang, G.; Zhou, J.; Chen, Y.
2006-01-01T23:59:59.000Z
the energy saving potential and economical benefits of microturbine and its combined production system, a simple mathematical model of microturbine is proposed. Part-load characteristics of main components are also considered for analyzing the unit...
THE STATE OF THE ART OF NUMERICAL MODELING OF THERMOHYDROLOGIC FLOW IN FRACTURED ROCK MASSES
Wang, J.S.Y.
2013-01-01T23:59:59.000Z
studies related to the Cerro Prieto Field: Proceedings, 1stSymposium on the Cerro Prieto Geothermal Field, San Diego,modeling studies of the Cerro Prieto Reservoir--A progress
Lusk, Miriam Beatriz
2011-12-31T23:59:59.000Z
To capture the fracture process and non-linear behavior at the element and structural level of granular materials, concrete, under the presence of pre-existing imperfections, a constitutive model and a mesh free method is ...
Jia, Dongxing 1984-
2012-11-15T23:59:59.000Z
The Unidirectional Hybrid Wave Model (UHWM) predicts irregular wave kinematics and pressure accurately in comparison with its linear counterpart and modification, especially near the free surface. Hence, in using the ...
Numerical modeling of fluid flow and time-lapse seismics to monitor ...
santos
May 30, 2014 ... including the presence of shale seals and fractures and fractal variations of the ... In the Black-Oil model employed, brine is NOT present, OIL is.
Final Report: A Model Management System for Numerical Simulations of Subsurface Processes
Zachmann, David
2013-10-07T23:59:59.000Z
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.
and free surface models and a global heat transfer model, with moving boundaries. An axisymmetric fluid to determine flow field, after the phase boundaries have been determined, by the heat transfer model. A finite field, from which temperature gradients are determined. The heat transfer model is furthermore expanded
Synthesis of Numerical Methods for Modeling Wave Energy Converter-Point Absorbers: Preprint
Li, Y.; Yu, Y. H.
2012-05-01T23:59:59.000Z
During the past few decades, wave energy has received significant attention among all ocean energy formats. Industry has proposed hundreds of prototypes such as an oscillating water column, a point absorber, an overtopping system, and a bottom-hinged system. In particular, many researchers have focused on modeling the floating-point absorber as the technology to extract wave energy. Several modeling methods have been used such as the analytical method, the boundary-integral equation method, the Navier-Stokes equations method, and the empirical method. However, no standardized method has been decided. To assist the development of wave energy conversion technologies, this report reviews the methods for modeling the floating-point absorber.
None
2005-07-01T23:59:59.000Z
This work was carried out to understand the behavior of the solid and gas phases in a CFB riser. Only the riser is modeled as a straight pipe. A model with linear algebraic approximation to solids viscosity of the form, {musubs} = 5.34{epsisubs}, ({espisubs} is the solids volume fraction) with an appropriate boundary condition at the wall obtained by approximate momentum balance solution at the wall to acount for the solids recirculation is tested against experimental results. The work done was to predict the flow patterns in the CFB risers from available experimental data, including data from a 7.5-cm-ID CFB riser at the Illinois Institute of Technology and data from a 20.0-cm-ID CFB riser at the Particulate Solid Research, Inc., facility. This research aims at modeling the removal of hydrogen sulfide from hot coal gas using zinc oxide as the sorbent in a circulating fluidized bed and in the process indentifying the parameters that affect the performance of the sulfidation reactor. Two different gas-solid reaction models, the unreacted shrinking core (USC) and the grain model were applied to take into account chemical reaction resistances. Also two different approaches were used to affect the hydrodynamics of the process streams. The first model takes into account the effect of micro-scale particle clustering by adjusting the gas-particle drag law and the second one assumes a turbulent core with pseudo-steady state boundary condition at the wall. A comparison is made with experimental results.
Modeling and numerical techniques for high-speed digital simulation of nuclear power plants
Wulff, W.; Cheng, H.S.; Mallen, A.N.
1987-01-01T23:59:59.000Z
Conventional computing methods are contrasted with newly developed high-speed and low-cost computing techniques for simulating normal and accidental transients in nuclear power plants. Six principles are formulated for cost-effective high-fidelity simulation with emphasis on modeling of transient two-phase flow coolant dynamics in nuclear reactors. Available computing architectures are characterized. It is shown that the combination of the newly developed modeling and computing principles with the use of existing special-purpose peripheral processors is capable of achieving low-cost and high-speed simulation with high-fidelity and outstanding user convenience, suitable for detailed reactor plant response analyses.
Numerical Modeling of the Transient Thermal Interference of Vertical U-Tube Haet Exchangers
Muraya, Norman K.
Significance of Terms . . . . . . . . . . . . . . . . . . . . . . . . Computation of Matrices . Time-stepping Scheme Theory . Summary . 37 39 42 49 52 52 52 53 IV VALIDATION OF HEAT TRANSFER MODELS , 54 IV.1 Overall U-value , 54 IV.2 Single Source Heating.......................................... 69 Xlll CHAPTER Page V APPLICATIONS OF THE HEAT TRANSFER MODELS 0 0 0 0 0 0 0 o. iO Vol Equivalent Diameter 0 0 0 0 0 0 0 0 0 0 • 0 0 • 0 • 0 • 0 • 0 • 0 0 0 0 0 0 0 o' iO V.2 Short Circuiting and Equivalent Diameter .. 0 • 0 • 0 ••• 0 0 • 0 • 0 i6 V...
A nanoscale numerical model of calcium silicate hydrate P.C. Fonseca a,
Andrade, Jose
to pre- dict the bulk properties of cement and concrete, such as shrinkage, creep, permeability and the properties of structural concrete are not fully under- stood. Models are becoming increasingly important, and cracking. CSH is responsible for much of the cohesive proper- ties in concrete but the chemical origin
Experimental Validation of a Numerical Multizone Airflow and Pollutant Transfer Model
Paris-Sud XI, Université de
and long-term assessment of the performances of ventilation systems, the experimental house MARIA and ventilation systems are modeled in MATLAB/Simulink environment. This paper quickly describes the multi exhaust, balanced and natural ventilation systems. In addition, the virtual laboratory SIMBAD Building
Model Validation and Spatial Interpolation by Combining Observations with Outputs from Numerical
Washington at Seattle, University of
""r,c,rn The authors are for hel]JfuI #12;Abstract Constructing maps of pollution levels is vital for air quality concentrations. Key tlJords: air pollution, Ba~yesian inference, change of support, likelihood approaches, Matern Resolutions 2.5 Modeling a Nonstationary Covariance . 3 Estimation 3.1 Algorithm 4 Application: Air Pollution
Roberts, Malcolm
analysis, and high-performance computing. I have experience working in the areas of compu- tational fluid mechanics, subgrid models for turbulence, and flows in complex geometries. 1 Introduction Consider, which is the product of the characteristic length and velocity of the flow divided by the kinematic
A numerical procedure to model and monitor CO2 sequestration in ...
santos
Sep 7, 2012 ... The model of the formation is based on the porosity and clay .... The porosity and dry bulk and shear modulus of the formation were obtained using a pore ... The viscosity, density and bulk modulus of CO2 needed for the flow ...
Numerical Modeling of Nonlinear Coupling between Lines/Beams with Multiple Floating Bodies
Yang, Chan K.
2010-07-14T23:59:59.000Z
to be capable of modeling the tendon disconnection both at the top and the bottom connection as well as the down stroke behavior for the pinned bottom joint. The performance of the tie-down clamp of derrick is also investigated by using six degrees of freedom...
Gotika, Priyanka
2012-02-14T23:59:59.000Z
-SMOOTH DYNAMICAL SYSTEMS : : : : : : : : : : : : : : : : : : : : : 7 A. Classical simple harmonic oscillator . . . . . . . . . . . . . 7 B. Spring-mass-dash-pot systems . . . . . . . . . . . . . . . . 7 C. Coulomb friction model... . . . . . . . . . 22 B. Lumped parameter system with a Bingham-type dash- pot of non-linear monotonic characteristics . . . . . . . . . 24 C. Lumped parameter system with a Bingham-type dash- pot of non-linear non-monotonic characteristics . . . . . . 27 D. Lumped...
Observational and Numerical Modeling Studies of Turbulence on the Texas-Louisiana Continental Shelf
Zhang, Zheng
2013-05-24T23:59:59.000Z
calculated from velocity measurements in the bot- tom boundary layer (BBL), using the Kolmogorov’s -5/3 law in the inertial subrange of energy spectra of vertical velocity fluctuations in each burst measurement. Four second-moment turbulence closure models...
A multiple layer numerical model of the formation of the low-level jet
Shen, Tsu-Cheng
1980-01-01T23:59:59.000Z
310 320 e(kj Fig. 3. The initial vertical distribution of potential tempera- ture in the model . Table 3. The vertical distribution of the variables used for setting the initial conditions. LEVEL k (m s ) kh(m s ) u(ms ) v(ms ) w(ms ) e (k) 00 p...
Numerical study of three-dimensional PIC for the surface plasmon excitation based on Drude model
Liu, La-Qun; Wang, Hui-Hui; Liu, Da-Gang
2015-01-01T23:59:59.000Z
This paper explores the time-domain equations of noble metals, in which Drude model is adopted to describe the dielectric constant, to implement three-dimensional particle-in-cell (PIC) simulations for the surface plasmon excitation with the finite-difference time-domain method (FDTD). A three-dimensional model for an electron bunch movement near the metal film is constructed, and particle-in-cell (PIC) simulations are carried out with various metal films of different thicknesses. The frequency of surface plasmon obtained from PIC simulation is agreed with that from theory. Furthermore, the surface plasmon wave properties of excitation and propagation with the metal film is summarized by PIC results.
Laboratory experiments and numerical modeling of wave attenuation through artificial vegetation
Augustin, Lauren Nicole
2009-05-15T23:59:59.000Z
). The first hydrodynamic model developed by Price et al. (1968) simulated the effects of seaweed as a high viscous layer. Mork (1996) extended the idea of the high viscous layer and developed a theory for kelp plants that took into account not only viscous... has been validated by artificial laboratory kelp experiments for the species Laminaria Hyperborea, and is assumed appropriate for representing wave transformation and damping over submerged vegetative fields of variable depths. Mork (1996) studied...
Exact theory and numeric results for short pulse ionization of simple model atom in one dimension
Rokhlenko, Alexander
2015-01-01T23:59:59.000Z
Our exact theory for continuous harmonic perturbation of a one dimensional model atom by parametric variations of its potential is generalized for the cases when a) the atom is exposed to short pulses of an external harmonic electric field and b) the forcing is represented by short bursts of different shape changing the strength of the binding potential. This work is motivated not only by the wide use of laser pulses for atomic ionization, but also by our earlier study of the same model which successfully described the ionization dynamics in all orders, i.e. the multi-photon processes, though being treated by the non-relativistic Schr\\"odinger equation. In particular, it was shown that the bound atom cannot survive the excitation of its potential caused by any non-zero frequency and amplitude of the continuous harmonic forcing. Our present analysis found important laws of the atomic ionization by short pulses, in particular the efficiency of ionizing this model system and presumably real ones as well.
Dubus, Guillaume; Fromang, Sébastien
2015-01-01T23:59:59.000Z
Detailed modeling of the high-energy emission from gamma-ray binaries has been propounded as a path to pulsar wind physics. Fulfilling this ambition requires a coherent model of the flow and its emission in the region where the pulsar wind interacts with the stellar wind of its companion. We developed a code that follows the evolution and emission of electrons in the shocked pulsar wind based on inputs from a relativistic hydrodynamical simulation. The code is used to model the well-documented spectral energy distribution and orbital modulations from LS 5039. The pulsar wind is fully confined by a bow shock and a back shock. The particles are distributed into a narrow Maxwellian, emitting mostly GeV photons, and a power law radiating very efficiently over a broad energy range from X-rays to TeV gamma rays. Most of the emission arises from the apex of the bow shock. Doppler boosting shapes the X-ray and VHE lightcurves, constraining the system inclination to $i\\approx 35^{\\rm o}$. There is a tension between th...
Volker Sick; Dennis N. Assanis
2002-11-27T23:59:59.000Z
Through the combination of advanced imaging laser diagnostics with multi-dimensional computer models, a new understanding of the performance of direct-injection gasoline engines is pursuit. The work focuses on the fuel injection process, the breakup of the liquid into a fine spray and the mixing of the fuel with the in-cylinder gases. Non-intrusive laser diagnostics will be used to measure the spatial distribution of droplets and vaporized fuel with very high temporal resolution. These data along with temperature measurements will be used to validate a new spray breakup model for gasoline direct-injection. Experimental data on near wall fuel distributions will be used for comparison with a model that predicts the spray-wall interaction and the dynamics of the liquid film on the surface. Quantitative measurements of local nitric oxide concentrations inside the combustion chamber will provide a critical test for a numerical simulation of the nitric oxide formation process. This model is based on a modified flamelet approach and will be used to study the effects of exhaust gas recirculation.
Numerical modeling of two-phase behavior in the PEFC gas diffusion layer
Mukherjee, Partha Pa223876 [Los Alamos National Laboratory; Kang, Qinjun [Los Alamos National Laboratory; Mukundan, Rangachary [Los Alamos National Laboratory; Borup, Rod L [Los Alamos National Laboratory
2009-01-01T23:59:59.000Z
A critical performance limitation in the polymer electrolye fuel cell (PEFC) is attributed to the mass transport loss originating from suboptimal liquid water transport and flooding phenomena. Liquid water can block the porous pathways in the fibrous gas diffusion layer (GDL) and the catalyst layer (CL), thus hindering oxygen transport from the flow field to the electrochemically actives sites in the catalyst layer. In this paper, the study of the two phase behavior and the durability implications due to the wetting characteristics in the carbon paper GDL are presented using a pore-scale modeling framework.
Development and validation of a vertically two-dimensional mesoscale numerical model
Walters, Michael Kent
1985-01-01T23:59:59.000Z
) ds ' ? gHpls1 h s + gHps g the first law of thermodynamics, and the boundary condition s = 0 at the top of the model. The derived equations used are the following: ~ = ? V Vp ? s W + ~ (~) a d as Tp dt 1 D ? V p( Q ) C T s = --(~) 1 7dt 1 Q... are evenly spaced at an interval of 5 x/2, while the vertical grid points are evenly spaced at an increment of Q s/2. . K4&H&&%& . W&H&&& ~ ~ ~ ~ ~ ~ ~ ~ I 12 4 x is equal to 45 km, while 6 s is equivalent to 250 m in the simulations...
Numerical modeling of boundary layer flow under shoaling and breaking waves
Pattipawaej, Olga Catherina
1998-01-01T23:59:59.000Z
velocities and shear stresses in the bottom boundary layer under breaking waves. The experiment wss conducted in the Precision Wave Tank located in the Ocean Engineering Laboratory at the University of Delaware. The flume was 33 m long, 0. 6 m wide, and 1... modeling of surf zone hydrodynamics. " Rep. No. CACR-95-97, Center for Applied Coastal Research, University Delaware, Newark, Delaware. Cox, D. T. , Kobayashi, N. , and Okayasu, A. (1996). "Bottom shear stress in the surf zone. " J. Geophys. Res. , 101(C...
Mathematical and numerical models to achieve high speed with special-purpose parallel processors
Cheng, H.S.; Wulff, W.; Mallen, A.N.
1986-07-01T23:59:59.000Z
One simulation facility that has been developed is the BNL Plant Analyzer, currently set up for BWR plant simulations at up to seven times faster than real-time process speeds. The principal hardware components of the BNL Plant Analyzer are two units of special-purpose parallel processors, the AD10 of Applied Dynamics International and a PDP-11/34 host computer. The AD10 is specifically designed for time-critical system simulations, utilizing the modern parallel processing technology with pipeline architecture. The simulator employs advanced modeling techniques and efficient integration techniques in conjunction with the parallel processors to achieve high speed performance.
Parallel 3D Finite Element Numerical Modelling of DC Electron Guns
Prudencio, E.; Candel, A.; Ge, L.; Kabel, A.; Ko, K.; Lee, L.; Li, Z.; Ng, C.; Schussman, G.; /SLAC
2008-02-04T23:59:59.000Z
In this paper we present Gun3P, a parallel 3D finite element application that the Advanced Computations Department at the Stanford Linear Accelerator Center is developing for the analysis of beam formation in DC guns and beam transport in klystrons. Gun3P is targeted specially to complex geometries that cannot be described by 2D models and cannot be easily handled by finite difference discretizations. Its parallel capability allows simulations with more accuracy and less processing time than packages currently available. We present simulation results for the L-band Sheet Beam Klystron DC gun, in which case Gun3P is able to reduce simulation time from days to some hours.
Numerical study of a binary Yukawa model in regimes characteristic of globular proteins in solutions
Giacometti, Achille; Gazzillo, Domenico; Pastore, Giorgio; Das, Tushar Kanti [Istituto Nazionale per la Fisica della Materia and Dipartimento di Chimica Fisica, Universita di Venezia, S. Marta DD 2137, I-30123 Venice (Italy); Dipartimento di Fisica Teorica, Universita di Trieste and INFM-DEMOCRITOS, National Simulation Center, Strada Costiera 11, Miramare, I-34100 Trieste (Italy); ICTP, Diploma Course, Strada Costiera 11, Miramare P.O Box 586, I-34100 Trieste (Italy)
2005-03-01T23:59:59.000Z
The main goal of this paper is to assess the limits of validity, in the regime of low concentration and strong Coulomb coupling (high molecular charges), of a simple perturbative approximation to the radial distribution functions (RDF's), based upon a low-density expansion of the potential of mean force and proposed to describe protein-protein interactions in a recent small-angle-scattering (SAS) experimental study. A highly simplified Yukawa (screened Coulomb) model of monomers and dimers of a charged globular protein ({beta}-lactoglobulin) in solution is considered. We test the accuracy of the RDF approximation, as a necessary complementary part of the previous experimental investigation, by comparison with the fluid structure predicted by approximate integral equations and exact Monte Carlo (MC) simulations. In the MC calculations, an Ewald construction for Yukawa potentials has been used to take into account the long-range part of the interactions in the weakly screened cases. Our results confirm that the perturbative first-order approximation is valid for this system even at strong Coulomb coupling, provided that the screening is not too weak (i.e., for Debye length smaller than monomer radius). A comparison of the MC results with integral equation calculations shows that both the hypernetted-chain (HNC) and Percus-Yevick closures have a satisfactory behavior under these regimes, with the HNC being superior throughout. The relevance of our findings for interpreting SAS results is also discussed.
Chiswell, S
2009-01-11T23:59:59.000Z
Assimilation of radar velocity and precipitation fields into high-resolution model simulations can improve precipitation forecasts with decreased 'spin-up' time and improve short-term simulation of boundary layer winds (Benjamin, 2004 & 2007; Xiao, 2008) which is critical to improving plume transport forecasts. Accurate description of wind and turbulence fields is essential to useful atmospheric transport and dispersion results, and any improvement in the accuracy of these fields will make consequence assessment more valuable during both routine operation as well as potential emergency situations. During 2008, the United States National Weather Service (NWS) radars implemented a significant upgrade which increased the real-time level II data resolution to 8 times their previous 'legacy' resolution, from 1 km range gate and 1.0 degree azimuthal resolution to 'super resolution' 250 m range gate and 0.5 degree azimuthal resolution (Fig 1). These radar observations provide reflectivity, velocity and returned power spectra measurements at a range of up to 300 km (460 km for reflectivity) at a frequency of 4-5 minutes and yield up to 13.5 million point observations per level in super-resolution mode. The migration of National Weather Service (NWS) WSR-88D radars to super resolution is expected to improve warning lead times by detecting small scale features sooner with increased reliability; however, current operational mesoscale model domains utilize grid spacing several times larger than the legacy data resolution, and therefore the added resolution of radar data is not fully exploited. The assimilation of super resolution reflectivity and velocity data into high resolution numerical weather model forecasts where grid spacing is comparable to the radar data resolution is investigated here to determine the impact of the improved data resolution on model predictions.
Rockhold, Mark L.; White, Mark D.; Freeman, Eugene J.
2004-10-12T23:59:59.000Z
This letter report documents initial numerical analyses conducted by PNNL to provide support for a feasibility study on decommissioning of the canyon buildings at Hanford. The 221-U facility is the first of the major canyon buildings to be decommissioned. The specific objective of this modeling effort was to provide estimates of potential rates of migration of residual contaminants out of the 221-U facility during the first 40 years after decommissioning. If minimal contaminant migration is predicted to occur from the facility during this time period, then the structure may be deemed to provide a level of groundwater protection that is essentially equivalent to the liner and leachate collection systems that are required at conventional landfills. The STOMP code was used to simulate transport of selected radionuclides out of a canyon building, representative of the 221-U facility after decommissioning, for a period of 40 years. Simulation results indicate that none of the selected radionuclides that were modeled migrated beyond the concrete structure of the facility during the 40-year period of interest. Jacques (2001) identified other potential contaminants in the 221-U facility that were not modeled, however, including kerosene, phenol, and various metals. Modeling of these contaminants was beyond the scope of this preliminary effort due to increased complexity. Simulation results indicate that contaminant release from the canyon buildings will be diffusion controlled at early times. Advection is expected to become much more important at later times, after contaminants have diffused out of the facility and into the surrounding soil environment. After contaminants have diffused out of the facility, surface infiltration covers will become very important for mitigating further transport of contaminants in the underlying vadose zone and groundwater.
Mickens, R.E.
1997-12-12T23:59:59.000Z
The major thrust of this proposal was to continue our investigations of so-called non-standard finite-difference schemes as formulated by other authors. These schemes do not follow the standard rules used to model continuous differential equations by discrete difference equations. The two major aspects of this procedure consist of generalizing the definition of the discrete derivative and using a nonlocal model (on the computational grid or lattice) for nonlinear terms that may occur in the differential equations. Our aim was to investigate the construction of nonstandard finite-difference schemes for several classes of ordinary and partial differential equations. These equations are simple enough to be tractable, yet, have enough complexity to be both mathematically and scientifically interesting. It should be noted that all of these equations differential equations model some physical phenomena under an appropriate set of experimental conditions. The major goal of the project was to better understand the process of constructing finite-difference models for differential equations. In particular, it demonstrates the value of using nonstandard finite-difference procedures. A secondary goal was to construct and study a variety of analytical techniques that can be used to investigate the mathematical properties of the obtained difference equations. These mathematical procedures are of interest in their own right and should be a valuable contribution to the mathematics research literature in difference equations. All of the results obtained from the research done under this project have been published in the relevant research/technical journals or submitted for publication. Our expectation is that these results will lead to improved finite difference schemes for the numerical integration of both ordinary and partial differential equations. Section G of the Appendix gives a concise summary of the major results obtained under funding by the grant.
11. NUMERICAL TECHNIQUES 1 Numerical identification of effective multipole
Boyer, Edmond
11. NUMERICAL TECHNIQUES 1 Numerical identification of effective multipole moments of polarizable of the induced multipole moments. A general multipole theory is available in the literature, however, only linear multipole model is usually exploited when determining numerically these effective moments. Since this axial
Mohaghegh, Shahab
with a large number of producers, second, to a CO2 sequestration project in Australia, and finally to a numerical simulation study of potential carbon storage site in the United States. The numerical reservoir
Luo, Meng
Advanced High Strength Steels (AHSS) are increasingly used in automotive industry due to their superior strength and substantial weight advantage. However, their compromised ductility gives rise to numerous manufacturing ...
Houze, Jr., Robert A. [University of Washington Dept. of Atmospheric Sciences
2013-11-13T23:59:59.000Z
We examined cloud radar data in monsoon climates, using cloud radars at Darwin in the Australian monsoon, on a ship in the Bay of Bengal in the South Asian monsoon, and at Niamey in the West African monsoon. We followed on with a more in-depth study of the continental MCSs over West Africa. We investigated whether the West African anvil clouds connected with squall line MCSs passing over the Niamey ARM site could be simulated in a numerical model by comparing the observed anvil clouds to anvil structures generated by the Weather Research and Forecasting (WRF) mesoscale model at high resolution using six different ice-phase microphysical schemes. We carried out further simulations with a cloud-resolving model forced by sounding network budgets over the Niamey region and over the northern Australian region. We have devoted some of the effort of this project to examining how well satellite data can determine the global breadth of the anvil cloud measurements obtained at the ARM ground sites. We next considered whether satellite data could be objectively analyzed to so that their large global measurement sets can be systematically related to the ARM measurements. Further differences were detailed between the land and ocean MCS anvil clouds by examining the interior structure of the anvils with the satellite-detected the CloudSat Cloud Profiling Radar (CPR). The satellite survey of anvil clouds in the Indo-Pacific region was continued to determine the role of MCSs in producing the cloud pattern associated with the MJO.
Prickett, T.A.
1980-04-01T23:59:59.000Z
Specifications are given which are necessary to develop a three-dimensional numerical model capable of simulating regional mass transport of radionuclides from a deep waste repository. The model to be developed will include all of the significant mass transport processes including flow, chemical, and thermal advection, mechanical dispersion, molecular diffusion, ion exchange reactions, and radioactive decay. The model specifications also include that density and viscosity fluid properties be functions of pressure, temperature, and concentration and take into account fluid and geologic heterogenieties by allowing possible assignment of individual values to every block of the model. The model specifications furthermore include the repository shape, input/output information, boundary conditions, and the need for documentation and a user's manual. Model code validation can be accomplished with the included known analytical or laboratory solutions. It is recommended that an existing finite-difference model (developed by INTERCOMP and INTERA, Inc.) be used as a starting point either as an acceptable basic code for modification or as a pattern for the development of a completely different numerical scheme. A ten-step plan is given to outline the general procedure for development of the code.
Guerin, Gilles
migrations. The comparison of successive 3D surveys indicates changes in seismic attributes within reservoirs: Bottom Hole Temperatures (BHT) and 3D seismic interpretation. More than 600 BHT from 200 wells allowed us from the main seismic reflectors interpreted from 3D surveys, we perform a 3D numerical simulation
of formation heterogeneity on convective mixing. Finite volume based numerical code was developed to capture of Canadian Petroleum Technology, 44(10). - Neufeld, J. A., Hesse M. A., Riaz, A., Hallworth M. A., Tchelepi, H. A., and Huppert H. E., 2010, Convective dissolution of carbon dioxide in saline aquifers, Geophys
Kirby, James T.
Palma, Canary Islands): Tsunami source and near field effects S. M. Abadie,1 J. C. Harris,2 S. T. Grilli of the Cumbre Vieja Volcano (CVV; La Palma, Canary Island, Spain) through numerical simulations performed in two of such wave trains on La Palma and other Canary Islands are assessed in detail in the paper. Citation: Abadie
Tao, Wei-Kuo; Houze, Robert, A., Jr.; Zeng, Xiping
2013-03-14T23:59:59.000Z
This three-year project, in cooperation with Professor Bob Houze at University of Washington, has been successfully finished as planned. Both ARM (the Atmospheric Radiation Measurement Program) data and cloud-resolving model (CRM) simulations were used to identify the water budgets of clouds observed in two international field campaigns. The research results achieved shed light on several key processes of clouds in climate change (or general circulation models), which are summarized below. 1. Revealed the effect of mineral dust on mesoscale convective systems (MCSs) Two international field campaigns near a desert and a tropical coast provided unique data to drive and evaluate CRM simulations, which are TWP-ICE (the Tropical Warm Pool International Cloud Experiment) and AMMA (the African Monsoon Multidisciplinary Analysis). Studies of the two campaign data were contrasted, revealing that much mineral dust can bring about large MCSs via ice nucleation and clouds. This result was reported as a PI presentation in the 3rd ASR Science Team meeting held in Arlington, Virginia in March 2012. A paper on the studies was published in the Journal of the Atmospheric Sciences (Zeng et al. 2013). 2. Identified the effect of convective downdrafts on ice crystal concentration Using the large-scale forcing data from TWP-ICE, ARM-SGP (the Southern Great Plains) and other field campaigns, Goddard CRM simulations were carried out in comparison with radar and satellite observations. The comparison between model and observations revealed that convective downdrafts could increase ice crystal concentration by up to three or four orders, which is a key to quantitatively represent the indirect effects of ice nuclei, a kind of aerosol, on clouds and radiation in the Tropics. This result was published in the Journal of the Atmospheric Sciences (Zeng et al. 2011) and summarized in the DOE/ASR Research Highlights Summaries (see http://www.arm.gov/science/highlights/RMjY5/view). 3. Used radar observations to evaluate model simulations In cooperation with Profs. Bob Houze at University of Washington and Steven Rutledge at Colorado State University, numerical model results were evaluated with observations from W- and C-band radars and CloudSat/TRMM satellites. These studies exhibited some shortcomings of current numerical models, such as too little of thin anvil clouds, directing the future improvement of cloud microphysics parameterization in CRMs. Two papers of Powell et al (2012) and Zeng et al. (2013), summarizing these studies, were published in the Journal of the Atmospheric Sciences. 4. Analyzed the water budgets of MCSs Using ARM data from TWP-ICE, ARM-SGP and other field campaigns, the Goddard CRM simulations were carried out to analyze the water budgets of clouds from TWP-ICE and AMMA. The simulations generated a set of datasets on clouds and radiation, which are available http://cloud.gsfc.nasa.gov/. The cloud datasets were available for modelers and other researchers aiming to improve the representation of cloud processes in multi-scale modeling frameworks, GCMs and climate models. Special datasets, such as 3D cloud distributions every six minutes for TWP-ICE, were requested and generated for ARM/ASR investigators. Data server records show that 86,206 datasets were downloaded by 120 users between April of 2010 and January of 2012. 5. MMF simulations The Goddard MMF (multi-scale modeling framework) has been improved by coupling with the Goddard Land Information System (LIS) and the Goddard Earth Observing System Model, Version 5 (GOES5). It has also been optimized on NASA HEC supercomputers and can be run over 4000 CPUs. The improved MMF with high horizontal resolution (1 x 1 degree) is currently being applied to cases covering 2005 and 2006. The results show that the spatial distribution pattern of precipitation rate is well simulated by the MMF through comparisons with satellite retrievals from the CMOPRH and GPCP data sets. In addition, the MMF results were compared with three reanalyses (MERRA, ERA-Interim and CFSR). Although the MMF tends
Sengupta, Arkaprabha
2010-01-01T23:59:59.000Z
Phase transformation in Nitinol polycrystals 3.1 Martensiticvector pairs (m ? , b ? ) for Nitinol . . . . . . . . . . .4 Constitutive model of Nitinol 4.1 Constitutive models for
D'Agnese, F.A.; O'Brien, G.M.; Faunt, C.C.; Belcher, W.R.; San Juan, Carma
2002-11-22T23:59:59.000Z
In the early 1990's, two numerical models of the Death Valley regional ground-water flow system were developed by the U.S. Department of Energy. In general, the two models were based on the same basic hydrogeologic data set. In 1998, the U.S. Department of Energy requested that the U.S. Geological Survey develop and maintain a ground-water flow model of the Death Valley region in support of U.S. Department of Energy programs at the Nevada Test Site. The purpose of developing this ''second-generation'' regional model was to enhance the knowledge and understanding of the ground-water flow system as new information and tools are developed. The U.S. Geological Survey also was encouraged by the U.S. Department of Energy to cooperate to the fullest extent with other Federal, State, and local entities in the region to take advantage of the benefits of their knowledge and expertise. The short-term objective of the Death Valley regional ground-water flow system project was to develop a steady-stat e representation of the predevelopment conditions of the ground-water flow system utilizing the two geologic interpretations used to develop the previous numerical models. The long-term objective of this project was to construct and calibrate a transient model that simulates the ground-water conditions of the study area over the historical record that utilizes a newly interpreted hydrogeologic conceptual model. This report describes the result of the predevelopment steady-state model construction and calibration.
Yu-Shu Wu
2006-02-28T23:59:59.000Z
A three-dimensional site-scale numerical model has been developed to simulate water and gas flow, heat transfer, and radionuclide transport in the unsaturated zone of Yucca Mountain, Nevada, the American underground repository site for high level radioactive waste. The modeling approach is based on a mathematical formulation of coupled multiphase fluid and heat flow and tracer transport through porous and fractured rock. This model is intended for use in predicting current and future conditions in the unsaturated zone, so as to aid in assessing the system performance of the repository. In particular, an integrated modeling methodology is discussed for integrating a wide variety of moisture, pneumatic, thermal, and isotopic geochemical data into comprehensive modeling analyses. The reliability and accuracy of the model predictions were the subject of a comprehensive model calibration study, in which the model was calibrated against measured data, including liquid saturation, water potential, and temperature. This study indicates that the model is able to reproduce the overall system behavior at Yucca Mountain with respect to moisture profiles, pneumatic pressure and chloride concentration variations in different geological units, and ambient geothermal conditions.
Dadi, Sireesh Kumar
2011-10-21T23:59:59.000Z
: the Neuman model, which assumes instantaneous drainage from the unsaturated zone; the Moench model, which considered gradual drainage from the unsaturated zone using a series of exponential terms in the water table boundary condition; and the Mathias...
Schuh, Harald
with a dynamic Earth system model (completed in October 2004) Variations of Earth rotation are caused Earth system model DyMEG has been developed. It is based on the balance of angular momentum
Cirpka, Olaf Arie
. It takes into account the two phases CO2 and brine and the components CO2 and water which can dissolve the conceptual model for a non-isothermal composi- tional CO2-water (brine) model based on the simulator MUFTE. MODEL For the description of the flow and transport processes of carbon dioxide and brine in a rock
Numerical studies of a one-dimensional three-spin spin-glass model with long-range interactions
Larson, Derek; Katzgraber, Helmut G.; Moore, M. A.; Young, A. P.
2010-01-01T23:59:59.000Z
We study a p-spin spin-glass model to understand if the finite-temperature glass transition found in the mean-field regime of p-spin models, and used to model the behavior of structural glasses, persists in the nonmean-field ...
The Numerical Simulation of Turbulence
W. Schmidt
2007-12-06T23:59:59.000Z
In this contribution, I give an overview of the various approaches toward the numerical modelling of turbulence, particularly, in the interstellar medium. The discussion is placed in a physical context, i. e. computational problems are motivated from basic physical considerations. Presenting selected examples for solutions to these problems, I introduce the basic ideas of the most commonly used numerical methods.
NUMERICAL ANALYSIS KENDALL E. ATKINSON
Atkinson, Kendall
of mathematics and computer science that creates, analyzes, and implements algorithms for solving nu- merically mathematical models in science and engineering, and numerical analysis of increasing sophistication has been of numerical analysis varies from quite theoretical mathematical studies (e.g. see [5]) to computer science
ATS 680 A6: Applied Numerical Weather Prediction MW, 1:00-1:50 PM, ACRC Room 212B
, Parameterization Schemes: Keys to Understanding Numerical Weather Prediction Models, Cambridge University Press in class. Numerical model The primary numerical model that will be u
Grossman, A.; Errico, R.M.
1999-11-29T23:59:59.000Z
The vertical motion of constant density atmospheric balloons has been considered via an equation of motion for the vertical displacement of a balloon, due to vertical air motion, which can be numerically solved for balloon positions. Initial calculations are made for a constant density atmosphere. Various vertical wind models with relatively large amplitudes are applied to the model to determine how tightly the balloons are coupled to the reference level and the time scale for the balloons to change to the wind driven reference altitude. A surface launch of a balloon to a 6 km reference altitude is modeled using a detailed atmospheric pressure-density-temperature profile in the equation of motion. The results show the balloons to be relatively tightly coupled ({approx} 50-100 m) to the reference altitude.
Buscaglia, Gustavo C.
2012-01-01T23:59:59.000Z
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING Int. J. Numer. Meth. Biomed iterative coupling of dimen- sionally heterogeneous flow models in computational hemodynamics. INTRODUCTION In recent years, several efforts have been directed at integrating different mathematical models
Shekhar, Ravi
2009-05-15T23:59:59.000Z
and amplitude variation with offset (AVO) results for our example model predicts that CO2 is easier to detect than brine in the fractured reservoirs. The effects of geochemical processes on seismics are simulated by time-lapse modeling for t = 1000 years. My...
Class Generation for Numerical Wind Atlases
Class Generation for Numerical Wind Atlases RisÃ¸ National Laboratory Wind Energy Department and The Technical University of Denmark Informatics and Mathematical Modelling Department Nicholas J. Cutler s000144 Constructing a Numerical Wind Atlas 5 2.1 Introduction
Numerical modeling of a 2K J-T heat exchanger used in Fermilab Vertical Test Stand VTS-1
Gupta, Prabhat Kumar [Raja Ramanna Centre for Advanced Technology (RRCAT), Indore (MP), India; Rabehl, Roger [FNAL
2014-07-01T23:59:59.000Z
Fermilab Vertical Test Stand-1 (VTS-1) is in operation since 2007 for testing the superconducting RF cavities at 2 K. This test stand has single layer coiled finned tubes heat exchanger before J-T valve. A finite difference based thermal model has been developed in Engineering Equation Solver (EES) to study its thermal performance during filling and refilling to maintain the constant liquid level of test stand. The model is also useful to predict its performance under other various operating conditions and will be useful to design the similar kind of heat exchanger for future needs. Present paper discusses the different operational modes of this heat exchanger and its thermal characteristics under these operational modes. Results of this model have also been compared with the experimental data gathered from the VTS-1 heat exchanger and they are in good agreement with the present model.
Bader, Aaron Craig
2012-01-01T23:59:59.000Z
In this thesis we discuss measurements and modeling of minority heated fast-ion distributions in the Ion Cyclotron Range of Frequencies (ICRF) on the Alcator C-Mod tokamak. Analysis of fast-ions >100Te is important for ...
Kandelous, Maziar M.; Šim?nek, Ji?í
2010-01-01T23:59:59.000Z
and clay (33.5%), and water contents for pressures of -33times, antecedent water contents, and emitter locations, andthe saturated soil water content as one of the model
Lee, C.; Schwab, D. J.; Beletsky, D.; Stroud, J.; Lesht, B.; PNNL; NOAA; Univ. of Michigan; Univ. of Pennsylvania
2007-02-17T23:59:59.000Z
A two-dimensional sediment transport model capable of simulating sediment resuspension of mixed (cohesive plus noncohesive) sediment is developed and applied to quantitatively simulate the March 1998 resuspension events in southern Lake Michigan. Some characteristics of the model are the capability to incorporate several floc size classes, a physically based settling velocity formula, bed armoring, and sediment availability limitation. Important resuspension parameters were estimated from field and laboratory measurement data. The model reproduced the resuspension plume (observed by the SeaWIFS satellite and field instruments) and recently measured sedimentation rate distribution (using radiotracer techniques) fairly well. Model results were verified with field measurements of suspended sediment concentration and settling flux (by ADCPs and sediment traps). Both wave conditions and sediment bed properties (critical shear stress, fine sediment fraction, and limited sediment availability or source) are the critical factors that determine the concentration distribution and width of the resuspension plume. The modeled sedimentation pattern shows preferential accumulation of sediment on the eastern side of the lake, which agrees with the observed sedimentation pattern despite a predominance of particle sources from the western shoreline. The main physical mechanisms determining the sedimentation pattern are (1) the two counter-rotating circulation gyres producing offshore mass transport along the southeastern coast during northerly wind and (2) the settling velocity of sediment flocs which controls the deposition location.
Lee, Cheegwan; Schwab, David J.; Beletsky, Dmitry; Stroud, Jonathan; Lesht, B. M.
2007-02-17T23:59:59.000Z
A two-dimensional sediment transport model capable of simulating sediment resuspension of mixed (cohesive+noncohesive) sediment is developed and applied to quantitatively simulate the March 1998 resuspension events in southern Lake Michigan. Some characteristics of the model are the capability to incorporate several floc size classes, a physically-based settling velocity formula, bed armoring, and sediment availability limitation. Important resuspension parameters were estimated from field and laboratory measurement data. The model reproduced the resuspension plume (observed by the SeaWIFS satellite and field instruments) and recently measured sedimentation rate distribution (using radiotracer techniques) fairly well. Model results were verified with field measurements of suspended sediment concentration and settling flux (by ADCPs and sediment traps). Both wave conditions and sediment bed properties (critical shear stress, fine sediment fraction, and limited sediment availability or source) are the critical factors that determine the concentration distribution and width of the resuspension plume. The modeled sedimentation pattern shows preferential accumulation of sediment on the eastern side of the lake, which agrees with the observed sedimentation pattern despite a predominance of particle sources from the western shoreline. The main physical mechanisms determining the sedimentation pattern are 1) the two counter-rotating circulation gyres producing offshore mass transport along the southeastern coast during northerly wind and 2) the settling velocity of sediment flocs which controls the deposition location.
Modelling and numerical approximation of a 2.5D set of equations for mesoscale atmospheric processes
Kalise, Dante
2011-01-01T23:59:59.000Z
The set of 3D inviscid primitive equations for the atmosphere is dimensionally reduced by a Discontinuous Galerkin discretization in one horizontal direction. The resulting model is a 2D system of balance laws where with a source term depending on the layering procedure and the choice of coupling fluxes, which is established in terms of upwind considerations. The "2.5D" system is discretized via a WENO-TVD scheme based in a flux limiter centered approach. We study four tests cases related to atmospheric phenomena to analyze the physical validity of the model.
NUMERICAL ANALYSIS: This refers to the analysis of mathematical problems by numerical means, es-
Atkinson, Kendall
NUMERICAL ANALYSIS: This refers to the analysis of mathematical problems by numerical means, es- pecially mathematical problems arising from models based on calculus. Effective numerical analysis requires such as rootfinding and numerical integration; but we will also look at the structure of computers and the impli
Mueller, Karl
for the mechanics and kinematics of compressive wedges Phaedra Upton,1,2 Karl Mueller,3 and Yue-Gau Chen4 Received develop three-dimensional mechanical models of a compressive wedge and investigate how the form and kinematics of the outboard wedge are affected by variation in initial topography, material properties
Scott Geraedts; Olexei Motrunich
2014-08-05T23:59:59.000Z
We study a topological phase of interacting bosons in (3+1) dimensions which is protected by charge conservation and time-reversal symmetry. We present an explicit lattice model which realizes this phase and which can be studied in sign-free Monte Carlo simulations. The idea behind our model is to bind bosons to topological defects called hedgehogs. We determine the phase diagram of the model and identify a phase where such bound states are proliferated. In this phase we observe a Witten effect in the bulk whereby an external monopole binds half of the elementary boson charge, which confirms that it is a bosonic topological insulator. We also study the boundary between the topological insulator and a trivial insulator. We find a surface phase diagram which includes exotic superfluids, a topologically ordered phase, and a phase with a Hall effect quantized to one-half of the value possible in a purely two-dimensional system. We also present models that realize symmetry-enriched topologically-ordered phases by binding multiple hedgehogs to each boson; these phases show charge fractionalization and intrinsic topological order as well as a fractional Witten effect.
Preprint of the paper "A Numerical Refraction-Diffraction Model in the Ray-front Coordinate System"
Colominas, Ignasi
theory was adapted to water waves3 , allowing engineers and scientists to assess the diffractive effects" (CD-ROM), Parte IV: "Fluid Mechanics", Sección 6: "Shallow Water and Water Waves". S.R. Idelsohn, E, Ray Methods. Abstract. A model for the propagation of water waves is presented, which retrieves
Pruess, Karsten; Zhang, Keni
2008-11-17T23:59:59.000Z
For purposes of geologic storage, CO2 would be injected into saline formations at supercritical temperature and pressure conditions, and would form a separate phase that is immiscible with the aqueous phase (brine). At typical subsurface temperature and pressure conditions, supercritical CO2 (scCO2) has lower density than the aqueous phase and would experience an upward buoyancy force. Accordingly, the CO2 is expected to accumulate beneath the caprock at the top of the permeable interval, and could escape from the storage formation wherever (sub-)vertical pathways are available, such as fractures or faults through the caprock, or improperly abandoned wells. Over time, an increasing fraction of CO2 may dissolve in the aqueous phase, and eventually some of the aqueous CO2 may react with rock minerals to form poorly soluble carbonates. Dissolution into the aqueous phase and eventual sequestration as carbonates are highly desirable processes as they would increase permanence and security of storage. Dissolution of CO2 will establish phase equilibrium locally between the overlying CO2 plume and the aqueous phase beneath. If the aqueous phase were immobile, CO2 dissolution would be limited by the rate at which molecular diffusion can remove dissolved CO2 from the interface between CO2-rich and aqueous phases. This is a slow process. However, dissolution of CO2 is accompanied by a small increase in the density of the aqueous phase, creating a negative buoyancy force that can give rise to downward convection of CO2-rich brine, which in turn can greatly accelerate CO2 dissolution. This study explores the process of dissolution-diffusion-convection (DDC), using high-resolution numerical simulation. We find that geometric features of convection patterns are very sensitive to small changes in problem specifications, reflecting self-enhancing feedbacks and the chaotic nature of the process. Total CO2 dissolution rates on the other hand are found to be quite robust against modest changes in problem parameters, and are essentially constant as long as no dissolved CO2 reaches the lower boundary of the system.
Keisuke Nishida; Masaki Shimizu; Daikou Shiota; Hiroyuki Takasaki; Tetsuya Magara; Kazunari Shibata
2008-09-04T23:59:59.000Z
The plasmoid-induced-reconnection model explaining solar flares based on bursty reconnection produced by an ejecting plasmoid suggests a possible relation between the ejection velocity of a plasmoid and the rate of magnetic reconnection. In this study, we focus on the quantitative description of this relation. We performed magnetohydrodynamic (MHD) simulations of solar flares by changing the values of resistivity and the plasmoid velocity. The plasmoid velocity has been changed by applying an additional force to the plasmoid to see how the plasmoid velocity affects the reconnection rate. An important result is that the reconnection rate has a positive correlation with the plasmoid velocity, which is consistent with the plasmoid-induced-reconnection model for solar flares. We also discuss an observational result supporting this positive correlation.
Ali, Melkamu; Ye, Sheng; Li, Hongyi; Huang, Maoyi; Leung, Lai-Yung R.; Fiori, Aldo; Sivapalan, Murugesu
2014-07-19T23:59:59.000Z
Subsurface stormflow is an important component of the rainfall-runoff response, especially in steep forested regions. However; its contribution is poorly represented in current generation of land surface hydrological models (LSMs) and catchment-scale rainfall-runoff models. The lack of physical basis of common parameterizations precludes a priori estimation (i.e. without calibration), which is a major drawback for prediction in ungauged basins, or for use in global models. This paper is aimed at deriving physically based parameterizations of the storage-discharge relationship relating to subsurface flow. These parameterizations are derived through a two-step up-scaling procedure: firstly, through simulations with a physically based (Darcian) subsurface flow model for idealized three dimensional rectangular hillslopes, accounting for within-hillslope random heterogeneity of soil hydraulic properties, and secondly, through subsequent up-scaling to the catchment scale by accounting for between-hillslope and within-catchment heterogeneity of topographic features (e.g., slope). These theoretical simulation results produced parameterizations of the storage-discharge relationship in terms of soil hydraulic properties, topographic slope and their heterogeneities, which were consistent with results of previous studies. Yet, regionalization of the resulting storage-discharge relations across 50 actual catchments in eastern United States, and a comparison of the regionalized results with equivalent empirical results obtained on the basis of analysis of observed streamflow recession curves, revealed a systematic inconsistency. It was found that the difference between the theoretical and empirically derived results could be explained, to first order, by climate in the form of climatic aridity index. This suggests a possible codependence of climate, soils, vegetation and topographic properties, and suggests that subsurface flow parameterization needed for ungauged locations must account for both the physics of flow in heterogeneous landscapes, and the co-dependence of soil and topographic properties with climate, including possibly the mediating role of vegetation.
Victoria-Rueda, Carlos Humberto
1981-01-01T23:59:59.000Z
). The amount of organic carbon yield by the leaves of Thalassia compares well, and in many cases surpasses, average values reported for tropical rain forest and reef ecosystems and intensively cultivated agroecosystems such as hybrid corn, rice... course of my graduate studies. I extend my grat1tude to Dr, Rezneat Darnell, co-chairman of the committee, who initiated my interest in the top1c of ecosystem modeling. His constructive crit1cism and contr1bution to this work is fully appreciated. I...
Pal, Sujay
2015-01-01T23:59:59.000Z
In this thesis, we theoretically predict the normal characteristics of Very Low Frequency (3~30 kHz) radio wave propagation through Earth-ionosphere waveguide corresponding to normal behavior of the D-region ionosphere. We took the VLF narrow band data from the receivers of Indian Centre for Space Physics (ICSP) to validate our models. Detection of sudden ionospheric disturbances (SIDs) are common to all the measurements. We apply our theoretical models to infer the D-region characteristics and to reproduce the observed VLF signal behavior corresponding to such SIDs. We develop a code based on ray theory to simulate the diurnal behavior of VLF signals over short propagation paths (2000~3000 km). The diurnal variation from this code are comparable to the variation obtained from a more general Long Wave Propagation Capability (LWPC) code which is based on mode theory approach. We simulate the observational results obtained during the Total Solar Eclipse of July 22, 2009 in India. We also report and simulate a h...
Yang, L.H.; Ding, Y.M. [China University of Mining & Technology, Xuzhou (China). College of Resources and Geoscience
2009-07-01T23:59:59.000Z
The aim of this article is to discuss the distribution law of the gasification agent concentration in a deep-going way during underground coal gasification and the new method of solving the problem for the convection diffusion of the gas. In this paper, the basic features of convection diffusion for the gas produced in underground coal gasification are studied. On the basis of the model experiment, through the analysis of the distribution and patterns of variation for the fluid concentration field in the process of the combustion and gasification of the coal seams within the gasifier, the 3-D non-linear unstable mathematical models on the convection diffusion for oxygen are established. In order to curb such pseudo-physical effects as numerical oscillation and surfeit which frequently occurred in the solution of the complex mathematical models, the novel finite unit algorithm, the upstream weighted multi-cell balance method is advanced in this article, and its main derivation process is introduced.
Rockhold, Mark L.; Bacon, Diana H.; Freedman, Vicky L.; Lindberg, Michael J.; Clayton, Ray E.
2012-03-19T23:59:59.000Z
To characterize the extent of contamination under the 324 Building, a pit was excavated on the north side of the building in 2010 by Washington Closure Hanford LLC (WCH). Horizontal closed-end steel access pipes were installed under the foundation of the building from this pit and were used for measuring temperatures and exposure rates under the B-Cell. The deployed sensors measured elevated temperatures of up to 61 C (142 F) and exposure rates of up to 8,900 R/hr. WCH suspended deactivation of the facility because it recognized that building safety systems and additional characterization data might be needed for remediation of the contaminated material. The characterization work included additional field sampling, laboratory measurements, and numerical flow and transport modeling. Laboratory measurements of sediment physical, hydraulic, and geochemical properties were performed by Pacific Northwest National Laboratory (PNNL) and others. Geochemical modeling and subsurface flow and transport modeling also were performed by PNNL to evaluate the possible extent of contamination in the unsaturated sand and gravel sediments underlying the building. Historical records suggest that the concentrated 137Cs- and 90Sr-bearing liquid wastes that were spilled in B-Cell were likely from a glass-waste repository testing program associated with the Federal Republic of Germany (FRG). Incomplete estimates of the aqueous chemical composition (no anion data provided) of the FRG waste solutions were entered into a geochemical speciation model and were charge balanced with nitrate to estimate waste composition. Additional geochemical modeling was performed to evaluate reactions of the waste stream with the concrete foundation of the building prior to the stream entering the subsurface.
Movshovitz, Naor; Asphaug, Erik; Korycansky, Donald, E-mail: nmovshov@ucsc.edu [Department of Earth and Planetary Sciences, University of California, Santa Cruz, 1156 High Street, Santa Cruz, CA 95064 (United States)
2012-11-10T23:59:59.000Z
We advance the modeling of rubble-pile solid bodies by re-examining the tidal breakup of comet Shoemaker-Levy 9, an event that occurred during a 1.33 R encounter with Jupiter in 1992 July. Tidal disruption of the comet nucleus led to a chain of sub-nuclei {approx}100-1000 m diameter; these went on to collide with the planet two years later. They were intensively studied prior to and during the collisions, making SL9 the best natural benchmark for physical models of small-body disruption. For the first time in the study of this event, we use numerical codes treating rubble piles as collections of polyhedra. This introduces forces of dilatation and friction, and inelastic response. As in our previous studies we conclude that the progenitor must have been a rubble pile, and we obtain approximately the same pre-breakup diameter ({approx}1.5 km) in our best fits to the data. We find that the inclusion of realistic fragment shapes leads to grain locking and dilatancy, so that even in the absence of friction or other dissipation we find that disruption is overall more difficult than in our spheres-based simulations. We constrain the comet's bulk density at {rho}{sub bulk} {approx} 300-400 kg m{sup -3}, half that of our spheres-based predictions and consistent with recent estimates derived from spacecraft observations.
Geometric Numerical Methods for Numerical Weather Prediction
Langdon, Stephen
-Mesh (HPM) Method · Label space is discretised into N particles with coordinates on the momentum phase space and Sij = (1 - ^2xx)-1. Geometric Numerical Methods for Numerical Weather Prediction p. 8/28 #12;HPM Equations of shallow water motions · The canonical HPM equations of 1D shallow water motion on TS1 are P
Numerical Modeling | Open Energy Information
AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on DeliciousPlasmaP a gHigh Plains WindInformationNVN-079666New look toNotes from theFractures below a Borehole
Lateral boundary errors in regional numerical weather
?umer, Slobodan
Lateral boundary errors in regional numerical weather prediction models Author: Ana Car Advisor weather services for short- range forecasts. These models are covering smaller areas with higher reso Introduction Equations for numerical weather prediction (NWP) are mathematical represen- ation of physical
High performance computing and numerical Volker Springel
Masci, Frank
High performance computing and numerical modelling Volker Springel Lecture Notes 43rd Saas-mail: volker.springel@h-its.org 1 arXiv:1412.5187v1[astro-ph.GA]16Dec2014 #12;#12;Contents High performance computing and numerical modelling . . . . . . . . . . . . . . . 1 Volker Springel 1 Preamble
Souza, Paulo Victor Santos
2015-01-01T23:59:59.000Z
Flow of viscous fluids are not usually discussed in detail in general and basic courses of physics. This is due in part to the fact that the Navier-Stokes equation has analytical solution only for a few restricted cases, while more sophisticated problems can only be solved by numerical methods. In this text, we present a computer simulation of wind tunnel, i.e., we present a set of programs to solve the Navier-Stokes equation for an arbitrary object inserted in a wind tunnel. The tunnel enables us to visualize the formation of vortices behind object, the so-called von K\\'arm\\'an vortices, and calculate the drag force on the object. We believe that this numerical wind tunnel can support the teacher and allow a more elaborate discussion of viscous flow. The potential of the tunnel is exemplified by the study of the drag on a simplified model of wing whose angle of attack can be controlled. A link to download the programs that make up the tunnel appears at the end.
Eindhoven, Technische Universiteit
Numerical Simulation of PulseTube Refrigerators: 1D model I.A. Lyulina 1 , R.M.M. Mattheij 1 , A of a pulsetube refrigerator. Conservation equations describing compressible gas flow in the tube are solved, calculate the average enthalpy flow and estimate the refrigeration power. Keywords: pulsetube refrigerator
Numerical simulation of hydraulic fracturing
Warner, Joseph Barnes
1987-01-01T23:59:59.000Z
~ared that the results of such treatments were not always adequately described by the two-dimensional models. With recent advances in hydraulic fracturing and computing technology, attempts have been made to formulate more realistic fracture models. These three...NUMERICAL SIMULATION OF HYDRAULIC FRACTURING A Thesis by JOSEPH BARNES WARNER Submitted to the Graduate College of Texas A&M University in partial fulfillment of the requirements for the degree of MASTER OF SCIENCE May 1987 Maj or Subj ect...
Knight, Matthew M.; Schleicher, David G.; Schwieterman, Edward W.; Christensen, Samantha R. [Lowell Observatory, 1400 W. Mars Hill Road, Flagstaff, AZ 86001 (United States); Farnham, Tony L., E-mail: knight@lowell.edu [Department of Astronomy, University of Maryland, College Park, MD 20742-2421 (United States)
2012-11-01T23:59:59.000Z
We report on photometry and imaging of Comet 10P/Tempel 2 obtained at Lowell Observatory from 1983 through 2011. We measured a nucleus rotation period of 8.950 {+-} 0.002 hr from 16 nights of imaging acquired between 2010 September and 2011 January. This rotation period is longer than the period we previously measured in 1999, which was itself longer than the period measured in 1988, and demonstrates that Tempel 2 is continuing to spin down, presumably due to torques caused by asymmetric outgassing. A nearly linear jet was observed which varied little during a rotation cycle in both R and CN images acquired during the 1999 and 2010 apparitions. We measured the projected direction of this jet throughout the two apparitions and, under the assumption that the source region of the jet was near the comet's pole, determined a rotational pole direction of R.A./decl. = 151 Degree-Sign /+59 Degree-Sign from CN measurements and R.A./decl. = 173 Degree-Sign /+57 Degree-Sign from dust measurements (we estimate a circular uncertainty of 3 Degree-Sign for CN and 4 Degree-Sign for dust). Different combinations of effects likely bias both gas and dust solutions and we elected to average these solutions for a final pole direction of R.A./decl. = 162 Degree-Sign {+-} 11 Degree-Sign /+58 Degree-Sign {+-} 1 Degree-Sign . Photoelectric photometry was acquired on 3 nights in 1983, 2 nights in 1988, 19 nights in 1999/2000, and 10 nights in 2010/2011. The activity exhibited a steep 'turn-on' {approx}3 months prior to perihelion (the exact timing of which varies) and a relatively smooth decline after perihelion. The activity during the 1999 and 2010 apparitions was similar; limited data in 1983 and 1988 (along with IUE data from the literature) were systematically higher and the difference cannot be explained entirely by the smaller perihelion distance. We measured a 'typical' composition, in agreement with previous investigators. Monte Carlo numerical modeling with our pole solution best replicated the observed coma morphology for a source region located near a comet latitude of +80 Degree-Sign and having a radius of {approx}10 Degree-Sign . Our model reproduced the seasonal changes in activity, suggesting that the majority of Tempel 2's activity originates from a small active region located near the pole. We also find that a cosine-squared solar angle function gives the best fit as compared to a standard cosine function.
2005-10-20T23:59:59.000Z
Page 1. LAB #8. Numerical Methods. Goal: The purpose of this lab is to explain how computers numerically ... Also you will examine what .... (7) Now consider the differential equation ... 3-exp(2*y)+sqrt(t)/y; (Don't forget the “;” at the end.).
Numerical simulation of micro-fluidic passive and active mixers
Kumar, Saurabh
2002-01-01T23:59:59.000Z
Numerical simulations of mixing using passive and active techniques are performed. For passive mixing, numerical modeling of a micro-fluidic device, build by Holden and Cremer, was performed. The micro-fluidic device consists of a Y...
Numerical Simulation of Laminar Reacting Flows with Complex Chemistry
Bell, John B.
Numerical Simulation of Laminar Reacting Flows with Complex Chemistry M S Day and J B Bell Lawrence: 47.40.Fw, 82.40.Py Submitted to: Combust. Theory Modelling #12;Numerical Simulation of Laminar
Meng, Da; Zheng, Bin; Lin, Guang; Sushko, Maria L.
2014-08-29T23:59:59.000Z
We have developed efficient numerical algorithms for the solution of 3D steady-state Poisson-Nernst-Planck equations (PNP) with excess chemical potentials described by the classical density functional theory (cDFT). The coupled PNP equations are discretized by finite difference scheme and solved iteratively by Gummel method with relaxation. The Nernst-Planck equations are transformed into Laplace equations through the Slotboom transformation. Algebraic multigrid method is then applied to efficiently solve the Poisson equation and the transformed Nernst-Planck equations. A novel strategy for calculating excess chemical potentials through fast Fourier transforms is proposed which reduces computational complexity from O(N2) to O(NlogN) where N is the number of grid points. Integrals involving Dirac delta function are evaluated directly by coordinate transformation which yields more accurate result compared to applying numerical quadrature to an approximated delta function. Numerical results for ion and electron transport in solid electrolyte for Li ion batteries are shown to be in good agreement with the experimental data and the results from previous studies.
Multigrid methods for improving the variational data assimilation in numerical weather prediction
Kwak, Do Young
Multigrid methods for improving the variational data assimilation in numerical weather prediction conditions are needed to solve numerical weather prediction models: initial condition and boundary condition: numerical weather prediction, variational data assimilation, minimization procedure, multigrid methods, cell
Disruptive Innovation in Numerical Hydrodynamics
Waltz, Jacob I. [Los Alamos National Laboratory
2012-09-06T23:59:59.000Z
We propose the research and development of a high-fidelity hydrodynamic algorithm for tetrahedral meshes that will lead to a disruptive innovation in the numerical modeling of Laboratory problems. Our proposed innovation has the potential to reduce turnaround time by orders of magnitude relative to Advanced Simulation and Computing (ASC) codes; reduce simulation setup costs by millions of dollars per year; and effectively leverage Graphics Processing Unit (GPU) and future Exascale computing hardware. If successful, this work will lead to a dramatic leap forward in the Laboratory's quest for a predictive simulation capability.
Introduction to exact numerical computation
Escardó, Martín
Introduction to exact numerical computation Notes 2 Floating-point computation 4 3 Exact numerical (called the operational semantics) for computing a syntactical representative of the mathematical entity
Threedimensional numerical simulation for various geometries
Herbin, Raphaèle
modelling and numerical simulation of natural gasfed solid oxide cells (Solid Oxide Fuel Cell, SOFC) at a stationary regime. The principle of a Solid Oxide Fuel Cell (SOFC) is based on the conversion of the chemical is taken into account in the present model. The SOFC systems seem to be of great interest for use
Xu, Tianfu; Apps, John A.; Pruess, Karsten; Yamamoto, Hajime
2008-01-01T23:59:59.000Z
Australia. Am. Assoc. Petrol. Geol. Bull. 76(8), 1224-1249.modeling. Mineral. and Petrol. 59, 121-140. Gunter, W.D. ,
Witherspoon, P.A.
2010-01-01T23:59:59.000Z
n The variable is the hydraulic head Energy T r a n s p o rhydraulic fractur ing process requires the modeling of the energy
Gupta, A.; Moridis, G.J.; Kneafsey, T.J.; Sloan, Jr., E.D.
2009-08-15T23:59:59.000Z
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.
Pohlmann Karl,Ye Ming
2012-03-01T23:59:59.000Z
Models of groundwater flow for the Yucca Flat area of the Nevada National Security Site (NNSS) are under development by the U.S. Department of Energy (DOE) for corrective action investigations of the Yucca Flat-Climax Mine Corrective Action Unit (CAU). One important aspect of these models is the quantity of inter-basin groundwater flow from regional systems to the north. This component of flow, together with its uncertainty, must be properly accounted for in the CAU flow models to provide a defensible regional framework for calculations of radionuclide transport that will support determinations of the Yucca Flat-Climax Mine contaminant boundary. Because characterizing flow boundary conditions in northern Yucca Flat requires evaluation to a higher level of detail than the scale of the Yucca Flat-Climax Mine CAU model can efficiently provide, a study more focused on this aspect of the model was required.
McConnell, Joshua B
2000-01-01T23:59:59.000Z
failure in the borehole wall. Derivation of energy, mass and momentum conservation equations was performed to aid in the identification of applicable loads acting on the melted material that will form the liner. A finite difference model was coded...
Feng, Yang
2012-07-16T23:59:59.000Z
The hypoxic area in the Gulf of Mexico is the second largest in the world, which has received extensive scientific study and management interest. Previous modeling studies have concluded that the increased hypoxic area in ...
Numerical Integration Gordon K. Smyth
Smyth, Gordon K.
Numerical Integration Gordon K. Smyth in Encyclopedia of Biostatistics (ISBN 0471 975761) Edited by Peter Armitage and Theodore Colton John Wiley & Sons, Ltd, Chichester, 1998 #12;Numerical Integration Numerical integration is the study of how the numerical value of an integral can be found. Also called
Prowell, I.; Robertson, A.; Jonkman, J.; Stewart, G. M.; Goupee, A. J.
2013-01-01T23:59:59.000Z
Realizing the critical importance the role physical experimental tests play in understanding the dynamics of floating offshore wind turbines, the DeepCwind consortium conducted a one-fiftieth-scale model test program where several floating wind platforms were subjected to a variety of wind and wave loading condition at the Maritime Research Institute Netherlands wave basin. This paper describes the observed behavior of a tension-leg platform, one of three platforms tested, and the systematic effort to predict the measured response with the FAST simulation tool using a model primarily based on consensus geometric and mass properties of the test specimen.
G. Litak; M. Wendeker; M. Krupa; J. Czarnigowski
2004-05-23T23:59:59.000Z
We examine a simple, fuel-air, model of combustion in a spark ignition (si) engine with indirect injection. In our two fluid model, variations of fuel mass burned in cycle sequences appear due to stochastic fluctuations of a fuel feed amount. We have shown that a small amplitude of these fluctuations affects considerably the stability of a combustion process strongly depending on the quality of air-fuel mixture. The largest influence was found in the limit of a lean combustion. The possible effect of nonlinearities in the combustion process has been also discussed.
Bonèa, Janez
, such as the self-consistent Born approximation SCBA ,26 self- consistent perturbational approach SCPA ,5 the Green's function, defined in imaginary time. The SCBA and SCPA calculations are likewise limited results8,9,18 are not necessarily compa- rable with SCBA, SCPA, and ED results of the t-J model
Mohaghegh, Shahab
of SPE copyright. Abstract Latest advances in shale gas reservoir simulation and modeling have made it possible to optimize and enhance the production from organic rich shale gas reservoirs. Reservoir simulator is no longer used with a simple description of the complex shale gas reservoirs, but with multiple, equally
Eric M. Leonard; Mitchell A. Plummer; Paul E. Carrara
2014-04-01T23:59:59.000Z
Well-preserved moraines from the penultimate, or Bull Lake, glaciation of Snowmass Creek Valley in the Elk Range of Colorado present an opportunity to examine the character of the high-altitude climate in the Rocky Mountains during Marine Oxygen Isotope Stage 6. This study employs a 2-D coupled mass/energy balance and flow model to assess the magnitudes of temperature and precipitation change that could have sustained the glacier in mass-balance equilibrium at its maximum extent during the Bull Lake glaciation. Variable substrate effects on glacier flow and ice thickness make the modeling somewhat more complex than in geologically simpler settings. Model results indicate that a temperature depression of about 6.7°C compared with the present (1971–2000 AD) would have been necessary to sustain the Snowmass Creek glacier in mass-balance equilibrium during the Bull Lake glaciation, assuming no change in precipitation amount or seasonality. A 50% increase or decrease from modern precipitation would have been coupled with 5.2°C and 9.1°C Bull Lake temperature depressions respectively. Uncertainty in these modeled temperature depressions is about 1°C.
Hanna, S.R.; Chang, J.C.
1997-01-01T23:59:59.000Z
The HGSYSTEM/UF{sub 6} model was developed for use in preparing Safety Analysis Reports (SARs) by estimating the consequences of possible accidental releases of UF{sub 6} to the atmosphere at the gaseous diffusion plants (GDPs) located in Portsmouth, Ohio, and Paducah, Kentucky. Although the latter report carries a 1996 date, the work that is described was completed in late 1994. When that report was written, the primary release scenarios of interest were thought to be gas pipeline and liquid tank ruptures over open terrain away from the influence of buildings. However, upon further analysis of possible release scenarios, the developers of the SARs decided it was necessary to also consider accidental releases within buildings. Consequently, during the fall and winter of 1995-96, modules were added to HGSYSTEM/UF{sub 6} to account for flow and dispersion around buildings. The original HGSYSTEM/UF{sub 6} model also contained a preliminary method for accounting for the possible lift-off of ground-based buoyant plumes. An improved model and a new set of wind tunnel data for buoyant plumes trapped in building recirculation cavities have become available that appear to be useful for revising the lift-off algorithm and modifying it for use in recirculation cavities. This improved lift-off model has been incorporated in the updated modules for dispersion around buildings.
Dimov, Ivan
investigations related to air pollution in Bulgaria ([6, 7]), Denmark ([8-11]), England ([12]), Europe ([2, 4 been used in some inter-comparisons of European large-scale air pollution models ([20, 21]). Recently Publishers CHAPTER 4 Part C: Sensitivity of European Pollution Levels to Changes of Human-Made Emissions
NUMERICAL SIMULATIONS OF TRANSVERSE COMPRESSION AND DENSIFICATION IN WOOD
Nairn, John A.
NUMERICAL SIMULATIONS OF TRANSVERSE COMPRESSION AND DENSIFICATION IN WOOD John A. Nairn1 Professor- terials is a useful tool for stress analysis and for failure modeling. Although FEA of wood as an anisotropic continuum is used, numerical modeling of realistic wood structures, including details of wood
Numerical simulation of tsunami waves generated by deformable submarine landslides
Kirby, James T.
Numerical simulation of tsunami waves generated by deformable submarine landslides Gangfeng Ma a wave model Tsunami wave Numerical modeling a b s t r a c t This paper presents a new submarine of landslide motion and associated tsunami wave generation on parameters including sediment settling velocity
Numerical Weather Forecasting at the Savannah River Site
Buckley, R.L. [Westinghouse Savannah River Company, AIKEN, SC (United States)
1998-08-01T23:59:59.000Z
This paper discusses the use of an advanced three-dimensional prognostic numerical model to provide space and time-dependent meteorological data for use in the WIND System dispersion models.
Numerical simulations of mine penetration in soft sediments
Yao, Zhigang
2003-01-01T23:59:59.000Z
This study is concerned with developing a predictive model for depth of penetration of mines in soft sediments. The primary objective of this numerical study is to model soil resistance as a function of penetration depth using the finite element...
Numerical simulation of oxidation and its effect on the crack growth resistance of titanium alloys
Triharjanto, Robertus Heru
1997-01-01T23:59:59.000Z
of titanium alloys used as metal matrix in MMC is imperative. This research effort develops a numerical modeling procedure to simulate the oxidation propagating from the crack surfaces and evaluate the effect on the mechanical behavior. The numerical modeling...
Numerical Integration Gordon K. Smyth
Smyth, Gordon K.
Numerical Integration Gordon K. Smyth May 1997 Numerical integration is the study of how the numerical value of an integral can be found. Also called quadrature, which refers to finding a square whose \\Lambda . Of central interest is the process of approximating a definite integral from values of the in
Numerical simulation of an electroweak oscillon
Graham, N. [Department of Physics, Middlebury College, Middlebury, Vermont 05753 (United States)
2007-10-15T23:59:59.000Z
Numerical simulations of the bosonic sector of the SU(2)xU(1) electroweak standard model in 3+1 dimensions have demonstrated the existence of an oscillon--an extremely long-lived, localized, oscillatory solution to the equations of motion--when the Higgs mass is equal to twice the W{sup {+-}} boson mass. It contains total energy roughly 30 TeV localized in a region of radius 0.05 fm. A detailed description of these numerical results is presented.
Status of NINJA: the Numerical INJection Analysis project
Benjamin Aylott; John G. Baker; William D. Boggs; Michael Boyle; Patrick R. Brady; Duncan A. Brown; Bernd Brügmann; Luisa T. Buchman; Alessandra Buonanno; Laura Cadonati; Jordan Camp; Manuela Campanelli; Joan Centrella; Shourov Chatterjis; Nelson Christensen; Tony Chu; Peter Diener; Nils Dorband; Zachariah B. Etienne; Joshua Faber; Stephen Fairhurst; Benjamin Farr; Sebastian Fischetti; Gianluca Guidi; Lisa M. Goggin; Mark Hannam; Frank Herrmann; Ian Hinder; Sascha Husa; Vicky Kalogera; Drew Keppel; Lawrence E. Kidder; Bernard J. Kelly; Badri Krishnan; Pablo Laguna; Carlos O. Lousto; Ilya Mandel; Pedro Marronetti; Richard Matzner; Sean T. McWilliams; Keith D. Matthews; R. Adam Mercer; Satyanarayan R. P. Mohapatra; Abdul H. Mroué; Hiroyuki Nakano; Evan Ochsner; Yi Pan; Larne Pekowsky; Harald P. Pfeiffer; Denis Pollney; Frans Pretorius; Vivien Raymond; Christian Reisswig; Luciano Rezzolla; Oliver Rinne; Craig Robinson; Christian Röver; Lucía Santamaría; Bangalore Sathyaprakash; Mark A. Scheel; Erik Schnetter; Jennifer Seiler; Stuart L. Shapiro; Deirdre Shoemaker; Ulrich Sperhake; Alexander Stroeer; Riccardo Sturani; Wolfgang Tichy; Yuk Tung Liu; Marc van der Sluys; James R. van Meter; Ruslan Vaulin; Alberto Vecchio; John Veitch; Andrea Viceré; John T. Whelan; Yosef Zlochower
2009-05-26T23:59:59.000Z
The 2008 NRDA conference introduced the Numerical INJection Analysis project (NINJA), a new collaborative effort between the numerical relativity community and the data analysis community. NINJA focuses on modeling and searching for gravitational wave signatures from the coalescence of binary system of compact objects. We review the scope of this collaboration and the components of the first NINJA project, where numerical relativity groups shared waveforms and data analysis teams applied various techniques to detect them when embedded in colored Gaussian noise.
CSU ATS703 Fall 2012 Numerical Weather Prediction
CSU ATS703 Fall 2012 Numerical Weather Prediction ATS703 is based on the course notes and papers method. A crucial element of accurate weather prediction is initialization, which is briefly discussed in Chapter 11. In the next decade, numerical weather prediction will expe- rience a revolution in model
Applications of Numerical Optimal Control to Nonlinear Hybrid Systems
Zefran, Milo?
Applications of Numerical Optimal Control to Nonlinear Hybrid Systems Shangming Wei a , Kasemsak the system has no mem- ory, the embedding technique can be used to reduce the hybrid optimal control problem of the controllers. Key words: Hybrid systems, optimal control, numerical optimization, model predictive control
MULTISCALE NUMERICAL STUDY OF TURBULENT FLOW AND BUBBLE ENTRAINMENT
Kirby, James T.
MULTISCALE NUMERICAL STUDY OF TURBULENT FLOW AND BUBBLE ENTRAINMENT IN THE SURF ZONE BY GANGFENG MA . . . . . . . . . . . . . . . . . . . . . 3 1.1.2 Numerical Investigations . . . . . . . . . . . . . . . . . . . . . 5 1.2 Bubble Entrainment . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15 2.4 Bubble Entrainment Model . . . . . . . . . . . . . . . . . . . . . . . 18 2.5 Bubble Breakup
Numerical Simulation of a Natural Circulation Steam Generator
Weinmüller, Ewa B.
Numerical Simulation of a Natural Circulation Steam Generator W. Linzer \\Lambda , K. Ponweiser circulation steam generator. We focus on a model with a simple geometry consisting of two vertical pipes properties of water and steam. We present a numerical algorithm based on an explicit upwind discretization
Numerical simulation of sandstone reservoir models
Gross, Stephen Joseph
1983-01-01T23:59:59.000Z
0 0. 1 0. 2 Or3 0. 4 O. S 0. 6 Or7 0. 8 Or9 I GAS SATURATIOMr FRACTION OF PORE VOLUME Figure 4 ? Relative permeability of gas-oil system for Layer 1 Table 1 - Layer Properties 'v ~La er ~md ~md 1 500 250 0. 22 0. 20 0. 20 0. 05 2 200 100 0. 21... ? - 2o. 'a STB/o """ 50. 0 STB/0 r J't. , ~~'$ ' v 'i ti, 0 2 4 6 8 10 12 14 16 18 OIL RfCIMRYi 2BIP Figure 12 - Gas-oil ratio behavior of so'lution gas drive runs for Case 4, a homogeneous interval 3500 A G E 3000 2saa 2000 lsaa ? CASE 1...
Numerical and experimental modelling of microwave applicators
Dibben, David
1995-10-24T23:59:59.000Z
area of use is for the processing of food, the tempering of fruit, meat, fish and dairy produce, for the heating and cooking of meals and for sterilisation and pasteurisation. The domestic microwave oven is now a common household appliance throughout...
numerical modeling | OpenEI Community
AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on DeliciousPlasmaP a g eWorksVillage of Brewster, OhioLonghui Hydropowercitiesimages Homeload datamulticolornet
Category:Numerical Modeling | Open Energy Information
AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on DeliciousPlasmaP a gHigh4-FD-aBeijingCalifornia/IncentivesCatahoula
OpenEI Community - numerical modeling
AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on DeliciousPlasmaP a gHigh Plains WindInformationNVN-079666NewOlsen Wind FarmPublic Artnotes-call
Modeling, Analysis and Numerical Approximations of Image ...
... fields, which have a broad range of applications in science, engineering and industry. Not long ago, statistical and ad hoc methods had been main tools for studying and analyzing image processing and computer vision problems. ... interesting and challenging mathematical problems in the areas of Calculus of Variation, ...
Numerical study of high heat ux pool boiling heat transfer Ying He a,*, Masahiro Shoji b
Maruyama, Shigeo
Numerical study of high heat ¯ux pool boiling heat transfer Ying He a,*, Masahiro Shoji b , Shigeo simulation model of boiling heat transfer is proposed based on a numerical macrolayer model [S. Maruyama, M. Shoji, S. Shimizu, A numerical simulation of transition boiling heat transfer, in: Proceedings
Numerical investigation of the heating process inside an industrial furnace
Wolper, Pierre
Numerical investigation of the heating process inside an industrial furnace Proposition: Combined furnace taking into account convective, conductive and radiative heat transfer. The model: Catalysis, Energy Materials, Performance Materials and Recycling. Each business area is divided into market
An investigation of the numerical treatment of condensation
Sasson, Joseph
The simulation of complete condensation continues to challenge the numerical methods currently used for multi-phase flow modeling; especially at low pressures, the change of phase process from a two-phase mixture to liquid ...
Numerical Simulation of Laminar Reacting Flows with Complex Chemistry
Bell, John B.
Numerical Simulation of Laminar Reacting Flows with Complex Chemistry M S Day and J B Bell Lawrence Simulation of Laminar Reacting Flows 2 1. Introduction Detailed modelling of time-dependent reacting ows
Mathiesen, Patrick; Kleissl, Jan
2011-01-01T23:59:59.000Z
transport and numerical weather modeling. J. Applied cross correlations. Weather and Forecasting, 8:4, 401?of radiation for numerical weather prediction and climate
Simple intrinsic defects in GaAs : numerical supplement.
Schultz, Peter Andrew
2012-04-01T23:59:59.000Z
This Report presents numerical tables summarizing properties of intrinsic defects in gallium arsenide, GaAs, as computed by density functional theory. This Report serves as a numerical supplement to the results published in: P.A. Schultz and O.A. von Lilienfeld, 'Simple intrinsic defects in GaAs', Modelling Simul. Mater. Sci Eng., Vol. 17, 084007 (2009), and intended for use as reference tables for a defect physics package in device models. The numerical results for density functional theory calculations of properties of simple intrinsic defects in gallium arsenide are presented.
On Numerical Methods for Hyperbolic Conservation Laws and Related Equations
Bürger, Raimund
A classical kinematical model of sedimentation of small equal-sized particles dispersed in a viscous fluidOn Numerical Methods for Hyperbolic Conservation Laws and Related Equations Modelling Sedimentation with nonlocal flux, systems of nonlinear conservation modelling the sedimentation of polydisperse suspensions
A numerical study of pseudoscalar inflation
Cheng, Shu-Lin; Ng, Kin-Wang
2015-01-01T23:59:59.000Z
A numerical study of a pseudoscalar inflation having an axion-photon-like coupling is performed by solving numerically the coupled differential equations of motion for inflaton and photon mode functions from the onset of inflation to the end of reheating. The backreaction due to particle production is also included self-consistently. We find that this particular inflation model realizes the idea of a warm inflation in which a steady thermal bath is established by the particle production. In most cases this thermal bath exceeds the amount of radiation released in the reheating process. In the strong coupling regime, the transition from the inflationary to the radiation-dominated phase does not involve either a preheating or reheating process. In addition, energy density peaks produced near the end of inflation may lead to the formation of primordial black holes.
A numerical study of pseudoscalar inflation
Shu-Lin Cheng; Wolung Lee; Kin-Wang Ng
2015-08-02T23:59:59.000Z
A numerical study of a pseudoscalar inflation having an axion-photon-like coupling is performed by solving numerically the coupled differential equations of motion for inflaton and photon mode functions from the onset of inflation to the end of reheating. The backreaction due to particle production is also included self-consistently. We find that this particular inflation model realizes the idea of a warm inflation in which a steady thermal bath is established by the particle production. In most cases this thermal bath exceeds the amount of radiation released in the reheating process. In the strong coupling regime, the transition from the inflationary to the radiation-dominated phase does not involve either a preheating or reheating process. In addition, energy density peaks produced near the end of inflation may lead to the formation of primordial black holes.
An Interactive Course in Numerical Methods for the Earth Susan E. Allen
Stockie, John
are typically in physics, mathematics, chemistry, biology or geography, and exposure to numerical computing science are becoming more and more dependent on numerical simulation, modelling, and computationally courses in numerical methods using MatlabTM as the programming language, or a graduate-level mathematics
Merging of coronal and heliospheric numerical two-dimensional D. Odstrcil,1,5
California at Berkeley, University of
spatial and temporal scales between the Sun and Earth. Specialized physically based numerical models have been developed to address particular aspects of the entire system. However, an integrated modeling physically based numerical models have been developed to address particular aspects of the entire system
RELAP-7 Numerical Stabilization: Entropy Viscosity Method
R. A. Berry; M. O. Delchini; J. Ragusa
2014-06-01T23:59:59.000Z
The RELAP-7 code is the next generation nuclear reactor system safety analysis code being developed at the Idaho National Laboratory (INL). The code is based on the INL's modern scientific software development framework, MOOSE (Multi-Physics Object Oriented Simulation Environment). The overall design goal of RELAP-7 is to take advantage of the previous thirty years of advancements in computer architecture, software design, numerical integration methods, and physical models. The end result will be a reactor systems analysis capability that retains and improves upon RELAP5's capability and extends the analysis capability for all reactor system simulation scenarios. RELAP-7 utilizes a single phase and a novel seven-equation two-phase flow models as described in the RELAP-7 Theory Manual (INL/EXT-14-31366). The basic equation systems are hyperbolic, which generally require some type of stabilization (or artificial viscosity) to capture nonlinear discontinuities and to suppress advection-caused oscillations. This report documents one of the available options for this stabilization in RELAP-7 -- a new and novel approach known as the entropy viscosity method. Because the code is an ongoing development effort in which the physical sub models, numerics, and coding are evolving, so too must the specific details of the entropy viscosity stabilization method. Here the fundamentals of the method in their current state are presented.
Benjamin Aylott; John G. Baker; William D. Boggs; Michael Boyle; Patrick R. Brady; Duncan A. Brown; Bernd Brügmann; Luisa T. Buchman; Alessandra Buonanno; Laura Cadonati; Jordan Camp; Manuela Campanelli; Joan Centrella; Shourov Chatterji; Nelson Christensen; Tony Chu; Peter Diener; Nils Dorband; Zachariah B. Etienne; Joshua Faber; Stephen Fairhurst; Benjamin Farr; Sebastian Fischetti; Gianluca Guidi; Lisa M. Goggin; Mark Hannam; Frank Herrmann; Ian Hinder; Sascha Husa; Vicky Kalogera; Drew Keppel; Lawrence E. Kidder; Bernard J. Kelly; Badri Krishnan; Pablo Laguna; Carlos O. Lousto; Ilya Mandel; Pedro Marronetti; Richard Matzner; Sean T. McWilliams; Keith D. Matthews; R. Adam Mercer; Satyanarayan R. P. Mohapatra; Abdul H. Mroué; Hiroyuki Nakano; Evan Ochsner; Yi Pan; Larne Pekowsky; Harald P. Pfeiffer; Denis Pollney; Frans Pretorius; Vivien Raymond; Christian Reisswig; Luciano Rezzolla; Oliver Rinne; Craig Robinson; Christian Röver; Lucía Santamaría; Bangalore Sathyaprakash; Mark A. Scheel; Erik Schnetter; Jennifer Seiler; Stuart L. Shapiro; Deirdre Shoemaker; Ulrich Sperhake; Alexander Stroeer; Riccardo Sturani; Wolfgang Tichy; Yuk Tung Liu; Marc van der Sluys; James R. van Meter; Ruslan Vaulin; Alberto Vecchio; John Veitch; Andrea Viceré; John T. Whelan; Yosef Zlochower
2009-07-09T23:59:59.000Z
The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational-wave data analysis communities. The purpose of NINJA is to study the sensitivity of existing gravitational-wave search algorithms using numerically generated waveforms and to foster closer collaboration between the numerical relativity and data analysis communities. We describe the results of the first NINJA analysis which focused on gravitational waveforms from binary black hole coalescence. Ten numerical relativity groups contributed numerical data which were used to generate a set of gravitational-wave signals. These signals were injected into a simulated data set, designed to mimic the response of the Initial LIGO and Virgo gravitational-wave detectors. Nine groups analysed this data using search and parameter-estimation pipelines. Matched filter algorithms, un-modelled-burst searches and Bayesian parameter-estimation and model-selection algorithms were applied to the data. We report the efficiency of these search methods in detecting the numerical waveforms and measuring their parameters. We describe preliminary comparisons between the different search methods and suggest improvements for future NINJA analyses.
Numerical Approximation of Vortex Density Evolution in a Superconductor.
Styles, Vanessa
Numerical Approximation of Vortex Density Evolution in a Superconductor. C.M. Elliott & V. Styles Abstract A #12;nite volume/element approximation of a mean #12;eld model of superconducting vortices in one approximations of a two-dimensional version of the mean #12;eld model of superconducting vortices considered
Simulation of the secondary settling process with reliable numerical methods
Bürger, Raimund
are placed in alphabetical order Abstract A consistent model for the settling-compression-dispersion process The sedimentation process in the secondary settling tank (SST) is still a challenge in modelling the full Sedimentation, clarifier, conservation law, partial differential equation, numerical flux INTRODUCTION
IVCNZ 2002 SUBMISSION 1 Pose Estimation by Applied Numerical Techniques
McCane, Brendan
process of how to deform the model in order to bring it into agreement with the input image. An energy. Keywords--- Numerical Optimisation, Pose Estimation, AnalysisÂbyÂSynthesis, Optimisation Algorithms (HCI). Traditionally, researchers have split approaches into an appearanceÂbased approach and a model
Numerical simulation of turbulent jet primary breakup in Diesel engines
Helluy, Philippe
Numerical simulation of turbulent jet primary breakup in Diesel engines Peng Zeng1 Marcus Herrmann" IRMA Strasbourg, 23.Jan.2008 #12;Introduction DNS of Primary Breakup in Diesel Injection Phase Transition Modeling Turbulence Modeling Summary Outline 1 Introduction 2 DNS of Primary Breakup in Diesel
Numerical Calculation of Granular Entropy
Asenjo, Daniel; Paillusson, Fabien; Frenkel, Daan
2014-03-05T23:59:59.000Z
We present numerical simulations that allow us to compute the number of ways in which N particles can pack into a given volume V. Our technique modifies the method of Xu, Frenkel, and Liu [Phys. Rev. Lett. 106, 245502 (2011)] and outperforms...
Planning numerical approximations Richard Power
Williams, Sandra
Planning numerical approximations Richard Power Sandra Williams 21st September 2009 #12;Table proportions (e.g., more than a quarter, 25.9 per cent) Proportions are a convenient well-defined subproblem Common in factual discourse (e.g., newspaper articles) Important for generating from data (but neglected
Numerical Simulation of Transpiration Cooling
to facilitate such numerical simulations for a carbon/carbon material mounted in the side wall of a hot gas channel that are able to capture a spatially varying interplay between the hot gas flow and the coolant itself. This calls for new cooling technologies, for example, in the combustion chamber of the rocket
Numerical uncertainty in computational engineering and physics
Hemez, Francois M [Los Alamos National Laboratory
2009-01-01T23:59:59.000Z
Obtaining a solution that approximates ordinary or partial differential equations on a computational mesh or grid does not necessarily mean that the solution is accurate or even 'correct'. Unfortunately assessing the quality of discrete solutions by questioning the role played by spatial and temporal discretizations generally comes as a distant third to test-analysis comparison and model calibration. This publication is contributed to raise awareness of the fact that discrete solutions introduce numerical uncertainty. This uncertainty may, in some cases, overwhelm in complexity and magnitude other sources of uncertainty that include experimental variability, parametric uncertainty and modeling assumptions. The concepts of consistency, convergence and truncation error are overviewed to explain the articulation between the exact solution of continuous equations, the solution of modified equations and discrete solutions computed by a code. The current state-of-the-practice of code and solution verification activities is discussed. An example in the discipline of hydro-dynamics illustrates the significant effect that meshing can have on the quality of code predictions. A simple method is proposed to derive bounds of solution uncertainty in cases where the exact solution of the continuous equations, or its modified equations, is unknown. It is argued that numerical uncertainty originating from mesh discretization should always be quantified and accounted for in the overall uncertainty 'budget' that supports decision-making for applications in computational physics and engineering.
Numerical semigroups Easy bounds on ng
Elizalde, Sergi
Numerical semigroups Easy bounds on ng Improved bounds on ng The number of numerical semigroups #12;Numerical semigroups Easy bounds on ng Improved bounds on ng Motivation Definitions The tree on ng Improved bounds on ng Motivation Definitions The tree T of numerical semigroups Succession rules
Graduate Survey of Numerical Methods Background material
Corless, Robert M.
Numerical Methods, Numerical Analysis, Scientific Computing, Com- putational Mathematics, Computational. Numerical Meth- ods are what one uses to solve a problem from continuous mathematics (vide Nick Trefethen), and in what sense. Scientific Computing is the use of numerical methods to solve problems of scientific
Numerical bifurcation analysis of piecewise smooth systems
( ) ( ) ( ) ( ), , outoutinin inin xffxff xhhxgg xxxx == == #12;Numerical bifurcation analysis of piecewise smooth systems INRIA
Dynamics of Dry Friction: A Numerical Investigation
Y. F. Lim; Kan Chen
1998-03-18T23:59:59.000Z
We perform extended numerical simulation of the dynamics of dry friction, based on a model derived from the phenomenological description proposed by T. Baumberger et al.. In the case of small deviation from the steady sliding motion, the model is shown to be equivalent to the state- and rate-dependent friction law which was first introduced by Rice and Ruina on the basis of experiments on rocks. We obtain the dynamical phase diagram that agrees well with the experimental results on the paper-on-paper systems. In particular, the bifurcation between stick-slip and steady sliding are shown to change from a direct (supercritical) Hopf type to an inverted (subcritical) one as the driving velocity increases, in agreement with the experiments.
Numerical Simulations of Gamma-Ray Burst Explosions
Lazzati, Davide; López-Cámara, Diego
2015-01-01T23:59:59.000Z
Gamma-ray bursts are a complex, non-linear system that evolves very rapidly through stages of vastly different conditions. They evolve from scales of few hundred kilometers where they are very dense and hot to cold and tenuous on scales of parsecs. As such, our understanding of such a phenomenon can truly increase by combining theoretical and numerical studies adopting different numerical techniques to face different problems and deal with diverse conditions. In this review, we will describe the tremendous advancement in our comprehension of the bursts phenomenology through numerical modeling. Though we will discuss studies mainly based on jet dynamics across the progenitor star and the interstellar medium, we will also touch upon other problems such as the jet launching, its acceleration, and the radiation mechanisms. Finally, we will describe how combining numerical results with observations from Swift and other instruments resulted in true understanding of the bursts phenomenon and the challenges still lyi...
Numerical Simulation of Groundwater Withdrawal at the Nevada Test Site
Carroll, Rosemary; Giroux, Brian; Pohll, Greg; Hershey, Ronald; Russell, Charles; Howcroft, William
2004-01-28T23:59:59.000Z
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 Methods for a Nonlinear BVP Arising in Physical Oceanography
Riccardo Fazio; Alessandra Jannelli
2013-10-08T23:59:59.000Z
In this paper we report and compare the numerical results for an ocean circulation model obtained by the classical truncated boundary formulation, the free boundary approach and a quasi-uniform grid treatment of the problem. We apply a shooting method to the truncated boundary formulation and finite difference methods to both the free boundary approach and the quasi-uniform grid treatment. Using the shooting method, supplemented by the Newton's iterations, we show that the ocean circulation model cannot be considered as a simple test case. In fact, for this method we are forced to use as initial iterate a value close to the correct missing initial condition in order to be able to get a convergent numerical solution. The reported numerical results allow us to point out how the finite difference method with a quasi-uniform grid is the less demanding approach and that the free boundary approach provides a more reliable formulation than the classical truncated boundary formulation.
NUMERICAL SIMULATION OF NATURAL GAS-SWIRL BURNER
Ala Qubbaj
2005-03-01T23:59:59.000Z
A numerical simulation of a turbulent natural gas jet diffusion flame at a Reynolds number of 9000 in a swirling air stream is presented. The numerical computations were carried out using the commercially available software package CFDRC. The instantaneous chemistry model was used as the reaction model. The thermal, composition, flow (velocity), as well as stream function fields for both the baseline and air-swirling flames were numerically simulated in the near-burner region, where most of the mixing and reactions occur. The results were useful to interpret the effects of swirl in enhancing the mixing rates in the combustion zone as well as in stabilizing the flame. The results showed the generation of two recirculating regimes induced by the swirling air stream, which account for such effects. The present investigation will be used as a benchmark study of swirl flow combustion analysis as a step in developing an enhanced swirl-cascade burner technology.
Numerical methods in stochastic processes
Corbett, Patrick James
1969-01-01T23:59:59.000Z
) for a large class of stochastic processes. The numerical solution of this problem was approached in the following manner. First of a1 1 equation (1. 1) was replaced by a di fference equation. Then the diff'erence ecuation wa. , solved for P(x, t...'s which would be used for i = n-1 would be those corresponding to r = 5. QP Using (2. 9) as an approximation for ? ', equation (2. 22) results 3x' in the difference equation (2. 10) where Z. is the expression on the right side of ecuation (2. 9...
Numerical Simulations of Bouncing Jets
Bonito, Andrea; Lee, Sanghyun
2015-01-01T23:59:59.000Z
Bouncing jets are fascinating phenomenons occurring under certain conditions when a jet impinges on a free surface. This effect is observed when the fluid is Newtonian and the jet falls in a bath undergoing a solid motion. It occurs also for non-Newtonian fluids when the jets falls in a vessel at rest containing the same fluid. We investigate numerically the impact of the experimental setting and the rheological properties of the fluid on the onset of the bouncing phenomenon. Our investigations show that the occurrence of a thin lubricating layer of air separating the jet and the rest of the liquid is a key factor for the bouncing of the jet to happen. The numerical technique that is used consists of a projection method for the Navier-Stokes system coupled with a level set formulation for the representation of the interface. The space approximation is done with adaptive finite elements. Adaptive refinement is shown to be very important to capture the thin layer of air that is responsible for the bouncing.
Capturing Incomplete Information in Resource Allocation Problems through Numerical
Powell, Warren B.
Capturing Incomplete Information in Resource Allocation Problems through Numerical Patterns Arun information in resource allocation models. We use a popular goodness-of-fit measure known as the Cramer-Von Mises metric based on the empirical distribution function (EDF) as the foundation of our approach. We
Breakdown of analyticity: From rigorous results to numerics
Doedel, Eusebius
work with Rafael de la Llave #12;Computing the breakdown Quasi-Periodic Solutions Outline QuasiBreakdown of analyticity: From rigorous results to numerics Renato Calleja Applied Mathematics-Periodic Solutions Twist Maps Models arising in Statistical Mechanics Computation of the Breakdown Small divisors
Numerical Convergence Properties of Option Pricing PDEs with Uncertain Volatility
Forsyth, Peter A.
is then simply the cost of this hedge. The principal source of risk is the price of the underlying assetNumerical Convergence Properties of Option Pricing PDEs with Uncertain Volatility D. M. Pooley , P, 2001 Abstract The pricing equations derived from uncertain volatility models in finance are often cast
On the numerical solution of some problems of environmental pollution
Ehrhardt, Matthias
world and their importance will even increase in the future. High pollution of air, water and soil mayOn the numerical solution of some problems of environmental pollution Quang A Dang1 , Matthias to the mathematical formula- tion of air pollution models. Let G be a cylindrical domain in the three
Simple intrinsic defects in InAs : numerical predictions.
Schultz, Peter Andrew
2013-03-01T23:59:59.000Z
This Report presents numerical tables summarizing properties of intrinsic defects in indium arsenide, InAs, as computed by density functional theory using semi-local density functionals, intended for use as reference tables for a defect physics package in device models.
A Numerical Investigation of Wettability Alteration during Immiscible CO2
Hossain, M. Enamul
A Numerical Investigation of Wettability Alteration during Immiscible CO2 Flooding Process was devoted to laboratory experiments including core and micro-model flooding to investigate wettability that the inclusion of wettability alteration is very influential parameter in enhancing oil recovery and sweeping
RisR1252(EN) The Numerical Wind Atlas
RisøR1252(EN) The Numerical Wind Atlas -- the KAMM/WAsP Method Helmut P. Frank, Ole Rathmann The method of combining the Karlsruhe Atmospheric Mesoscale Model, KAMM, with the Wind Atlas Analysis and Application Program, WAsP, to make local predictions of the wind resource is presented. It combines
Numerical and experimental investigation of nonlinear properties of rubber absorber
Paris-Sud XI, Université de
Numerical and experimental investigation of nonlinear properties of rubber absorber in rail, Shanghai, 201804 Abstract: A dynamic nonlinear model of rubber absorber in railway fastening is proposed in this paper based on a superposition principal demonstrating that the restoring force of rubber components
Numerical simulation of electrokinetically driven micro flows
Hahm, Jungyoon
2005-11-01T23:59:59.000Z
Spectral element based numerical solvers are developed to simulate electrokinetically driven flows for micro-fluidic applications. Based on these numerical solvers, basic phenomena and devices for electrokinetic applications in micro and nano flows...
An implicit numerical algorithm general relativistic hydrodynamics
A. Hujeirat
2008-01-09T23:59:59.000Z
An implicit numerical algorithm general relativistic hydrodynamics This article has been replaced by arXiv:0801.1017
From Numerical Analysis to Computational Science
Li, Tiejun
. Numerical computing has, of course, been part of mathematics for a very long time. Al- gorithms by the namesFrom Numerical Analysis to Computational Science Bj¨orn Engquist · Gene Golub 1. Introduction The modern development of numerical computing is driven by the rapid in- crease in computer performance
Simulating Reionization in Numerical Cosmology
Aaron Sokasian; Tom Abel; Lars E. Hernquist
2001-05-10T23:59:59.000Z
The incorporation of radiative transfer effects into cosmological hydrodynamical simulations is essential for understanding how the intergalactic medium (IGM) makes the transition from a neutral medium to one that is almost fully ionized. Here, we present an approximate numerical method designed to study in a statistical sense how a cosmological density field is ionized by a set of discrete point sources. A diffuse background radiation field is also computed self-consistently in our procedure. The method requires relatively few time steps and can be employed with simulations having high resolution. We describe the details of the algorithm and provide a description of how the method can be applied to the output from a pre-existing cosmological simulation to study the systematic reionization of a particular ionic species. As a first application, we compute the reionization of He II by quasars in the redshift range 3 to 6.
Turbulent Combustion in Type Ia Supernova Models
F. K. Roepke; W. Hillebrandt
2006-09-15T23:59:59.000Z
We review the astrophysical modeling of type Ia supernova explosions and describe numerical methods to implement numerical simulations of these events. Some results of such simulations are discussed.
Farquharson, Colin G.
Comparison of integral equation and physical scale modelling of the electromagnetic response history of EM numerical modelling in geophysics. Â· Another integral equation modelling program;Introduction: a brief history Â· Two main approaches to numerical modelling: integral equation; finite
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. 2005; 12:683
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. 2005; 12:683 Published Numerical Linear Algebra and its Applications The fourth workshop of the ERCIM Working Group on `Matrix Computations and Statistics' and the First International workshop on `Numerical Linear Algebra and its
Area E Numerics and Scientific Computing Foundation in Numerical Mathematics F4E1
Sturm, Karl-Theodor
Area E Numerics and Scientific Computing Foundation in Numerical Mathematics F4E1: · Scientific Selected Topics in Scientific Computing · V5E3 Advanced Topics in Numerical Methods in Science matrices) · Computational Finance (e.g. option pricing, fast numerical methods) · Visualization Methods (e
Non-Gaussian numerical errors versus mass hierarchy
Y. Meurice; M. B. Oktay
2000-05-12T23:59:59.000Z
We probe the numerical errors made in renormalization group calculations by varying slightly the rescaling factor of the fields and rescaling back in order to get the same (if there were no round-off errors) zero momentum 2-point function (magnetic susceptibility). The actual calculations were performed with Dyson's hierarchical model and a simplified version of it. We compare the distributions of numerical values obtained from a large sample of rescaling factors with the (Gaussian by design) distribution of a random number generator and find significant departures from the Gaussian behavior. In addition, the average value differ (robustly) from the exact answer by a quantity which is of the same order as the standard deviation. We provide a simple model in which the errors made at shorter distance have a larger weight than those made at larger distance. This model explains in part the non-Gaussian features and why the central-limit theorem does not apply.
Direct numerical simulation of turbulent reacting flows
Chen, J.H. [Sandia National Laboratories, Livermore, CA (United States)
1993-12-01T23:59:59.000Z
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 POOL BOILING FOR STEADY STATE AND TRANSIENT HEATING
Maruyama, Shigeo
1 NUMERICAL SIMULATION OF POOL BOILING FOR STEADY STATE AND TRANSIENT HEATING Ying He, Masahiro role in nucleate and transition boiling heat transfer at high heat flux. Many experiments have been in the numerical simulation of boiling heat transfer. In this study, based on the macrolayer evaporation model
Impact of Turbulence Closures and Numerical Errors for the Optimization of Flow Control Devices
Paris-Sud XI, Université de
Impact of Turbulence Closures and Numerical Errors for the Optimization of Flow Control Devices J the use of a Kriging-based global optimization method to determine optimal control parameters conduct an optimization process and measure the impact of numerical and modeling errors on the optimal
Numerical study on optimal Stirling engine regenerator matrix designs taking into account matrix design that improves the efficiency of a Stirling engine has been developed in a numerical study of the existing SM5 Stirling engine. A new, detailed, one-dimensional Stirling engine model that delivers results
NUMERICAL STUDY OF A TURBULENT HYDRAULIC JUMP Qun Zhao 1 Shubhra K. Misra1
Zhao, Qun
. Hydraulic jumps are commonly used as energy dissipators and they have been studied intensively by hydraulicNUMERICAL STUDY OF A TURBULENT HYDRAULIC JUMP Qun Zhao 1 Shubhra K. Misra1 Ib A. Svendsen 1 (Member of a turbulent hydraulic jump. The numerical model is based on RIPPLE (Kothe et al., 1994) with two turbulence
Physical mechanisms involved in grooved flat heat pipes: experimental and numerical analyses
Paris-Sud XI, Université de
1 Physical mechanisms involved in grooved flat heat pipes: experimental and numerical analyses.lefevre@insa-lyon.fr) Abstract An experimental database, obtained with flat plate heat pipes (FPHP) with longitudinal grooves to results of numerical models. Physical mechanisms involved in grooved heat pipes are discussed, including
MATHEMATICAL THEORY AND NUMERICAL METHODS FOR ...
2015-02-05T23:59:59.000Z
Weizhu Bao. Department of Mathematics and Center for Computational Science and Engineering ... Numerical methods for computing dynamics of GPE. 36. 4.1.
Adjoint Sensitivity Analysis for Numerical Weather Prediction
Alexandru Cioaca
2011-09-02T23:59:59.000Z
Sep 2, 2011 ... Adjoint Sensitivity Analysis for Numerical Weather Prediction: ... regional generation cost in the state of Illinois with respect to wind speed and ...
Paris-Sud XI, Université de
-induced snow transport for a high5 alpine area in the Berchtesgaden National Park (Germany), and we discuss sublimation from the surface15 of the trees, and a model for the simulation of wind-induced snow transport heat of fusion (melting or freezing) and the latent heat of vaporization (evaporation or condensation
Boyer, Edmond
uncontrollably if icing occurs on the horizontal or vertical stabilizer. Obstruction of pitot tube will causeSimHydro 2014:Modelling of rapid transitory flows,11-13 June 2014, Sophia Antipolis Marechal E at low temperatures. If the flow rate is increased, sudden releases of large quantities of ice may occur
Numerical Weather Forecasting at the Savannah River Site
Buckley, R.L. [Westinghouse Savannah River Company, AIKEN, SC (United States)
1998-11-01T23:59:59.000Z
Weather forecasts at the Savannah River Site (SRS) are important for applications to emergency response. The fate of accidentally-released radiological materials and toxic chemicals can be determined by providing wind and turbulence input to atmospheric transport models. This operation has been routinely performed at SRS using the WIND System, a system of computer models and monitors which collect data from towers situated throughout the SRS. However, the information provided to these models is spatially homogeneous (in one or two dimensions) with an elementary forecasting capability. This paper discusses the use of an advanced three-dimensional prognostic numerical model to provide space and time-dependent meteorological data for use in the WIND System dispersion models. The extensive meteorological data collection at SRS serves as a ground truth for further model development as well as for use in other applications.
Numerical solution of plasma fluid equations using locally refined grids
Colella, P., LLNL
1997-01-26T23:59:59.000Z
This paper describes a numerical method for the solution of plasma fluid equations on block-structured, locally refined grids. The plasma under consideration is typical of those used for the processing of semiconductors. The governing equations consist of a drift-diffusion model of the electrons and an isothermal model of the ions coupled by Poisson's equation. A discretization of the equations is given for a uniform spatial grid, and a time-split integration scheme is developed. The algorithm is then extended to accommodate locally refined grids. This extension involves the advancement of the discrete system on a hierarchy of levels, each of which represents a degree of refinement, together with synchronization steps to ensure consistency across levels. A brief discussion of a software implementation is followed by a presentation of numerical results.
Numerical analysis of a microwave torch with axial gas injection
Gritsinin, S. I.; Davydov, A. M.; Kossyi, I. A., E-mail: kossyi@fpl.gpi.ru [Russian Academy of Sciences, Prokhorov General Physics Institute (Russian Federation); Kulumbaev, E. B. [National Research University Belgorod State University (Russian Federation)] [National Research University Belgorod State University (Russian Federation); Lelevkin, V. M. [Kyrgyz-Russian Slavic University (Kyrgyzstan)] [Kyrgyz-Russian Slavic University (Kyrgyzstan)
2013-07-15T23:59:59.000Z
The characteristics of a microwave discharge in an argon jet injected axially into a coaxial channel with a shortened inner electrode are numerically analyzed using a self-consistent equilibrium gas-dynamic model. The specific features of the excitation and maintenance of the microwave discharge are determined, and the dependences of the discharge characteristics on the supplied electromagnetic power and gas flow rate are obtained. The calculated results are compared with experimental data.
Pham, Van Sang
We present a systematic, multiscale, fully detailed numerical modeling for dynamics of fluid flow and ion transport covering Ohmic, limiting, and overlimiting current regimes in conductance of ion-selective membrane. By ...
t -software package for numerical simulations of radioactive contaminant transport in groundwater
Frolkovic, Peter
in groundwater Peter Frolkovic Michael Lampe Gabriel Wittum September 16, 2005 Abstract The software package r3t or as numerical solution of some groundwater flow model, e.g., the density driven flow problem. The matrix Di = Di
On the numerical simulation of the instationary free fall of a solid in a fluid.
from meteorology, sedimentology and aerospace engineering to biology. We present a new numerical method], sedimentology [1], aerospace engineering [21] and biology (e.g. models for animal flights [2]). Recently
Three-dimensional numerical manifold method simulations for blocky rock analysis
Shentu, Longfei
2011-01-01T23:59:59.000Z
After decades of development, people realize that there are wider and more various applications of numerical modeling and analysis. However, current feasible software tools cannot satisfy engineering and commercial goals. ...
Nonclassicality thresholds for multiqubit states: Numerical analysis
Gruca, Jacek; Zukowski, Marek [Institute of Theoretical Physics and Astrophysics, University of Gdansk, PL-80-952 Gdansk (Poland); Laskowski, Wieslaw [Institute of Theoretical Physics and Astrophysics, University of Gdansk, PL-80-952 Gdansk (Poland); Fakultaet fuer Physik, Ludwig-Maximilians Universitaet Muenchen, D-80799 Muenchen (Germany); Max-Planck-Institut fuer Quantenoptik, D-85748 Garching (Germany); Kiesel, Nikolai [Faculty of Physics, University of Vienna, A-1090 Vienna (Austria); Wieczorek, Witlef; Weinfurter, Harald [Fakultaet fuer Physik, Ludwig-Maximilians Universitaet Muenchen, D-80799 Muenchen (Germany); Max-Planck-Institut fuer Quantenoptik, D-85748 Garching (Germany); Schmid, Christian [Fakultaet fuer Physik, Ludwig-Maximilians Universitaet Muenchen, D-80799 Muenchen (Germany); Max-Planck-Institut fuer Quantenoptik, D-85748 Garching (Germany); European Organisation for Astronomical Research in the Southern Hemisphere, D-85748 Garching (Germany)
2010-07-15T23:59:59.000Z
States that strongly violate Bell's inequalities are required in many quantum-informational protocols as, for example, in cryptography, secret sharing, and the reduction of communication complexity. We investigate families of such states with a numerical method which allows us to reveal nonclassicality even without direct knowledge of Bell's inequalities for the given problem. An extensive set of numerical results is presented and discussed.
Numerical Study of a Turbulent Hydraulic Jump
Zhao, Qun
Numerical Study of a Turbulent Hydraulic Jump Qun Zhao, Shubhra Misra, Ib. A. Svendsen and James T of a Turbulent Hydraulic Jump p.1/14 #12;Objective Our ultimate goal is to study the breaking waves. Numerical Study of a Turbulent Hydraulic Jump p.2/14 #12;A moving bore Qiantang Bore China (Courtesy of Dr J
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. 2011; 18:961980
De Sterck, Hans
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. 2011; 18:961980 Published aggregation method of [1], and the Markov chain algebraic multigrid (MCAMG) method of [3] using the OTF
Numerical Methods of Computational Electromagnetics for Complex Inhomogeneous Systems
Cai, Wei
2014-05-15T23:59:59.000Z
Understanding electromagnetic phenomena is the key in many scientific investigation and engineering designs such as solar cell designs, studying biological ion channels for diseases, and creating clean fusion energies, among other things. The objectives of the project are to develop high order numerical methods to simulate evanescent electromagnetic waves occurring in plasmon solar cells and biological ion-channels, where local field enhancement within random media in the former and long range electrostatic interactions in the latter are of major challenges for accurate and efficient numerical computations. We have accomplished these objectives by developing high order numerical methods for solving Maxwell equations such as high order finite element basis for discontinuous Galerkin methods, well-conditioned Nedelec edge element method, divergence free finite element basis for MHD, and fast integral equation methods for layered media. These methods can be used to model the complex local field enhancement in plasmon solar cells. On the other hand, to treat long range electrostatic interaction in ion channels, we have developed image charge based method for a hybrid model in combining atomistic electrostatics and continuum Poisson-Boltzmann electrostatics. Such a hybrid model will speed up the molecular dynamics simulation of transport in biological ion-channels.
West, Matthew
momentum and energy conservation characteristics, as several numerical examples demonstrate. Copyright 2005INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2005; 64Center for Advanced Computing Research, California Institute of Technology, Mail Code 158-79, Pasadena
NUMERICAL SIMULATIONS OF CHROMOSPHERIC MICROFLARES
Jiang, R. L.; Fang, C.; Chen, P. F., E-mail: fangc@nju.edu.c [Department of Astronomy, Nanjing University, Nanjing 210093 (China)
2010-02-20T23:59:59.000Z
With gravity, ionization, and radiation being considered, we perform 2.5 dimensional (2.5D) compressible resistive magnetohydrodynamic (MHD) simulations of chromospheric magnetic reconnection using the CIP-MOCCT scheme. The temperature distribution of the quiet-Sun atmospheric model VALC and the helium abundance (10%) are adopted. Our 2.5D MHD simulation reproduces qualitatively the temperature enhancement observed in chromospheric microflares. The temperature enhancement DELTAT is demonstrated to be sensitive to the background magnetic field, whereas the total evolution time DELTAt is sensitive to the magnitude of the anomalous resistivity. Moreover, we found a scaling law, which is described as DELTAT/DELTAt {approx} n{sub H} {sup -1.5} B {sup 2.1}eta{sub 0} {sup 0.88}. Our results also indicate that the velocity of the upward jet is much greater than that of the downward jet, and the X-point may move up or down.
Connecting Numerical Relativity and Data Analysis of Gravitational Wave Detectors
Deirdre Shoemaker; Karan Jani; Lionel London; Larne Pekowsky
2015-03-09T23:59:59.000Z
Gravitational waves deliver information in exquisite detail about astrophysical phenomena, among them the collision of two black holes, a system completely invisible to the eyes of electromagnetic telescopes. Models that predict gravitational wave signals from likely sources are crucial for the success of this endeavor. Modeling binary black hole sources of gravitational radiation requires solving the Eintein equations of General Relativity using powerful computer hardware and sophisticated numerical algorithms. This proceeding presents where we are in understanding ground-based gravitational waves resulting from the merger of black holes and the implications of these sources for the advent of gravitational-wave astronomy.
Numerical and experimental analyses of resin infusion manufacturing processes of composite materials
Wang, Peng; Molimard, Jérôme; Vautrin, Alain; Minni, Jean-Christophe; 10.1177/0021998311421990
2012-01-01T23:59:59.000Z
Liquid resin infusion (LRI) processes are promising manufacturing routes to produce large, thick, or complex structural parts. They are based on the resin flow induced, across its thickness, by a pressure applied onto a preform/resin stacking. However, both thickness and fiber volume fraction of the final piece are not well controlled since they result from complex mechanisms which drive the transient mechanical equilibrium leading to the final geometrical configuration. In order to optimize both design and manufacturing parameters, but also to monitor the LRI process, an isothermal numerical model has been developed which describes the mechanical interaction between the deformations of the porous medium and the resin flow during infusion.1, 2 With this numerical model, it is possible to investigate the LRI process of classical industrial part shapes. To validate the numerical model, first in 2D, and to improve the knowledge of the LRI process, this study details a comparison between numerical simulations and...
The Numerical Simulation Of A Transitional Flow In The VKI-GENOA Turbine Cascade
Yershov, Sergiy; Yakovlev, Viktor; Gryzun, Maria
2015-01-01T23:59:59.000Z
This study presents a numerical simulation of a 3D viscous flow in the VKI-Genoa cascade that takes into account the laminar-turbulent transition. The numerical simulation is performed using the Reynolds-averaged Navier-Stokes equations and the two-equation k-omega SST turbulence model. The algebraic Production Term Modification model is used for modeling the laminar-turbulent transition. Computations of both fully turbulent and transitional flows are carried out. The contours of the Mach number, the turbulence kinetic energy, the entropy function, as well as limiting streamlines are presented. The analysis of the numerical results demonstrates the influence of the laminar-turbulent transition on the secondary flow pattern. The comparison between the present computational results and the existing experimental and numerical data shows that the proposed approach reflects sufficiently the physics of the laminar-turbulent transition in turbine cascades.
Numerical weather forecasting at the Savannah River Site
Buckley, R.L. [Westinghouse Savannah River Site, Aiken, SC (United States)
1998-12-31T23:59:59.000Z
Weather forecasts at the Savannah River Site (SRS) are important for applications to emergency response. The fate of accidentally released radiological materials and toxic chemicals can be determined by providing wind and turbulence input to atmospheric transport models. This operation has been routinely performed at SRS using the WIND system, a system of computer models and monitors that collects data from towers situated throughout the SRS. However, the information provided to these models is spatially homogeneous (in one or two dimensions) with an elementary forecasting capability. This paper discusses the use of an advanced three-dimensional prognostic numerical model to provide space- and time-dependent meteorological data for use in the WIND system dispersion models. The extensive meteorological data collection at SRS serves as a ground truth for further model development as well as for use in other applications. A prognostic mesoscale model, the regional atmospheric modeling system (RAMS), is used to provide these forecasts. Use of RAMS allows for incorporation of mesoscale features such as the sea breeze, which has been shown to affect local weather conditions. This paper discusses the mesoscale model and its configuration for the operational simulation, as well as an application using a dispersion model at the SRS.
Multiscale Problems: Numerical Analysis and Scientific Computing
Wirosoetisno, Djoko
/Transport Equations Alternative (Wind, Solar, Wave, . . . ) Weather and Climate Prediction Wave Energy Fuel CellsMultiscale Problems: Numerical Analysis and Scientific Computing with Applications in Energy inverse problems, data assimilation, stochastic differential equations Networks, compressed sensing
A destructuration theory and its application to SANICLAY model
Taiebat, Mahdi; Dafalias, Yannis; Peek, Ralf
2010-01-01T23:59:59.000Z
AND ANALYTICAL METHODS IN GEOMECHANICS Int. J. Numer. Anal.a¨ a Numerical Models in Geomechanics, Proceedings of NUMOGAnalytical Methods in Geomechanics 2006; 30(12):1231–1257.
Numerical determination of the material properties of porous dust cakes
D. Paszun; C. Dominik
2008-02-13T23:59:59.000Z
The formation of planetesimals requires the growth of dust particles through collisions. Micron-sized particles must grow by many orders of magnitude in mass. In order to understand and model the processes during this growth, the mechanical properties, and the interaction cross sections of aggregates with surrounding gas must be well understood. Recent advances in experimental (laboratory) studies now provide the background for pushing numerical aggregate models onto a new level. We present the calibration of a previously tested model of aggregate dynamics. We use plastic deformation of surface asperities as the physical model to bring critical velocities for sticking into accordance with experimental results. The modified code is then used to compute compression strength and the velocity of sound in the aggregate at different densities. We compare these predictions with experimental results and conclude that the new code is capable of studying the properties of small aggregates.
Spherically symmetric cosmological spacetimes with dust and radiation — numerical implementation
Lim, Woei Chet [Department of Mathematics, University of Waikato, Private Bag 3105, Hamilton 3240 (New Zealand); Regis, Marco [Dipartimento di Fisica, Università di Torino and INFN, Torino (Italy); Clarkson, Chris, E-mail: wclim@waikato.ac.nz, E-mail: regis@to.infn.it, E-mail: chris.clarkson@gmail.com [Astrophysics, Cosmology and Gravity Centre, and Department of Mathematics and Applied Mathematics, University of Cape Town, Rondebosch 7701, Cape Town (South Africa)
2013-10-01T23:59:59.000Z
We present new numerical cosmological solutions of the Einstein Field Equations. The spacetime is spherically symmetric with a source of dust and radiation approximated as a perfect fluid. The dust and radiation are necessarily non-comoving due to the inhomogeneity of the spacetime. Such a model can be used to investigate non-linear general relativistic effects present during decoupling or big-bang nucleosynthesis, as well as for investigating void models of dark energy with isocurvature degrees of freedom. We describe the full evolution of the spacetime as well as the redshift and luminosity distance for a central observer. After demonstrating accuracy of the code, we consider a few example models, and demonstrate the sensitivity of the late time model to the degree of inhomogeneity of the initial radiation contrast.
Numerical Estimation of Frictional Torques with Rate and State Friction
Arun K. Singh; T. N. Singh
2015-01-20T23:59:59.000Z
In this paper, numerical estimation of frictional torques is carried out of a rotary elastic disc on a hard and rough surface under different rotating conditions. A one dimensional spring- mass rotary system is numerically solved under the quasistatic condition with the rate and state dependent friction model. It is established that torque of frictional strength as well as torque of steady dynamic stress increases with radius and found to be maximum at the periphery of the disc. Torque corresponding to frictional strength estimated using the analytical solution matches closely with the simulation only in the case of high stiffness of the connecting spring. In steady relaxation simulation, a steadily rotating disc is suddenly stopped and relaxational angular velocity and corresponding frictional torque decreases with both steady angular velocity and stiffness of the connecting spring in the velocity strengthening regime. In velocity weakening regime, in contrast, torque of relaxation stress deceases but relaxation velocity increases. The reason for the contradiction is explained.
Numerical method for shear bands in ductile metal with inclusions
Plohr, Jee Yeon N [Los Alamos National Laboratory; Plohr, Bradley J [Los Alamos National Laboratory
2010-01-01T23:59:59.000Z
A numerical method for mesoscale simulation of high strain-rate loading of ductile metal containing inclusions is described. Because of small-scale inhomogeneities, such a composite material is prone to localized shear deformation (adiabatic shear bands). The modeling framework is the Generalized Method of Cells of Paley and Aboudi [Mech. Materials, vol. 14, pp. /27-139, 1992], which ensures that the micromechanical response of the material is reflected in the behavior of the composite at the mesoscale. To calculate the effective plastic strain rate when shear bands are present, the analytic and numerical analysis of shear bands by Glimm, Plohr, and Sharp [Mech. Materials, vol. 24, pp. 31-41, 1996] is adapted and extended.
Accurate complex scaling of three dimensional numerical potentials
Cerioni, Alessandro [European Synchrotron Radiation Facility, 6 rue Horowitz, BP220 38043 Grenoble Cedex 9 (France); Genovese, Luigi; Duchemin, Ivan; Deutsch, Thierry [Laboratoire de simulation atomistique (L-Sim), SP2M, UMR-E CEA/UJF-Grenoble 1, INAC, Grenoble F-38054 (France)
2013-05-28T23:59:59.000Z
The complex scaling method, which consists in continuing spatial coordinates into the complex plane, is a well-established method that allows to compute resonant eigenfunctions of the time-independent Schroedinger operator. Whenever it is desirable to apply the complex scaling to investigate resonances in physical systems defined on numerical discrete grids, the most direct approach relies on the application of a similarity transformation to the original, unscaled Hamiltonian. We show that such an approach can be conveniently implemented in the Daubechies wavelet basis set, featuring a very promising level of generality, high accuracy, and no need for artificial convergence parameters. Complex scaling of three dimensional numerical potentials can be efficiently and accurately performed. By carrying out an illustrative resonant state computation in the case of a one-dimensional model potential, we then show that our wavelet-based approach may disclose new exciting opportunities in the field of computational non-Hermitian quantum mechanics.
Toward standard testbeds for numerical relativity
Miguel Alcubierre; Gabrielle Allen; Carles Bona; David Fiske; Tom Goodale; F. Siddharta Guzman; Ian Hawke; Scott H. Hawley; Sascha Husa; Michael Koppitz; Christiane Lechner; Denis Pollney; David Rideout; Marcelo Salgado; Erik Schnetter; Edward Seidel; Hisa-aki Shinkai; Bela Szilagyi; Deirdre Shoemaker; Ryoji Takahashi; Jeffrey Winicour
2003-05-06T23:59:59.000Z
In recent years, many different numerical evolution schemes for Einstein's equations have been proposed to address stability and accuracy problems that have plagued the numerical relativity community for decades. Some of these approaches have been tested on different spacetimes, and conclusions have been drawn based on these tests. However, differences in results originate from many sources, including not only formulations of the equations, but also gauges, boundary conditions, numerical methods, and so on. We propose to build up a suite of standardized testbeds for comparing approaches to the numerical evolution of Einstein's equations that are designed to both probe their strengths and weaknesses and to separate out different effects, and their causes, seen in the results. We discuss general design principles of suitable testbeds, and we present an initial round of simple tests with periodic boundary conditions. This is a pivotal first step toward building a suite of testbeds to serve the numerical relativists and researchers from related fields who wish to assess the capabilities of numerical relativity codes. We present some examples of how these tests can be quite effective in revealing various limitations of different approaches, and illustrating their differences. The tests are presently limited to vacuum spacetimes, can be run on modest computational resources, and can be used with many different approaches used in the relativity community.
Rutqvist, Jonny; Majer, Ernie; Oldenburg, Curt; Peterson, John; Vasco, Don
2006-01-01T23:59:59.000Z
coupled reservoir geomechanical numerical modeling, (2) datacoupled reservoir geomechanical numerical modeling, (2) datareservoir geomechanical analysis will be corroborated with, and constrained by, data
On the Numerical Dispersion of Electromagnetic Particle-In-Cell Code : Finite Grid Instability
Meyers, Michael David [Los Alamos National Lab. (LANL), Los Alamos, NM (United States); Univ. of California, Los Angeles, CA (United States) Dept. of Physics and Astronomy; Huang, Chengkun [Los Alamos National Lab. (LANL), Los Alamos, NM (United States); Zeng, Yong [Los Alamos National Lab. (LANL), Los Alamos, NM (United States); Yi, Sunghwan [Los Alamos National Lab. (LANL), Los Alamos, NM (United States); Albright, Brian James [Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
2014-07-15T23:59:59.000Z
The Particle-In-Cell (PIC) method is widely used in relativistic particle beam and laser plasma modeling. However, the PIC method exhibits numerical instabilities that can render unphysical simulation results or even destroy the simulation. For electromagnetic relativistic beam and plasma modeling, the most relevant numerical instabilities are the finite grid instability and the numerical Cherenkov instability. We review the numerical dispersion relation of the electromagnetic PIC algorithm to analyze the origin of these instabilities. We rigorously derive the faithful 3D numerical dispersion of the PIC algorithm, and then specialize to the Yee FDTD scheme. In particular, we account for the manner in which the PIC algorithm updates and samples the fields and distribution function. Temporal and spatial phase factors from solving Maxwell's equations on the Yee grid with the leapfrog scheme are also explicitly accounted for. Numerical solutions to the electrostatic-like modes in the 1D dispersion relation for a cold drifting plasma are obtained for parameters of interest. In the succeeding analysis, we investigate how the finite grid instability arises from the interaction of the numerical 1D modes admitted in the system and their aliases. The most significant interaction is due critically to the correct representation of the operators in the dispersion relation. We obtain a simple analytic expression for the peak growth rate due to this interaction.
Provably Efficient Algorithms for Numerical Tensor Edgar Solomonik
California at Berkeley, University of
Provably Efficient Algorithms for Numerical Tensor Algebra Edgar Solomonik Electrical Engineering for Numerical Tensor Algebra by Edgar Solomonik A dissertation submitted in partial satisfaction-Gordon Fall 2014 #12;Provably Efficient Algorithms for Numerical Tensor Algebra Copyright 2014 by Edgar
Numerical Simulation on Laser Fusion in China
Zhu Shaoping; Pei Wenbing; Xu Yan; Gu Peijun; Lan Ke; Ye Wenhua; Wu Junfeng; Li Jinghong; Gao Yaoming; Zheng Chunyang; Li Shuanggui; Mo Zeyao; Yan Jun [Institute of Applied Physics and Computational Mathematics, Beijing 100088 (China); Zhang Weiyan [National High-Tech Inertial Confinement Fusion Committee of China, Beijing 100088 (China)
2009-05-02T23:59:59.000Z
Numerical simulation is a powerful tool to get insight into the physics of laser fusion. Much effort has been devoted to develop the numerical simulation code series named LARED in China. The code series LARED are composed of six parts and enable us to have the simulation capability for the key processes in laser fusion. In recent years, a number of numerical simulations using LARED have been carried out and the simulation is checked by experiments done at the laser facility SG-II and SG-III prototype. In the present talk, some details of LARED code series will be introduced, and some simulation results, especially recent work on the opacities, will be shown.
KWIC Index for Numerical Linear Algebra
Carpenter, J.A.
1983-07-01T23:59:59.000Z
This report is a sequel to ORNL/CSD-106 in the ongoing supplements to Professor A.S. Householder's KWIC Index for Numerical Algebra. Beginning with the previous supplement, the subject has been restricted to Numerical Linear Algebra, roughly characterized by the American Mathematical Society's classification sections 15 and 65F but with little coverage of infinite matrices, matrices over fields of characteristics other than zero, operator theory, optimization and those parts of matrix theory primarily combinatorial in nature. Some consideration is given to the uses of graph theory in Numerical Linear Algebra, particularly with respect to algorithms for sparse matrix computations. The period covered by this report is roughly the calendar year 1982 as measured by the appearance of the articles in the American Mathematical Society's Contents of Mathematical Publications lagging actual appearance dates by up to nearly half a year. The review citations are limited to the Mathematical Reviews (MR).
Analytical-Numerical Modeling Of Komatiite Lava Emplacement And...
thicknesses over felsic tuffaceous substrates with various degrees of consolidation and water contents. Thermo-mechanical erosion becomes more effective for substrates that are...
Numerical Modeling of HCCI Combustion | Department of Energy
Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site
AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on DeliciousPlasma | Department ofEnergy 9 LehmanDepartment of Energy Nuclear Safety Regulatory83-2007iEnergyHCCI
Numerical Modeling of PCCI Combustion | Department of Energy
Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site
AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on DeliciousPlasma | Department ofEnergy 9 LehmanDepartment of Energy Nuclear Safety
Numerical modeling of magnetohydrodynamic activity in the Swarthmore Spheromak Experiment
Brown, Michael R.
simulation data will be used as input for a particle orbit and energization code. © 2001 American Institute of Physics. DOI: 10.1063/1.1362294 I. OVERVIEW There are several unresolved mysteries related to solar and astrophysical magnetism. First is the origin of solar and astrophysical magnetic fields: the dynamo problem
Numerical modeling of seawater intrusion into endorheic hydrological systems
Lyakhovsky, Vladimir
concentration to brine through evaporation, the dense brine starts to flow below the encroaching seawater body. Below it an interface is recognized between the encroaching seawater body and a brine body underneath
A Numeric Predictive Failure Model for Percussive Excavation
Green, Alex Nicholas
2011-01-01T23:59:59.000Z
Densities, 500, 1250, and 1750 BPM, 70 mm Depth, 5 mm/s, ODensities, 0,500, and 1750 BPM, 70 mm Depth, 5 mm/s, 90 Ofrom 4.16 Hz to 29.16 Hz (250 BPM – 1750 BPM). In order to
Challenges for Numerical Modeling of Enhanced Geothermal Systems
J. P. Fairley; S. E. Ingebritsen; R. K. Podgorney
2010-08-01T23:59:59.000Z
A recent guest editorial by Wood (2009) pointed out the potential of enhanced geothermal systems (EGS) as a future source of “green” energy and suggested that EGS offers research opportunities for hydrogeologists seeking to become involved in the world’s energy future. Although EGS may have a bright future as a sustainable, low-carbon emission energy source, significant technical challenges must be overcome before this promising energy resource can be commercially viable. Because pilot EGS projects in the United States face very different economic constraints than current European projects, there is a real need to make technological advances to improve the return on capital investment. In this article, we amplify on Wood’s excellent editorial by describing some of the challenges that exist for the simulation of EGS.
Bidimensional Numerical Model for Polyurethane Smoldering in a Fixed Bed
Ghabi, Chekib; Rein, Guillermo; Ben Ticha, Hmaied; Sassi, Mohamed
2005-01-01T23:59:59.000Z
FOR POLYUERETHANE SMOLDERING IN A FIXED BED Chekib GHABI 1 ,experiments. In this paper, a fixed bed of polyurethane foam
Geology 775 NUMERICAL MODELING FOR EARTH SCIENCE APPLICATIONS
and engineers. Brian D. Hahn, 2002. [On reserve] Matlab, Second Edition: A Practical Introduction to Programming with Applications. Amos Gilat, 2010. [The 2008 edition is pretty reasonable on Amazon] Course Goals and Learning to address them 3. Write computer programs that use finite difference and finite element methods to solve
Geology 575 NUMERICAL MODELING FOR EARTH SCIENCE APPLICATIONS
and engineers. Brian D. Hahn, 2002. [On reserve] Matlab, Second Edition: A Practical Introduction to Programming with Applications. Amos Gilat, 2010. [The 2008 edition is pretty reasonable on Amazon] Course Goals and Learning to address them 3. Write computer programs that use finite difference and finite element methods to solve
Numerical modelling of current transfer in nonlinear anisotropic conductive media
Baranowski, Robert Paul
on the nature of current transport. The main motivation for this work was the desire for a better understanding of the conceptually difficult behaviour of current transport in superconducting bodies and examines current transfer quantitatively for a number...
Simplified Numerical Models for Complex Air Supply Jelena Srebric*
Chen, Qingyan "Yan"
greatly improve our knowledge of the performance of ventilation systems. Unfortunately, a major bottleneck
Modeling and numerical approximation of two-phase ...
2011-02-09T23:59:59.000Z
Feb 9, 2011 ... phase incompressible flows with matching or non-matching density. ... The interfacial dynamics of two-phase immiscible fluids have attracted.
1-Dimensional Numerical Model of Thermal Conduction and Vapor Diffusion
SchÃ¶rghofer, Norbert
developed by Samar Khatiwala, 2001 extended to variable thermal properties and irregular grid by Norbert Sch for c. Upper boundary condition: a) Radiation Q + k T z z=0 = T4 z=0 Q is the incoming solar flux of Water Vapor with Phase Transitions developed by Norbert SchÂ¨orghofer, 2003Â2004 3 phases: vapor, free
Numerical model for steel catenary riser on seafloor support
You, Jung Hwan
2007-04-25T23:59:59.000Z
on catenary riser and stresses in riser pipe at the touch down point (TDP) was conducted over 3 months at Watchet Harbor in the west of England by the STRIDE ? JIP, 2H Offshore Engineering Ltd in 2000 (Willis and West, 2001). The purpose of the full... published data and data from the pipe and soil interaction experiments conducted within the STRIDE and CARISIMA JIP?s. They describe an example of the development of a pipe and soil interaction curve with an unloading and reloading cycle, as presented...
Numerical and analytical modeling of sanding onset prediction
Yi, Xianjie
2004-09-30T23:59:59.000Z
To provide technical support for sand control decision-making, it is necessary to predict the production condition at which sand production occurs. Sanding onset prediction involves simulating the stress state on the surface of an oil/gas producing...
A UNIFIED NUMERICAL MODEL FOR SATURATED-UNSATURATED GROUNDWATER FLOW
Narasimhan, T.N.
2011-01-01T23:59:59.000Z
saturation) 1/L c c Compression Index; slope of the best-line is called the compression index (C ) in the case of thestraightline is the compression index, C of c Constant air
A UNIFIED NUMERICAL MODEL FOR SATURATED-UNSATURATED GROUNDWATER FLOW
Narasimhan, T.N.
2011-01-01T23:59:59.000Z
Saturated-Unsaturated Groundwater Flow Ph.D. Dissertation in~ " Fundamental principles of groundwater flow uv e in Flowunsaturated flow in a groundwater basi.n 11 9 Hater
Mathematical Theory and Numerical Methods for Biomolecular Modeling /
White, Michael Robert
2013-01-01T23:59:59.000Z
50] B. Li. , X. -L. Cheng, and Z. -F. Zhang. Dielectric127:084503, 2007. [14] L. -T. Cheng, Z. Wang, P. Setny, J.63] P. Setny, Z. Wang, L. -T. Cheng, B. Li, J. A. McCammon,
Numerical modeling of methane venting from lake sediments
Scandella, Benjamin P. (Benjamin Paul)
2010-01-01T23:59:59.000Z
The dynamics of methane transport in lake sediments control the release of methane into the water column above, and the portion that reaches the atmosphere may contribute significantly to the greenhouse effect. The observed ...
Numerical Modeling of Transient Basin and Range Extensional Geothermal...
km relief fromthe valley floor) separated by a thick sequence (about4 km) of clastic sediments derived from the adjacentranges, and a relatively permeable, high angle faultthat...
Hydrodynamic and numerical modeling of a spherical homogeneous.pdf
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE:1 First Use of Energy for All Purposes (Fuel and Nonfuel), 2002; Level:5 TablesExports(Journal Article)41clothThe Bonneville PowerTariff Pages default Sign In About | Careers | ContactStatistical Self-Similarity in Peak
Analytical-Numerical Modeling Of Komatiite Lava Emplacement And Thermal
AFDC Printable Version Share this resource Send a link to EERE: Alternative Fuels Data Center Home Page to someone by E-mail Share EERE: Alternative Fuels Data Center Home Page on Facebook Tweet about EERE: Alternative Fuels Data Center Home Page on Twitter Bookmark EERE: Alternative Fuels Data Center Home Page on Google Bookmark EERE: Alternative Fuels Data Center Home Page on DeliciousPlasmaP a gHigh4-FD-a <Alaska UndergroundInformation An EmpiricalGeothermal-Energy| Open1986) |
Modeling of fluids and waves with analytics and numerics
Liang, Xiangdong, Ph. D. Massachusetts Institute of Technology
2013-01-01T23:59:59.000Z
Capillary instability (Plateau-Rayleigh instability) has been playing an important role in experimental work such as multimaterial fiber drawing and multilayer particle fabrication. Motivated by complex multi-fluid geometries ...