Multiscale vision model for event detection and reconstruction in two-photon imaging data

DEFF Research Database (Denmark)

Brazhe, Alexey; Mathiesen, Claus; Lind, Barbara Lykke

2014-01-01

on a modified multiscale vision model, an object detection framework based on the thresholding of wavelet coefficients and hierarchical trees of significant coefficients followed by nonlinear iterative partial object reconstruction, for the analysis of two-photon calcium imaging data. The framework is discussed...... of the multiscale vision model is similar in the denoising, but provides a better segmenation of the image into meaningful objects, whereas other methods need to be combined with dedicated thresholding and segmentation utilities....

Computer Aided Multi-scale Design of SiC-Si3N4 Nanoceramic Composites for High-Temperature Structural Applications

Energy Technology Data Exchange (ETDEWEB)

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

Multiscale geometric modeling of macromolecules I: Cartesian representation

Science.gov (United States)

Xia, Kelin; Feng, Xin; Chen, Zhan; Tong, Yiying; Wei, Guo-Wei

2014-01-01

This paper focuses on the geometric modeling and computational algorithm development of biomolecular structures from two data sources: Protein Data Bank (PDB) and Electron Microscopy Data Bank (EMDB) in the Eulerian (or Cartesian) representation. Molecular surface (MS) contains non-smooth geometric singularities, such as cusps, tips and self-intersecting facets, which often lead to computational instabilities in molecular simulations, and violate the physical principle of surface free energy minimization. Variational multiscale surface definitions are proposed based on geometric flows and solvation analysis of biomolecular systems. Our approach leads to geometric and potential driven Laplace-Beltrami flows for biomolecular surface evolution and formation. The resulting surfaces are free of geometric singularities and minimize the total free energy of the biomolecular system. High order partial differential equation (PDE)-based nonlinear filters are employed for EMDB data processing. We show the efficacy of this approach in feature-preserving noise reduction. After the construction of protein multiresolution surfaces, we explore the analysis and characterization of surface morphology by using a variety of curvature definitions. Apart from the classical Gaussian curvature and mean curvature, maximum curvature, minimum curvature, shape index, and curvedness are also applied to macromolecular surface analysis for the first time. Our curvature analysis is uniquely coupled to the analysis of electrostatic surface potential, which is a by-product of our variational multiscale solvation models. As an expository investigation, we particularly emphasize the numerical algorithms and computational protocols for practical applications of the above multiscale geometric models. Such information may otherwise be scattered over the vast literature on this topic. Based on the curvature and electrostatic analysis from our multiresolution surfaces, we introduce a new concept, the

Multiscale geometric modeling of macromolecules I: Cartesian representation

Energy Technology Data Exchange (ETDEWEB)

Xia, Kelin [Department of Mathematics, Michigan State University, MI 48824 (United States); Feng, Xin [Department of Computer Science and Engineering, Michigan State University, MI 48824 (United States); Chen, Zhan [Department of Mathematics, Michigan State University, MI 48824 (United States); Tong, Yiying [Department of Computer Science and Engineering, Michigan State University, MI 48824 (United States); Wei, Guo-Wei, E-mail: wei@math.msu.edu [Department of Mathematics, Michigan State University, MI 48824 (United States); Department of Biochemistry and Molecular Biology, Michigan State University, MI 48824 (United States)

2014-01-15

This paper focuses on the geometric modeling and computational algorithm development of biomolecular structures from two data sources: Protein Data Bank (PDB) and Electron Microscopy Data Bank (EMDB) in the Eulerian (or Cartesian) representation. Molecular surface (MS) contains non-smooth geometric singularities, such as cusps, tips and self-intersecting facets, which often lead to computational instabilities in molecular simulations, and violate the physical principle of surface free energy minimization. Variational multiscale surface definitions are proposed based on geometric flows and solvation analysis of biomolecular systems. Our approach leads to geometric and potential driven Laplace–Beltrami flows for biomolecular surface evolution and formation. The resulting surfaces are free of geometric singularities and minimize the total free energy of the biomolecular system. High order partial differential equation (PDE)-based nonlinear filters are employed for EMDB data processing. We show the efficacy of this approach in feature-preserving noise reduction. After the construction of protein multiresolution surfaces, we explore the analysis and characterization of surface morphology by using a variety of curvature definitions. Apart from the classical Gaussian curvature and mean curvature, maximum curvature, minimum curvature, shape index, and curvedness are also applied to macromolecular surface analysis for the first time. Our curvature analysis is uniquely coupled to the analysis of electrostatic surface potential, which is a by-product of our variational multiscale solvation models. As an expository investigation, we particularly emphasize the numerical algorithms and computational protocols for practical applications of the above multiscale geometric models. Such information may otherwise be scattered over the vast literature on this topic. Based on the curvature and electrostatic analysis from our multiresolution surfaces, we introduce a new concept, the

Use of ARM Data to address the Climate Change Further Development and Applications of A Multi-scale Modeling Framework

Energy Technology Data Exchange (ETDEWEB)

David A. Randall; Marat Khairoutdinov

2007-12-14

The Colorado State University (CSU) Multi-scale Modeling Framework (MMF) is a new type of general circulation model (GCM) that replaces the conventional parameterizations of convection, clouds and boundary layer with a cloud-resolving model (CRM) embedded into each grid column. The MMF that we have been working with is a “super-parameterized” version of the Community Atmosphere Model (CAM). As reported in the publications listed below, we have done extensive work with the model. We have explored the MMF’s performance in several studies, including an AMIP run and a CAPT test, and we have applied the MMF to an analysis of climate sensitivity.

Multi-scale graph-cut algorithm for efficient water-fat separation.

Science.gov (United States)

Berglund, Johan; Skorpil, Mikael

2017-09-01

To improve the accuracy and robustness to noise in water-fat separation by unifying the multiscale and graph cut based approaches to B 0 -correction. A previously proposed water-fat separation algorithm that corrects for B 0 field inhomogeneity in 3D by a single quadratic pseudo-Boolean optimization (QPBO) graph cut was incorporated into a multi-scale framework, where field map solutions are propagated from coarse to fine scales for voxels that are not resolved by the graph cut. The accuracy of the single-scale and multi-scale QPBO algorithms was evaluated against benchmark reference datasets. The robustness to noise was evaluated by adding noise to the input data prior to water-fat separation. Both algorithms achieved the highest accuracy when compared with seven previously published methods, while computation times were acceptable for implementation in clinical routine. The multi-scale algorithm was more robust to noise than the single-scale algorithm, while causing only a small increase (+10%) of the reconstruction time. The proposed 3D multi-scale QPBO algorithm offers accurate water-fat separation, robustness to noise, and fast reconstruction. The software implementation is freely available to the research community. Magn Reson Med 78:941-949, 2017. © 2016 International Society for Magnetic Resonance in Medicine. © 2016 International Society for Magnetic Resonance in Medicine.

Multiscale Drivers of Global Environmental Health

Science.gov (United States)

Desai, Manish Anil

In this dissertation, I motivate, develop, and demonstrate three such approaches for investigating multiscale drivers of global environmental health: (1) a metric for analyzing contributions and responses to climate change from global to sectoral scales, (2) a framework for unraveling the influence of environmental change on infectious diseases at regional to local scales, and (3) a model for informing the design and evaluation of clean cooking interventions at community to household scales. The full utility of climate debt as an analytical perspective will remain untapped without tools that can be manipulated by a wide range of analysts, including global environmental health researchers. Chapter 2 explains how international natural debt (IND) apportions global radiative forcing from fossil fuel carbon dioxide and methane, the two most significant climate altering pollutants, to individual entities -- primarily countries but also subnational states and economic sectors, with even finer scales possible -- as a function of unique trajectories of historical emissions, taking into account the quite different radiative efficiencies and atmospheric lifetimes of each pollutant. Owing to its straightforward and transparent derivation, IND can readily operationalize climate debt to consider issues of equity and efficiency and drive scenario exercises that explore the response to climate change at multiple scales. Collectively, the analyses presented in this chapter demonstrate how IND can inform a range of key question on climate change mitigation at multiple scales, compelling environmental health towards an appraisal of the causes and not just the consequences of climate change. The environmental change and infectious disease (EnvID) conceptual framework of Chapter 3 builds on a rich history of prior efforts in epidemiologic theory, environmental science, and mathematical modeling by: (1) articulating a flexible and logical system specification; (2) incorporating

Anatomy and Physiology of Multiscale Modeling and Simulation in Systems Medicine.

Science.gov (United States)

Mizeranschi, Alexandru; Groen, Derek; Borgdorff, Joris; Hoekstra, Alfons G; Chopard, Bastien; Dubitzky, Werner

2016-01-01

Systems medicine is the application of systems biology concepts, methods, and tools to medical research and practice. It aims to integrate data and knowledge from different disciplines into biomedical models and simulations for the understanding, prevention, cure, and management of complex diseases. Complex diseases arise from the interactions among disease-influencing factors across multiple levels of biological organization from the environment to molecules. To tackle the enormous challenges posed by complex diseases, we need a modeling and simulation framework capable of capturing and integrating information originating from multiple spatiotemporal and organizational scales. Multiscale modeling and simulation in systems medicine is an emerging methodology and discipline that has already demonstrated its potential in becoming this framework. The aim of this chapter is to present some of the main concepts, requirements, and challenges of multiscale modeling and simulation in systems medicine.

Multiscale modeling of pedestrian dynamics

CERN Document Server

Cristiani, Emiliano; Tosin, Andrea

2014-01-01

This book presents mathematical models and numerical simulations of crowd dynamics. The core topic is the development of a new multiscale paradigm, which bridges the microscopic and macroscopic scales taking the most from each of them for capturing the relevant clues of complexity of crowds. The background idea is indeed that most of the complex trends exhibited by crowds are due to an intrinsic interplay between individual and collective behaviors. The modeling approach promoted in this book pursues actively this intuition and profits from it for designing general mathematical structures susceptible of application also in fields different from the inspiring original one. The book considers also the two most traditional points of view: the microscopic one, in which pedestrians are tracked individually, and the macroscopic one, in which pedestrians are assimilated to a continuum. Selected existing models are critically analyzed. The work is addressed to researchers and graduate students.

Multiscale analysis and computation for flows in heterogeneous media

Energy Technology Data Exchange (ETDEWEB)

Efendiev, Yalchin [Texas A & M Univ., College Station, TX (United States); Hou, T. Y. [California Inst. of Technology (CalTech), Pasadena, CA (United States); Durlofsky, L. J. [Stanford Univ., CA (United States); Tchelepi, H. [Stanford Univ., CA (United States)

2016-08-04

Our work in this project is aimed at making fundamental advances in multiscale methods for flow and transport in highly heterogeneous porous media. The main thrust of this research is to develop a systematic multiscale analysis and efficient coarse-scale models that can capture global effects and extend existing multiscale approaches to problems with additional physics and uncertainties. A key emphasis is on problems without an apparent scale separation. Multiscale solution methods are currently under active investigation for the simulation of subsurface flow in heterogeneous formations. These procedures capture the effects of fine-scale permeability variations through the calculation of specialized coarse-scale basis functions. Most of the multiscale techniques presented to date employ localization approximations in the calculation of these basis functions. For some highly correlated (e.g., channelized) formations, however, global effects are important and these may need to be incorporated into the multiscale basis functions. Other challenging issues facing multiscale simulations are the extension of existing multiscale techniques to problems with additional physics, such as compressibility, capillary effects, etc. In our project, we explore the improvement of multiscale methods through the incorporation of additional (single-phase flow) information and the development of a general multiscale framework for flows in the presence of uncertainties, compressible flow and heterogeneous transport, and geomechanics. We have considered (1) adaptive local-global multiscale methods, (2) multiscale methods for the transport equation, (3) operator-based multiscale methods and solvers, (4) multiscale methods in the presence of uncertainties and applications, (5) multiscale finite element methods for high contrast porous media and their generalizations, and (6) multiscale methods for geomechanics. Below, we present a brief overview of each of these contributions.

Multi-scale modeling of dispersed gas-liquid two-phase flow

NARCIS (Netherlands)

Deen, N.G.; Sint Annaland, van M.; Kuipers, J.A.M.

2004-01-01

In this work the concept of multi-scale modeling is demonstrated. The idea of this approach is to use different levels of modeling, each developed to study phenomena at a certain length scale. Information obtained at the level of small length scales can be used to provide closure information at the

Mechanical Properties of Graphene Nanoplatelet/Carbon Fiber/Epoxy Hybrid Composites: Multiscale Modeling and Experiments

Science.gov (United States)

Hadden, C. M.; Klimek-McDonald, D. R.; Pineda, E. J.; King, J. A.; Reichanadter, A. M.; Miskioglu, I.; Gowtham, S.; Odegard, G. M.

2015-01-01

Because of the relatively high specific mechanical properties of carbon fiber/epoxy composite materials, they are often used as structural components in aerospace applications. Graphene nanoplatelets (GNPs) can be added to the epoxy matrix to improve the overall mechanical properties of the composite. The resulting GNP/carbon fiber/epoxy hybrid composites have been studied using multiscale modeling to determine the influence of GNP volume fraction, epoxy crosslink density, and GNP dispersion on the mechanical performance. The hierarchical multiscale modeling approach developed herein includes Molecular Dynamics (MD) and micromechanical modeling, and it is validated with experimental testing of the same hybrid composite material system. The results indicate that the multiscale modeling approach is accurate and provides physical insight into the composite mechanical behavior. Also, the results quantify the substantial impact of GNP volume fraction and dispersion on the transverse mechanical properties of the hybrid composite, while the effect on the axial properties is shown to be insignificant.

Mechanical Properties of Graphene Nanoplatelet Carbon Fiber Epoxy Hybrid Composites: Multiscale Modeling and Experiments

Science.gov (United States)

Hadden, Cameron M.; Klimek-McDonald, Danielle R.; Pineda, Evan J.; King, Julie A.; Reichanadter, Alex M.; Miskioglu, Ibrahim; Gowtham, S.; Odegard, Gregory M.

2015-01-01

Because of the relatively high specific mechanical properties of carbon fiber/epoxy composite materials, they are often used as structural components in aerospace applications. Graphene nanoplatelets (GNPs) can be added to the epoxy matrix to improve the overall mechanical properties of the composite. The resulting GNP/carbon fiber/epoxy hybrid composites have been studied using multiscale modeling to determine the influence of GNP volume fraction, epoxy crosslink density, and GNP dispersion on the mechanical performance. The hierarchical multiscale modeling approach developed herein includes Molecular Dynamics (MD) and micromechanical modeling, and it is validated with experimental testing of the same hybrid composite material system. The results indicate that the multiscale modeling approach is accurate and provides physical insight into the composite mechanical behavior. Also, the results quantify the substantial impact of GNP volume fraction and dispersion on the transverse mechanical properties of the hybrid composite, while the effect on the axial properties is shown to be insignificant.

Hybrid continuum–molecular modelling of multiscale internal gas flows

International Nuclear Information System (INIS)

Patronis, Alexander; Lockerby, Duncan A.; Borg, Matthew K.; Reese, Jason M.

2013-01-01

We develop and apply an efficient multiscale method for simulating a large class of low-speed internal rarefied gas flows. The method is an extension of the hybrid atomistic–continuum approach proposed by Borg et al. (2013) [28] for the simulation of micro/nano flows of high-aspect ratio. The major new extensions are: (1) incorporation of fluid compressibility; (2) implementation using the direct simulation Monte Carlo (DSMC) method for dilute rarefied gas flows, and (3) application to a broader range of geometries, including periodic, non-periodic, pressure-driven, gravity-driven and shear-driven internal flows. The multiscale method is applied to micro-scale gas flows through a periodic converging–diverging channel (driven by an external acceleration) and a non-periodic channel with a bend (driven by a pressure difference), as well as the flow between two eccentric cylinders (with the inner rotating relative to the outer). In all these cases there exists a wide variation of Knudsen number within the geometries, as well as substantial compressibility despite the Mach number being very low. For validation purposes, our multiscale simulation results are compared to those obtained from full-scale DSMC simulations: very close agreement is obtained in all cases for all flow variables considered. Our multiscale simulation is an order of magnitude more computationally efficient than the full-scale DSMC for the first and second test cases, and two orders of magnitude more efficient for the third case

Multi-scale associations between vegetation cover and woodland bird communities across a large agricultural region.

Directory of Open Access Journals (Sweden)

Karen Ikin

Full Text Available Improving biodiversity conservation in fragmented agricultural landscapes has become an important global issue. Vegetation at the patch and landscape-scale is important for species occupancy and diversity, yet few previous studies have explored multi-scale associations between vegetation and community assemblages. Here, we investigated how patch and landscape-scale vegetation cover structure woodland bird communities. We asked: (1 How is the bird community associated with the vegetation structure of woodland patches and the amount of vegetation cover in the surrounding landscape? (2 Do species of conservation concern respond to woodland vegetation structure and surrounding vegetation cover differently to other species in the community? And (3 Can the relationships between the bird community and the woodland vegetation structure and surrounding vegetation cover be explained by the ecological traits of the species comprising the bird community? We studied 103 woodland patches (0.5 - 53.8 ha over two time periods across a large (6,800 km(2 agricultural region in southeastern Australia. We found that both patch vegetation and surrounding woody vegetation cover were important for structuring the bird community, and that these relationships were consistent over time. In particular, the occurrence of mistletoe within the patches and high values of woody vegetation cover within 1,000 ha and 10,000 ha were important, especially for bird species of conservation concern. We found that the majority of these species displayed similar, positive responses to patch and landscape vegetation attributes. We also found that these relationships were related to the foraging and nesting traits of the bird community. Our findings suggest that management strategies to increase both remnant vegetation quality and the cover of surrounding woody vegetation in fragmented agricultural landscapes may lead to improved conservation of bird communities.

Implementation of dust emission and chemistry into the Community Multiscale Air Quality modeling system and initial application to an Asian dust storm episode

Directory of Open Access Journals (Sweden)

K. Wang

2012-11-01

Full Text Available The US Environmental Protection Agency's (EPA Community Multiscale Air Quality (CMAQ modeling system version 4.7 is further developed to enhance its capability in simulating the photochemical cycles in the presence of dust particles. The new model treatments implemented in CMAQ v4.7 in this work include two online dust emission schemes (i.e., the Zender and Westphal schemes, nine dust-related heterogeneous reactions, an updated aerosol inorganic thermodynamic module ISORROPIA II with an explicit treatment of crustal species, and the interface between ISORROPIA II and the new dust treatments. The resulting improved CMAQ (referred to as CMAQ-Dust, offline-coupled with the Weather Research and Forecast model (WRF, is applied to the April 2001 dust storm episode over the trans-Pacific domain to examine the impact of new model treatments and understand associated uncertainties. WRF/CMAQ-Dust produces reasonable spatial distribution of dust emissions and captures the dust outbreak events, with the total dust emissions of ~111 and 223 Tg when using the Zender scheme with an erodible fraction of 0.5 and 1.0, respectively. The model system can reproduce well observed meteorological and chemical concentrations, with significant improvements for suspended particulate matter (PM, PM with aerodynamic diameter of 10 μm, and aerosol optical depth than the default CMAQ v4.7. The sensitivity studies show that the inclusion of crustal species reduces the concentration of PM with aerodynamic diameter of 2.5 μm (PM_2.5 over polluted areas. The heterogeneous chemistry occurring on dust particles acts as a sink for some species (e.g., as a lower limit estimate, reducing O₃ by up to 3.8 ppb (~9% and SO₂ by up to 0.3 ppb (~27% and as a source for some others (e.g., increasing fine-mode SO₄²⁻ by up to 1.1 μg m⁻³ (~12% and PM_2.5 by up to 1.4 μg m⁻³ (~3% over the domain. The

Development of a Quasi-3D Multiscale Modeling Framework: Motivation, basic algorithm and preliminary results

Directory of Open Access Journals (Sweden)

Joon-Hee Jung

2010-11-01

Full Text Available A new framework for modeling the atmosphere, which we call the quasi-3D (Q3D multi-scale modeling framework (MMF, is developed with the objective of including cloud-scale three-dimensional effects in a GCM without necessarily using a global cloud-resolving model (CRM. It combines a GCM with a Q3D CRM that has the horizontal domain consisting of two perpendicular sets of channels, each of which contains a locally 3D grid-point array. For computing efficiency, the widths of the channels are chosen to be narrow. Thus, it is crucial to select a proper lateral boundary condition to realistically simulate the statistics of cloud and cloud-associated processes. Among the various possibilities, a periodic lateral boundary condition is chosen for the deviations from background fields that are obtained by interpolations from the GCM grid points. Since the deviations tend to vanish as the GCM grid size approaches that of the CRM, the whole system of the Q3D MMF can converge to a fully 3D global CRM. Consequently, the horizontal resolution of the GCM can be freely chosen depending on the objective of application, without changing the formulation of model physics. To evaluate the newly developed Q3D CRM in an efficient way, idealized experiments have been performed using a small horizontal domain. In these tests, the Q3D CRM uses only one pair of perpendicular channels with only two grid points across each channel. Comparing the simulation results with those of a fully 3D CRM, it is concluded that the Q3D CRM can reproduce most of the important statistics of the 3D solutions, including the vertical distributions of cloud water and precipitants, vertical transports of potential temperature and water vapor, and the variances and covariances of dynamical variables. The main improvement from a corresponding 2D simulation appears in the surface fluxes and the vorticity transports that cause the mean wind to change. A comparison with a simulation using a coarse

Sensitivity of the Community Multiscale Air Quality (CMAQ model v4.7 results for the eastern United States to MM5 and WRF meteorological drivers

Directory of Open Access Journals (Sweden)

K. W. Appel

2010-02-01

Full Text Available This paper presents a comparison of the operational performances of two Community Multiscale Air Quality (CMAQ model v4.7 simulations that utilize input data from the 5th-generation Mesoscale Model (MM5 and the Weather Research and Forecasting (WRF meteorological models. Two sets of CMAQ model simulations were performed for January and August 2006. One set utilized MM5 meteorology (MM5-CMAQ and the other utilized WRF meteorology (WRF-CMAQ, while all other model inputs and options were kept the same. For January, predicted ozone (O₃ mixing ratios were higher in the Southeast and lower Mid-west regions in the WRF-CMAQ simulation, resulting in slightly higher bias and error as compared to the MM5-CMAQ simulations. The higher predicted O₃ mixing ratios are attributed to less dry deposition of O₃ in the WRF-CMAQ simulation due to differences in the calculation of the vegetation fraction between the MM5 and WRF models. The WRF-CMAQ results showed better performance for particulate sulfate (SO₄²⁻, similar performance for nitrate (NO₃⁻, and slightly worse performance for nitric acid (HNO₃, total carbon (TC and total fine particulate (PM_2.5 mass than the corresponding MM5-CMAQ results. For August, predictions of O₃ were notably higher in the WRF-CMAQ simulation, particularly in the southern United States, resulting in increased model bias. Concentrations of predicted particulate SO₄²⁻ were lower in the region surrounding the Ohio Valley and higher along the Gulf of Mexico in the WRF-CMAQ simulation, contributing to poorer model performance. The primary causes of the differences in the MM5-CMAQ and WRF-CMAQ simulations appear to be due to differences in the calculation of wind speed, planetary boundary layer height, cloud cover and the friction velocity (u_∗ in the MM5 and WRF model simulations, while

Multi-Scale Modelling of the Gamma Radiolysis of Nitrate Solutions

OpenAIRE

Horne, Gregory; Donoclift, Thomas; Sims, Howard E.; M. Orr, Robin; Pimblott, Simon

2016-01-01

A multi-scale modelling approach has been developed for the extended timescale long-term radiolysis of aqueous systems. The approach uses a combination of stochastic track structure and track chemistry as well as deterministic homogeneous chemistry techniques and involves four key stages; radiation track structure simulation, the subsequent physicochemical processes, nonhomogeneous diffusion-reaction kinetic evolution, and homogeneous bulk chemistry modelling. The first three components model...

Multiscale modeling of θ' precipitation in Al-Cu binary alloys

International Nuclear Information System (INIS)

Vaithyanathan, V.; Wolverton, C.; Chen, L.Q.

2004-01-01

We present a multiscale model for studying the growth and coarsening of θ' precipitates in Al-Cu alloys. Our approach utilizes a novel combination of the mesoscale phase-field method with atomistic approaches such as first-principles total energy and linear response calculations, as well as a mixed-space cluster expansion coupled with Monte Carlo simulations. We give quantitative first-principles predictions of: (i) bulk energetics of the Al-Cu solid solution and θ ' precipitate phases, (ii) interfacial energies of the coherent and semi-coherent θ ' /Al interfaces, and (iii) stress-free misfit strains and coherency strain energies of the θ ' /Al system. These first-principles data comprise all the necessary energetic information to construct our phase-field model of microstructural evolution. Using our multiscale approach, we elucidate the effects of various energetic contributions on the equilibrium shape of θ ' precipitates, finding that both the elastic energy and interfacial energy anisotropy contributions play critical roles in determining the aspect ratio of θ ' precipitates. Additionally, we have performed a quantitative study of the morphology of two-dimensional multi-precipitate microstructures during growth and coarsening, and compared the calculated results with experimentally observed morphologies. Our multiscale first-principles/phase-field method is completely general and should therefore be applicable to a wide variety of problems in microstructural evolution

Multi-scale simulation of droplet-droplet interactions and coalescence

CSIR Research Space (South Africa)

Musehane, Ndivhuwo M

2016-10-01

Full Text Available Conference on Computational and Applied Mechanics Potchefstroom 3–5 October 2016 Multi-scale simulation of droplet-droplet interactions and coalescence 1,2Ndivhuwo M. Musehane?, 1Oliver F. Oxtoby and 2Daya B. Reddy 1. Aeronautic Systems, Council... topology changes that result when droplets interact. This work endeavours to eliminate the need to use empirical correlations based on phenomenological models by developing a multi-scale model that predicts the outcome of a collision between droplets from...

Cloud Feedbacks on Greenhouse Warming in a Multi-Scale Modeling Framework with a Higher-Order Turbulence Closure

Science.gov (United States)

Cheng, Anning; Xu, Kuan-Man

2015-01-01

Five-year simulation experiments with a multi-scale modeling Framework (MMF) with a advanced intermediately prognostic higher-order turbulence closure (IPHOC) in its cloud resolving model (CRM) component, also known as SPCAM-IPHOC (super parameterized Community Atmospheric Model), are performed to understand the fast tropical (30S-30N) cloud response to an instantaneous doubling of CO2 concentration with SST held fixed at present-day values. SPCAM-IPHOC has substantially improved the low-level representation compared with SPCAM. It is expected that the cloud responses to greenhouse warming in SPCAM-IPHOC is more realistic. The change of rising motion, surface precipitation, cloud cover, and shortwave and longwave cloud radiative forcing in SPCAM-IPHOC from the greenhouse warming will be presented in the presentation.

Development of a Sampling-Based Global Sensitivity Analysis Workflow for Multiscale Computational Cancer Models

Science.gov (United States)

Wang, Zhihui; Deisboeck, Thomas S.; Cristini, Vittorio

2014-01-01

There are two challenges that researchers face when performing global sensitivity analysis (GSA) on multiscale in silico cancer models. The first is increased computational intensity, since a multiscale cancer model generally takes longer to run than does a scale-specific model. The second problem is the lack of a best GSA method that fits all types of models, which implies that multiple methods and their sequence need to be taken into account. In this article, we therefore propose a sampling-based GSA workflow consisting of three phases – pre-analysis, analysis, and post-analysis – by integrating Monte Carlo and resampling methods with the repeated use of analysis of variance (ANOVA); we then exemplify this workflow using a two-dimensional multiscale lung cancer model. By accounting for all parameter rankings produced by multiple GSA methods, a summarized ranking is created at the end of the workflow based on the weighted mean of the rankings for each input parameter. For the cancer model investigated here, this analysis reveals that ERK, a downstream molecule of the EGFR signaling pathway, has the most important impact on regulating both the tumor volume and expansion rate in the algorithm used. PMID:25257020

Using the Community Multiscale Air Quality (CMAQ) model to estimate public health impacts of PM2.5 from individual power plants.

Science.gov (United States)

Buonocore, Jonathan J; Dong, Xinyi; Spengler, John D; Fu, Joshua S; Levy, Jonathan I

2014-07-01

We estimated PM2.5-related public health impacts/ton emitted of primary PM2.5, SO2, and NOx for a set of power plants in the Mid-Atlantic and Lower Great Lakes regions of the United States, selected to include varying emission profiles and broad geographic representation. We then developed a regression model explaining variability in impacts per ton emitted using the population distributions around each plant. We linked outputs from the Community Multiscale Air Quality (CMAQ) model v 4.7.1 with census data and concentration-response functions for PM2.5-related mortality, and monetized health estimates using the value-of-statistical-life. The median impacts for the final set of plants were $130,000/ton for primary PM2.5 (range: $22,000-230,000), $28,000/ton for SO2 (range: $19,000-33,000), and $16,000/ton for NOx (range: $7100-26,000). Impacts of NOx were a median of 34% (range: 20%-75%) from ammonium nitrate and 66% (range: 25%-79%) from ammonium sulfate. The latter pathway is likely from NOx enhancing atmospheric oxidative capacity and amplifying sulfate formation, and is often excluded. Our regression models explained most of the variation in impact/ton estimates using basic population covariates, and can aid in estimating impacts averted from interventions such as pollution controls, alternative energy installations, or demand-side management. Copyright © 2014 Elsevier Ltd. All rights reserved.

Multiscale Simulation Framework for Coupled Fluid Flow and Mechanical Deformation

Energy Technology Data Exchange (ETDEWEB)

Hou, Thomas [California Inst. of Technology (CalTech), Pasadena, CA (United States); Efendiev, Yalchin [Stanford Univ., CA (United States); Tchelepi, Hamdi [Texas A & M Univ., College Station, TX (United States); Durlofsky, Louis [Stanford Univ., CA (United States)

2016-05-24

Our work in this project is aimed at making fundamental advances in multiscale methods for flow and transport in highly heterogeneous porous media. The main thrust of this research is to develop a systematic multiscale analysis and efficient coarse-scale models that can capture global effects and extend existing multiscale approaches to problems with additional physics and uncertainties. A key emphasis is on problems without an apparent scale separation. Multiscale solution methods are currently under active investigation for the simulation of subsurface flow in heterogeneous formations. These procedures capture the effects of fine-scale permeability variations through the calculation of specialized coarse-scale basis functions. Most of the multiscale techniques presented to date employ localization approximations in the calculation of these basis functions. For some highly correlated (e.g., channelized) formations, however, global effects are important and these may need to be incorporated into the multiscale basis functions. Other challenging issues facing multiscale simulations are the extension of existing multiscale techniques to problems with additional physics, such as compressibility, capillary effects, etc. In our project, we explore the improvement of multiscale methods through the incorporation of additional (single-phase flow) information and the development of a general multiscale framework for flows in the presence of uncertainties, compressible flow and heterogeneous transport, and geomechanics. We have considered (1) adaptive local-global multiscale methods, (2) multiscale methods for the transport equation, (3) operator-based multiscale methods and solvers, (4) multiscale methods in the presence of uncertainties and applications, (5) multiscale finite element methods for high contrast porous media and their generalizations, and (6) multiscale methods for geomechanics.

Evaluating 20th Century precipitation characteristics between multi-scale atmospheric models with different land-atmosphere coupling

Science.gov (United States)

Phillips, M.; Denning, A. S.; Randall, D. A.; Branson, M.

2016-12-01

Multi-scale models of the atmosphere provide an opportunity to investigate processes that are unresolved by traditional Global Climate Models while at the same time remaining viable in terms of computational resources for climate-length time scales. The MMF represents a shift away from large horizontal grid spacing in traditional GCMs that leads to overabundant light precipitation and lack of heavy events, toward a model where precipitation intensity is allowed to vary over a much wider range of values. Resolving atmospheric motions on the scale of 4 km makes it possible to recover features of precipitation, such as intense downpours, that were previously only obtained by computationally expensive regional simulations. These heavy precipitation events may have little impact on large-scale moisture and energy budgets, but are outstanding in terms of interaction with the land surface and potential impact on human life. Three versions of the Community Earth System Model were used in this study; the standard CESM, the multi-scale `Super-Parameterized' CESM where large-scale parameterizations have been replaced with a 2D cloud-permitting model, and a multi-instance land version of the SP-CESM where each column of the 2D CRM is allowed to interact with an individual land unit. These simulations were carried out using prescribed Sea Surface Temperatures for the period from 1979-2006 with daily precipitation saved for all 28 years. Comparisons of the statistical properties of precipitation between model architectures and against observations from rain gauges were made, with specific focus on detection and evaluation of extreme precipitation events.

Multiscale sampling model for motion integration.

Science.gov (United States)

Sherbakov, Lena; Yazdanbakhsh, Arash

2013-09-30

Biologically plausible strategies for visual scene integration across spatial and temporal domains continues to be a challenging topic. The fundamental question we address is whether classical problems in motion integration, such as the aperture problem, can be solved in a model that samples the visual scene at multiple spatial and temporal scales in parallel. We hypothesize that fast interareal connections that allow feedback of information between cortical layers are the key processes that disambiguate motion direction. We developed a neural model showing how the aperture problem can be solved using different spatial sampling scales between LGN, V1 layer 4, V1 layer 6, and area MT. Our results suggest that multiscale sampling, rather than feedback explicitly, is the key process that gives rise to end-stopped cells in V1 and enables area MT to solve the aperture problem without the need for calculating intersecting constraints or crafting intricate patterns of spatiotemporal receptive fields. Furthermore, the model explains why end-stopped cells no longer emerge in the absence of V1 layer 6 activity (Bolz & Gilbert, 1986), why V1 layer 4 cells are significantly more end-stopped than V1 layer 6 cells (Pack, Livingstone, Duffy, & Born, 2003), and how it is possible to have a solution to the aperture problem in area MT with no solution in V1 in the presence of driving feedback. In summary, while much research in the field focuses on how a laminar architecture can give rise to complicated spatiotemporal receptive fields to solve problems in the motion domain, we show that one can reframe motion integration as an emergent property of multiscale sampling achieved concurrently within lamina and across multiple visual areas.

Computational design and multiscale modeling of a nanoactuator using DNA actuation

International Nuclear Information System (INIS)

Hamdi, Mustapha

2009-01-01

Developments in the field of nano-biodevices coupling nanostructures and biological components are of great interest in medical nanorobotics. As the fundamentals of bio/non-bio interaction processes are still poorly understood in the design of these devices, design tools and multiscale dynamics modeling approaches are necessary at the fabrication pre-project stage. This paper proposes a new concept of optimized carbon nanotube based servomotor design for drug delivery and biomolecular transport applications. The design of an encapsulated DNA-multi-walled carbon nanotube actuator is prototyped using multiscale modeling. The system is parametrized by using a quantum level approach and characterized by using a molecular dynamics simulation. Based on the analysis of the simulation results, a servo nanoactuator using ionic current feedback is simulated and analyzed for application as a drug delivery carrier.

Simulation of left atrial function using a multi-scale model of the cardiovascular system.

Directory of Open Access Journals (Sweden)

Antoine Pironet

Full Text Available During a full cardiac cycle, the left atrium successively behaves as a reservoir, a conduit and a pump. This complex behavior makes it unrealistic to apply the time-varying elastance theory to characterize the left atrium, first, because this theory has known limitations, and second, because it is still uncertain whether the load independence hypothesis holds. In this study, we aim to bypass this uncertainty by relying on another kind of mathematical model of the cardiac chambers. In the present work, we describe both the left atrium and the left ventricle with a multi-scale model. The multi-scale property of this model comes from the fact that pressure inside a cardiac chamber is derived from a model of the sarcomere behavior. Macroscopic model parameters are identified from reference dog hemodynamic data. The multi-scale model of the cardiovascular system including the left atrium is then simulated to show that the physiological roles of the left atrium are correctly reproduced. This include a biphasic pressure wave and an eight-shaped pressure-volume loop. We also test the validity of our model in non basal conditions by reproducing a preload reduction experiment by inferior vena cava occlusion with the model. We compute the variation of eight indices before and after this experiment and obtain the same variation as experimentally observed for seven out of the eight indices. In summary, the multi-scale mathematical model presented in this work is able to correctly account for the three roles of the left atrium and also exhibits a realistic left atrial pressure-volume loop. Furthermore, the model has been previously presented and validated for the left ventricle. This makes it a proper alternative to the time-varying elastance theory if the focus is set on precisely representing the left atrial and left ventricular behaviors.

Multiscale Simulation of Breaking Wave Impacts

DEFF Research Database (Denmark)

Lindberg, Ole

compare reasonably well. The incompressible and inviscid ALE-WLS model is coupled with the potential flow model of Engsig-Karup et al. [2009], to perform multiscale calculation of breaking wave impacts on a vertical breakwater. The potential flow model provides accurate calculation of the wave...... with a potential flow model to provide multiscale calculation of forces from breaking wave impacts on structures....

Multi-scale habitat selection modeling: A review and outlook

Science.gov (United States)

Kevin McGarigal; Ho Yi Wan; Kathy A. Zeller; Brad C. Timm; Samuel A. Cushman

2016-01-01

Scale is the lens that focuses ecological relationships. Organisms select habitat at multiple hierarchical levels and at different spatial and/or temporal scales within each level. Failure to properly address scale dependence can result in incorrect inferences in multi-scale habitat selection modeling studies.

Conformal-Based Surface Morphing and Multi-Scale Representation

Directory of Open Access Journals (Sweden)

Ka Chun Lam

2014-05-01

Full Text Available This paper presents two algorithms, based on conformal geometry, for the multi-scale representations of geometric shapes and surface morphing. A multi-scale surface representation aims to describe a 3D shape at different levels of geometric detail, which allows analyzing or editing surfaces at the global or local scales effectively. Surface morphing refers to the process of interpolating between two geometric shapes, which has been widely applied to estimate or analyze deformations in computer graphics, computer vision and medical imaging. In this work, we propose two geometric models for surface morphing and multi-scale representation for 3D surfaces. The basic idea is to represent a 3D surface by its mean curvature function, H, and conformal factor function λ, which uniquely determine the geometry of the surface according to Riemann surface theory. Once we have the (λ, H parameterization of the surface, post-processing of the surface can be done directly on the conformal parameter domain. In particular, the problem of multi-scale representations of shapes can be reduced to the signal filtering on the λ and H parameters. On the other hand, the surface morphing problem can be transformed to an interpolation process of two sets of (λ, H parameters. We test the proposed algorithms on 3D human face data and MRI-derived brain surfaces. Experimental results show that our proposed methods can effectively obtain multi-scale surface representations and give natural surface morphing results.

Multiscale Modeling of Carbon/Phenolic Composite Thermal Protection Materials: Atomistic to Effective Properties

Science.gov (United States)

Arnold, Steven M.; Murthy, Pappu L.; Bednarcyk, Brett A.; Lawson, John W.; Monk, Joshua D.; Bauschlicher, Charles W., Jr.

2016-01-01

Next generation ablative thermal protection systems are expected to consist of 3D woven composite architectures. It is well known that composites can be tailored to achieve desired mechanical and thermal properties in various directions and thus can be made fit-for-purpose if the proper combination of constituent materials and microstructures can be realized. In the present work, the first, multiscale, atomistically-informed, computational analysis of mechanical and thermal properties of a present day - Carbon/Phenolic composite Thermal Protection System (TPS) material is conducted. Model results are compared to measured in-plane and out-of-plane mechanical and thermal properties to validate the computational approach. Results indicate that given sufficient microstructural fidelity, along with lowerscale, constituent properties derived from molecular dynamics simulations, accurate composite level (effective) thermo-elastic properties can be obtained. This suggests that next generation TPS properties can be accurately estimated via atomistically informed multiscale analysis.

Multi-scale modeling of the thermo-mechanical behavior of particle-based composites

International Nuclear Information System (INIS)

Di Paola, F.

2010-01-01

The aim of this work was to perform numerical simulations of the thermal and mechanical behavior of a particle-based nuclear fuel. This is a refractory composite material made of UO 2 spherical particles which are coated with two layers of pyrocarbon and embedded in a graphite matrix at a high volume fraction (45%). The objective was to develop a multi-scale modeling of this composite material which can estimate its mean behavior as well as the heterogeneity of the local mechanical variables. The first part of this work was dedicated to the modeling of the microstructure in 3D. To do this, we developed tools to generate random distributions of spheres, meshes and to characterize the morphology of the microstructure towards the finite element code Cast3M. A hundred of numerical samples of the composite were created. The second part was devoted to the characterization of the thermo-elastic behavior by the finite element modeling of the samples. We studied the influence of different modeling parameters, one of them is the boundary conditions. We proposed a method to vanish the boundary conditions effects from the computed solution by analyzing it on an internal sub-volume of the sample obtained by erosion. Then, we determined the effective properties (elastic moduli, thermal conductivity and thermal expansion) and the stress distribution within the matrix. Finally, in the third part we proposed a multi-scale modeling to determine the mean values and the variance and covariance of the local mechanical variables for any macroscopic load. This statistical approach have been used to estimate the intra-phase distribution of these variables in the composite material. (author) [fr

Multi-scale modeling of the thermo-mechanical behavior of particle-based composites

International Nuclear Information System (INIS)

Di Paola, F.

2010-11-01

The aim of this work was to perform numerical simulations of the thermal and mechanical behavior of a particle-based nuclear fuel. This is a refractory composite material made of UO 2 spherical particles which are coated with two layers of pyrocarbon and embedded in a graphite matrix at a high volume fraction (45 %). The objective was to develop a multi-scale modeling of this composite material which can estimate its mean behavior as well as the heterogeneity of the local mechanical variables. The first part of this work was dedicated to the modeling of the microstructure in 3D. To do this, we developed tools to generate random distributions of spheres, meshes and to characterize the morphology of the microstructure towards the finite element code Cast3M. A hundred of numerical samples of the composite were created. The second part was devoted to the characterization of the thermo-elastic behavior by the finite element modeling of the samples. We studied the influence of different modeling parameters, one of them is the boundary conditions. We proposed a method to vanish the boundary conditions effects from the computed solution by analyzing it on an internal sub-volume of the sample obtained by erosion. Then, we determined the effective properties (elastic moduli, thermal conductivity and thermal expansion) and the stress distribution within the matrix. Finally, in the third part we proposed a multi-scale modeling to determine the mean values and the variance and covariance of the local mechanical variables for any macroscopic load. This statistical approach have been used to estimate the intra-phase distribution of these variables in the composite material. (author)

The Multi-Scale Model Approach to Thermohydrology at Yucca Mountain

International Nuclear Information System (INIS)

Glascoe, L; Buscheck, T A; Gansemer, J; Sun, Y

2002-01-01

The Multi-Scale Thermo-Hydrologic (MSTH) process model is a modeling abstraction of them1 hydrology (TH) of the potential Yucca Mountain repository at multiple spatial scales. The MSTH model as described herein was used for the Supplemental Science and Performance Analyses (BSC, 2001) and is documented in detail in CRWMS M and O (2000) and Glascoe et al. (2002). The model has been validated to a nested grid model in Buscheck et al. (In Review). The MSTH approach is necessary for modeling thermal hydrology at Yucca Mountain for two reasons: (1) varying levels of detail are necessary at different spatial scales to capture important TH processes and (2) a fully-coupled TH model of the repository which includes the necessary spatial detail is computationally prohibitive. The MSTH model consists of six ''submodels'' which are combined in a manner to reduce the complexity of modeling where appropriate. The coupling of these models allows for appropriate consideration of mountain-scale thermal hydrology along with the thermal hydrology of drift-scale discrete waste packages of varying heat load. Two stages are involved in the MSTH approach, first, the execution of submodels, and second, the assembly of submodels using the Multi-scale Thermohydrology Abstraction Code (MSTHAC). MSTHAC assembles the submodels in a five-step process culminating in the TH model output of discrete waste packages including a mountain-scale influence

Polynomial Chaos Characterization of Uncertainty in Multiscale Models and Behavior of Carbon Reinforced Composites

Energy Technology Data Exchange (ETDEWEB)

Mehrez, Loujaine [University of Southern California; Ghanem, Roger [University of Southern California; Aitharaju, Venkat [General Motors; Rodgers, William [General Motors

2017-10-23

Design of non-crimp fabric (NCF) composites entails major challenges pertaining to (1) the complex fine-scale morphology of the constituents, (2) the manufacturing-produced inconsistency of this morphology spatially, and thus (3) the ability to build reliable, robust, and efficient computational surrogate models to account for this complex nature. Traditional approaches to construct computational surrogate models have been to average over the fluctuations of the material properties at different scale lengths. This fails to account for the fine-scale features and fluctuations in morphology, material properties of the constituents, as well as fine-scale phenomena such as damage and cracks. In addition, it fails to accurately predict the scatter in macroscopic properties, which is vital to the design process and behavior prediction. In this work, funded in part by the Department of Energy, we present an approach for addressing these challenges by relying on polynomial chaos representations of both input parameters and material properties at different scales. Moreover, we emphasize the efficiency and robustness of integrating the polynomial chaos expansion with multiscale tools to perform multiscale assimilation, characterization, propagation, and prediction, all of which are necessary to construct the data-driven surrogate models required to design under the uncertainty of composites. These data-driven constructions provide an accurate map from parameters (and their uncertainties) at all scales and the system-level behavior relevant for design. While this perspective is quite general and applicable to all multiscale systems, NCF composites present a particular hierarchy of scales that permits the efficient implementation of these concepts.

Multiscale computer modeling in biomechanics and biomedical engineering

CERN Document Server

2013-01-01

This book reviews the state-of-the-art in multiscale computer modeling, in terms of both accomplishments and challenges. The information in the book is particularly useful for biomedical engineers, medical physicists and researchers in systems biology, mathematical biology, micro-biomechanics and biomaterials who are interested in how to bridge between traditional biomedical engineering work at the organ and tissue scales, and the newer arenas of cellular and molecular bioengineering.

Comparison of single- and multi-scale models for the prediction of the Culicoides biting midge distribution in Germany

Directory of Open Access Journals (Sweden)

Renke Lühken

2016-05-01

Full Text Available This study analysed Culicoides presence-absence data from 46 sampling sites in Germany, where monitoring was carried out from April 2007 until May 2008. Culicoides presence-absence data were analysed in relation to land cover data, in order to study whether the prevalence of biting midges is correlated to land cover data with respect to the trapping sites. We differentiated eight scales, i.e. buffer zones with radii of 0.5, 1, 2, 3, 4, 5, 7.5 and 10 km, around each site, and chose several land cover variables. For each species, we built eight single-scale models (i.e. predictor variables from one of the eight scales for each model based on averaged, generalised linear models and two multiscale models (i.e. predictor variables from all of the eight scales based on averaged, generalised linear models and generalised linear models with random forest variable selection. There were no significant differences between performance indicators of models built with land cover data from different buffer zones around the trapping sites. However, the overall performance of multi-scale models was higher than the alternatives. Furthermore, these models mostly achieved the best performance for the different species using the index area under the receiver operating characteristic curve. However, as also presented in this study, the relevance of the different variables could significantly differ between various scales, including the number of species affected and the positive or negative direction. This is an even more severe problem if multi-scale models are concerned, in which one model can have the same variable at different scales but with different directions, i.e. negative and positive direction of the same variable at different scales. However, multi-scale modelling is a promising approach to model the distribution of Culicoides species, accounting much more for the ecology of biting midges, which uses different resources (breeding sites, hosts, etc. at

Multiscale model reduction for shale gas transport in fractured media

KAUST Repository

Akkutlu, I. Y.; Efendiev, Yalchin R.; Vasilyeva, Maria

2016-01-01

fracture distributions on an unstructured grid; (2) develop GMsFEM for nonlinear flows; and (3) develop online basis function strategies to adaptively improve the convergence. The number of multiscale basis functions in each coarse region represents

GEMMER: GEnome-wide tool for Multi-scale Modeling data Extraction and Representation for Saccharomyces cerevisiae.

Science.gov (United States)

Mondeel, Thierry D G A; Crémazy, Frédéric; Barberis, Matteo

2018-02-01

Multi-scale modeling of biological systems requires integration of various information about genes and proteins that are connected together in networks. Spatial, temporal and functional information is available; however, it is still a challenge to retrieve and explore this knowledge in an integrated, quick and user-friendly manner. We present GEMMER (GEnome-wide tool for Multi-scale Modelling data Extraction and Representation), a web-based data-integration tool that facilitates high quality visualization of physical, regulatory and genetic interactions between proteins/genes in Saccharomyces cerevisiae. GEMMER creates network visualizations that integrate information on function, temporal expression, localization and abundance from various existing databases. GEMMER supports modeling efforts by effortlessly gathering this information and providing convenient export options for images and their underlying data. GEMMER is freely available at http://gemmer.barberislab.com. Source code, written in Python, JavaScript library D3js, PHP and JSON, is freely available at https://github.com/barberislab/GEMMER. M.Barberis@uva.nl. Supplementary data are available at Bioinformatics online. © The Author(s) 2018. Published by Oxford University Press.

Multiscale modeling of materials. Materials Research Society symposium proceedings: Volume 538

International Nuclear Information System (INIS)

Bulatov, V.V.; Diaz de la Rubia, T.; Phillips, R.; Kaxiras, E.; Ghoniem, N.

1999-01-01

The symposium, Multiscale Modeling of Materials, was held at the 1998 MRS Fall Meeting in Boston, Massachusetts, November 30 to December 3. Though multiple scale models are not new the topic has recently taken on a new sense of urgency. This is in large part due to the recognition that brute force computational approaches often fall short of allowing for direct simulation of both the characteristic structures and temporal processes found in real materials. As a result, a number of hybrid approaches are now finding favor in which ideas borrowed from distinct disciplines or modeling paradigms are unified to produce more powerful techniques. Topics included are modeling dislocation properties and behavior, defect dynamics and microstructural evolution, crystal defects and interfaces, novel methods for materials modeling, and non-crystalline and nanocrystalline materials. Eighty papers have been processed separately for inclusion on the data base

Bridging scales through multiscale modeling: A case study on Protein Kinase A

Directory of Open Access Journals (Sweden)

Sophia P Hirakis

2015-09-01

Full Text Available The goal of multiscale modeling in biology is to use structurally based physico-chemical models to integrate across temporal and spatial scales of biology and thereby improve mechanistic understanding of, for example, how a single mutation can alter organism-scale phenotypes. This approach may also inform therapeutic strategies or identify candidate drug targets that might otherwise have been overlooked. However, in many cases, it remains unclear how best to synthesize information obtained from various scales and analysis approaches, such as atomistic molecular models, Markov state models (MSM, subcellular network models, and whole cell models. In this paper, we use protein kinase A (PKA activation as a case study to explore how computational methods that model different physical scales can complement each other and integrate into an improved multiscale representation of the biological mechanisms. Using measured crystal structures, we show how molecular dynamics (MD simulations coupled with atomic-scale MSMs can provide conformations for Brownian dynamics (BD simulations to feed transitional states and kinetic parameters into protein-scale MSMs. We discuss how milestoning can give reaction probabilities and forward-rate constants of cAMP association events by seamlessly integrating MD and BD simulation scales. These rate constants coupled with MSMs provide a robust representation of the free energy landscape, enabling access to kinetic and thermodynamic parameters unavailable from current experimental data. These approaches have helped to illuminate the cooperative nature of PKA activation in response to distinct cAMP binding events. Collectively, this approach exemplifies a general strategy for multiscale model development that is applicable to a wide range of biological problems.

Community Land Model Version 3.0 (CLM3.0) Developer's Guide

Energy Technology Data Exchange (ETDEWEB)

Hoffman, FM

2004-12-21

This document describes the guidelines adopted for software development of the Community Land Model (CLM) and serves as a reference to the entire code base of the released version of the model. The version of the code described here is Version 3.0 which was released in the summer of 2004. This document, the Community Land Model Version 3.0 (CLM3.0) User's Guide (Vertenstein et al., 2004), the Technical Description of the Community Land Model (CLM) (Oleson et al., 2004), and the Community Land Model's Dynamic Global Vegetation Model (CLM-DGVM): Technical Description and User's Guide (Levis et al., 2004) provide the developer, user, or researcher with details of implementation, instructions for using the model, a scientific description of the model, and a scientific description of the Dynamic Global Vegetation Model integrated with CLM respectively. The CLM is a single column (snow-soil-vegetation) biogeophysical model of the land surface which can be run serially (on a laptop or personal computer) or in parallel (using distributed or shared memory processors or both) on both vector and scalar computer architectures. Written in Fortran 90, CLM can be run offline (i.e., run in isolation using stored atmospheric forcing data), coupled to an atmospheric model (e.g., the Community Atmosphere Model (CAM)), or coupled to a climate system model (e.g., the Community Climate System Model Version 3 (CCSM3)) through a flux coupler (e.g., Coupler 6 (CPL6)). When coupled, CLM exchanges fluxes of energy, water, and momentum with the atmosphere. The horizontal land surface heterogeneity is represented by a nested subgrid hierarchy composed of gridcells, landunits, columns, and plant functional types (PFTs). This hierarchical representation is reflected in the data structures used by the model code. Biophysical processes are simulated for each subgrid unit (landunit, column, and PFT) independently, and prognostic variables are maintained for each subgrid unit

Foundations for a multiscale collaborative Earth model

KAUST Repository

Afanasiev, M.

2015-11-11

of the CSEM development, the broad global updates mostly act to remove artefacts from the assembly of the initial CSEM. During the future evolution of the CSEM, the reference data set will be used to account for the influence of small-scale refinements on large-scale global structure. The CSEM as a computational framework is intended to help bridging the gap between local, regional and global tomography, and to contribute to the development of a global multiscale Earth model. While the current construction serves as a first proof of concept, future refinements and additions will require community involvement, which is welcome at this stage already.

Multi-scale modeling with cellular automata: The complex automata approach

NARCIS (Netherlands)

Hoekstra, A.G.; Falcone, J.-L.; Caiazzo, A.; Chopard, B.

2008-01-01

Cellular Automata are commonly used to describe complex natural phenomena. In many cases it is required to capture the multi-scale nature of these phenomena. A single Cellular Automata model may not be able to efficiently simulate a wide range of spatial and temporal scales. It is our goal to

Multi-scale modelling of non-uniform consolidation of uncured toughened unidirectional prepregs

Science.gov (United States)

Sorba, G.; Binetruy, C.; Syerko, E.; Leygue, A.; Comas-Cardona, S.; Belnoue, J. P.-H.; Nixon-Pearson, O. J.; Ivanov, D. S.; Hallett, S. R.; Advani, S. G.

2018-05-01

Consolidation is a crucial step in manufacturing of composite parts with prepregs because its role is to eliminate inter- and intra-ply gaps and porosity. Some thermoset prepreg systems are toughened with thermoplastic particles. Depending on their size, thermoplastic particles can be either located in between plies or distributed within the inter-fibre regions. When subjected to transverse compaction, resin will bleed out of low-viscosity unidirectional prepregs along the fibre direction, whereas one would expect transverse squeeze flow to dominate for higher viscosity prepregs. Recent experimental work showed that the consolidation of uncured toughened prepregs involves complex flow and deformation mechanisms where both bleeding and squeeze flow patterns are observed [1]. Micrographs of compacted and cured samples confirm these features as shown in Fig.1. A phenomenological model was proposed [2] where bleeding flow and squeeze flow are combined. A criterion for the transition from shear flow to resin bleeding was also proposed. However, the micrographs also reveal a resin rich layer between plies which may be contributing to the complex flow mechanisms during the consolidation process. In an effort to provide additional insight into these complex mechanisms, this work focuses on the 3D numerical modelling of the compaction of uncured toughened prepregs in the cross-ply configuration described in [1]. A transversely isotropic fluid model is used to describe the flow behaviour of the plies coupled with interplay resin flow of an isotropic fluid. The multi-scale flow model used is based on [3, 4]. A numerical parametric study is carried out where the resin viscosity, permeability and inter-ply thickness are varied to identify the role of important variables. The squeezing flow and the bleeding flow are compared for a range of process parameters to investigate the coupling and competition between the two flow mechanisms. Figure 4 shows the predicted displacement of

Multiscale Modeling of Grain Boundary Segregation and Embrittlement in Tungsten for Mechanistic Design of Alloys for Coal Fired Plants

Energy Technology Data Exchange (ETDEWEB)

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.

Multiscale Mathematics for Biomass Conversion to Renewable Hydrogen

Energy Technology Data Exchange (ETDEWEB)

Plechac, Petr [Univ. of Tennessee, Knoxville, TN (United States). Dept. of Mathematics; Univ. of Delaware, Newark, DE (United States). Dept. of Mathematics; Vlachos, Dionisios [Univ. of Delaware, Newark, DE (United States). Dept. of Chemical and Biomolecular Engineering; Katsoulakis, Markos [Univ. of Massachusetts, Amherst, MA (United States). Dept. of Mathematics

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.

Multi-Scale Modeling of the Gamma Radiolysis of Nitrate Solutions.

Science.gov (United States)

Horne, Gregory P; Donoclift, Thomas A; Sims, Howard E; Orr, Robin M; Pimblott, Simon M

2016-11-17

A multiscale modeling approach has been developed for the extended time scale long-term radiolysis of aqueous systems. The approach uses a combination of stochastic track structure and track chemistry as well as deterministic homogeneous chemistry techniques and involves four key stages: radiation track structure simulation, the subsequent physicochemical processes, nonhomogeneous diffusion-reaction kinetic evolution, and homogeneous bulk chemistry modeling. The first three components model the physical and chemical evolution of an isolated radiation chemical track and provide radiolysis yields, within the extremely low dose isolated track paradigm, as the input parameters for a bulk deterministic chemistry model. This approach to radiation chemical modeling has been tested by comparison with the experimentally observed yield of nitrite from the gamma radiolysis of sodium nitrate solutions. This is a complex radiation chemical system which is strongly dependent on secondary reaction processes. The concentration of nitrite is not just dependent upon the evolution of radiation track chemistry and the scavenging of the hydrated electron and its precursors but also on the subsequent reactions of the products of these scavenging reactions with other water radiolysis products. Without the inclusion of intratrack chemistry, the deterministic component of the multiscale model is unable to correctly predict experimental data, highlighting the importance of intratrack radiation chemistry in the chemical evolution of the irradiated system.

A phenomenological variational multiscale constitutive model for intergranular failure in nanocrystalline materials

KAUST Repository

Siddiq, A.; El Sayed, Tamer S.

2013-01-01

We present a variational multiscale constitutive model that accounts for intergranular failure in nanocrystalline fcc metals due to void growth and coalescence in the grain boundary region. Following previous work by the authors, a nanocrystalline

Evaluation of a seven-year air quality simulation using the Weather Research and Forecasting (WRF)/Community Multiscale Air Quality (CMAQ) models in the eastern United States.

Science.gov (United States)

Zhang, Hongliang; Chen, Gang; Hu, Jianlin; Chen, Shu-Hua; Wiedinmyer, Christine; Kleeman, Michael; Ying, Qi

2014-03-01

The performance of the Weather Research and Forecasting (WRF)/Community Multi-scale Air Quality (CMAQ) system in the eastern United States is analyzed based on results from a seven-year modeling study with a 4-km spatial resolution. For 2-m temperature, the monthly averaged mean bias (MB) and gross error (GE) values are generally within the recommended performance criteria, although temperature is over-predicted with MB values up to 2K. Water vapor at 2-m is well-predicted but significant biases (>2 g kg(-1)) were observed in wintertime. Predictions for wind speed are satisfactory but biased towards over-prediction with 0nitrate and sulfate concentrations are also well reproduced. The other unresolved PM2.5 components (OTHER) are significantly overestimated by more than a factor of two. No conclusive explanations can be made regarding the possible cause of this universal overestimation, which warrants a follow-up study to better understand this problem. Copyright © 2013 Elsevier B.V. All rights reserved.

A MULTISCALE FRAMEWORK FOR THE STOCHASTIC ASSIMILATION AND MODELING OF UNCERTAINTY ASSOCIATED NCF COMPOSITE MATERIALS

Energy Technology Data Exchange (ETDEWEB)

Mehrez, Loujaine [University of Southern California; Ghanem, Roger [University of Southern California; McAuliffe, Colin [Altair Engineering, Inc.; Aitharaju, Venkat [General Motors; Rodgers, William [General Motors

2016-06-06

multiscale framework to construct stochastic macroscopic constitutive material models is proposed. A spectral projection approach, specifically polynomial chaos expansion, has been used to construct explicit functional relationships between the homogenized properties and input parameters from finer scales. A homogenization engine embedded in Multiscale Designer, software for composite materials, has been used for the upscaling process. The framework is demonstrated using non-crimp fabric composite materials by constructing probabilistic models of the homogenized properties of a non-crimp fabric laminate in terms of the input parameters together with the homogenized properties from finer scales.

Multiscale modeling for the prediction of casting defects in investment cast aluminum alloys

International Nuclear Information System (INIS)

Hamilton, R.W.; See, D.; Butler, S.; Lee, P.D.

2003-01-01

Macroscopic modeling of heat transfer and fluid flow is now routinely used for the prediction of macroscopic defects in castings, while microscopic models are used to investigate the effects of alloy changes on typical microstructures. By combining these two levels of modeling it is possible to simulate the casting process over a wider range of spatial and temporal scales. This paper presents a multiscale model where micromodels for dendrite arm spacing and microporosity are incorporated into a macromodel of heat transfer and in order to predict the as cast microstructure and prevalence of microscopic defects, specifically porosity. The approach is applied to aluminum alloy (L169) investment castings. The models are compared with results obtained by optical image analysis of prepared slices, and X-ray tomography of volume samples from the experiments. Multiscale modeling is shown to provide the designer with a useful tool to improve the properties of the final casting by testing how altering the casting process affects the final microstructure including porosity

A multi-scale modeling of surface effect via the modified boundary Cauchy-Born model

Energy Technology Data Exchange (ETDEWEB)

Khoei, A.R., E-mail: arkhoei@sharif.edu; Aramoon, A.

2012-10-01

In this paper, a new multi-scale approach is presented based on the modified boundary Cauchy-Born (MBCB) technique to model the surface effects of nano-structures. The salient point of the MBCB model is the definition of radial quadrature used in the surface elements which is an indicator of material behavior. The characteristics of quadrature are derived by interpolating data from atoms laid in a circular support around the quadrature, in a least-square scene. The total-Lagrangian formulation is derived for the equivalent continua by employing the Cauchy-Born hypothesis for calculating the strain energy density function of the continua. The numerical results of the proposed method are compared with direct atomistic and finite element simulation results to indicate that the proposed technique provides promising results for modeling surface effects of nano-structures. - Highlights: Black-Right-Pointing-Pointer A multi-scale approach is presented to model the surface effects in nano-structures. Black-Right-Pointing-Pointer The total-Lagrangian formulation is derived by employing the Cauchy-Born hypothesis. Black-Right-Pointing-Pointer The radial quadrature is used to model the material behavior in surface elements. Black-Right-Pointing-Pointer The quadrature characteristics are derived using the data at the atomistic level.

A multi-scale computational model of the effects of TMS on motor cortex [version 3; referees: 2 approved

Directory of Open Access Journals (Sweden)

Hyeon Seo

2017-05-01

Full Text Available The detailed biophysical mechanisms through which transcranial magnetic stimulation (TMS activates cortical circuits are still not fully understood. Here we present a multi-scale computational model to describe and explain the activation of different pyramidal cell types in motor cortex due to TMS. Our model determines precise electric fields based on an individual head model derived from magnetic resonance imaging and calculates how these electric fields activate morphologically detailed models of different neuron types. We predict neural activation patterns for different coil orientations consistent with experimental findings. Beyond this, our model allows us to calculate activation thresholds for individual neurons and precise initiation sites of individual action potentials on the neurons’ complex morphologies. Specifically, our model predicts that cortical layer 3 pyramidal neurons are generally easier to stimulate than layer 5 pyramidal neurons, thereby explaining the lower stimulation thresholds observed for I-waves compared to D-waves. It also shows differences in the regions of activated cortical layer 5 and layer 3 pyramidal cells depending on coil orientation. Finally, it predicts that under standard stimulation conditions, action potentials are mostly generated at the axon initial segment of cortical pyramidal cells, with a much less important activation site being the part of a layer 5 pyramidal cell axon where it crosses the boundary between grey matter and white matter. In conclusion, our computational model offers a detailed account of the mechanisms through which TMS activates different cortical pyramidal cell types, paving the way for more targeted application of TMS based on individual brain morphology in clinical and basic research settings.

Multiscale Molecular Dynamics Model for Heterogeneous Charged Systems

Science.gov (United States)

Stanton, L. G.; Glosli, J. N.; Murillo, M. S.

2018-04-01

Modeling matter across large length scales and timescales using molecular dynamics simulations poses significant challenges. These challenges are typically addressed through the use of precomputed pair potentials that depend on thermodynamic properties like temperature and density; however, many scenarios of interest involve spatiotemporal variations in these properties, and such variations can violate assumptions made in constructing these potentials, thus precluding their use. In particular, when a system is strongly heterogeneous, most of the usual simplifying assumptions (e.g., spherical potentials) do not apply. Here, we present a multiscale approach to orbital-free density functional theory molecular dynamics (OFDFT-MD) simulations that bridges atomic, interionic, and continuum length scales to allow for variations in hydrodynamic quantities in a consistent way. Our multiscale approach enables simulations on the order of micron length scales and 10's of picosecond timescales, which exceeds current OFDFT-MD simulations by many orders of magnitude. This new capability is then used to study the heterogeneous, nonequilibrium dynamics of a heated interface characteristic of an inertial-confinement-fusion capsule containing a plastic ablator near a fuel layer composed of deuterium-tritium ice. At these scales, fundamental assumptions of continuum models are explored; features such as the separation of the momentum fields among the species and strong hydrogen jetting from the plastic into the fuel region are observed, which had previously not been seen in hydrodynamic simulations.

Multi-scale Material Parameter Identification Using LS-DYNA® and LS-OPT®

Energy Technology Data Exchange (ETDEWEB)

Stander, Nielen; Basudhar, Anirban; Basu, Ushnish; Gandikota, Imtiaz; Savic, Vesna; Sun, Xin; Choi, Kyoo Sil; Hu, Xiaohua; Pourboghrat, F.; Park, Taejoon; Mapar, Aboozar; Kumar, Shavan; Ghassemi-Armaki, Hassan; Abu-Farha, Fadi

2015-09-14

Ever-tightening regulations on fuel economy, and the likely future regulation of carbon emissions, demand persistent innovation in vehicle design to reduce vehicle mass. Classical methods for computational mass reduction include sizing, shape and topology optimization. One of the few remaining options for weight reduction can be found in materials engineering and material design optimization. Apart from considering different types of materials, by adding material diversity and composite materials, an appealing option in automotive design is to engineer steel alloys for the purpose of reducing plate thickness while retaining sufficient strength and ductility required for durability and safety. A project to develop computational material models for advanced high strength steel is currently being executed under the auspices of the United States Automotive Materials Partnership (USAMP) funded by the US Department of Energy. Under this program, new Third Generation Advanced High Strength Steel (i.e., 3GAHSS) are being designed, tested and integrated with the remaining design variables of a benchmark vehicle Finite Element model. The objectives of the project are to integrate atomistic, microstructural, forming and performance models to create an integrated computational materials engineering (ICME) toolkit for 3GAHSS. The mechanical properties of Advanced High Strength Steels (AHSS) are controlled by many factors, including phase composition and distribution in the overall microstructure, volume fraction, size and morphology of phase constituents as well as stability of the metastable retained austenite phase. The complex phase transformation and deformation mechanisms in these steels make the well-established traditional techniques obsolete, and a multi-scale microstructure-based modeling approach following the ICME [0]strategy was therefore chosen in this project. Multi-scale modeling as a major area of research and development is an outgrowth of the Comprehensive

Coarse-graining and hybrid methods for efficient simulation of stochastic multi-scale models of tumour growth

International Nuclear Information System (INIS)

Cruz, Roberto de la; Guerrero, Pilar; Calvo, Juan; Alarcón, Tomás

2017-01-01

The development of hybrid methodologies is of current interest in both multi-scale modelling and stochastic reaction–diffusion systems regarding their applications to biology. We formulate a hybrid method for stochastic multi-scale models of cells populations that extends the remit of existing hybrid methods for reaction–diffusion systems. Such method is developed for a stochastic multi-scale model of tumour growth, i.e. population-dynamical models which account for the effects of intrinsic noise affecting both the number of cells and the intracellular dynamics. In order to formulate this method, we develop a coarse-grained approximation for both the full stochastic model and its mean-field limit. Such approximation involves averaging out the age-structure (which accounts for the multi-scale nature of the model) by assuming that the age distribution of the population settles onto equilibrium very fast. We then couple the coarse-grained mean-field model to the full stochastic multi-scale model. By doing so, within the mean-field region, we are neglecting noise in both cell numbers (population) and their birth rates (structure). This implies that, in addition to the issues that arise in stochastic-reaction diffusion systems, we need to account for the age-structure of the population when attempting to couple both descriptions. We exploit our coarse-graining model so that, within the mean-field region, the age-distribution is in equilibrium and we know its explicit form. This allows us to couple both domains consistently, as upon transference of cells from the mean-field to the stochastic region, we sample the equilibrium age distribution. Furthermore, our method allows us to investigate the effects of intracellular noise, i.e. fluctuations of the birth rate, on collective properties such as travelling wave velocity. We show that the combination of population and birth-rate noise gives rise to large fluctuations of the birth rate in the region at the leading edge

A multiscale modeling approach for biomolecular systems

Energy Technology Data Exchange (ETDEWEB)

Bowling, Alan, E-mail: bowling@uta.edu; Haghshenas-Jaryani, Mahdi, E-mail: mahdi.haghshenasjaryani@mavs.uta.edu [The University of Texas at Arlington, Department of Mechanical and Aerospace Engineering (United States)

2015-04-15

This paper presents a new multiscale molecular dynamic model for investigating the effects of external interactions, such as contact and impact, during stepping and docking of motor proteins and other biomolecular systems. The model retains the mass properties ensuring that the result satisfies Newton’s second law. This idea is presented using a simple particle model to facilitate discussion of the rigid body model; however, the particle model does provide insights into particle dynamics at the nanoscale. The resulting three-dimensional model predicts a significant decrease in the effect of the random forces associated with Brownian motion. This conclusion runs contrary to the widely accepted notion that the motor protein’s movements are primarily the result of thermal effects. This work focuses on the mechanical aspects of protein locomotion; the effect ATP hydrolysis is estimated as internal forces acting on the mechanical model. In addition, the proposed model can be numerically integrated in a reasonable amount of time. Herein, the differences between the motion predicted by the old and new modeling approaches are compared using a simplified model of myosin V.

Multiscale Modeling of Wear Degradation

KAUST Repository

Moraes, Alvaro; Ruggeri, Fabrizio; Tempone, Raul; Vilanova, Pedro

2016-01-01

Cylinder liners of diesel engines used for marine propulsion are naturally subjected to a wear process, and may fail when their wear exceeds a specified limit. Since failures often represent high economical costs, it is utterly important to predict and avoid them. In this work [4], we model the wear process using a pure jump process. Therefore, the inference goal here is to estimate: the number of possible jumps, its sizes, the coefficients and the shapes of the jump intensities. We propose a multiscale approach for the inference problem that can be seen as an indirect inference scheme. We found that using a Gaussian approximation based on moment expansions, it is possible to accurately estimate the jump intensities and the jump amplitudes. We obtained results equivalent to the state of the art but using a simpler and less expensive approach.

Multiscale Modeling of Wear Degradation

KAUST Repository

Moraes, Alvaro; Ruggeri, Fabrizio; Tempone, Raul; Vilanova, Pedro

2015-01-01

Cylinder liners of diesel engines used for marine propulsion are naturally subjected to a wear process, and may fail when their wear exceeds a specified limit. Since failures often represent high economical costs, it is utterly important to predict and avoid them. In this work [4], we model the wear process using a pure jump process. Therefore, the inference goal here is to estimate: the number of possible jumps, its sizes, the coefficients and the shapes of the jump intensities. We propose a multiscale approach for the inference problem that can be seen as an indirect inference scheme. We found that using a Gaussian approximation based on moment expansions, it is possible to accurately estimate the jump intensities and the jump amplitudes. We obtained results equivalent to the state of the art but using a simpler and less expensive approach.

Multiscale Modeling of Wear Degradation

KAUST Repository

Moraes, Alvaro

2015-01-07

Cylinder liners of diesel engines used for marine propulsion are naturally subjected to a wear process, and may fail when their wear exceeds a specified limit. Since failures often represent high economical costs, it is utterly important to predict and avoid them. In this work [4], we model the wear process using a pure jump process. Therefore, the inference goal here is to estimate: the number of possible jumps, its sizes, the coefficients and the shapes of the jump intensities. We propose a multiscale approach for the inference problem that can be seen as an indirect inference scheme. We found that using a Gaussian approximation based on moment expansions, it is possible to accurately estimate the jump intensities and the jump amplitudes. We obtained results equivalent to the state of the art but using a simpler and less expensive approach.

Multiscale Modeling of Wear Degradation

KAUST Repository

Moraes, Alvaro

2014-01-06

Cylinder liners of diesel engines used for marine propulsion are naturally subjected to a wear process, and may fail when their wear exceeds a specified limit. Since failures often represent high economical costs, it is utterly important to predict and avoid them. In this work [4], we model the wear process using a pure jump process. Therefore, the inference goal here is to estimate: the number of possible jumps, its sizes, the coefficients and the shapes of the jump intensities. We propose a multiscale approach for the inference problem that can be seen as an indirect inference scheme. We found that using a Gaussian approximation based on moment expansions, it is possible to accurately estimate the jump intensities and the jump amplitudes. We obtained results equivalent to the state of the art but using a simpler and less expensive approach.

Multiscale Modeling of Wear Degradation

KAUST Repository

Moraes, Alvaro

2016-01-06

Cylinder liners of diesel engines used for marine propulsion are naturally subjected to a wear process, and may fail when their wear exceeds a specified limit. Since failures often represent high economical costs, it is utterly important to predict and avoid them. In this work [4], we model the wear process using a pure jump process. Therefore, the inference goal here is to estimate: the number of possible jumps, its sizes, the coefficients and the shapes of the jump intensities. We propose a multiscale approach for the inference problem that can be seen as an indirect inference scheme. We found that using a Gaussian approximation based on moment expansions, it is possible to accurately estimate the jump intensities and the jump amplitudes. We obtained results equivalent to the state of the art but using a simpler and less expensive approach.

Improved convergence of gradient-based reconstruction using multi-scale models

International Nuclear Information System (INIS)

Cunningham, G.S.; Hanson, K.M.; Koyfman, I.

1996-01-01

Geometric models have received increasing attention in medical imaging for tasks such as segmentation, reconstruction, restoration, and registration. In order to determine the best configuration of the geometric model in the context of any of these tasks, one needs to perform a difficult global optimization of an energy function that may have many local minima. Explicit models of geometry, also called deformable models, snakes, or active contours, have been used extensively to solve image segmentation problems in a non-Bayesian framework. Researchers have seen empirically that multi-scale analysis is useful for convergence to a configuration that is near the global minimum. In this type of analysis, the image data are convolved with blur functions of increasing resolution, and an optimal configuration of the snake is found for each blurred image. The configuration obtained using the highest resolution blur is used as the solution to the global optimization problem. In this article, the authors use explicit models of geometry for a variety of Bayesian estimation problems, including image segmentation, reconstruction and restoration. The authors introduce a multi-scale approach that blurs the geometric model, rather than the image data, and show that this approach turns a global, highly nonquadratic optimization into a sequence of local, approximately quadratic problems that converge to the global minimum. The result is a deterministic, robust, and efficient optimization strategy applicable to a wide variety of Bayesian estimation problems in which geometric models of images are an important component

The aerosol distribution in Europe derived with the Community Multiscale Air Quality (CMAQ model: comparison to near surface in situ and sunphotometer measurements

Directory of Open Access Journals (Sweden)

V. Matthias

2008-09-01

Full Text Available The aerosol distribution in Europe was simulated with the Community Multiscale Air Quality (CMAQ model system version 4.5 for the years 2000 and 2001. The results were compared with daily averages of PM₁₀ measurements taken in the framework of EMEP and with aerosol optical depth (AOD values measured within AERONET. The modelled total aerosol mass is typically about 30–60% lower than the corresponding measurements. However a comparison of the chemical composition of the aerosol revealed a considerably better agreement between the modelled and the measured aerosol components for ammonium, nitrate and sulfate, which are on average only 15–20% underestimated. Sligthly worse agreement was determined for sea salt, that was only avaliable at two sites. The largest discrepancies result from the aerosol mass which was not chemically specified by the measurements. The agreement between measurements and model is better in winter than in summer. The modelled organic aerosol mass is higher in summer than in winter but it is significantly underestimated by the model. This could be one of the main reasons for the discrepancies between measurements and model results. The other is that primary coarse particles are underestimated in the emissions. The probability distribution function of the PM₁₀ measurements follows a log-normal distribution at most sites. The model is only able to reproduce this distribution function at non-coastal low altitude stations. The AOD derived from the model results is 20–70% lower than the values observed within AERONET. This is mainly attributed to the missing aerosol mass in the model. The day-to-day variability of the AOD and the log-normal distribution functions are quite well reproduced by the model. The seasonality on the other hand is underestimated by the model results because better agreement is achieved in winter.

Evaluation of the Community Multiscale Air Quality Model for Simulating Winter Ozone Formation in the Uinta Basin with Intensive Oil and Gas Production

Science.gov (United States)

Matichuk, R.; Tonnesen, G.; Luecken, D.; Roselle, S. J.; Napelenok, S. L.; Baker, K. R.; Gilliam, R. C.; Misenis, C.; Murphy, B.; Schwede, D. B.

2015-12-01

The western United States is an important source of domestic energy resources. One of the primary environmental impacts associated with oil and natural gas production is related to air emission releases of a number of air pollutants. Some of these pollutants are important precursors to the formation of ground-level ozone. To better understand ozone impacts and other air quality issues, photochemical air quality models are used to simulate the changes in pollutant concentrations in the atmosphere on local, regional, and national spatial scales. These models are important for air quality management because they assist in identifying source contributions to air quality problems and designing effective strategies to reduce harmful air pollutants. The success of predicting oil and natural gas air quality impacts depends on the accuracy of the input information, including emissions inventories, meteorological information, and boundary conditions. The treatment of chemical and physical processes within these models is equally important. However, given the limited amount of data collected for oil and natural gas production emissions in the past and the complex terrain and meteorological conditions in western states, the ability of these models to accurately predict pollution concentrations from these sources is uncertain. Therefore, this presentation will focus on understanding the Community Multiscale Air Quality (CMAQ) model's ability to predict air quality impacts associated with oil and natural gas production and its sensitivity to input uncertainties. The results will focus on winter ozone issues in the Uinta Basin, Utah and identify the factors contributing to model performance issues. The results of this study will help support future air quality model development, policy and regulatory decisions for the oil and gas sector.

Coarse-graining and hybrid methods for efficient simulation of stochastic multi-scale models of tumour growth

Science.gov (United States)

de la Cruz, Roberto; Guerrero, Pilar; Calvo, Juan; Alarcón, Tomás

2017-12-01

The development of hybrid methodologies is of current interest in both multi-scale modelling and stochastic reaction-diffusion systems regarding their applications to biology. We formulate a hybrid method for stochastic multi-scale models of cells populations that extends the remit of existing hybrid methods for reaction-diffusion systems. Such method is developed for a stochastic multi-scale model of tumour growth, i.e. population-dynamical models which account for the effects of intrinsic noise affecting both the number of cells and the intracellular dynamics. In order to formulate this method, we develop a coarse-grained approximation for both the full stochastic model and its mean-field limit. Such approximation involves averaging out the age-structure (which accounts for the multi-scale nature of the model) by assuming that the age distribution of the population settles onto equilibrium very fast. We then couple the coarse-grained mean-field model to the full stochastic multi-scale model. By doing so, within the mean-field region, we are neglecting noise in both cell numbers (population) and their birth rates (structure). This implies that, in addition to the issues that arise in stochastic-reaction diffusion systems, we need to account for the age-structure of the population when attempting to couple both descriptions. We exploit our coarse-graining model so that, within the mean-field region, the age-distribution is in equilibrium and we know its explicit form. This allows us to couple both domains consistently, as upon transference of cells from the mean-field to the stochastic region, we sample the equilibrium age distribution. Furthermore, our method allows us to investigate the effects of intracellular noise, i.e. fluctuations of the birth rate, on collective properties such as travelling wave velocity. We show that the combination of population and birth-rate noise gives rise to large fluctuations of the birth rate in the region at the leading edge of

Peridynamic Multiscale Finite Element Methods

Energy Technology Data Exchange (ETDEWEB)

Costa, Timothy [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Bond, Stephen D. [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Littlewood, David John [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Moore, Stan Gerald [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)

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

Development and application of a multiscale model for the magnetic fusion edge plasma region

International Nuclear Information System (INIS)

Hasenbeck, Felix Martin Michael

2016-01-01

remains diffusive. An analytic expression for the diffusion coefficient D is found, being consistent with results from percolation theory. The multiscale model and the coupling procedure are benchmarked with a one-dimensional test problem which consists of coupling the local version of the drift fluid code ATTEMPT to a 1D macroscale code to determine the time-dependent evolution of the flux surface averaged density in radial direction in the tokamak edge region. The reference simulation is given by a simulation of the non-local version of ATTEMPT, accounting for both the mesoscale and the macroscale evolution. Results of the coupled code simulations show that the macroscale evolution of the density and the radial particle flux can be reproduced with typical uncertainties of 6 and 22%, respectively. Time savings with respect to the non-local simulations are of a factor of ten for a typical macroscale evolution time of 10 milliseconds while a speedup factor of the order of 50 is achievable for an edge region with a radial extent of ∝ 30 cm and typical radial profile lengths of ∝5 cm. The multiscale model for two-dimensional, stationary problems is realized by coupling of the B2 code and the local version of the ATTEMPT code. The results of the corresponding coupled code simulations for experiments at the tokamak TEX TOR provide plasma profiles in agreement with experimental measurements with uncertainties regarding the electron density and electron temperature measured at the outer midplane around 10 to 25%. Poloidally and radially dependent profiles of the radial particle diffusion coefficients D, self-consistently determined by B2-ATTEMPT, have typical values of 0.3 to 0.9 m"2s"-"1 and are within a 10 to 30% range of effective diffusion coefficients employed in B2-EIRENE simulations to fit simulation results to measured density profiles. The poloidal dependence of D as given by the B2-ATTEMPT simulations accounts for the experimentally confirmed ballooning character

Development and application of a multiscale model for the magnetic fusion edge plasma region

Energy Technology Data Exchange (ETDEWEB)

Hasenbeck, Felix Martin Michael

2016-07-01

remains diffusive. An analytic expression for the diffusion coefficient D is found, being consistent with results from percolation theory. The multiscale model and the coupling procedure are benchmarked with a one-dimensional test problem which consists of coupling the local version of the drift fluid code ATTEMPT to a 1D macroscale code to determine the time-dependent evolution of the flux surface averaged density in radial direction in the tokamak edge region. The reference simulation is given by a simulation of the non-local version of ATTEMPT, accounting for both the mesoscale and the macroscale evolution. Results of the coupled code simulations show that the macroscale evolution of the density and the radial particle flux can be reproduced with typical uncertainties of 6 and 22%, respectively. Time savings with respect to the non-local simulations are of a factor of ten for a typical macroscale evolution time of 10 milliseconds while a speedup factor of the order of 50 is achievable for an edge region with a radial extent of ∝ 30 cm and typical radial profile lengths of ∝5 cm. The multiscale model for two-dimensional, stationary problems is realized by coupling of the B2 code and the local version of the ATTEMPT code. The results of the corresponding coupled code simulations for experiments at the tokamak TEX TOR provide plasma profiles in agreement with experimental measurements with uncertainties regarding the electron density and electron temperature measured at the outer midplane around 10 to 25%. Poloidally and radially dependent profiles of the radial particle diffusion coefficients D, self-consistently determined by B2-ATTEMPT, have typical values of 0.3 to 0.9 m{sup 2}s{sup -1} and are within a 10 to 30% range of effective diffusion coefficients employed in B2-EIRENE simulations to fit simulation results to measured density profiles. The poloidal dependence of D as given by the B2-ATTEMPT simulations accounts for the experimentally confirmed ballooning

Magnetic hysteresis at the domain scale of a multi-scale material model for magneto-elastic behaviour

Energy Technology Data Exchange (ETDEWEB)

Vanoost, D., E-mail: dries.vanoost@kuleuven-kulak.be [KU Leuven Technology Campus Ostend, ReMI Research Group, Oostende B-8400 (Belgium); KU Leuven Kulak, Wave Propagation and Signal Processing Research Group, Kortrijk B-8500 (Belgium); Steentjes, S. [Institute of Electrical Machines, RWTH Aachen University, Aachen D-52062 (Germany); Peuteman, J. [KU Leuven Technology Campus Ostend, ReMI Research Group, Oostende B-8400 (Belgium); KU Leuven, Department of Electrical Engineering, Electrical Energy and Computer Architecture, Heverlee B-3001 (Belgium); Gielen, G. [KU Leuven, Department of Electrical Engineering, Microelectronics and Sensors, Heverlee B-3001 (Belgium); De Gersem, H. [KU Leuven Kulak, Wave Propagation and Signal Processing Research Group, Kortrijk B-8500 (Belgium); TU Darmstadt, Institut für Theorie Elektromagnetischer Felder, Darmstadt D-64289 (Germany); Pissoort, D. [KU Leuven Technology Campus Ostend, ReMI Research Group, Oostende B-8400 (Belgium); KU Leuven, Department of Electrical Engineering, Microelectronics and Sensors, Heverlee B-3001 (Belgium); Hameyer, K. [Institute of Electrical Machines, RWTH Aachen University, Aachen D-52062 (Germany)

2016-09-15

This paper proposes a multi-scale energy-based material model for poly-crystalline materials. Describing the behaviour of poly-crystalline materials at three spatial scales of dominating physical mechanisms allows accounting for the heterogeneity and multi-axiality of the material behaviour. The three spatial scales are the poly-crystalline, grain and domain scale. Together with appropriate scale transitions rules and models for local magnetic behaviour at each scale, the model is able to describe the magneto-elastic behaviour (magnetostriction and hysteresis) at the macroscale, although the data input is merely based on a set of physical constants. Introducing a new energy density function that describes the demagnetisation field, the anhysteretic multi-scale energy-based material model is extended to the hysteretic case. The hysteresis behaviour is included at the domain scale according to the micro-magnetic domain theory while preserving a valid description for the magneto-elastic coupling. The model is verified using existing measurement data for different mechanical stress levels. - Highlights: • A ferromagnetic hysteretic energy-based multi-scale material model is proposed. • The hysteresis is obtained by new proposed hysteresis energy density function. • Avoids tedious parameter identification.

Quantification of structural uncertainties in multi-scale models; case study of the Lublin Basin, Poland

Science.gov (United States)

Małolepszy, Zbigniew; Szynkaruk, Ewa

2015-04-01

The multiscale static modeling of regional structure of the Lublin Basin is carried on in the Polish Geological Institute, in accordance with principles of integrated 3D geological modelling. The model is based on all available geospatial data from Polish digital databases and analogue archives. Mapped regional structure covers the area of 260x80 km located between Warsaw and Polish-Ukrainian border, along NW-SE-trending margin of the East European Craton. Within the basin, the Paleozoic beds with coalbearing Carboniferous and older formations containing hydrocarbons and unconventional prospects are covered unconformably by Permo-Mesozoic and younger rocks. Vertical extent of the regional model is set from topographic surface to 6000 m ssl and at the bottom includes some Proterozoic crystalline formations of the craton. The project focuses on internal consistency of the models built at different scales - from basin (small) scale to field-scale (large-scale). The models, nested in the common structural framework, are being constructed with regional geological knowledge, ensuring smooth transition in the 3D model resolution and amount of geological detail. Major challenge of the multiscale approach to subsurface modelling is the assessment and consistent quantification of various types of geological uncertainties tied to those various scale sub-models. Decreasing amount of information with depth and, particularly, very limited data collected below exploration targets, as well as accuracy and quality of data, all have the most critical impact on the modelled structure. In deeper levels of the Lublin Basin model, seismic interpretation of 2D surveys is sparsely tied to well data. Therefore time-to-depth conversion carries one of the major uncertainties in the modeling of structures, especially below 3000 m ssl. Furthermore, as all models at different scales are based on the same dataset, we must deal with different levels of generalization of geological structures. The

Revisiting of Multiscale Static Analysis of Notched Laminates Using the Generalized Method of Cells

Science.gov (United States)

Naghipour Ghezeljeh, Paria; Arnold, Steven M.; Pineda, Evan J.

2016-01-01

Composite material systems generally exhibit a range of behavior on different length scales (from constituent level to macro); therefore, a multiscale framework is beneficial for the design and engineering of these material systems. The complex nature of the observed composite failure during experiments suggests the need for a three-dimensional (3D) multiscale model to attain a reliable prediction. However, the size of a multiscale three-dimensional finite element model can become prohibitively large and computationally costly. Two-dimensional (2D) models are preferred due to computational efficiency, especially if many different configurations have to be analyzed for an in-depth damage tolerance and durability design study. In this study, various 2D and 3D multiscale analyses will be employed to conduct a detailed investigation into the tensile failure of a given multidirectional, notched carbon fiber reinforced polymer laminate. Threedimensional finite element analysis is typically considered more accurate than a 2D finite element model, as compared with experiments. Nevertheless, in the absence of adequate mesh refinement, large differences may be observed between a 2D and 3D analysis, especially for a shear-dominated layup. This observed difference has not been widely addressed in previous literature and is the main focus of this paper.

Cellular potts models multiscale extensions and biological applications

CERN Document Server

Scianna, Marco

2013-01-01

A flexible, cell-level, and lattice-based technique, the cellular Potts model accurately describes the phenomenological mechanisms involved in many biological processes. Cellular Potts Models: Multiscale Extensions and Biological Applications gives an interdisciplinary, accessible treatment of these models, from the original methodologies to the latest developments. The book first explains the biophysical bases, main merits, and limitations of the cellular Potts model. It then proposes several innovative extensions, focusing on ways to integrate and interface the basic cellular Potts model at the mesoscopic scale with approaches that accurately model microscopic dynamics. These extensions are designed to create a nested and hybrid environment, where the evolution of a biological system is realistically driven by the constant interplay and flux of information between the different levels of description. Through several biological examples, the authors demonstrate a qualitative and quantitative agreement with t...

SU-F-18C-15: Model-Based Multiscale Noise Reduction On Low Dose Cone Beam Projection

International Nuclear Information System (INIS)

Yao, W; Farr, J

2014-01-01

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

A multi-resolution assessment of the Community Multiscale Air Quality (CMAQ model v4.7 wet deposition estimates for 2002–2006

Directory of Open Access Journals (Sweden)

K. W. Appel

2011-05-01

Full Text Available This paper examines the operational performance of the Community Multiscale Air Quality (CMAQ model simulations for 2002–2006 using both 36-km and 12-km horizontal grid spacing, with a primary focus on the performance of the CMAQ model in predicting wet deposition of sulfate (SO₄⁼, ammonium (NH₄⁺ and nitrate (NO₃⁻. Performance of the wet deposition estimates from the model is determined by comparing CMAQ predicted concentrations to concentrations measured by the National Acid Deposition Program (NADP, specifically the National Trends Network (NTN. For SO₄⁼ wet deposition, the CMAQ model estimates were generally comparable between the 36-km and 12-km simulations for the eastern US, with the 12-km simulation giving slightly higher estimates of SO₄⁼ wet deposition than the 36-km simulation on average. The result is a slightly larger normalized mean bias (NMB for the 12-km simulation; however both simulations had annual biases that were less than ±15 % for each of the five years. The model estimated SO₄⁼ wet deposition values improved when they were adjusted to account for biases in the model estimated precipitation. The CMAQ model underestimates NH₄⁺ wet deposition over the eastern US, with a slightly larger underestimation in the 36-km simulation. The largest underestimations occur in the winter and spring periods, while the summer and fall have slightly smaller underestimations of NH₄⁺ wet deposition. The underestimation in NH₄⁺ wet deposition is likely due in part to the poor temporal and spatial representation of ammonia (NH₃ emissions, particularly those emissions associated with fertilizer applications and NH₃ bi-directional exchange. The model performance for estimates of NO₃⁻ wet deposition are

Multiscale virtual particle based elastic network model (MVP-ENM) for normal mode analysis of large-sized biomolecules.

Science.gov (United States)

Xia, Kelin

2017-12-20

In this paper, a multiscale virtual particle based elastic network model (MVP-ENM) is proposed for the normal mode analysis of large-sized biomolecules. The multiscale virtual particle (MVP) model is proposed for the discretization of biomolecular density data. With this model, large-sized biomolecular structures can be coarse-grained into virtual particles such that a balance between model accuracy and computational cost can be achieved. An elastic network is constructed by assuming "connections" between virtual particles. The connection is described by a special harmonic potential function, which considers the influence from both the mass distributions and distance relations of the virtual particles. Two independent models, i.e., the multiscale virtual particle based Gaussian network model (MVP-GNM) and the multiscale virtual particle based anisotropic network model (MVP-ANM), are proposed. It has been found that in the Debye-Waller factor (B-factor) prediction, the results from our MVP-GNM with a high resolution are as good as the ones from GNM. Even with low resolutions, our MVP-GNM can still capture the global behavior of the B-factor very well with mismatches predominantly from the regions with large B-factor values. Further, it has been demonstrated that the low-frequency eigenmodes from our MVP-ANM are highly consistent with the ones from ANM even with very low resolutions and a coarse grid. Finally, the great advantage of MVP-ANM model for large-sized biomolecules has been demonstrated by using two poliovirus virus structures. The paper ends with a conclusion.

Land-Atmosphere Coupling in the Multi-Scale Modelling Framework

Science.gov (United States)

Kraus, P. M.; Denning, S.

2015-12-01

The Multi-Scale Modeling Framework (MMF), in which cloud-resolving models (CRMs) are embedded within general circulation model (GCM) gridcells to serve as the model's cloud parameterization, has offered a number of benefits to GCM simulations. The coupling of these cloud-resolving models directly to land surface model instances, rather than passing averaged atmospheric variables to a single instance of a land surface model, the logical next step in model development, has recently been accomplished. This new configuration offers conspicuous improvements to estimates of precipitation and canopy through-fall, but overall the model exhibits warm surface temperature biases and low productivity.This work presents modifications to a land-surface model that take advantage of the new multi-scale modeling framework, and accommodate the change in spatial scale from a typical GCM range of ~200 km to the CRM grid-scale of 4 km.A parameterization is introduced to apportion modeled surface radiation into direct-beam and diffuse components. The diffuse component is then distributed among the land-surface model instances within each GCM cell domain. This substantially reduces the number excessively low light values provided to the land-surface model when cloudy conditions are modeled in the CRM, associated with its 1-D radiation scheme. The small spatial scale of the CRM, ~4 km, as compared with the typical ~200 km GCM scale, provides much more realistic estimates of precipitation intensity, this permits the elimination of a model parameterization of canopy through-fall. However, runoff at such scales can no longer be considered as an immediate flow to the ocean. Allowing sub-surface water flow between land-surface instances within the GCM domain affords better realism and also reduces temperature and productivity biases.The MMF affords a number of opportunities to land-surface modelers, providing both the advantages of direct simulation at the 4 km scale and a much reduced

Application of Mortar Coupling in Multiscale Modelling of Coupled Flow, Transport, and Biofilm Growth in Porous Media

Science.gov (United States)

Laleian, A.; Valocchi, A. J.; Werth, C. J.

2017-12-01

Multiscale models of reactive transport in porous media are capable of capturing complex pore-scale processes while leveraging the efficiency of continuum-scale models. In particular, porosity changes caused by biofilm development yield complex feedbacks between transport and reaction that are difficult to quantify at the continuum scale. Pore-scale models, needed to accurately resolve these dynamics, are often impractical for applications due to their computational cost. To address this challenge, we are developing a multiscale model of biofilm growth in which non-overlapping regions at pore and continuum spatial scales are coupled with a mortar method providing continuity at interfaces. We explore two decompositions of coupled pore-scale and continuum-scale regions to study biofilm growth in a transverse mixing zone. In the first decomposition, all reaction is confined to a pore-scale region extending the transverse mixing zone length. Only solute transport occurs in the surrounding continuum-scale regions. Relative to a fully pore-scale result, we find the multiscale model with this decomposition has a reduced run time and consistent result in terms of biofilm growth and solute utilization. In the second decomposition, reaction occurs in both an up-gradient pore-scale region and a down-gradient continuum-scale region. To quantify clogging, the continuum-scale model implements empirical relations between porosity and continuum-scale parameters, such as permeability and the transverse dispersion coefficient. Solutes are sufficiently mixed at the end of the pore-scale region, such that the initial reaction rate is accurately computed using averaged concentrations in the continuum-scale region. Relative to a fully pore-scale result, we find accuracy of biomass growth in the multiscale model with this decomposition improves as the interface between pore-scale and continuum-scale regions moves downgradient where transverse mixing is more fully developed. Also, this

Multi-time, multi-scale correlation functions in turbulence and in turbulent models

NARCIS (Netherlands)

Biferale, L.; Boffetta, G.; Celani, A.; Toschi, F.

1999-01-01

A multifractal-like representation for multi-time, multi-scale velocity correlation in turbulence and dynamical turbulent models is proposed. The importance of subleading contributions to time correlations is highlighted. The fulfillment of the dynamical constraints due to the equations of motion is

Multi-scale dynamic modeling of atmospheric pollution in urban environment

International Nuclear Information System (INIS)

Thouron, Laetitia

2017-01-01

Urban air pollution has been identified as an important cause of health impacts, including premature deaths. In particular, ambient concentrations of gaseous pollutants such as nitrogen dioxide (NO 2 ) and particulate matter (PM10 and PM2.5) are regulated, which means that emission reduction strategies must be put in place to reduce these concentrations in places where the corresponding regulations are not respected. Besides, air pollution can contribute to the contamination of other media, for example through the contribution of atmospheric deposition to runoff contamination. The multifactorial and multi-scale aspects of urban make the pollution sources difficult to identify. Indeed, the urban environment is a heterogeneous space characterized by complex architectural structures (old buildings alongside a more modern building, residential, commercial, industrial zones, roads, etc.), non-uniform atmospheric pollutant emissions and therefore the population exposure to pollution is variable in space and time. The modeling of urban air pollution aims to understand the origin of pollutants, their spatial extent and their concentration/deposition levels. Some pollutants have long residence times and can stay several weeks in the atmosphere (PM2.5) and therefore be transported over long distances, while others are more local (NO x in the vicinity of traffic). The spatial distribution of a pollutant will therefore depend on several factors, and in particular on the surfaces encountered. Air quality depends strongly on weather, buildings (canyon-street) and emissions. The aim of this thesis is to address some of these aspects by modeling: (1) urban background pollution with a transport-chemical model (Polyphemus / POLAIR3D), which makes it possible to estimate atmospheric pollutants by type of urban surfaces (roofs, walls and roadways), (2) street-level pollution by explicitly integrating the effects of the building in a three-dimensional way with a multi-scale model of

Multiscale Models for the Two-Stream Instability

Science.gov (United States)

Joseph, Ilon; Dimits, Andris; Banks, Jeffrey; Berger, Richard; Brunner, Stephan; Chapman, Thomas

2017-10-01

Interpenetrating streams of plasma found in many important scenarios in nature and in the laboratory can develop kinetic two-stream instabilities that exchange momentum and energy between the streams. A quasilinear model for the electrostatic two-stream instability is under development as a component of a multiscale model that couples fluid simulations to kinetic theory. Parameters of the model will be validated with comparison to full kinetic simulations using LOKI and efficient strategies for numerical solution of the quasilinear model and for coupling to the fluid model will be discussed. Extending the kinetic models into the collisional regime requires an efficient treatment of the collision operator. Useful reductions of the collision operator relative to the full multi-species Landau-Fokker-Plank operator are being explored. These are further motivated both by careful consideration of the parameter orderings relevant to two-stream scenarios and by the particular 2D+2V phase space used in the LOKI code. Prepared for US DOE by LLNL under Contract DE-AC52-07NA27344 and LDRD project 17- ERD-081.

Review on modeling development for multiscale chemical reactions coupled transport phenomena in solid oxide fuel cells

International Nuclear Information System (INIS)

Andersson, Martin; Yuan, Jinliang; Sunden, Bengt

2010-01-01

A literature study is performed to compile the state-of-the-art, as well as future potential, in SOFC modeling. Principles behind various transport processes such as mass, heat, momentum and charge as well as for electrochemical and internal reforming reactions are described. A deeper investigation is made to find out potentials and challenges using a multiscale approach to model solid oxide fuel cells (SOFCs) and combine the accuracy at microscale with the calculation speed at macroscale to design SOFCs, based on a clear understanding of transport phenomena, chemical reactions and functional requirements. Suitable methods are studied to model SOFCs covering various length scales. Coupling methods between different approaches and length scales by multiscale models are outlined. Multiscale modeling increases the understanding for detailed transport phenomena, and can be used to make a correct decision on the specific design and control of operating conditions. It is expected that the development and production costs will be decreased and the energy efficiency be increased (reducing running cost) as the understanding of complex physical phenomena increases. It is concluded that the connection between numerical modeling and experiments is too rare and also that material parameters in most cases are valid only for standard materials and not for the actual SOFC component microstructures.

Cost and time-effective method for multi-scale measures of rugosity, fractal dimension, and vector dispersion from coral reef 3D models.

Directory of Open Access Journals (Sweden)

G C Young

Full Text Available We present a method to construct and analyse 3D models of underwater scenes using a single cost-effective camera on a standard laptop with (a free or low-cost software, (b no computer programming ability, and (c minimal man hours for both filming and analysis. This study focuses on four key structural complexity metrics: point-to-point distances, linear rugosity (R, fractal dimension (D, and vector dispersion (1/k. We present the first assessment of accuracy and precision of structure-from-motion (SfM 3D models from an uncalibrated GoPro™ camera at a small scale (4 m2 and show that they can provide meaningful, ecologically relevant results. Models had root mean square errors of 1.48 cm in X-Y and 1.35 in Z, and accuracies of 86.8% (R, 99.6% (D at scales 30-60 cm, 93.6% (D at scales 1-5 cm, and 86.9 (1/k. Values of R were compared to in-situ chain-and-tape measurements, while values of D and 1/k were compared with ground truths from 3D printed objects modelled underwater. All metrics varied less than 3% between independently rendered models. We thereby improve and rigorously validate a tool for ecologists to non-invasively quantify coral reef structural complexity with a variety of multi-scale metrics.

Multiscale Persistent Functions for Biomolecular Structure Characterization

Energy Technology Data Exchange (ETDEWEB)

Xia, Kelin [Nanyang Technological University (Singapore). Division of Mathematical Sciences, School of Physical, Mathematical Sciences and School of Biological Sciences; Li, Zhiming [Central China Normal University, Wuhan (China). Key Laboratory of Quark and Lepton Physics (MOE) and Institute of Particle Physics; Mu, Lin [Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States). Computer Science and Mathematics Division

2017-11-02

Here in this paper, we introduce multiscale persistent functions for biomolecular structure characterization. The essential idea is to combine our multiscale rigidity functions (MRFs) with persistent homology analysis, so as to construct a series of multiscale persistent functions, particularly multiscale persistent entropies, for structure characterization. To clarify the fundamental idea of our method, the multiscale persistent entropy (MPE) model is discussed in great detail. Mathematically, unlike the previous persistent entropy (Chintakunta et al. in Pattern Recognit 48(2):391–401, 2015; Merelli et al. in Entropy 17(10):6872–6892, 2015; Rucco et al. in: Proceedings of ECCS 2014, Springer, pp 117–128, 2016), a special resolution parameter is incorporated into our model. Various scales can be achieved by tuning its value. Physically, our MPE can be used in conformational entropy evaluation. More specifically, it is found that our method incorporates in it a natural classification scheme. This is achieved through a density filtration of an MRF built from angular distributions. To further validate our model, a systematical comparison with the traditional entropy evaluation model is done. Additionally, it is found that our model is able to preserve the intrinsic topological features of biomolecular data much better than traditional approaches, particularly for resolutions in the intermediate range. Moreover, by comparing with traditional entropies from various grid sizes, bond angle-based methods and a persistent homology-based support vector machine method (Cang et al. in Mol Based Math Biol 3:140–162, 2015), we find that our MPE method gives the best results in terms of average true positive rate in a classic protein structure classification test. More interestingly, all-alpha and all-beta protein classes can be clearly separated from each other with zero error only in our model. Finally, a special protein structure index (PSI) is proposed, for the first

Modeling ozone and aerosol formation and transport in the pacific northwest with the community Multi-Scale Air Quality (CMAQ) modeling system.

Science.gov (United States)

O'Neill, Susan M; Lamb, Brian K; Chen, Jack; Claiborn, Candis; Finn, Dennis; Otterson, Sally; Figueroa, Cristiana; Bowman, Clint; Boyer, Mike; Wilson, Rob; Arnold, Jeff; Aalbers, Steven; Stocum, Jeffrey; Swab, Christopher; Stoll, Matt; Dubois, Mike; Anderson, Mary

2006-02-15

The Community Multi-Scale Air Quality (CMAQ) modeling system was used to investigate ozone and aerosol concentrations in the Pacific Northwest (PNW) during hot summertime conditions during July 1-15, 1996. Two emission inventories (El) were developed: emissions for the first El were based upon the National Emission Trend 1996 (NET96) database and the BEIS2 biogenic emission model, and emissions for the second El were developed through a "bottom up" approach that included biogenic emissions obtained from the GLOBEIS model. The two simulations showed that elevated PM2.5 concentrations occurred near and downwind of the Interstate-5 corridor along the foothills of the Cascade Mountains and in forested areas of central Idaho. The relative contributions of organic and inorganic aerosols varied by region, but generally organic aerosols constituted the largest fraction of PM2.5. In wilderness areas near the 1-5 corridor, organic carbon from anthropogenic sources contributed approximately 50% of the total organic carbon with the remainder from biogenic precursors, while in wilderness areas in Idaho, biogenic organic carbon accounted for 80% of the total organic aerosol. Regional analysis of the secondary organic aerosol formation in the Columbia River Gorge, Central Idaho, and the Olympics/Puget Sound showed that the production rate of secondary organic carbon depends on local terpene concentrations and the local oxidizing capacity of the atmosphere, which was strongly influenced by anthropogenic emissions. Comparison with observations from 12 IMPROVE sites and 21 ozone monitoring sites showed that results from the two El simulations generally bracketed the average observed PM parameters and that errors calculated for the model results were within acceptable bounds. Analysis across all statistical parameters indicated that the NW-AIRQUEST El solution performed better at predicting PM2.5, PM1, and beta(ext) even though organic carbon PM was over-predicted, and the NET96 El

Multiscale modeling of alloy solidification using a database approach

Science.gov (United States)

Tan, Lijian; Zabaras, Nicholas

2007-11-01

A two-scale model based on a database approach is presented to investigate alloy solidification. Appropriate assumptions are introduced to describe the behavior of macroscopic temperature, macroscopic concentration, liquid volume fraction and microstructure features. These assumptions lead to a macroscale model with two unknown functions: liquid volume fraction and microstructure features. These functions are computed using information from microscale solutions of selected problems. This work addresses the selection of sample problems relevant to the interested problem and the utilization of data from the microscale solution of the selected sample problems. A computationally efficient model, which is different from the microscale and macroscale models, is utilized to find relevant sample problems. In this work, the computationally efficient model is a sharp interface solidification model of a pure material. Similarities between the sample problems and the problem of interest are explored by assuming that the liquid volume fraction and microstructure features are functions of solution features extracted from the solution of the computationally efficient model. The solution features of the computationally efficient model are selected as the interface velocity and thermal gradient in the liquid at the time the sharp solid-liquid interface passes through. An analytical solution of the computationally efficient model is utilized to select sample problems relevant to solution features obtained at any location of the domain of the problem of interest. The microscale solution of selected sample problems is then utilized to evaluate the two unknown functions (liquid volume fraction and microstructure features) in the macroscale model. The temperature solution of the macroscale model is further used to improve the estimation of the liquid volume fraction and microstructure features. Interpolation is utilized in the feature space to greatly reduce the number of required

Performance and diagnostic evaluation of ozone predictions by the Eta-Community Multiscale Air Quality Forecast System during the 2002 New England Air Quality Study.

Science.gov (United States)

Yu, Shaocai; Mathur, Rohit; Kang, Daiwen; Schere, Kenneth; Eder, Brian; Pleim, Jonathan

2006-10-01

A real-time air quality forecasting system (Eta-Community Multiscale Air Quality [CMAQ] model suite) has been developed by linking the National Centers for Environmental Estimation Eta model to the U.S. Environmental Protection Agency (EPA) CMAQ model. This work presents results from the application of the Eta-CMAQ modeling system for forecasting ozone (O3) over the Northeastern United States during the 2002 New England Air Quality Study (NEAQS). Spatial and temporal performance of the Eta-CMAQ model for O3 was evaluated by comparison with observations from the EPA Air Quality System (AQS) network. This study also examines the ability of the model to simulate the processes governing the distributions of tropospheric O3 on the basis of the intensive datasets obtained at the four Atmospheric Investigation, Regional Modeling, Analysis, and Estimation (AIRMAP) and Harvard Forest (HF) surface sites. The episode analysis reveals that the model captured the buildup of O3 concentrations over the northeastern domain from August 11 and reproduced the spatial distributions of observed O3 very well for the daytime (8:00 p.m.) of both August 8 and 12 with most of normalized mean bias (NMB) within +/- 20%. The model reproduced 53.3% of the observed hourly O3 within a factor of 1.5 with NMB of 29.7% and normalized mean error of 46.9% at the 342 AQS sites. The comparison of modeled and observed lidar O3 vertical profiles shows that whereas the model reproduced the observed vertical structure, it tended to overestimate at higher altitude. The model reproduced 64-77% of observed NO2 photolysis rate values within a factor of 1.5 at the AIRMAP sites. At the HF site, comparison of modeled and observed O3/nitrogen oxide (NOx) ratios suggests that the site is mainly under strongly NOx-sensitive conditions (>53%). It was found that the modeled lower limits of the O3 production efficiency values (inferred from O3-CO correlation) are close to the observations.

Multi-scale modelling and numerical simulation of electronic kinetic transport

International Nuclear Information System (INIS)

Duclous, R.

2009-11-01

This research thesis which is at the interface between numerical analysis, plasma physics and applied mathematics, deals with the kinetic modelling and numerical simulations of the electron energy transport and deposition in laser-produced plasmas, having in view the processes of fuel assembly to temperature and density conditions necessary to ignite fusion reactions. After a brief review of the processes at play in the collisional kinetic theory of plasmas, with a focus on basic models and methods to implement, couple and validate them, the author focuses on the collective aspect related to the free-streaming electron transport equation in the non-relativistic limit as well as in the relativistic regime. He discusses the numerical development and analysis of the scheme for the Vlasov-Maxwell system, and the selection of a validation procedure and numerical tests. Then, he investigates more specific aspects of the collective transport: the multi-specie transport, submitted to phase-space discontinuities. Dealing with the multi-scale physics of electron transport with collision source terms, he validates the accuracy of a fast Monte Carlo multi-grid solver for the Fokker-Planck-Landau electron-electron collision operator. He reports realistic simulations for the kinetic electron transport in the frame of the shock ignition scheme, the development and validation of a reduced electron transport angular model. He finally explores the relative importance of the processes involving electron-electron collisions at high energy by means a multi-scale reduced model with relativistic Boltzmann terms

Multiscale crystal defect dynamics: A coarse-grained lattice defect model based on crystal microstructure

Science.gov (United States)

Lyu, Dandan; Li, Shaofan

2017-10-01

Crystal defects have microstructure, and this microstructure should be related to the microstructure of the original crystal. Hence each type of crystals may have similar defects due to the same failure mechanism originated from the same microstructure, if they are under the same loading conditions. In this work, we propose a multiscale crystal defect dynamics (MCDD) model that models defects by considering its intrinsic microstructure derived from the microstructure or material genome of the original perfect crystal. The main novelties of present work are: (1) the discrete exterior calculus and algebraic topology theory are used to construct a scale-up (coarse-grained) dual lattice model for crystal defects, which may represent all possible defect modes inside a crystal; (2) a higher order Cauchy-Born rule (up to the fourth order) is adopted to construct atomistic-informed constitutive relations for various defect process zones, and (3) an hierarchical strain gradient theory based finite element formulation is developed to support an hierarchical multiscale cohesive (process) zone model for various defects in a unified formulation. The efficiency of MCDD computational algorithm allows us to simulate dynamic defect evolution at large scale while taking into account atomistic interaction. The MCDD model has been validated by comparing of the results of MCDD simulations with that of molecular dynamics (MD) in the cases of nanoindentation and uniaxial tension. Numerical simulations have shown that MCDD model can predict dislocation nucleation induced instability and inelastic deformation, and thus it may provide an alternative solution to study crystal plasticity.

Sustainable Manufacturing via Multi-Scale, Physics-Based Process Modeling and Manufacturing-Informed Design

Energy Technology Data Exchange (ETDEWEB)

None

2017-04-01

This factsheet describes a project that developed and demonstrated a new manufacturing-informed design framework that utilizes advanced multi-scale, physics-based process modeling to dramatically improve manufacturing productivity and quality in machining operations while reducing the cost of machined components.

A multiscale and multiphysics numerical framework for modelling of hygrothermal ageing in laminated composites

NARCIS (Netherlands)

Barcelos Carneiro M Rocha, Iuri; van der Meer, F.P.; Nijssen, RPL; Sluijs, Bert

2017-01-01

In this work, a numerical framework for modelling of hygrothermal ageing in laminated composites is proposed. The model consists of a macroscopic diffusion analysis based on Fick's second law coupled with a multiscale FE² stress analysis in order to take microscopic degradation

Evaluating and Improving Cloud Processes in the Multi-Scale Modeling Framework

Energy Technology Data Exchange (ETDEWEB)

Ackerman, Thomas P. [Univ. of Washington, Seattle, WA (United States)

2015-03-01

The research performed under this grant was intended to improve the embedded cloud model in the Multi-scale Modeling Framework (MMF) for convective clouds by using a 2-moment microphysics scheme rather than the single moment scheme used in all the MMF runs to date. The technical report and associated documents describe the results of testing the cloud resolving model with fixed boundary conditions and evaluation of model results with data. The overarching conclusion is that such model evaluations are problematic because errors in the forcing fields control the results so strongly that variations in parameterization values cannot be usefully constrained

EPOS Multi-Scale Laboratory platform: a long-term reference tool for experimental Earth Sciences

Science.gov (United States)

Trippanera, Daniele; Tesei, Telemaco; Funiciello, Francesca; Sagnotti, Leonardo; Scarlato, Piergiorgio; Rosenau, Matthias; Elger, Kirsten; Ulbricht, Damian; Lange, Otto; Calignano, Elisa; Spiers, Chris; Drury, Martin; Willingshofer, Ernst; Winkler, Aldo

2017-04-01

With continuous progress on scientific research, a large amount of datasets has been and will be produced. The data access and sharing along with their storage and homogenization within a unique and coherent framework is a new challenge for the whole scientific community. This is particularly emphasized for geo-scientific laboratories, encompassing the most diverse Earth Science disciplines and typology of data. To this aim the "Multiscale Laboratories" Work Package (WP16), operating in the framework of the European Plate Observing System (EPOS), is developing a virtual platform of geo-scientific data and services for the worldwide community of laboratories. This long-term project aims at merging the top class multidisciplinary laboratories in Geoscience into a coherent and collaborative network, facilitating the standardization of virtual access to data, data products and software. This will help our community to evolve beyond the stage in which most of data produced by the different laboratories are available only within the related scholarly publications (often as print-version only) or they remain unpublished and inaccessible on local devices. The EPOS multi-scale laboratory platform will provide the possibility to easily share and discover data by means of open access, DOI-referenced, online data publication including long-term storage, managing and curation services and to set up a cohesive community of laboratories. The WP16 is starting with three pilot cases laboratories: (1) rock physics, (2) palaeomagnetic, and (3) analogue modelling. As a proof of concept, first analogue modelling datasets have been published via GFZ Data Services (http://doidb.wdc-terra.org/search/public/ui?&sort=updated+desc&q=epos). The datasets include rock analogue material properties (e.g. friction data, rheology data, SEM imagery), as well as supplementary figures, images and movies from experiments on tectonic processes. A metadata catalogue tailored to the specific communities

Predicting FLDs Using a Multiscale Modeling Scheme

Science.gov (United States)

Wu, Z.; Loy, C.; Wang, E.; Hegadekatte, V.

2017-09-01

The measurement of a single forming limit diagram (FLD) requires significant resources and is time consuming. We have developed a multiscale modeling scheme to predict FLDs using a combination of limited laboratory testing, crystal plasticity (VPSC) modeling, and dual sequential-stage finite element (ABAQUS/Explicit) modeling with the Marciniak-Kuczynski (M-K) criterion to determine the limit strain. We have established a means to work around existing limitations in ABAQUS/Explicit by using an anisotropic yield locus (e.g., BBC2008) in combination with the M-K criterion. We further apply a VPSC model to reduce the number of laboratory tests required to characterize the anisotropic yield locus. In the present work, we show that the predicted FLD is in excellent agreement with the measured FLD for AA5182 in the O temper. Instead of 13 different tests as for a traditional FLD determination within Novelis, our technique uses just four measurements: tensile properties in three orientations; plane strain tension; biaxial bulge; and the sheet crystallographic texture. The turnaround time is consequently far less than for the traditional laboratory measurement of the FLD.

Multiscale modeling of polyisoprene on graphite

International Nuclear Information System (INIS)

Pandey, Yogendra Narayan; Brayton, Alexander; Doxastakis, Manolis; Burkhart, Craig; Papakonstantopoulos, George J.

2014-01-01

The local dynamics and the conformational properties of polyisoprene next to a smooth graphite surface constructed by graphene layers are studied by a multiscale methodology. First, fully atomistic molecular dynamics simulations of oligomers next to the surface are performed. Subsequently, Monte Carlo simulations of a systematically derived coarse-grained model generate numerous uncorrelated structures for polymer systems. A new reverse backmapping strategy is presented that reintroduces atomistic detail. Finally, multiple extensive fully atomistic simulations with large systems of long macromolecules are employed to examine local dynamics in proximity to graphite. Polyisoprene repeat units arrange close to a parallel configuration with chains exhibiting a distribution of contact lengths. Efficient Monte Carlo algorithms with the coarse-grain model are capable of sampling these distributions for any molecular weight in quantitative agreement with predictions from atomistic models. Furthermore, molecular dynamics simulations with well-equilibrated systems at all length-scales support an increased dynamic heterogeneity that is emerging from both intermolecular interactions with the flat surface and intramolecular cooperativity. This study provides a detailed comprehensive picture of polyisoprene on a flat surface and consists of an effort to characterize such systems in atomistic detail

Experimentally-based multiscale model of the elastic moduli of bovine trabecular bone and its constituents

Energy Technology Data Exchange (ETDEWEB)

Hamed, Elham [University of Illinois at Urbana-Champaign, Department of Mechanical Science and Engineering, 1206 West Green Street, Urbana, IL 61801 (United States); Novitskaya, Ekaterina, E-mail: eevdokim@ucsd.edu [University of California, San Diego, Department of Mechanical and Aerospace Engineering, Materials Science and Engineering Program, 9500 Gilman Dr., La Jolla, CA 92093 (United States); Li, Jun; Jasiuk, Iwona [University of Illinois at Urbana-Champaign, Department of Mechanical Science and Engineering, 1206 West Green Street, Urbana, IL 61801 (United States); McKittrick, Joanna [University of California, San Diego, Department of Mechanical and Aerospace Engineering, Materials Science and Engineering Program, 9500 Gilman Dr., La Jolla, CA 92093 (United States)

2015-09-01

The elastic moduli of trabecular bone were modeled using an analytical multiscale approach. Trabecular bone was represented as a porous nanocomposite material with a hierarchical structure spanning from the collagen–mineral level to the trabecular architecture level. In parallel, compression testing was done on bovine femoral trabecular bone samples in two anatomical directions, parallel to the femoral neck axis and perpendicular to it, and the measured elastic moduli were compared with the corresponding theoretical results. To gain insights on the interaction of collagen and minerals at the nanoscale, bone samples were deproteinized or demineralized. After such processing, the treated samples remained as self-standing structures and were tested in compression. Micro-computed tomography was used to characterize the hierarchical structure of these three bone types and to quantify the amount of bone porosity. The obtained experimental data served as inputs to the multiscale model and guided us to represent bone as an interpenetrating composite material. Good agreement was found between the theory and experiments for the elastic moduli of the untreated, deproteinized, and demineralized trabecular bone. - Highlights: • A multiscale model was used to predict the elastic moduli of trabecular bone. • Samples included demineralized, deproteinized and untreated bone. • The model portrays bone as a porous, interpenetrating two phase composite. • The experimental elastic moduli for trabecular bone fell between theoretical bounds.

Experimentally-based multiscale model of the elastic moduli of bovine trabecular bone and its constituents

International Nuclear Information System (INIS)

Hamed, Elham; Novitskaya, Ekaterina; Li, Jun; Jasiuk, Iwona; McKittrick, Joanna

2015-01-01

The elastic moduli of trabecular bone were modeled using an analytical multiscale approach. Trabecular bone was represented as a porous nanocomposite material with a hierarchical structure spanning from the collagen–mineral level to the trabecular architecture level. In parallel, compression testing was done on bovine femoral trabecular bone samples in two anatomical directions, parallel to the femoral neck axis and perpendicular to it, and the measured elastic moduli were compared with the corresponding theoretical results. To gain insights on the interaction of collagen and minerals at the nanoscale, bone samples were deproteinized or demineralized. After such processing, the treated samples remained as self-standing structures and were tested in compression. Micro-computed tomography was used to characterize the hierarchical structure of these three bone types and to quantify the amount of bone porosity. The obtained experimental data served as inputs to the multiscale model and guided us to represent bone as an interpenetrating composite material. Good agreement was found between the theory and experiments for the elastic moduli of the untreated, deproteinized, and demineralized trabecular bone. - Highlights: • A multiscale model was used to predict the elastic moduli of trabecular bone. • Samples included demineralized, deproteinized and untreated bone. • The model portrays bone as a porous, interpenetrating two phase composite. • The experimental elastic moduli for trabecular bone fell between theoretical bounds

A multi-scale energy demand model suggests sharing market risks with intelligent energy cooperatives

NARCIS (Netherlands)

G. Methenitis (Georgios); M. Kaisers (Michael); J.A. La Poutré (Han)

2015-01-01

textabstractIn this paper, we propose a multi-scale model of energy demand that is consistent with observations at a macro scale, in our use-case standard load profiles for (residential) electric loads. We employ the model to study incentives to assume the risk of volatile market prices for

Multiscale Modeling of Mesoscale and Interfacial Phenomena

Science.gov (United States)

Petsev, Nikolai Dimitrov

we provide a novel and general framework for multiscale modeling of systems featuring one or more dissolved species. This makes it possible to retain molecular detail for parts of the problem that require it while using a simple, continuum description for parts where high detail is unnecessary, reducing the number of degrees of freedom (i.e. number of particles) dramatically. This opens the possibility for modeling ion transport in biological processes and biomolecule assembly in ionic solution, as well as electrokinetic phenomena at interfaces such as corrosion. The number of particles in the system is further reduced through an integrated boundary approach, which we apply to colloidal suspensions. In this thesis, we describe this general framework for multiscale modeling single- and multicomponent systems, provide several simple equilibrium and non-equilibrium case studies, and discuss future applications.

A Multi-Scale Energy Food Systems Modeling Framework For Climate Adaptation

Science.gov (United States)

Siddiqui, S.; Bakker, C.; Zaitchik, B. F.; Hobbs, B. F.; Broaddus, E.; Neff, R.; Haskett, J.; Parker, C.

2016-12-01

Our goal is to understand coupled system dynamics across scales in a manner that allows us to quantify the sensitivity of critical human outcomes (nutritional satisfaction, household economic well-being) to development strategies and to climate or market induced shocks in sub-Saharan Africa. We adopt both bottom-up and top-down multi-scale modeling approaches focusing our efforts on food, energy, water (FEW) dynamics to define, parameterize, and evaluate modeled processes nationally as well as across climate zones and communities. Our framework comprises three complementary modeling techniques spanning local, sub-national and national scales to capture interdependencies between sectors, across time scales, and on multiple levels of geographic aggregation. At the center is a multi-player micro-economic (MME) partial equilibrium model for the production, consumption, storage, and transportation of food, energy, and fuels, which is the focus of this presentation. We show why such models can be very useful for linking and integrating across time and spatial scales, as well as a wide variety of models including an agent-based model applied to rural villages and larger population centers, an optimization-based electricity infrastructure model at a regional scale, and a computable general equilibrium model, which is applied to understand FEW resources and economic patterns at national scale. The MME is based on aggregating individual optimization problems for relevant players in an energy, electricity, or food market and captures important food supply chain components of trade and food distribution accounting for infrastructure and geography. Second, our model considers food access and utilization by modeling food waste and disaggregating consumption by income and age. Third, the model is set up to evaluate the effects of seasonality and system shocks on supply, demand, infrastructure, and transportation in both energy and food.

3D multiscale crack propagation using the XFEM applied to a gas turbine blade

Science.gov (United States)

Holl, Matthias; Rogge, Timo; Loehnert, Stefan; Wriggers, Peter; Rolfes, Raimund

2014-01-01

This work presents a new multiscale technique to investigate advancing cracks in three dimensional space. This fully adaptive multiscale technique is designed to take into account cracks of different length scales efficiently, by enabling fine scale domains locally in regions of interest, i.e. where stress concentrations and high stress gradients occur. Due to crack propagation, these regions change during the simulation process. Cracks are modeled using the extended finite element method, such that an accurate and powerful numerical tool is achieved. Restricting ourselves to linear elastic fracture mechanics, the -integral yields an accurate solution of the stress intensity factors, and with the criterion of maximum hoop stress, a precise direction of growth. If necessary, the on the finest scale computed crack surface is finally transferred to the corresponding scale. In a final step, the model is applied to a quadrature point of a gas turbine blade, to compute crack growth on the microscale of a real structure.

Physics-based hybrid method for multiscale transport in porous media

Science.gov (United States)

Yousefzadeh, Mehrdad; Battiato, Ilenia

2017-09-01

Despite advancements in the development of multiscale models for flow and reactive transport in porous media, the accurate, efficient and physics-based coupling of multiple scales in hybrid models remains a major theoretical and computational challenge. Improving the predictivity of macroscale predictions by means of multiscale algorithms relative to classical at-scale models is the primary motivation for the development of multiscale simulators. Yet, very few are the quantitative studies that explicitly address the predictive capability of multiscale coupling algorithms as it is still generally not possible to have a priori estimates of the errors that are present when complex flow processes are modeled. We develop a nonintrusive pore-/continuum-scale hybrid model whose coupling error is bounded by the upscaling error, i.e. we build a predictive tightly coupled multiscale scheme. This is accomplished by slightly enlarging the subdomain where continuum-scale equations are locally invalid and analytically defining physics-based coupling conditions at the interfaces separating the two computational sub-domains, while enforcing state variable and flux continuity. The proposed multiscale coupling approach retains the advantages of domain decomposition approaches, including the use of existing solvers for each subdomain, while it gains flexibility in the choice of the numerical discretization method and maintains the coupling errors bounded by the upscaling error. We implement the coupling in finite volumes and test the proposed method by modeling flow and transport through a reactive channel and past an array of heterogeneously reactive cylinders.

Multiscale Modeling of UHTC: Thermal Conductivity

Science.gov (United States)

Lawson, John W.; Murry, Daw; Squire, Thomas; Bauschlicher, Charles W.

2012-01-01

We are developing a multiscale framework in computational modeling for the ultra high temperature ceramics (UHTC) ZrB2 and HfB2. These materials are characterized by high melting point, good strength, and reasonable oxidation resistance. They are candidate materials for a number of applications in extreme environments including sharp leading edges of hypersonic aircraft. In particular, we used a combination of ab initio methods, atomistic simulations and continuum computations to obtain insights into fundamental properties of these materials. Ab initio methods were used to compute basic structural, mechanical and thermal properties. From these results, a database was constructed to fit a Tersoff style interatomic potential suitable for atomistic simulations. These potentials were used to evaluate the lattice thermal conductivity of single crystals and the thermal resistance of simple grain boundaries. Finite element method (FEM) computations using atomistic results as inputs were performed with meshes constructed on SEM images thereby modeling the realistic microstructure. These continuum computations showed the reduction in thermal conductivity due to the grain boundary network.

Multiscale methods in computational fluid and solid mechanics

NARCIS (Netherlands)

Borst, de R.; Hulshoff, S.J.; Lenz, S.; Munts, E.A.; Brummelen, van E.H.; Wall, W.; Wesseling, P.; Onate, E.; Periaux, J.

2006-01-01

First, an attempt is made towards gaining a more systematic understanding of recent progress in multiscale modelling in computational solid and fluid mechanics. Sub- sequently, the discussion is focused on variational multiscale methods for the compressible and incompressible Navier-Stokes

Multiscale thermohydrologic model: addressing variability and uncertainty at Yucca Mountain

International Nuclear Information System (INIS)

Buscheck, T; Rosenberg, N D; Gansemer, J D; Sun, Y

2000-01-01

Performance assessment and design evaluation require a modeling tool that simultaneously accounts for processes occurring at a scale of a few tens of centimeters around individual waste packages and emplacement drifts, and also on behavior at the scale of the mountain. Many processes and features must be considered, including non-isothermal, multiphase-flow in rock of variable saturation and thermal radiation in open cavities. Also, given the nature of the fractured rock at Yucca Mountain, a dual-permeability approach is needed to represent permeability. A monolithic numerical model with all these features requires too large a computational cost to be an effective simulation tool, one that is used to examine sensitivity to key model assumptions and parameters. We have developed a multi-scale modeling approach that effectively simulates 3D discrete-heat-source, mountain-scale thermohydrologic behavior at Yucca Mountain and captures the natural variability of the site consistent with what we know from site characterization and waste-package-to-waste-package variability in heat output. We describe this approach and present results examining the role of infiltration flux, the most important natural-system parameter with respect to how thermohydrologic behavior influences the performance of the repository

Stochastic scalar mixing models accounting for turbulent frequency multiscale fluctuations

International Nuclear Information System (INIS)

Soulard, Olivier; Sabel'nikov, Vladimir; Gorokhovski, Michael

2004-01-01

Two new scalar micromixing models accounting for a turbulent frequency scale distribution are investigated. These models were derived by Sabel'nikov and Gorokhovski [Second International Symposium on Turbulence and Shear FLow Phenomena, Royal Institute of technology (KTH), Stockholm, Sweden, June 27-29, 2001] using a multiscale extension of the classical interaction by exchange with the mean (IEM) and Langevin models. They are, respectively, called Extended IEM (EIEM) and Extended Langevin (ELM) models. The EIEM and ELM models are tested against DNS results in the case of the decay of a homogeneous scalar field in homogeneous turbulence. This comparison leads to a reformulation of the law governing the mixing frequency distribution. Finally, the asymptotic behaviour of the modeled PDF is discussed

Training Systems Modelers through the Development of a Multi-scale Chagas Disease Risk Model

Science.gov (United States)

Hanley, J.; Stevens-Goodnight, S.; Kulkarni, S.; Bustamante, D.; Fytilis, N.; Goff, P.; Monroy, C.; Morrissey, L. A.; Orantes, L.; Stevens, L.; Dorn, P.; Lucero, D.; Rios, J.; Rizzo, D. M.

2012-12-01

The goal of our NSF-sponsored Division of Behavioral and Cognitive Sciences grant is to create a multidisciplinary approach to develop spatially explicit models of vector-borne disease risk using Chagas disease as our model. Chagas disease is a parasitic disease endemic to Latin America that afflicts an estimated 10 million people. The causative agent (Trypanosoma cruzi) is most commonly transmitted to humans by blood feeding triatomine insect vectors. Our objectives are: (1) advance knowledge on the multiple interacting factors affecting the transmission of Chagas disease, and (2) provide next generation genomic and spatial analysis tools applicable to the study of other vector-borne diseases worldwide. This funding is a collaborative effort between the RSENR (UVM), the School of Engineering (UVM), the Department of Biology (UVM), the Department of Biological Sciences (Loyola (New Orleans)) and the Laboratory of Applied Entomology and Parasitology (Universidad de San Carlos). Throughout this five-year study, multi-educational groups (i.e., high school, undergraduate, graduate, and postdoctoral) will be trained in systems modeling. This systems approach challenges students to incorporate environmental, social, and economic as well as technical aspects and enables modelers to simulate and visualize topics that would either be too expensive, complex or difficult to study directly (Yasar and Landau 2003). We launch this research by developing a set of multi-scale, epidemiological models of Chagas disease risk using STELLA® software v.9.1.3 (isee systems, inc., Lebanon, NH). We use this particular system dynamics software as a starting point because of its simple graphical user interface (e.g., behavior-over-time graphs, stock/flow diagrams, and causal loops). To date, high school and undergraduate students have created a set of multi-scale (i.e., homestead, village, and regional) disease models. Modeling the system at multiple spatial scales forces recognition that

Multiscale analysis of structure development in expanded starch snacks

Science.gov (United States)

van der Sman, R. G. M.; Broeze, J.

2014-11-01

In this paper we perform a multiscale analysis of the food structuring process of the expansion of starchy snack foods like keropok, which obtains a solid foam structure. In particular, we want to investigate the validity of the hypothesis of Kokini and coworkers, that expansion is optimal at the moisture content, where the glass transition and the boiling line intersect. In our analysis we make use of several tools, (1) time scale analysis from the field of physical transport phenomena, (2) the scale separation map (SSM) developed within a multiscale simulation framework of complex automata, (3) the supplemented state diagram (SSD), depicting phase transition and glass transition lines, and (4) a multiscale simulation model for the bubble expansion. Results of the time scale analysis are plotted in the SSD, and give insight into the dominant physical processes involved in expansion. Furthermore, the results of the time scale analysis are used to construct the SSM, which has aided us in the construction of the multiscale simulation model. Simulation results are plotted in the SSD. This clearly shows that the hypothesis of Kokini is qualitatively true, but has to be refined. Our results show that bubble expansion is optimal for moisture content, where the boiling line for gas pressure of 4 bars intersects the isoviscosity line of the critical viscosity 106 Pa.s, which runs parallel to the glass transition line.

Multiscale integration of -omic, imaging, and clinical data in biomedical informatics.

Science.gov (United States)

Phan, John H; Quo, Chang F; Cheng, Chihwen; Wang, May Dongmei

2012-01-01

This paper reviews challenges and opportunities in multiscale data integration for biomedical informatics. Biomedical data can come from different biological origins, data acquisition technologies, and clinical applications. Integrating such data across multiple scales (e.g., molecular, cellular/tissue, and patient) can lead to more informed decisions for personalized, predictive, and preventive medicine. However, data heterogeneity, community standards in data acquisition, and computational complexity are big challenges for such decision making. This review describes genomic and proteomic (i.e., molecular), histopathological imaging (i.e., cellular/tissue), and clinical (i.e., patient) data; it includes case studies for single-scale (e.g., combining genomic or histopathological image data), multiscale (e.g., combining histopathological image and clinical data), and multiscale and multiplatform (e.g., the Human Protein Atlas and The Cancer Genome Atlas) data integration. Numerous opportunities exist in biomedical informatics research focusing on integration of multiscale and multiplatform data.

Contextual Multi-Scale Region Convolutional 3D Network for Activity Detection

KAUST Repository

Bai, Yancheng

2018-01-28

Activity detection is a fundamental problem in computer vision. Detecting activities of different temporal scales is particularly challenging. In this paper, we propose the contextual multi-scale region convolutional 3D network (CMS-RC3D) for activity detection. To deal with the inherent temporal scale variability of activity instances, the temporal feature pyramid is used to represent activities of different temporal scales. On each level of the temporal feature pyramid, an activity proposal detector and an activity classifier are learned to detect activities of specific temporal scales. Temporal contextual information is fused into activity classifiers for better recognition. More importantly, the entire model at all levels can be trained end-to-end. Our CMS-RC3D detector can deal with activities at all temporal scale ranges with only a single pass through the backbone network. We test our detector on two public activity detection benchmarks, THUMOS14 and ActivityNet. Extensive experiments show that the proposed CMS-RC3D detector outperforms state-of-the-art methods on THUMOS14 by a substantial margin and achieves comparable results on ActivityNet despite using a shallow feature extractor.

Contextual Multi-Scale Region Convolutional 3D Network for Activity Detection

KAUST Repository

Bai, Yancheng; Xu, Huijuan; Saenko, Kate; Ghanem, Bernard

2018-01-01

Activity detection is a fundamental problem in computer vision. Detecting activities of different temporal scales is particularly challenging. In this paper, we propose the contextual multi-scale region convolutional 3D network (CMS-RC3D) for activity detection. To deal with the inherent temporal scale variability of activity instances, the temporal feature pyramid is used to represent activities of different temporal scales. On each level of the temporal feature pyramid, an activity proposal detector and an activity classifier are learned to detect activities of specific temporal scales. Temporal contextual information is fused into activity classifiers for better recognition. More importantly, the entire model at all levels can be trained end-to-end. Our CMS-RC3D detector can deal with activities at all temporal scale ranges with only a single pass through the backbone network. We test our detector on two public activity detection benchmarks, THUMOS14 and ActivityNet. Extensive experiments show that the proposed CMS-RC3D detector outperforms state-of-the-art methods on THUMOS14 by a substantial margin and achieves comparable results on ActivityNet despite using a shallow feature extractor.

Multiscale and Multiphysics Modeling of Additive Manufacturing of Advanced Materials

Science.gov (United States)

Liou, Frank; Newkirk, Joseph; Fan, Zhiqiang; Sparks, Todd; Chen, Xueyang; Fletcher, Kenneth; Zhang, Jingwei; Zhang, Yunlu; Kumar, Kannan Suresh; Karnati, Sreekar

2015-01-01

The objective of this proposed project is to research and develop a prediction tool for advanced additive manufacturing (AAM) processes for advanced materials and develop experimental methods to provide fundamental properties and establish validation data. Aircraft structures and engines demand materials that are stronger, useable at much higher temperatures, provide less acoustic transmission, and enable more aeroelastic tailoring than those currently used. Significant improvements in properties can only be achieved by processing the materials under nonequilibrium conditions, such as AAM processes. AAM processes encompass a class of processes that use a focused heat source to create a melt pool on a substrate. Examples include Electron Beam Freeform Fabrication and Direct Metal Deposition. These types of additive processes enable fabrication of parts directly from CAD drawings. To achieve the desired material properties and geometries of the final structure, assessing the impact of process parameters and predicting optimized conditions with numerical modeling as an effective prediction tool is necessary. The targets for the processing are multiple and at different spatial scales, and the physical phenomena associated occur in multiphysics and multiscale. In this project, the research work has been developed to model AAM processes in a multiscale and multiphysics approach. A macroscale model was developed to investigate the residual stresses and distortion in AAM processes. A sequentially coupled, thermomechanical, finite element model was developed and validated experimentally. The results showed the temperature distribution, residual stress, and deformation within the formed deposits and substrates. A mesoscale model was developed to include heat transfer, phase change with mushy zone, incompressible free surface flow, solute redistribution, and surface tension. Because of excessive computing time needed, a parallel computing approach was also tested. In addition

Multiscale empirical interpolation for solving nonlinear PDEs

KAUST Repository

Calo, Victor M.

2014-12-01

In this paper, we propose a multiscale empirical interpolation method for solving nonlinear multiscale partial differential equations. The proposed method combines empirical interpolation techniques and local multiscale methods, such as the Generalized Multiscale Finite Element Method (GMsFEM). To solve nonlinear equations, the GMsFEM is used to represent the solution on a coarse grid with multiscale basis functions computed offline. Computing the GMsFEM solution involves calculating the system residuals and Jacobians on the fine grid. We use empirical interpolation concepts to evaluate these residuals and Jacobians of the multiscale system with a computational cost which is proportional to the size of the coarse-scale problem rather than the fully-resolved fine scale one. The empirical interpolation method uses basis functions which are built by sampling the nonlinear function we want to approximate a limited number of times. The coefficients needed for this approximation are computed in the offline stage by inverting an inexpensive linear system. The proposed multiscale empirical interpolation techniques: (1) divide computing the nonlinear function into coarse regions; (2) evaluate contributions of nonlinear functions in each coarse region taking advantage of a reduced-order representation of the solution; and (3) introduce multiscale proper-orthogonal-decomposition techniques to find appropriate interpolation vectors. We demonstrate the effectiveness of the proposed methods on several nonlinear multiscale PDEs that are solved with Newton\\'s methods and fully-implicit time marching schemes. Our numerical results show that the proposed methods provide a robust framework for solving nonlinear multiscale PDEs on a coarse grid with bounded error and significant computational cost reduction.

Multiscale Modeling of a Conditionally Disordered pH-Sensing Chaperone

OpenAIRE

Ahlstrom, Logan S.; Law, Sean M.; Dickson, Alex; Brooks, Charles L.

2015-01-01

The pH-sensing chaperone HdeA promotes the survival of enteropathogenic bacteria during transit through the harshly acidic environment of the mammalian stomach. At low pH, HdeA transitions from an inactive, folded, dimer to chaperone-active, disordered, monomers to protect against the acid-induced aggregation of periplasmic proteins. Toward achieving a detailed mechanistic understanding of the pH response of HdeA, we develop a multiscale modeling approach to capture its pH-dependent thermodyn...

A Posteriori Analysis of Adaptive Multiscale Operator Decomposition Methods for Multiphysics Problems

Energy Technology Data Exchange (ETDEWEB)

Donald Estep; Michael Holst; Simon Tavener

2010-02-08

This project was concerned with the accurate computational error estimation for numerical solutions of multiphysics, multiscale systems that couple different physical processes acting across a large range of scales relevant to the interests of the DOE. Multiscale, multiphysics models are characterized by intimate interactions between different physics across a wide range of scales. This poses significant computational challenges addressed by the proposal, including: (1) Accurate and efficient computation; (2) Complex stability; and (3) Linking different physics. The research in this project focused on Multiscale Operator Decomposition methods for solving multiphysics problems. The general approach is to decompose a multiphysics problem into components involving simpler physics over a relatively limited range of scales, and then to seek the solution of the entire system through some sort of iterative procedure involving solutions of the individual components. MOD is a very widely used technique for solving multiphysics, multiscale problems; it is heavily used throughout the DOE computational landscape. This project made a major advance in the analysis of the solution of multiscale, multiphysics problems.

A Novel Multiscale Ensemble Carbon Price Prediction Model Integrating Empirical Mode Decomposition, Genetic Algorithm and Artificial Neural Network

Directory of Open Access Journals (Sweden)

Bangzhu Zhu

2012-02-01

Full Text Available Due to the movement and complexity of the carbon market, traditional monoscale forecasting approaches often fail to capture its nonstationary and nonlinear properties and accurately describe its moving tendencies. In this study, a multiscale ensemble forecasting model integrating empirical mode decomposition (EMD, genetic algorithm (GA and artificial neural network (ANN is proposed to forecast carbon price. Firstly, the proposed model uses EMD to decompose carbon price data into several intrinsic mode functions (IMFs and one residue. Then, the IMFs and residue are composed into a high frequency component, a low frequency component and a trend component which have similar frequency characteristics, simple components and strong regularity using the fine-to-coarse reconstruction algorithm. Finally, those three components are predicted using an ANN trained by GA, i.e., a GAANN model, and the final forecasting results can be obtained by the sum of these three forecasting results. For verification and testing, two main carbon future prices with different maturity in the European Climate Exchange (ECX are used to test the effectiveness of the proposed multiscale ensemble forecasting model. Empirical results obtained demonstrate that the proposed multiscale ensemble forecasting model can outperform the single random walk (RW, ARIMA, ANN and GAANN models without EMD preprocessing and the ensemble ARIMA model with EMD preprocessing.

Multiscale Cues Drive Collective Cell Migration

Science.gov (United States)

Nam, Ki-Hwan; Kim, Peter; Wood, David K.; Kwon, Sunghoon; Provenzano, Paolo P.; Kim, Deok-Ho

2016-07-01

To investigate complex biophysical relationships driving directed cell migration, we developed a biomimetic platform that allows perturbation of microscale geometric constraints with concomitant nanoscale contact guidance architectures. This permits us to elucidate the influence, and parse out the relative contribution, of multiscale features, and define how these physical inputs are jointly processed with oncogenic signaling. We demonstrate that collective cell migration is profoundly enhanced by the addition of contract guidance cues when not otherwise constrained. However, while nanoscale cues promoted migration in all cases, microscale directed migration cues are dominant as the geometric constraint narrows, a behavior that is well explained by stochastic diffusion anisotropy modeling. Further, oncogene activation (i.e. mutant PIK3CA) resulted in profoundly increased migration where extracellular multiscale directed migration cues and intrinsic signaling synergistically conspire to greatly outperform normal cells or any extracellular guidance cues in isolation.

Multiscale modeling of three-dimensional genome

Science.gov (United States)

Zhang, Bin; Wolynes, Peter

The genome, the blueprint of life, contains nearly all the information needed to build and maintain an entire organism. A comprehensive understanding of the genome is of paramount interest to human health and will advance progress in many areas, including life sciences, medicine, and biotechnology. The overarching goal of my research is to understand the structure-dynamics-function relationships of the human genome. In this talk, I will be presenting our efforts in moving towards that goal, with a particular emphasis on studying the three-dimensional organization, the structure of the genome with multi-scale approaches. Specifically, I will discuss the reconstruction of genome structures at both interphase and metaphase by making use of data from chromosome conformation capture experiments. Computationally modeling of chromatin fiber at atomistic level from first principles will also be presented as our effort for studying the genome structure from bottom up.

Industrial process system assessment: bridging process engineering and life cycle assessment through multiscale modeling.

Science.gov (United States)

The Industrial Process System Assessment (IPSA) methodology is a multiple step allocation approach for connecting information from the production line level up to the facility level and vice versa using a multiscale model of process systems. The allocation procedure assigns inpu...

Three-Phase Carbon Fiber Amine Functionalized Carbon Nanotubes Epoxy Composite: Processing, Characterisation, and Multiscale Modeling

Directory of Open Access Journals (Sweden)

Kamal Sharma

2014-01-01

Full Text Available The present paper discusses the key issues of carbon nanotube (CNT dispersion and effect of functionalisation on the mechanical properties of multiscale carbon epoxy composites. In this study, CNTs were added in epoxy matrix and further reinforced with carbon fibres. Predetermined amounts of optimally amine functionalised CNTs were dispersed in epoxy matrix, and unidirectional carbon fiber laminates were produced. The effect of the presence of CNTs (1.0 wt% in the resin was reflected by pronounced increase in Young’s modulus, inter-laminar shear strength, and flexural modulus by 51.46%, 39.62%, and 38.04%, respectively. However, 1.5 wt% CNT loading in epoxy resin decreased the overall properties of the three-phase composites. A combination of Halpin-Tsai equations and micromechanics modeling approach was also used to evaluate the mechanical properties of multiscale composites and the differences between the predicted and experimental values are reported. These multiscale composites are likely to be used for potential missile and aerospace structural applications.

Information Management Workflow and Tools Enabling Multiscale Modeling Within ICME Paradigm

Science.gov (United States)

Arnold, Steven M.; Bednarcyk, Brett A.; Austin, Nic; Terentjev, Igor; Cebon, Dave; Marsden, Will

2016-01-01

With the increased emphasis on reducing the cost and time to market of new materials, the need for analytical tools that enable the virtual design and optimization of materials throughout their processing - internal structure - property - performance envelope, along with the capturing and storing of the associated material and model information across its lifecycle, has become critical. This need is also fueled by the demands for higher efficiency in material testing; consistency, quality and traceability of data; product design; engineering analysis; as well as control of access to proprietary or sensitive information. Fortunately, material information management systems and physics-based multiscale modeling methods have kept pace with the growing user demands. Herein, recent efforts to establish workflow for and demonstrate a unique set of web application tools for linking NASA GRC's Integrated Computational Materials Engineering (ICME) Granta MI database schema and NASA GRC's Integrated multiscale Micromechanics Analysis Code (ImMAC) software toolset are presented. The goal is to enable seamless coupling between both test data and simulation data, which is captured and tracked automatically within Granta MI®, with full model pedigree information. These tools, and this type of linkage, are foundational to realizing the full potential of ICME, in which materials processing, microstructure, properties, and performance are coupled to enable application-driven design and optimization of materials and structures.

A Multiscale Finite Element Model Validation Method of Composite Cable-Stayed Bridge Based on Structural Health Monitoring System

Directory of Open Access Journals (Sweden)

Rumian Zhong

2015-01-01

Full Text Available A two-step response surface method for multiscale finite element model (FEM updating and validation is presented with respect to Guanhe Bridge, a composite cable-stayed bridge in the National Highway number G15, in China. Firstly, the state equations of both multiscale and single-scale FEM are established based on the basic equation in structural dynamic mechanics to update the multiscale coupling parameters and structural parameters. Secondly, based on the measured data from the structural health monitoring (SHM system, a Monte Carlo simulation is employed to analyze the uncertainty quantification and transmission, where the uncertainties of the multiscale FEM and measured data were considered. The results indicate that the relative errors between the calculated and measured frequencies are less than 2%, and the overlap ratio indexes of each modal frequency are larger than 80% without the average absolute value of relative errors. These demonstrate that the proposed method can be applied to validate the multiscale FEM, and the validated FEM can reflect the current conditions of the real bridge; thus it can be used as the basis for bridge health monitoring, damage prognosis (DP, and safety prognosis (SP.

A multi-scale model for correlation in B cell VDJ usage of zebrafish

International Nuclear Information System (INIS)

Pan, Keyao; Deem, Michael W

2011-01-01

The zebrafish (Danio rerio) is one of the model animals used for the study of immunology because the dynamics in the adaptive immune system of zebrafish are similar to that in higher animals. In this work, we built a multi-scale model to simulate the dynamics of B cells in the primary and secondary immune responses of zebrafish. We use this model to explain the reported correlation between VDJ usage of B cell repertoires in individual zebrafish. We use a delay ordinary differential equation (ODE) system to model the immune responses in the 6-month lifespan of a zebrafish. This mean field theory gives the number of high-affinity B cells as a function of time during an infection. The sequences of those B cells are then taken from a distribution calculated by a 'microscopic' random energy model. This generalized NK model shows that mature B cells specific to one antigen largely possess a single VDJ recombination. The model allows first-principle calculation of the probability, p, that two zebrafish responding to the same antigen will select the same VDJ recombination. This probability p increases with the B cell population size and the B cell selection intensity. The probability p decreases with the B cell hypermutation rate. The multi-scale model predicts correlations in the immune system of the zebrafish that are highly similar to that from experiment

A multi-scale model for correlation in B cell VDJ usage of zebrafish

Science.gov (United States)

Pan, Keyao; Deem, Michael W.

2011-10-01

The zebrafish (Danio rerio) is one of the model animals used for the study of immunology because the dynamics in the adaptive immune system of zebrafish are similar to that in higher animals. In this work, we built a multi-scale model to simulate the dynamics of B cells in the primary and secondary immune responses of zebrafish. We use this model to explain the reported correlation between VDJ usage of B cell repertoires in individual zebrafish. We use a delay ordinary differential equation (ODE) system to model the immune responses in the 6-month lifespan of a zebrafish. This mean field theory gives the number of high-affinity B cells as a function of time during an infection. The sequences of those B cells are then taken from a distribution calculated by a 'microscopic' random energy model. This generalized NK model shows that mature B cells specific to one antigen largely possess a single VDJ recombination. The model allows first-principle calculation of the probability, p, that two zebrafish responding to the same antigen will select the same VDJ recombination. This probability p increases with the B cell population size and the B cell selection intensity. The probability p decreases with the B cell hypermutation rate. The multi-scale model predicts correlations in the immune system of the zebrafish that are highly similar to that from experiment.

Microstructure representations for sound absorbing fibrous media: 3D and 2D multiscale modelling and experiments

Science.gov (United States)

Zieliński, Tomasz G.

2017-11-01

The paper proposes and investigates computationally-efficient microstructure representations for sound absorbing fibrous media. Three-dimensional volume elements involving non-trivial periodic arrangements of straight fibres are examined as well as simple two-dimensional cells. It has been found that a simple 2D quasi-representative cell can provide similar predictions as a volume element which is in general much more geometrically accurate for typical fibrous materials. The multiscale modelling allowed to determine the effective speeds and damping of acoustic waves propagating in such media, which brings up a discussion on the correlation between the speed, penetration range and attenuation of sound waves. Original experiments on manufactured copper-wire samples are presented and the microstructure-based calculations of acoustic absorption are compared with the corresponding experimental results. In fact, the comparison suggested the microstructure modifications leading to representations with non-uniformly distributed fibres.

A mathematical framework for multiscale science and engineering: the variational multiscale method and interscale transfer operators

International Nuclear Information System (INIS)

Shadid, John Nicolas; Lehoucq, Richard B.; Christon, Mark Allen; Slepoy, Alexander; Bochev, Pavel Blagoveston; Collis, Samuel Scott; Wagner, Gregory John

2004-01-01

Existing approaches in multiscale science and engineering have evolved from a range of ideas and solutions that are reflective of their original problem domains. As a result, research in multiscale science has followed widely diverse and disjoint paths, which presents a barrier to cross pollination of ideas and application of methods outside their application domains. The status of the research environment calls for an abstract mathematical framework that can provide a common language to formulate and analyze multiscale problems across a range of scientific and engineering disciplines. In such a framework, critical common issues arising in multiscale problems can be identified, explored and characterized in an abstract setting. This type of overarching approach would allow categorization and clarification of existing models and approximations in a landscape of seemingly disjoint, mutually exclusive and ad hoc methods. More importantly, such an approach can provide context for both the development of new techniques and their critical examination. As with any new mathematical framework, it is necessary to demonstrate its viability on problems of practical importance. At Sandia, lab-centric, prototype application problems in fluid mechanics, reacting flows, magnetohydrodynamics (MHD), shock hydrodynamics and materials science span an important subset of DOE Office of Science applications and form an ideal proving ground for new approaches in multiscale science.

Inter-subject FDG PET Brain Networks Exhibit Multi-scale Community Structure with Different Normalization Techniques.

Science.gov (United States)

Sperry, Megan M; Kartha, Sonia; Granquist, Eric J; Winkelstein, Beth A

2018-07-01

Inter-subject networks are used to model correlations between brain regions and are particularly useful for metabolic imaging techniques, like 18F-2-deoxy-2-(18F)fluoro-D-glucose (FDG) positron emission tomography (PET). Since FDG PET typically produces a single image, correlations cannot be calculated over time. Little focus has been placed on the basic properties of inter-subject networks and if they are affected by group size and image normalization. FDG PET images were acquired from rats (n = 18), normalized by whole brain, visual cortex, or cerebellar FDG uptake, and used to construct correlation matrices. Group size effects on network stability were investigated by systematically adding rats and evaluating local network connectivity (node strength and clustering coefficient). Modularity and community structure were also evaluated in the differently normalized networks to assess meso-scale network relationships. Local network properties are stable regardless of normalization region for groups of at least 10. Whole brain-normalized networks are more modular than visual cortex- or cerebellum-normalized network (p network resolutions where modularity differs most between brain and randomized networks. Hierarchical analysis reveals consistent modules at different scales and clustering of spatially-proximate brain regions. Findings suggest inter-subject FDG PET networks are stable for reasonable group sizes and exhibit multi-scale modularity.

A Physiologically Based, Multi-Scale Model of Skeletal Muscle Structure and Function

Science.gov (United States)

Röhrle, O.; Davidson, J. B.; Pullan, A. J.

2012-01-01

Models of skeletal muscle can be classified as phenomenological or biophysical. Phenomenological models predict the muscle’s response to a specified input based on experimental measurements. Prominent phenomenological models are the Hill-type muscle models, which have been incorporated into rigid-body modeling frameworks, and three-dimensional continuum-mechanical models. Biophysically based models attempt to predict the muscle’s response as emerging from the underlying physiology of the system. In this contribution, the conventional biophysically based modeling methodology is extended to include several structural and functional characteristics of skeletal muscle. The result is a physiologically based, multi-scale skeletal muscle finite element model that is capable of representing detailed, geometrical descriptions of skeletal muscle fibers and their grouping. Together with a well-established model of motor-unit recruitment, the electro-physiological behavior of single muscle fibers within motor units is computed and linked to a continuum-mechanical constitutive law. The bridging between the cellular level and the organ level has been achieved via a multi-scale constitutive law and homogenization. The effect of homogenization has been investigated by varying the number of embedded skeletal muscle fibers and/or motor units and computing the resulting exerted muscle forces while applying the same excitatory input. All simulations were conducted using an anatomically realistic finite element model of the tibialis anterior muscle. Given the fact that the underlying electro-physiological cellular muscle model is capable of modeling metabolic fatigue effects such as potassium accumulation in the T-tubular space and inorganic phosphate build-up, the proposed framework provides a novel simulation-based way to investigate muscle behavior ranging from motor-unit recruitment to force generation and fatigue. PMID:22993509

A physiologically based, multi-scale model of skeletal muscle structure and function

Directory of Open Access Journals (Sweden)

Oliver eRöhrle

2012-09-01

Full Text Available Models of skeletal muscle can be classified as phenomenological or biophysical. Phenomenological models predict the muscle's response to a specified input based on experimental measurements. Prominent phenomenological models are the Hill-type muscle models, which have been incorporated into rigid-body modelling frameworks, and three-dimensional continuum-mechanical models. Biophysically based models attempt to predict the muscle's response as emerging from the underlying physiology of the system. In this contribution, the conventional biophysically based modelling methodology is extended to include several structural and functional characteristics of skeletal muscle. The result is a physiologically based, multi-scale skeletal muscle finite element model that is capable of representing detailed, geometrical descriptions of skeletal muscle fibres and their grouping. Together with a well-established model of motor unit recruitment, the electro-physiological behaviour of single muscle fibres within motor units is computed and linked to a continuum-mechanical constitutive law. The bridging between the cellular level and the organ level has been achieved via a multi-scale constitutive law and homogenisation. The effect of homogenisation has been investigated by varying the number of embedded skeletal muscle fibres and/or motor units and computing the resulting exerted muscle forces while applying the same excitatory input. All simulations were conducted using an anatomically realistic finite element model of the Tibialis Anterior muscle. Given the fact that the underlying electro-physiological cellular muscle model is capable of modelling metabolic fatigue effects such as potassium accumulation in the T-tubular space and inorganic phosphate build-up, the proposed framework provides a novel simulation-based way to investigate muscle behaviour ranging from motor unit recruitment to force generation and fatigue.

Multiscale landscape genomic models to detect signatures of selection in the alpine plant Biscutella laevigata.

Science.gov (United States)

Leempoel, Kevin; Parisod, Christian; Geiser, Céline; Joost, Stéphane

2018-02-01

Plant species are known to adapt locally to their environment, particularly in mountainous areas where conditions can vary drastically over short distances. The climate of such landscapes being largely influenced by topography, using fine-scale models to evaluate environmental heterogeneity may help detecting adaptation to micro-habitats. Here, we applied a multiscale landscape genomic approach to detect evidence of local adaptation in the alpine plant Biscutella laevigata . The two gene pools identified, experiencing limited gene flow along a 1-km ridge, were different in regard to several habitat features derived from a very high resolution (VHR) digital elevation model (DEM). A correlative approach detected signatures of selection along environmental gradients such as altitude, wind exposure, and solar radiation, indicating adaptive pressures likely driven by fine-scale topography. Using a large panel of DEM-derived variables as ecologically relevant proxies, our results highlighted the critical role of spatial resolution. These high-resolution multiscale variables indeed indicate that the robustness of associations between genetic loci and environmental features depends on spatial parameters that are poorly documented. We argue that the scale issue is critical in landscape genomics and that multiscale ecological variables are key to improve our understanding of local adaptation in highly heterogeneous landscapes.

Multiscale digital Arabidopsis predicts individual organ and whole-organism growth.

Science.gov (United States)

Chew, Yin Hoon; Wenden, Bénédicte; Flis, Anna; Mengin, Virginie; Taylor, Jasper; Davey, Christopher L; Tindal, Christopher; Thomas, Howard; Ougham, Helen J; de Reffye, Philippe; Stitt, Mark; Williams, Mathew; Muetzelfeldt, Robert; Halliday, Karen J; Millar, Andrew J

2014-09-30

Understanding how dynamic molecular networks affect whole-organism physiology, analogous to mapping genotype to phenotype, remains a key challenge in biology. Quantitative models that represent processes at multiple scales and link understanding from several research domains can help to tackle this problem. Such integrated models are more common in crop science and ecophysiology than in the research communities that elucidate molecular networks. Several laboratories have modeled particular aspects of growth in Arabidopsis thaliana, but it was unclear whether these existing models could productively be combined. We test this approach by constructing a multiscale model of Arabidopsis rosette growth. Four existing models were integrated with minimal parameter modification (leaf water content and one flowering parameter used measured data). The resulting framework model links genetic regulation and biochemical dynamics to events at the organ and whole-plant levels, helping to understand the combined effects of endogenous and environmental regulators on Arabidopsis growth. The framework model was validated and tested with metabolic, physiological, and biomass data from two laboratories, for five photoperiods, three accessions, and a transgenic line, highlighting the plasticity of plant growth strategies. The model was extended to include stochastic development. Model simulations gave insight into the developmental control of leaf production and provided a quantitative explanation for the pleiotropic developmental phenotype caused by overexpression of miR156, which was an open question. Modular, multiscale models, assembling knowledge from systems biology to ecophysiology, will help to understand and to engineer plant behavior from the genome to the field.

Multi-Scale Parameter Identification of Lithium-Ion Battery Electric Models Using a PSO-LM Algorithm

Directory of Open Access Journals (Sweden)

Wen-Jing Shen

2017-03-01

Full Text Available This paper proposes a multi-scale parameter identification algorithm for the lithium-ion battery (LIB electric model by using a combination of particle swarm optimization (PSO and Levenberg-Marquardt (LM algorithms. Two-dimensional Poisson equations with unknown parameters are used to describe the potential and current density distribution (PDD of the positive and negative electrodes in the LIB electric model. The model parameters are difficult to determine in the simulation due to the nonlinear complexity of the model. In the proposed identification algorithm, PSO is used for the coarse-scale parameter identification and the LM algorithm is applied for the fine-scale parameter identification. The experiment results show that the multi-scale identification not only improves the convergence rate and effectively escapes from the stagnation of PSO, but also overcomes the local minimum entrapment drawback of the LM algorithm. The terminal voltage curves from the PDD model with the identified parameter values are in good agreement with those from the experiments at different discharge/charge rates.

Residual-driven online generalized multiscale finite element methods

KAUST Repository

Chung, Eric T.

2015-09-08

The construction of local reduced-order models via multiscale basis functions has been an area of active research. In this paper, we propose online multiscale basis functions which are constructed using the offline space and the current residual. Online multiscale basis functions are constructed adaptively in some selected regions based on our error indicators. We derive an error estimator which shows that one needs to have an offline space with certain properties to guarantee that additional online multiscale basis function will decrease the error. This error decrease is independent of physical parameters, such as the contrast and multiple scales in the problem. The offline spaces are constructed using Generalized Multiscale Finite Element Methods (GMsFEM). We show that if one chooses a sufficient number of offline basis functions, one can guarantee that additional online multiscale basis functions will reduce the error independent of contrast. We note that the construction of online basis functions is motivated by the fact that the offline space construction does not take into account distant effects. Using the residual information, we can incorporate the distant information provided the offline approximation satisfies certain properties. In the paper, theoretical and numerical results are presented. Our numerical results show that if the offline space is sufficiently large (in terms of the dimension) such that the coarse space contains all multiscale spectral basis functions that correspond to small eigenvalues, then the error reduction by adding online multiscale basis function is independent of the contrast. We discuss various ways computing online multiscale basis functions which include a use of small dimensional offline spaces.

High resolution multi-scale air quality modelling for all streets in Denmark

DEFF Research Database (Denmark)

Jensen, Steen Solvang; Ketzel, Matthias; Becker, Thomas

2017-01-01

The annual concentrations of NO2, PM2.5 and PM10 in 2012 have for the first time been modelled for all 2.4 million addresses in Denmark based on a multi-scale air quality modelling approach. All addresses include residential, industrial, institutional, shop, school, restaurant addresses etc...... concentrations of NO2 for the five available street monitoring stations are within −27% to +12%. The model results were also verified with comparisons with previous model results for NO2 at 98 selected streets in Copenhagen and 31 streets in Aalborg. The verification showed good correlation in Copenhagen (r2 = 0...

Development of hybrid 3-D hydrological modeling for the NCAR Community Earth System Model (CESM)

Energy Technology Data Exchange (ETDEWEB)

Zeng, Xubin [Univ. of Arizona, Tucson, AZ (United States); Troch, Peter [Univ. of Arizona, Tucson, AZ (United States); Pelletier, Jon [Univ. of Arizona, Tucson, AZ (United States); Niu, Guo-Yue [Univ. of Arizona, Tucson, AZ (United States); Gochis, David [NCAR Research Applications Lab., Boulder, CO (United States)

2015-11-15

This is the Final Report of our four-year (3-year plus one-year no cost extension) collaborative project between the University of Arizona (UA) and the National Center for Atmospheric Research (NCAR). The overall objective of our project is to develop and evaluate the first hybrid 3-D hydrological model with a horizontal grid spacing of 1 km for the NCAR Community Earth System Model (CESM).

Multiscale modeling of dislocation processes in BCC tantalum: bridging atomistic and mesoscale simulations

International Nuclear Information System (INIS)

Yang, L H; Tang, M; Moriarty, J A

2001-01-01

Plastic deformation in bcc metals at low temperatures and high-strain rates is controlled by the motion of a/2 screw dislocations, and understanding the fundamental atomistic processes of this motion is essential to develop predictive multiscale models of crystal plasticity. The multiscale modeling approach presented here for bcc Ta is based on information passing, where results of simulations at the atomic scale are used in simulations of plastic deformation at mesoscopic length scales via dislocation dynamics (DD). The relevant core properties of a/2 screw dislocations in Ta have been obtained using quantum-based interatomic potentials derived from model generalized pseudopotential theory and an ab-initio data base together with an accurate Green's-function simulation method that implements flexible boundary conditions. In particular, the stress-dependent activation enthalpy for the lowest-energy kink-pair mechanism has been calculated and fitted to a revealing analytic form. This is the critical quantity determining dislocation mobility in the DD simulations, and the present activation enthalpy is found to be in good agreement with the previous empirical form used to explain the temperature dependence of the yield stress

Shaken but not stirred: Multiscale habitat suitability modeling of sympatric marten species (Martes martes and Martes foina) in the northern Iberian Peninsula

Science.gov (United States)

Maria Vergara; Samuel A. Cushman; Fermin Urra; Aritz Ruiz-Gonzalez

2016-01-01

Multispecies and multiscale habitat suitability models (HSM) are important to identify the environmental variables and scales influencing habitat selection and facilitate the comparison of closely related species with different ecological requirements. Objectives This study explores the multiscale relationships of habitat suitability for the pine (Martes...

Multiscale forward electromagnetic model of uterine contractions during pregnancy

International Nuclear Information System (INIS)

La Rosa, Patricio S; Eswaran, Hari; Preissl, Hubert; Nehorai, Arye

2012-01-01

Analyzing and monitoring uterine contractions during pregnancy is relevant to the field of reproductive health assessment. Its clinical importance is grounded in the need to reliably predict the onset of labor at term and pre-term. Preterm births can cause health problems or even be fatal for the fetus. Currently, there are no objective methods for consistently predicting the onset of labor based on sensing of the mechanical or electrophysiological aspects of uterine contractions. Therefore, modeling uterine contractions could help to better interpret such measurements and to develop more accurate methods for predicting labor. In this work, we develop a multiscale forward electromagnetic model of myometrial contractions during pregnancy. In particular, we introduce a model of myometrial current source densities and compute its magnetic field and action potential at the abdominal surface, using Maxwell’s equations and a four-compartment volume conductor geometry. To model the current source density at the myometrium we use a bidomain approach. We consider a modified version of the Fitzhugh-Nagumo (FHN) equation for modeling ionic currents in each myocyte, assuming a plateau-type transmembrane potential, and we incorporate the anisotropic nature of the uterus by designing conductivity-tensor fields. We illustrate our modeling approach considering a spherical uterus and one pacemaker located in the fundus. We obtained a travelling transmembrane potential depolarizing from −56 mV to −16 mV and an average potential in the plateau area of −25 mV with a duration, before hyperpolarization, of 35 s, which is a good approximation with respect to the average recorded transmembrane potentials at term reported in the technical literature. Similarly, the percentage of myometrial cells contracting as a function of time had the same symmetric properties and duration as the intrauterine pressure waveforms of a pregnant human myometrium at term. We introduced a multiscale

Multiscale Phase Inversion of Seismic Data

KAUST Repository

Fu, Lei

2017-12-02

We present a scheme for multiscale phase inversion (MPI) of seismic data that is less sensitive to the unmodeled physics of wave propagation and a poor starting model than standard full waveform inversion (FWI). To avoid cycle-skipping, the multiscale strategy temporally integrates the traces several times, i.e. high-order integration, to produce low-boost seismograms that are used as input data for the initial iterations of MPI. As the iterations proceed, higher frequencies in the data are boosted by using integrated traces of lower order as the input data. The input data are also filtered into different narrow frequency bands for the MPI implementation. At low frequencies, we show that MPI with windowed reflections approximates wave equation inversion of the reflection traveltimes, except no traveltime picking is needed. Numerical results with synthetic acoustic data show that MPI is more robust than conventional multiscale FWI when the initial model is far from the true model. Results from synthetic viscoacoustic and elastic data show that MPI is less sensitive than FWI to some of the unmodeled physics. Inversion of marine data shows that MPI is more robust and produces modestly more accurate results than FWI for this data set.

Multi-scale modeling of urban air pollution: development and application of a Street-in-Grid model (v1.0) by coupling MUNICH (v1.0) and Polair3D (v1.8.1)

OpenAIRE

Y. Kim; Y. Wu; C. Seigneur; Y. Roustan

2018-01-01

A new multi-scale model of urban air pollution is presented. This model combines a chemistry–transport model (CTM) that includes a comprehensive treatment of atmospheric chemistry and transport on spatial scales down to 1 km and a street-network model that describes the atmospheric concentrations of pollutants in an urban street network. The street-network model is the Model of Urban Network of Intersecting Canyons and Highways (MUNICH), which consists of two main components...

Multi-scale damage modelling in a ceramic matrix composite using a finite-element microstructure meshfree methodology

Science.gov (United States)

2016-01-01

The problem of multi-scale modelling of damage development in a SiC ceramic fibre-reinforced SiC matrix ceramic composite tube is addressed, with the objective of demonstrating the ability of the finite-element microstructure meshfree (FEMME) model to introduce important aspects of the microstructure into a larger scale model of the component. These are particularly the location, orientation and geometry of significant porosity and the load-carrying capability and quasi-brittle failure behaviour of the fibre tows. The FEMME model uses finite-element and cellular automata layers, connected by a meshfree layer, to efficiently couple the damage in the microstructure with the strain field at the component level. Comparison is made with experimental observations of damage development in an axially loaded composite tube, studied by X-ray computed tomography and digital volume correlation. Recommendations are made for further development of the model to achieve greater fidelity to the microstructure. This article is part of the themed issue ‘Multiscale modelling of the structural integrity of composite materials’. PMID:27242308

Modelling future impacts of air pollution using the multi-scale UK Integrated Assessment Model (UKIAM).

Science.gov (United States)

Oxley, Tim; Dore, Anthony J; ApSimon, Helen; Hall, Jane; Kryza, Maciej

2013-11-01

Integrated assessment modelling has evolved to support policy development in relation to air pollutants and greenhouse gases by providing integrated simulation tools able to produce quick and realistic representations of emission scenarios and their environmental impacts without the need to re-run complex atmospheric dispersion models. The UK Integrated Assessment Model (UKIAM) has been developed to investigate strategies for reducing UK emissions by bringing together information on projected UK emissions of SO2, NOx, NH3, PM10 and PM2.5, atmospheric dispersion, criteria for protection of ecosystems, urban air quality and human health, and data on potential abatement measures to reduce emissions, which may subsequently be linked to associated analyses of costs and benefits. We describe the multi-scale model structure ranging from continental to roadside, UK emission sources, atmospheric dispersion of emissions, implementation of abatement measures, integration with European-scale modelling, and environmental impacts. The model generates outputs from a national perspective which are used to evaluate alternative strategies in relation to emissions, deposition patterns, air quality metrics and ecosystem critical load exceedance. We present a selection of scenarios in relation to the 2020 Business-As-Usual projections and identify potential further reductions beyond those currently being planned. © 2013.

Final Technical Report "Multiscale Simulation Algorithms for Biochemical Systems"

Energy Technology Data Exchange (ETDEWEB)

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.

Process analysis and sensitivity study of regional ozone formation over the Pearl River Delta, China, during the PRIDE-PRD2004 campaign using the Community Multiscale Air Quality modeling system

Directory of Open Access Journals (Sweden)

X. Wang

2010-05-01

Full Text Available In this study, the Community Multiscale Air Quality (CMAQ modeling system is used to simulate the ozone (O₃ episodes during the Program of Regional Integrated Experiments of Air Quality over the Pearl River Delta, China, in October 2004 (PRIDE-PRD2004. The simulation suggests that O₃ pollution is a regional phenomenon in the Pearl River Delta (PRD. Elevated O₃ levels often occurred in the southwestern inland PRD, Pearl River estuary (PRE, and southern coastal areas during the 1-month field campaign. Three evolution patterns of simulated surface O₃ are summarized based on different near-ground flow conditions. More than 75% of days featured interactions between weak synoptic forcing and local sea-land circulation. Integrated process rate (IPR analysis shows that photochemical production is a dominant contributor to O₃ enhancement from 09:00 to 15:00 local standard time in the atmospheric boundary layer over most areas with elevated O₃ occurrence in the mid-afternoon. The simulated ozone production efficiency is 2–8 O₃ molecules per NO_x molecule oxidized in areas with high O₃ chemical production. Precursors of O₃ originating from different source regions in the central PRD are mixed during the course of transport to downwind rural areas during nighttime and early morning, where they then contribute to the daytime O₃ photochemical production. The sea-land circulation plays an important role on the regional O₃ formation and distribution over PRD. Sensitivity studies suggest that O₃ formation is volatile-organic-compound-limited in the central inland PRD, PRE, and surrounding coastal areas with less chemical aging (NO_x/NO_y>0.6, but is NO_x-limited in the rural southwestern PRD with aged air (NO_x/NO_y<0.3.

The Lagrangian particle dispersion model FLEXPART-WRF VERSION 3.1

Energy Technology Data Exchange (ETDEWEB)

Brioude, J.; Arnold, D.; Stohl, A.; Cassiani, M.; Morton, Don; Seibert, P.; Angevine, W. M.; Evan, S.; Dingwell, A.; Fast, Jerome D.; Easter, Richard C.; Pisso, I.; Bukhart, J.; Wotawa, G.

2013-11-01

The Lagrangian particle dispersion model FLEXPART was originally designed for cal- culating long-range and mesoscale dispersion of air pollutants from point sources, such as after an accident in a nuclear power plant. In the meantime FLEXPART has evolved into a comprehensive tool for atmospheric transport modeling and analysis at different scales. This multiscale need from the modeler community has encouraged new developments in FLEXPART. In this document, we present a version that works with the Weather Research and Forecasting (WRF) mesoscale meteoro- logical model. Simple procedures on how to run FLEXPART-WRF are presented along with special options and features that differ from its predecessor versions. In addition, test case data, the source code and visualization tools are provided to the reader as supplementary material.

On a multiscale approach for filter efficiency simulations

KAUST Repository

Iliev, Oleg

2014-07-01

Filtration in general, and the dead end depth filtration of solid particles out of fluid in particular, is intrinsic multiscale problem. The deposition (capturing of particles) essentially depends on local velocity, on microgeometry (pore scale geometry) of the filtering medium and on the diameter distribution of the particles. The deposited (captured) particles change the microstructure of the porous media what leads to change of permeability. The changed permeability directly influences the velocity field and pressure distribution inside the filter element. To close the loop, we mention that the velocity influences the transport and deposition of particles. In certain cases one can evaluate the filtration efficiency considering only microscale or only macroscale models, but in general an accurate prediction of the filtration efficiency requires multiscale models and algorithms. This paper discusses the single scale and the multiscale models, and presents a fractional time step discretization algorithm for the multiscale problem. The velocity within the filter element is computed at macroscale, and is used as input for the solution of microscale problems at selected locations of the porous medium. The microscale problem is solved with respect to transport and capturing of individual particles, and its solution is postprocessed to provide permeability values for macroscale computations. Results from computational experiments with an oil filter are presented and discussed.

Multi-scale Mexican spotted owl (Strix occidentalis lucida) nest/roost habitat selection in Arizona and a comparison with single-scale modeling results

Science.gov (United States)

Brad C. Timm; Kevin McGarigal; Samuel A. Cushman; Joseph L. Ganey

2016-01-01

Efficacy of future habitat selection studies will benefit by taking a multi-scale approach. In addition to potentially providing increased explanatory power and predictive capacity, multi-scale habitat models enhance our understanding of the scales at which species respond to their environment, which is critical knowledge required to implement effective...

A variational multiscale constitutive model for nanocrystalline materials

KAUST Repository

Gurses, Ercan

2011-03-01

This paper presents a variational multi-scale constitutive model in the finite deformation regime capable of capturing the mechanical behavior of nanocrystalline (nc) fcc metals. The nc-material is modeled as a two-phase material consisting of a grain interior phase and a grain boundary effected zone (GBAZ). A rate-independent isotropic porous plasticity model is employed to describe the GBAZ, whereas a crystal-plasticity model which accounts for the transition from partial dislocation to full dislocation mediated plasticity is employed for the grain interior. The constitutive models of both phases are formulated in a small strain framework and extended to finite deformation by use of logarithmic and exponential mappings. Assuming the rule of mixtures, the overall behavior of a given grain is obtained via volume averaging. The scale transition from a single grain to a polycrystal is achieved by Taylor-type homogenization where a log-normal grain size distribution is assumed. It is shown that the proposed model is able to capture the inverse HallPetch effect, i.e., loss of strength with grain size refinement. Finally, the predictive capability of the model is validated against experimental results on nanocrystalline copper and nickel. © 2010 Elsevier Ltd. All rights reserved.

Modeling Urban Collaborative Growth Dynamics Using a Multiscale Simulation Model for the Wuhan Urban Agglomeration Area, China

Directory of Open Access Journals (Sweden)

Yan Yu

2018-05-01

Full Text Available Urban agglomeration has become the predominant form of urbanization in China. In this process, spatial interaction evidently played a significant role in promoting the collaborative development of these correlated cities. The traditional urban model’s focus on individual cities should be transformed to an urban system model. In this study, a multi-scale simulation model has been proposed to simulate the agglomeration development process of the Wuhan urban agglomeration area by embedding the multi-scale spatial interaction into the transition rule system of cellular automata (CA. A system dynamic model was used to predict the demand for new urban land at an aggregated urban agglomeration area scale. A data field approach was adopted to measuring the interaction of intercity at city scale. Neighborhood interaction was interpreted with a logistic regression method at the land parcel scale. Land use data from 1995, 2005, and 2015 were used to calibrate and evaluate the model. The simulation results show that there has been continuing urban growth in the Wuhan urban agglomeration area from 1995 to 2020. Although extension-sprawl was the predominant pattern of urban spatial expansion, the trend of extensive growth to intensive growth is clear during the entire period. The spatial interaction among these cities has been reinforced, which guided the collaborative development and formed the regional urban system network.

Multi-scale modeling of inter-granular fracture in UO2

Energy Technology Data Exchange (ETDEWEB)

Chakraborty, Pritam [Idaho National Lab. (INL), Idaho Falls, ID (United States); Zhang, Yongfeng [Idaho National Lab. (INL), Idaho Falls, ID (United States); Tonks, Michael R. [Idaho National Lab. (INL), Idaho Falls, ID (United States); Biner, S. Bulent [Idaho National Lab. (INL), Idaho Falls, ID (United States)

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.

Multiscale modeling of radiation damage in Fe-based alloys in the fusion environment

International Nuclear Information System (INIS)

Wirth, B.D.; Odette, G.R.; Marian, J.; Ventelon, L.; Young-Vandersall, J.A.; Zepeda-Ruiz, L.A.

2004-01-01

Ferritic alloys represent a technologically important class of candidate materials for fusion first wall and blanket structures. A detailed understanding of the mechanisms of defect accumulation and microstructure evolution, and the corresponding effects on mechanical properties is required to predict their in-service structural performance limits. The physical processes involved in radiation damage, and its effects on mechanical properties, are inherently multiscale and hierarchical, spanning length and time scales from the atomic nucleus to meters and picosecond to decades. In this paper, we present a multiscale modeling methodology to describe radiation effects within the fusion energy environment. Selected results from atomic scale investigation are presented, focusing on (i) the mechanisms of self-interstitial dislocation loop formation with Burgers vector of a in iron relative to vanadium, (ii) helium transport and (iii) the interaction between helium and small self-interstitial clusters in iron, and (iv) dislocation-helium bubble interactions in fcc aluminum

Generalized Multiscale Finite-Element Method (GMsFEM) for elastic wave propagation in heterogeneous, anisotropic media

KAUST Repository

Gao, Kai

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 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 boundaries and the 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.

The Community WRF-Hydro Modeling System Version 4 Updates: Merging Toward Capabilities of the National Water Model

Science.gov (United States)

McAllister, M.; Gochis, D.; Dugger, A. L.; Karsten, L. R.; McCreight, J. L.; Pan, L.; Rafieeinasab, A.; Read, L. K.; Sampson, K. M.; Yu, W.

2017-12-01

The community WRF-Hydro modeling system is publicly available and provides researchers and operational forecasters a flexible and extensible capability for performing multi-scale, multi-physics options for hydrologic modeling that can be run independent or fully-interactive with the WRF atmospheric model. The core WRF-Hydro physics model contains very high-resolution descriptions of terrestrial hydrologic process representations such as land-atmosphere exchanges of energy and moisture, snowpack evolution, infiltration, terrain routing, channel routing, basic reservoir representation and hydrologic data assimilation. Complementing the core physics components of WRF-Hydro are an ecosystem of pre- and post-processing tools that facilitate the preparation of terrain and meteorological input data, an open-source hydrologic model evaluation toolset (Rwrfhydro), hydrologic data assimilation capabilities with DART and advanced model visualization capabilities. The National Center for Atmospheric Research (NCAR), through collaborative support from the National Science Foundation and other funding partners, provides community support for the entire WRF-Hydro system through a variety of mechanisms. This presentation summarizes the enhanced user support capabilities that are being developed for the community WRF-Hydro modeling system. These products and services include a new website, open-source code repositories, documentation and user guides, test cases, online training materials, live, hands-on training sessions, an email list serve, and individual user support via email through a new help desk ticketing system. The WRF-Hydro modeling system and supporting tools which now include re-gridding scripts and model calibration have recently been updated to Version 4 and are merging toward capabilities of the National Water Model.

Multiscale patterns and drivers of arbuscular mycorrhizal fungal communities in the roots and root-associated soil of a wild perennial herb.

Science.gov (United States)

Rasmussen, Pil U; Hugerth, Luisa W; Blanchet, F Guillaume; Andersson, Anders F; Lindahl, Björn D; Tack, Ayco J M

2018-03-24

Arbuscular mycorrhizal (AM) fungi form diverse communities and are known to influence above-ground community dynamics and biodiversity. However, the multiscale patterns and drivers of AM fungal composition and diversity are still poorly understood. We sequenced DNA markers from roots and root-associated soil from Plantago lanceolata plants collected across multiple spatial scales to allow comparison of AM fungal communities among neighbouring plants, plant subpopulations, nearby plant populations, and regions. We also measured soil nutrients, temperature, humidity, and community composition of neighbouring plants and nonAM root-associated fungi. AM fungal communities were already highly dissimilar among neighbouring plants (c. 30 cm apart), albeit with a high variation in the degree of similarity at this small spatial scale. AM fungal communities were increasingly, and more consistently, dissimilar at larger spatial scales. Spatial structure and environmental drivers explained a similar percentage of the variation, from 7% to 25%. A large fraction of the variation remained unexplained, which may be a result of unmeasured environmental variables, species interactions and stochastic processes. We conclude that AM fungal communities are highly variable among nearby plants. AM fungi may therefore play a major role in maintaining small-scale variation in community dynamics and biodiversity. © 2018 The Authors New Phytologist © 2018 New Phytologist Trust.

Mapping rice ecosystem dynamics and greenhouse gas emissions using multiscale imagery and biogeochemical models

Science.gov (United States)

Salas, W.; Torbick, N.

2017-12-01

Rice greenhouse gas (GHG) emissions in production hot spots have been mapped using multiscale satellite imagery and a processed-based biogeochemical model. The multiscale Synthetic Aperture Radar (SAR) and optical imagery were co-processed and fed into a machine leanring framework to map paddy attributes that are tuned using field observations and surveys. Geospatial maps of rice extent, crop calendar, hydroperiod, and cropping intensity were then used to parameterize the DeNitrification-DeComposition (DNDC) model to estimate emissions. Results, in the Red River Detla for example, show total methane emissions at 345.4 million kgCH4-C equivalent to 11.5 million tonnes CO2e (carbon dioxide equivalent). We further assessed the role of Alternative Wetting and Drying and the impact on GHG and yield across production hot spots with uncertainty estimates. The approach described in this research provides a framework for using SAR to derive maps of rice and landscape characteristics to drive process models like DNDC. These types of tools and approaches will support the next generation of Monitoring, Reporting, and Verification (MRV) to combat climate change and support ecosystem service markets.

Technical Work Plan for: Additional Multiscale Thermohydrologic Modeling

International Nuclear Information System (INIS)

B. Kirstein

2006-01-01

will be evaluated and justified. Some of this evaluation will be conducted in conjunction with the post-model development validation activity involving comparisons of predicted TH conditions with measured TH conditions in the DST. The expected result is that, consistent with what was found in Revision 03 of Multiscale Thermohydrologic Model (BSC 2005 [DIRS 173944], Section 6.3.9), near-field/in-drift TH behavior is insensitive to a wide range of host-rock hydrologic property values. It is the intention of the work described in this TWP to propagate the new infiltration fluxes from the replacement infiltration model, by using the percolation fluxes from the revised site-scale unsaturated zone (UZ) flow model that has applied those new infiltration fluxes. The percolation flux distributions will be obtained from the updated site-scale UZ flow model, which has applied updated infiltration flux maps. Another objective of the work scope is to develop, implement, and validate a revised TH submodel-construction approach. This revised approach utilizes interpolation among a set of generic LDTH submodels that are run for a range of percolation flux histories that cover a sufficiently broad range of infiltration flux uncertainty, as well as for four host-rock units (two lithophysal units and two nonlithophysal units), and for three thermal property sets (low, mean, and high). A key motivation for this revised LDTH submodel-construction approach is to enable the MSTHM to be more flexible in addressing a broad range of infiltration flux cases. This approach allows the generic LDTH submodel simulations to be conducted prior to receiving percolation flux maps

Study on high density multi-scale calculation technique

International Nuclear Information System (INIS)

Sekiguchi, S.; Tanaka, Y.; Nakada, H.; Nishikawa, T.; Yamamoto, N.; Yokokawa, M.

2004-01-01

To understand degradation of nuclear materials under irradiation, it is essential to know as much about each phenomenon observed from multi-scale points of view; they are micro-scale in atomic-level, macro-level in structural scale and intermediate level. In this study for application to meso-scale materials (100A ∼ 2μm), computer technology approaching from micro- and macro-scales was developed including modeling and computer application using computational science and technology method. And environmental condition of grid technology for multi-scale calculation was prepared. The software and MD (molecular dynamics) stencil for verifying the multi-scale calculation were improved and their movement was confirmed. (A. Hishinuma)

Multi-scale modeling of spin transport in organic semiconductors

Science.gov (United States)

Hemmatiyan, Shayan; Souza, Amaury; Kordt, Pascal; McNellis, Erik; Andrienko, Denis; Sinova, Jairo

In this work, we present our theoretical framework to simulate simultaneously spin and charge transport in amorphous organic semiconductors. By combining several techniques e.g. molecular dynamics, density functional theory and kinetic Monte Carlo, we are be able to study spin transport in the presence of anisotropy, thermal effects, magnetic and electric field effects in a realistic morphologies of amorphous organic systems. We apply our multi-scale approach to investigate the spin transport in amorphous Alq3 (Tris(8-hydroxyquinolinato)aluminum) and address the underlying spin relaxation mechanism in this system as a function of temperature, bias voltage, magnetic field and sample thickness.

A Multiscale Closed-Loop Cardiovascular Model, with Applications to Heart Pacing and Hemorrhage

Science.gov (United States)

Canuto, Daniel; Eldredge, Jeff; Chong, Kwitae; Benharash, Peyman; Dutson, Erik

2017-11-01

A computational tool is developed for simulating the dynamic response of the human cardiovascular system to various stressors and injuries. The tool couples zero-dimensional models of the heart, pulmonary vasculature, and peripheral vasculature to one-dimensional models of the major systemic arteries. To simulate autonomic response, this multiscale circulatory model is integrated with a feedback model of the baroreflex, allowing control of heart rate, cardiac contractility, and peripheral impedance. The performance of the tool is demonstrated in two scenarios: increasing heart rate by stimulating the sympathetic nervous system, and an acute 10 percent hemorrhage from the left femoral artery.

Multiscale Modeling of Wear Degradation in Cylinder Liners

KAUST Repository

Moraes, Alvaro; Ruggeri, Fabrizio; Tempone, Raul; Vilanova, Pedro

2014-01-01

both to predict and to avoid them. To achieve this, a monitoring system of the wear level should be implemented to decrease the risk of failure. In this work, we take a first step into the development of a multiscale indirect inference methodology

ESCOMPTE 2001: multi-scale modelling and experimental validation

Science.gov (United States)

Cousin, F.; Tulet, P.; Rosset, R.

2003-04-01

ESCOMPTE is a European pollution field experiment located in the Marseille / Fos-Berre area in the summer 2001.This Mediterranean area, with frequent pollution peaks, is characterized by a complex topography subject to sea breeze regimes, together with intense localized urban, industrial and biogenic sources. Four POI have been selected, the most significant being POI2a / b, a 6-day pollution episode extensively documented for dynamics, radiation, gas phase and aerosols, with surface measurements (including measurements at sea in the gulf of Genoa, on board instrumented ferries between Marseille and Corsica), 7 aircrafts, lidar, radar and constant-level flight balloon soundings. The two-way mesoscale model MESO-NH-C (MNH-C), with horizontal resolutions of 9 and 3 km and high vertical resolution (up to 40 levels in the first 2 km), embedded in the global CTM Mocage, has been run for all POIs, with a focus here on POI2b (June 24-27,2001), a typical high pollution episode. The multi-scale modelling system MNH-C+MOCAGE allows to simulate local and regional pollution issued from emission sources in the Marseille / Fos-Berre area as well as from remote sources (e.g. the Po Valley and / or western Mediterranean sources) and their associated transboundary pollution fluxes. Detailed dynamical, chemical and aerosol (both modal and sectional spectra with organics and inorganics) simulations generally favorably compare to surface(continental and on ships), lidar and along-flight aircraft measurements.

Effect of Mesoscale and Multiscale Modeling on the Performance of Kevlar Woven Fabric Subjected to Ballistic Impact: A Numerical Study

Science.gov (United States)

Jia, Xin; Huang, Zhengxiang; Zu, Xudong; Gu, Xiaohui; Xiao, Qiangqiang

2013-12-01

In this study, an optimal finite element model of Kevlar woven fabric that is more computational efficient compared with existing models was developed to simulate ballistic impact onto fabric. Kevlar woven fabric was modeled to yarn level architecture by using the hybrid elements analysis (HEA), which uses solid elements in modeling the yarns at the impact region and uses shell elements in modeling the yarns away from the impact region. Three HEA configurations were constructed, in which the solid element region was set as about one, two, and three times that of the projectile's diameter with impact velocities of 30 m/s (non-perforation case) and 200 m/s (perforation case) to determine the optimal ratio between the solid element region and the shell element region. To further reduce computational time and to maintain the necessary accuracy, three multiscale models were presented also. These multiscale models combine the local region with the yarn level architecture by using the HEA approach and the global region with homogenous level architecture. The effect of the varying ratios of the local and global area on the ballistic performance of fabric was discussed. The deformation and damage mechanisms of fabric were analyzed and compared among numerical models. Simulation results indicate that the multiscale model based on HEA accurately reproduces the baseline results and obviously decreases computational time.

Towards the design of 3D multiscale instructive tissue engineering constructs: Current approaches and trends.

Science.gov (United States)

Oliveira, Sara M; Reis, Rui L; Mano, João F

2015-11-01

The design of 3D constructs with adequate properties to instruct and guide cells both in vitro and in vivo is one of the major focuses of tissue engineering. Successful tissue regeneration depends on the favorable crosstalk between the supporting structure, the cells and the host tissue so that a balanced matrix production and degradation are achieved. Herein, the major occurring events and players in normal and regenerative tissue are overviewed. These have been inspiring the selection or synthesis of instructive cues to include into the 3D constructs. We further highlight the importance of a multiscale perception of the range of features that can be included on the biomimetic structures. Lastly, we focus on the current and developing tissue-engineering approaches for the preparation of such 3D constructs: top-down, bottom-up and integrative. Bottom-up and integrative approaches present a higher potential for the design of tissue engineering devices with multiscale features and higher biochemical control than top-down strategies, and are the main focus of this review. Copyright © 2015 Elsevier Inc. All rights reserved.

Deductive multiscale simulation using order parameters

Science.gov (United States)

Ortoleva, Peter J.

2017-05-16

Illustrative embodiments of systems and methods for the deductive multiscale simulation of macromolecules are disclosed. In one illustrative embodiment, a deductive multiscale simulation method may include (i) constructing a set of order parameters that model one or more structural characteristics of a macromolecule, (ii) simulating an ensemble of atomistic configurations for the macromolecule using instantaneous values of the set of order parameters, (iii) simulating thermal-average forces and diffusivities for the ensemble of atomistic configurations, and (iv) evolving the set of order parameters via Langevin dynamics using the thermal-average forces and diffusivities.

A Liver-Centric Multiscale Modeling Framework for Xenobiotics.

Directory of Open Access Journals (Sweden)

James P Sluka

Full Text Available We describe a multi-scale, liver-centric in silico modeling framework for acetaminophen pharmacology and metabolism. We focus on a computational model to characterize whole body uptake and clearance, liver transport and phase I and phase II metabolism. We do this by incorporating sub-models that span three scales; Physiologically Based Pharmacokinetic (PBPK modeling of acetaminophen uptake and distribution at the whole body level, cell and blood flow modeling at the tissue/organ level and metabolism at the sub-cellular level. We have used standard modeling modalities at each of the three scales. In particular, we have used the Systems Biology Markup Language (SBML to create both the whole-body and sub-cellular scales. Our modeling approach allows us to run the individual sub-models separately and allows us to easily exchange models at a particular scale without the need to extensively rework the sub-models at other scales. In addition, the use of SBML greatly facilitates the inclusion of biological annotations directly in the model code. The model was calibrated using human in vivo data for acetaminophen and its sulfate and glucuronate metabolites. We then carried out extensive parameter sensitivity studies including the pairwise interaction of parameters. We also simulated population variation of exposure and sensitivity to acetaminophen. Our modeling framework can be extended to the prediction of liver toxicity following acetaminophen overdose, or used as a general purpose pharmacokinetic model for xenobiotics.

Multiscale modelling of solidification microstructures, including microsegregation and microporosity, in an Al-Si-Cu alloy

International Nuclear Information System (INIS)

Lee, P.D.; Chirazi, A.; Atwood, R.C.; Wang, W.

2004-01-01

Phase transition phenomena in metallic alloys involve complex physical processes occurring over a wide range of temporal, spatial and energy scales. Multiscale modelling is a powerful methodology for understanding these complex systems. In this paper, a multiscale model of grain and pore formation is presented during solidification. At the microscale, a combined stochastic-deterministic approach based on the cellular automata method is used to solve multicomponent diffusion in a three-phase system (liquid, solid and gas), simulating the nucleation and growth of both grains and pores. The impingement of the growing pores upon the developing solid is also solved to predict the tortuous shape of the porosity, a critical factor for fatigue properties. The micromodel is coupled with a finite element method (FEM) solution of the macroscale heat transfer and fluid flow in industrial castings through the temperature and pressure fields. The result model was used to investigate the influence of local solidification time, hydrogen content, local metallostatic pressure and alloy composition upon the predicted grain structure and pore morphology. Comparison of the model predictions to both laboratory and industrial scale castings are presented

Multiscale computing in the exascale era

NARCIS (Netherlands)

Alowayyed, S.; Groen, D.; Coveney, P.V.; Hoekstra, A.G.

We expect that multiscale simulations will be one of the main high performance computing workloads in the exascale era. We propose multiscale computing patterns as a generic vehicle to realise load balanced, fault tolerant and energy aware high performance multiscale computing. Multiscale computing

Multiscale Modeling of Fracture Processes in Cementitious Materials

NARCIS (Netherlands)

Qian, Z.

2012-01-01

Concrete is a composite construction material, which is composed primarily of coarse aggregates, sands and cement paste. The fracture processes in concrete are complicated, because of the multiscale and multiphase nature of the material. In the past decades, comprehensive effort has been put to

Resolving Nuclear Reactor Lifetime Extension Questions: A Combined Multiscale Modeling and Positron Characterization approach

International Nuclear Information System (INIS)

Wirth, B; Asoka-Kumar, P; Denison, A; Glade, S; Howell, R; Marian, J; Odette, G; Sterne, P

2004-01-01

The objective of this work is to determine the chemical composition of nanometer precipitates responsible for irradiation hardening and embrittlement of reactor pressure vessel steels, which threaten to limit the operating lifetime of nuclear power plants worldwide. The scientific approach incorporates computational multiscale modeling of radiation damage and microstructural evolution in Fe-Cu-Ni-Mn alloys, and experimental characterization by positron annihilation spectroscopy and small angle neutron scattering. The modeling and experimental results are

Biology meets Physics: Reductionism and Multi-scale Modeling of Morphogenesis

DEFF Research Database (Denmark)

Green, Sara; Batterman, Robert

2017-01-01

A common reductionist assumption is that macro-scale behaviors can be described "bottom-up" if only sufficient details about lower-scale processes are available. The view that an "ideal" or "fundamental" physics would be sufficient to explain all macro-scale phenomena has been met with criticism ...... modeling in developmental biology. In such contexts, the relation between models at different scales and from different disciplines is neither reductive nor completely autonomous, but interdependent....... from philosophers of biology. Specifically, scholars have pointed to the impossibility of deducing biological explanations from physical ones, and to the irreducible nature of distinctively biological processes such as gene regulation and evolution. This paper takes a step back in asking whether bottom......-up modeling is feasible even when modeling simple physical systems across scales. By comparing examples of multi-scale modeling in physics and biology, we argue that the “tyranny of scales” problem present a challenge to reductive explanations in both physics and biology. The problem refers to the scale...

Integrating Cellular Metabolism into a Multiscale Whole-Body Model

Science.gov (United States)

Krauss, Markus; Schaller, Stephan; Borchers, Steffen; Findeisen, Rolf; Lippert, Jörg; Kuepfer, Lars

2012-01-01

Cellular metabolism continuously processes an enormous range of external compounds into endogenous metabolites and is as such a key element in human physiology. The multifaceted physiological role of the metabolic network fulfilling the catalytic conversions can only be fully understood from a whole-body perspective where the causal interplay of the metabolic states of individual cells, the surrounding tissue and the whole organism are simultaneously considered. We here present an approach relying on dynamic flux balance analysis that allows the integration of metabolic networks at the cellular scale into standardized physiologically-based pharmacokinetic models at the whole-body level. To evaluate our approach we integrated a genome-scale network reconstruction of a human hepatocyte into the liver tissue of a physiologically-based pharmacokinetic model of a human adult. The resulting multiscale model was used to investigate hyperuricemia therapy, ammonia detoxification and paracetamol-induced toxication at a systems level. The specific models simultaneously integrate multiple layers of biological organization and offer mechanistic insights into pathology and medication. The approach presented may in future support a mechanistic understanding in diagnostics and drug development. PMID:23133351

Bayesian data assimilation for stochastic multiscale models of transport in porous media.

Energy Technology Data Exchange (ETDEWEB)

Marzouk, Youssef M. (Massachusetts Institute of Technology, Cambridge, MA); van Bloemen Waanders, Bart Gustaaf (Sandia National Laboratories, Albuquerque NM); Parno, Matthew (Massachusetts Institute of Technology, Cambridge, MA); Ray, Jaideep; Lefantzi, Sophia; Salazar, Luke (Sandia National Laboratories, Albuquerque NM); McKenna, Sean Andrew (Sandia National Laboratories, Albuquerque NM); Klise, Katherine A. (Sandia National Laboratories, Albuquerque NM)

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

Collaborating for Multi-Scale Chemical Science

Energy Technology Data Exchange (ETDEWEB)

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.

Implementation of Grid-computing Framework for Simulation in Multi-scale Structural Analysis

Directory of Open Access Journals (Sweden)

Data Iranata

2010-05-01

Full Text Available A new grid-computing framework for simulation in multi-scale structural analysis is presented. Two levels of parallel processing will be involved in this framework: multiple local distributed computing environments connected by local network to form a grid-based cluster-to-cluster distributed computing environment. To successfully perform the simulation, a large-scale structural system task is decomposed into the simulations of a simplified global model and several detailed component models using various scales. These correlated multi-scale structural system tasks are distributed among clusters and connected together in a multi-level hierarchy and then coordinated over the internet. The software framework for supporting the multi-scale structural simulation approach is also presented. The program architecture design allows the integration of several multi-scale models as clients and servers under a single platform. To check its feasibility, a prototype software system has been designed and implemented to perform the proposed concept. The simulation results show that the software framework can increase the speedup performance of the structural analysis. Based on this result, the proposed grid-computing framework is suitable to perform the simulation of the multi-scale structural analysis.

Multi-scale functional and taxonomic β-diversity of the macroinvertebrate communities in a Mediterranean coastal lagoon

Directory of Open Access Journals (Sweden)

D. CABANA

2017-03-01

Full Text Available Benthic macroinvertebrate communities form the basis of the intricate lagoonal food web. Understanding their functional and taxonomic response, from a β-diversity perspective, is essential to disclose underlying patterns with potential applicability in conservation and management actions. Within the central lagoon of Messolonghi we studied the main environmental components structuring the macroinvertebrate community. We analyzed the β-taxonomic and β-functional diversity across the main habitats and seasons, over a year time frame. Our results outline habitat type and vegetation biomass as the major factors structuring the communities. We found environmental variability to have a positive correlation with functional β-diversity, however no correlation was found with taxonomic β-diversity. Across the seasons an asynchronous response of the functional and taxonomic β-diversity was identified. The taxonomic composition displayed significant heterogeneity during the driest period and the functional during the rainy season. Across the habitats the unvegetated presented higher taxonomic homogeneity and functionally heterogeneity, contrary the vegetated habitats present higher taxonomic variability and functional homogeneity. Across the seasons and habitats a pattern of functional redundancy and taxonomic replacement was identified. Besides high functional turnover versus low taxonomic turnover was documented in an anthropogenic organically enriched habitat We conclude that habitats display independent functional and taxonomic seasonal patterns, thus different processes may contribute to their variability. The framework presented here highlights the importance of studying both β-diversity components framed in a multiscale approach to better understand ecological processes and variability patterns. These results are important to understand macroinvertebrate community assembly processes and are valuable for conservation purposes.

Multi-Scale Validation of a Nanodiamond Drug Delivery System and Multi-Scale Engineering Education

Science.gov (United States)

Schwalbe, Michelle Kristin

2010-01-01

This dissertation has two primary concerns: (i) evaluating the uncertainty and prediction capabilities of a nanodiamond drug delivery model using Bayesian calibration and bias correction, and (ii) determining conceptual difficulties of multi-scale analysis from an engineering education perspective. A Bayesian uncertainty quantification scheme…

MULTISCALE THERMOHYDROLOGIC MODEL

International Nuclear Information System (INIS)

T. Buscheck

2005-01-01

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

MULTISCALE THERMOHYDROLOGIC MODEL

Energy Technology Data Exchange (ETDEWEB)

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

Key Factors Influencing the Energy Absorption of Dual-Phase Steels: Multiscale Material Model Approach and Microstructural Optimization

Science.gov (United States)

Belgasam, Tarek M.; Zbib, Hussein M.

2018-06-01

The increase in use of dual-phase (DP) steel grades by vehicle manufacturers to enhance crash resistance and reduce body car weight requires the development of a clear understanding of the effect of various microstructural parameters on the energy absorption in these materials. Accordingly, DP steelmakers are interested in predicting the effect of various microscopic factors as well as optimizing microstructural properties for application in crash-relevant components of vehicle bodies. This study presents a microstructure-based approach using a multiscale material and structure model. In this approach, Digimat and LS-DYNA software were coupled and employed to provide a full micro-macro multiscale material model, which is then used to simulate tensile tests. Microstructures with varied ferrite grain sizes, martensite volume fractions, and carbon content in DP steels were studied. The impact of these microstructural features at different strain rates on energy absorption characteristics of DP steels is investigated numerically using an elasto-viscoplastic constitutive model. The model is implemented in a multiscale finite-element framework. A comprehensive statistical parametric study using response surface methodology is performed to determine the optimum microstructural features for a required tensile toughness at different strain rates. The simulation results are validated using experimental data found in the literature. The developed methodology proved to be effective for investigating the influence and interaction of key microscopic properties on the energy absorption characteristics of DP steels. Furthermore, it is shown that this method can be used to identify optimum microstructural conditions at different strain-rate conditions.

Toward the multiscale nature of stress corrosion cracking

Directory of Open Access Journals (Sweden)

Xiaolong Liu

2018-02-01

Full Text Available This article reviews the multiscale nature of stress corrosion cracking (SCC observed by high-resolution characterizations in austenite stainless steels and Ni-base superalloys in light water reactors (including boiling water reactors, pressurized water reactors, and supercritical water reactors with related opinions. A new statistical summary and comparison of observed degradation phenomena at different length scales is included. The intrinsic causes of this multiscale nature of SCC are discussed based on existing evidence and related opinions, ranging from materials theory to practical processing technologies. Questions of interest are then discussed to improve bottom-up understanding of the intrinsic causes. Last, a multiscale modeling and simulation methodology is proposed as a promising interdisciplinary solution to understand the intrinsic causes of the multiscale nature of SCC in light water reactors, based on a review of related supporting application evidence.

Multiscale modeling of high contrast brinkman equations with applications to deformable porous media

KAUST Repository

Brown, Donald

2013-06-18

Simulating porous media flows has a wide range of applications. Often, these applications involve many scales and multi-physical processes. A useful tool in the analysis of such problems in that of homogenization as an averaged description is derived circumventing the need for complicated simulation of the fine scale features. In this work, we recall recent developments of homogenization techniques in the application of flows in deformable porous media. In addition, homogenization of media with high-contrast. In particular, we recall the main ideas of the homogenization of slowly varying Stokes flow and summarize the results of [4]. We also present the ideas for extending these techniques to high-contrast deformable media [3]. These ideas are connected by the modeling of multiscale fluid-structure interaction problems. © 2013 American Society of Civil Engineers.

A multi-scale correlative investigation of ductile fracture

International Nuclear Information System (INIS)

Daly, M.; Burnett, T.L.; Pickering, E.J.; Tuck, O.C.G.; Léonard, F.; Kelley, R.; Withers, P.J.; Sherry, A.H.

2017-01-01

The use of novel multi-scale correlative methods, which involve the coordinated characterisation of matter across a range of length scales, are becoming of increasing value to materials scientists. Here, we describe for the first time how a multi-scale correlative approach can be used to investigate the nature of ductile fracture in metals. Specimens of a nuclear pressure vessel steel, SA508 Grade 3, are examined following ductile fracture using medium and high-resolution 3D X-ray computed tomography (CT) analyses, and a site-specific analysis using a dual beam plasma focused ion beam scanning electron microscope (PFIB-SEM). The methods are employed sequentially to characterise damage by void nucleation and growth in one volume of interest, allowing for the imaging of voids that ranged in size from less than 100 nm to over 100 μm. This enables the examination of voids initiated at carbide particles to be detected, as well as the large voids initiated at inclusions. We demonstrate that this multi-scale correlative approach is a powerful tool, which not only enhances our understanding of ductile failure through detailed characterisation of microstructure, but also provides quantitative information about the size, volume fractions and spatial distributions of voids that can be used to inform models of failure. It is found that the vast majority of large voids nucleated at MnS inclusions, and that the volume of a void varied according to the volume of its initiating inclusion raised to the power 3/2. The most severe voiding was concentrated within 500 μm of the fracture surface, but measurable damage was found to extend to a depth of at least 3 mm. Microvoids associated with carbides (carbide-initiated voids) were found to be concentrated around larger inclusion-initiated voids at depths of at least 400 μm. Methods for quantifying X-ray CT void data are discussed, and a procedure for using this data to calibrate parameters in the Gurson-Tvergaard Needleman (GTN

Multiscale Modeling of Wear Degradation in Cylinder Liners

KAUST Repository

Moraes, Alvaro

2014-03-20

Every mechanical system is naturally subjected to some kind of wear process that, at some point, will cause failure in the system if no monitoring or treatment process is applied. Since failures often lead to high economical costs, it is essential both to predict and to avoid them. To achieve this, a monitoring system of the wear level should be implemented to decrease the risk of failure. In this work, we take a first step into the development of a multiscale indirect inference methodology for state-dependent Markovian pure jump processes. This allows us to model the evolution of the wear level and to identify when the system reaches some critical level that triggers a maintenance response. Since the likelihood function of a discretely observed pure jump process does not have an expression that is simple enough for standard nonsampling optimization methods, we approximate this likelihood by expressions from upscaled models of the data. We use the Master Equation (ME) to assess the goodness-of-fit and to compute the distribution of the hitting time to the critical level.

Linear theory for filtering nonlinear multiscale systems with model error.

Science.gov (United States)

Berry, Tyrus; Harlim, John

2014-07-08

In this paper, we study filtering of multiscale dynamical systems with model error arising from limitations in resolving the smaller scale processes. In particular, the analysis assumes the availability of continuous-time noisy observations of all components of the slow variables. Mathematically, this paper presents new results on higher order asymptotic expansion of the first two moments of a conditional measure. In particular, we are interested in the application of filtering multiscale problems in which the conditional distribution is defined over the slow variables, given noisy observation of the slow variables alone. From the mathematical analysis, we learn that for a continuous time linear model with Gaussian noise, there exists a unique choice of parameters in a linear reduced model for the slow variables which gives the optimal filtering when only the slow variables are observed. Moreover, these parameters simultaneously give the optimal equilibrium statistical estimates of the underlying system, and as a consequence they can be estimated offline from the equilibrium statistics of the true signal. By examining a nonlinear test model, we show that the linear theory extends in this non-Gaussian, nonlinear configuration as long as we know the optimal stochastic parametrization and the correct observation model. However, when the stochastic parametrization model is inappropriate, parameters chosen for good filter performance may give poor equilibrium statistical estimates and vice versa; this finding is based on analytical and numerical results on our nonlinear test model and the two-layer Lorenz-96 model. Finally, even when the correct stochastic ansatz is given, it is imperative to estimate the parameters simultaneously and to account for the nonlinear feedback of the stochastic parameters into the reduced filter estimates. In numerical experiments on the two-layer Lorenz-96 model, we find that the parameters estimated online , as part of a filtering

Developing Flexible, Integrated Hydrologic Modeling Systems for Multiscale Analysis in the Midwest and Great Lakes Region

Science.gov (United States)

Hamlet, A. F.; Chiu, C. M.; Sharma, A.; Byun, K.; Hanson, Z.

2016-12-01

Physically based hydrologic modeling of surface and groundwater resources that can be flexibly and efficiently applied to support water resources policy/planning/management decisions at a wide range of spatial and temporal scales are greatly needed in the Midwest, where stakeholder access to such tools is currently a fundamental barrier to basic climate change assessment and adaptation efforts, and also the co-production of useful products to support detailed decision making. Based on earlier pilot studies in the Pacific Northwest Region, we are currently assembling a suite of end-to-end tools and resources to support various kinds of water resources planning and management applications across the region. One of the key aspects of these integrated tools is that the user community can access gridded products at any point along the end-to-end chain of models, looking backwards in time about 100 years (1915-2015), and forwards in time about 85 years using CMIP5 climate model projections. The integrated model is composed of historical and projected future meteorological data based on station observations and statistical and dynamically downscaled climate model output respectively. These gridded meteorological data sets serve as forcing data for the macro-scale VIC hydrologic model implemented over the Midwest at 1/16 degree resolution. High-resolution climate model (4km WRF) output provides inputs for the analyses of urban impacts, hydrologic extremes, agricultural impacts, and impacts to the Great Lakes. Groundwater recharge estimated by the surface water model provides input data for fine-scale and macro-scale groundwater models needed for specific applications. To highlight the multi-scale use of the integrated models in support of co-production of scientific information for decision making, we briefly describe three current case studies addressing different spatial scales of analysis: 1) Effects of climate change on the water balance of the Great Lakes, 2) Future

Fusing Multiscale Charts into 3D ENC Systems Based on Underwater Topography and Remote Sensing Image

Directory of Open Access Journals (Sweden)

Tao Liu

2015-01-01

Full Text Available The purpose of this study is to propose an approach to fuse multiscale charts into three-dimensional (3D electronic navigational chart (ENC systems based on underwater topography and remote sensing image. This is the first time that the fusion of multiscale standard ENCs in the 3D ENC system has been studied. First, a view-dependent visualization technology is presented for the determination of the display condition of a chart. Second, a map sheet processing method is described for dealing with the map sheet splice problem. A process order called “3D order” is designed to adapt to the characteristics of the chart. A map sheet clipping process is described to deal with the overlap between the adjacent map sheets. And our strategy for map sheet splice is proposed. Third, the rendering method for ENC objects in the 3D ENC system is introduced. Fourth, our picking-up method for ENC objects is proposed. Finally, we implement the above methods in our system: automotive intelligent chart (AIC 3D electronic chart display and information systems (ECDIS. And our method can handle the fusion problem well.

Model reduction of multiscale chemical langevin equations: a numerical case study.

Science.gov (United States)

Sotiropoulos, Vassilios; Contou-Carrere, Marie-Nathalie; Daoutidis, Prodromos; Kaznessis, Yiannis N

2009-01-01

Two very important characteristics of biological reaction networks need to be considered carefully when modeling these systems. First, models must account for the inherent probabilistic nature of systems far from the thermodynamic limit. Often, biological systems cannot be modeled with traditional continuous-deterministic models. Second, models must take into consideration the disparate spectrum of time scales observed in biological phenomena, such as slow transcription events and fast dimerization reactions. In the last decade, significant efforts have been expended on the development of stochastic chemical kinetics models to capture the dynamics of biomolecular systems, and on the development of robust multiscale algorithms, able to handle stiffness. In this paper, the focus is on the dynamics of reaction sets governed by stiff chemical Langevin equations, i.e., stiff stochastic differential equations. These are particularly challenging systems to model, requiring prohibitively small integration step sizes. We describe and illustrate the application of a semianalytical reduction framework for chemical Langevin equations that results in significant gains in computational cost.

Multi-scale Modeling of the Impact Response of a Strain Rate Sensitive High-Manganese Austenitic Steel

Directory of Open Access Journals (Sweden)

Orkun eÖnal

2014-09-01

Full Text Available A multi-scale modeling approach was applied to predict the impact response of a strain rate sensitive high-manganese austenitic steel. The roles of texture, geometry and strain rate sensitivity were successfully taken into account all at once by coupling crystal plasticity and finite element (FE analysis. Specifically, crystal plasticity was utilized to obtain the multi-axial flow rule at different strain rates based on the experimental deformation response under uniaxial tensile loading. The equivalent stress – equivalent strain response was then incorporated into the FE model for the sake of a more representative hardening rule under impact loading. The current results demonstrate that reliable predictions can be obtained by proper coupling of crystal plasticity and FE analysis even if the experimental flow rule of the material is acquired under uniaxial loading and at moderate strain rates that are significantly slower than those attained during impact loading. Furthermore, the current findings also demonstrate the need for an experiment-based multi-scale modeling approach for the sake of reliable predictions of the impact response.

Modelling an industrial anaerobic granular reactor using a multi-scale approach

DEFF Research Database (Denmark)

Feldman, Hannah; Flores Alsina, Xavier; Ramin, Pedram

2017-01-01

The objective of this paper is to show the results of an industrial project dealing with modelling of anaerobic digesters. A multi-scale mathematical approach is developed to describe reactor hydrodynamics, granule growth/distribution and microbial competition/inhibition for substrate/space within...... the biofilm. The main biochemical and physico-chemical processes in the model are based on the Anaerobic Digestion Model No 1 (ADM1) extended with the fate of phosphorus (P), sulfur (S) and ethanol (Et-OH). Wastewater dynamic conditions are reproduced and data frequency increased using the Benchmark...... simulations show the effects on the overall process performance when operational (pH) and loading (S:COD) conditions are modified. Lastly, the effect of intra-granular precipitation on the overall organic/inorganic distribution is assessed at: 1) different times; and, 2) reactor heights. Finally...

Multiscale Retinex

Directory of Open Access Journals (Sweden)

Ana Belén Petro

2014-04-01

Full Text Available While the retinex theory aimed at explaining human color perception, its derivations have led to efficient algorithms enhancing local image contrast, thus permitting among other features, to "see in the shadows". Among these derived algorithms, Multiscale Retinex is probably the most successful center-surround image filter. In this paper, we offer an analysis and implementation of Multiscale Retinex. We point out and resolve some ambiguities of the method. In particular, we show that the important color correction final step of the method can be seriously improved. This analysis permits to come up with an automatic implementation of Multiscale Retinex which is as faithful as possible to the one described in the original paper. Overall, this implementation delivers excellent results and confirms the validity of Multiscale Retinex for image color restoration and contrast enhancement. Nevertheless, while the method parameters can be fixed, we show that a crucial choice must be left to the user, depending on the lightning condition of the image: the method must either be applied to each color independently if a color balance is required, or to the luminance only if the goal is to achieve local contrast enhancement. Thus, we propose two slightly different algorithms to deal with both cases.

Multi-scale finite element modeling allows the mechanics of amphibian neurulation to be elucidated

International Nuclear Information System (INIS)

Chen Xiaoguang; Wayne Brodland, G

2008-01-01

The novel multi-scale computational approach introduced here makes possible a new means for testing hypotheses about the forces that drive specific morphogenetic movements. A 3D model based on this approach is used to investigate neurulation in the axolotl (Ambystoma mexicanum), a type of amphibian. The model is based on geometric data from 3D surface reconstructions of live embryos and from serial sections. Tissue properties are described by a system of cell-based constitutive equations, and parameters in the equations are determined from physical tests. The model includes the effects of Shroom-activated neural ridge reshaping and lamellipodium-driven convergent extension. A typical whole-embryo model consists of 10 239 elements and to run its 100 incremental time steps requires 2 days. The model shows that a normal phenotype does not result if lamellipodium forces are uniform across the width of the neural plate; but it can result if the lamellipodium forces decrease from a maximum value at the mid-sagittal plane to zero at the plate edge. Even the seemingly simple motions of neurulation are found to contain important features that would remain hidden, they were not studied using an advanced computational model. The present model operates in a setting where data are extremely sparse and an important outcome of the study is a better understanding of the role of computational models in such environments

Multi-scale finite element modeling allows the mechanics of amphibian neurulation to be elucidated

Science.gov (United States)

Chen, Xiaoguang; Brodland, G. Wayne

2008-03-01

The novel multi-scale computational approach introduced here makes possible a new means for testing hypotheses about the forces that drive specific morphogenetic movements. A 3D model based on this approach is used to investigate neurulation in the axolotl (Ambystoma mexicanum), a type of amphibian. The model is based on geometric data from 3D surface reconstructions of live embryos and from serial sections. Tissue properties are described by a system of cell-based constitutive equations, and parameters in the equations are determined from physical tests. The model includes the effects of Shroom-activated neural ridge reshaping and lamellipodium-driven convergent extension. A typical whole-embryo model consists of 10 239 elements and to run its 100 incremental time steps requires 2 days. The model shows that a normal phenotype does not result if lamellipodium forces are uniform across the width of the neural plate; but it can result if the lamellipodium forces decrease from a maximum value at the mid-sagittal plane to zero at the plate edge. Even the seemingly simple motions of neurulation are found to contain important features that would remain hidden, they were not studied using an advanced computational model. The present model operates in a setting where data are extremely sparse and an important outcome of the study is a better understanding of the role of computational models in such environments.

A multiscale method for modeling high-aspect-ratio micro/nano flows

Science.gov (United States)

Lockerby, Duncan; Borg, Matthew; Reese, Jason

2012-11-01

In this paper we present a new multiscale scheme for simulating micro/nano flows of high aspect ratio in the flow direction, e.g. within long ducts, tubes, or channels, of varying section. The scheme consists of applying a simple hydrodynamic description over the entire domain, and allocating micro sub-domains in very small ``slices'' of the channel. Every micro element is a molecular dynamics simulation (or other appropriate model, e.g., a direct simulation Monte Carlo method for micro-channel gas flows) over the local height of the channel/tube. The number of micro elements as well as their streamwise position is chosen to resolve the geometrical features of the macro channel. While there is no direct communication between individual micro elements, coupling occurs via an iterative imposition of mass and momentum-flux conservation on the macro scale. The greater the streamwise scale of the geometry, the more significant is the computational speed-up when compared to a full MD simulation. We test our new multiscale method on the case of a converging/diverging nanochannel conveying a simple Lennard-Jones liquid. We validate the results from our simulations by comparing them to a full MD simulation of the same test case. Supported by EPSRC Programme Grant, EP/I011927/1.

Multiscale Modeling of Point and Line Defects in Cubic Lattices

National Research Council Canada - National Science Library

Chung, P. W; Clayton, J. D

2007-01-01

.... This multiscale theory explicitly captures heterogeneity in microscopic atomic motion in crystalline materials, attributed, for example, to the presence of various point and line lattice defects...

Information theory and stochastics for multiscale nonlinear systems

CERN Document Server

Majda, Andrew J; Grote, Marcus J

2005-01-01

This book introduces mathematicians to the fascinating emerging mathematical interplay between ideas from stochastics and information theory and important practical issues in studying complex multiscale nonlinear systems. It emphasizes the serendipity between modern applied mathematics and applications where rigorous analysis, the development of qualitative and/or asymptotic models, and numerical modeling all interact to explain complex phenomena. After a brief introduction to the emerging issues in multiscale modeling, the book has three main chapters. The first chapter is an introduction to information theory with novel applications to statistical mechanics, predictability, and Jupiter's Red Spot for geophysical flows. The second chapter discusses new mathematical issues regarding fluctuation-dissipation theorems for complex nonlinear systems including information flow, various approximations, and illustrates applications to various mathematical models. The third chapter discusses stochastic modeling of com...

Multiscale mechanistic modeling in pharmaceutical research and development.

Science.gov (United States)

Kuepfer, Lars; Lippert, Jörg; Eissing, Thomas

2012-01-01

Discontinuation of drug development projects due to lack of efficacy or adverse events is one of the main cost drivers in pharmaceutical research and development (R&D). Investments have to be written-off and contribute to the total costs of a successful drug candidate receiving marketing authorization and allowing return on invest. A vital risk for pharmaceutical innovator companies is late stage clinical failure since costs for individual clinical trials may exceed the one billion Euro threshold. To guide investment decisions and to safeguard maximum medical benefit and safety for patients recruited in clinical trials, it is therefore essential to understand the clinical consequences of all information and data generated. The complexity of the physiological and pathophysiological processes and the sheer amount of information available overcharge the mental capacity of any human being and prevent a prediction of the success in clinical development. A rigorous integration of knowledge, assumption, and experimental data into computational models promises a significant improvement of the rationalization of decision making in pharmaceutical industry. We here give an overview of the current status of modeling and simulation in pharmaceutical R&D and outline the perspectives of more recent developments in mechanistic modeling. Specific modeling approaches for different biological scales ranging from intracellular processes to whole organism physiology are introduced and an example for integrative multiscale modeling of therapeutic efficiency in clinical oncology trials is showcased.

Generalized Multiscale Finite-Element Method (GMsFEM) for elastic wave propagation in heterogeneous, anisotropic media

Energy Technology Data Exchange (ETDEWEB)

Gao, Kai, E-mail: kaigao87@gmail.com [Department of Geology and Geophysics, Texas A& M University, College Station, TX 77843 (United States); Fu, Shubin, E-mail: shubinfu89@gmail.com [Department of Mathematics, Texas A& M University, College Station, TX 77843 (United States); Gibson, Richard L., E-mail: gibson@tamu.edu [Department of Geology and Geophysics, Texas A& M University, College Station, TX 77843 (United States); Chung, Eric T., E-mail: tschung@math.cuhk.edu.hk [Department of Mathematics, The Chinese University of Hong Kong, Shatin, NT (Hong Kong); Efendiev, Yalchin, E-mail: efendiev@math.tamu.edu [Department of Mathematics, Texas A& M University, College Station, TX 77843 (United States); Numerical Porous Media SRI Center (NumPor), King Abdullah University of Science and Technology, Thuwal (Saudi Arabia)

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.

Generalized Multiscale Finite-Element Method (GMsFEM) for elastic wave propagation in heterogeneous, anisotropic media

International Nuclear Information System (INIS)

Gao, Kai; Fu, Shubin; Gibson, Richard L.; Chung, Eric T.; Efendiev, Yalchin

2015-01-01

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

High-Efficiency Multiscale Modeling of Cell Deformations in Confined Microenvironments in Microcirculation and Microfluidic Devices

Science.gov (United States)

Lu, Huijie; Peng, Zhangli

2017-11-01

Our goal is to develop a high-efficiency multiscale modeling method to predict the stress and deformation of cells during the interactions with their microenvironments in microcirculation and microfluidic devices, including red blood cells (RBCs) and circulating tumor cells (CTCs). There are more than 1 billion people in the world suffering from RBC diseases, e.g. anemia, sickle cell diseases, and malaria. The mechanical properties of RBCs are changed in these diseases due to molecular structure alternations, which is not only important for understanding the disease pathology but also provides an opportunity for diagnostics. On the other hand, the mechanical properties of cancer cells are also altered compared to healthy cells. This can lead to acquired ability to cross the narrow capillary networks and endothelial gaps, which is crucial for metastasis, the leading cause of cancer mortality. Therefore, it is important to predict the deformation and stress of RBCs and CTCs in microcirculations. We are developing a high-efficiency multiscale model of cell-fluid interaction to study these two topics.

Microstructure-based multiscale modeling of elevated temperature deformation in aluminum alloys

International Nuclear Information System (INIS)

Krajewski, Paul E.; Hector, Louis G.; Du Ningning; Bower, Allan F.

2010-01-01

A multiscale model for predicting elevated temperature deformation in Al-Mg alloys is presented. Constitutive models are generated from a theoretical methodology and used to investigate the effects of grain size on formability. Flow data are computed with a polycrystalline, microstructure-based model which accounts for grain boundary sliding, stress-induced diffusion, and dislocation creep. Favorable agreement is found between the computed flow data and elevated temperature tensile measurements. A creep constitutive model is then fit to the computed flow data and used in finite-element simulations of two simple gas pressure forming processes, where favorable results are observed. These results are fully consistent with gas pressure forming experiments, and suggest a greater role for constitutive models, derived largely from theoretical methodologies, in the design of Al alloys with enhanced elevated temperature formability. The methodology detailed herein provides a framework for incorporation of results from atomistic-scale models of dislocation creep and diffusion.

Nuclear Energy Advanced Modeling and Simulation (NEAMS) Accident Tolerant Fuels High Impact Problem: Coordinate Multiscale U₃Si₂ Modeling

Energy Technology Data Exchange (ETDEWEB)

Gamble, K. A. [Idaho National Lab. (INL), Idaho Falls, ID (United States); Hales, J. D. [Idaho National Lab. (INL), Idaho Falls, ID (United States); Miao, Y. [Argonne National Lab. (ANL), Argonne, IL (United States); Andersson, D. [Los Alamos National Lab. (LANL), Los Alamos, NM (United States); Zhang, Y. [Idaho National Lab. (INL), Idaho Falls, ID (United States)

2017-07-26

Since the events at the Fukushima-Daiichi nuclear power plant in March 2011 significant research has unfolded at national laboratories, universities and other institutions into alternative materials that have potential enhanced accident tolerance when compared to traditional \\uo~fuel zircaloy clad fuel rods. One of the potential replacement fuels is uranium silicide (\\usi) for its higher thermal conductivity and uranium density. The lower melting temperature is of potential concern during postulated accident conditions. Another disadvantage for \\usi~ is the lack of experimental data under power reactor conditions. Due to the aggressive development schedule for inserting some of the potential materials into lead test assemblies or rods by 2022~\\cite{bragg-sitton_2014} multiscale multiphysics modeling approaches have been used to provide insight into these materials. \\\\ \

Multi-scale modeling of urban air pollution: development and application of a Street-in-Grid model (v1.0) by coupling MUNICH (v1.0) and Polair3D (v1.8.1)

Science.gov (United States)

Kim, Youngseob; Wu, You; Seigneur, Christian; Roustan, Yelva

2018-02-01

A new multi-scale model of urban air pollution is presented. This model combines a chemistry-transport model (CTM) that includes a comprehensive treatment of atmospheric chemistry and transport on spatial scales down to 1 km and a street-network model that describes the atmospheric concentrations of pollutants in an urban street network. The street-network model is the Model of Urban Network of Intersecting Canyons and Highways (MUNICH), which consists of two main components: a street-canyon component and a street-intersection component. MUNICH is coupled to the Polair3D CTM of the Polyphemus air quality modeling platform to constitute the Street-in-Grid (SinG) model. MUNICH is used to simulate the concentrations of the chemical species in the urban canopy, which is located in the lowest layer of Polair3D, and the simulation of pollutant concentrations above rooftops is performed with Polair3D. Interactions between MUNICH and Polair3D occur at roof level and depend on a vertical mass transfer coefficient that is a function of atmospheric turbulence. SinG is used to simulate the concentrations of nitrogen oxides (NOx) and ozone (O3) in a Paris suburb. Simulated concentrations are compared to NOx concentrations measured at two monitoring stations within a street canyon. SinG shows better performance than MUNICH for nitrogen dioxide (NO2) concentrations. However, both SinG and MUNICH underestimate NOx. For the case study considered, the model performance for NOx concentrations is not sensitive to using a complex chemistry model in MUNICH and the Leighton NO-NO2-O3 set of reactions is sufficient.

Multi-scale modeling of urban air pollution: development and application of a Street-in-Grid model (v1.0 by coupling MUNICH (v1.0 and Polair3D (v1.8.1

Directory of Open Access Journals (Sweden)

Y. Kim

2018-02-01

Full Text Available A new multi-scale model of urban air pollution is presented. This model combines a chemistry–transport model (CTM that includes a comprehensive treatment of atmospheric chemistry and transport on spatial scales down to 1 km and a street-network model that describes the atmospheric concentrations of pollutants in an urban street network. The street-network model is the Model of Urban Network of Intersecting Canyons and Highways (MUNICH, which consists of two main components: a street-canyon component and a street-intersection component. MUNICH is coupled to the Polair3D CTM of the Polyphemus air quality modeling platform to constitute the Street-in-Grid (SinG model. MUNICH is used to simulate the concentrations of the chemical species in the urban canopy, which is located in the lowest layer of Polair3D, and the simulation of pollutant concentrations above rooftops is performed with Polair3D. Interactions between MUNICH and Polair3D occur at roof level and depend on a vertical mass transfer coefficient that is a function of atmospheric turbulence. SinG is used to simulate the concentrations of nitrogen oxides (NOx and ozone (O3 in a Paris suburb. Simulated concentrations are compared to NOx concentrations measured at two monitoring stations within a street canyon. SinG shows better performance than MUNICH for nitrogen dioxide (NO2 concentrations. However, both SinG and MUNICH underestimate NOx. For the case study considered, the model performance for NOx concentrations is not sensitive to using a complex chemistry model in MUNICH and the Leighton NO–NO2–O3 set of reactions is sufficient.

Multi-scale modeling for prediction of distributed cellular properties in response to substrate spatial gradients in a continuously run microreactor

DEFF Research Database (Denmark)

Lencastre Fernandes, Rita; Krühne, Ulrich; Nopens, Ingmar

2012-01-01

microbioreactor is simulated. A multiscale model consisting of the coupling of a population balance model, a kinetic model and a flow model was developed in order to predict simultaneously local concentrations of substrate (glucose), product (ethanol) and biomass, as well as the local cell size distributions....

A generalized multiscale finite element method for elastic wave propagation in fractured media

KAUST Repository

Chung, Eric T.

2016-02-26

In this paper, we consider elastic wave propagation in fractured media applying a linear-slip model to represent the effects of fractures on the wavefield. Fractured media, typically, are highly heterogeneous due to multiple length scales. Direct numerical simulations for wave propagation in highly heterogeneous fractured media can be computationally expensive and require some type of model reduction. We develop a multiscale model reduction technique that captures the complex nature of the media (heterogeneities and fractures) in the coarse scale system. The proposed method is based on the generalized multiscale finite element method, where the multiscale basis functions are constructed to capture the fine-scale information of the heterogeneous, fractured media and effectively reduce the degrees of freedom. These multiscale basis functions are coupled via the interior penalty discontinuous Galerkin method, which provides a block-diagonal mass matrix. The latter is needed for fast computation in an explicit time discretization, which is used in our simulations. Numerical results are presented to show the performance of the presented multiscale method for fractured media. We consider several cases where fractured media contain fractures of multiple lengths. Our numerical results show that the proposed reduced-order models can provide accurate approximations for the fine-scale solution.

A generalized multiscale finite element method for elastic wave propagation in fractured media

KAUST Repository

Chung, Eric T.; Efendiev, Yalchin R.; Gibson, Richard L.; Vasilyeva, Maria

2016-01-01

In this paper, we consider elastic wave propagation in fractured media applying a linear-slip model to represent the effects of fractures on the wavefield. Fractured media, typically, are highly heterogeneous due to multiple length scales. Direct numerical simulations for wave propagation in highly heterogeneous fractured media can be computationally expensive and require some type of model reduction. We develop a multiscale model reduction technique that captures the complex nature of the media (heterogeneities and fractures) in the coarse scale system. The proposed method is based on the generalized multiscale finite element method, where the multiscale basis functions are constructed to capture the fine-scale information of the heterogeneous, fractured media and effectively reduce the degrees of freedom. These multiscale basis functions are coupled via the interior penalty discontinuous Galerkin method, which provides a block-diagonal mass matrix. The latter is needed for fast computation in an explicit time discretization, which is used in our simulations. Numerical results are presented to show the performance of the presented multiscale method for fractured media. We consider several cases where fractured media contain fractures of multiple lengths. Our numerical results show that the proposed reduced-order models can provide accurate approximations for the fine-scale solution.

A multiphysics and multiscale model for low frequency electromagnetic direct-chill casting

International Nuclear Information System (INIS)

Košnik, N; Guštin, A Z; Mavrič, B; Šarler, B

2016-01-01

Simulation and control of macrosegregation, deformation and grain size in low frequency electromagnetic (EM) direct-chill casting (LFEMC) is important for downstream processing. Respectively, a multiphysics and multiscale model is developed for solution of Lorentz force, temperature, velocity, concentration, deformation and grain structure of LFEMC processed aluminum alloys, with focus on axisymmetric billets. The mixture equations with lever rule, linearized phase diagram, and stationary thermoelastic solid phase are assumed, together with EM induction equation for the field imposed by the coil. Explicit diffuse approximate meshless solution procedure [1] is used for solving the EM field, and the explicit local radial basis function collocation method [2] is used for solving the coupled transport phenomena and thermomechanics fields. Pressure-velocity coupling is performed by the fractional step method [3]. The point automata method with modified KGT model is used to estimate the grain structure [4] in a post-processing mode. Thermal, mechanical, EM and grain structure outcomes of the model are demonstrated. A systematic study of the complicated influences of the process parameters can be investigated by the model, including intensity and frequency of the electromagnetic field. The meshless solution framework, with the implemented simplest physical models, will be further extended by including more sophisticated microsegregation and grain structure models, as well as a more realistic solid and solid-liquid phase rheology. (paper)

Fully automated segmentation of oncological PET volumes using a combined multiscale and statistical model

International Nuclear Information System (INIS)

Montgomery, David W. G.; Amira, Abbes; Zaidi, Habib

2007-01-01

The widespread application of positron emission tomography (PET) in clinical oncology has driven this imaging technology into a number of new research and clinical arenas. Increasing numbers of patient scans have led to an urgent need for efficient data handling and the development of new image analysis techniques to aid clinicians in the diagnosis of disease and planning of treatment. Automatic quantitative assessment of metabolic PET data is attractive and will certainly revolutionize the practice of functional imaging since it can lower variability across institutions and may enhance the consistency of image interpretation independent of reader experience. In this paper, a novel automated system for the segmentation of oncological PET data aiming at providing an accurate quantitative analysis tool is proposed. The initial step involves expectation maximization (EM)-based mixture modeling using a k-means clustering procedure, which varies voxel order for initialization. A multiscale Markov model is then used to refine this segmentation by modeling spatial correlations between neighboring image voxels. An experimental study using an anthropomorphic thorax phantom was conducted for quantitative evaluation of the performance of the proposed segmentation algorithm. The comparison of actual tumor volumes to the volumes calculated using different segmentation methodologies including standard k-means, spatial domain Markov Random Field Model (MRFM), and the new multiscale MRFM proposed in this paper showed that the latter dramatically reduces the relative error to less than 8% for small lesions (7 mm radii) and less than 3.5% for larger lesions (9 mm radii). The analysis of the resulting segmentations of clinical oncologic PET data seems to confirm that this methodology shows promise and can successfully segment patient lesions. For problematic images, this technique enables the identification of tumors situated very close to nearby high normal physiologic uptake. The

Integrating cellular metabolism into a multiscale whole-body model.

Directory of Open Access Journals (Sweden)

Markus Krauss

Full Text Available Cellular metabolism continuously processes an enormous range of external compounds into endogenous metabolites and is as such a key element in human physiology. The multifaceted physiological role of the metabolic network fulfilling the catalytic conversions can only be fully understood from a whole-body perspective where the causal interplay of the metabolic states of individual cells, the surrounding tissue and the whole organism are simultaneously considered. We here present an approach relying on dynamic flux balance analysis that allows the integration of metabolic networks at the cellular scale into standardized physiologically-based pharmacokinetic models at the whole-body level. To evaluate our approach we integrated a genome-scale network reconstruction of a human hepatocyte into the liver tissue of a physiologically-based pharmacokinetic model of a human adult. The resulting multiscale model was used to investigate hyperuricemia therapy, ammonia detoxification and paracetamol-induced toxication at a systems level. The specific models simultaneously integrate multiple layers of biological organization and offer mechanistic insights into pathology and medication. The approach presented may in future support a mechanistic understanding in diagnostics and drug development.

Multiscale Modeling of Carbon Fiber Reinforced Polymer (CFRP) for Integrated Computational Materials Engineering Process

Energy Technology Data Exchange (ETDEWEB)

Gao, Jiaying; Liang, Biao; Zhang, Weizhao; Liu, Zeliang; Cheng, Puikei; Bostanabad, Ramin; Cao, Jian; Chen, Wei; Liu, Wing Kam; Su, Xuming; Zeng, Danielle; Zhao, John

2017-10-23

In this work, a multiscale modeling framework for CFRP is introduced to study hierarchical structure of CFRP. Four distinct scales are defined: nanoscale, microscale, mesoscale, and macroscale. Information at lower scales can be passed to higher scale, which is beneficial for studying effect of constituents on macroscale part’s mechanical property. This bottom-up modeling approach enables better understanding of CFRP from finest details. Current study focuses on microscale and mesoscale. Representative volume element is used at microscale and mesoscale to model material’s properties. At microscale, unidirection CFRP (UD) RVE is used to study properties of UD. The UD RVE can be modeled with different volumetric fraction to encounter non-uniform fiber distribution in CFRP part. Such consideration is important in modeling uncertainties at microscale level. Currently, we identified volumetric fraction as the only uncertainty parameters in UD RVE. To measure effective material properties of UD RVE, periodic boundary conditions (PBC) are applied to UD RVE to ensure convergence of obtained properties. Properties of UD is directly used at mesoscale woven RVE modeling, where each yarn is assumed to have same properties as UD. Within woven RVE, there can be many potential uncertainties parameters to consider for a physical modeling of CFRP. Currently, we will consider fiber misalignment within yarn and angle between wrap and weft yarns. PBC is applied to woven RVE to calculate its effective material properties. The effect of uncertainties are investigated quantitatively by Gaussian process. Preliminary results of UD and Woven study are analyzed for efficacy of the RVE modeling. This work is considered as the foundation for future multiscale modeling framework development for ICME project.

Multi-scale Modeling of Plasticity in Tantalum.

Energy Technology Data Exchange (ETDEWEB)

Lim, Hojun [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Battaile, Corbett Chandler. [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Carroll, Jay [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Buchheit, Thomas E. [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Boyce, Brad [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Weinberger, Christopher [Drexel Univ., Philadelphia, PA (United States)

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

A Multiscale Time-Splitting Discrete Fracture Model of Nanoparticles Transport in Fractured Porous Media

KAUST Repository

El-Amin, Mohamed F.; Kou, Jisheng; Sun, Shuyu

2017-01-01

Recently, applications of nanoparticles have been considered in many branches of petroleum engineering, especially, enhanced oil recovery. The current paper is devoted to investigate the problem of nanoparticles transport in fractured porous media, numerically. We employed the discrete-fracture model (DFM) to represent the flow and transport in the fractured formations. The system of the governing equations consists of the mass conservation law, Darcy's law, nanoparticles concentration in water, deposited nanoparticles concentration on the pore-wall, and entrapped nanoparticles concentration in the pore-throat. The variation of porosity and permeability due to the nanoparticles deposition/entrapment on/in the pores is also considered. We employ the multiscale time-splitting strategy to control different time-step sizes for different physics, such as pressure and concentration. The cell-centered finite difference (CCFD) method is used for the spatial discretization. Numerical examples are provided to demonstrate the efficiency of the proposed multiscale time splitting approach.

A Multiscale Time-Splitting Discrete Fracture Model of Nanoparticles Transport in Fractured Porous Media

KAUST Repository

El-Amin, Mohamed F.

2017-06-06

Recently, applications of nanoparticles have been considered in many branches of petroleum engineering, especially, enhanced oil recovery. The current paper is devoted to investigate the problem of nanoparticles transport in fractured porous media, numerically. We employed the discrete-fracture model (DFM) to represent the flow and transport in the fractured formations. The system of the governing equations consists of the mass conservation law, Darcy\\'s law, nanoparticles concentration in water, deposited nanoparticles concentration on the pore-wall, and entrapped nanoparticles concentration in the pore-throat. The variation of porosity and permeability due to the nanoparticles deposition/entrapment on/in the pores is also considered. We employ the multiscale time-splitting strategy to control different time-step sizes for different physics, such as pressure and concentration. The cell-centered finite difference (CCFD) method is used for the spatial discretization. Numerical examples are provided to demonstrate the efficiency of the proposed multiscale time splitting approach.

Multiscale Thermohydrologic Model Supporting the License Application for the Yucca Mountain Repository

International Nuclear Information System (INIS)

T.A. Buscheck; Y. Sun; Y. Hao

2006-01-01

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

Integrated Multiscale Latent Variable Regression and Application to Distillation Columns

Directory of Open Access Journals (Sweden)

Muddu Madakyaru

2013-01-01

Full Text Available Proper control of distillation columns requires estimating some key variables that are challenging to measure online (such as compositions, which are usually estimated using inferential models. Commonly used inferential models include latent variable regression (LVR techniques, such as principal component regression (PCR, partial least squares (PLS, and regularized canonical correlation analysis (RCCA. Unfortunately, measured practical data are usually contaminated with errors, which degrade the prediction abilities of inferential models. Therefore, noisy measurements need to be filtered to enhance the prediction accuracy of these models. Multiscale filtering has been shown to be a powerful feature extraction tool. In this work, the advantages of multiscale filtering are utilized to enhance the prediction accuracy of LVR models by developing an integrated multiscale LVR (IMSLVR modeling algorithm that integrates modeling and feature extraction. The idea behind the IMSLVR modeling algorithm is to filter the process data at different decomposition levels, model the filtered data from each level, and then select the LVR model that optimizes a model selection criterion. The performance of the developed IMSLVR algorithm is illustrated using three examples, one using synthetic data, one using simulated distillation column data, and one using experimental packed bed distillation column data. All examples clearly demonstrate the effectiveness of the IMSLVR algorithm over the conventional methods.

Multiscale modelling for better hygrothermal prediction of porous building materials

Directory of Open Access Journals (Sweden)

Belarbi Rafik

2018-01-01

Full Text Available The aim of this work is to understand the influence of the microstructuralgeometric parameters of porous building materials on the mechanisms of coupled heat, air and moisture transfers, in order to predict behavior of the building to control and improve it in its durability. For this a multi-scale approach is implemented. It consists of mastering the dominant physical phenomena and their interactions on the microscopic scale. Followed by a dual-scale modelling, microscopic-macroscopic, of coupled heat, air and moisture transfers that takes into account the intrinsic properties and microstructural topology of the material using X-ray tomography combined with the correlation of 3D images were undertaken. In fact, the hygromorphicbehavior under hydric solicitations was considered. In this context, a model of coupled heat, air and moisture transfer in porous building materials was developed using the periodic homogenization technique. These informations were subsequently implemented in a dynamic computation simulation that model the hygrothermalbehaviourof material at the scale of the envelopes and indoor air quality of building. Results reveals that is essential to consider the local behaviors of materials, but also to be able to measure and quantify the evolution of its properties on a macroscopic scale from the youngest age of the material. In addition, comparisons between experimental and numerical temperature and relative humidity profilesin multilayers wall and in building envelopes were undertaken. Good agreements were observed.

Distributed multiscale computing with MUSCLE 2, the Multiscale Coupling Library and Environment

NARCIS (Netherlands)

Borgdorff, J.; Mamonski, M.; Bosak, B.; Kurowski, K.; Ben Belgacem, M.; Chopard, B.; Groen, D.; Coveney, P.V.; Hoekstra, A.G.

2014-01-01

We present the Multiscale Coupling Library and Environment: MUSCLE 2. This multiscale component-based execution environment has a simple to use Java, C++, C, Python and Fortran API, compatible with MPI, OpenMP and threading codes. We demonstrate its local and distributed computing capabilities and