Modeling Crack Propagation in Polycrystalline Microstructure Using Variational Multiscale Method
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Sun, S.; Sundararaghavan, V.
2016-01-01
Crack propagation in a polycrystalline microstructure is analyzed using a novel multiscale model. The model includes an explicit microstructural representation at critical regions (stress concentrators such as notches and cracks) and a reduced order model that statistically captures the microstructure at regions far away from stress concentrations. Crack propagation is modeled in these critical regions using the variational multiscale method. In this approach, a discontinuous displacement field is added to elements that exceed the critical values of normal or tangential tractions during loading. Compared to traditional cohesive zone modeling approaches, the method does not require the use of any specialmore » interface elements in the microstructure and thus can model arbitrary crack paths. The capability of the method in predicting both intergranular and transgranular failure modes in an elastoplastic polycrystal is demonstrated under tensile and three-point bending loads.« less
Multiscale Subsurface Biogeochemical Modeling
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Biogeochemical Modeling Multiscale Subsurface Biogeochemical Modeling ScheibeSmaller.jpg Simulation of flow inside an experimental packed bed, performed on Franklin Key...
Peridynamic Multiscale Finite Element Methods
Costa, Timothy; Bond, Stephen D.; Littlewood, David John; Moore, Stan Gerald
2015-12-01
The problem of computing quantum-accurate design-scale solutions to mechanics problems is rich with applications and serves as the background to modern multiscale science research. The prob- lem can be broken into component problems comprised of communicating across adjacent scales, which when strung together create a pipeline for information to travel from quantum scales to design scales. Traditionally, this involves connections between a) quantum electronic structure calculations and molecular dynamics and between b) molecular dynamics and local partial differ- ential equation models at the design scale. The second step, b), is particularly challenging since the appropriate scales of molecular dynamic and local partial differential equation models do not overlap. The peridynamic model for continuum mechanics provides an advantage in this endeavor, as the basic equations of peridynamics are valid at a wide range of scales limiting from the classical partial differential equation models valid at the design scale to the scale of molecular dynamics. In this work we focus on the development of multiscale finite element methods for the peridynamic model, in an effort to create a mathematically consistent channel for microscale information to travel from the upper limits of the molecular dynamics scale to the design scale. In particular, we first develop a Nonlocal Multiscale Finite Element Method which solves the peridynamic model at multiple scales to include microscale information at the coarse-scale. We then consider a method that solves a fine-scale peridynamic model to build element-support basis functions for a coarse- scale local partial differential equation model, called the Mixed Locality Multiscale Finite Element Method. Given decades of research and development into finite element codes for the local partial differential equation models of continuum mechanics there is a strong desire to couple local and nonlocal models to leverage the speed and state of the
Multiscale Thermohydrologic Model
T. Buscheck
2004-10-12
The purpose of the multiscale thermohydrologic model (MSTHM) is to predict the possible range of thermal-hydrologic conditions, resulting from uncertainty and variability, in the repository emplacement drifts, including the invert, and in the adjoining host rock for the repository at Yucca Mountain. Thus, the goal is to predict the range of possible thermal-hydrologic conditions across the repository; this is quite different from predicting a single expected thermal-hydrologic response. The MSTHM calculates the following thermal-hydrologic parameters: temperature, relative humidity, liquid-phase saturation, evaporation rate, air-mass fraction, gas-phase pressure, capillary pressure, and liquid- and gas-phase fluxes (Table 1-1). These thermal-hydrologic parameters are required to support ''Total System Performance Assessment (TSPA) Model/Analysis for the License Application'' (BSC 2004 [DIRS 168504]). The thermal-hydrologic parameters are determined as a function of position along each of the emplacement drifts and as a function of waste package type. These parameters are determined at various reference locations within the emplacement drifts, including the waste package and drip-shield surfaces and in the invert. The parameters are also determined at various defined locations in the adjoining host rock. The MSTHM uses data obtained from the data tracking numbers (DTNs) listed in Table 4.1-1. The majority of those DTNs were generated from the following analyses and model reports: (1) ''UZ Flow Model and Submodels'' (BSC 2004 [DIRS 169861]); (2) ''Development of Numerical Grids for UZ Flow and Transport Modeling'' (BSC 2004); (3) ''Calibrated Properties Model'' (BSC 2004 [DIRS 169857]); (4) ''Thermal Conductivity of the Potential Repository Horizon'' (BSC 2004 [DIRS 169854]); (5) ''Thermal Conductivity of the Non-Repository Lithostratigraphic Layers'' (BSC 2004 [DIRS 170033]); (6) ''Ventilation Model and Analysis Report'' (BSC 2004 [DIRS 169862]); (7) ''Heat Capacity
Wagner, Gregory John; Collis, Samuel Scott; Templeton, Jeremy Alan; Lehoucq, Richard B.; Parks, Michael L.; Jones, Reese E.; Silling, Stewart Andrew; Scovazzi, Guglielmo; Bochev, Pavel B.
2007-10-01
This report is a collection of documents written as part of the Laboratory Directed Research and Development (LDRD) project A Mathematical Framework for Multiscale Science and Engineering: The Variational Multiscale Method and Interscale Transfer Operators. We present developments in two categories of multiscale mathematics and analysis. The first, continuum-to-continuum (CtC) multiscale, includes problems that allow application of the same continuum model at all scales with the primary barrier to simulation being computing resources. The second, atomistic-to-continuum (AtC) multiscale, represents applications where detailed physics at the atomistic or molecular level must be simulated to resolve the small scales, but the effect on and coupling to the continuum level is frequently unclear.
MULTISCALE THERMOHYDROLOGIC MODEL
T. Buscheck
2005-07-07
The intended purpose of the multiscale thermohydrologic model (MSTHM) is to predict the possible range of thermal-hydrologic conditions, resulting from uncertainty and variability, in the repository emplacement drifts, including the invert, and in the adjoining host rock for the repository at Yucca Mountain. The goal of the MSTHM is to predict a reasonable range of possible thermal-hydrologic conditions within the emplacement drift. To be reasonable, this range includes the influence of waste-package-to-waste-package heat output variability relevant to the license application design, as well as the influence of uncertainty and variability in the geologic and hydrologic conditions relevant to predicting the thermal-hydrologic response in emplacement drifts. This goal is quite different from the goal of a model to predict a single expected thermal-hydrologic response. As a result, the development and validation of the MSTHM and the associated analyses using this model are focused on the goal of predicting a reasonable range of thermal-hydrologic conditions resulting from parametric uncertainty and waste-package-to-waste-package heat-output variability. Thermal-hydrologic conditions within emplacement drifts depend primarily on thermal-hydrologic conditions in the host rock at the drift wall and on the temperature difference between the drift wall and the drip-shield and waste-package surfaces. Thus, the ability to predict a reasonable range of relevant in-drift MSTHM output parameters (e.g., temperature and relative humidity) is based on valid predictions of thermal-hydrologic processes in the host rock, as well as valid predictions of heat-transfer processes between the drift wall and the drip-shield and waste-package surfaces. Because the invert contains crushed gravel derived from the host rock, the invert is, in effect, an extension of the host rock, with thermal and hydrologic properties that have been modified by virtue of the crushing (and the resulting
X. Frank Xu
2010-03-30
Multiscale modeling of stochastic systems, or uncertainty quantization of multiscale modeling is becoming an emerging research frontier, with rapidly growing engineering applications in nanotechnology, biotechnology, advanced materials, and geo-systems, etc. While tremendous efforts have been devoted to either stochastic methods or multiscale methods, little combined work had been done on integration of multiscale and stochastic methods, and there was no method formally available to tackle multiscale problems involving uncertainties. By developing an innovative Multiscale Stochastic Finite Element Method (MSFEM), this research has made a ground-breaking contribution to the emerging field of Multiscale Stochastic Modeling (MSM) (Fig 1). The theory of MSFEM basically decomposes a boundary value problem of random microstructure into a slow scale deterministic problem and a fast scale stochastic one. The slow scale problem corresponds to common engineering modeling practices where fine-scale microstructure is approximated by certain effective constitutive constants, which can be solved by using standard numerical solvers. The fast scale problem evaluates fluctuations of local quantities due to random microstructure, which is important for scale-coupling systems and particularly those involving failure mechanisms. The Green-function-based fast-scale solver developed in this research overcomes the curse-of-dimensionality commonly met in conventional approaches, by proposing a random field-based orthogonal expansion approach. The MSFEM formulated in this project paves the way to deliver the first computational tool/software on uncertainty quantification of multiscale systems. The applications of MSFEM on engineering problems will directly enhance our modeling capability on materials science (composite materials, nanostructures), geophysics (porous media, earthquake), biological systems (biological tissues, bones, protein folding). Continuous development of MSFEM will
Multi-Scale Multi-physics Methods Development for the Calculation...
Office of Scientific and Technical Information (OSTI)
by developing multi-scale, multi-physics methods and implementing them within the ... suitable scales for a physical and mathematical model and then deriving and applying ...
Towards a Multiscale Approach to Cybersecurity Modeling
Hogan, Emilie A.; Hui, Peter SY; Choudhury, Sutanay; Halappanavar, Mahantesh; Oler, Kiri J.; Joslyn, Cliff A.
2013-11-12
We propose a multiscale approach to modeling cyber networks, with the goal of capturing a view of the network and overall situational awareness with respect to a few key properties--- connectivity, distance, and centrality--- for a system under an active attack. We focus on theoretical and algorithmic foundations of multiscale graphs, coming from an algorithmic perspective, with the goal of modeling cyber system defense as a specific use case scenario. We first define a notion of \\emph{multiscale} graphs, in contrast with their well-studied single-scale counterparts. We develop multiscale analogs of paths and distance metrics. As a simple, motivating example of a common metric, we present a multiscale analog of the all-pairs shortest-path problem, along with a multiscale analog of a well-known algorithm which solves it. From a cyber defense perspective, this metric might be used to model the distance from an attacker's position in the network to a sensitive machine. In addition, we investigate probabilistic models of connectivity. These models exploit the hierarchy to quantify the likelihood that sensitive targets might be reachable from compromised nodes. We believe that our novel multiscale approach to modeling cyber-physical systems will advance several aspects of cyber defense, specifically allowing for a more efficient and agile approach to defending these systems.
Chakraborty, Pritam; Zhang, Yongfeng; Tonks, Michael R.
2015-12-07
In this study, the fracture behavior of brittle materials is strongly influenced by their underlying microstructure that needs explicit consideration for accurate prediction of fracture properties and the associated scatter. In this work, a hierarchical multi-scale approach is pursued to model microstructure sensitive brittle fracture. A quantitative phase-field based fracture model is utilized to capture the complex crack growth behavior in the microstructure and the related parameters are calibrated from lower length scale atomistic simulations instead of engineering scale experimental data. The workability of this approach is demonstrated by performing porosity dependent intergranular fracture simulations in UO_{2} and comparing the predictions with experiments.
A multiscale two-point flux-approximation method
Myner, Olav Lie, Knut-Andreas
2014-10-15
A large number of multiscale finite-volume methods have been developed over the past decade to compute conservative approximations to multiphase flow problems in heterogeneous porous media. In particular, several iterative and algebraic multiscale frameworks that seek to reduce the fine-scale residual towards machine precision have been presented. Common for all such methods is that they rely on a compatible primaldual coarse partition, which makes it challenging to extend them to stratigraphic and unstructured grids. Herein, we propose a general idea for how one can formulate multiscale finite-volume methods using only a primal coarse partition. To this end, we use two key ingredients that are computed numerically: (i) elementary functions that correspond to flow solutions used in transmissibility upscaling, and (ii) partition-of-unity functions used to combine elementary functions into basis functions. We exemplify the idea by deriving a multiscale two-point flux-approximation (MsTPFA) method, which is robust with regards to strong heterogeneities in the permeability field and can easily handle general grids with unstructured fine- and coarse-scale connections. The method can easily be adapted to arbitrary levels of coarsening, and can be used both as a standalone solver and as a preconditioner. Several numerical experiments are presented to demonstrate that the MsTPFA method can be used to solve elliptic pressure problems on a wide variety of geological models in a robust and efficient manner.
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Chakraborty, Pritam; Zhang, Yongfeng; Tonks, Michael R.
2015-12-07
In this study, the fracture behavior of brittle materials is strongly influenced by their underlying microstructure that needs explicit consideration for accurate prediction of fracture properties and the associated scatter. In this work, a hierarchical multi-scale approach is pursued to model microstructure sensitive brittle fracture. A quantitative phase-field based fracture model is utilized to capture the complex crack growth behavior in the microstructure and the related parameters are calibrated from lower length scale atomistic simulations instead of engineering scale experimental data. The workability of this approach is demonstrated by performing porosity dependent intergranular fracture simulations in UO2 and comparing themore » predictions with experiments.« less
Final Report for Integrated Multiscale Modeling of Molecular Computing Devices
Glotzer, Sharon C.
2013-08-28
In collaboration with researchers at Vanderbilt University, North Carolina State University, Princeton and Oakridge National Laboratory we developed multiscale modeling and simulation methods capable of modeling the synthesis, assembly, and operation of molecular electronics devices. Our role in this project included the development of coarse-grained molecular and mesoscale models and simulation methods capable of simulating the assembly of millions of organic conducting molecules and other molecular components into nanowires, crossbars, and other organized patterns.
Multiscale Concrete Modeling of Aging Degradation
Hammi, Yousseff; Gullett, Philipp; Horstemeyer, Mark F.
2015-07-31
In this work a numerical finite element framework is implemented to enable the integration of coupled multiscale and multiphysics transport processes. A User Element subroutine (UEL) in Abaqus is used to simultaneously solve stress equilibrium, heat conduction, and multiple diffusion equations for 2D and 3D linear and quadratic elements. Transport processes in concrete structures and their degradation mechanisms are presented along with the discretization of the governing equations. The multiphysics modeling framework is theoretically extended to the linear elastic fracture mechanics (LEFM) by introducing the eXtended Finite Element Method (XFEM) and based on the XFEM user element implementation of Giner et al. [2009]. A damage model that takes into account the damage contribution from the different degradation mechanisms is theoretically developed. The total contribution of damage is forwarded to a Multi-Stage Fatigue (MSF) model to enable the assessment of the fatigue life and the deterioration of reinforced concrete structures in a nuclear power plant. Finally, two examples are presented to illustrate the developed multiphysics user element implementation and the XFEM implementation of Giner et al. [2009].
Moist multi-scale models for the hurricane embryo
Majda, Andrew J.; Xing, Yulong; Mohammadian, Majid
2010-01-01
Determining the finite-amplitude preconditioned states in the hurricane embryo, which lead to tropical cyclogenesis, is a central issue in contemporary meteorology. In the embryo there is competition between different preconditioning mechanisms involving hydrodynamics and moist thermodynamics, which can lead to cyclogenesis. Here systematic asymptotic methods from applied mathematics are utilized to develop new simplified moist multi-scale models starting from the moist anelastic equations. Three interesting multi-scale models emerge in the analysis. The balanced mesoscale vortex (BMV) dynamics and the microscale balanced hot tower (BHT) dynamics involve simplified balanced equations without gravity waves for vertical vorticity amplification due to moist heat sources and incorporate nonlinear advective fluxes across scales. The BMV model is the central one for tropical cyclogenesis in the embryo. The moist mesoscale wave (MMW) dynamics involves simplified equations for mesoscale moisture fluctuations, as well as linear hydrostatic waves driven by heat sources from moisture and eddy flux divergences. A simplified cloud physics model for deep convection is introduced here and used to study moist axisymmetric plumes in the BHT model. A simple application in periodic geometry involving the effects of mesoscale vertical shear and moist microscale hot towers on vortex amplification is developed here to illustrate features of the coupled multi-scale models. These results illustrate the use of these models in isolating key mechanisms in the embryo in a simplified content.
A multilevel multiscale mimetic method for an anisotropic infiltration problem
Lipnikov, Konstantin; Moulton, David; Svyatskiy, Daniil
2009-01-01
Modeling of multiphase flow and transport in highly heterogeneous porous media must capture a broad range of coupled spatial and temporal scales. Recently, a hierarchical approach dubbed the Multilevel Multiscale Mimetic (M3) method, was developed to simulate two-phase flow in porous media. The M{sup 3} method is locally mass conserving at all levels in its hierarchy, it supports unstructured polygonal grids and full tensor permeabilities, and it can achieve large coarsening factors. In this work we consider infiltration of water into a two-dimensional layered medium. The grid is aligned with the layers but not the coordinate axes. We demonstrate that with an efficient temporal updating strategy for the coarsening parameters, fine-scale accuracy of prominent features in the flow is maintained by the M{sup 3} method.
Towards a Multiscale Approach to Cybersecurity Modeling (Conference...
Office of Scientific and Technical Information (OSTI)
We propose a multiscale approach to modeling cyber networks, with the goal of capturing a view of the network and overall situational awareness with respect to a few key ...
Multiscale Modeling of Energy Storage Materials | Argonne Leadership...
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Lithium ions are shown as spheres. Gregory A. Voth Multiscale Modeling of Energy Storage ... listed the transformation of the energy system of the country as a central goal to ...
Multiscale Modeling of Energy Storage Materials | Argonne Leadership
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Computing Facility Multiscale Modeling of Energy Storage Materials PI Name: Gregory A. Voth PI Email: gavoth@uchicago.edu Institution: University of Chicago and Argonne National Laboratory Allocation Program: INCITE Allocation Hours at ALCF: 25,000,000 Year: 2012 Research Domain: Materials Science The leadership-class computing resources provided by the INCITE program will be used for the multiscale modeling of charge transport processes in materials relevant to fuel cell and battery
A bidirectional coupling procedure applied to multiscale respiratory modeling
Kuprat, A.P.; Kabilan, S.; Carson, J.P.; Corley, R.A.; Einstein, D.R.
2013-07-01
In this study, we present a novel multiscale computational framework for efficiently linking multiple lower-dimensional models describing the distal lung mechanics to imaging-based 3D computational fluid dynamics (CFDs) models of the upper pulmonary airways in order to incorporate physiologically appropriate outlet boundary conditions. The framework is an extension of the modified Newton’s method with nonlinear Krylov accelerator developed by Carlson and Miller [1], Miller [2] and Scott and Fenves [3]. Our extensions include the retention of subspace information over multiple timesteps, and a special correction at the end of a timestep that allows for corrections to be accepted with verified low residual with as little as a single residual evaluation per timestep on average. In the case of a single residual evaluation per timestep, the method has zero additional computational cost compared to uncoupled or unidirectionally coupled simulations. We expect these enhancements to be generally applicable to other multiscale coupling applications where timestepping occurs. In addition we have developed a “pressure-drop” residual which allows for stable coupling of flows between a 3D incompressible CFD application and another (lower-dimensional) fluid system. We expect this residual to also be useful for coupling non-respiratory incompressible fluid applications, such as multiscale simulations involving blood flow. The lower-dimensional models that are considered in this study are sets of simple ordinary differential equations (ODEs) representing the compliant mechanics of symmetric human pulmonary airway trees. To validate the method, we compare the predictions of hybrid CFD–ODE models against an ODE-only model of pulmonary airflow in an idealized geometry. Subsequently, we couple multiple sets of ODEs describing the distal lung to an imaging-based human lung geometry. Boundary conditions in these models consist of atmospheric pressure at the mouth and intrapleural
A Bidirectional Coupling Procedure Applied to Multiscale Respiratory Modeling
Kuprat, Andrew P.; Kabilan, Senthil; Carson, James P.; Corley, Richard A.; Einstein, Daniel R.
2013-07-01
In this study, we present a novel multiscale computational framework for efficiently linking multiple lower-dimensional models describing the distal lung mechanics to imaging-based 3D computational fluid dynamics (CFD) models of the upper pulmonary airways in order to incorporate physiologically appropriate outlet boundary conditions. The framework is an extension of the Modified Newton’s Method with nonlinear Krylov accelerator developed by Carlson and Miller [1, 2, 3]. Our extensions include the retention of subspace information over multiple timesteps, and a special correction at the end of a timestep that allows for corrections to be accepted with verified low residual with as little as a single residual evaluation per timestep on average. In the case of a single residual evaluation per timestep, the method has zero additional computational cost compared to uncoupled or unidirectionally coupled simulations. We expect these enhancements to be generally applicable to other multiscale coupling applications where timestepping occurs. In addition we have developed a “pressure-drop” residual which allows for stable coupling of flows between a 3D incompressible CFD application and another (lower-dimensional) fluid system. We expect this residual to also be useful for coupling non-respiratory incompressible fluid applications, such as multiscale simulations involving blood flow. The lower-dimensional models that are considered in this study are sets of simple ordinary differential equations (ODEs) representing the compliant mechanics of symmetric human pulmonary airway trees. To validate the method, we compare the predictions of hybrid CFD-ODE models against an ODE-only model of pulmonary airflow in an idealized geometry. Subsequently, we couple multiple sets of ODEs describing the distal lung to an imaging-based human lung geometry. Boundary conditions in these models consist of atmospheric pressure at the mouth and intrapleural pressure applied to the multiple
Gao, Kai; Fu, Shubin; Gibson, Richard L.; Chung, Eric T.; Efendiev, Yalchin
2015-08-15
It is important to develop fast yet accurate numerical methods for seismic wave propagation to characterize complex geological structures and oil and gas reservoirs. However, the computational cost of conventional numerical modeling methods, such as finite-difference method and finite-element method, becomes prohibitively expensive when applied to very large models. We propose a Generalized Multiscale Finite-Element Method (GMsFEM) for elastic wave propagation in heterogeneous, anisotropic media, where we construct basis functions from multiple local problems for both the boundaries and interior of a coarse node support or coarse element. The application of multiscale basis functions can capture the fine scale medium property variations, and allows us to greatly reduce the degrees of freedom that are required to implement the modeling compared with conventional finite-element method for wave equation, while restricting the error to low values. We formulate the continuous Galerkin and discontinuous Galerkin formulation of the multiscale method, both of which have pros and cons. Applications of the multiscale method to three heterogeneous models show that our multiscale method can effectively model the elastic wave propagation in anisotropic media with a significant reduction in the degrees of freedom in the modeling system.
Gao, Kai; Fu, Shubin; Gibson, Richard L.; Chung, Eric T.; Efendiev, Yalchin
2015-04-14
It is important to develop fast yet accurate numerical methods for seismic wave propagation to characterize complex geological structures and oil and gas reservoirs. However, the computational cost of conventional numerical modeling methods, such as finite-difference method and finite-element method, becomes prohibitively expensive when applied to very large models. We propose a Generalized Multiscale Finite-Element Method (GMsFEM) for elastic wave propagation in heterogeneous, anisotropic media, where we construct basis functions from multiple local problems for both the boundaries and interior of a coarse node support or coarse element. The application of multiscale basis functions can capture the fine scale medium property variations, and allows us to greatly reduce the degrees of freedom that are required to implement the modeling compared with conventional finite-element method for wave equation, while restricting the error to low values. We formulate the continuous Galerkin and discontinuous Galerkin formulation of the multiscale method, both of which have pros and cons. Applications of the multiscale method to three heterogeneous models show that our multiscale method can effectively model the elastic wave propagation in anisotropic media with a significant reduction in the degrees of freedom in the modeling system.
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Gao, Kai; Fu, Shubin; Gibson, Richard L.; Chung, Eric T.; Efendiev, Yalchin
2015-04-14
It is important to develop fast yet accurate numerical methods for seismic wave propagation to characterize complex geological structures and oil and gas reservoirs. However, the computational cost of conventional numerical modeling methods, such as finite-difference method and finite-element method, becomes prohibitively expensive when applied to very large models. We propose a Generalized Multiscale Finite-Element Method (GMsFEM) for elastic wave propagation in heterogeneous, anisotropic media, where we construct basis functions from multiple local problems for both the boundaries and interior of a coarse node support or coarse element. The application of multiscale basis functions can capture the fine scale mediummore » property variations, and allows us to greatly reduce the degrees of freedom that are required to implement the modeling compared with conventional finite-element method for wave equation, while restricting the error to low values. We formulate the continuous Galerkin and discontinuous Galerkin formulation of the multiscale method, both of which have pros and cons. Applications of the multiscale method to three heterogeneous models show that our multiscale method can effectively model the elastic wave propagation in anisotropic media with a significant reduction in the degrees of freedom in the modeling system.« less
Smith, Kandler; Graf, Peter; Jun, Myungsoo; Yang, Chuanbo; Li, Genong; Li, Shaoping; Hochman, Amit; Tselepidakis, Dimitrios
2015-06-09
This presentation provides an update on improvements in computational efficiency in a nonlinear multiscale battery model for computer aided engineering.
Fast multiscale Gaussian beam methods for wave equations in bounded convex domains
Bao, Gang; Department of Mathematics, Michigan State University, East Lansing, MI 48824 ; Lai, Jun; Qian, Jianliang
2014-03-15
Motivated by fast multiscale Gaussian wavepacket transforms and multiscale Gaussian beam methods which were originally designed for pure initial-value problems of wave equations, we develop fast multiscale Gaussian beam methods for initial boundary value problems of wave equations in bounded convex domains in the high frequency regime. To compute the wave propagation in bounded convex domains, we have to take into account reflecting multiscale Gaussian beams, which are accomplished by enforcing reflecting boundary conditions during beam propagation and carrying out suitable reflecting beam summation. To propagate multiscale beams efficiently, we prove that the ratio of the squared magnitude of beam amplitude and the beam width is roughly conserved, and accordingly we propose an effective indicator to identify significant beams. We also prove that the resulting multiscale Gaussian beam methods converge asymptotically. Numerical examples demonstrate the accuracy and efficiency of the method.
Mathematical and Numerical Analyses of Peridynamics for Multiscale Materials Modeling
Du, Qiang
2014-11-12
The rational design of materials, the development of accurate and efficient material simulation algorithms, and the determination of the response of materials to environments and loads occurring in practice all require an understanding of mechanics at disparate spatial and temporal scales. The project addresses mathematical and numerical analyses for material problems for which relevant scales range from those usually treated by molecular dynamics all the way up to those most often treated by classical elasticity. The prevalent approach towards developing a multiscale material model couples two or more well known models, e.g., molecular dynamics and classical elasticity, each of which is useful at a different scale, creating a multiscale multi-model. However, the challenges behind such a coupling are formidable and largely arise because the atomistic and continuum models employ nonlocal and local models of force, respectively. The project focuses on a multiscale analysis of the peridynamics materials model. Peridynamics can be used as a transition between molecular dynamics and classical elasticity so that the difficulties encountered when directly coupling those two models are mitigated. In addition, in some situations, peridynamics can be used all by itself as a material model that accurately and efficiently captures the behavior of materials over a wide range of spatial and temporal scales. Peridynamics is well suited to these purposes because it employs a nonlocal model of force, analogous to that of molecular dynamics; furthermore, at sufficiently large length scales and assuming smooth deformation, peridynamics can be approximated by classical elasticity. The project will extend the emerging mathematical and numerical analysis of peridynamics. One goal is to develop a peridynamics-enabled multiscale multi-model that potentially provides a new and more extensive mathematical basis for coupling classical elasticity and molecular dynamics, thus enabling next
Multi-scale Modeling of Plasticity in Tantalum.
Lim, Hojun; Battaile, Corbett Chandler.; Carroll, Jay; Buchheit, Thomas E.; Boyce, Brad; Weinberger, Christopher
2015-12-01
In this report, we present a multi-scale computational model to simulate plastic deformation of tantalum and validating experiments. In atomistic/ dislocation level, dislocation kink- pair theory is used to formulate temperature and strain rate dependent constitutive equations. The kink-pair theory is calibrated to available data from single crystal experiments to produce accurate and convenient constitutive laws. The model is then implemented into a BCC crystal plasticity finite element method (CP-FEM) model to predict temperature and strain rate dependent yield stresses of single and polycrystalline tantalum and compared with existing experimental data from the literature. Furthermore, classical continuum constitutive models describing temperature and strain rate dependent flow behaviors are fit to the yield stresses obtained from the CP-FEM polycrystal predictions. The model is then used to conduct hydro- dynamic simulations of Taylor cylinder impact test and compared with experiments. In order to validate the proposed tantalum CP-FEM model with experiments, we introduce a method for quantitative comparison of CP-FEM models with various experimental techniques. To mitigate the effects of unknown subsurface microstructure, tantalum tensile specimens with a pseudo-two-dimensional grain structure and grain sizes on the order of millimeters are used. A technique combining an electron back scatter diffraction (EBSD) and high resolution digital image correlation (HR-DIC) is used to measure the texture and sub-grain strain fields upon uniaxial tensile loading at various applied strains. Deformed specimens are also analyzed with optical profilometry measurements to obtain out-of- plane strain fields. These high resolution measurements are directly compared with large-scale CP-FEM predictions. This computational method directly links fundamental dislocation physics to plastic deformations in the grain-scale and to the engineering-scale applications. Furthermore, direct
Multiscale Modeling of Hemodynamic Disorders" | Argonne Leadership
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Computing Facility Multiscale Modeling of Hemodynamic Disorders" Authors: Fedosov, D., Pivkin, I., Pan, W., Dao, M., Caswell, B., Karniadakis, G.E. This book offers a mathematical update of the state of the art of the research in the field of mathematical and numerical models of the circulatory system. It is structured into different chapters, written by outstanding experts in the field. Many fundamental issues are considered, such as: the mathematical representation of vascular
Multiscale Modeling of Malaria | Argonne Leadership Computing Facility
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Malaria Authors: Karniadakis, G.E., Parasitic infectious diseases like malaria and certain hereditary hematologic disorders are often associated with major changes in the shape and viscoelastic properties of red blood cells. Such changes can disrupt blood flow and, possibly, brain perfusion, as in the case of cerebral malaria. In recent work on stochastic multiscale models-in conjunction with large-scale parallel computing-we were able to quantify, for the first time, the main biophysical
Multiscale Modeling of Sickle Anemia Blood Blow by Dissipative Partice
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Dynamics | Argonne Leadership Computing Facility Sickle Anemia Blood Blow by Dissipative Partice Dynamics Authors: Huan, L., Caswell, B., Karniadakis, G A multi-scale model for sickle red blood cell is developed based on Dissipative Particle Dynamics (DPD). Different cell morphologies (sickle, granular, elongated shapes) typically observed in in vitro and in vivo are constructed and the deviations from the biconcave shape is quantified by the Asphericity and Elliptical shape factors. The
Kim, G.; Pesaran, A.; Smith, K.; Graf, P.; Jun, M.; Yang, C.; Li, G.; Li, S.; Hochman, A.; Tselepidakis, D.; White, J.
2014-06-01
This presentation discusses the significant enhancement of computational efficiency in nonlinear multiscale battery model for computer aided engineering in current research at NREL.
Multi-Scale Multi-Dimensional Model for Better Cell Design and Management (Presentation)
Kim, G.-H.; Smith, K.
2008-09-01
Describes NREL's R&D to develop a multi-scale model to assist in designing better, more reliable lithium-ion battery cells for advanced vehicles.
Multiscale modeling for fluid transport in nanosystems.
Lee, Jonathan W.; Jones, Reese E.; Mandadapu, Kranthi Kiran; Templeton, Jeremy Alan; Zimmerman, Jonathan A.
2013-09-01
Atomistic-scale behavior drives performance in many micro- and nano-fluidic systems, such as mircrofludic mixers and electrical energy storage devices. Bringing this information into the traditionally continuum models used for engineering analysis has proved challenging. This work describes one such approach to address this issue by developing atomistic-to-continuum multi scale and multi physics methods to enable molecular dynamics (MD) representations of atoms to incorporated into continuum simulations. Coupling is achieved by imposing constraints based on fluxes of conserved quantities between the two regions described by one of these models. The impact of electric fields and surface charges are also critical, hence, methodologies to extend finite-element (FE) MD electric field solvers have been derived to account for these effects. Finally, the continuum description can have inconsistencies with the coarse-grained MD dynamics, so FE equations based on MD statistics were derived to facilitate the multi scale coupling. Examples are shown relevant to nanofluidic systems, such as pore flow, Couette flow, and electric double layer.
Multi-Scale Multi-Dimensional Ion Battery Performance Model
Energy Science and Technology Software Center (OSTI)
2007-05-07
The Multi-Scale Multi-Dimensional (MSMD) Lithium Ion Battery Model allows for computer prediction and engineering optimization of thermal, electrical, and electrochemical performance of lithium ion cells with realistic geometries. The model introduces separate simulation domains for different scale physics, achieving much higher computational efficiency compared to the single domain approach. It solves a one dimensional electrochemistry model in a micro sub-grid system, and captures the impacts of macro-scale battery design factors on cell performance and materialmore » usage by solving cell-level electron and heat transports in a macro grid system.« less
Fluid simulations with atomistic resolution: a hybrid multiscale method with field-wise coupling
Borg, Matthew K. [Department of Mechanical and Aerospace Engineering, University of Strathclyde, Glasgow G1 1XJ (United Kingdom)] [Department of Mechanical and Aerospace Engineering, University of Strathclyde, Glasgow G1 1XJ (United Kingdom); Lockerby, Duncan A., E-mail: duncan.lockerby@warwick.ac.uk [School of Engineering, University of Warwick, Coventry CV4 7AL (United Kingdom); Reese, Jason M., E-mail: jason.reese@strath.ac.uk [Department of Mechanical and Aerospace Engineering, University of Strathclyde, Glasgow G1 1XJ (United Kingdom)
2013-12-15
We present a new hybrid method for simulating dense fluid systems that exhibit multiscale behaviour, in particular, systems in which a NavierStokes model may not be valid in parts of the computational domain. We apply molecular dynamics as a local microscopic refinement for correcting the NavierStokes constitutive approximation in the bulk of the domain, as well as providing a direct measurement of velocity slip at bounding surfaces. Our hybrid approach differs from existing techniques, such as the heterogeneous multiscale method (HMM), in some fundamental respects. In our method, the individual molecular solvers, which provide information to the macro model, are not coupled with the continuum grid at nodes (i.e. point-wise coupling), instead coupling occurs over distributed heterogeneous fields (here referred to as field-wise coupling). This affords two major advantages. Whereas point-wise coupled HMM is limited to regions of flow that are highly scale-separated in all spatial directions (i.e. where the state of non-equilibrium in the fluid can be adequately described by a single strain tensor and temperature gradient vector), our field-wise coupled HMM has no such limitations and so can be applied to flows with arbitrarily-varying degrees of scale separation (e.g. flow from a large reservoir into a nano-channel). The second major advantage is that the position of molecular elements does not need to be collocated with nodes of the continuum grid, which means that the resolution of the microscopic correction can be adjusted independently of the resolution of the continuum model. This in turn means the computational cost and accuracy of the molecular correction can be independently controlled and optimised. The macroscopic constraints on the individual molecular solvers are artificial body-force distributions, used in conjunction with standard periodicity. We test our hybrid method on the Poiseuille flow problem for both Newtonian (Lennard-Jones) and non
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Gao, Kai; Chung, Eric T.; Gibson, Richard L.; Fu, Shubin; Efendiev, Yalchin
2015-06-05
The development of reliable methods for upscaling fine scale models of elastic media has long been an important topic for rock physics and applied seismology. Several effective medium theories have been developed to provide elastic parameters for materials such as finely layered media or randomly oriented or aligned fractures. In such cases, the analytic solutions for upscaled properties can be used for accurate prediction of wave propagation. However, such theories cannot be applied directly to homogenize elastic media with more complex, arbitrary spatial heterogeneity. We therefore propose a numerical homogenization algorithm based on multiscale finite element methods for simulating elasticmore » wave propagation in heterogeneous, anisotropic elastic media. Specifically, our method used multiscale basis functions obtained from a local linear elasticity problem with appropriately defined boundary conditions. Homogenized, effective medium parameters were then computed using these basis functions, and the approach applied a numerical discretization that is similar to the rotated staggered-grid finite difference scheme. Comparisons of the results from our method and from conventional, analytical approaches for finely layered media showed that the homogenization reliably estimated elastic parameters for this simple geometry. Additional tests examined anisotropic models with arbitrary spatial heterogeneity where the average size of the heterogeneities ranged from several centimeters to several meters, and the ratio between the dominant wavelength and the average size of the arbitrary heterogeneities ranged from 10 to 100. Comparisons to finite-difference simulations proved that the numerical homogenization was equally accurate for these complex cases.« less
Gao, Kai; Chung, Eric T.; Gibson, Richard L.; Fu, Shubin; Efendiev, Yalchin
2015-06-05
The development of reliable methods for upscaling fine scale models of elastic media has long been an important topic for rock physics and applied seismology. Several effective medium theories have been developed to provide elastic parameters for materials such as finely layered media or randomly oriented or aligned fractures. In such cases, the analytic solutions for upscaled properties can be used for accurate prediction of wave propagation. However, such theories cannot be applied directly to homogenize elastic media with more complex, arbitrary spatial heterogeneity. We therefore propose a numerical homogenization algorithm based on multiscale finite element methods for simulating elastic wave propagation in heterogeneous, anisotropic elastic media. Specifically, our method used multiscale basis functions obtained from a local linear elasticity problem with appropriately defined boundary conditions. Homogenized, effective medium parameters were then computed using these basis functions, and the approach applied a numerical discretization that is similar to the rotated staggered-grid finite difference scheme. Comparisons of the results from our method and from conventional, analytical approaches for finely layered media showed that the homogenization reliably estimated elastic parameters for this simple geometry. Additional tests examined anisotropic models with arbitrary spatial heterogeneity where the average size of the heterogeneities ranged from several centimeters to several meters, and the ratio between the dominant wavelength and the average size of the arbitrary heterogeneities ranged from 10 to 100. Comparisons to finite-difference simulations proved that the numerical homogenization was equally accurate for these complex cases.
Integrated Multiscale Modeling of Molecular Computing Devices
Weinan E
2012-03-29
The main bottleneck in modeling transport in molecular devices is to develop the correct formulation of the problem and efficient algorithms for analyzing the electronic structure and dynamics using, for example, the time-dependent density functional theory. We have divided this task into several steps. The first step is to developing the right mathematical formulation and numerical algorithms for analyzing the electronic structure using density functional theory. The second step is to study time-dependent density functional theory, particularly the far-field boundary conditions. The third step is to study electronic transport in molecular devices. We are now at the end of the first step. Under DOE support, we have made subtantial progress in developing linear scaling and sub-linear scaling algorithms for electronic structure analysis. Although there has been a huge amount of effort in the past on developing linear scaling algorithms, most of the algorithms developed suffer from the lack of robustness and controllable accuracy. We have made the following progress: (1) We have analyzed thoroughly the localization properties of the wave-functions. We have developed a clear understanding of the physical as well as mathematical origin of the decay properties. One important conclusion is that even for metals, one can choose wavefunctions that decay faster than any algebraic power. (2) We have developed algorithms that make use of these localization properties. Our algorithms are based on non-orthogonal formulations of the density functional theory. Our key contribution is to add a localization step into the algorithm. The addition of this localization step makes the algorithm quite robust and much more accurate. Moreover, we can control the accuracy of these algorithms by changing the numerical parameters. (3) We have considerably improved the Fermi operator expansion (FOE) approach. Through pole expansion, we have developed the optimal scaling FOE algorithm.
Multiscale Design of Advanced Materials based on Hybrid Ab Initio and Quasicontinuum Methods
Luskin, Mitchell
2014-03-12
This project united researchers from mathematics, chemistry, computer science, and engineering for the development of new multiscale methods for the design of materials. Our approach was highly interdisciplinary, but it had two unifying themes: first, we utilized modern mathematical ideas about change-of-scale and state-of-the-art numerical analysis to develop computational methods and codes to solve real multiscale problems of DOE interest; and, second, we took very seriously the need for quantum mechanics-based atomistic forces, and based our methods on fast solvers of chemically accurate methods.
A MULTISCALE, CELL-BASED FRAMEWORK FOR MODELING CANCER DEVELOPMENT
JIANG, YI
2007-01-16
Cancer remains to be one of the leading causes of death due to diseases. We use a systems approach that combines mathematical modeling, numerical simulation, in vivo and in vitro experiments, to develop a predictive model that medical researchers can use to study and treat cancerous tumors. The multiscale, cell-based model includes intracellular regulations, cellular level dynamics and intercellular interactions, and extracellular level chemical dynamics. The intracellular level protein regulations and signaling pathways are described by Boolean networks. The cellular level growth and division dynamics, cellular adhesion and interaction with the extracellular matrix is described by a lattice Monte Carlo model (the Cellular Potts Model). The extracellular dynamics of the signaling molecules and metabolites are described by a system of reaction-diffusion equations. All three levels of the model are integrated through a hybrid parallel scheme into a high-performance simulation tool. The simulation results reproduce experimental data in both avasular tumors and tumor angiogenesis. By combining the model with experimental data to construct biologically accurate simulations of tumors and their vascular systems, this model will enable medical researchers to gain a deeper understanding of the cellular and molecular interactions associated with cancer progression and treatment.
A Mathematical Analysis of Atomistic-to-Continuum (AtC) Multiscale Coupling Methods
Gunzburger, Max
2013-11-13
We have worked on several projects aimed at improving the efficiency and understanding of multiscale methods, especially those applicable to problems involving atomistic-to-continuum coupling. Activities include blending methods for AtC coupling and efficient quasi-continuum methods for problems with long-range interactions.
Watt-Sun: A Multi-Scale, Multi-Model, Machine-Learning Solar Forecasting
Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site
Technology | Department of Energy Watt-Sun: A Multi-Scale, Multi-Model, Machine-Learning Solar Forecasting Technology Watt-Sun: A Multi-Scale, Multi-Model, Machine-Learning Solar Forecasting Technology IBM logo.png As part of this project, new solar forecasting technology will be developed that leverages big data processing, deep machine learning, and cloud modeling integrated in a universal platform with an open architecture. Similar to the Watson computer system, this proposed technology
SU-F-18C-15: Model-Based Multiscale Noise Reduction On Low Dose Cone Beam Projection
Yao, W; Farr, J
2014-06-15
Purpose: To improve image quality of low dose cone beam CT for patient positioning in radiation therapy. Methods: In low dose cone beam CT (CBCT) imaging systems, Poisson process governs the randomness of photon fluence at x-ray source and the detector because of the independent binomial process of photon absorption in medium. On a CBCT projection, the variance of fluence consists of the variance of noiseless imaging structure and that of Poisson noise, which is proportional to the mean (noiseless) of the fluence at the detector. This requires multiscale filters to smoothen noise while keeping the structure information of the imaged object. We used a mathematical model of Poisson process to design multiscale filters and established the balance of noise correction and structure blurring. The algorithm was checked with low dose kilo-voltage CBCT projections acquired from a Varian OBI system. Results: From the investigation of low dose CBCT of a Catphan phantom and patients, it showed that our model-based multiscale technique could efficiently reduce noise and meanwhile keep the fine structure of the imaged object. After the image processing, the number of visible line pairs in Catphan phantom scanned with 4 ms pulse time was similar to that scanned with 32 ms, and soft tissue structure from simulated 4 ms patient head-and-neck images was also comparable with scanned 20 ms ones. Compared with fixed-scale technique, the image quality from multiscale one was improved. Conclusion: Use of projection-specific multiscale filters can reach better balance on noise reduction and structure information loss. The image quality of low dose CBCT can be improved by using multiscale filters.
Multiscale Modeling of the Orthotropic Behaviour of PA6-6 overmoulded Composites using MMI Approach
Bikard, Jerome; Robert, Gilles; Moulinjeune, Olivier [RHODIA ENGINEERING PLASTICS, Technyl Application Center Avenue Ramboz, BP 64, 69192 Saint FONS CEDEX (France)
2011-05-04
In this study the MMI ConfidentDesign multiscale approach (consisting in a non-linear multiscale simulation based on DIGIMAT registered including the injection modeling of the filled polymer and a multiscale mechanical model using the fiber orientation tensor resulting from the injection) has been combined with an orthotropic damageable elastic simulation. The anisotropic properties (including rupture criterion) are estimated and a multiscale simulation including the heterogeneous material properties issued from injection process is done. The impact of fiber ratios is then investigated. The structural simulation predicts stresses localized close to the punch, as well in injected PA66 than in composite part. Greater the fiber volume ratio, greater the modulus and more brittle the composite.
Components for Atomistic-to-Continuum Multiscale Modeling of Flow in Micro- and Nanofluidic Systems
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Adalsteinsson, Helgi; Debusschere, Bert J.; Long, Kevin R.; Najm, Habib N.
2008-01-01
Micro- and nanofluidics pose a series of significant challenges for science-based modeling. Key among those are the wide separation of length- and timescales between interface phenomena and bulk flow and the spatially heterogeneous solution properties near solid-liquid interfaces. It is not uncommon for characteristic scales in these systems to span nine orders of magnitude from the atomic motions in particle dynamics up to evolution of mass transport at the macroscale level, making explicit particle models intractable for all but the simplest systems. Recently, atomistic-to-continuum (A2C) multiscale simulations have gained a lot of interest as an approach to rigorously handle particle-levelmore » dynamics while also tracking evolution of large-scale macroscale behavior. While these methods are clearly not applicable to all classes of simulations, they are finding traction in systems in which tight-binding, and physically important, dynamics at system interfaces have complex effects on the slower-evolving large-scale evolution of the surrounding medium. These conditions allow decomposition of the simulation into discrete domains, either spatially or temporally. In this paper, we describe how features of domain decomposed simulation systems can be harnessed to yield flexible and efficient software for multiscale simulations of electric field-driven micro- and nanofluidics.« less
Bayesian data assimilation for stochastic multiscale models of transport in porous media.
Marzouk, Youssef M.; van Bloemen Waanders, Bart Gustaaf; Parno, Matthew; Ray, Jaideep; Lefantzi, Sophia; Salazar, Luke; McKenna, Sean Andrew; Klise, Katherine A.
2011-10-01
We investigate Bayesian techniques that can be used to reconstruct field variables from partial observations. In particular, we target fields that exhibit spatial structures with a large spectrum of lengthscales. Contemporary methods typically describe the field on a grid and estimate structures which can be resolved by it. In contrast, we address the reconstruction of grid-resolved structures as well as estimation of statistical summaries of subgrid structures, which are smaller than the grid resolution. We perform this in two different ways (a) via a physical (phenomenological), parameterized subgrid model that summarizes the impact of the unresolved scales at the coarse level and (b) via multiscale finite elements, where specially designed prolongation and restriction operators establish the interscale link between the same problem defined on a coarse and fine mesh. The estimation problem is posed as a Bayesian inverse problem. Dimensionality reduction is performed by projecting the field to be inferred on a suitable orthogonal basis set, viz. the Karhunen-Loeve expansion of a multiGaussian. We first demonstrate our techniques on the reconstruction of a binary medium consisting of a matrix with embedded inclusions, which are too small to be grid-resolved. The reconstruction is performed using an adaptive Markov chain Monte Carlo method. We find that the posterior distributions of the inferred parameters are approximately Gaussian. We exploit this finding to reconstruct a permeability field with long, but narrow embedded fractures (which are too fine to be grid-resolved) using scalable ensemble Kalman filters; this also allows us to address larger grids. Ensemble Kalman filtering is then used to estimate the values of hydraulic conductivity and specific yield in a model of the High Plains Aquifer in Kansas. Strong conditioning of the spatial structure of the parameters and the non-linear aspects of the water table aquifer create difficulty for the ensemble Kalman
Sondak, David; Oberai, Assad A.
2012-10-15
Novel large eddy simulation (LES) models are developed for incompressible magnetohydrodynamics (MHD). These models include the application of the variational multiscale formulation of LES to the equations of incompressible MHD. Additionally, a new residual-based eddy viscosity model is introduced for MHD. A mixed LES model that combines the strengths of both of these models is also derived. The new models result in a consistent numerical method that is relatively simple to implement. The need for a dynamic procedure in determining model coefficients is no longer required. The new LES models are tested on a decaying Taylor-Green vortex generalized to MHD and benchmarked against classical LES turbulence models. The LES simulations are run in a periodic box of size [-{pi}, {pi}]{sup 3} with 32 modes in each direction and are compared to a direct numerical simulation (DNS) with 512 modes in each direction. The new models are able to account for the essential MHD physics which is demonstrated via comparisons of energy spectra. We also compare the performance of our models to a DNS simulation by Pouquet et al.['The dynamics of unforced turbulence at high Reynolds number for Taylor-Green vortices generalized to MHD,' Geophys. Astrophys. Fluid Dyn. 104, 115-134 (2010)], for which the ratio of DNS modes to LES modes is 262:144.
A multiscale method for the analysis of defect behavior in MO during electron irradiation
Rest, J.; Insepov, Z.; Ye, B.; Yun, D.
2014-10-01
In order to overcome a lack of experimental information on values for key materials properties and kinetic coefficients, a multiscale modeling approach is applied to defect behavior in irradiated Mo where key materials properties, such as point defect (vacancy and interstitial) migration enthalpies as well as kinetic factors such as dimer formation, defect recombination, and self interstitial–interstitial loop interaction coefficients, are obtained by molecular dynamics calculations and implemented into rate-theory simulations of defect behavior. The multiscale methodology is validated against interstitial loop growth data obtained from electron irradiation of pure Mo. It is shown that the observed linear behavior of the loop diameter vs. the square root of irradiation time is a direct consequence of the 1D migration of self-interstitial atoms.
Mathematical and Numerical Analyses of Peridynamics for Multiscale Materials Modeling
Gunzburger, Max
2015-02-17
We have treated the modeling, analysis, numerical analysis, and algorithmic development for nonlocal models of diffusion and mechanics. Variational formulations were developed and finite element methods were developed based on those formulations for both steady state and time dependent problems. Obstacle problems and optimization problems for the nonlocal models were also treated and connections made with fractional derivative models.
A second gradient theoretical framework for hierarchical multiscale modeling of materials
Luscher, Darby J; Bronkhorst, Curt A; Mc Dowell, David L
2009-01-01
A theoretical framework for the hierarchical multiscale modeling of inelastic response of heterogeneous materials has been presented. Within this multiscale framework, the second gradient is used as a non local kinematic link between the response of a material point at the coarse scale and the response of a neighborhood of material points at the fine scale. Kinematic consistency between these scales results in specific requirements for constraints on the fluctuation field. The wryness tensor serves as a second-order measure of strain. The nature of the second-order strain induces anti-symmetry in the first order stress at the coarse scale. The multiscale ISV constitutive theory is couched in the coarse scale intermediate configuration, from which an important new concept in scale transitions emerges, namely scale invariance of dissipation. Finally, a strategy for developing meaningful kinematic ISVs and the proper free energy functions and evolution kinetics is presented.
A voxel-based multiscale model to simulate the radiation response of hypoxic tumors
Espinoza, I.; Peschke, P.; Karger, C. P.
2015-01-15
Purpose: In radiotherapy, it is important to predict the response of tumors to irradiation prior to the treatment. This is especially important for hypoxic tumors, which are known to be highly radioresistant. Mathematical modeling based on the dose distribution, biological parameters, and medical images may help to improve this prediction and to optimize the treatment plan. Methods: A voxel-based multiscale tumor response model for simulating the radiation response of hypoxic tumors was developed. It considers viable and dead tumor cells, capillary and normal cells, as well as the most relevant biological processes such as (i) proliferation of tumor cells, (ii) hypoxia-induced angiogenesis, (iii) spatial exchange of cells leading to tumor growth, (iv) oxygen-dependent cell survival after irradiation, (v) resorption of dead cells, and (vi) spatial exchange of cells leading to tumor shrinkage. Oxygenation is described on a microscopic scale using a previously published tumor oxygenation model, which calculates the oxygen distribution for each voxel using the vascular fraction as the most important input parameter. To demonstrate the capabilities of the model, the dependence of the oxygen distribution on tumor growth and radiation-induced shrinkage is investigated. In addition, the impact of three different reoxygenation processes is compared and tumor control probability (TCP) curves for a squamous cells carcinoma of the head and neck (HNSSC) are simulated under normoxic and hypoxic conditions. Results: The model describes the spatiotemporal behavior of the tumor on three different scales: (i) on the macroscopic scale, it describes tumor growth and shrinkage during radiation treatment, (ii) on a mesoscopic scale, it provides the cell density and vascular fraction for each voxel, and (iii) on the microscopic scale, the oxygen distribution may be obtained in terms of oxygen histograms. With increasing tumor size, the simulated tumors develop a hypoxic core. Within the
Collaborating for Multi-Scale Chemical Science
William H. Green
2006-07-14
Advanced model reduction methods were developed and integrated into the CMCS multiscale chemical science simulation software. The new technologies were used to simulate HCCI engines and burner flames with exceptional fidelity.
Frequency-domain multiscale quantum mechanics/electromagnetics simulation method
Meng, Lingyi; Yin, Zhenyu; Yam, ChiYung E-mail: ghc@everest.hku.hk; Koo, SiuKong; Chen, GuanHua E-mail: ghc@everest.hku.hk; Chen, Quan; Wong, Ngai
2013-12-28
A frequency-domain quantum mechanics and electromagnetics (QM/EM) method is developed. Compared with the time-domain QM/EM method [Meng et al., J. Chem. Theory Comput. 8, 11901199 (2012)], the newly developed frequency-domain QM/EM method could effectively capture the dynamic properties of electronic devices over a broader range of operating frequencies. The system is divided into QM and EM regions and solved in a self-consistent manner via updating the boundary conditions at the QM and EM interface. The calculated potential distributions and current densities at the interface are taken as the boundary conditions for the QM and EM calculations, respectively, which facilitate the information exchange between the QM and EM calculations and ensure that the potential, charge, and current distributions are continuous across the QM/EM interface. Via Fourier transformation, the dynamic admittance calculated from the time-domain and frequency-domain QM/EM methods is compared for a carbon nanotube based molecular device.
A Multiscale Modeling Approach to Analyze Filament-Wound Composite Pressure Vessels
Nguyen, Ba Nghiep; Simmons, Kevin L.
2013-07-22
A multiscale modeling approach to analyze filament-wound composite pressure vessels is developed in this article. The approach, which extends the Nguyen et al. model [J. Comp. Mater. 43 (2009) 217] developed for discontinuous fiber composites to continuous fiber ones, spans three modeling scales. The microscale considers the unidirectional elastic fibers embedded in an elastic-plastic matrix obeying the Ramberg-Osgood relation and J2 deformation theory of plasticity. The mesoscale behavior representing the composite lamina is obtained through an incremental Mori-Tanaka type model and the Eshelby equivalent inclusion method [Proc. Roy. Soc. Lond. A241 (1957) 376]. The implementation of the micro-meso constitutive relations in the ABAQUS finite element package (via user subroutines) allows the analysis of a filament-wound composite pressure vessel (macroscale) to be performed. Failure of the composite lamina is predicted by a criterion that accounts for the strengths of the fibers and of the matrix as well as of their interface. The developed approach is demonstrated in the analysis of a filament-wound pressure vessel to study the effect of the lamina thickness on the burst pressure. The predictions are favorably compared to the numerical and experimental results by Lifshitz and Dayan [Comp. Struct. 32 (1995) 313].
Horstemeyer, Mark R.; Chaudhuri, Santanu
2015-09-30
A multiscale modeling Internal State Variable (ISV) constitutive model was developed that captures the fundamental structure-property relationships. The macroscale ISV model used lower length scale simulations (Butler-Volmer and Electronics Structures results) in order to inform the ISVs at the macroscale. The chemomechanical ISV model was calibrated and validated from experiments with magnesium (Mg) alloys that were investigated under corrosive environments coupled with experimental electrochemical studies. Because the ISV chemomechanical model is physically based, it can be used for other material systems to predict corrosion behavior. As such, others can use the chemomechanical model for analyzing corrosion effects on their designs.
Multiscale Modeling of the Deformation of Advanced Ferritic Steels for Generation IV Nuclear Energy
Nasr M. Ghoniem; Nick Kioussis
2009-04-18
The objective of this project is to use the multi-scale modeling of materials (MMM) approach to develop an improved understanding of the effects of neutron irradiation on the mechanical properties of high-temperature structural materials that are being developed or proposed for Gen IV applications. In particular, the research focuses on advanced ferritic/ martensitic steels to enable operation up to 650-700°C, compared to the current 550°C limit on high-temperature steels.
Sondak, D.; Shadid, J. N.; Oberai, A. A.; Pawlowski, R. P.; Cyr, E. C.; Smith, T. M.
2015-04-29
New large eddy simulation (LES) turbulence models for incompressible magnetohydrodynamics (MHD) derived from the variational multiscale (VMS) formulation for finite element simulations are introduced. The new models include the variational multiscale formulation, a residual-based eddy viscosity model, and a mixed model that combines both of these component models. Each model contains terms that are proportional to the residual of the incompressible MHD equations and is therefore numerically consistent. Moreover, each model is also dynamic, in that its effect vanishes when this residual is small. The new models are tested on the decaying MHD Taylor Green vortex at low and high Reynolds numbers. The evaluation of the models is based on comparisons with available data from direct numerical simulations (DNS) of the time evolution of energies as well as energy spectra at various discrete times. Thus a numerical study, on a sequence of meshes, is presented that demonstrates that the large eddy simulation approaches the DNS solution for these quantities with spatial mesh refinement.
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Sondak, D.; Shadid, J. N.; Oberai, A. A.; Pawlowski, R. P.; Cyr, E. C.; Smith, T. M.
2015-04-29
New large eddy simulation (LES) turbulence models for incompressible magnetohydrodynamics (MHD) derived from the variational multiscale (VMS) formulation for finite element simulations are introduced. The new models include the variational multiscale formulation, a residual-based eddy viscosity model, and a mixed model that combines both of these component models. Each model contains terms that are proportional to the residual of the incompressible MHD equations and is therefore numerically consistent. Moreover, each model is also dynamic, in that its effect vanishes when this residual is small. The new models are tested on the decaying MHD Taylor Green vortex at low and highmore » Reynolds numbers. The evaluation of the models is based on comparisons with available data from direct numerical simulations (DNS) of the time evolution of energies as well as energy spectra at various discrete times. Thus a numerical study, on a sequence of meshes, is presented that demonstrates that the large eddy simulation approaches the DNS solution for these quantities with spatial mesh refinement.« less
Carbon Capture Simulation Initiative: A Case Study in Multi-Scale Modeling and New Challenges
Miller, David; Syamlal, Madhava; Mebane, David; Storlie, Curtis; Bhattacharyya, Debangsu; Sahinidis, Nikolaos V.; Agarwal, Deborah A.; Tong, Charles; Zitney, Stephen E.; Sarkar, Avik; Sun, Xin; Sundaresan, Sankaran; Ryan, Emily M.; Engel, David W.; Dale, Crystal
2014-04-01
Advanced multi-scale modeling and simulation has the potential to dramatically reduce development time, resulting in considerable cost savings. The Carbon Capture Simulation Initiative is a partnership among national laboratories, industry and universities that is developing and deploying a suite of multi-scale modeling and simulation tools including basic data submodels, steady-state and dynamic process models, process optimization and uncertainty quantification tools, an advanced dynamic process control framework, high-resolution filtered computational-fluid-dynamic (CFD) submodels, validated high-fidelity device-scale CFD models with quantified uncertainty, and a risk analysis framework. These tools and models enable basic data submodels, including thermodynamics and kinetics, to be used within detailed process models to synthesize and optimize a process. The resulting process informs the development of process control systems and more detailed simulations of potential equipment to improve the design and reduce scale-up risk. Quantification and propagation of uncertainty across scales is an essential part of these tools and models.
Carbon Capture Simulation Initiative: A Case Study in Multi-Scale Modeling and New Challenges
Miller, David C; Syamlal, Madhava; Zitney, Stephen E.
2014-06-07
Abstract: Advanced multi-scale modeling and simulation has the potential to dramatically reduce development time, resulting in considerable cost savings. The Carbon Capture Simulation Initiative is a partnership among national laboratories, industry and universities that is developing and deploying a suite of multi-scale modeling and simulation tools including basic data submodels, steady-state and dynamic process models, process optimization and uncertainty quantification tools, an advanced dynamic process control framework, high-resolution filtered computational-fluid-dynamic (CFD) submodels, validated high-fidelity device-scale CFD models with quantified uncertainty, and a risk analysis framework. These tools and models enable basic data submodels, including thermodynamics and kinetics, to be used within detailed process models to synthesize and optimize a process. The resulting process informs the development of process control systems and more detailed simulations of potential equipment to improve the design and reduce scale-up risk. Quantification and propagation of uncertainty across scales is an essential part of these tools and models.
Sustainable Manufacturing via Multi-Scale, Physics-Based Process...
Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site
Sustainable Manufacturing via Multi-Scale, Physics-Based Process Modeling and ... design framework enabled by multi-scale, physics-based process models. ...
Multiscale modeling of blood-plasma separation in bifurcations | Argonne
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Leadership Computing Facility modeling of blood-plasma separation in bifurcations Authors: Li, X. J., A. S. Popel, G. E. Karniadakis Motion of a suspension of red blood cells (RBCs) flowing in a Y-shaped bifurcating microfluidic channel is investigated using a low-dimensional RBC validated 3D model based on dissipative particle dynamics. No-slip wall boundary and adaptive boundary conditions were implemented to model hydrodynamic flow within a specific wall structure of diverging
Assessment of Multi-Scale T/H Codes and Models for DNB CP
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Assessment of Multi-Scale Thermal-Hydraulic Codes and Models for DNB Challenge Problem Applications L3.AMA.CP.P8.01 Yixing Sung, Jin Yan, Liping Cao, Vefa N. Kucukboyaci, Emre Tatli Westinghouse Electric Company LLC Mark A. Christon, Jozsef Bakosi Los Alamos National Laboratories Robert K. Salko Oak Ridge National Laboratories Hongbin Zhang Idaho National Laboratory March 31, 2014 CASL-U-2014-0032-000 L3.AMA.CP.P8.01 ii CASL-U-2014-0032-000 REVISION LOG Revision Date Affected Pages Revision
Safer Batteries through Coupled Multiscale Modeling (ICCS 2015)
Turner, John A; Allu, Srikanth; Berrill, Mark A; Elwasif, Wael R; Kalnaus, Sergiy; Kumar, Abhishek; Lebrun-Grandie, Damien T; Pannala, Dr. Sreekanth; Simunovic, Srdjan
2015-01-01
Batteries are highly complex electrochemical systems, with performance and safety governed by coupled nonlinear electrochemical-electrical-thermal-mechanical processes over a range of spatiotemporal scales. We describe a new, open source computational environment for battery simulation known as VIBE - the Virtual Integrated Battery Environment. VIBE includes homogenized and pseudo-2D electrochemistry models such as those by Newman-Tiedemann-Gu (NTG) and Doyle- Fuller-Newman (DFN, a.k.a. DualFoil) as well as a new advanced capability known as AMPERES (Advanced MultiPhysics for Electrochemical and Renewable Energy Storage). AMPERES provides a 3D model for electrochemistry and full coupling with 3D electrical and thermal models on the same grid. VIBE/AMPERES has been used to create three-dimensional battery cell and pack models that explicitly simulate all the battery components (current collectors, electrodes, and separator). The models are used to predict battery performance under normal operations and to study thermal and mechanical response under adverse conditions.
Predictive Maturity of Multi-Scale Simulation Models for Fuel Performance
Atamturktur, Sez; Unal, Cetin; Hemez, Francois; Williams, Brian; Tome, Carlos
2015-03-16
The project proposed to provide a Predictive Maturity Framework with its companion metrics that (1) introduce a formalized, quantitative means to communicate information between interested parties, (2) provide scientifically dependable means to claim completion of Validation and Uncertainty Quantification (VU) activities, and (3) guide the decision makers in the allocation of Nuclear Energy’s resources for code development and physical experiments. The project team proposed to develop this framework based on two complimentary criteria: (1) the extent of experimental evidence available for the calibration of simulation models and (2) the sophistication of the physics incorporated in simulation models. The proposed framework is capable of quantifying the interaction between the required number of physical experiments and degree of physics sophistication. The project team has developed this framework and implemented it with a multi-scale model for simulating creep of a core reactor cladding. The multi-scale model is composed of the viscoplastic self-consistent (VPSC) code at the meso-scale, which represents the visco-plastic behavior and changing properties of a highly anisotropic material and a Finite Element (FE) code at the macro-scale to represent the elastic behavior and apply the loading. The framework developed takes advantage of the transparency provided by partitioned analysis, where independent constituent codes are coupled in an iterative manner. This transparency allows model developers to better understand and remedy the source of biases and uncertainties, whether they stem from the constituents or the coupling interface by exploiting separate-effect experiments conducted within the constituent domain and integral-effect experiments conducted within the full-system domain. The project team has implemented this procedure with the multi- scale VPSC-FE model and demonstrated its ability to improve the predictive capability of the model. Within this
Wang, Yuan; Wang, Minghuai; Zhang, Renyi; Ghan, Steven J.; Lin, Yun; Hu, Jiaxi; Pan, Bowen; Levy, Misti; Jiang, Jonathan; Molina, Mario J.
2014-05-13
Atmospheric aerosols impact weather and global general circulation by modifying cloud and precipitation processes, but the magnitude of cloud adjustment by aerosols remains poorly quantified and represents the largest uncertainty in estimated forcing of climate change. Here we assess the impacts of anthropogenic aerosols on the Pacific storm track using a multi-scale global aerosol-climate model (GCM). Simulations of two aerosol scenarios corresponding to the present day and pre-industrial conditions reveal long-range transport of anthropogenic aerosols across the north Pacific and large resulting changes in the aerosol optical depth, cloud droplet number concentration, and cloud and ice water paths. Shortwave and longwave cloud radiative forcing at the top of atmosphere are changed by - 2.5 and + 1.3 W m-2, respectively, by emission changes from pre-industrial to present day, and an increased cloud-top height indicates invigorated mid-latitude cyclones. The overall increased precipitation and poleward heat transport reflect intensification of the Pacific storm track by anthropogenic aerosols. Hence, this work provides for the first time a global perspective of the impacts of Asian pollution outflows from GCMs. Furthermore, our results suggest that the multi-scale modeling framework is essential in producing the aerosol invigoration effect of deep convective clouds on the global scale.
Multiscale model of metal alloy oxidation at grain boundaries
Sushko, Maria L. Alexandrov, Vitaly; Schreiber, Daniel K.; Rosso, Kevin M.; Bruemmer, Stephen M.
2015-06-07
High temperature intergranular oxidation and corrosion of metal alloys is one of the primary causes of materials degradation in nuclear systems. In order to gain insights into grain boundary oxidation processes, a mesoscale metal alloy oxidation model is established by combining quantum Density Functional Theory (DFT) and mesoscopic Poisson-Nernst-Planck/classical DFT with predictions focused on Ni alloyed with either Cr or Al. Analysis of species and fluxes at steady-state conditions indicates that the oxidation process involves vacancy-mediated transport of Ni and the minor alloying element to the oxidation front and the formation of stable metal oxides. The simulations further demonstrate that the mechanism of oxidation for Ni-5Cr and Ni-4Al is qualitatively different. Intergranular oxidation of Ni-5Cr involves the selective oxidation of the minor element and not matrix Ni, due to slower diffusion of Ni relative to Cr in the alloy and due to the significantly smaller energy gain upon the formation of nickel oxide compared to that of Cr{sub 2}O{sub 3}. This essentially one-component oxidation process results in continuous oxide formation and a monotonic Cr vacancy distribution ahead of the oxidation front, peaking at alloy/oxide interface. In contrast, Ni and Al are both oxidized in Ni-4Al forming a mixed spinel NiAl{sub 2}O{sub 4}. Different diffusivities of Ni and Al give rise to a complex elemental distribution in the vicinity of the oxidation front. Slower diffusing Ni accumulates in the oxide and metal within 3 nm of the interface, while Al penetrates deeper into the oxide phase. Ni and Al are both depleted from the region 3–10 nm ahead of the oxidation front creating voids. The oxide microstructure is also different. Cr{sub 2}O{sub 3} has a plate-like structure with 1.2–1.7 nm wide pores running along the grain boundary, while NiAl{sub 2}O{sub 4} has 1.5 nm wide pores in the direction parallel to the grain boundary and 0.6 nm pores in the perpendicular
Tourret, Damien; Clarke, Amy J.; Imhoff, Seth D.; Gibbs, Paul J.; Gibbs, John W.; Karma, Alain
2015-05-27
We present a three-dimensional extension of the multiscale dendritic needle network (DNN) model. This approach enables quantitative simulations of the unsteady dynamics of complex hierarchical networks in spatially extended dendritic arrays. We apply the model to directional solidification of Al-9.8 wt.%Si alloy and directly compare the model predictions with measurements from experiments with in situ x-ray imaging. The focus is on the dynamical selection of primary spacings over a range of growth velocities, and the influence of sample geometry on the selection of spacings. Simulation results show good agreement with experiments. The computationally efficient DNN model opens new avenues for investigating the dynamics of large dendritic arrays at scales relevant to solidification experiments and processes.
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Tourret, Damien; Clarke, Amy J.; Imhoff, Seth D.; Gibbs, Paul J.; Gibbs, John W.; Karma, Alain
2015-05-27
We present a three-dimensional extension of the multiscale dendritic needle network (DNN) model. This approach enables quantitative simulations of the unsteady dynamics of complex hierarchical networks in spatially extended dendritic arrays. We apply the model to directional solidification of Al-9.8 wt.%Si alloy and directly compare the model predictions with measurements from experiments with in situ x-ray imaging. The focus is on the dynamical selection of primary spacings over a range of growth velocities, and the influence of sample geometry on the selection of spacings. Simulation results show good agreement with experiments. The computationally efficient DNN model opens new avenues formore » investigating the dynamics of large dendritic arrays at scales relevant to solidification experiments and processes.« less
Multi-scale modeling of inter-granular fracture in UO2
Chakraborty, Pritam; Zhang, Yongfeng; Tonks, Michael R.; Biner, S. Bulent
2015-03-01
A hierarchical multi-scale approach is pursued in this work to investigate the influence of porosity, pore and grain size on the intergranular brittle fracture in UO2. In this approach, molecular dynamics simulations are performed to obtain the fracture properties for different grain boundary types. A phase-field model is then utilized to perform intergranular fracture simulations of representative microstructures with different porosities, pore and grain sizes. In these simulations the grain boundary fracture properties obtained from molecular dynamics simulations are used. The responses from the phase-field fracture simulations are then fitted with a stress-based brittle fracture model usable at the engineering scale. This approach encapsulates three different length and time scales, and allows the development of microstructurally informed engineering scale model from properties evaluated at the atomistic scale.
T.A> Buscheck; Y. Sun; Y. Hao
2006-03-28
The MultiScale ThermoHydrologic Model (MSTHM) predicts thermal-hydrologic (TH) conditions within emplacement tunnels (drifts) and in the adjoining host rock at Yucca Mountain, Nevada, which is the proposed site for a radioactive waste repository in the US. Because these predictions are used in the performance assessment of the Yucca Mountain repository, they must address the influence of variability and uncertainty of the engineered- and natural-system parameters that significantly influence those predictions. Parameter-sensitivity studies show that the MSTHM predictions adequately propagate the influence of parametric variability and uncertainty. Model-validation studies show that the influence of conceptual-model uncertainty on the MSTHM predictions is insignificant compared to that of parametric uncertainty, which is propagated through the MSTHM.
Multiscale Multiphysics Lithium-Ion Battery Model with Multidomain Modular Framework
Kim, G. H.
2013-01-01
Lithium-ion batteries (LIBs) powering recent wave of personal ubiquitous electronics are also believed to be a key enabler of electrification of vehicle powertrain on the path toward sustainable transportation future. Over the past several years, National Renewable Energy Laboratory (NREL) has developed the Multi-Scale Multi-Domain (MSMD) model framework, which is an expandable platform and a generic modularized flexible framework resolving interactions among multiple physics occurring in varied length and time scales in LIB[1]. NREL has continued to enhance the functionality of the framework and to develop constituent models in the context of the MSMD framework responding to U.S. Department of Energy's CAEBAT program objectives. This talk will introduce recent advancements in NREL's LIB modeling research in regards of scale-bridging, multi-physics integration, and numerical scheme developments.
Uncertainty Quantification and Management for Multi-scale Nuclear Materials Modeling
McDowell, David; Deo, Chaitanya; Zhu, Ting; Wang, Yan
2015-10-21
Understanding and improving microstructural mechanical stability in metals and alloys is central to the development of high strength and high ductility materials for cladding and cores structures in advanced fast reactors. Design and enhancement of radiation-induced damage tolerant alloys are facilitated by better understanding the connection of various unit processes to collective responses in a multiscale model chain, including: dislocation nucleation, absorption and desorption at interfaces; vacancy production, radiation-induced segregation of Cr and Ni at defect clusters (point defect sinks) in BCC Fe-Cr ferritic/martensitic steels; investigation of interaction of interstitials and vacancies with impurities (V, Nb, Ta, Mo, W, Al, Si, P, S); time evolution of swelling (cluster growth) phenomena of irradiated materials; and energetics and kinetics of dislocation bypass of defects formed by interstitial clustering and formation of prismatic loops, informing statistical models of continuum character with regard to processes of dislocation glide, vacancy agglomeration and swelling, climb and cross slip.
Pesaran, A.; Kim, G. H.; Smith, K.; Santhanagopalan, S.; Lee, K. J.
2012-05-01
This 2012 Annual Merit Review presentation gives an overview of the Computer-Aided Engineering of Batteries (CAEBAT) project and introduces the Multi-Scale, Multi-Dimensional model for modeling lithium-ion batteries for electric vehicles.
Anh Bui; Nam Dinh; Brian Williams
2013-09-01
In addition to validation data plan, development of advanced techniques for calibration and validation of complex multiscale, multiphysics nuclear reactor simulation codes are a main objective of the CASL VUQ plan. Advanced modeling of LWR systems normally involves a range of physico-chemical models describing multiple interacting phenomena, such as thermal hydraulics, reactor physics, coolant chemistry, etc., which occur over a wide range of spatial and temporal scales. To a large extent, the accuracy of (and uncertainty in) overall model predictions is determined by the correctness of various sub-models, which are not conservation-laws based, but empirically derived from measurement data. Such sub-models normally require extensive calibration before the models can be applied to analysis of real reactor problems. This work demonstrates a case study of calibration of a common model of subcooled flow boiling, which is an important multiscale, multiphysics phenomenon in LWR thermal hydraulics. The calibration process is based on a new strategy of model-data integration, in which, all sub-models are simultaneously analyzed and calibrated using multiple sets of data of different types. Specifically, both data on large-scale distributions of void fraction and fluid temperature and data on small-scale physics of wall evaporation were simultaneously used in this works calibration. In a departure from traditional (or common-sense) practice of tuning/calibrating complex models, a modern calibration technique based on statistical modeling and Bayesian inference was employed, which allowed simultaneous calibration of multiple sub-models (and related parameters) using different datasets. Quality of data (relevancy, scalability, and uncertainty) could be taken into consideration in the calibration process. This work presents a step forward in the development and realization of the CIPS Validation Data Plan at the Consortium for Advanced Simulation of LWRs to enable
Luo, Jian; Tomar, Vikas; Zhou, Naixie; Lee, Hongsuk
2013-06-30
Based on a recent discovery of premelting-like grain boundary segregation in refractory metals occurring at high temperatures and/or high alloying levels, this project investigated grain boundary segregation and embrittlement in tungsten (W) based alloys. Specifically, new interfacial thermodynamic models have been developed and quantified to predict high-temperature grain boundary segregation in the W-Ni binary alloy and W-Ni-Fe, W-Ni-Ti, W-Ni-Co, W-Ni-Cr, W-Ni-Zr and W-Ni-Nb ternary alloys. The thermodynamic modeling results have been experimentally validated for selected systems. Furthermore, multiscale modeling has been conducted at continuum, atomistic and quantum-mechanical levels to link grain boundary segregation with embrittlement. In summary, this 3-year project has successfully developed a theoretical framework in combination with a multiscale modeling strategy for predicting grain boundary segregation and embrittlement in W based alloys.
Discharge Performance of Li-O_{2} Batteries Using a Multiscale Modeling Approach
Bao, Jie; Xu, Wu; Bhattacharya, Priyanka; Stewart, Mark L.; Zhang, Jiguang; Pan, Wenxiao
2015-06-10
To study the discharge performance of Li–O_{2} batteries, we propose a multiscale modeling framework that links models in an upscaling fashion from the nanoscale to mesoscale and finally to the device scale. We have effectively reconstructed the microstructure of a Li–O_{2} air electrode in silico, conserving the porosity, surface-to-volume ratio, and pore size distribution of the real air electrode structure. The mechanism of rate-dependent morphology of Li_{2}O_{2} growth is incorporated into the mesoscale model. The correlation between the active-surface-to-volume ratio and averaged Li_{2}O_{2} concentration is derived to link different scales. The proposed approach’s accuracy is first demonstrated by comparing the predicted discharge curves of Li–O_{2} batteries with experimental results at the high current density. Next, the validated modeling approach effectively captures the significant improvement in discharge capacity due to the formation of Li_{2}O_{2} particles. Finally, it predicts the discharge capacities of Li–O_{2} batteries with different air electrode microstructure designs and operating conditions.
Investigating ice nucleation in cirrus clouds with an aerosol-enabled Multiscale Modeling Framework
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Zhang, Chengzhu; Wang, Minghuai; Morrison, H.; Somerville, Richard C.; Zhang, Kai; Liu, Xiaohong; Li, J-L F.
2014-11-06
In this study, an aerosol-dependent ice nucleation scheme [Liu and Penner, 2005] has been implemented in an aerosol-enabled multi-scale modeling framework (PNNL MMF) to study ice formation in upper troposphere cirrus clouds through both homogeneous and heterogeneous nucleation. The MMF model represents cloud scale processes by embedding a cloud-resolving model (CRM) within each vertical column of a GCM grid. By explicitly linking ice nucleation to aerosol number concentration, CRM-scale temperature, relative humidity and vertical velocity, the new MMF model simulates the persistent high ice supersaturation and low ice number concentration (10 to 100/L) at cirrus temperatures. The low ice numbermore » is attributed to the dominance of heterogeneous nucleation in ice formation. The new model simulates the observed shift of the ice supersaturation PDF towards higher values at low temperatures following homogeneous nucleation threshold. The MMF models predict a higher frequency of midlatitude supersaturation in the Southern hemisphere and winter hemisphere, which is consistent with previous satellite and in-situ observations. It is shown that compared to a conventional GCM, the MMF is a more powerful model to emulate parameters that evolve over short time scales such as supersaturation. Sensitivity tests suggest that the simulated global distribution of ice clouds is sensitive to the ice nucleation schemes and the distribution of sulfate and dust aerosols. Simulations are also performed to test empirical parameters related to auto-conversion of ice crystals to snow. Results show that with a value of 250 μm for the critical diameter, Dcs, that distinguishes ice crystals from snow, the model can produce good agreement to the satellite retrieved products in terms of cloud ice water path and ice water content, while the total ice water is not sensitive to the specification of Dcs value.« less
Investigating ice nucleation in cirrus clouds with an aerosol-enabled Multiscale Modeling Framework
Zhang, Chengzhu; Wang, Minghuai; Morrison, H.; Somerville, Richard C.; Zhang, Kai; Liu, Xiaohong; Li, J-L F.
2014-12-01
In this study, an aerosol-dependent ice nucleation scheme [Liu and Penner, 2005] has been implemented in an aerosol-enabled multi-scale modeling framework (PNNL MMF) to study ice formation in upper troposphere cirrus clouds through both homogeneous and heterogeneous nucleation. The MMF model represents cloud scale processes by embedding a cloud-resolving model (CRM) within each vertical column of a GCM grid. By explicitly linking ice nucleation to aerosol number concentration, CRM-scale temperature, relative humidity and vertical velocity, the new MMF model simulates the persistent high ice supersaturation and low ice number concentration (10 to 100/L) at cirrus temperatures. The low ice number is attributed to the dominance of heterogeneous nucleation in ice formation. The new model simulates the observed shift of the ice supersaturation PDF towards higher values at low temperatures following homogeneous nucleation threshold. The MMF models predict a higher frequency of midlatitude supersaturation in the Southern hemisphere and winter hemisphere, which is consistent with previous satellite and in-situ observations. It is shown that compared to a conventional GCM, the MMF is a more powerful model to emulate parameters that evolve over short time scales such as supersaturation. Sensitivity tests suggest that the simulated global distribution of ice clouds is sensitive to the ice nucleation schemes and the distribution of sulfate and dust aerosols. Simulations are also performed to test empirical parameters related to auto-conversion of ice crystals to snow. Results show that with a value of 250 ?m for the critical diameter, Dcs, that distinguishes ice crystals from snow, the model can produce good agreement to the satellite retrieved products in terms of cloud ice water path and ice water content, while the total ice water is not sensitive to the specification of Dcs value.
Investigating ice nucleation in cirrus clouds with an aerosol-enabled Multiscale Modeling Framework
Zhang, Chengzhu; Wang, Minghuai; Morrison, H.; Somerville, Richard C.; Zhang, Kai; Liu, Xiaohong; Li, J-L F.
2014-11-06
In this study, an aerosol-dependent ice nucleation scheme [Liu and Penner, 2005] has been implemented in an aerosol-enabled multi-scale modeling framework (PNNL MMF) to study ice formation in upper troposphere cirrus clouds through both homogeneous and heterogeneous nucleation. The MMF model represents cloud scale processes by embedding a cloud-resolving model (CRM) within each vertical column of a GCM grid. By explicitly linking ice nucleation to aerosol number concentration, CRM-scale temperature, relative humidity and vertical velocity, the new MMF model simulates the persistent high ice supersaturation and low ice number concentration (10 to 100/L) at cirrus temperatures. The low ice number is attributed to the dominance of heterogeneous nucleation in ice formation. The new model simulates the observed shift of the ice supersaturation PDF towards higher values at low temperatures following homogeneous nucleation threshold. The MMF models predict a higher frequency of midlatitude supersaturation in the Southern hemisphere and winter hemisphere, which is consistent with previous satellite and in-situ observations. It is shown that compared to a conventional GCM, the MMF is a more powerful model to emulate parameters that evolve over short time scales such as supersaturation. Sensitivity tests suggest that the simulated global distribution of ice clouds is sensitive to the ice nucleation schemes and the distribution of sulfate and dust aerosols. Simulations are also performed to test empirical parameters related to auto-conversion of ice crystals to snow. Results show that with a value of 250 ?m for the critical diameter, Dcs, that distinguishes ice crystals from snow, the model can produce good agreement to the satellite retrieved products in terms of cloud ice water path and ice water content, while the total ice water is not sensitive to the specification of Dcs value.
Toward Multi-scale Modeling and simulation of conduction in heterogeneous materials.
Lechman, Jeremy B.; Battaile, Corbett Chandler.; Bolintineanu, Dan; Cooper, Marcia A.; Erikson, William W.; Foiles, Stephen M.; Kay, Jeffrey J; Phinney, Leslie M.; Piekos, Edward S.; Specht, Paul Elliott; Wixom, Ryan R.; Yarrington, Cole
2015-01-01
This report summarizes a project in which the authors sought to develop and deploy: (i) experimental techniques to elucidate the complex, multiscale nature of thermal transport in particle-based materials; and (ii) modeling approaches to address current challenges in predicting performace variability of materials (e.g., identifying and characterizing physical- chemical processes and their couplings across multiple length and time scales, modeling infor- mation transfer between scales, and statically and dynamically resolving material structure and its evolution during manufacturing and device performance). Experimentally, several capabilities were sucessfully advanced. As discussed in Chapter 2 a flash diffusivity capabil- ity for measuring homogeneous thermal conductivity of pyrotechnic powders (and beyond) was advanced; leading to enhanced characterization of pyrotechnic materials and properties impacting component development. Chapter 4 describes sucess for the first time, although preliminary, in resolving thermal fields at speeds and spatial scales relevant to energetic components. Chapter 7 summarizes the first ever (as far as the authors know) application of TDTR to actual pyrotechnic materials. This is the first attempt to actually characterize these materials at the interfacial scale. On the modeling side, new capabilities in image processing of experimental microstructures and direct numerical simulation on complicated structures were advanced (see Chapters 3 and 5). In addition, modeling work described in Chapter 8 led to improved prediction of interface thermal conductance from first principles calculations. Toward the second point, for a model system of packed particles, significant headway was made in implementing numerical algorithms and collecting data to justify the approach in terms of highlighting the phenomena at play and pointing the way forward in de- veloping and informing the kind of modeling approach oringinally envisioned (see Chapter 6
A Unified Multi-Scale Model for Pore-Scale Flow Simulations in Soils
Yang, Xiaofan; Liu, Chongxuan; Shang, Jianying; Fang, Yilin; Bailey, Vanessa L.
2014-01-30
Pore-scale simulations have received increasing interest in subsurface sciences to provide mechanistic insights into the macroscopic phenomena of water flow and reactive transport processes. The application of the pore scale simulations to soils and sediments is, however, challenged because of the characterization limitation that often only allows partial resolution of pore structure and geometry. A significant proportion of the pore space in soils and sediments is below the spatial resolution, forming a mixed media of pore and porous domains. Here we reported a unified multi-scale model (UMSM) that can be used to simulate water flow and transport in mixed media of pore and porous domains under both saturated and unsaturated conditions. The approach modifies the classic Navier-Stokes equation by adding a Darcy term to describe fluid momentum and uses a generalized mass balance equation for saturated and unsaturated conditions. By properly defining physical parameters, the UMSM can be applied in both pore and porous domains. This paper describes the set of equations for the UMSM, a series of validation cases under saturated or unsaturated conditions, and a real soil case for the application of the approach.
MULTI-SCALE MODELING AND APPROXIMATION ASSISTED OPTIMIZATION OF BARE TUBE HEAT EXCHANGERS
Bacellar, Daniel; Ling, Jiazhen; Aute, Vikrant; Radermacher, Reinhard; Abdelaziz, Omar
2014-01-01
Air-to-refrigerant heat exchangers are very common in air-conditioning, heat pump and refrigeration applications. In these heat exchangers, there is a great benefit in terms of size, weight, refrigerant charge and heat transfer coefficient, by moving from conventional channel sizes (~ 9mm) to smaller channel sizes (< 5mm). This work investigates new designs for air-to-refrigerant heat exchangers with tube outer diameter ranging from 0.5 to 2.0mm. The goal of this research is to develop and optimize the design of these heat exchangers and compare their performance with existing state of the art designs. The air-side performance of various tube bundle configurations are analyzed using a Parallel Parameterized CFD (PPCFD) technique. PPCFD allows for fast-parametric CFD analyses of various geometries with topology change. Approximation techniques drastically reduce the number of CFD evaluations required during optimization. Maximum Entropy Design method is used for sampling and Kriging method is used for metamodeling. Metamodels are developed for the air-side heat transfer coefficients and pressure drop as a function of tube-bundle dimensions and air velocity. The metamodels are then integrated with an air-to-refrigerant heat exchanger design code. This integration allows a multi-scale analysis of air-side performance heat exchangers including air-to-refrigerant heat transfer and phase change. Overall optimization is carried out using a multi-objective genetic algorithm. The optimal designs found can exhibit 50 percent size reduction, 75 percent decrease in air side pressure drop and doubled air heat transfer coefficients compared to a high performance compact micro channel heat exchanger with same capacity and flow rates.
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Wei, Yaxing; Liu, Shishi; Huntzinger, Deborah N.; Michalak, Anna M.; Viovy, Nicolas; Post, Wilfred M.; Schwalm, Christopher R.; Schaeffer, Kevin; Jacobson, Andrew R.; Lu, Chaoqun; et al
2014-12-05
Ecosystems are important and dynamic components of the global carbon cycle, and terrestrial biospheric models (TBMs) are crucial tools in further understanding of how terrestrial carbon is stored and exchanged with the atmosphere across a variety of spatial and temporal scales. Improving TBM skills, and quantifying and reducing their estimation uncertainties, pose significant challenges. The Multi-scale Synthesis and Terrestrial Model Intercomparison Project (MsTMIP) is a formal multi-scale and multi-model intercomparison effort set up to tackle these challenges. The MsTMIP protocol prescribes standardized environmental driver data that are shared among model teams to facilitate model model and model observation comparisons. Inmore » this article, we describe the global and North American environmental driver data sets prepared for the MsTMIP activity to both support their use in MsTMIP and make these data, along with the processes used in selecting/processing these data, accessible to a broader audience. Based on project needs and lessons learned from past model intercomparison activities, we compiled climate, atmospheric CO2 concentrations, nitrogen deposition, land use and land cover change (LULCC), C3 / C4 grasses fractions, major crops, phenology and soil data into a standard format for global (0.5⁰ x 0.5⁰ resolution) and regional (North American: 0.25⁰ x 0.25⁰ resolution) simulations. In order to meet the needs of MsTMIP, improvements were made to several of the original environmental data sets, by improving the quality, and/or changing their spatial and temporal coverage, and resolution. The resulting standardized model driver data sets are being used by over 20 different models participating in MsTMIP. Lastly, the data are archived at the Oak Ridge National Laboratory Distributed Active Archive Center (ORNL DAAC, http://daac.ornl.gov) to provide long-term data management and distribution.« less
Huntzinger, D.N.; Schwalm, C.; Michalak, A.M; Schaefer, K.; King, A.W.; Wei, Y.; Jacobson, A.; Liu, S.; Cook, R.; Post, W.M.; Berthier, G.; Hayes, D.; Huang, M.; Ito, A.; Lei, H.; Lu, C.; Mao, J.; Peng, C.H.; Peng, S.; Poulter, B.; Riccuito, D.; Shi, X.; Tian, H.; Wang, W.; Zeng, N.; Zhao, F.; Zhu, Q.
2013-01-01
Terrestrial biosphere models (TBMs) have become an integral tool for extrapolating local observations and understanding of land-atmosphere carbon exchange to larger regions. The North American Carbon Program (NACP) Multi-scale synthesis and Terrestrial Model Intercomparison Project (MsTMIP) is a formal model intercomparison and evaluation effort focused on improving the diagnosis and attribution of carbon exchange at regional and global scales. MsTMIP builds upon current and past synthesis activities, and has a unique framework designed to isolate, interpret, and inform understanding of how model structural differences impact estimates of carbon uptake and release. Here we provide an overview of the MsTMIP effort and describe how the MsTMIP experimental design enables the assessment and quantification of TBM structural uncertainty. Model structure refers to the types of processes considered (e.g. nutrient cycling, disturbance, lateral transport of carbon), and how these processes are represented (e.g. photosynthetic formulation, temperature sensitivity, respiration) in the models. By prescribing a common experimental protocol with standard spin-up procedures and driver data sets, we isolate any biases and variability in TBM estimates of regional and global carbon budgets resulting from differences in the models themselves (i.e. model structure) and model-specific parameter values. An initial intercomparison of model structural differences is represented using hierarchical cluster diagrams (a.k.a. dendrograms), which highlight similarities and differences in how models account for carbon cycle, vegetation, energy, and nitrogen cycle dynamics. We show that, despite the standardized protocol used to derive initial conditions, models show a high degree of variation for GPP, total living biomass, and total soil carbon, underscoring the influence of differences in model structure and parameterization on model estimates.
Multiscale modeling and characterization for performance and safety of lithium-ion batteries
Pannala, Sreekanth; Turner, John A.; Allu, Srikanth; Elwasif, Wael R.; Kalnaus, Sergiy; Simunovic, Srdjan; Kumar, Abhishek; Billings, Jay Jay; Wang, Hsin; Nanda, Jagjit
2015-08-19
Lithium-ion batteries are highly complex electrochemical systems whose performance and safety are governed by coupled nonlinear electrochemical-electrical-thermal-mechanical processes over a range of spatiotemporal scales. In this paper we describe a new, open source computational framework for Lithium-ion battery simulations that is designed to support a variety of model types and formulations. This framework has been used to create three-dimensional cell and battery pack models that explicitly simulate all the battery components (current collectors, electrodes, and separator). The models are used to predict battery performance under normal operations and to study thermal and mechanical safety aspects under adverse conditions. The model development and validation are supported by experimental methods such as IR-imaging, X-ray tomography and micro-Raman mapping.
Multiscale modeling and characterization for performance and safety of lithium-ion batteries
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Pannala, Sreekanth; Turner, John A.; Allu, Srikanth; Elwasif, Wael R.; Kalnaus, Sergiy; Simunovic, Srdjan; Kumar, Abhishek; Billings, Jay Jay; Wang, Hsin; Nanda, Jagjit
2015-08-19
Lithium-ion batteries are highly complex electrochemical systems whose performance and safety are governed by coupled nonlinear electrochemical-electrical-thermal-mechanical processes over a range of spatiotemporal scales. In this paper we describe a new, open source computational framework for Lithium-ion battery simulations that is designed to support a variety of model types and formulations. This framework has been used to create three-dimensional cell and battery pack models that explicitly simulate all the battery components (current collectors, electrodes, and separator). The models are used to predict battery performance under normal operations and to study thermal and mechanical safety aspects under adverse conditions. The modelmore » development and validation are supported by experimental methods such as IR-imaging, X-ray tomography and micro-Raman mapping.« less
Welz, Oliver; Burke, Michael P.; Antonov, Ivan O.; Goldsmith, C. Franklin; Savee, John David; Osborn, David L.; Taatjes, Craig A.; Klippenstein, Stephen J.; Sheps, Leonid
2015-04-10
We studied low-temperature propane oxidation at P = 4 Torr and T = 530, 600, and 670 K by time-resolved multiplexed photoionization mass spectrometry (MPIMS), which probes the reactants, intermediates, and products with isomeric selectivity using tunable synchrotron vacuum UV ionizing radiation. The oxidation is initiated by pulsed laser photolysis of oxalyl chloride, (COCl)_{2}, at 248 nm, which rapidly generates a ~1:1 mixture of 1-propyl (n-propyl) and 2-propyl (i-propyl) radicals via the fast Cl + propane reaction. At all three temperatures, the major stable product species is propene, formed in the propyl + O_{2} reactions by direct HO_{2} elimination from both n- and i-propyl peroxy radicals. The experimentally derived propene yields relative to the initial concentration of Cl atoms are (20 4)% at 530 K, (55 11)% at 600 K, and (86 17)% at 670 K at a reaction time of 20 ms. The lower yield of propene at low temperature reflects substantial formation of propyl peroxy radicals, which do not completely decompose on the experimental time scale. In addition, C_{3}H_{6}O isomers methyloxirane, oxetane, acetone, and propanal are detected as minor products. Our measured yields of oxetane and methyloxirane, which are coproducts of OH radicals, suggest a revision of the OH formation pathways in models of low-temperature propane oxidation. The experimental results are modeled and interpreted using a multiscale informatics approach, presented in detail in a separate publication (Burke, M. P.; Goldsmith, C. F.; Klippenstein, S. J.; Welz, O.; Huang H.; Antonov I. O.; Savee J. D.; Osborn D. L.; Zdor, J.; Taatjes, C. A.; Sheps, L. Multiscale Informatics for Low-Temperature Propane Oxidation: Further Complexities in Studies of Complex Reactions. J. Phys. Chem A. 2015, DOI: 10.1021/acs.jpca.5b01003). Additionally, we found that the model predicts the time profiles and yields of the experimentally observed primary products well, and
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Welz, Oliver; Burke, Michael P.; Antonov, Ivan O.; Goldsmith, C. Franklin; Savee, John David; Osborn, David L.; Taatjes, Craig A.; Klippenstein, Stephen J.; Sheps, Leonid
2015-04-10
We studied low-temperature propane oxidation at P = 4 Torr and T = 530, 600, and 670 K by time-resolved multiplexed photoionization mass spectrometry (MPIMS), which probes the reactants, intermediates, and products with isomeric selectivity using tunable synchrotron vacuum UV ionizing radiation. The oxidation is initiated by pulsed laser photolysis of oxalyl chloride, (COCl)2, at 248 nm, which rapidly generates a ~1:1 mixture of 1-propyl (n-propyl) and 2-propyl (i-propyl) radicals via the fast Cl + propane reaction. At all three temperatures, the major stable product species is propene, formed in the propyl + O2 reactions by direct HO2 elimination frommore » both n- and i-propyl peroxy radicals. The experimentally derived propene yields relative to the initial concentration of Cl atoms are (20 ± 4)% at 530 K, (55 ± 11)% at 600 K, and (86 ± 17)% at 670 K at a reaction time of 20 ms. The lower yield of propene at low temperature reflects substantial formation of propyl peroxy radicals, which do not completely decompose on the experimental time scale. In addition, C3H6O isomers methyloxirane, oxetane, acetone, and propanal are detected as minor products. Our measured yields of oxetane and methyloxirane, which are coproducts of OH radicals, suggest a revision of the OH formation pathways in models of low-temperature propane oxidation. The experimental results are modeled and interpreted using a multiscale informatics approach, presented in detail in a separate publication (Burke, M. P.; Goldsmith, C. F.; Klippenstein, S. J.; Welz, O.; Huang H.; Antonov I. O.; Savee J. D.; Osborn D. L.; Zádor, J.; Taatjes, C. A.; Sheps, L. Multiscale Informatics for Low-Temperature Propane Oxidation: Further Complexities in Studies of Complex Reactions. J. Phys. Chem A. 2015, DOI: 10.1021/acs.jpca.5b01003). Additionally, we found that the model predicts the time profiles and yields of the experimentally observed primary products well, and shows satisfactory agreement for products
Welz, Oliver; Burke, Michael P.; Antonov, Ivan O.; Goldsmith, C. Franklin; Savee, John David; Osborn, David L.; Taatjes, Craig A.; Klippenstein, Stephen J.; Sheps, Leonid
2015-04-10
We studied low-temperature propane oxidation at P = 4 Torr and T = 530, 600, and 670 K by time-resolved multiplexed photoionization mass spectrometry (MPIMS), which probes the reactants, intermediates, and products with isomeric selectivity using tunable synchrotron vacuum UV ionizing radiation. The oxidation is initiated by pulsed laser photolysis of oxalyl chloride, (COCl)_{2}, at 248 nm, which rapidly generates a ~1:1 mixture of 1-propyl (n-propyl) and 2-propyl (i-propyl) radicals via the fast Cl + propane reaction. At all three temperatures, the major stable product species is propene, formed in the propyl + O_{2} reactions by direct HO_{2} elimination from both n- and i-propyl peroxy radicals. The experimentally derived propene yields relative to the initial concentration of Cl atoms are (20 ± 4)% at 530 K, (55 ± 11)% at 600 K, and (86 ± 17)% at 670 K at a reaction time of 20 ms. The lower yield of propene at low temperature reflects substantial formation of propyl peroxy radicals, which do not completely decompose on the experimental time scale. In addition, C_{3}H_{6}O isomers methyloxirane, oxetane, acetone, and propanal are detected as minor products. Our measured yields of oxetane and methyloxirane, which are coproducts of OH radicals, suggest a revision of the OH formation pathways in models of low-temperature propane oxidation. The experimental results are modeled and interpreted using a multiscale informatics approach, presented in detail in a separate publication (Burke, M. P.; Goldsmith, C. F.; Klippenstein, S. J.; Welz, O.; Huang H.; Antonov I. O.; Savee J. D.; Osborn D. L.; Zádor, J.; Taatjes, C. A.; Sheps, L. Multiscale Informatics for Low-Temperature Propane Oxidation: Further Complexities in Studies of Complex Reactions. J. Phys. Chem A. 2015, DOI: 10.1021/acs.jpca.5b01003). Additionally, we found that the model predicts the time profiles and yields of the experimentally observed primary products well
Wang, Minghuai; Larson, Vincent E.; Ghan, Steven J.; Ovchinnikov, Mikhail; Schanen, D.; Xiao, Heng; Liu, Xiaohong; Rasch, Philip J.; Guo, Zhun
2015-06-01
In this study, a higher-order turbulence closure scheme, called Cloud Layers Unified by Binormals (CLUBB), is implemented into a Multi-scale Modeling Framework (MMF) model to improve low cloud simulations. The performance of CLUBB in MMF simulations with two different microphysics configurations (one-moment cloud microphysics without aerosol treatment and two-moment cloud microphysics coupled with aerosol treatment) is evaluated against observations and further compared with results from the Community Atmosphere Model, Version 5 (CAM5) with conventional cloud parameterizations. CLUBB is found to improve low cloud simulations in the MMF, and the improvement is particularly evident in the stratocumulus-to-cumulus transition regions. Compared to the single-moment cloud microphysics, CLUBB with two-moment microphysics produces clouds that are closer to the coast, and agrees better with observations. In the stratocumulus-to cumulus transition regions, CLUBB with two-moment cloud microphysics produces shortwave cloud forcing in better agreement with observations, while CLUBB with single moment cloud microphysics overestimates shortwave cloud forcing. CLUBB is further found to produce quantitatively similar improvements in the MMF and CAM5, with slightly better performance in the MMF simulations (e.g., MMF with CLUBB generally produces low clouds that are closer to the coast than CAM5 with CLUBB). Improved low cloud simulations in MMF make it an even more attractive tool for studying aerosol-cloud-precipitation interactions.
A multiscale MDCT image-based breathing lung model with time-varying regional ventilation
Yin, Youbing, E-mail: youbing-yin@uiowa.edu [Department of Mechanical and Industrial Engineering, The University of Iowa, Iowa City, IA 52242 (United States) [Department of Mechanical and Industrial Engineering, The University of Iowa, Iowa City, IA 52242 (United States); IIHR-Hydroscience and Engineering, The University of Iowa, Iowa City, IA 52242 (United States); Department of Radiology, The University of Iowa, Iowa City, IA 52242 (United States); Choi, Jiwoong, E-mail: jiwoong-choi@uiowa.edu [Department of Mechanical and Industrial Engineering, The University of Iowa, Iowa City, IA 52242 (United States) [Department of Mechanical and Industrial Engineering, The University of Iowa, Iowa City, IA 52242 (United States); IIHR-Hydroscience and Engineering, The University of Iowa, Iowa City, IA 52242 (United States); Hoffman, Eric A., E-mail: eric-hoffman@uiowa.edu [Department of Radiology, The University of Iowa, Iowa City, IA 52242 (United States); Department of Biomedical Engineering, The University of Iowa, Iowa City, IA 52242 (United States); Department of Internal Medicine, The University of Iowa, Iowa City, IA 52242 (United States); Tawhai, Merryn H., E-mail: m.tawhai@auckland.ac.nz [Auckland Bioengineering Institute, The University of Auckland, Auckland (New Zealand); Lin, Ching-Long, E-mail: ching-long-lin@uiowa.edu [Department of Mechanical and Industrial Engineering, The University of Iowa, Iowa City, IA 52242 (United States) [Department of Mechanical and Industrial Engineering, The University of Iowa, Iowa City, IA 52242 (United States); IIHR-Hydroscience and Engineering, The University of Iowa, Iowa City, IA 52242 (United States)
2013-07-01
A novel algorithm is presented that links local structural variables (regional ventilation and deforming central airways) to global function (total lung volume) in the lung over three imaged lung volumes, to derive a breathing lung model for computational fluid dynamics simulation. The algorithm constitutes the core of an integrative, image-based computational framework for subject-specific simulation of the breathing lung. For the first time, the algorithm is applied to three multi-detector row computed tomography (MDCT) volumetric lung images of the same individual. A key technique in linking global and local variables over multiple images is an in-house mass-preserving image registration method. Throughout breathing cycles, cubic interpolation is employed to ensure C{sub 1} continuity in constructing time-varying regional ventilation at the whole lung level, flow rate fractions exiting the terminal airways, and airway deformation. The imaged exit airway flow rate fractions are derived from regional ventilation with the aid of a three-dimensional (3D) and one-dimensional (1D) coupled airway tree that connects the airways to the alveolar tissue. An in-house parallel large-eddy simulation (LES) technique is adopted to capture turbulent-transitional-laminar flows in both normal and deep breathing conditions. The results obtained by the proposed algorithm when using three lung volume images are compared with those using only one or two volume images. The three-volume-based lung model produces physiologically-consistent time-varying pressure and ventilation distribution. The one-volume-based lung model under-predicts pressure drop and yields un-physiological lobar ventilation. The two-volume-based model can account for airway deformation and non-uniform regional ventilation to some extent, but does not capture the non-linear features of the lung.
Grell, Georg; Fast, Jerome D.; Gustafson, William I.; Peckham, Steven E.; McKeen, Stuart A.; Salzmann, Marc; Freitas, Saulo
2010-01-01
This is a conference proceeding that is now being put together as a book. This is chapter 2 of the book: "INTEGRATED SYSTEMS OF MESO-METEOROLOGICAL AND CHEMICAL TRANSPORT MODELS" published by Springer. The chapter title is "On-line Chemistry within WRF: Description and Evaluation of a State-of-the-Art Multiscale Air Quality and Weather Prediction Model." The original conference was the COST-728/NetFAM workshop on Integrated systems of meso-meteorological and chemical transport models, Danish Meteorological Institute, Copenhagen, May 21-23, 2007.
Freed, Alan D.; Einstein, Daniel R.; Carson, James P.; Jacob, Rick E.
2012-03-01
In the first year of this contractual effort a hypo-elastic constitutive model was developed and shown to have great potential in modeling the elastic response of parenchyma. This model resides at the macroscopic level of the continuum. In this, the second year of our support, an isotropic dodecahedron is employed as an alveolar model. This is a microscopic model for parenchyma. A hopeful outcome is that the linkage between these two scales of modeling will be a source of insight and inspiration that will aid us in the final year's activity: creating a viscoelastic model for parenchyma.
Multiscale Simulation Framework for Coupled Fluid Flow and Mechanical Deformation
Tchelepi, Hamdi
2014-11-14
A multiscale linear-solver framework for the pressure equation associated with flow in highly heterogeneous porous formations was developed. The multiscale based approach is cast in a general algebraic form, which facilitates integration of the new scalable linear solver in existing flow simulators. The Algebraic Multiscale Solver (AMS) is employed as a preconditioner within a multi-stage strategy. The formulations investigated include the standard MultiScale Finite-Element (MSFE) andMultiScale Finite-Volume (MSFV) methods. The local-stage solvers include incomplete factorization and the so-called Correction Functions (CF) associated with the MSFV approach. Extensive testing of AMS, as an iterative linear solver, indicate excellent convergence rates and computational scalability. AMS compares favorably with advanced Algebraic MultiGrid (AMG) solvers for highly detailed three-dimensional heterogeneous models. Moreover, AMS is expected to be especially beneficial in solving time-dependent problems of coupled multiphase flow and transport in large-scale subsurface formations.
Behafarid, F.; Shaver, D. R.; Bolotnov, I. A.; Jansen, K. E.; Antal, S. P.; Podowski, M. Z.
2012-07-01
The required technological and safety standards for future Gen IV Reactors can only be achieved if advanced simulation capabilities become available, which combine high performance computing with the necessary level of modeling detail and high accuracy of predictions. The purpose of this paper is to present new results of multi-scale three-dimensional (3D) simulations of the inter-related phenomena, which occur as a result of fuel element heat-up and cladding failure, including the injection of a jet of gaseous fission products into a partially blocked Sodium Fast Reactor (SFR) coolant channel, and gas/molten sodium transport along the coolant channels. The computational approach to the analysis of the overall accident scenario is based on using two different inter-communicating computational multiphase fluid dynamics (CMFD) codes: a CFD code, PHASTA, and a RANS code, NPHASE-CMFD. Using the geometry and time history of cladding failure and the gas injection rate, direct numerical simulations (DNS), combined with the Level Set method, of two-phase turbulent flow have been performed by the PHASTA code. The model allows one to track the evolution of gas/liquid interfaces at a centimeter scale. The simulated phenomena include the formation and breakup of the jet of fission products injected into the liquid sodium coolant. The PHASTA outflow has been averaged over time to obtain mean phasic velocities and volumetric concentrations, as well as the liquid turbulent kinetic energy and turbulence dissipation rate, all of which have served as the input to the core-scale simulations using the NPHASE-CMFD code. A sliding window time averaging has been used to capture mean flow parameters for transient cases. The results presented in the paper include testing and validation of the proposed models, as well the predictions of fission-gas/liquid-sodium transport along a multi-rod fuel assembly of SFR during a partial loss-of-flow accident. (authors)
Modeling, Analysis and Simulation of Multiscale Preferential Flow - 8/05-8/10 - Final Report
Ralph Showalter; Malgorzata Peszynska
2012-07-03
The research agenda of this project are: (1) Modeling of preferential transport from mesoscale to macroscale; (2) Modeling of fast flow in narrow fractures in porous media; (3) Pseudo-parabolic Models of Dynamic Capillary Pressure; (4) Adaptive computational upscaling of flow with inertia from porescale to mesoscale; (5) Adaptive modeling of nonlinear coupled systems; and (6) Adaptive modeling and a-posteriori estimators for coupled systems with heterogeneous data.
Multi-scale Modelling of bcc-Fe Based Alloys for Nuclear Applications
Malerba, Lorenzo
2008-07-01
, advanced techniques to fit interatomic potentials consistent with thermodynamics are proposed and the results of their application to the mentioned alloys are presented. Next, the development of advanced methods, based on the use of artificial intelligence, to improve both the physical reliability and the computational efficiency of kinetic Monte Carlo codes for the study of point-defect clustering and phase changes beyond the scale of MD, is reported. These recent progresses bear the promise of being able, in the near future, of producing reliable tools for the description of the microstructure evolution of realistic model alloys under irradiation. (author)
Costigan, Keeley Rochelle; Sauer, Jeremy A.; Dubey, Manvendra Krishna
2015-07-10
This report discusses the ghgas IC project which when applied, allows for an evaluation of LANL's HIGRAD model which can be used to create atmospheric simulations.
Evaluation of the Multi-scale Modeling Framework Using Data from...
Office of Scientific and Technical Information (OSTI)
Unfortunately, the traditional parametric approach of diagnosing cloud and radiation properties for gridcells that are tens to hundreds kilometers across from large-scale model ...
The Radiative Properties of Small Clouds: Multi-Scale Observations and Modeling
Feingold, Graham; McComiskey, Allison
2013-09-25
Warm, liquid clouds and their representation in climate models continue to represent one of the most significant unknowns in climate sensitivity and climate change. Our project combines ARM observations, LES modeling, and satellite imagery to characterize shallow clouds and the role of aerosol in modifying their radiative effects.
Multiscale modeling and characterization for performance and safety of lithium-ion batteries
Pannala, S. Turner, J. A.; Allu, S.; Elwasif, W. R.; Kalnaus, S.; Simunovic, S.; Kumar, A.; Billings, J. J.; Wang, H.; Nanda, J.
2015-08-21
Lithium-ion batteries are highly complex electrochemical systems whose performance and safety are governed by coupled nonlinear electrochemical-electrical-thermal-mechanical processes over a range of spatiotemporal scales. Gaining an understanding of the role of these processes as well as development of predictive capabilities for design of better performing batteries requires synergy between theory, modeling, and simulation, and fundamental experimental work to support the models. This paper presents the overview of the work performed by the authors aligned with both experimental and computational efforts. In this paper, we describe a new, open source computational environment for battery simulations with an initial focus on lithium-ion systems but designed to support a variety of model types and formulations. This system has been used to create a three-dimensional cell and battery pack models that explicitly simulate all the battery components (current collectors, electrodes, and separator). The models are used to predict battery performance under normal operations and to study thermal and mechanical safety aspects under adverse conditions. This paper also provides an overview of the experimental techniques to obtain crucial validation data to benchmark the simulations at various scales for performance as well as abuse. We detail some initial validation using characterization experiments such as infrared and neutron imaging and micro-Raman mapping. In addition, we identify opportunities for future integration of theory, modeling, and experiments.
Multiscale Universal Interface: A concurrent framework for coupling heterogeneous solvers
Tang, Yu-Hang; Kudo, Shuhei; Bian, Xin; Li, Zhen; Karniadakis, George Em
2015-09-15
Graphical abstract: - Abstract: Concurrently coupled numerical simulations using heterogeneous solvers are powerful tools for modeling multiscale phenomena. However, major modifications to existing codes are often required to enable such simulations, posing significant difficulties in practice. In this paper we present a C++ library, i.e. the Multiscale Universal Interface (MUI), which is capable of facilitating the coupling effort for a wide range of multiscale simulations. The library adopts a header-only form with minimal external dependency and hence can be easily dropped into existing codes. A data sampler concept is introduced, combined with a hybrid dynamic/static typing mechanism, to create an easily customizable framework for solver-independent data interpretation. The library integrates MPI MPMD support and an asynchronous communication protocol to handle inter-solver information exchange irrespective of the solvers' own MPI awareness. Template metaprogramming is heavily employed to simultaneously improve runtime performance and code flexibility. We validated the library by solving three different multiscale problems, which also serve to demonstrate the flexibility of the framework in handling heterogeneous models and solvers. In the first example, a Couette flow was simulated using two concurrently coupled Smoothed Particle Hydrodynamics (SPH) simulations of different spatial resolutions. In the second example, we coupled the deterministic SPH method with the stochastic Dissipative Particle Dynamics (DPD) method to study the effect of surface grafting on the hydrodynamics properties on the surface. In the third example, we consider conjugate heat transfer between a solid domain and a fluid domain by coupling the particle-based energy-conserving DPD (eDPD) method with the Finite Element Method (FEM)
Broader source: Energy.gov [DOE]
Micro-structural modeling tools for metals are being developed and used to demonstrate a design framework to improve the understanding of dynamic response and statistical variability. This project will enable design engineers to evaluate the effects of design changes and material selection; anticipate quality and cost prior to implementation on the factory floor; and enable low-waste, low-cost manufacturing. Third Wave Systems, Inc. - Minneapolis, MN
Understanding Creep Mechanisms in Graphite with Experiments, Multiscale Simulations, and Modeling
Eapen, Jacob; Murty, Korukonda; Burchell, Timothy
2014-06-02
Disordering mechanisms in graphite have a long history with conflicting viewpoints. Using Raman and x-ray photon spectroscopy, electron microscopy, x-ray diffraction experiments and atomistic modeling and simulations, the current project has developed a fundamental understanding of early-to-late state radiation damage mechanisms in nuclear reactor grade graphite (NBG-18 and PCEA). We show that the topological defects in graphite play an important role under neutron and ion irradiation.
SUSTAINABLE MANUFACTURING VIA MULTI-SCALE PHYSICS-BASED PROCESS MODELING
Broader source: Energy.gov (indexed) [DOE]
AND MANUFACTURING-INFORMED DESIGN | Department of Energy Third Wave Systems, Inc. - Minneapolis, MN Micro-structural modeling tools for metals are being developed and used to demonstrate a design framework to improve the understanding of dynamic response and statistical variability. This project will enable design engineers to evaluate the effects of design changes and material selection; anticipate quality and cost prior to implementation on the factory floor; and enable low-waste, low-cost
Multiscale modeling of thermal conductivity of high burnup structures in UO_{2} fuels
Bai, Xian -Ming; Tonks, Michael R.; Zhang, Yongfeng; Hales, Jason D.
2015-12-22
The high burnup structure forming at the rim region in UO_{2} based nuclear fuel pellets has interesting physical properties such as improved thermal conductivity, even though it contains a high density of grain boundaries and micron-size gas bubbles. To understand this counterintuitive phenomenon, mesoscale heat conduction simulations with inputs from atomistic simulations and experiments were conducted to study the thermal conductivities of a small-grain high burnup microstructure and two large-grain unrestructured microstructures. We concluded that the phonon scattering effects caused by small point defects such as dispersed Xe atoms in the grain interior must be included in order to correctly predict the thermal transport properties of these microstructures. In extreme cases, even a small concentration of dispersed Xe atoms such as 10^{-5} can result in a lower thermal conductivity in the large-grain unrestructured microstructures than in the small-grain high burnup structure. The high-density grain boundaries in a high burnup structure act as defect sinks and can reduce the concentration of point defects in its grain interior and improve its thermal conductivity in comparison with its large-grain counterparts. Furthermore, an analytical model was developed to describe the thermal conductivity at different concentrations of dispersed Xe, bubble porosities, and grain sizes. Upon calibration, the model is robust and agrees well with independent heat conduction modeling over a wide range of microstructural parameters.
Multiscale modeling of thermal conductivity of high burnup structures in UO2 fuels
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Bai, Xian -Ming; Tonks, Michael R.; Zhang, Yongfeng; Hales, Jason D.
2015-12-22
The high burnup structure forming at the rim region in UO2 based nuclear fuel pellets has interesting physical properties such as improved thermal conductivity, even though it contains a high density of grain boundaries and micron-size gas bubbles. To understand this counterintuitive phenomenon, mesoscale heat conduction simulations with inputs from atomistic simulations and experiments were conducted to study the thermal conductivities of a small-grain high burnup microstructure and two large-grain unrestructured microstructures. We concluded that the phonon scattering effects caused by small point defects such as dispersed Xe atoms in the grain interior must be included in order to correctlymore » predict the thermal transport properties of these microstructures. In extreme cases, even a small concentration of dispersed Xe atoms such as 10-5 can result in a lower thermal conductivity in the large-grain unrestructured microstructures than in the small-grain high burnup structure. The high-density grain boundaries in a high burnup structure act as defect sinks and can reduce the concentration of point defects in its grain interior and improve its thermal conductivity in comparison with its large-grain counterparts. Furthermore, an analytical model was developed to describe the thermal conductivity at different concentrations of dispersed Xe, bubble porosities, and grain sizes. Upon calibration, the model is robust and agrees well with independent heat conduction modeling over a wide range of microstructural parameters.« less
Swift, D. C.; Paisley, Dennis L.; Kyrala, George A.; Hauer, Allan
2002-01-01
Ab initio quantum mechanics was used to construct a thermodynamically complete and rigorous equation of state for beryllium in the hexagonal and body-centred cubic structures, and to predict elastic constants as a function of compression. The equation of state agreed well with Hugoniot data and previously-published equations of state, but the temperatures were significantly different. The hexagonal/bcc phase boundary agreed reasonably well with published data, suggesting that the temperatures in our new equation of state were accurate. Shock waves were induced in single crystals and polycrystalline foils of beryllium, by direct illumination using the TRIDENT laser at Los Alamos. The velocity history at the surface of the sample was measured using a line-imaging VISAR, and transient X-ray diffraction (TXD) records were obtained with a plasma backlighter and X-ray streak cameras. The VISAR records exhibited elastic precursors, plastic waves, phase changes and spall. Dual TXD records were taken, in Bragg and Laue orientations. The Bragg lines moved in response to compression in the uniaxial direction. Because direct laser drive was used, the results had to be interpreted with the aid of radiation hydrodynamics simulations to predict the loading history for each laser pulse. In the experiments where there was evidence of polymorphism in the VISAR record, additional lines appeared in the Bragg and Laue records. The corresponding pressures were consistent with the phase boundary predicted by the quantum mechanical equation of state for beryllium. A model of the response of a single crystal of beryllium to shock loading is being developed using these new theoretical and experimental results. This model will be used in meso-scale studies of the response of the microstructure, allowing us to develop a more accurate representation of the behaviour of polycrystalline beryllium.
Uncertainty quantification and multiscale mathematics. (Conference...
Office of Scientific and Technical Information (OSTI)
quantification and multiscale mathematics. Citation Details In-Document Search Title: Uncertainty quantification and multiscale mathematics. Authors: Trucano, Timothy Guy ...
Uncertainty quantification and multiscale mathematics. (Conference...
Office of Scientific and Technical Information (OSTI)
Uncertainty quantification and multiscale mathematics. Citation Details In-Document Search Title: Uncertainty quantification and multiscale mathematics. No abstract prepared. ...
Plimpton, Steven J.; Sershen, Cheryl L.; May, Elebeoba E.
2015-01-01
This paper describes a method for incorporating a diffusion field modeling oxygen usage and dispersion in a multi-scale model of Mycobacterium tuberculosis (Mtb) infection mediated granuloma formation. We implemented this method over a floating-point field to model oxygen dynamics in host tissue during chronic phase response and Mtb persistence. The method avoids the requirement of satisfying the Courant-Friedrichs-Lewy (CFL) condition, which is necessary in implementing the explicit version of the finite-difference method, but imposes an impractical bound on the time step. Instead, diffusion is modeled by a matrix-based, steady state approximate solution to the diffusion equation. Moreover, presented in figure 1 is the evolution of the diffusion profiles of a containment granuloma over time.
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Plimpton, Steven J.; Sershen, Cheryl L.; May, Elebeoba E.
2015-01-01
This paper describes a method for incorporating a diffusion field modeling oxygen usage and dispersion in a multi-scale model of Mycobacterium tuberculosis (Mtb) infection mediated granuloma formation. We implemented this method over a floating-point field to model oxygen dynamics in host tissue during chronic phase response and Mtb persistence. The method avoids the requirement of satisfying the Courant-Friedrichs-Lewy (CFL) condition, which is necessary in implementing the explicit version of the finite-difference method, but imposes an impractical bound on the time step. Instead, diffusion is modeled by a matrix-based, steady state approximate solution to the diffusion equation. Moreover, presented in figuremore » 1 is the evolution of the diffusion profiles of a containment granuloma over time.« less
Dynamic Multiscale Averaging (DMA) of Turbulent Flow
Richard W. Johnson
2012-09-01
A new approach called dynamic multiscale averaging (DMA) for computing the effects of turbulent flow is described. The new method encompasses multiple applications of temporal and spatial averaging, that is, multiscale operations. Initially, a direct numerical simulation (DNS) is performed for a relatively short time; it is envisioned that this short time should be long enough to capture several fluctuating time periods of the smallest scales. The flow field variables are subject to running time averaging during the DNS. After the relatively short time, the time-averaged variables are volume averaged onto a coarser grid. Both time and volume averaging of the describing equations generate correlations in the averaged equations. These correlations are computed from the flow field and added as source terms to the computation on the next coarser mesh. They represent coupling between the two adjacent scales. Since they are computed directly from first principles, there is no modeling involved. However, there is approximation involved in the coupling correlations as the flow field has been computed for only a relatively short time. After the time and spatial averaging operations are applied at a given stage, new computations are performed on the next coarser mesh using a larger time step. The process continues until the coarsest scale needed is reached. New correlations are created for each averaging procedure. The number of averaging operations needed is expected to be problem dependent. The new DMA approach is applied to a relatively low Reynolds number flow in a square duct segment. Time-averaged stream-wise velocity and vorticity contours from the DMA approach appear to be very similar to a full DNS for a similar flow reported in the literature. Expected symmetry for the final results is produced for the DMA method. The results obtained indicate that DMA holds significant potential in being able to accurately compute turbulent flow without modeling for practical
Weston, David; Hanson, Paul J; Norby, Richard J; Tuskan, Gerald A; Wullschleger, Stan D
2012-01-01
Network analysis is now a common statistical tool for molecular biologists. Network algorithms are readily used to model gene, protein and metabolic correlations providing insight into pathways driving biological phenomenon. One output from such an analysis is a candidate gene list that can be responsible, in part, for the biological process of interest. The question remains, however, as to whether molecular network analysis can be used to inform process models at higher levels of biological organization. In our previous work, transcriptional networks derived from three plant species were constructed, interrogated for orthology and then correlated to photosynthetic inhibition at elevated temperature. One unique aspect of that study was the link from co-expression networks to net photosynthesis. In this addendum, we propose a conceptual model where traditional network analysis can be linked to whole-plant models thereby informing predictions on key processes such as photosynthesis, nutrient uptake and assimilation, and C partitioning.
A Hybrid Multiscale Framework for Subsurface Flow and Transport Simulations
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Scheibe, Timothy D.; Yang, Xiaofan; Chen, Xingyuan; Hammond, Glenn E.
2015-06-01
Extensive research efforts have been invested in reducing model errors to improve the predictive ability of biogeochemical earth and environmental system simulators, with applications ranging from contaminant transport and remediation to impacts of biogeochemical elemental cycling (e.g., carbon and nitrogen) on local ecosystems and regional to global climate. While the bulk of this research has focused on improving model parameterizations in the face of observational limitations, the more challenging type of model error/uncertainty to identify and quantify is model structural error which arises from incorrect mathematical representations of (or failure to consider) important physical, chemical, or biological processes, properties, ormore » system states in model formulations. While improved process understanding can be achieved through scientific study, such understanding is usually developed at small scales. Process-based numerical models are typically designed for a particular characteristic length and time scale. For application-relevant scales, it is generally necessary to introduce approximations and empirical parameterizations to describe complex systems or processes. This single-scale approach has been the best available to date because of limited understanding of process coupling combined with practical limitations on system characterization and computation. While computational power is increasing significantly and our understanding of biological and environmental processes at fundamental scales is accelerating, using this information to advance our knowledge of the larger system behavior requires the development of multiscale simulators. Accordingly there has been much recent interest in novel multiscale methods in which microscale and macroscale models are explicitly coupled in a single hybrid multiscale simulation. A limited number of hybrid multiscale simulations have been developed for biogeochemical earth systems, but they mostly utilize application
A Hybrid Multiscale Framework for Subsurface Flow and Transport Simulations
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Scheibe, Timothy D.; Yang, Xiaofan; Chen, Xingyuan; Hammond, Glenn E.
2015-06-01
Extensive research efforts have been invested in reducing model errors to improve the predictive ability of biogeochemical earth and environmental system simulators, with applications ranging from contaminant transport and remediation to impacts of biogeochemical elemental cycling (e.g., carbon and nitrogen) on local ecosystems and regional to global climate. While the bulk of this research has focused on improving model parameterizations in the face of observational limitations, the more challenging type of model error/uncertainty to identify and quantify is model structural error which arises from incorrect mathematical representations of (or failure to consider) important physical, chemical, or biological processes, properties, ormoresystem states in model formulations. While improved process understanding can be achieved through scientific study, such understanding is usually developed at small scales. Process-based numerical models are typically designed for a particular characteristic length and time scale. For application-relevant scales, it is generally necessary to introduce approximations and empirical parameterizations to describe complex systems or processes. This single-scale approach has been the best available to date because of limited understanding of process coupling combined with practical limitations on system characterization and computation. While computational power is increasing significantly and our understanding of biological and environmental processes at fundamental scales is accelerating, using this information to advance our knowledge of the larger system behavior requires the development of multiscale simulators. Accordingly there has been much recent interest in novel multiscale methods in which microscale and macroscale models are explicitly coupled in a single hybrid multiscale simulation. A limited number of hybrid multiscale simulations have been developed for biogeochemical earth systems, but they mostly utilize application
A Hybrid Multiscale Framework for Subsurface Flow and Transport Simulations
Scheibe, Timothy D.; Yang, Xiaofan; Chen, Xingyuan; Hammond, Glenn E.
2015-06-01
Extensive research efforts have been invested in reducing model errors to improve the predictive ability of biogeochemical earth and environmental system simulators, with applications ranging from contaminant transport and remediation to impacts of biogeochemical elemental cycling (e.g., carbon and nitrogen) on local ecosystems and regional to global climate. While the bulk of this research has focused on improving model parameterizations in the face of observational limitations, the more challenging type of model error/uncertainty to identify and quantify is model structural error which arises from incorrect mathematical representations of (or failure to consider) important physical, chemical, or biological processes, properties, or system states in model formulations. While improved process understanding can be achieved through scientific study, such understanding is usually developed at small scales. Process-based numerical models are typically designed for a particular characteristic length and time scale. For application-relevant scales, it is generally necessary to introduce approximations and empirical parameterizations to describe complex systems or processes. This single-scale approach has been the best available to date because of limited understanding of process coupling combined with practical limitations on system characterization and computation. While computational power is increasing significantly and our understanding of biological and environmental processes at fundamental scales is accelerating, using this information to advance our knowledge of the larger system behavior requires the development of multiscale simulators. Accordingly there has been much recent interest in novel multiscale methods in which microscale and macroscale models are explicitly coupled in a single hybrid multiscale simulation. A limited number of hybrid multiscale simulations have been developed for biogeochemical earth systems, but they mostly utilize application
Graph modeling systems and methods
Neergaard, Mike
2015-10-13
An apparatus and a method for vulnerability and reliability modeling are provided. The method generally includes constructing a graph model of a physical network using a computer, the graph model including a plurality of terminating vertices to represent nodes in the physical network, a plurality of edges to represent transmission paths in the physical network, and a non-terminating vertex to represent a non-nodal vulnerability along a transmission path in the physical network. The method additionally includes evaluating the vulnerability and reliability of the physical network using the constructed graph model, wherein the vulnerability and reliability evaluation includes a determination of whether each terminating and non-terminating vertex represents a critical point of failure. The method can be utilized to evaluate wide variety of networks, including power grid infrastructures, communication network topologies, and fluid distribution systems.
Niyogi, Devdutta S.
2013-06-07
The CLASIC experiment was conducted over the US southern great plains (SGP) in June 2007 with an objective to lead an enhanced understanding of the cumulus convection particularly as it relates to land surface conditions. This project was design to help assist with understanding the overall improvement of land atmosphere convection initiation representation of which is important for global and regional models. The study helped address one of the critical documented deficiency in the models central to the ARM objectives for cumulus convection initiation and particularly under summer time conditions. This project was guided by the scientific question building on the CLASIC theme questions: What is the effect of improved land surface representation on the ability of coupled models to simulate cumulus and convection initiation? The focus was on the US Southern Great Plains region. Since the CLASIC period was anomalously wet the strategy has been to use other periods and domains to develop the comparative assessment for the CLASIC data period, and to understand the mechanisms of the anomalous wet conditions on the tropical systems and convection over land. The data periods include the IHOP 2002 field experiment that was over roughly same domain as the CLASIC in the SGP, and some of the DOE funded Ameriflux datasets.
Multiscale Molecular Simulations at the Petascale (Parallelization of
Office of Scientific and Technical Information (OSTI)
Reactive Force Field Model for Blue Gene/Q): ALCF-2 Early Science Program Technical Report (Technical Report) | SciTech Connect Multiscale Molecular Simulations at the Petascale (Parallelization of Reactive Force Field Model for Blue Gene/Q): ALCF-2 Early Science Program Technical Report Citation Details In-Document Search Title: Multiscale Molecular Simulations at the Petascale (Parallelization of Reactive Force Field Model for Blue Gene/Q): ALCF-2 Early Science Program Technical Report
Yu, Sungduk; Pritchard, Michael S.
2015-12-17
The effect of global climate model (GCM) time step—which also controls how frequently global and embedded cloud resolving scales are coupled—is examined in the Superparameterized Community Atmosphere Model ver 3.0. Systematic bias reductions of time-mean shortwave cloud forcing (~10 W/m^{2}) and longwave cloud forcing (~5 W/m^{2}) occur as scale coupling frequency increases, but with systematically increasing rainfall variance and extremes throughout the tropics. An overarching change in the vertical structure of deep tropical convection, favoring more bottom-heavy deep convection as a global model time step is reduced may help orchestrate these responses. The weak temperature gradient approximation is more faithfully satisfied when a high scale coupling frequency (a short global model time step) is used. These findings are distinct from the global model time step sensitivities of conventionally parameterized GCMs and have implications for understanding emergent behaviors of multiscale deep convective organization in superparameterized GCMs. Lastly, the results may also be useful for helping to tune them.
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Yu, Sungduk; Pritchard, Michael S.
2015-12-17
The effect of global climate model (GCM) time step—which also controls how frequently global and embedded cloud resolving scales are coupled—is examined in the Superparameterized Community Atmosphere Model ver 3.0. Systematic bias reductions of time-mean shortwave cloud forcing (~10 W/m2) and longwave cloud forcing (~5 W/m2) occur as scale coupling frequency increases, but with systematically increasing rainfall variance and extremes throughout the tropics. An overarching change in the vertical structure of deep tropical convection, favoring more bottom-heavy deep convection as a global model time step is reduced may help orchestrate these responses. The weak temperature gradient approximation is more faithfullymore » satisfied when a high scale coupling frequency (a short global model time step) is used. These findings are distinct from the global model time step sensitivities of conventionally parameterized GCMs and have implications for understanding emergent behaviors of multiscale deep convective organization in superparameterized GCMs. Lastly, the results may also be useful for helping to tune them.« less
Parallel Paradigm for Ultraparallel Multi-Scale Brain Blood Flow
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Simulations | Argonne Leadership Computing Facility Parallel Paradigm for Ultraparallel Multi-Scale Brain Blood Flow Simulations Authors: Grinberg, L., Karniadakis, G.E. In this paper we present one approach in building a scalable solver NekTarG for solution of multi-scale and large size problems [1]. NekTarG has been designed for multi-scale blood modeling. The macro-vascular scales describing the flow dynamics in large vessels are coupled to the mesovascular scales unfolding dynamics of
Multi-Scale Simulations Solve a Plasma Turbulence Mystery
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Multi-Scale Simulations Solve a Plasma Turbulence Mystery Multi-Scale Simulations Solve a Plasma Turbulence Mystery Coupled Model Reproduces Experimental Electron Heat Losses March 7, 2016 Contact: Kathy Kincade, kkincade@lbl.gov, +1 510 495 2124 turb cross High-res image of the inside of the Alcator C-Mod tokamak, with a representative cross-section of a plasma. The inset shows the approximate domain for one of the multi-scale simulations and a graphic of the plasma turbulence in the
Multiscale reactive molecular dynamics | Argonne Leadership Computing
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Facility reactive molecular dynamics Authors: Chris KnighT, Gerrick E. Lindberg, Gregory A. Voth Many processes important to chemistry, materials science, and biology cannot be described without considering electronic and nuclear-level dynamics and their coupling to slower, cooperative motions of the system. These inherently multiscale problems require computationally efficient and accurate methods to converge statistical properties. In this paper, a method is presented that uses data
Andrade, JosÃÂ© E; Rudnicki, John W
2012-12-14
In this project, a predictive multiscale framework will be developed to simulate the strong coupling between solid deformations and fluid diffusion in porous rocks. We intend to improve macroscale modeling by incorporating fundamental physical modeling at the microscale in a computationally efficient way. This is an essential step toward further developments in multiphysics modeling, linking hydraulic, thermal, chemical, and geomechanical processes. This research will focus on areas where severe deformations are observed, such as deformation bands, where classical phenomenology breaks down. Multiscale geometric complexities and key geomechanical and hydraulic attributes of deformation bands (e.g., grain sliding and crushing, and pore collapse, causing interstitial fluid expulsion under saturated conditions), can significantly affect the constitutive response of the skeleton and the intrinsic permeability. Discrete mechanics (DEM) and the lattice Boltzmann method (LBM) will be used to probe the microstructure---under the current state---to extract the evolution of macroscopic constitutive parameters and the permeability tensor. These evolving macroscopic constitutive parameters are then directly used in continuum scale predictions using the finite element method (FEM) accounting for the coupled solid deformation and fluid diffusion. A particularly valuable aspect of this research is the thorough quantitative verification and validation program at different scales. The multiscale homogenization framework will be validated using X-ray computed tomography and 3D digital image correlation in situ at the Advanced Photon Source in Argonne National Laboratories. Also, the hierarchical computations at the specimen level will be validated using the aforementioned techniques in samples of sandstone undergoing deformation bands.
A Multiscale Bidirectional Coupling Framework
Kabilan, Senthil; Kuprat, Andrew P.; Hlastala, Michael P.; Corley, Richard A.; Einstein, Daniel R.
2011-12-01
The lung is geometrically articulated across multiple scales from the trachea to the alveoli. A major computational challenge is to tightly link ODEs that describe lower scales to 3D finite element or finite volume models of airway mechanics using iterative communication between scales. In this study, we developed a novel multiscale computational framework for bidirectionally coupling 3D CFD models and systems of lower order ODEs. To validate the coupling framework, a four and eight generation Weibel lung model was constructed. For the coupled CFD-ODE simulations, the lung models were truncated at different generations and a RL circuit represented the truncated portion. The flow characteristics from the coupled models were compared to untruncated full 3D CFD models at peak inhalation and peak exhalation. Results showed that at no time or simulation was the difference in mass flux and/or pressure at a given location between uncoupled and coupled models was greater than 2.43%. The flow characteristics at prime locations for the coupled models showed good agreement to uncoupled models. Remarkably, due to reuse of the Krylov subspace, the cost of the ODE coupling is not much greater than uncoupled full 3D-CFD computations with simple prescribed pressure values at the outlets.
Costigan, Keeley Rochelle; Dubey, Manvendra Krishna
2015-07-10
Atmospheric models are compared in collaboration with LANL and the University of Michigan to understand emissions and the condition of the atmosphere from a model perspective.
MULTISCALE MATHEMATICS FOR BIOMASS CONVERSION TO RENEWABLE HYDROGEN
Vlachos, Dionisios; Plechac, Petr; Katsoulakis, Markos
2013-09-05
The overall objective of this project is to develop multiscale models for understanding and eventually designing complex processes for renewables. To the best of our knowledge, our work is the first attempt at modeling complex reacting systems, whose performance relies on underlying multiscale mathematics. Our specific application lies at the heart of biofuels initiatives of DOE and entails modeling of catalytic systems, to enable economic, environmentally benign, and efficient conversion of biomass into either hydrogen or valuable chemicals. Specific goals include: (i) Development of rigorous spatio-temporal coarse-grained kinetic Monte Carlo (KMC) mathematics and simulation for microscopic processes encountered in biomass transformation. (ii) Development of hybrid multiscale simulation that links stochastic simulation to a deterministic partial differential equation (PDE) model for an entire reactor. (iii) Development of hybrid multiscale simulation that links KMC simulation with quantum density functional theory (DFT) calculations. (iv) Development of parallelization of models of (i)-(iii) to take advantage of Petaflop computing and enable real world applications of complex, multiscale models. In this NCE period, we continued addressing these objectives and completed the proposed work. Main initiatives, key results, and activities are outlined.
Kim, G. H.; Smith, K.
2009-05-01
Addresses battery requirements for electric vehicles using a model that evaluates physical-chemical processes in lithium-ion batteries, from atomic variations to vehicle interface controls.
Multi-Scale Multi-Dimensional Li-Ion Battery Model for Better Design and Management (Presentation)
Kim, G.-H.; Smith, K.
2008-10-01
The developed model used is to provide a better understanding and help answer engineering questions about improving the design, operational strategy, management, and safety of cells.
Multilingual interfaces for parallel coupling in multiphysics and multiscale systems.
Ong, E. T.; Larson, J. W.; Norris, B.; Jacob, R. L.; Tobis, M.; Steder, M.; Mathematics and Computer Science; Univ. of Wisconsin; Australian National Univ.; Univ. of Chicago
2007-01-01
Multiphysics and multiscale simulation systems are emerging as a new grand challenge in computational science, largely because of increased computing power provided by the distributed-memory parallel programming model on commodity clusters. These systems often present a parallel coupling problem in their intercomponent data exchanges. Another potential problem in these coupled systems is language interoperability between their various constituent codes. In anticipation of combined parallel coupling/language interoperability challenges, we have created a set of interlanguage bindings for a successful parallel coupling library, the Model Coupling Toolkit. We describe the method used for automatically generating the bindings using the Babel language interoperability tool, and illustrate with short examples how MCT can be used from the C++ and Python languages. We report preliminary performance reports for the MCT interpolation benchmark. We conclude with a discussion of the significance of this work to the rapid prototyping of large parallel coupled systems.
Freed, Alan D.; Einstein, Daniel R.
2011-04-14
An isotropic constitutive model for the parenchyma of lung has been derived from the theory of hypo-elasticity. The intent is to use it to represent the mechanical response of this soft tissue in sophisticated, computational, fluid-dynamic models of the lung. This demands that the continuum model be accurate, yet simple and effcient. An objective algorithm for its numeric integration is provided. The response of the model is determined for several boundary-value problems whose experiments are used for material characterization. The effective elastic, bulk, and shear moduli, and Poissons ratio, as tangent functions, are also derived. The model is characterized against published experimental data for lung. A bridge between this continuum model and a dodecahedral model of alveolar geometry is investigated, with preliminary findings being reported.
Method and apparatus for modeling interactions
Xavier, Patrick G.
2002-01-01
The present invention provides a method and apparatus for modeling interactions that overcomes drawbacks. The method of the present invention comprises representing two bodies undergoing translations by two swept volume representations. Interactions such as nearest approach and collision can be modeled based on the swept body representations. The present invention is more robust and allows faster modeling than previous methods.
Multi-scale First-Principles Modeling of Three-Phase System of Polymer Electrolyte Membrane Fuel Cel
Brunello, Giuseppe; Choi, Ji; Harvey, David; Jang, Seung
2012-07-01
The three-phase system consisting of Nafion, graphite and platinum in the presence of water is studied using molecule dynamics simulation. The force fields describing the molecular interaction between the components in the system are developed to reproduce the energies calculated from density functional theory modeling. The configuration of such complicated three-phase system is predicted through MD simulations. The nanophase-segregation and transport properties are investigated from the equilibrium state. The coverage of the electrolyte on the platinum surface and the dissolution of oxygen are analyzed.
Zhang, Xuesong; Sahajpal, Ritvik; Manowitz, D.; Zhao, Kaiguang; LeDuc, Stephen D.; Xu, Min; Xiong, Wei; Zhang, Aiping; Izaurralde, Roberto C.; Thomson, Allison M.; West, Tristram O.; Post, W. M.
2014-05-01
The development of effective measures to stabilize atmospheric CO2 concentration and mitigate negative impacts of climate change requires accurate quantification of the spatial variation and magnitude of the terrestrial carbon (C) flux. However, the spatial pattern and strength of terrestrial C sinks and sources remain uncertain. In this study, we designed a spatially-explicit agroecosystem modeling system by integrating the Environmental Policy Integrated Climate (EPIC) model with multiple sources of geospatial and surveyed datasets (including crop type map, elevation, climate forcing, fertilizer application, tillage type and distribution, and crop planting and harvesting date), and applied it to examine the sensitivity of cropland C flux simulations to two widely used soil databases (i.e. State Soil Geographic-STATSGO of a scale of 1:250,000 and Soil Survey Geographic-SSURGO of a scale of 1:24,000) in Iowa, USA. To efficiently execute numerous EPIC runs resulting from the use of high resolution spatial data (56m), we developed a parallelized version of EPIC. Both STATSGO and SSURGO led to similar simulations of crop yields and Net Ecosystem Production (NEP) estimates at the State level. However, substantial differences were observed at the county and sub-county (grid) levels. In general, the fine resolution SSURGO data outperformed the coarse resolution STATSGO data for county-scale crop-yield simulation, and within STATSGO, the area-weighted approach provided more accurate results. Further analysis showed that spatial distribution and magnitude of simulated NEP were more sensitive to the resolution difference between SSURGO and STATSGO at the county or grid scale. For over 60% of the cropland areas in Iowa, the deviations between STATSGO- and SSURGO-derived NEP were larger than 1MgCha(-1)yr(-1), or about half of the average cropland NEP, highlighting the significant uncertainty in spatial distribution and magnitude of simulated C fluxes resulting from
A Many-Task Parallel Approach for Multiscale Simulations of Subsurface Flow and Reactive Transport
Scheibe, Timothy D.; Yang, Xiaofan; Schuchardt, Karen L.; Agarwal, Khushbu; Chase, Jared M.; Palmer, Bruce J.; Tartakovsky, Alexandre M.
2014-12-16
Continuum-scale models have long been used to study subsurface flow, transport, and reactions but lack the ability to resolve processes that are governed by pore-scale mixing. Recently, pore-scale models, which explicitly resolve individual pores and soil grains, have been developed to more accurately model pore-scale phenomena, particularly reaction processes that are controlled by local mixing. However, pore-scale models are prohibitively expensive for modeling application-scale domains. This motivates the use of a hybrid multiscale approach in which continuum- and pore-scale codes are coupled either hierarchically or concurrently within an overall simulation domain (time and space). This approach is naturally suited to an adaptive, loosely-coupled many-task methodology with three potential levels of concurrency. Each individual code (pore- and continuum-scale) can be implemented in parallel; multiple semi-independent instances of the pore-scale code are required at each time step providing a second level of concurrency; and Monte Carlo simulations of the overall system to represent uncertainty in material property distributions provide a third level of concurrency. We have developed a hybrid multiscale model of a mixing-controlled reaction in a porous medium wherein the reaction occurs only over a limited portion of the domain. Loose, minimally-invasive coupling of pre-existing parallel continuum- and pore-scale codes has been accomplished by an adaptive script-based workflow implemented in the Swift workflow system. We describe here the methods used to create the model system, adaptively control multiple coupled instances of pore- and continuum-scale simulations, and maximize the scalability of the overall system. We present results of numerical experiments conducted on NERSC supercomputing systems; our results demonstrate that loose many-task coupling provides a scalable solution for multiscale subsurface simulations with minimal overhead.
Multiscale simulations of blood-flow: from a platelet to an artery |
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Argonne Leadership Computing Facility simulations of blood-flow: from a platelet to an artery Authors: L. Grinberg, M. Deng, H. Lei, J. A. Insley, G. E. Karniadakis We review our recent advances on multiscale modeling of blood flow including blood rheology. We focus on the objectives, methods, computational complexity and overall methodology for simulations at the level of glycocalyx (< 1 μm), blood cells (2--8μm) and up to larger arteries (O(cm)). The main findings of our research and
Hybrid multiscale simulation of a mixing-controlled reaction
Scheibe, Timothy D.; Schuchardt, Karen L.; Agarwal, Khushbu; Chase, Jared M.; Yang, Xiaofan; Palmer, Bruce J.; Tartakovsky, Alexandre M.; Elsethagen, Todd O.; Redden, George D.
2015-09-01
Continuum-scale models have been used to study subsurface flow, transport, and reactions for many years but lack the capability to resolve fine-grained processes. Recently, pore-scale models, which operate at scales of individual soil grains, have been developed to more accurately model and study pore-scale phenomena, such as mineral precipitation and dissolution reactions, microbially-mediated surface reactions, and other complex processes. However, these highly-resolved models are prohibitively expensive for modeling domains of sizes relevant to practical problems. To broaden the utility of pore-scale models for larger domains, we developed a hybrid multiscale model that initially simulates the full domain at the continuum scale and applies a pore-scale model only to areas of high reactivity. Since the location and number of pore-scale model regions in the model varies as the reactions proceed, an adaptive script defines the number and location of pore regions within each continuum iteration and initializes pore-scale simulations from macroscale information. Another script communicates information from the pore-scale simulation results back to the continuum scale. These components provide loose coupling between the pore- and continuum-scale codes into a single hybrid multiscale model implemented within the SWIFT workflow environment. In this paper, we consider an irreversible homogenous bimolecular reaction (two solutes reacting to form a third solute) in a 2D test problem. This paper is focused on the approach used for multiscale coupling between pore- and continuum-scale models, application to a realistic test problem, and implications of the results for predictive simulation of mixing-controlled reactions in porous media. Our results and analysis demonstrate that loose coupling provides a feasible, efficient and scalable approach for multiscale subsurface simulations.
Spectral method for a kinetic swarming model
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Gamba, Irene M.; Haack, Jeffrey R.; Motsch, Sebastien
2015-04-28
Here we present the first numerical method for a kinetic description of the Vicsek swarming model. The kinetic model poses a unique challenge, as there is a distribution dependent collision invariant to satisfy when computing the interaction term. We use a spectral representation linked with a discrete constrained optimization to compute these interactions. To test the numerical scheme we investigate the kinetic model at different scales and compare the solution with the microscopic and macroscopic descriptions of the Vicsek model. Lastly, we observe that the kinetic model captures key features such as vortex formation and traveling waves.
Other: Multiscale Simulation of Blood Flow in Brain Arteries...
Office of Scientific and Technical Information (OSTI)
Multiscale Simulation of Blood Flow in Brain Arteries with an Aneurysm Citation Details Title: Multiscale Simulation of Blood Flow in Brain Arteries with an Aneurysm
Multiscale Blood Flow Simulations | Argonne Leadership Computing Facility
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
a continuum model Brain blood flow simulation with NekTar; a continuum model. Leopold Grinberg, George Em. Karniadakis, Brown University; Vitali Morozov, Joseph A. Insley, Michael E. Papka, Kalyan Kumaran; Argonne National Laboratory; Dmitry A Fedosov, Forschungszentrum Juelich Multiscale Blood Flow Simulations PI Name: George Karniadakis PI Email: gk@dam.brown.edu Institution: Brown University Allocation Program: INCITE Allocation Hours at ALCF: 71 Million Year: 2013 Research Domain: Biological
"Multiscale Capabilities for Exploring Transport Phenomena in...
Office of Scientific and Technical Information (OSTI)
in Batteries": Ab Initio Calculations on Defective LiFePO4 Citation Details In-Document Search Title: "Multiscale Capabilities for Exploring Transport Phenomena in Batteries": Ab ...
Multiscale Simulations of Human Pathologies | Argonne Leadership...
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Inset shows the time evolution of thrombus formation. George Karniadakis, Brown University Multiscale Simulations of Human Pathologies PI Name: George Karniadakis PI Email: ...
Multiscale Simulations of Human Pathologies | Argonne Leadership...
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Multiscale Simulations of Human Pathologies PI Name: George Karniadakis PI Email: ... of Sickle Hemoglobin Predicting Human Blood Viscosity in Silico Tightly Coupled ...
Multiscale Simulations of Human Pathologies | Argonne Leadership...
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Karniadakis, Paris Perdikaris, and Yue Yu, Brown University; Leopold Grinberg, IBM T. J. Watson Research Center and Brown University Multiscale Simulations of Human Pathologies PI ...
Computational Physics and Methods
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
2 Computational Physics and Methods Performing innovative simulations of physics phenomena on tomorrow's scientific computing platforms Growth and emissivity of young galaxy hosting a supermassive black hole as calculated in cosmological code ENZO and post-processed with radiative transfer code AURORA. image showing detailed turbulence simulation, Rayleigh-Taylor Turbulence imaging: the largest turbulence simulations to date Advanced multi-scale modeling Turbulence datasets Density iso-surfaces
Davtyan, Aram; Dama, James F.; Voth, Gregory A.; Andersen, Hans C.
2015-04-21
Coarse-grained (CG) models of molecular systems, with fewer mechanical degrees of freedom than an all-atom model, are used extensively in chemical physics. It is generally accepted that a coarse-grained model that accurately describes equilibrium structural properties (as a result of having a well constructed CG potential energy function) does not necessarily exhibit appropriate dynamical behavior when simulated using conservative Hamiltonian dynamics for the CG degrees of freedom on the CG potential energy surface. Attempts to develop accurate CG dynamic models usually focus on replacing Hamiltonian motion by stochastic but Markovian dynamics on that surface, such as Langevin or Brownian dynamics. However, depending on the nature of the system and the extent of the coarse-graining, a Markovian dynamics for the CG degrees of freedom may not be appropriate. In this paper, we consider the problem of constructing dynamic CG models within the context of the Multi-Scale Coarse-graining (MS-CG) method of Voth and coworkers. We propose a method of converting a MS-CG model into a dynamic CG model by adding degrees of freedom to it in the form of a small number of fictitious particles that interact with the CG degrees of freedom in simple ways and that are subject to Langevin forces. The dynamic models are members of a class of nonlinear systems interacting with special heat baths that were studied by Zwanzig [J. Stat. Phys. 9, 215 (1973)]. The properties of the fictitious particles can be inferred from analysis of the dynamics of all-atom simulations of the system of interest. This is analogous to the fact that the MS-CG method generates the CG potential from analysis of equilibrium structures observed in all-atom simulation data. The dynamic models generate a non-Markovian dynamics for the CG degrees of freedom, but they can be easily simulated using standard molecular dynamics programs. We present tests of this method on a series of simple examples that demonstrate that
Multiscale characterization and analysis of shapes
Prasad, Lakshman; Rao, Ramana
2002-01-01
An adaptive multiscale method approximates shapes with continuous or uniformly and densely sampled contours, with the purpose of sparsely and nonuniformly discretizing the boundaries of shapes at any prescribed resolution, while at the same time retaining the salient shape features at that resolution. In another aspect, a fundamental geometric filtering scheme using the Constrained Delaunay Triangulation (CDT) of polygonized shapes creates an efficient parsing of shapes into components that have semantic significance dependent only on the shapes' structure and not on their representations per se. A shape skeletonization process generalizes to sparsely discretized shapes, with the additional benefit of prunability to filter out irrelevant and morphologically insignificant features. The skeletal representation of characters of varying thickness and the elimination of insignificant and noisy spurs and branches from the skeleton greatly increases the robustness, reliability and recognition rates of character recognition algorithms.
The Adaptive Multi-scale Simulation Infrastructure
Tobin, William R.
2015-09-01
The Adaptive Multi-scale Simulation Infrastructure (AMSI) is a set of libraries and tools developed to support the development, implementation, and execution of general multimodel simulations. Using a minimal set of simulation meta-data AMSI allows for minimally intrusive work to adapt existent single-scale simulations for use in multi-scale simulations. Support for dynamic runtime operations such as single- and multi-scale adaptive properties is a key focus of AMSI. Particular focus has been spent on the development on scale-sensitive load balancing operations to allow single-scale simulations incorporated into a multi-scale simulation using AMSI to use standard load-balancing operations without affecting the integrity of the overall multi-scale simulation.
Method and apparatus for modeling interactions
Xavier, Patrick G.
2000-08-08
A method and apparatus for modeling interactions between bodies. The method comprises representing two bodies undergoing translations and rotations by two hierarchical swept volume representations. Interactions such as nearest approach and collision can be modeled based on the swept body representations. The present invention can serve as a practical tool in motion planning, CAD systems, simulation systems, safety analysis, and applications that require modeling time-based interactions. A body can be represented in the present invention by a union of convex polygons and convex polyhedra. As used generally herein, polyhedron includes polygon, and polyhedra includes polygons. The body undergoing translation can be represented by a swept body representation, where the swept body representation comprises a hierarchical bounding volume representation whose leaves each contain a representation of the region swept by a section of the body during the translation, and where the union of the regions is a superset of the region swept by the surface of the body during translation. Interactions between two bodies thus represented can be modeled by modeling interactions between the convex hulls of the finite sets of discrete points in the swept body representations.
Multiscale MonteCarlo equilibration: Pure Yang-Mills theory
Endres, Michael G.; Brower, Richard C.; Orginos, Kostas; Detmold, William; Pochinsky, Andrew V.
2015-12-29
In this study, we present a multiscale thermalization algorithm for lattice gauge theory, which enables efficient parallel generation of uncorrelated gauge field configurations. The algorithm combines standard Monte Carlo techniques with ideas drawn from real space renormalization group and multigrid methods. We demonstrate the viability of the algorithm for pure Yang-Mills gauge theory for both heat bath and hybrid Monte Carlo evolution, and show that it ameliorates the problem of topological freezing up to controllable lattice spacing artifacts.
MULTISCALE DYNAMICS OF SOLAR MAGNETIC STRUCTURES
Uritsky, Vadim M.; Davila, Joseph M.
2012-03-20
Multiscale topological complexity of the solar magnetic field is among the primary factors controlling energy release in the corona, including associated processes in the photospheric and chromospheric boundaries. We present a new approach for analyzing multiscale behavior of the photospheric magnetic flux underlying these dynamics as depicted by a sequence of high-resolution solar magnetograms. The approach involves two basic processing steps: (1) identification of timing and location of magnetic flux origin and demise events (as defined by DeForest et al.) by tracking spatiotemporal evolution of unipolar and bipolar photospheric regions, and (2) analysis of collective behavior of the detected magnetic events using a generalized version of the Grassberger-Procaccia correlation integral algorithm. The scale-free nature of the developed algorithms makes it possible to characterize the dynamics of the photospheric network across a wide range of distances and relaxation times. Three types of photospheric conditions are considered to test the method: a quiet photosphere, a solar active region (NOAA 10365) in a quiescent non-flaring state, and the same active region during a period of M-class flares. The results obtained show (1) the presence of a topologically complex asymmetrically fragmented magnetic network in the quiet photosphere driven by meso- and supergranulation, (2) the formation of non-potential magnetic structures with complex polarity separation lines inside the active region, and (3) statistical signatures of canceling bipolar magnetic structures coinciding with flaring activity in the active region. Each of these effects can represent an unstable magnetic configuration acting as an energy source for coronal dissipation and heating.
Energy Science and Technology Software Center (OSTI)
2009-08-01
The code to be released is a new addition to the LAMMPS molecular dynamics code. LAMMPS is developed and maintained by Sandia, is publicly available, and is used widely by both natioanl laboratories and academics. The new addition to be released enables LAMMPS to perform molecular dynamics simulations of shock waves using the Multi-scale Shock Simulation Technique (MSST) which we have developed and has been previously published. This technique enables molecular dynamics simulations of shockmore » waves in materials for orders of magnitude longer timescales than the direct, commonly employed approach.« less
Adaptive model training system and method
Bickford, Randall L; Palnitkar, Rahul M
2014-11-18
An adaptive model training system and method for filtering asset operating data values acquired from a monitored asset for selectively choosing asset operating data values that meet at least one predefined criterion of good data quality while rejecting asset operating data values that fail to meet at least the one predefined criterion of good data quality; and recalibrating a previously trained or calibrated model having a learned scope of normal operation of the asset by utilizing the asset operating data values that meet at least the one predefined criterion of good data quality for adjusting the learned scope of normal operation of the asset for defining a recalibrated model having the adjusted learned scope of normal operation of the asset.
Adaptive model training system and method
Bickford, Randall L; Palnitkar, Rahul M; Lee, Vo
2014-04-15
An adaptive model training system and method for filtering asset operating data values acquired from a monitored asset for selectively choosing asset operating data values that meet at least one predefined criterion of good data quality while rejecting asset operating data values that fail to meet at least the one predefined criterion of good data quality; and recalibrating a previously trained or calibrated model having a learned scope of normal operation of the asset by utilizing the asset operating data values that meet at least the one predefined criterion of good data quality for adjusting the learned scope of normal operation of the asset for defining a recalibrated model having the adjusted learned scope of normal operation of the asset.
Final Technical Report "Multiscale Simulation Algorithms for Biochemical Systems"
Petzold, Linda R.
2012-10-25
Biochemical systems are inherently multiscale and stochastic. In microscopic systems formed by living cells, the small numbers of reactant molecules can result in dynamical behavior that is discrete and stochastic rather than continuous and deterministic. An analysis tool that respects these dynamical characteristics is the stochastic simulation algorithm (SSA, Gillespie, 1976), a numerical simulation procedure that is essentially exact for chemical systems that are spatially homogeneous or well stirred. Despite recent improvements, as a procedure that simulates every reaction event, the SSA is necessarily inefficient for most realistic problems. There are two main reasons for this, both arising from the multiscale nature of the underlying problem: (1) stiffness, i.e. the presence of multiple timescales, the fastest of which are stable; and (2) the need to include in the simulation both species that are present in relatively small quantities and should be modeled by a discrete stochastic process, and species that are present in larger quantities and are more efficiently modeled by a deterministic differential equation (or at some scale in between). This project has focused on the development of fast and adaptive algorithms, and the fun- damental theory upon which they must be based, for the multiscale simulation of biochemical systems. Areas addressed by this project include: (1) Theoretical and practical foundations for ac- celerated discrete stochastic simulation (tau-leaping); (2) Dealing with stiffness (fast reactions) in an efficient and well-justified manner in discrete stochastic simulation; (3) Development of adaptive multiscale algorithms for spatially homogeneous discrete stochastic simulation; (4) Development of high-performance SSA algorithms.
Method of and apparatus for modeling interactions
Budge, Kent G.
2004-01-13
A method and apparatus for modeling interactions can accurately model tribological and other properties and accommodate topological disruptions. Two portions of a problem space are represented, a first with a Lagrangian mesh and a second with an ALE mesh. The ALE and Lagrangian meshes are constructed so that each node on the surface of the Lagrangian mesh is in a known correspondence with adjacent nodes in the ALE mesh. The interaction can be predicted for a time interval. Material flow within the ALE mesh can accurately model complex interactions such as bifurcation. After prediction, nodes in the ALE mesh in correspondence with nodes on the surface of the Lagrangian mesh can be mapped so that they are once again adjacent to their corresponding Lagrangian mesh nodes. The ALE mesh can then be smoothed to reduce mesh distortion that might reduce the accuracy or efficiency of subsequent prediction steps. The process, from prediction through mapping and smoothing, can be repeated until a terminal condition is reached.
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Lim, H.; Hale, L. M.; Zimmerman, J. A.; Battaile, C. C.; Weinberger, C. R.
2015-01-05
In this study, we develop an atomistically informed crystal plasticity finite element (CP-FE) model for body-centered-cubic (BCC) α-Fe that incorporates non-Schmid stress dependent slip with temperature and strain rate effects. Based on recent insights obtained from atomistic simulations, we propose a new constitutive model that combines a generalized non-Schmid yield law with aspects from a line tension (LT) model for describing activation enthalpy required for the motion of dislocation kinks. Atomistic calculations are conducted to quantify the non-Schmid effects while both experimental data and atomistic simulations are used to assess the temperature and strain rate effects. The parameterized constitutive equationmore » is implemented into a BCC CP-FE model to simulate plastic deformation of single and polycrystalline Fe which is compared with experimental data from the literature. This direct comparison demonstrates that the atomistically informed model accurately captures the effects of crystal orientation, temperature and strain rate on the flow behavior of siangle crystal Fe. Furthermore, our proposed CP-FE model exhibits temperature and strain rate dependent flow and yield surfaces in polycrystalline Fe that deviate from conventional CP-FE models based on Schmid's law.« less
Tome, Carlos N; Caro, J A; Lebensohn, R A; Unal, Cetin; Arsenlis, A; Marian, J; Pasamehmetoglu, K
2010-01-01
Advancing the performance of Light Water Reactors, Advanced Nuclear Fuel Cycles, and Advanced Reactors, such as the Next Generation Nuclear Power Plants, requires enhancing our fundamental understanding of fuel and materials behavior under irradiation. The capability to accurately model the nuclear fuel systems to develop predictive tools is critical. Not only are fabrication and performance models needed to understand specific aspects of the nuclear fuel, fully coupled fuel simulation codes are required to achieve licensing of specific nuclear fuel designs for operation. The backbone of these codes, models, and simulations is a fundamental understanding and predictive capability for simulating the phase and microstructural behavior of the nuclear fuel system materials and matrices. In this paper we review the current status of the advanced modeling and simulation of nuclear reactor cladding, with emphasis on what is available and what is to be developed in each scale of the project, how we propose to pass information from one scale to the next, and what experimental information is required for benchmarking and advancing the modeling at each scale level.
Adomian Decomposition Method for Quark Gluon Plasma Model
Constantinescu, Radu; Ionescu, Carmen; Stoicescu, Mihai
2011-10-03
The paper investigates the possibility of obtaining analytical solutions for the Quark Gluon Plasma model using the Adomian decomposition method.
Yang, Judith C.
2015-01-09
The purpose of this grant is to develop the multi-scale theoretical methods to describe the nanoscale oxidation of metal thin films, as the PI (Yang) extensive previous experience in the experimental elucidation of the initial stages of Cu oxidation by primarily in situ transmission electron microscopy methods. Through the use and development of computational tools at varying length (and time) scales, from atomistic quantum mechanical calculation, force field mesoscale simulations, to large scale Kinetic Monte Carlo (KMC) modeling, the fundamental underpinings of the initial stages of Cu oxidation have been elucidated. The development of computational modeling tools allows for accelerated materials discovery. The theoretical tools developed from this program impact a wide range of technologies that depend on surface reactions, including corrosion, catalysis, and nanomaterials fabrication.
Liu, Dajiang [Ames Laboratory; Evans, James W. [Ames Laboratory
2013-12-01
A realistic molecular-level description of catalytic reactions on single-crystal metal surfaces can be provided by stochastic multisite lattice-gas (msLG) models. This approach has general applicability, although in this report, we will focus on the example of CO-oxidation on the unreconstructed fcc metal (100) or M(100) surfaces of common catalyst metals M = Pd, Rh, Pt and Ir (i.e., avoiding regimes where Pt and Ir reconstruct). These models can capture the thermodynamics and kinetics of adsorbed layers for the individual reactants species, such as CO/M(100) and O/M(100), as well as the interaction and reaction between different reactant species in mixed adlayers, such as (CO + O)/M(100). The msLG models allow population of any of hollow, bridge, and top sites. This enables a more flexible and realistic description of adsorption and adlayer ordering, as well as of reaction configurations and configuration-dependent barriers. Adspecies adsorption and interaction energies, as well as barriers for various processes, constitute key model input. The choice of these energies is guided by experimental observations, as well as by extensive Density Functional Theory analysis. Model behavior is assessed via Kinetic Monte Carlo (KMC) simulation. We also address the simulation challenges and theoretical ramifications associated with very rapid diffusion and local equilibration of reactant adspecies such as CO. These msLG models are applied to describe adsorption, ordering, and temperature programmed desorption (TPD) for individual CO/M(100) and O/M(100) reactant adlayers. In addition, they are also applied to predict mixed (CO + O)/M(100) adlayer structure on the nanoscale, the complete bifurcation diagram for reactive steady-states under continuous flow conditions, temperature programmed reaction (TPR) spectra, and titration reactions for the CO-oxidation reaction. Extensive and reasonably successful comparison of model predictions is made with experimental data. Furthermore
Vikas Tomer; John Renaud
2010-08-31
It is estimated that by using better and improved high temperature structural materials, the power generation efficiency of the power plants can be increased by 15% resulting in significant cost savings. One such promising material system for future high-temperature structural applications in power plants is Silicon Carbide-Silicon Nitride (SiC-Si{sub 3}N{sub 4}) nanoceramic matrix composites. The described research work focuses on multiscale simulation-based design of these SiC-Si{sub 3}N{sub 4} nanoceramic matrix composites. There were two primary objectives of the research: (1) Development of a multiscale simulation tool and corresponding multiscale analyses of the high-temperature creep and fracture resistance properties of the SiC-Si{sub 3}N{sub 4} nanocomposites at nano-, meso- and continuum length- and timescales; and (2) Development of a simulation-based robust design optimization methodology for application to the multiscale simulations to predict the range of the most suitable phase morphologies for the desired high-temperature properties of the SiC-Si{sub 3}N{sub 4} nanocomposites. The multiscale simulation tool is based on a combination of molecular dynamics (MD), cohesive finite element method (CFEM), and continuum level modeling for characterizing time-dependent material deformation behavior. The material simulation tool is incorporated in a variable fidelity model management based design optimization framework. Material modeling includes development of an experimental verification framework. Using material models based on multiscaling, it was found using molecular simulations that clustering of the SiC particles near Si{sub 3}N{sub 4} grain boundaries leads to significant nanocomposite strengthening and significant rise in fracture resistance. It was found that a control of grain boundary thicknesses by dispersing non-stoichiometric carbide or nitride phases can lead to reduction in strength however significant rise in fracture strength. The
Collaboratory for Multiscale Chemical Science (CMCS)
Allison, Thomas C
2012-07-03
This document provides details of the contributions made by NIST to the Collaboratory for Multiscale Chemical Science (CMCS) project. In particular, efforts related to the provision of data (and software in support of that data) relevant to the combustion pilot project are described.
Practical Use of Computationally Frugal Model Analysis Methods
Hill, Mary C.; Kavetski, Dmitri; Clark, Martyn; Ye, Ming; Arabi, Mazdak; Lu, Dan; Foglia, Laura; Mehl, Steffen
2015-03-21
Computationally frugal methods of model analysis can provide substantial benefits when developing models of groundwater and other environmental systems. Model analysis includes ways to evaluate model adequacy and to perform sensitivity and uncertainty analysis. Frugal methods typically require 10s of parallelizable model runs; their convenience allows for other uses of the computational effort. We suggest that model analysis be posed as a set of questions used to organize methods that range from frugal to expensive (requiring 10,000 model runs or more). This encourages focus on method utility, even when methods have starkly different theoretical backgrounds. We note that many frugal methods are more useful when unrealistic process-model nonlinearities are reduced. Inexpensive diagnostics are identified for determining when frugal methods are advantageous. Examples from the literature are used to demonstrate local methods and the diagnostics. We suggest that the greater use of computationally frugal model analysis methods would allow questions such as those posed in this work to be addressed more routinely, allowing the environmental sciences community to obtain greater scientific insight from the many ongoing and future modeling efforts
Practical Use of Computationally Frugal Model Analysis Methods
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Hill, Mary C.; Kavetski, Dmitri; Clark, Martyn; Ye, Ming; Arabi, Mazdak; Lu, Dan; Foglia, Laura; Mehl, Steffen
2015-03-21
Computationally frugal methods of model analysis can provide substantial benefits when developing models of groundwater and other environmental systems. Model analysis includes ways to evaluate model adequacy and to perform sensitivity and uncertainty analysis. Frugal methods typically require 10s of parallelizable model runs; their convenience allows for other uses of the computational effort. We suggest that model analysis be posed as a set of questions used to organize methods that range from frugal to expensive (requiring 10,000 model runs or more). This encourages focus on method utility, even when methods have starkly different theoretical backgrounds. We note that many frugalmore » methods are more useful when unrealistic process-model nonlinearities are reduced. Inexpensive diagnostics are identified for determining when frugal methods are advantageous. Examples from the literature are used to demonstrate local methods and the diagnostics. We suggest that the greater use of computationally frugal model analysis methods would allow questions such as those posed in this work to be addressed more routinely, allowing the environmental sciences community to obtain greater scientific insight from the many ongoing and future modeling efforts« less
Geoelectrical Measurement of Multi-Scale Mass Transfer Parameters
Day-Lewis, Frederick; Singha, Kamini; Haggerty, Roy; Johnson, Tim; Binley, Andrew; Lane, John
2014-01-16
-part research plan involving (1) development of computer codes and techniques to estimate mass-transfer parameters from time-lapse electrical data; (2) bench-scale experiments on synthetic materials and materials from cores from the Hanford 300 Area; and (3) field demonstration experiments at the DOE’s Hanford 300 Area. In a synergistic add-on to our workplan, we analyzed data from field experiments performed at the DOE Naturita Site under a separate DOE SBR grant, on which PI Day-Lewis served as co-PI. Techniques developed for application to Hanford datasets also were applied to data from Naturita. 1. Introduction The Department of Energy (DOE) faces enormous scientific and engineering challenges associated with the remediation of legacy contamination at former nuclear weapons production facilities. Selection, design and optimization of appropriate site remedies (e.g., pump-and-treat, biostimulation, or monitored natural attenuation) requires reliable predictive models of radionuclide fate and transport; however, our current modeling capabilities are limited by an incomplete understanding of multi-scale mass transfer—its rates, scales, and the heterogeneity of controlling parameters. At many DOE sites, long “tailing” behavior, concentration rebound, and slower-than-expected cleanup are observed; these observations are all consistent with multi-scale mass transfer [Haggerty and Gorelick, 1995; Haggerty et al., 2000; 2004], which renders pump-and-treat remediation and biotransformation inefficient and slow [Haggerty and Gorelick, 1994; Harvey et al., 1994; Wilson, 1997]. Despite the importance of mass transfer, there are significant uncertainties associated with controlling parameters, and the prevalence of mass transfer remains a point of debate [e.g., Hill et al., 2006; Molz et al., 2006] for lack of experimental methods to verify and measure it in situ or independently of tracer breakthrough. There is a critical need for new field-experimental techniques to
Weather Research and Forecasting Model with the Immersed Boundary Method
Energy Science and Technology Software Center (OSTI)
2012-05-01
The Weather Research and Forecasting (WRF) Model with the immersed boundary method is an extension of the open-source WRF Model available for wwww.wrf-model.org. The new code modifies the gridding procedure and boundary conditions in the WRF model to improve WRF's ability to simutate the atmosphere in environments with steep terrain and additionally at high-resolutions.
Multiscale modeling and characterization for performance and...
Office of Scientific and Technical Information (OSTI)
Lithium-ion batteries are highly complex electrochemical systems whose performance and safety are governed by coupled nonlinear electrochemical-electrical-thermal-mechanical ...
A New Computational Paradigm in Multiscale Simulations: Application to
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Brain Blood Flow | Argonne Leadership Computing Facility A New Computational Paradigm in Multiscale Simulations: Application to Brain Blood Flow Authors: Grinberg, L., Insley, J.A., Morozov, V., Papka, M.E., Karniadakis, G.E., Fedosov, D., Kumaran, K. Interfacing atomistic-based with continuum-based simulation codes is now required in many multiscale physical and biological systems. We present the computational advances that have enabled the first multiscale simulation on 190,740 processors
Systems and methods for modeling and analyzing networks
Hill, Colin C; Church, Bruce W; McDonagh, Paul D; Khalil, Iya G; Neyarapally, Thomas A; Pitluk, Zachary W
2013-10-29
The systems and methods described herein utilize a probabilistic modeling framework for reverse engineering an ensemble of causal models, from data and then forward simulating the ensemble of models to analyze and predict the behavior of the network. In certain embodiments, the systems and methods described herein include data-driven techniques for developing causal models for biological networks. Causal network models include computational representations of the causal relationships between independent variables such as a compound of interest and dependent variables such as measured DNA alterations, changes in mRNA, protein, and metabolites to phenotypic readouts of efficacy and toxicity.
A New Computational Paradigm in Multiscale Simulations: Application...
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
A New Computational Paradigm in Multiscale Simulations: Application to Brain Blood Flow ... We study blood flow in a patient-specific cerebrovasculature with a brain aneurysm, and ...
Multi-Scale Initial Conditions For Cosmological Simulations ...
Office of Scientific and Technical Information (OSTI)
Journal Article: Multi-Scale Initial Conditions For Cosmological Simulations Citation ... OSTI Identifier: 1028674 Report Number(s): SLAC-PUB-14463 Journal ID: ISSN 0035-8711; ...
Theory & Modeling | Argonne National Laboratory
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
interactions with nanostructures Methods and software development, including multiscale approaches to assembly Group Lead Stephen Gray People Maria K. Y. Chan Larry Curtiss...
An adaptive wavelet stochastic collocation method for irregular...
Office of Scientific and Technical Information (OSTI)
adaptive (MdMrA) sparse grid stochastic collocation method, that utilizes hierarchical multiscale piecewise Riesz basis functions constructed from interpolating wavelets. ...
Multiscale Computation. Needs and Opportunities for BER Science
Scheibe, Timothy D.; Smith, Jeremy C.
2015-01-01
The Environmental Molecular Sciences Laboratory (EMSL), a scientific user facility managed by Pacific Northwest National Laboratory for the U.S. Department of Energy, Office of Biological and Environmental Research (BER), conducted a one-day workshop on August 26, 2014 on the topic of “Multiscale Computation: Needs and Opportunities for BER Science.” Twenty invited participants, from various computational disciplines within the BER program research areas, were charged with the following objectives; Identify BER-relevant models and their potential cross-scale linkages that could be exploited to better connect molecular-scale research to BER research at larger scales and; Identify critical science directions that will motivate EMSL decisions regarding future computational (hardware and software) architectures.
Optical Measurement Methods used in Calibration and Validation of Modeled
Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site
Injection Spray Characteristics | Department of Energy Optical Measurement Methods used in Calibration and Validation of Modeled Injection Spray Characteristics Optical Measurement Methods used in Calibration and Validation of Modeled Injection Spray Characteristics Poster presented at the 16th Directions in Engine-Efficiency and Emissions Research (DEER) Conference in Detroit, MI, September 27-30, 2010. p-07_klyza.pdf (363.89 KB) More Documents & Publications Effect of Ambient Pressure
Modeling Microbes: New Methods for Integrated Metabolic and Regulatory
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Network Reconstruction | Argonne National Laboratory Modeling Microbes: New Methods for Integrated Metabolic and Regulatory Network Reconstruction September 8, 2016 10:30AM to 11:30PM Presenter Jose Faria Location Building 240, Room 4301 Type Seminar Series MCS Seminar Abstract: The reconstruction of genome-scale metabolic models (GEMs) from genome functional annotations is, nowadays, a routine practice in systems biology research. The models have been successfully used to predict organisms'
Self-Consistent Multiscale Theory of Internal Wave, Mean-Flow Interactions
Holm, D.D.; Aceves, A.; Allen, J.S.; Alber, M.; Camassa, R.; Cendra, H.; Chen, S.; Duan, J.; Fabijonas, B.; Foias, C.; Fringer, O.; Gent, P.R.; Jordan, R.; Kouranbaeva, S.; Kovacic, G.; Levermore, C.D.; Lythe, G.; Lifschitz, A.; Marsden, J.E.; Margolin, L.; Newberger, P.; Olson, E.; Ratiu, T.; Shkoller, S.; Timofeyev, I.; Titi, E.S.; Wynn, S.
1999-06-03
This is the final report of a three-year, Laboratory Directed Research and Development (LDRD) project at Los Alamos National Laboratory (LANL). The research reported here produced new effective ways to solve multiscale problems in nonlinear fluid dynamics, such as turbulent flow and global ocean circulation. This was accomplished by first developing new methods for averaging over random or rapidly varying phases in nonlinear systems at multiple scales. We then used these methods to derive new equations for analyzing the mean behavior of fluctuation processes coupled self consistently to nonlinear fluid dynamics. This project extends a technology base relevant to a variety of multiscale problems in fluid dynamics of interest to the Laboratory and applies this technology to those problems. The project's theoretical and mathematical developments also help advance our understanding of the scientific principles underlying the control of complex behavior in fluid dynamical systems with strong spatial and temporal internal variability.
Byrne, Jason P.; Morgan, Huw; Habbal, Shadia R.; Gallagher, Peter T.
2012-06-20
Studying coronal mass ejections (CMEs) in coronagraph data can be challenging due to their diffuse structure and transient nature, and user-specific biases may be introduced through visual inspection of the images. The large amount of data available from the Solar and Heliospheric Observatory (SOHO), Solar TErrestrial RElations Observatory (STEREO), and future coronagraph missions also makes manual cataloging of CMEs tedious, and so a robust method of detection and analysis is required. This has led to the development of automated CME detection and cataloging packages such as CACTus, SEEDS, and ARTEMIS. Here, we present the development of a new CORIMP (coronal image processing) CME detection and tracking technique that overcomes many of the drawbacks of current catalogs. It works by first employing the dynamic CME separation technique outlined in a companion paper, and then characterizing CME structure via a multiscale edge-detection algorithm. The detections are chained through time to determine the CME kinematics and morphological changes as it propagates across the plane of sky. The effectiveness of the method is demonstrated by its application to a selection of SOHO/LASCO and STEREO/SECCHI images, as well as to synthetic coronagraph images created from a model corona with a variety of CMEs. The algorithms described in this article are being applied to the whole LASCO and SECCHI data sets, and a catalog of results will soon be available to the public.
Probability of detection models for eddy current NDE methods
Rajesh, S.N.
1993-04-30
The development of probability of detection (POD) models for a variety of nondestructive evaluation (NDE) methods is motivated by a desire to quantify the variability introduced during the process of testing. Sources of variability involved in eddy current methods of NDE include those caused by variations in liftoff, material properties, probe canting angle, scan format, surface roughness and measurement noise. This thesis presents a comprehensive POD model for eddy current NDE. Eddy current methods of nondestructive testing are used widely in industry to inspect a variety of nonferromagnetic and ferromagnetic materials. The development of a comprehensive POD model is therefore of significant importance. The model incorporates several sources of variability characterized by a multivariate Gaussian distribution and employs finite element analysis to predict the signal distribution. The method of mixtures is then used for estimating optimal threshold values. The research demonstrates the use of a finite element model within a probabilistic framework to the spread in the measured signal for eddy current nondestructive methods. Using the signal distributions for various flaw sizes the POD curves for varying defect parameters have been computed. In contrast to experimental POD models, the cost of generating such curves is very low and complex defect shapes can be handled very easily. The results are also operator independent.
Multiscale Monte Carlo equilibration: Pure Yang-Mills theory
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Endres, Michael G.; Brower, Richard C.; Orginos, Kostas; Detmold, William; Pochinsky, Andrew V.
2015-12-29
In this study, we present a multiscale thermalization algorithm for lattice gauge theory, which enables efficient parallel generation of uncorrelated gauge field configurations. The algorithm combines standard Monte Carlo techniques with ideas drawn from real space renormalization group and multigrid methods. We demonstrate the viability of the algorithm for pure Yang-Mills gauge theory for both heat bath and hybrid Monte Carlo evolution, and show that it ameliorates the problem of topological freezing up to controllable lattice spacing artifacts.
Material point method modeling in oil and gas reservoirs
Vanderheyden, William Brian; Zhang, Duan
2016-06-28
A computer system and method of simulating the behavior of an oil and gas reservoir including changes in the margins of frangible solids. A system of equations including state equations such as momentum, and conservation laws such as mass conservation and volume fraction continuity, are defined and discretized for at least two phases in a modeled volume, one of which corresponds to frangible material. A material point model technique for numerically solving the system of discretized equations, to derive fluid flow at each of a plurality of mesh nodes in the modeled volume, and the velocity of at each of a plurality of particles representing the frangible material in the modeled volume. A time-splitting technique improves the computational efficiency of the simulation while maintaining accuracy on the deformation scale. The method can be applied to derive accurate upscaled model equations for larger volume scale simulations.
Localized Scale Coupling and New Educational Paradigms in Multiscale Mathematics and Science
LEAL, L. GARY
2013-06-30
One of the most challenging multi-scale simulation problems in the area of multi-phase materials is to develop effective computational techniques for the prediction of coalescence and related phenomena involving rupture of a thin liquid film due to the onset of instability driven by van der Waals or other micro-scale attractive forces. Accurate modeling of this process is critical to prediction of the outcome of milling processes for immiscible polymer blends, one of the most important routes to new advanced polymeric materials. In typical situations, the blend evolves into an ?emulsion? of dispersed phase drops in a continuous matrix fluid. Coalescence is then a critical factor in determining the size distribution of the dispersed phase, but is extremely difficult to predict from first principles. The thin film separating two drops may only achieve rupture at dimensions of approximately 10 nm while the drop sizes are 0(10 ?m). It is essential to achieve very accurate solutions for the flow and for the interface shape at both the macroscale of the full drops, and within the thin film (where the destabilizing disjoining pressure due to van der Waals forces is proportional approximately to the inverse third power of the local film thickness, h-3). Furthermore, the fluids of interest are polymeric (through Newtonian) and the classical continuum description begins to fail as the film thins ? requiring incorporation of molecular effects, such as a hybrid code that incorporates a version of coarse grain molecular dynamics within the thin film coupled with a classical continuum description elsewhere in the flow domain. Finally, the presence of surface active additions, either surfactants (in the form of di-block copolymers) or surface-functionalized micro- or nano-scale particles, adds an additional level of complexity, requiring development of a distinct numerical method to predict the nonuniform concentration gradients of these additives that are responsible for
High-Order/Low-Order methods for ocean modeling
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Newman, Christopher; Womeldorff, Geoff; Chacón, Luis; Knoll, Dana A.
2015-06-01
We examine a High Order/Low Order (HOLO) approach for a z-level ocean model and show that the traditional semi-implicit and split-explicit methods, as well as a recent preconditioning strategy, can easily be cast in the framework of HOLO methods. The HOLO formulation admits an implicit-explicit method that is algorithmically scalable and second-order accurate, allowing timesteps much larger than the barotropic time scale. We demonstrate how HOLO approaches, in particular the implicit-explicit method, can provide a solid route for ocean simulation to heterogeneous computing and exascale environments.
Curve fitting methods for solar radiation data modeling
Karim, Samsul Ariffin Abdul E-mail: balbir@petronas.com.my; Singh, Balbir Singh Mahinder E-mail: balbir@petronas.com.my
2014-10-24
This paper studies the use of several type of curve fitting method to smooth the global solar radiation data. After the data have been fitted by using curve fitting method, the mathematical model of global solar radiation will be developed. The error measurement was calculated by using goodness-fit statistics such as root mean square error (RMSE) and the value of R{sup 2}. The best fitting methods will be used as a starting point for the construction of mathematical modeling of solar radiation received in Universiti Teknologi PETRONAS (UTP) Malaysia. Numerical results indicated that Gaussian fitting and sine fitting (both with two terms) gives better results as compare with the other fitting methods.
A meshless method for modeling convective heat transfer
Carrington, David B
2010-01-01
A meshless method is used in a projection-based approach to solve the primitive equations for fluid flow with heat transfer. The method is easy to implement in a MATLAB format. Radial basis functions are used to solve two benchmark test cases: natural convection in a square enclosure and flow with forced convection over a backward facing step. The results are compared with two popular and widely used commercial codes: COMSOL, a finite element model, and FLUENT, a finite volume-based model.
Spectral characteristics of background error covariance and multiscale data assimilation
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Li, Zhijin; Cheng, Xiaoping; Gustafson, Jr., William I.; Vogelmann, Andrew M.
2016-05-17
The steady increase of the spatial resolutions of numerical atmospheric and oceanic circulation models has occurred over the past decades. Horizontal grid spacing down to the order of 1 km is now often used to resolve cloud systems in the atmosphere and sub-mesoscale circulation systems in the ocean. These fine resolution models encompass a wide range of temporal and spatial scales, across which dynamical and statistical properties vary. In particular, dynamic flow systems at small scales can be spatially localized and temporarily intermittent. Difficulties of current data assimilation algorithms for such fine resolution models are numerically and theoretically examined. Ourmore » analysis shows that the background error correlation length scale is larger than 75 km for streamfunctions and is larger than 25 km for water vapor mixing ratios, even for a 2-km resolution model. A theoretical analysis suggests that such correlation length scales prevent the currently used data assimilation schemes from constraining spatial scales smaller than 150 km for streamfunctions and 50 km for water vapor mixing ratios. Moreover, our results highlight the need to fundamentally modify currently used data assimilation algorithms for assimilating high-resolution observations into the aforementioned fine resolution models. Lastly, within the framework of four-dimensional variational data assimilation, a multiscale methodology based on scale decomposition is suggested and challenges are discussed.« less
Multilevel method for modeling large-scale networks.
Safro, I. M.
2012-02-24
Understanding the behavior of real complex networks is of great theoretical and practical significance. It includes developing accurate artificial models whose topological properties are similar to the real networks, generating the artificial networks at different scales under special conditions, investigating a network dynamics, reconstructing missing data, predicting network response, detecting anomalies and other tasks. Network generation, reconstruction, and prediction of its future topology are central issues of this field. In this project, we address the questions related to the understanding of the network modeling, investigating its structure and properties, and generating artificial networks. Most of the modern network generation methods are based either on various random graph models (reinforced by a set of properties such as power law distribution of node degrees, graph diameter, and number of triangles) or on the principle of replicating an existing model with elements of randomization such as R-MAT generator and Kronecker product modeling. Hierarchical models operate at different levels of network hierarchy but with the same finest elements of the network. However, in many cases the methods that include randomization and replication elements on the finest relationships between network nodes and modeling that addresses the problem of preserving a set of simplified properties do not fit accurately enough the real networks. Among the unsatisfactory features are numerically inadequate results, non-stability of algorithms on real (artificial) data, that have been tested on artificial (real) data, and incorrect behavior at different scales. One reason is that randomization and replication of existing structures can create conflicts between fine and coarse scales of the real network geometry. Moreover, the randomization and satisfying of some attribute at the same time can abolish those topological attributes that have been undefined or hidden from
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Li, Zhijin; Cheng, Xiaoping; Gustafson, William I.; Vogelmann, Andrew M.
2016-05-17
The steady increase of the spatial resolutions of numerical atmospheric and oceanic circulation models has occurred over the past decades. Horizontal grid spacing down to the order of 1 km is now often used to resolve cloud systems in the atmosphere and sub-mesoscale circulation systems in the ocean. These fine resolution models encompass a wide range of temporal and spatial scales, across which dynamical and statistical properties vary. In particular, dynamic flow systems at small scales can be spatially localized and temporarily intermittent. Difficulties of current data assimilation algorithms for such fine resolution models are numerically and theoretically examined. Ourmore » analysis shows that the background error correlation length scale is larger than 75 km for streamfunctions and is larger than 25 km for water vapor mixing ratios, even for a 2-km resolution model. A theoretical analysis suggests that such correlation length scales prevent the currently used data assimilation schemes from constraining spatial scales smaller than 150 km for streamfunctions and 50 km for water vapor mixing ratios. Moreover, our results highlight the need to fundamentally modify currently used data assimilation algorithms for assimilating high-resolution observations into the aforementioned fine resolution models. Within the framework of four-dimensional variational data assimilation, a multiscale methodology based on scale decomposition is suggested and challenges are discussed.« less
(Environmental and geophysical modeling, fracture mechanics, and boundary element methods)
Gray, L.J.
1990-11-09
Technical discussions at the various sites visited centered on application of boundary integral methods for environmental modeling, seismic analysis, and computational fracture mechanics in composite and smart'' materials. The traveler also attended the International Association for Boundary Element Methods Conference at Rome, Italy. While many aspects of boundary element theory and applications were discussed in the papers, the dominant topic was the analysis and application of hypersingular equations. This has been the focus of recent work by the author, and thus the conference was highly relevant to research at ORNL.
Bayesian methods for characterizing unknown parameters of material models
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Emery, J. M.; Grigoriu, M. D.; Field Jr., R. V.
2016-02-04
A Bayesian framework is developed for characterizing the unknown parameters of probabilistic models for material properties. In this framework, the unknown parameters are viewed as random and described by their posterior distributions obtained from prior information and measurements of quantities of interest that are observable and depend on the unknown parameters. The proposed Bayesian method is applied to characterize an unknown spatial correlation of the conductivity field in the definition of a stochastic transport equation and to solve this equation by Monte Carlo simulation and stochastic reduced order models (SROMs). As a result, the Bayesian method is also employed tomore » characterize unknown parameters of material properties for laser welds from measurements of peak forces sustained by these welds.« less
Arctic sea ice modeling with the material-point method.
Peterson, Kara J.; Bochev, Pavel Blagoveston
2010-04-01
Arctic sea ice plays an important role in global climate by reflecting solar radiation and insulating the ocean from the atmosphere. Due to feedback effects, the Arctic sea ice cover is changing rapidly. To accurately model this change, high-resolution calculations must incorporate: (1) annual cycle of growth and melt due to radiative forcing; (2) mechanical deformation due to surface winds, ocean currents and Coriolis forces; and (3) localized effects of leads and ridges. We have demonstrated a new mathematical algorithm for solving the sea ice governing equations using the material-point method with an elastic-decohesive constitutive model. An initial comparison with the LANL CICE code indicates that the ice edge is sharper using Materials-Point Method (MPM), but that many of the overall features are similar.
A New Method of Comparing Forcing Agents in Climate Models
Kravitz, Benjamin S.; MacMartin, Douglas; Rasch, Philip J.; Jarvis, Andrew
2015-10-14
We describe a new method of comparing different climate forcing agents (e.g., CO2, CH4, and solar irradiance) that avoids many of the ambiguities introduced by temperature-related climate feedbacks. This is achieved by introducing an explicit feedback loop external to the climate model that adjusts one forcing agent to balance another while keeping global mean surface temperature constant. Compared to current approaches, this method has two main advantages: (i) the need to define radiative forcing is bypassed and (ii) by maintaining roughly constant global mean temperature, the effects of state dependence on internal feedback strengths are minimized. We demonstrate this approach for several different forcing agents and derive the relationships between these forcing agents in two climate models; comparisons between forcing agents are highly linear in concordance with predicted functional forms. Transitivity of the relationships between the forcing agents appears to hold within a wide range of forcing. The relationships between the forcing agents obtained from this method are consistent across both models but differ from relationships that would be obtained from calculations of radiative forcing, highlighting the importance of controlling for surface temperature feedback effects when separating radiative forcing and climate response.
Editorial: Mathematical Methods and Modeling in Machine Fault Diagnosis
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Yan, Ruqiang; Chen, Xuefeng; Li, Weihua; Sheng, Shuangwen
2014-12-18
Modern mathematics has commonly been utilized as an effective tool to model mechanical equipment so that their dynamic characteristics can be studied analytically. This will help identify potential failures of mechanical equipment by observing change in the equipment’s dynamic parameters. On the other hand, dynamic signals are also important and provide reliable information about the equipment’s working status. Modern mathematics has also provided us with a systematic way to design and implement various signal processing methods, which are used to analyze these dynamic signals, and to enhance intrinsic signal components that are directly related to machine failures. This special issuemore » is aimed at stimulating not only new insights on mathematical methods for modeling but also recently developed signal processing methods, such as sparse decomposition with potential applications in machine fault diagnosis. Finally, the papers included in this special issue provide a glimpse into some of the research and applications in the field of machine fault diagnosis through applications of the modern mathematical methods.« less
Editorial: Mathematical Methods and Modeling in Machine Fault Diagnosis
Yan, Ruqiang; Chen, Xuefeng; Li, Weihua; Sheng, Shuangwen
2014-12-18
Modern mathematics has commonly been utilized as an effective tool to model mechanical equipment so that their dynamic characteristics can be studied analytically. This will help identify potential failures of mechanical equipment by observing change in the equipment’s dynamic parameters. On the other hand, dynamic signals are also important and provide reliable information about the equipment’s working status. Modern mathematics has also provided us with a systematic way to design and implement various signal processing methods, which are used to analyze these dynamic signals, and to enhance intrinsic signal components that are directly related to machine failures. This special issue is aimed at stimulating not only new insights on mathematical methods for modeling but also recently developed signal processing methods, such as sparse decomposition with potential applications in machine fault diagnosis. Finally, the papers included in this special issue provide a glimpse into some of the research and applications in the field of machine fault diagnosis through applications of the modern mathematical methods.
Progress in Fast, Accurate Multi-scale Climate Simulations
Collins, William D; Johansen, Hans; Evans, Katherine J; Woodward, Carol S.; Caldwell, Peter
2015-01-01
We present a survey of physical and computational techniques that have the potential to con- tribute to the next generation of high-fidelity, multi-scale climate simulations. Examples of the climate science problems that can be investigated with more depth include the capture of remote forcings of localized hydrological extreme events, an accurate representation of cloud features over a range of spatial and temporal scales, and parallel, large ensembles of simulations to more effectively explore model sensitivities and uncertainties. Numerical techniques, such as adaptive mesh refinement, implicit time integration, and separate treatment of fast physical time scales are enabling improved accuracy and fidelity in simulation of dynamics and allow more complete representations of climate features at the global scale. At the same time, part- nerships with computer science teams have focused on taking advantage of evolving computer architectures, such as many-core processors and GPUs, so that these approaches which were previously considered prohibitively costly have become both more efficient and scalable. In combination, progress in these three critical areas is poised to transform climate modeling in the coming decades.
Progress in fast, accurate multi-scale climate simulations
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Collins, W. D.; Johansen, H.; Evans, K. J.; Woodward, C. S.; Caldwell, P. M.
2015-06-01
We present a survey of physical and computational techniques that have the potential to contribute to the next generation of high-fidelity, multi-scale climate simulations. Examples of the climate science problems that can be investigated with more depth with these computational improvements include the capture of remote forcings of localized hydrological extreme events, an accurate representation of cloud features over a range of spatial and temporal scales, and parallel, large ensembles of simulations to more effectively explore model sensitivities and uncertainties. Numerical techniques, such as adaptive mesh refinement, implicit time integration, and separate treatment of fast physical time scales are enablingmore » improved accuracy and fidelity in simulation of dynamics and allowing more complete representations of climate features at the global scale. At the same time, partnerships with computer science teams have focused on taking advantage of evolving computer architectures such as many-core processors and GPUs. As a result, approaches which were previously considered prohibitively costly have become both more efficient and scalable. In combination, progress in these three critical areas is poised to transform climate modeling in the coming decades.« less
Progress in fast, accurate multi-scale climate simulations
Collins, W. D.; Johansen, H.; Evans, K. J.; Woodward, C. S.; Caldwell, P. M.
2015-06-01
We present a survey of physical and computational techniques that have the potential to contribute to the next generation of high-fidelity, multi-scale climate simulations. Examples of the climate science problems that can be investigated with more depth with these computational improvements include the capture of remote forcings of localized hydrological extreme events, an accurate representation of cloud features over a range of spatial and temporal scales, and parallel, large ensembles of simulations to more effectively explore model sensitivities and uncertainties. Numerical techniques, such as adaptive mesh refinement, implicit time integration, and separate treatment of fast physical time scales are enabling improved accuracy and fidelity in simulation of dynamics and allowing more complete representations of climate features at the global scale. At the same time, partnerships with computer science teams have focused on taking advantage of evolving computer architectures such as many-core processors and GPUs. As a result, approaches which were previously considered prohibitively costly have become both more efficient and scalable. In combination, progress in these three critical areas is poised to transform climate modeling in the coming decades.
An Efficient Implementation of Multiscale Simulation Software PNP-cDFT
Meng, Da; Lin, Guang; Sushko, Maria L.
2012-07-23
An efficient implementation of PNP-cDFT, a multiscale method for computing the chemical potentials of charged species is designed and evaluated. Spatial decomposition of the multi particle system is employed in the parallelization of classical density functional theory (cDFT) algorithm. Furthermore, a truncated strategy is used to reduce the computational complexity of cDFT algorithm. The simulation results show that the parallel implementation has close to linear scalability in parallel computing environments for both 1D and 3D systems. It also shows that the truncated versions of cDFT improve the efficiency of the methods substantially.
Thermodynamic Development of Corrosion Rate Modeling in Iron...
Office of Scientific and Technical Information (OSTI)
... multiscale neutronics research which was performed ... deterministic transport method was not available, this ... transient of a flow reversal and is qualitative in nature. ...
Bridging the PSI Knowledge Gap: A Multi-Scale Approach
Wirth, Brian D
2015-01-08
Plasma-surface interactions (PSI) pose an immense scientific hurdle in magnetic confinement fusion and our present understanding of PSI in confinement environments is highly inadequate; indeed, a recent Fusion Energy Sciences Advisory Committee report found that 4 out of the 5 top five fusion knowledge gaps were related to PSI. The time is appropriate to develop a concentrated and synergistic science effort that would expand, exploit and integrate the wealth of laboratory ion-beam and plasma research, as well as exciting new computational tools, towards the goal of bridging the PSI knowledge gap. This effort would broadly advance plasma and material sciences, while providing critical knowledge towards progress in fusion PSI. This project involves the development of a Science Center focused on a new approach to PSI science; an approach that both exploits access to state-of-the-art PSI experiments and modeling, as well as confinement devices. The organizing principle is to develop synergistic experimental and modeling tools that treat the truly coupled multi-scale aspect of the PSI issues in confinement devices. This is motivated by the simple observation that while typical lab experiments and models allow independent manipulation of controlling variables, the confinement PSI environment is essentially self-determined with few outside controls. This means that processes that may be treated independently in laboratory experiments, because they involve vastly different physical and time scales, will now affect one another in the confinement environment. Also, lab experiments cannot simultaneously match all exposure conditions found in confinement devices typically forcing a linear extrapolation of lab results. At the same time programmatic limitations prevent confinement experiments alone from answering many key PSI questions. The resolution to this problem is to usefully exploit access to PSI science in lab devices, while retooling our thinking from a linear and de
Next Generation Multi-Scale Quantum Simulation Software for Strongly Correlated Materials
Jarrell, Mark
2014-11-18
The goal of this project was to develop a new formalism for the correlated electron problem, which we call, the Multi Scale Many Body formalism. This report will focus on the work done at the Louisiana State University (LSU) since the mid term report. The LSU group moved from the University of Cincinnati (UC) to LSU in the summer of 2008. In the last full year at UC, only half of the funds were received and it took nearly two years for the funds to be transferred from UC to LSU . This effectively shut down the research at LSU until the transfer was completed in 2011, there were also two no-cost extensions of the grant until August of this year. The grant ended for the other SciDAC partners at Davis and ORNL in 2011. Since the mid term report, the LSU group has published 19 papers [P1-P19] acknowledging this SciDAC, which are listed below. In addition, numerous invited talked acknowledged the SciDAC. Below, we will summarize the work at LSU since the mid-term report and mainly since funding resumed. The projects include the further development of multi-scale methods for correlated systems (1), the study of quantum criticality at finite doping in the Hubbard model (2), the description of a promising new method to study Anderson localization with a million-fold reduction of computational complexity!, the description of other projects (4), and (5) a workshop to close out the project that brought together exascale program developers (Stellar, MPI, OpenMP,...) with applications developers.
Methods for Developing Emissions Scenarios for Integrated Assessment Models
Prinn, Ronald; Webster, Mort
2007-08-20
The overall objective of this research was to contribute data and methods to support the future development of new emissions scenarios for integrated assessment of climate change. Specifically, this research had two main objectives: 1. Use historical data on economic growth and energy efficiency changes, and develop probability density functions (PDFs) for the appropriate parameters for two or three commonly used integrated assessment models. 2. Using the parameter distributions developed through the first task and previous work, we will develop methods of designing multi-gas emission scenarios that usefully span the joint uncertainty space in a small number of scenarios. Results on the autonomous energy efficiency improvement (AEEI) parameter are summarized, an uncertainty analysis of elasticities of substitution is described, and the probabilistic emissions scenario approach is presented.
Modeling granular phosphor screens by Monte Carlo methods
Liaparinos, Panagiotis F.; Kandarakis, Ioannis S.; Cavouras, Dionisis A.; Delis, Harry B.; Panayiotakis, George S.
2006-12-15
The intrinsic phosphor properties are of significant importance for the performance of phosphor screens used in medical imaging systems. In previous analytical-theoretical and Monte Carlo studies on granular phosphor materials, values of optical properties, and light interaction cross sections were found by fitting to experimental data. These values were then employed for the assessment of phosphor screen imaging performance. However, it was found that, depending on the experimental technique and fitting methodology, the optical parameters of a specific phosphor material varied within a wide range of values, i.e., variations of light scattering with respect to light absorption coefficients were often observed for the same phosphor material. In this study, x-ray and light transport within granular phosphor materials was studied by developing a computational model using Monte Carlo methods. The model was based on the intrinsic physical characteristics of the phosphor. Input values required to feed the model can be easily obtained from tabulated data. The complex refractive index was introduced and microscopic probabilities for light interactions were produced, using Mie scattering theory. Model validation was carried out by comparing model results on x-ray and light parameters (x-ray absorption, statistical fluctuations in the x-ray to light conversion process, number of emitted light photons, output light spatial distribution) with previous published experimental data on Gd{sub 2}O{sub 2}S:Tb phosphor material (Kodak Min-R screen). Results showed the dependence of the modulation transfer function (MTF) on phosphor grain size and material packing density. It was predicted that granular Gd{sub 2}O{sub 2}S:Tb screens of high packing density and small grain size may exhibit considerably better resolution and light emission properties than the conventional Gd{sub 2}O{sub 2}S:Tb screens, under similar conditions (x-ray incident energy, screen thickness)
A robust absorbing layer method for anisotropic seismic wave modeling
Mtivier, L.; Brossier, R.; Labb, S.; Operto, S.; Virieux, J.
2014-12-15
When applied to wave propagation modeling in anisotropic media, Perfectly Matched Layers (PML) exhibit instabilities. Incoming waves are amplified instead of being absorbed. Overcoming this difficulty is crucial as in many seismic imaging applications, accounting accurately for the subsurface anisotropy is mandatory. In this study, we present the SMART layer method as an alternative to PML approach. This method is based on the decomposition of the wavefield into components propagating inward and outward the domain of interest. Only outgoing components are damped. We show that for elastic and acoustic wave propagation in Transverse Isotropic media, the SMART layer is unconditionally dissipative: no amplification of the wavefield is possible. The SMART layers are not perfectly matched, therefore less accurate than conventional PML. However, a reasonable increase of the layer size yields an accuracy similar to PML. Finally, we illustrate that the selective damping strategy on which is based the SMART method can prevent the generation of spurious S-waves by embedding the source in a small zone where only S-waves are damped.
Neural node network and model, and method of teaching same
Parlos, Alexander G.; Atiya, Amir F.; Fernandez, Benito; Tsai, Wei K.; Chong, Kil T.
1995-01-01
The present invention is a fully connected feed forward network that includes at least one hidden layer 16. The hidden layer 16 includes nodes 20 in which the output of the node is fed back to that node as an input with a unit delay produced by a delay device 24 occurring in the feedback path 22 (local feedback). Each node within each layer also receives a delayed output (crosstalk) produced by a delay unit 36 from all the other nodes within the same layer 16. The node performs a transfer function operation based on the inputs from the previous layer and the delayed outputs. The network can be implemented as analog or digital or within a general purpose processor. Two teaching methods can be used: (1) back propagation of weight calculation that includes the local feedback and the crosstalk or (2) more preferably a feed forward gradient decent which immediately follows the output computations and which also includes the local feedback and the crosstalk. Subsequent to the gradient propagation, the weights can be normalized, thereby preventing convergence to a local optimum. Education of the network can be incremental both on and off-line. An educated network is suitable for modeling and controlling dynamic nonlinear systems and time series systems and predicting the outputs as well as hidden states and parameters. The educated network can also be further educated during on-line processing.
Neural node network and model, and method of teaching same
Parlos, A.G.; Atiya, A.F.; Fernandez, B.; Tsai, W.K.; Chong, K.T.
1995-12-26
The present invention is a fully connected feed forward network that includes at least one hidden layer. The hidden layer includes nodes in which the output of the node is fed back to that node as an input with a unit delay produced by a delay device occurring in the feedback path (local feedback). Each node within each layer also receives a delayed output (crosstalk) produced by a delay unit from all the other nodes within the same layer. The node performs a transfer function operation based on the inputs from the previous layer and the delayed outputs. The network can be implemented as analog or digital or within a general purpose processor. Two teaching methods can be used: (1) back propagation of weight calculation that includes the local feedback and the crosstalk or (2) more preferably a feed forward gradient decent which immediately follows the output computations and which also includes the local feedback and the crosstalk. Subsequent to the gradient propagation, the weights can be normalized, thereby preventing convergence to a local optimum. Education of the network can be incremental both on and off-line. An educated network is suitable for modeling and controlling dynamic nonlinear systems and time series systems and predicting the outputs as well as hidden states and parameters. The educated network can also be further educated during on-line processing. 21 figs.
Multi-Scale Initial Conditions For Cosmological Simulations
Hahn, Oliver; Abel, Tom; /KIPAC, Menlo Park /ZAH, Heidelberg /HITS, Heidelberg
2011-11-04
We discuss a new algorithm to generate multi-scale initial conditions with multiple levels of refinements for cosmological 'zoom-in' simulations. The method uses an adaptive convolution of Gaussian white noise with a real-space transfer function kernel together with an adaptive multi-grid Poisson solver to generate displacements and velocities following first- (1LPT) or second-order Lagrangian perturbation theory (2LPT). The new algorithm achieves rms relative errors of the order of 10{sup -4} for displacements and velocities in the refinement region and thus improves in terms of errors by about two orders of magnitude over previous approaches. In addition, errors are localized at coarse-fine boundaries and do not suffer from Fourier-space-induced interference ringing. An optional hybrid multi-grid and Fast Fourier Transform (FFT) based scheme is introduced which has identical Fourier-space behaviour as traditional approaches. Using a suite of re-simulations of a galaxy cluster halo our real-space-based approach is found to reproduce correlation functions, density profiles, key halo properties and subhalo abundances with per cent level accuracy. Finally, we generalize our approach for two-component baryon and dark-matter simulations and demonstrate that the power spectrum evolution is in excellent agreement with linear perturbation theory. For initial baryon density fields, it is suggested to use the local Lagrangian approximation in order to generate a density field for mesh-based codes that is consistent with the Lagrangian perturbation theory instead of the current practice of using the Eulerian linearly scaled densities.
Assessment of Molecular Modeling & Simulation
2002-01-03
This report reviews the development and applications of molecular and materials modeling in Europe and Japan in comparison to those in the United States. Topics covered include computational quantum chemistry, molecular simulations by molecular dynamics and Monte Carlo methods, mesoscale modeling of material domains, molecular-structure/macroscale property correlations like QSARs and QSPRs, and related information technologies like informatics and special-purpose molecular-modeling computers. The panel's findings include the following: The United States leads this field in many scientific areas. However, Canada has particular strengths in DFT methods and homogeneous catalysis; Europe in heterogeneous catalysis, mesoscale, and materials modeling; and Japan in materials modeling and special-purpose computing. Major government-industry initiatives are underway in Europe and Japan, notably in multi-scale materials modeling and in development of chemistry-capable ab-initio molecular dynamics codes.
Symmetry Methods for a Geophysical Mass Flow Model
Torrisi, Mariano; Tracina, Rita
2011-09-14
In the framework of symmetry analysis, the class of 2 x 2 PDE systems to whom belong the Savage and Hutter model and the Iverson model is considered. New classes of exact solutions are found.
A Perspective on Coupled Multiscale Simulation and Validation in Nuclear Materials
M. P. Short; D. Gaston; C. R. Stanek; S. Yip
2014-01-01
The field of nuclear materials encompasses numerous opportunities to address and ultimately solve longstanding industrial problems by improving the fundamental understanding of materials through the integration of experiments with multiscale modeling and high-performance simulation. A particularly noteworthy example is an ongoing study of axial power distortions in a nuclear reactor induced by corrosion deposits, known as CRUD (Chalk River unidentified deposits). We describe how progress is being made toward achieving scientific advances and technological solutions on two fronts. Specifically, the study of thermal conductivity of CRUD phases has augmented missing data as well as revealed new mechanisms. Additionally, the development of a multiscale simulation framework shows potential for the validation of a new capability to predict the power distribution of a reactor, in effect direct evidence of technological impact. The material- and system-level challenges identified in the study of CRUD are similar to other well-known vexing problems in nuclear materials, such as irradiation accelerated corrosion, stress corrosion cracking, and void swelling; they all involve connecting materials science fundamentals at the atomistic- and mesoscales to technology challenges at the macroscale.
Trabanino, Rene J; Vaidehi, Nagarajan; Hall, Spencer E; Goddard, William A; Floriano, Wely
2013-02-05
The invention provides computer-implemented methods and apparatus implementing a hierarchical protocol using multiscale molecular dynamics and molecular modeling methods to predict the presence of transmembrane regions in proteins, such as G-Protein Coupled Receptors (GPCR), and protein structural models generated according to the protocol. The protocol features a coarse grain sampling method, such as hydrophobicity analysis, to provide a fast and accurate procedure for predicting transmembrane regions. Methods and apparatus of the invention are useful to screen protein or polynucleotide databases for encoded proteins with transmembrane regions, such as GPCRs.
Piri, Mohammad
2014-03-31
Under this project, a multidisciplinary team of researchers at the University of Wyoming combined state-of-the-art experimental studies, numerical pore- and reservoir-scale modeling, and high performance computing to investigate trapping mechanisms relevant to geologic storage of mixed scCO{sub 2} in deep saline aquifers. The research included investigations in three fundamental areas: (i) the experimental determination of two-‐phase flow relative permeability functions, relative permeability hysteresis, and residual trapping under reservoir conditions for mixed scCO{sub 2}-‐brine systems; (ii) improved understanding of permanent trapping mechanisms; (iii) scientifically correct, fine grid numerical simulations of CO{sub 2} storage in deep saline aquifers taking into account the underlying rock heterogeneity. The specific activities included: (1) Measurement of reservoir-‐conditions drainage and imbibition relative permeabilities, irreducible brine and residual mixed scCO{sub 2} saturations, and relative permeability scanning curves (hysteresis) in rock samples from RSU; (2) Characterization of wettability through measurements of contact angles and interfacial tensions under reservoir conditions; (3) Development of physically-‐based dynamic core-‐scale pore network model; (4) Development of new, improved high-‐ performance modules for the UW-‐team simulator to provide new capabilities to the existing model to include hysteresis in the relative permeability functions, geomechanical deformation and an equilibrium calculation (Both pore-‐ and core-‐scale models were rigorously validated against well-‐characterized core-‐ flooding experiments); and (5) An analysis of long term permanent trapping of mixed scCO{sub 2} through high-‐resolution numerical experiments and analytical solutions. The analysis takes into account formation heterogeneity, capillary trapping, and relative permeability hysteresis.
Multiscale Toxicology- Building the Next Generation Tools for Toxicology
Retterer, S. T.; Holsapple, M. P.
2013-10-31
A Cooperative Research and Development Agreement (CRADA) was established between Battelle Memorial Institute (BMI), Pacific Northwest National Laboratory (PNNL), Oak Ridge National Laboratory (ORNL), Brookhaven National Laboratory (BNL), Lawrence Livermore National Laboratory (LLNL) with the goal of combining the analytical and synthetic strengths of the National Laboratories with BMI's expertise in basic and translational medical research to develop a collaborative pipeline and suite of high throughput and imaging technologies that could be used to provide a more comprehensive understanding of material and drug toxicology in humans. The Multi-Scale Toxicity Initiative (MSTI), consisting of the team members above, was established to coordinate cellular scale, high-throughput in vitro testing, computational modeling and whole animal in vivo toxicology studies between MSTI team members. Development of a common, well-characterized set of materials for testing was identified as a crucial need for the initiative. Two research tracks were established by BMI during the course of the CRADA. The first research track focused on the development of tools and techniques for understanding the toxicity of nanomaterials, specifically inorganic nanoparticles (NPs). ORNL"s work focused primarily on the synthesis, functionalization and characterization of a common set of NPs for dissemination to the participating laboratories. These particles were synthesized to retain the same surface characteristics and size, but to allow visualization using the variety of imaging technologies present across the team. Characterization included the quantitative analysis of physical and chemical properties of the materials as well as the preliminary assessment of NP toxicity using commercially available toxicity screens and emerging optical imaging strategies. Additional efforts examined the development of high-throughput microfluidic and imaging assays for measuring NP uptake, localization, and
Multi-scale framework for the accelerated design of high-efficiency...
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Multi-scale framework for the accelerated design of high-efficiency organic photovoltaic cells Organic and hybrid organicinorganic solar cells (OSC) offer a promising low-cost...
Putting ''place'' a multiscale context: perspectices from the sustainability sciences
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Wilbanks, Thomas J.
2015-05-04
This paper summarizes a number of perspectives that have emerged from the sustainability sciences in recent decades (NRC, 1999; Kates et al., 2001; NRC, 2006; Kates, 2010) that shed light on the role of place in multi-scale sustainability science and vice-versa, ranging from the importance of the co-production of knowledge for sustainable development to threats to a ''sense of place'' from global environmental and economic changes.
Computational Design of Novel Multiscale Concrete Rheometers | Argonne
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Leadership Computing Facility Suspended particles in a rheometer This simulation image shows suspended particles in a rheometer for NIST's proposed mortar SRM. The spheres, which are color coded by their starting location in the rheometer, are suspended in a cement paste with properties derived from NIST's cement paste SRM. Nicos Martys and Steven G. Satterfield, National Institute of Standards and Technology Computational Design of Novel Multiscale Concrete Rheometers PI Name: William
Multi-Scale Characterization of Improved Algae Strains
Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site
Slide 1 DOE Bioenergy Technologies Office (BETO) 2015 Project Peer Review Multi-Scale Characterization of Improved Algae Strains March 23, 2015 Algae Technology Area Review Dr. Taraka Dale Los Alamos National Laboratory This presentation does not contain any proprietary, confidential, or otherwise restricted information LA-UR-15-21927 Operated by Los Alamos National Security, LLC for the U.S. Department of Energy's NNSA UNCLASSIFIED Goal Statement The overall goal of this project is to develop a
Lifetime statistics of quantum chaos studied by a multiscale analysis
Di Falco, A.; Krauss, T. F. [School of Physics and Astronomy, University of St. Andrews, North Haugh, St. Andrews, KY16 9SS (United Kingdom); Fratalocchi, A. [PRIMALIGHT, Faculty of Electrical Engineering, Applied Mathematics and Computational Science, King Abdullah University of Science and Technology (KAUST), Thuwal 23955-6900 (Saudi Arabia)
2012-04-30
In a series of pump and probe experiments, we study the lifetime statistics of a quantum chaotic resonator when the number of open channels is greater than one. Our design embeds a stadium billiard into a two dimensional photonic crystal realized on a silicon-on-insulator substrate. We calculate resonances through a multiscale procedure that combines energy landscape analysis and wavelet transforms. Experimental data is found to follow the universal predictions arising from random matrix theory with an excellent level of agreement.
Method of modeling transmissions for real-time simulation
Hebbale, Kumaraswamy V.
2012-09-25
A transmission modeling system includes an in-gear module that determines an in-gear acceleration when a vehicle is in gear. A shift module determines a shift acceleration based on a clutch torque when the vehicle is shifting between gears. A shaft acceleration determination module determines a shaft acceleration based on at least one of the in-gear acceleration and the shift acceleration.
System and method for modeling and analyzing complex scenarios
Shevitz, Daniel Wolf
2013-04-09
An embodiment of the present invention includes a method for analyzing and solving possibility tree. A possibility tree having a plurality of programmable nodes is constructed and solved with a solver module executed by a processor element. The solver module executes the programming of said nodes, and tracks the state of at least a variable through a branch. When a variable of said branch is out of tolerance with a parameter, the solver disables remaining nodes of the branch and marks the branch as an invalid solution. The valid solutions are then aggregated and displayed as valid tree solutions.
Review of Wind Energy Forecasting Methods for Modeling Ramping Events
Wharton, S; Lundquist, J K; Marjanovic, N; Williams, J L; Rhodes, M; Chow, T K; Maxwell, R
2011-03-28
Tall onshore wind turbines, with hub heights between 80 m and 100 m, can extract large amounts of energy from the atmosphere since they generally encounter higher wind speeds, but they face challenges given the complexity of boundary layer flows. This complexity of the lowest layers of the atmosphere, where wind turbines reside, has made conventional modeling efforts less than ideal. To meet the nation's goal of increasing wind power into the U.S. electrical grid, the accuracy of wind power forecasts must be improved. In this report, the Lawrence Livermore National Laboratory, in collaboration with the University of Colorado at Boulder, University of California at Berkeley, and Colorado School of Mines, evaluates innovative approaches to forecasting sudden changes in wind speed or 'ramping events' at an onshore, multimegawatt wind farm. The forecast simulations are compared to observations of wind speed and direction from tall meteorological towers and a remote-sensing Sound Detection and Ranging (SODAR) instrument. Ramping events, i.e., sudden increases or decreases in wind speed and hence, power generated by a turbine, are especially problematic for wind farm operators. Sudden changes in wind speed or direction can lead to large power generation differences across a wind farm and are very difficult to predict with current forecasting tools. Here, we quantify the ability of three models, mesoscale WRF, WRF-LES, and PF.WRF, which vary in sophistication and required user expertise, to predict three ramping events at a North American wind farm.
Search Method for Real-time Knowledge Discovery Modeled on the...
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Search Method for Real-time Knowledge Discovery Modeled on the Human Brain Oak Ridge ... information processing properties of the human brain to computational knowledge discovery. ...
MREG V1.1 : a multi-scale image registration algorithm for SAR applications.
Eichel, Paul H.
2013-08-01
MREG V1.1 is the sixth generation SAR image registration algorithm developed by the Signal Processing&Technology Department for Synthetic Aperture Radar applications. Like its predecessor algorithm REGI, it employs a powerful iterative multi-scale paradigm to achieve the competing goals of sub-pixel registration accuracy and the ability to handle large initial offsets. Since it is not model based, it allows for high fidelity tracking of spatially varying terrain-induced misregistration. Since it does not rely on image domain phase, it is equally adept at coherent and noncoherent image registration. This document provides a brief history of the registration processors developed by Dept. 5962 leading up to MREG V1.1, a full description of the signal processing steps involved in the algorithm, and a user's manual with application specific recommendations for CCD, TwoColor MultiView, and SAR stereoscopy.
MULTI-SCALE MORPHOLOGICAL ANALYSIS OF SDSS DR5 SURVEY USING THE METRIC SPACE TECHNIQUE
Wu Yongfeng; Batuski, David J.; Khalil, Andre
2009-12-20
Following the novel development and adaptation of the Metric Space Technique (MST), a multi-scale morphological analysis of the Sloan Digital Sky Survey (SDSS) Data Release 5 (DR5) was performed. The technique was adapted to perform a space-scale morphological analysis by filtering the galaxy point distributions with a smoothing Gaussian function, thus giving quantitative structural information on all size scales between 5 and 250 Mpc. The analysis was performed on a dozen slices of a volume of space containing many newly measured galaxies from the SDSS DR5 survey. Using the MST, observational data were compared to galaxy samples taken from N-body simulations with current best estimates of cosmological parameters and from random catalogs. By using the maximal ranking method among MST output functions, we also develop a way to quantify the overall similarity of the observed samples with the simulated samples.
Yortsos, Y.C.
2001-05-29
This report is an investigation of various multi-phase and multiscale transport and reaction processes associated with heavy oil recovery. The thrust areas of the project include the following: Internal drives, vapor-liquid flows, combustion and reaction processes, fluid displacements and the effect of instabilities and heterogeneities and the flow of fluids with yield stress. These find respective applications in foamy oils, the evolution of dissolved gas, internal steam drives, the mechanics of concurrent and countercurrent vapor-liquid flows, associated with thermal methods and steam injection, such as SAGD, the in-situ combustion, the upscaling of displacements in heterogeneous media and the flow of foams, Bingham plastics and heavy oils in porous media and the development of wormholes during cold production.
Yortsos, Yanis C.
2001-08-07
This project is an investigation of various multi-phase and multiscale transport and reaction processes associated with heavy oil recovery. The thrust areas of the project include the following: Internal drives, vapor-liquid flows, combustion and reaction processes, fluid displacements and the effect of instabilities and heterogeneities and the flow of fluids with yield stress. These find respective applications in foamy oils, the evolution of dissolved gas, internal steam drives, the mechanics of concurrent and countercurrent vapor-liquid flows, associated with thermal methods and steam injection, such as SAGD, the in-situ combustion, the upscaling of displacements in heterogeneous media and the flow of foams, Bingham plastics and heavy oils in porous media and the development of wormholes during cold production.
Robertson, Eric P; Christiansen, Richard L.
2007-05-29
A method of optically determining a change in magnitude of at least one dimensional characteristic of a sample in response to a selected chamber environment. A magnitude of at least one dimension of the at least one sample may be optically determined subsequent to altering the at least one environmental condition within the chamber. A maximum change in dimension of the at least one sample may be predicted. A dimensional measurement apparatus for indicating a change in at least one dimension of at least one sample. The dimensional measurement apparatus may include a housing with a chamber configured for accommodating pressure changes and an optical perception device for measuring a dimension of at least one sample disposed in the chamber. Methods of simulating injection of a gas into a subterranean formation, injecting gas into a subterranean formation, and producing methane from a coal bed are also disclosed.
Robertson, Eric P; Christiansen, Richard L.
2007-10-23
A method of optically determining a change in magnitude of at least one dimensional characteristic of a sample in response to a selected chamber environment. A magnitude of at least one dimension of the at least one sample may be optically determined subsequent to altering the at least one environmental condition within the chamber. A maximum change in dimension of the at least one sample may be predicted. A dimensional measurement apparatus for indicating a change in at least one dimension of at least one sample. The dimensional measurement apparatus may include a housing with a chamber configured for accommodating pressure changes and an optical perception device for measuring a dimension of at least one sample disposed in the chamber. Methods of simulating injection of a gas into a subterranean formation, injecting gas into a subterranean formation, and producing methane from a coal bed are also disclosed.
An improved multiscale model for dilute turbulent gas particle...
Office of Scientific and Technical Information (OSTI)
as the scale-up in the design of circulating fluidized combustor and coal gasifications. ... Country of Publication: United States Language: English Subject: 01 COAL, LIGNITE, AND ...
Multiscale model of heat dissipation mechanisms during field...
Office of Scientific and Technical Information (OSTI)
OSTI Identifier: 22492702 Resource Type: Journal Article Resource Relation: Journal Name: Applied Physics Letters; Journal Volume: 108; Journal Issue: 3; Other Information: (c) ...
Utilizing CLASIC observations and multiscale models to study...
Office of Scientific and Technical Information (OSTI)
... Close Cite: Bibtex Format Close 0 pages in this document matching the terms "" Search For Terms: Enter terms in the toolbar above to search the full text of this document for ...
Richardson, John G.
2009-11-17
An impedance estimation method includes measuring three or more impedances of an object having a periphery using three or more probes coupled to the periphery. The three or more impedance measurements are made at a first frequency. Three or more additional impedance measurements of the object are made using the three or more probes. The three or more additional impedance measurements are made at a second frequency different from the first frequency. An impedance of the object at a point within the periphery is estimated based on the impedance measurements and the additional impedance measurements.
Radiation Damage in Nuclear Fuel for Advanced Burner Reactors: Modeling and Experimental Validation
Jensen, Niels Gronbech; Asta, Mark; Ozolins, Nigel Browning'Vidvuds; de Walle, Axel van; Wolverton, Christopher
2011-12-29
The consortium has completed its existence and we are here highlighting work and accomplishments. As outlined in the proposal, the objective of the work was to advance the theoretical understanding of advanced nuclear fuel materials (oxides) toward a comprehensive modeling strategy that incorporates the different relevant scales involved in radiation damage in oxide fuels. Approaching this we set out to investigate and develop a set of directions: 1) Fission fragment and ion trajectory studies through advanced molecular dynamics methods that allow for statistical multi-scale simulations. This work also includes an investigation of appropriate interatomic force fields useful for the energetic multi-scale phenomena of high energy collisions; 2) Studies of defect and gas bubble formation through electronic structure and Monte Carlo simulations; and 3) an experimental component for the characterization of materials such that comparisons can be obtained between theory and experiment.
A multilingual programming model for coupled systems.
Ong, E. T.; Larson, J. W.; Norris, B.; Tobis, M.; Steder, M.; Jacob, R. L.; Mathematics and Computer Science; Univ. of Wisconsin; Univ. of Chicago; The Australian National Univ.
2008-01-01
Multiphysics and multiscale simulation systems share a common software requirement-infrastructure to implement data exchanges between their constituent parts-often called the coupling problem. On distributed-memory parallel platforms, the coupling problem is complicated by the need to describe, transfer, and transform distributed data, known as the parallel coupling problem. Parallel coupling is emerging as a new grand challenge in computational science as scientists attempt to build multiscale and multiphysics systems on parallel platforms. An additional coupling problem in these systems is language interoperability between their constituent codes. We have created a multilingual parallel coupling programming model based on a successful open-source parallel coupling library, the Model Coupling Toolkit (MCT). This programming model's capabilities reach beyond MCT's native Fortran implementation to include bindings for the C++ and Python programming languages. We describe the method used to generate the interlanguage bindings. This approach enables an object-based programming model for implementing parallel couplings in non-Fortran coupled systems and in systems with language heterogeneity. We describe the C++ and Python versions of the MCT programming model and provide short examples. We report preliminary performance results for the MCT interpolation benchmark. We describe a major Python application that uses the MCT Python bindings, a Python implementation of the control and coupling infrastructure for the community climate system model. We conclude with a discussion of the significance of this work to productivity computing in multidisciplinary computational science.
Three-Dimensional Lithium-Ion Battery Model (Presentation)
Kim, G. H.; Smith, K.
2008-05-01
Nonuniform battery physics can cause unexpected performance and life degradations in lithium-ion batteries; a three-dimensional cell performance model was developed by integrating an electrode-scale submodel using a multiscale modeling scheme.
Multiscale Simulation of Moist Global Atmospheric Flows
Grabowski, Wojciech W.; Smolarkiewicz, P. K.
2015-04-13
The overarching goal of this award was to include phase changes of the water substance and accompanying latent heating and precipitation processes into the all-scale nonhydrostatic atmospheric dynamics EUlerian/LAGrangian (EULAG) model. The model includes fluid flow solver that is based on either an unabbreviated set of the governing equations (i.e., compressible dynamics) or a simplified set of equations without sound waves (i.e., sound-proof, either anelastic or pseudo-incompressible). The latter set has been used in small-scale dynamics for decades, but its application to the all-scale dynamics (from small-scale to planetary) has never been studied in practical implementations. The highlight of the project is the development of the moist implicit compressible model that can be run by applying time steps, as long as the anelastic model is limited only by the computational stability of the fluid flow and not by the speed of sound waves that limit the stability of explicit compressible models. Applying various versions of the EULAG model within the same numerical framework allows for an unprecedented comparison of solutions obtained with various sets of the governing equations and straightforward evaluation of the impact of various physical parameterizations on the model solutions. The main outcomes of this study are reported in three papers, two published and one currently under review. These papers include comparisons between model solutions for idealized moist problems across the range of scales from small to planetary. These tests include: moist thermals rising in the stable-stratified environment (following Grabowski and Clark, J. Atmos. Sci. 1991) and in the moist-neutral environment (after Bryan and Fritsch, Mon. Wea. Rev. 2002), moist flows over a mesoscale topography (as in Grabowski and Smolarkiewicz, Mon. Wea. Rev. 2002), deep convection in a sheared environment (following Weisman and Klemp, Mon. Wea. Rev. 1982), moist extension of the baroclinic wave on
Multi-scale statistical analysis of coronal solar activity
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Gamborino, Diana; del-Castillo-Negrete, Diego; Martinell, Julio J.
2016-07-08
Multi-filter images from the solar corona are used to obtain temperature maps that are analyzed using techniques based on proper orthogonal decomposition (POD) in order to extract dynamical and structural information at various scales. Exploring active regions before and after a solar flare and comparing them with quiet regions, we show that the multi-scale behavior presents distinct statistical properties for each case that can be used to characterize the level of activity in a region. Information about the nature of heat transport is also to be extracted from the analysis.
Synthesis of Numerical Methods for Modeling Wave Energy Converter-Point Absorbers: Preprint
Li, Y.; Yu, Y. H.
2012-05-01
During the past few decades, wave energy has received significant attention among all ocean energy formats. Industry has proposed hundreds of prototypes such as an oscillating water column, a point absorber, an overtopping system, and a bottom-hinged system. In particular, many researchers have focused on modeling the floating-point absorber as the technology to extract wave energy. Several modeling methods have been used such as the analytical method, the boundary-integral equation method, the Navier-Stokes equations method, and the empirical method. However, no standardized method has been decided. To assist the development of wave energy conversion technologies, this report reviews the methods for modeling the floating-point absorber.
Proposed SPAR Modeling Method for Quantifying Time Dependent Station Blackout Cut Sets
John A. Schroeder
2010-06-01
Abstract: The U.S. Nuclear Regulatory Commissions (USNRCs) Standardized Plant Analysis Risk (SPAR) models and industry risk models take similar approaches to analyzing the risk associated with loss of offsite power and station blackout (LOOP/SBO) events at nuclear reactor plants. In both SPAR models and industry models, core damage risk resulting from a LOOP/SBO event is analyzed using a combination of event trees and fault trees that produce cut sets that are, in turn, quantified to obtain a numerical estimate of the resulting core damage risk. A proposed SPAR method for quantifying the time-dependent cut sets is sometimes referred to as a convolution method. The SPAR method reflects assumptions about the timing of emergency diesel failures, the timing of subsequent attempts at emergency diesel repair, and the timing of core damage that may be different than those often used in industry models. This paper describes the proposed SPAR method.
Methods to Register Models and Input/Output Parameters for Integrated Modeling
Droppo, James G.; Whelan, Gene; Tryby, Michael E.; Pelton, Mitchell A.; Taira, Randal Y.; Dorow, Kevin E.
2010-07-10
Significant resources can be required when constructing integrated modeling systems. In a typical application, components (e.g., models and databases) created by different developers are assimilated, requiring the framework’s functionality to bridge the gap between the user’s knowledge of the components being linked. The framework, therefore, needs the capability to assimilate a wide range of model-specific input/output requirements as well as their associated assumptions and constraints. The process of assimilating such disparate components into an integrated modeling framework varies in complexity and difficulty. Several factors influence the relative ease of assimilating components, including, but not limited to, familiarity with the components being assimilated, familiarity with the framework and its tools that support the assimilation process, level of documentation associated with the components and the framework, and design structure of the components and framework. This initial effort reviews different approaches for assimilating models and their model-specific input/output requirements: 1) modifying component models to directly communicate with the framework (i.e., through an Application Programming Interface), 2) developing model-specific external wrappers such that no component model modifications are required, 3) using parsing tools to visually map pre-existing input/output files, and 4) describing and linking models as dynamic link libraries. Most of these approaches are illustrated using the widely distributed modeling system called Framework for Risk Analysis in Multimedia Environmental Systems (FRAMES). The review concludes that each has its strengths and weakness, the factors that determine which approaches work best in a given application.
Sereda, Yuriy V.; Ortoleva, Peter J.
2014-04-07
A closed kinetic equation for the single-particle density of a viscous simple liquid is derived using a variational method for the Liouville equation and a coarse-grained mean-field (CGMF) ansatz. The CGMF ansatz is based on the notion that during the characteristic time of deformation a given particle interacts with many others so that it experiences an average interaction. A trial function for the N-particle probability density is constructed using a multiscale perturbation method and the CGMF ansatz is applied to it. The multiscale perturbation scheme is based on the ratio of the average nearest-neighbor atom distance to the total size of the assembly. A constraint on the initial condition is discovered which guarantees that the kinetic equation is mass-conserving and closed in the single-particle density. The kinetic equation has much of the character of the Vlasov equation except that true viscous, and not Landau, damping is accounted for. The theory captures condensation kinetics and takes much of the character of the Gross-Pitaevskii equation in the weak-gradient short-range force limit.
Physics-based statistical model and simulation method of RF propagation in urban environments
Pao, Hsueh-Yuan; Dvorak, Steven L.
2010-09-14
A physics-based statistical model and simulation/modeling method and system of electromagnetic wave propagation (wireless communication) in urban environments. In particular, the model is a computationally efficient close-formed parametric model of RF propagation in an urban environment which is extracted from a physics-based statistical wireless channel simulation method and system. The simulation divides the complex urban environment into a network of interconnected urban canyon waveguides which can be analyzed individually; calculates spectral coefficients of modal fields in the waveguides excited by the propagation using a database of statistical impedance boundary conditions which incorporates the complexity of building walls in the propagation model; determines statistical parameters of the calculated modal fields; and determines a parametric propagation model based on the statistical parameters of the calculated modal fields from which predictions of communications capability may be made.
Mixed-RKDG Finite Element Methods for the 2-D Hydrodynamic Model for Semiconductor Device Simulation
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Chen, Zhangxin; Cockburn, Bernardo; Jerome, Joseph W.; Shu, Chi-Wang
1995-01-01
In this paper we introduce a new method for numerically solving the equations of the hydrodynamic model for semiconductor devices in two space dimensions. The method combines a standard mixed finite element method, used to obtain directly an approximation to the electric field, with the so-called Runge-Kutta Discontinuous Galerkin (RKDG) method, originally devised for numerically solving multi-dimensional hyperbolic systems of conservation laws, which is applied here to the convective part of the equations. Numerical simulations showing the performance of the new method are displayed, and the results compared with those obtained by using Essentially Nonoscillatory (ENO) finite difference schemes. Frommore » the perspective of device modeling, these methods are robust, since they are capable of encompassing broad parameter ranges, including those for which shock formation is possible. The simulations presented here are for Gallium Arsenide at room temperature, but we have tested them much more generally with considerable success.« less
Gao Yajun
2008-08-15
A previously established Hauser-Ernst-type extended double-complex linear system is slightly modified and used to develop an inverse scattering method for the stationary axisymmetric general symplectic gravity model. The reduction procedures in this inverse scattering method are found to be fairly simple, which makes the inverse scattering method applied fine and effective. As an application, a concrete family of soliton double solutions for the considered theory is obtained.
A Comparison of Multiscale Variations of Decade-long Cloud Fractions...
Office of Scientific and Technical Information (OSTI)
A Comparison of Multiscale Variations of Decade-long Cloud Fractions from Six Different Platforms over the Southern Great Plains in the United States Citation Details In-Document ...